求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
$ b2 g' F  u$ v. _( q; I# f8 G  b6 lk=0.5244;
R=0.15;
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
% ]' V8 g$ l$ i1 H4 t4 E* XE=[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];
: c4 M0 s; a. c& J6 uE=E/10000;
c=0.0075;
( s) T' o9 v, R2 x4 [" qoptions = optimset('Algorithm','active-set'); % run active-set algorithm
aeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
. U& L& R2 I, y% {lb=[0,0,0,0,0]';
4 D! o$ G% E# ?. Z' oub=[1,1,1,1,1]';
& F! }7 d4 G; Z5 n$ t9 e8 W+ bx0=[0.02,0.03,0.05,0.11,0.12]';
; q2 x5 m& _3 `, W[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)

function f = myfun( x,t,E,u,c )    %目标函数
3 u$ T( u" n3 G%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
5 \2 W; q9 j" O% `' r) s. U% ff=t*sqrt(x'*E*x)-(u'-c)*x;
7 e/ s4 ~5 w0 W% e$ U2 C: Aend

' ^! _3 D* ^& {; j2 }( a/ q4 Zfunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
8 w* |/ K* ~0 X; ~. {%   Detailed explanation goes here
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
end
4 T) _8 U2 v) [! [

4 G9 D1 \0 k$ t' R+ `  p- }  e! |运行结果：
Solver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
' e6 ~+ D1 r% \: |% |, Doptions.MaxFunEvals = 500 (the default value).

x =
   -0.0175
    0.4875
   -0.0175
    0.3074
    0.2328
( \/ m# u  ^# n
fval =
   -0.0390
  D0 ]$ U: {: N' t6 w
( O% t; w- W# Q( N4 ^8 oexitflag =
     0
4 w) e5 u- ?, ?6 s) o! y4 u
: }2 a9 B3 v9 o, V4 |# z# L5 Koutput =
' u! ~3 o* j0 B# a( d         iterations: 51
, Q. ]. J& H1 ~2 g          funcCount: 505
       lssteplength: 0.1250
8 ?4 s% Q0 r) V9 K5 V- q  Q8 T9 W           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
/ v9 u: V$ f& @# S% z      firstorderopt: 2.4732
/ |: \, B+ j4 `    constrviolation: 0.0509
$ Q& N( g; O$ B1 }& m) z            message: [1x145 char]
4 o4 y2 U: K; A. `% E9 _9 |各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12
, |- @* K- m5 X3 U4 X! m

( \4 L% E0 e" ]8 }& X2 b- e: u/ c    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!