求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
% G7 T" l' }+ `/ \1 c; nR=0.15;
- d# R; A  u& m' O) l! X" X& r% Ru=[0.1221,0.1786,0.1452,0.1054,0.1132]';
$ \8 k! N3 @8 z9 P& W5 ~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];
3 y; P5 e- C6 f) K: `E=E/10000;
c=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
aeq=[1+c,1+c,1+c,1+c,1+c];
4 z- W+ q' `# W! h7 y0 g$ U1 Cbeq=1;
& W0 P" q" ^8 u, Y' Z5 p( `& Y/ }' nlb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
- `( Z3 _! J' c( ?[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)

function f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
6 W0 k8 w; M& O0 l' l: |/ [! d/ t" e%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
end
- ]( `& D4 `& E; h9 r7 @
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
" M4 A5 Y/ V+ w0 s4 ~3 w%UNTITLED4 Summary of this function goes here
- `2 ~  i7 r  b9 `' z%   Detailed explanation goes here
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
end

# U* d8 ]$ Q' J
5 C  ]1 \" Z: F1 N4 o+ p/ a) O" ^( ?运行结果：
% d& m5 e: k9 J! O( CSolver stopped prematurely.
! r- L% m4 R) e4 s+ B- l/ i! |fmincon stopped because it exceeded the function evaluation limit,
9 d( u8 V9 ?; A5 h; L' n7 M# l. yoptions.MaxFunEvals = 500 (the default value).
. _2 u& l7 F+ U" J1 |- O2 `
* `' @2 u/ f( d( ?( Sx =
   -0.0175
    0.4875
   -0.0175
    0.3074
) U  Z0 b9 n8 b0 i: T8 M( l    0.2328
: g- u+ l, n- N1 i% ]6 ?8 H
fval =
   -0.0390

exitflag =
     0

output =
1 |6 N% O4 z& t9 i! Q         iterations: 51
          funcCount: 505
& J3 T& Y/ U$ q! [- s       lssteplength: 0.1250
           stepsize: 0.0812
  |+ i0 ]6 |4 S6 _9 Q5 a7 d) X          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
! x0 x! ~4 G" U8 d1 y- @      firstorderopt: 2.4732
    constrviolation: 0.0509
) F6 d3 x( f  V/ G4 t7 M2 N            message: [1x145 char]
! k% c$ Z2 z) `4 |) w3 n各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

/ b  \3 P) H1 K. c
' t/ G( G3 V' r; q$ l( D% l8 l    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!