数学建模社区-数学中国

标题: 请高手指点 [打印本页]

作者: zhanfeng 时间: 2004-12-17 21:58
标题: 请高手指点

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

function f=fh3(n)+ |$ D) U$ L+ i# A- C g2=2^0.5-1; ( v8 b( P! y( p( j& a& C4 T& Vmu0=4*pi*10^(-7); : t8 X' y! Z' F& l7 {. G Qgaf=7.8*10^3; 1 a y; u0 }% b2 I Qgacu=8.9*10^3;- G; N% D( C' ] gat=4.5*10^3;1 q/ N. D% ` ^3 N. B gaa=2.7*10^3;6 w# J6 d6 ~* W p=2.7*10^4;. b5 V; A5 p3 [8 @ q=3*10^(-5); 9 f) G+ B( M& \6 `9 c& [rou=0.7; % F- i: s D' O9 v6 a. Dkb=0.9; + v) o2 e4 X, e1 f%j0=8*10^6;9 O. }; @& ]3 N3 B7 T1 t j0=2*10^6;6 `! i x% S% v" D' a9 a kd=3;6 R0 O& N& A/ f, d! z bm=2;( ~7 W" u4 v8 o, h; D ht=2*10^(-3); % c8 s1 m3 T; N K# l; T8 Bmu1=10^8/1.75; s4 H$ v* n: ]6 Vkt1=2*pi*gat*ht; ' u- Q& J1 A3 F8 Gkt2=2*pi*gaa*ht;; n9 H& j' W: [; Z# Y: m5 [" X kr=1.02;( E$ F a- w& T ku=(mu0/2/bm)^0.5;# i. S( T5 A2 G8 G3 M: i kx=mu0^0.5*kb*kd^2*j0;- k! }: g% U L, K kp=(2*q/pi/p)^0.5; ( w2 E& K3 O' Y6 j5 t$ qkg=mu0^0.5*kb*kd*j0*g2; 6 E7 O, Y5 K; K2 S/ A8 crx=ku*kx-kp*kg; 4 j! @4 E( _2 `2 V2 l0 Kra=ku*kg;; ~3 ~0 o9 b I3 q; ^, L* v rc=kp*kx; * q/ l- I9 E: |# f" Dkj=mu0*kb*kd*j0; % m q) R' g) R0 r! ?% Ckj=kj/kx; ( f3 w6 E B% \h0=3*ht; - B4 Z: s: {/ h8 Ch1=2*(1+1/rou)*q/pi; * t: m% v1 T$ O8 qh2=2*(kd-1); 3 ]. [/ X. p5 C+ R3 D/ Gh3=1/bm*(mu0*p*q/(2*pi))^0.5;1 ^8 c7 e& V. q" v7 t* V h4=1/kj*(mu0*p*q/(2*pi))^0.5; ) D [# P; R& n( Y2 E4 l8 Uhg=2*pi*kr*ht; % t" z. } r* P1 w& b7 Q& kw1=gaf/bm*(pi*p*q*mu0/2)^0.5; ' J2 T% S4 ~8 Y" d" ~* ]* `0 Qw1p=8*g2*kd*w1; ' H( |. L& e" v+ \8 R4 `0 v5 P* Hw2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; ( N+ J& }, l/ V* c. C( j1 {w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;7 H2 j, X R: H! ~ D w3pp=hg*h3; ' p p1 N3 ?: w: Vw3=w3p+w3pp;! M, N' C$ u$ u9 ?4 X w4=2*kd*gacu/kj*(2*pi*p*q*mu0); 4 ~8 Q) b/ s! B8 L( P- ow5=5*pi*ht*gat; - s# ^# b, y) B9 X6 t9 r- zw6=hg*h0*gat; : T7 h h# w# \# U- X% gw7=hg*h1*gat; & n5 U: A8 o: z% n' k; qw8=hg*h0*gat; 6 `# J! \5 g" b' C5 ow9=hg*h4*gat;! t& Y; d3 j+ Y5 N: l t 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);, R" l. Q. i/ f, G

作者: zhanfeng 时间: 2004-12-17 22:01

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

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

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

请高手指点！！！

作者: zhanfeng 时间: 2004-12-17 22:03

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

作者: zhanfeng 时间: 2004-12-17 22:05

我也是应人之邀啊！

邱先谢谢各位了！

作者: zhanfeng 时间: 2004-12-18 20:56

高手请进！

谢谢！！！

作者: ilikenba 时间: 2004-12-19 23:34

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

作者: zhanfeng 时间: 2004-12-20 17:11

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

谢谢！

作者: Allos 时间: 2004-12-22 21:28
大不了GA

