QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5666|回复: 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之间。

) l: O9 H. ^, q0 z, O

function f=fh3(n)% }2 p( b- x7 b+ `5 T: U& `- k g2=2^0.5-1;" I7 B# Q; ]* m6 g2 l3 `2 d mu0=4*pi*10^(-7);5 ^( Y) o f9 }) B" ]% Y gaf=7.8*10^3; . S- u: i0 f0 }6 H4 _gacu=8.9*10^3;+ S# g) |8 X; K! y. p7 a gat=4.5*10^3;: B% l/ l7 Z+ K+ W. M6 S9 t3 r gaa=2.7*10^3;/ Z' `! B; v3 S* w1 ?+ t3 U p=2.7*10^4; . ]. M2 ]( R: H4 h+ a m6 I: pq=3*10^(-5); : L5 Y) v& ?* r7 {4 Z. g |rou=0.7;2 i% e, J9 C% P+ M& [# D- e& k* k B g kb=0.9; ; v5 k, C# b0 K/ D2 w% F%j0=8*10^6; ) j! o5 e3 w; s4 M2 rj0=2*10^6;) m6 x" B# S7 g: u P( _- }3 l4 p" T kd=3;9 t, }) u9 S: {6 {( ]8 N bm=2;; \1 T- ?1 ]5 n8 Q' t' C8 n5 {9 Z ht=2*10^(-3); 0 Z \; D5 z4 Mmu1=10^8/1.75; 1 q6 `7 z9 v+ K0 o$ Wkt1=2*pi*gat*ht;( ~% H. `' ]7 l; u+ ~ kt2=2*pi*gaa*ht;3 t$ o5 o( J/ B* h4 a$ F kr=1.02; & j8 b5 J8 P& o5 Fku=(mu0/2/bm)^0.5;6 \0 C0 g. }- \: K& v6 l kx=mu0^0.5*kb*kd^2*j0;% \- l, q& w; n! C/ ~6 V' E kp=(2*q/pi/p)^0.5;# @) ^; A$ k" o3 ^: ] kg=mu0^0.5*kb*kd*j0*g2;+ ~5 V4 w) X. C, h- Z$ F/ v rx=ku*kx-kp*kg; 5 x0 g" A/ w) I* _; f0 xra=ku*kg; $ ~/ K1 e2 F, q1 `' w! I! A/ vrc=kp*kx; 0 q- g2 V' g. @; S, ]# Z% [/ Okj=mu0*kb*kd*j0; 7 k7 q2 z; r% o' C$ tkj=kj/kx; / }/ g( J' Z, Y% Q7 x5 W# ^! K# X6 Ph0=3*ht; 7 {7 c- Z' g2 M' W" G( Ih1=2*(1+1/rou)*q/pi;( Y9 }" x% V9 @ h2=2*(kd-1);% T ^8 O+ ^+ N- d: D; u$ S9 h% ~ h3=1/bm*(mu0*p*q/(2*pi))^0.5;. S: R$ q# }% E- p3 Y h4=1/kj*(mu0*p*q/(2*pi))^0.5; N4 O+ l5 z! y1 i7 |4 Fhg=2*pi*kr*ht;6 i5 S8 z$ E8 e4 Z w1=gaf/bm*(pi*p*q*mu0/2)^0.5;- g5 U0 \* F2 `4 D& W w1p=8*g2*kd*w1;, H M, k' l. ?$ S' u8 [ w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;. t6 C ^, {3 }5 l! ] Z5 O7 ]% z0 } w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;( o8 r/ k# q) G+ J% T, k3 X4 `# I w3pp=hg*h3; 8 d: o2 T8 _* L6 v, _# }w3=w3p+w3pp;8 i g4 O: }0 Q/ l! c; g- N w4=2*kd*gacu/kj*(2*pi*p*q*mu0);5 U/ n& e1 ?8 T. l w5=5*pi*ht*gat; ) X: f0 C' F1 `4 @w6=hg*h0*gat;* x2 `$ l5 J- G; [ w7=hg*h1*gat;& N& y6 L+ l0 B* F. C* h5 M% I w8=hg*h0*gat;1 }- P+ ]: \; w! o! w w9=hg*h4*gat; ! @1 x6 C3 \1 S- {& af=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);; O9 M0 m: X+ V( ?8 z9 K! P

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-4-20 05:03 , Processed in 0.588248 second(s), 99 queries .

    回顶部