求解含参方程组

domenica

求解图里的含参方程组。
参数：i(虚数）、k、e（自然对数）、L1和L2。
! s' v1 n. ^' m  G/ E: y需要求出a、c1、c2、d1、d2和g含有上面参数的表达式，一共6个方程，6个未知数，可以求解，但是应该怎么才能求出含参的表达式呢？
) s( a8 p2 N/ u6 ]

mytomorrow

建议先做数学推导，依次找到未知变量间的关系，直到找到一个方程Y(可能是超越的)，解出该方程Y即可依次代入其他关系方程，一个一个求出其他未知量。方程Y可以使用solve():
1 i, k- Q0 w+ v: ]4 ~clear;clc;
7 p% E0 |0 x/ K4 vsyms x y
y=cos(x)+x;
! T: O. G! Z* N: A; g3 y4 X; Vx0=vpa(solve(y))%cos(x)+x=0

baqiao

参考如下代码片段，全局搜索一组方程的解：
y(x)=a*exp(-c*(x-b)^2)+d;
5 C1 ~" ~+ O! k; s6 x: [2 V! oyp(a,b,c,d,x)=y(x);
6 C9 D, ^3 y3 t: R, G, y  |dy(x)=diff(y,x); dyp(a,b,c,d,x)=dy(x);
func1(a,b,c,d)=y(p1(1))-p1(2);
; R1 e6 \/ i, D. z3 bfunc2(a,b,c,d)=y(p2(1))-p2(2);
, Z+ |' R- D/ u) b7 @/ t2 a* ufunc3(a,b,c,d)=dy(p1(1))-k1;
func4(a,b,c,d)=dy(p2(1))-k2;
* |; c& J4 [1 y6 Vfuncs=[func1(a,b,c,d);func2(a,b,c,d);func3(a,b,c,d);func4(a,b,c,d)];
vs=[a,b,c,d];
! z8 o" V; N) I. @% uvs0=[h,L,1.5*L/2,h];
0 D4 ~2 x6 x- }. pfuncHd=matlabFunction(funcs,'vars',{[a,b,c,d]});
n=4;m=20;
4 F3 |' F( Z6 V1 o) r[vss,fval,exitFlag] = GlobalSolve(funcHd,n);

NingW

楼上都太厉害了啊