Tema 5 Grado en Ingeniería Mecánica

ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

CONOCIMIENTOS PREVIOS

Para poder seguir adecuadamente este tema, se requiere que el alumno repase y ponga al día sus
conocimientos en los siguientes contenidos:
• Cálculo diferencial de funciones de una variable.
• Cálculo integral de funciones de una variable.
• Derivación de la función compuesta e implícita.
• Representación de curvas planas.

OBJETIVOS ESPECÍFICOS

Los objetivos específicos de este tema son:

1. Poder definir los conceptos básicos del tema.
2. Saber si una función es solución de una e.d.o. de segundo orden . Dadas una familia
biparamétrica de curvas y una e.d.o. de segundo orden saber comprobar que la familia es la
solución general de la ecuación. Para una familia biparamétrica de curvas, saber encontrar
la ecuación cuya solución general es la familia dada.
3. Conocida la solución general de una e.d.o., poder encontrar la particular correspondiente a
unas condiciones iniciales o de frontera. Dado un problema de valor inicial, saber analizar si
existe solución y ésta es única.
4. Poder utilizar el método de reducción de orden para encontrar una segunda solución de una
ecuación lineal homogénea de segundo orden, conocida una primera solución.
5. Determinar si un conjunto de soluciones de una ecuación homogénea es un sistema
fundamental.
6. En ecuaciones lineales no homogéneas, conocer la construcción de la solución general a
partir de una particular y de la general de la homogénea asociada. Saber utilizar el principio
de superposición, para simplificar la resolución.
7. Saber encontrar la solución general de una e.d.o. lineal homogénea de segundo orden de
coeficientes constantes, manejando la relación entre las raíces de la ecuación característica
y el sistema fundamental de soluciones.
8. Discriminar cuándo se pueden utilizar los métodos de coeficientes indeterminados y
variación de constantes para la búsqueda de soluciones particulares y saber emplearlos.
9. Saber utilizar el cambio de variable en ecuaciones diferenciales de orden superior a uno.
10. Poder transformar un problema de valor inicial sobre una e.d.o. lineal en un sistema de
ecuaciones de primer orden para su posterior resolución analítica o numérica.
11. Saber obtener la solución analítica de un sistema de ecuaciones diferenciales lineales de
coeficientes constantes de primer orden.
2 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

ECUACIONES DIFERENCIALES LINEALES. GENERALIDADES

1 Introducción

Una aplicación clásica de las ecuaciones diferenciales se presenta en el estudio de un circuito

eléctrico que consiste en resistores, inductores y capacitores, al cual se aplica una fuerza
electromotriz. En este caso una aplicación de las leyes de Kirchoff conduce a la ecuación

d 2q dq 1
L 2
R  q  E (t )
dt dt C
donde L es la inductancia, R la resistencia, C la capacitancia, E (t ) la fuerza electromotriz, q(t )
la carga y t el tiempo. Esta ecuación diferencial constituye un ejemplo de las ecuaciones
diferenciales que se estudiarán en este tema: las ecuaciones lineales de segundo orden.

2 Definición ecuación diferencial lineal de orden n

Ecuación diferencial lineal de orden n.- Es una ecuación de la forma:

n 
y  a1 (x )y (n 1)    an 1 (x )y   an (x )y  b(x )

- Si b(x )  0 la ecuación diferencial lineal de orden n se llama homogénea

n 
y  a1 (x )y (n 1)    an 1 (x )y   an (x )y  0

- Si b(x )  0 , la ecuación lineal se llama no homogénea o completa.

Ecuación diferencial lineal de orden n de coeficientes constantes.- Es una ecuación de la

forma:
n 
y  a1y (n 1)    an 1y   an y  b(x ), ai  

Centraremos nuestro estudio en las ecuaciones lineales de segundo orden, cuya forma canónica es:
y   p(x )y   q(x )y  r (x ) (1)

Los resultados obtenidos para estas ecuaciones se generalizan fácilmente a ecuaciones lineales de
orden superior a 2.

3 Soluciones

Solución general.- Una familia biparamétrica de funciones (x , y,C 1, C2 )  0 es solución

general de la ecuación (1) si cada función de la familia satisface la ecuación y si dados unos
valores yo e y1 , existen C 1 y C 2 tales que la correspondiente función cumpla que
y xo   yo y que y ' x o   y1 .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
3

Solución particular.- Es cada una de las funciones que se obtiene de la solución general
al dar valores a los parámetros.

Solución singular.- Es una solución que no puede extraerse de la solución general.

Conocida la solución general, puede obtenerse la ecuación diferencial de la familia por derivación y
eliminación de parámetros.

C2
Ejemplo: Dada la solución general y  C 1x 2  , encontrar la EDO cuya solución es la familia de
x
curvas biparamétricas.
Para ello eliminamos los dos parámetros del sistema formado por

C2 
y  C 1x 2  (I ) 
x  (I )  x II  y  xy '  3C 1x 2
C2 
y '  2C 1x  (II )   C2
 (II )  x III
x2
2C 2    y ' xy ''  3
x3
y ''  2C 1  (III )
x3 

Sustituyendo en (I)

y  xy ' y ' x  x 2y ''

y   3y  y  x 2y ''  2y  x 2y ''  0
3 3

4 Problemas asociados

Problema de valores iniciales (PVI).- Es el formado por la ecuación diferencial junto con los
valores que debe de tomar la solución y su derivada en un mismo punto x  x o :
y(x 0 )  y 0 , y (x 0 )  y1

Problema de valores en la frontera.- Es el formado por la ecuación diferencial junto con

los valores fijados para la solución en dos puntos distintos x  x o y x  x 1 :
y(x 0 )  y 0 , y(x 1 )  y1

Problema mixto.- Es el formado por la ecuación diferencial junto con los valores fijados de
la solución y su derivada en dos puntos distintos x  x o y x  x 1 :
1y(x 0 )  2y (x 0 )  0, 1y(x 1 )  2y (x 1 )  0, i , i  

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

4 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

5 Existencia y unicidad de un PVI

Las condiciones de existencia y unicidad de la solución del problema de valor inicial vienen dadas por
el siguiente teorema:

TEOREMA DE EXISTENCIA Y UNICIDAD.- Si p(x ) , q(x ) y r (x ) son funciones continuas en

un intervalo a, b  , si x o es un punto cualquiera de a, b  y si y 0 e y1 son números
cualesquiera, entonces la ecuación y   p(x )y   q(x )y  r (x ) tiene solución única, y(x ) ,
en este intervalo, cumpliendo y(x 0 )  y 0 , y (x 0 )  y1 .

Ejemplo: Se considera el problema de valor inicial siguiente

x 2y '' y  x


y(1)  2 y ' 1  0

Vamos a ver que se puede asegurar que solo hay una curva solución de este problema de valor inicial.

y x
Solución La ecuación diferencial se puede escribir y '' 2
 2 donde
x x

1 1
p x   0 q x   r x  
x2 x
Estas funciones son continuas para todo x real distinto de cero. Por lo tanto, son continuas en un
entorno del punto x o  1 . Aplicando el teorema de existencia y unicidad se puede probar que existe
una única solución pasando por el punto (1,2) y con pendiente 0 en dicho punto.
Buscando la solución en Matlab y representando la familia de curvas solución general y la solución
que satisface el problema de valor inicial.
y   p(x )y   q(x )y  r (x )


y(x 0 )  y 0 , y (x 0 )  y1

syms y(x) C1 C2
ecuacion=diff(y,x,2)==-y/x^2+1;
sol(x,C1,C2)=dsolve(ecuacion)
hold on
for k1=0:2
for k2=-1:2;
fplot(sol(x,k1,k2),[0,2])
end
end
cond1=y(1)==2;der(x)=diff(y,1);cond2=der(1)==0;
sol1(x)=dsolve(ecuacion,cond1,cond2)
fplot(sol1(x),[0 2], 'Linewidth',4)
hold off

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
5

-1

-2

0 0.5 1 1.5 2

6 Solución general de una ecuación lineal de segundo orden

TEOREMA.- Si se denota por y p cualquier solución particular de la ecuación completa y por

yh la solución general de la homogénea asociada, entonces la solución general de la
completa es: y  y p  yh

Así pues, a la hora de resolver una ecuación diferencial lineal debemos resolver dos problemas:
• encontrar la ecuación general de la homogénea asociada y
• encontrar una particular de la completa.

ECUACIONES DIFERENCIALES LINEALES HOMOGÉNEAS

Consideramos en este apartado las ecuaciones diferenciales de la forma

y   p(x )y   q(x )y  0

7 Soluciones

Comenzaremos observando que cualquier combinación lineal de dos soluciones de la ecuación

homogénea, es también solución. En efecto, si y1 (x ) e y2 (x ) son soluciones de

y   p(x )y   q(x )y  0

la función y(x )  C 1y1 (x )  C 2y2 (x ) lo es también (para comprobarlo basta derivar y(x ) y sustituir
en la ecuación).

El hecho de que en la función y(x ) así construida estemos incluyendo dos constantes C 1 y C 2 , debe
hacernos pensar que si sabemos encontrar y1 (x ) e y2 (x ) tendremos ya la solución general de la

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

6 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

homogénea: estamos incluyendo los dos grados de libertad que debe presentar toda solución
general de una ecuación lineal de segundo orden. Por supuesto, esto será así siempre que y1 (x ) e
y2 (x ) no sean proporcionales, es decir siempre que sean linealmente independientes.

Funciones linealmente independientes.- El conjunto de funciones

y (x ), y (x ), , y (x )
1 2 n

son linealmente independientes en el intervalo a, b  , si la ecuación

1y1 (x )  2y2 (x )    n yn (x )  0
admite como única solución 1  2    n  0 en a, b  .


Sistema fundamental de soluciones.- La pareja de soluciones y1 (x ), y2 (x ) , forman un 
sistema fundamental de soluciones de la ecuación y   p(x )y   q(x )y  0 , en a, b  si
son linealmente independientes en a, b  .

 
En el caso de dos funciones, y1 (x ), y2 (x ) , serán linealmente dependientes si existe una constante

k tal que y1(x )  k y2 (x ) y linealmente independientes si el cociente y1(x ) / y2 (x ) es una función

no constante de x . En la práctica es frecuente estudiar la dependencia lineal a través del
Wronskiano.

Wronskiano de un conjunto de funciones.- Se llama Wronskiano del conjunto de

 
funciones y1 (x ), y2 (x ), , yn (x ) al determinante funcional:

y1 y2  yn
y y2  yn
W (y1, y2 , , yn )  1
   
