求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
8 @4 o, y6 Q, B, |k=0.5244;
R=0.15;
* ^9 o  M# H2 p# w) lu=[0.1221,0.1786,0.1452,0.1054,0.1132]';
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];
E=E/10000;
9 h' C. ?* V! }' sc=0.0075;
. B+ E% U8 i4 c1 Toptions = optimset('Algorithm','active-set'); % run active-set algorithm
aeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
, g) |* a5 W0 {4 b/ f$ v0 Plb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
) R5 C! L& y' hx0=[0.02,0.03,0.05,0.11,0.12]';
1 Q. Z) `0 {3 L, g% Y$ f[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)

function f = myfun( x,t,E,u,c )    %目标函数
+ a( {6 U3 K; F& b6 J%UNTITLED3 Summary of this function goes here
6 `, |: Z; n' F%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
end

* C# R( t! c( P( rfunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
* F) I) t( V" R7 K%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
- c! T8 l) g; Q( g; b; }' h9 ]! \C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
  s9 [8 R' Z" [) e) J& ]Ceq=[];
end
$ y% G3 [/ u: R: M# i4 j, h
; V3 m7 ^3 t  y1 n& o
运行结果：
) ^6 Q2 I3 F; x- W. s; o$ y, ZSolver stopped prematurely.
$ R+ x& |4 o+ Y$ K; O2 e1 afmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).

x =
   -0.0175
    0.4875
   -0.0175
    0.3074
    0.2328
, l$ `1 V: c- E
9 f: Y' a0 T2 x% ^) Jfval =
   -0.0390

exitflag =
/ v! n  i4 A$ g. l( B: H+ \7 E7 t2 c     0

& u& v9 h. D; ?$ l! Aoutput =
( \( b  |, b6 Q5 Z" b1 M( t9 S9 L         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
1 z& @( {& l" {           stepsize: 0.0812
0 a9 B* a* d3 j+ F7 P8 |          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
    constrviolation: 0.0509
            message: [1x145 char]
. z* H2 U/ r2 R8 J7 k各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

难哦。。。。。。。。。。。。。。。。。研究研究。。。。。。。。。。。

zfq12

难哦。。。。。。。。。。。。。。。。。研究研究。。。。。。。。。。。

zhangal

不会啊，帮不上》》》》》》》》》》》》》》》》》》》》》》》》》》》》

亚瑟王的冬瓜

回复 3# zfq12
0 s9 j* z1 u4 ^% ?6 |
9 A- G3 n& V! C% P) d
    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

我把maxfunevals 改为1000000，maxiter改为100000，仍是过早结束

韩国福娃

不错不错，值得学习啊!