求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
2 f$ u/ [( T. p6 ?: _  C' @k=0.5244;
" K* H* E. A) E" a; D* _R=0.15;
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
5 E1 p5 F8 g" I& }, KE=[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];
E=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;
7 v8 n- T  a2 [3 P7 h: Z4 Tlb=[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)
% o6 M9 M4 U) p/ x
function f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
! Z3 w2 n3 _  F3 W* k! j3 L, w1 Y& d%   Detailed explanation goes here
) o. T  ?, ]1 Q! p3 I6 r% Rf=t*sqrt(x'*E*x)-(u'-c)*x;
end
- B, @! P: T) {: E
- C0 [. @5 c; L9 m9 [! J! }function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
# r, K: F" x, B9 }* {1 L%   Detailed explanation goes here
) r$ _: e  o1 P1 U6 _C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
end

/ Q6 c% T) l$ K6 H4 U; M
运行结果：
+ W/ I9 O& @/ F6 V6 f' P7 I! OSolver stopped prematurely.
: ]3 ~! ~% {% K5 {/ p( Z0 I" k4 Cfmincon stopped because it exceeded the function evaluation limit,
1 Z# e! e7 E: U8 y4 T0 Voptions.MaxFunEvals = 500 (the default value).

x =
  N) j2 V/ R7 _" @   -0.0175
    0.4875
( Y8 w3 u) U" Q   -0.0175
    0.3074
0 u1 h: f6 ^( H3 M& b/ F6 y    0.2328

& k& J- g1 t- \& T: {/ Ofval =
6 P+ q1 k4 }$ y9 T# _8 d+ Q0 ], A   -0.0390
& z; ]) v, A" \
exitflag =
     0
2 `/ \0 w# H; o
9 v4 a& z: I7 e7 Q" I& x. doutput =
         iterations: 51
          funcCount: 505
; o7 ^5 U8 L$ z       lssteplength: 0.1250
           stepsize: 0.0812
( j2 v6 c' J7 Y! a5 V          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
$ n- l  [9 x$ Q" e7 b      firstorderopt: 2.4732
    constrviolation: 0.0509
! v6 r$ M) r  a  K& S0 K            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12
5 @9 I1 n" z- F; T( |, S& ^
6 f* d6 v4 A, T0 Y( d
    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!