求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
, K* I3 L+ `; ?7 g& I; q5 v) Ck=0.5244;
5 m2 L7 T6 L9 l* ]) _R=0.15;
u=[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];
6 [7 \5 {$ a2 t  h: V1 xE=E/10000;
c=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
aeq=[1+c,1+c,1+c,1+c,1+c];
! g; a0 w4 m- L3 p) U) a4 abeq=1;
lb=[0,0,0,0,0]';
' J- R# P  D2 _2 a& w. G$ pub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
[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)
# ~' V8 @! ?5 Z" g/ ?7 b' ]
1 @; g; ^8 E4 V! r9 Rfunction f = myfun( x,t,E,u,c )    %目标函数
$ r8 R. M& E2 E9 Y6 ~2 ]8 ]%UNTITLED3 Summary of this function goes here
3 L: M: x2 T8 ^) F- b* g%   Detailed explanation goes here
# ^) K8 l4 G5 A) ]f=t*sqrt(x'*E*x)-(u'-c)*x;
end
5 u5 ?% Z4 B$ P# ^5 Y. x; L
, A+ [2 D+ @, v0 Wfunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
* e9 `# X( r  o8 k2 g+ {! R# D%   Detailed explanation goes here
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
. O4 f6 `0 {' Z7 [+ NCeq=[];
end


5 m& ~; e! c% V$ E/ a9 _运行结果：
" m7 D% [. \) k+ H& q% u1 eSolver stopped prematurely.
/ o' w5 R! `+ @( Z# ifmincon stopped because it exceeded the function evaluation limit,
: R2 D! z% t" Woptions.MaxFunEvals = 500 (the default value).
* z6 w, ]: b0 }, z! }' p
x =
   -0.0175
    0.4875
' T( P1 e% ~% e( n- p* X: S   -0.0175
    0.3074
+ z/ `! y" I3 ~. z    0.2328
7 l5 u) ^" i/ V3 D! f% o
fval =
   -0.0390
$ s6 o* ^8 }& P. A
exitflag =
     0
7 [9 r/ Q8 f; {# y8 H' p* P, Z
output =
- p$ F% Y1 }* M+ o         iterations: 51
          funcCount: 505
" T6 W: v2 _5 E8 R  M  c; ^       lssteplength: 0.1250
3 b' m& o! k+ h7 T! W           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
! L5 F/ R0 |; F$ s+ Y      firstorderopt: 2.4732
    constrviolation: 0.0509
; l) M1 }; ]0 K, e* ^) t8 ]            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12
: u! D3 |# Z9 N, q
; y" f. z6 o4 R+ n3 K! h( Q- X
  L* [' T7 L# R( g& F; y    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!