找回密码
 注册
关于网站域名变更的通知
查看: 487|回复: 3
打印 上一主题 下一主题

求解含参方程组

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-7-13 14:27 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
求解图里的含参方程组。6 N7 l* z+ b. \7 b$ N6 |: g+ ~
参数:i(虚数)、k、e(自然对数)、L1和L2。
# X" g3 b& a) G& Y7 g! C需要求出a、c1、c2、d1、d2和g含有上面参数的表达式,一共6个方程,6个未知数,可以求解,但是应该怎么才能求出含参的表达式呢?
  }# w! U: I: N& R+ B0 W/ e0 `
: O0 A; g1 l$ w1 p) Q 2 R6 L6 |, F2 _# g3 {0 l( r

该用户从未签到

2#
发表于 2020-7-13 15:16 | 只看该作者
建议先做数学推导,依次找到未知变量间的关系,直到找到一个方程Y(可能是超越的),解出该方程Y即可依次代入其他关系方程,一个一个求出其他未知量。方程Y可以使用solve():
) `  Y, D: A) w# i# rclear;clc;% N# E$ T  ]2 L, \$ S# f0 F. z6 t9 }
syms x y1 t7 u7 J( I! S  b3 L# b7 e$ ~; w
y=cos(x)+x;
  k' Q- H+ O4 t( ^6 sx0=vpa(solve(y))%cos(x)+x=0

该用户从未签到

3#
发表于 2020-7-13 15:17 | 只看该作者
参考如下代码片段,全局搜索一组方程的解:  E* _; \) }0 {0 e2 y' @7 n+ }" Q7 k/ \
y(x)=a*exp(-c*(x-b)^2)+d;$ o( K6 U; K9 ]  _
yp(a,b,c,d,x)=y(x);
& j/ p2 @$ _, y- Tdy(x)=diff(y,x); dyp(a,b,c,d,x)=dy(x);
- n( z3 ~0 L  z5 Hfunc1(a,b,c,d)=y(p1(1))-p1(2);, X9 K, Y# {! G2 w
func2(a,b,c,d)=y(p2(1))-p2(2);0 J8 }4 ?, P. X/ ]. n
func3(a,b,c,d)=dy(p1(1))-k1;
  v8 W/ v) H( p% O) E! zfunc4(a,b,c,d)=dy(p2(1))-k2;- n7 C/ ]% A" x. R  }
funcs=[func1(a,b,c,d);func2(a,b,c,d);func3(a,b,c,d);func4(a,b,c,d)];8 r  F5 {: u2 f, t/ }
vs=[a,b,c,d];
+ D0 L$ {0 y! a: ^& Y8 b# }vs0=[h,L,1.5*L/2,h];
) b5 h+ Q8 O0 v; U! `/ I6 `5 E/ BfuncHd=matlabFunction(funcs,'vars',{[a,b,c,d]});- R9 g( r$ x1 L) n! c
n=4;m=20;7 _- a3 c$ t. _% l% X6 d5 u) }
[vss,fval,exitFlag] = GlobalSolve(funcHd,n);

该用户从未签到

4#
发表于 2020-7-13 15:17 | 只看该作者
楼上都太厉害了啊
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-5 11:59 , Processed in 0.156250 second(s), 26 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表