求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
R=0.15;
4 C% {; [2 ?9 n8 N! c" @2 Ku=[0.1221,0.1786,0.1452,0.1054,0.1132]';
. @, w# }2 O7 S' Z# z9 C: oE=[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];
/ ~' h8 Z$ ]2 |! ~2 DE=E/10000;
c=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
aeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
' H5 G( j( B9 @+ E/ nlb=[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)

function f = myfun( x,t,E,u,c )    %目标函数
" O; e! o1 Q2 @4 m9 D1 d  H8 r%UNTITLED3 Summary of this function goes here
, P% h$ P; M; ^%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
1 w3 `5 g2 ]" p) D' vend
, n9 e- S5 ~( ^* L
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%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);
( B4 y* r7 m& R$ ]& s$ X7 ZCeq=[];
end
9 a1 _* c, e# d* q

# z3 J' U% r* f+ d/ y) s1 L运行结果：
* B& p0 R" e! F5 fSolver stopped prematurely.
; X) c( P5 \  i3 M, sfmincon stopped because it exceeded the function evaluation limit,
$ c* S1 N# P! o' ^% ?/ a$ W# Koptions.MaxFunEvals = 500 (the default value).

x =
   -0.0175
2 z3 `5 P% y* _4 S- [9 V  i8 h: ?    0.4875
   -0.0175
- T2 r, k; r; q9 C    0.3074
% U* W) d  u( r) T    0.2328

+ ]& H, Y" L$ e1 Tfval =
   -0.0390

exitflag =
     0

output =
         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
           stepsize: 0.0812
: S+ y! B. J/ N; ^0 x. x* E( X; _          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
. p+ T$ x; W( m1 m) U# L9 ?# t4 z    constrviolation: 0.0509
            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12


' @5 ~& B0 @1 }& |! G4 A    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!