ING. EDGAR H.

MEDINA TAPIA MAT 1207

Tema 2.- Ecuaciones Diferenciales Ordinarias Lineales de Orden n

2.1.- Preliminares

Se denomina EDO lineal de orden 𝑛𝑛 con Coeficientes Variables 𝐴𝐴𝑖𝑖 (𝑥𝑥) 𝑐𝑐𝑐𝑐𝑐𝑐 𝑖𝑖 = 0, 1,2, ⋯ , 𝑛𝑛; a toda ecuación de la
forma:

𝐴𝐴𝑛𝑛 (𝑥𝑥)𝑦𝑦 (𝑛𝑛) + 𝐴𝐴𝑛𝑛−1 (𝑥𝑥)𝑦𝑦 (𝑛𝑛−1) + ⋯ + 𝐴𝐴1 (𝑥𝑥)𝑦𝑦 ′ + 𝐴𝐴0 (𝑥𝑥)𝑦𝑦 = 𝑅𝑅(𝑥𝑥) (2.1.1)

𝑑𝑑 𝑘𝑘
Donde 𝑦𝑦 (𝑘𝑘) = 𝑑𝑑𝑥𝑥 𝑘𝑘 𝑘𝑘 = 0, 1,2, ⋯ , 𝑛𝑛 𝑦𝑦 𝑅𝑅(𝑥𝑥) es una función de 𝑥𝑥

Si 𝑅𝑅(𝑥𝑥) = 0 ; la EDO recibe el nombre de Reducida o Incompleta (homogénea) de orden 𝑛𝑛

Además, si las 𝐴𝐴𝑖𝑖 (𝑥𝑥) = 𝑎𝑎𝑖𝑖 son constantes la Edo es lineal con Coeficientes Constantes

Representación de una EDO Lineal por medio de los Operadores 𝑫𝑫𝒌𝒌y 𝑳𝑳(𝑫𝑫) .- Si se emplea el operador
𝑑𝑑 𝑘𝑘
derivada 𝐷𝐷𝑘𝑘 = 𝑑𝑑𝑑𝑑 𝑘𝑘 ; 𝑘𝑘 = 0,1,2, ⋯ , 𝑛𝑛, la EDO de enésimo orden puede ser escrita como

𝐴𝐴𝑛𝑛 (𝑥𝑥)𝐷𝐷𝑛𝑛 𝑦𝑦 + 𝐴𝐴𝑛𝑛−1 (𝑥𝑥)𝐷𝐷𝑛𝑛−1 𝑦𝑦 + ⋯ + 𝐴𝐴1 (𝑥𝑥)𝐷𝐷𝐷𝐷 + 𝐴𝐴0 (𝑥𝑥)𝑦𝑦 = 𝑅𝑅(𝑥𝑥) (2.1.2)

𝐿𝐿(𝐷𝐷)[𝑦𝑦] = [𝐴𝐴𝑛𝑛 (𝑥𝑥)𝐷𝐷𝑛𝑛 + 𝐴𝐴𝑛𝑛−1 (𝑥𝑥)𝐷𝐷𝑛𝑛−1 + ⋯ + 𝐴𝐴1 (𝑥𝑥)𝐷𝐷 + 𝐴𝐴0 (𝑥𝑥)]{𝑦𝑦} = 𝑅𝑅(𝑥𝑥) (2.1.3)

Donde 𝐿𝐿(𝐷𝐷) = 𝐴𝐴𝑛𝑛 (𝑥𝑥)𝐷𝐷𝑛𝑛 + 𝐴𝐴𝑛𝑛−1 (𝑥𝑥)𝐷𝐷𝑛𝑛−1 + ⋯ + 𝐴𝐴1 (𝑥𝑥)𝐷𝐷 + 𝐴𝐴0 (𝑥𝑥) es el operador diferencial polinomial de orden 𝑛𝑛

Solución de una EDO Lineal

Toda función 𝑦𝑦 = 𝑢𝑢(𝑥𝑥) continua y diferenciable en la región de integración, será solución de la ecuación lineal, si
esta función y sus derivadas 𝑢𝑢′ , 𝑢𝑢′′ , ⋯ , 𝑢𝑢(𝑛𝑛)transforman la ecuación dada en una identidad.

Propiedades de las Soluciones

1. Si 𝑦𝑦1 es una solución de la ecuación reducida 𝐿𝐿(𝐷𝐷)[𝑦𝑦] = 0 , entonces 𝑐𝑐1 𝑦𝑦1 es también solución.
2. Principio de Superposición. - Si se tiene un conjunto de 𝑛𝑛 soluciones {𝑦𝑦1, 𝑦𝑦2, ⋯ , 𝑦𝑦𝑛𝑛 } de la ecuación
reducida 𝐿𝐿(𝐷𝐷)[𝑦𝑦] = 0 , todas Linealmente independientes entre sí, entonces para toda 𝑐𝑐𝑖𝑖 𝑐𝑐𝑐𝑐𝑐𝑐 𝑖𝑖 =
1,2, ⋯ 𝑛𝑛; la combinación lineal
𝑦𝑦𝑐𝑐 = 𝑐𝑐1 𝑦𝑦1 + 𝑐𝑐2 𝑦𝑦2 + ⋯ + 𝑐𝑐𝑛𝑛 𝑦𝑦𝑛𝑛
Es también una solución y se denomina función o solución complementaria.
3. La solución general de la ecuación completa 𝐿𝐿(𝐷𝐷)[𝑦𝑦] = 𝑅𝑅(𝑥𝑥) es igual a
𝑦𝑦 = 𝑦𝑦𝑐𝑐 + 𝑦𝑦𝑝𝑝
Donde 𝑦𝑦𝑝𝑝 es la solución particular de 𝐿𝐿(𝐷𝐷)[𝑦𝑦] = 𝑅𝑅(𝑥𝑥)
4. El Wronskiano. - Se denomina el Wronskiano de un conjunto de 𝑛𝑛 soluciones {𝑦𝑦1 , 𝑦𝑦2 , ⋯ , 𝑦𝑦𝑛𝑛 } de la
ecuación reducida 𝐿𝐿(𝐷𝐷)[𝑦𝑦] = 0 a la determinante:

