|
|
参考如下代码片段，全局搜索一组方程的解：
. P9 h$ P" ]' Z9 by(x)=a*exp(-c*(x-b)^2)+d;- a, m9 l+ d5 }. P/ G* v1 `( q
yp(a,b,c,d,x)=y(x);! o8 l9 W8 V9 b9 q+ W" {
dy(x)=diff(y,x); dyp(a,b,c,d,x)=dy(x);
, I3 z$ I' s% u: Y+ h, o% cfunc1(a,b,c,d)=y(p1(1))-p1(2);
8 t9 U0 ~8 K1 P# N6 Ffunc2(a,b,c,d)=y(p2(1))-p2(2);& H3 n( R' J5 ~- s# z4 \
func3(a,b,c,d)=dy(p1(1))-k1;- a, ^9 p$ W/ X2 K2 Z
func4(a,b,c,d)=dy(p2(1))-k2;& G5 T) e" o5 z0 m2 J' X+ T. w. u
funcs=[func1(a,b,c,d);func2(a,b,c,d);func3(a,b,c,d);func4(a,b,c,d)];4 s! e8 L2 B- N3 g- A- s+ |7 e. }
vs=[a,b,c,d];, j% M8 p5 u/ s
vs0=[h,L,1.5*L/2,h];, C( c' z% p4 F4 V+ b0 \& d8 s
funcHd=matlabFunction(funcs,'vars',{[a,b,c,d]});
) O$ B; D7 O z, z( A6 n* l5 Wn=4;m=20;
+ U, ~1 |1 I( A; F[vss,fval,exitFlag] = GlobalSolve(funcHd,n); |
|
/1 
发表于 2020-7-13 14:27
收藏
淘帖
支持！
反对！