求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
% ?2 X+ ]" J# ]4 uk=0.5244;
; {+ D2 d1 A, b$ K( }8 bR=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;
6 R7 D6 ?8 T& u- Xc=0.0075;
" _/ O- t4 r5 M& b5 _options = optimset('Algorithm','active-set'); % run active-set algorithm
# @, w: E, g1 ?; Z  b4 s/ w6 `4 Naeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
. c- l2 U/ I, R[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)

* M* R% G! O8 L/ }. K" Rfunction f = myfun( x,t,E,u,c )    %目标函数
7 a5 S: Y3 U: {2 b/ Q& M%UNTITLED3 Summary of this function goes here
. A' M) x1 w& ~7 u; m0 [  h2 x%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
end

1 o' i% |( l- Y. Rfunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
7 Q# l) R5 i2 k% P: E6 U3 q! J3 Y* b%UNTITLED4 Summary of this function goes here
* ^, L7 |) R8 d8 v+ H8 ~  k%   Detailed explanation goes here
6 }# b3 H2 p4 t4 x# FC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
: d/ g) T+ U/ g- d/ c9 _+ ^Ceq=[];
" Q8 B* ^, }2 m8 A  wend

$ `1 Q- P3 L9 a! g
运行结果：
Solver stopped prematurely.
$ I# }% }( @) k, P/ Z8 J. ufmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).

, ~, H! T  Z$ `8 Wx =
   -0.0175
    0.4875
   -0.0175
    0.3074
$ D; m, i# d6 r' m* p    0.2328
6 c* [1 P/ f  @
fval =
* o6 [/ q( Q& n, R0 r0 @: d   -0.0390
- m) U* k4 J+ X( e, r2 \: f
- q# A/ z4 K! o7 g( U4 L. e. ?exitflag =
+ U$ c) W' M/ b7 y" s. G# ^     0

! F6 ?5 |/ d6 k: U9 ], g7 t4 R5 V: n" {output =
         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
4 `5 R/ i+ @3 E2 K1 R# K      firstorderopt: 2.4732
    constrviolation: 0.0509
            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12


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

火箭发动机

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

韩国福娃

不错不错，值得学习啊!