QQ登录

只需要一步,快速开始

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

0 m; ]2 K' z a

function f=fh3(n) , ^* ]. e; A) n( g' d* eg2=2^0.5-1; 8 h% c, x; q( `6 N# hmu0=4*pi*10^(-7);# l/ c4 g- h I* _ gaf=7.8*10^3; Z* R% F+ f& ]5 N6 ^gacu=8.9*10^3; ) X0 k7 ~# I% G' [/ V4 ugat=4.5*10^3; 1 t9 _( R5 j6 ?; B$ k Mgaa=2.7*10^3; - `8 {0 S. I/ @% Fp=2.7*10^4;2 |9 g4 B- n8 ~% x; c q=3*10^(-5); 8 @3 P0 k+ F# L F! K% arou=0.7; 5 `* i* _2 Z, B! |kb=0.9; . e# H" E2 K- c- V) ~%j0=8*10^6; @4 d* N, o: T% ~ j0=2*10^6;5 J7 [+ t% O6 H kd=3; % d0 e" t7 ]6 O- f V: u( ybm=2; * F& Z4 W$ H6 }: `ht=2*10^(-3);) E7 u7 T' R! l mu1=10^8/1.75;: R7 [6 @+ w# w, N kt1=2*pi*gat*ht; 6 N# o# B: e3 mkt2=2*pi*gaa*ht; $ D! e6 \9 @ bkr=1.02;" _, P1 J# o9 \9 Q- I% o ku=(mu0/2/bm)^0.5;$ |/ n, U: R' u* Q! c, p% D0 F( M kx=mu0^0.5*kb*kd^2*j0; , k2 a) k6 C) [$ s m5 P- \9 _9 skp=(2*q/pi/p)^0.5; . Z3 B7 j8 S- o( h! t& ^, k1 _, Ukg=mu0^0.5*kb*kd*j0*g2; $ P8 I1 c# P* U, z# _$ P7 i$ mrx=ku*kx-kp*kg;0 F+ a3 n% G/ p4 _6 ^ ra=ku*kg;& h% ~) _7 `4 p) {- q2 ^ rc=kp*kx; $ |. D& c6 B8 E }& [; s! Jkj=mu0*kb*kd*j0;: y: U! k3 z& I7 s/ s0 }6 S: G1 l kj=kj/kx;/ V) C! S: L1 g. I4 o& a/ d h0=3*ht; 6 w9 s+ N& u# ch1=2*(1+1/rou)*q/pi;5 T, C' W5 D$ m- N h2=2*(kd-1);+ }0 |8 J1 Z& M h3=1/bm*(mu0*p*q/(2*pi))^0.5;& p: {) r# k% u1 Z h4=1/kj*(mu0*p*q/(2*pi))^0.5;0 B! @ O; }& w! h Z hg=2*pi*kr*ht; 0 }( t4 S; _, p8 q. sw1=gaf/bm*(pi*p*q*mu0/2)^0.5; 3 Y* ~9 Y9 y' m6 sw1p=8*g2*kd*w1; 5 {' }7 F0 e/ s6 ?* V2 Nw2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;+ }3 p1 @( H/ T* | \3 k( U: r w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;% n/ I" u, y7 G+ t w3pp=hg*h3;' n% m8 ?' |- q w3=w3p+w3pp; ' r$ x- d( K- ~3 ~( \, Q: i" Kw4=2*kd*gacu/kj*(2*pi*p*q*mu0); ( y4 d3 [! Y0 F& z4 J/ Y" Y" R- U" Lw5=5*pi*ht*gat;7 o* S% q' {. R4 U2 c w6=hg*h0*gat;5 t! A% O+ L% K+ [/ \ w7=hg*h1*gat; ( v4 C; E& _ s: F+ _0 @; Jw8=hg*h0*gat; 5 R1 S6 S5 b8 ?* _4 fw9=hg*h4*gat;3 X; [+ F5 ~/ x+ x! V 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);- x' Y) H1 q. o4 ^

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 04:02 , Processed in 3.152441 second(s), 99 queries .

    回顶部