求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
R=0.15;
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
  q  d9 o/ B3 u9 jE=[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;
! j1 i: _' D4 x7 C! L+ [0 t3 Ioptions = optimset('Algorithm','active-set'); % run active-set algorithm
# G# y2 h$ U0 ~+ _# y5 I+ A! v2 oaeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
' L( h+ o2 Y4 O& z0 Kub=[1,1,1,1,1]';
x0=[0.02,0.03,0.05,0.11,0.12]';
6 Y* L7 i0 V# g, t[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)

/ y1 L6 r7 M( a" O: Tfunction f = myfun( x,t,E,u,c )    %目标函数
2 a& [" U; v  S, n3 H6 S%UNTITLED3 Summary of this function goes here
4 J! q$ @# V$ T3 b$ x; K%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
2 Q4 @/ ]1 U7 k. _. W% L" ~3 pend
7 v8 U, O7 \! R: R1 X- f5 d
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
: _$ F0 K% H2 S%UNTITLED4 Summary of this function goes here
( u( V; G' |% A- P/ C7 G% h%   Detailed explanation goes here
9 \1 O7 U; V. h/ ~- |( `C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
- D7 L7 \0 D6 Q0 M) E9 qend


运行结果：
. Z( U2 I; p) ]/ r7 h! fSolver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
# l" h* p5 E+ v% R. ~' Poptions.MaxFunEvals = 500 (the default value).

x =
# v2 a6 Z+ m7 f' [% J   -0.0175
    0.4875
   -0.0175
/ R( v# i, Q  O+ h9 z; ]  S    0.3074
    0.2328
! i! G3 ~! S" L5 y! F2 Y1 m) |
( I; y5 T. C/ X% G7 gfval =
8 D$ o3 E( d) L   -0.0390

exitflag =
9 m; W% f/ @2 ~% Y, P     0
8 W! I. `9 \8 b
output =
9 A/ Y3 Z3 u6 T; J) e# F$ X0 R         iterations: 51
          funcCount: 505
0 F, q1 [7 n3 g4 t# z9 A       lssteplength: 0.1250
$ V1 h/ Q8 Z- y. C3 o           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
2 Q9 W4 o, ~. Y    constrviolation: 0.0509
! i, k" a/ K/ x# p6 R: ^$ h2 M            message: [1x145 char]
9 a! I, N' x9 ~$ }# C* O( P各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

韩国福娃

不错不错，值得学习啊!

火箭发动机

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

亚瑟王的冬瓜

回复 3# zfq12

. H0 k: ^! b4 Z+ o- c3 A( b
    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

zhangal

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

zfq12

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

zfq12

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