求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
k=0.5244;
R=0.15;
+ M7 E% V; c8 @% t: w- lu=[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];
% R4 ^- R- j4 e& \; N) Y( xE=E/10000;
c=0.0075;
; T9 K: V; b/ K' `0 g; Coptions = optimset('Algorithm','active-set'); % run active-set algorithm
" E8 d2 f! f# J8 |aeq=[1+c,1+c,1+c,1+c,1+c];
beq=1;
lb=[0,0,0,0,0]';
; M) {; a' F" K; ]6 _4 x0 n6 x4 S- Kub=[1,1,1,1,1]';
/ j+ \: Z+ H, _; ^9 Y3 e( Kx0=[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)
; m$ u+ N1 M$ R" _& H$ D7 E
function 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;
" G- k2 e2 d3 d% D  l+ q: Q1 ?end

function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
& o- |$ r2 Z$ p/ V5 V- D%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
% T5 J" d% V' _( w$ ?/ t/ IC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
* ?/ \# j* h# B+ ~6 j5 H: P& D2 mCeq=[];
) u. _: v4 B3 [5 L& o% hend
$ e3 f) V8 L* b. P- ]9 l3 J
4 V1 c8 B- @* k0 {% {) U
0 F) r* Y. ^' m运行结果：
Solver stopped prematurely.
  b+ k5 V4 e" i3 z) ofmincon stopped because it exceeded the function evaluation limit,
; }& |& M4 [9 ?7 c3 Z( Z8 boptions.MaxFunEvals = 500 (the default value).
( P9 {: [+ t. ]/ s: u, Q
x =
   -0.0175
; z% |) y4 T0 z. b, a7 I9 z! J    0.4875
- E* {5 {/ r6 N   -0.0175
    0.3074
' Y/ J% u* V- ]    0.2328

fval =
   -0.0390
1 _) K" U" n6 U% c
exitflag =
; E, D9 C) n9 a+ ~( z  O     0

8 s- P0 H, S! K0 houtput =
         iterations: 51
2 o/ u) M6 O. a          funcCount: 505
( Z) k9 J: d- S+ r  D9 C: U) k5 h: \       lssteplength: 0.1250
           stepsize: 0.0812
8 [' r" o% [. L1 }$ O% v7 W7 g          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
1 v; _  @- K7 T9 \# ^( s      firstorderopt: 2.4732
    constrviolation: 0.0509
            message: [1x145 char]
各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

1 P% U6 b; @( O
" U, z2 N3 |3 t/ ?' B8 k    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!