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 !
Partager