% Campiona il segnale x(t)=sin(t*2*pi*3) con frequenza di campionamento
% a scelta ottenendo il segnale campionato xn.
% Mostra l'effetto dell'alias in frequenza quando non si rispetta il teorema 
% del campionamento che dice che fs>6.
% Si provi a vedere cosa succede se si utilizza fs=24, fs=12, fs=4, fs=3.


clear all, close all, clc

%% Parametri
BeginTime = 0;
EndTime  = 1;
dt = 0.001;
TrueF = 3;

%% Esecuzione Programma
t=(BeginTime:dt:EndTime); % Array dei tempi
k=length(t); % vedi significato inglese
x(1:k)=sin(2*pi*TrueF*t); % Costruisco un vettore con al suo interno i campioni della sinusoide
figure(1), % Apro la figura 1
plot(t,x)
grid on
xlabel('time')
title('sinusoid frequency = 3Hz')
hold on
fs=input('frequenza di campionamento fs? '); % Stampo a video e aspetto la risposta, che salvo in fs
T=1/fs;
nT=(BeginTime:T:EndTime); % Vettore Tempo Campionato
kk=length(nT);
xn=sin(2*pi*TrueF*nT);
figure(1)
stem(nT,xn,'r','filled')
hold off