求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
R=0.15;
- P8 T; x' D0 m' |& L" Qu=[0.1221,0.1786,0.1452,0.1054,0.1132]';
# I* v$ {9 I! C' ~9 _$ `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];
% Z4 [9 \7 D3 g0 _& wE=E/10000;
c=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
  e. [! Y7 c: Naeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
, w0 g4 t4 M6 A/ gub=[1,1,1,1,1]';
# Y* J0 h3 k+ S# L% rx0=[0.02,0.03,0.05,0.11,0.12]';
+ C9 b/ y2 ?2 i# _, S[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)

  ?0 ^) b8 p1 f) q- C& h! ]8 Bfunction f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
0 h# G/ M5 _: p, p* E7 R( i%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
  n/ O9 E, \, J& F+ n) Iend
* m2 M6 q- Q- o/ e" x" W9 N& x
0 @1 X  B! t6 K, s( m; sfunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
0 Q% U; A) O# ?- Q2 S& @%UNTITLED4 Summary of this function goes here
0 |- m* B$ ?: \- x8 T%   Detailed explanation goes here
* w9 B; _3 D  xC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
/ h* @9 a' L5 V& r- y# P2 }Ceq=[];
1 ], Z8 j  l# z/ Gend
7 v, j3 j7 D2 [# g4 l  R
& v* U% R: w, }5 ~
. I/ ]( Q" x" f9 x# w! `; E运行结果：
Solver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).

x =
) H" B4 s. c, c+ j" S: f   -0.0175
    0.4875
. p: o2 J& N* P* ^$ s   -0.0175
    0.3074
( S# w7 x1 ?4 r. A% e3 _1 C    0.2328
% x( O6 A' A+ C# y$ h7 a1 a
1 J* _2 K4 p+ R" Ifval =
   -0.0390

exitflag =
     0
5 A  D1 q# {: O, A8 F  M
output =
9 f9 p) o- y( T  p$ Z4 k# s4 w! ^4 F: c         iterations: 51
          funcCount: 505
) X$ A+ x* q6 {6 P$ c; I4 S       lssteplength: 0.1250
( b* [6 |$ {; M8 t4 o1 N% K3 E           stepsize: 0.0812
' _) w1 R" T9 _* R5 s$ A          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
' Z8 T3 ]4 ]4 f3 I9 A      firstorderopt: 2.4732
    constrviolation: 0.0509
" O9 w8 U! s/ H1 ^            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12
/ k: v: z0 _. B2 d( J) l6 s

. ~1 O. N7 _* Q5 [5 X+ c    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!