ex020900.m

Digital Signal Processing Using MATLAB V.4随书代码

% Example 2.9

%

% x(n)=[3,11,7,0,-1,4,2]; nx = [-3:3]

% y(n)=x(n-2)+w(n)

% ryx = cross(y,x)

%

% noise sequence 1

x = [3, 11, 7, 0, -1, 4, 2]; nx=[-3:3]; % given signal x(n)

[y,ny] = sigshift(x,nx,2);              % obtain x(n-2)

w = randn(1,length(y)); nw = ny;        % generate w(n)

[y,ny] = sigadd(y,ny,w,nw);             % obtain y(n) = x(n-2) + w(n)

[x,nx] = sigfold(x,nx);                 % obtain x(-n)

[rxy,nrxy] = conv_m(y,ny,x,nx);         % cross-corrlation

subplot(1,1,1)

subplot(2,1,1);stem(nrxy,rxy)

axis([-4,8,-50,250]);xlabel('lag variable l')

ylabel('rxy');title('Crosscorrelation: noise sequence 1')

gtext('Maximum')

%

% noise sequence 2

x = [3, 11, 7, 0, -1, 4, 2]; nx=[-3:3]; % given signal x(n)

[y,ny] = sigshift(x,nx,2);              % obtain x(n-2)

w = randn(1,length(y)); nw = ny;        % generate w(n)

[y,ny] = sigadd(y,ny,w,nw);             % obtain y(n) = x(n-2) + w(n)

[x,nx] = sigfold(x,nx);                 % obtain x(-n)

[rxy,nrxy] = conv_m(y,ny,x,nx);         % cross-corrlation

subplot(2,1,2);stem(nrxy,rxy)

gtext('Maximum')

axis([-4,8,-50,250]);xlabel('lag variable l')

ylabel('rxy');title('Crosscorrelation: noise sequence 2')