QQ登录

只需要一步,快速开始

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

$ {1 m+ s' L0 p% s% r

function f=fh3(n) 2 z7 X( I/ q$ [g2=2^0.5-1;6 c( e9 m- }/ [4 U mu0=4*pi*10^(-7); : C6 [# t! L7 F! }: ^gaf=7.8*10^3;. _9 d0 Z- D O4 B: \8 Y8 t" p% b3 \ gacu=8.9*10^3;9 C3 P% |% x% ]" i gat=4.5*10^3; % [# X( s% ^! y9 pgaa=2.7*10^3; 0 y+ |: P* G$ S1 \p=2.7*10^4; & }0 K' u- f* Z( ?+ u5 aq=3*10^(-5); & H! R- [: \$ |# _4 Krou=0.7;1 i8 L& Y: F' r- F; W" q kb=0.9;& p/ a0 b! q; Z1 {8 b, q %j0=8*10^6;1 A+ ]6 t* z& X3 H/ \8 M2 n( R j0=2*10^6; % j! q- }2 }( ~$ r p2 w3 A- @kd=3;, c2 S5 ?2 q: Y- n- L5 c bm=2;9 f% Q" t4 a8 O1 P* J5 i# I ht=2*10^(-3);! P3 N: I7 m3 V/ v6 i* [ mu1=10^8/1.75;- W! M1 m) u- k kt1=2*pi*gat*ht; " ]) O0 m6 M/ k" Y5 x; x6 _kt2=2*pi*gaa*ht;2 t0 Y( g- k* n$ H% _( v kr=1.02;: ]# v+ g% l& N6 y ku=(mu0/2/bm)^0.5;7 M% z+ j2 E4 J5 p1 T9 i3 c, q: j4 V kx=mu0^0.5*kb*kd^2*j0; 1 L d( }" l( l" I( Rkp=(2*q/pi/p)^0.5; & _' M* D8 p/ _% K4 D1 e0 \. o4 [. kkg=mu0^0.5*kb*kd*j0*g2;% o% u, k& t$ y6 z8 _$ e rx=ku*kx-kp*kg;" J# {2 P& K, ^. x5 o/ ^" G) @ ra=ku*kg; / A4 ]# z2 f* ^& g0 U* prc=kp*kx;* J7 Q) r1 K$ I, |2 A+ Q E kj=mu0*kb*kd*j0;7 o+ C" b) i2 F: B# J5 {% c- v kj=kj/kx;6 n& [6 u$ U: u! l2 `7 B h0=3*ht; |/ h- u1 z1 Z1 n3 {( I/ @h1=2*(1+1/rou)*q/pi;2 n* P5 A5 A/ I1 P8 ^ h2=2*(kd-1); + O4 e- A/ _5 a& t4 Uh3=1/bm*(mu0*p*q/(2*pi))^0.5; ; }9 ~2 q$ ?( `8 V9 gh4=1/kj*(mu0*p*q/(2*pi))^0.5;3 _. P2 g, z" z1 C4 z6 E! V hg=2*pi*kr*ht; . o: l2 T- N0 R! d/ bw1=gaf/bm*(pi*p*q*mu0/2)^0.5; ! [ s1 w) \3 xw1p=8*g2*kd*w1; ( @. ]: x% {. i3 c( Kw2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;! M4 ?# l8 o& y) _ w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;2 Z) U" z7 g V! H- W) r( Q w3pp=hg*h3; - s$ S7 k' E3 Z0 O2 Ww3=w3p+w3pp;# v+ A" ]( }) a$ F( n w4=2*kd*gacu/kj*(2*pi*p*q*mu0); 9 j' H- f- L- c5 V- P+ w- Dw5=5*pi*ht*gat;1 u- v% [) d' I* p w6=hg*h0*gat;- ?$ J6 s7 A) y2 ^: Y3 r' U w7=hg*h1*gat;. N z4 I4 E* M( ]( b- \! m, w w8=hg*h0*gat;/ e# Q# F# S3 G w9=hg*h4*gat; ' U D% c8 Q. a. _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);5 I; f4 ^0 T' x- a/ N# B

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 15:29 , Processed in 0.373645 second(s), 99 queries .

    回顶部