求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
  @" g, Y( a: U! Y. q, `: L9 Tk=0.5244;
R=0.15;
: P/ K# x1 Q, U6 X& u0 yu=[0.1221,0.1786,0.1452,0.1054,0.1132]';
6 M  e5 _+ m+ k6 m# 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;
( G$ w, W. m5 v3 r4 zc=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
aeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
* @8 p  p4 p, u9 gub=[1,1,1,1,1]';
) m" p) H, f9 J9 }x0=[0.02,0.03,0.05,0.11,0.12]';
7 n+ }- ^+ t/ _4 s+ Q( z5 Q9 W5 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)
1 W) x% [% Q0 b4 }1 V
function f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
1 y. m! `& R- ^6 V%   Detailed explanation goes here
  ~5 h' ^; p& k+ f/ |6 `% h6 R) Ef=t*sqrt(x'*E*x)-(u'-c)*x;
end
  Z. _1 z& F: O* u( c6 ~
7 m% }; ?1 P& c# L; R# Ffunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
& M9 Z: W3 t4 D! ^* d1 O- E+ L& r%UNTITLED4 Summary of this function goes here
! f- h. t' u) ]9 t2 v0 j8 _4 s%   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
9 ~, P/ t) P% t: E
! d* h. |8 H( X- o' k1 ]( n
运行结果：
Solver stopped prematurely.
0 S7 ]+ E) r0 y/ @fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).
: d  o7 ?% {! k$ M' P% |
! v; _! k  A: [# Vx =
   -0.0175
$ l5 C  S0 [# Y0 Q    0.4875
   -0.0175
' |& j- R" W  j3 ^. H* H5 u% ~    0.3074
7 O0 X8 A9 y$ F( x. J( j: y    0.2328

fval =
2 j1 d0 S( ]4 E$ u6 L; A   -0.0390
" G$ t. l$ u5 r  S% U
, Z3 B  O6 Z% @( \exitflag =
     0

output =
         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
; c& A/ p8 l' h. q" l0 ^      firstorderopt: 2.4732
2 t" H. ?& x2 \3 l& r' }3 J% O    constrviolation: 0.0509
            message: [1x145 char]
, h: N6 d$ l, \7 k8 c各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12
9 h% H' V/ J" l) O# L/ F$ d

    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!