求助: fmincon的参数问题

亚瑟王的冬瓜

t=2.3263;
3 @$ p6 W, n( x$ B' bk=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];
# Q2 T7 ]/ G; [) s( [) x; \E=E/10000;
7 o- x: d4 x) o1 [c=0.0075;
, a$ D* k# f2 l. @$ h7 G. T" Qoptions = optimset('Algorithm','active-set'); % run active-set algorithm
" Z  Z" n' b  paeq=[1+c,1+c,1+c,1+c,1+c];
% I0 n8 c' S: y# Zbeq=1;
lb=[0,0,0,0,0]';
ub=[1,1,1,1,1]';
x0=[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)

( q) R& B7 ^( E; b$ ^( lfunction f = myfun( x,t,E,u,c )    %目标函数
%UNTITLED3 Summary of this function goes here
/ v' _- S% J: \% F0 z%   Detailed explanation goes here
2 }" x  k4 `9 Z0 zf=t*sqrt(x'*E*x)-(u'-c)*x;
/ E+ f/ N; L8 lend

function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
%UNTITLED4 Summary of this function goes here
%   Detailed explanation goes here
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
Ceq=[];
! X% p  C: k. P! G4 {% vend
6 a5 N; F. k! E" T' D: G

运行结果：
Solver stopped prematurely.
" P' s) s6 ?2 Z% M4 ~/ kfmincon stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 500 (the default value).

x =
% P. X" R: H( H% y, n4 ?+ ~   -0.0175
, [+ q) I; V! p% C4 V- F/ g6 c    0.4875
' ~7 B5 {8 ~6 I. h9 d, t   -0.0175
' t5 d9 B  n& F% ?0 w0 }    0.3074
0 b: ?6 X) i6 C    0.2328
) g' N& C. N- \$ c" W
' p, [: ?( S& j- gfval =
: S+ A: A! J1 r( r( Y# E   -0.0390

exitflag =
     0
) G8 A4 e5 {. A% U7 R
output =
         iterations: 51
          funcCount: 505
       lssteplength: 0.1250
, b9 K$ r( _; G9 t9 I% [# D           stepsize: 0.0812
. m8 F0 Q; r! l. l8 Q2 e% H: ^          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
; U  t- U, `  g. l/ x  X      firstorderopt: 2.4732
9 o" x8 e" f9 M) g    constrviolation: 0.0509
, m, |! N) c- A- d7 ]% I, W7 e            message: [1x145 char]
7 _( M8 n3 a" Y  S8 ]& `' {* |各位大虾，程序如上所述，为什么函数提前就结束了，是不是要改默认参数设置，小弟也试过，但是捣鼓不出来，希望各位高人不吝赐教，小弟先拜谢了  Orz

zfq12

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

zfq12

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

zhangal

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

亚瑟王的冬瓜

回复 3# zfq12

; |4 D, x5 s1 X8 ?: u" T* f7 B
& z$ i9 r3 b: b( u) _" K" [2 r    我也试过把lb,ub上下界去掉，结果程序跑了一个上午都还在busy。。。不知何故啊

火箭发动机

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

韩国福娃

不错不错，值得学习啊!