求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
! t& L  G" ~# nk=0.5244;
R=0.15;
# h$ X" \  i) W4 ^' u* o5 ku=[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;
c=0.0075;
/ e, C9 g6 Y: h# toptions = optimset('Algorithm','active-set'); % run active-set algorithm
1 X# m  T) u; k2 F6 f- Q; naeq=[1+c,1+c,1+c,1+c,1+c];
) `: S- `& v3 o3 Q1 V; A$ k* sbeq=1;
) x& D  v7 y: b' j( u# h% ulb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
5 N% F1 ~$ D6 L[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 `6 ^5 L' \( L8 ^& O+ S- Z, x
/ I4 r$ y2 I8 S% ]/ O! B8 a! sfunction f = myfun( x,t,E,u,c )    %目标函数
* I8 Q: n9 c# g( q%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
end
9 W7 G: [' l( {; i
3 w$ Q) R0 @# U3 \# ~0 X( b: K" Rfunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
8 ]% X4 J/ k" `%UNTITLED4 Summary of this function goes here
4 A/ Q/ B0 V! x4 s0 z%   Detailed explanation goes here
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
" g/ i9 Z/ ^1 I8 m+ l) NCeq=[];
end

  b' f! }! y$ ^9 G% ~
运行结果：
+ V# R) S1 p& h! ?0 h3 cSolver stopped prematurely.
8 B- W7 M$ `; [' g7 M1 ~fmincon stopped because it exceeded the function evaluation limit,
& Z% v' s2 b1 y+ ]$ u0 J0 woptions.MaxFunEvals = 500 (the default value).
' Q# w, \2 X. ^- C2 e8 [2 u- F/ X
" f% |/ v" U* q, I$ Y5 `x =
   -0.0175
    0.4875
) f& l& g' w6 S3 l9 Y1 v   -0.0175
8 K2 }8 f! `3 O    0.3074
: d5 k: i% u0 {/ w    0.2328

. K/ Q- Q  F6 D5 j# a& Jfval =
) f' y/ `. ?; m# J' [9 l   -0.0390
: F9 {" ?! a% i
exitflag =
     0

, @5 T# P: H+ T# n; ^2 Zoutput =
+ c! @, v1 R- O% Z+ w6 o         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
$ N4 G! F0 f/ M- I! H8 ^           stepsize: 0.0812
. O* v$ M9 [. }% k* r9 m          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
& _. @' [, J1 V7 e! J5 {    constrviolation: 0.0509
            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12


8 j" |! `2 I7 w9 C3 z1 n% I$ d    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!