sacf_d.m

来自「计量工具箱」· M 代码 · 共 38 行

% PURPOSE: demo of sacf() 
%          sample autocorrelation function estimates
% 
%---------------------------------------------------
% USAGE: sacf_d
%---------------------------------------------------


nobs = 110;

% generate arma model data

e = randn(nobs,1);

y = zeros(nobs,1);

for i=3:nobs;
y(i,1) =  10 +  0.8*y(i-1,1) + e(i,1) + 0.25*e(i-1,1);
end;

yt = y(11:nobs,1);

m = 30;

acf = sacf(y,m);
pause;

pacf = spacf(y,m);
pause;

tt=1:m;
ttp = tt';

fprintf(1,'Lag     acf        pacf \n');
for i=1:m
fprintf(1,'%2d %10.4f %10.4f \n',ttp(i,1),acf(i,1),pacf(i,1));
end;