求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
$ w0 A  ]/ S# wk=0.5244;
/ L/ f# N9 t/ hR=0.15;
5 I8 D. _. p+ E% w8 v" O0 j5 Vu=[0.1221,0.1786,0.1452,0.1054,0.1132]';
( t4 W9 A9 W4 K- s- p4 Q$ }# ^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];
E=E/10000;
* N3 O; n! l  ?- q, ]9 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;
4 W$ V; |) \  q4 Zlb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
0 L' I. ~9 z  W, \x0=[0.02,0.03,0.05,0.11,0.12]';
& k! W% O) z! B7 n[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)
* l2 Z, t& w) v) F; I
function f = myfun( x,t,E,u,c )    %目标函数
( @1 S2 M( A8 z7 k6 w%UNTITLED3 Summary of this function goes here
! h  H# r# r6 {3 S) b( l%   Detailed explanation goes here
5 M+ L, D  U$ t4 j" X4 uf=t*sqrt(x'*E*x)-(u'-c)*x;
end
0 Q& i4 c: e' w6 b! R$ P
. o. f. e0 S( [& q8 Ffunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
, ?$ S9 F, X3 n( F' a( U* O# k%   Detailed explanation goes here
2 c9 x! G% r. OC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
* p1 w  W$ D# i! L, S' ^Ceq=[];
end


运行结果：
Solver stopped prematurely.
# c+ H+ ]+ v7 n1 P1 ~2 Ufmincon stopped because it exceeded the function evaluation limit,
7 |: r1 v# j' F  X7 m( n2 toptions.MaxFunEvals = 500 (the default value).

x =
   -0.0175
; D# b. u  O; Q) E& Y/ d# C    0.4875
   -0.0175
6 a; C/ l! N  W3 t    0.3074
    0.2328

fval =
) y. a3 W8 T( a4 T  i   -0.0390

' H2 g3 e- {) H  c: zexitflag =
9 e# c2 M0 E' q0 \5 ?     0
! r; c* x6 Z  [
output =
) V* v7 P0 z' ^& ]0 K& f         iterations: 51
9 ]) c& Y+ x; n" R: b; e          funcCount: 505
       lssteplength: 0.1250
           stepsize: 0.0812
# `* l2 ?6 v. T9 W! Q" }- ^  g/ w          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
    constrviolation: 0.0509
3 I9 X) E) |3 ^. H2 g5 b5 e% t+ W            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12
' h1 l, j1 [* y
' n# |, d8 R8 z/ p
    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!