求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
1 n* H5 a( Z) R% _5 ek=0.5244;
) e1 x! @: ]; o* L% X# hR=0.15;
6 ~& |) d# b# Du=[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;
0 C1 Y( }: {6 `. x* `c=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
  H3 b) Y8 Y6 }3 Maeq=[1+c,1+c,1+c,1+c,1+c];
2 o; U. i7 O( y) X3 l! Vbeq=1;
5 B7 Y+ b. i% Klb=[0,0,0,0,0]';
' Q: I" @9 C3 d1 R9 C6 G! mub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
: c2 j  N' P/ o7 ?[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)
/ e3 f, i. ~+ I
function f = myfun( x,t,E,u,c )    %目标函数
% l' d% A7 V; I$ o( }%UNTITLED3 Summary of this function goes here
9 W4 V0 \* b$ o" q, j%   Detailed explanation goes here
" [1 d4 i  T6 k0 H: e" d; T; \f=t*sqrt(x'*E*x)-(u'-c)*x;
# `. B- T( C. R7 V7 M0 S; G+ lend

$ q! C2 \8 K4 v% _function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
9 N* B/ E& d9 {& q2 r5 e%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
  Q% T, X. ?( E. HC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
end

+ a2 R- I0 D; D. O; u& z, E
) {% F0 N6 Q: g9 z0 H+ p& ?运行结果：
3 d8 x5 {- m: L  U% F$ JSolver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
" b+ @9 N- O& y2 o4 E4 y1 q3 {  Coptions.MaxFunEvals = 500 (the default value).
) \) p7 n% A9 X- X
* w+ S1 p3 F7 Y9 R$ T1 D- Y+ \x =
# }: L' v! Z  u1 H. U; O/ O2 F   -0.0175
) }3 u: R- z9 I) i. Y6 y    0.4875
8 r2 W$ ~6 `; h1 A7 o$ r$ t) d   -0.0175
    0.3074
    0.2328

" R; @4 _9 y5 @6 w* o6 u2 ~! r; m3 Ifval =
   -0.0390

- ^% K# r+ ~# m5 p9 @- Hexitflag =
" p7 B, _  |7 l7 }. \* t7 @7 T0 P     0

output =
- _- `4 g; [: x         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
5 N# T3 J2 T4 @  [* U      firstorderopt: 2.4732
' \) o9 }' a( v3 d8 }( U    constrviolation: 0.0509
            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

/ R, |; z; A3 l% k+ b# ]- G
0 i: w( m3 S/ h1 q, K1 c    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!