数学建模社区-数学中国

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

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

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

function f=fh3(n) 2 K4 g. m* J3 ~! ]g2=2^0.5-1; : ^8 V# p# ]7 U4 {! Xmu0=4*pi*10^(-7); _- B5 e: C: J( p5 Z6 T- S! A! Wgaf=7.8*10^3;0 J; K6 H5 v6 i9 w+ ] gacu=8.9*10^3;- {) U) b3 x5 a1 b X) z/ F gat=4.5*10^3; * }, d3 j/ `1 @: e- ]; `( `* Dgaa=2.7*10^3;+ r1 {( }8 @# `' O5 x1 i' D: \ p=2.7*10^4;$ z5 f4 O, @% d) ?' @, ^3 r q=3*10^(-5);( |) Q# m( `! O8 v+ c7 p4 z* | rou=0.7; 6 v9 A4 Y& G( G* c0 ]8 Z' c9 ?kb=0.9;' ?4 Q* H* q% P; V' G3 W* ^ %j0=8*10^6; 4 [. B. R% X+ p9 [3 L6 fj0=2*10^6; * R+ a, B+ \( f$ skd=3; 3 J( ^4 E7 U: H. B* a4 e7 mbm=2;; \0 o. d1 ? l% r% h9 Z* J ht=2*10^(-3); 3 o* i7 w# w- A) z2 g6 ]. O$ gmu1=10^8/1.75;( U* _$ I: o6 K) E+ C% F# I2 h kt1=2*pi*gat*ht; ; m9 j" v! {. okt2=2*pi*gaa*ht; & h5 ]8 o2 O; h( Ckr=1.02; ( M5 o8 p. |7 Y# zku=(mu0/2/bm)^0.5; 3 d8 J: N7 @; }, S, W" b6 qkx=mu0^0.5*kb*kd^2*j0;( [+ h! `: X. t) ]/ W1 n kp=(2*q/pi/p)^0.5; % b4 j( i2 b4 m7 Q a9 ukg=mu0^0.5*kb*kd*j0*g2; 8 C2 T7 R& s- w* I$ Wrx=ku*kx-kp*kg; & q9 ]. B, m# X( p: \ra=ku*kg; 8 }, `4 r$ J4 @( Wrc=kp*kx; # R( }; g( ?" d( Xkj=mu0*kb*kd*j0;1 b- A/ s! V/ [+ Q% M' T% V kj=kj/kx; ; z0 B5 T, [6 y' h8 J' `8 zh0=3*ht;* M# |* y7 @4 a" x h1=2*(1+1/rou)*q/pi;5 M: o! l5 o2 r' c8 ^ h2=2*(kd-1);6 h7 q4 Z2 E$ R4 \6 x, ^" o h3=1/bm*(mu0*p*q/(2*pi))^0.5; & I% W# J w, jh4=1/kj*(mu0*p*q/(2*pi))^0.5;1 l+ D, C5 o' f7 B& d; D hg=2*pi*kr*ht; w& c" k9 r# S w1=gaf/bm*(pi*p*q*mu0/2)^0.5;( j1 V' Y' d- n$ v6 p q w1p=8*g2*kd*w1; ! j% A7 H) Y9 d6 p8 M( R- mw2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; & X' h f i4 a0 i+ w6 }4 ~w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;7 [& z/ p+ c2 n3 u w3pp=hg*h3; R( m4 p9 q+ D3 _( R4 y w3=w3p+w3pp;8 Z# t: T$ C% a& d w4=2*kd*gacu/kj*(2*pi*p*q*mu0); y! r O' j$ r' Fw5=5*pi*ht*gat;" B7 Z# \4 K% p% H8 a$ x w6=hg*h0*gat;6 }6 D1 d3 U, p' l& I" \ w7=hg*h1*gat;8 z4 ^7 _+ r) F2 J. W4 X0 e w8=hg*h0*gat;! D3 V+ u( V; x. [1 ^4 C4 c w9=hg*h4*gat; ' q! R. W- D' z+ e- pf=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);" c3 A6 |/ E _0 j6 z% t) f2 q9 T- v9 b

