求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
3 y- l6 q; ?0 e9 s: Qk=0.5244;
R=0.15;
1 X0 P* D7 s9 B# N9 \u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
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];
E=E/10000;
c=0.0075;
* H; U7 @* y) {( V! S9 v& moptions = optimset('Algorithm','active-set'); % run active-set algorithm
& z2 i4 N5 ~( {aeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
3 r6 w5 D/ ]9 c$ h8 ^. f/ e! dub=[1,1,1,1,1]';
/ {. C* n% R6 d& ?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)

7 \6 ~/ M; [5 u4 Cfunction f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
$ l% i1 |6 x6 M/ u' i' b: n%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
9 w0 i5 G1 y8 Iend

function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
6 N9 ^% z- A% S" A8 H%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);
Ceq=[];
- [! W) _6 M. X( ?+ _/ lend


运行结果：
* O7 |$ R; `2 {+ lSolver stopped prematurely.
5 i8 }# I& ~4 ]6 d: f" j, jfmincon stopped because it exceeded the function evaluation limit,
$ T+ j8 e: R9 q  f; x: toptions.MaxFunEvals = 500 (the default value).
$ o# C4 \# M4 r: o5 g. y1 Y
x =
6 `4 _1 }2 o9 A   -0.0175
    0.4875
9 z9 ?( r% A! }8 s! m   -0.0175
    0.3074
% T3 p" ^' @/ I% n: c& ~: A4 C    0.2328

fval =
. {/ J2 O: @$ E   -0.0390

exitflag =
     0
; O5 m& l: K9 E" @7 R8 Q5 c
output =
) s% _8 V# u- u5 O         iterations: 51
& R1 g- ^' o# y3 j          funcCount: 505
2 w3 A  n9 k) S$ x5 C       lssteplength: 0.1250
           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
$ Y) D/ z  B. \; s  L    constrviolation: 0.0509
3 X' X, A/ C' [( e9 }            message: [1x145 char]
7 C' p* N$ Z1 B9 O1 J7 t. Z各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

+ @" ]3 f+ t7 B3 |9 P, ~
    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!