voir la Feuille de calcul:

frmax=NbPoint/t[NbPoint]	fréquence max , fonction des paramètres choisis
	NbPoint est choisi = 4096 pour la FFT
	t[NbPoint]=Durée Totale sera adaptée à la fréquence fm pour "voir" quelques périodes
m=200	taux de modulation (x 100)
fp=1000	porteuse
fm=50	information
x=T	variable temporelle
y=sin(fp*T*6.283)*(1.0000001+(m/100)*cos(fm*T*6.283))	signal modulé, ici il y a surmodulation (m/100=2) - fenètre 1
z=y*sin(fp*T*6.283)	apres la multiplication - fenètre 2
zfiltre=PBas(z,100)	après le filtre RC pour 100Hz = signal détecté - fenètre 4
fr=t*frmax/t[NbPoint]	variable fr=fréquence
Sfftz=FFTM(z)	FFT du signal après la multiplication - fenètre3
Sfftzf=FFTM(zfiltre)	FFT du signal après le filtre - fenètre 5
Il faut ajuster les abscisses des fenètres 3 et 5 au maximum à 3*fp ( ou 2*fp+fm)