QQ登录

只需要一步,快速开始

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

, C' n" h. ]# g: C4 e

function f=fh3(n)7 t" p% a) ?3 |+ G; u g2=2^0.5-1;3 S* V2 F4 F: j# I8 k5 k, N mu0=4*pi*10^(-7); 8 T0 o8 G& [, Pgaf=7.8*10^3; ( k6 R5 r+ {+ l* f) j/ q& s3 Ygacu=8.9*10^3; $ C8 Z( b# M0 v* o* q& g' L8 y1 _gat=4.5*10^3;5 D* k6 r' N; D" Q gaa=2.7*10^3; # ^- S: ~! M E7 ?" P( h1 ^p=2.7*10^4; 6 K# J& S+ g: Fq=3*10^(-5); $ i" |* u0 l1 `, [4 D) d2 drou=0.7;! z. e# k& r; W) C. t4 ^! J+ X kb=0.9; 1 | w" _8 _1 d% X1 p) e%j0=8*10^6; - W5 K( U3 B9 m4 k; h Vj0=2*10^6; " [# U! j, i1 Gkd=3;% \1 _- K( Z) U7 _, P; N bm=2; / v7 @: k7 n4 q; t. Aht=2*10^(-3); & b/ M% ^! }1 |$ Umu1=10^8/1.75; 6 M7 x) `" J) u$ p) _kt1=2*pi*gat*ht; g+ Y3 e: d) Y( o1 S) w: `0 s5 x) a, vkt2=2*pi*gaa*ht;: j. _2 j. B* o9 y kr=1.02; # G2 R, ^6 p2 {2 l* ?/ ], eku=(mu0/2/bm)^0.5; # b2 E7 J3 A, E, B4 mkx=mu0^0.5*kb*kd^2*j0;$ x5 o3 Y0 `, S! x: J' `- E* |# X kp=(2*q/pi/p)^0.5;4 S9 m7 I$ ?0 g1 P6 I kg=mu0^0.5*kb*kd*j0*g2;& C0 h2 W# z1 F1 m; Y rx=ku*kx-kp*kg; & v ?4 R) E6 T+ Y4 Ura=ku*kg;; k9 |( w4 f% S# H' N rc=kp*kx;7 X0 y* ~* W6 X2 `7 f kj=mu0*kb*kd*j0;/ g# U$ {( f+ W, E kj=kj/kx;3 x* x* a" m# J5 S. ~ o' Q: ~0 z( K9 U h0=3*ht;% E2 v; i+ H1 v2 z) J6 g$ z h1=2*(1+1/rou)*q/pi; 0 L+ h. g0 D8 v+ E1 C+ ?% K Ah2=2*(kd-1); : C1 k& ^0 L# Kh3=1/bm*(mu0*p*q/(2*pi))^0.5; * c7 o' ~3 p. H, _+ uh4=1/kj*(mu0*p*q/(2*pi))^0.5; $ s9 B! t6 h2 e e/ p Phg=2*pi*kr*ht;* e9 c: M7 I/ O. }" e3 Z* j Z w1=gaf/bm*(pi*p*q*mu0/2)^0.5;" z; y, Z; b. Q! b& I. s F w1p=8*g2*kd*w1; ; G% l! G, m6 P# Sw2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; & Y; L# p/ `6 h) dw3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;0 g4 q3 ]# ]/ s( V. ^' O7 } w3pp=hg*h3;+ ~9 [# N' {4 Q9 d( k w3=w3p+w3pp; , X) a1 X1 Z' T; |, ~- b+ ww4=2*kd*gacu/kj*(2*pi*p*q*mu0); : E0 ~2 S$ `1 Z+ {0 E2 s1 l3 Bw5=5*pi*ht*gat; A5 J( s) i B2 I# e4 d w6=hg*h0*gat; ( Y7 y% U3 R7 ?& [5 |w7=hg*h1*gat;2 y7 }- F; K+ E; }9 }$ o+ N3 B w8=hg*h0*gat;# |$ \9 X D. E% e& Y, n w9=hg*h4*gat; . {# o( o3 B1 t8 [0 r* D; q$ xf=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);" b$ ]( L4 W0 a" l% h( d2 t. 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-11 22:20 , Processed in 0.485756 second(s), 100 queries .

    回顶部