• function s(n)
%
% SF Calcule la série de Fourier correspondant au signal carré
% et affiche la somme de la série jusqu'à l'harmonique de
% rang n
%
% n : nombre d'harmoniques
%
close all;
m=floor(n/2);
Ts=0.22e-3; % Période d'échantillonnage
t0=24.e-3; % Période du signal carré
N=500; % nombre de points
p=(-m-1):m ;
cn=((-1).^p)./(2*p+1)/pi;
cc=(2*pi*i*Ts)*(2*p+1)/t0;
for k=0:N-1
vn=exp(cc*k);
s(k+1)=cn*vn';
end
t=(0:N-1)*Ts;
f=(0:N-1)/N/Ts;
subplot(211); plot(t,real(s)); grid; title('signal')
subplot(212); plot(f,abs(fft(s)/N)); grid; title('signal')