ccs.m

来自「对平稳时间序列中的参数估计方法（包括矩估计」· M 代码 · 共 19 行

function [Y,B]=ccs(p,N,A,X)
% 得出残差数据
% Y为残差值
% B为得到残差值用到的自回归参数值
% p为自回归阶数
% N为样本数据个数
% A为自回归参数值
% X为样本数据
T=zeros(N,N);
for j=1:N
    for r=1:p
        T(j,j)=1;
        if j-r>=1
            T(j,j-r)=-A(r);
        end
    end
end
Y=T*X;
B=A;