QQ登录

只需要一步,快速开始

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

2 I, s5 x6 d8 a C& ~% M

function f=fh3(n)$ m8 d* S- _2 l1 D* k% C7 _# l g2=2^0.5-1;4 |! x% `: L8 K6 F4 y1 B mu0=4*pi*10^(-7);, _+ V _2 O' S gaf=7.8*10^3; # }+ G3 g- A+ z. Q9 ]9 ~! i- Igacu=8.9*10^3; e5 i' x% G4 h; G2 jgat=4.5*10^3;6 e& l1 p) Z: [% C+ s) h6 \2 x" A gaa=2.7*10^3;* E7 e8 @' o3 t3 ^ p=2.7*10^4; + Z; V9 e/ U) R, E; Y2 d4 Iq=3*10^(-5); `4 l& a/ P% e; @% Wrou=0.7;5 H: a. s: x/ A1 `7 P1 S7 ]$ k kb=0.9;3 M! l3 s' f& E: C! N9 R- e( A9 c %j0=8*10^6; 6 h2 p$ j( P2 c2 k" {% Ej0=2*10^6;7 [9 j# L4 _% V+ q" M2 ~4 f# e kd=3; 5 Z: n' w3 f: U6 G/ C' Xbm=2; " O! |8 J; ~3 [. J: O, @ht=2*10^(-3);5 a# u5 i; y) Q3 R- R8 i7 @" j8 w+ ^ mu1=10^8/1.75; 8 {9 U$ i5 f j3 c0 Z* D* E2 Bkt1=2*pi*gat*ht; 1 ^' T' [7 }2 `* ]; m. P; e/ k4 ykt2=2*pi*gaa*ht; 4 g* o( }9 ^' Y& W- C2 kkr=1.02; " k4 P; D9 s! h$ O. g2 m3 mku=(mu0/2/bm)^0.5; ) H5 }5 e, `. h, g$ z1 E+ B5 wkx=mu0^0.5*kb*kd^2*j0;' n; M: I& U2 t* f5 O S kp=(2*q/pi/p)^0.5;9 h( \6 l" d4 ^) O kg=mu0^0.5*kb*kd*j0*g2;% ]0 {, U7 _$ ]6 B3 o rx=ku*kx-kp*kg;; ]: O: ^- I: w9 ]/ N% o# l% I$ L& g' C ra=ku*kg; 4 i+ H7 E' L* Y0 T2 Irc=kp*kx;. A3 }9 [0 u% \2 ?% k7 O4 p9 } kj=mu0*kb*kd*j0; # s" [7 A) f5 I, J0 Ikj=kj/kx; 8 ^3 y1 ^. _/ }) ]5 L: j5 e( J: A( B, mh0=3*ht;2 p; E$ m+ b. V) e) ] h1=2*(1+1/rou)*q/pi; , k' j- A: ~* m/ jh2=2*(kd-1);3 N a0 o( e6 F5 a1 o; Y h3=1/bm*(mu0*p*q/(2*pi))^0.5;, ~9 }# b& F% B* _2 u6 b h4=1/kj*(mu0*p*q/(2*pi))^0.5; \ y v; |; Chg=2*pi*kr*ht;5 D% y2 \# f% m w1=gaf/bm*(pi*p*q*mu0/2)^0.5;7 j' o! _% j( K/ M4 p3 [: Z5 e w1p=8*g2*kd*w1;2 t2 @/ V+ ~# ?. E, P0 E w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; e5 h7 W0 A7 Yw3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;+ w% R& f2 T. D5 M9 q5 T9 y5 t6 u2 l w3pp=hg*h3;4 `( C, ^* g. d4 r' y, M5 O w3=w3p+w3pp;! m: G4 g% w! P8 k( I w4=2*kd*gacu/kj*(2*pi*p*q*mu0); + B% x+ e* }/ x" |7 ^" iw5=5*pi*ht*gat;8 I" T u4 [4 t7 R# s4 o- o w6=hg*h0*gat;' d j. k9 w. r, X$ K5 L5 S$ T$ x* F w7=hg*h1*gat;4 H* w7 Y5 n/ I v$ G* y3 ] w8=hg*h0*gat;0 B0 F( T7 I" ?' t/ j w9=hg*h4*gat;' r0 K9 v, }4 s9 O; 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); % F. ]1 w! y% y2 F

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-17 22:02 , Processed in 0.439774 second(s), 99 queries .

    回顶部