Ecuaciones Diferenciales de Segundo Orden

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 29

Ecuaciones Diferenciales

Ordinarias de Segundo
Orden
Con Coeficientes Constantes

Integrantes:
- Cholán
- Gálvez Silva, Oscar
- Mendoza Vásquez, Luis
- Ludeña
- Ludeña

Profesor: Sevillano Castro, Rodolfo


Ecuaciones Diferenciales Ordinarias

Indice

I. Introducción……………………………………….02

II. Objetivos…………………………………………..03

III. Fundamento Teórico……………………………....04

III.I Ecuaciones Diferenciales de segundo orden


Homogéneas………………………………………….04

III.II Ecuaciones Diferenciales de segundo orden No


Homogéneas………………………………………….10

IV. Ejercicios de aplicación……………………………19

1
Ecuaciones Diferenciales Ordinarias

I. INTRODUCCIÓN

2
Ecuaciones Diferenciales Ordinarias

Ejemplo.- Si f1 ( x )  cos x ; f 2 ( x )  sin x entonces

II. OBJETIVOS.
Principal
o Estudiar las ecuaciones diferenciales ordinarias lineales de
segundo orden con coeficientes constantes homogéneo y no
homogéneo.
Secundario
o Aprender a resolver ejercicios y problemas usando las
ecuaciones diferenciales
o Aprender los métodos existentes para resolver las ecuaciones
diferenciales ordinarias de segundo orden homogénea y no
homogénea.
o Aprender las aplicaciones de este tipo de ecuaciones
diferenciales.

3
Ecuaciones Diferenciales Ordinarias

III. FUNDAMENTO TEÓRICO


 DEFINICIÓN:

Una ecuación diferencial es una ecuación en la que intervienen


derivadas de una o más funciones desconocidas. Dependiendo del
número de variables independientes respecto de las que se deriva, las
ecuaciones diferenciales se dividen en:

o Ecuaciones diferenciales ordinarias: aquellas que contienen


derivadas respecto a una sola variable independiente.
o Ecuaciones en derivadas parciales: aquellas que contienen
derivadas respecto a dos o más variables.

III.I ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN


HOMOGENEAS:

Una ecuación diferencial de Segundo Orden con coeficientes


constantes homogénea es de la forma:
ay ´ ´ +by ´ +cy =0
La función "y", solución general de la ecuación diferencial anterior, es
de la forma y ( x ) =k e rx. Donde "k" es una constante que da la
generalidad de la solución.
Entonces el objetivo ahora será hallar el valor de r.
Bien, de la solución general tenemos: y ´=kr erx
y ´ ´=k r 2 e rx
Reemplazando en ay ´ ´ +by ´ ´ +cy =0 tenemos:
ak r e +bkr e rx +ck erx =0
2 rx

k e rx [ ar 2 +br + c ] =0
Ahora bien, k≠ 0porque si no tuviéramos la solución trivial y como
también e rx ≠ 0, entonces a r 2 +br + c=0. A esta expresión se la
denomina Ecuación Auxiliar y es útil para hallar r.
Observe que la ecuación auxiliar es una ecuación cuadrática cuyas
raíces se las puede determinar empleando la formula general
−b ± √b 2−4 ac
r 1 , r 2=
2a

AQUÍ SE PRESENTAN TRES CASOS:

 CASO I
Discriminante positivo[ b 2−4 ac> 0 ]. Entonces r 1y r 2son raíces reales y
diferentes. En este caso se dice que existen dos soluciones
fundamentales
y 1 ( x )=k 1 erx
y 2 ( x ) =k 2 erx

La solución General estaría dada por la combinación lineal de las


soluciones fundamentales

4
Ecuaciones Diferenciales Ordinarias

y ( x ) =k 1 e rx +k 2 e rx

 CASO II
Discriminante cero[ b 2−4 ac=0 ] . Entonces r 1 y r 2 son raíces reales e
iguales.

 CASO III
Discriminante negativo[b2−4 ac <0]. Entonces r 1= λ+ μiy r 2= λ−μi son
raíces complejas conjugadas
Reemplazando en y ( x ) =c 1 e rx + c 2 e rx tenemos:
y ( x ) =C1 e ( λ+ μi ) x +C 2 e ( λ−μi ) x
y ( x ) =C1 e λx e μix +C 2 e λx e− μix
y ( x ) =e λx [ C1 e μix +C 2 e−μix ]
Como e iμx=cos μx +i sin μx y e−iμx =cos μx+ isin μx
Reemplazando tenemos:
y (x )=e λx [ C1 ( cos μx +isin μx )+ C2 ( cos μx−isin μx ) ]
y (x )=e λx ¿

Por lo tanto la solución sería


y (x )=e λx [ k 1 sin ( μx ) + k 2 cos ( μx ) ]

DEPENDENCIA LINEAL. WRONSKIANO:

El método usado para obtener todas las soluciones de las ecuaciones


