QQ登录

只需要一步,快速开始

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

; Q7 y& P @( P3 ` u/ r

function f=fh3(n)( X1 T; @ y) Z1 w% ? g2=2^0.5-1; 0 i) G0 ]+ g- ~mu0=4*pi*10^(-7);' G* M6 S) K/ Q" x2 U4 H gaf=7.8*10^3; ( o0 R* z5 }9 ~gacu=8.9*10^3; 4 h. S6 `( M0 x6 Z2 W4 V. Cgat=4.5*10^3;$ f7 \8 M6 z& \/ H8 X( e gaa=2.7*10^3;( f" y9 r8 U4 x& C; ?' p p=2.7*10^4;" {( Y6 T$ H% u7 H q=3*10^(-5);8 l/ d: w7 |4 [6 {6 C: T. z rou=0.7; 0 z! d4 a# C9 C" v8 ]6 bkb=0.9; % c* y( ` g9 p%j0=8*10^6;/ [7 _4 q6 c5 \/ G, s8 b+ x { j0=2*10^6;9 j: b0 q8 q! p) N7 T kd=3;$ D. ]2 l# _2 t7 t7 y3 k bm=2; 2 R( Q6 A6 D% J% b. p3 |' i: `9 R% W9 dht=2*10^(-3); 7 M; K& v" l' j: _ ?mu1=10^8/1.75; " N" O# H7 s0 q6 O- ?* c+ Skt1=2*pi*gat*ht; 0 x! i1 W3 c& T8 {' Qkt2=2*pi*gaa*ht;! J5 Q* b/ e& B& Q6 K5 g1 i" a kr=1.02;: L& x0 Y) `: C5 i0 T5 Y* G ku=(mu0/2/bm)^0.5;8 |! j1 ]6 t+ s D0 O& Q2 R$ l+ Q kx=mu0^0.5*kb*kd^2*j0; 5 G& v8 @+ ?: l7 s0 R9 c8 h2 Mkp=(2*q/pi/p)^0.5;" O4 Q: L- H" G1 h) l kg=mu0^0.5*kb*kd*j0*g2;+ Y3 P+ d# n8 T# Z: z) U rx=ku*kx-kp*kg;( Z0 F) M( `4 x* @! q0 v- [/ V ra=ku*kg;/ c/ u! g' ~; j rc=kp*kx;/ R4 c: U2 w8 Q kj=mu0*kb*kd*j0;% w: O7 S) B: Z" X! Z4 X kj=kj/kx; " I6 D+ A4 D& M7 ~, j& N. a; Nh0=3*ht; + @9 I! n( y7 a5 f5 _' o8 }2 u2 w& rh1=2*(1+1/rou)*q/pi; 2 X5 s& [: u M! H3 |; Ah2=2*(kd-1); # o% @! {, U, n4 j1 d% E+ ^h3=1/bm*(mu0*p*q/(2*pi))^0.5;, @" B! P& Z4 W' F8 i" `5 v h4=1/kj*(mu0*p*q/(2*pi))^0.5;; ]# ]6 U# B2 D7 E: Z# x hg=2*pi*kr*ht; ' W; F7 A& }1 G3 g5 ew1=gaf/bm*(pi*p*q*mu0/2)^0.5; 4 [4 o T1 A4 D& yw1p=8*g2*kd*w1; . D% @; I% p4 L8 h1 W4 o# [w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;: r9 M' n, |. }% n3 Z) x w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;& d/ X; ^5 `( n, w, H& W7 Q w3pp=hg*h3; 0 H3 ~" O/ Y8 E8 i" t# N- \7 rw3=w3p+w3pp; 2 S4 d5 A+ m8 `w4=2*kd*gacu/kj*(2*pi*p*q*mu0);. x, ?/ V' Y% }2 J w5=5*pi*ht*gat; 4 L2 y- N" s" U$ D3 A ^# C9 Jw6=hg*h0*gat;4 n" @. {( f9 U w7=hg*h1*gat;: H9 ?/ L$ t/ s& c6 s w8=hg*h0*gat; ; G% b, Y, J2 o. P% @6 ]( [w9=hg*h4*gat;1 M6 y9 ?& v3 d5 P0 u0 ? 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); * }1 u J2 g5 E

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-9-13 20:12 , Processed in 1.399153 second(s), 98 queries .

    回顶部