y1(n 1 y2(n 1  yn(n 1


Ejemplo: Calcular el wronskiano del conjunto de funciones cos x , sen x  
cos x  sen x 
W (cos x , sen x )   cos2 x   sen 2 x   1
sen x  cos x 

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
7

TEOREMA (Propiedad del Wronskiano de soluciones).- Sean y1 (x ) e y2 (x ) dos soluciones

de la ecuación homogénea y   p(x )y   q(x )y  0 en el intervalo a, b  , siendo p x  y
q x  continuas en a, b  .
• Estas soluciones son linealmente dependientes si y solo si el wronskiano es
idénticamente nulo en a, b 
• Estas soluciones son linealmente independientes si y solo si W (y1, y2 )  0 para todo x
en a, b  .

Ejemplo: Consideremos la EDO lineal homogénea y '' y  0 . Es fácil ver que y1  cos x  ,
y2  sen x  son soluciones de esa ecuación diferencial. Como hemos visto en el ejemplo anterior


W cos x , sen x   1  0 
Por lo tanto, son soluciones linealmente independientes en  de la EDO y '' y  0 .

TEOREMA (Solución general de la e. d. lineal homogénea de segundo orden).- Sean

y1 (x ) e y2 (x ) dos soluciones linealmente independientes de la ecuación homogénea
y   p(x )y   q(x )y  0 en el intervalo a, b  , siendo p x  y q x  continuas en a, b  .
Entonces la solución general de esa ecuación en a, b  es y(x )  C 1y1 (x )  C 2y2 (x )

Esto significa que el problema de encontrar la solución general de una ecuación lineal homogénea
de segundo orden está resuelto si se han encontrado dos soluciones no proporcionales. Los
resultados de los teoremas anteriores, son generalizables al caso de una ecuación lineal homogénea
de orden n .

Ejemplo: Teniendo en cuenta que y1  cos x  , y2  sen x  son soluciones linealmente

independientes de y '' y  0 , la solución general de esta ecuación diferencial lineal homogénea es
y(x )  C 1 cos x   C 2sen x  .

-1

-2

-3

-4
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

8 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

Ejemplo: Demostrar que e x como e x son soluciones de la EDO y ' '− y = 0 . Demuestre que esas
funciones son linealmente independientes y halle la solución general de la ecuación diferencial.
Solución

El determinante Wronskiano del conjunto e x , e x   es

ex e x
 
W e x , e x 
ex e x
   
 e x e x  e x e x  2  0

Como el Wronskiano es distinto de cero el conjunto e x , e x   es linealmente independiente y

forman un conjunto fundamental de soluciones. La solución general es y  C 1e x  C 2e x .

8 Reducción de orden

Por reducción de orden entendemos el proceso por el cual, dada una solución y1 (x ) de la ecuación

y   p(x )y   q(x )y  0

se encuentra otra solución y2 (x ) , de manera que y1 (x ) e y2 (x ) son independientes. En lo que sigue,

supondremos que las funciones p(x ) y q(x ) son continuas.

Si y1 (x ) es una solución de una ecuación lineal homogénea de orden 2, entonces la sustitución

y2 (x )  v(x )y1 (x )

la reduce a una de primer orden. Descripción del método:

1. Se escribe y2 (x ) = v(x )y1 (x ) , siendo el objetivo encontrar v(x ) .

2. Se deriva dos veces y2 (x ) = v(x )y1 (x ) y se impone que se verifique la ecuación homogénea:
y   p(x )y   q(x )y  0

3. Se agrupan los términos en v (x ) (sólo hay uno), en v (x ) y en v(x ) .

4. El término en v(x ) debe anularse por ser y1 (x ) solución, luego resulta una ecuación de
primer orden 1 en v (x ) , de la cual se obtiene esta función.

5. Se integra v (x ) para encontrar v(x ) y con ella la solución y2 (x ) requerida.

Ejemplo. Suponiendo que y1 (x )  x 2 es solución de la ecuación diferencial x 2y '' 7xy ' 20y  0
, encontrar la solución general en el intervalo 0,a  . Es fácil ver que y1 (x )  x 2 cumple
x 2y1 '' 7xy1 ' 20y1  0 , ya que y '1 (x )  2x 3 , y ''1 (x )  6x 4 .

1De aquí el nombre del método: para encontrar la solución de una de segundo se resuelve una de primer orden.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
9

Para obtener una solución linealmente independiente consideramos una solución de la forma

y2 (x )  x 2v x  . Se tendrá

y2 (x )  x 2v x  y '2 (x )  2x 3 v x   x 2v ' x 

y ''2 (x )  6x 4 v x   2x 3v ' x   2x 3v ' x   x 2v '' x 

Sustituyendo

x 2y2 '' 7xy2 ' 20y2  0

x 2 6x 4 v x   4x 3v ' x   x 2v '' x   7x 2x 3v x   x 2v ' x   20x 2v x   0
 
 

v '' x   v ' x  4x 1  7x 1   v x  6x 2  14x 2  20x 2   0

  
 
0

11 v '' x  11
v '' x   v ' x   0    log v ' x   11 log x 
x v ' x  x

x 12
 v ' x   x 11  v x  
12

x 12 x 10
Por lo tanto, y2 (x )  x 2  . La solución general entonces de la ecuación diferencial es
12 12

y(x )  C 1x 2  C 2x 10

ECUACIONES DIFERENCIALES LINEALES NO HOMOGÉNEAS

Como ya hemos visto, son ecuaciones de la forma

y   p(x )y   q(x )y  r (x )

En esta sección describiremos el método de variación de constantes para determinar una solución
particular de esta ecuación. Previamente se puede aplicar el principio de superposición para
desglosar la resolución de una ecuación de este tipo en procesos más sencillos.

8 Principio de superposición

PRINCIPIO DE SUPERPOSICIÓN.- Si y1 (x ) es solución de

y (x )  p(x )y (x )  q(x )y(x ) = r1 (x )
e y2 (x ) es solución de
y (x )  p(x )y (x )  q(x )y(x ) = r2 (x )
entonces 1y1 (x )  2y2 (x ) es solución de
y (x )  p(x )y (x )  q(x )y(x ) = 1r1 (x )  2r2 (x )

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

10 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

Ejemplo: Se tiene que

• y1  4x 2 es solución particular de y '' 3y ' 4y  16x 2  24x  8

• y2  e 2x es solución particular de y '' 3y ' 4y  2e 2x

Por lo tanto, y 3  y1  y2  4x 2  e 2x es solución de y '' 3y ' 4y  16x 2  24x  8  2e 2x

9 Método de variación de constantes

Con este método podremos encontrar una solución particular para una ecuación lineal no
homogénea, siempre que se conozca el sistema fundamental de soluciones de la ecuación
homogénea asociada.
• La idea, que da nombre al método, es sustituir las constantes de la familia solución general
de la homogénea C 1y1 (x )  C 2y2 (x ) por funciones, para hallar una solución particular de la
forma y p (x ) = C 1 (x )y1 (x )  C 2 (x )y2 (x )

• Puesto que deben hallarse dos incógnitas, necesitaremos una condición añadida, que
además permita simplificar los cálculos; esta condición es

C 1(x )y1 (x )  C 2(x )y2 (x ) = 0 (1)

• El primer paso del desarrollo del método consiste en derivar dos veces y p (x ) e imponer que
verifique la ecuación: y   p(x )y   q(x )y  r (x ) ; de ahí se obtendrá la ecuación en C 1(x )
y C 2(x ) siguiente

C 1(x )y1(x )  C 2(x )y2(x ) = r (x ) (2)

• Después, se resuelve el sistema formado por las ecuaciones (1) y (2) en las incógnitas C 1(x )
y C 2(x ) ; este sistema siempre tiene solución única puesto que las soluciones y1 e y2 son
no proporcionales.

• Por último se integran C 1(x ) y C 2(x ) para obtener la solución y p (x ) buscada.

1
Ejemplo: Consideremos la ecuación diferencial y '' y  x
. En la sesión 2 de la semana
e 1
anterior se demostró que la solución general de la homogénea asociada es y  C 1e x  C 2e x .
Vamos a encontrar la solución particular utilizando el método de variación de las constantes.
Para encontrar una solución particular de la completa se considera una solución particular de la
forma

y p  C 1 x e x  C 2 x e x considerando que se cumple C '1 x e x  C '2 x e x  0

Sustituyendo en la ecuación diferencial, se tendrá que cumplir

1
y p '' y p  x
e 1

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
11

Como

y 'p  C '1 x e x  C 1 x e x  C '2 x e x  C 2 x e x  C 1 x e x  C 2 x e x

y ''p  C 1 ' x e x  C 1 x e x  C '2 x e x  C 2 x e x

al sustituir en el ecuación diferencial

C 1 ' x e x  C 1 x e x  C '2 x e x  C 2 x e x  C 1 x e x  C 2 x e x  x

1
C 1 ' x e x  C '2 x e x  x
e 1

Para encontrar se resuelve el sistema

C ' x e x  C ' x e x  0
 1 2

C ' x e x  C ' x e x  1
 1 2
ex  1

Este sistema se resuelve utilizando la regla de Cramer suponiendo C '1 x  y C '2 x  como
incógnitas. Se tendrá,

0 e x
1
e x
x
e 1 e x e x 1 e x x
C '1 x     C 1 x    2 e dx 
 
log e x  1  
e x e x 2 ex  1  x
1  Hacer
cambio
2 2 2
ex e x e x t

ex 0
1
ex
C '2 x  
x
e 1 ex
 C 2 x  
ex 
log e x  1 
e x
e x


2 e x  1   2 e x
1 
dx  
2
x x
e e

Sustituyendo en y p  C 1 x e x  C 2 x e x se obtendrá la solución particular

1 x
e x x  x log e  1 x  
e x  e x x x 1

 2
x

y   log e  1 
2

 e 
2  2
e  log e x  1
2
 e 
2 2
 

ECUACIONES DIFERENCIALES LINEALES DE COEFICIENTES CONSTANTES DE ORDEN 2

Estas ecuaciones son un caso particular de las de coeficientes variables, sin embargo, tienen su
propio método de resolución, tanto para obtener la ecuación general de la homogénea asociada,
como para localizar una solución particular de la completa.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

12 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

10 Ecuaciones lineales de coeficientes constantes

Estas ecuaciones son de la forma: y (x )  py (x )  qy(x ) = R x  con p, q  

11 Ecuaciones lineales de coeficientes constantes homogéneas

En esta sesión resolveremos la ecuación diferencial homogénea y en la siguiente cómo obtener la

solución particular de la completa.
La forma general de las ecuaciones de este tipo es
y (x )  py (x )  qy(x ) = 0 con p, q   (3)

Introducción:
[Link]
equations/linear-homogeneous-2nd-order/v/2nd-order-linear-homogeneous-differential-
equations-2

Definición (Ecuación característica).- Se llama ecuación característica de la ecuación

y (x )  py (x )  qy(x ) = 0 , a la ecuación algebraica asociada: r 2  pr  q = 0

Conociendo las raíces de esta ecuación podemos formar la solución general de (3). Según que
p 2  4q sea positivo, nulo o negativo, nos encontraremos con tres situaciones diferentes:

CASO 1.- La ecuación característica tiene dos raíces reales distintas, r1 y r2 :

rx r x
Entonces el sistema fundamental de soluciones es y1 (x ) = e 1 , y2 (x ) = e 2

rx r x
y la solución general será: y(x ) = C 1e 1  C 2e 2

CASO 2.- La ecuación característica tiene una raíz real doble, r :

Entonces el sistema fundamental de soluciones 2 es y1 (x ) = e rx , y2 (x ) = xe rx

y la solución general será: y(x ) = e rx (C 1  C 2x )

CASO 3.- La ecuación característica tiene dos raíces complejas conjugadas, r1 = a  bi y

r2 = a  bi :

2La segunda solución se encuentra por el método de reducción de orden.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
13

Entonces el sistema fundamental de soluciones es:

y1 (x ) = eax cos bx , y2 (x ) = eax sen bx

y la solución general será: y(x ) = eax (C 1 cos bx  C 2 sen bx )

Ejemplo: Supongamos la ecuación diferencial y '' 5y ' 6y  0 . Para llegar al polinomio

característico, basta considerar que se buscan soluciones de la forma e rx . Si e rx es solución de la
ecuación diferencial se tendrá que cumplir

r 2e rx  5re x  6e rx  0 rx r 2  5r  6  0  r  2r  3  0  r1  2, r2  3

e 0

Entonces e 2x , e 3x son solución de la ecuación diferencial de coeficientes constantes. Además, se

puede comprobar que son linealmente independientes, por lo que la solución general de la ecuación
diferencial es

y  C 1e 2x  C 2e 3x

Nota: Estamos en el caso 1, donde el polinomio característico tiene dos raíces reales distintas.
Ejemplo: Supongamos la ecuación diferencial y '' 2y ' y  0 . Para llegar al polinomio
característico, basta considerar que se buscan soluciones de la forma e rx . Si e rx es solución de la
ecuación diferencial se tendrá que cumplir

r 2e rx  2re x  e rx  0 rx r 2  2r  1  0  r  1  0  r  1 doble 

2

e 0

Entonces y1  e x sería solución de la ecuación diferencial de coeficientes constantes. Como

necesitamos otra solución linealmente independiente consideramos y2  e x v x  . Es fácil ver que
utilizando el método de reducción de orden v x   x , y por tanto, y2  xe x .

La solución general de la ecuación diferencial será entonces

y  C 1e x  C 2xe x

Nota: Estamos en el caso 2, donde el polinomio característico tiene una raíz real doble.
Ejemplo: Supongamos la ecuación diferencial y '' 4y ' 5y  0 . Para llegar al polinomio
característico, basta considerar que se buscan soluciones de la forma e rx . Si e rx es solución de la
ecuación diferencial se tendrá que cumplir

4  16  20
r 2e rx  4re x  5e rx  0 rx r 2  4r  5  0  r   2  i
e 0 2
2i x 2i x
Entonces y1  e , y1  e puede demostrarse que son linealmente independientes y que
la solución general será una combinación lineal de ellas. Teniendo en cuenta que
2i x
y1  e  e 2x cos x  isenx 

2i x
y2  e  e 2x cos x  isenx 

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

14 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

 
También sería un sistema fundamental e 2x cos x , e 2x senx , lo que permite no trabajar con números
complejos.
La solución general de la ecuación diferencial será entonces

y  C 1e 2x cos x  C 2e 2x senx

Nota: Estamos en el caso 3, donde el polinomio característico tiene dos raíces complejas conjugadas.

12 Ecuaciones lineales de coeficientes constantes no homogéneas

La forma general de las ecuaciones de este tipo es

y (x )  py (x )  qy(x ) = R(x )

Sabemos que la solución de esta ecuación se obtiene mediante la suma y p  yh donde yh es la

solución general de la homogénea asociada, que ya podemos calcular, e y p es una solución

particular de la completa. Esta solución particular puede calcularse a partir de la general de la
homogénea, por el método de coeficientes indeterminados o por el método de variación de
constantes descrito en el apartado de las ecuaciones lineales de coeficientes variables.

MÉTODO DE COEFICIENTES INDETERMINADOS

Este método se aplica únicamente cuando los coeficientes son constantes; su utilización permite
construir una solución particular de la ecuación

y (x )  py (x )  qy(x ) = R(x )

cuando r (x ) es combinación lineal o producto de funciones polinómicas, exponenciales, cosenos y

senos.
La idea del método es buscar una solución particular del mismo tipo que el término R(x ) , puesto
que la parte izquierda de la ecuación es únicamente una combinación lineal de y(x ) y sus derivadas.

Es imprescindible haber resuelto la ecuación característica de la homogénea asociada y disponer ya

del sistema fundamental de soluciones de esa ecuación.

El primer paso del método es el diseño de una propuesta de solución particular, y p (x ) , del mismo
tipo que la función R(x ) (Ver cuadro). Esta función y p (x ) dependerá de unos coeficientes aún por
determinar 3.

El segundo paso es el cálculo de esos coeficientes, para lo cual hay que derivar dos veces y p (x ) e
imponer que cumpla la ecuación; se obtendrá un sistema lineal de ecuaciones de la misma dimensión
que el número de coeficientes a determinar.

3
De aquí el nombre de este método.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
15

En el diseño de y p (x ) es muy importante tener en cuenta que no debe contener ningún sumando
que sea solución de la homogénea asociada.
En virtud del principio de superposición, a una suma de funciones R(x ) de las del cuadro, le
corresponderá una suma de las soluciones propuestas para cada sumando.
Lo mismo puede indicarse si el segundo término de la ecuación es producto de funciones R(x ) de
las del cuadro, teniendo en cuenta que habrá que añadir un único factor x s .

R(x ) y p (x )

ae kx Ax se kx

b0  b1x    bm x m x s (B0  B1x    Bm x m )

a 0 cos kx  a1 sen kx x s (A0 cos kx  A1 sen kx )

en todos los casos s es el menor entero no negativo tal que

ningún sumando de y p (x ) sea solución de la homogénea
asociada.

Cuadro 1.- Diferentes soluciones particulares según diferentes valores de R(x ) .

Hemos visto que en el caso de las ecuaciones de coeficientes constantes la búsqueda de la solución
particular puede hacerse por dos métodos: la ventaja del de variación de constantes respecto al de
coeficientes indeterminados es que no restringe el tipo de R(x ) ; su desventaja es que conlleva
integración, lo que en muchos casos complica su aplicación respecto al de coeficientes
indeterminados, en el que no se calculan primitivas.

APLICACIONES: VIBRACIONES EN SISTEMAS MECÁNICOS

Una masa m está sujeta a una pared mediante un resorte o muelle. Este muelle ejerce sobre el objeto
una fuerza de recuperación, Fm. Las vibraciones del objeto se producen a consecuencia de la pérdida
del equilibrio del sistema; las fuerzas presentes en el sistema tenderán a recuperar el equilibrio.
Dependiendo de cuáles sean esas fuerzas se presentan distintos casos.

13 Vibraciones armónicas simples no amortiguadas

Supongamos que la superficie es lo suficientemente lisa como para despreciar la fuerza de

rozamiento.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

16 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

La única fuerza presente en el sistema es la del muelle 4: Fm  k x t  donde k es una constante

positiva que depende de la rigidez del muelle.
La segunda ley de Newton establece que la fuerza total que actúa sobre un cuerpo es igual al
producto de su masa por su aceleración, esto es

mx '' t   Fm

El problema de valor inicial que modela las vibraciones del objeto es

k
x '' t   x t   0, x 0  xo , x ' 0  0
m

14 Vibraciones armónicas amortiguadas

En el sistema anterior consideremos que se agrega un sistema mecánico (amortiguador) que tiene
el efecto de reducir la velocidad cuando el sistema se encuentra vibrando. El amortiguador ejerce
una fuerza dependiendo de la velocidad de la masa.

Por simplicidad supongamos que esta fuerza en magnitud es proporcional a la rapidez, entonces la
fuerza que ejerce el amortiguador es

Fc  c  x ' t  c  0

Observad que el signo negativo indica que la fuerza de amortiguación va en sentido contrario a la
del cuerpo. La fuerza total ejercida por la masa es Fm  Fc .

La ecuación para este sistema 5 es

m x '' t   c x ' t   k x t   0

La ecuación característica de la ecuación diferencial es

m r2  c r  k  0

Las dos soluciones de esta ecuación son:

c  c 2  4mk c  c 2  4mk
r1  r2 
2m 2m

El signo del radicando c 2  4mk determina el tipo de movimiento del sistema:

4 La fuerza que actúa sobre un objeto causada por un resorte tiene una magnitud proporcioanl a la elongación del
resorte a partir de su longitud natural y de sentido opuesto a la elongación (ley de Hooke)
5 La misma ecuación diferencial modela el sistema masa-resorte colocado verticalmente

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
17

1. Movimiento sobreamoriguado c 2  4mk  0 . En este caso

x t   c1e 1  c2e 2
rt rt

Las dos funciones exponenciales son decrecientes y el sistema tiende rápidamente a su

posición de equilibrio.

2. Movimiento críticamente amortiguado c  4mk . En este caso

c

x t   c1  c2t e 1  c1  c2t e
rt t
2m

La función posición contiene un término exponencial decreciente multiplicada por una

función lineal en el tiempo. Se espera que la posición decrezca hacia la posición de equilibrio
sin vibrar, la manera qn la que lo haga dependerá de las condiciones iniciales.

3. Movimiento subamortiguado c 2  4mk  0 . En este caso

    

c
  4mk  c 2  4mk  c 2 
x t   e
t
2m
c1 cos  t   c2 sen  t 
  2m   2m 


El sistema oscilará alrededor de la posición de equilibrio. Como el factor exponencial es

decreciente, se espera que la amplitud de vibración sea cada vez más pequeña.

15 Vibraciones forzadas

Además de las fuerzas de recuperación y de rozamiento sobre el sistema puede actuar una fuerza
externa dada por Fe  f t  . Se tendrá entonces

mx '' t   Fm  Fc  Fe

es decir,

mx '' t   c x ' t   k x t   f t 

ECUACIONES DIFERENCIALES LINEALES DE COEFICIENTES CONSTANTES HOMOGÉNEAS DE

ORDEN N

16 Soluciones

En general, para una ecuación lineal de orden n de la forma

n 
y  a1y (n 1)    an 1y   an y  0

la ecuación característica, es la siguiente ecuación de grado n :

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

18 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

r n  a1r n 1    an 1r  an  0

Además, la solución general de la ecuación se escribe como una combinación lineal de n soluciones
linealmente independientes, es decir, en base a un sistema fundamental de n soluciones.
Generalizando lo expuesto para orden 2, veamos qué relación guardan las raíces de la ecuación
característica con las n soluciones que forman el sistema fundamental:

1. A cada raíz real r de multiplicidad s , le corresponden las soluciones independientes

e rx , xe rx , x 2e rx , , x s 1e rx

2. A cada par de raíces complejas conjugadas a  bi de multiplicidad s , le corresponden

las soluciones independientes

eax cos bx , eax sen bx , xeax cos bx , xeax sen bx , , x s 1eax cos bx , x s 1eax sen bx

SISTEMAS DE ECUACIONES DIFERENCIALES ORDINARIAS

17 Definiciones básicas

Un sistema de orden 2 con condiciones iniciales tiene la forma

x ' (t )  f t, x , x  x 1 to   y1
 1 1 1 2
 '
x 2 (t )  f2 t, x 1, x 2  x 2 to   y2

En el caso de que f1 y f2 sean lineales entonces el sistema se denominará lineal
x ' (t )  a t  x t   a t  x t   g t  x 1 to   y1
 1 11 1 12 2 1
 '
x 2 (t )  a21 t  x 1 t   a22 t  x 2 t   g 2 t  x 2 to   y2

Matricialmente,
x ' (t ) a t  a t  x (t ) g (t )
 1   11  1   1 
 '     
12
 
x 2 (t ) a21 t  a22 t  x 2 (t ) g 2 (t )

18 Resolución de sistemas lineales

El método que desarrollaremos para la resolución analítica de un sistema de ecuaciones lineales de

coeficientes constantes consiste en eliminar una de las variables junto con su derivada y resolver la
ecuación diferencial de segundo orden que resulta del proceso de eliminación. Sustituyendo en el
sistema la solución obtenida, se calcula la función que falta.

Denotamos por D la derivación respecto de x , con D 2 la derivación respecto de x dos veces y, en

general, por D n la derivación respecto de x , n veces.
Así,

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
19


y   py   qy  D 2y  pDy  qy  D 2  pD  q y 
Este operador permite utilizar técnicas algebraicas en la eliminación.

x (t )  3x (t )  4y(t )

Por ejemplo, sea el sistema 
y (t )  4x (t )  7y(t )

Reescribimos el sistema utilizando el operador D
(D  3)x  4y  0



 4x  (D  7)y  0

En el sistema algebraico obtenido, se elimina una de las variables. Si multiplicamos la primera

ecuación por (D  7) y la segunda por −4 , se elimina la variable y , obteniéndose la ecuación de
segundo orden

x (t )  4x (t )  5x  0

Esta ecuación se resolverá mediante la ecuación característica. Una vez obtenida y t  , se substituye
en la primera ecuación del sistema para calcular y t  .

Ejercicios propuestos

1 3
Halla, mediante eliminación de los En los siguientes problemas de valor
parámetros, la e.d.o. de segundo orden cuya inicial, determinar si se puede aplicar el
solución general es la familia teorema de existencia y unicidad de las
C ecuaciones lineales y en caso afirmativo definir
y = C 1x 2  2 el intervalo donde se garantiza la existencia de
x
solución única.
Solución: x y   2y = 0
2
a. y   yy   x 2  1 ;
y(0)  1, y (0)  1
2
Utilizar el teorema de existencia para b. (1  x )y   xy   2y  sen x ;
discutir la existencia y unicidad de la solución de y(0)  1, y (0)  1
los siguientes problemas de valor inicial:

x 2y   y  cos x Solución: a) No se puede aplicar el teorema de
a)   existencia porque no es un problema lineal; b)

 y(1)  y 0 , y (1)  y1 La solución es única y está definida en (,1)


 x
