TECNÓLOGICO NACIONAL DE MÉXICO

CAMPÚS TECNÓLOGICO DE CERRO AZUL

ASIGNATURA:

ECUACIONES DIFERENCIALES

DOCENTE:

ING. ARRIETA OSORIO MARIANA

TRABAJO:

TEMA 4 SISTEMAS DE ECUACIONES

DIFERENCIALES LINEALES: INVESTIGACIÓN

NO. DE CONTROL:

20500471

ALUMNO:

JIMENEZ ROSALES JORGE ALBERTO

CARRERA:

INGENIERÍA EN SISTEMAS COMPUTACIONALES

SEMESTRE:

PERIODO:

ENERO-JUNIO

CERRO AZUL, VER. A 31 DE MAYO DE 2022

 ÍNDICE
INTRODUCCIÓN………………………………………………………………………...…..4

4. SISTEMAS DE ECUACIONES DIFERENCIALES LINEALES………………………..5

4.1 TEORÍA PRELIMINAR………………………………………………………..…5

4.1.1 SISTEMAS DE ECUACIONES DIFERENCIALES LINEALES…….6

4.1.2 SISTEMAS DE ECUACIONES DIFERENCIALES LINEALES

HOMOGÉNEOS.……………………………………………………………13

4.1.3 SOLUCIÓN GENERAL Y SOLUCIÓN PARTICULAR DE

SISTEMAS DE ECUACIONES DIFERENCIALES LINEALES…………22

4.2 MÉTODOS DE SOLUCIÓN PARA SISTEMAS DE ECUACIONES

DIFERENCIALES LINEALES……………………………………………………...24

4.3 MÉTODO DE LOS OPERADORES…………………………………………..25

4.4 UTILIZANDO LA TRANSFORMADA DE LAPLACE……………………......30

4.5 APLICACIONES……………………………………………………………......32

CONCLUSION………………………………………………………………………………34

FUENTES DE INFORMACION……………………………………………………………35
 INTRODUCCIÓN

Las ecuaciones diferenciales conforman, por una sección, una de las gigantes
subramas del Estudio matemático; con relevantes contactos con otras ramas de las
Matemáticas, como la Geometría diferencial, la Teoría de variable compleja, la Mejora
y el Cálculo de variaciones. Sin embargo, las ecuaciones diferenciales son un
instrumento omnipresente en Física e Ingeniería a partir de que Galileo y Newton
fundaron la Física actualizada. Actualmente además poseen aplicaciones importantes
en Química, Biología y Ciencias sociales. Las ecuaciones lineales son relevantes (en
Matemáticas, Física e Ingeniería), ya que, o bien corresponden con la naturaleza de
los inconvenientes, o bien conforman una primera aproximación a modelos no lineales.
En los últimos 100 años fueron desarrollándose lentamente modelos no lineales de
ecuaciones diferenciales ordinarias, apoyándose primero en la exploración cualitativo
y luego además en los pc y los programas informáticos de cálculo científico. Sin
embargo, los modelos lineales siguen siendo primordiales:.
 4. SISTEMAS DE ECUACIONES
DIFERENCIALES LINEALES

4.1 TEORÍA PRELIMINAR

Esta unidad se dedica al estudio de sistemas de
Ecuaciones Diferenciales de Primer Orden que son
casos especiales de sistemas que tienen la forma
normal. Es decir, las ecuaciones que se vieron en el
tema 1 que fue Ecuaciones Diferenciales Ordinarias
de Primer Orden se aplicaran en este nuevo tema,
pero ahora en base a un sistema de Ecuaciones
Diferenciales que como bien se sabe es la resolución de un conjunto de ecuaciones y
esas ecuaciones tendrán el mismo valor de resultados.

𝑑𝑥1𝑑𝑡= 𝑔1(𝑡,𝑥1,𝑥2,𝑥𝑛)

𝑑𝑥2𝑑𝑡= 𝑔2(𝑡,𝑥1,𝑥2,𝑥𝑛)

𝑑𝑥𝑛𝑑𝑡= 𝑔𝑛 𝑡,1,𝑥2,𝑥𝑛 .

Un sistema tal como el ejemplo anterior de n ecuaciones diferenciales de primer orden

se llama sistema de primer orden. Una ecuación diferencial de primer orden es una
ecuación diferencial ordinaria donde intervienen derivadas de primer orden respecto a
una variable independiente. Estas ecuaciones, junto con su condicional inicial, se
pueden encontrar expresadas en la forma explícita:

