求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
R=0.15;
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
2 i  o8 {7 p# y% Y+ SE=[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];
. s* k; J6 L' n4 o0 NE=E/10000;
c=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
' M2 j9 o" ]0 p6 b5 e( T5 Qaeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
' R" r  ^& v+ ~! E9 Rlb=[0,0,0,0,0]';
; B  L5 C- G. }- K! Dub=[1,1,1,1,1]';
  d, K3 f. u( lx0=[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)
4 l# c5 H/ N2 l. g7 M
function f = myfun( x,t,E,u,c )    %目标函数
$ c3 i; o5 W2 s- G5 n" w, x% @%UNTITLED3 Summary of this function goes here
. I) M! z+ J8 C, z! `# h" O( Z1 ~%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
' U& N6 g8 ^- h% v: q6 k, [end
4 A+ t8 Z) |7 c  f, `
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
2 u0 y# z, R" s! d%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);
( a& [8 P8 t( V4 QCeq=[];
end


$ \/ j  r& z; Y# G$ n. G7 ?运行结果：
Solver stopped prematurely.
# e% D% c) K9 h7 D- }fmincon stopped because it exceeded the function evaluation limit,
$ [9 j- {: Q1 s! ]; _# H) a5 Uoptions.MaxFunEvals = 500 (the default value).
! @5 R  q0 w8 m
! i5 Z" ]% C; D# @0 x5 dx =
   -0.0175
    0.4875
   -0.0175
( \/ D6 s5 L6 D    0.3074
- _' S3 {$ n: Q9 R/ k' U4 U+ u/ z    0.2328
5 B% j: r- ^+ F& g
fval =
   -0.0390
3 H* g$ u/ y# ?2 [/ D! l( |) e4 x
exitflag =
& {9 K8 p9 R& w0 e! b     0
2 ^1 i) f& C0 ]1 @; ?
0 R/ K, O$ `# L2 T1 ^0 coutput =
         iterations: 51
          funcCount: 505
0 W/ s1 \/ F4 e# B- s       lssteplength: 0.1250
2 n4 F$ k8 \9 B2 A6 W3 |           stepsize: 0.0812
6 l2 K7 @+ j* L% O0 a# l          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
% P- M9 u7 z) k% v    constrviolation: 0.0509
            message: [1x145 char]
# Z3 \" m8 \5 G& u: y各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12
1 \0 a/ B5 G5 {9 Z+ ^

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

火箭发动机

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

韩国福娃

不错不错，值得学习啊!