求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
% _" ~3 T% K7 x" E* k3 u3 z8 _k=0.5244;
  @$ H( P; z; A& Z  s* j& E' G( NR=0.15;
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
4 I) Q8 I( X* a1 p) TE=[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];
" t% e0 r4 x% `$ k6 d; zE=E/10000;
; }, g) O9 s( g& nc=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
8 s. ]( ^& a  ?( m, ^/ l) jaeq=[1+c,1+c,1+c,1+c,1+c];
+ g' A5 g* }8 lbeq=1;
) {3 `% ^5 O7 \% w8 glb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
6 R& P; _( l" Q" P* r[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)
- d, O$ V/ @0 l$ q
( T7 u7 F0 `& F7 N5 vfunction f = myfun( x,t,E,u,c )    %目标函数
7 k( k  I+ |2 z" E( v8 @%UNTITLED3 Summary of this function goes here
/ r' e# y5 y( g2 b; z* J/ m8 y%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
; Q0 f/ r7 ~4 N" e; send
0 M0 O9 d# s) f! l9 B1 \6 s
& s* _$ P$ `3 @function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
( e9 ]7 r8 a: @) b$ d' q%   Detailed explanation goes here
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
, ~# V/ t7 _0 o1 f$ ]end


$ @# R2 p( d3 {+ J* B6 S" g+ p+ b运行结果：
Solver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
$ R( I7 ]9 b* T" g. G6 _; }0 G* hoptions.MaxFunEvals = 500 (the default value).

x =
% K  K5 ~1 x& Z$ u( o0 p" q8 ^- R' ~   -0.0175
    0.4875
7 Q: M7 X( }* l  M, @- o+ x. U   -0.0175
; `" s2 |$ g+ c+ Z    0.3074
    0.2328

fval =
   -0.0390
9 j; D& r# l6 d, Z+ |6 \8 Q
exitflag =
     0
. o3 ~& u7 y* Y. n( b
# I0 Z7 B$ j: k$ _% N7 a2 S/ U( ?output =
         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
1 l6 M; l0 |% ?3 \/ Y5 m! O$ z' X           stepsize: 0.0812
& `0 A' B% {6 S  \, t          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
    constrviolation: 0.0509
# ]( C1 f# }! E! h& l; X            message: [1x145 char]
" p; J% m# K$ Z. o各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

$ f$ ?; S$ C# Z
8 E+ Z  {/ T0 L+ e" S    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!