e y   1 y   y  log x
b)  x 3 4 1

y(1)  y 0 , y (1)  y1 a. Demostrar que y1 (x ) 
x
e

donde y 0 e y1 son constantes reales. y2 (x )  x 3 son soluciones de la ecuación lineal

Solución: a) Tiene solución única definida en homogénea x 2y   xy   3y  0 en el
0,  ; b) Tiene solución única definida en intervalo 0,  .
0, 3 . b. Probar que son soluciones linealmente
independientes en 0,  .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

20 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

c. Escribir la solución general.

8
1 Partimos de nuevo de la familia de
Solución: c) y  C 1  C 2x 3
x curvas y(x ) = C 1e x  C 2e x que es la solución
general de y  = y , no sólo porque cada una de
5 las curvas es solución de la ecuación (como
Determinar si las siguientes funciones comprobaste antes), sino también porque para
pueden ser wronskianos en 1,1 de dos cada par de condiciones iniciales y(0) = y 0 ,
soluciones de alguna ecuación lineal
y (0) = y1 existen unos valores adecuados de
homogénea de segundo orden con p(x ) y q(x )
continuas. C 1 y C 2 de forma que la curva que los lleva
a) w(x )  3e 2x b) w(x )  x 2 cumple esas condiciones.
a. Representa con Matlab, en el intervalo
1  
c) w(x )  d) w(x )  0 1,1 , una muestra de soluciones de
x 1
e) ¿Y la función w(x )  3(x  1)2 en el
y  = y , primero fijando un punto de paso,
P (0, 2) , y dejando variar la pendiente y
intervalo 0, 2 ?
después fijando una pendiente, y (0)  2 ,
Solución: a) Sí, b) No, c) Sí, d) Sí, e) No.
y dejando variar el punto de paso.
b. Adapta el fichero anterior para que repita el
6 mismo proceso para la familia de curvas
Considerar dos funciones linealmente
ax
y(x ) = e x (C 1  C 2x ) . Para ello debes
independientes de la forma e y encontrar la
primero encontrar las fórmulas que dan C 1
ecuación diferencial homogénea cuya solución
general sea una combinación lineal de ambas. y C 2 en función de y(0) e y (0) .
c. Repite el apartado anterior para la familia
7 de curvas y(x ) = e x (C 1 cos x  C 2 sen x )
La expresión
y(x ) = C 1e x  C 2e x
se llama familia biparamétrica de curvas porque 9
Los coeficientes de la familia de curvas
cada par de valores reales (C 1,C 2 ) proporciona y(x ) = C 1e x  C 2e x
una curva diferente. también pueden ajustarse para que la función
a. Comprueba que todas las curvas de la pase por dos puntos dados. Éstas son las
familia anterior verifican la relación condiciones de frontera.
y  = y .
a. Representa, en el intervalo  0, 2 , una
b. Representa con Matlab, en el intervalo
2, 2 , una muestra de curvas de esa muestra de cinco curvas de la familia
  pasando por el punto P (0, 2) , es decir
