EDA365电子论坛网

标题: 求解含参方程组 [打印本页]
作者: domenica    时间: 2020-7-13 14:27
标题: 求解含参方程组
求解图里的含参方程组。# _0 x6 }2 v. X' f3 G9 |
参数：i(虚数）、k、e（自然对数）、L1和L2。
$ K% n1 E( W/ w8 K需要求出a、c1、c2、d1、d2和g含有上面参数的表达式，一共6个方程，6个未知数，可以求解，但是应该怎么才能求出含参的表达式呢？; `1 @5 i/ T+ p& |2 t

4 G% Z  L4 w+ z" Q& D! B: @4 I+ L% P
( j% H9 |7 ]$ w. x' ]& O) Y
作者: mytomorrow    时间: 2020-7-13 15:16
建议先做数学推导，依次找到未知变量间的关系，直到找到一个方程Y(可能是超越的)，解出该方程Y即可依次代入其他关系方程，一个一个求出其他未知量。方程Y可以使用solve():
  I* n1 a6 J8 G% e9 ^clear;clc;
( }9 \, S8 B# ksyms x y
, K- v, e8 Q* }, B* hy=cos(x)+x;1 X5 j4 {/ }/ g# w+ [0 t
x0=vpa(solve(y))%cos(x)+x=0
作者: baqiao    时间: 2020-7-13 15:17
参考如下代码片段，全局搜索一组方程的解：
8 R& f3 P. g( p9 jy(x)=a*exp(-c*(x-b)^2)+d;8 d' }. g9 G( ~6 L/ ]0 o; T6 D5 w) W
yp(a,b,c,d,x)=y(x);) {; n) A/ \6 A; ]6 |$ y
dy(x)=diff(y,x); dyp(a,b,c,d,x)=dy(x);
  i' H7 ^. u( W8 I! ^9 ]; c- ?" R2 gfunc1(a,b,c,d)=y(p1(1))-p1(2);9 B3 s0 L0 v: k% |  _1 _/ _
func2(a,b,c,d)=y(p2(1))-p2(2);
' e2 q" k9 D) z7 }3 b) Y0 ^func3(a,b,c,d)=dy(p1(1))-k1;
# |  v5 o) {! yfunc4(a,b,c,d)=dy(p2(1))-k2;
: S- r5 Q& Q! a- [6 Efuncs=[func1(a,b,c,d);func2(a,b,c,d);func3(a,b,c,d);func4(a,b,c,d)];7 ^/ E3 C1 c4 a! H9 n0 ~: O
vs=[a,b,c,d];
/ b' b: v( {$ ?* r& ]; Hvs0=[h,L,1.5*L/2,h];
, r. l8 y+ N$ k) ^/ W* HfuncHd=matlabFunction(funcs,'vars',{[a,b,c,d]});+ g) B8 h6 }$ _
n=4;m=20;9 c- P$ z; u) A' I4 Q& I; G
[vss,fval,exitFlag] = GlobalSolve(funcHd,n);
作者: NingW    时间: 2020-7-13 15:17
楼上都太厉害了啊



欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/) Powered by Discuz! X3.2