QQ登录

只需要一步,快速开始

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

2 y) h0 r. `+ g% n5 L. N# g

function f=fh3(n). y! w! w- y6 o g2=2^0.5-1;, Q& H j/ u- i9 H' ?1 ~: h1 t mu0=4*pi*10^(-7); 3 n5 o; {7 U5 w* Igaf=7.8*10^3; 7 D1 p) O7 C# y: n5 W- j1 o/ q6 Jgacu=8.9*10^3; 3 {# P. S9 w. x4 K# Xgat=4.5*10^3; ) [, Q& D( N' F, b: qgaa=2.7*10^3; 9 B# i) `6 |7 L+ `+ s' bp=2.7*10^4; Q& R- k& q/ @) q) H4 C q=3*10^(-5);) R8 [" m" H. r; B u rou=0.7; 4 e3 c9 L( B( ]- @kb=0.9; * h, Y2 n4 z; ]( t& N% F r%j0=8*10^6; 9 S M! _: @0 l0 C% Hj0=2*10^6;7 e# h* C! y! F x kd=3; " `5 |5 X! s8 \& q+ Abm=2; ; |7 b' U" G2 v- Fht=2*10^(-3); 2 Z2 a6 }: z8 N; c6 ^. Pmu1=10^8/1.75; % `2 n3 Y7 O$ d# hkt1=2*pi*gat*ht; 4 O, k7 t8 Z& P* Akt2=2*pi*gaa*ht; 6 o2 N5 R- n& w* l9 Pkr=1.02;7 t. o, I" h& d) @7 e+ o ku=(mu0/2/bm)^0.5;8 C5 V: v u) b. D kx=mu0^0.5*kb*kd^2*j0; ' N g Z4 Y) Q- q, z6 f8 Ekp=(2*q/pi/p)^0.5;5 K1 V7 n; g7 Z( _0 ~7 H% P2 d kg=mu0^0.5*kb*kd*j0*g2; . F' X* q1 u% N0 urx=ku*kx-kp*kg; . B1 B- \+ w* _2 Pra=ku*kg;3 L) `& t( y8 F @" X rc=kp*kx;) \( R; Y+ y( f* j: d( \' Z' a kj=mu0*kb*kd*j0;7 R1 z5 I9 O. M4 W9 {0 g. d kj=kj/kx; ) d2 q1 f8 S6 Z4 `h0=3*ht; ; ^$ P' n8 y4 W- zh1=2*(1+1/rou)*q/pi; 4 E: ?% k2 p6 r1 ^8 ?5 r# X) w% Xh2=2*(kd-1); : `1 W/ i( @! ?! o& v( ?h3=1/bm*(mu0*p*q/(2*pi))^0.5; 4 {8 k4 K# {! fh4=1/kj*(mu0*p*q/(2*pi))^0.5;. T( g5 T3 D# z; u ~ hg=2*pi*kr*ht; + E& p! ^2 _6 R, q7 Kw1=gaf/bm*(pi*p*q*mu0/2)^0.5;2 @: V/ b/ n0 }* c! C; x! ] w1p=8*g2*kd*w1;8 K% N. p. F; ^: {# T( Q w2=2*kd^2*gaf*(2*pi*mu0*p*q)^0.5/kj; / e& d( K% B6 | ~2 ?8 Y+ aw3p=2*g2*gacu/kj*(2*pi*p*q*mu0)^0.5;: k* R2 G `7 p" W& f; K5 @ w3pp=hg*h3; 1 C$ e0 q1 z* ]- Q: U' Hw3=w3p+w3pp; ; o# L4 ~, ? P8 i [9 G) Zw4=2*kd*gacu/kj*(2*pi*p*q*mu0);* l) T& S- R5 f9 h- T w5=5*pi*ht*gat; 4 [2 V1 i3 o- ]w6=hg*h0*gat; 1 R/ G5 i8 ], D1 e8 q; |6 aw7=hg*h1*gat;( {8 G3 @- x5 @, g: `- r w8=hg*h0*gat;2 S Y$ B: @" W, t4 @: `% ` w9=hg*h4*gat;/ H; p0 m( v; s" w1 z3 S7 s 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); ' Z* F% N+ h9 s5 d

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-13 05:03 , Processed in 0.489617 second(s), 99 queries .

    回顶部