求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
$ f: X3 y, S: T. R9 n0 c. A7 lR=0.15;
u=[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];
5 L( M" a4 X4 t+ V- cE=E/10000;
c=0.0075;
' R! {9 p( e- I2 `8 L/ coptions = optimset('Algorithm','active-set'); % run active-set algorithm
aeq=[1+c,1+c,1+c,1+c,1+c];
" h6 \5 w9 L$ |8 j& m' k. Ebeq=1;
lb=[0,0,0,0,0]';
! v8 S3 w  E" p/ ?! Eub=[1,1,1,1,1]';
+ K4 a& C' t8 p7 _- Y1 _x0=[0.02,0.03,0.05,0.11,0.12]';
$ D; a, ~3 h, m* [# a[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 V/ k) I" c# q( g$ L  `
+ m3 ?4 F! n- k# Lfunction f = myfun( x,t,E,u,c )    %目标函数
. F5 j$ i$ W: P% {3 D& s4 Y/ b%UNTITLED3 Summary of this function goes here
, E& T  c2 M5 f%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
end
' C& R4 y, D* R" i" C( e- s7 h( M
7 l2 X* a& O# c4 gfunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
8 H: Z! v! h4 B, t1 c: iC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
end

/ O8 r# E- A* f/ ^
' @: P3 O5 J. A0 N7 u# M5 c! I运行结果：
# ~8 l; n$ l/ m) R; J( n7 S% wSolver stopped prematurely.
4 ]. l/ [" v3 Jfmincon stopped because it exceeded the function evaluation limit,
4 m: |8 }' F4 y  W3 {4 q" `options.MaxFunEvals = 500 (the default value).

x =
   -0.0175
    0.4875
; t" U; R8 s. P' n  m   -0.0175
% w6 i/ n8 S; C    0.3074
    0.2328

( c2 x- }0 O+ K8 I/ cfval =
   -0.0390
# g" Q) P* y1 W) `
exitflag =
     0

output =
         iterations: 51
) m" x# l5 P2 w/ H* U          funcCount: 505
       lssteplength: 0.1250
* R; Z! g0 I' D' F7 A: h( m           stepsize: 0.0812
  {( i3 s" Y( w. z$ ]          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
7 E" e( X' H0 m1 s* h7 ^    constrviolation: 0.0509
0 w4 s% ~2 q& G2 O3 A            message: [1x145 char]
0 W) R1 ]" {: A: A: b0 d. i; P' V各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

3 e- F' ?8 P$ ]! I( V
8 q, I- K6 Y* v: f% ~% w6 v' ?4 d    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!