Matlab的fmincon一个问题

uelophim

求教大神：
6 o4 m8 e3 a8 h% y& \2 T. ]  [function f=nsllh(x)
9 P9 _( t, Y: i  L; L) iclose=xlsread('HW_3_Part2Data','#3','D31257');
nvix=xlsread('HW_3_Part2Data','#3','H4:H1257');
alpha=zeros(1,1);beta=zeros(1,1);theta=zeros(1,1);gamma=zeros(1,1);
x=[alpha,beta,theta,gamma];
1 H% y. H( H4 A/ S1 g) hr=zeros(1254,1);
for i=1:1254
+ v) C/ R% g6 ^5 E/ K    r(i+1,1)=log(close(i+1,1)/close(i,1));
" k6 t8 w/ r) h6 jend
0 o, c5 c3 Q( q0 c- Gsigma=var(r);
omega=sigma*(1-alpha*(1+theta^2)-beta);
. k6 r, e6 r6 n. m7 W# @! \0 Mnpsigma=ones(1254,1);
$ t& H. L5 L- z8 Q$ F+ [npsigma(1,1)=sigma;
for i=1:1253
( ]2 ?# ]4 D) A2 ?npsigma(i+1,1)=omega+alpha*(r(i,1)-theta*npsigma(i,1)^0.5)^2+beta*npsigma(i,1)+gamma*nvix(i,1)^2/252;
+ W. N! u2 P- h3 q  dend
alpha=x(1);beta=x(2);theta=x(3);gamma=x(4);
for i=1:1254
/ g5 Y6 N1 C$ u% _# `5 O& n* fnllh(i,1)=0.5*log(2*3.1415926)+0.5*log(npsigma(i,1))+0.5*(r(i,1)^2/npsigma(i,1));
end
f=sum(nllh);
$ Y% h2 B5 k* m7 t4 L9 N
5 E" h9 H. i3 ?4 v另外nonlcon的code是：
. n, G, Q5 V' }function [c,ceq]=mycon(x)
alpha=zeros(1,1);beta=zeros(1,1);theta=zeros(1,1);
5 R' Y% s( @( d1 |  J  ]x=[alpha,beta,theta];
c=alpha*(1+theta^2)+beta-1;
( y0 S9 E2 X1 `0 Nceq=[];
4 }1 R/ ]4 I0 u
run出来显示：
% ^9 z/ S) Q5 q, I# Hx0=[0.04;0.5;2;0.07];
) G2 x* ~5 {, _2 \Optimization terminated: first-order optimality measure less
! L# h* v0 D" `* L6 X5 a5 O7 Bthan options.TolFun and maximum constraint violation is less
$ S) _# d, ?' p! V! p7 ?3 A% ]- V7 nthan options.TolCon.
No active inequalities.

x =
3 o$ V/ G( n! H
5 J4 v- J2 C/ s2 [0 L$ C  I  g- p    0.0400
1 C2 z' A( q) E4 L( Q    0.5000
# F- y9 \: Y3 f/ p! {    2.0000
& i! }( n% v$ E$ R. S* V5 I5 P    0.0700
( e0 W5 l, t" q+ Q* s; t2 x4 T

3 C3 d; |9 w0 @# Z2 v  j' Wfval =
) Q  t/ J) A  ]$ R# ^$ g, C
-3.6963e+003
1 }6 o& A- T( ?5 v
+ o  d1 u( W9 y6 \
4 f. H% c/ W6 n0 U: Zexitflag =
% w" K% G3 D8 E2 ~0 m9 a* t1 A
7 o8 v4 V0 m! w( I' t0 h     1
* R. Y/ ]! g/ k+ _+ x/ x
) t  D$ v: R  p* C+ E- f" X: Q8 N( {0 F: k
output =

       iterations: 1
) S( w  _- B+ y5 k7 C        funcCount: 10
- {) ]4 t2 \# ?. d         stepsize: 1
        algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
    firstorderopt: 0
     cgiterations: []
          message: [1x144 char]

7 d/ m9 a# K! b, f' X9 s4 }, p
( @, P; w) H) a" p! H( [$ [( E8 s9 A  rlamda =

         lower: [4x1 double]
         upper: [4x1 double]
' ]8 B0 P) s# E9 U5 I         eqlin: [1x0 double]
      eqnonlin: [1x0 double]
       ineqlin: [1x0 double]
    ineqnonlin: 0
7 I8 F; R! O5 }7 Y: V- B# N- _
, V) e, W( Z0 l& q. ]& Q为什么会这样呢？fval已经改变了但是x不会变？而且iteration只有1

" W! ?0 t! K% H5 x1 E

大小的小

帮你顶一下

shuddkk

1、没有你的xls数据文件，没有进行测试；
2、你在哪儿使用了fmincon，没加上啊；
5 M$ t$ X, s7 T3、调整初值，再测试；
4、设置fmincon的options进行测试；

小白的白

来学习一下

uelophim

shuddkk 发表于 2020-8-5 13:17
% ?, U  c8 z/ v5 L$ ?9 k' J- y1、没有你的xls数据文件，没有进行测试；
2 M. C) t9 _$ c5 h4 [+ R2、你在哪儿使用了fmincon，没加上啊；
3、调整初值，再测试 ...

! T# M3 C% F: b* r7 Z" A; C好的，我先去试试，谢谢