%《數字信號處理教程——MATLAB釋義與實現》第二章例2.3.6程序hc236
% 信號序列的周期延拓
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
x = [3, 5, -7, 2, -1, -3, 2]; nx=[-3:3];    % 給定信號x(n)
[y,ny] = seqshift(x,nx,2);                  % 得到遲延信號 x(n-2)
w = randn(1,length(y)); nw = ny;            % 產生 w(n)
[y,ny] = seqadd(y,ny,w,nw);                 % 得到遲延加擾信號 y(n) 
[x1,nx1] = seqfold(x,nx);                   % 得到x1(n)= x(-n)
[rxy,nrxy] = convwthn(y,ny,x1,nx1);         % 把y(n)與x1(n)卷積求互相關
stem(nrxy,rxy)
set(gcf,'color','w')                        % 置圖形背景色為白
line([-4,8],[0,0])
title('延遲加噪序列與原序列的互相關')
xlabel('ny'),ylabel('rxy')