求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
R=0.15;
1 e; l- E+ @4 \. r5 ~6 Z- mu=[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];
E=E/10000;
, n2 x" e* }  _4 }& k  i/ K; @c=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
$ x* F  [% m* q* `, {aeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
/ \& U- C- A7 i5 j  dx0=[0.02,0.03,0.05,0.11,0.12]';
[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)
9 Q6 D  B; H( j, k& J
function f = myfun( x,t,E,u,c )    %目标函数
) @; E4 q$ M! F) a% ^; `6 L%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
. k, p9 [' _0 ?* h% l$ {) xf=t*sqrt(x'*E*x)-(u'-c)*x;
end
8 z% L; h. s, h- ]$ z# O5 {
$ ~/ q: f8 J! z" F, u6 b' K2 ?, {function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
% U4 @" N) Y" B7 u%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
$ M. F& Q4 b8 [7 EC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
! y: J; v3 u( u& j* c8 }Ceq=[];
/ N9 I# g1 G$ cend
. }1 Y) h/ A; x# K9 _, F
( b0 x4 x, F0 l
2 A/ J& P# T9 j运行结果：
Solver stopped prematurely.
+ _1 w& Z9 x+ Lfmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).
# O6 b7 J$ s2 {2 H
x =
# t  V( C& L8 D5 C4 H7 [   -0.0175
6 ^" F# d5 e% n$ Z2 q0 H7 O    0.4875
   -0.0175
    0.3074
2 x! @- w8 q) m7 U( y1 D- P    0.2328

fval =
6 o, V% v5 ^% m# d! q# e0 D   -0.0390
# z! c, l: @( ]! _6 ^3 U5 s
9 j- o" @; O7 [. Iexitflag =
( a5 z- R! ~3 \( u     0

output =
         iterations: 51
          funcCount: 505
6 Y# _- H5 `3 b- Y" r  J) r# h       lssteplength: 0.1250
3 c; ^4 Y3 ^) u) {1 v9 L, f6 d           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
0 O5 r% r, m; C( ~) s, D    constrviolation: 0.0509
" m" O8 Q' F' n3 S: D            message: [1x145 char]
  d7 m8 l  K+ U1 w  B9 W. Z7 j各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12
4 @$ k$ a6 p, j" y* Z
! b6 j8 N- ~5 z/ M- h' n8 h
1 @. `8 \5 \" d1 l' V5 ^    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!