% Pas temporel égal à 0.0001
dt=1e-4;
% Fréquence d'échantillonnage Fs=10000Hz
Fs=1/dt;
%
% Création d'un vecteur temps contenant 10001 valeurs
% 
temps=0:dt:1;
% Un signal utile modulant possédant quatre composante fréquentielle entre
% 100 Hz et 400 Hz, et avec des amplitudes maximales respectivement égales
% à 1, 0.8, 0.5 et 0.3
% Ce signal possède 10001 échantillons, car il utilise le vecteur temps
%
xbf=[1 .8 .5 .3]*cos(2*pi*[100 200 300 400]'*temps); 
% Affichage du signal xbf sur 200 points uniquement
figure;
plot(temps(1:200),xbf(1:200));
title('Signal Original dans le domaine temporel')

Nw=length(temps);
df=Fs/Nw;
freq=-(Fs-df)/2:df:(Fs-df)/2;
% Porteuse de fréquence 1000 Hz
x0=cos(2*pi*1000*temps);
% Modulation sans porteuse
x=x0.*xbf;
%affiche le spectre
figure;
Nw=length(temps);
df=Fs/Nw;
%
% Calcul du module du spectre du signal utile modulant xbf en utilisant la FFT
%
X=abs(fft(xbf,Nw))/Nw;
X=fftshift(2*X);
%
% Calcul du module du spectre du signal modulé x en utilisant la FFT
%
S=abs(fft(x,Nw))/Nw;
S=fftshift(2*S);
%
% Traçage des deux spectres, du signal modulant et du signal modulé, sur un
% même repère
%
%
plot(freq,S,'b',freq,X/2,'r');
legend('spectre S(f)','spectre X_bf')
grid;
axis([0 2000 0 .5]);
title('Translation de fréquence')
xlabel('fréquence (Hz)')
ylabel('échelle linéaire')