求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
6 q1 L, _! R  }+ E  D4 Lk=0.5244;
1 F" o8 t$ z1 j4 u: \& j% eR=0.15;
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
( A+ D# p1 `  J* IE=[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;
c=0.0075;
2 k3 Y$ @3 A6 Q/ Z* Z7 g0 X" yoptions = optimset('Algorithm','active-set'); % run active-set algorithm
, h! u) W0 F$ }. C8 Paeq=[1+c,1+c,1+c,1+c,1+c];
8 B  }) d% k3 U% ]- ]2 b+ o" ubeq=1;
lb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
# C5 F- @% K# G. B* ^[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)
# p2 a: s- t- u$ w
function f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
  }' I# k+ b2 i4 ^f=t*sqrt(x'*E*x)-(u'-c)*x;
/ b9 J, X9 U! }  p6 ^end
2 M4 S% X. P( @6 r
% u5 F% m3 g5 R4 g& g+ s  kfunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
0 ^4 ]% E" F( `5 R. y! M& z' R%UNTITLED4 Summary of this function goes here
1 K# Y# o4 S! F4 E" t%   Detailed explanation goes here
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
/ F& V* b  J' z9 ~5 ACeq=[];
end
9 Z+ [/ ?; D6 ]9 m9 i8 \7 j, i
4 b; ^2 M- x5 \8 j' J  D& s
运行结果：
Solver stopped prematurely.
7 V% G2 d4 T& T! nfmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).
- \0 \& m8 p8 r( t4 A4 T" W$ [
x =
   -0.0175
    0.4875
$ c! [3 ^, X% ?   -0.0175
$ u4 _8 R, O; J& b    0.3074
    0.2328

fval =
   -0.0390

9 ?, l$ M- U1 [+ q: |' l9 ]exitflag =
% L% y9 R1 B; u. R, k' R0 J! Z     0
* N7 f5 a+ N$ W# O5 W* A
output =
         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
           stepsize: 0.0812
/ s& j5 |8 t3 T/ z9 X          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
' }6 |0 I6 Y( l    constrviolation: 0.0509
            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12


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

火箭发动机

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

韩国福娃

不错不错，值得学习啊!