example4_x.m

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

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