FR 2846768 A1 20040507 FR 0213682 A 20021031 L'invention a pour objet un procédé de simulation numérique d'une interaction non linéaire entre une source excitatrice et une onde dans un résonateur et peut s'appliquer, en particulier, à la synthèse numérique, en temps réel, d'un phénomène oscillant tel que le son émis par un instrument de musique fonctionnant en oscillations entretenues, comme un instrument à vent. On sait que la simulation numérique d'un phénomène auto-oscillant peut s'effectuer par discrétisation, dans le domaine temporel, d'équations constituant le modèle du phénomène tel qu'il se traduit physiquement. Dans ce cas, cependant, la synthèse sonore, en temps réel, est difficile à réaliser et c'est pourquoi, depuis quelques années, on a développé des méthodes du genre guide d'ondes numériques , basées sur un formalisme traitement du signal de la propagation des ondes dans les deux directions du corps de l'instrument. Habituellement, de telles méthodes utilisent la décomposition de la vibration au sein du résonateur sous forme d'une somme des ondes allerretour se propageant, respectivement, de la source à l'extrémité du résonateur et de cette extrémité à la source. Une telle méthode ne permet pas d'exprimer simplement le couplage non linéaire qui existe entre la source excitatrice et le résonateur de l'instrument et limite la paramétrisation physique des algorithmes de synthèse. De plus, une telle méthode limite le choix de la géométrie du corps du résonateur à des portions cylindriques de tubes et rend difficile la prise en compte de la conicité du résonateur, par exemple pour des instruments du type saxophone. L'invention a pour objet de résoudre ces inconvénients et de lever ces limitations grâce à une nouvelle méthode de simulation et de synthèse en temps réel d'un phénomène oscillant, applicable spécialement mais non limitativement aux instruments à vent. En particulier, l'invention a pour objet un procédé de simulation permettant de prendre en compte les processus physiques régissant le fonctionnement d'un instrument réel et dont la mise en #uvre numérique peut être particulièrement simple. En outre, à partir d'une méthode de base applicable à des instruments du type clarinette, à résonateur cylindrique, l'invention peut être adaptée à la simulation d'autres types d'instruments. Par ailleurs, l'invention ne se limite pas à la simulation d'instruments de musique mais peut s'appliquer, d'une façon générale à la synthèse numérique, en temps réel de phénomènes oscillant de toutes sortes. D'une façon générale, l'invention a donc pour objet la simulation d'une interaction non linéaire entre une source excitatrice et une onde dans un résonateur, au moyen d'outils de calcul de signaux numériques, soit à partir d'équations dont la solution correspond à la manifestation physique d'un phénomène à simuler, soit à partir d'une mesure de cette manifestation physique. Conformément à l'invention, le phénomène à simuler se traduisant, à chaque instant et en un point donné du résonateur, par une relation d'impédance entre deux variables représentatives de l'effet et de la cause dudit phénomène, on transcrit directement l'équation de l'impédance sous forme d'un modèle faisant appel à des guides d'onde pour réaliser une interaction non linéaire entre les deux variables de la relation d'impédance. A cet effet, le modèle comporte, d'une part, au moins une partie linéaire représentant directement l'impédance dite d'entrée du résonateur, c'est-à-dire l'impédance au point où se manifeste l'interaction non linéaire et, d'autre part, une partie non linéaire modélisant le rôle de la source excitatrice du phénomène à simuler. En particulier, pour la synthèse numérique, en temps réel, d'un phénomène oscillant, l'invention permet, à partir d'un système d'équations entre au moins deux variables représentatives du comportement du résonateur, d'établir une expression de l'impédance d'entrée du résonateur en terme de guide d'ondes linéaires, sans décomposition en ondes allerretour, de façon à réaliser au moins une partie linéaire du modèle qui peut être couplée à une boucle non linéaire faisant intervenir l'évolution de la nonlinéarité telle qu'elle s'exprime entre les deux variables de la relation d'impédance d'entrée du résonateur. De façon particulièrement avantageuse, cette partie linéaire du modèle constitue la somme de deux guides d'ondes élémentaires réalisant une fonction de transfert entre les deux variables de la relation d'impédance. Selon une autre caractéristique particulièrement avantageuse, le modèle est piloté par au moins deux paramètres représentatifs de l'interaction physique non linéaire entre la source et le résonateur, au moyen d'une boucle reliant la sortie à l'entrée de la partie linéaire et comprenant une fonction non linéaire jouant le rôle de source excitatrice pour le résonateur. Comme on l'a indiqué, l'invention s'applique, en particulier, à la synthèse en temps réel des sons produits par un instrument à vent. Dans ce cas, les deux variables de la relation d'impédance sont la pression et le débit acoustiques à l'entrée du résonateur. Dans le cas d'un résonateur cylindrique ayant une extrémité ouverte, il est particulièrement avantageux de réaliser la partie linéaire du modèle de transcription numérique de l'équation de l'impédance sous forme d'une somme de deux guides d'ondes élémentaires ayant pour source excitatrice le débit à l'entrée du résonateur et réalisant la fonction de transfert : dans laquelle : - # est la pulsation de l'onde - Ze(co) est l'impédance d'entrée du résonateur, - Pe(co) et Ue(#) sont les transformées de Fourier des valeurs adimensionnées de la pression et du débit à l'entrée du résonateur, - k(co) est une fonction de la pulsation de l'onde qui dépend du phénomène à simuler, - L est la longueur du résonateur. Selon une autre caractéristique, chacun des deux guides d'ondes fait intervenir un filtre ayant pour fonction de transfert : - F(c)2 - - exp(-2ik(w)L) et représentant un trajet aller-retour d'une onde avec changement de signe à l'extrémité ouverte du résonateur, chaque guide d'onde correspondant à un terme de l'équation de l'impédance. Un tel modèle peut avantageusement être piloté par la longueur du résonateur et au moins deux paramètres représentatifs de l'interaction physique non linéaire entre la pression et le débit à l'entrée du résonateur, au moyen d'une boucle reliant la sortie à l'entrée de la partie linéaire et comprenant une fonction non linéaire jouant le rôle de source excitatrice pour le résonateur. En particulier, pour la synthèse en temps réel des sons à simuler, on réalise une formulation, dans le domaine temporel, de la réponse impulsionnelle du résonateur, en approximant les pertes représentées par le filtre au moyen d'un filtre numérique approché. L'invention couvre d'autres caractéristiques essentielles mentionnées dans les revendications et portant, en particulier, sur les équations utilisées par l'outil de calcul des signaux numériques et qui conduisent à des modèles à guides d'ondes dépendant du phénomène à simuler. En effet, selon une caractéristique essentielle de l'invention, la méthode proposée pour la simulation d'un phénomène simple comme la propagation d'une onde dans un résonateur cylindrique, peut être adaptée de multiples façons pour la simulation de phénomènes plus complexes et, en particulier, de divers types d'instruments. Dans la description qui va suivre, on va donc exposer en détail la méthode de simulation, les équations utilisées et le modèle à mettre en #uvre pour la synthèse du son d'un instrument à résonateur acoustique cylindrique à anche, du type clarinette, ainsi que certaines adaptations pour la simulation d'autres types d'instruments. La figure 1 donne deux diagrammes représentant respectivement, à gauche la fonction de transfert, en Hertz, d'un modèle d'anche à un mode et à droite, la réponse impulsionnelle en fonction des échantillons, avec une fréquence d'échantillonnage de 44 100 Hertz. La figure 2 est un schéma de calcul par combinaison de guides d'ondes, représentant l'impédance d'entrée d'un résonateur cylindrique. La figure 3 donne deux diagrammes représentant respectivement, pour un résonateur cylindrique, en haut l'impédance d'entrée en fonction de la fréquence indiquée en Hertz et, en bas, la réponse impulsionnelle en fonction du temps, en secondes. La figure 4 est un schéma de calcul d'un modèle de simulation d'un instrument à anche à résonateur cylindrique. La figure 5 donne deux diagrammes analogues à la figure 3, représentant respectivement, pour un modèle de résonateur calculé selon l'invention, en haut l'impédance d'entrée approchée et en bas la réponse impulsionnelle approchée. La figure 6 donne deux diagrammes analogues à la figure 1, représentant respectivement, pour un modèle d'anche calculé selon l'invention, à gauche la fonction de transfert et à droite la réponse impulsionnelle. La figure 7a montre les variations, en fonction du temps indiqué en secondes, de la pression acoustique interne au niveau de l'embouchure d'un résonateur cylindrique. Les figures 7b et 7c sont des agrandissements des transitoires d'attaque et d'extinction. La figure 8 donne deux diagrammes représentant respectivement, à gauche la fonction de transfert et à droite la réponse impulsionnelle, pour un modèle d'anche multimode calculé selon l'invention. La figure 9 donne deux diagrammes représentant le spectre de la pression acoustique externe, respectivement en haut pour une anche à un seul mode et en bas pour une anche à modes multiples. La figure 10 est un schéma de calcul représentant l'impédance d'un résonateur cylindrique avec impédance terminale. La figure 11est un schéma de calcul représentant l'impédance d'un résonateur conique. La figure 12 est un schéma de calcul représentant l'impédance d'un résonateur simplifié pour cuivres. La figure 13 est un schéma général de calcul représentant l'impédance d'une combinaison en parallèle de résonateurs cylindriques. Comme on le sait, les phénomènes physiques mis en jeu lors de la production du son de la clarinette, s'expriment d'une part, par une équation de propagation linéaire des ondes dans le tuyau avec perte, et d'autre part, par une équation non linéaire reliant le débit avec la pression et le déplacement de l'anche au niveau de l'embouchure de l'instrument. Un modèle de simulation du son comprend donc une partie linéaire du modèle correspondant au résonateur de l'instrument qui, dans le cas de la clarinette est constitué d'un tube cylindrique. Pour cette géométrie, en supposant que le rayon du tube est grand devant les épaisseurs des couches limites, la pression acoustique à l'intérieur du tube est régie par une équation de la forme : où a = Rc 2 3 1 + (cp -1)$:J., R étant le rayon du tube, c'est-à-dire Rcy2 acv 7mm dans le cas de la clarinette. Les valeurs des constantes physiques, en unités mKs, sont : c=340, 1,=4.10-8, lt=5.6.10-8, Cp = 1.4. Cv En recherchant les solutions du type exp (i(#t-k(#)x)), # étant la pulsation de l'onde, on peut écrire : En remplçant#1+x par la valeur approchée 1 + x/2 lorsque x est petit, l'expression approchée classique de k(co), que nous utiliserons par la suite, devient : On sait que, si l'on considère un tuyau de longueur infinie que l'on suppose excité en x=0 et t=0 par une impulsion Dirac #(#)#(t), en tout point x>0, la pression acoustique propagée à partir de cette source s'écrit sous la forme d'une somme continue de toutes les ondes susceptibles de se propager dans le tuyau : p(x,t) = Jexp(-ik(a>)x)exp(ia>t)d# qui apparaît comme la Transformée de Fourier inverse de la valeur exp (-ik(#)). La fonction de transfert d'un tuyau de longueur L, qui constitue le filtre du modèle guide d'onde représentant la propagation, la dissipation et la dispersion est donc : La dissipation représentée par le module de F(co) et la dispersion représentée par la phase de F(co) sont donc proportionnelles à ##, tandis que le retard de propagation est donné par L . La longueur du tuyau sera c donc le paramètre de contrôle de la hauteur et son rayon contrôle les pertes. On sait, d'autre part, que les transformées de Fourier des pressions et débits adimensionnés à l'entrée (Pe(), Ue()) et en sortie (Ps(co), US()), du résonateur, sont liés par le système d'équations : Pe(co) = cos(k()L)Ps(#)+ /sin(k((D)L)Us(cù) Ue(m) = /sin(k(M)L)Ps(G))+ cos(k(#)L)Us() De façon classique, afin de modéliser la pression acoustique interne, le rayonnement peut être négligé. L'extrémité ouverte de l'instrument est donc parfaitement réfléchissante, ce qui entraîne que Ps(#)=0. Ceci permet d'exprimer la relation entre pression et débit à l'entrée du résonateur : P,.(co) = i tan(k((D)L)Ue(0)) = Ze (ro)Ue(ro) (4) où Ze(#) = i tan(k(co)L) est l'impédance d'entrée normalisée du résonateur. Dans le cas d'un modèle classique d'anche ou de lèvres à un mode, le déplacement adimensionné x (t) de l'anche par rapport à son point d'équilibre, et la pression acoustique pe(t) qui le produit sont liés par l'équation : dans laquelle #r=2#fr correspond à la fréquence de résonance fr, par exemple 2500 Hz et qr est le facteur de qualité de l'anche, par exemple 0,2. En écrivant, l'équation (5) dans le domaine de Fourier, on obtient la fonction de transfert de l'anche : dont la réponse impulsionnelle est donnée par : Comme on le verra plus loin, il est important de noter que x(0)=0. Par ailleurs, dans le cas d'un instrument à anche de type clarinette ou à embouchure de type trompette, la pression acoustique pe(t) et le débit acoustique ue(t) (adimensionnés) à l'entrée du résonateur sont reliés de façon non-linéaire par l'équation : Dans le cas d'un instrument à anche, le paramètre est caractéristique de l'embouchure et prend en compte la position des lèvres et le rapport de section entre le bec et le résonateur. Des valeurs classiques peuvent être de 0,2 à 0,6. Le paramètre y est le rapport entre la pression à l'intérieur de la bouche d'un instrumentiste et la pression de plaquage statique de l'anche. Pour un tuyau sans perte, il va de - pour la mise en vibration à - pour la position d'anche battante. Les paramètres C, et y sont donc deux paramètres de jeu importants dans la mesure où ils représentent, respectivement, la manière dont l'instrumentiste pince l'anche et la pression du souffle dans l'instrument. En combinant l'équation de déplacement de l'anche, la relation d'impédance et la caractéristique non linéaire, il apparaît que la pression acoustique au niveau de l'embouchure est contrôlée par le système d'équations suivant : La figure 1 donne deux diagrammes indiquant respectivement, à gauche la fonction de transfert et à droite, la réponse impulsionnelle du modèle d'anche, pour une fréquence de résonance fr= 2 500Hz et un facteur de qualité qr=0,2. L'idée de l'invention est donc de trouver une formulation dans le domaine temporel de la relation d'impédance permettant de résoudre ce système de trois équations, en modélisant la relation d'impédance en termes de modèles de guides d'onde élémentaires. Pour modéliser l'impédance d'entrée du résonateur en terme de modèles de guides d'onde élémentaires, on écrit la transformée de Fourier de l'impédance Ze(#) sous la forme : Cette expression peut s'écrire sous la forme : La figure 2 représente un modèle de calcul par combinaison de guides d'ondes directement issu de cette dernière équation et dont la fonction de transfert est l'impédance d'entrée du résonateur. Il est constitué d'une somme de deux guides d'ondes élémentaires. L'élément supérieur correspond au premier terme de l'équation (12) tandis que l'élément inférieur correspond au deuxième. Le filtre dont la fonction de transfert est - F(co)2 = - exp(-2ik(co)L) représente un trajet aller-retour, avec changement de signe de la pression acoustique à l'extrémité ouverte. A titre d'exemple, la figure 3 donne deux diagrammes représentant respectivement, pour un tuyau de longueur L = 0.5m et de rayon R=7mm, en haut la variation de l'impédance d'entrée du résonateur en fonction de la fréquence et, en bas, la réponse impulsionnelle du modèle à guides d'onde correspondant, calculée par transformée de Fourier inverse de l'impédance. Le système des trois équations physiques couplées (9), (10), (11) permet d'introduire la non-linéarité sous forme d'une boucle reliant la sortie pe du résonateur à l'entrée ue. La figure 4 donne un schéma de calcul équivalent permettant, pour la simulation d'un instrument à anche ou embouchure, de coupler non linéairement le déplacement de l'anche et la pression acoustique avec le débit acoustique à l'entrée du résonateur, en calculant, à chaque instant échantillonné, la pression acoustique interne au niveau de l'embouchure. Le modèle est entièrement piloté par la longueur L du résonateur et au moins deux paramètres et y représentatifs de l'interaction physique non linéaire entre la source et le résonateur, au moyen d'une boucle reliant l'entrée à la sortie de la partie linéaire et comprenant une fonction nonlinéaire jouant le rôle de source excitatrice pour le résonateur. Comme le montre la Figure 4, la partie linéaire reprend le schéma de la Figure 2 et la fonction non-linéaire f est contrôlée par les deux paramètres et y permettant de simuler le jeu d'un instrumentiste, et a comme paramètres d'entrée, dans le cas d'une clarinette, la pression à l'embouchure et le déplacement x (t) de l'anche par rapport à son point d'équilibre, calculé en fonction de la pression à l'embouchure, par un modèle d'anche (m). Aucun signal d'entrée n'est nécessaire, dans la mesure où y est directement proportionnel à la pression dans la bouche de l'instrumentiste. C'est donc la non-linéarité elle-même, et son évolution imposée par l'instrumentiste, qui joue le rôle de source excitatrice, conformément au modèle physique. Pour la synthèse, en temps réel, des sons à simuler, le modèle nécessite un échantillonnage numérique et, pour cela, on réalise une formulation, dans le domaine temporel, de la réponse impulsionnelle du résonateur, correspondant à la transformée de Fourier inverse de l'impédance. Cette formulation dans le domaine temporel permet de calculer la pression pe(t) à l'embouchure en fonction du débit ue(t) mais il est nécessaire, pour cela, d'approximer les pertes représentées par le filtre F(co) au moyen d'un filtre numérique approché. On va donc utiliser une approximation de la fonction de transfert du filtre, - F(co)2 = - exp(-2ik(co)L) , dont les coefficients sont déterminés à partir de variables physiques telles que la longueur du résonateur et son rayon, de façon à pouvoir apporter les modifications nécessaires en fonction de la géométrie du résonateur. A cet effet, on exprime analytiquement les coefficients du filtre numérique comme des fonctions de paramètres physiques. En pratique, il est particulièrement avantageux d'utiliser un filtre à un pô;e en exprimant la fonction de transfert approchée sous la forme : dans laquelle tu = # , fe étant la fréquence d'échantillonnage, et @e D = feL est le retard pur correspondant à un trajet aller ou retour des ondes c dans le résonateur. Les paramètres bo et ai sont exprimés en fonction des paramètres physiques de telle sorte que IF(w)2 = F(rs) - pour deux valeurs données de #. La première valeur #1 retenue est celle de la fréquence fondamentale de jeu. Ceci permet d'assurer un temps de décroissance de la fréquence fondamentale de la réponse impulsionnelle du modèle à guides d'ondes utilisant le filtre approché, identique à celui du modèle à guides d'ondes utilisant le filtre exact. La deuxième valeur #2 retenue est celle d'un harmonique choisi de façon à obtenir une décroissance globale identique des réponses impulsionnelles des guides d'onde, respectivement, exact et approché. Le choix de cette deuxième valeur #2 est donc plus libre. Elle correspond, par exemple, au second pic de résonance dans le cas de la clarinette mais, dans certains cas, comme on le verra plus loin dans le cas de la trompette, il peut être préférable de choisir un harmonique de rang supérieur. Le système d'équations à résoudre est donné par : En négligeant la dispersion introduite par la partie non-linéaire de la -H phase de F(c), les fréquences des harmoniques sont : tuk = Lfe ou 2k-1est le rang de l'harmonique. En notant : Ci = Cos(Ti3), C2 = cos('UJ2),F1 = F(o), )2 12@ F2 = IF(wZ )2 2, A = Fc,, A2 = F2c2, les coefficients ai et bo sont donnés par : A partir de l'expression de l'impédance d'entrée du résonateur : et en notant z=exp(im), on tire directement : d'où l'on tire l'équation aux différences : Pe (n) = ue (n) - a1ue (n -1) - boue (n - 2D) + a1pe (n -1) - boPe (n - 2D) (16) La figure 5 donne deux diagrammes montrant respectivement, en haut la variation de l'impédance d'entrée du résonateur approché en fonction de la fréquence et, en bas, la réponse impulsionnelle du modèle à guides d'onde correspondant, calculée à partir de l'équation aux différences pour un tuyau cylindrique de longueur L=0.5m et de rayon R=7mm. On constate que le modèle ainsi établi donne une réponse impulsionnelle très proche de celle du résonateur, représentée sur la figure 1. De même que pour le filtre représentant les pertes, la relation entre la pression acoustique et le déplacement de l'anche doit être discrétisée dans le domaine temporel. Or, la réponse impulsionnelle de l'anche est une sinusoïde amortie de façon exponentielle qui satisfait à la condition x(0)=0, comme on l'a indiqué plus haut. Il est donc possible de construire un filtre numérique pour lequel le déplacement de l'anche à l'instant tn = ../fe est fe fonction de la pression acoustique à l'instant tn-1= n-1/fe et non pas tn. Ceci permet de respecter la propriété x(0)=0 du système continu lorsque l'anche est soumise à une excitation Dirac. Pour satisfaire cette condition, au lieu d'utiliser la transformation bilinéaire pour approximer les termes i# et -#2, on utilise, selon l'invention, les expressions i# # fe/2(z - z-1) et -#2#fe2(z - 2 + z-1), qui correspondent à une différentiation numérique centrée exacte au second ordre. Avec ces approximations, la fonction numérique de transfert de l'anche est donnée par : d'où l'on tire l'équation aux différences : x(n) =0xpe(n)+b1aPe(n-1)+a1ax(n-1)+a2aX(n-2) (17) dans laquelle les coefficients b1a, a1a et a2a sont définis par : La figure 6 donne, pour un tel modèle d'anche approché, deux diagrammes représentant, à gauche la fonction de transfert et à droite la réponse impulsionnelle, avec une fréquence d'échantillonnage fe=44 100 Hertz, les valeurs des paramètres étant fr=2500Hz et qr=0,2. On voit que les diagrammes obtenus sont très proches de ceux de la figure 1. Ceci étant établi, on va maintenant exposer une méthode de résolution explicite, selon l'invention, pour le couplage des équations aux différences avec la caractéristique non linéaire. Les formulations échantillonnées des réponses impulsionnelles du déplacement de l'anche et de l'impédance permettent, en effet, d'écrire l'équivalent échantillonné du système d'équation (9,10, 11) indiqué plus haut, sous la forme : x(n)=b1ape(n-1)+a1ax(n-1)+a2ax(n-2) (18) pe(n)=ue(n)-a, ue (p- )-boue(n-2)+a, Pe(n- )-bope(n-2) (19) Ce système d'équations est implicite, du fait que le calcul de pe(n) par l'équation d'impédance nécessite la connaissance de ue(n) et que cette valeur est elle-même obtenue à partir de l'équation non linéaire et nécessite la connaissance de pe(n). Cependant, comme on l'a indiqué plus haut, le calcul de x (n) nécessite pas la connaissance de pe(n) mais celle de pe(n-1) qui est connu à l'instant n. Ceci permet, selon l'invention, de résoudre simplement et exactement le système couplé. Pour cela, les termes des équations 19 et 20 ci-dessus, qui ne dépendent pas de l'échantillon de temps n, peuvent, en effet, être regroupés dans les expressions : V=-ai ue(n-1 )-boUe(n-2 D)+ai Pe(n-1 )-boPe(n-2D) W = 1/2 (1 - sign(y - x(n) - 1)) x Ç(1 - y + x(n)) ce qui permet d'écrire : pe (n) =ue(n)+V Pour généraliser la méthode, il est intéressant d'associer ue (n) à un coefficient bco = 1 dans le cas d'un tuyau cylindrique. Les deux équations 19 et 20 ci-dessus peuvent alors être écrites sous la forme : Pe(n)=bocoue(n)+V Etant donné que le terme! (1 - sign(y - x(n) - 1)) annule W lorsque (1-y+x(n)) est négatif, W reste toujours positif. Si l'on considère successivement les deux cas : y-pe(n)>0 et y-pe(n) De la sorte, le calcul de la pression acoustique et du débit au niveau de l'embouchure à un instant échantillonné n, peut s'effectuer en utilisant de façon séquentielle les équations suivantes : x(n)=biaPe(n-1)+a1ax(n-1)+a2aX(n-2) (21) V=-aiue(n-1 )-boue(n-2D)+aipe(n-1 )-bope(n-2D) (22) W = 1/2 (1 - sign(y - x(n) - 1))Ç(1 - y + x(n)) (23) Pe(n)=bcoue(n)+V (25) Ainsi, l'invention permet de résoudre dans le domaine temporel le système d'équations régissant la modélisation physique de l'instrument, à partir d'une formulation échantillonnée équivalente de la réponse impulsionnelle du déplacement de l'anche, de la relation d'impédance et de la caractéristique non linéaire, qui se traduit par le système d'équations 18, 19, 20, dans lequel : - l'équation (18) est une transcription numérique du modèle (m) de la figure 4, - l'équation (19) est une transcription numérique du modèle d'impédance de la figure 2, - l'équation (20) est une transcription numérique de la caractéristique non linéaire reliant le déplacement de l'anche et la pression acoustique avec le débit acoustique. En regroupant les termes qui ne dépendent pas de l'échantillon temporel n, le procédé selon l'invention permet, en effet, de déterminer le débit et la pression à l'entrée du résonateur par un calcul séquentiel des équations 21 à 25, et de résoudre, dans le domaine temporel, le système d'équations 9,10, 11 régissant la modélisation physique d'un instrument à anche du type clarinette, afin de synthétiser les sons produits par un tel instrument. A titre d'exemple, la figure 7 montre la variation de la pression acoustique interne au niveau de l'embouchure, calculée par un tel modèle non linéaire faisant appel à des guides d'ondes, pour un tuyau de longueur L=0,5m et de rayon R=7mm, les valeurs des paramètres étant y=0,4, =0,4, fr=2205 Hz, qr=0,3. On observe trois phases : le transitoire d'attaque correspondant à une augmentation brutale de y et , le régime permanent pendant lequel y et diminuent progressivement, de façon linéaire, jusqu'au seuil d'oscillation, et le transitoire d'extinction. En pratique, la mise en #uvre numérique d'un tel modèle à guides d'ondes non linéaire, peut se faire en utilisant des éléments se trouvant dans le commerce. Par exemple, on peut effectuer une mise en #uvre numérique en langage C sous forme d'un objet externe clarinette pour l'environnement connu sous la marque Max-MSP, piloté à partir de commandes MIDI fournies par un contrôleur Yamaha WX5#. Ce contrôleur mesure la pression des lèvres sur l'anche, qui contrôle le paramètre , et la pression du souffle, qui contrôle le paramètre y. Ces informations reçues au format MIDI (donc entre 0 et 127) sont renormalisées pour correspondre à l'échelle des paramètres physiques. L'accord du guide d'onde est réalisé à partir de l'information MIDI pitch contrôlée à partir du doigté qui détermine la longueur L du tuyau. Cependant, comme dans un instrument réel, la hauteur change en fonction des paramètres physiques tels que y, , #r et qr. Or, en réalité, un instrument de musique n'est pas parfaitement accordé pour tous les doigtés. L'utilisation d'un filtre passe-tout pour la mise en #uvre de la partie fractionnaire du retard D n'est donc pas indispensable. En pratique, il apparaît que les sensations de jeu d'un tel instrument virtuel sont tout à fait comparables à celles d'un instrument réel. Cependant, le procédé qui vient d'être décrit pour la simulation du son d'un instrument à anche du type clarinette, peut encore être perfectionné. On sait, en effet, que la pression acoustique au niveau de l'embouchure n'est pas la variable représentative du son perçu. Il est donc intéressant de calculer la pression externe qui, pour un tuyau cylindrique, peut s'exprimer comme la dérivée en fonction du temps, du débit sortant : Pext(t)= dus (t) En négligeant encore le rayonnement, ce qui entraîne ps(t)=O, dt il vient : Pe(ro) = i sin(k(co)L Us(o Ue(co) = cos(k(#)L)Us(#) d'où l'on tire : Us(co) =exp( -ik(ro)L)(P e(ro )+Ue(ro Du point de vue perception, le terme exp(-ik(co)L) est négligeable. L'expression ci-dessus peut donc être simplifiée et devient : Pext(t)=-(Pe(t) + ue(t (26) Ainsi, d'un point de vue numérique, le calcul, à chaque instant échantillonné (n), de la pression externe pext(n), se réduit à une simple différence entre la somme de la pression interne et du débit, entre l'instant (n) et l'instant (n-1). D'autre part, des mesures physiques ont montré que les vibrations d'une anche étaient plus complexes qu'une simple sinusoïde amortie. Il est donc apparu que le modèle simple décrit plus haut pouvait être perfectionné de façon à améliorer la qualité du son produit, tel qu'il est perçu. Pour cela, on considère un modèle très simplifié d'anche sous forme d'une corde encastrée libre. Ainsi, le modèle d'impédance d'un tuyau cylindrique décrit plus haut, qui génère principalement des harmoniques impaires, sert de base pour la réalisation d'un modèle d'anche à mode multiple satisfaisant la même condition x(0)=0, ce qui permet de conserver le schéma de résolution numérique des équations 21 à 25 indiquées plus haut. Etant donné que la valeur de l'impédance est réelle pour tous les pics de résonance, sa réponse impulsionnelle est une somme de fonctions cosinus. D'autre part, on a vu que la réponse impulsionnelle d'un modèle d'anche à mode unique est une fonction sinusoïdale. On va donc utiliser, dans la définition du modèle, la fonction de transfert entre une sinusoïde amortie et une cosinusoïde amortie. En conservant les mêmes notations, la transformée de Fourier de x(t)= 4 w r expC- 2 rorqrt) sin # 4 - qr2wtJ q2 2 est donnée par : X(co) = 2 2 ro2 r .. De même, la transformée de Fourier de y(t) = . 200, z expC- 2 COrqrt) cos # 4 - qCOrt] est donnée par : V4-qr2 12 fonction transfert entre y(t) x(t) est alors : X(ro) 4 - q2 ru La fonction de transfert entre y(t) et x(t) est alors : # # = # # - . Y(co) rorqr + 2iw On peut ainsi écrire le modèle de transmittance sous la forme : Pour déterminer les trois inconnues #a, C, B, on s'impose trois conditions. La première est de conserver la fréquence du premier pic, qui est le pic maximal. On choisit pour cela, coa= #r en supposant que le décalage de fréquence résultant du facteur de qualité qr est négligeable. La seconde condition, satisfaite par le modèle d'anche à mode unique, consiste à imposer une valeur unitaire du module de transmittance pour #=0 de façon à maintenir X(#) # Pe(#) aux basses fréquences. La troisième condition est une valeur imposée de 1 pour le module qr de la transmittance à la fréquence #r, afin de conserver la hauteur du pic du modèle d'anche à mode unique. Grâce à ces trois conditions, la fonction de transfert ainsi réalisée reproduit les caractéristiques principales du modèle d'anche à mode unique. Les deux premières conditions conduisent au système d'équations : en posant : A, = 2 , A2= 2 corqr et A3=A2Aqr, les coefficients q? + B et C résolvant le système sont donnés par : De même que pour le cas de l'anche à mode unique, on va construire un modèle numérique dans lequel le déplacement de l'anche à l'instant échantillonné tn=n/f@est fonction de la pression acoustique, non pas à l'instant Te tn mais à l'instant tn-1= n - 1/f@. Ceci est possible du fait que la réponse Te impulsionnelle du modèle est la somme de fonctions sinusoïdales amorties. Pour remplir cette condition, on utilise une approximation de i# sous la forme i# #fe (z - 1) = fe(exp(i#) - 1). De plus, pour ajouter au modèle un paramètre de contrôle supplémentaire, on remplace le coefficient B par un filtre tel que B ba/aaexp(-i#). Ainsi, on peut ajuster l'amortissement des harmoniques 1 - aa exp(-i#) en fonction de l'amortissement de la fondamentale. Pour conserver la caractéristique X(0)=1, les paramètres ba et aa sont liés par l'équation B= ba 1 -a. Le terme exp (-i##/#a) est remplacé par son équivalent échantillonné : z-Da=exp(-i#Da) avec le retard Da défini par: Da=E #fe dans lequel E indique la partie entière. Si l'on note : P = - ffirqr ' la fonction de transfert numérique s'écrit : ce qui conduit à l'équation aux différences : x(n)=baPe(n-)+ba2Pe(n-2)+baolPeO-a-1) +a,lx(n-1) +aa2x(n-2)+aaDX(n-Da)+aaDix(n-Da-1) (29) dans laquelle les coefficients aa1, aa2, aaD2, aaD1 sont définis par : et les coefficients ba1, ba2, baD1 par : On remarque que l'équation 29 ainsi établie permet de déterminer le déplacement adimensionné x (n) de l'anche à l'instant échantillonné (n), à partir des instants précédents. Il apparaît donc que le schéma de calcul numérique du son dans le cas d'un modèle d'anche multimode peut être le même que celui de l'anche à mode unique en remplaçant l'équation 21 par l'équation 29 qui est une autre transcription numérique du modèle (m) de la figure 4. A titre d'exemple, la figure 8, qui est analogue aux figures 1 et 6, donne, pour un modèle approché d'anche multimode calculé selon l'invention, deux diagrammes représentant respectivement, en trait plein, à gauche la fonction de transfert et à droite la réponse impulsionnelle avec une fréquence d'échantillonnage fe=44100 Hertz, les paramètres ayant des valeurs fr=1837,5 Hz, qr=0,2, aa=0. Sur les mêmes diagrammes, on a superposé, en traits interrompus, la fonction de transfert et la réponse impulsionnelle du modèle de anche unique telles que représentées sur la figure 6. Selon un autre développement de l'invention, on peut aussi améliorer le modèle de son de clarinette de façon à le rendre plus naturel en lui incorporant un certain bruit, le système étant ainsi plus réaliste. Etant donné que le bruit est créé par une turbulence au niveau de l'anche avant le début du tuyau, le bruit est ajouté à x(t). Il apparaît d'autre part que, en pratique, le niveau de bruit dépend de la pression du souffle alors que sa couleur dépend de la pression des lèvres sur l'anche. En effet, d'un point de vue physique, plus l'anche est pressée, plus faible est l'ouverture entre l'anche et le tuyau et plus importante est la turbulence. On utilisera donc un simple modèle de bruit dont le niveau est piloté par y et la brillance pilotée par #. Pour la mise en #uvre numérique d'un tel modèle de bruit, on utilisera un filtrage passe-bas d'un bruit blanc. La fonction de transfert de ce filtre est donnée par Tb(z) = bb(1 - ab)/1 - abz-1, le coefficient bb étant piloté par [gamma], et le coefficient ab piloté par . Les lois de variation de bb et aa peuvent être déterminées de telle sorte que le son simulé par le modèle soit aussi réaliste que possible. Les deux diagrammes de la figure 9 montrent à titre d'exemple, la variation du module du spectre de la pression acoustique externe correspondant au son produit par le modèle, respectivement sur le diagramme du haut pour une anche à mode unique et sur le diagramme du bas pour une anche à mode multiple avec bruit additionnel, les paramètres de simulation étant les suivants : fr=2205 Hz, qr=0.25, aa=0, y=0. 44, =0.4, L=0,48 m, R=7.10-3 m Comme on l'a indiqué, le procédé selon l'invention, tel qu'il vient d'être décrit en détail, a pour objet la simulation de sons produits par un instrument de musique à anche et à résonateur cylindrique, du type clarinette. Mais l'invention ne se limite pas à une telle application et peut, au contraire, faire l'objet de nombreux développements. En effet, à partir du modèle physique non linéaire, faisant appel à des guides d'onde et schématisé sur la figure 4, ainsi que de sa transcription numérique selon la séquence des équations 21 à 25, il va être possible, de simuler le fonctionnement d'un résonateur ayant une géométrie quelconque, en modifiant le modèle d'impédance pour résonateur cylindrique et les équations aux différences associées, pour les remplacer par des impédances et des équations aux différences plus complexes. Tout en conservant les propriétés du modèle qui vient d'être décrit, d'autres modèles d'impédance plus complexes peuvent, en effet, être bâtis en combinant certains éléments d'impédance en parallèle ou en série et en proposant des approximations numériques permettant un usage explicite des variables physiques et un contrôle plus flexible de l'instrument numérique. A titre d'exemple, on va maintenant décrire certains développements du modèle de base pour résonateur cylindrique, en se référant aux figures 10 à 13 qui représentent des schémas de calcul équivalents faisant appel à des guides d'ondes et correspondant à des résonateurs ayant des géométries diverses. D'une façon générale, sur ces schémas, l'opérateur C(#) représente l'impédance d'entrée et C-1(#) l'admittance d'entrée d'un résonateur cylindrique, le modèle numérique correspondant à C-1(#) étant obtenu en changeant simplement le signe du coefficient bo. Un premier perfectionnement du modèle de base qui vient d'être décrit en se référant aux figures 2 et 4, va permettre, en faisant appel, de façon analogue, a des guides d'onde, de réaliser un modèle physique pour résonateur cylindrique avec impédance terminale. Un tel élément permettra, par exemple, de relier entre elles des portions de résonateurs cylindriques ayant des longueurs et des sections différentes, de façon à simuler l'impédance d'entrée d'un conduit de section variable, ou bien de prendre en compte l'impédance de rayonnement. Pour cela, on considère le formalisme de la ligne de transmission liant la pression et le débit acoustiques à l'entrée du résonateur (Pe(#), Ue(#)) et à son extrémité ouverte (Ps(#), Us(#)). Si l'on note Zc = #c/#R2 l'impédance caractéristique on a : Pe(co) = cos(k(o)L)Ps(G)) + iZcsin(k(ro)L)Us(ro) Ue(ro) = # sin(k(m)L)Ps(m) + cos(k()L)US() En notant Vs = Ps (ro) l'impédance de sortie et R(co) = # - , Us (ro) Zc + Zs (ro) on peut écrire l'impédance d'entrée de deux façons différentes : L'équation 31 montre donc que l'impédance d'un résonateur cylindrique avec impédance terminale peut être obtenue à partir de l'impédance d'un résonateur cylindrique sans impédance terminale, en remplaçant : exp(-2ik(ro)L) par R(w)exp(-2ik()L). La Figure 10 donne un schéma équivalent de calcul faisant appel à des guides d'onde, pour la mise en #uvre de l'équation 30, permettant de calculer l'impédance d'un résonateur cylindrique avec impédance terminale. Un tel modèle permet de générer en cascade l'impédance d'entrée d'un conduit ayant une géométrie quelconque et peut même s'appliquer à la simulation du conduit vocal. Mais l'invention peut faire l'objet de bien d'autres développements et permet en particulier, à partir du modèle de base physique pour résonateur cylindrique schématisé sur la figure 4, de construire des modèles spécifiques pour la simulation de divers instruments de musique. C'est ainsi que dans un premier développement du modèle de base, on va construire un modèle pour résonateur conique, utilisable, par exemple, pour la simulation d'un saxophone. Si l'on appelle L la longueur du tuyau, R son rayon d'entrée, # son ouverture, la distance xe entre le sommet du cône et l'entrée est : L'impédance d'entrée, relative à l'impédance caractéristique Z. = #c est alors donnée par l'expression : #R2 qui peut s'écrire sous la forme : La Figure 11 donne un schéma équivalent de calcul faisant appel à des guides d'onde, dans lequel l'élément noté D est l'opérateur différentiel D=ico et l'élément noté C-1(#) correspond au schéma de la figure 2 en remplaçant -exp (-2ik(co)L) par +exp (-2ik(co)L). Les équations aux différences correspondantes peuvent être établies en suivant un processus analogue à celui qui a été décrit plus haut. Ainsi, en utilisant la transformée bilinéaire pour approximer i#, la fonction de transfert numérique d'un résonateur conique est donnée par : En notant : Gp = 1+ 2fe 1 Xe et Gm = 1- 2fe 1 Xe , la fonction de transfert C C se réduit à : d'où l'on tire l'équation aux différences : pe(n) = bcoue(n) +bciue(n-1)+bc2Ue(n-2)+bcDUe(n-2D)+bCDiue(n-2D-1) + ac1pe(n-1)+ac2Pe(n-2)+acDPe(n-2D)+acDiPe(n-2D-1) (33) dans laquelle les coefficients bco, bac1, bc2, bco et bcD1 sont définis par: et les coefficients ac1, ac2, aco et acD1 sont définis par : De façon analogue, l'invention peut s'appliquer au cas de résonateurs courts qui apparaissent, par exemple, dans l'embouchure d'un cuivre ou d'un trou de registre. Pour cela, on admet que le rayon du résonateur court est suffisamment grand pour conserver le modèle de perte utilisé jusqu'à présent. Cette approximation d'un résonateur court se fera en approximant l'impédance Zi(#)= i tan (k(#)I) pour des valeurs faibles de k(co)l. On arrive ainsi à l'expression l,(ro)= i tan (k((o)i) = G(m) + iroH(ro) (34) Dans un autre développement du modèle de base pour résonateur cylindrique, l'invention permet également de simuler l'ensemble d'un résonateur de type cuivre. Dans ce cas, on modélisera l'embouchure par un résonateur de Helmholtz comprenant une cavité hémisphérique couplée avec un tuyau court cylindrique et un résonateur principal à tuyau conique. L'impédance d'entrée de l'ensemble du résonateur peut s'exprimer par : dans laquelle V= 4 #R3b est le volume de la cavité hémisphérique, L1 6 est la longueur du tuyau court, L2 est la longueur du tuyau conique, Z1 et Z2 sont les impédances caractéristiques des deux tuyaux qui dépendent de leurs rayons, k1(#) et k2(#) tiennent compte des pertes et du rayon R1 et R2 de chaque tuyau. Cette équation peut être écrite sous la forme : On peut ainsi établir le schéma équivalent de calcul représenté sur la figure 12, dans lequel l'opérateur noté C1(#) correspond au schéma de la figure 4 et l'opérateur noté S(co) correspond à l'impédance d'entrée du tuyau conique et au schéma de la figure 11. Des méthodes de résolution analogues à celles qui ont été décrites précédemment permettent d'exprimer le modèle équivalent numérique et les équations aux différences correspondantes. L'impédance du tuyau court peut être modélisée en utilisant l'approximation exprimée par l'équation (34) précédente. L'impédance du tuyau conique est représentée par le modèle correspondant à l'équation aux différences (33). L'admittance de la cavité ico V/#c2 est approximée par la transformation bilinéaire d z - 1 , où d=2fe. #c2 z + 1 En considérant l'association de la cavité hémisphérique et du tuyau court comme un résonateur de Helmholtz ayant pour fréquence de résonance (0, c #- , on peut utiliser cette fréquence coh pour approximer G(#) par G(#h) et H(co) par H(COh). Les deux fréquences utilisées pour le calcul des # # ~# u c(127rL+9x, +16L) calcul des coefficients ai and bo sont w, c(12L + 9ZXe + 1$L) qui 0)1 - 4L(4L + 3nxe + 4xe) QUOI correspond au premier pic d'impédance du tuyau conique, et #2 = #h. De plus, on utilisera Zn = #c/S@ = Z2 pour normaliser l'impédance d'entrée. S2 L'impédance numérique d'un résonateur de type cuivre est ainsi donnée par l'expression : qui peut être simplifiée sous la forme : d'où l'on tire l'équation aux différences : Les coefficients découlent d'un calcul direct à partir de l'équation (35). L'invention peut encore être appliquée à la modélisation d'un résonateur cylindrique avec trous de registre. On utilisera, à cet effet, des éléments faisant appel à des guides d'ondes et correspondant respectivement à un modèle physique de tuyau cylindrique avec impédance terminale représentant un tuyau de longueur L1 entre l'embouchure et le trou de registre, un modèle de tuyau court qui représente le trou de registre de longueur ht et le modèle de base pour tuyau cylindrique représentant un tuyau de longueur L2 entre le trou de registre et l'extrémité ouverte. L'impédance terminale de la première partie du tuyau peut s'écrire : En écrivant cette expression sous la forme : on peut établir un schéma général de calcul par guides d'onde représenté sur la Figure13, qui, dans le cas qui vient d'être décrit, est utilisé pour calculer l'impédance terminale par une combinaison en parallèle des impédances Z2C2(co) et ZtCt(co). Dans le cas limite d'un trou de registre fermé, avec deux parties du tuyau ayant le même rayon, on pose ZtCt(#)=# et Z2=Zc d'où il vient Zs(co)=iZc tan(k(ro)L2) et Zc(to)=iZc tan(k(ro)(L1+L2)). L'impédance totale d'entrée du tuyau peut alors s'exprimer par : En utilisant i#=fe(1-Z-1) dans l'impédance du tuyau court, comme précédemment, il en découle directement l'équation aux différences. Ainsi, le modèle de simulation du résonateur cylindrique du type clarinette, obtenu par transposition directe des équations simplifiées du comportement physique de l'instrument, peut être adapté à la simulation d'instruments à résonateur non cylindrique, tel que le saxophone, la trompette ou autres instruments à tuyau de cuivre. Mais l'invention ne se limite pas à un tel mode de réalisation et aux adaptations qui viennent d'être décrites car, sans s'écarter du cadre de protection défini par les revendications, elle peut s'appliquer à la simulation d'autres types d'instruments, par exemple à corde frottée comme le violon ou frappée comme le piano, les variables de la relation d'impédance étant alors la force exercée en un point sur la corde par l'archet ou le marteau et la vitesse de la corde en ce point. En effet, dans le cas d'une corde, la relation d'impédance mécanique entre la force exercée sur la corde en un point quelconque et la vitesse en résultant, permet d'exprimer l'admittance mécanique de la corde en ce point sous forme d'une combinaison en parallèle des admittances de chaque portion de corde selon le schéma de la figure 13, dans lequel les éléments notés C2 et Ct représentent l'admittance de chaque portion de corde et correspondent au schéma de la figure 2, le filtre F(#) du modèle de base pouvant être exprimé, soit à partir d'un modèle de propagation d'ondes de flexion dans une corde munie de raideur, soit à partir de mesures expérimentales. En effet, l'admittance mécanique de chaque portion joue le rôle d'une impédance acoustique. En outre, la méthode selon l'invention pourrait aussi s'appliquer à la simulation d'autres phénomènes oscillants, grâce à une adaptation de certaines équations aux différences et un choix de coefficients tenant compte des caractéristiques physiques des phénomènes à simuler. REVENDICATIONS 1. Procédé de simulation numérique d'une interaction non linéaire entre une source excitatrice et une onde dans un résonateur, au moyen d'outils de calcul de signaux numériques à partir d'équations dont la solution correspond à la manifestation physique d'un phénomène à simuler qui se traduit, à chaque instant et en chaque point du résonateur, par une relation d'impédance entre deux variables représentatives de l'effet et de la cause dudit phénomène à simuler, procédé dans lequel on transcrit directement l'équation de l'impédance sous forme d'un modèle faisant appel à des guides d'ondes pour réaliser une interaction non linéaire entre les deux variables de la relation d'impédance. 2. Procédé selon la revendication 1, de simulation d'un phénomène oscillant, caractérisé par le fait que le modèle comporte, d'une part au moins une partie linéaire représentant directement l'impédance d'entrée du résonateur et, d'autre part, une partie non linéaire modélisant le rôle de la source excitatrice du phénomène à simuler. 3. Procédé de simulation selon l'une des revendications 1 et 2, pour la synthèse numérique, en temps réel, d'un phénomène oscillant, caractérisé par le fait que, à partir d'un système d'équations entre au moins deux variables représentatives du comportement d'une onde dans le résonateur, on établit une expression de l'impédance d'entrée du résonateur en terme de guide d'onde linéaire, sans décomposition en onde aller-retour, de façon à réaliser au moins une partie linéaire du modèle. 4. Procédé selon la revendication 3, caractérisé par le fait que la partie linéaire du modèle est couplée à une partie non linéaire faisant intervenir l'évolution de la non-linéarité telle qu'elle s'exprime entre les deux variables de la relation d'impédance d'entrée du résonateur. 5. Procédé selon la revendication 4, caractérisé par le fait que la partie linéaire du modèle de simulation numérique de l'équation de l'impédance fait appel à deux guides d'ondes élémentaires réalisant une fonction de transfert entre les deux variables de la relation d'impédance. 6. Procédé selon la revendication 5, caractérisé par le fait que la partie linéaire à deux guides d'onde du modèle est couplée à une boucle reliant la sortie à l'entrée de ladite partie linéaire et comprenant une fonction faisant intervenir la non-linéarité telle qu'elle s'exprime physiquement. 7. Procédé selon la revendication 6, caractérisé par le fait que le modèle est piloté par au moins deux paramètres représentatifs de l'interaction physique non-linéaire entre la source et le résonateur, au moyen d'une boucle reliant la sortie à l'entrée de la partie linéaire et comprenant une fonction non-linéaire jouant le rôle de source excitatrice pour le résonateur. 8. Procédé selon l'une des revendications précédentes, pour la synthèse en temps réel de sons produits par un instrument à vent, caractérisé par le fait que les deux variables de la relation d'impédance sont la pression et le débit acoustiques à l'entrée du résonateur. 9. Procédé selon la revendication 8, caractérisé par le fait que, pour un résonateur cylindrique ayant une extrémité ouverte, la partie linéaire du modèle de transcription numérique de l'équation de l'impédance constitue la somme de deux guides d'ondes élémentaires ayant pour source excitatrice le débit à l'entrée du résonateur et réalise la fonction de transfert : dans laquelle : - # est la pulsation de l'onde - Ze(co) est l'impédance d'entrée du résonateur, - Pe(co) et Ue(#) sont les transformées de Fourier des valeurs adimensionnées de la pression et du débit à l'entrée du résonateur, - k(co) est une fonction de la pulsation de l'onde qui dépend du phénomène à simuler, - L est la longueur du résonateur. 10. Procédé selon la revendication 9, caractérisé par le fait que chacun des deux guides d'ondes fait intervenir un filtre ayant pour fonction de transfert : - F(#)2 = - exp(-2ik(#)L) et représentant un trajet aller-retour d'une onde avec changement de signe à l'extrémité ouverte du résonateur, chaque guide d'onde correspondant à un terme de l'équation de l'impédance. 11. Procédé selon la revendication 10, caractérisé par le fait que le modèle est piloté par la longueur du résonateur et au moins deux paramètres représentatifs de l'interaction physique non linéaire entre la pression et le débit à l'entrée du résonateur, au moyen d'une boucle reliant la sortie à l'entrée de la partie linéaire et comprenant une fonction non linéaire jouant le rôle de source excitatrice pour le résonateur. 12. Procédé selon la revendication 11, caractérisé par le fait que la fonction non linéaire a comme paramètres d'entrée la pression et le déplacement de l'organe de formation des vibrations et est contrôlée par au moins deux paramètres de simulation du jeu d'un instrumentiste. 13. Procédé selon la revendication 12, caractérisé par le fait que les paramètres de jeu pour le contrôle de la fonction non linéaire sont : - un paramètre # caractéristique de l'embouchure et de l'action de l'instrumentiste sur l'organe de formation des vibrations, - un paramètre y représentatif de la pression appliquée sur l'organe de formation des vibrations. 14. Procédé selon l'une des revendications 10 à 13, caractérisé par le fait que, pour la synthèse en temps réel des sons à simuler, on réalise une formulation dans le domaine temporel de la réponse impulsionnelle de l'impédance du résonateur, en approximant les pertes représentées par le filtre au moyen d'un filtre numérique approché. 15. Procédé selon la revendication 14, caractérisé par le fait que l'on utilise un filtre numérique à un pôle en exprimant la fonction de transfert approchée sous la forme : dans laquelle : - # = #/f, fe étant la fréquence d'échantillonnage, @e - D feL est le retard pur correspondant à un trajet aller ou retour c de l'onde dans le résonateur, les coefficients bo et a1 sont exprimés en fonction des paramètres physiques de telle sorte que F(co) ,2 = ##(#)#2pour une valeur #1 de la pulsation correspondant à la fréquence fondamentale de jeu et une autre valeur #2 correspondant à un harmonique. 16. Procédé selon la revendication 15, caractérisé par le fait que les coefficients bo et a1 sont obtenus par résolution du système d'équation avec IF(ú))212 = exp(-2ac 2L), lesdits coefficients étant donnés par les formules : dans lesquelles 17. Procédé selon la revendication 16, de simulation d'un instrument à anche à résonateur cylindrique, à partir d'une modélisation physique régie par le système d'équation : dans lequel (Or est la fréquence de résonance et qr est le facteur de qualité de l'anche, caractérisé par le fait que ledit système d'équations est résolu dans le domaine temporel à partir d'une formulation échantillonnée équivalente de la réponse impulsionnelle du déplacement de l'anche et de la relation d'impédance qui se traduit par le système d'équations : x(n)=b1ape(n-1)+aiax(n-1)+a2aX(n-2) (18) Pe(n)=ue(n)-a1u0 (n-1)-boue(n-2D)+aipe(n-1)bope(n-2D) (19) ue(n) = 2 (1- sign(y - x(n) - 1))sign(Y - pe (n))n)- y + x(n ly - p-, (n)l (20) lesdites équations étant utilisées de manière séquentielle en regroupant les termes ne dépendant pas de l'échantillon temporel n, de façon à calculer successivement : 18. Procédé selon la revendication 17, pour une simulation plus réaliste du son produit, caractérisé par le fait que, en négligeant le rayonnement, la pression externe s'exprime comme la dérivée en fonction du temps, du débit sortant, sous la forme et est calculée, à chaque instant échantillonné (n), par différence entre les sommes de la pression interne pe et du débit ue, respectivement à l'instant (n) et à l'instant (n-1). 19. Procédé selon la revendication 17 de simulation d'un instrument à anche multimode, caractérisé par le fait que le calcul de la pression acoustique et du débit au niveau de l'embouchure s'effectue par résolution séquentielle d'un système d'équations dans lequel le déplacement de l'anche à chaque instant (n) est de la forme : x(n)=baipe(n-1)+ba2pe(n-2)+baoipe(r'-Da-1)+aaiX(r)-1)+aa2x(n-2)+aaDX(n-Da)+aaDix(n-Da-1) les coefficients aa1, aa2, aaD2, aaD1 étant définis par : et les coefficients ba1, ba2, baD1 par : les équations suivantes étant les mêmes que pour une anche à mode unique. 20. Procédé selon l'une des revendications 17 à 19, caractérisé par le fait que l'on réalise un modèle pour résonateur cylindrique avec impédance terminale à partir du modèle de base correspondant à un résonateur cylindrique et constituant la somme de deux guides d'onde faisant intervenir chacun un filtre ayant pour fonction de transfert -F(#)2=-exp(-2ik(#)L), en remplaçant l'expression exp(-2ik()L) par l'expression R(c)exp(-2ik(w)L), dans laquelle R(co) = e 5 zc + Zs (w) Zc étant l'impédance caractéristique P 2 etZs l'impédance de sortie S nR u (Co) 21. Procédé selon l'une des revendications 17 à 20, caractérisé par le fait que, à partir du modèle d'impédance pour un instrument à résonateur cylindrique et des équations aux différences associées, on construit d'autres modèles d'impédance plus complexes pour la simulation de phénomènes oscillants produits par un résonateur de forme quelconque en combinant des éléments d'impédance en parallèle ou en série et en utilisant des approximations numériques pour un usage explicite des variables physiques intervenant dans la production desdits phénomènes oscillants et un contrôle plus flexible du résultat de la simulation. 22. Procédé selon l'une des revendications 17 à 21, caractérisé par le fait que, à partir du modèle de base pour résonateur cylindrique dans lequel la réponse impulsionnelle du déplacement de l'anche se traduit par un système d'équations aux différences donnant, à chaque instant (n), le déplacement x (n) la pression pe(n) et le débit ue(n) à l'entrée du résonateur, on construit un modèle pour résonateur conique dans lequel l'équation de la pression est de la forme : pe(n) = bcoue(n) +bciUe(n-1)+bc2Ue(n-2)+bCDUe(n-2D)+bCDiUe(n-2D-1) + acipe(n-1)+ac2pe(n-2)+acDpe(n-2D)+acDipe(n-2D-1) (33) dans laquelle les coefficients bco, bc1, bc2, bco et bcD1 sont définis par : et les coefficients aci, ac2, acD et acD1 sont définis par : 23. Procédé selon l'une des revendications 17 à 21, caractérisé par le fait que, à partir du modèle de base pour résonateur cylindrique, on construit un modèle pour un résonateur court ayant une longueur I, en approximant l'impédance selon l'expression : Z@(#)= i tan (k(co)l) =- G(co) + i#H(#) (34) dans laquelle G(co) = 1 - exp - ac 2 et (co)=1(1 - G(co . 1+exp -ac 2 24. Procédé selon l'ensemble des revendications 21,22, 23, pour la simulation d'un instrument du type cuivre, caractérisé par le fait que l'on modélise l'embouchure par un résonateur de Helmholtz comprenant une cavité hémisphérique couplée avec un tuyau court cylindrique et un résonateur principal à tuyau conique, l'impédance d'entrée de l'ensemble du résonateur pouvant s'exprimer par l'expression : dans laquelle V=4/6 #Rb3 est le volume de la cavité hémisphérique, L1 est la longueur du tuyau court, L2 est la longueur du tuyau conique, Z1 et Z2 sont les impédances caractéristiques des deux tuyaux qui dépendent de leurs rayons, k1(#) et k2(co) tiennent compte des pertes et du rayon R1 et R2 de chaque tuyau, et que, à partir du modèle de base pour résonateur cylindrique, de ses extensions au tuyau conique et au tuyau court, on construit un modèle pour résonateur de type cuivre en exprimant la pression à l'embouchure par l'équation aux différences : 25. Procédé selon la revendication 21, pour la simulation d'un phénomène oscillant dans lequel les deux variables physiques de la relation d'impédance sont la force exercée en un point d'un système mécanique tel qu'une corde engendrant des vibrations et la vitesse en ce point, caractérisé par le fait que l'on exprime l'admittance en ce point sous la forme d'une combinaison en parallèle des admittances de chaque portion de corde, de part et d'autre dudit point, chaque admittance mécanique étant obtenue à partir du modèle de base décrivant l'impédance acoustique d'un résonateur à tuyau cylindrique, en exprimant la vitesse au point considéré de la corde en fonction de la force exercée en ce point, le filtre F(co) du modèle de base pouvant être exprimé, soit à partir d'un modèle de propagation d'ondes de flexion dans une corde munie de raideur, soit à partir de mesures expérimentales.