《DSP using MATLAB》Problem 2.10

代码：

%% ------------------------------------------------------------------------
%%            Output Info about this m-file
fprintf('
***********************************************************
');
fprintf('        <DSP using MATLAB> Problem 2.10 

');

banner();
%% ------------------------------------------------------------------------

nx = [0:1:199]; 
 x = cos(0.2*pi*nx) + 0.5*cos(0.6*pi*nx); 

alpha = 0.1; k = 50;

[xk, nxk] = sigshift(x, nx, k);
[y, ny] = sigadd(x, nx, alpha*xk, nxk);

[x1, nx1] = sigfold(x, nx);
[y1, ny1] = sigfold(y, ny);

[Rxy, nxy] = conv_m(x1, nx1, y, ny);
[Rxx, nxx] = conv_m(x1, nx1, x, nx);

figure('NumberTitle', 'off', 'Name', 'Problem 2.10 original sequence')
set(gcf,'Color',[1,1,1])                  % 改变坐标外围背景颜色

subplot(2,1,1); stem(nx, x); title('x sequence');
xlabel('nx'); ylabel('x(n)') ;
grid on
subplot(2,1,2); stem(ny, y); title('y sequence');
xlabel('ny'); ylabel('y(n)');
grid on;


figure('NumberTitle', 'off', 'Name', 'Problem 2.10 correlation sequence')
set(gcf,'Color',[1,1,1])                  % 改变坐标外围背景颜色

subplot(2,1,1); stem(nxy, Rxy); title('Crosscorrelation Rxy(l) sequence');
xlabel('nxy'); ylabel('Rxy(nxy)') ;
grid on
subplot(2,1,2); stem(nxx, Rxx); title('Autocorrelation Rxx(l) sequence');
xlabel('nxx'); ylabel('Rxx(nxx)');
grid on;

　　运行结果：