求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
& G: t- K# W1 d$ K  hk=0.5244;
R=0.15;
, x0 Q, G0 F, j! S8 cu=[0.1221,0.1786,0.1452,0.1054,0.1132]';
, I$ v1 [# e6 @! O3 AE=[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 w% o+ l1 QE=E/10000;
2 a  z) `' @. E/ G0 F# |5 H/ G+ L) ec=0.0075;
, q8 R4 k+ ?$ c# q$ F+ Woptions = optimset('Algorithm','active-set'); % run active-set algorithm
  ?! p- [8 n. d6 ^$ D, [aeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
* V4 U* D% B# F0 W4 nlb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
- `6 {' S3 Z3 cx0=[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)

function f = myfun( x,t,E,u,c )    %目标函数
/ N5 z$ C/ Q( P/ O( A%UNTITLED3 Summary of this function goes here
6 b& U  O/ o7 [$ H: Y) Z%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
1 x; A0 C2 u, W' k6 Wend

, u1 s5 b% G8 ]) W9 }9 nfunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
% v: g- R( @$ a8 `$ w%   Detailed explanation goes here
5 H5 x* b5 r" jC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
end
( B2 O" U& w2 X3 n( [  _9 N

运行结果：
3 o5 I' ]8 X4 S/ }2 tSolver stopped prematurely.
3 I4 f: [4 P8 N4 M, l- Jfmincon stopped because it exceeded the function evaluation limit,
' _: p! M* a3 Aoptions.MaxFunEvals = 500 (the default value).
5 Q7 U; n* d* J2 k! v" Y2 [
  ]" M% g, H' [) B, X$ |' v/ Bx =
   -0.0175
    0.4875
7 n, _/ k# J7 s# z( d' E2 \   -0.0175
    0.3074
+ I) k. ^# C2 p) C8 V, g- V    0.2328
+ J7 z  N" W) w& S; {( M
5 _0 I& w) g; Xfval =
9 U' j0 T! K& y+ k4 F0 t  C8 a   -0.0390
4 L$ p0 T! d$ y  T4 }! }0 ?
exitflag =
0 e% P7 D0 I3 v2 ?     0
+ h0 k+ o4 G8 W2 ^, P& f
output =
+ @6 C! n/ z+ ?: c         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
* i# U5 O  ]' r0 M. ^: e) k7 W           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
6 ^% x8 }, I: u- i      firstorderopt: 2.4732
, H& Y3 G4 {% }) v1 _% h' ]$ ^    constrviolation: 0.0509
6 X5 n: }2 V' W            message: [1x145 char]
' v9 n8 I' s* a7 U各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12
2 l/ [, R) ~' D8 A

& @6 S: S5 E. F+ v  o    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!