求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
( {& P4 e$ u' UR=0.15;
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
' s& V7 J% R' i" n% {: pE=[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];
, \+ C3 [: R' k. RE=E/10000;
$ x. n( L) ], T: R9 ]c=0.0075;
. N4 p$ `. D' G% m% ~, Hoptions = optimset('Algorithm','active-set'); % run active-set algorithm
aeq=[1+c,1+c,1+c,1+c,1+c];
! e) m0 D8 E* l: l1 Kbeq=1;
lb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
+ R5 V9 F: Z" I' y4 ~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)

- b0 ~8 q: R( Bfunction f = myfun( x,t,E,u,c )    %目标函数
  E% ^* p% n7 H+ t6 ^2 H%UNTITLED3 Summary of this function goes here
1 w5 D8 d" }. k%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
: G. w  _, A& ~! q# Wend
3 }4 t2 C) T0 Y! O$ G8 l, R) d
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
; K2 _; s; Q/ L9 V%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
' g6 t# p; H  f2 CC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
8 o. @) L/ ^) x& {Ceq=[];
end


$ J! y" l  H2 H- V  R  S  B运行结果：
Solver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).
( D. i! N# K5 l
x =
0 n) y8 e* G& l1 {% R- D# T   -0.0175
& G- c/ a- ^3 H9 k5 ~    0.4875
   -0.0175
    0.3074
    0.2328
2 o* o7 q! b* d5 a
fval =
7 s: S/ L: I+ i, ]! ^9 E- m   -0.0390
3 j) O& k4 P  v( h3 X$ S8 t
7 o2 s3 {% c7 v: sexitflag =
- ]( ^- X/ E1 x1 k& F3 q! w     0

output =
, O6 v! Q; [( i' K" m! Z4 Y         iterations: 51
          funcCount: 505
# [( y  R: G. F4 K8 |3 o6 z       lssteplength: 0.1250
           stepsize: 0.0812
# }$ H. }) a2 F: o9 @9 \7 B2 [- Z          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
! ]0 v1 F  }. E# V      firstorderopt: 2.4732
    constrviolation: 0.0509
3 E% V: o0 L1 T            message: [1x145 char]
6 H9 G- j6 [; J6 k, F各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12


$ @; C* Y$ Q$ h& q* _5 a    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!