求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
( ]7 r& O; B) \9 ]" P3 [) k8 IR=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];
' O0 V( o' G0 {3 o; p% r( ^5 e8 D+ zE=E/10000;
* M- g" o" l4 q2 S' l9 P' d- }c=0.0075;
options = optimset('Algorithm','active-set'); % run active-set algorithm
6 ?8 Q% J6 k! }8 Waeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
9 {$ p# v9 M/ jlb=[0,0,0,0,0]';
* [. ^) o: h2 I* \2 Dub=[1,1,1,1,1]';
  I7 q* o2 T& z, w- I$ K- V: L3 \. L6 bx0=[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)
/ v. q  l) ]9 i
- ~2 r/ z8 R0 ~" Z( n0 nfunction f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
8 \' q/ J5 v8 M# w: h%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
$ f5 G9 x9 r8 ]end

function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
  b) z0 L1 U+ Q+ [%UNTITLED4 Summary of this function goes here
6 D) U7 a# L. z$ a5 K3 E) 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=[];
( ^8 J4 q6 j" g* h# ~end

7 n$ V! K+ R$ h/ x3 t+ B
运行结果：
Solver stopped prematurely.
fmincon stopped because it exceeded the function evaluation limit,
  v1 h+ ^; g1 i1 M( Voptions.MaxFunEvals = 500 (the default value).
' B1 }; B0 s' J: s) B# i
/ `* y* t# ]* ~; N- Hx =
' ]" H9 p% T9 Q  n5 S3 e" s   -0.0175
1 B9 `. i4 u3 g    0.4875
- b& P% c' R! U+ v( v2 M( ?$ M5 T   -0.0175
+ e2 ~: {/ G+ }* m    0.3074
! g$ l  _$ `  M$ l! C2 n    0.2328
9 I& s  |+ |7 j
fval =
4 [! N2 R9 F, H( o( F& v& A   -0.0390

exitflag =
' v7 Q4 K& }* e  x& D     0

output =
         iterations: 51
          funcCount: 505
7 Q. r, m1 R+ E       lssteplength: 0.1250
           stepsize: 0.0812
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
3 _- g  C% }: V% i; ?    constrviolation: 0.0509
# A- `8 d( }! e3 {            message: [1x145 char]
  J4 O6 k/ v( A$ O各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

韩国福娃

不错不错，值得学习啊!

火箭发动机

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

亚瑟王的冬瓜

回复 3# zfq12
0 v6 G, `! }( c- e3 f7 x' L

0 J& V2 P$ B0 {- S6 Z" C    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

zhangal

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

zfq12

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

zfq12

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