La ecuación de onda en una aplicación no uniforme de la cadena es : $$ u_{tt} = c(x)^2 u_{xx} $$ $$ u(x,0) = f(x) = e^\frac{(x-\mu)^2}{2 \sigma^2} , \:\:u(0,t) = 0\:,\:\:u(L,t) = 0, \:\: u_{t}(x,0) = -cf'(x) $$ con $c(x) = c_{1}$$0 \le x \le \frac{L}{2} $, e $c(x) = c_{2}$$ \frac{L}{2} < x \le L$. Donde $L$ es la longitud de la cadena.
He intentado usar el método de diferencias finitas para animar la solución y el resultado es bastante agradable, cambiando el habitual esquema de la diferencia de $$ u^{m+1} = (2I - (c\triangle t)^2A)u^{m} - u^{m-1} $$ a $$ u^{m+1} = (2I - C(\triangle t^2)A)u^{m} - u^{m-1} $$ donde $u(x_{i},t_{m}) \approx u_{i}^{m} $, $x_{i} = i\triangle x$, $t_{m}= m \triangle t$.
$ i = 0,1,2,3,...,l$ donde $x_{l} = L$ y
$ m = 0,1,2,3,.....$
también : $ u^{m} = \left[ \begin{array}{c} u_{1}^{m} \\ u_{2}^{m} \\ . \\ . \\ u_{l-2}^{m} \\ u_{l-1}^{m} \end{array} \right] $ , $A = \frac{1}{(\triángulo x)^2} \left( \begin{array}{cccccc} 2 & -1 & 0 & 0 & 0 & 0 \\ -1 & 2 & -1 & 0 & 0 & 0 \\ 0 & . & . & . & 0 & 0 \\ 0 & 0 & . & . & . & 0 \\ 0 & 0 & 0 & -1 & 2 & -1 \\ 0 & 0 & 0 & 0 & -1 & 2 \\ \end{array} \right) $ , $ C = \left( \begin{array}{cccccccc} c_{1}^2 & 0 & 0 & 0 & 0 & 0 & 0 & 0\\ 0 & c_{1}^2 & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & . & 0 & 0 & 0 & 0 & 0\\ 0 & 0 & 0 & c_{1}^2 & 0 & 0 & 0 & 0\\ 0 & 0 & 0 & 0 & c_{2}^2 & 0 & 0 & 0\\ 0 & 0 & 0 & 0 & 0 & . & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & c_{2}^2 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 & c_{2}^2 \\ \end{array} \right) $
$A $ $C$ matriz cuadrada del mismo tamaño (que es $l-2 $).
El resultado se ve bien (http://www.youtube.com/watch?v=9lDRms5i0Hc). Cuando la onda pasa a través de la articulación, algunos son reflejadas y algunos de transmisión sexual, pero mi profesor dijo que la solución numérica está mal y debo añadir otra condición para $u$ en la articulación.
Yo quería demostrar mediante la comparación de mi solución numérica con la solución exacta para la no uniforme de la cadena de la ecuación de onda, alguien puede compartir su pensamiento para la solución exacta, y el esquema numérico que he hecho? Gracias