求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
0 |! j% v, V7 B( `: JR=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;
. h5 y& c% M  g; I4 Zc=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
* q: M; F2 r0 {" N9 W3 Q( a5 ?aeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
2 Y8 U7 v  b* `, l( qlb=[0,0,0,0,0]';
( M6 k- k6 `% Oub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
1 R2 n0 _5 k' r- `2 V" f& ^[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)

% u3 n% q* ^0 }function f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
6 Z* i/ [& n5 E$ K0 ?%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
end
) [! c( N3 j4 O
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
8 S( c* s1 {: f4 N# l: |Ceq=[];
end


运行结果：
* m* S  S& W. U# M3 `1 HSolver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).

9 _: \# ^7 I& ]$ ]7 Gx =
1 _2 N# y6 p& `! z   -0.0175
9 ~# K6 ~, p* L! v7 a. n2 f/ \    0.4875
( Z, Q) W; Y* A8 m; Y: S   -0.0175
4 |9 j9 q* T0 d/ e( R$ t" u    0.3074
( m3 H: e  `  O5 g& ~    0.2328
# Q' |! {& t1 ]4 A6 z
0 v; }& l7 \* C$ ]. Pfval =
   -0.0390

exitflag =
/ Y( \1 q9 }/ t1 }. Z8 N# C     0

output =
         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
/ ?9 d. q9 _$ `0 ~- b0 [: {           stepsize: 0.0812
) B7 W# ~8 X3 f4 B  u. p/ O/ `6 K          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
1 R) p$ |6 i" Y7 G3 F+ v* G      firstorderopt: 2.4732
    constrviolation: 0.0509
            message: [1x145 char]
, `. a# n, M9 @4 K各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

韩国福娃

不错不错，值得学习啊!

火箭发动机

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

亚瑟王的冬瓜

回复 3# zfq12
% |7 a$ k# m' r9 w. Q' y

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

zhangal

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

zfq12

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

zfq12

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