求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
9 F" s( j  i) ]/ @( _4 Vk=0.5244;
R=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];
E=E/10000;
c=0.0075;
2 Q- N6 D1 Y' |0 N8 boptions = optimset('Algorithm','active-set'); % run active-set algorithm
aeq=[1+c,1+c,1+c,1+c,1+c];
! k, V2 q5 d0 j. A% V; i5 _) ]beq=1;
. F; h' y& N0 W/ \0 h0 E5 ?lb=[0,0,0,0,0]';
$ q) i5 I: F& h! o8 rub=[1,1,1,1,1]';
; o2 F/ L8 H# }& M1 f& z  ~x0=[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)
- Y: W8 j: ]% i
function f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
, h4 g- t( b& W% x; q%   Detailed explanation goes here
1 K$ m" a/ P( _5 n3 K( Q+ y2 rf=t*sqrt(x'*E*x)-(u'-c)*x;
7 ?2 i% ~+ n. d9 bend
, M: |/ t8 ~4 E  ~* a7 r6 `
* z- d! b$ ^: [1 j( q' D2 Ifunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
# X# J% [1 m( b# {+ uC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
1 g. Y2 i0 Q$ }  jend


运行结果：
Solver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).

& G+ D9 h) x; j# g2 ~0 N3 wx =
   -0.0175
; R  @2 V- Y: o/ ?# P& ]1 Q    0.4875
7 w2 p  s3 V9 z& g1 l   -0.0175
8 E. ^% N  E1 N$ R: Z; {5 h7 C    0.3074
6 T3 M  l+ I. U0 l# N5 ~: N! ]    0.2328
( D" f5 U' V2 g. F  T
( y& S0 J  X& r* T# w: z* o1 Ufval =
   -0.0390
+ W% V- R- v4 _' H- A
2 C" u8 R; X$ pexitflag =
$ M# N+ _/ X+ ]     0
4 h' z2 A0 ]1 d( O, Y" Z, M; J
output =
         iterations: 51
* Q. i- g1 b$ w+ v8 B          funcCount: 505
       lssteplength: 0.1250
           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
    constrviolation: 0.0509
            message: [1x145 char]
; y2 f( u9 {  O  e各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12


    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!