diferenciales lineales se basa en el hecho de que, al igual que los vectores
ordinarios del plano o del espacio, la solución de la ecuación homogénea:
d2 y dy
2
+ P ( x ) +Q ( x ) y=0….. (1)
dx dx
Se pueden considerar como vectores en un cierto espacio vectorial ℇ de
funciones. En efecto, como los vectores del plano, las funciones se
pueden sumar entre sí y multiplicar por números, reales o complejos:
( f + g )( x )=f ( x )+ g ( x ) . y ( λf ) ( x ) =λf ( x) y, en el caso particular de las
soluciones de (1), tenemos las siguientes propiedades:
 Proposición 1:
La suma de dos soluciones de la ecuación (1) en un cierto
intervalo es también solución en el mismo intervalo. Lo mismo
ocurre con el producto de una solución de (1) por cualquier
número, real o complejo.
 Proposición 2:
El conjunto de todas las soluciones de la ecuación diferencial (1)
en un cierto intervalo tiene estructura de espacio vectorial. Cada
solución es pues un vector de dicho espacio.

 Proposición 3:

5
Ecuaciones Diferenciales Ordinarias

Dos soluciones de la ecuación diferencial (1) en un intervalo I son


linealmente independientes en I si el wronskiano de dichas
funciones es distinto de cero en un punto cualquiera de I, en cuyo
caso es distinto de cero en todos los puntos de I.
Esta proposición es equivalente a la siguiente:

 Proposición 4:
Dos soluciones de la ecuación diferencial (1) en un intervalo I son
linealmente dependientes en I si el wronskiano de dichas funciones
se anula en un punto cualquiera de I, en cuyo caso se anula en
todos los puntos de I.
 Proposición 5:
Este resultado caracteriza a todas las soluciones de la ecuación
diferencial lineal de segundo orden completo:
d2 y dy
2
+ P ( x ) +Q ( x ) y=R (x)
dx dx
Se puede escribir como la suma de la solución general de la
ecuación homogénea (1) y de una solución particular de la
completa.

 Proposición 6:

Sea y 1 ( x)una solución, distinta de la trivial (solución nula), de la


ecuación diferencial lineal homogénea (1) en un intervalo [a,b].
Entonces existe una segunda solución y 2 ( x ) =v ( x )
y 1 ( x ) de ( 1 ) en el mismointervalo , donde
1
v ( x )=∫ 2 e−∫ P ( x ) dx dx ,
y 1 (x)
Y las dos soluciones y 1 e y 2son linealmente independientes en [a,b]

EJEMPLOS DE ECUACIONES DIFERENCIALES DE SEGUNDO


ORDEN HOMOGÉNEAS:

o Ejemplo 1:
y ´ ´ + y ´ + y=0

Solución:
Polinomio característico a la ecuación diferencial es:
−1 √ 3
y ´ ´ + y ´ + y=0 es P ( r )=r 2 +r +1=0 , de donde r 1= + i ,
2 2
−1 √ 3
r 2= − i
2 2
−x −x
3 3
Luego el sistema de solución es:e 2
cos √ x, e 2
sen √ xy la solución
2 2
general es:

6
Ecuaciones Diferenciales Ordinarias

−x −x
3 3
C1e 2
cos √ x+ C2 e 2
sen √ x
2 2

o Ejemplo 2:
y ´ ´ +6 y ´ +13 y=0 ; y ( 0 )=1 ; y ´ ( 0 )=1
Solución
Ecuación auxiliar sería:
−6 ± √36−4 (1 ) ( 13 )
r 1 , r 2=
2 ( 1)
−6 ± √−16
r 1 , r 2= √ −1=i
2 ( 1)
Hallando las raíces tenemos:
−6 ± √−16 √ −1
r 1 , r 2=
2 (1)
−6 ± 4 i
r 1 , r 2=
2 (1 )
r 1=−3+2i ∨r 2=−3−2 i
En este caso λ=−3 y o μ=2, por tanto la solución general sería:
y (0)=e−3( 0) [ k 1 sin ( 2 ( 0 ) ) +k 2 cos ( 2 ( 0 ) ) ]
Como y (0)¿ 1
Entonces:
y (0)=e−3( 0) [ k 1 sin ( 2 ( 0 ) ) +k 2 cos ( 2 ( 0 ) ) ]
1=( 1 ) [ k 1 ( 0 ) +k 2 ( 1 ) ]
1=k 2
Como y ´ ( 0 )=1
Entonces:
y ´ ( 0 )=e−3 ( x ) [ 2 k 1 cos ( 2 x )−2 k 2 sen ( 2 x ) ]−3 e−3 x [ k 1 sen (2 x ) +k 2 cos ( 2 x ) ]

y ´ ( 0 )=e−3 ( 0) [ 2k 1 cos ( 0 )−2 k 2 sen ( 0 ) ]−3 e−30 [ k 1 sen ( 0 ) +k 2 cos ( 0 ) ]


1=1 k 1 −3 k 2
Resolviendo simultáneamente:
1 3
+ k =k
2 2 2 1
1 3
+ 1=k 1
2 2
k 1=2
Por lo tanto la solución sería:
y ( x ) =e−3 x [ 2sin ( 2 x )+ cos ( 2 x ) ]