𝑦𝑦1 𝑦𝑦2 ⋯ 𝑦𝑦𝑛𝑛

𝑦𝑦1 ′ 𝑦𝑦2 ′ ⋯ 𝑦𝑦𝑛𝑛 ′
𝑊𝑊[𝑦𝑦1 , 𝑦𝑦2 , ⋯ , 𝑦𝑦𝑛𝑛 ]=� �
⋮ ⋮ ⋱ ⋮
𝑦𝑦1 (𝑛𝑛−1) 𝑦𝑦2 (𝑛𝑛−1) ⋯ 𝑦𝑦𝑛𝑛 (𝑛𝑛−1)

13
ING. EDGAR H. MEDINA TAPIA MAT 1207

Teorema. - Si las funciones 𝐴𝐴𝑖𝑖 (𝑥𝑥) de la ecuación lineal de orden n son continuas en la región de integración,
entonces las 𝑛𝑛 soluciones 𝑦𝑦1, 𝑦𝑦2 , ⋯ , 𝑦𝑦𝑛𝑛 son linealmente independientes (L.I.) si y solo sí su Wronskiano es distinto
de cero.

Observación. - Una EDO lineal de orden n 𝐿𝐿(𝐷𝐷)[𝑦𝑦] = 0 tiene exactamente n soluciones L.I.

2.2.- EDO Lineal Reducida de Orden n con Coeficientes Constantes

Sea la ecuación lineal:

𝐿𝐿(𝐷𝐷)[𝑦𝑦] = [𝑎𝑎𝑛𝑛 𝐷𝐷𝑛𝑛 + 𝑎𝑎𝑛𝑛−1 𝐷𝐷𝑛𝑛−1 + ⋯ + 𝑎𝑎1 𝐷𝐷 + 𝑎𝑎0 ]{𝑦𝑦} = 0 (2.2.1)

Donde las 𝑎𝑎𝑖𝑖 son constantes.

Una solución no trivial 𝑦𝑦 ≠ 0 de la ecuación reducida tiene la forma 𝑦𝑦 = 𝑒𝑒 𝜆𝜆𝜆𝜆 .

Ahora bien, cuando se sustituye esta solución en la ecuación lineal, se obtiene:

𝐿𝐿(𝐷𝐷)�𝑒𝑒 𝜆𝜆𝜆𝜆 � = 𝑎𝑎𝑛𝑛 𝜆𝜆𝑛𝑛 𝑒𝑒 𝜆𝜆𝜆𝜆 + 𝑎𝑎𝑛𝑛−1 𝜆𝜆𝑛𝑛−1 𝑒𝑒 𝜆𝜆𝜆𝜆 + ⋯ + 𝑎𝑎1 𝜆𝜆𝑒𝑒 𝜆𝜆𝜆𝜆 + 𝑎𝑎0 𝑒𝑒 𝜆𝜆𝜆𝜆 = 0

Factorizando 𝑒𝑒 𝜆𝜆𝜆𝜆 , resulta (𝑎𝑎𝑛𝑛 𝜆𝜆𝑛𝑛 + 𝑎𝑎𝑛𝑛−1 𝜆𝜆𝑛𝑛−1 + ⋯ + 𝑎𝑎1 𝜆𝜆 + 𝑎𝑎0 )𝑒𝑒 𝜆𝜆𝜆𝜆 = 0

Como 𝑦𝑦 = 𝑒𝑒 𝜆𝜆𝜆𝜆 ≠ 0 entonces, se tiene una ecuación polinomial de grado 𝑛𝑛, denominada ecuación característica.

𝑎𝑎𝑛𝑛 𝜆𝜆𝑛𝑛 + 𝑎𝑎𝑛𝑛−1 𝜆𝜆𝑛𝑛−1 + ⋯ + 𝑎𝑎1 𝜆𝜆 + 𝑎𝑎0 = 0

El tipo de soluciones que se obtengan de esta ecuación, determina también el tipo de soluciones de la ecuación
diferencial reducida 𝐿𝐿(𝐷𝐷)[𝑦𝑦] = 0, por esta razón se tienen los siguientes casos:

1º La Ecuación Característica, tiene Raíces Reales y Distintas

Sean 𝜆𝜆1 , 𝜆𝜆2 , ⋯ , 𝜆𝜆𝑛𝑛 las raíces de la ecuación característica, todas ellas reales y distintas; en estas condiciones, las
n soluciones de la ecuación reducida son 𝑦𝑦1 = 𝑒𝑒 𝜆𝜆1𝑥𝑥 ; 𝑦𝑦2 = 𝑒𝑒 𝜆𝜆2𝑥𝑥 ; ⋯ ; 𝑦𝑦𝑛𝑛 = 𝑒𝑒 𝜆𝜆𝑛𝑛 𝑥𝑥 , luego, la función o solución
complementaria será igual:

𝑦𝑦𝑐𝑐 = 𝑐𝑐1 𝑒𝑒 𝜆𝜆1𝑥𝑥 + 𝑐𝑐2 𝑒𝑒 𝜆𝜆2𝑥𝑥 + ⋯ + 𝑐𝑐𝑛𝑛 𝑒𝑒 𝜆𝜆𝑛𝑛𝑥𝑥

2º La Ecuación Característica, tiene Raíces Complejas y Conjugadas

A cada par conjugado de raíces complejas de la ecuación característica 𝜆𝜆1 = 𝛼𝛼 + 𝑖𝑖𝑖𝑖 y 𝜆𝜆2 = 𝛼𝛼 − 𝑖𝑖𝑖𝑖, le
corresponde también un par de soluciones para la EDO lineal de las formas 𝑦𝑦1 = 𝑒𝑒 (𝛼𝛼+𝑖𝑖𝑖𝑖) 𝑥𝑥 y 𝑦𝑦2 = 𝑒𝑒 (𝛼𝛼−𝑖𝑖𝑖𝑖) 𝑥𝑥 , la
combinación de estas soluciones es también una solución, es decir

𝑦𝑦 = 𝑘𝑘1 𝑦𝑦1 + 𝑘𝑘2 𝑦𝑦2 = 𝑒𝑒 𝛼𝛼𝛼𝛼 �𝑘𝑘1 𝑒𝑒 𝑖𝑖𝑖𝑖𝑖𝑖 + 𝑘𝑘2 𝑒𝑒 −𝑖𝑖𝑖𝑖𝑖𝑖 �

