%calcolo di alcuni segnali notevoli
% rihiede la funzione delta_unitaria.m (nella stessa cartella di lavoro)
close all
clear all;
fs = 1000;  % Number of samples per second
tMin=-1;
tMax=1;
t = tMin:1/fs:tMax;                      % Time vector
Tau=0.5
A=1
% rettangolare
x1=A*(abs(t)<Tau/2)+0.5*(abs(t)==Tau/2);
figure(1)
plot(t,x1)         % Plotting the pulse
axis([tMin tMax min(x1) max(x1)*1.1])   
xlabel('Time (sec)');
ylabel('Amplitude');

%triangolare
x2=A*(1-abs(t/(Tau/2))).*(abs(t)<Tau/2);
figure(2)
plot(t,x2)         % Plotting the pulse
axis([tMin tMax  min(x2) max(x2)]*1.1)   
xlabel('Time (sec)');
ylabel('Amplitude');

% impulso unitario
x3=zeros(1,length(t));
[~,k]=min(abs(t));
%x3(k)=1/(t(2)-t(1));  % dirac
x3(k)=1;   %unitario
figure(3)
plot(t,x3)         % Plotting the pulse
axis([tMin tMax  min(x3) max(x3)*1.1])   
xlabel('Time (sec)');
ylabel('Amplitude');

%treno di impulsi unitari
T=0.25;
k=T;
x4=delta_unitaria(t); %delta in 0
  while k<=max(t)
      x4=x4+delta_unitaria(t-k)+delta_unitaria(t+k);
      k=T+k;
  end
figure(4)
plot(t,x4)         % Plotting the pulse
axis([tMin tMax  min(x4) max(x4)]*1.1)   
xlabel('Time (sec)');
ylabel('Amplitude');

%gradino
x5=zeros(1,length(t));
x5=(t>=0);
figure(5)
plot(t,x5)         % Plotting the pulse
axis([tMin tMax  min(x5) max(x5)*1.1])   
xlabel('Time (sec)');
ylabel('Amplitude');

%semirampa
x6=zeros(1,length(t));
x6=(-t+1).*(t/2>0);
figure(6)
plot(t,x6)         % Plotting the pulse
axis([tMin tMax  min(x6) max(x6)*1.1])   
xlabel('Time (sec)');
ylabel('Amplitude');

% exp
x7=A*((t)>=0).*exp(-6*t);
figure(7)
plot(t,x7)         % Plotting the pulse
axis([tMin tMax min(x1) max(x1)*1.1])   
xlabel('Time (sec)');
ylabel('Amplitude');