求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
0 c8 x5 _& d5 f% _5 C# ?/ e! BR=0.15;
# P) H" b' T6 x, {6 i( \u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
" D& x9 ]6 E8 u7 ^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];
" C5 d/ ]0 W/ iE=E/10000;
& f4 j1 K; S6 P) Hc=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
: x; k6 i( Q2 S* P" e! Paeq=[1+c,1+c,1+c,1+c,1+c];
+ g4 I  L( k6 J# z& W" y8 f5 \beq=1;
7 V, h! O6 k$ ]lb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
5 p( B+ k! d" r* fx0=[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)

" O" N( Z1 G' ~/ J( O9 X- {9 Ffunction f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
+ _, A8 y  a* A" y" r) H9 F%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
4 O' ~9 r% z9 B2 H2 s- G' s3 C1 ]7 ]! hend
2 r+ U9 M2 S, A7 J/ T( N
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
: X, X8 q+ g$ Z2 r# jC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
& d: X$ W: D% J6 z, WCeq=[];
; P" Y1 w' X' z+ \' Oend

1 C( `4 G# \6 l) h+ L
2 t9 O2 l" k. O; C; F运行结果：
% b, J4 S" b8 T/ sSolver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).

+ v5 [: D" I! J9 _, b9 S$ ~0 _x =
   -0.0175
; E4 S: ~8 T+ |    0.4875
& g, r( ?6 S8 z3 L   -0.0175
! K; H) V0 r' b9 Q    0.3074
    0.2328
6 {# w! G4 P8 b& x. }
fval =
. g; @! p8 O; k" |   -0.0390

/ r3 S- Q/ q! B* }. fexitflag =
     0

output =
         iterations: 51
          funcCount: 505
5 t* l) A( y6 l1 p: o. j       lssteplength: 0.1250
) K+ Y5 q/ z: a1 p$ b  n2 H           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
+ J" |  v2 a6 P5 k$ x& I) q      firstorderopt: 2.4732
2 j) ?6 q7 u  u    constrviolation: 0.0509
            message: [1x145 char]
7 l5 L, {' k$ {) O$ Y) q各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

/ a; p: {" p% v# p) g. }3 M" t
' e( \- a+ B5 h    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!