Simulation corde de piano sous MATLAB

Bonjour à tous,

J'essaye de simuler une corde de piano avec raideur et pertes par viscoelasticité et rayonnement.
J'ai un schéma de la forme u(i, n+1) = f(u(i, n-1),u(i-2,n), u(i-1,n), u(i,n), u(i+1,n), u(i+2,n))
Du coup je me suis fait 5 listes, une pour chaque terme (de dimension Nbr d'éléments +4).
Ensuite je fais une double boucle, la première est sur l'espace, la seconde sur le temps (celle sur l'espace est contenue dans celle sur le temps pour etre le plus clair possible). Au début de la boucle d'espace je "réattribue" les différentes listes au niveau suivant (je ne sais pas si c'est clair, grosso modo L(i-2) <- L(i-1) etc ...) et au début de la boucle de temps je fais de même avec les listes correspondant au pas de temps.
Cette méthode vous parait-elle correcte ou complètement fausse dès le départ ?
Ensuite je ne sais pas comment afficher l'évolution temporelle de ma corde, comment pourrais-je m'y prendre ?

Merci d'avance pour votre aide !

Salut!
La vibration d'une corde obéit à l'équation de d'Alembert. C'est une équation aux dérivées partielle du second ordre. Elle relie une dérivée seconde par rapport à x et la dérivée seconde par rapport à t. La meilleure stratégie consiste à appliquer la méthode des dérivées finies par rapport à x et la méthode de Rung-Kutta par rapport à t.
Jean-Marc Blanc