example1_8.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 12 行

%例1-8程序,example1-8
x=[3,11,7,0,-1,4,2];nx=[-3:3]; %给出x(n)
[y,ny]=sigshift(x,nx,2);       %得到x(n-2)
w=randn(1,length(y)); nw=ny;   %随机产生w(n)
[y,ny]=sigadd(y,ny,w,nw);      %得到y(n)
[x,nx]=sigfold(x,nx);          %得到x(-n)
[rxy,nr]=conv_m(y,ny,x,nx);    %卷积运算
%画图
stem(nr,rxy)
axis([-4,8,-50,250]);xlabel('时间');ylabel('幅度');title('互相关');
m1=[-4,8];m2=[0,0];
line(m1,m2);