求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
. `) m) C  z7 W4 b/ G7 N) x6 y9 \% Ck=0.5244;
R=0.15;
. H# c% x% L+ F) m& ?u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
" q$ w! v& K3 a9 i, mE=[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];
; P; O* c% v+ k- ME=E/10000;
# S' Y$ E- ~6 J* a& Xc=0.0075;
# x" y# O' J6 Q( \6 Zoptions = optimset('Algorithm','active-set'); % run active-set algorithm
4 `2 A) A$ V( O9 }, D8 H& x* Xaeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
6 S. R! n, r4 M9 @' Sx0=[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)
) j$ Q1 p- d) D0 O& D! ]
function f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
# t9 [6 I$ {8 D; u- O6 i' ^5 X%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
end

function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
+ j9 O8 H8 S$ r1 B: Q8 r%   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
& z7 x" H4 S) K- t' U2 E* E
' R0 g; D# e( r/ D
& T' W1 F$ a( Z! C( `运行结果：
& B& h7 B$ b- o0 M8 b* zSolver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
& _5 Z! @5 \2 a- j/ m* T8 Zoptions.MaxFunEvals = 500 (the default value).

0 _  }& |( t' Z6 x# q, \x =
   -0.0175
% t5 j, q+ s$ h- Y" g    0.4875
   -0.0175
! v$ P% D* L! g  @3 R    0.3074
. n  S* m& _+ R% ?, \: P+ N  H* g    0.2328
+ N1 y5 [$ r- V- m1 H  i7 w; n
: e: Q3 N* O, Tfval =
   -0.0390

exitflag =
% F( |  e! o4 @     0

0 J( p+ r8 G* U6 \! Joutput =
         iterations: 51
9 L8 x* ^& Y  E. `  C4 K' B  r          funcCount: 505
       lssteplength: 0.1250
           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
2 y5 t$ O: \8 ?2 {% j, y- q    constrviolation: 0.0509
            message: [1x145 char]
" U+ z' ~: i. \+ n) J各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12
: `2 a! v  l, J  F

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

火箭发动机

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

韩国福娃

不错不错，值得学习啊!