familia y comprueba gráficamente que son
cumpliendo y(0) = 2 , y luego cinco curvas
cóncavas en el semiplano y > 0 y
pasando por Q(2, 2) , es decir, cumpliendo
convexas en el semiplano y < 0 .
¿Presentan puntos de inflexión? Toma para que y(2) = 2 . Adapta ese fichero para que
C 1 y C 2 los valores 2, 1, 0,1, 2 . repita el mismo proceso para la familia de
curvas
c. Modifica el programa anterior para que
y(x ) = C 1 cos x  C 2 sen x
dibuje una muestra de curvas de la familia
y(x ) = C 1 cos x  C 2 sen x y los puntos P (0, 2) y Q( / 2,1) . Para
d. omprueba que todas las curvas de esta ello debes primero encontrar las fórmulas
familia son soluciones de y  = y que dan C 1 y C 2 en función de y(0) e
y( / 2) .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
21

10 x 1
y p (x ) = x 2  x log | |
Sabiendo que la función dada, y1 (x ) , es x 1
solución de la ecuación, encuentra la solución x2  1
 log | x 2  1 |
general de las siguientes ecuaciones: 2
a. (x  1)y   xy   y = 0 , y1 (x ) = x
b. x 2y   2xy   2y = 0 , y1 (x ) = x 14
Efectúa el cambio de variable x = e t
c. x y   xy   4y = 0 , y1 (x ) = x
2 2
para transformar las ecuaciones siguientes en
Comprueba los resultados con Matlab, otras de coeficientes constantes y resuelve
utilizando dsolve. éstas últimas.
Solución: a) y(x ) = C 1x  C 2e x a. x 2y   3xy   4y = 0
b. x 2y   4xy   6y = x 3
b) y(x ) = C 1x  C 2x 2
Comprueba los resultados con Matlab,
c) y(x ) = C 1x 2  C 2x 2 utilizando dsolve.
Solución: a) y(x ) = x 2 (C 1  C 2 log x )
11 b) y(x ) = x 3 (C 1  log x )  C 2x 2
Dada la ecuación diferencial
1
x 2y '' 4xy ' 6y  . Se pide
x
15
a. Demostrar que y1  x 2 , y2  x 3 es un Halla la solución particular de los
sistema fundamental de soluciones de la siguientes problemas de valor inicial
ecuación homogénea.
a. y   10y   25y = 0 , y(0) = 3 ,
b. Encontrar la solución general de la ecuación
diferencial. y (0) = 13
1 b. y   6y   13y = 0 , y(0) = 2 ,
Solucion b) y1   C 1x 2  C 2x 3
12x y (0) = 0
c. y   y   6y = 0 , y(0) = 7 ,
12 y (0) = 1
Dada la ecuación diferencial
x 2y '' xy ' y  4x log x . Se pide Comprueba los resultados con Matlab,
utilizando dsolve.
a. Demostrar que y1  x , y2  x log x es un
Solución: a) y(x ) = e 5x (3  2x )
sistema fundamental de soluciones de la
ecuación homogénea. b) y(x ) = e 3x (2 cos 2x  3 sen 2x )
b. Encontrar la solución general de la ecuación
diferencial. c) y(x ) = 4e 2x  3e 3x
2
Solucion b) y1  x  log 3 x  C 1x  C 2x log x
3 16
Encuentra la solución general de cada
una de las siguientes ecuaciones lineales
13 a. y   y  = 1 / x
Para la siguiente ecuación, encuentra la
solución particular que se obtiene mediante el b. y   2y   y = x 2e x cos x
método de variación de constantes c. y   y = (e x  1)1
(x 2  1)y   2xy   2y = x 2  1 d. y   y = 2x sen x
sabiendo que la solución general de la
homogénea asociada es
Comprueba los resultados con Matlab,
yh (x ) = C 1x  C 2 (1  x 2 ) utilizando dsolve.
Comprueba los resultados con Matlab, Solución:
utilizando dsolve. ex
Solución: a) y(x ) = C 1  C 2e x  log x  e x  dx
x
b) y(x ) = e x [C 1  C 2x 

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

22 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

 4x sen x  (6  x 2 ) cos x posición y la velocidad de la masa en el tiempo,

la frecuencia de oscilación, la amplitud, el
c) y(x ) = C 1e x  C 2e x  ángulo de fase y las energías cinética y potencial
en el tiempo t.
1 1 El problema que modela esta situación es:
 (e x  e x ) log(e x  1)  (xe x  1)
