Al igual que con cualquier ODE de segundo orden no homogéneo, primero encontraremos la solución a la ecuación homogénea (llame a esto [matemáticas] y_h [/ matemáticas]), y luego encontraremos una solución particular (llame a esto [matemáticas] y_p [/ math]) a la ecuación no homogénea. Agregar estos dos produce la solución general de la ODE.
Paso uno: solución homogénea
Entonces, en este caso, nuestra ecuación característica / auxiliar / complementaria es [matemática] r ^ 2 + 1 = 0 [/ matemática], y sus raíces son [matemática] r = \ pm i [/ matemática], así que desde la solución a una EDO cuya ecuación característica tiene raíces complejas [matemáticas] \ alpha \ pm \ beta i [/ matemáticas] toma la forma [matemáticas] C_1e ^ {\ alpha x} \ cos (\ beta x) + C_2e ^ {\ alpha x } \ sin (\ beta x) [/ math],
[matemáticas] y_h (x) = C_1 \ cos (x) + C_2 \ sin (x) [/ matemáticas]
- ¿Cuál es la solución de la ecuación diferencial d ^ 2y / dx ^ 2 + y = secx?
- ¿Cuál será el grado y el orden de la ecuación diferencial cuando las potencias fraccionarias no se pueden eliminar como 2 (dy / dx) ^ (1/5) +5 (dy / dx) ^ (1/2) = (dy / dx) ^ (2/3)?
- ¿Cuál es la solución de la ecuación diferencial 3xydx + (1 + x ^ 2) dy = 0?
- ¿Cómo resuelvo [matemáticas] \ frac {| x + y |} {1+ | x + y |} <\ frac {| x |} {1+ | x |} + \ frac {| y |} {1 + | y |} [/ matemáticas]?
- ¿Las soluciones diferenciales parciales homogéneas lineales tienen una solución única, o dos o más?
Paso dos: solución particular
Ahora encontraremos una solución particular, y podemos hacerlo utilizando la variación de parámetros. Generalmente, la variación de parámetros es un método bastante poderoso porque produce una integral que generalmente es computable en términos de funciones elementales, e incluso si no lo es, casi siempre se puede expresar como una especie de función integral especial.
Recuerde que la solución particular, [math] y_p [/ math], toma la forma [math] v_1y_1 + v_2y_2 [/ math] donde
[matemáticas] \ displaystyle v_1 = \ int {\ dfrac {-y_2 f (x)} {aW [y_1, y_2]} \ dx} [/ math]
[matemáticas] \ displaystyle v_2 = \ int {\ dfrac {y_1 f (x)} {aW [y_1, y_2]} \ dx} [/ matemáticas]
[math] a [/ math] es el coeficiente de la segunda derivada en el ODE original, [math] f (x) [/ math] es el término forzado (el RHS del ODE original), [math] y_1 [/ math] y [math] y_2 [/ math] son las dos soluciones linealmente independientes, y [math] W [y_1, y_2] [/ math] es el Wronskian, que viene dado por
[matemáticas] \ begin {bmatrix} y_1 & y_2 \\ y_1 ‘& y_2’ \ end {bmatrix} = y_1 y_2 ‘- y_2 y_1’ [/ math]
Calculemos primero el Wronskian:
[matemáticas] \ begin {bmatrix} \ cos (x) & \ sin (x) \\ – \ sin (x) & \ cos (x) \ end {bmatrix} = \ cos ^ 2 (x) + \ sin ^ 2 (x) = 1 [/ matemáticas]
Ahora conectaremos todo a la fórmula integral para [math] v_1 [/ math]:
[matemáticas] \ displaystyle v_1 = \ int {-2 \ sin (x) e ^ x \ dx} [/ matemáticas]
Integrando por partes con [math] u = \ sin (x) [/ math] y [math] dv = e ^ x [/ math] produce
[matemáticas] \ displaystyle v_1 = -2 \ sin (x) e ^ x + 2 \ int {\ cos (x) e ^ x \ dx} [/ matemáticas]
La segunda integral se puede integrar por partes nuevamente con [math] s = \ cos (x) [/ math] y [math] dt = e ^ x [/ math]:
[matemáticas] \ displaystyle v_1 = -2 \ sin (x) e ^ x + 2 \ cos (x) e ^ x – 2 \ int {\ sin (x) e ^ x \ dx} [/ matemáticas]
Tenga en cuenta que la integral que tenemos ahora es exactamente la misma que [math] v_1 [/ math], entonces:
[matemáticas] v_1 = -2 \ sin (x) e ^ x + 2 \ cos (x) e ^ x – v_1 [/ matemáticas]
[matemáticas] \ implica v_1 = e ^ x (\ cos (x) – \ sin (x)) [/ matemáticas]
A continuación, intentaremos encontrar [math] v_2 [/ math]:
[matemáticas] \ displaystyle v_2 = \ int {2 \ cos (x) e ^ x \ dx} [/ matemáticas]
Podemos usar un método muy similar para resolver esta integral (integrar por partes dos veces hasta que la integral se vea como la integral original, luego reorganizar), y resulta ser:
[matemáticas] v_2 = e ^ x (\ sin (x) + \ cos (x)) [/ matemáticas]
Finalmente, dado que ahora sabemos [matemáticas] v_1 [/ matemáticas] y [matemáticas] v_2 [/ matemáticas], podemos escribir la solución particular, [matemáticas] y_p [/ matemáticas]:
[matemáticas] y_p = e ^ x \ cos (x) (\ cos (x) – \ sin (x)) + e ^ x \ sin (x) (\ sin (x) + \ cos (x)) [/ matemáticas]
Sin embargo, podemos distribuir para simplificar:
[matemáticas] y_p = e ^ x (\ cos ^ 2 (x) – \ cos (x) \ sin (x) + \ sin ^ 2 (x) + \ sin (x) \ cos (x)) [/ matemáticas ]
[matemáticas] \ implica y_p = e ^ x [/ matemáticas]
Paso tres: principio de superposición
Ahora, agregaremos las soluciones particulares y homogéneas para obtener la solución general:
[matemáticas] y (x) = y_h + y_p = \ en caja {C_1 \ cos (x) + C_2 \ sin (x) + e ^ x} [/ matemáticas]