¿Por qué todas las soluciones de este sistema de ecuaciones diferenciales de péndulo son una combinación lineal de las dos soluciones dadas?

Actualmente estoy tratando de hacer un informe de laboratorio para un experimento de péndulos acoplados en el que encontramos el siguiente sistema lineal de ecuaciones diferenciales de segundo orden (que describe la posición en función del tiempo de las dos masas):

d 2 d t 2 ( X 1 X 2 ) = ( ω 0 2 k k k ω 0 2 k ) ( X 1 X 2 )

Y las notas de laboratorio establecen que todas las soluciones son combinaciones lineales de las siguientes dos soluciones

( X 1 X 2 ) = A s y metro ( 1 1 ) porque ( ω s y metro t + ϕ s y metro )

( X 1 X 2 ) = A a s y metro ( 1 1 ) porque ( ω a s y metro t + ϕ a s y metro )

dónde ω s y metro = ω 0 y ω a s y metro = ω 0 2 + 2 k

Puedo entender intuitivamente por qué las dos soluciones abarcan todas las soluciones, pero realmente no puedo entenderlo matemáticamente. Creo que mi confusión proviene del hecho de que abarcan un espacio de funciones en lugar de vectores finitos. ¿Puede alguien explicar completamente de una manera simple por qué las dos soluciones abarcan todas las soluciones?

¿Puede encontrar que este documento sobre modos normales que describe las coordenadas normales así como los modos normales le ayudará?
En términos generales, eso se debe a que puede ver la matriz correspondiente al sistema de ecuaciones como una aplicación lineal (al igual que en Álgebra lineal). Como tienes dos soluciones independientes lineales y la dimensión del espacio (de funciones que son soluciones) es dos, tienes una base X 1 , X 2 del espacio, y por lo tanto puede abarcar todas las soluciones de aquellos.
¿Cómo sabes exactamente que la dimensión del espacio de las funciones de solución es dos? Por ejemplo, ¿cómo sé que una función continua arbitraria, como una curva gaussiana o una exponencial decreciente, no puede ser una solución? (por supuesto que entiendo por qué estos dos ejemplos no pueden funcionar físicamente, pero quiero demostrar que ninguna otra solución puede funcionar).

Respuestas (1)

Recuerde del álgebra lineal que dado un operador A actuando sobre vectores v de R norte , la ecuacion

A v = 0
define el núcleo o espacio nulo del operador A . Este es el subespacio de todos los vectores que son asignados a 0 por dicho operador y su dimensión corresponde a norte Rango ( A ) . Este es el caso de los espacios de dimensión finita, pero se generaliza a los espacios vectoriales de funciones y se aplican casi las mismas cosas.

Así que aquí tienes un operador diferencial

A = ( d d t 2 + ω 0 2 + k 2 k k d d t 2 + ω 0 2 + k 2 )
que está actuando sobre algún espacio espacio de funciones V (No entraré en detalle en qué espacio es este pero, suponga funciones "agradables"). Como puede ver, el operador en realidad está actuando sobre el producto directo de dos de esos espacios, V 2 = V × V , es decir, actúa sobre dos copias de algún espacio de funciones, por lo que los vectores se ven como v = ( F 1 ( t ) , F 2 ( t ) ) y estamos resolviendo para
(1) A v = 0
como antes. Ahora, bajo la misma filosofía del caso habitual del álgebra lineal, lo más que puede pasar es que el operador envíe todo a cero, en cuyo caso el espacio nulo tendrá dimensión 2 norte dónde norte sería ahora la dimensión de V y podría ser infinito. Sin embargo, de la teoría de las ecuaciones diferenciales sabemos que las ODE de segundo orden deberán suministrarse con dos condiciones de contorno, y esto es exactamente lo que está relacionado con la libertad de "direcciones" que tiene en este espacio nulo. Entonces, si piensa en las soluciones como direcciones en un espacio vectorial, esto le dice que hay como máximo dos direcciones en las que puede "moverse" o cambiar las condiciones sin dejar de satisfacer (1) . Esto le pasaría a cada copia por separado si no hubiera términos fuera de la diagonal, pero como aquí no son cero, lo que hacen es acoplar las ecuaciones y así forzar una relación particular entre los componentes, dada por los autovalores. Al fin y al cabo, lo importante es darse cuenta de que todo se reduce a la dimensión de un espacio nulo de un operador en un espacio de funciones y esto está limitado por el orden de las ODE involucradas y el tamaño del sistema.