求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
; M+ K( x. p' p. H3 w( Z1 {: Vk=0.5244;
, o- Z' I) L4 `2 OR=0.15;
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
. R* j5 C2 Q' ]8 @/ d; W5 LE=[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;
* R+ |  b( J2 Sc=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
- B& m6 b" m. j7 F  c- Daeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
; r; p, f5 S0 d) S. }x0=[0.02,0.03,0.05,0.11,0.12]';
2 F8 [( q" k" X+ b' 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)

. @0 M  l0 ~! m1 `" l/ d( t* Zfunction f = myfun( x,t,E,u,c )    %目标函数
3 h' r2 {8 b; w2 [%UNTITLED3 Summary of this function goes here
, z1 \8 L6 u& k5 v5 o" I& C%   Detailed explanation goes here
5 J: r8 O& W. i& t; Z" }f=t*sqrt(x'*E*x)-(u'-c)*x;
end

, Z( A/ V3 O7 ?& }3 @2 Z, ]0 ~8 v. s$ {function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
end
+ V  [+ V0 W; h9 `2 i- l

- V: y  A+ y9 f$ A  P运行结果：
Solver stopped prematurely.
% t. x7 L' i3 Q( ^fmincon stopped because it exceeded the function evaluation limit,
0 A9 T7 w7 k& |& zoptions.MaxFunEvals = 500 (the default value).

0 }' d9 ^& J4 ?7 P, Yx =
   -0.0175
    0.4875
   -0.0175
+ B4 u( r8 U# K" Y5 J    0.3074
    0.2328
3 l+ @, q, s: [' X0 f  O1 J5 O
fval =
) m# `  S1 Y# Q1 M5 k   -0.0390
$ @4 Q' c1 m, Z& P( z
exitflag =
* \1 m/ ]6 Q9 u6 S1 ^     0
8 {! W+ |# D* J: ?
' p/ s  ~3 ?$ J) L  E. B* ^" e* Moutput =
1 F' Y3 \# ?. W( \) x. [         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
$ q8 B- r2 n: S$ u6 `/ _           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
/ b# e- H0 v3 g2 x5 {/ q: v/ f      firstorderopt: 2.4732
    constrviolation: 0.0509
0 H2 _! t+ Y$ O& `3 {; w            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

, y. [! Q+ \. r+ S" g
    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!