QQ登录

只需要一步,快速开始

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

* ` w+ ^' Z3 r: I( C4 G: C

function f=fh3(n) " d9 o5 r, V# o, D4 ? g& u9 ^ |g2=2^0.5-1; Y5 I1 w6 a8 o+ J- x" ~ mu0=4*pi*10^(-7); # s# e2 O6 z. }/ _' [+ n7 Q1 lgaf=7.8*10^3; # S* Z/ Z& l6 ?+ R5 K' Wgacu=8.9*10^3; ( M- G6 n+ W3 z( o J+ A2 bgat=4.5*10^3;3 @" s& L( A' S# ~6 f; t+ a, v0 r gaa=2.7*10^3; 4 m/ k3 ~# @$ P! {* Op=2.7*10^4;2 A+ O) i9 K: s8 i, | q=3*10^(-5); 2 L0 w$ X) c3 Erou=0.7;3 E) W' P" b) [9 O* m; w' v" C kb=0.9;) `5 t6 p9 q' `- j6 K3 N %j0=8*10^6; $ _- I" l$ D6 H5 m: a2 N; d2 Q6 J' Cj0=2*10^6; 4 Y- W- d$ D' \4 S# ^kd=3; : T# b6 s( I3 G& k) ] ebm=2; , n( |0 r3 w3 G0 ~) \. ^- Dht=2*10^(-3); " f2 u% `) u( A1 _mu1=10^8/1.75;( g, {4 G/ E" X, r$ v kt1=2*pi*gat*ht;3 y+ F5 b" S2 Z6 `* q+ W kt2=2*pi*gaa*ht;% p$ P4 H0 F, u' B9 G+ d" ?/ d W3 t kr=1.02; & U e& \2 w- M0 k" m* ~ cku=(mu0/2/bm)^0.5; " W: j$ z/ N% {+ W! [( \: Okx=mu0^0.5*kb*kd^2*j0;# b+ T3 o. k$ H; J6 U# t: C6 n& ] kp=(2*q/pi/p)^0.5;. ? t4 I% J& ~6 P; A kg=mu0^0.5*kb*kd*j0*g2; - ?9 v; E8 O5 |9 hrx=ku*kx-kp*kg;/ g! y# h8 K+ R9 A) |/ Q ra=ku*kg;3 ~) r# r) `4 e( c9 M rc=kp*kx; + }5 Y+ }1 `7 V9 v5 D% j* Zkj=mu0*kb*kd*j0;1 c2 x( w }, o kj=kj/kx; ' | w% ?9 b0 Yh0=3*ht; 7 @: p+ X6 b3 F& r8 b) jh1=2*(1+1/rou)*q/pi; - [3 K# R' H. jh2=2*(kd-1);5 Z( {* D" g: Z0 g+ p' ] h3=1/bm*(mu0*p*q/(2*pi))^0.5;& |; Q" Q6 @' e& ~5 Z2 V: w9 W h4=1/kj*(mu0*p*q/(2*pi))^0.5;7 W$ E7 J$ }3 o hg=2*pi*kr*ht; . W' ]2 B1 x% N# o0 p7 l, B" dw1=gaf/bm*(pi*p*q*mu0/2)^0.5; c) y9 i& W. A7 Iw1p=8*g2*kd*w1;; q$ S# o) Q9 _, Y. F6 d! z w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; ( r8 Y! m- [; i+ G/ Ew3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;8 ?0 N, Q$ d9 x4 t8 F1 z# f H( u# ^ w3pp=hg*h3;9 V' z: S" g. U& @' _3 F% | w3=w3p+w3pp;8 A7 u! ?0 o6 C7 x- u5 n6 d) Z w4=2*kd*gacu/kj*(2*pi*p*q*mu0);0 q0 Y# y9 B; y6 u! X {" i, p w5=5*pi*ht*gat; ( d% [3 Q% J, T( B2 d* [: Zw6=hg*h0*gat;" t ?- u) H& K/ L5 u+ b0 J w7=hg*h1*gat;& H: V' N. E" S: f5 u) H w8=hg*h0*gat; 5 a2 a& e( c% Z$ iw9=hg*h4*gat; ) Q) f8 S: s+ w2 \" uf=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); 0 b( u( P+ Z! K7 \3 X7 r: {

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 22:58 , Processed in 2.149058 second(s), 104 queries .

    回顶部