求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
6 h' D7 T8 R" nk=0.5244;
R=0.15;
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
+ O( j) L* p# e" ]7 n5 T3 ]- TE=[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];
4 t" F& @( `9 R+ B7 KE=E/10000;
c=0.0075;
: O& k7 N' m5 D3 v6 C% E: E. k: @options = optimset('Algorithm','active-set'); % run active-set algorithm
9 W3 H# D5 R7 H+ d8 u: ]) ]aeq=[1+c,1+c,1+c,1+c,1+c];
& @( x% e3 Y, ?  N* J; g) d: ybeq=1;
9 a$ U# o4 z7 ~' f  Ilb=[0,0,0,0,0]';
( I* Q  |" L! v) Bub=[1,1,1,1,1]';
; m/ k. v9 f; G4 S/ m: Mx0=[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)
5 d2 }# v. ^; a' f9 |! H' y# r
function f = myfun( x,t,E,u,c )    %目标函数
, R/ ^/ c( h# s2 O6 l6 G%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
end
" q' n7 L0 t. e$ U2 Z/ ]8 R
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
! j, s2 X4 @3 V/ N* G! v%UNTITLED4 Summary of this function goes here
/ e/ {# Z9 `/ v%   Detailed explanation goes here
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
  S" B, N& n: k1 N* K* V' PCeq=[];
end

. A3 p8 S& A2 p
5 l& K0 g5 }8 v* a运行结果：
: d7 F8 u& k" s9 ~) K! u4 n" ySolver stopped prematurely.
3 @6 F2 [! G; G$ J! K! Ofmincon stopped because it exceeded the function evaluation limit,
; j. c* r; W2 r3 ^: }. ?; D7 |options.MaxFunEvals = 500 (the default value).
& Y0 s. v' s& m9 V6 i
# G  r, p. @9 w4 z& E! a$ ax =
   -0.0175
* A; Q( |8 |# V, V8 i: P3 c    0.4875
! _0 o' K" @( |3 s9 ]   -0.0175
4 Q/ w7 b$ J8 I4 T9 s: [, N1 t    0.3074
2 O) M7 J1 h2 a    0.2328
( D, ^0 H5 n! \8 D5 q
( B# w! z( L" S* Efval =
   -0.0390
, [! i* Y% @+ {; l
% F- T6 D' {9 P+ [3 \" H# s- q5 Gexitflag =
: I6 ^1 J+ E& o% `# s- f1 ~* Z     0
! I" h( y+ J: V  a$ s
output =
         iterations: 51
$ [3 K# K, q$ w. ^  P          funcCount: 505
       lssteplength: 0.1250
           stepsize: 0.0812
7 O! Q* N! W7 b6 G          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
7 q* U; P5 @  y. r4 y    constrviolation: 0.0509
            message: [1x145 char]
2 D2 k. L7 q: I% P$ N. K7 k各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12


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

火箭发动机

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

韩国福娃

不错不错，值得学习啊!