QQ登录

只需要一步,快速开始

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

5 B1 g J" [( Z6 a

function f=fh3(n), F0 T+ h1 N! G! J. a! ?" ` g2=2^0.5-1; & b+ z% G% O1 [7 ~! lmu0=4*pi*10^(-7);/ p0 n" Y- h- ^* |9 V- { gaf=7.8*10^3; ; u0 y# s' L+ @* Y/ Y' {4 V# cgacu=8.9*10^3;: r6 f( h/ R: p' o5 J( [ gat=4.5*10^3; : L7 l. A0 z! n# X* |' zgaa=2.7*10^3; / a: }3 B5 n( ~5 `$ c/ L8 Zp=2.7*10^4; " {: a9 v9 |4 g+ d' r) B8 Aq=3*10^(-5);+ B6 s( F6 t2 V" V rou=0.7;* f9 L+ H$ d2 o' I, A% \ kb=0.9; 6 q+ t! N7 W" C) x) o%j0=8*10^6; ! `7 ^( n+ x( W6 A \- Yj0=2*10^6;" D6 r: u) Z$ w3 J# U kd=3; % T3 d6 \8 }# g# f9 Pbm=2;/ n- _/ f7 ]9 |: U/ u; _ ht=2*10^(-3); 4 W# O' z* [( q: z. omu1=10^8/1.75; 9 t# V) B1 ~# b: } x3 Ukt1=2*pi*gat*ht; 9 Q4 b/ c& d: I8 w1 B. c3 bkt2=2*pi*gaa*ht; & V: j3 c" P* I y4 ikr=1.02;, b; v1 u% r1 J% W" X/ Y# m* H ku=(mu0/2/bm)^0.5; ; N1 F8 e5 w0 k+ |4 V! w2 Zkx=mu0^0.5*kb*kd^2*j0;- I. ~7 W7 W2 M: C1 T kp=(2*q/pi/p)^0.5;7 X/ u* U {: ]8 c5 O# _3 L kg=mu0^0.5*kb*kd*j0*g2;! m4 M- E* P. \9 W9 l) t! } rx=ku*kx-kp*kg; 2 z( M# i- n) h2 Fra=ku*kg; ( i0 |) N0 H. X& m9 xrc=kp*kx; 1 O: ]. U, E3 A. Dkj=mu0*kb*kd*j0;: h- A! F7 C0 ~$ x kj=kj/kx; 6 {" ^1 \8 w1 p. d( bh0=3*ht; 6 i1 ~) W l2 C7 N" Bh1=2*(1+1/rou)*q/pi; ; G& E: X0 D9 ~- C' }' O7 eh2=2*(kd-1);- |8 P0 f* ?+ u; O0 G h3=1/bm*(mu0*p*q/(2*pi))^0.5;) ?; e. j: o4 F0 b, A S h4=1/kj*(mu0*p*q/(2*pi))^0.5; * f: n/ B) q `1 i, T! c5 shg=2*pi*kr*ht;' ]* |! F. X1 g0 {, x- Y& ]# T w1=gaf/bm*(pi*p*q*mu0/2)^0.5;) J: i. E% W7 c$ t w1p=8*g2*kd*w1; % V! Y. a" g# ]. |. p* J! \w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; + r7 [7 w+ \5 g5 g$ {4 E9 _w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;+ S- s( o$ r5 A& C; z w3pp=hg*h3;0 t: W% f: m/ S/ k# m* l9 e# r w3=w3p+w3pp;+ J: i& E( X( H) O0 ] w4=2*kd*gacu/kj*(2*pi*p*q*mu0); : v( A9 D8 R" }: m3 O H, Gw5=5*pi*ht*gat;; k- X; r" ^) j7 `* A1 Y% e; T w6=hg*h0*gat; 5 X# d- j% ?% \- Hw7=hg*h1*gat;2 x4 \. T# K! Z9 j6 j9 u+ s$ G/ c w8=hg*h0*gat;. {& w' E6 O- H/ R* {' ^0 e w9=hg*h4*gat; 4 V6 N" z4 r" G' i6 b; T7 ef=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);$ Z# V' G6 K! B# K# o g* s @

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-17 21:16 , Processed in 0.475719 second(s), 99 queries .

    回顶部