求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
R=0.15;
- o, P+ B5 e. G8 A! X! ?u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
$ T5 K3 N, v- ~; A, AE=[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;
options = optimset('Algorithm','active-set'); % run active-set algorithm
* w' G3 ~0 Z# d/ c5 x2 Paeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
& s" g$ N* s! w+ d& v/ ^6 jub=[1,1,1,1,1]';
; r4 V8 Z1 a. ~. J/ L- @- F  f; ux0=[0.02,0.03,0.05,0.11,0.12]';
  P7 P+ |- V6 G6 k, t% V4 Q* l[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)

* e0 a' B! a- y5 e; Bfunction f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
f=t*sqrt(x'*E*x)-(u'-c)*x;
end
: H$ h" \! k7 H- j" O% t
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
& K8 Z; n* U" @7 t%UNTITLED4 Summary of this function goes here
4 h3 P$ ~" w) P- I7 \# x. q%   Detailed explanation goes here
* ]1 x7 `" x" b: f- T( G4 n$ g! QC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
  |* N* Y0 |( [* @% S' B+ zend

: g# S! H9 m7 v* m3 Y; g
运行结果：
( d3 d5 y3 F6 p1 b* b0 y; ^Solver stopped prematurely.
; u( V6 [$ ^" ]- q; kfmincon stopped because it exceeded the function evaluation limit,
0 v/ f% `0 Z$ E0 z3 j# I) Hoptions.MaxFunEvals = 500 (the default value).
* E0 a' \# u( Y8 y3 A
% z9 Q2 P- }5 s' o& ^x =
# e* T" }5 B8 R* f, e3 y   -0.0175
    0.4875
4 v" ~( m# y# o8 n$ `   -0.0175
    0.3074
    0.2328
/ R4 T1 s2 o0 f* \' F
fval =
6 L, G9 t2 I0 g   -0.0390
  Z$ A* O1 s" k6 T$ R
exitflag =
6 j/ \8 w7 ?, h     0

! N5 C3 Y6 f3 I4 U# Y6 {/ k5 moutput =
         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
) M& f# ]( W% E0 G           stepsize: 0.0812
7 D- g! M, B4 g8 ~( K+ j          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
      firstorderopt: 2.4732
    constrviolation: 0.0509
            message: [1x145 char]
# |( y9 A  c  `4 L2 l2 m各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

5 L, e+ p/ P' x) i) X7 p" C9 g% G( Z
6 o) a& W+ ?6 _# G8 d( P4 r    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!