La respuesta de Nikhil Patel cubre las más convenientes. Esta respuesta es simplemente en aras de la integridad, ya que parece ser un depósito de algún tipo para resolver DE.
Resolveré el ODE de segundo orden reduciéndolo a un sistema de dos ODE de primer orden. Esto se puede hacer mediante la sustitución adecuada de la siguiente manera.
[matemáticas] y ” + 4y = 0 [/ matemáticas]
Poner [matemáticas] y_1 = y [/ matemáticas] y [matemáticas] y_2 = y ‘= y’ _ {1} [/ matemáticas]
- ¿Debo tomar álgebra lineal o ecuaciones diferenciales como un CS mayor?
- Cómo construir un lagrangiano a partir de una ecuación de campo dada (además de resolver la ecuación diferencial EL)
- Cómo evaluar [math] \ int _ {- \ infty} ^ \ infty e ^ {- ax ^ 2} \ cos (kx) dx [/ math]
- Cómo resolver el ODE [matemáticas] \ frac {d \ theta} {dt} = – \ sqrt {2K (\ cos (\ theta) – \ cos (\ theta_0))} [/ matemáticas] donde [matemáticas] K [ / math] es una constante positiva, y [math] \ cos (\ theta) \ geq \ cos (\ theta_0) [/ math]
- ¿Es útil conocer ecuaciones diferenciales parciales como científico de datos? Si es así, ¿cuáles son algunas de sus aplicaciones?
Por lo tanto, el sistema de dos EDO de primer orden se convierte en
[matemáticas] y ‘_ {1} = y_ {2} [/ matemáticas]
[matemáticas] y ‘_ {2} = -4y_ {1} [/ matemáticas]
que se puede escribir en forma de matriz como
[matemáticas] \ begin {bmatrix} y’_1 \\ y’_2 \ end {bmatrix} = [/ math] [math] \ left [\ begin {array} {lr} 0 & 1 \\ -4 & 0 \ end {array} \ right] [/ math] [math] \ begin {bmatrix} y_1 \\ y_2 \ end {bmatrix} [/ math]
[matemáticas] \ implica \ textbf {Y ‘} = A \ textbf {Y} [/ matemáticas]
donde [math] \ textbf {Y} = [/ math] [math] \ begin {bmatrix} y_1 \\ y_2 \ end {bmatrix} [/ math] y
[matemática] A = [/ matemática] [matemática] \ left [\ begin {array} {lr} 0 & 1 \\ -4 & 0 \ end {array} \ right] [/ math].
Sustituyendo [math] \ textbf {Y} = \ textbf {X} e ^ {\ lambda t} [/ math] que al diferenciar con respecto a ‘t’ se convierte en [math] \ textbf {Y ‘} = \ lambda \ textbf {X} e ^ {\ lambda t} = A \; [/ math] [math] \ textbf {X} e ^ {\ lambda t}, [/ math] al cancelar el término exponencial de ambos lados,
[matemáticas] A \ textbf {X} = \ lambda \ textbf {X} [/ matemáticas].
Este es el problema del valor propio que se resuelve para los valores propios y los vectores propios correspondientes.
[matemáticas] | A- \ lambda I | = 0 [/ matemáticas] da [matemáticas] {\ lambda} ^ 2 + 4 = 0 [/ matemáticas]
[matemáticas] \ lambda_1 = 2i; \; \; \; \ lambda_2 = -2i [/ math] son los valores propios.
Los vectores propios correspondientes son [math] \ textbf {X} _ {1} = \ begin {bmatrix} 1 \\ 2i \ end {bmatrix} [/ math] y [math] \ textbf {X} _ {2} = \ begin {bmatrix} 1 \\ -2i \ end {bmatrix} [/ math]
La solución general se da como [math] \ textbf {Y} = c_ {1} \ textbf {X} _ {1} e ^ {2it} + [/ math] [math] c_ {2} \ textbf {X} _ {2} e ^ {- 2it} [/ math] donde [math] c_1 [/ math] y [math] c_2 [/ math] son constantes arbitrarias. [matemáticas] [/ matemáticas]
En términos de sus componentes,
[matemáticas] y_1 = c_1 e ^ {2it} + c_2 e ^ {- 2it} \; \; \; \; [/matemáticas]
[matemáticas] y_2 = 2ic_1e ^ {2it} -2ic_2 e ^ {- 2it} (= y ‘_ {1}) [/ matemáticas]
Como sustituimos [matemáticas] \; \; \; \; y = y_1 [/ matemáticas]
La solución es por lo tanto
[matemáticas] \ displaystyle y = [/ matemáticas] [matemáticas] c_1e ^ {2it} + c_2e ^ {- 2it} [/ matemáticas]
Lo anterior se puede reducir a la combinación seno-coseno familiar usando la fórmula de Euler.
Las trayectorias para el sistema son elipses, con el punto central del sistema en el origen [matemática] (y_1 = 0, y_2 = 0) [/ matemática].
[Imagen: diagrama de fase generado en Matlab]
También puede usar Matlab para generar soluciones.
función dy = secondOrderLinear (t, y)
Función de espacio de estado%
dy = ceros (2,1);
dy (1) = y (2);
dy (2) = – 4 * y (1);
fin
Script% Solver
cierra todo; limpiar todo; clc;
[t, y] = ode45 (@secondOrderLinear, [0 100], [0 1]);
% de velocidad inicial que se supone que es 1 arbitrariamente
plot (t, y (:, 1)); cuadrícula activa;
[Imagen: Gráfico del tiempo de desplazamiento]