𝐝𝐲
= 𝐟 (𝐱, 𝐲)
{𝐝𝐱
𝐲(𝐱 𝟎 ) = 𝐲𝟎

O como su forma implícita:

𝐝𝐲
𝐟 (𝐱, 𝐲, ) = 𝟎 𝐜𝐨𝐧 𝐲(𝐱 𝟎 ) = 𝐲𝟎
𝐝𝐱
4.1.1 SISTEMAS DE ECUACIONES
DIFERENCIALES LINEALES

Un sistema de ecuaciones diferenciales es un conjunto de una o más ecuaciones en

las que aparecen una o más funciones incógnitas, pero
todas ellas dependiendo de una sola variable
independiente. El orden de un sistema de ecuaciones
diferenciales es el orden de la derivada de mayor orden
que aparece en el sistema.

La forma general de un sistema de dos ecuaciones y de

primer orden es:

𝐱′ = 𝐟(𝐭, 𝐱, 𝐲)
𝐲 ′ = 𝐠(𝐭, 𝐱, 𝐲)

Donde f y g son funciones de las tres variables 𝒙, 𝒚 y 𝒕 (variable independiente del

sistema).
Una solución del sistema en el intervalo (𝒂, 𝒃) es un par de funciones 𝒙(𝒕), 𝒚(𝒕) que
satisfacen:
𝐱’(𝐭) = 𝐟(𝐭, 𝐱(𝐭), 𝐲(𝐭))
𝐲’(𝐭) = 𝐠(𝐭, 𝐱(𝐭), 𝐲(𝐭))
Idénticamente para todo 𝒕 ∈ (𝒂, 𝒃).

EJEMPLO: Probar que las funciones 𝒙(𝒕) = 𝒆𝒕 , 𝒚(𝒕) = 𝒆−𝒕 son soluciones del
sistema en toda la recta real.

𝐱’ = 𝐱 𝟐 𝐲
𝐲’ = −𝐱𝐲 𝟐
En efecto, por una parte:
𝐱′(𝐭) = 𝐞𝐭 𝐲 𝐱(𝐭)𝟐 𝐲(𝐭) = 𝐞𝟐 𝐞−𝐭 = 𝐞𝐭 . Así pues 𝒙(𝒕) = 𝒙(𝒕)𝟐 𝒚(𝒕)
Y se satisface la primera ecuación idénticamente. Además

𝒚′(𝒕) = −𝒆 − 𝒕 𝒚 − 𝒙(𝒕)𝒚(𝒕)𝟐 = −𝒆𝒕 𝒆 − 𝟐𝒕 = −𝒆−𝒕 .

Entonces:
𝐲′(𝐭) = −𝐱(𝐭)𝐲(𝐭)𝟐

Con lo que se satisface la segunda ecuación.

En consecuencia, las funciones

x(𝐭) = 𝐞𝐭 , 𝐲(𝐭) = 𝐞−𝐭
Son soluciones del sistema.

Los sistemas no tienen por qué ser de dos 2 ecuaciones. En general, un sistema de n
ecuaciones diferenciales con n funciones incógnitas 𝐱𝟏(𝐭), . . . , 𝐱 𝐧 (𝐭) en la variable
independiente t, tendrá la siguiente forma:

𝐱’𝟏 = 𝐟𝟏(𝐭, 𝐱 𝟏 , 𝐱 𝟐 , . . . , 𝐱 𝐧 )
𝐱’𝟐 = 𝐟𝟐(𝐭, 𝐱 𝟏 , 𝐱 𝟐 , . . . , 𝐱 𝐧 )
.
∶
𝐱’𝐧 = 𝐟𝐧(𝐭, 𝐱 𝟏 , 𝐱 𝟐 , . . . , 𝐱 𝐧 )

Donde 𝐟𝟏, 𝐟𝟐, . . . , 𝐟𝐧 son funciones de las 𝐧 + 𝟏 variables 𝐭, 𝐱 𝟏 , 𝐱 𝟐 , . . . , 𝐱 𝐧 . Siempre

exigiremos que el número de ecuaciones y de incógnitas sea el mismo. Y a este
número común le llamaremos la dimensión del sistema.

Si cada una de las funciones 𝐟𝟏, . . . , 𝐟𝐧 en el sistema son lineales, entonces el sistema
se dice que es lineal.
La forma general de un sistema lineal de n ecuaciones de primer orden es:

𝐱’𝟏 = 𝐚𝟏𝟏(𝐭)𝐱 𝟏 + 𝐚𝟏𝟐(𝐭)𝐱 𝟐 + · · · 𝐚𝟏𝐧(𝐭)𝐱 𝐧 + 𝐛𝟏(𝐭)

𝐱’𝟐 = 𝐚𝟐𝟏(𝐭)𝐱 𝟏 + 𝐚𝟐𝟐(𝐭)𝐱 𝟐 + · · · 𝐚𝟐𝐧(𝐭)𝐱 𝐧 + 𝐛𝟐(𝐭)
.
:
𝐱’𝐧 = 𝐚𝐧𝟏(𝐭)𝐱 𝟏 + 𝐚𝐧𝟐(𝐭)𝐱 𝟐 + · · · 𝐚𝐧𝐧(𝐭)𝐱 𝐧 + 𝐛𝐧(𝐭)

Así los sistemas:

𝐱 𝟏 ’ = 𝟑𝐱 𝟏 − 𝟓𝐱 𝟐
𝐱 𝟐 ’ = −𝟐𝐱 𝟏

𝐱’ = −𝟐𝐭𝐲
𝐲’ = 𝟑𝐞𝐭𝐱 + 𝟒 𝐜𝐨𝐬 𝐭

𝟓
𝐮𝟏 ’ = −𝐜𝐨𝐬(𝐭)𝐮𝟏 − 𝐮𝟐
𝐭
𝐮𝟐 ’ = 𝐮𝟏 ’ − 𝐬𝐞𝐧(𝐭)𝐮𝟐 + √𝐭

Son lineales, mientras que los sistemas:

𝐱’ = 𝟑𝐱𝐲 − 𝟓𝐲
𝐲’ = −𝟐𝐱

𝐱 𝟏 ’ = −𝟐𝐱 𝟐
𝐱 𝟐 ’ = 𝟑𝐱 𝟐 𝟏 + 𝟒

𝐮𝟏 ’ = −𝐜𝐨𝐬(𝐭𝐮𝟏)
𝐮𝟐 ’ = 𝐮𝟏 − 𝐬𝐢𝐧(𝐭)

No lo son.
Si cada una de las funciones 𝐛𝟏(𝐭), 𝐛𝟐(𝐭), . . . , 𝐛𝐧(𝐭) son idénticamente cero, entonces
el sistema se dice que es homogéneo y en caso contrario, no homogéneo.
Los sistemas lineales son los más simples entre todos los sistemas de primer orden,
pero incluso éstos son muy difíciles de resolver analíticamente. Suficientemente
difíciles son los sistemas lineales con coeficientes constantes; es decir, aquellos en
los que las funciones 𝐚𝐢𝐣(𝐭) son funciones constantes. Éstos son los únicos que
podremos resolver analíticamente.

Lo mismo que sucede con los sistemas de ecuaciones algebraicas lineales, la notación
matricial es útil a fin de simplificar la exposición y destacar las propiedades de estos
sistemas.
En este punto si no se está suficientemente familiarizado con la teoría básica de
matrices, conviene estudiar o repasar sobre este tema para que no haya ningún
problema de entendimiento sobre este subtema y como hemos hecho en la sección
anterior para sistemas en general, llamaremos 𝐱(𝐭) al vector (función vectorial):

𝐱 𝟏 (𝐭)
𝐱 (𝐭)
𝐱(𝐭) = ( 𝟐 )
⁞
𝐱 𝐧 (𝐭)
De forma que:

𝐱′𝟏 (𝐭)
𝐱′ (𝐭)
x’(𝐭) = ( 𝟐 )
⁞
𝐱′𝐧 (𝐭)
Es el vector derivado de 𝐱(𝐭). Y si ponemos

𝒂𝟏𝟏 (𝒕) 𝒂𝟏𝟐 (𝒕) …. 𝒂𝟏𝒏 (𝒕) 𝒃𝟏 (𝒕)

𝒂𝟐𝟏 (𝒕) 𝒂𝟐𝟐 (𝒕) …. 𝒂𝟐𝒏 (𝒕) 𝒃 (𝒕)
𝑨(𝒕) = ⁞ ⁞ …. ⁞ y 𝒃(𝒕) = ( 𝟐 )
⁞
𝒂𝒏𝟏 (𝒕) 𝒂𝒏𝟐 (𝒕) …. 𝒂𝒏𝒏 (𝒕) 𝒃𝒏 (𝒕)
( )
Entonces el sistema se puede escribir abreviadamente:

𝐱 ′ (𝐭) = 𝐀(𝐭)𝐱(𝐭) + 𝐛(𝐭)

Y si es homogéneo (El vector 𝐛(𝐭) = 𝟎) entonces:

𝐱 ′ (𝐭) = 𝐀(𝐭)𝐱(𝐭)

Por ejemplo, el sistema:

𝐱’𝟏 = 𝐱𝟏 + 𝟐𝐱𝟐 + 𝐜𝐨𝐬 𝐭
𝐱’𝟐 = 𝟐𝐱𝟏 + 𝐱𝟐 + 𝐭 𝟐

Se puede escribir en notación matricial como:

𝟏𝟐 𝐜𝐨𝐬 𝒕
𝒙′ = ( )𝒙 + ( 𝟐 )
𝟐 𝟏 𝒕

De forma que la matriz del sistema es:

𝟏𝟐
𝐀(𝐭) = ( )
𝟐 𝟏

(Constante, en este caso) y el vector de los términos independientes es:

𝐜𝐨𝐬 𝒕
𝒃(𝒕) = ( 𝟐 )
𝒕

Se trata, entonces, de un sistema no homogéneo, cuya parte homogénea es:

𝟏𝟐
𝐱′ = ( )𝐱
𝟐 𝟏
Se puede comprobar, mediante sustitución, que los dos siguientes conjuntos de
funciones son soluciones del sistema homogéneo:

𝐱 𝟏 = 𝐞−𝐭 , 𝐱 𝟐 (𝐭) = −𝐞−𝐭

Estos dos conjuntos de soluciones se pueden escribir de forma más compacta usando
la notación vectorial:

−𝐭 𝟑𝐭
𝐱(𝐭) = ( 𝐞 −𝐭 ) , 𝐲(𝐭) = (𝐞𝟑𝐭 )
−𝐞 𝐞

En efecto, sustituyendo directamente en el sistema tendríamos:

−𝒕 𝟏𝟐 −𝒕 −𝒕
𝒙′(𝒕) = (−𝒆−𝒕 ) y 𝑨𝒙(𝒕) = ( ) ( 𝒆 −𝒕 ) =(−𝒆−𝒕 )
𝒆 𝟐 𝟏 −𝒆 𝒆

Y lo mismo para 𝐲(𝐭). Diremos, simplemente, que los vectores 𝐱(𝐭) y 𝐲(𝐭) son
soluciones del sistema inicial.

El método se basa en la eliminación que se utiliza para la resolución de sistemas de

ecuaciones algebraicas. En el caso de sistemas de ecuaciones diferenciales lineales,
el método de eliminación reduce el sistema a una sola ecuación diferencial de orden
n con coeficientes constantes en términos de una de las variables.

Para aplicar el método es necesario expresar el sistema en términos del operador

diferencial D. Veamos algunos ejemplos de cómo expresar un sistema lineal en
términos de “D”.
EJEMPLO 2. Exprese en notación matricial los sistemas siguientes:

a) (𝐃𝟐 + 𝐃)[𝐱] + (𝐃 + 𝟏)[ 𝐲] = 𝟎

𝐃𝟑 [𝐱] + (𝐃𝟐 + 𝟏)[ 𝐲] = 𝟎

La parte izquierda del sistema puede escribirse como el producto matricial de la matriz
operacional por el vector solución:

(𝐃𝟐 + 𝐃) (𝐃 + 𝟏) 𝐱 𝟎
[ ] [ ] = [ ]
𝐃𝟑 (𝐃𝟐 + 𝟏) 𝐲 𝟎

Observe que la matriz del sistema contiene operadores diferenciales. Por esta razón
se le llama matriz operacional, ya que puede interpretarse como un operador que
opera sobre un vector cuyos componentes son funciones.

b) 𝑫[𝒙] + 𝒚 + 𝒛 = 𝟎
−𝒙 + 𝑫[𝒚] + 𝒛 = 𝟎 →
𝒙 + 𝒚 + 𝑫[𝒁] = 𝟎

