% Program 2_8
% Computation of Autocorrelation of a
% Noise Corrupted Sinusoidal Sequence
%
N = 96;
n = 1:N;
x = cos(pi*0.25*n); % Generate the sinusoidal sequence
d = rand(1,N) - 0.5; % Generate the noise sequence
y = x + d; % Generate the noise-corrupted sinusoidal sequence
r = conv(d, fliplr(d)); % Compute the correlation sequence
k = -28:28;
stem(k, r(68:124));
xlabel('Lag index'); ylabel('Amplitude');