% experiment 1
%
fs=375; % fs>fN=2fB
T=1/fs;
T0=1/10; % Omega_0=MCF(200pi, 40pi, 60pi)
t=0:T:2*T0-T; % With fs=375 Hz, the period is N=75 samples covers two periods 
              % of the continuous time signal.
xp=3*cos(200*pi*t+2)-sin(40*pi*t)+2*cos(60*pi*t-pi/2);
plot(t,xp);title('a periodic signal');xlabel('t');
Xp=fft(xp);
fr=fs/length(t);
f=0:fr:fs-fr;
figure;stem(f,abs(Xp));title('Fourier representation of the signal');
xlabel('Hz');