求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
R=0.15;
) ?" m* L& J: m0 N5 p! j& Y0 Vu=[0.1221,0.1786,0.1452,0.1054,0.1132]';
2 z. M" `/ T( R3 B' N( z" p* zE=[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;
& p9 Q8 }' K. V, E% R1 Jc=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
aeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
, p, ^. w, `- z' P6 {ub=[1,1,1,1,1]';
x0=[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 )    %目标函数
4 C$ A' Z0 b4 j3 k; p/ ^%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
& ?0 y4 c+ @3 \f=t*sqrt(x'*E*x)-(u'-c)*x;
1 P, ?3 X2 B$ D; [7 b2 G8 [end

" A: Z) F8 f, V# W/ g. {7 }3 R' {function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
* R- b' q9 Q5 p% [2 n1 Q: }%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
% `+ L" j! g/ B; v; O0 B8 s3 ~( U1 N5 ~C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
0 {1 j; [( }, c) u3 x: W) s! kCeq=[];
end

, P: g- J% [3 A0 D8 V
运行结果：
4 H4 F* h/ }# P0 SSolver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).
" C8 _# s' F% M0 y  r) a9 s. r
+ K2 k9 R% |; G$ j+ D# ]x =
3 q- j- M$ f5 B* d( P5 c   -0.0175
  b8 g5 ]$ T1 D% d  s    0.4875
0 R; u; P5 {/ O9 I7 l* W   -0.0175
1 k& G6 a( g" n* L' q    0.3074
0 u2 h  s4 W% `" @  v8 Z6 t    0.2328

fval =
   -0.0390

exitflag =
     0
0 l" P/ ]. ~3 P7 h/ T5 x
$ V8 z3 l0 t# f4 [; K, x& Uoutput =
         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
8 I! K$ q7 ~* b& r  i  I# k           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
! z* @  O2 \/ U: Q      firstorderopt: 2.4732
    constrviolation: 0.0509
" k) ]) r+ z3 m2 Z            message: [1x145 char]
/ }1 v) ]+ B" q) Q7 {9 o* V& o8 N各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

( F3 m- q- J7 E6 y5 ?% s$ B  \
( f/ y5 x- k# ~6 R$ Y- z7 h    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!