求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
8 e% m5 w/ Z# x+ dR=0.15;
' e' k' B  g$ q5 Z3 su=[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];
# B, u1 B( y; ~) f5 WE=E/10000;
% J6 ]( D; M% ?5 Fc=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
  Q& \( ^3 Q9 K3 _! ]; A* xaeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
  C2 S  Y* y& C2 k7 W3 e5 o[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)
# r+ X7 ?5 [$ K& Z- F% z3 v* J
function f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
6 {1 p# {6 [# Tf=t*sqrt(x'*E*x)-(u'-c)*x;
5 h% h7 Q  U- V% Send
) f5 ~2 y1 s9 L: J! l/ M' F
. O" ?1 p5 K. v; |4 R- E3 i7 ufunction [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
) Z! ]1 M( O% S8 D' @%   Detailed explanation goes here
4 w/ w# p5 }& P! C; sC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
/ ]- X. R, u% D3 ]) h+ C" aCeq=[];
; r+ T* n1 y  _0 X2 Xend
$ b; [  [7 X1 s- L: r4 x# R

运行结果：
Solver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
2 Y. @6 `) e- j! W9 ]  q8 e& ~options.MaxFunEvals = 500 (the default value).

x =
6 Q' E' g0 n6 |0 A   -0.0175
    0.4875
   -0.0175
9 i: n: \9 ^4 z! R/ h  `    0.3074
    0.2328

, Z2 j' s9 `, R: rfval =
   -0.0390
  k/ ~+ V7 R0 d. {. P
- k3 H! V3 a- b+ t& g2 K3 a, sexitflag =
! p6 N: j/ S) A) o" [     0
) d8 {' U1 m  R6 G: I
# I! S5 }* V: ]3 g& u( P% i. Xoutput =
0 X7 h/ N  V1 Q% s" R% s         iterations: 51
4 a! n+ m4 E2 I/ F6 \          funcCount: 505
       lssteplength: 0.1250
           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
3 m1 X/ `+ m! N! O4 {      firstorderopt: 2.4732
    constrviolation: 0.0509
            message: [1x145 char]
3 d( j9 E5 m9 J各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

韩国福娃

不错不错，值得学习啊!

火箭发动机

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

亚瑟王的冬瓜

回复 3# zfq12


/ E: h$ J* G) d; G% Z    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

zhangal

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

zfq12

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

zfq12

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