数学建模社区-数学中国

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

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

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

function f=fh3(n) $ e- ~( V6 ] wg2=2^0.5-1;8 U# r1 d+ k8 e' u& m$ O' u mu0=4*pi*10^(-7); 8 c& R0 Y$ h9 S: ]/ {* bgaf=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); / |# A& e, @/ t' ~7 J9 R' @9 ^rou=0.7; # Y' `8 G1 @' ~' q+ skb=0.9; ! m5 e6 |# N5 s+ {%j0=8*10^6; - J/ S. n# t9 R3 U9 ?' O# N- c3 Ij0=2*10^6; - ~- [/ b3 g/ ekd=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; 7 v- s7 c2 a- ?3 r9 y/ C9 |0 Q* pkt1=2*pi*gat*ht; 1 m8 p0 S5 j. T0 fkt2=2*pi*gaa*ht;9 h/ B7 e5 A$ m9 t+ l kr=1.02; . q' l1 e. C3 E, sku=(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; % X4 F: }7 j3 K3 Trx=ku*kx-kp*kg;1 _6 ~( n, F0 ` ra=ku*kg; " Q( t7 u E3 M( p+ ~) Arc=kp*kx; & N. O$ V5 @6 T) U; M3 P* xkj=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); 2 s" A0 `: I- r6 Oh3=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; 6 k2 ]8 b" r- x0 t" A; X* phg=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; 1 P1 w9 K" R2 I0 S3 |w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; ) N& H( d2 z+ @# A' y: e( vw3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5; 4 m6 A4 d8 S7 \, L/ O9 _w3pp=hg*h3; 9 U# D- `; A$ q- p/ D4 fw3=w3p+w3pp; " K3 k$ s- W: [, p/ J/ rw4=2*kd*gacu/kj*(2*pi*p*q*mu0); : h7 Q/ t0 v& |, ~w5=5*pi*ht*gat; B0 T x S% n" iw6=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; 5 |# C+ j: N; F- ?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

作者: 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; . R& G$ R$ @# _6 `. C' Zmu0=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; ) O# k: H }# j* O9 xgat=4.5*10^3; + A3 o/ Z: B" T5 \$ }+ bgaa=2.7*10^3;# V" g1 j- [* L p=2.7*10^4; " S2 D( _' |- f8 q: o: x' ^q=3*10^(-5); 0 K5 R0 F: ?; B6 k0 G- Trou=0.7;; c! C9 p7 C1 D% g! E: ]& S5 B8 W kb=0.9; + Q8 W% R7 {: H4 s/ i%j0=8*10^6; 6 f0 c: f; n6 f8 D& k; Z+ r3 cj0=2*10^6; 4 g& ]+ m/ H5 l! `. jkd=3; # r& }9 j% Y7 z$ N/ mbm=2;* _6 Q# j2 ?9 { ht=2*10^(-3); * t- m) h6 ^2 A" C$ Y+ i# _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; / p$ L1 B5 t0 |6 ukr=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; 7 H2 g% @5 W. Ckp=(2*q/pi/p)^0.5; % y) m& g7 \1 p6 |8 e; |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; 2 s. l1 Y1 a8 ^$ @3 {0 ?$ {5 s0 C. @rc=kp*kx;, j3 o& @4 y& R( U kj=mu0*kb*kd*j0; % W6 f7 i8 m2 B; m& Mkj=kj/kx; ) }- V" o8 E! M, ?& ?' ^h0=3*ht; + k) e& v" B6 B8 Q2 J- _- d5 th1=2*(1+1/rou)*q/pi;3 [5 z1 p7 L) h7 R( V$ q h2=2*(kd-1); : x& y/ a3 L: gh3=1/bm*(mu0*p*q/(2*pi))^0.5; ! L# L' y1 w5 eh4=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; $ r( Y- @: {0 }) hw3=w3p+w3pp;% [! d) G+ `6 l" d3 C- n8 F( O) } w4=2*kd*gacu/kj*(2*pi*p*q*mu0); ( m* D% Y8 S9 b/ v% o9 `3 ^7 Jw5=5*pi*ht*gat; " ^3 V$ P* }* N: rw6=hg*h0*gat; 5 q% A9 M, ^) q9 H8 ^- n4 \w7=hg*h1*gat;! S% L+ ^' t- S4 k5 I w8=hg*h0*gat; 6 z) _) g3 h8 j, t! z' [# J) zw9=hg*h4*gat; & C& h+ E' T7 s: G! tmin=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/)