求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
- D; d) M! z# I# w. d  |6 ~/ pk=0.5244;
- a) }3 y* G" W* _( sR=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;
7 O/ b7 `$ N% c$ i1 mc=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
) B. |, S8 ~% @8 S% n1 Q, s0 t) Laeq=[1+c,1+c,1+c,1+c,1+c];
: n) `" ?6 j7 dbeq=1;
lb=[0,0,0,0,0]';
  T" Y6 c; ^+ }9 l8 Z& Pub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
& e* e2 m2 m% ?# Q6 |8 k: m  n[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 )    %目标函数
& ^7 Q( e5 S2 L2 Q+ Z, F, c  D. s%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
: k5 ?9 Y+ T: ^, x; nf=t*sqrt(x'*E*x)-(u'-c)*x;
0 ^6 r3 l8 G4 I) ?& kend

9 W1 w0 F" q# W4 @# ^9 j; ffunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
5 B8 Q: p: w& E%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
* G9 W) v% S8 A- ~% a7 JCeq=[];
" n) p9 m+ W: ^: S, Z# vend
' p8 b7 ]  b# v( ^+ N1 T" ^2 G0 t

8 `* e& d* R) Z8 c' O# ]) B1 F4 o- \运行结果：
Solver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).
# q& o! ^0 @0 L
x =
1 b9 R) w% s- _6 J3 d   -0.0175
    0.4875
   -0.0175
    0.3074
    0.2328
& d9 l# n. i  F2 H/ V) Y3 A" P
fval =
   -0.0390

exitflag =
     0
5 L6 Z  ]3 L0 A' r& v  \6 m9 n6 I
! g& w- I( S7 B3 Y. Routput =
+ C' l% r* j9 k7 \% f& c% r0 y         iterations: 51
          funcCount: 505
1 W: [5 `. C8 _' G       lssteplength: 0.1250
5 h8 k7 @( h8 s5 y  E$ u           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
    constrviolation: 0.0509
9 Y, E5 q* Y& A& o# U/ i            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

2 O7 j1 V" g& E3 ]9 ^$ c
    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!