close all;
clear all;
clc;
fs = 1e3;
f = 2e2;
D = 0.3;
StopTime = 1;
t = (0 : 1/fs : StopTime-1/fs)';
N = size(t,1);
w = 0.2;
yp = pulstran(t,D,@rectpuls,w);
x = (10)*sin(2*pi*f*t);
display(size(yp));
display(size(x));
z = yp.*x;
noisy_z = z + 0.5*randn(size(z));
y = xcorr(noisy_z,'biased');
display(size(y));
subplot(2,1,1);
plot(t*(1/(f*StopTime)),z);
title('Pulsed sine signal');
subplot(2,1,2);
plot(t*(1/(f*StopTime)),noisy_z);
title('Noisy signal');
l = length(t);
X1 = fft(y,l);
X2 = X1(1:l/2);
X = abs(X2)/N;
dF = fs/N;
f1 = 0:dF:fs/2-dF;
display(size(f1));
display(size(X));
figure;
plot(f1,(10*log10((X).^2*1e3)));
xlabel('Frequency (in hertz)');
ylabel('Power (in dbm)');
title('Autocorrelation spectrum');
grid on;