Si se aplican las relaciones de Euler 𝑒𝑒 ±𝑖𝑖𝑖𝑖𝑖𝑖 = cos 𝛽𝛽𝛽𝛽 ± 𝑖𝑖 sin 𝛽𝛽𝛽𝛽 , resulta:

𝑦𝑦 = 𝑒𝑒 𝛼𝛼𝛼𝛼 [𝑘𝑘1 (cos 𝛽𝛽𝛽𝛽 + 𝑖𝑖 sin 𝛽𝛽𝛽𝛽) + 𝑘𝑘2 (cos 𝛽𝛽𝛽𝛽 − 𝑖𝑖 sin 𝛽𝛽𝛽𝛽)]

Reorganizando este resultado, se tiene

𝑦𝑦 = 𝑒𝑒 𝛼𝛼𝛼𝛼 [(𝑘𝑘1 + 𝑘𝑘2 ) cos 𝛽𝛽𝛽𝛽 + 𝑖𝑖 (𝑘𝑘1 − 𝑘𝑘2 ) sin 𝛽𝛽𝛽𝛽]

14
ING. EDGAR H. MEDINA TAPIA MAT 1207

Denominando 𝑐𝑐1 = (𝑘𝑘1 + 𝑘𝑘2 ) ; 𝑐𝑐2 = 𝑖𝑖(𝑘𝑘1 − 𝑘𝑘2 ) se obtiene finalmente la solución:

𝑦𝑦 = 𝑒𝑒 𝛼𝛼𝛼𝛼 (𝑐𝑐1 cos 𝛽𝛽𝛽𝛽 + 𝑐𝑐2 sin 𝛽𝛽𝛽𝛽)

3º La Ecuación Característica, tiene Raíces Múltiples o Repetidas

Si 𝜆𝜆 = 𝜆𝜆𝑚𝑚 es una raíz de multiplicidad 𝑚𝑚 de la ecuación característica y si 𝜆𝜆𝑚𝑚+1 , ⋯ , 𝜆𝜆𝑛𝑛 son las demás raíces
todas distintas de 𝜆𝜆𝑚𝑚 ; entonces para una solución de la forma 𝑦𝑦 = 𝑒𝑒 𝜆𝜆𝜆𝜆 , se tiene

𝐿𝐿(𝐷𝐷)�𝑒𝑒 𝜆𝜆𝜆𝜆 � = (𝑎𝑎𝑛𝑛 𝜆𝜆𝑛𝑛 + 𝑎𝑎𝑛𝑛−1 𝜆𝜆𝑛𝑛−1 + ⋯ + 𝑎𝑎1 𝜆𝜆 + 𝑎𝑎0 )𝑒𝑒 𝜆𝜆𝜆𝜆 = (𝜆𝜆 − 𝜆𝜆𝑚𝑚 )𝑚𝑚 ℎ(𝜆𝜆)𝑒𝑒 𝜆𝜆𝜆𝜆

Derivando esta relación parcialmente respecto de 𝜆𝜆 , resulta

𝜕𝜕 𝜕𝜕
�𝐿𝐿(𝐷𝐷)�𝑒𝑒 𝜆𝜆𝜆𝜆 �� = 𝑚𝑚(𝜆𝜆 − 𝜆𝜆𝑚𝑚 )𝑚𝑚−1 ℎ(𝜆𝜆)𝑒𝑒 𝜆𝜆𝜆𝜆 + (𝜆𝜆 − 𝜆𝜆𝑚𝑚 )𝑚𝑚 �ℎ(𝜆𝜆)𝑒𝑒 𝜆𝜆𝜆𝜆 �
𝜕𝜕𝜆𝜆 𝜕𝜕𝜆𝜆

Recordando que: 𝜕𝜕 𝜕𝜕
�𝐿𝐿(𝐷𝐷)�𝑒𝑒 𝜆𝜆𝜆𝜆 �� = 𝐿𝐿(𝐷𝐷) � �𝑒𝑒 𝜆𝜆𝜆𝜆 �� = 𝐿𝐿(𝐷𝐷)�𝑥𝑥𝑒𝑒 𝜆𝜆𝜆𝜆 �
𝜕𝜕𝜆𝜆 𝜕𝜕𝜆𝜆
Y sustituyendo 𝜆𝜆 = 𝜆𝜆𝑚𝑚 , se obtiene 𝐿𝐿(𝐷𝐷)�𝑥𝑥𝑒𝑒 𝜆𝜆𝑚𝑚𝑥𝑥 � = 0

Esto significa que 𝑥𝑥𝑒𝑒 𝜆𝜆𝑚𝑚𝑥𝑥 es también una solución. Si se repite esta operación (𝑚𝑚 − 1) veces se obtienen las
soluciones 𝑥𝑥𝑒𝑒 𝜆𝜆𝑚𝑚𝑥𝑥 , 𝑥𝑥 2 𝑒𝑒 𝜆𝜆𝑚𝑚 𝑥𝑥 , ⋯ , 𝑥𝑥 𝑚𝑚−1 𝑒𝑒 𝜆𝜆𝑚𝑚𝑥𝑥 , todas linealmente independientes. Por lo tanto, la solución de la
ecuación para una raíz múltiple 𝜆𝜆 = 𝜆𝜆𝑚𝑚 está dada por la fórmula

𝑦𝑦 = (𝑐𝑐1 + 𝑐𝑐2 𝑥𝑥 + 𝑐𝑐3 𝑥𝑥 2 + ⋯ + 𝑐𝑐𝑚𝑚 𝑥𝑥 𝑚𝑚−1 )𝑒𝑒 𝜆𝜆𝑚𝑚𝑥𝑥

2.3.- EDO Lineal Completa de Orden n con Coeficientes Constantes

Sea la ecuación:

𝐿𝐿(𝐷𝐷)[𝑦𝑦] = [𝑎𝑎𝑛𝑛 𝐷𝐷𝑛𝑛 + 𝑎𝑎𝑛𝑛−1 𝐷𝐷𝑛𝑛−1 + ⋯ + 𝑎𝑎1 𝐷𝐷 + 𝑎𝑎0 ]{𝑦𝑦} = 𝑅𝑅(𝑥𝑥) (2.3.1.1)
Donde las 𝑎𝑎𝑖𝑖 son constantes.

