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;