求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
2 V% I4 q4 O  q3 Q. k: ~& vk=0.5244;
* V7 [- N+ {+ Z8 Y3 m, g* BR=0.15;
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
7 ]$ U3 l  n+ [# J# S* C4 oE=[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];
3 q1 c; [7 z, z6 ?7 f. CE=E/10000;
: B- z- ^( }, P9 }  n' l6 T8 O" ]c=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]';
* |7 r. ]+ u/ o5 `3 i: J" H: q/ xub=[1,1,1,1,1]';
( w7 I! v. Y* ux0=[0.02,0.03,0.05,0.11,0.12]';
% x' s  y2 i; V! w[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 )    %目标函数
%UNTITLED3 Summary of this function goes here
; }' Z/ h$ ]* V. t, D1 m4 n% S( k%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
end

function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
) c, m! M3 Q5 q* I- J2 r  x%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);
& L4 D  }% b8 [8 ]Ceq=[];
end


运行结果：
- c7 r& ]8 _& {" j2 [* JSolver stopped prematurely.
$ O2 \! M9 Q8 i- A% I" O, Jfmincon stopped because it exceeded the function evaluation limit,
9 w" e" a% [; K  X0 Hoptions.MaxFunEvals = 500 (the default value).
+ ~' e4 u: l7 ^4 ^% a) D
! ^8 `7 g2 k7 ~( }x =
. O0 j7 g6 q0 x# }& W  k5 U   -0.0175
    0.4875
   -0.0175
! W# b9 }6 u+ y. v" V5 T4 o    0.3074
5 l# H% P8 u/ P' o  r/ m    0.2328

* f& T6 g8 S3 }: Qfval =
   -0.0390

" r* V: ^$ h+ gexitflag =
     0
* [. \2 r1 e4 k! p) Q
- e. u5 _' _$ W0 B! t. z  a+ Koutput =
         iterations: 51
8 s& k; L# E" Q4 ]8 M: T- i          funcCount: 505
       lssteplength: 0.1250
9 V2 f3 \, a! p; S           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
! W: _3 c( J) X8 Z. v! U" ~      firstorderopt: 2.4732
    constrviolation: 0.0509
/ k( m: c9 p, |- r3 y6 H9 l, _- l            message: [1x145 char]
- @2 ~8 J, g* ~* x! I2 I9 p各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12
  l4 _* p8 ]% V# z# U
. @: X: ~1 l5 W# N0 k" F
0 V* i, o3 i1 L) c" \9 B$ G    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!