El cálculo de los modos propios de una membrana semicircular se reduce al siguiente problema de valores propios
donde la región de interés es un semicírculo definido por y φ ∈ [ 0 , π ] .
Es apropiado trabajar en coordenadas cilíndricas, donde el laplaciano se escribe como
Las condiciones de contorno fijan el valor de en el límite del semicírculo, donde u = 0 .
Primero, hacemos una discretización de con u i j = u ( r i , φ j ) , donde r i = ( i + 1yφj=(j+1i,j=0...N-1yhr=1/N,hr=π/N. Esta es unamallacentrada.
Luego usamos una aproximación de diferencia finita para el laplaciano y obtenemos
o
Debido a que nuestra malla está centrada, tenemos que hacer el siguiente reemplazo en la ecuación anterior: . Este reemplazo también nos ayuda a deshacernos de la singularidad de coordenadas parai=0.
Las condiciones de contorno en y r = 0 , 1 puede ser todo manejado con el mismo truco , donde se parte en el límite
u i , j + 1 = - u i , j u i - 1 , j = - u i , j u i + 1 , j = - u i , j .
La matriz es una matriz real asimétrica y los valores propios y los vectores propios se pueden obtener con una rutina dgeev
de LAPACK.
Las soluciones analíticas pueden obtenerse fácilmente mediante el método de separación de variables.
Son
Aquí está la gráfica de la solución analítica para la primera función propia:
La siguiente gráfica muestra la comparación de resultados numéricos para tres discretizaciones diferentes, en la medida en que mis recursos computacionales me permiten ir.