7
Ecuaciones Diferenciales Ordinarias

o Ejemplo 3: Hallar la solución general de 4 y''  4 y'  17 y con las


condiciones iniciales y(0) = -1 e y’(0) = 2

o Ejemplo 4: Resolver la ecuación diferencial y” +6y’ +9y =0. Con las


condiciones iniciales dadas y(0) = 2 e y’(0) = 1

8
Ecuaciones Diferenciales Ordinarias

o Ejemplo 5: Resolver la siguiente ecuación diferencial y” + 6’ +


12y=0

o Ejemplo 6: Solucionar la siguiente ecuación diferencial y” – 5y’ =0


donde y(0) = 1 y’(0) = 2

9
Ecuaciones Diferenciales Ordinarias

III.II ECUACIONES DIFERENCIALES NO HOMOGENEAS DE SEGUNDO


ORDEN:

Tiene la siguiente forma


2
d y dy
a 2 2 +a1 +a 0 y=F (x )
d x dx

Donde F(x) es continua en algún intervalo [a, b] y a 2,a 1,a 0 son constantes con
a 2≠0

Para ello será conveniente usar OPERADORES tales como:


D2 D ASI
d2 y
D2 y ;es lo mismo que: 2
d x
2
d y
D y es lo mismo que 2
d x
ENTONCES
d2 y dy
a 2 2 +a1 +a 0 y=F ( x ) = a 2D2 y +a 1D y +a 0 y =F (x)
d x dx

( a 2D2 +a 1D +a 0 ¿ y =F (x)

10
Ecuaciones Diferenciales Ordinarias

Ø(D) y =F
EJEMPLO
La ecuación diferencial Ø(D) y =F, donde F=4e x -cos(x)
Y Ø(D)= xD 2 -2D +3x2
Es una manera abreviada de escribir:
d2 y dy 2
x 2 −2 + 3 x y=¿4e x -cos(x)
dx dx
OPERADOR ANULADOR
Del cálculo diferencial se sabe que se puede derivar varias veces un polinomio
hasta que la derivada n-ésima sea cero a partir de algún entero n
Sea: y= f(x), una función que tiene al menos n derivada. Si
(a nDn +a n−1 Dn−1 +…+a 1D1 +a 0 )f(x)=0
Entonces se dice que el operador diferencial
a nDn +a n−1 Dn−1 +…+a 1D1 +a 0 anula a f
SE DESTACAN LOS SIGUIENTES OPERADORES ANULADORES
Dn anula a la funciones 1, x, x2,…,xn-1
Y=c0+ c1x+ c2x2+…+ cn-1xn-1
Ejemplo
D5(x4-3x2+2x+1)=0
El operador diferencial (D-ᾳ)n anula a:
e ᾳx , xe ᾳx , …, x n−1 e ᾳx
Y por lo tanto anula a cualquier combinación lineal de las mismas:
Y= c0e ᾳx + c1xe ᾳx + c2x2e ᾳx +…+ cn-1xn-1e ᾳx
Ejemplo
(D-3)e3x =0
Operador anulador asociado con polinomios característicos de raíces complejas:
El operador diferencial
[D2-2αD + (α2+β2)]
Anula a cada una de las funciones:
e ∝ xcos βx, xe ∝ xcos βx,x2e ∝ xcos βx,…., xn-1e ∝ xcos βx
e ∝ xsen βx, xe ∝ xsen βx,x2e ∝ xsen βx,…., xn-1e ∝ xsen βx
Ejemplo
Encontrar el operador anulador para e 2 x (sen3x-cos3x)
Solución
La expresión e 2 x (sen3x-cos3x) se puede escribir como:
e 2 x sen3x-e 2 x cos3x)
Aplicando el operador anulador a cada término se tiene:
(D2-4D+13)e 2 x cos3x=0 Y (D2-4D+13)e 2 x sen3x=0

OBTENCION DE LA SOLUCION GENERAL DE LA ECUACION


DIFERENCIAL LINEAL DE SEGUNDO ORDEN NO HOMOGENEA
a 2 y " +a 1 y ' + a 0 y = F(x) la solución general es:
Y=C1y1(x)+ C2y2(x) + yp
En donde:
Yc= C1y1(x)+ C2y2(x)
Yc=solución complementaria
YP=solución particular
CÁLCULO DE LA SOLUCION PARTICULAR
Seguimos el siguiente procedimiento:

11
Ecuaciones Diferenciales Ordinarias

1. Dividir toda la ecuación por a 2 para convertirla en:

a1 a F (x )
y" + y' + + 0 y =
a2 a2 a2
2. Hallar la solución complementaria

Yc= C1y1(x)+ C2y2(x)


3. Usar el operador anulador para la función F(x), aquí es donde se aplica el método de
los coeficientes indeterminados
4. Si no es posible hallar un operador anulador para F(x), entonces hallamos el
wronskiano

W( y 1(x), y 2(x))
y 2 f ( x)
5. Halar u1'(x)=
W ¿¿

