数学建模社区-数学中国

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

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

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

: r# u+ L0 ?' x7 e. e5 O8 H

function f=fh3(n)- J1 a7 \: s5 U' O. q0 @ g2=2^0.5-1; - n/ T" [% i4 ]% N7 q9 ?6 Umu0=4*pi*10^(-7);! d$ C3 g0 R6 X4 q gaf=7.8*10^3;5 N) s V' Q j, {& a gacu=8.9*10^3; ' j% X) P! c2 E* Y& ugat=4.5*10^3; 0 O7 v( t6 O3 b( ~6 |3 {6 ^gaa=2.7*10^3; 0 l7 Z5 I/ N: N5 E; I6 h' u2 s: mp=2.7*10^4;' ]# P5 S* a: P1 s+ @0 ?& f q=3*10^(-5); " l8 _: T2 C' d9 F: Erou=0.7;6 B5 i8 A" B; e& A kb=0.9; ( J, Y8 q" T9 E- M- G: o& B%j0=8*10^6; 2 q1 B8 V G. j. r0 D% j0 H [j0=2*10^6; , F( u- n+ k( w9 i2 ]* |0 M. mkd=3; ' {9 F/ Z# z9 q0 cbm=2; # E- j4 G# f4 }- b! [ht=2*10^(-3); * _) P$ m( s2 l) x2 b' T* Tmu1=10^8/1.75; $ [- ?( _7 L" }kt1=2*pi*gat*ht; 7 c6 p% }+ S5 e9 h& W/ Xkt2=2*pi*gaa*ht;- \$ m6 s0 ^9 z6 z t2 v kr=1.02; ) ^0 d8 |' I D& \, b. D# bku=(mu0/2/bm)^0.5; ' D2 X' H2 J( x$ P) f& ]kx=mu0^0.5*kb*kd^2*j0; " A# Z0 ^7 M% g+ h9 S. |$ P0 M3 g! {3 Lkp=(2*q/pi/p)^0.5; % I5 W% G2 S X7 ]2 Kkg=mu0^0.5*kb*kd*j0*g2;2 t2 S4 t" ^/ S' C rx=ku*kx-kp*kg;9 d/ e3 {0 l- x6 a+ h ra=ku*kg;. M A& ]" R' i1 j: b% v- @' W rc=kp*kx; 3 i& a4 u' o: R- nkj=mu0*kb*kd*j0; 4 C7 O! }* P0 `% s6 tkj=kj/kx;- j2 A0 h, _; U; s$ r h0=3*ht;- U6 i$ i& Q5 H K j& F3 S8 { h1=2*(1+1/rou)*q/pi;* U. j4 S% _$ s h2=2*(kd-1); ) |8 e, @% q0 B, t) z7 `h3=1/bm*(mu0*p*q/(2*pi))^0.5;$ Q& L- b( J* b6 E1 U h4=1/kj*(mu0*p*q/(2*pi))^0.5;* l9 I; \6 V6 ~* i$ U hg=2*pi*kr*ht;$ Q2 A$ g; E# T$ s& t: P w1=gaf/bm*(pi*p*q*mu0/2)^0.5; 1 ]( ]& m) ^& |3 @w1p=8*g2*kd*w1;; O. \0 A4 i G( c, K& z w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; 0 ]* q: z# ~) M* \; v/ gw3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5; 1 r, i3 W8 i) i" F& ?' w, Q' Ow3pp=hg*h3; # M7 \: o) j& A/ s; y1 i: lw3=w3p+w3pp; ' Q! g" j4 a, } `" s+ r) uw4=2*kd*gacu/kj*(2*pi*p*q*mu0);8 K+ r+ _8 K# r3 Q" O% O w5=5*pi*ht*gat;) U/ C$ T- d( V5 m2 n* ?) K w6=hg*h0*gat;0 u4 I$ T3 }1 ^/ d! U T w7=hg*h1*gat;) `% l/ ~4 W9 t7 l w8=hg*h0*gat; 9 I5 d& O! ^0 B- U4 Y' xw9=hg*h4*gat;% s) q- b' X8 d( S8 I2 K) G7 E% x 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); , Z" f5 b9 s) A: _# b+ I


作者: 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
不会吧!如果导数不连续和不可导是差不多的性质!你看看是不是半光滑的!. Q; P- }1 r0 e
[此贴子已经被作者于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;8 { T; }2 Y& l' |! M- J mu0=4*pi*10^(-7);! `* s3 t2 v' [* j4 m8 W( b! r( q gaf=7.8*10^3;; _+ z/ P! q- a1 u4 ~6 F: @+ j gacu=8.9*10^3; # u g% X: K3 B0 ~& \gat=4.5*10^3;6 |/ T) i2 B# R: y$ J gaa=2.7*10^3;2 p5 t! [& x. C/ t7 Y. ] p=2.7*10^4;1 ]) O1 @: E4 E" w q=3*10^(-5); $ S9 X% o' U) I- v& N- ~( D; Nrou=0.7;$ J% M! R, e0 Y" A7 l0 _) c kb=0.9; ! `- q0 t: l) v2 P; R%j0=8*10^6; " {# k, o- H. ~( d( ~/ F, P! Hj0=2*10^6;/ L$ M$ f; ^6 P d+ }$ j1 [3 ?- C kd=3; 7 m: b7 Z7 B8 d. Ubm=2; r N2 n3 L! O" s' y9 q: F9 B ht=2*10^(-3);& g+ q2 r% i. X- S9 g: o mu1=10^8/1.75;$ P# J. D" u" J+ Z9 h kt1=2*pi*gat*ht; " ?: I. ?( w8 I3 l* Qkt2=2*pi*gaa*ht; * R/ W0 w w6 J3 i0 O6 |1 qkr=1.02; 2 y+ ]. i& J* `( W4 a- Xku=(mu0/2/bm)^0.5;& u' G' i! y7 [( O" J* N9 h kx=mu0^0.5*kb*kd^2*j0;# z& w+ g) f. K3 I, | kp=(2*q/pi/p)^0.5;. k. ^ `0 l. T. }& f kg=mu0^0.5*kb*kd*j0*g2;# S7 m1 \# F4 [: G p rx=ku*kx-kp*kg;( |# P4 @2 J: X, E ra=ku*kg;0 a8 \0 K$ c7 v! X) {7 Z rc=kp*kx;. Q- [/ h+ l4 u) ^; ~% F# j) \$ [ kj=mu0*kb*kd*j0;" Y/ _3 K& X9 \ kj=kj/kx; ( b. e% j3 R1 g- f( lh0=3*ht;; c! t- J0 Y E$ ?7 _0 |3 O" E8 J6 E) P h1=2*(1+1/rou)*q/pi; ' c3 D c* G* Uh2=2*(kd-1); # @7 x4 V6 h: K% E% `h3=1/bm*(mu0*p*q/(2*pi))^0.5; % W7 Z0 ?) x5 x' `$ m2 i; Th4=1/kj*(mu0*p*q/(2*pi))^0.5; v. V8 E' J6 ]* ^& ] hg=2*pi*kr*ht;$ o# O/ Z& A2 X4 }: A! y w1=gaf/bm*(pi*p*q*mu0/2)^0.5; " j. `# w; m! A1 ]w1p=8*g2*kd*w1; 0 |( r# X. v a' |! W, bw2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; 8 C7 \0 p; k4 z9 _" O5 Bw3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5; 1 U4 r4 W& x2 }5 Ww3pp=hg*h3;. Z3 v# [4 c5 T# | w3=w3p+w3pp;) A3 \: G9 ^1 o; R w4=2*kd*gacu/kj*(2*pi*p*q*mu0);. X' f- C4 g& I* c' t w5=5*pi*ht*gat; 4 c1 P( k- u1 Rw6=hg*h0*gat; 8 d7 O% z# C- ?1 {- I3 l8 {, Lw7=hg*h1*gat;" R9 b" Y2 O5 E9 \- ? w8=hg*h0*gat;& x9 C! G# u' l w9=hg*h4*gat; + Q, x( R8 D7 {* bmin=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;

这是后来修改后的,请高手指点。$ c3 @# V7 s9 ?3 s9 z( x






欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5