La solución general de la ecuación completa está dada por combinación 𝑦𝑦 = 𝑦𝑦𝑐𝑐 + 𝑦𝑦𝑝𝑝 donde 𝑦𝑦𝑐𝑐 es la función
complementaria, 𝑦𝑦𝑝𝑝 es la solución particular.

Teorema. - Sean 𝑦𝑦𝑝𝑝1 ; 𝑦𝑦𝑝𝑝2 ; ⋯ ; 𝑦𝑦𝑝𝑝𝑝𝑝 soluciones particulares de la EDO lineal completa de orden n correspondientes
a las funciones 𝑅𝑅1 ; 𝑅𝑅2 ; ⋯ ; 𝑅𝑅𝑘𝑘 ; entonces la combinación

𝑦𝑦𝑝𝑝 = 𝑦𝑦𝑝𝑝1 + 𝑦𝑦𝑝𝑝2 + ⋯ + 𝑦𝑦𝑝𝑝𝑝𝑝

Es la solución particular de la ecuación completa 𝐿𝐿(𝐷𝐷)[𝑦𝑦] = 𝑅𝑅1 (𝑥𝑥) + 𝑅𝑅2 (𝑥𝑥) + ⋯ + 𝑅𝑅𝑛𝑛 (𝑥𝑥)

La solución particular 𝑦𝑦𝑝𝑝 de la EDO completa 𝐿𝐿(𝐷𝐷)[𝑦𝑦] = 𝑅𝑅(𝑥𝑥) puede ser obtenida por medio de la aplicación de
los siguientes métodos:

15
ING. EDGAR H. MEDINA TAPIA MAT 1207

2.3.1.- El Método de los Coeficientes Indeterminados

Sea la ecuación diferencial lineal de orden n

𝐿𝐿(𝐷𝐷)[𝑦𝑦] = [𝑎𝑎𝑛𝑛 𝐷𝐷𝑛𝑛 + 𝑎𝑎𝑛𝑛−1 𝐷𝐷𝑛𝑛−1 + ⋯ + 𝑎𝑎1 𝐷𝐷 + 𝑎𝑎0 ]{𝑦𝑦} = 𝑅𝑅(𝑥𝑥) (2.3.1.2)
Este método es aplicable únicamente cuando la función 𝑅𝑅(𝑥𝑥) tiene en forma parcial o total la siguiente
configuración

𝑅𝑅(𝑥𝑥) = 𝑒𝑒 𝛼𝛼𝛼𝛼 [𝑃𝑃𝑟𝑟 (𝑥𝑥) 𝑐𝑐𝑐𝑐𝑐𝑐 𝛽𝛽𝛽𝛽 + 𝑄𝑄𝑠𝑠 (𝑥𝑥) 𝑠𝑠𝑠𝑠𝑠𝑠 𝛽𝛽𝛽𝛽]

En este caso 𝑃𝑃𝑟𝑟 (𝑥𝑥) y𝑄𝑄𝑠𝑠 (𝑥𝑥) son polinomios de grados 𝑟𝑟 y 𝑠𝑠 respectivamente. El método plantea una solución
particular 𝑦𝑦𝑝𝑝 de la forma

𝑦𝑦𝑝𝑝 = 𝑥𝑥 𝑘𝑘 𝑅𝑅�(𝑥𝑥) = 𝑥𝑥 𝑘𝑘 𝑒𝑒 𝛼𝛼𝛼𝛼 �𝑃𝑃�𝑡𝑡 (𝑥𝑥) 𝑐𝑐𝑐𝑐𝑐𝑐 𝛽𝛽𝛽𝛽 + 𝑄𝑄�𝑡𝑡 (𝑥𝑥) 𝑠𝑠𝑠𝑠𝑠𝑠 𝛽𝛽𝛽𝛽�

Dónde: 𝑃𝑃�𝑡𝑡 (𝑥𝑥) , 𝑄𝑄�𝑡𝑡 (𝑥𝑥)Ambos son polinomios de grado 𝑡𝑡; siendo 𝑡𝑡 = 𝑚𝑚á𝑥𝑥(𝑟𝑟, 𝑠𝑠) el máximo valor entre los grados de
los polinomios 𝑃𝑃�𝑡𝑡 (𝑥𝑥) y 𝑄𝑄�𝑡𝑡. El parámetro 𝑘𝑘 Es el menor entero positivo que elimina términos semejantes entre la
solución particular y la función complementaria y coincide con la multiplicidad de la solución complementaria.

Los coeficientes de los polinomios 𝑃𝑃�𝑡𝑡 (𝑥𝑥) y 𝑄𝑄�𝑡𝑡 se obtienen cuando se sustituye la solución particular planteada 𝑦𝑦𝑝𝑝
en la ecuación diferencial.

2.3.2.- El método del Operador Inverso

Sea la ecuación diferencial lineal de orden n con coeficientes constantes

𝐿𝐿(𝐷𝐷)[𝑦𝑦] = [𝑎𝑎𝑛𝑛 𝐷𝐷𝑛𝑛 + 𝑎𝑎𝑛𝑛−1 𝐷𝐷𝑛𝑛−1 + ⋯ + 𝑎𝑎1 𝐷𝐷 + 𝑎𝑎0 ]{𝑦𝑦} = 𝑅𝑅(𝑥𝑥) (2.3.2.1)

𝑑𝑑 𝑘𝑘
Donde 𝐷𝐷𝑘𝑘 = 𝑑𝑑𝑑𝑑 𝑘𝑘 ; 𝑘𝑘 = 0,1,2, ⋯ , 𝑛𝑛 . En estas condiciones, la solución particular de la EDO formalmente puede
ser escrita como

1
𝑦𝑦𝑝𝑝 = {𝑅𝑅(𝑥𝑥)} (2.3.2.2)
𝐿𝐿(𝐷𝐷)
1
Donde 𝐿𝐿(𝐷𝐷) es el Operador Inverso.

