QQ登录

只需要一步,快速开始

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

% [- [- H4 G& u4 q" C

function f=fh3(n)6 X' w" w5 r+ |% m: W g2=2^0.5-1;" @, W, ~$ _1 c, E mu0=4*pi*10^(-7);& u" u& ~3 `( G* B- n% d gaf=7.8*10^3;5 O8 m8 U3 W! W, }/ k _ gacu=8.9*10^3; 5 g. ~! _2 |7 w6 R4 r8 kgat=4.5*10^3;" E% g8 Y! O% m8 B gaa=2.7*10^3;3 ?5 W/ D+ F4 x" G. j% t p=2.7*10^4;; r9 S" T) }! f2 B% [/ u q=3*10^(-5); - k+ y j8 K/ _2 Hrou=0.7; * b& C; M8 c8 l% wkb=0.9; 3 w; H3 B5 T J1 W2 v%j0=8*10^6; C1 A) D/ e' o0 P j0=2*10^6;6 i% f+ \+ B2 I4 O8 e! u% } kd=3;: d! @3 {+ j0 d3 X1 R bm=2;- X0 \7 n1 |5 v1 V! O; x( k ht=2*10^(-3);5 _" T w; B4 v; D; v( T& U mu1=10^8/1.75;! I# ~/ h7 y! X4 h$ ?* S5 v kt1=2*pi*gat*ht;9 x% C' G& Y3 \* A2 ]2 b% Q. Y8 M& A kt2=2*pi*gaa*ht; ( p( [6 `/ p6 x! ^kr=1.02;# `; h+ ?# x6 M' I8 E ku=(mu0/2/bm)^0.5; 2 d: i4 j' x. Z; Kkx=mu0^0.5*kb*kd^2*j0; 8 g: _% p0 h8 kkp=(2*q/pi/p)^0.5; * q) _1 q( a/ W+ h( I* z5 o. Gkg=mu0^0.5*kb*kd*j0*g2; " u: q) w0 y9 I4 j) O _' Lrx=ku*kx-kp*kg;3 b, H+ A# X* r1 z6 h. q1 g( O ra=ku*kg;9 s4 {4 p4 C: m3 D0 {9 y) U7 u rc=kp*kx;% K% Z6 G' Q k kj=mu0*kb*kd*j0; ! n' j$ Q9 G. B7 B) lkj=kj/kx; 8 J. o l6 U. j2 R2 A0 @h0=3*ht;+ i/ d8 F H* W; k- O h1=2*(1+1/rou)*q/pi; ; B; U) M4 a5 ^4 Ph2=2*(kd-1); ' @% z7 Y9 w3 l! R K9 a: Sh3=1/bm*(mu0*p*q/(2*pi))^0.5;- B' j4 b9 z7 P h4=1/kj*(mu0*p*q/(2*pi))^0.5; 4 c2 e0 f* q4 r% Shg=2*pi*kr*ht;8 a3 I% i! X' y3 B6 u8 l( o w1=gaf/bm*(pi*p*q*mu0/2)^0.5; 1 [! Z, H: s2 P& b1 w' }' C( Jw1p=8*g2*kd*w1;, O0 s1 {$ ~' u8 e! j w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; ! z' R5 t9 F2 y* K( Rw3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5; ! }2 L+ w5 t4 }& j- ]& S2 f1 jw3pp=hg*h3;- }, G. j# d1 c' P- c w3=w3p+w3pp;/ _2 e* T. d9 ?, }' h: s w4=2*kd*gacu/kj*(2*pi*p*q*mu0);! c$ O9 f6 H5 e1 x/ v' D w5=5*pi*ht*gat; * D1 J6 d: p# N1 N3 Z2 Bw6=hg*h0*gat;: t& h- D1 i' o# l# u w7=hg*h1*gat;0 t* Y G& F) s' A9 y w8=hg*h0*gat; ; h$ r, [+ z( K1 v" R: @4 Ew9=hg*h4*gat;# D1 O3 b$ D! f; P6 D0 ^1 i) } A) 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);5 E0 ]/ ]# B$ p' m7 w. x

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-10 16:49 , Processed in 0.482516 second(s), 98 queries .

    回顶部