请高手指点

zhanfeng

求fh3函数在[0.001 0.01]和[0.003 0.1]之间的最小值，估计在1至5之间。

function f=fh3(n) " H3 k3 {* W5 Bg2=2^0.5-1; mu0=4*pi*10^(-7); ' {. _/ p: Y6 B7 Xgaf=7.8*10^3; gacu=8.9*10^3; Y8 q- U1 j& x, ngat=4.5*10^3; gaa=2.7*10^3; + P7 \4 h& S; a3 ~8 D e' Np=2.7*10^4; , J$ E8 N j8 E5 C- [3 zq=3*10^(-5); rou=0.7; : G! a' r( I8 U0 T+ Bkb=0.9; %j0=8*10^6; 3 z8 r6 N$ I# B0 ?6 t; [& I$ l8 Ij0=2*10^6; kd=3; , J& U6 A" e6 i p1 `; d* q% ^" Q9 O, S$ fbm=2; ht=2*10^(-3); , l/ B5 [# A% \& |0 Gmu1=10^8/1.75; ; Q, Q3 B; W. U$ N; e. p( \1 i1 nkt1=2*pi*gat*ht; 6 ?7 H9 H& D- `4 V; ikt2=2*pi*gaa*ht; kr=1.02; ku=(mu0/2/bm)^0.5; kx=mu0^0.5*kb*kd^2*j0; kp=(2*q/pi/p)^0.5; kg=mu0^0.5*kb*kd*j0*g2; ( }& }$ o% g# e, q; Xrx=ku*kx-kp*kg; 5 U" m6 N, C; |% o: ]ra=ku*kg; # V: X" k/ C6 R3 w& ?! lrc=kp*kx; kj=mu0*kb*kd*j0; + n: e( x; U+ ikj=kj/kx; h0=3*ht; $ Z; n5 u8 K1 L \7 C* N) ^h1=2*(1+1/rou)*q/pi; ! U* ^2 P3 k7 q4 E! ^h2=2*(kd-1); h3=1/bm*(mu0*p*q/(2*pi))^0.5; h4=1/kj*(mu0*p*q/(2*pi))^0.5; % Y: H% \4 y. ], `) q: xhg=2*pi*kr*ht; 3 G! x$ a# Q; Q& B$ b+ E; ~w1=gaf/bm*(pi*p*q*mu0/2)^0.5; w1p=8*g2*kd*w1; w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; 4 ~2 |7 Q# Z$ d/ h7 R" p9 v0 w# Z# i3 R& Ew3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5; w3pp=hg*h3; w3=w3p+w3pp; ( X8 h# K7 L3 ?: @w4=2*kd*gacu/kj*(2*pi*p*q*mu0); 1 L; N, `6 v9 q) \w5=5*pi*ht*gat; 9 h- j5 f& `3 e. ^. ?w6=hg*h0*gat; 5 i7 S. c d! O- a4 |! p( qw7=hg*h1*gat; 1 f% d; |- b7 ^w8=hg*h0*gat; & O8 J7 n' A/ ^1 ^( m! ew9=hg*h4*gat; f=w1p*r*x^0.5+w2*x^1.5/(g2*r-kd*x)+w3*r*x^(-0.5)+w4*x^0.5+w5*r^2+w6*r+w7/r+w8*r*x+w9*r*x^(-0.5)/(g2*r-kd*x); : F% x: g5 s+ B1 O- u* k* R

zhanfeng

[0.001 0.01],[0.003 0.1]为向量

我算了很多次，没有理想的结果，

fmins,fminsearch,fmincon我都用过，或许是没有用恰当。

请高手指点！！！

zhanfeng

kj=kj/kx;是多余的，忘记删去了

zhanfeng

我也是应人之邀啊！

邱先谢谢各位了！

zhanfeng

高手请进！

谢谢！！！

ilikenba

我没有太看懂你的目标函数，估计应该是一个非线性的约束最优化问题吧！这种问题最好先研究一下目标函数的性质！如果是光滑凸的函数，那样的话先用L乘子化成无约束的然后再用牛顿法求解就行了！如果函数十分复杂，那你可能得到了局部最优解，这是需要作适当的扰动才能跳出！

zhanfeng

我试试看，目标函数是f。

谢谢！

Allos

大不了GA

zhanfeng

哦

zhanfeng

目标函数有开平方项，且在分母上，相当复杂。

并且目标函数的最大特点是对R的一阶偏导不连续。

敬请高手指点！