求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
R=0.15;
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
  S: i8 N$ F) J6 z  ZE=[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];
4 F% v$ W1 Q7 }: ]9 M7 T# k' t$ K( dE=E/10000;
3 k7 R, Y4 O- |8 K1 @' l4 Cc=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
. N0 `  w& k7 T/ D, {4 v" @aeq=[1+c,1+c,1+c,1+c,1+c];
( O5 m$ x0 ?$ M9 a( y7 u0 j: L' R% Jbeq=1;
$ {3 A! }' o2 @! D: g: elb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
7 T% V" K* l5 ~3 i6 {# `[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)
* v& n0 \! R, o! t8 j
function f = myfun( x,t,E,u,c )    %目标函数
5 C. U: g  X, @1 x6 }6 M%UNTITLED3 Summary of this function goes here
8 k6 L) {0 y7 m, |% t* f0 f%   Detailed explanation goes here
  O4 f+ i/ g" H) p- K4 if=t*sqrt(x'*E*x)-(u'-c)*x;
end
" c/ m, _" i1 y9 Y
  e* c7 a  n+ Q5 P' Jfunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
6 j8 d' F! N1 }1 l4 Q$ \%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);
  E/ |/ L9 r% a5 ^! @Ceq=[];
5 M5 k' m0 u- L; R5 Q: ~end
' D# j( d9 S+ j

运行结果：
Solver stopped prematurely.
( s$ A; ]- d# O  N: n$ r& ]3 Gfmincon stopped because it exceeded the function evaluation limit,
  h- H% }3 S; k; R- I1 O8 Boptions.MaxFunEvals = 500 (the default value).
2 g5 f( _0 r) J$ j: X9 j
x =
' r! G! {" u7 k   -0.0175
    0.4875
1 \% H1 [7 }; [- U5 l& l' Z   -0.0175
5 |/ y0 b, t4 ^7 _* G2 K1 y    0.3074
$ G9 n' B4 T( L8 v* y/ }2 R6 V1 q    0.2328

6 Q5 `+ w) |. x2 |fval =
  G+ y  f8 g" [" H- \9 P   -0.0390
9 j, Y: d9 E" [+ e
' p" g. m7 D  V# J+ Qexitflag =
     0
/ n. y2 C: f+ [; ?% A) J& v
output =
1 {4 w' f' _' C3 `. S) B         iterations: 51
4 W; d' T+ m1 F% Y* T$ w8 T6 J( w          funcCount: 505
       lssteplength: 0.1250
8 R7 x' F3 q0 W/ K* v           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
7 r% c& o% S( X- ~& r- ?      firstorderopt: 2.4732
( x. I/ q: @$ ]# C& x4 F% c/ W    constrviolation: 0.0509
            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

韩国福娃

不错不错，值得学习啊!

火箭发动机

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

亚瑟王的冬瓜

回复 3# zfq12
5 K! q5 T: [  L8 ?% K3 Q' A  O/ p
! N0 @1 I+ H! f, e* `8 Z3 H  ~
    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

zhangal

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

zfq12

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

zfq12

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