QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5313|回复: 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 e2 W6 ? N, R" w/ b3 y( P7 J: p

function f=fh3(n)2 _3 ^2 K6 `. A% t9 o# |2 E g2=2^0.5-1;" _, a+ i/ \6 a) p n% { mu0=4*pi*10^(-7);: }, Z" B7 B6 `. N+ J6 g: i6 R gaf=7.8*10^3; 0 Y- v' B( W7 B- ]) s$ m4 vgacu=8.9*10^3;8 k+ Z/ j) D) r& r& `; Y; M) N5 N gat=4.5*10^3;2 g2 X% k" _4 E gaa=2.7*10^3;8 B) a" ?" @# T* g# L9 z p=2.7*10^4;5 ^/ \* H2 u3 E. s$ m, \8 w4 R q=3*10^(-5); 1 F6 A/ M% O" P* [; erou=0.7; 0 y+ g" d* M5 c/ skb=0.9; # B6 r$ _7 [% k4 F- @%j0=8*10^6;7 C( s! o9 w$ Y9 r0 l" e" m j0=2*10^6; ! B! V, `, b y: qkd=3; 9 n7 i, l, _, [" A; L& Pbm=2; ) j6 M; R. A( w5 uht=2*10^(-3);6 L" ]& O8 u* X5 \0 p" E7 y4 \ mu1=10^8/1.75;, }: E, a. D ]. U4 s kt1=2*pi*gat*ht;' X& Z. u! \1 c: [ kt2=2*pi*gaa*ht; . v) E7 N9 }* _, V! M1 Xkr=1.02; * I( C& b5 E/ t( _. b8 ]ku=(mu0/2/bm)^0.5;0 j: b, Z! ~4 ^" M% H% z& u kx=mu0^0.5*kb*kd^2*j0; @7 A8 l* ^- Z5 ]/ Jkp=(2*q/pi/p)^0.5; $ M# @, ], J) a4 z! gkg=mu0^0.5*kb*kd*j0*g2;$ L& x! A- X3 B( f rx=ku*kx-kp*kg; ( A6 V* p* ]5 ?ra=ku*kg;0 C' d# }: Q X$ @4 s9 ~ y rc=kp*kx; l, b/ \) ^$ a2 ], A: ?" F+ H1 X8 \ kj=mu0*kb*kd*j0; 3 y) Z# r: O# t8 Z! r0 X8 Q3 Zkj=kj/kx;2 j# s9 p( C, T. m. I h0=3*ht; 9 P; D, o" f( j3 A4 l- G: |h1=2*(1+1/rou)*q/pi; 9 l3 m, l3 x" Y. yh2=2*(kd-1); 9 u/ B9 Y1 R+ F. jh3=1/bm*(mu0*p*q/(2*pi))^0.5;5 X, z( M- G* N. |3 D' t! I5 H h4=1/kj*(mu0*p*q/(2*pi))^0.5;7 y$ h4 m6 X6 y hg=2*pi*kr*ht; 0 f- A# i3 d3 b- A7 D$ Rw1=gaf/bm*(pi*p*q*mu0/2)^0.5;( X. l0 |1 c8 E w1p=8*g2*kd*w1;8 W! k, f8 g; |( i) c; }2 {* w2 B* h w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;% E3 X( L: \# X4 }) ` w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;9 J! N, I9 Y6 W+ [ w3pp=hg*h3; ' O- @$ w$ _6 j8 _! ]w3=w3p+w3pp;# ^& v0 @1 m5 O' G w4=2*kd*gacu/kj*(2*pi*p*q*mu0); 6 A6 _9 w. l: J6 F( |+ Uw5=5*pi*ht*gat; ! Y/ t0 c3 t8 `; Y4 Jw6=hg*h0*gat;9 m8 I" P8 u6 l; T w7=hg*h1*gat;1 g$ X6 {2 |! W1 y( [ w8=hg*h0*gat;$ V1 K8 ~$ l: P w9=hg*h4*gat;6 K- I" X4 g5 P" ?- I: O7 _( p' K 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);! W# d1 G8 v) z( Q

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, 2025-8-2 11:45 , Processed in 0.634882 second(s), 98 queries .

    回顶部