2 2 10x '' 10x  0, x 0  0.02, x ' 0  0 .
d) y(x ) = C 1 cos x  C 2 sen x   
Solución: x t   0.02 cos t  0.02 sen t   .
x 1 x 2   2 
sen x     cos x
2  4 2  w  1rad / seg. , T  2seg. , A  0.02 m. , la
1 1
frecuencia de oscilación f   osc/seg
17 T 2
Encuentra la solución general de 
cada una de las siguientes ecuaciones (Hertz), el ángulo de fase es   ,
2
lineales
1
a. y ' ' 3y ' 4y  x  2 Ec  mv 2  0.002 sen2 t ,
2
b. y ' ' y '  3x 2  x 1
E p  kx 2  0.002 cos2 t
c. y ' ' 2y ' 5y  2e x  3 sen x  3x 2  x 2

d. y ' ' 7y ' 12y  e 2x (x 3  5x 2 )

20
Solución Un resorte cuelga verticalmente de un
techo, el resorte sufre una deformación en su
x 11
a. y(x ) = C 1e x  C 2e 4x   longitud de un centímetro cuando se coloca una
4 16 masa de 1.5 kg y después el sistema queda en
7x 2 equilibrio. Posteriormente se elonga una
b. y(x ) = C 1  C 2e x  x 3   7x
2 cantidad de 1.5 cm. y se suelta a partir del
reposo. Determinar la constante del resorte, la
1
c. y(x ) = e x C 1 cos 2x  C 2 sen 2x   e x  posición y la velocidad de la masa en el tiempo
4 t  0 . ¿Cuál es la frecuencia de oscilaciones de
3 3 3 7 16 la masa y la amplitud del movimiento?
 cos x  sen x  x 2  x 
10 5 5 25 125 Nota: En la posición de equilibrio, la constante
x 3 x 2 9 25  mg
d. y(x ) = C 1e 3x  C 2e 4x  e 2x    x   del resorte es: k  . Considerando como
 2 4 4 8  l
origen de coordenadas a la posición de
equilibrio y la dirección positiva del eje vertical
18
Halla la solución general de la siguiente hacia abajo, entonces, la posición de la masa
ecuación, utilizando el método de coeficientes viene dada por la solución del problema de
indeterminados para encontrar una particular valor inicial
de la ecuación completa mx '' kx  0, x 0  xo , x ' 0  vo
y   9y =sen 4x
Solución:
x t   0.015 cos  980 t 
Nota: Debes escribir sen 4 x como combinación
de cosenos.
v t   0.015 980 sen  980 t 
Solución: 980
y(x ) = C 1 cos 3x  C 2 sen 3x  w  980 rad / seg , f   4.98
2
1 1 1 oscilaciones en un segundo. La amplitud es
  cos 2x  cos 4x
24 10 56 A  0.015m
19
Considerar una masa de 10 Kg. Que está
21
unida a una pared una distancia de 0.02 m. y se Un péndulo de masa m  2kg y de
suelta a partir del reposo, determinar la longitud 2.45 m. está suspendido sobre un

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
23

marco horizontal. El péndulo se levanta un Solución:

ángulo de 10º y se suelta con una velocidad  1
angular de -0.4 rad/seg. Se pide:  t   cos 2t  sen 2t  A sen wt  
18 5
a. ¿cuántos ciclos (oscilaciones) completos
2 2
habrá completado el péndulo después de    1
10 seg.? siendo A       ,
18   5 
b. ¿en qué momento la masa pasa por la
posición de equilibrio con velocidad  / 18
tg   ,   2.4242 , w  2 , T   ,
angular positiva por tercera vez? 0.2
c. ¿En qué instantes la masa alcanza sus en 10 segundos se realizarán 10 /   3 ciclos,
desplazamientos extremos?
 10  0.01
d. ¿Cuál es la posición de la masa del péndulo
en cualquier tiempo?
e. ¿Cuál es la posición de la masa a los 10 22
seg.? Encontrar la solución general de los
siguientes sistemas de ecuaciones:
x   3y
a.   
y  2x  y


 x   y  t2
b. 

 y  x  1


x   y   2x  0
c.   
x  y   x  y  sen t

Nota: La ecuación diferencial del moviemiento y x
Solución: e.d.o. y   ;
g y x
es  '' sen   0 . Suponiendo ángulos
l 1 y
pequeños se puede suponer sen    , 2
 
log x 2  y 2  arctg  C
x
obteniendo la ecuación diferencial
g
 ''   0 .
l

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

24 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

Test de autoevaluación

1 3
Si en la ecuación xy   y   9x 2 se Dado el problema de valores en la
t
realiza el cambio de variable x  e , se obtiene frontera,
la ecuación: y   16y  0, y(0)  0, y( / 2)  1
A) y   9e 3t podemos afirmar que:
B) y   9e 2t A) Una solución de este problema es
y  0.
C) y   2y   9e 3t
B) Este problema no tiene solución.
D) Ninguna de las anteriores.
C) Este problema tiene infinitas
soluciones.
2 D) Ninguna de las anteriores.
Se considera la ecuación diferencial
y   y   0 . Justificar cuál de los conjuntos
4
y (x ), y (x ) son un sistema de soluciones
1 2 Dada la ecuación diferencial,
linealmente independientes: y   3y   2y  4x 2
A) y1 (x )  0 , y2 (x )  e x podemos afirmar que una ecuación particular es
de la forma:
B) y1 (x )  xe x , y2 (x )  e x
A) y p  Ax 2
1 x
C) y1 (x )  1  2e x , y2 (x )  e B) y p  x 2e x .
4
D) Ninguna de las anteriores. C) y p  ax 2  bx  c
D) Ninguna de las anteriores.

1 2 3 4
A C B C

Ejercicios resueltos
DEFINICIONES BÁSICAS

1
Encontrar mediante eliminación de constantes, la ecuación diferencial de la familia
y = C 1x  C 2x 2

Solución
Derivando la expresión dos veces, se tiene un sistema que nos permite despejar las constantes

 C  1 y '' 

y '  C 1  2C 2x 
 2 

 2 
y ''  2C 2  

 C 1  y ' y '' x 
 

Así sustituyendo en la expresión inicial
1 
y  y ' y '' x  x   y '' x 2  x 2y '' 2xy ' 2y  0
 2 

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
25

se obtiene la ecuación diferencial de la cual es solución la familia biparamétrica dada.

2
Demostrar que y1 x   0 , y2 x   x 2 sen x  son soluciones de la ecuación differencial

 
x 2y '' 4xy ' x 2  6  0 . ¿Contradice el teorema de existencia y unicddad de soluciones?

Solución

Es sencillo probar que y1 x   0 es solución de la ecuación, así como ver que

y1  0   y1 '  0   0 .

Para comprobar que y2 x   x 2 sen x  lo es, basta derivar

y '2 x   x 2 cos x   3x sen x 

 
y ''2 x   4x cos x   2  x 3 sen x 

 
y susituir y2 , y2' , y2'' en x 2y '' 4xy ' x 2  6  0 para ver que se anula. Fácilmente se obtiene

que y2' 0  y2'' 0  0 .

Para poder aplicar el teorema de existencia y unicidad, escribimos la ecuación de la forma.

4  6
y '' y ' 1  2  y  0
x  x 

4
Observando esta expresión, se deduce que el teorema no puede aplicarse porque p x   y
x
6
q x   1  no están definidas para x=0. Por lo tanto, no hay contradicción en la existencia de
x2
dos solucionjes pasando por el mismo punto con la misma pendiente en ese punto.

RESOLUCIÓN CON MATLAB

dsolve(ecuación,condiciones)
Devuelve la solución de la ecuación diferencial con las condiciones iniciales o de frontera
especificadas.
syms y(x)
der(x)=diff(y);der2(x)=diff(y,x,2);
ecuacion=der2(x)-x*der(x)+y==0;
cond1=der(0)==1;cond2=der2(0)==0;
sol(x)=dsolve(ecuacion,cond1,cond2)
fplot(sol(x),[-3,3])

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

26 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

3
Halla y representa en 0.5, 0.5 la solución de los siguientes problemas de valor inicial,
 
utilizando el comando dsolve:

a) y   10y   25y = 0 , y(0) = 3 , y (0) = 13

b) y   6y   13y = 0 , y(0) = 2 , y (0) = 0

c) y   y   6y = 0 , y(0) = 7 , y (0) = 1

Solución a)

syms y(x)
eqn= diff(y,x,2)-10*diff(y,x)+25*y==0;
Dy=diff(y,x);
cond=[y(0)==3,Dy(0)==13]
sol(x)=dsolve(eqn,cond)
fplot(sol(x),[-0.5,0.5])

20

15

10

-0.5 0 0.5

EDOS LINEALES

4
Determinar si la familia y  C 1x  C 2x 2 es solución general de la ecuación
x 2y '' 2xy   2y  0 en cualquier intervalo que no contenga al cero.

Solución
Para que una combinación lineal de y1 x   x e y2 x   x 2 sea la solución general de la
ecuación diferencial se tendrá que demostrar que son soluciones linealmente independientes.

• Soluciones. Se puede comprobar fácilmente que y1 x   x e y2 x   x 2 verfican la

ecuación, es decir, se cumple:

x 2y1  2xy1  2y1  0 x 2y2  2xy2  2y2  0

• Linealmente independientes. Calculando el wronskiano

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
27

x x2
W y1, y2    x 2  0 para todo x  a, b  , 0  a, b 
1 2x
se deduce que y son linealmente independientes en cualquier intervalo que no contenga
al cero.

Constituyen por tanto un sistema fundamental de soluciones.

5
Encontrar la solución general de la ecuación x 2y '' x  2 xy ' x  2 y  0 sabiendo
que y1 x   x es solución de la misma.

Solución
Se debe buscar una solución independiente con y1 x  de la forma y2 x   v x  y1 x  .
Derivando e imponiendo que y2 x  verifique la ecuación, se obtiene la expresión

   
x 2 vy1''  2v 'y1  v '' y1  x  2 x vy1'  v ' y1  x  2 vy1  0
donde se ha omitido, para mayor claridad, la dependencia de la x en las distintas funciones.
Agrupando los términos en v, v’ y v’’

    
v '' x 2y1  v ' 2x 2y1'  x x  2 y1  v x 2y1''  x  2 xy1'  x  2 y1  0 
De esta forma observamos que el coeficiente env es nuelo por ser y1 x  solución de la ecuación
x 2y '' x  2 xy ' x  2 y  0 .

