QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5802|回复: 13
打印 上一主题 下一主题

请高手指点

[复制链接]
字体大小: 正常 放大
zhanfeng        

8

主题

2

听众

152

积分

升级  26%

该用户从未签到

跳转到指定楼层
1#
发表于 2004-12-17 21:58 |只看该作者 |倒序浏览
|招呼Ta 关注Ta

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

: G* q, c/ h; U1 @5 F! O5 |

function f=fh3(n) 1 k% u+ u2 D2 B2 b8 y- p9 r7 j5 rg2=2^0.5-1;7 O1 w+ T% |; l0 k mu0=4*pi*10^(-7); 6 n) |7 W. x, h& xgaf=7.8*10^3; 3 O! x& |( o$ W" y) y( P2 igacu=8.9*10^3;' Q- n6 r! M" W( B: P gat=4.5*10^3; " B r# T- l G, @& ~2 h$ cgaa=2.7*10^3; , k! R0 H+ B4 Wp=2.7*10^4;2 P+ }3 O. b" ~( G q=3*10^(-5); 5 n8 M6 w5 h- ^. _) vrou=0.7;" V9 `# n# S5 e) l0 y kb=0.9; & F( J) r! X5 A* g" Q" C6 A/ c%j0=8*10^6; # W8 ^. r+ D+ N8 [/ e4 }$ o/ Oj0=2*10^6; 8 ` E0 Q- v& ?$ ]' u( {5 h& @kd=3;6 [* J: v6 J) P; Z5 G' `0 V, M: F bm=2;7 o( E) F7 o6 u# `6 T ht=2*10^(-3); + E5 v' z0 V1 w6 Jmu1=10^8/1.75;6 {7 l6 ~$ n6 k& ~2 J/ F kt1=2*pi*gat*ht;. }8 p" r- q' g* z3 @ kt2=2*pi*gaa*ht;* s( ?4 [& |8 u+ t, J kr=1.02;+ H9 k# E# T5 |. I ku=(mu0/2/bm)^0.5; ( r' H! u3 |# q) k6 [' q7 ]kx=mu0^0.5*kb*kd^2*j0; % r* ]" { ~* z# Ykp=(2*q/pi/p)^0.5;" l3 D( p5 ^! J5 A. B w4 A# u; e t kg=mu0^0.5*kb*kd*j0*g2; 8 J/ |* @7 _, K q+ d' ^" rrx=ku*kx-kp*kg;# u" m* G. { _- f3 H4 q ra=ku*kg;' J9 ]' j9 ~* y8 S D rc=kp*kx;7 N1 _$ N. j- O6 ? kj=mu0*kb*kd*j0;6 ^4 h/ J6 v% ^1 o. d9 x/ Y7 x$ q kj=kj/kx;) h/ F- Y3 c) B; ?( ]6 r" n h0=3*ht;! ^8 U3 ?4 b- w; B& Q) m* F h1=2*(1+1/rou)*q/pi;7 V4 h# m& I2 C1 ~ h2=2*(kd-1); " R3 Y# Y8 S& ~h3=1/bm*(mu0*p*q/(2*pi))^0.5;* Y! c, u) R' y- N3 M! q h4=1/kj*(mu0*p*q/(2*pi))^0.5; % Y. S" a% p+ Bhg=2*pi*kr*ht; , ]; D( Z& E/ ^4 ~w1=gaf/bm*(pi*p*q*mu0/2)^0.5; 2 D) J/ H3 J4 ]8 G# B" Ow1p=8*g2*kd*w1; 0 S( b: G) @" Xw2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; V. }- m* L! e+ m- Fw3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5; ) f# D2 Z+ K& }! aw3pp=hg*h3; ; o9 y: i9 C1 L% w. \/ F8 x- fw3=w3p+w3pp; 5 d- d" J3 C9 x3 O1 J8 e0 {2 V) kw4=2*kd*gacu/kj*(2*pi*p*q*mu0); : W1 W7 s8 I! i. ?* fw5=5*pi*ht*gat;6 {3 J! W0 L$ Z w6=hg*h0*gat; * p0 `9 F2 p. qw7=hg*h1*gat; # s6 F8 d7 E7 L8 N" X8 Yw8=hg*h0*gat;2 g1 i# k# n, U9 T1 P* P2 H8 G6 q w9=hg*h4*gat;0 z. C" T8 A5 x* l$ x0 Z 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);) \, f: _' T+ P" P ~" x

zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
www.madio.net 相见恨晚! 众里寻她千百度,蓦然回首,那人却在灯火阑珊处!
zhanfeng        

8

主题

2

听众

152

积分

升级  26%

该用户从未签到

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

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

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

请高手指点!!!

www.madio.net 相见恨晚! 众里寻她千百度,蓦然回首,那人却在灯火阑珊处!
回复

使用道具 举报

zhanfeng        

8

主题

2

听众

152

积分

升级  26%

该用户从未签到

回复

使用道具 举报

zhanfeng        

8

主题

2

听众

152

积分

升级  26%

该用户从未签到

回复

使用道具 举报

zhanfeng        

8

主题

2

听众

152

积分

升级  26%

该用户从未签到

回复

使用道具 举报

ilikenba 实名认证       

1万

主题

49

听众

2万

积分

  • TA的每日心情
    奋斗
    2024-6-23 05:14
  • 签到天数: 1043 天

    [LV.10]以坛为家III

    社区QQ达人 新人进步奖 优秀斑竹奖 发帖功臣

    群组万里江山

    群组sas讨论小组

    群组长盛证券理财有限公司

    群组C 语言讨论组

    群组Matlab讨论组

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

    回复

    使用道具 举报

    zhanfeng        

    8

    主题

    2

    听众

    152

    积分

    升级  26%

    该用户从未签到

    回复

    使用道具 举报

    Allos        

    13

    主题

    2

    听众

    103

    积分

    升级  1.5%

    该用户从未签到

    回复

    使用道具 举报

    zhanfeng        

    8

    主题

    2

    听众

    152

    积分

    升级  26%

    该用户从未签到

    回复

    使用道具 举报

    zhanfeng        

    8

    主题

    2

    听众

    152

    积分

    升级  26%

    该用户从未签到

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

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

    敬请高手指点!

    www.madio.net 相见恨晚! 众里寻她千百度,蓦然回首,那人却在灯火阑珊处!
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-9-13 20:18 , Processed in 0.471194 second(s), 98 queries .

    回顶部