求fh3函数在[0.001 0.01]和[0.003 0.1]之间的最小值,估计在1至5之间。
function f=fh3(n) g2=2^0.5-1;8 U# r1 d+ k8 e' u& m$ O' u mu0=4*pi*10^(-7); gaf=7.8*10^3;) W- G$ B$ b& V( e0 T- L$ } gacu=8.9*10^3;9 e8 A- _ Z7 }4 V8 |- j% x gat=4.5*10^3;. e; d: a7 n& J gaa=2.7*10^3;0 I" J" x3 G. o& a! d p=2.7*10^4;' v, B8 j1 V8 ^ q=3*10^(-5); rou=0.7; kb=0.9; %j0=8*10^6; j0=2*10^6; kd=3;$ l5 b H% U" b; X bm=2; [8 a: Q: [9 g4 r' t/ z3 ` ht=2*10^(-3);6 R) t7 D0 b( ?+ r4 Y5 X mu1=10^8/1.75; kt1=2*pi*gat*ht; kt2=2*pi*gaa*ht;9 h/ B7 e5 A$ m9 t+ l kr=1.02; ku=(mu0/2/bm)^0.5;& G6 P/ W4 c: C1 e8 B kx=mu0^0.5*kb*kd^2*j0;& x2 O. g5 X/ ?9 t2 G kp=(2*q/pi/p)^0.5;2 k; N; |# R6 B; K8 i kg=mu0^0.5*kb*kd*j0*g2; rx=ku*kx-kp*kg;1 _6 ~( n, F0 ` ra=ku*kg; rc=kp*kx; kj=mu0*kb*kd*j0;0 Q; L0 [# t& `2 r1 F9 {9 M3 _+ { kj=kj/kx;7 C* n+ K9 n) q4 X P h0=3*ht; x7 } v h+ g; p Z0 X2 B h1=2*(1+1/rou)*q/pi;4 j" M; O1 e$ I' d4 U/ M9 r* l h2=2*(kd-1); h3=1/bm*(mu0*p*q/(2*pi))^0.5;0 a6 \+ n. E# d! w! ^" w0 D" n h4=1/kj*(mu0*p*q/(2*pi))^0.5; hg=2*pi*kr*ht;3 S9 B: |! F9 Z( J/ e w1=gaf/bm*(pi*p*q*mu0/2)^0.5;% r9 T" s- X5 o" F2 q 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;# S% W. R) p! _- j( _ w7=hg*h1*gat;6 f. l. S3 l: l: H8 k4 ]6 f w8=hg*h0*gat;1 W" ^' T8 Y% s3 p) F- D) @ 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);: a$ g1 z: g7 I/ M7 ], F: r
[0.001 0.01],[0.003 0.1]为向量
我算了很多次,没有理想的结果,
fmins,fminsearch,fmincon我都用过,或许是没有用恰当。
请高手指点!!!
kj=kj/kx;是多余的,忘记删去了
我也是应人之邀啊!
邱先谢谢各位了!
高手请进!
谢谢!!!
我没有太看懂你的目标函数,估计应该是一个非线性的约束最优化问题吧!这种问题最好先研究一下目标函数的性质!如果是光滑凸的函数,那样的话先用L乘子化成无约束的然后再用牛顿法求解就行了!如果函数十分复杂,那你可能得到了局部最优解,这是需要作适当的扰动才能跳出!
我试试看,目标函数是f。
谢谢!
哦
目标函数有开平方项,且在分母上,相当复杂。
并且目标函数的最大特点是对R的一阶偏导不连续。
敬请高手指点!
建议楼主能够把问题用解析方式表达出来,这样更便于讨论。而不是放一段程序。
g2=2^0.5-1; mu0=4*pi*10^(-7);3 A1 C3 S W4 O& R gaf=7.8*10^3;8 ~, @# _' H, o7 r: | gacu=8.9*10^3; gat=4.5*10^3; gaa=2.7*10^3;# V" g1 j- [* L p=2.7*10^4; q=3*10^(-5); rou=0.7;; c! C9 p7 C1 D% g! E: ]& S5 B8 W kb=0.9; %j0=8*10^6; j0=2*10^6; kd=3; bm=2;* _6 Q# j2 ?9 { ht=2*10^(-3); mu1=10^8/1.75;; J0 e5 J' F b" f$ s$ a& p kt1=2*pi*gat*ht;& m5 ~. c$ V+ z) [: l# \ kt2=2*pi*gaa*ht; kr=1.02;- y$ }. _: p# K; Y ku=(mu0/2/bm)^0.5;3 w2 k: |9 C7 i2 u8 o0 Y: N0 _ kx=mu0^0.5*kb*kd^2*j0; kp=(2*q/pi/p)^0.5; kg=mu0^0.5*kb*kd*j0*g2;% y% b8 N% R& ]- d9 n0 P" U rx=ku*kx-kp*kg;) r* `) A7 C# y7 K. w ~ ra=ku*kg; rc=kp*kx;, j3 o& @4 y& R( U kj=mu0*kb*kd*j0; kj=kj/kx; h0=3*ht; h1=2*(1+1/rou)*q/pi;3 [5 z1 p7 L) h7 R( V$ q h2=2*(kd-1); h3=1/bm*(mu0*p*q/(2*pi))^0.5; h4=1/kj*(mu0*p*q/(2*pi))^0.5;$ d8 }. O0 S0 G" H$ K hg=2*pi*kr*ht;( f; c7 P. @- x9 J( D& {6 |) i w1=gaf/bm*(pi*p*q*mu0/2)^0.5;: R# W" ]$ D m/ E$ i4 X& ^- N U8 j) N w1p=8*g2*kd*w1;+ T" a5 A/ U4 d- } w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; _+ |4 x( ? L" S1 h w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;+ M5 b1 x; n$ ~ w3pp=hg*h3; w3=w3p+w3pp;% [! d) G+ `6 l" d3 C- n8 F( O) } w4=2*kd*gacu/kj*(2*pi*p*q*mu0); w5=5*pi*ht*gat; w6=hg*h0*gat; w7=hg*h1*gat;! S% L+ ^' t- S4 k5 I w8=hg*h0*gat; w9=hg*h4*gat; min=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);
Subject T
h0+2*hp/rou+h2*X+2*hd+hz-0.25<=0;
4*kx*j0/mu1/kd*(pi*p*q/2/mu0)^0.5*(g2*R-kd*X)^(-0.5)-10<=0;
这是后来修改后的,请高手指点。/ H5 W3 \/ M( l
| 欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) | Powered by Discuz! X2.5 |