QQ登录

只需要一步,快速开始

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

' c) r/ n0 `7 e% O. ]

function f=fh3(n) 5 Z5 O- N) z3 H( c0 Z' Og2=2^0.5-1;2 S4 I2 S2 ?' u3 J mu0=4*pi*10^(-7); 9 V2 w; s) b2 ggaf=7.8*10^3;" Y. U8 @' l4 g7 `% O* Y2 [ gacu=8.9*10^3; ; g0 Q N' A9 Q: v9 K5 K" M# a2 n: h7 Qgat=4.5*10^3;; M p$ `; N) g4 s gaa=2.7*10^3;/ f3 V1 |! j, G+ L m) ] p=2.7*10^4;, Z/ g. u9 G9 h- _; b7 \7 J G q=3*10^(-5); ) \1 D: s: o# t/ p) N$ @rou=0.7; . b4 ~6 _4 _9 f: W8 c9 Q, ikb=0.9;: F6 D6 b8 D" f0 O3 a! q4 O %j0=8*10^6; , \9 k- K- v$ i+ N7 Y) _" Tj0=2*10^6;- s# B4 |& ~6 f kd=3; 1 H9 A8 r2 q$ y# \bm=2;! _8 A+ f6 { h: h; v ht=2*10^(-3);) w& d I, j9 O4 I9 {3 m8 z, g mu1=10^8/1.75; & \: E6 T, P; s; |% Okt1=2*pi*gat*ht;) H. S7 O7 l0 C0 ?' `* Y) b" U$ y6 Q kt2=2*pi*gaa*ht;4 l1 ^+ `) B7 X) c! Z kr=1.02;9 E/ R+ v6 G! M" J' k% G& D ku=(mu0/2/bm)^0.5;: b2 q* g v$ [' @$ L( @ kx=mu0^0.5*kb*kd^2*j0; 1 H' ?/ `4 k9 B5 D7 g* G# o4 W# {kp=(2*q/pi/p)^0.5; 3 E0 |* d2 F9 r7 ?8 `kg=mu0^0.5*kb*kd*j0*g2;: k6 B) W. x. W4 m* x7 E rx=ku*kx-kp*kg; ' z" y& D/ x- d7 a$ [. c5 I5 Xra=ku*kg;0 ]& l; q$ Y& B' D9 ^6 X8 O" c, E rc=kp*kx; 0 l; X- h7 V0 p$ g1 i! d) U6 _kj=mu0*kb*kd*j0;, s1 A) g' O# E/ G! `6 q8 M kj=kj/kx;* R& i' G8 P; r h0=3*ht; ) ^" B2 e5 @! w3 w4 C0 ?h1=2*(1+1/rou)*q/pi; ' d: G$ w( j% A2 }6 Sh2=2*(kd-1); 2 F; G% g# {7 g* ~ S/ z/ R% ah3=1/bm*(mu0*p*q/(2*pi))^0.5;4 }, g7 Z7 u6 o3 Q h4=1/kj*(mu0*p*q/(2*pi))^0.5;- K \: ]8 N( i) l2 z hg=2*pi*kr*ht; 0 @$ z: Y. N: v9 K9 ^w1=gaf/bm*(pi*p*q*mu0/2)^0.5; - A7 s; B3 t" r! P& T9 R$ a1 qw1p=8*g2*kd*w1;" I W1 ?) V2 d& V w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;' ?/ l, _$ z4 ]# D w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;+ e" H* N5 r* M' P* y w3pp=hg*h3; % R, g" T' I7 x& gw3=w3p+w3pp; 1 S! X1 t1 I6 I0 t# M) R7 H. B4 }w4=2*kd*gacu/kj*(2*pi*p*q*mu0);4 t8 u0 \0 |. R! l9 K! F+ l w5=5*pi*ht*gat; ; o. w( C( M3 ?1 s0 x1 X0 Uw6=hg*h0*gat;. |$ `; z# Q" H" o/ ^ w7=hg*h1*gat; - H: B1 p. E& F' t& ww8=hg*h0*gat; - g, W: ` h% O* @w9=hg*h4*gat; - r/ y8 V1 l: h3 \& h9 {0 ~1 p% sf=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);5 W+ p0 k( G# O2 H4 J! _1 ?

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-18 01:17 , Processed in 1.040359 second(s), 98 queries .

    回顶部