Propiedades del Operador 𝑳𝑳(𝑫𝑫). - Sean 𝐿𝐿𝑛𝑛 (𝐷𝐷) 𝑦𝑦 𝐿𝐿𝑚𝑚 (𝐷𝐷) dos operadores polinomiales diferenciales de órdenes
𝑛𝑛 𝑦𝑦 𝑚𝑚 respectivamente con 𝑎𝑎 𝑦𝑦 𝑏𝑏 constantes, entonces

1.- 𝐿𝐿𝑛𝑛 (𝐷𝐷)[𝑎𝑎𝑎𝑎(𝑥𝑥) + 𝑏𝑏𝑏𝑏(𝑥𝑥)] = 𝑎𝑎𝐿𝐿𝑛𝑛 (𝐷𝐷){𝑓𝑓(𝑥𝑥)} + 𝑏𝑏𝐿𝐿𝑛𝑛 (𝐷𝐷){𝑔𝑔(𝑥𝑥)}

2.- [𝐿𝐿𝑚𝑚 (𝐷𝐷) + 𝐿𝐿𝑛𝑛 (𝐷𝐷)]{𝑓𝑓(𝑥𝑥)} = [𝐿𝐿𝑛𝑛 (𝐷𝐷) + 𝐿𝐿𝑚𝑚 (𝐷𝐷)]{𝑓𝑓(𝑥𝑥)} = 𝐿𝐿𝑚𝑚 (𝐷𝐷){𝑓𝑓(𝑥𝑥)} + 𝐿𝐿𝑛𝑛 (𝐷𝐷){𝑓𝑓(𝑥𝑥)}

3.- [𝐿𝐿𝑛𝑛 (𝐷𝐷)𝐿𝐿𝑚𝑚 (𝐷𝐷)]{𝑓𝑓(𝑥𝑥)} = 𝐿𝐿𝑛𝑛 (𝐷𝐷)[𝐿𝐿𝑚𝑚 (𝐷𝐷){𝑓𝑓(𝑥𝑥)}] = 𝐿𝐿𝑚𝑚 (𝐷𝐷)[𝐿𝐿𝑛𝑛 (𝐷𝐷){𝑓𝑓(𝑥𝑥)}]

4.- 𝐿𝐿(𝐷𝐷){𝑒𝑒 𝛼𝛼𝛼𝛼 𝑓𝑓(𝑥𝑥)} = 𝑒𝑒 𝛼𝛼𝛼𝛼 𝐿𝐿(𝐷𝐷 + 𝛼𝛼){𝑓𝑓(𝑥𝑥)}

16
ING. EDGAR H. MEDINA TAPIA MAT 1207

𝟏𝟏 1 1
Propiedades del Operador Inverso .- Sean 𝐿𝐿 𝑦𝑦 dos operadores inversos, entonces,
𝑳𝑳(𝑫𝑫) 𝑚𝑚 (𝐷𝐷) 𝐿𝐿𝑛𝑛 (𝐷𝐷)

1 1 1 1
1. - �
𝐿𝐿𝑚𝑚(𝐷𝐷) 𝐿𝐿𝑛𝑛 (𝐷𝐷)
{𝑅𝑅(𝑥𝑥)}� = �
𝐿𝐿𝑚𝑚(𝐷𝐷) 𝐿𝐿𝑛𝑛 (𝐷𝐷)
{𝑅𝑅(𝑥𝑥)}�

1 1 1
2. - 𝐿𝐿𝑚𝑚 (𝐷𝐷)
{𝑅𝑅(𝑥𝑥)} =
𝐿𝐿𝑛𝑛 (𝐷𝐷)𝐿𝐿𝑛𝑛 (𝐷𝐷)
[𝐿𝐿𝑛𝑛 (𝐷𝐷){𝑅𝑅(𝑥𝑥)}] = 𝐿𝐿𝑛𝑛 (𝐷𝐷) �
𝐿𝐿𝑚𝑚(𝐷𝐷)𝐿𝐿𝑛𝑛 (𝐷𝐷)
{𝑅𝑅(𝑥𝑥)}�

1 1
3. - 𝐿𝐿 {𝑒𝑒 𝛼𝛼𝛼𝛼 𝐺𝐺(𝑥𝑥)} = 𝑒𝑒 𝛼𝛼𝛼𝛼 {𝐺𝐺(𝑥𝑥)} 𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒
𝑛𝑛 (𝐷𝐷) 𝐿𝐿𝑛𝑛 (𝐷𝐷+𝛼𝛼)

1 1 𝐿𝐿′ (𝐷𝐷)
4. - 𝐿𝐿(𝐷𝐷)
{𝑥𝑥𝑥𝑥(𝑥𝑥)} = 𝑥𝑥
𝐿𝐿(𝐷𝐷)
{𝐺𝐺(𝑥𝑥)} −
[𝐿𝐿(𝐷𝐷)]2
{𝐺𝐺(𝑥𝑥)} 𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 𝑑𝑑𝑑𝑑 𝑥𝑥

2.3.2.1. - Integración Mediante el Operador Inverso por definición

Primer Caso

1 1
𝑦𝑦𝑝𝑝 = {𝑅𝑅(𝑥𝑥)} = 𝑘𝑘 {𝑅𝑅(𝑥𝑥)}
𝐿𝐿(𝐷𝐷) 𝐷𝐷
1 𝑑𝑑𝑑𝑑
Para 𝒌𝒌 = 𝟏𝟏.- 𝑦𝑦 = 𝐷𝐷 {𝑅𝑅(𝑥𝑥)} es la solución particular de la EDO 𝑑𝑑𝑑𝑑
= 𝑅𝑅(𝑥𝑥) expresada también como 𝑦𝑦 = ∫ 𝑅𝑅(𝑥𝑥)𝑑𝑑𝑑𝑑
por consiguiente, se tiene la siguiente equivalencia

1
𝑦𝑦𝑝𝑝 = {𝑅𝑅(𝑥𝑥)} = � 𝑅𝑅(𝑥𝑥)𝑑𝑑𝑑𝑑 (2.3.2.3)
𝐷𝐷

Para 𝒌𝒌 = 𝒏𝒏 > 𝟏𝟏. - Se realiza la operación anterior en forma repetida, obteniendo la equivalencia

