求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
* x4 ]* S5 g$ j$ v7 rk=0.5244;
9 e1 A7 w2 g% N! [8 e3 ~5 e9 ?* e- q: x  mR=0.15;
* `( c! Y% I/ ^2 H9 {2 t3 Su=[0.1221,0.1786,0.1452,0.1054,0.1132]';
' r$ F2 m! f' ]- M6 VE=[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];
6 u. g% w; s/ S9 P' v4 oE=E/10000;
; J* Y/ m  H( \" x0 F, sc=0.0075;
6 q6 g. h7 d) @" E4 c. ]+ ~9 ]2 S9 Soptions = optimset('Algorithm','active-set'); % run active-set algorithm
5 O3 o* X7 w0 B3 B2 V5 maeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
3 ^  w# y5 [4 Elb=[0,0,0,0,0]';
7 P/ f% O' F: n4 {% t; x  q: k* qub=[1,1,1,1,1]';
: E3 r8 @4 D  Y. ?* zx0=[0.02,0.03,0.05,0.11,0.12]';
/ D' L* g- O5 q/ [[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)
# Q; M8 y; C. Z
# B& {: _' h  U; M- T) N+ U  H" a9 Bfunction f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
%   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
* T( @$ x4 K8 @7 ]4 G$ LC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
end
3 U8 {4 t& g% ]  A% W
  E9 @/ k1 [1 h0 ^( f+ }4 d, V
运行结果：
! [- p2 u& w* N1 e; [  S# FSolver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).

x =
* V  x2 j0 q4 Y/ r# o$ E1 M   -0.0175
    0.4875
3 j7 r; M! k: u; O/ R# ], E   -0.0175
/ S1 Z; d9 Q3 V  m$ H    0.3074
3 f& S% U* F* {9 S    0.2328
+ h2 U, e" S8 n7 b8 S1 z
' N' k  L; W; b& M' `8 ~fval =
   -0.0390
, R% h4 d$ X9 p& z6 Y5 W! Y8 p
: F4 f& d) Q: ]& Gexitflag =
. C" x: r6 x/ ?- m8 L1 ]% q     0

0 p" m- |. Q# Z1 E9 @output =
" b; A9 ]" b, B* x; i0 @" ?         iterations: 51
          funcCount: 505
4 }; x- ]9 x0 |$ }( N       lssteplength: 0.1250
           stepsize: 0.0812
# E9 m. p& L4 _. y- G          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
7 h% J8 Z6 m( a$ G! E      firstorderopt: 2.4732
- ~) S- t6 O: {7 T/ h5 z    constrviolation: 0.0509
            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

; h5 b. [' z9 D
    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!