% PART 1
% Càlcul de la TF d'un pols rectangular discret
fm=50;
% Freqüència de mostreig
T=90e-3;
% Durada del pols analògic en segons
% Generació del pols discret
L=T*fm;
% Longitud del pols discret
n = 0:L;
% Vector d'índexs temporals del senyal pols discret xd
xd=1.*(n>=0)&(n<L);
% Càlcul dels N punts d'observació de la TF del pols discret
N=1001;
F=linspace(-2,2,N);

Xd = xy_trF(xd, n, F);
% !!!!!Introduiu aqui el codi necessari per generar el vector Xd

% Dibuix transformada
figure
subplot(3,1,1);
stem(n,xd);
title('Pols digital causal');
xlabel('n');
subplot(3,1,2);
plot(F,abs(Xd));
title('Mòdul TF');
xlabel('F');
subplot(3,1,3);
plot(F,angle(Xd));
title('Fase TF');
xlabel('F');
% PART 2
% Comparació de TF de pols rectangular analògic i discret
% Càlcul dels N punts d'observació de la TF del pols analògic
N=1001;
f=linspace(-2*fm,2*fm,N);
Xa=T*exp(-j*2*pi*f*T/2).*sinc(T*f);
% Dibuix comparat
figure
subplot(4,1,1);
stem(n,xd);
title('Pols digital causal');
xlabel('n');
subplot(4,1,2); plot(f,abs(Xa));
title('Mòdul TF senyal analògic |Xa(f)|');
xlabel('f(Hz)');

subplot(4,1,3);
plot(F,abs(Xd),'r');
title('Mòdul TF senyal digital |Xd(F)|');
xlabel('F');
subplot(4,1,4);plot(fm*F,abs(Xd)/fm,'r',f,abs(Xa),'b');
title ('Comparació de |Xd(f/fm)|/fm (vermell) i |Xa(f)|(blau)');
xlabel('f(Hz)');