QQ登录

只需要一步,快速开始

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

' \8 {! Y4 D1 r: `1 A- E

function f=fh3(n)$ u' e9 {' u7 j7 u g2=2^0.5-1;& Q" s0 n' b6 k S1 y$ p mu0=4*pi*10^(-7); ; E7 x" J) ~/ `- k1 Jgaf=7.8*10^3; 9 u' n) m ^; F @gacu=8.9*10^3; * x# x- ?& q5 Wgat=4.5*10^3; ; @" G8 i8 Z7 H$ f( o) P% d) x6 h5 Wgaa=2.7*10^3; * `4 G7 | U% {* F. M$ W2 |: s0 Tp=2.7*10^4;0 M( Q6 p6 Z4 u% N. ^3 I q=3*10^(-5);0 o% z) |/ X& X% w8 x7 H+ L% n0 C rou=0.7; + a3 k6 G3 ]- L- pkb=0.9;2 S+ s7 H* R" G5 f; G8 m0 P" q; T %j0=8*10^6;- z( ?7 W# [7 ^( o9 o4 j: t j0=2*10^6; 7 e' Z* B3 ?+ a1 `% f# M& r+ w2 xkd=3; + ]% S- ]. z* E$ {1 `, `. a4 lbm=2; V0 u* I3 d. b6 X2 \" }1 G& }" p ht=2*10^(-3); , S2 t+ a6 ?3 |( F, @: ~# P+ x5 Umu1=10^8/1.75; + Q) f/ C) f% l/ D! }4 ]kt1=2*pi*gat*ht;, j- b0 A, y* C kt2=2*pi*gaa*ht; 8 H \. {1 H' l+ q$ _' a. K8 bkr=1.02;' e0 U1 ?2 B% k |" F6 U- I ku=(mu0/2/bm)^0.5;+ _$ a# m$ R- o! G kx=mu0^0.5*kb*kd^2*j0; + [& r2 r, Q# h n7 b9 n n4 ]0 ^: S1 I; c5 \kp=(2*q/pi/p)^0.5;2 [0 G+ e8 K5 T3 }" i& m kg=mu0^0.5*kb*kd*j0*g2; % L) R) s" K! R* Jrx=ku*kx-kp*kg; ' S" N0 t, G$ z3 Dra=ku*kg;( A/ Z: o& r8 r, M5 T* F rc=kp*kx;- D- ]: i. g# r5 W8 d; h kj=mu0*kb*kd*j0;# o& M2 j* i4 O( n kj=kj/kx;' N' a9 C5 A2 G" F6 B h0=3*ht; # h7 B7 d" h& U4 |7 eh1=2*(1+1/rou)*q/pi;9 l- r) y P; Z- s! Q! {8 u h2=2*(kd-1); . {# K" p5 O* mh3=1/bm*(mu0*p*q/(2*pi))^0.5;, T( n+ }6 H$ i5 }0 v h4=1/kj*(mu0*p*q/(2*pi))^0.5; ; L! h! H5 C, s8 @2 D8 H: ^hg=2*pi*kr*ht; % |+ p6 Z3 s) [0 {+ U f. aw1=gaf/bm*(pi*p*q*mu0/2)^0.5; # W/ c% E. O7 Tw1p=8*g2*kd*w1;( y* R/ ]& r5 b. ~$ F" F6 y5 W/ j, G2 ` w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;. O7 y# z4 }8 O4 s$ x w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5; 3 P( |1 W B8 u+ G) E; g5 g! sw3pp=hg*h3; K* o. }) w- }8 |4 D$ iw3=w3p+w3pp;3 M1 g& k9 ]8 V; g- | w4=2*kd*gacu/kj*(2*pi*p*q*mu0); , ]' M0 x3 {6 \w5=5*pi*ht*gat;7 V/ u" t" ?1 u% {! U w6=hg*h0*gat;* h7 a: A9 S2 x- ~6 k& ^4 Q" @; B w7=hg*h1*gat;8 g5 J& }0 l$ H w8=hg*h0*gat;& P: E: T$ G! _% d w9=hg*h4*gat;) q0 v: m) s3 D$ m 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); * c) C6 S: j! _

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-1 22:13 , Processed in 1.978248 second(s), 99 queries .

    回顶部