1
𝑦𝑦𝑝𝑝 = {𝑅𝑅(𝑥𝑥)} = � ⋯ � 𝑅𝑅(𝑥𝑥) 𝑑𝑑𝑥𝑥 𝑛𝑛 (2.3.2.4)
𝐷𝐷𝑛𝑛
Segundo Caso

1 1
𝑦𝑦𝑝𝑝 = {𝑅𝑅(𝑥𝑥)} = {𝑅𝑅(𝑥𝑥)}
𝐿𝐿(𝐷𝐷) (𝐷𝐷 − 𝜆𝜆1 )(𝐷𝐷 − 𝜆𝜆2 ) ⋯ (𝐷𝐷 − 𝜆𝜆𝑛𝑛 )

1 𝑑𝑑𝑑𝑑
Para 𝒌𝒌 = 𝟏𝟏.- 𝑦𝑦 = (𝐷𝐷−𝜆𝜆 ) {𝑅𝑅(𝑥𝑥)} es la solución particular de la ecuación 𝑑𝑑𝑑𝑑
− 𝜆𝜆1 𝑥𝑥 = 𝑅𝑅(𝑥𝑥) escrita también como
1
𝑦𝑦 = 𝑒𝑒 𝜆𝜆1𝑥𝑥 ∫ 𝑒𝑒 −𝜆𝜆1𝑥𝑥 𝑅𝑅(𝑥𝑥)𝑑𝑑𝑑𝑑; como consecuencia se verifica

1
𝑦𝑦𝑝𝑝 = {𝑅𝑅(𝑥𝑥)} = 𝑒𝑒 𝜆𝜆1𝑥𝑥 � 𝑒𝑒 −𝜆𝜆1 𝑥𝑥 𝑅𝑅(𝑥𝑥)𝑑𝑑𝑑𝑑 (2.3.2.5)
(𝐷𝐷 − 𝜆𝜆1 )

Para 𝒌𝒌 = 𝒏𝒏 > 𝟏𝟏. - La solución particular, puede ser obtenida mediante los siguientes procedimientos

1º Procedimiento. – Se aplica en forma repetida la equivalencia del caso 𝑘𝑘 = 1, esto es

1
𝑦𝑦𝑝𝑝 = {𝑅𝑅(𝑥𝑥)} = 𝑒𝑒 𝜆𝜆𝑛𝑛 𝑥𝑥 � 𝑒𝑒 (𝜆𝜆𝑛𝑛−1−𝜆𝜆𝑛𝑛 )𝑥𝑥 ⋯ 𝑒𝑒 𝜆𝜆2𝑥𝑥 � 𝑒𝑒 (𝜆𝜆1−𝜆𝜆2)𝑥𝑥 � 𝑒𝑒 −𝜆𝜆1𝑥𝑥 𝑅𝑅(𝑥𝑥)𝑑𝑑𝑥𝑥 𝑛𝑛 (2.3.2.5)
(𝐷𝐷 − 𝜆𝜆𝑛𝑛 ) ⋯ (𝐷𝐷 − 𝜆𝜆2 )(𝐷𝐷 − 𝜆𝜆1 )

17
ING. EDGAR H. MEDINA TAPIA MAT 1207

𝟏𝟏
2º Procedimiento. – Se desarrolla el operador 𝑳𝑳(𝑫𝑫)
en fracciones parciales y luego aplicar a cada fracción la
equivalencia para el caso 𝑘𝑘 = 1, es decir

1 1 𝐴𝐴1 𝐴𝐴2 𝐴𝐴𝑛𝑛

= = + + ⋯+
𝐿𝐿(𝐷𝐷) (𝐷𝐷 − 𝜆𝜆1 )(𝐷𝐷 − 𝜆𝜆2 ) ⋯ (𝐷𝐷 − 𝜆𝜆𝑛𝑛 ) 𝐷𝐷 − 𝜆𝜆1 𝐷𝐷 − 𝜆𝜆2 𝐷𝐷 − 𝜆𝜆𝑛𝑛

1 1
𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷, 𝐴𝐴𝑘𝑘 = � =
𝐿𝐿′(𝐷𝐷) 𝐷𝐷=𝜆𝜆𝑘𝑘 𝐿𝐿′(𝜆𝜆𝑘𝑘 )

Luego, la solución particular 𝑦𝑦𝑝𝑝 será igual a

𝑛𝑛 𝑛𝑛
1 𝐴𝐴𝑘𝑘
𝑦𝑦𝑝𝑝 = {𝑅𝑅(𝑥𝑥)} = � {𝑅𝑅(𝑥𝑥)} = � 𝐴𝐴𝑘𝑘 𝑒𝑒 𝜆𝜆𝑘𝑘𝑥𝑥 � 𝑒𝑒 −𝜆𝜆𝑘𝑘𝑥𝑥 𝑅𝑅(𝑥𝑥)𝑑𝑑𝑑𝑑 (2.3.2.5)
𝐿𝐿(𝐷𝐷) (𝐷𝐷 − 𝜆𝜆𝑘𝑘 )
𝑘𝑘=1 𝑘𝑘=1

2.3.2.2.- Integración por los Métodos Abreviados del Operador Inverso

La aplicación del método del Operador Inverso se simplifica en los siguientes casos:

(1) Si 𝑅𝑅(𝑥𝑥) = 𝐴𝐴𝐴𝐴 𝛼𝛼𝛼𝛼

La solución particular para [𝐿𝐿(𝐷𝐷)]{𝑦𝑦} = 𝐴𝐴𝐴𝐴 𝛼𝛼𝛼𝛼 se obtiene por medio de las siguientes fórmulas

Para 𝐿𝐿(𝛼𝛼) ≠ 0

1 𝐴𝐴 𝛼𝛼𝛼𝛼
𝑦𝑦𝑝𝑝 = {𝐴𝐴𝑒𝑒 𝛼𝛼𝛼𝛼 } = 𝑒𝑒 (2.3.2.6)
𝐿𝐿(𝐷𝐷) 𝐿𝐿(𝛼𝛼)

Para 𝐿𝐿(𝛼𝛼) = 𝐿𝐿′ (𝛼𝛼) = ⋯ = 𝐿𝐿(𝑘𝑘−1) (𝛼𝛼) = 0 ; 𝐿𝐿(𝑘𝑘) (𝛼𝛼) ≠ 0

