|
|
参考如下代码片段,全局搜索一组方程的解:
, {- Q" c. h& Ty(x)=a*exp(-c*(x-b)^2)+d;' }! o6 G& p5 q- L: X
yp(a,b,c,d,x)=y(x);/ S( E' L8 z; c5 P* U' J
dy(x)=diff(y,x); dyp(a,b,c,d,x)=dy(x); e# M; M6 }5 n, o4 q/ m7 P! Z( { p
func1(a,b,c,d)=y(p1(1))-p1(2);
* e5 e) n2 G$ a t/ \" N: gfunc2(a,b,c,d)=y(p2(1))-p2(2);
. o) n. K( q7 Q; Mfunc3(a,b,c,d)=dy(p1(1))-k1;
/ n' Z. n6 d3 u, \+ ufunc4(a,b,c,d)=dy(p2(1))-k2;0 F: i: L! H/ K9 Y
funcs=[func1(a,b,c,d);func2(a,b,c,d);func3(a,b,c,d);func4(a,b,c,d)];2 [6 n0 P6 z( F1 ~: N
vs=[a,b,c,d];
2 Q5 b7 E- u$ j" vvs0=[h,L,1.5*L/2,h];
, l1 h7 c# o6 ]- t- LfuncHd=matlabFunction(funcs,'vars',{[a,b,c,d]});
* U5 w0 b- J8 N2 Yn=4;m=20;* C; g+ E# p7 `# V
[vss,fval,exitFlag] = GlobalSolve(funcHd,n); |
|