求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
5 i' S  l9 b4 V  Q6 B! h/ U7 Nk=0.5244;
R=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];
8 T* V2 |) C  r( y! GE=E/10000;
3 O) ^9 v8 k; |- o( j& ac=0.0075;
' A7 ]1 I- F- l4 |options = optimset('Algorithm','active-set'); % run active-set algorithm
* g3 B- t8 k4 f8 k* Qaeq=[1+c,1+c,1+c,1+c,1+c];
  g/ o( ^/ B% z* [beq=1;
lb=[0,0,0,0,0]';
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)
  x& H' i9 p9 l+ q3 f6 X
1 V, r% W. r1 h) a" z3 wfunction f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
, D% f* Y9 |1 E1 ef=t*sqrt(x'*E*x)-(u'-c)*x;
end
* E& l/ h8 Q' \; }6 s: H
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
- ~! w5 a3 _# p, I9 `%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);
Ceq=[];
2 }( V1 [7 O- o0 H$ ]. |end

0 V3 M, e  m( G8 e  s4 r9 W' l
运行结果：
Solver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).

3 ]! I5 o" k- s. z; dx =
   -0.0175
" Y  f. u# d; a4 k8 Z    0.4875
   -0.0175
0 M: P% n" [; F' A0 P4 Z    0.3074
    0.2328
5 u( Y, h* W! ~" U
fval =
6 Z9 z3 v& \  h# C2 @   -0.0390

* ~" t: T2 s! R2 e) Gexitflag =
8 J; a1 ?9 ?; C$ V' D& H! r     0

0 j- ~+ @* a2 b) k* f) L7 qoutput =
         iterations: 51
9 g6 K' }3 |6 e* M, b          funcCount: 505
" a$ D9 m/ Z# ]4 Z% Y( S! Y       lssteplength: 0.1250
           stepsize: 0.0812
0 ^! y/ y& ~3 r4 f! z7 A* |          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
0 c8 N3 X$ z6 z3 O* L* g      firstorderopt: 2.4732
    constrviolation: 0.0509
  o7 D! {8 x% O+ Y. O            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12


* U1 i7 @0 B# _" g8 W- Z* V5 D/ T    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!