QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 9319|回复: 7
打印 上一主题 下一主题

求助: fmincon的参数问题

[复制链接]
字体大小: 正常 放大

3

主题

3

听众

20

积分

升级  15.79%

该用户从未签到

自我介绍
对新事物永远热爱
跳转到指定楼层
1#
发表于 2010-5-20 10:10 |只看该作者 |正序浏览
|招呼Ta 关注Ta
t=2.3263;6 F# I. V" N5 h/ R! q+ d
k=0.5244;
( ]7 r& O; B) \9 ]" P3 [) k8 IR=0.15;# D; B) V) V5 t# N4 v9 w$ q' o6 O
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';6 l# H( d+ S" q+ Y
E=[1234,864,-617,-246,987;864,642,586,-198,-815;-617,586,1433,469,-963;-246,-198,469,346,25;987,-815,-963,25,1827];
' O0 V( o' G0 {3 o; p% r( ^5 e8 D+ zE=E/10000;
* M- g" o" l4 q2 S' l9 P' d- }c=0.0075;  a- ?5 ?# O8 G% ^
options = optimset('Algorithm','active-set'); % run active-set algorithm
6 ?8 Q% J6 k! }8 Waeq=[1+c,1+c,1+c,1+c,1+c];, U0 q% J1 r5 C; i
beq=1;
9 {$ p# v9 M/ jlb=[0,0,0,0,0]';
* [. ^) o: h2 I* \2 Dub=[1,1,1,1,1]';
  I7 q* o2 T& z, w- I$ K- V: L3 \. L6 bx0=[0.02,0.03,0.05,0.11,0.12]';$ Z3 b5 Z+ ~' {1 t' |
[x,fval,exitflag,output] = fmincon(@(x) myfun( x,t,E,u,c),x0,[],[],aeq,beq,lb,ub,@(x) mycon( x,t,k,E,u,c,R),options)
/ v. q  l) ]9 i
- ~2 r/ z8 R0 ~" Z( n0 nfunction f = myfun( x,t,E,u,c )    %目标函数( o/ C# s0 J2 V" e
%UNTITLED3 Summary of this function goes here
8 \' q/ J5 v8 M# w: h%   Detailed explanation goes here0 `0 t. m2 v7 y6 |+ ]
f=t*sqrt(x'*E*x)-(u'-c)*x;
$ f5 G9 x9 r8 ]end3 A: Y8 h- l3 ~% P+ D% v! N
% \( X6 p0 M: t, D( E2 y! n4 m
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
  b) z0 L1 U+ Q+ [%UNTITLED4 Summary of this function goes here
6 D) U7 a# L. z$ a5 K3 E) q%   Detailed explanation goes here3 n! N# Z: g- D( i
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);5 l" @9 L' }5 m2 b, t
Ceq=[];
( ^8 J4 q6 j" g* h# ~end$ h9 S4 R: P* p$ ]$ t! g( ]

7 n$ V! K+ R$ h/ x3 t+ B6 K1 ^# r! h' H3 ^! t' G% V& g
运行结果:9 z' S, D( T1 N
Solver stopped prematurely.) h3 X/ V# [2 r  J- }( u: f1 m: g
fmincon stopped because it exceeded the function evaluation limit,
  v1 h+ ^; g1 i1 M( Voptions.MaxFunEvals = 500 (the default value).
' B1 }; B0 s' J: s) B# i
/ `* y* t# ]* ~; N- Hx =
' ]" H9 p% T9 Q  n5 S3 e" s   -0.0175
1 B9 `. i4 u3 g    0.4875
- b& P% c' R! U+ v( v2 M( ?$ M5 T   -0.0175
+ e2 ~: {/ G+ }* m    0.3074
! g$ l  _$ `  M$ l! C2 n    0.2328
9 I& s  |+ |7 j- G& A; c: _( \% |
fval =
4 [! N2 R9 F, H( o( F& v& A   -0.03900 y5 r( L5 l+ w" r3 r8 N
4 I. k7 \5 e7 H% y" \. z7 G6 N+ W$ s
exitflag =
' v7 Q4 K& }* e  x& D     0! |9 E. K) H- [1 g
$ a' n6 V% K- w' r1 c" T  ^
output = % n. Z$ o+ C* a7 z3 I3 T5 c
         iterations: 51+ A; A; H' G2 q# ~- }9 Z
          funcCount: 505
7 Q. r, m1 R+ E       lssteplength: 0.12502 X; N8 X& _" x" Q( @7 M2 t) d
           stepsize: 0.08124 C7 H8 I3 ^2 I; F$ S
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'! L4 W2 [8 E/ |
      firstorderopt: 2.4732
3 _- g  C% }: V% i; ?    constrviolation: 0.0509
# A- `8 d( }! e3 {            message: [1x145 char]
  J4 O6 k/ v( A$ O各位大虾,程序如上所述,为什么函数提前就结束了,是不是要改默认参数设置,小弟也试过,但是捣鼓不出来,希望各位高人不吝赐教,小弟先拜谢了  Orz
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

0

主题

3

听众

7

积分

升级  2.11%

该用户从未签到

自我介绍
888888
回复

使用道具 举报

alair002        
头像被屏蔽

1

主题

4

听众

328

积分

升级  9.33%

  • TA的每日心情
    擦汗
    2012-2-6 07:40
  • 签到天数: 6 天

    [LV.2]偶尔看看I

    提示: 作者被禁止或删除 内容自动屏蔽
    回复

    使用道具 举报

    0

    主题

    0

    听众

    2

    积分

    升级  40%

    该用户从未签到

    回复

    使用道具 举报

    3

    主题

    3

    听众

    20

    积分

    升级  15.79%

    该用户从未签到

    自我介绍
    对新事物永远热爱
    回复 3# zfq12
    0 v6 G, `! }( c- e3 f7 x' L  }& F6 r7 `# I3 P! W4 B, N$ S$ X

    0 J& V2 P$ B0 {- S6 Z" C    我也试过把lb,ub上下界去掉,结果程序跑了一个上午都还在busy。。。不知何故啊
    回复

    使用道具 举报

    zhangal 实名认证       

    5

    主题

    3

    听众

    45

    积分

    升级  42.11%

    该用户从未签到

    自我介绍
    冰雪聪明 热情大方 诚心诚意
    不会啊,帮不上》》》》》》》》》》》》》》》》》》》》》》》》》》》》
    回复

    使用道具 举报

    zfq12 实名认证       

    2

    主题

    6

    听众

    981

    积分

    升级  95.25%

    该用户从未签到

    自我介绍
    我很喜欢数学,希望在数学中国能结识更多的朋友。。。。。

    新人进步奖

    群组西南大学建模组

    群组数学建摸协会

    群组Matlab讨论组

    回复

    使用道具 举报

    zfq12 实名认证       

    2

    主题

    6

    听众

    981

    积分

    升级  95.25%

    该用户从未签到

    自我介绍
    我很喜欢数学,希望在数学中国能结识更多的朋友。。。。。

    新人进步奖

    群组西南大学建模组

    群组数学建摸协会

    群组Matlab讨论组

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-19 19:04 , Processed in 0.512299 second(s), 94 queries .

    回顶部