Procédé de simulation numérique d’un son d’un instrument de musique comprenant une anche, un canal d’anche et un résonateur, ce procédé étant mis en œuvre par un ordinateur et comprenant les étapes suivantes : - réception d’un signal représentatif d’une pression appliquée en entrée de l’instrument et destinée à se propager dans l’instrument, - réception d’un signal représentatif d’un doigté appliqué à l’instrument, - utilisation d’un modèle de l’anche sous forme d’un modèle masse-ressort à deux paramètres, pulsation de résonance ω r et amortissement q r , pour déterminer un déplacement de l’anche, - utilisation d’un modèle de canal d’anche pour déterminer un débit acoustique à partir du déplacement de l’anche, - utilisation d’un modèle du résonateur pour réaliser, à partir du débit acoustique et de la pression, une décomposition modale dans laquelle la pression à l’entrée du résonateur est décomposée en une somme de N m pressions modales p n , et - utilisation d’un modèle de son rayonné pour générer un son à partir des pressions modales et du débit acoustique. Procédé selon la revendication 1, caractérisé en ce que le son étant fonction d’une dérivée de la pression et d’une dérivée du débit acoustique, on affecte des coefficients multiplicateurs r n aux dérivées des pressions modales pour définir une enveloppe spectrale du son attendu et défini par le doigté. Procédé selon la revendication 1 ou 2, caractérisé en ce que le son étant fonction d’une dérivée de la pression et d’une dérivée du débit acoustique, on affecte un coefficient multiplicateur r u à la dérivée du débit acoustique pour définir une enveloppe spectrale du son attendu et défini par le doigté. Procédé selon la revendication 2 ou 3, caractérisé en ce que les coefficients multiplicateurs sont fonction du doigté, des paramètres de la pression en amont de l’anche p m et de l’ouverture de l’anche. Procédé selon les revendications 2 et 3, caractérisé en ce que lors d’une transition entre doigtés, le coefficient multiplicateur r u est mis à zéro, et les coefficients multiplicateurs r n sont mis à zéro sauf le premier. Procédé selon l’une quelconque des revendications précédentes, caractérisé en ce qu’il comprend une réception d’un signal représentatif d’une action sur l’ouverture au repos de l’anche ; et en ce que N m est compris entre 2 et 30. Procédé selon l’une quelconque des revendications précédentes, caractérisé en ce que la décomposition modale comprend une étape de calcul d’une composante strictement causale des pressions modales. Procédé selon la revendication 7, caractérisé en ce que la composante strictement causale des pressions modales comprend une composante retardée du débit acoustique à laquelle est ajoutée une correction. Procédé selon la revendication 8, caractérisé en ce que la correction comprend l’application d’un bruit de turbulence sous forme d’un bruit blanc ou d’un bruit coloré. Procédé selon la revendication 8 ou 9, caractérisé en ce que la correction comprend un sinus dont la fréquence est fonction du doigté, et dont la phase est incrémentée à chaque pas de temps en fonction du doigté. Procédé selon la revendication 10, caractérisé en ce que le sinus est multiplié par une partie positive de la pression en amont de l’anche, cette partie positive étant filtrée par un filtre passe-haut. Procédé selon l’une quelconque des revendications précédentes, caractérisé en ce qu’une méthode de l’invariance de la réponse impulsionnelle est appliquée sur les paramètres du modèle masse-ressort pour déterminer un filtre numérique liant la position de l’anche sans limitation x t à l’écart de pression entre l’amont et l’aval de l’anche. Procédé selon la revendication 12, caractérisé en ce que le filtre numérique comprend une correction de façon à ce que l’ouverture de l’anche ne conserve que la partie positive de la position de l’anche. Procédé selon la revendication 13, caractérisé en ce que la correction est définie par l’équation suivante : , étant l’ouverture d’anche corrigée, étant la position de l’anche sans limitation. Procédé, caractérisé en ce qu’on ajoute un paramètre de régularisation η dans une équation de l’ouverture d’anche selon la revendication 13 ou dans l’équation d’ouverture d’anche corrigée selon la revendication 14. Procédé selon la revendication 15, caractérisé en ce que la régularisation consiste à modifier la valeur absolue dans l’équation de l’ouverture d’anche corrigée telle que : Procédé selon la revendication 15 ou 16, caractérisé en ce que pour simuler la production d’un second registre, le paramètre de régularisation η est compris entre 0,015 et 0,05. Procédé selon la revendication 15 ou 16, caractérisé en ce que pour simuler la production d’un registre altissimo, le paramètre de régularisation η est compris entre 0,04 et 0,2. Procédé selon l’une quelconque des revendications précédentes, caractérisé en ce que pour simuler la production d’un second registre, la pulsation de résonance ω r et l’amortissement q r sont déterminés tels que : ω r est compris entre 0,9 ω note et 1,4 ω note , ω note étant la pulsation de la note attendue, définie par le doigté, q r est compris entre 0,02 et 0,4. Procédé selon l’une quelconque des revendications 1 à 18, caractérisé en ce que pour simuler la production d’un registre altissimo, la pulsation de résonance ω r et l’amortissement q r sont déterminés tels que : ω r est compris entre 0,8 ω note et 1,5 ω note , ω note étant la pulsation de la note attendue, définie par le doigté, q r est compris entre 0,0005 et 0,04. Procédé selon l’une quelconque des revendications précédentes, caractérisé en ce que les modèles de l’anche, du canal de l’anche, du résonateur et du son rayonné sont régis par des équations comprenant des coefficients sur une partie desquels un filtre passe-bas d’ordre 1 est appliqué pour rendre moins brutales des transitions entre doigtés ; des coefficients liés à la pression dans la décomposition modale ne subissent pas ce filtre passe-bas d’ordre 1. Produit programme d'ordinateur comprenant des instructions qui, lorsque le programme est exécuté par un ordinateur, conduisent celui-ci à mettre en œuvre le procédé selon l’une quelconque des revendications précédentes. Système de traitement de données comprenant un processeur adapté au procédé selon l’une quelconque des revendications 1 à 21.