求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
, t& S8 a! t0 A# {! \k=0.5244;
R=0.15;
5 J8 C8 o( a1 V! x3 u- x- Q. Pu=[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];
+ X( D# ?8 ]3 IE=E/10000;
- ~* m* S$ Y! rc=0.0075;
1 q- g- C& t1 c) Noptions = optimset('Algorithm','active-set'); % run active-set algorithm
4 I7 F1 q# L1 J1 raeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
/ O4 f1 H  {' a5 q3 T7 Ox0=[0.02,0.03,0.05,0.11,0.12]';
" s: r# |+ A- h* y7 x# G# _[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)
" J" i) L1 y) Q, l# r0 o
$ Z  q( F% K" f, zfunction f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
3 w# {2 ]: l" l5 [end

4 T  {6 T+ C( S! j9 A4 mfunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
2 p# V7 ]% b* g  E%UNTITLED4 Summary of this function goes here
* i0 K( Q) [8 [9 {3 B5 t%   Detailed explanation goes here
9 L4 q+ ^. w! X4 CC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
end
- y5 ^: |1 p- v- M. O- W% h" b5 U6 s

0 g: y7 X, W6 \( O- C4 d/ Z; ]% w% v运行结果：
* @& C+ b# W/ N7 s! oSolver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).
1 L/ {1 i  @% i6 i1 y/ m# \7 ~
' U" u: `' ]( Dx =
( \- W. b$ i+ P* G) @" y3 s   -0.0175
! V. B& N# D* B$ Q- g    0.4875
   -0.0175
6 E+ e, {/ a( e; w6 E    0.3074
3 `% G7 s% {- Y* t' P" H/ C    0.2328
& |3 G' `' E' K; ^
fval =
* ~5 X2 ]6 b" o; `- V& c( w   -0.0390
* h5 [& w; q! b* C2 H! t+ _' F5 c+ H) C# a
exitflag =
     0

& h* a" r2 Q' B  s, X* x  }output =
         iterations: 51
          funcCount: 505
: C; o, I% ]# s# t       lssteplength: 0.1250
           stepsize: 0.0812
1 q/ y5 T4 p' p/ ^7 {          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
' c( U5 U4 @6 u9 R4 M2 z. ~    constrviolation: 0.0509
6 S6 Z1 d- c+ [4 E8 m& U2 s            message: [1x145 char]
/ I0 S: T0 N! d  ]0 @) R; f各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12


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

火箭发动机

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

韩国福娃

不错不错，值得学习啊!