y 1 f ( x)
u2'(x)=
W ¿¿
F (x)
Donde: f (x)=
a2
6. Determinar u1(x)=∫ u 1 '( x ¿) ¿dx y u2(x)=∫ u 2 '( x¿) ¿dx
7. La solución particular será de la forma YP=u1(x) y 1(x)+ u2(x) y 2(x)

Otro concepto

Una ecuación diferencial lineal de segundo no homogénea de coeficiente


constante, es de la forma: y” + f(x) y’ + g(x) y = r(x), donde f(x) y g(x) son
constantes (I)
La diferencia con las anteriores ecuaciones estudiadas estriba en que está
igualada a una función de la variable independiente x. Esto nos sugiere una
relación entre:
y” + f(x)y’ + g(x)y = 0 y y” + f(x)y’ + g(x)y =
r(x)
Llamaremos yh a la solución general de la ecuación homogénea correspondiente
y yp a una solución particular de la no homogénea que podamos encontrar de
alguna manera; entonces se puede establecer el siguiente teorema.
Teorema. Si yh es la solución general de y” + f(x)y’ + g(x)y = 0 y y p es
cualquier solución particular de (1)
Y = yh + yp es la solución general de (1).
Demostración: supongamos que y = yh + yp es solución de (1):
y’ = yh‘ + yp ‘
y” = yh” + yp” , sustituyendo en (1):
yh” + yp” + f(x)( yh‘ + yp ‘ ) + g(x)( yh + yp) = r(x)
agrupando: (yh” + f(x)yh+g(x)yh´) + (yp” + f(x)yp+g(x)yp´) = r(x)
Como yh es solución de la homogénea, el primer paréntesis se hace cero, y como
ypes la solución de la no homogénea (1), el segundo paréntesis se convierte en
r(x),
0+ r(x) = r(x)
Por tanto, si satisface y = yh + yp a la ecuación (1).

12
Ecuaciones Diferenciales Ordinarias

Como la solución yh por los métodos anteriores, el problema se reduce entonces


a encontrar la solución para yp resolver las ecuaciones no homogéneas los
métodos para encontrar yp son: coeficientes indeterminados y variación de
parámetros.
El método de variación de parámetros, llamado también método general supone
el cambio de constantes c1 y c2 de la solución yp, por funciones de x. El método
de coeficientes indeterminados es más sencillo y se usa para ciertos tipos de
función r(x).
METODO DE COEFICIENTES INDETERMINADOS PARA OBTENER
yp

Se usa para tres formas de r(x): r(x) = polinomio, r(x) = exponencial, r(x) =
función trigonométrica o combinaciones de ella, que puede resumirse, en forma
general, de la siguiente manera:
r(x) = eαx [Pm(x) cosβx + Qn(x) senβx]
donde λ =α±i β es raíz de la ecuación auxiliar y P m(x) y Qn(x) son polinomios
de grado m y n, respectivamente. Se busca una solución particular y p de la
forma:
yp = xz eαx [pk(x)cosβx + qk(x)senβx]
Donde k = máx(m,n), pk(x) y qk(x) son polinomios de grado k, cuyos
coeficientes están indeterminados, y z es la multiplicidad de la raízλ = α±i β de
la ecuación auxiliar. La forma yp se puede resumir en el siguiente cuando

Forma de r(x) Raíces de la ecuación Forma de yp para


auxiliar k = máx(m,n)
x
α no es raíz Pm(x)e α
x
Pm(x)e α Α es raíz repetida z veces x

(de orden z) xnpm(x)e α

Forma de r(x) Raíces de la ecuación Forma de yp para


auxiliar k = máx(m,n)
λt ≠ 0, i = 1, 2, … z Pm(x)
Pm(x)
Alguna λt = 0 xzpm(x)

Forma de r(x) Raíces de la ecuación Forma de yp para


auxiliar k = máx(m,n)
Pm(x)cosβx + Qn(x)senβx ± iβ no son raíces ph((x)cosβx + qk(x)senβx)
± iβ son raíces de orden z xz(pk(x)cosβx +
qh(x)senβx)

Forma de r(x) Raíces de la ecuación Forma de yp para


auxiliar k = máx(m,n)
x x

e α [Pm(x)cosβx+Qnsenβx] α± iβ no son raíces α


e (ph(x)cosβx+qksenβx)
α± iβ son raíces de x
xze α (ph(x)cosβx+qksenβx)
orden z

13
Ecuaciones Diferenciales Ordinarias

 EJEMPLOS DE ECUACIONES DIFERENCIALES DE SEGUNDO


ORDEN NO HOMOGÉNEAS:

