U11.Ecuaciones Diferenciales 2da Parte
U11.Ecuaciones Diferenciales 2da Parte
U11.Ecuaciones Diferenciales 2da Parte
Ecuaciones
Diferenciales
2°PARTE
Unidad 11
Autores: Elena Arlauskas, Julio Ferreiro y Gabriela Righetti 1
Análisis Matemático II Ecuaciones Diferenciales
Función homogénea
Necesitamos definir función homogénea para reconocer luego las ED homogéneas. Una función
de dos variables dada por f(x; y) es homogénea de grado m con respecto a sus dos variables si, para
toda constante 0 , se verifica que f (x; y) m f ( x; y ) .El exponente m de indica el grado de
homogeneidad de la función f.
f (x; y) x 4 4 x 2 y 2 4 x 4 4 x 2 y 2 4 f ( x; y)
3xy
b) g ( x; y ) es homogénea de grado cero porque se verifica que
x y2
2
3xy 3 xy
g (x; y ) 0 0 g ( x, y )
x y
2 2
x y
2 2
𝑦
Para resolver realizamos la sustitución z= para que la ecuación diferencial se transforma en una
𝑥
ecuación a variables separables.
Ejemplo: y ’ xy = x 2 + 2 y 2
x2 y2 x y
y ' y '
xy y x
Haciendo la sustitución y = z.x , y reemplazando y’ se tiene:
1 1 dz 1 dx dx
z 'x z z z ' x x z.dz . z.dz
z z dx z x x
1 2 y2
Resulta: z ln Cx . Reemplazando: 2
ln CX y 2 2 x 2 ln Cx S.G
2 2x
Se llaman así las ecuaciones del tipo ay´´ + by´+ cy = k(x) con b y c constantes y k(x) una función
continua en un cierto intervalo I del eje x. Para este tipo de EDO vale un teorema de existencia y
unicidad en todo el intervalo I de números reales, cuyo enunciado es:
Teorema de existencia y unicidad:
Si x0 es un punto del intervalo de números reales I, cualesquiera sean las condiciones iniciales
y(x0)= u0, y´(x0) = u1 (con u0 y u1 números reales arbitrarios), existe en todo I una y sólo una
solución de ay´´ + by´+ cy = k(x) que verifica tales condiciones iniciales.
Definiciones: Si k(x) = 0 para todo x del intervalo I, se dice que la ecuación es homogénea.
Si k(x) ≠ 0 se dice que la ecuación es no homogénea.
Son las ecuaciones del tipo ay´´+ by´+ cy = 0 con a, b y c constantes (a≠0), para las cuales vale el
siguiente
Teorema: Si y = y1(x) y y = y2(x) son soluciones de ay´´ + by´+ cy = 0 para todo x real, entonces
la combinación lineal y = C1 y1(x) + C2 y2(x) es también una solución, en todo el eje real, para
cualesquiera constantes C1 y C2.
Demostración: Como y1 (x) e y2 (x) son soluciones de ay´´ + by´+ cy = 0 se verifica que
ay1´´(x) + b y1´(x) + c y1 (x) = 0 y ay2´´(x) + b y2´ (x) + c y2 (x) = 0
Si se multiplica la primera ecuación por C1, la segunda por C2 y se suman, el resultado es:
a[C1 y1´´(x) + C2 y2´´(x)] + b [C1 y1´(x) + C2 y2´(x)] + c [C1 y1 (x) + C2 y2 (x)] = 0 por lo tanto, es
cierto que y = C1 y1 (x) + C2 y2 (x) es también solución.
Por tanto, para encontrar la SG basta obtener dos de estas funciones y1 e y2, como veremos más
adelante.
Nota: Es posible demostrar que, en la SG, están comprendidas todas las soluciones de dicha EDO,
o sea, que ay´´ +b y´+ c y = 0 no tiene soluciones singulares.
Solución general
Supongamos que de la EDO homogénea ay´´ +b y´+ c y = 0, se conocen dos soluciones y1(x)
e y2 (x). Vimos que la combinación lineal de ellas y=C1 y1(x)+C2 y2 (x) proporciona, para
cualquier valor de las constantes, una nueva solución. ¿Es esta combinación lineal la solución
general de la EDO?
Para que así sea, es necesario y suficiente que, fijadas las condiciones iniciales y(x0)= u0, y´(x0) = u1
(donde u0 y u1 son números reales arbitrarios), puedan determinarse los valores de las constantes C1
y C2 de manera de obtener una solución particular de la EDO que verifique las condiciones iniciales
impuestas. O sea, es necesario y suficiente que el sistema
u 0 C1 y1 ( x0 ) C 2 y 2 ( x0 )
admita una y sólo una solución C1 y C2.
u1 C1 y1 ( x0 ) C 2 y 2 ( x0 )
y1 ( x0 ) y 2 ( x0 )
Para que esto sea posible, el determinante de los coeficientes debe ser 0.
y1 ( x0 ) y 2 ( x0 )
Determinante wronskiano
y1 ( x) y 2 ( x)
Definición: Un determinante del tipo W ( x) se llama determinante wronskiano
y1 ( x) y 2 ( x)
Entonces, podemos decir que, para que y = C1 y1(x) + C2 y2 (x) exprese la SG de la EDO, es
necesario y suficiente que el wronskiano de las dos funciones y1(x), y2 (x) no se anule en el punto x0;
pero como x0 es un punto arbitrario cualquiera del eje real, resulta el siguiente:
Definición: Se dice que dos funciones y1(x), y2 (x) son linealmente dependientes si difieren por un
factor de proporcionalidad constante para todo valor de x de un intervalo real I.
y1 ( x)
Es decir, si resulta C para todo valor de x de un intervalo real I, las funciones son
y 2 ( x)
linealmente dependientes. En caso contrario, se dicen linealmente independientes.
Es posible demostrar que si el wronskiano de dos funciones derivables no se anula en algún punto
de un intervalo abierto (W(x0) 0), entonces las funciones son linealmente independientes en
dicho intervalo. W ( x0 ) 0 func. lin. indep. (1)
Si las funciones son linealmente independientes y, además, son soluciones de ay´´ +b y´+ c y =
0, entonces se cumple W(x0) 0 para todo x0 de I.
func. lin. indep. soluciones de la ED W ( x0 ) 0 (2)
Ejemplo 3:
sen x 1
a) sen x y -2 sen x no son linealmente independientes pues por lo cual resulta
2 sen x 2
sen x 2 sen x
W ( x) 0 para todo x . Aquí se cumple el contrarrecíproco de (1)
cos x 2 cos x
x e3x
b) e3x y x e3x son linealmente independientes pues x , que no es constante en todo un
e3x
intervalo I de números reales. Las funciones son además soluciones de y´´ - 6y + 9y = 0
(verifíquenlo)
e3x xe3 x
W ( x) e6x nunca se anula . Aquí se cumple (2)
3e 3 x e 3xe
3x 3x
La independencia lineal de las soluciones para las EDO del tipo ay´´ +b y´+ c y = 0 y el hecho de
que W(x) 0 para todos los reales, son condiciones equivalentes.
Queda el problema de hallar las dos soluciones y1(x), y2 (x) para expresar la SG de la EDO lineal
homogénea de 2º orden con coeficientes constantes.
Autores: Elena Arlauskas, Julio Ferreiro y Gabriela Righetti 5
Análisis Matemático II Ecuaciones Diferenciales
Actividad 2
Analizar si son linealmente independientes los siguientes pares de funciones:
b) y1 x e y 2 x e mx
mx
a) y1 e mx y 2 e px c) y1 e mx cos px y 2 e mx sin px
Rta: a) l.i. b) no son l.i. c) l.i.
am² + b m + c = 0
La última ecuación es muy importante para encontrar las soluciones de ay´´ + by´+ cy = 0, y se le
da el nombre especial de ecuación auxiliar o ecuación característica de la ED.
Las raíces de la ecuación auxiliar pueden encontrarse factorizando, o bien, aplicando la fórmula
para resolver ecuaciones de segundo grado:
- b± b 2 - 4ac
m1,2 =
2a
De modo que, si b² - 4ac es positivo, cero o negativo, la ecuación auxiliar tendrá respectivamente:
a) dos raíces reales y distintas m1 y m2
b) una raíz real doble m
c) dos raíces complejas conjugadas
e m1 x e m2 x
e m1 m2 x m2 m1 e m1 m2 x 0
1 1
para las cuales W ( x)
m1e m1 x m2 e m2 x m1 m2
Entonces, la SG es y C1e 1 C 2 e 2
mx mx
b) Si la ecuación auxiliar tiene una raíz m, doble, para hallar la SG deberemos hallar otra función
- b± b 2 - 4ca
Efectivamente: Sabemos que las raíces de am² + bm +c = 0 son m1,2 = .
2
Derivamos: y 2 e mx mx e mx y 2 m e mx m e mx m 2 x e mx 2 m e mx m 2 x e mx
Reemplazamos en la ED y reacomodamos factorizando:
a (2 m e mx + m² x emx) + b (m x emx + emx) + c x emx =
(am² + b m + c) x emx + (2am+b) emx =
= 0 x emx + 0 emx = 0 porque ambas expresiones entre paréntesis son nulas. Recordemos que
para que la ecuación auxiliar tenga raíz doble debe ser b² - 4ac = 0, de donde se obtiene
b
m= - o bien 2am + b = 0.
2a
Al satisfacerse la identidad, queda probado que la y2 = x e mx es también solución.
Observemos que las soluciones y1 e y2 son linealmente independientes. Entonces, la SG en el caso
c) Si la ecuación auxiliar de ay´´ + by´+ cy = 0 tiene raíces complejas s ± ti. Entonces la solución
general de la ecuación diferencial puede escribirse como y = e sx ( C1 cos (tx) + C2 sen (tx)).
ei t x ei t x
ye sx
y esx cos(t x)
2
Por lo tanto, y = e s t cos (t x) es una solución particular de ay´´ + by´+ cy = 0
Tomando C = - C = i/2 se llega a la solución particular y = e sx sen tx.
Es decir: Si la ecuación auxiliar de am² + bm + c = 0 tiene dos raíces complejas s ± ti , entonces la
solución general de ay´´ + by´ +cy = 0 es:
Actividad 3
1) Verificar que la función nula y = 0 es solución (trivial) de cualquier EDO homogénea.
2) Encontrar la SG o la SP, según se pida:
a) y’’ – y = 0
b) y”- 2y’+3y = 0
c) y”-2y’+ y = 0
d) 8 y’’ -2 y’ – 3y = 0 con y=f(x) que satisfaga las condiciones f (0) = 3 , f ’(0) = -1
Rtas: a) y=C1 ex + C2 e-x b) y = ex (C1 cos( 2 x)+ C2 sen( 2 x)) c) y=C1 ex+ C2 x ex d)
1 3
13 x 2 x
y e 2 e4
5 5
Nosotros nos limitaremos al caso en que el segundo miembro k(x) sea una función del tipo
(1) k(x) = Pn(x) e rx cos (px) o bien del tipo
(2) k(x) = Pn(x) e rx sen (px)
donde r y p son números reales y Pn (x)= bn xn +...+ b1 x+ b0 es un polinomio de grado n 0.
Pueden aparecer casos particulares, según se anulen, o no, los parámetros a y p o el cuando el grado
del polinomio Pn es cero:
Si p = 0, la k(x) de (1) se reduce a la forma
(3) Pn(x) e rx
Si p = 0, r = 0, la k(x) de (1) se reduce a la forma
(4) Pn(x)
Si p = 0, r 0 , Pn(x) = b0 , la k(x) de (1) se reduce a la forma
(5) b0 e rx
Si p 0 , r = 0 , Pn(x) = b0 , la k(x) de (1) se reduce a la forma
(6) b0 cos (px)
Situaciones similares ocurren con la k(x) = Pn(x) e rx sen (px)
Puede demostrarse (aunque no lo haremos) que, tanto en el caso que k(x) sea como en (1) o como en
(2), la ED no homogénea admite una solución particular del tipo:
(7) yp=xs(α+βx+….+δxn)erxcos(px)+( γ+λx+…+θxn)erxsen(px)
Donde el exponente s es el menor entero no negativo tal que ningún término de yp sea solución de
la ED homogénea asociada ay´´ + b y´ + c y = 0.
Las constantes , ,..., , , ,..., deberán determinarse exigiendo que (7) sea solución de la ED
ay´´ + b y´ + c y = k(x). Tales constantes se denominan coeficientes indeterminados.
Se hace notar que la expresión se simplifica bastante en los casos particulares mencionados antes:
en el caso (3) , la (7) se reduce a
(8) y p x s x ... x n e ax
(9) y p x s x ... x n
en el caso (5) la (7) se reduce a
(10) y p x s e ax
Observación: la forma (11) también será la yp propuesta para cuando k(x) = b0 sen (px)
Supongamos ahora que el 2º miembro de la ED no homogénea y´´ + by´ + c y = k(x) sea la suma
de dos funciones k1(x) y k2(x), o sea:
ay´´ + by´ + c y = k1(x)+ k2(x).
Para determinar una solución particular yp vale el siguiente principio de superposición de
soluciones:
Teorema:
Considerando separadamente las ecuaciones ay´´ + by´ + c y = k1(x) y ay´´ + by´ + c y = k2(x).
Si y1 es SP de ay´´ + by´ + c y = k1(x) y2 es SP de ay´´ + by´ + c y = k2(x), entonces y1 +
y2 es SP de ay´´ + by´ + c y= k1(x)+ k2(x)
Resumen:
El método se puede aplicar sin problemas cuando k(x) es:
1) Polinómica
2) Exponencial
3) Combinación de senos y cosenos
4) Producto de exponencial y polinomio
5) Combinación lineal de cualquiera de los casos anteriores.
n
1) k(x) = bi x i
i 0
2) k(x) = B. e x
n
4) k(x) = bi xi e x
i 0
es raíz simple n
x. Ai xi e x Los coeficientes a
i 0 determinar son los Ai.
es raíz doble n
x 2 Ai xi e x
i 0
5) Si Q(x) es una combinación lineal de los casos anteriores la yp que se propone es una
combinación lineal de las que se propondrían aisladamente.
imponer una solución del tipo y p x s x e x donde s es el menor entero no negativo tal que
Ahora se deben calcular las constantes y imponiendo que yp satisfaga la ED dada. Para ello,
escribimos las derivadas
y p e x x e x
y p 2 e x x e x
y luego sustituimos en la ED
Autores: Elena Arlauskas, Julio Ferreiro y Gabriela Righetti 12
Análisis Matemático II Ecuaciones Diferenciales
2 e x x e x e x x e x 2 x e x ( x 2) e x
Dividiendo por ex y agrupando términos queda
2 x 2 x 2
Por el principio de identidad de polinomios debe ser
2 2 2 1
1 5
Con lo cual
2 4
5 1
Y por lo tanto resulta y p x e x de donde la solución general está dada por:
4 2
5 1
y y h y p C1e x C2 e 2 x x e x
4 2
donde s es el menor entero no negativo tal que y p x s e x no coincida con ningún término de
Ahora se debe calcular la constante imponiendo que yp satisfaga la ED dada. Para esto, hallamos
las derivadas y p e x x e x x e x
y p e x x e x
Sustituyendo en la ED:
e x x e x x e x 2x e x e x
1
Dividiendo por e-x y aplicando el principio de identidad de polinomios, resulta , o sea que
3
1 1
y p x e x ; y la SG es y y h y p C1 e x C 2 e 2 x x e x
3 3
forma (9)
y p x s x ... x n
donde s es el menor entero no negativo tal que todos los términos de yp sean diferentes de los
términos de yh.. Como esta condición se cumple para s = 0, y como el polinomio de k(x) es de
primer grado, resulta y p = a + bx .
2 x x .
2 1 1 1
Por el principio de identidad de los polinomios debe ser ;
2 0 2 4
1 1
Con lo cual, la SG es y y h y p C1 e x C 2 e 2 x x
2 4
forma (8)
y p x s x ... x n e ax
donde s es el menor entero no negativo tal que todos los términos de yp sean diferentes de los
términos de yh.. Como esta condición se cumple para s = 2, y como el polinomio de k(x) es de
primer grado, resulta y p x 2 x e 3 x x 2 x 3 e 3 x .
Ahora se deben calcular las constantes y imponiendo que yp satisfaga la ED dada.
Para ello, escribimos las derivadas
y p 2x 3x 2 e 3x 3 x 2 x 3 e 3x
y p 2 6x e 3x 6 2x 3x 2 e 3x
9 x 2 x 3 e 3x .
6 1 1
Y, aplicando el principio de identidad de los polinomios, debe ser: 0;
2 0 6
1 3 3x
De lo que resulta que la SG es y y h y p C1 e 3 x C 2 x e 3 x x e
6
donde s es el menor entero no negativo tal que todos los términos de yp sean diferentes de los
términos de yh.. Como esta condición se cumple para s = 0, resulta
y p xcos x xsenx
forma y p x s cos(2 x) sen(2 x) donde s es el menor entero no negativo tal que todos
los términos de yp sean diferentes de los términos de yh.. Como esta condición se cumple para s = 1,
resulta y p x cos(2 x) sen(2 x)
donde s es el menor entero no negativo tal que todos los términos de yp sean diferentes de los
términos de yh.. Como esta condición se cumple para s = 0, resulta y p cos x sen x
Ahora se deben calcular las constantes , imponiendo que yp satisfaga la ED dada. Para ello,
calculamos las derivadas de yp, sustituimos en la ED y simplificamos la expresión que obtiene en el
primer miembro. Igualando los respectivos coeficientes de sen x y cos x de ambos miembros, se
2 1 1
tiene: ; 0
2 0 2
1
De lo que resulta que la SG es y y h y p C1 e x C 2 e x sen x
2
Ejemplo 15:
Dada la ED y´´ - y’ = 2 + e3x resulta yh = C1 ex + C2 e-x (comprobarlo).
Para hallar la yp se tendrá en cuenta el teorema anterior. Consideremos separadamente las
ecuaciones y´´ - y’ = 2 y´´ - y’ = e3x. A ellas les corresponden, respectivamente,
1 1
2 . Luego, la SG es y y h y p C1 e x C 2 e x 2 e 3 x
8 8
La fuerza recuperadora es conservativa, por lo que tiene asociado una energía potencial
1 2.
V ( x) kx
2
diferencial
donde los puntos indican derivación respecto del tiempo, y es la frecuencia natural de
vibración. La solución general a esta ecuación se puede escribir de la forma
donde .
La solución general a esta ecuación depende de la relación entre y .
Tenemos tres casos:
Oscilador infraamortiguado:
Este es el caso . La solución es de la forma (2)
donde
Oscilador crítico:
En este caso . La solución general es (3)
Oscilador sobreamortiguado:
Aquí se tiene y la nueva solución general es (4) donde
En este caso, obtenemos una solución secular, es decir, cuya amplitud aumenta en el tiempo hasta
hacerse muy grande. Físicamente, esta solución no tiene sentido, ya que tarde o temprano el
rozamiento, que siempre existe pero que en este caso despreciamos, entrará en juego impidiendo
que la amplitud de oscilación crezca indefinidamente.
homogénea (oscilador sin forzar), dada por las ecuaciones (2), (3) y (4) ,según la relación entre y
, y además
Como vemos, la solución particular, proporcional a , es la única que importa para tiempos
grandes, ya que todas las soluciones de la ecuación homogénea decaen exponencialmente. Así,
pues, tenemos un estado estacionario, correspondiente a oscilaciones de amplitud .
En este caso, la solución es válida para todas las frecuencias de la fuerza de forzamiento. Sin
embargo, vemos que la amplitud de respuesta es máxima para una frecuencia: la de resonancia,
Bibliografía utilizada