作者: 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;: g: f# T4 P) u/ l mu0=4*pi*10^(-7);6 R. S2 O1 g0 l gaf=7.8*10^3;! A- y* v, i' _" ^ C gacu=8.9*10^3; ) I+ ~. p$ ]$ Y0 T, K7 ^gat=4.5*10^3; 4 R1 @" I5 _8 O: B( `gaa=2.7*10^3;8 D R7 R. @9 O1 t3 e# S3 z$ [6 E p=2.7*10^4; 9 s+ f/ R# k4 R d) `" qq=3*10^(-5); - E' s% k' T; Trou=0.7;4 \6 h! D* K- u5 t kb=0.9; ! Q- r" L+ [& N2 Z! W, V m: w0 p3 K& Y%j0=8*10^6;9 Q' M( U% v" Q% L6 f1 E j0=2*10^6;; K s& u- ~3 w) x/ D. g) I kd=3;6 V$ w* a# ^$ O9 b2 I" _ bm=2;$ d& n n: V4 N2 { ht=2*10^(-3);) q+ O6 Y7 I; B mu1=10^8/1.75; - [4 B+ Z9 z. P& p( r; ekt1=2*pi*gat*ht;" m2 l8 [# v8 H% b8 H" a kt2=2*pi*gaa*ht; " Q5 S; U2 p2 V9 Y4 A5 `* Gkr=1.02;" Q! W7 p3 G3 a/ c3 S' K ku=(mu0/2/bm)^0.5;( O9 X0 l. S5 h* y+ e7 A kx=mu0^0.5*kb*kd^2*j0;7 n7 p' a% I Y. |6 Z m! n kp=(2*q/pi/p)^0.5;7 ^8 H4 P4 w" ^3 u8 ^: n( U kg=mu0^0.5*kb*kd*j0*g2;* H, [5 l! L" f3 o3 w( ^ rx=ku*kx-kp*kg;1 |3 ]; Y7 t' p! C+ k& P ra=ku*kg;6 r; s+ [# y# B" k1 @ rc=kp*kx;5 G! s' n$ O; ^* B+ h/ r7 G kj=mu0*kb*kd*j0;0 G9 x5 y0 u' n: }5 J" X8 _7 B kj=kj/kx; $ ^% C0 l( w' j3 {7 fh0=3*ht;; f! i4 V$ z- h2 o! \ h1=2*(1+1/rou)*q/pi;5 ]+ N5 }7 A5 t8 m h2=2*(kd-1);( N# d0 m Q- j* Q5 o( ^% [ h3=1/bm*(mu0*p*q/(2*pi))^0.5;% c$ ^4 ^' Y% `$ I8 R! W h4=1/kj*(mu0*p*q/(2*pi))^0.5;: Q9 ?& q; Y0 t- Q( M# a' {- o8 z hg=2*pi*kr*ht; " p0 q6 i1 ?: W0 w3 Gw1=gaf/bm*(pi*p*q*mu0/2)^0.5; . ]0 h/ O9 Z, y( i9 [( w' Iw1p=8*g2*kd*w1; 8 P$ ^3 x, }& g+ aw2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;8 v, G+ K! O: x5 P* a w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;% C1 j( l9 U6 F w3pp=hg*h3; 2 \( C4 S5 [' ?% F( o* V7 k5 Ew3=w3p+w3pp;; L; F0 r( ?. N4 v4 Q w4=2*kd*gacu/kj*(2*pi*p*q*mu0); , O# o6 g, _! F) f4 e* rw5=5*pi*ht*gat; \/ v, z1 k% c9 Q1 [w6=hg*h0*gat; * v0 V, L. v7 a5 n& ^2 k# e6 E1 Ew7=hg*h1*gat; ) ~0 k0 Q. T8 Zw8=hg*h0*gat;+ _ Q: L* g; F1 e$ A# T w9=hg*h4*gat;4 a3 h4 G9 q4 l" J8 B" z 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;

这是后来修改后的，请高手指点。; ?3 I. Z5 n* H5 }2 B

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