ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
SISTEMAS DE ECUACIONES DIFERENCIALES LINEALES
Hasta ahora hemos estudiado ecuaciones diferenciales que sólo contienen dos
variables, aquí consideraremos ecuaciones en las que intervienen más de dos
variables. Si una de las variables es independiente, las ecuaciones son Ecuaciones
diferenciales Ordinarias; si son independientes más variables, las ecuaciones se
llaman Ecuaciones Diferenciales entre Derivadas Parciales, como mencionamos en
el Módulo 1.
Para resolver un sistema de n-ecuaciones diferenciales ordinarias con n-variables
dependientes, consiste en obtener un conjunto en el que se puedan eliminar todas
menos una de las variables dependientes, por ejemplo “x”. La ecuación resultante de
la eliminación se resuelve entonces para esta variable “x”. Cada una de las variables
dependientes se obtiene de forma análoga.
El método de eliminación sistemática para resolver sistemas de ecuaciones lineales
con coeficientes constantes está basado en el principio algebraico de eliminación de
variables. El análogo de multiplicar una ecuación algebraica por una constante es
operar en una EDO con alguna combinación de derivadas. El proceso de eliminación
se acelera al escribir nuevamente cada ecuación de un sistema empleando una
notación de operador diferencial.
1
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
Método de los Operadores
Un “operador” lo podemos definir como un objeto matemático que toma una función
y nos devuelve otra. Por ejemplo, la operación de derivar hace precisamente esto:
𝑑𝑓
Dada la función 𝑦 = 𝑓(𝑥) nos devuelve otra función 𝑦 ′ = 𝑑𝑥 que llamamos derivada
de la función original. Lo que hacemos a continuación es definir el operador asociado
a la derivación.
Definición:
Dada una ecuación Diferencial Ordinaria 𝐹(𝑥, 𝑦, 𝑦 ′ , … , 𝑦 (𝑛) ) = 0, definimos su
operador diferencial P(D) sustituyendo cada derivada 𝑦 (𝑛) , por su operador 𝐷𝑛 .
El operador “D” aplicado a la función “y”; 𝐷𝑦 devuelve la derivada de la función:
𝐷𝑦 = 𝑦′
Aplicando el operador “D”, repetidas veces obtenemos las derivadas sucesivas de la
función:
𝐷(𝐷𝑦) = 𝐷2 𝑦 = 𝑦′′
𝐷(𝐷2 𝑦) = 𝐷3 𝑦 = 𝑦′′′
𝐷𝑛 𝑦 = 𝑦 (𝑛)
2
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
Ejemplo 1
Escriba con el operador “D” la siguiente Ecuación Diferencial
𝑥´´ + 2𝑥´ + 𝑦´´ = 𝑥 + 3𝑦 + 𝑠𝑒𝑛𝑡
Solución
𝑥´´ + 2𝑥´ − 𝑥 − 𝑦´´ − 3𝑦 = 𝑠𝑒𝑛𝑡
(𝐷2 + 2𝐷 − 1)𝑥 − (𝐷2 + 3)𝑦 = 𝑠𝑒𝑛𝑡
Propiedades de los operadores diferenciales
✓ La diferenciación es lineal 𝐷(𝑓 + 𝑔) = (𝐷𝑓) + (𝐷𝑔) ; 𝐷(𝑎𝑓) = 𝑎(𝐷𝑓)
en donde f y g son funciones y a es una constante.
✓ Cualquier polinomio en D con funciones como coeficientes es también un
operador diferencial. También se pueden componer operadores diferenciales con
la regla 𝐷1 °𝐷2 (𝑓) = 𝐷1 (𝐷2 (𝑓))
✓ Esta última propiedad dota al conjunto de los operadores lineales, sobre un cierto
espacio de funciones reales, de estructura de espacio vectorial sobre y
de módulo izquierdo sobre el mismo conjunto de funciones. Eso último implica
a su vez que el conjunto de operadores constituyen un álgebra asociativa.
✓ Se requiere de algo de cuidado: primero, cualesquiera coeficientes de función en
el operador D2 deben ser diferenciables tantas veces como requiera la aplicación
de D1. Para obtener un anillo de dichos operadores se debe suponer que se
utilizan derivadas de todos los órdenes. Segundo, este anillo no debe
3
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
ser conmutativo: un operador gD no es el mismo en general que Dg. De hecho
se tiene por ejemplo la relación básica en mecánica cuántica: Dx − xD = 1.
Recordemos que la derivación es una operación lineal:
′
(𝑐1 𝑓(𝑥) + 𝑐2 𝑔(𝑥)) = 𝑐1 𝑓 ′ (𝑥) + 𝑐2 𝑔′(𝑥) por lo tanto el operador “D” también
será lineal.
Si queremos resolver un Sistema de Ecuaciones Diferenciales Ordinarias con
coeficientes constantes, utilizando el Método de los Operadores, escribimos las
Ecuaciones Diferenciales del sistema con la notación “D”, y note que hay una
sorprendente semejanza entre este procedimiento y el método de Resolución de
un Sistema de Ecuaciones Algebraicas. Esto se debe a que es posible tratar el
operador “D” como una variable.
Ejemplo 2
Resuelva el siguiente sistema por el método de los operadores.
𝑑𝑥
= 3𝑦 𝐷𝑥 − 3𝑦 = 0
𝑑𝑡
{𝑑𝑦 entonces {
= 2𝑥 −2𝑥 + 𝐷𝑦 = 0
𝑑𝑡
Reduciendo algebraicamente el sistema para la variable “y”
𝐷 Dx − 3𝑦 = 0
{
3 −2𝑥 + Dy = 0
𝐷2 𝑥 − 3𝐷𝑦 = 0
−6𝑥 + 3𝐷𝑦 = 0
4
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
Tenemos (𝐷2 − 6)𝑥 = 0.
𝐷2 𝑥 − 6𝑥 = 0
Que es una ecuación, lineal homogénea con coeficientes constantes y la
resolvemos:
Utilizando la auxiliar 𝑚2 − 6 = 0 → 𝑚1,2 = ±√6
Así 𝑥(𝑡) = 𝑐1 𝑒 −√6𝑡 + 𝑐2 𝑒 √6𝑡
Ahora conseguimos 𝑦(𝑡) reemplazando en cualquier ecuación del sistema,
de forma conveniente; usamos la primera ecuación del sistema, pero
debemos derivar 𝑥(𝑡)
Entonces; 𝑥′(𝑡) = −√6𝑐1 𝑒 −√6𝑡 + 𝑐2 √6𝑒 √6𝑡
Sustituyendo en Dx − 3𝑦 = 0
−√6𝑐1 𝑒 −√6𝑡 + 𝑐2 √6𝑒 √6𝑡 − 3𝑦 = 0
−√6𝑐1 𝑒 −√6𝑡 + 𝑐2 √6𝑒 √6𝑡
𝑦(𝑡) =
3
Por lo tanto
𝑥(𝑡) = 𝑐1 𝑒 −√6𝑡 + 𝑐2 𝑒 √6𝑡
{ −√6𝑐1 𝑒 −√6𝑡 +𝑐2 √6𝑒 √6𝑡 son las soluciones del sistema
𝑦(𝑡) = 3
5
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
Ejemplo 3
Resuelva el sistema de Ecuaciones Diferenciales
𝑑𝑥 𝑑𝑦
2 − 5𝑥 + = 𝑒𝑡
{ 𝑑𝑡 𝑑𝑡
𝑑𝑥 𝑑𝑦
−𝑥 + = 5𝑒 𝑡
𝑑𝑡 𝑑𝑡
Solución
2𝐷𝑥 − 5𝑥 + 𝐷𝑦 = 𝑒 𝑡 −2𝐷𝑥 + 5𝑥 − 𝐷𝑦 = −𝑒 𝑡
{ ~ {
𝐷𝑥 − 𝑥 + 𝐷𝑦 = 5𝑒 𝑡 𝐷𝑥 − 𝑥 + 𝐷𝑦 = 5𝑒 𝑡
−𝐷𝑥 + 4𝑥 = 4𝑒 𝑡
(𝐷 − 4)𝑥 = −4𝑒 𝑡
Utilizando la ecuación auxiliar en la ecuación diferencial lineal de primer
orden 𝑚 − 4 = 0 → 𝑚 = 4 Entonces 𝑥𝑐 (𝑡) = 𝑐1 𝑒 4𝑡
Por medio d ellos coeficientes indeterminados
Supongamos 𝑥𝑝 = 𝐴𝑒 𝑡
𝑥′𝑝 = 𝐴𝑒 𝑡
4
Reemplazando 𝐴𝑒 𝑡 − 4(𝐴𝑒 𝑡 ) = −4𝑒 𝑡 → 𝐴=3
4
De allí que 𝑥𝑝 = 3 𝑒 𝑡
4
Por lo tanto 𝑥(𝑡) = 𝑐1 𝑒 4𝑡 + 3 𝑒 𝑡
6
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
Ahora para encontrar 𝑦(𝑡), reemplazamos en la segunda ecuación del
sistema
4
Considerando la derivada de 𝑥(𝑡) = 𝑐1 𝑒 4𝑡 + 3 𝑒 𝑡 , que es
4
𝑥′(𝑡) = 4𝑐1 𝑒 4𝑡 + 3 𝑒 𝑡 ,
Ahora sustituyendo en 𝐷𝑥 − 𝑥 + 𝐷𝑦 = 5𝑒 𝑡
4 4
(4𝑐1 𝑒 4𝑡 + 3 𝑒 𝑡 ) − (𝑐1 𝑒 4𝑡 + 3 𝑒 𝑡 ) + 𝐷𝑦 = 5𝑒 𝑡
𝐷𝑦 = 5𝑒 𝑡 − 3𝑐1 𝑒 4𝑡
Para conseguir 𝑦(𝑡) integramos ∫ 𝐷𝑦 = ∫(5𝑒 𝑡 − 3𝑐1 𝑒 4𝑡 )𝑑𝑡
3
𝑦(𝑡) = 5𝑒 𝑡 − 4 𝑐1 𝑒 4𝑡 + 𝑐2
4
𝑥(𝑡) = 𝑐1 𝑒 4𝑡 + 3 𝑒 𝑡
Así la solución del sistema es { 3
𝑦(𝑡) = 5𝑒 𝑡 − 4 𝑐1 𝑒 4𝑡 + 𝑐2
7
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
EJERCICIOS PROPUESTOS
8
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
APLICACIONES DE LAS ECUACIONES DIFERENCIALES
LINEALES DE SEGUNDO ORDEN CON COEFICIENTES
CONSTANTES
En este módulo se presentarán dos tópicos en los cuales las ecuaciones diferenciales
ordinarias de segundo orden desempeñan un papel vital para su modelación, a saber,
vibraciones mecánicas y circuitos eléctricos. Un principio fundamental de la física
establece que los sistemas físicos tienden a estar en una posición de mínima energía
potencial denominada posición de equilibrio y si, por alguna razón, el sistema es
forzado a salir de ese equilibrio, entonces tenderá a regresar a él.
Ecuación Diferencial de las vibraciones de una masa en un resorte.
Las vibraciones aparecen cuando se aplica una pequeña fuerza a un sistema físico que
se encuentra inicialmente en un estado de equilibrio estable. Cuando esta fuerza
desaparece, el sistema tiende a regresar a su posición de equilibrio. Para entender el
proceso físico que ocurre, recordemos que un sistema físico está en una posición de
equilibrio estable cuando se encuentra en un mínimo de energía potencial.
Para que abandone esa posición es necesario proporcionarle energía mediante la
acción de una fuerza. Cuando se deja de aplicar la fuerza, el sistema ha adquirido
energía potencial que, al intentar retornar a la posición de equilibrio, se transforma en
energía cinética. Es decir, cuando pasa la posición de equilibrio tendrá energía
cinética y no se detendrá; continuará su movimiento transformando ahora, hasta que
desaparezca, su energía cinética en potencial. Esta transferencia entre energía cinética
9
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
y potencial se repetirá indefinidamente a menos que algún mecanismo permita la
disipación de energía mecánica.
Movimiento Libre No Amortiguado
Para iniciar el estudio de las vibraciones mecánicas, analicemos una situación
cotidiana y simple. Consideremos un cuerpo de masa m que está unido a una pared
por medio de un resorte de constante k (sistema masa-resorte) el cual se encuentra
sobre una mesa horizontal. Por simplicidad supongamos también que no existe
fricción entre el cuerpo y la mesa y que el sistema se encuentra inicialmente en
equilibrio. De repente, el resorte se comprime (o se elonga) una distancia pequeña 𝑥0 ,
medida desde la posición de equilibrio y se le aplica una velocidad 𝑣0 . Desde ese
momento, el resorte ejerce una fuerza sobre la masa que tiende a regresarla a su
posición de equilibrio inicial. En general, esta fuerza depende de la distancia
comprimida (o elongada) del resorte. Si la compresión (o elongación) es pequeña, se
puede suponer que la fuerza es directamente proporcional a dicha deformación y que
siempre apunta hacia la posición de equilibrio o en sentido contrario a la deformación.
Dicha suposición se conoce como ley de Hooke para resortes lineales. Es decir, la
fuerza 𝐹𝑅 que en todo momento ejerce el resorte sobre la masa está dada por
𝐹𝑅 = −𝑘𝑥 , donde x es la deformación y k > 0 es la constante del resorte. Por otra
parte, y de acuerdo con la segunda ley de Newton, la suma de todas la fuerzas que se
aplican a un cuerpo produce un cambio a su movimiento que se rige por la ecuación
𝒅𝟐 𝒙
𝑭 = 𝒎𝒂 = 𝒎 𝒅𝒕𝟐
10
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
Igualando estos dos resultados, se obtiene el PVI que modela el sistema masa-resorte:
𝒅𝟐 𝒙
𝒎 𝒅𝒕𝟐 = −𝒌𝒙 ; con las condiciones iniciales 𝑥(0) = 𝑥0 ; 𝑣(0) = 𝑣0 , o
equivalentemente:
𝒅𝟐 𝒙
𝒎 𝒅𝒕𝟐 + 𝒌𝒙 = 𝟎 ; con las condiciones iniciales 𝑥(0) = 𝑥0 ; 𝑥′(0) = 𝑣0 ,
El modelo encontrado es una ecuación diferencial de segundo orden con coeficientes
constantes. Para resolverla proponemos una ecuación diferencial de segundo orden
𝒅𝟐 𝒙
con coeficientes constantes. De allí que 𝒎 𝒅𝒕𝟐 + 𝒌𝒙 = 𝟎
𝒅𝟐 𝒙 𝒌
+ 𝒙=𝟎
𝒅𝒕𝟐 𝒎
𝑘
Sea 𝑚 = 𝑤 2 una constante
𝒅𝟐 𝒙
Entonces + 𝒘𝟐 𝒙 = 𝟎 es la ecuación a resolver
𝒅𝒕𝟐
Usando la ecuación auxiliar 𝑚2 + 𝑤 2 = 0 → 𝑚1,2 = ±𝑖
(tenga cuidado ya que esta m, es la considerada en el módulo 2, para las raíces
de la ecuación auxiliar, no se trata de la masa).
Al colocar m (masa) en una posición inicial 𝑥0 e imprimirle una velocidad
inicial 𝑣0 , m tiende a oscilar en torno a su posición de equilibrio. En la
siguiente figura se muestra un esquema del resorte vertical
11
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
Ejemplo 1
Una fuerza de 400 Newtons estira un resorte 2m. Una masa de 50 Kg se
sujeta al extremo del resorte e inicialmente se libera desde la posición de
equilibrio con velocidad ascendente de 10m/seg. Encuentre la ecuación de
movimiento.
Solución
400
Sea 𝐹𝑤 = 𝑘𝑥 → 𝑘= = 200
2
x(0)=0 , x’(0)= -10m/seg
𝑑2 𝑥
Luego 50 𝑑𝑡 2 + 200𝑥 = 0
𝑑2 𝑥
+ 4𝑥 = 0
𝑑𝑡 2
Usando la ecuación auxiliar 𝑚2 + 4 = 0 → 𝑚 = ±2𝑖
Así 𝑥(𝑡) = 𝑐1 𝐶𝑜𝑠2𝑡 + 𝑐2 𝑆𝑒𝑛2𝑡
Como x(0)=0 entonces 𝑐1 = 0
Como x’(0)= -10 entonces 𝑐2 = −5
De tal manera que 𝒙(𝒕) = −𝟓𝑺𝒆𝒏𝟐𝒕
12
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
Movimiento Libre Amortiguado
El concepto de movimiento armónico simple es un tanto irreal, dado que el
movimiento descrito asume que no hay fuerzas de retardo que actúen sobre la masa
puesta en movimiento. A menos que la masa esté suspendida en un vacío perfecto,
habrá al menos una fuerza resistente debida al medio ambiente. En el estudio de la
mecánica, las fuerzas amortiguadoras que actúan sobre un cuerpo se consideran
proporcionales a la potencia de la velocidad instantánea. Cuando no se aplican otras
fuerzas externas al sistema, de la segunda ley de Newton se deduce que
𝑑2𝑥 𝒅𝒙
𝑚 = −𝑘𝑥 − 𝜷
𝑑𝑡 2 𝒅𝒕
𝑑2𝑥 𝒅𝒙
𝑚 2
+𝜷 + 𝑘𝑥 = 0
𝑑𝑡 𝒅𝒕
Que también es una ecuación diferencial lineal de segundo orden con coeficientes
constantes.
13
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
Ejemplo 2
Una masa de 5 kg se une a un resorte de constante k=5N/m y a un amortiguador de
constante 𝛽 = 26𝑁 ∙ 𝑠/𝑚. La masa se suelta del punto 𝑥0 = −0.1m con velocidad
𝑣0 = 1.94 m/s; determine:
1. La posición, velocidad y aceleración de la masa en el tiempo 𝑡 ≥ 0.
2. El tiempo en que la masa cruza por la posición de equilibrio y la velocidad en ese
instante.
3. El tiempo en que la velocidad de la masa es 0 m/s, así como la posición y
aceleración en ese instante.
Solución
La ecuación con respecto a la posición de equilibrio, está dada por la solución del PVI
𝑑2 𝑥 𝑑𝑥
5 𝑑𝑡 2 + 26 𝑑𝑡 + 5𝑥 = 0 𝑥(0) = −0.1m ; 𝑥′(0) = 1.94 m/s
1
Utilizando la ecuación auxiliar 5𝑚2 + 26𝑚 + 5 = 0 → 𝑚1 = −5 , 𝑚2 = − 5
1
De allí que 𝑥(𝑡) = 𝑐1 𝑒 −5𝑡 + 𝑐2 𝑒 −5𝑡
Como 𝑥(0) = −0.1m entonces −0.1 = 𝑐1 + 𝑐2
1
1
Como 𝑥′(0) = 1.94 m/s reemplazamos en 𝑥′(𝑡) = −5𝑐1 𝑒 −5𝑡 − 5 𝑐2 𝑒 −5𝑡
1
Entonces 1.94 = −5𝑐1 − 5 𝑐2
−0.1 = 𝑐1 + 𝑐2
Resolviendo el sistema de ecuaciones obtenido { 1
1.94 = −5𝑐1 − 5 𝑐2
14
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
2 3
Tenemos que 𝑐1 = − 5 ; 𝑐2 = 10
1) Por lo tanto :
1
2 3
𝑥(𝑡) = − 5 𝑒 −5𝑡 + 10 𝑒 −5𝑡 es la posición de la masa en el tiempo 𝑡 ≥ 0
1
3
𝑥′(𝑡) = 2𝑒 −5𝑡 − 50 𝑒 −5𝑡 es la velocidad de la masa en el tiempo 𝑡 ≥ 0
1
3
𝑥 ′′ (𝑡) = −10𝑒 −5𝑡 + 𝑒 −5 𝑡 es la aceleración de la masa en el tiempo 𝑡 ≥ 0
250
2) La masa cruza por la posición de equilibrio cuando x=0
2 3 −1𝑡 2 −5𝑡 3 −1𝑡
0 = − 𝑒 −5𝑡 + 𝑒 5 →= 𝑒 = 𝑒 5
5 10 5 10
1
𝑒 −5𝑡 3 3
→ 1 =4 → 𝑒 (5−5)𝑡 = 4
− 𝑡
𝑒 5
24 3
Luego aplicando ln a ambos miembros - 𝑡 = 𝑙𝑛 (4) → 𝑡 ≈ 0.0599 𝑠
5
La masa pasa por la posición de equilibrio a los 0.0599 segundos
1
3
Además, 𝑥 ′ (𝑡) = 2𝑒 −5𝑡 − 50 𝑒 −5𝑡 →
3 −1(0.0599)
𝑣(0.0599𝑠) = 2𝑒 −5(0.0599) − 𝑒 5
50
𝑣 ≈ 1.4228 𝑚/𝑠 representa la velocidad en el tiempo t=0.0599 seg.
3. La máxima separación de la posición de equilibrio se obtiene cuando la velocidad
1
3
es cero; esto ocurre cuando 𝑣(𝑡) = 0 → 0 = 2𝑒 −5𝑡 − 50 𝑒 −5𝑡 → 𝑡 ≈ 0.7305𝑠
En este tiempo la posición está dada por
1
2 3
𝑥(0.7305𝑠) = − 5 𝑒 −5(0.7305𝑠) + 10 𝑒 −5(0.7305𝑠)
𝑥(0.7305𝑠) ≈ 0.2489𝑚
15
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
1
3
Y la aceleración por: 𝑥 ′′ (0.7305𝑠) = −10𝑒 −5(0.7305𝑠) + 250 𝑒 −5(0.7305𝑠)
𝑥 ′′ (0.7305𝑠) = 𝑎 ≈ −0.2489𝑚/𝑠 2
Movimiento Forzado
La fuerza de excitación desempeña un papel diferente al de las otras fuerzas internas
del sistema, pues a veces provoca una reducción de la velocidad y en otras provoca
un aumento. Es decir, la fuerza de excitación puede reducir o aumentar la energía
cinética del sistema. Cuando la fuerza de excitación sea distinta de cero, diremos que
el sistema masa-resorte-amortiguador está forzado.
Suponga que tomamos en consideración una fuerza externa f(t) que actúa sobre una
masa vibratoria en un resorte. La inclusión de f(t) en la formulación de la segunda ley
de Newton de la ecuación diferencial del movimiento forzado:
𝑑2𝑥 𝑑𝑥
𝑚 2
= −𝑘𝑥 − 𝛽 + 𝑓(𝑡)
𝑑𝑡 𝑑𝑡
𝑑2𝑥 𝑑𝑥
𝑚 2 +𝛽 + 𝑘𝑥 = 𝐹(𝑡)
𝑑𝑡 𝑑𝑡
Ejemplo 4
Una masa de 1 slug, se sujeta a un resorte cuya constante es de 5lbs/pie. En un inicio,
la masa se libera desde 1pie por debajo de la posición de equilibrio a velocidad
descendente de 5pies/s, y el movimiento subsiguiente se desarrolla en un medio que
ofrece una fuerza de amortiguamiento numéricamente igual al doble de la velocidad
instantánea. Encuentre la ecuación de movimiento si la masa es conducida por una
fuerza externa igual a 𝑓(𝑡) = 12𝐶𝑜𝑠2𝑡 + 3𝑆𝑒𝑛2𝑡
16
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
Solución
Usamos el modelo
𝑑2𝑥 𝑑𝑥
𝑚 2 +𝛽 + 𝑘𝑥 = 𝐹(𝑡)
𝑑𝑡 𝑑𝑡
Además; 𝑚𝑎𝑠𝑎 = 𝑚 = 1𝑠𝑙𝑢𝑔 , 𝛽 = 2 , 𝑘 = 5 , 𝑥(0) = 1 , 𝑥 ′ (0) = 5
𝑑2 𝑥 𝑑𝑥
+ 2 + 5𝑥 = 12𝐶𝑜𝑠2𝑡 + 3𝑆𝑒𝑛2𝑡
𝑑𝑡 2 𝑑𝑡
Resolvemos la ecuación homogénea asociada
𝑑2𝑥 𝑑𝑥
2
+2 + 5𝑥 = 0
𝑑𝑡 𝑑𝑡
Entonces 𝑚2 + 2𝑚 + 5 = 0 → 𝑚1,2 = −1 ± 2𝑖
Luego 𝑥(𝑡) = 𝑒 −𝑡 (𝑐1 𝐶𝑜𝑠2𝑡 + 𝑐2 𝑆𝑒𝑛2𝑡)
Por otro lado, utilizando coeficientes indeterminados
consideremos 𝑥𝑝 = 𝐴𝐶𝑜𝑠2𝑡 + 𝐵𝑆𝑒𝑛2𝑡 (note que no hay superposición de soluciones)
𝑥′𝑝 = −2𝐴𝑆𝑒𝑛2𝑡 + 2𝐵𝐶𝑜𝑠2𝑡
𝑥′′𝑝 = −4𝐴𝐶𝑜𝑠2𝑡 − 4𝐵𝑆𝑒𝑛2𝑡
(−4𝐴𝐶𝑜𝑠2𝑡 − 4𝐵𝑆𝑒𝑛2𝑡) + 2(−2𝐴𝑆𝑒𝑛2𝑡 + 2𝐵𝐶𝑜𝑠2𝑡) + 5(𝐴𝐶𝑜𝑠2𝑡 + 𝐵𝑆𝑒𝑛2𝑡) =
12𝐶𝑜𝑠2𝑡 + 3𝑆𝑒𝑛2𝑡
Igualando coeficientes encontramos que 𝐴 = 0 , 𝐵 = 3
De allí que 𝑥𝑝 = 3𝑆𝑒𝑛2𝑡
Por lo tanto 𝑥(𝑡) = 𝑒 −𝑡 (𝑐1 𝐶𝑜𝑠2𝑡 + 𝑐2 𝑆𝑒𝑛2𝑡) + 3𝑆𝑒𝑛2𝑡
Ahora como 𝑥(0) = 1 , 𝑑𝑒 𝑥(𝑡) = 𝑒 −𝑡 (𝑐1 𝐶𝑜𝑠2𝑡 + 𝑐2 𝑆𝑒𝑛2𝑡) + 3𝑆𝑒𝑛2𝑡 →
𝑐1 = 1
17
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
𝑥 ′ (𝑡) = −𝑒 −𝑡 (𝑐1 𝐶𝑜𝑠2𝑡 + 𝑐2 𝑆𝑒𝑛2𝑡) + 𝑒 −𝑡 (−2𝑐1 𝑆𝑒𝑛2𝑡 + 2𝑐2 𝐶𝑜𝑠2𝑡) + 6𝐶𝑜𝑠2𝑡
𝑝𝑎𝑟𝑎 𝑥 ′ (0) = 5 → 5 = −𝑐1 + 2𝑐2 + 6 → 𝑐2 = 0
Así 𝑥(𝑡) = 𝑒 −𝑡 𝐶𝑜𝑠2𝑡 + 3𝑆𝑒𝑛2𝑡
Sistemas Análogos
Circuitos Eléctricos en Serie
Si i(t) denota la corriente en el circuito eléctrico LRC en serie, entonces las caídas
de voltaje a través del inductor, del resistor, y del capacitor
son como se ilustra. Por la segunda ley de Kirchhoff, la suma de esos voltajes es
igual al voltaje E(t) aplicado al circuito, es decir:
𝑑𝑖 1
𝐿 + 𝑅𝑖 + 𝑞 = 𝐸(𝑡)
𝑑𝑡 𝐶
Pero la carga q(t) presente en el capacitor está relacionada con la corriente i(t) por
𝑑𝑞
𝑖= , y por lo tanto, se convierte en la ecuación diferencial lineal de segundo
𝑑𝑡
orden:
𝑑2𝑞 𝑑𝑞 1
𝐿 +𝑅 + 𝑞 = 𝐸(𝑡)
𝑑𝑡 𝑑𝑡 𝐶
18
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
Ejemplo 5
Encuentre una carga en el capacitor de un circuito LRC en serie cuando la
inductancia equivale a ¼ Hz, la resistencia 20Ώ, la capacitancia 1/300F y la
E(t)=0, además q(0)=4C y i(0)=0
1 𝑑2𝑞 𝑑𝑞 1
2
+ 20 + 𝑞=0
4 𝑑𝑡 𝑑𝑡 300
𝑑2 𝑞 𝑑𝑞
+ 80 + 1200𝑞 = 0
𝑑𝑡 2 𝑑𝑡
𝑞(𝑡) = 𝑐1 ℮−60𝑡 + 𝑐2 ℮−20𝑡
𝑑𝑎𝑑𝑎𝑠 𝑙𝑎 𝑐𝑜𝑛𝑑𝑖𝑐𝑖𝑜𝑛𝑒𝑠 𝑑𝑒 𝑣𝑎𝑙𝑜𝑟 𝑖𝑛𝑖𝑐𝑖𝑎𝑙
𝑐1 = −2
𝑐2 = 6
𝑞(𝑡) = −2℮−60𝑡 + 6℮−20𝑡
El siguiente es un ejemplo que ilustra, un sistema de Ecuaciones Diferenciales
Lineales con Coeficientes Constantes, mismo que estudiamos en la sección 4.6 del
módulo 4.
19
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
Ejemplo 6
Resuelva el siguiente circuito, utilizando ecuaciones diferenciales.
Determine 𝐼1 , 𝐼2 , 𝐼3 . Suponga que todas las corrientes iniciales se anulan
𝑑𝑖1 1
Considerando el modelo matemático 𝐿 + 𝑅𝑖1 + 𝐶 𝑞 = 𝐸(𝑡)
𝑑𝑡
𝑑𝑖1
Sea Malla 1: 0.005 + 10𝑖2 = 50
𝑑𝑡
𝑑𝑖3
Malla 2: 20𝑖3 + 0.01 − 10𝑖2 = 0
𝑑𝑡
𝑖1 = 𝑖2 + 𝑖3
𝑑𝑖1
Entonces + 2000𝑖2 = 10000
𝑑𝑡
𝑑𝑖3
2000𝑖3 + − 1000𝑖2 = 0
𝑑𝑡
𝑑𝑖1 𝑑𝑖2 𝑑𝑖3
Dado que 𝑖1 = 𝑖2 + 𝑖3 entonces = +
𝑑𝑡 𝑑𝑡 𝑑𝑡
𝑑𝑖2 𝑑𝑖3
+ + 2000𝑖2 = 10000
𝑑𝑡 𝑑𝑡
Así sustituyendo { 𝑑𝑖3
2000𝑖3 + − 1000𝑖2 = 0
𝑑𝑡
Utilizando Operadores “D”
𝐷𝑖2 + 𝐷𝑖3 + 2000𝑖2 = 10000
{
2000𝑖3 + 𝐷𝑖3 − 1000𝑖2 = 0 (∗)
𝐷𝑖3 + (𝐷 + 2000)𝑖2 = 10000
𝑎𝑠𝑜𝑐𝑖𝑎𝑛𝑑𝑜 {
(𝐷 + 2000)𝑖3 − 1000𝑖2 = 0
20
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
1000𝐷𝑖3 = 10 000 000
{
(𝐷 + 2000)2 𝑖3 = 0
(𝐷2 + 4000𝐷 + 20002 )𝑖3 + 1000𝐷𝑖3 = 10 000 000
(𝐷2 + 5000𝐷 + 20002 )𝑖3 = 10 000 000
Usando la ecuación auxiliar
𝑚2 + 5000𝑚 + 4 000 000 = 0 → 𝑚1 = −1000 , 𝑚2 = −4000
Por lo tanto 𝑖3𝑐 (𝑡) = 𝑐1 𝑒 −1000𝑡 + 𝑐2 𝑒 −4000𝑡
Ahora supongamos 𝑖3𝑝 (𝑡) = 𝐴 → 𝑖′3𝑝 (𝑡) = 0 → 𝑖′′3𝑝 (𝑡) = 0
5
Sustituyendo 0 + 5000(0) + 20002 𝐴 = 10 000 000 → 𝐴= 2
5
Entonces 𝑖3𝑝 (𝑡) = 2
𝟓
Por lo tanto 𝒊𝟑 (𝒕) = 𝒄𝟏 𝒆−𝟏𝟎𝟎𝟎𝒕 + 𝒄𝟐 𝒆−𝟒𝟎𝟎𝟎𝒕 + 𝟐
𝑑𝑖3
Por otro lado = −1000𝑐1 𝑒 −1000𝑡 − 4000𝑐2 𝑒 −4000𝑡
𝑑𝑡
Reemplazando en (*)
5
2000 (𝑐1 𝑒 −1000𝑡 + 𝑐2 𝑒 −4000𝑡 + 2 ) − 1000𝑐1 𝑒 −1000𝑡 −
4000𝑐2 𝑒 −4000𝑡 − 1000𝑖2 = 0
Reduciendo 𝒊𝟐 (𝒕) = 𝒄𝟏 𝒆−𝟏𝟎𝟎𝟎𝒕 − 𝟐𝒄𝟐 𝒆−𝟒𝟎𝟎𝟎𝒕 + 𝟓
𝟓
Como 𝒊𝟑 (𝟎) = 𝟎 → 𝒄𝟏 + 𝒄𝟐 = − 𝟐
5 10
Como 𝒊𝟐 (𝟎) = 𝟎 → 𝒄𝟏 − 𝟐𝒄𝟐 = −𝟓 → 𝒄𝟐 = 6 ; 𝑐1 = − 3
𝟏𝟎 𝟓 𝟓
Así, 𝒊𝟑 (𝒕) = − 𝒆−𝟏𝟎𝟎𝟎𝒕 + 𝟔 𝒆−𝟒𝟎𝟎𝟎𝒕 + 𝟐
𝟑
21
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
𝟏𝟎 𝟓
𝒊𝟐 (𝒕) = − 𝒆−𝟏𝟎𝟎𝟎𝒕 − 𝟑 𝒆−𝟒𝟎𝟎𝟎𝒕 + 𝟓
𝟑
Como 𝑖1 = 𝑖2 + 𝑖3 entonces
𝟐𝟎 −𝟏𝟎𝟎𝟎𝒕 𝟓 −𝟒𝟎𝟎𝟎𝒕 𝟏𝟓
𝒊𝟏 (𝒕) = − 𝒆 − 𝒆 +
𝟑 𝟔 𝟐
22
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
PROBLEMAS PROPUESTOS
23
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
24
ECUACIONES DIFERENCIALES ORDINARIAS
Profesora Aura Del C. Gómez
CIRCUITOS
25
