N° 794 BULLETIN DE L'UNION DES PHYSICIENS 929
Modulation de fréquence et synthèse F.M.
par Daniel BEAUFILS
Institut National de Recherche Pédagogique
Département Technologies Nouvelles
91, rue Gabriel Péri - 92120 Montrouge
e-mail : beaufils@[Link]
RÉSUMÉ
Dans un précédent article nous avons présenté une introduction possible du
principe des modulations d’amplitude et de fréquence fondée sur leurs utilisations
dans le domaine sonore et musical. Nous proposons ici quelques compléments
spécifiques des signaux obtenus par modulation de fréquence : la génération et l’étude
spectrale des signaux sont effectuées en référence au principe de la «synthèse F.M.»
qui est utilisée en synthèse sonore (sons harmoniques et inharmoniques) et est à la
base des synthétiseurs.
Dans un précédent article [1] nous avons présenté une introduction possible du
principe des modulations d’amplitude et de fréquence fondée sur leurs utilisations dans
le domaine sonore et musical. Nous proposons ici quelques compléments concernant
le cas particulier de la synthèse par modulation de fréquence. Cette «synthèse F.M.»,
ainsi nommée en création musicale, a été inventée en 1967 dans le département de
Musique de Standford par John CHOWNING. Les éléments que nous donnons ci-dessous,
s’appuient sur quelques ouvrages dont en particulier celui que J. CHOWNING lui-même
a écrit avec D. BRISTOW1 [2] dans le cadre de leur contribution à la conception des
synthétiseurs Yamaha©.
1. L’ouvrage de J. CHOWNING et D. BRISTOW est intéressant dans la mesure où il s’adresse à
des non-spécialistes et contient toutes les explications, depuis la définition de la notion de
fréquence jusqu’aux tableaux des fonctions de Bessel ; de nombreux exercices d’applica-
tion sont proposés, mais il faut connaître le DX7 (ou équivalent) pour suivre les consi-
gnes ; par ailleurs, il est malheureusement épuisé et ne semble pas avoir fait l’objet d’une
référence d’édition particulière.
Vol. 91 - Mai 1997 D. BEAUFILS
930 BULLETIN DE L'UNION DES PHYSICIENS
LA MODULATION DE FRÉQUENCE
Le principe théorique est la modulation de la fréquence d’un signal (dit «onde
porteuse» dans le cas des télécommunications) par un signal de plus basse fréquence
(signal modulant). La relation mathématique représentant l’onde modulée s’écrit
généralement :
c a
s = s0 .sin 2 p fp t + [Link] 2 p fm t fh (1)
Cette expression, où β est appelé indice de modulation, se comprend en calculant
la fréquence «instantanée» comme la dérivée temporelle de la phase ; on obtient bien
en effet la relation donnant la modulation de la fréquence de base par la fonction de
fréquence fm(2)
a
f = fp + fm . [Link] 2 p fm t f (2)
Une telle relation peut être utilisée pour générer, non seulement la fonction
mathématique et son graphe, mais également un signal acoustique si les fréquences de
la porteuse et de la modulante sont du domaine audible.
a f
Figure 1 : Exemple d’un signal modulé en fréquence (indice égal à 4) :
sin 2 * pi * f * t+4 *sin ( 2 * pi * f * t / 10 ) ce signal est obtenu à partir de l’expression mathématique (édi-
teur syntaxique du logiciel GoldWaves© et peut être écouté si l’on dispose d’un ordinateur muni d’une
«carte son» munie d’enceintes acoustiques.
2. Cette relation permet de relier l’indice de modulation à l’excursion temporelle en fréquence :
b = Df/fm.
Modulation de fréquence et synthèse F.M. B.U.P. n° 794
BULLETIN DE L'UNION DES PHYSICIENS 931
LE SPECTRE D’UN SIGNAL MODULÉ EN FRÉQUENCE
Le calcul du spectre de Fourier d’un signal défini par la relation (1) présente quelques
difficultés mathématiques ; celui-ci conduit en effet à des fonctions intégrales : les
fonctions de Bessel3 (notées Jk(x)). De façon plus précise, le spectre est donné par le
développement :
c a
sin 2 p fp t + [Link] 2 p fm t = fh
af c h
J 0 b .sin 2 p fp t +
J a b [Link] c 2 p f t + 2 p f t h – J a b [Link] c 2 p f t – 2 p f t h +
1 p m 1 p m
J a b [Link] 2 p f t + 4 p f t h – J a b [Link] 2 p f t – 4 p f t h +
2 p m 2 p m
J a b [Link] 2 p f t + 6 p f t h – J a b [Link] 2 p f t – 6 p f t h+K
3 p m 3 p m
Soit, sous forme d’un tableau de «coefficients d’amplitude» :
Fréquence (k = 0, 1, ...) «Amplitude»
fp + k. fm J k ( b)
fp – k. fm ( – 1)k . J k ( b )
Tableau 1
On obtient donc un spectre centré sur la fréquence de la porteuse et comportant
des composantes espacées d’une valeur égale à la fréquence modulante. Mais il
convient de faire deux remarques. La première concerne l’étendue du spectre qui,
contrairement à la modulation d’amplitude, n’est pas limité : le nombre des fréquences
latérales (inférieures et supérieures) est infini. Toutefois, compte tenu du comportement
des fonctions de Bessel, on peut définir une largeur «effective» : pour une valeur
donnée de l’indice de modulation les amplitudes des composantes sont en effet
négligeables au-delà d’une certaine valeur de k. Dans la pratique, on considère les
valeurs |k| inférieures ou égale à b + 2.
La seconde remarque concerne les amplitudes relatives : l’importance de la
fréquence centrale (celle de la porteuse) dépend fortement de l’indice de modulation
: on voit facilement que pour un indice égal à 2.4, l’amplitude de cette composante est
3. Voir compléments en annexe.
Vol. 91 - Mai 1997 D. BEAUFILS
932 BULLETIN DE L'UNION DES PHYSICIENS
quasiment nulle, et que l’essentiel du spectre est donné par deux ou trois fréquences
latérales.
Précisons sur un exemple. Considérons pour cela une fréquence porteuse de
1000 Hz et une fréquence de modulation de 100 Hz. Nous allons calculer le spectre
pour trois valeurs de l’indice de modulation : b = 1, 2.4 et 4. Pour cela, nous allons
dresser le tableau 2 des valeurs des amplitudes des principales composantes.
«Amplitude» «Amplitude»
b k de 0 à + N J k(b )
(fréq. sup.) (fréq. inf.)
J0(1) 0.76 0.76
J1(1) 0.44 – 0.44
1 k de 0 à + 3
J2(1) 0.11 0.11
J3(1) 0.02 – 0.02
J0(2.4) 0 0
J1(2.4) 0.78 – 0.78
2.4 k de 0 à + 4 J2(2.4) 0.45 0.45
J3(2.4) 0.13 – 0.13
J4(2.4) 0.05 0.05
J0(4) – 0.39 – 0.39
J1(4) – 0.07 0.07
J2(4) 0.36 0.36
4 k de 0 à + 6 J3(4) 0.43 – 0.43
J4(4) 0.28 0.28
J5(4) 0.13 – 0.13
J6(4) 0.05 0.05
Tableau 2
Dans ce tableau, nous avons donné des amplitudes avec le signe attribué par le
calcul. Ce signe correspond à un déphasage de p dont nous ne tenons pas compte pour
obtenir les représentations spectrales simples correspondantes.
Figure 2
Modulation de fréquence et synthèse F.M. B.U.P. n° 794
BULLETIN DE L'UNION DES PHYSICIENS 933
Ci-dessous (figure 3), nous donnons le résultat de l’analyse spectrale par F.F.T.
d’un son créé à partir de l’expression mathématique (indice de modulation égal à 4 ;
voir figure 1).
Figure 3 : Calcul de la F.F.T. et représentation graphique (logiciel Virtual Waves©).
LA BASE DE LA SYNTHÈSE DE SONS MUSICAUX
La structure des spectres des signaux obtenus par modulation de fréquence montre
immédiatement l’intérêt pour la synthèse de sons riches : contrairement à la synthèse
additive, un nombre élevé d’harmoniques peut être facilement obtenu avec un seul
générateur. Pour obtenir un tel son il suffit de choisir une fréquence de porteuse qui
soit un multiple de la fréquence de modulation ; le spectre est alors semblable à celui
d’un son d’instrument : fondamentale (éventuellement absente [3]) et harmoniques
(fréquences multiples de l’intervalle).
Mais bien évidemment, le cas général où aucune relation particulière n’existe entre
la porteuse et la modulante peut être exploité pour la création de sons inharmoniques4.
Les sonorités ainsi obtenues sont très riches, nouvelles, parfois inouïes et ont été une
partie du succès de la synthèse F.M. Bien évidemment, l’un des exemples de synthèse
aujourd’hui banal est celui des sons de cloche. Ceux-ci peuvent en effet être très
éloignés d’un son harmonique, et la synthèse F.M. permet de les générer très
facilement ; il convient seulement, pour donner la «vie» au son de jouer sur la variation
de l’amplitude au cours du temps : amplitude de la porteuse qui tend vers zéro, mais
4. On retrouve alors le fait que le signal résultant n’est périodique (donc décomposable en
série de Fourier) que si les fréquences sont commensurables.
Vol. 91 - Mai 1997 D. BEAUFILS
934 BULLETIN DE L'UNION DES PHYSICIENS
aussi amplitude de la modulante qui peut évoluer et donc entraîner une variation de
timbre5.
L’étendue des possibles est encore plus grande, puisque rien n’interdit d’utiliser
de très forts indices de modulation, de même que des fréquences de modulation plus
élevées que celle de base. Bien évidemment, nous quittons ici le domaine intéressant
la communication6 puisque le résultat repose en particulier sur les effets de combinaison
de composantes avec des fréquences «négatives». Pour calculer le spectre il faut en
général bien tenir compte des signes des coefficients d’amplitudes.
Prenons un nouvel exemple : celui d’une modulation avec un indice égal à 1.5,
mais avec une fréquence de base de 100 Hz et une modulation à 400 Hz. Les
coefficients sont donnés par les valeurs des fonctions de Bessel Jk(1.5), avec k = 0, 1,
2, 3.
f = fp + fm Coefficient f = fp – fm Coefficient
Jk(b) Jk(b)
(Hz) d’amplitude (Hz) d’amplitude
100 J0(1.5) = 0.51
500 J1(1.5) = 0.56 – 1100 J– 1(1.5) = – 0.56
900 J2(1.5) = 0.24 – 700 J– 2(1.5) = 0.24
1300 J3(1.5) = 0.05 – 300 J– 3(1.5) = – 0.056
Tableau 3
Dans ce tableau 3, nous avons fait apparaître non seulement l’amplitude avec un
signe, comme précédemment, mais également les composantes correspondant à des
«fréquences négatives».
5. Contrairement à ce qui est sous-entendu (quand ce n’est pas écrit !) dans de nombreux
manuels scolaires, l’évolution temporelle du spectre est essentielle à la perception du tim-
bre ; elle est même inhérente à celle de «son musical» ! On notera à ce propos la remar-
que sur la nécessité de prendre en compte le régime transitoire des oscillateurs pour une
étude sérieuse de l’acoustique (voir document d’accompagnement des programmes de Ter-
minale S, version de janvier 1995, pages 72-73).
6. Il convient donc d’éviter ici le terme de «porteuse».
Modulation de fréquence et synthèse F.M. B.U.P. n° 794
BULLETIN DE L'UNION DES PHYSICIENS 935
Pour obtenir le spectre d’am- Figure 4
plitude (qui donnera le spec-
tre en énergie), il faut alors
tenir compte de ces signes
qui correspondent à des dé-
phasages de p.
De sorte, que pour “rame-
ner” les composantes de fré- Figure 5
quence négative dans la par-
tie “réelle” du spectre, il faut
effectuer le changement de
signe correspondant, addi-
tionner algébriquement les
composantes qui peuvent se Figure 6
recouvrir (ce n’est pas le cas
ici) avant de faire la repré-
sentation en valeurs abso-
lues.
Ci-dessous, le spectre obtenu (figure 7) par calcul numérique après la synthèse
d’un tel signal (Logiciel Virtual Waves ©).
Figure 7
On voit donc ici un exemple d’obtention d’un son harmonique ne comportant que
des harmoniques impaires. Mais on comprend alors également qu’un «désaccord» (en
utilisant une fréquence de base différente) entraînera la disparition du caractère
Vol. 91 - Mai 1997 D. BEAUFILS
936 BULLETIN DE L'UNION DES PHYSICIENS
harmonique et donc un changement du timbre, et pourra dans certains cas faire apparaître
des battements. Les spectres ci-dessous (figures 8 à 10) ont ceux de tels signaux.
Même modulation mais avec une fréquence «porteuse» de 150 Hz.
Figure 8
Même modulation mais avec une fréquence «porteuse» de 190 Hz.
Figure 9
Forme d’onde du dernier signal synthétisé : on y voit à la fois la variation de la
fréquence et l’effet sur l’amplitude (perçu comme des battements).
Figure 10
Modulation de fréquence et synthèse F.M. B.U.P. n° 794
BULLETIN DE L'UNION DES PHYSICIENS 937
LA MODULATION PAR UN SIGNAL COMPLEXE : LA «SYNTHÈSE F.M.»
Dans les paragraphes précédents nous n’avons présenté que le cas particulier d’une
modulation par un signal simple : un signal harmonique ne comportant donc qu’une
seule fréquence. Bien souvent les explications s’arrêtent là, alors que bien évidemment,
la modulation n’a d’intérêt que pour des signaux modulant complexes ! Nous avons
détaillé l’exemple de la transmission d’une parole dans le cas de la modulation
d’amplitude7, et nous voulons ici présenter le cas de la modulation de fréquence pour
la synthèse de sons musicaux.
Nous faisons ici référence aux «algorithmes»8 implantés dans les synthétiseurs
type DX7© (Yamaha). Un module (ou algorithme) de synthèse F.M. est une association
particulière de six oscillateurs, dont un au moins est modulé en fréquence. L’objectif
n’est évidemment pas d’en faire une présentation détaillée (il y a trente-deux
algorithmes classiques) mais de donner quelques indications sur leur nature et sur les
utilisations simples. Les exemples ci-dessous sont les algorithmes numérotés 3, 5 et
18 chez Yamaha :
Dans l’algorithme a, par exemple, les trois oscilla-
teurs 1, 3 et 5 sont associés en synthèse additive, et
sont respectivement modulés par les oscillateurs 2, (a)
4 et 6. La boucle sur l’oscillateur 6 indique qu’une (c)
rétroaction est possible, permettant ainsi de géné-
rer des signaux très complexes.
Dans l’algorithme b, les sorties de deux oscilla-
teurs (1 et 4), modulés par des oscillateurs (2 et 5)
eux mêmes modulés (par 3 et 6), sont ajoutées en
sortie. (b)
Dans l’algorithme c, la sortie est celle d’un seul
oscillateur (1) mais modulé par trois oscillateurs Figure 11
complexes (2, 3 et 4).
Pour chacun des six oscillateurs, il est bien évidemment possible de choisir la
fréquence, l’amplitude, l’indice de modulation ou de rétroaction, ainsi que la durée
totale et l’enveloppe temporelle. Comme exemple simple, il peut être intéressant de
montrer l’évolution spectrale liée à la simple variation (ici diminution) de l’amplitude
du signal modulant.
7. Article précédemment cité.
8. Le terme est ici mis entre guillemets puisqu’il s’agit de la terminologie Yamaha ; il fait
référence aux méthodes implantées dans ce synthétiseur numérique.
Vol. 91 - Mai 1997 D. BEAUFILS
938 BULLETIN DE L'UNION DES PHYSICIENS
Figure 12 : Module de synthèse, ana- Figure 13 : Signal de référence :
lyseur, et sortie. Pour cet exemple on oscillateur n° 1 : 1000 Hz, le signal monte très rapidement, est
n’utilisera qu’une partie «d’algo- maintenu durant près d’une seconde puis s’éteint rapidement.
rithme» pour moduler un oscillateur
(n° 1) par un second oscillateur (n° 2)
(Logiciel Virtual Waves©).
Figure 14 : Signal modulant : Figure 15 : Spectre tridimensionnel
oscillateur n° 2 (400 Hz) d’amplitude constante. (Logiciel Virtual Waves©).
Dans le spectre du signal obtenu, on voit bien l’effet de la modulation de
fréquence : plusieurs fréquences latérales, faible amplitude relative de la fréquence
«porteuse». Dans l’exemple qui suit, seule l’amplitude du signal de modulation a été
modifiée par rapport au cas précédent : elle diminue au cours du temps. Dès lors, l’effet
est nettement visible sur le spectre correspondant : à l’extinction des fréquences
latérales, correspond une augmentation progressive du signal «porteur». En terme
Modulation de fréquence et synthèse F.M. B.U.P. n° 794
BULLETIN DE L'UNION DES PHYSICIENS 939
«sonore», on voit donc apparaître une enveloppe croissant progressivement pour la
fréquence centrale et on entend un son initialement riche et dont le timbre évolue vers
un son dit «pur».
Figure 16 : Signal modulant : Figure 17 : L’augmentation progressive de
forte amplitude initiale et diminution rapide. l’amplitude de la «porteuse» n’est pas dans le signal
n° 1 : elle est due à la variation de la modulation de
fréquence.
On comprend alors que, devant l’étendue des possibles ouverte par le jeu des
paramétrages, s’arrête ici le domaine des explications physico-mathématiques et
commence le rôle de la pratique, de l’intuition et de l’inspiration propre à la création
musicale.
Spectre d’un son complexe synthétisé en
un seul calcul de «synthèse F.M.» : type
son de cloche
(logiciel Virtual Waves©).
Figure 18
Notons enfin que, tous les logiciels dits «audionumériques» permettent la
sauvegarde des résultats sous forme de fichiers de type «Wave» (.WAV pour compatible
PC) avec la fréquence d’échantillonnage de son choix (généralement 22 050 Hz ou
Vol. 91 - Mai 1997 D. BEAUFILS
940 BULLETIN DE L'UNION DES PHYSICIENS
44 100 Hz). Ils peuvent alors être reproduits par tout logiciel gérant la sortie audio
d’une carte sonore.
BIBLIOGRAPHIE
[1] D. BEAUFILS : «Modulation d’amplitude et modulation de fréquence : les avantages
du domaine sonore», B.U.P. n° 793, avril 1997, p. 741-749.
[2] J. CHOWNING et D. BRISTOW : «F.M. ; théorie et applications (par des musiciens
et pour des musiciens)», Yamaha (1985 ?).
[3] J.-P. LE CARDONNEL : «A propos des sons périodiques sans fondamental», B.U.P.
n° 767, juillet-août-septembre 1994, pp. 1397-1402.
[4] F. BROWN : «La musique par ordinateur», collection Que-Sais-Je ?, n° 2011, PUF,
1982.
[5] C. CANCE : «Analyse spectrale par transformée de Fourier», in actes de
l’Université d’été Outils informatiques d’investigation scientifique pour l’ensei-
gnement des sciences physiques, Paris : UdP-INRP, 1995, pp. 117-127.
[6] J. CHOWNING : «The synthesis of complex audio spectra by means of frequency
modulation», in Fondation of computer music, C. Roads et J. Strawn Eds.,
Cambridge, Massaschussets : MIT Press, 1985.
[7] J. ESQUIEU : «Analyse spectrale numérique», B.U.P. n° 754, avril 1993, pp. 775-
797.
[8] N. FOURNEL : «Le son par l’image», Ordinateur et Musique, n° 19, 1996,
pp. 10-13.
[9] R. KERAVEC et R. MOREAU : «Moduler en fréquence, wobuler : morceaux
choisis», B.U.P. n° 772, mars 1995, pp. 441-458.
[10] H. REINHARD : «Équations différentielles», Dunod Université, 1989, p. 452.
[11] S. SAUNDERS : «Improved F.M. audio synthesis methods for real-time digital
music generation», in Fondation of computer music, C. Roads et J. Strawn Eds.,
Cambridge, Massaschussets : MIT Press, 1985.
[12] R. TOURNIER : «Modulations», B.U.P. n° 766, juin 1994, pp. 1151-1164.
[13] B. VELAY : «Simulation de modulations de fréquence avec le logiciel SPICE»,
B.U.P. n° 772, mars 1995, pp. 459-470.
Modulation de fréquence et synthèse F.M. B.U.P. n° 794
BULLETIN DE L'UNION DES PHYSICIENS 941
LOGICIELS CITÉS
GOLDWAVES : logiciel en shareware, livré en démonstration sur certaines disquettes et
accessible par téléchargement à [Link]
VIRTUAL WAVES : logiciel français (auteur Nicolas Fournel), commercialisé par
Synoptic (Montreuil). Une présentation détaillée de la première version a été faite dans
les numéros de la revue Ordinateur et Musique (1996). Version de démonstration
téléchargeable à : [Link]
REGRESSI : logiciel distribué par Micrélec. Il permet notamment d’obtenir facilement
le tracé de fonctions définies par une intégrale (et non par une primitive).
NOTE
Les illustrations utilisées sont des copies d’écrans partielles des fenêtres des
logiciels : celles-ci contiennent généralement de nombreux paramètres de configura-
tion qui auraient surchargé la présentation sans apporter d’information utile à notre
propos.
Vol. 91 - Mai 1997 D. BEAUFILS
942 BULLETIN DE L'UNION DES PHYSICIENS
Annexe
Quelques compléments sur les fonctions de Bessel
DÉFINITION DES FONCTIONS DE BESSEL
Les fonctions de Bessel se rencontrent peu souvent, hormis le cas bien connu de
la fonction J1 qui intervient dans le calcul de la figure de diffraction par une pupille
circulaire. Dans les ouvrages de mathématiques (voir [10] par exemple), ces fonctions
sont introduites comme les solutions de l’équation différentielle :
c
x 2 y" + xy' + x 2 – n2 y = 0 h
Celles qui nous intéressent ici sont en fait les fonctions dites de Bessel de première
espèce, et seulement celles pour lesquelles n est un entier (n). Ces fonctions sont
notées Jn. Dans ce cas, l’une des propriétés est celle implicitement utilisée dans l’article
:
J – n ( x) = (– 1)n J n ( x )
L’expression intégrale des fonctions Jn peut s’écrire de différentes façons :
Jn =
1
p z p
0
cos ([Link]( t ) – n.t ). dt
ou : J2n =
1
p zp
0
cos ([Link]( t )).cos (2 n.t ). dt
J 2 n+1 =
1
p z
0
p
sin ([Link]( t )).sin ((2 n +1). t ). dt
Ces expressions permettent de voir le lien avec la décomposition de la fonction
modulée en fréquence donnée dans l’article.
GRAPHES DE QUELQUES FONCTIONS DE BESSEL
Ci-après nous donnons, afin de permettre une approche semi-quantitative, le
graphe des fonctions de Bessel J0, J1, J2 et J3. Ces graphes ont été obtenus par calcul
d’intégrale suivant l’une des formules ci-dessus avec le logiciel Regressi© : la variable
Modulation de fréquence et synthèse F.M. B.U.P. n° 794
BULLETIN DE L'UNION DES PHYSICIENS 943
x ayant été créée avec ses valeurs successives (commande Remplir), une variable telle
que J1 peut être définie comme fonction par :
a f
(1 / p) * INTG 0, p, t, cos (x * sin( t ) – t )
Figure 19 : Graphe de la fonction J0.
Figure 20 : Graphe de la fonction J1.
Vol. 91 - Mai 1997 D. BEAUFILS
944 BULLETIN DE L'UNION DES PHYSICIENS
Figure 21 : Graphe de la fonction J2.
Figure 22 : Graphe de la fonction J3
Ces fonctions sont visualisables à l’aide de logiciels tels que VirtualWaves cité
dans le texte. Ci-dessous, le spectre en trois dimensions d’un signal dont on fait varier
le taux de modulation de 0 à 400 % de façon linéaire dans le temps. Cela revient à
faire varier à la manière d’une vobulation, l’indice de modulation de 0 (à l’instant
Modulation de fréquence et synthèse F.M. B.U.P. n° 794
BULLETIN DE L'UNION DES PHYSICIENS 945
initial) à 4 (date finale choisie égale à 400 ms). Les courbes d’enveloppe temporelle
qui apparaissent pour la fréquence centrale et les bandes latérales ne sont autre que les
fonctions de Bessel ci-dessus (du moins leur valeur absolue, puisqu’il s’agit
uniquement des amplitudes) : on voit en particulier l’annulation de la fréquence
centrale pour t = 240 ms environ.
Figure 23 : Calcul et représentation de 0 à 240 %.
Figure 24 : Calcul et représentation complète de 0 à 400 %.
Vol. 91 - Mai 1997 D. BEAUFILS