% esegui il seguente script per generare il grafico di 
% funzioni pseudo-periodiche del tipo y= exp(-t/tau).*sin(omega*t)
% con diversi valori di tau e di omega

clear all
close all
clc

% definisci tau ed omega
tau=5;
omega=3;

% definisci il vettore dei tempi con un passo opportuno (per visualizzare
% correttamente le funzioni exp e sin)

t_0=0; %estremo inf
t_f=20; % estremo sup
t_s=0.01; % passo
% definisci il vettore dei tempi t
t=t_0:t_s:t_f;


y_exp=exp(-t/tau);
y_sin=sin(omega*t);
% definisci y
y=y_exp.*y_sin;

figure
hold on
xlabel('time')
ylabel('y')
plot(t,y,'-r','linewidth',2)

% definisci funz. exp. con due diverse costanti di tempo
tau_a=3;
tau_b=6;
y_exp_a=exp(-t/tau_a);
y_exp_b=exp(-t/tau_b);

% definisci funz. sin. con due pulsazioni diverse
omega_1=1;
omega_2=2;
y_sin_1=sin(omega_1*t);
y_sin_2=sin(omega_2*t);


figure
subplot(2,2,1)
% plot le funz. exp.
hold on
grid on
plot(t, y_exp_a,'-b');
plot(t, y_exp_b,'-r');
ylabel('y')
%xlabel('t')
legend('y_1=e^{-t/3}','y_2=e^{-t/6}')

% plot le funz. sin
subplot(2,2,2)
hold on
grid on
plot(t, y_sin_1,'-k');
plot(t, y_sin_2,'-g');
legend('y_1=sin(t)','y_2=sin(2t)')
%xlabel('t')
%ylabel('y')

% plot exp(-t/3), exp(-t/6), sin(t), exp(-t/3).*sin(t)
% e exp(-t/6).*sin(t)
subplot(2,2,3)
hold on
grid on
y_pseudo_1a=y_exp_a.*y_sin_1;
y_pseudo_1b=y_exp_b.*y_sin_1;

plot(t, y_exp_a,'--b');
plot(t, y_exp_b,'--r');
plot(t, y_sin_1,'--k');

plot(t,y_pseudo_1a,'-b','linewidth',3)
plot(t,y_pseudo_1b,'-r','linewidth',3)

legend('y_1=e^{-t/3}','y_2=e^{-t/6}','y_3=sin(t)','y_4=y_1 y_3','y_5=y_2 y_3');
xlabel('t')
ylabel('y')


% plot exp(-t/3), exp(-t/6), sin(2t), exp(-t/3).*sin(2t)
% e exp(-t/6).*sin(2t)
subplot(2,2,4)
hold on
grid on
y_pseudo_2a=y_exp_a.*y_sin_2;
y_pseudo_2b=y_exp_b.*y_sin_2;

plot(t, y_exp_a,'--b');
plot(t, y_exp_b,'--r');
plot(t, y_sin_2,'--g');

plot(t,y_pseudo_2a,'-c','linewidth',3)
plot(t,y_pseudo_2b,'-m','linewidth',3)
xlabel('t')
%ylabel('y')

legend('y_1=e^{-t/3}','y_2=e^{-t/6}','y_3=sin(2t)','y_4=y_1 y_3','y_5=y_2 y_3');