La expresión se reduce a


v '' x 3  v ' 2x 2  x 2 x  2  0  v '' v '  0

que es una ecuación lineal de primer orden, cuya función incógnita es v ' x  . Resolviendo esta
última ecuación, sabemos que

v ' x   ke x  v x   Ke x  C

De esta expresión se puede obtener la forma general de y2 x  ,

y2 x   Kxe x  Cx

Dado que basta encontrar una función independiente con y1 x  , se hará K=1 y C=0 para concluir
que podemos tomar como expresión de la solución general de la ecuación

y x   C 1x  C 2xe x

6 6
Dada la ecuación diferencial x 2y '' 8xy ' 12y  . Se pide
x2

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

28 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

a) Demostrar que y1  x 3 , y2  x 4 es un sistema fundamental de soluciones de la ecuación

homogénea.
b) Encontrar la solución general de la ecuación diferencial.

Solución a)

Basta ver que y1  x 3 , y2  x 4 son soluciones de la ecuación diferencial homogénea

(cumplen la ecuación diferencial homogénea) y son linealmente independientes.

Solución b)

8 12 6
Como la ecuación diferencial es lineal y '' y ' 2 y  4 la solución general será
x x x

yG  yH  y p siendo yH solución general del homogéneo e y p solución particular de la

completa.

Por el primer apartado se tiene que yH  C 1x 3  C 2x 4

Para obtener yp se considera y p  C 1 x  x 3  C 2 x  x 4 con la condición

C '1 x  x 3  C '2 x  x 4  0 . Sustituyendo en la ecuación diferencial se deberá cumplir:

C '1 x  x  C '2 x  x  0

 3 4
C ' x  x 3  C ' x  x 4  0
 
 1 2
 6 

3 C '1 x  x  4C '2 x  x  4
4 5 
3 C '1 x   4C '2 x  x 1  6


 x 


0 x 4
6 x 1 6x 4
C 1 x    6  C 1 x   6x
x 3 x 4 x 4
3 4x 1

x 3 0
3 6 6x 3
C 2 x     6x  C 2 x   3x 2
x 3
x 4
x 4
3 4x 1

Se tiene que y p  6x  x 3  3x 2  x 4  3x 2 .

Por lo tanto, la solución general de la ecuación diferencial es yG  C 1x 3  C 2x 4  3x 2

7
Se consideran y1(x ) e y2 (x ) soluciones de la ecuación y '' x 2y '  1 . Determinar,

 
justificadamente, si y1 x  2y2 (x ) e y1 x  2y2 (x ) son también soluciones de la ecuación
diferencial.
Solución

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
29

y  y1 x   2y2 (x ) sí es solución ya que

y '' x 2y '  y1 '' x 2y1 ' 2y2 '' 2x 2y2 '  1  2  1

y  y1 x   2y2 (x ) no es solución ya que

y '' x 2y '  y1 '' x 2y1 ' 2y2 '' 2x 2y2 '  1  2  1

8
Encontrar la solución general para la siguiente ecuación

x   
2
2
 1 y '' 2xy ' 2y  x 2  1

sabiendo que y1 x   x es solución de la ecuación lineal homogénea asociada.

Solución
Se recurrirá a la reducción de orden para encontrar otra independiente con ella. Se considera

y2 x   xv x   y2' x   xv ' x   v x  y2''  2v ' x   xv '' x 

Por lo tanto,

v '' 2  1 
  v x   K x   C 
v' 
x x2 1   x 

de donde

y2 x   x 2  1  
yH x   c1x  c2 x 2  1 
Una vez obtenida esta solución general de la homogénea, se buscará una particular de la
completa. Se hará mediante variación de parámetros:

 
y p x   c1 x  x  c2 x  x 2  1 con c1' x  x  c2' x  x 2  1  0  
Derivando y sustituyendo en la ecuación,

y p' x   c1 x  x  2xc2 x 

  c1' x   2xc2' x   x 2  1
y p'' x   c1' x   2xc2' x  

Así, se dispone del sistema de dos ecuaciones con dos incógnitas siguiente

c1' x   2xc2' x   x 2  1 

 c1' x   x 2  1
 

c1' x  x  c2' x  x 2  1  0
 c2' x   x 

Integrando se tiene

x3 

c1 x     x
 x2 2

x 2
3 


 y p x  
6
x 3  
c1 x   

2 



Prof. Elena E. Alvarez Saiz Universidad de Cantabria

30 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

Retomando la solución general de la homogénea, se concluye finalmente que la solución general

de la completa es

x2 2

y x   c1x  c2 x 2  1   6

x 3 
9
Justificar la certeza o falsedad de las siguientes afirmaciones:
2
(a) La siguiente ecuación x 2y '' xy ' 5x  e x es una ecuación diferencial lineal
homogénea.
x
(b) La ecuación y1(x )  e es solución de la ecuación xy   (x  2)y   2y  0 . Si
y2 (x )
y2 (x ) es otra solución de esta ecuación, la función v(x )  cumple que es
y1(x )
solución de la ecuación xv   (x  2)v   0

(c) Si y1(x ) e y2 (x ) son solución de la ecuación y '' x 2y '  1 entonces y1 x  y2 (x ) 

también es solución de la ecuación diferencial.
Solución a)
Se trata de una ecuación diferencial lineal por ser de la forma

y '' p x  y ' q x  y  r x 
2 2
1 e x  5x 1 e x  5x
y '' y '  p x   q x   0 r x  
x x2 x x2
Pero, no es homogénea ya que r x  0 
Solución b)
x
Es fácil ver que y1(x )  e es solución: xe x  (x  2)e x  2e x  0 .

 
Considerando y2 (x )  v x y1 x y teniendo en cuenta que

y2 (x )  v x e x

y '2 (x )  v ' x e x  v x e x

y2'' (x )  v '' x e x  2v ' x e x  v x e x

Como y2 (x ) es solución de la ecuación diferencial, se tendrá que cumplir

xy2'' (x )  x  2 y2'' (x )  2y2 x   0

Sustituyendo

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
31

   
x v '' e x  2v ' e x  ve x  x  2 v ' e x  ve x  2ve x  0

Dividiendo por e x y operando

xv '' 2xv ' xv   xv ' 2v ' xv  2v   2v  0

x v '' x  2 v '  0 ¨

La afirmación es cierta
Solución c)


Si y1 x es solución de la ecuación y '' x 2y '  1 se cumplirá

y1''  x 2y1'  1


Si y2 x es solución de la ecuación y '' x 2y '  1 se cumplirá

y2''  x 2y2'  1

Por lo tanto,

y 1
 y2  '' x 2 y1  y2  '  y1''  x 2y1''  y2''  x 2y1''  1  1  0

Al no cumplir la ecuación diferencial no es solución.

10
Dada la ecuación diferencial x 2y '' xy ' y  0 , se pide:

(a) Comprobar con Matlab que y1  x log x  es solución de la ecuación diferencial.

(b) Determinar a mano la solución general de la ecuación diferencial. Puedes utilizar Matlab
para realizar las integrales.

(c) Calcular con Matlab la solución general y representar la solución que verifica y 1  2 ,

y ' 1  0 .

Solución a)
El código Matlab para ver que es solución puede ser el siguiente que comprueba que cumple la
ecuación diferencial
syms y(x)
y=x*cos(log(x));
dy=diff(y,x);dy2=diff(dy,x);
simplify(x^2*dy2-x*dy+2*y)

Solución b)
Como la ecuación diferencial es lineal homogénea utilizaremos el método de reducción de orden.
Dada y1  x log x  buscamos y2 x   y1 x  v x  linealmente independiente con ella. Se
tendrá que

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

32 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

y2' x   y1' x  v x   y1 x  v ' x 

y2'' x   y1'' x  v x   2y1' x  v ' x   y1 x  v '' x 

Sustituyendo en la ecuación diferencial x 2y2''  xy2'  2y2  0 , es decir,

x 2 y1'' x  v x   2y1' x  v ' x   y1 x  v '' x   x y ' x  v x   y x  v ' x   2y x  v x   0

 1 
  1 1

Simplificando, utilizando que y1 es solución

x 2 2y1' x  v ' x   y1 x  v '' x   x y1 x  v ' x   0

   

v '' x  x 2y1 x   v ' x  2x 2y1' x   xy1 x   0

 

v '' x  2x 2y1' x   xy1 x  2xy1' x   y1 x 

 
v ' x  x 2y1 x  x y1 x 

Utilizando Matlab

v '' x 

  
cos log x   2 sin log x  
v ' x  x cos log x 

Integrando con Matlab

 
 1

log v ' x   2 log cos log x   log(x )  log 
 x cos2 log x 

 

1
v ' x  
x cos log x 
2

u'
Esta integral es inmediata, teniendo en cuenta que  du  tg u  , luego
cos2 u 

 
v x   tg cos log x  
Por lo tanto,

   
y2 x   x cos log x  tg cos log x   x sen log x    
La solución general de la ecuación diferencial es

 
yG x   C 1 x cos log x   C 2x sen log x    C 1,C 2  

Código Matlab
syms y(x)
y1=x*cos(log(x));
dy=diff(y1,x);dy2=diff(dy,x);

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
33

aux=simplify(-(2*x*dy-y1)/(x*y1))
aux1=int(aux)
v1=simplify(exp(aux1))
%Comprobación
y2=x*sin(log(x))
simplify(x^2*diff(y2,x,2)-x*diff(y2,1)+2*y2)

Solución c)
El código para calcular la solución con Matlab es
syms y(x)
dy=diff(y,x);dy2=diff(dy,x);
eqn=x^2*dy2-x*dy+2*y==0
dsolve(eqn)

Para encontrar la solución al problema de valor inicial dado se añadirá al código anterior
cond=[y(1)==1,dy(1)==2];
solu=dsolve(eqn,cond)

Para representar la solución anterior, teniendo en cuenta que la solución está definida para x>0,
se considera un intervalo dentro de este dominio
fplot(solu,[0.2,4])

11
Sabiendo que y1(x )  e x es solución de la siguiente ecuación
xy   (x  1)y   y  0
Encontrar su solución general.
Solución

Si y1  e x , otra solución de la ecuación diferencial homogénea xy   (x  1)y   y  0 que

sea linealmente independiente con ella será de la forma y2  v x e x . Deberá cumplir entonces

xy2''  (x  1)y2'  y2  0

y2'  v ' x e x  v x e x y2''  v '' x e x  2v ' x e x  v x e x

