请高手指点
<P>求fh3函数在和之间的最小值，估计在1至5之间。</P><P>function f=fh3(n)
g2=2^0.5-1;
mu0=4*pi*10^(-7);
gaf=7.8*10^3;
gacu=8.9*10^3;
gat=4.5*10^3;
gaa=2.7*10^3;
p=2.7*10^4;
q=3*10^(-5);
rou=0.7;
kb=0.9;
%j0=8*10^6;
j0=2*10^6;
kd=3;
bm=2;
ht=2*10^(-3);
mu1=10^8/1.75;
kt1=2*pi*gat*ht;
kt2=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;
rx=ku*kx-kp*kg;
ra=ku*kg;
rc=kp*kx;
kj=mu0*kb*kd*j0;
kj=kj/kx;
h0=3*ht;
h1=2*(1+1/rou)*q/pi;
h2=2*(kd-1);
h3=1/bm*(mu0*p*q/(2*pi))^0.5;
h4=1/kj*(mu0*p*q/(2*pi))^0.5;
hg=2*pi*kr*ht;
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;
w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;
w3pp=hg*h3;
w3=w3p+w3pp;
w4=2*kd*gacu/kj*(2*pi*p*q*mu0);
w5=5*pi*ht*gat;
w6=hg*h0*gat;
w7=hg*h1*gat;
w8=hg*h0*gat;
w9=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);
</P> <P>,为向量</P><P>我算了很多次，没有理想的结果，</P><P>fmins,fminsearch,fmincon我都用过，或许是没有用恰当。</P><P>请高手指点！！！</P> <P>kj=kj/kx;是多余的，忘记删去了</P> <P>我也是应人之邀啊！</P><P>邱先谢谢各位了！</P> <P>高手请进！</P><P>谢谢！！！</P> <P>我没有太看懂你的目标函数，估计应该是一个非线性的约束最优化问题吧！这种问题最好先研究一下目标函数的性质！如果是光滑凸的函数，那样的话先用L乘子化成无约束的然后再用牛顿法求解就行了！如果函数十分复杂，那你可能得到了局部最优解，这是需要作适当的扰动才能跳出！</P> <P>我试试看，目标函数是f。</P><P>谢谢！</P> 大不了GA <P>哦</P> <P>目标函数有开平方项，且在分母上，相当复杂。</P><P>并且目标函数的最大特点是对R的一阶偏导不连续。</P><P>敬请高手指点！</P>
页:
[1]
2