Tom Colinot, Buet Crampon Réverbération à plaque Acoustique Musicale
Exercice : Réverb à plaque avec Magpie
27 septembre 2024
[Jae and Smith, 1983]
Algorithme de Karplus-Strong, version standard [Karplus and Strong, 1983]
1. Soit l'équation aux dérivées partielles régissant la position h(x, t) d'une corde sans raideur :
∂2h 2
2∂ h
− c = 0.
∂t2 ∂x2
2. Si l'onde progressive vaut h+ (0, t) à l'abscisse x = 0, que vaut-elle en x = ℓ ? Si l'onde
régressive vaut h− (ℓ, t) en x = ℓ, que vaut-elle en x = 0 ?
3. On considère la condition aux limites suivante en x = ℓ :
∂h ∂h
−T (ℓ, t) = a (ℓ, t)
∂x ∂t
Si une onde pure progressive h+ = H + ej(ωt−kx) donnée arrive sur cette terminaison, que
vaut l'onde régressive h− = H − ej(ωt+kx) qui s'y ré échit ? On rappelle que h = h+ + h− .
4. La méthode de synthèse dite de Karplus-Strong décompose la corde en un retard pur suivi
d'un ltre passe-bas. Identi ez ces deux éléments dans la fonction de transfert précédente.
5. Le ltre passe-bas proposé par Strong, de fonction de transfert
1 1 −
Ha (z) = + z 1. (1)
2 2
ne s'appuie pas sur des considérations liées à la physique de la corde ou de son extrémité,
mais plutôt sur des arguments purement numériques de coût de calcul. Par rapport à un
ltre plus standard du type
1
Ha′ (z) = , (2)
1 + a1 z−1 + a2 z −2
en quoi le ltre Ha permet-il de limiter le coût de calcul numérique ?
6. La fonction de transfert discrète associée au retard pur est la suivante :
Hb (z) = z −N (3)
où N est un nombre d'échantillons. Que doit valoir N en fonction de la célérité c, la longueur
ℓ et la fréquence d'échantillonnage Fs pour que le retard corresponde à un aller-retour sur
la corde ?
7. Quel type de nombre est N ? Quelle limitation cela implique-t-il sur la fréquence fondamen-
tale de la corde ? Indication : on pourra se baser sur la formule de fréquence fondamentale
de la corde f0 = c/(2ℓ).
8. Pour surmonter ce problème, Jae [Jae and Smith, 1983] propose un ltre passe-tout per-
mettant d'introduire un retard de propagation non-entier sans aecter l'amplitude de l'onde.
9. L'excitation proposée initialement par Karplus et Strong se base sur une séquence aléatoire
un prenant les valeurs 1 ou -1 avec une probabilité de 0.5 chacune. On dé nit la fonction
d'autocorrélation du processus un comme l'espérance suivante
Ru [k] = E(u[n]u[n − k]) (4)
Que vaut Ru [k] pour k = 0 ? Pour k ̸= 0 ? Quelle fonction reconnait-on ?
10. Donner la densité spectrale de puissance du signal u, dé nie comme la transformée de Fourier
de Ru
Su = F(Ru ) (5)
Master 2 Acoustique & Musicologie 1 27 septembre 2024
Tom Colinot, Buet Crampon Réverbération à plaque Acoustique Musicale
11. En réalité, le signal d'excitation x[n] est à support ni et vaut
x[n] = rectN [n]u[n] (6)
avec rectN la fenêtre rectangulaire qui vaut 1 pour n entre 0 et N, et 0 ailleurs. Réécrire la
fonction
Sx = F(Rx ) (7)
en fonction de Su et Srect .
Sx = F(E(u[n]rectN [n]rectN [n − k]u[n − k])) (8)
Sx = F(E(u[n]u[n − k])E(rectN [n]rectN [n − k])) (9)
Sx = Su ∗ Srect (10)
où ∗ désigne le produit de convolution.
2
12. On admet que ( ?) Srect [k] = sin(πk/N
πk/N
)
. Que vaut alors la densité spectrale de puissance
du signal d'excitation Sx ?
13. Jae [Jae and Smith, 1983] ajoute également un ltre passe-bas au signal d'excitation.
Quel paramètre de l'excitation peut être modélisé par cet ajout ?
14. A n de prendre en compte l'abscisse du point d'excitation, [Jae and Smith, 1983] propose
d'utiliser un ltre dit 'en peigne' de fonction de transfert
He (z) = 1 − z µN , (11)
avec µ un paramètre inférieur à 1. Que vaut la réponse en fréquence d'un tel ltre ? En
particulier, quels sont les fréquences où cette réponse s'annule ? Ont-elles un rapport avec
les harmoniques de la corde ? Indication : on rappelle que z = exp(−jωTs ).
Références
[Jae and Smith, 1983] Jae, D. A. and Smith, J. O. (1983). Extensions of the karplus-strong
plucked-string algorithm. Computer Music Journal, 7(2) :5669.
[Karplus and Strong, 1983] Karplus, K. and Strong, A. (1983). Digital synthesis of plucked-string
and drum timbres. Computer Music Journal, 7(2) :4355.
Master 2 Acoustique & Musicologie 2 27 septembre 2024