es decir,
   
x v '' x  e x  2v ' x  e x  v x  e x   (x  1) v ' x  e x  v x  e x   v x  e x  0
   
Simplificando

xv '' x   2xv ' x   (x  1)v ' x   0

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

34 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

xv '' x   (x  1)v ' x   0

v '' x  1x 1
  1 log v ' x   log x  x
v ' x  x x

v ' x   xe x v x    xe
x
dx  xe x  e x

Por lo tanto,

 
y2  xe x  e x e x  x  1

La solución general es

yGH  C 1e x  C 2 x  1

12
Verificar que x 2 e x 1 son solución de la ecuacion diferencial x 2y '' 2y  0 . Calcular la
solución general de x 2y '' 2y  0 .
Solución
Se comprueba que cumplen la ecuación
y1  x 2 y1'  2x y1''  2  x 2 y1''  3 y1  0
y 2  x 1 y2'  x 2 y2''  2x 3  x 2 y2''  3 y2  0

Como además son linealmente independientes forman un sistema fundamentalmente de

soluciones:
y1 y2 x2 x 1
W y1, y2     1  2  3  0
y1' y2' 2x x 2

La solución de la edo es yGH  C 1x 2  C 2x 1 C 1, C 2  

EDOS COEFICIENTES CONSTANTES

13
Conociendo las raíces de sus ecuaciones características, formar las ecuaciones
diferenciales lineales homogéneas correspondientes.

(a) r1  1, r2  2 (b) r1  1  r2 (c) r1  3  2i r2  3  2i

Solución
(a) Ecuación característica: r 2  3r  2  0
Ecuación diferencial: y '' 3y ' 2y  0

(b) Ecuación característica: r 2  2r  1  0

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
35

Ecuación diferencial: y '' 2y ' y  0

(c) Ecuación característica: r 2  6r  13  0

Ecuación diferencial: y '' 6y ' 13y  0