o Ejemplo1:
Sea y } +a {y} ^ {'} +by=R(x ¿
Solución:
La solución general es de la forma y(t) = y c + yp
La solución complementaria satisface la ecuación homogénea
y }c +4 {y} rsub {c} rsup {'} +9 {y} rsub {c} = ¿
La ecuación auxiliar esr 2 + 4 r +9 . Hallando las raíces tenemos

−4 ± √ 16−4 (9)
r 1 , r 2=
2
−4 ± √−20
r 1 , r 2=
2
−4 ± √ 20 √−1
r 1 , r 2=
2
−4 ± √ ( 5 ) 4 √−1
r 1 , r 2=
2
−4 ± 2 √ 5 i
r 1 , r 2=
2
−4+2 √ 5i
r 1= → r 1=−2+ √ 5 i
2
−4−2 √ 5 i
r 2= →r 2=−2− √5 i
2

Por tanto:

y c ( x )=e2 x [ k 1 sen ( √ 5 x ) +k 2 cos ( √ 5 x ) ]

o Ejemplo 2:
dy 2 dy 2
2
+3 +2 y=4 x
dx dx
Solución:
Sea p ( m )=m 2+ 3 m=0 la ecuación característica donde m 1=−1; m 2=−2
luego la solución general de la ecuación diferencial homogénea o solución
complementaria es:
y c =c1 e +¿c e ¿
−x −2x
2

Usando el operador diferencial se tiene ( D¿¿ 2+3 D+2) y=x 2 ¿ como


3 3 2
D 3 x2 =0 se tiene D (D¿¿ 2+3 D+2) y=4 D x =0 ¿. Luego la ecuación
auxiliar asociada es
m3 (m¿¿ 2+3 m+2)=0¿ Por lo tanto la solución general es:
y=c1 e− x + c2 e−2 x +c 3 +c 4 x+ c 5 x 2
solucion deSolución particular de la ecuación no homogénea.
la ecuacion
homogenia

14
Ecuaciones Diferenciales Ordinarias

A los coeficientes c 3 , c 4 , c 5 se le denomine indeterminados y se procede a


determinarlos en la siguiente formula:
Se supone que la solución particular tiene la forma y p=c3 + c 4 x +c 5 x2
como debe ser solución de la ecuación diferencial dada, se tiene:
y 'p' +3 yp+ 2 yp=2 c 3 +3 c 4 +2 c 5+ ( 2 c 4 +6 c 5 ) x +2 c 5 x 2=4 x 2
igualando coefientes:
2 c 3+ 3 c 4 +2 c5 =02 c 4 +6 c 5=02 c 5=4
Resulta c 3=0, c 4 =−6 y c 5=−2
Entonces y p=7−6 x +2 x2
Luego la solución general viene expresado por y= y c+¿ y ¿p

Entonces y=c1 e− x + c2 e−2 x +7−6 x+2 x 2

o Ejemplo 3:
y ' ' −5 y ' +6 y =x e 3 x

Solución.

Sea m 2−5 m+6=0, la ecuación característica donde m1=2 , m2=3 luego la


solución general de la ecuación diferencial homogénea o solución
complementaria es:

y c =c1 e 2 x + c 2 e 3 x

Ahora aplicamos el operador anulador a nuestra ecuación se tiene

( D−3 )2( D−3)( D−2) y=( D−3 )2 x e3 x =0


Luego la solución general es:

y g=c 1 e2 x +c 2 e 3 x + c 3 x e3 x + c 4 x 2 e 3 x

15
Ecuaciones Diferenciales Ordinarias

Donde una solución particular es:

y p=c3 x e3 x + c 4 x 2 e3 x

Luego reemplazando en la ecuación original para hallar la constante se


tiene

(e 3 x ( c 3 x+ c 4 x 2 ) )-(5( {c} rsub {3} x+ {c} rsub {4} {x} ^ {2} ))'+( {e} ^ {3x} ( {c} rsub {3} x

De la cual obtenemos: c 3=−1 / 4


c 4 =1/2
Así
−1 3 x x 2 3 x
y p= xe + e
4 2
Finalmente:
3x 3x 1 3 x x2 3 x
y=c1 e +c 2 e − x e + e
4 2

o Ejemplo 4:
y -4y=x {e} ^ {x} +cos2

Solución:
Sea m2−4=0 la ecuación característica donde m 1=2 , m 2=−2 luego la
solución general de la ecuación diferencial homogénea a solución
complementaria es:
y=c1 e 2 x +c 2 e−2 x

Observamos que un operador anulador para x e x es ( D−1 )2 y para cos2x es


2
D2 +4 luego el operador anulador para ( x e x + 2 x ¿ será ( D−1 ) ( D2 +4 )
Ahora aplicando el operador anulador a nuestra ecuación se tiene:
( D−1 )2 ( D 2 +4 ) ( D2−4) y=( D−1 )2 ( D2 +4 ) ( x e x + cos 2 x ¿=0
Por lo tanto la solución general es:
y c =c1 e 2 x + c 2 e−2 x + y =c 3 e x +c 4 x e x + y=c 5 cos 2 x+ c6 sen 2 x

16
Ecuaciones Diferenciales Ordinarias

Donde una solución particular es:


y p=c3 e x +c 4 x e x + c 5 cos 2 x +c 6 sen 2 x
Luego reemplazando en la ecuación original para hallar las constantes de
la solución particular se tiene:
(c 3 e x + c 4 x e x +c 5 cos 2 x+ c6 sen 2 x ¿ -4(
c 3 e x + c 4 x e x +c 5 cos 2 x+ c6 sen 2 x ¿=x e x +cos 2 x

De donde
c 3=−2/9
c 4 =−1/3
c 5=−1 /8
c 6=0
Así la solución particular es
−2 x 1 1
y p= e − c 4 x e x − cos 2 x
9 3 8
Finalmente la solución general es:
−2 x 1 1
y g=c 1 e2 x +c 2 e−2 x − e − c 4 x e x − cos 2 x
9 3 8

o Ejemplo 5:
Y” – 5y’ +6y = xe3x

17
Ecuaciones Diferenciales Ordinarias

o Ejemplo 6:

Sea y”+4 y’+9y= x2 +3x Hallar la solución General

SOLUCIÓN:
la solución general de la forma y(t) = yc + yp
Primero hallemos yc
La solución complementaria satisface la ecuación
homogénea y”c +4y’c +9yc =0
La ecuación auxiliar es r2 + 4r + 9 = 0 . Hallando las raíces
tenemos.
r1,r2¿−4 ±
√(16−¿ 4 (9)) ¿
2
−4 ± √−20
r1,r2¿
2
−4 ± √ 20 √−1
r1,r2¿
2
−4 ± √ ( 5 ) 4 √ −1
r1,r2 ¿
2
−4 ± √ 20 √−1
r1,r2¿
2

18
Ecuaciones Diferenciales Ordinarias

−4 ± 2 √5 i
r1,r2¿
2
−4 ± √ 20 √−1
r1 ¿ »r1¿−2+ √ 5 i
2
−4 ± √ 20 √−1
r2¿ »r2¿−2−√ 5 i
2
Por tanto yc (x) = e-2x[k1sen(√5x) + k2 cos (√5x)]
Segundo, hallamos yp

Como g(x)=x2 +3x (polinomio de grado 2) entonces la solución


particular es de la forma yp(x)= Ax2+Bx+C (polinomio
generalizado de grado 2). Luego debemos determinar los
coeficientes A,B y C.

La solución particular debe satisfacer a la ecuación no


homogénea; es decir, yp” +4yp’+9yp = x2 + 3x
Hallemos la primera y segunda derivada para yp(x)= Ax2+Bx+C
yp’= 2Ax+b
yp”= 2A

Remplazando y agrupando

2A + 8Ax + 4b + 9Ax2 + bx + c = x2 + 3x
9Ax + (8A + 9b)x + (2A + 4b + 9c) = x2 + 3x + 0
2

Si dos polinomios son iguales sus coeficientes deben ser iguales


9 A= 1
{
Entonces= 8 A+ 9 B=3
2 A +4 B+9 C=0

Resolviendo el sistema simultáneo tenemos:

1 19 94
A¿ ,B¿ y C=
9 81 729
1 19 94
Por, tanto yp(x) = x2 + x-
9 81 729

Finalmente la solución general sería:


1 19 94
y(x)= e-2x[k1sen(√5x) + k2 cos (√5x)] + x2 + x-
9 81 729

19
Ecuaciones Diferenciales Ordinarias

IV. EJERCICIOS DE APLICACIÓN

1. El sistema de la figura consta de una masa, dos muelles y un


Amortiguador de características: m =20kg ; k =50 N /m ; k= 70 N /m ; c=
80 N s /m
Determinar la Ecuación diferencial del movimiento y su solución general

RESOLUCIÓN
a) La ecuación diferencial del movimiento de una vibración libre amortiguada
es:
mx ´´ +cx ´ +k p x=0
La constante equivalente de los muelles en paralelo es k p=k 1 +k 2=120 N /m
Por lo que:
20 x ´´ +80 x ´ + 120 x=0

Simplificando queda:
x ´ ´ + 4 x ´ +6 x=0

Cuya ecuación característica es:


r2 +4r +6=0 donde r= - 2± √ 2t

la solución general es :
x= A e 2 t sen ( √ 2t +φ ´ )

20
Ecuaciones Diferenciales Ordinarias

2. Determine la ecuación del movimiento y la frecuencia natural de


vibración del sistema masa resorte mostrado en la figura

Aplicando la segunda ley de newton ∑ F=ma


Para oscilaciones verticales las fuerzas que actúan son la fuerza del resorte
K(δ est + x) y el peso mg de la masa por tanto la ecuación del movimiento es:
mx´´ =- K(δ est + x) +mg

Donde x´´ es la segunda derivada respecto al tiempo y δ est es la deflexión


estática debido al peso de la masa que actúa sobre el resorte entonces mg= Kδ
est y la ecuación del movimiento se conviene en:

mx´´+ Kx = 0
que es la ecuación del MAS la solución mas general de esta ecuación es :

x= A sen √ k /m t + B cos√ k /m t
o x = C cos (√ k /m + ∅)

donde A, B y C son constantes que dependen de de las condiciones iniciales


x(0)
y x´(0) en la solución general deben aparecer dos constantes puesto que se
trata de una ecuación diferencial de segundo orden para un desplazaminto
inicial tenemos que A=0 Y B = X(0) y por tanto :

X= X(0) cos (√ k /m t)
Físicamente esta ecuación representa una ecuación libre sin amortiguamiento
en la cual se produce un ciclo cuando √ k /m t varía 360° por tanto:

Periodo T = y frecuencia natural f =1/T
√ k /m
Donde √ k /m =Wn rad/s es la frecuencia natural angular.

21
Ecuaciones Diferenciales Ordinarias

3. Un manómetro utilizado en un laboratorio de mecánica de fluidos tiene


una sección transversal de área A .Si la columna de líquido de longitud
L y densidad ρ se pone en movimiento como se muestra en la figura
encuentre la frecuencia de movimiento resultante.

Aplicando la segunda ley de Newton:


∑ F=mx´ ´ o -2AX0 ρ=LAρX´´

Dónde:
X´´+ (2g/L)X=0 y Wn=√ 2 g /L rad/s

Método de la energía:

1 1
EC = mx’2 = (LAρ ) x’2
2 2
1 1
Ep= kx2 = (2Aρg/x) x2 = Ax2 ρg
2 2
d
y (Ec+Ep) =LAρ X’X’’+ 2AρgXX’=0
dt
Dónde:
X’’+ (2g/L) = 0 y Wn=√ 2 g /L rad/s

22
Ecuaciones Diferenciales Ordinarias

4. Del siguiente sistema hallar la ecuación característica donde K=1250N/m y


D=2000N-s/m; M=10Kg y F(t)=1N

d2 x dx d2 x dx
m 2
+D + kx=F 10 2
+2000 + 1250 x=1
dt dt dt dt

2 2
2 d x dx 10 d x 2000 d x 1250 1
T 2
+ 2 ԐT + x=K 2
+ + x=
dt dt 1250 dt 1250 dt 1250 1250

−3 d2 x dx −4
8∗10 2
+1.6 + x=8∗10 T 2=8∗10−3 →T =√ 8∗10−3 →T =0.089
dt dt

1.6
2 ԐT =1.6 → Ԑ = → Ԑ =8.988 Ԑ >1: sobre amortiguad .
2∗0.089
−Ԑ ± √ Ԑ 2−1
P1,2=
T
−8.988 ± √ 8.9882−1
P1,2=
0.089
P1=−0.627˄ P2=−199.4
k P2 P ∗t k P1 P ∗t
X ( t )= e + 1
e +k 2

P 1−P2 P2 −P 1

8∗10−4∗−194.4 −0.627∗t (
X ( t )= e + 8∗10−4 )∗(−0.627) e−199.4∗t + 8∗10−4
−0.627−(−199.4)
X ( t )=−8∗10−4∗e−0.627∗t +2.52∗10−6 e−199.4∗t +8∗10− 4

23
Ecuaciones Diferenciales Ordinarias

5. Una masa de 2lb, un resorte se estira 0.5ft en t=0, luego se suelta la masa y se
' −4
estira el resorte 2/3ft debajo de la posición inicial. Hacia arriba: X (t )= ft .
3

d2 x k
+ x=0
dt 2 m
mg=ks
1
k =2→ k =4
2
w=mg
1
2=m∗32→ m=
16
4
x+ x=0 → x +64 x=0
1
16
m +64=0 → m 2−( 8i )=0
2

( m−8 i ) ( m+ 8i )=0
m 1=0+ 8i ˄m 2=0−8 i
x ( t )=c1 cos 8 t +c 2 seno 8 t
x ( 0 )=c 1 cos 0+c 2 seno 0
2
=c
3 1
x ' ( t)=−8 c1 seno 8 t+8 c 2 cos 8 t
x ' ( 0)=−8 c 1 seno 0+ 8 c 2 cos 0
−4 −1
=8 c 2 → c 2=
3 6
2 1
x ( t )= cos 8 t− seno 8 t
3 6

24
Ecuaciones Diferenciales Ordinarias

6. Resolver e interpretar el problema de valor inicial:


d²x/dt² + 16 x = 0
x(0) = 10, dx/dt = 0
t=0
Solución:
Una formulación equivalente del problema es: se estira hacia abajo de un cuerpo
que pende de un resorte hasta que esté 10 unidades bajo la posición de equilibrio
y luego se le retiene hasta t = 0; se le suelta a continuación de manera que parta
de un estado de reposo. Aplicando las condiciones iniciales a la solución:
x (t) = C1 cos 4t + C2 sen 4t.
Resulta x (0) = 10 = C1 . 1 + C2 . 0
de modo que C1 = 10 y por lo tanto
x (t) = 10 cos 4t + C2 sen 4t.
dx/dt = 40 sen 4t + 4C2 cos 4t
dx/dt = 0 = 4C2 . 1
ecuación del movimiento x (t) = 10 cos 4t + 0
está claro que la solución indica que el sistema permanece en movimiento una
vez puesto en movimiento y la masa va y viene 10 unidades a cada lado de la
posición de equilibrio x=o.

7. Una masa de 8 lb de peso estira 2ft un resorte. Si una fuerza de amortiguamiento


numéricamente igual a 2 veces la velocidad instantánea actúa sobre el
contrapeso, deduzca la ecuación del movimiento si la masa se suelta de la
posición de equilibrio con una velocidad hacia arriba de 3 ft/s.

Solución:
De acuerdo con la ley de Hooke, 8 =K(2) da k=4 lb/ft. Entonces W=mg da
m=(8/32) =1/4 slug. Entonces la ecuación diferencial del movimiento es:
1 d2 x dx d2 x dx
2
=−4 x− o sea 2 +8 + 16 x =0 … …(S)
4 dt dt dt dt
La ecuación auxiliar de (17) es m2+8m+16=(m+4)2=0, de forma que m1=m2=-
4.luego el sistema es críticamente armónico y
x ( t )=c1 e−4 t + c 2 te −4 t
Aplicar las condiciones iniciales X(0)=0 y X’(0)=-3 veos, a su vez que c1=0 y
c2=-3 así la ecuación del movimiento es
x ( t )=−3t e−4 t

25
Ecuaciones Diferenciales Ordinarias

Para graficar x(t). de x’(t)=−3 t e− 4 t (1−4 t) tenemos que x’(t)=0 cuando t=1/4.
1 1 −1
El desplazamiento extremo correspondiente es x () ()
4
=−3
4
e =−0.276.
En la figura vemos que podemos interpretar este valor como el punto en que el
contrapaso alcanza una altura máxima de 0.276ft sobre su posición de equilibrio.

8. Un objeto que pesa 16 lb se une a un resorte de 5 ft de longitud. En la posición


de equilibrio, el resorte mide 8.2 ft. Si el peso se eleva y se suelta del reposo en
un punto a 2ft arriba de la posición de equilibrio, determine los desplazamientos,
x(t). Concidere que el medio que rodea al sistema ofrece una resistencia al
movimiento numéricamente igual a la velocidad instantánea
Solución:

S=8.2−5=3.2 ft

m=Ks → 16=K (3.2) → K =5 lb/ ft

16 1
W = mg → m= = slug
32 2
La ecuación diferencial es:
1 d2 x dx
2
=−5− ó
2 dt dt
2
d x dx
2
+2 +10 x =0
dt dt

m 2 +2 m+10=0 → m 1=−1+3 i → m 2=−1−3 i, lo cual implica el sistemas


subamortiguado y que:
x ( t )=e−t (c 1 cos 3 t+c 2 seno 3 t )
x ( 0 )=c 1 cos 0+c 2 seno 0
c 1=−2
x ' (t)=0
−2
c 2=
3
−4 −1
=8 c 2 → c 2=
3 6
2 1
x ( t )= cos 8 t− seno 8 t
3 6

26
Ecuaciones Diferenciales Ordinarias

Entonces la ecuación de movimiento e


−2
x ( t )=e−t (−2 cos 3 t+ seno 3 t)
3
9. Interpretar y resolver el problema de valor inicial
1 d2 x dx 1 '
2
+1.2 +2 x=5 cos 4 t ; x ( 0 )= , x ( 0 )=0
5 dt dt 2
Solución:
1 d2 x 1
2
+1.2+2 x=5 cos 4 t , x ( 0 )= , x ´ ( 0 ) =0
5 dt 2

1
m= ft → K=2 lb /ft
2
π
Movimiento amortiguado:β=1.2 → T= s( se inicia cuando T =0)
2
f ( t )=5 cos 4 t
Multiplicamos por 5 la ecuación diferencial:
d2 x dx
2
+6 +10 x=0
dt dt
Resolvemos con los métodos acostumbrados
m2 +6 m+ 10=0 → m1=−3+i→ m 2=−3−i
Entonces
x ( t )=e−3 t (c 1 cost+ c2 sent )
Aplicamos x p ( t )= Acos 4 t+ Bsen 4 t ,

x ´ p ( t )=−4 Asen 4 t +4 Bcos 4 t , x ´ ´ p =−16 Acos 4 t −16 Bseno 4 t


De modo que
x ´ ´ p+ 6 x ´ p +10 x p=(−6 A+ 24 B ) cos 4 t + (−24 A−6 B ) seno 4 t=25 cos 4 t
El sistema resultante de ecuaciones
-6A+24B=25, -24A-6B=0
−25 50
A= B=
102 51
−3 t 25 51
x ( t )=e ( c1 cost +c 2 sent )− cos 4 t + sen 4 t
102 50
Cuando t=0 en la ecuación de arriba obtenemos:
38
c 1=
51
Si diferenciamos la expresión y hacemos t=0 obtenemos:
−86
c 2=
51
Por consiguiente, la ecuación de movimiento es:
38 86 25 51
(
x ( t )=e−3 t
51 51 )
cost− sent −
102
cos 4 t+ sen 4 t
50

27
Ecuaciones Diferenciales Ordinarias

Bibliografía
 http://www.elprepo.com/2009/05/libro-pre-calculo-moises-villena.html
 Guía UPN Cálculo

28

También podría gustarte