% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Convoluzione di segnali continui a supporto limitato %%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Calcoliamo con Matlab una versione approssimata della
% convoluzione dei segnali continui
% f(t) = A_f * rect ((t- t_f )/ D_f )
% g(t) = A_g * rect ((t- t_g )/ D_g )

clc
clear all
close all

%% Segnali di interesse
% Parametri che definiscono i segnali
Af = 1;
Ag = 1;
tf = 1;
tg = 1;
Df = 1;
Dg = 1;
% Approssimazione discreta
T=0.001;
t=0:T:2;
f=Af*rectpuls(t-tf,Df);
g=Ag*rectpuls(t-tg,Dg);
% Attenzione a moltiplicare per il passo di campionamento
conv_fg=T*conv(f,g);
t_conv=2*t(1):T:2*t(end);

%% Grafici
subplot(221)
plot(t,f,'LineWidth',2);
xlabel('t')
ylabel('f(t)')
axis([xlim -.2 1.2])
grid on

subplot(222)
plot(t,g,'LineWidth' ,2);
xlabel('t')
ylabel('g(t)')
axis([xlim -.2 1.2])
grid on

subplot(212)
plot(t_conv,conv_fg,'b','LineWidth',2)
xlabel('t')
ylabel('f*g(t)')
axis([xlim -.2 1.2])
grid on