%
% Script di esempio di campionamento di una funzione
% con campionatore 'reale' e campionatore 'ideale' con
% delta di Dirac
%

tb=0:0.1:10; % base dei tempi
Nc = length(tb);
xt=sin(tb); % funzione da campionare

figure, plot(tb,xt), % Sinusoide continua
title('x(t) = sin(t)'), xlabel('Tempo'), ylabel('x(t)') 

T=1; % Tempo di campionamento
k=0;
deltak = zeros(size(xt)); % Delta di Dirac

for t=1:Nc,
    if(tb(t) == k*T),
        Xkt(k+1)=sin(k*T); % effettuo il campionamento
        dT(k+1) = k;       % indice del campionamento
        deltak(t) = 1;     % Valore 1 della delta di Dirac (t=kT)
    k=k+1;
    end
end

figure, plot(tb,xt,'b--',dT,Xkt,'ro'), xlabel('Tempo')
ylabel('x(t) e x(kT)'), title('Sinusoide e campioni')
legend('x(t)','x(kT)')

figure, plot(dT,Xkt,'ro'), xlabel('Tempo')
ylabel('x(kT)'), title('Campioni della sinusoide')

xtdk = xt.*deltak; % Convoluzione

figure
subplot(211), plot(tb,xt,'b--'), 
xlabel('Tempo'), title('Sinusoide'), ylabel('x(t)')
subplot(212), plot(tb,deltak,'r:')
xlabel('Tempo'), title('Delta di Dirac'), ylabel('\delta(t)')

figure
plot(tb,xtdk,'m:',dT,Xkt,'ro')
xlabel('Tempo'), title('Funzione Campionata con \delta(t)'), ylabel('x(t)*\delta(t-kT)')

figure
subplot(211), plot(dT,Xkt,'ro')
xlabel('Tempo'), title('Funzione Campionata Reale'), ylabel('x(kT)')
subplot(212), plot(tb,xtdk,'m:',dT,Xkt,'ro')
xlabel('Tempo'), title('Funzione Campionata con \delta(t)'), ylabel('x(t)*\delta(t-kT)')