QQ登录

只需要一步,快速开始

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

9 c: S7 \( i Y" L

function f=fh3(n) * m7 ]- j! ?; r/ t' H+ c: wg2=2^0.5-1; ' b3 D; r% f8 T; A3 ^mu0=4*pi*10^(-7); 4 T: X: l/ i+ C6 ugaf=7.8*10^3; & |3 ?$ v0 w, N: t. b2 Pgacu=8.9*10^3; 5 }& R2 p( v: V/ E" H' m- j: Ngat=4.5*10^3;/ w9 D4 y! H9 M1 T; |, N4 k, s gaa=2.7*10^3;+ n7 f+ j) S) Z p=2.7*10^4; 8 M/ t) x1 l5 R | m5 D+ v, Uq=3*10^(-5);2 h( x$ Y- x/ }; | rou=0.7;, v E( L+ B @4 N kb=0.9;5 R( ?( R1 K G1 Z; [$ d0 P X %j0=8*10^6; - J& N; j- F4 i0 jj0=2*10^6; " N+ s, k1 l; {) [- y9 |2 {kd=3; ( b4 ]! Y6 h% U: T5 pbm=2; % H" H2 Z5 I9 P. p. x/ } C0 ~ht=2*10^(-3);/ j! L( Y. F9 r! k7 k mu1=10^8/1.75;" N/ G7 @7 ?) c8 C9 g( G# | kt1=2*pi*gat*ht; ?8 m$ |0 m- o5 O/ Z: | kt2=2*pi*gaa*ht;) u% `5 y0 B* u( V8 g9 D! @ kr=1.02; / F1 M! J2 H: ?% Dku=(mu0/2/bm)^0.5; ! `, w$ C+ s$ o1 P& Lkx=mu0^0.5*kb*kd^2*j0;5 n5 W% `/ B: X kp=(2*q/pi/p)^0.5; ; Y6 @' l7 y; ]/ X' Fkg=mu0^0.5*kb*kd*j0*g2; 6 k; U1 [1 l! a0 H* C- ?rx=ku*kx-kp*kg; ! w0 O+ _* e1 h# E8 t3 Nra=ku*kg;% I; H( P4 z7 {. R rc=kp*kx;3 W" c7 d, I% O W5 ^0 [ kj=mu0*kb*kd*j0; # B5 k6 D$ A+ X7 Mkj=kj/kx;! }" E. C+ }4 {: L+ D% ^ h0=3*ht; 4 i2 m: T7 `/ yh1=2*(1+1/rou)*q/pi; i7 G% N" x- {8 p$ D h2=2*(kd-1); H- e# u1 b8 e$ Z: Q h3=1/bm*(mu0*p*q/(2*pi))^0.5;0 p% O; _ T ~- I+ ~- ] h4=1/kj*(mu0*p*q/(2*pi))^0.5;0 |# \( u) A4 a3 t Y hg=2*pi*kr*ht; % I' y1 Q( r' K/ d4 P7 X# W, gw1=gaf/bm*(pi*p*q*mu0/2)^0.5; " Z; l9 Y$ x2 H$ Mw1p=8*g2*kd*w1;; n4 e- i* P. E% P' @7 N* L, _6 x w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj;5 D% @4 U( w; j$ C& a6 B; n; y. Z w3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5; ? v- b( W- S* } w3pp=hg*h3; ; z, f6 \& E {w3=w3p+w3pp; ' d* `- j( T: Xw4=2*kd*gacu/kj*(2*pi*p*q*mu0); 9 y s; J) n# T; tw5=5*pi*ht*gat;: ?0 \% ^) z6 I8 X- ]: |) X w6=hg*h0*gat;7 W( w- L4 _: T* P9 l# ^0 T w7=hg*h1*gat;! y: e. K- K ]/ z& c w8=hg*h0*gat;9 {# V8 e/ _, h! I: \ w9=hg*h4*gat;% |% A$ o/ x, @; s; o 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);! e2 ]6 p! n4 o: a' u

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-10 18:57 , Processed in 0.432728 second(s), 98 queries .

    回顶部