求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
7 k. \9 K5 E- `' ck=0.5244;
R=0.15;
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;
7 b7 B- X  t* b7 Hc=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
# q  k) @5 U. g+ V6 _aeq=[1+c,1+c,1+c,1+c,1+c];
+ ]# Q0 C- m( Y$ {' {beq=1;
lb=[0,0,0,0,0]';
4 O$ T& U+ S2 L, mub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
; ?: u3 |" Q# h( 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)

, u' c, B) L" gfunction f = myfun( x,t,E,u,c )    %目标函数
7 S$ c( ^* Q! g%UNTITLED3 Summary of this function goes here
$ I( `1 j1 l$ T% N%   Detailed explanation goes here
9 p, K, H' T5 F# T, m# @/ ff=t*sqrt(x'*E*x)-(u'-c)*x;
end

function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
* |2 c3 x0 N/ i) R% `0 z%UNTITLED4 Summary of this function goes here
2 n" o  s$ w+ g  j& m# _, z2 z%   Detailed explanation goes here
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
7 E2 w! k2 z: t6 h2 \  oend
' i! w# o  h% o  G0 ^0 u0 U
8 T8 p& @! K; q% V  ^
运行结果：
# _6 f8 }5 J# R. _Solver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).

$ U/ M: a" Z, V, z0 u/ O+ \x =
   -0.0175
! B  ^" L* O. T! Z: L. i    0.4875
- V0 z$ ]4 f) [1 N$ p   -0.0175
    0.3074
+ _9 w1 w7 Q4 H' l  f3 p+ r  q    0.2328
) H% I2 r9 |, [4 f4 b$ Y
& V0 Q  C8 c# g: Efval =
0 h" y% k% F9 X$ U* R0 t   -0.0390

exitflag =
0 Z/ t: \- y. ]3 _1 ~# e     0
7 E$ g2 ?4 R0 x* }! a
1 n: v, u$ t, i. x6 z$ Coutput =
1 {# \9 e% N& ^" B: M+ M         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
6 y0 {# q+ s# g& y6 z  [& ^7 }' b    constrviolation: 0.0509
            message: [1x145 char]
0 ~4 l- D. w* D" g+ e7 \各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12
* i1 m1 |0 h- L

    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!