求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
7 D( E- |6 J/ N5 Rk=0.5244;
R=0.15;
& t' M0 Z3 ]1 {8 U8 m# ?u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
! `; c* E' X$ {4 }2 }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];
; S6 u* ~2 `* B1 N" T6 M9 B4 gE=E/10000;
c=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
9 k' l+ E- \+ [4 eaeq=[1+c,1+c,1+c,1+c,1+c];
, a5 t/ R# w) \beq=1;
lb=[0,0,0,0,0]';
- K& x5 H8 |$ E5 ^2 w8 g" J/ Fub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
1 i3 c/ w% l9 N6 W0 H[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 z$ G# O/ O$ y" o- p%UNTITLED3 Summary of this function goes here
; B8 j) M/ g( K: v, `, w$ y%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
end

function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
" m# i4 h8 X9 z( WC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
5 D7 q- s9 r: HCeq=[];
end


运行结果：
Solver stopped prematurely.
/ ?2 N% {3 J) M" K5 Z! Jfmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).
- q, c% d7 @$ H+ o
x =
' _: e8 J6 Z3 I& {   -0.0175
    0.4875
% ?7 T: s" ~" w" H2 t' F. U7 S* b   -0.0175
$ z6 ?# @. ^9 z# v6 Q, P" P    0.3074
    0.2328

fval =
: Q, y5 B% i  r3 }+ D2 {: q   -0.0390
, R8 V  W, Z) u3 o2 B! a
+ ?9 H, K5 }3 N" B( u+ @( i/ H* Yexitflag =
% z1 S( R0 @, x     0

output =
         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
' ~9 H: w" E$ W4 c) [  T6 `9 V    constrviolation: 0.0509
            message: [1x145 char]
7 Z9 g0 l* }/ j4 K' O: n. B2 e& \各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12


* ~. \% S! c0 q7 `# h6 s    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!