Rechercher : dans
Par :

[MATLAB] Modélisation vibrations tambour

Dernière réponse le 31 oct 2008 à 10:19:15 snakehill, le 30 oct 2008 à 22:09:02 
 Signaler ce message aux modérateurs

Bonsoir,
je suis peut être en vacances, mais je m'acharne et bloque sans trouver de solution à un problème de calcul de valeurs propres de matrice.
Voici le topo c'est une modélisation de vibration sur la surface d'un tambour sous MATLAB.
On a créé une matrice 20x15 que l'on rempli de zeros. Jusque là tout va bien.
J'ai développé une fonction qui utiliser la méthode des puissances itérées afin de trouver les valeurs propres d'une matrice numériquement, et en utilisant une autre fonction , celle de jacobi , on trouve le vecteur propre correspondant.

Seul petit hic , l'équation différentielle du tambour est
laplacien(u) = lambda*u <=> d²(u)/dx² + d²(u)/dy² = lambda*u

qui peut s'écrire (grace a de grossières approximation) :
u(x-1,y)+u(x+1,y)+u(x,y-1)+u(x,y+1)-4*u(x,y)=lambda*u
<=>u(x-1,y)+u(x+1,y)+u(x,y-1)+u(x,y+1)-(4+lambda)*u(x,y)=0

Cette équation , pour chaque point sera mise en matrice le lambda étant la valeur propre de la matrice , mais de quelle matrice ?

C'est assez confus dans ma tête et je n'en vois vraiment pas le bout ...

Configuration: Windows XP
Firefox 3.0.3

Meilleures réponses pour « [MATLAB] Modélisation vibrations tambour » dans :

1

 Sacabouffe, le 31 oct 2008 à 10:19:15
  • +1

Salut
La matrice du laplacien discret.
Si t'écris u comme un vecteur, ton équation s'écrit A u = λ uA est la matrice du Δ discret.
Ciao Thought I heard a rumbling, calling to my name
Two hundred million guns are loaded, Satan cries "Take aim!"

Répondre à Sacabouffe