𝑫 𝟏 𝟏 𝒙 𝟎
[ −𝟏 𝑫 𝟏] [𝒚] = [𝟎]
𝟏 𝟏 𝑫 𝒛 𝟎

Al escribir los sistemas en forma matricial en términos de D, notamos que la notación

matricial para el caso de sistemas escritos en la forma normal es distinta.
Esto se debe a que para ese caso particular existe notación especial que es la que se
estudió en el tema anterior.
Note que todos los operadores diferenciales que aparecen en las matrices de los
ejemplos anteriores son de coeficientes constantes.
Ahora que sabemos cómo expresar un sistema lineal con coeficientes constantes en
términos del operador diferencial D, y en forma matricial, comenzaremos a resolver
sistemas de ecuaciones diferenciales. Primero veremos el caso homogéneo.
4.1.2 SISTEMAS DE ECUACIONES
DIFERENCIALES LINEALES
HOMOGÉNEOS

Nos centramos, por simplicidad, en los sistemas lineales

de dimensión, aunque las ideas que vamos a desarrollar
son válidas para sistemas de cualquier dimensión.
Los resultados generales los enunciaremos para sistemas
de dimensión n.

Los sistemas lineales de dimensión dos son de la forma:

𝐗 ′𝟏 = 𝐚𝟏𝟏 (𝐭)𝐱 𝟏 + 𝐚𝟏𝟐 (𝐭)𝐱 𝟐 + 𝐛𝟏 (𝐭) 𝐱 ′ = 𝐀(𝐭)𝐱 + 𝐛(𝐭),

𝐗 ′𝟐 = 𝐚𝟐𝟏 (𝐭)𝐱 𝟏 + 𝐚𝟐𝟐 (𝐭)𝐱 𝟐 + 𝐛𝟐 (𝐭)

Siendo

𝐚𝟏𝟏 (𝐭) 𝐚𝟏𝟐 (𝐭) 𝐛𝟏 (𝐭)

𝐀(𝐭) = ( ) 𝐲 𝐛(𝐭) = ( )
𝐚𝟐𝟏 (𝐭) 𝐚𝟐𝟐 (𝐭) 𝐛𝟐 (𝐭)

A estos sistemas también se les llama sistemas planos (o planares, según autores).
La primera propiedad que observamos de los sistemas lineales es que las soluciones
de estos sistemas cumplen el principio de superposición. Es decir, si

𝐮 (𝐭) 𝐮 (𝐭)
𝐮(𝐭) = ( 𝟏 ) y 𝐮(𝐭) = ( 𝟏 )
𝐮𝟐 (𝐭) 𝐮𝟐 (𝐭)

Son soluciones del sistema anterior a y b son números enteros reales entonces
𝐚𝐮(𝐭) + 𝐛𝐯(𝐭) también es solución del sistema anterior. Hay una forma en
matemáticas de expresar esta idea:

Las soluciones del sistema lineal 𝐱’ = 𝐀(𝐭)𝐱 forman un espacio vectorial.

Se trata, en realidad, de un subespacio vectorial del espacio vectorial de las funciones

continuas, pero esto no tiene, aquí, ninguna importancia.
Veamos que para los sistemas planes, en efecto funciona el principio de
superposición:

Sea

𝐰 ′ (𝐭) = 𝐚𝐮(𝐭) + 𝐛𝐯(𝐭), y veamos que 𝐰 ′ (𝐭) = 𝐀(𝐭)𝐰(𝐭):

