求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
9 O# C' ^( ~, y. Fk=0.5244;
R=0.15;
) E. E* G6 I3 f0 Zu=[0.1221,0.1786,0.1452,0.1054,0.1132]';
; E# e* ]" f$ ME=[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];
+ y6 _2 g- j6 T0 ~! XE=E/10000;
c=0.0075;
) z7 L0 h; x/ O& P* N+ ^options = optimset('Algorithm','active-set'); % run active-set algorithm
" S, S' F1 h4 b2 n4 F& M/ Z; _* Naeq=[1+c,1+c,1+c,1+c,1+c];
9 B. y0 @& B: Dbeq=1;
* l4 }  l% d2 F/ j5 jlb=[0,0,0,0,0]';
$ ]9 q2 L3 @% _; jub=[1,1,1,1,1]';
! W; n9 O' p6 gx0=[0.02,0.03,0.05,0.11,0.12]';
( ?- ?. V; `9 r2 i# O/ t[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)

4 K  p% ?; \2 W9 O1 D4 ^  ]function f = myfun( x,t,E,u,c )    %目标函数
  B+ p. G8 |& R# J( ?( `: |  c: s%UNTITLED3 Summary of this function goes here
! p( F1 ^8 ]7 I3 b. {* v' U%   Detailed explanation goes here
  y7 ?* V+ r, E2 h, r9 N* ^f=t*sqrt(x'*E*x)-(u'-c)*x;
end
: V0 r& C4 v2 @& I" V1 |/ u
9 D1 J% g5 D8 N$ |function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
6 T+ V- {# r" G: z* I7 S# h%UNTITLED4 Summary of this function goes here
2 i& o  v9 I3 \, ^, V%   Detailed explanation goes here
# E5 \; D, x! q# E/ y) A! W. w4 uC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
8 ?% R  m7 f' b, P) aCeq=[];
end


运行结果：
) ?4 j* x- i6 g3 q0 v6 U6 P5 SSolver stopped prematurely.
- a: {- |5 J9 o# Kfmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).

x =
   -0.0175
: n# N  y/ T* O/ S( ?    0.4875
   -0.0175
; Q2 ~8 y: |# @  V# m1 n' H: L5 j    0.3074
* y' p% m5 {# Y, g# [. i    0.2328

! T  N$ H4 M+ }. gfval =
6 T. j8 F6 E! c# N+ T: `   -0.0390
/ n) L) ]+ @5 s& F. H8 \
- }& ]+ Q; P: `* k) S4 Rexitflag =
     0

output =
) [  ~' V- @2 Z3 D0 ?) A6 l         iterations: 51
. g) P% v4 b& g6 [          funcCount: 505
, t2 i& k8 B6 V       lssteplength: 0.1250
  X* E5 M4 m/ m& a           stepsize: 0.0812
0 `5 y& x, M) R9 W$ E          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
/ |1 w' c# j  w    constrviolation: 0.0509
% V5 l" K$ G1 T2 n, |0 F            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

' q. T8 |2 f; {- V4 @1 g
3 s$ p6 a& e7 V5 l    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!