作者: zhanfeng 时间: 2004-12-25 11:48

哦

作者: zhanfeng 时间: 2004-12-28 21:35

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

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

敬请高手指点！

作者: ilikenba 时间: 2004-12-28 21:44
不会吧！如果导数不连续和不可导是差不多的性质！你看看是不是半光滑的！

[此贴子已经被作者于2005-1-1 23:03:14编辑过]

作者: sunchuren 时间: 2004-12-31 09:33

建议楼主能够把问题用解析方式表达出来，这样更便于讨论。而不是放一段程序。

作者: zhanfeng 时间: 2005-1-6 20:48
好建议！

作者: zhanfeng 时间: 2005-1-6 20:54

g2=2^0.5-1; $ E+ @, c8 S! [0 w8 w# ]% M$ X. Lmu0=4*pi*10^(-7);8 E2 {. z! h9 C; C gaf=7.8*10^3;2 k) [( [! I& p o1 g3 Z6 K+ B gacu=8.9*10^3;2 x6 \1 h& |! u: u gat=4.5*10^3; + m2 c9 V% O1 [1 o" H# A8 |5 A- N2 Lgaa=2.7*10^3; F* V8 o; H0 sp=2.7*10^4;' |/ Y6 o0 K7 @1 ^ q=3*10^(-5);; A/ e8 N/ g+ a% T$ w rou=0.7; ) ]" C$ t0 }; Q; Pkb=0.9;3 { Z9 N4 G! I- O %j0=8*10^6; " T! Z4 p& R) Ej0=2*10^6; % F7 G- Y6 s$ m+ Z$ _! k5 Q5 dkd=3;/ E( p K' W) x bm=2;% P v1 c b5 j2 \; f0 }5 e5 Y ht=2*10^(-3);8 [; k3 D4 `9 D; `1 I6 m8 M" g- y6 e mu1=10^8/1.75; % [) i% @7 @* y2 L$ tkt1=2*pi*gat*ht; / b1 m2 v+ I' [8 P* J- Kkt2=2*pi*gaa*ht;' K I+ m7 t: |; { kr=1.02;9 F) _2 _1 A3 O2 J/ \3 D( ] ku=(mu0/2/bm)^0.5;( U% \/ @) v' ]/ Y1 G kx=mu0^0.5*kb*kd^2*j0; 1 Y1 j6 p0 f/ k! g1 qkp=(2*q/pi/p)^0.5;) [* _5 m( z' p/ f0 l6 ~5 B7 p kg=mu0^0.5*kb*kd*j0*g2; ( w: T$ P/ W2 \9 P. |9 vrx=ku*kx-kp*kg; + B P1 R4 t% p7 kra=ku*kg;" r2 I$ X% {8 V% f- @6 R rc=kp*kx; 3 }; }' r/ R) R/ O, Bkj=mu0*kb*kd*j0; 5 E; O" N8 g5 |; D! @0 A) @kj=kj/kx;- w$ |& p) Q% G h0=3*ht; " z6 Z& G7 O/ o; D M% Nh1=2*(1+1/rou)*q/pi;, F0 e: ^ e! }9 l- ] h2=2*(kd-1); & d' c0 R: i5 f- uh3=1/bm*(mu0*p*q/(2*pi))^0.5;& k: k$ ^0 j+ x h4=1/kj*(mu0*p*q/(2*pi))^0.5;; u1 K1 n6 b- ^ h c hg=2*pi*kr*ht; 6 z4 }- {4 n& s# a/ i* O4 J3 ]7 aw1=gaf/bm*(pi*p*q*mu0/2)^0.5; O; A7 P6 a7 i( T+ o w1p=8*g2*kd*w1;- D- D# E- U/ K1 i. [ w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;7 ~/ [% f0 b. F+ ?* O w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5; - w p, W. y y. e3 M* ~: @7 T( pw3pp=hg*h3; * P3 |, V. H- d: D& t) _( Iw3=w3p+w3pp; % X6 [1 V+ Y% mw4=2*kd*gacu/kj*(2*pi*p*q*mu0);* n" e. s: H6 ]" l, g4 a2 L w5=5*pi*ht*gat;- k, f8 @7 e/ j* d( y6 B6 F" L w6=hg*h0*gat;+ ?) y5 k& w+ ^ w7=hg*h1*gat;! i- P& `1 _# Q% O' H9 W% `; H w8=hg*h0*gat;- \6 ~; k% n8 m! ?1 }/ _4 Z w9=hg*h4*gat;5 l" R# \5 ]5 p1 } 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;

这是后来修改后的，请高手指点。& w9 @" J! k Z3 O7 V+ ^3 S Q

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)