1 𝐴𝐴𝑥𝑥 𝑘𝑘
𝑦𝑦𝑝𝑝 = {𝐴𝐴𝑒𝑒 𝛼𝛼𝛼𝛼 } = (𝑘𝑘) 𝑒𝑒 𝛼𝛼𝛼𝛼 (2.3.2.7)
𝐿𝐿(𝐷𝐷) 𝐿𝐿 (𝛼𝛼)

(2) 𝑆𝑆𝑆𝑆 𝑅𝑅(𝑥𝑥) = 𝑃𝑃𝑘𝑘 (𝑥𝑥) es un polinomio de grado 𝑘𝑘

En este caso se aplicará una de las siguientes equivalencias

𝑘𝑘 𝑘𝑘
1 1
= � 𝐷𝐷ℓ = 1 + 𝐷𝐷 + 𝐷𝐷2 + ⋯ + 𝐷𝐷𝑘𝑘 + ⋯ , = �(−1)ℓ 𝐷𝐷ℓ = 1 − 𝐷𝐷 + 𝐷𝐷2 + ⋯ − 𝐷𝐷𝑘𝑘 + ⋯
1 − 𝐷𝐷 1 + 𝐷𝐷
ℓ=0 ℓ=0

1
O alternativamente se realizará la división 𝐿𝐿(𝐷𝐷)
obteniendo la solución particular para 𝐿𝐿(𝐷𝐷){𝑦𝑦} = 𝑃𝑃𝑘𝑘 (𝑥𝑥) de la
forma
1
𝑦𝑦𝑝𝑝 = {𝑃𝑃 (𝑥𝑥)} = [𝑏𝑏0 + 𝑏𝑏1𝐷𝐷 + 𝑏𝑏2 𝐷𝐷2 + ⋯ + 𝑏𝑏𝑘𝑘 𝐷𝐷𝑘𝑘 ]{𝑃𝑃𝑘𝑘 (𝑥𝑥)} (2.3.2.8)
𝐿𝐿(𝐷𝐷) 𝑘𝑘
Donde,
𝐷𝐷𝑘𝑘 𝑥𝑥 𝑛𝑛 = 0 ; 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑞𝑞𝑞𝑞𝑞𝑞 𝑘𝑘 ≥ 𝑛𝑛

18
ING. EDGAR H. MEDINA TAPIA MAT 1207

(3) 𝑆𝑆𝑆𝑆 𝑅𝑅(𝑥𝑥) = 𝐴𝐴 𝑐𝑐𝑐𝑐𝑐𝑐 𝛽𝛽𝛽𝛽 + 𝐵𝐵 𝑠𝑠𝑠𝑠𝑠𝑠 𝛽𝛽𝛽𝛽

La solución particular en este caso se obtiene mediante las fórmulas

Para 𝐿𝐿(−𝛽𝛽2 ) ≠ 0

1 1
𝑦𝑦𝑝𝑝 = {𝐴𝐴 𝑐𝑐𝑐𝑐𝑐𝑐 𝛽𝛽𝛽𝛽 + 𝐵𝐵 𝑠𝑠𝑠𝑠𝑠𝑠 𝛽𝛽𝛽𝛽} = (𝐴𝐴 𝑐𝑐𝑐𝑐𝑐𝑐 𝛽𝛽𝛽𝛽 + 𝐵𝐵 𝑠𝑠𝑠𝑠𝑠𝑠 𝛽𝛽𝛽𝛽) (2.3.2.9)
𝐿𝐿(𝐷𝐷2 ) 𝐿𝐿(−𝛽𝛽2 )

Para 𝐿𝐿(−𝛽𝛽2 ) = 𝐿𝐿′ (−𝛽𝛽2 ) = ⋯ = 𝐿𝐿(𝑘𝑘−1) (−𝛽𝛽2 ) = 0 ; 𝐿𝐿(𝑘𝑘) (−𝛽𝛽2 ) ≠ 0

1 𝑥𝑥 𝑘𝑘
𝑦𝑦𝑝𝑝 = {𝐴𝐴 𝑐𝑐𝑐𝑐𝑐𝑐 𝛽𝛽𝛽𝛽 + 𝐵𝐵 𝑠𝑠𝑠𝑠𝑠𝑠 𝛽𝛽𝛽𝛽} = (𝐴𝐴 𝑐𝑐𝑐𝑐𝑐𝑐 𝛽𝛽𝛽𝛽 + 𝐵𝐵 𝑠𝑠𝑠𝑠𝑠𝑠 𝛽𝛽𝛽𝛽) (2.3.2.10)
𝐿𝐿(𝐷𝐷2 ) 𝐿𝐿(𝑘𝑘) (−𝛽𝛽2 )

(4) Si 𝑅𝑅(𝑥𝑥) = 𝑃𝑃𝑘𝑘 (𝑥𝑥)𝑒𝑒 𝛼𝛼𝛼𝛼 [ 𝐴𝐴 𝑐𝑐𝑐𝑐𝑐𝑐 𝛽𝛽𝛽𝛽 + 𝐵𝐵 𝑠𝑠𝑠𝑠𝑠𝑠 𝛽𝛽𝛽𝛽]

La solución particular 𝑦𝑦𝑃𝑃 se obtiene usando las fórmulas de equivalencia de Euler para las funciones de
𝑐𝑐𝑐𝑐𝑐𝑐 𝛽𝛽𝛽𝛽 y 𝑠𝑠𝑠𝑠𝑠𝑠 𝛽𝛽𝛽𝛽; luego, se empleará la propiedad del desplazamiento del operador inverso y finalmente se
aplicarán los procedimientos del segundo caso para integrar polinomios

