hc236.m

来自「数字信号处理教程：MATLAB释义与实现 书中实例码源」· M 代码 · 共 16 行

%《数字信号处理教程——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')