QQ登录

只需要一步,快速开始

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

. [ j( ]: J4 g

function f=fh3(n)0 _, Q1 R+ Q6 p/ n1 t; r g2=2^0.5-1; % f* h' k8 d6 O% b$ K. `; @1 \mu0=4*pi*10^(-7); $ B! I: H$ ~1 U e! h7 g9 v* Y/ _gaf=7.8*10^3;9 w5 f9 [3 O! I gacu=8.9*10^3;- f$ p# W! ^+ P! R3 p y" @ gat=4.5*10^3;3 R6 l, h' ]% B+ R, }3 n gaa=2.7*10^3; 4 Y! O0 b9 j/ F4 Pp=2.7*10^4; ) M- ]: x1 Y/ uq=3*10^(-5); 0 u* w5 ^1 |& H2 ~. ~rou=0.7; + `+ g- Z' m$ k9 b* @* W* Okb=0.9;) z- \- N7 f7 ^" X5 V %j0=8*10^6; . E. @! T! C) W oj0=2*10^6; 0 {' X9 J. F% e9 X5 |: Ukd=3; ( N. H8 Y2 t7 F* s7 U. Dbm=2;9 U* A& I. V, A9 s) Q2 v& A" Q( O ht=2*10^(-3); 7 l7 B: \ M8 B% o9 Hmu1=10^8/1.75;5 l8 Y/ X; K/ F kt1=2*pi*gat*ht; ( y1 ?! n" p$ j$ fkt2=2*pi*gaa*ht; $ w9 Q7 T: F- u) N5 W' P0 t \) i; @' z0 zkr=1.02;* I4 g( a2 x8 ~* ?4 f7 l ku=(mu0/2/bm)^0.5;$ F/ _0 z5 M7 Q kx=mu0^0.5*kb*kd^2*j0;5 u, ^0 h/ H: ?' q% A; S* X kp=(2*q/pi/p)^0.5;/ K4 |" F% C/ t: Y3 w9 e kg=mu0^0.5*kb*kd*j0*g2; ; h6 H* r" E) }# ^" Mrx=ku*kx-kp*kg; 6 l) j/ G% u' i/ {. C! Cra=ku*kg; , q* }8 g% U4 a. erc=kp*kx; 5 z2 g; _( s7 `" c' c9 lkj=mu0*kb*kd*j0; 3 ^/ d6 |0 Y& Lkj=kj/kx; 6 j) L% g# ] Lh0=3*ht; , g4 |; B! l& y) p0 @; uh1=2*(1+1/rou)*q/pi; G8 s8 s# N. m, j, [1 E+ F h2=2*(kd-1); $ g) @. a! Z, d2 [# G* s0 P8 ^h3=1/bm*(mu0*p*q/(2*pi))^0.5;" a6 x9 \- t7 n$ b: c3 Z h4=1/kj*(mu0*p*q/(2*pi))^0.5; 7 n6 y( F) J) Yhg=2*pi*kr*ht;7 r0 b" n8 D: j3 m) t, K w1=gaf/bm*(pi*p*q*mu0/2)^0.5;- B( O2 ?- k/ a, h8 Q( } w1p=8*g2*kd*w1;1 K8 G2 @! X4 @& y; |8 N! n- Y w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; + `* ~. {2 B2 u) L" i7 O3 I: Lw3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;, l* W& w$ H) g6 N+ |& Z& t w3pp=hg*h3;' K: U, U V2 x0 e0 i w3=w3p+w3pp; : w; r5 j& {8 v; Iw4=2*kd*gacu/kj*(2*pi*p*q*mu0); ' L- v6 H! O7 P% pw5=5*pi*ht*gat; 5 H8 `0 c" j9 p+ q, C/ q& x* hw6=hg*h0*gat;6 _1 l7 ^. g, z( k$ T. m w7=hg*h1*gat;, ~+ y8 B0 W0 ] w8=hg*h0*gat;) [) |& \% S0 R y' q) ~2 k5 u w9=hg*h4*gat;: W* t( }) |/ Y6 S) @, { 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);2 I$ ~% A2 x- S" r

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-6-11 16:47 , Processed in 0.509057 second(s), 99 queries .

    回顶部