1
𝑦𝑦𝑝𝑝 = {𝑃𝑃 (𝑥𝑥)𝑒𝑒 𝛼𝛼𝛼𝛼 [𝐴𝐴 𝑐𝑐𝑐𝑐𝑐𝑐 𝛽𝛽𝛽𝛽 + 𝐵𝐵 𝑠𝑠𝑠𝑠𝑠𝑠 𝛽𝛽𝛽𝛽]}
𝐿𝐿(𝐷𝐷) 𝑘𝑘
1 1
= (𝐴𝐴 − 𝑖𝑖𝑖𝑖)𝑒𝑒 (𝛼𝛼+𝑖𝑖𝑖𝑖)𝑥𝑥 {𝑃𝑃 (𝑥𝑥)}
2 𝐿𝐿(𝐷𝐷 + 𝛼𝛼 + 𝑖𝑖𝑖𝑖) 𝑘𝑘 (2.3.2.10)
1 1
+ (𝐴𝐴 + 𝑖𝑖𝑖𝑖)𝑒𝑒 (𝛼𝛼−𝑖𝑖𝑖𝑖)𝑥𝑥 {𝑃𝑃 (𝑥𝑥)}
2 𝐿𝐿(𝐷𝐷 + 𝛼𝛼 − 𝑖𝑖𝑖𝑖) 𝑘𝑘

2.3.3.- El Método de Variación de Parámetros

Sea la ecuación diferencial lineal de orden n

𝐿𝐿(𝐷𝐷)[𝑦𝑦] = [𝑎𝑎𝑛𝑛 𝐷𝐷𝑛𝑛 + 𝑎𝑎𝑛𝑛−1 𝐷𝐷𝑛𝑛−1 + ⋯ + 𝑎𝑎1 𝐷𝐷 + 𝑎𝑎0 ]{𝑦𝑦} = 𝑅𝑅(𝑥𝑥) (2.3.3.1)

Cuya función o solución complementaria está dada por la ecuación

𝑦𝑦𝑐𝑐 = 𝑐𝑐1 𝑦𝑦1 + 𝑐𝑐2 𝑦𝑦2 + ⋯ + 𝑐𝑐𝑛𝑛 𝑦𝑦𝑛𝑛 (2.3.3.2)

El método de variación de parámetros plantea una solución particular a partir de la solución complementaria,
esto es

𝑦𝑦𝑝𝑝 = 𝑢𝑢1 𝑦𝑦1 + 𝑢𝑢2 𝑦𝑦2 + ⋯ + 𝑢𝑢𝑛𝑛 𝑦𝑦𝑛𝑛 (2.3.3.3)

En este caso, las 𝑢𝑢𝑘𝑘 = 𝑢𝑢𝑘𝑘 (𝑥𝑥) son funciones que han sustituido a las constantes de la solución complementaria. Las
funciones 𝑢𝑢𝑘𝑘 (𝑥𝑥) se determinan a través del siguiente procedimiento

1º Se genera el siguiente sistema de ecuaciones

𝑢𝑢′1 𝑦𝑦1 + 𝑢𝑢′2 𝑦𝑦2 + ⋯ + 𝑢𝑢′𝑛𝑛 𝑦𝑦𝑛𝑛 = 0

𝑢𝑢′1 𝑦𝑦′1 + 𝑢𝑢′2 𝑦𝑦′2 + ⋯ + 𝑢𝑢′𝑛𝑛 𝑦𝑦′𝑛𝑛 = 0 (2.3.3.4)
⋮
(𝑛𝑛−1) (𝑛𝑛−1) (𝑛𝑛−1)
𝑢𝑢′1 𝑦𝑦1 + 𝑢𝑢′2 𝑦𝑦2 + ⋯ + 𝑢𝑢′𝑛𝑛 𝑦𝑦𝑛𝑛 = 𝑅𝑅(𝑥𝑥)

19
ING. EDGAR H. MEDINA TAPIA MAT 1207

2º Se resuelve este sistema para obtener las funciones 𝑢𝑢′𝑘𝑘 (𝑥𝑥); en este caso, si se emplea el método de las
determinantes, se tienen
𝑦𝑦1 𝑦𝑦2 ⋯ 𝑦𝑦𝑘𝑘 ⋯ 𝑦𝑦𝑛𝑛 𝑦𝑦1 𝑦𝑦2 ⋯ 0 ⋯ 𝑦𝑦𝑛𝑛
𝑦𝑦′1 𝑦𝑦′2 ⋯ 𝑦𝑦′𝑘𝑘 ⋯ 𝑦𝑦′𝑛𝑛 𝑦𝑦′1 𝑦𝑦′2 ⋯ 0 ⋯ 𝑦𝑦′𝑛𝑛
𝑊𝑊 = � ⋮ ⋮ ⋱ ⋮ ⋱ ⋮ � 𝑊𝑊𝑘𝑘 = � ⋮ ⋮ ⋱ ⋮ ⋱ ⋮ � (2.3.3.5)
(𝑛𝑛−1) (𝑛𝑛−1) ⋯ (𝑛𝑛−1) ⋯ (𝑛𝑛−1) (𝑛𝑛−1) (𝑛𝑛−1) ⋯
𝑦𝑦1 𝑦𝑦2 𝑦𝑦𝑘𝑘 𝑦𝑦𝑛𝑛 𝑦𝑦1 𝑦𝑦2 𝑅𝑅(𝑥𝑥)⋯𝑦𝑦𝑛𝑛(𝑛𝑛−1)

Por lo tanto, las funciones 𝑢𝑢′𝑘𝑘 (𝑥𝑥) se obtienen de

𝑊𝑊𝑘𝑘
𝑢𝑢′𝑘𝑘 = ; 𝑘𝑘 = 1,2, ⋯ 𝑛𝑛 (2.3.3.6)
𝑊𝑊

Luego, las funciones 𝑢𝑢𝑘𝑘 (𝑥𝑥) son iguales a

𝑊𝑊𝑘𝑘
𝑢𝑢𝑘𝑘 = ∫ 𝑊𝑊
𝑑𝑑𝑑𝑑 ; 𝑘𝑘 = 1,2, ⋯ 𝑛𝑛 (2.3.3.7)

Finalmente, la solución particular toma la forma

𝑊𝑊1 𝑊𝑊2 𝑊𝑊𝑛𝑛

𝑦𝑦𝑝𝑝 = 𝑦𝑦1 � 𝑑𝑑𝑑𝑑 + 𝑦𝑦2 � 𝑑𝑑𝑑𝑑 + ⋯ + 𝑦𝑦𝑛𝑛 � 𝑑𝑑𝑑𝑑
𝑊𝑊 𝑊𝑊 𝑊𝑊 (2.3.3.8)

20