𝐰 ′ (𝐭) = 𝐚𝐮′ (𝐭) + 𝐛𝐯 ′ (𝐭) = 𝐚𝐀(𝐭)𝐮(𝐭) + 𝐛𝐀(𝐭)𝐮(𝐭) = 𝐀(𝐭(𝐚𝐮(𝐭) + 𝐛𝐯(𝐭)) = 𝐀(𝐭)𝐰(𝐭)

La función vectorial 𝐰(𝐭) = 𝐚𝐮(𝐭) + 𝐛𝐯(𝐭) se dice que es una combinación lineal de
las funciones 𝒖(𝒕) y 𝒗(𝒕). Así pues, el principio de superposición dice que cualquier
combinación lineal de soluciones de un sistema lineal homogéneo es también una
solución del sistema.

El principio de superposición es una propiedad de los sistemas lineales que no es

verdadera en general para sistemas no lineales. Por ejemplo, si consideramos el
sistema:

𝐱 𝟏′ = 𝐱 𝟏𝟐
{ }
𝐱 𝟐′ = 𝐱 𝟏 𝐱 𝟐

Entonces
𝟏
−
𝐮(𝐭) = ( 𝐭 )
𝟎
Es una solución del sistema por que
𝟏 𝟐 𝟏
− 𝐮 𝟏 −
𝐮′(𝐭) = ( 𝐭 𝟐 ) 𝐲 ( ) = ( 𝐭𝟐)
𝐮𝟏 𝐮𝟐
𝟎 𝟎
Con lo que

𝐮𝟐
𝐮′(𝐭) = ( 𝟏 )
𝐮𝟏 𝐮𝟐
EJEMPLO: Dado el sistema
𝟏 𝟐
𝐱′ = ( )𝐱
𝟐 𝟏
Comprobar que
−𝐭 𝟑𝐭
𝐱 𝟏 (𝐭) = ( 𝐞 −𝐭 ) 𝐲 𝐱 𝟐 (𝐭) = (𝐞𝟑𝐭 )
−𝐞 𝐞
Son soluciones del sistema y probar que cualquier otra solución es combinación lineal
de estas dos.

En realidad, ya hemos visto que x1 (t) y x2 (t) son soluciones. Veamos ahora que si
x(t) es una solución del sistema, entonces x(t) es una combinación lineal de x1 (t) y x2
(t). Es decir, vamos a ver que existen unos números c1 y c2 tales que x(t)=c1x1(t)+
c2x2(t) para todo t. Encontrar estos números c1 y c2 para un valor de t concreto es muy
fácil. Por ejemplo, si supiéramos el valor de x(t) en t = 0, digamos
𝟑
𝐱(𝟎) = ( ) ,
−𝟏
Podríamos encontrar enseguida valores para c1 y c2 de modo que
𝐱(𝟎) = 𝐜𝟏𝐱𝟏(𝟎) + 𝐜𝟐𝐱𝟐(𝟎). En efecto, como
𝟏 𝟏
𝐱 𝟏 (𝟎) = ( ) 𝐲 𝐱 𝟐 (𝟎) = ( )
−𝟏 𝟏
Tendríamos que
𝟑 𝟏 𝟏 𝐜 +𝐜
𝐱(𝟎) = 𝐜𝟏 𝐱 𝟏 (𝟎) + 𝐜𝟏 𝐱 𝟏 (𝟎) ⟺ ( ) = 𝐜𝟏 ( ) + 𝐜𝟐 ( ) = −𝐜𝟏 + 𝐜𝟐
−𝟏 −𝟏 𝟏 𝟏 𝟐
𝟑 = 𝐜𝟏 + 𝐜𝟐
⟺{
−𝟏 = −𝐜𝟏 + 𝐂𝟐

Sólo hay que resolver este sencillo sistema de dos ecuaciones con dos incógnitas.
Pero este sistema no siempre tiene solución, ello depende de que la matriz de los
coeficientes tenga determinante distinto de cero. Éste es el caso en este ejemplo:
𝟏 𝟏
𝐝𝐞𝐭 ( ) = 𝟐 distinto de 𝟎
−𝟏 𝟏
De modo que el sistema tiene solución y además es única. En concreto c1=2 y c2=1
Debemos observar dos cosas:
𝐱 𝟏𝟎
Si el valor de x en t = 0 fuera otro, digamos 𝐱(𝟎) = ( 𝟎 ) entonces el sistema
𝐱𝟐
lineal a resolver sería
𝐱 𝟏𝟎 = 𝐜𝟏 + 𝐜𝟐
{ 𝟎
𝐱 𝟐 = −𝐜𝟏 + 𝐜𝟐
Que también tendría solución porque la existencia o no de soluciones de este sistema
depende de los términos independientes si no de la matriz de los coeficientes.

La matriz de los coeficientes depende exclusivamente de las funciones solución

dadas: 𝐱 𝟏 (𝐭) 𝐲 𝐱 𝟐 (𝐭) .Es, en efecto, la matriz cuyas columnas son 𝐱 𝟏 (𝟎) 𝐲 𝐱 𝟐 (𝟎):
𝟏 𝟏
(𝐱 𝟏 (𝟎) 𝐱 𝟐 (𝟎)) = ( )
−𝟏 𝟏

Resumiendo, hemos visto hasta ahora que si la matriz (𝐱 𝟏 (𝟎) 𝐱 𝟐 (𝟎)) tiene
determinante distinto de cero o lo que es lo mismo, si los vectores 𝐱 𝟏 (𝟎) 𝐲 𝐱 𝟐 (𝟎) son
linealmente independientes, entonces se pueden encontrar números c1 y c2 de modo
que 𝐱(𝟎) = 𝐜𝟏 𝐱𝟏 (𝟎) + 𝐜𝟐𝐱𝟐(𝟎).
Pero lo que nosotros pretendemos es mucho más ambicioso: pretendemos encontrar
unos números c1 y c2 tales que 𝐱(𝐭) = 𝐜𝟏𝐱𝟏(𝐭) + 𝐜𝟐𝐱𝟐(𝐭) para todo t y no sólo para
t= 0. Claro que deben existir tales números, éstos deberían ser hallados al resolver el
sistema 𝐱(𝟎) = 𝐜𝟏𝐱𝟏(𝟎) + 𝐜𝟐𝐱𝟐(𝟎) entonces 𝐱(𝐭) = 𝐜𝟏𝐱𝟏(𝐭) + 𝐜𝟐𝐱𝟐(𝐭) para todo t.
Esto parece una tarea difícil pero no lo es tanto por que hay una propiedad
fundamental de x(t) que todavía no hemos utilizado:
Es solución del sistema 𝐱’ = 𝐀(𝐭)𝐱. Además es la solución de este sistema que
𝐱 𝟏𝟎
cumple la condición inicial 𝐱(𝟎) = .
𝐱 𝟐𝟎
Ahora bien, si ponemos 𝐲(𝐭) = 𝐜𝟏𝐱𝟏(𝐭) + + 𝐜𝟐𝐱𝟐(𝐭) con c1 y c2 los números que
solucionan el sistema, tenemos que, por el principio de superposición 𝑦(𝑡) es solución
del sistema 𝐱’ = 𝐀(𝐭)𝐱 y cumple que 𝐲(𝟎) = + 𝐜𝟏𝐱𝟏(𝟎) + 𝐜𝟐𝐱𝟐(𝟎) = 𝐱(𝟎).
Es decir, 𝐲(𝐭) es una solución del sistema, que cumple la misma condición inicial que
𝐱(𝐭). El teorema de unicidad nos dice que esto sólo es posible si 𝐱(𝐭) = 𝐲(𝐭). Como
𝐲(𝐭) = 𝐜𝟏𝐱𝟏(𝐭) + 𝐜𝟐𝐱𝟐(𝐭) concluimos que 𝐱(𝐭) = 𝐜𝟏𝐱𝟏(𝐭) + 𝐜𝟐𝐱𝟐(𝐭) para todo t.
Si 𝐱𝟏(𝐭), 𝐱𝟐(𝐭), … . , 𝐱𝐧(𝐭) son exactamente n soluciones del sistema lineal homogéneo
𝐱 ′ = 𝐀 (𝐭)𝐱 tales que cualquier otra solución del sistema, 𝐱(𝐭), se puede poner como
combinación lineal de ellas, se dice que forman un sistema fundamental de soluciones
del sistema.
En el ejemplo anterior, las soluciones
𝟏 𝟐
𝐱𝟏 = ( ) 𝐱.
𝟐 𝟏

A la vista de este ejemplo y todo su desarrollo podemos enunciar el siguiente

resultado:

PROPOSICIÓN: Si 𝐱 𝟏 (𝐭), 𝐱 𝟐 (𝐭), … , 𝐱 𝐧 (𝐭) son n soluciones de un sistema lineal

homogéneo 𝐱 ′ = 𝐀(𝐭)𝐱 definido en el intervalo (𝐚, 𝐛), entonces forman un
sistema fundamental de soluciones si y solo si los vectores
𝐱 𝟏 (𝐭 𝟎 ), 𝐱 𝟐 (𝐭 𝟎 ), … , 𝐱 𝐧 (𝐭 𝟎 ) son linealmente independientes para algún 𝐭 𝟎 ∈ (𝐚, 𝐛).
Es decir, si y sólo si 𝐝𝐞𝐭(𝐱 𝟏 (𝐭 𝟎 )𝐱 𝟐 (𝐭 𝟎 ) … 𝐱 𝐧 (𝐭 𝟎 )) ≠ 𝟎 para algún 𝐭 𝟎 ∈ (𝐚, 𝐛).

En realidad, en el ejemplo de más arriba ya se han dado las ideas básicas para obtener
una demostración formal de esta Proposición para el caso 𝐧 = 𝟐. El caso general es
similar.
Veremos más adelante que todo sistema lineal homogéneo de dimensión 𝐧 admite un
sistema fundamental de 𝐧 soluciones. Este hecho, se puede enunciar
matemáticamente de la siguiente forma

TEOREMA: El conjunto de soluciones de un sistema n-dimensional lineal

homogéneo de primer orden 𝐱 𝟎 = 𝐀(𝐭)𝐱 forma un espacio vectorial de
dimensión 𝐧.
Si {𝐱 𝟏 (𝐭), 𝐱 𝟐 (𝐭), … , 𝐱 𝐧 (𝐭)} es un tal sistema, a la matriz

𝐗(𝐭) = (𝐱 𝟏 (𝐭)𝐱 𝟐 (𝐭) … 𝐱 𝐧 (𝐭))

Cuyas columnas son los vectores solución, se le llama matriz fundamental de

soluciones.
Podemos observar que

𝐗 ′ (𝐭) = (𝐱 ′ 𝟏 (𝐭) 𝐱 ′ 𝟐 (𝐭) ⋯ 𝐱 ′ 𝐧 (𝐭)) = (𝐀(𝐭)𝐱 𝟏 (𝐭) 𝐀(𝐭)𝐱 𝟐 (𝐭) ⋯ 𝐀(𝐭)𝐱 𝐧 (𝐭)) = 𝐀(𝐭)𝐗(𝐭) ,

Donde hemos hecho uso de la siguiente propiedad del producto de matrices: si

expresamos la matriz B en función de sus columnas: 𝑩 = (𝒃𝟏 𝒃𝟐 ⋯ 𝒃𝒏 ), entonces las
columnas de la matriz producto 𝐀𝐁 son 𝐀𝐛𝟏 , 𝐀𝐛𝟐 , … , 𝐀𝐛𝐧 ; es decir,

𝐀𝐁 = (𝐀𝐛𝟏 𝐀𝐛𝟐 ⋯ 𝐀𝐛𝐧 ).

Como conclusión de todo este proceso tenemos el siguiente teorema

TEOREMA: Una matriz 𝐗(𝐭) = (𝐱 𝟏 (𝐭) 𝐱𝟐 (𝐭) ⋯ 𝐱 𝐧 (𝐭))es una matriz fundamental
de soluciones del sistema lineal homogéneo 𝐱 ′ = 𝐀(𝐭)𝐱 de dimensión 𝐧 en el
intervalo (𝐚, 𝐛) si y solo si cumple las dos siguientes propiedades:

i. 𝑿′ (𝒕) = 𝑨(𝒕)𝑿(𝒕)
ii. 𝐝𝐞𝐭 𝑿(𝒕𝟎 ) ≠ 𝟎 Para algún 𝒕𝟎 ∈ (𝒂, 𝒃).

La primera condición indica que los vectores 𝐱 𝟏 (𝐭), 𝐱 𝟐 (𝐭), … , 𝐱 𝐧 (𝐭) (las columnas de
𝐗(𝐭)) son solución del sistema 𝐱 ′ = 𝐀(𝐭)𝐱. Y la segunda significa que estos vectores
son linealmente independientes en (𝐚, 𝐛). Es decir, el conjunto {𝐱 𝟏 (𝐭), 𝐱 𝟐 (𝐭), … , 𝐱 𝟐 (𝐭)}
forma un sistema fundamental de soluciones del sistema 𝐱 ′ = 𝐀(𝐭)𝐱.
Este Teorema nos indica el camino que debemos seguir para obtener la solución
general del sistema homogéneo 𝐱 ′ = 𝐀(𝐭)𝐱:

1. Encontrar n soluciones del sistema 𝐱 𝟏 (𝐭), 𝐱 𝟐 (𝐭), … , 𝐱 𝟐 (𝐭). Equivalentemente,

encontrar una matriz de n soluciones.
2. Comprobar que son linealmente independientes; es decir, que si (𝐚, 𝐛) es el
intervalo en el que está definido el sistema entonces
𝐝𝐞𝐭(𝒙𝟏 (𝒕𝟎 )𝒙𝟐 (𝒕𝟎 ) … 𝒙𝒏 (𝒕𝟎 )) = 𝐝𝐞𝐭 𝑿(𝒕𝟎 ) ≠ 𝟎 para algún 𝐭 𝟎 ∈ (𝐚, 𝐛).
3. Escribir la solución general del sistema como combinación lineal de las
soluciones encontradas. Es decir, la solución general del sistema sería:

𝐱(𝐭) = 𝐜𝟏 𝐱 𝟏 (𝐭) + 𝐜𝟐 𝐱 𝟐 (𝐭) + ⋯ + 𝐜𝐧 𝐱 𝐧 (𝐭)

Siendo 𝒄𝟏 , 𝒄𝟐 , … , 𝒄𝒏 constantes arbitrarias.

Si observamos que

𝐜𝟏
𝐜𝟐
𝐜𝟏 𝐱 𝟏 (𝐭) + 𝐜𝟐 𝐜𝟐 (𝐭) + ⋯ + 𝐜𝐧 𝐱 𝐧 (𝐭) = (𝐱 𝟏 (𝐭) 𝐱 𝟐 (𝐭) … 𝐱 𝟐 (𝐭)) ( ⋮ ) = 𝐗(𝐭)𝐜
𝐜𝐧

Con

𝐜𝟏
𝐜𝟐
𝐜 = ( ⋮ ),
𝐜𝟐

La solución general del sistema se podría escribir, en notación matricial, de la siguiente

forma:

𝐱(𝐭) = 𝐗(𝐭)𝐜.
 𝟎 𝟏
EJEMPLO: Dado el sistema 𝐱 ′ = ( ) 𝐱, comprobar que las funciones
−𝟐 𝟐
vectoriales

𝒆𝒕 𝑪𝒐𝒔 𝒕 𝒆𝒕 𝑺𝒆𝒏 𝒕
𝒙𝟏 (𝒕) = ( 𝒕 ) Y 𝒙𝟐 (𝒕) = ( 𝒕 )
𝒆 (𝑪𝒐𝒔 𝒕 − 𝑺𝒆𝒏 𝒕) 𝒆 (𝑪𝒐𝒔 𝒕 + 𝑺𝒆𝒏 𝒕)

Nótese que si pusiéramos

𝐞𝐭 𝐜𝐨𝐬𝐭 𝐭 𝐞𝐭 𝐬𝐞𝐧 𝐭 𝐜𝟏
𝐱(𝐓) = 𝐗(𝐭)𝐱 = ( 𝐭 (𝐜𝐨𝐬 𝐭 (𝐬𝐞𝐧 ) (𝐜 )
𝐞 𝐭 − 𝐬𝐞𝐧 𝐓) 𝐞 𝐭 + 𝐜𝐨𝐬 𝐭) 𝟐

𝐜𝟏 𝐞𝐭 𝐜𝐨𝐬 𝐭 + 𝐜𝟐 𝐞𝐭 𝐬𝐞𝐧 𝐭
=( )
𝐜𝟏 𝐞𝐭 (𝐜𝐨𝐬𝐭 𝐭 − 𝐬𝐞𝐧𝐭 𝐭) + 𝐜𝟐 𝐞𝐭 (𝐜𝐨𝐬 𝐭 + 𝐬𝐞𝐧 𝐭)

Forman un sistema fundamental de soluciones, escribir la solución general del sistema

𝟐
y hallar la única solución del sistema que cumple la condición inicial 𝐱(𝟎) = ( ).
𝟑
𝟎 𝟏
En primer lugar, hay que observar que la matriz del sistema 𝐀(𝐭) = ( ) es
−𝟐 𝟐
continua en toda la recta real, de modo que el intervalo de definición del sistema es
(−∞ + ∞). A continuación, tenemos que comprobar que 𝐱 𝟏 (𝐭) y 𝐱 𝟐 (𝐭) son soluciones
del sistema:

𝐞𝐭 (𝐜𝐨𝐬 𝐭 − 𝐬𝐞𝐧 𝐭) 𝟎 𝟏 𝐞𝐭 𝐜𝐨𝐬 𝐭

𝐱 ′ 𝟏 (𝐭) = ( ) 𝐀(𝐭)𝐱𝟏 (𝐭) = ( )( 𝐭 )
−𝟐𝐞𝐭 𝐬𝐞𝐧 𝐭 −𝟐 𝟐 𝐞 (𝐜𝐨𝐬 𝐭 − 𝐬𝐞𝐧 𝐭)
𝐞𝐭 (𝐜𝐨𝐬 𝐭 − 𝐬𝐞𝐧 𝐭)
=( ).
−𝟐𝐞𝐭 𝐬𝐞𝐧 𝐭

𝐞𝐭 (𝐬𝐞𝐧 𝐭 − 𝐜𝐨𝐬 𝐭) 𝟎 𝟏 𝐞𝐭 𝐬𝐞𝐧 𝐭

𝐱 ′ 𝟐 (𝐭) = ( ) 𝐀(𝐭)𝐱𝟐 (𝐭) = ( )( 𝐭 )
𝟐𝟐𝐞𝐭 𝐜𝐨𝐬 𝐭 −𝟐 𝟐 𝐞 (𝐜𝐨𝐬 𝐭 + 𝐬𝐞𝐧 𝐭)
𝐞𝐭 (𝐬𝐞𝐧 𝐭 − 𝐜𝐨𝐬 𝐭)
=( ).
𝟐𝐞𝐭 𝐬𝐞𝐧 𝐭

Calculamos

𝒆𝒕 𝒄𝒐𝒔 𝒕 𝒆𝒕 𝒔𝒆𝒏 𝒕
𝐝𝐞𝐭(𝒙𝟏 (𝒕)𝒙𝟐 (𝒕)) = ( 𝒕 ) = 𝒆𝒕 (𝑺𝒆𝒏𝟐 𝒕 + 𝑪𝒐𝒔𝟐 𝒕) = 𝒆𝒕
𝒆 (𝒄𝒐𝒔 𝒕 − 𝑺𝒆𝒏 𝒕) 𝒆𝒕 (𝒔𝒆𝒏 𝒕 + 𝒄𝒐𝒔 𝒕)
Así pues 𝐝𝐞𝐭(𝒙𝟏 (𝟎)𝒙𝟐 (𝟎)) = 𝒆𝟎 = 𝟏 ≠ 𝟎. esto significa que 𝐱 𝟏 (𝐭) y 𝐱 𝟐 (𝐭) son
linealmente independientes y forman un sistema fundamental de soluciones.
Equivalentemente, la matriz
𝐞𝐭 𝐜𝐨𝐬 𝐭 𝐞𝐭 𝐬𝐞𝐧 𝐭
𝐗(𝐭) = (𝐱 𝟏 (𝐭) 𝐱 𝟐 (𝐭)) = ( 𝐭 )
𝐞 (𝐜𝐨𝐬 𝐭 − 𝐬𝐞𝐧 𝐭) 𝐞𝐭 (𝐬𝐞𝐧 𝐭 + 𝐜𝐨𝐬 𝐭)

Es una matriz fundamental de soluciones. La solución general del sistema será

𝒆𝒕 𝒄𝒐𝒔 𝒕 𝒆𝒕 𝒔𝒆𝒏 𝒕
𝒙(𝑻) = 𝒄𝟏 𝒙𝟏 (𝒕) + 𝒄𝟐 𝒙𝟐 (𝒕) = ( 𝒕 ) + 𝒄𝟐 ( 𝒕 )
𝒆 (𝒄𝒐𝒔 𝒕 − 𝑺𝒆𝒏 𝒕) 𝒆 (𝒄𝒐𝒔 𝒕 + 𝒔𝒆𝒏 𝒕)
𝒄𝟏 𝒆𝒕 𝒄𝒐𝒔 𝒕 + 𝒄𝟐 𝒆𝒕 𝒔𝒆𝒏 𝒕
=( ).
𝒄𝟏 𝒆𝒕 (𝒄𝒐𝒔 𝒕 − 𝒔𝒆𝒏 𝒕) + 𝒄𝟐 𝒆𝒕 (𝒄𝒐𝒔 𝒕 + 𝒔𝒆𝒏 𝒕)
Nótese que si pusiéramos
𝐞𝐭 𝐜𝐨𝐬 𝐭 𝐞𝐭 𝐬𝐞𝐧 𝐭 𝐜𝟏
𝐱(𝐓) = 𝐗(𝐭)𝐜 = ( ) (𝐜 )
𝐞𝐭 (𝐜𝐨𝐬 𝐭 − 𝐬𝐞𝐧 𝐭) 𝐞 𝐭 (𝐬𝐞𝐧
𝐭 + 𝐜𝐨𝐬 𝐭) 𝟐

𝐜𝟏 𝐞𝐭 𝐜𝐨𝐬 𝐭 + 𝐜𝟐 𝐞𝐭 𝐬𝐞𝐧 𝐭
=( )
𝐜𝟏 𝐞𝐭 (𝐜𝐨𝐬 𝐭 − 𝐬𝐞𝐧 𝐭) + 𝐜𝟐 𝐞𝐭 (𝐜𝐨𝐬 𝐭 + 𝐬𝐞𝐧 𝐭)

Obtendríamos el mismo resultado. Finalmente, la solución que cumple la condición

𝟐
inicial 𝐱(𝐨) = ( ) será la que se obtiene imponiendo esta condición en la solución
𝟑
general:
𝟐 𝐜 ∙ 𝟏 + 𝐜𝟐 ∙ 𝟎
( )=( 𝟏 )
𝟑 𝐜𝟏 ∙ 𝟏 + 𝐜𝟐 ∙ 𝟏

Así 𝐜𝟏 = 𝟐 y 𝐜𝟐 = 𝟏 y la solución pedida será

𝐞𝐭 (𝟐 𝐜𝐨𝐬 𝐭 + 𝐬𝐞𝐧 𝐭)
𝐱(𝐭) = ( 𝐭 ).
𝐞 (𝟑 𝐜𝐨𝐬𝐭 𝐭 − 𝐬𝐞𝐧 𝐭)

𝐱 (𝐭)
O, escrito en función de las componentes 𝐱(𝐭) = ( 𝟏 ):
𝐱 𝟐 (𝐭)

𝐱 𝟏 (𝐭) = 𝐞𝐭 (𝟐 𝐜𝐨𝐬 𝐭 + 𝐬𝐞𝐧 𝐭)

𝐱 𝟐 (𝐭) = 𝐞𝐭 (𝟑 𝐬𝐨𝐬 𝐭 − 𝐬𝐞𝐧 𝐭).

En mi opinión para hallar la solución general de cualquier sistema lineal homogéneo

hay que diseñar un método para calcular n soluciones que sean linealmente
independientes, o equivalentemente, una matriz fundamental de soluciones.
4.1.3 SOLUCIÓN GENERAL Y SOLUCIÓN
PARTICULAR DE SISTEMAS DE
ECUACIONES DIFERENCIALES LINEALES

En general podemos decir que la solución de un sistema de

ecuación diferencial es llamada solución general si los
valores de las constantes no se obtienen en la solución final.
La misma solución puede convertirse en una solución
particular cuando tenemos el valor de las constantes
determinadas. Esto se hace en el caso que el sistema de
entrada de la ecuación diferencial sea un problema de valor
inicial con las condiciones iniciales establecidas para la
determinación de los términos constantes.

El ejemplo siguiente aclarará el procedimiento para resolver un sistema de ecuaciones

diferenciales lineales.

Determina el conjunto de ecuaciones como 𝒙𝑻(𝒕) = [(𝒙𝟏(𝒕), 𝒙𝟐(𝒕)] para el sistema

de ecuaciones 𝒅𝒙/ 𝒅𝒕 = 𝑨 ∗ 𝒙 con las condiciones iniciales establecidas como
𝒙(𝟎) = 𝒙𝟎 = (𝒙𝟎𝟏, 𝒙𝟎𝟐). El valor de la matriz A está dado como:

Entonces, el vector propio de la matriz es dado de la forma:

La matriz tiene un solo vector propio ya que ambos valores propios son los mismos.
Por lo tanto, la solución general del problema se da como:
Por consiguiente, un vector propio generalizado puede ser calculado como:

𝐯 = 𝐯𝐩 + 𝐬𝟏 ∗ 𝐯𝐡𝟏 𝐯𝟏 𝐯𝟐 − 𝟏 𝟎 + 𝐬𝟏 ∗ 𝟏 𝟏

La solución particular de este problema sería 𝒗𝑻(𝒑) = (−𝟏, 𝟎) = 𝒗𝟐, el cual es el

valor propio generalizado de esta matriz, junto con los valores propios repetidos y vh
es la solución homogénea dando el vector propio vector 1.

Y la solución general del problema es:

4.2 MÉTODOS DE SOLUCIÓN PARA
SISTEMAS DE ECUACIONES
DIFERENCIALES LINEALES
Un sistema de diferenciales lineales puede resolver las ecuaciones. Al igual que
existen varias técnicas para resolver una ecuación diferencial lineal, también las hay
para un sistema de ecuaciones diferenciales lineales. Como el método de eliminación
de Gauss, método separable y reducible etc. Sea un sistema de ecuaciones
diferenciales lineales representado como,

Entonces, la representación de la matriz equivalente de este sistema de ecuaciones

diferenciales lineales será,
 4.3 MÉTODO DE LOS
OPERADORES
MÉTODO DE ELIMINACIÓN

El método de eliminación sistemática para resolver sistemas de ecuaciones

diferenciales con coeficientes constantes se basa
en el principio algebraico de eliminación de
variables. Se verá que la operación análoga de
multiplicar una ecuación algebraica por una
constante es operar en una EDO (ecuación
diferencial ordinaria) con cierta combinación de
derivadas

ELIMINACIÓN SISTEMÁTICA: La eliminación de una incógnita en un sistema de

ecuaciones diferenciales lineales se agiliza al escribir una vez más cada ecuación del
sistema en notación de operador diferencial. Recuerde que una sola ecuación lineal.

𝒂𝒏 𝒚(𝒏) + 𝒂𝒏 − 𝟏 𝒚(𝒏 − 𝟏) + . . . + 𝒂𝟏 𝒚′ + 𝒂𝟎 𝒚 = 𝒈 (𝒕 )

Donde las ai , i = 0,1,..., n son constantes. Puede escribirse de nuevo como:

𝒂𝒏 𝑫𝒏 + 𝒂𝒏 − 𝟏 𝑫𝒏 − 𝟏 + . . . + 𝒂𝟏 𝑫 + 𝒂𝟎 𝒚 = 𝒈 (𝒕 )

Si el operador diferencial de n-ésimo orden 𝒂𝒏 𝑫 + 𝒂𝒏 − 𝟏 𝑫 + . . . + 𝒂𝟏 𝑫 + 𝒂𝟎 se

factoriza en operadores diferenciales de orden menor, entonces los factores
conmutan.
EJEMPLO: Escriba otra vez el siguiente sistema en términos del operador

𝑫𝒙 ′′ + 𝟐 𝒙′ + 𝒚′′ = 𝒙 + 𝟑 𝒚 + 𝒔𝒆𝒏𝒕
𝒙′ + 𝒚 ′ = 𝟒 𝒙 + 𝟐 𝒚 + 𝒆 − 𝒕

SOLUCIÓN: Primero se reúnen en un lado los términos con variables dependientes y

segundo se agrupan las mismas variables y tercero se reescribe cada una de las
ecuaciones del sistema en notación de operador diferencial.

(𝑫). 𝒙′′ + 𝟐𝒙′ − 𝒙 + 𝒚′′ − 𝟑𝒚 = 𝒔𝒆𝒏𝒕

𝒙′ − 𝟒 𝒙 + 𝒚 ′ − 𝟐 𝒚 = 𝒆 − 𝒕

(𝑫𝟐 + 𝟐𝑫 − 𝟏) 𝒙 + (𝑫𝟐 − 𝟑) 𝒚 = 𝒔𝒆𝒏𝒕

(𝑫 − 𝟒) 𝒙 + ( 𝑫 − 𝟐) 𝒚 = 𝒆 − 𝒕

EJEMPLO: Escriba otra vez el siguiente sistema en términos del operador

𝑫 𝒙′′′ − 𝟐 𝒙′′ + 𝒚′′ − 𝟐 𝒚′ = 𝒙′ − 𝟒 𝒙 + 𝟑 𝒚 + 𝒄𝒐𝒔𝒕

𝒙′′ + 𝒙′ − 𝟐 𝒚′′ = 𝟑 𝒙 + 𝟒 𝒚′ + 𝟓 𝒚 + 𝒔𝒆𝒏𝟐𝒕

SOLUCIÓN:
𝒙′′′ − 𝟐 𝒙′′ − 𝒙′ + 𝟒 𝒙 + 𝒚′′ − 𝟐 𝒚′ − 𝟑 𝒚 = 𝒄𝒐𝒔𝒕
𝒙′′ + 𝒙′ − 𝟑 𝒙 − 𝟐 𝒚′′ − 𝟒 𝒚′ − 𝟓 𝒚 = 𝒔𝒆𝒏𝟐𝒕

(𝑫𝟑 − 𝟐 𝑫𝟐 − 𝑫 + 𝟒) 𝒙 + ( 𝑫𝟐 − 𝟐𝑫 − 𝟑) 𝒚 = 𝒄𝒐𝒔 𝒕
(𝑫𝟐 + 𝑫 − 𝟑) 𝒙 − (𝟐𝑫𝟐 + 𝟒 𝑫 + 𝟓) 𝒚 = 𝒔𝒆𝒏𝟐𝒕
PROCEDIMIENTO PARA EL MÉTODO DE SOLUCIÓN DE SISTEMAS DE EDO
LINEALES POR ELIMINACIÓN SISTEMÁTICA PASO

PASO 1: Se cambia el sistema de la notación Leibniz a la notación Prima

PASO 2: Se reúnen en un lado los términos con variables dependientes y se
agrupan las mismas variables.
PASO 3: Se reescribe el sistema en términos del operador D. (Sistema
transformado)
PASO 4: Se elimina la variable x; multiplicando la ecuación (I) por el coeficiente
de x de la ecuación (II) y, multiplicando la ecuación (II) por el coeficiente de x
de la ecuación (I), tratando de que al multiplicarse ambas ecuaciones queden
con coeficientes de signo contrario para poder eliminarlos realizando la suma.
PASO 5: La ecuación factorizada obtenida en el paso 4 se resuelve utilizando
la ecuación característica (método de las m) para obtener las raíces.
PASO 6: Observando el tipo de raíces obtenidas en el paso 5 se decide la
solución complementaria y(t) que se tendrá.
PASO 7: Al realizar de nuevo los pasos 4, paso 5 y paso 6, se elimina la variable
y y se obtiene x(t).
PASO 8: Sustituyendo las ecuaciones obtenidas en los pasos 6 y paso7 (x (t),
y (t)) en la ecuación (I) del sistema transformado obtenido en el paso 3 se
calculan c3 y c4 en términos de c1 y c2 respectivamente.
PASO 9: se encuentra la solución del sistema original en términos de c1 y c2.

EJEMPLO
1: Por el método de eliminación sistemática resuelva el siguiente sistema de
ecuaciones lineales de primer orden.
𝒅𝒙
= 𝟑𝒚
𝒅𝒕

𝒅𝒚
= 𝟐𝒙
𝒅𝒕
PASO 1: Se cambia el sistema de la notación Leibniz a la notación

Prima x′ = 3y
▪ y′ = 2x

PASO 2: Se reúnen en un lado los términos con variables dependientes y se

agrupan las mismas variables.
▪ x′ − 3y = 0
▪ 2x − y′ = 0

PASO 3: Se reescribe el sistema en términos del operador D. (Sistema

transformado)
Dx− 3y = 0 (I)
2x − Dy = 0 (II)

PASO 4: Se elimina la variable x; multiplicando la ecuación (I) por el coeficiente

de x de la ecuación (II) y, multiplicando la ecuación (II) por el coeficiente de x
de la ecuación (I), tratando de que al multiplicarse ambas ecuaciones queden
con coeficientes de signo contrario para poder eliminarlos realizando la suma.
(−2) Dx − 3 y = 0
2 Dx + 6 y = 0

-2Dx + 6y = 0
2Dx − D2y = 0
______________
−D2y + 6y = 0

D2y − 6y = 0 → ( D2 − 6)y = 0 Ecuación factorizada

PASO 5: La ecuación factorizada obtenida en el paso (4) se resuelve utilizando

la ecuación característica (método de las m) para obtener las raíces.
m2 − 6 = 0 → m 2 = 6 → m1 = D1 = 6 , m2 = D2 = − 6
 PASO 6: Observando el tipo de raíces obtenidas en el paso (5) se decide la
solución complementaria y (t) que se tendrá. (en este caso la solución
complementaria corresponde al de raíces reales y distintas)

𝒚(𝒕) = 𝒄𝟏 𝒆𝒎𝟏 𝒕 + 𝒄𝟐 𝒆𝒎𝟐 𝒕

𝒚(𝒕) = 𝒄𝟏 𝒆√𝟔𝒕 + 𝒄𝟐 𝒆−√𝟔𝒕 (A)

PASO 7: Al realizar de nuevo los pasos (4), (5) y (6), se elimina la variable y y
se obtiene x (t).
𝒙(𝒕) = 𝒄𝟑 𝒆√𝟔𝒕 + 𝒄𝟒 𝒆−√𝟔𝒕 (B)

PASO 8: Sustituyendo las ecuaciones (A) y (B) en la ecuación (I) del sistema
transformado se calculan c3 y c4 en términos de c1 y c2 respectivamente.

√𝟔𝒄𝟑 𝒆√𝟔𝒕 − √𝟔𝒄𝟒 𝒆−√𝟔𝒕 − 𝟑(𝒄𝟏 𝒆√𝟔𝒕 − 𝒄𝟐 −𝒆√𝟔𝒕 ) = 𝟎

(√𝟔𝒄𝟑 − 𝟑𝒄𝟏 )𝒆√𝟔𝒕 + (−√𝟔𝒄𝟒 + 𝟑𝒄𝟐 )𝒆−√𝟔𝒕 = 𝟎

(√𝟔𝒄𝟑 − 𝟑𝒄𝟏 ) = 𝟎 (−√𝟔𝒄𝟒 + 𝟑𝒄𝟐 ) = 𝟎

𝟑 𝟑 √𝟔 𝟑√𝟔 √𝟔 √𝟔
√𝟔𝒄𝟑 = 𝟑𝒄𝟏 → 𝒄𝟏 → 𝒄𝟏 = 𝒄𝟏 = 𝒄 𝒄𝟑 = 𝒄
√𝟔 √𝟔 √𝟔 𝟔 𝟐 𝟏 𝟐 𝟏

𝟑 𝟑 −√𝟔 −𝟑√𝟔 −√𝟔

−√𝟔𝒄𝟒 = 𝟑𝒄𝟐 → 𝒄𝟐 → 𝒄𝟐 = 𝒄𝟐 = 𝒄 𝒄𝟒
−√𝟔 −√𝟔 −√𝟔 𝟔 𝟐 𝟐
√𝟔
=− 𝒄
𝟐 𝟐
PASO 9: se encuentra la solución del sistema original en términos de c1 y c2.
√𝟔 √𝟔
𝒙(𝒕) = 𝒄𝟏 𝒆√𝟔𝒕 − 𝒄𝟐 𝒆−√𝟔𝒕
𝟐 𝟐
𝒚(𝒕) = 𝒄𝟏 𝒆√𝟔𝒕 − 𝒄𝟐 𝒆−√𝟔𝒕
 4.4 UTILIZANDO LA
TRANSFORMADA DE LAPLACE
El procedimiento de transformada de Laplace puede aplicarse cuando se tiene un
sistema de ecuaciones diferenciales
con coeficientes constante. Se
resuelve el sistema en el dominio de
Laplace como un sistema de
ecuaciones algebraicas lineales.
Las soluciones en el dominio de
Laplace son luego invertidas al
dominio del tiempo.
EJEMPLO: SOLUCIÓN DE UN SISTEMA DE ECUACIONES DIFERENCIALES
Dado el siguiente sistema de ecuaciones diferenciales con condiciones iniciales
conocidas encuentra la solución para 𝑦1 y 𝑦2 en el dominio del tiempo.
𝒅𝒚𝟏 𝒅𝒚𝟐
− − 𝟐𝒚𝟏 + 𝟐𝒚𝟐 = 𝟏 − 𝟐𝒕
𝒅𝒕 𝒅𝒕
𝒅𝟐 𝒚 𝒅𝒚𝟐
+ 𝟐 + 𝒚𝟏 = 𝟎
𝒅𝒕𝟐 𝒅𝒕
𝒚𝟏 (𝟎) = 𝒚´𝟏 (𝟎) = 𝒚𝟐 (𝟎) = 𝟎
SOLUCIÓN
Aplicando la transformada de Laplace a las dos ecuaciones anteriores queda:
𝟏 𝟐
𝒚𝟏 (𝒔 − 𝟐) − 𝒚𝟐 (𝒔 − 𝟐) = −
𝒔 𝒔𝟐
𝒚𝟏 (𝒔𝟐 + 𝟏) − 𝒚𝟐 𝟐𝒔 = 𝟎
Despejando 𝑦2 de la ecuación 𝒚𝟏 (𝒔𝟐 + 𝟏) − 𝒚𝟐 𝟐𝒔 = 𝟎 y sustituyéndola en 𝒚𝟏 (𝒔 − 𝟐) −
𝟏 𝟐
𝒚𝟐 (𝒔 − 𝟐) = 𝒔 − 𝒔𝟐 queda:

(𝒔𝟐 + 𝟏)𝒚𝟏
𝒚𝟐 =
𝟐𝒔
(𝒔𝟐 + 𝟏)𝒚𝟏 𝟏 𝟐
𝒚𝟏 (𝒔 − 𝟐) − ( ) (𝒔 − 𝟐) = − 𝟐
𝟐𝒔 𝒔 𝒔
𝒔𝟑 − 𝟑𝒔 − 𝟐 (𝒔 + 𝟐)(𝒔 − 𝟏)𝟐 𝒔 − 𝟐
𝒚𝟏 ( )= = 𝟐
𝟐𝒔 𝟐𝒔 𝒔
Resolviendo para 𝒚𝟏
𝟐
𝒚𝟏 =
𝒔(𝒔 + 𝟏)𝟐

𝟐 (𝒔𝟐 +𝟏)𝒚𝟏
Sustituyendo 𝒚𝟏 = 𝒔(𝒔+𝟏)𝟐 en la ecuación 𝒚𝟐 = queda de la siguiente manera:
𝟐𝒔

(𝒔𝟐 + 𝟏) 𝟐 𝒔𝟐 + 𝟏
𝒚𝟐 = =− 𝟐
𝟐𝒔 𝒔(𝒔 + 𝟏)𝟐 𝒔 (𝒔 + 𝟏)𝟐
𝟐 (𝒔𝟐 +𝟏) 𝟐 𝒔𝟐 +𝟏
La inversión de las ecuaciones 𝒚𝟏 = 𝒔(𝒔+𝟏)𝟐 y 𝒚𝟐 = = − 𝒔𝟐 (𝒔+𝟏)𝟐
𝟐𝒔 𝒔(𝒔+𝟏)𝟐

nos da la solución del sistema de ecuaciones

𝒚𝟏 (𝒕) = 𝟐 − 𝟐(𝒕 + 𝟏)𝒆−𝒕
𝒚𝟐 = −𝒕 + 𝟐 − 𝟐𝒕𝒆−𝒕 − 𝟐𝒆−𝒕
 4.5 APLICACIONES
Los sistemas de ecuaciones
diferenciales lineales encuentran sus
aplicaciones en varios problemas que
surgen en el sistema del mundo real.
Uno de estos problemas puede ser
como el que se muestra a continuación:

Problemas eléctricos: Muchos de los

circuitos eléctricos pueden ser
reducidos para solucionar un sistema
de ecuaciones diferenciales. Sea un circuito eléctrico dado como:

TRAYECTORIAS ORTOGONALES

Como primera aplicación de algunos de los conceptos y procedimientos relativos a

ecuaciones diferenciales de primer orden que se han introducido hasta el momento,
se va a considerar el problema de hallar las trayectorias ortogonales a una familia de
curvas dada. El problema tiene su origen en algunas aplicaciones en las que aparecen
dos familias de curvas (planas) que son mutuamente ortogonales; esto es, cada curva
de una de las familias es ortogonal en todos sus puntos a cada curva de la otra familia.
En tales casos, se dice que cada una de ellas es una familia de curvas ortogonales a
la otra.
De esta manera, obtener las trayectorias ortogonales a una familia de curvas de las
cuales se conoce su expresión analítica 𝒇(𝒙; 𝒚; 𝑪) = 𝟎, requiere los siguientes pasos:

Diferenciar la ecuación f(𝒙; 𝒚; 𝑪) = 𝟎.

Eliminar C del sistema formado por la ecuación y su derivada, para así obtener
la ecuación diferencial de la familia.
Obtener la ecuación diferencial de la familia de trayectorias ortogonales
aplicando el resultado precedente.
Integrar la ecuación resultante.

MODELOS DE POBLACION

El análisis del problema de tasar el crecimiento de una población permite ilustrar las
diferentes etapas que se presentan en el estudio de un problema de Matemática
Aplicada.

Problema físico: Se trata de obtener una ley que permita predecir el

crecimiento de una población 𝒑(𝒕) que varía con el tiempo.
Formulación matemática: La primera dificultad que se encuentra al intentar
formular matemáticamente el problema es que una población solo toma valores
enteros y, por lo tanto, considerada como función del tiempo, es una función
discontinua. El problema se solventa observando que, si se trata de una
población elevada, la variación en una unidad es comparada con el número de
individuos que constituye el valor de la población y, por consiguiente, se puede
considerar que esta vara continua e incluso diferenciablemente con el tiempo.
Modelo matemático (Primera aproximación): En primera aproximación
podría hacerse la hipótesis de que \el número de individuos que nacen y que
mueren por unidad de tiempo es proporcional a la población existente".

DESINTEGRACION RADIACTIVA

Experimentalmente se ha comprobado que el ritmo de desintegración de los

elementos radiactivos es proporcional a la cantidad de elemento presente. Así, si Q(t)
designa dicha cantidad en cada instante de tiempo, se tiene que la ley de
desintegración es:

𝒅𝑸 𝒅𝒕 = 𝑲𝑸 (𝑲 < 𝟎)
 CONCLUSIÓN
En esta indagación vimos los Sistemas de ecuaciones diferenciales lineales, el cual
nos ayuda a entender mejor la aplicación a inconvenientes que implican bastante más
de una variable que es dependiente de procesos simultáneos. Un sistema de
ecuaciones diferenciales lineales es un grupo de numerosas ecuaciones diferentes
con algunas incógnitas. Además, hay los sistemas de ecuaciones diferenciales
lineales homogéneas, dichos sistemas se otorgan una vez que las funcionalidades
conocidas de las ecuaciones del sistema no se hallan en ellas. Hay 2 resoluciones de
sistemas de ecuaciones diferenciales lineales, una de ellas es la solución general la
cual se aplica si los valores de las constantes no se obtienen en la solución final y la
otra solución es la especial que es una vez que poseemos el costo de las constantes
determinadas. Estas resoluciones se realizan en la situación que el sistema de ingreso
de la ecuación diferencial presente un costo inicial con las condiciones iniciales
establecidas para que se determinen los términos constantes. Hay diversos
procedimientos de solución para sistemas de ecuaciones diferenciales lineales, como
lo son el procedimiento de supresión de Gauss, el procedimiento separable y reducible
entre otros. En esta averiguación nos enfocamos en el procedimiento de los
operadores, que se basa en transformar una funcionalidad en otra funcionalidad. Otra
de las formas de resolver los sistemas de ecuaciones diferenciales lineales es usando
transformadas de Laplace, que en mi criterio es el más difícil para la solución de los
sistemas de ecuaciones diferenciales lineales debido a que para eso poseemos que
usar en algunas ocasiones la regla de Cramer. Debido a esta indagación conocimos
como tienen la posibilidad de usar los sistemas de ecuaciones diferenciales lineales,
así como además como tienen la posibilidad de resolver y en que tienen la posibilidad
de utilizar, tales como en inconvenientes mecánicos de acoplamiento de resortes y en
inconvenientes eléctricos entre otros.
 FUENTES DE INFORMACIÓN
Zill, Dennis & Cullen, Michael (2008). Matemáticas Avanzadas Para Ingeniería
I Ecuaciones Diferenciales. Tercera Edición. Ed. McGraw-Hill.

Rainville, Earl (2006). Ecuaciones Diferenciales Elementales. Segunda Edición.

Ed. Trillas.

Ramiro Garcia Torres. (2015). Sistemas de ecuaciones diferenciales lineales.

16 junio 2016, de Instituto Politécnico Nacional Sitio web:
http://itpn.mx/recursosisc/4semestre/ecuacioneslineales/Unidad%20IV.pdf

Tania Del Ángel Delgado. (2012). UNIDAD 4 ECUACIONES DIFERENCIALES

4.1 TEORIA PRELIMINAR. 16 junio 2016, de Blogspot Sitio web:
http://tanniadelgadogonzalez.blogspot.mx/2012/06/unidad-4-ecuaciones-
diferenciales-41.html