|
参考如下代码片段,全局搜索一组方程的解:
7 |: a2 d _7 }y(x)=a*exp(-c*(x-b)^2)+d;! W5 ?9 F9 t; \7 c p. N
yp(a,b,c,d,x)=y(x);+ `- n. L0 D$ R2 |6 i( `: o+ m
dy(x)=diff(y,x); dyp(a,b,c,d,x)=dy(x);" E0 i4 D2 V: u6 [
func1(a,b,c,d)=y(p1(1))-p1(2);6 e r1 |* Z9 ]1 N+ O% q
func2(a,b,c,d)=y(p2(1))-p2(2);
# n3 k% ]4 v3 b3 [8 t1 Cfunc3(a,b,c,d)=dy(p1(1))-k1;) O+ J2 k: E$ u- S
func4(a,b,c,d)=dy(p2(1))-k2;" b& s* L$ h$ E" @5 `
funcs=[func1(a,b,c,d);func2(a,b,c,d);func3(a,b,c,d);func4(a,b,c,d)];3 Z# ]9 a* ?6 H- Y6 E% R
vs=[a,b,c,d];
# V" y2 v/ Q. L+ svs0=[h,L,1.5*L/2,h];6 {( `3 X7 f$ g1 e
funcHd=matlabFunction(funcs,'vars',{[a,b,c,d]});
1 s: I7 C: X( O! F! H4 B: pn=4;m=20;
; i. h0 x: \/ g+ w. W8 c- k[vss,fval,exitFlag] = GlobalSolve(funcHd,n); |
|