%%ex03

%a
figure(1)
f = 1;
T = 1;
Ta =  98 * T /100;
t = [2:Ta:(100*T-Ta)]';

x = sin(2*pi*t);

figure
Espetro(x,Ta);
title('Spectrum (a)')
[X,fx] = Espetro(x,Ta);



%b
f = 1; %gcd(5,6,7);
T = 1/f;
Ta = 50 * T/100;
t = [0:Ta:(50-Ta)]';

y = sin(10*pi*t)+cos(12*pi*t)+cos(14*pi*t-pi/4);

figure
Espetro(y,Ta);
title('Spectrum (b)')
[Y,fy] = Espetro(y,Ta);


%c
f = gcd(10,20);
T = 1/f;
Ta = 100 * T/100;
t = [0:Ta:(T*100)-Ta];

c = 10 + 14*cos(20*pi*t - pi/3) + 8*cos(40*pi*t + pi/2);


Espetro(c,Ta);
title('Spectrum (c)')
[C,fc] = Espetro(c,Ta);


% d
Fa = 1;
Ta =  1/Fa;
t=[0:Ta:5-Ta];
z=(square(2*pi*t)+1)/2;


Espetro(z,Ta);
title('Spectrum (d) (square wave [0,1])')
[Z,fz] = Espetro(z,Ta);



%e

%4

Reconstroi(X,fx);

Reconstroi(Y,fy);

Reconstroi(C,fc);

Reconstroi(Z,fz);