14
Encontrar la solución de la ecuación y ''' 3y '' 3y ' y  0 para la cual
y 0  1, y ' 0  2, y '' 0  3

Solución
Para encontrar la solución pedida, se calculan las raíces de la ecuación característica

r  1
3
r 3  3r 2  3r  1  0  0


La raíz es r=1 con multiplicidad 3, luego tres soluciones independentes son e x , xe x , x 2e x  y la
solución general es

y x   C 1e x  C 2xe x  C 3x 2e x

De todas las soluciones de la familia biparamétrica, se deben encontrar las que verifican las tres
condiciones pedidas.

y 0  0  C1  1  y x   e x  C 2xe x  C 3x 2e x

Derivando

y ' x   1  C 2 e x  C 2  2C 3  xe x  C 3x 2e x

Como

y ' 0  2  C2  1  y x   e x  xe x  C 3x 2e x

Derivando de nuevo

y '' x   3  2C 3 e x  1  4C 3  xe x  C 3x 2e x

Imponiendo la última condición, se tiene

y '' 0  3  C3  0  y x   e x  xe x

15
Justificar la certeza o falsedad de la siguiente afirmación:
t t
Las funciones y1  e cos t , y2  e sent son funciones linealmente independientes y
una combinación lineal de ellas es solución de una ecuación diferencial de orden 2 con
coeficientes constantes.
Solución


Considerando que las funciones dadas son de la forma eax cos bx , eax sen bx     con a y b
iguales a 1, basta considerar el polinomio que tiene como solución estas dos raíces complejas:
1  i es . El polinomio es el siguiente

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

36 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

r  1  i  r  1  i   r 2  2r  2
   
que será el polinomio característico de la ecuación diferencial y '' 2y ' 2y  0 que tiene a
esas dos funciones como soluciones linealmente independientes. La ecuación diferencial es de
orden 2 con coeficientes constantes.

16
Dada las funciones y1  e 2x , y2  xe 2x , justificar si constituyen un sistema
fundamental de soluciones de una ecuación diferencial lineal homogénea de coeficientes
constantes. En caso afirmativo, encontrar dicha ecuación diferencial.

Solución
Las dos funciones son linealmente independientes ya que

e 2x xe 2x
W y1, y2   2x 2x 2x
 e 4x  2xe 4x  2xe 4x  e 4x  0
2e e  2xe

Además, estas dos funciones son solución de la ecuación diferencial de coeficientes constantes
que tiene por polinomio característico
2
r  2  r 2  4r  4

que corresponde a la EDO de coeficientes constantes siguiente: y '' 4y ' 4y  0 . La solución

general de esta ecuación diferencial es yGH  C 1e 2x  C 2xe 2x .

17
Escribir, utilizando el método de coeficientes indeterminados, la expresión que tendría
una solución particular y p de la ecuación

y   2y   3y  g(x )
para los siguientes valores de g(x ) :

a) 7 cos 3x b) 5e 3x c) x 2 cos x d) 2xe x

Escribir la forma que tendría la solución particular (no se pide resolver).
Solución
Calculamos en primer lugar, la solución general de la ecuación homogénea asociada

r 2  2r  3  0

que tiene por raíces r  1 , r  3 . La solución general es yGH  C 1e x  C 2e 3x .

En función del término independiente de la ecuación diferencial, escribimos la solución particular

a) 7 cos 3x y p  A cos 3x  Bsen 3x

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
37

b) 5e 3x y p  Axe 3x

c) x 2 cos x  
y p  Ax 2  Bx  C cos x  sen x 

d) 2xe x y p  Ax  B  xe x

18
Encontrar la solución general de las ecuaciones

(a) y '' y ' 6y  20e 2x (b) y '' 10y ' 25y  14e 5x

(c) y '' y  2 cos x (d) y '' 2y ' 5y  25x 2  12

Solución a)
La ecuación característica tiene como raíces -2 y 3. Esto implica que la forma de la solución
particular de la completa es

y p x   Axe 2x

Derivando y sustituyendo en la ecuación se tiene

y p' x   A 2x  1e 2x 

2 x 
 y '' y ' 6y  5Ae 2x  20e 2x  A  4
y p x   A 4x  4e
''


Con lo que la solución general de la completa resulta

y x   C 1e 2x  C 2e 3x  4xe 2x

Solución b)

La ecuación característica es r  5  0 y la solución general de la homogenea es

2

yGH x   C 1e 5x  C 2xe 5x . Esto implica que la forma de la solución particular de la completa
es

y p x   Ax 2e 5x

Derivando y sustituyendo en la ecuación se tiene

 
y p' x   A 2x  5x 2 e 5x 

 y '' 10y ' 25y  2Ae 5x  14e 5x  A  7
5 x 
y p'' x   A 2  20x  25x e
2 

Con lo que la solución general de la completa resulta

y x   C 1e 5x  C 2xe 5x  7x 2e 5x

Solución c)
Las raíces de la ecuación característica son i , por lo que se debe buscar una solución de la
forma

y p x   Ax cos x   Bx sen x 

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

38 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

Derivando y sustituyendo en la ecuación se tiene

y p' x   A  Bx  cos x   B  Ax  sen x  



y p x   2B  Ax  cos x   2A  Bx  sen x 
''

Tras sustituir y simplificar, se tiene

y p''  y p' x   2B cos x  2A sen x   2 cos x   A  0 B 1

Por lo tanto, y p x   x sen x  , siendo entonces la solución general de la completa

y x   x sen x   C 1 cos x   C2 sen x 

Solución d)
La ecuación característica r 2  2r  5  0 , tiene como raíces 1  2i , 1  2i . La forma de la
solución general de la homogénea asociada es


yH x   e x C1 cos 2x   C 2 sen 2x 
Por ser el término independiente un polinomio de grado 2, y no ser el 0 raíz de la ecuación
característica, para la solución particular de la completa se debe ensayar un polinomio de grado
2

y p x   ax 2  bx  c

Derivando y sustituyendo en la ecuación completa se tiene

y p' x   A  Bx  cos x   B  Ax  sen x  



y p x   2B  Ax  cos x   2A  Bx  sen x 
''

Tras sustituir y simplificar, se tiene

y p' x   2ax  b 

  5ax  4a  5b  x  2a  2b  5c   25x  12
2 2

y p x   2a
''



a  5



 b  4


c  2



Por lo tanto,

 
y x   e x C 1 cos 2x   C 2 sen 2x   5x 2  4x  2

19
Encontrar la solución general de la ecuación y '' y ' 6y  e x hallando la solución
particular primero mediante el método de los coeficientes indeterminados y después por
variación de parámetros.
Solución
En primer lugar, se calculan las raíces de la ecuación característica

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
39

r2  r  6  0  r1  2 r2  3

La solución general de la homogénea asociada es

yH x   C 1e 2x  C 2e 3x

COEFICIENTES INDETERMINADOS

Dado que el término independiente es exponencial y que -1 no es raíz de la ecuación

característica, se debe buscar una solución particular de la forma

y p x   Ae x  y p' x   Ae x  y p'' x   Ae x

que sustituyendo en la ecuación proporciona

1 1 x
e x A  A  6A  e x  A  yp  e
4 4
VARIACIÓN DE PARÁMETROS

Se considera y p x   c1 x e 2x  c2 x e 3x , derivando

y p' x   c1' x e 2x  c2' x e 3x  2c1 x e 2x  3c2 x e 3x

Antes de volver a derivar hay que simplificar la primera derivada teniendo en cuenta la
condición adicional de que

c1' x e 2x  c2' x e 3x  0

Luego

y p' x   2c1 x e 2x  3c2 x e 3x

Derivando otra vez

y p'' x   2c1' x e 2x  3c2' x e 3x  4c1 x e 2x  9c2 x e 3x

y sustituyendo en la ecuación se tiene

2c1' x e 2x  3c2' x e 3x  e x

Por lo tanto, las funciones incógnitas se obtienen del sistema

c1' x e 2x  c2' x e 3x  0 


x 
2c1' x e 2x  3c2' x e 3x  e 

de donde

0 e 3x e 2 x 0
x 2 x
e 3e x
e x 2e e x e 4 x
c x  
'
1
 c x  
'
2

5e x 5 5e x 5
con lo que

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

40 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

e x e 4 x
c1 x   c2 x   
5 20

e x 2x e 4x 3x 1 x
y p x   e  e  e
5 20 4
La solución general de la completa es

1 x
y x   e  C 1e 2x  C 2e 3x
4

 
y x   e x C 1 cos 2x   C 2 sen 2x   5x 2  4x  2

20
Encontrar la solución particular de la ecuación y '' 4y  tg(2x ) .

Solución
En este caso el término independiente no es adecuado para utilizar el método de coeficientes
indeterminados por lo que se utiliza el método de variación de constantes. Las raíces de la
ecuación característica son 2i , por lo que se busca una solución particular de la forma

y p x   c1 x  cos 2x   c2 sen 2x 

Derivando y teniendo en cuenta la condición adicional

c1' x  cos 2x   c2' x  sen 2x   0

se tiene

y p' x   c1' x  cos 2x   c2' sen 2x   2c1 x  sen 2x   2c2 x  cos 2x  

y p' x   2c1 x  sen 2x   2c2 x  cos 2x 

Derivando de nuevo

y p'' x   2c1' x  sen 2x   2c2' x  cos 2x   4c1 x  cos 2x   4c2 x  sen 2x 

Sustituyendo en la ecuación

y p'' x   4y p x   2c1' x  sen 2x   2c2' x  cos 2x   tg 2x 

El sistema que hay que resolver es

1 
c x  cos 2x   c x  sen 2x   0
'
1
'
2



c1
'
x  
2
tg 2x  sen 2x  

 
2c1 x  sen 2x   2c2 x  cos 2x   tg 2x 
' '
1 
 c2' x   sen 2x  
2 
Integrando

1 sen 2x  1 1  cos 2x 

2 2
1
c1 x      
2  2  cos 2x  2  cos 2x 
tg 2x sen 2x dx  dx  dx 

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
41

  
1  1  sen 2x  1
  
1  1    2    
2 
  cos 2x  dx log 
  2 sen 2x
 cos 2x    cos 2x   
  
  

1 1
c2 x    sen 2x dx  cos 2x 
2 4
Finalmente resulta

1
y p x   c1 x  cos 2x   c2 x  sen 2x  
4  
cos 2x  log sec 2x   tg 2x  
 
21 e 2 x
Resolver la ecuación diferencial y '' 4y ' 4y  .
x2
Solución
Se trata de una ecuación diferencial de coeficientes constantes. El polinomio característico
asociado a la ecuación homogénea es
2
r 2  4r  4  r  2  0  yGH  C 1e 2x  C 2xe 2x

Para obtener la solución particular consideramos

y p  C 1 x e 2x  C 2 x  xe 2x cumpliendo C '1 x e 2x  C '2 x  xe 2x  0

Sustituyendo yp en la ecuación e imponiendo esta condición obtenemos

e 2 x

2C '1 x e 2x  C '2 x  e 2x  2xe 2x   x2

 
Calculamos C '1 x ,C '2 x del sistema

C '1 x e 2x  C '2 x  xe 2x  0 




2 x  
e 

2C '1 x e 2x  C '2 x  e 2x  2xe 2x  
x2 





C '1 x   C '2 x  x  0 


1 
2C '1 x   C '2 x 1  2x   
x 2 

 
Despejando de la primera ecuación C '1 x  C '2 x x , sustituyendo en la segunda

1 1 1
2xC '2 x   C '2 x 1  2x    C '2 x    C 2 x   
2 2 x
x x
1 1
C '1 x    x   C 1 x    log x
x2 x
La solución particular es

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

42 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

y p   log x e 2x  e 2x

La solución general es

yG  C 1e 2x  C 2xe 2x  log x e 2x

22 9x
Resolver la ecuación diferencial y '' 9y  .
e 3x
Solución

9x
La ecuación diferencial de coeficientes constantes y '' 9y  tiene como polinomio
e 3x
característico

r 2  9  r  3r  3  0

La solución general es yGH  C 1e 3x  C 2e 3x . Para encontrar la solución particular de la completa

consideramos

y p  x Ax  B e 3x 
y p  Ax 2  Bx e 3x 
y determinamos los coeficientes A y B sabiendo que es solución de la ecuación diferencial:
9x
y p''  9y p  3x  9xe 3x . Como se cumple
e


y p'  2Ax  B e 3x  3 Ax 2  Bx e 3x 

y p''  2Ae 3x  6 2Ax  B e 3x  9 Ax 2  Bx e 3x 
sustituyendo en la ecuación diferencial

  
2Ae 3x  6 2Ax  B e 3x  9 Ax 2  Bx e 3x  9 Ax 2  Bx e 3x  9xe 3x 
se tiene el siguiente sistema
12A  9
2A  6B  0
cuya solución es A  3 / 4, b  1 / 4 . Por lo tanto, la solución particular es

 3 1 
y p   x 2  x e 3x
 4 4 

y la solución general es
 3 1 
yG  C 1e 3x  C 2e 3x   x 2  x e 3x
 4 4 

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
43

23
Escribir la expresión de la solución particular de la ecuación diferencial
y '' 5y ' 4y  12x  5e 4x utilizando el método de coeficientes indeterminados. Nota: No se
pide calcular la solución particular.
Solución
Para encontrar la expresión de la solución particular, calculamos en primer lugar la solución de la
ecuación diferencial homogénea asociada y '' 5y ' 4y  0 . El polinomio característico es
r 2  5r  4  r  4r  1  0
La solución general de la homogénea asociada es
yGH x   C 1e x  C 2e 4x C 1, C 2  
En consecuencia, la solución particular será d ela forma

y p x   x Ax  B e 4x  Ax 2  Bx e 4x  A, B  
Aunque no pide resolverla para encontrar esta solución se tendrá en cuenta que

y p x   Ax 2  Bx e 4x 
  
y p' x   2Ax  B e 4x  4 Ax 2  Bx e 4x  4Ax 2  2A  4B  x  B e 4x 
  
y p'' x   8Ax  2A  4B  e 4x  4 4Ax 2  2A  4B  x  B e 4x  

 2 8Ax  8A  8B  x  A  4B e
2
 4x

 
Como y p''  5y p'  4y p  12x  5 e 4x , se tendrá
6Ax  2A  3B  e 4x  12x  5e 4x
 
Se debe resolver el siguiente sistema
6A  12 A  2
2A  3B  5  B  3
 
La solución es entonces y p  2x 2  3x e 4x y la solución general de la ecuación diferencial es


y  C 1e x  C 2e 4x  2x 2  3x e 4x  C 1, C 2  

24 et
Resolver la siguiente ecuación diferencial y '' 2y ' y 
t2
Solución
Como es una edo de segundo orden de coeficientes constantes, calculamos el polinomio
característico para obtener la solución general de la ecuación diferencial homogénea asociada.
Dado que

r 2  2r  1  0  r  1 doble
la solución general de la homogénea es
yGH  C 1et  C 2tet C 1, C 2  
Para obtener la solución particular de la completa utilizamos el método de variación de
constantes

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

44 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

y p  C 1 t et  C 2 t  tet C 1' t et  C 2' t  tet  0

Las funciones C1 y C2 deben cumplir

C 1' t et  C 2' t  tet  0 

t
e 

C 1' t et  C 2' t  et  tet   
t 2 
Resolviendo el sistema
0 tet
et e 2t
et 1  t 
2 1
C 1' t   t  t   C 1 t   log t
et tet e 1  t   te
2t 2t t
et et 1  t 

et 0
et e 2t
et
C 2' t   t2 2
 t 
1
 C 2 t   
1
et tet e 2t t2 t
et et 1  t 

La solución particular es
1
y p  log t et  tet  log t et  et
t
La solución de la edo dada es

y t   C 1et  C 2tet  log t et C 1, C 2  

Observar que la expresión et de la solución particular no se ha añadido al formar parte de la
expresión C 1et .

APLICACIONES

25
Aplicación. Oscilaciones mecánicas
Supongamos un muelle de masa m que colgamos de un resorte y que se agrega un sistema
mecánico (amortiguador) que tiene el efecto de reducir la velocidad cuando el sistema se
encuentra vibrando. Supongamos que f t  es una fuerza externa que actúa sobre el sistema.

Teniendo en cuenta la ley de Hooke, el movimiento del cuerpo viene dado

mx '' t   Fm  Fc  Fe
Donde Fm es la fuerza del muelle siendo k una constante que depende de la rigidez del muelle,
Fc es la fuerza que ejerce el amortiguador y Fe  f t  es la fuerza externa. La ecuación
diferencial de este sistema es

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
45

mx '' t   c x ' t   k x t   f t 

[Link]

Ejemplo: Un cuerpo oscilando sobre su posición de equilibrio (x=0) se mueve de acuerdo a la

ecuación diferencial x '' t   16x t   24 cos 4t  . Si en el instante inicial se encuentra en su
posición de equilibrio y con velocidad nula, ¿en qué posición se encontrará al cabo de 10
segundos? Representa la solución de la ecuación diferencial.
Solución con Matlab
syms x(t)
eqn= diff(x,t,2)+16*x==24*cos(4*t);
Dx=diff(x,t);
cond=[x(0)==0,Dx(0)==0]
sol(t)=dsolve(eqn,cond)
%En el instante t=10
sol(10)
ezplot(sol(t),[0,10])
hold on
plot(0,0,'o')
plot(10,sol(10),'*')
hold off

26
Aplicación. Oscilaciones mecánicas
Aplicación. Circuito eléctrico LRC
Una aplicación clásica de las ecuaciones diferenciales se presenta en el estudio de un circuito
eléctrico que consiste en resistores, inductores y capacitores, al cual se aplica una fuerza
electromotriz.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

46 T5  ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN Y SISTEMAS DE PRIMER ORDEN

En este caso una aplicación de las leyes de Kirchoff conduce a la ecuación

d 2q dq 1
L 2
R  q  E (t )
dt dt C
donde L es la inductancia, R la resistencia, C la capacitancia, E (t ) la fuerza electromotriz,
q(t ) la carga y t el tiempo. Esta ecuación diferencial constituye un ejemplo de las ecuaciones
diferenciales que se estudiarán en este tema: las ecuaciones lineales de segundo orden.
Teniendo en cuenta que la intensidad se define como la derivada de la carga, la ecuación
diferencial en términos de la intensidad es

d 2i di 1
L 2
 R  i  E '(t )
dt dt C
Ejemplo: Determina la carga q(t ) en un circuito en serie LRC , cuando:

L  0.25 henrios (h ) R  10 ohmios ()

C  0.001 faradios(f) E (t )  200 sen 40t voltios (V )

q(0)  3 culombios (C) i(0)  q (0)  0 amperios (A)

Solución con Matlab

syms q(t)
L=0.25;R=10;C=0.001;
eqn= L*diff(q,t,2)+R*diff(q,t)+q/C==200*sin(40*t);
Dq=diff(q,t);
cond=[q(0)==3,Dq(0)==0]
sol(t)=dsolve(eqn,cond)
ezplot(sol(t),[0,10])

 41 35  3 2
q(t )  e 20t  cos 60t  sen 60t   sen 40t  cos 40t
13 39  13 13

Puedes ver más ejercicios resueltos sobre ecuaciones diferenciales de segundo orden en
la página de Giematic UC
[Link]

Prof. Elena E. Alvarez Saiz Universidad de Cantabria