QQ登录

只需要一步,快速开始

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

: [. B1 }# @. X6 b

function f=fh3(n) j$ t+ x- E! i2 m* f3 i: U g2=2^0.5-1;8 p# e; `. q: Z mu0=4*pi*10^(-7);% {" \2 E. S" Q1 v gaf=7.8*10^3; : {4 ^$ _4 u! Y5 ?, C- _gacu=8.9*10^3;7 S1 i6 z C8 Z1 ?9 [: e gat=4.5*10^3;7 P3 N6 p# J" c$ v a gaa=2.7*10^3; 9 Y# A* B# c& k- H" ^$ U- Pp=2.7*10^4;, G n( \0 [" o6 K) Q+ J: n( Y q=3*10^(-5);. I4 H, l* E: @& d rou=0.7; . G1 N x- B1 k. C2 k" f# \: ekb=0.9; 0 O; c+ o2 T/ |% M) N z%j0=8*10^6; / P' r$ ~" ]2 d/ `7 ^& K" W( G$ }j0=2*10^6; 4 C' w0 t1 j4 Q! _kd=3;( b, W0 Y' D: P bm=2;" \" k4 f, z7 v' ?0 R6 c9 t# b2 ~, ~ ht=2*10^(-3);. r5 d+ a# y! B- ^. J |4 G mu1=10^8/1.75;: ]& ?; E7 m! L5 P- E% X kt1=2*pi*gat*ht; & E6 s5 g8 o* Y/ }kt2=2*pi*gaa*ht;0 Q6 d9 k, R8 p9 a kr=1.02; 1 W/ m& |/ Q4 L% a$ d0 bku=(mu0/2/bm)^0.5;( c9 f/ h5 c' S/ g# b6 D' x( N kx=mu0^0.5*kb*kd^2*j0;6 J. Z, D0 G, v* g. R/ {- { kp=(2*q/pi/p)^0.5;, E( D& }9 }- D* }# V kg=mu0^0.5*kb*kd*j0*g2;( ~8 A" E5 a; B8 N0 G rx=ku*kx-kp*kg; 1 H, \: E( D% @) Q3 K7 }+ nra=ku*kg;7 c: z' O7 H: I" W rc=kp*kx;) L5 f( ], K) y9 A2 ?9 E+ k kj=mu0*kb*kd*j0; 0 x+ ^' ?- g& y- W1 Skj=kj/kx; ' j# ^) b( s/ o) ]% p% c; e/ ]3 T; ^h0=3*ht; % W+ Y% t! m% m; d- E3 W' vh1=2*(1+1/rou)*q/pi; $ i: L1 }5 m2 Ch2=2*(kd-1); # A) E; P1 f4 {* }h3=1/bm*(mu0*p*q/(2*pi))^0.5;- w' w+ f' Y; ^ h4=1/kj*(mu0*p*q/(2*pi))^0.5;1 x1 P+ [9 E: | hg=2*pi*kr*ht; - [. ^8 O8 m% n) Uw1=gaf/bm*(pi*p*q*mu0/2)^0.5; ! E" g1 W; a5 h+ P# E! G0 `w1p=8*g2*kd*w1; 9 Y) p2 }0 G: k7 \( q) _& ~w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;' |0 d4 j p5 q: g. L w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;3 M3 M( k9 Z4 T0 w w3pp=hg*h3; - L% Y% W+ t6 e+ V) R4 I, e* qw3=w3p+w3pp;) D6 y( Y! {8 y$ c w4=2*kd*gacu/kj*(2*pi*p*q*mu0);2 F4 J7 n& Y5 s w5=5*pi*ht*gat;2 K `0 L2 V/ @5 v2 f% r w6=hg*h0*gat; 1 `! g/ O1 \* hw7=hg*h1*gat;$ A/ J: A, w9 [' Y1 F w8=hg*h0*gat;& b' e* k: V9 }" }4 L0 l& r; X, e0 [ w9=hg*h4*gat;4 f5 y! z1 R! w7 `/ _* 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" k- Q) Y8 i4 Z h' @8 M6 C# _5 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-18 10:49 , Processed in 1.972719 second(s), 98 queries .

    回顶部