juanji_zxiangguan.m

来自「用MATLAB实现的标准卷积的算法源程序」· M 代码 · 共 41 行

%comprasion between  signal shift 2 and shift 4 
x=[3,11,7,0,-1,4,2];nx=[-3:3]
[y,ny]=sigshift(x,nx,2);
%y=x;ny=nx+2;
w=randn(1,length(y));
%nw=ny;
nw=[-3:3]
[y,ny]=sigadd(y,ny,w,nw);
%y=y+w;
[x,nx]=sigfold(x,nx);
%nx=-nx;
%x=-x
[rxy,nrxy]=conv_m(y,ny,x,nx);
%subplot(1,1,1);
subplot(2,1,1);
stem(nrxy,rxy);
axis([-4,8,-50,250]);xlabel('lag variable 1')
ylabel('rxy');
title('crosscorrelation:noise sequence 1')



x=[3,11,7,0,-1,4,2];nx=[-3:3]
[y,ny]=sigshift(x,nx,4);
%y=x;ny=nx+2;
w=randn(1,length(y));
%nw=ny;
nw=[-3:3]
[y,ny]=sigadd(y,ny,w,nw);
%y=y+w;
[x,nx]=sigfold(x,nx);
%nx=-nx;
%x=-x
[rxy,nrxy]=conv_m(y,ny,x,nx);
%[rxy,nrxy]=conv_m(y,ny,x,nx);
%subplot(1,1,1);
subplot(2,1,2);
stem(nrxy,rxy);
axis([-4,8,-50,250]);xlabel('lag variable 1')
ylabel('rxy');
title('crosscorrelation:noise sequence 1')