求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
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];
( C& }  b" M* JE=E/10000;
% N: y; l/ k: |) o% O- @5 ~c=0.0075;
1 c  t3 s4 y' j9 i* }options = optimset('Algorithm','active-set'); % run active-set algorithm
aeq=[1+c,1+c,1+c,1+c,1+c];
! h( ]* S8 c8 X- p2 Z1 ^beq=1;
/ u' e9 M+ n* T5 [0 llb=[0,0,0,0,0]';
3 |0 O4 Y6 k& D8 u5 @; Fub=[1,1,1,1,1]';
4 t, S9 d* U* V" Hx0=[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)

function f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
3 h* }% j4 b' h# G%   Detailed explanation goes here
) E5 L; L2 S' A$ Y! U2 o3 nf=t*sqrt(x'*E*x)-(u'-c)*x;
end
! k9 @% Y* o3 X, ?. _- K
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
8 U# E; S8 d  Z3 S1 W+ x2 j8 Z%   Detailed explanation goes here
, {. ]2 d+ Q6 A: N1 n7 f( TC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
( h+ g! }/ m( m5 }+ G! iend
3 i% ?4 n+ H6 T8 B& c4 r
$ c, `, G- ?' g' h
运行结果：
" q: r% O4 C9 U; g  e) G2 j- I: cSolver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
( f4 |0 Z5 B; P6 L- E6 ?/ soptions.MaxFunEvals = 500 (the default value).
3 ^: J6 |# L; s3 ~7 t
+ N& l, g; ?& [. W! Y; D) U3 Ox =
- W5 ~* M9 n. f1 h7 L- d   -0.0175
    0.4875
   -0.0175
2 [: D* y5 v# R% K/ \    0.3074
0 \+ V' S3 ~6 _+ U+ c    0.2328
/ D! V# V9 f, o7 T1 l8 F
fval =
   -0.0390
1 y% i1 U& T" F! c0 d$ d0 U* s! L
+ P8 J- w3 t" x7 l3 _exitflag =
     0
  Z3 h1 [0 {7 y1 R
output =
         iterations: 51
8 m8 G* |7 |# z, h          funcCount: 505
; s9 d9 b- u7 w# z       lssteplength: 0.1250
+ u5 @0 V3 W( E5 Z           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
/ J) t) W/ N8 T3 O9 `% c, b/ z! Q      firstorderopt: 2.4732
) e. K- O$ x9 M3 h* ]1 S    constrviolation: 0.0509
/ j& i- X- p- `% T            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12


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

火箭发动机

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

韩国福娃

不错不错，值得学习啊!