Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Notas 2 Ecuaciones Diferenciales

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 296

Unidad 1

ECUACIONES DIFERENCIALES
ORDINARIAS DE PRIMERORDEN

MG NOLAN JARA JARA


MODELACIÓN
DINAMICA POBLACIONAL: La razón con la que dP
la población P crece es proporcional a la = kP
población P en el tiempo t. dt
DECAIMIENTO RADIOACTIVO: La razón con la
dA
que los núcleos de una sustancia se = kA;
desintegran es proporcional a la cantidad de dt
núcleos que quedan en el tiempo t.

LEY DE ENFRIAMIENTO/CALENTAMIENTO DE NEWTON: La


razón con la que cambia la temperatura T de un cuerpo es
= k(T −Tm )
proporcional a la diferencia entre la temperatura T del
dT
cuerpo y dt
la temperatura Tm del medio que lo rodea en el
tiempo t.
MEZCLAS: La razón con la que la cantidad de sal dA
cambia es igual a la razón de entrada menos la = R E − Rs
razón de salida en el tiempo t. dt
RE=(concentración de entrada de sal).(razón de entrada del fluido)
RS=(concentración de salida de sal).(razón de salida del fluido)
MG NOLAN JARA JARA 2
dx
= kx(n + 1 − x)
dt

MG NOLAN JARA JARA 3


EDO DE PRIMER ORDEN

Hallarla soluciondel problema de valor inicial :


1)y  = x, y(0) = 1
2) y = − y, y(0) = 1
1
3) y  = , y(0) = 4
y

MG NOLAN JARA JARA


INTRODUCCIÓN
Los términos ecuaciones y diferenciales
nos hacen pensar en la solución de ciertos
tipos de ecuaciones que contienen
derivadas ó diferenciales.

Ecuación con Ecuación con


derivadas diferenciales

dy
x + 2y = 0 xdx + (2xy −1)dy = 0
dx
MG NOLAN JARA JARA
ECUACIÓN DIFERENCIAL
Definición

Una ecuación diferencial (ED) es una ecuación que


relaciona a una función desconocida y una o más
derivadas de esta función con respecto a una o
más variables independientes.
Si la función desconocida depende de una sola
variable, la ecuación diferencial se llama ordinaria.
Si por el contrario dependiese de más de una
variable, se llama parcial .

MG NOLAN JARA JARA


EDO y EDP

EDO EDP

w w
dy
= x 2 y − cosy − = x2 + y
dx x y

y''' = y − y'
3

MG NOLAN JARA JARA


OBSERVACIÓN

Nuestra atención se centrará sobre ecuaciones


diferenciales ordinarias. Una ecuación diferencial
ordinaria es aquella que tiene a y como variable
dependiente y a x como variable independiente.

Se acostumbra expresar en la forma:

( , ( n)
F x, y, y y´´ ,..., y = 0 )
MG NOLAN JARA JARA
ORDEN

El orden de una ecuación diferencial es igual al de


la derivada de más alto orden que aparece (de
manera no trivial) en la ecuación.
Ejemplos:

3xy`` + 5(y`)3 = 6x 5y 4 es de 2do orden.

𝒚´´
xexy``` − + sen(xy)= 0 es de 3er orden.
𝒙

MG NOLAN JARA JARA


EDO LINEAL Una ecuación diferencial
ordinaria de orden n es lineal, si se puede
escribir de la forma:

(𝑛)
an (x)𝑦 + 𝑎𝑛−1(𝑥)𝑦 (𝑛−1) +…+𝑎0 𝑥 𝑦 = g(𝑥)

ak (x); k = 0,...n ; an (x)  0 funciones que dependen de x

Una ecuación diferencial ordinaria que no se


pueda expresar de esta forma es no-lineal.

MG NOLAN JARA JARA


SOLUCIÓN DE UNA EDO
Decimos que y =Ф(x) es una solución de la
ecuación diferencial:

( )
F x, y, y, y ,..., y(n ) = 0

en el intervalo I si: F ( x,, , ,...,  (n)


)= 0
para todo x є I.

MG NOLAN JARA JARA


EJEMPLOS

Indique si las funciones dadas son


soluciones de las EDO en cierto intervalo I.
1 4
función: y= x EDO : y'= xy1/ 2
16

función : y = 2x EDO : y'+3x = 5

1
función : y = EDO : xy'+ y =0
x

MG NOLAN JARA JARA


USO DE SOFTWARE

Utilice un software para comprobar si la siguiente función


es una familia de soluciones de la EDO.

función: y = Ae2 x + Bxe2x

EDO : 𝑦´´ − 4𝑦´ + 4𝑦 = 0

MG NOLAN JARA JARA


USO DEL DERIVE

Al sustituir la familia de funciones en la EDO,


se satisface la igualdad.
Luego, sí es una familia de soluciones de la EDO.
MG NOLAN JARA JARA
INTERVALO DE VALIDEZ

1 xy  + y = 0
y = f (x) = 𝑥
es solución de:

Como función Como solución:


I = 0,+
domf = R − 0
I = − ,0
MG NOLAN JARA JARA
INTERVALO DE VALIDEZ

1 
𝑦= 2 es solución general de: y + 2xy 2
=0
𝑥 +𝑐
Si y(0)=-1
1
𝑦 = 𝑓(𝑥) = 2
𝑥 −1

Como función: Como solución:


Dom(f )= 𝑥 ∈R −−1,1 I = 𝑥 ∈ −1,1
MG NOLAN JARA JARA
SOLUCIÓN IMPLÍCITA

La relación G(x,y)=0 se llama solución


implícita de una EDO en un intervalo I si
existe alguna función Ф que satisface tanto
la relación como la EDO en I.

MG NOLAN JARA JARA


EJEMPLOS

Indique si las relaciones dadas son


soluciones implícitas de las EDO en I.

x 2 + y 2 = 25 x
relación: EDO : y' =−
y

relación: 𝑥𝑦 2 − 𝑒 −𝑦 − 1 = 0 EDO : (2xy + e )y'+y


−y 2
=0

MG NOLAN JARA JARA


FAMILIA DE SOLUCIONES

Algunas veces, a una solución de una ecuación


diferencial se le llama integral de la ecuación y a su
gráfica curva integral o curva solución. Como la
solución general de una ecuación diferencial lineal
de orden n tiene n constantes, se acostumbra
llamarla familia n-paramétrica de soluciones y se
denota por:
G(x, y, c1 , c2,..., cn ) = 0
Esto quiere decir que una EDO tiene una cantidad
infinita de soluciones que corresponden a la
elección arbitraria de esos parámetros.
MG NOLAN JARA JARA
SOLUCIÓN GENERAL Y PARTICULAR

Si encontramos una familia n-paramétrica de


soluciones que contiene a TODAS las soluciones
de una EDO, llamaremos a esta familia solución
general de la EDO.
Dada una familia n-paramétrica de soluciones de
una EDO, una solución que se obtiene al dar
valores a los n parámetros se llama solución
particular.

Observación: Las EDO lineales siempre tienen


solución general.

MG NOLAN JARA JARA


EJEMPLOS
EDO : y = xy'+(y') 2

1)Verifique si la familia mostrada


es una familia uniparamétrica de
y = cx + c 2

soluciones de la EDO:
2) Verifique si las siguientes a) y = 3x + 9
funciones son soluciones de la b) y=0
EDO y clasifíquelas como
solución particular y/o trivial. c) y = −x 2 / 4
3) Diga si la familia uniparamétrica de 1) es o
no solución general de la EDO. ¿Por qué?

MG NOLAN JARA JARA


PROBLEMA DE VALOR INICIAL (PVI)

Un problema de valor inicial (o de Cauchy) consta


de una ecuación diferencial de orden n y de n
condiciones iniciales impuestas a la función
desconocida y a sus n-1 primeras derivadas en un
valor de la variable independiente.
Es decir:
dn y
dx n
= (
f x, y, y ,...,y (n−1)
)
y (x0 ) = y 0 ,y  (x0 ) = y1,..., y (n−1)
(x0 ) = y n−1

MG NOLAN JARA JARA


EJEMPLO

y´= 2x + 3y
y(1) = 2

MG NOLAN JARA JARA


MODELACIÓN (Paracaídas)

Por la 2da Ley de Newton:


2

 F = m dt2 = m dt
d y dv

dv
m = mg − kv
dt
mg = Peso
kv = Re sistencia Aire
kv mg F (resultan te) = ma

MG NOLAN JARA JARA


MODELACIÓN (varios)

Resolver los ejemplos del libro de Dennis G.


Zill páginas 20 a 26.

MG NOLAN JARA JARA


EDO DE PRIMER ORDEN
CLASIFICACIÓN

EDO VARIABLE SEPARABLE


EDO EXACTA
EDO HOMOGENEA
EDO LINEAL
EDO BERNOULLI
EDO RICATTI
EDO LAGRANGE
EDO CLAIROUT

MG NOLAN JARA JARA


EDO VARIABLE SEPARABLE

dy dy
=  F (x)dx
dy
1) = F (x)G( y)  = F (x)dx  
dx G( y) G( y)
2) P1(x)P2 ( y)dx + Q1 (x)Q2 ( y)dy = 0 
P1(x)P2 ( y)dx = −Q1 (x)Q2 ( y)dy 
P1(x) Q ( y)
dx = − 2 dy 
Q1 (x) P2 (y)
P1(x) Q ( y)
dx + 2 dy = 0 
Q1 (x) P2 (y)
P1(x) Q2 ( y)
 Q (x)  P (y) dy = c
dx +
1 2

MG NOLAN JARA JARA


27
EDO VARIABLE SEPARABLE
Resolver ejemplo 1 ó 2 de la pág. 46 de Zill.
1. (1+x)dy-ydx=0
2. yy`=-x,y(4)=-3
3. y`=(y+2)(y-2)
4. xdx+2y²dy=0
𝒅𝒚
5. 𝒆 𝟐𝒚
− 𝒚 𝒄𝒐𝒔𝒙 = 𝒆𝒚 𝒔𝒆𝒏 𝟐𝒙 ; 𝒚 𝟎 = 𝟎
𝒅𝒙
6. (xy²-y²+x-1)dx+(x²y-2xy+x²+2y-2x+2)dy=0

MG NOLAN JARA JARA


(xy²-y²+x-1)dx+(x²y-2xy+x²+2y-2x+2)dy=0

𝒅𝒚
7. 𝒚 𝒚
𝒆 𝒄𝒐𝒔𝒙 + 𝒄𝒐𝒔𝒙 + 𝒆 𝒔𝒆𝒏𝒙 + 𝒆 𝒚
= 𝟎; 𝒚 𝝅 = 𝟎
𝒅𝒙
𝟏−𝒔𝒆𝒏𝒙
Solución. 𝒚 = 𝒍𝒏 𝟏+𝒔𝒆𝒏𝒙

MG NOLAN JARA JARA


8. La rapidez con que se desintegran núcleos radiactivos es
proporcional al numero de núcleos que están presentes en
una muestra dada. La mitad del numero original de núcleos
radiactivos ha experimentado la desintegración en un periodo
de 1500 años. ¿ Que porcentaje de núcleos radiactivos
originales continuaran después de 4500 años?
¿ En cuantos años quedara solamente un decimo del numero
original de núcleos radiactivos?
Solución.
N(t)…cantidad de núcleos radiactivos después de t años
N´(t)… rapidez con que se desintegran núcleos radiactivos
N(0)=𝑵𝟎 … 𝒏𝒖𝒎𝒆𝒓𝒐 𝒐𝒓𝒊𝒈𝒊𝒏𝒂𝒍 𝒅𝒆 𝒏𝒖𝒄𝒍𝒆𝒐𝒔 𝒓𝒂𝒅𝒊𝒂𝒄𝒕𝒊𝒗𝒐𝒔
𝑵𝟎
N(1500)= 𝟐
𝑵
¿% N(4500)?; ¿t/N(t)= 𝟏𝟎𝟎?

MG NOLAN JARA JARA


𝑵𝟎
N´(t)=k N(t); N(0)=𝑵𝟎 ; N(1500)= 𝟐 ; ¿ 𝑵 𝒕 ?
N(t)=c𝒆𝒌𝒕 ; 𝑵𝟎 =N(0)=c 𝒆𝒌.𝟎 c→ 𝒄 = 𝑵𝟎 → N(t)=𝑵𝟎 𝒆𝒌𝒕

𝑵𝟎 𝟏 𝟏 𝟏
= N(1500)= 𝑵𝟎 𝒆𝒌.𝟏𝟓𝟎𝟎 →𝒆 𝒌.𝟏𝟓𝟎𝟎 = →k= ln
𝟐 𝟐 𝟏𝟓𝟎𝟎 𝟐
−𝒍𝒏𝟐
−𝒍𝒏𝟐
k=𝟏𝟓𝟎𝟎 → N(t)=𝑵𝟎 𝒆𝟏𝟓𝟎𝟎𝒕

−𝒍𝒏𝟐
(𝟒𝟓𝟎𝟎) 𝑵 𝑵𝟎
N(4500)= 𝑵𝟎 𝒆 𝟏𝟓𝟎𝟎 = 𝑵𝟎 𝒆−𝟑𝒍𝒏𝟐 = 𝒆𝒍𝒏𝟖
𝟎
= =12.5% 𝑵𝟎
𝟖

MG NOLAN JARA JARA


La EDO de la forma
y’= f ( ax + by + c)
se reduce a EDO de variable separable
haciendo
z = ax + by + c
Resolver:

1)y’= ax+by+c
2) y  − e x e y +1 = 0

MG NOLAN JARA JARA 31


2) y  − e x e y +1 = 0

MG NOLAN JARA JARA


FORMA DIFERENCIAL EXACTA
La forma diferencial
P (x, y)dx + Q(x, y)dy (1)
𝝏𝑷 𝝏𝑸
es exacta en un rectángulo R, si =
𝝏𝒚 𝝏𝒙

y existe una función f (x,y) que cumple:

𝝏𝒇 𝝏𝒇
=𝑷; =𝑸
𝝏𝒙 𝝏𝒚

para todo (x,y) en R.

a f (x,y) se le llama Función Potencial.

MG NOLAN JARA JARA


EDO EXACTA
Definición: Una EDO de la forma
P (x, y)dx + Q (x, y)dy = 0
se llama exacta si la forma diferencial (1) es
exacta.
Ejemplo:

(x − 2 y)dx − (2x + y)dy = 0


es exacta.

MG NOLAN JARA JARA


El diferencial de f(x,y) satisface:

df = P (x, y)dx + Q(x, y)dy


Ejemplo: Verificar en la forma diferencial
(x − 2y)dx − (2x + y)dy
2 2
x y
f (x, y) = − 2xy −
2 2
¿Es f única o no?

MG NOLAN JARA JARA


CONDICIÓN

En (1) Sean P, Q y sus derivadas parciales


de primer orden continuas en R, luego la
condición necesaria y suficiente para que la
forma diferencial sea exacta es:

P Q
=
y x

MG NOLAN JARA JARA


ECUACION DIFERENCIAL EXACTA
Sea la EDO P(x,y)dx + Q(x,y)dy = 0. Y si

P(x, y) Q(x, y)
= 
y x
La EDO P(x,y)dx + Q(x,y)dy = 0 es una Ecuación
diferencial exacta. Con solución general: f(x,y)=c

f (x, y) f (x, y)
P(x, y) = ; Q(x, y) =
x y
x y
f ( x, y) =  P (u, y)du +  Q (a, v)dv = c
a a

MG NOLAN JARA JARA


DEMOSTRACION. Si la EDO P(x, y)dx + Q(x, y)dy = 0 es exacta
f (x, y) f (x, y)
P(x, y) = ; Q(x, y) =
x y
x
f (x, y)
de P(x, y) =  f (x, y) =  P(u, y)du + h( y)… (*)
x a

f (x, y)  x   x

   P(u, y)du + h( y)  =
y
=
y  a  y
 P(u, y)du
a
+ h( y)… (*1)

f (x, y)  x
de Q(x, y) =
y
y (*1) Q(x, y) = 
y a
P(u, y)du + h( y)
x

 Q(a, y) = h( y)  h( y) =  Q (a, v)dv en (*)


a
x x

f (x, y) =  P(u, y)du +  Q(a, v)dv


a a

MG NOLAN JARA JARA


RESOLVER : 1) 3y + e + (3x + cos y) y = 0
x

2) ydx + (y 2
+ x)dy = 0
y2x − y
3)y =
x-x2 y
2xy + x2-x
4 )y  =
y-x2
y
5) dx + ( y3 + ln x)dy = 0
x
6)2xydx+ (x 2 −1) dy = 0
7)(e 2 y − y cos(xy))dx + (2 xe2 y − x cos(xy) + 2 y)dy = 0
xy2 − senxcos x
8) y = 2 , y(0) = 2
y(1− x )
MG NOLAN JARA JARA
REDUCIBLES A ECUACION DIFERENCIAL EXACTA
Se puede multiplicar una ecuación diferencial no
exacta por un factor u tal que se convierta en
exacta. Este factor se llama FACTOR INTEGRANTE.
uP(x, y)dx + uQ(x, y)dy = 0 es exacta 
(uP(x, y)) (uQ(x, y))
= 
y x
u P(x, y) u Q(x, y)
P(x, y) +u = Q(x, y) +u
y y x x
 P(x, y) Q(x, y)  u u
u −  = Q(x, y) − P(x, y)
 y x  x y

MG NOLAN JARA JARA


Si el factor u fuese sencillo, por
ejemplo u (x) o u (y)
podríamos simplificar la EDO
por ejemplo, si consideramos que
u = u (x):

MG NOLAN JARA JARA


u(x) Py − Q x
= = g(x)
u Q
 (l n u ) = g ( x )
'

 u(x) = e  g ( x )dx

 Py − Q x 
  dx
u(x) = e  
Q

MG NOLAN JARA JARA


Analogamente si consideramos u = u (y) :
u( y) Qx − Py
= = g( y)
u P
 (ln u) = g( y)
'

 u( y) = e  g ( y )dy

 Qx −Py 
 
 P
 dy

u( y)= e
Ejemplo:1)ydx + (x 2 y - x)dy = 0
 2x + y -1
2)y = y
x +y
MG NOLAN JARA JARA
3)ydx + ( x 2 y - x)dy = 0
2x2 + y
4) y  =
x - x2 y
2xy + y 2 - y y
5 )y  = ;4) d x + ( y 3 + ln x ) d x = 0
x +y x

MG NOLAN JARA JARA


FUNCIÓN HOMOGÉNEA
Definición: Una función f (x,y) se llama
homogénea de grado n si:

f (  x,  y) =  f ( x, y)
n

Ejemplo:

f ( x, y) = x − y + xy
2 2

MG NOLAN JARA JARA


EDO HOMOGENEAS
una función f(x, y) es HOMOGENEA de grado n si

f(  x,  y) = 
n f(x, y)

La EDO y’ = f(x,y) es HOMOGENEA si f es


homogénea de grado 0
La EDO P(x, y)dx + Q(x, y)dy = 0 Es HOMOGENEA
si P y Q son funciones homogéneas del mismo grado.
PASO DE HOMOGENEA A VARIABLES SEPARADAS
Las EDO HOMOGENEAS se integran con el cambio
y
z =  y( x) = xz( x)
x
y=zx→dy=zdx+xdz
MG NOLAN JARA JARA
RESOLVER :
3
2y
1)(y − x )dx −
2 2
dy = 0
x
2)xdy − ydx = x 2
+ y dx
2

3)4x − 3y + (2y − 3x)y = 0


4)xy = y + y −x 2 2

MG NOLAN JARA JARA 47


RESOLVER :
4
y
5)( x − y ) d x − 2 d y = 0
2 2
x
6) xdy + ydx = x 2
+ y dx2

7 ) 4 x + 3 y + (2 x − 3 y ) y  = 0
8) xy = y - y2 − x2
9 ) (y − x y ) = x 2 + y 2
2

MG NOLAN JARA JARA


y d x = 2 (x + y )d y
EJEMPLO:

dy y Hacemos: y = ux
=
dx 2(x + y) 
dy
=u+x
du
dx dx
du ux
u+x =
dx 2(x+ ux)
du u EDO de Variable
u+x = Separable
dx 2(1+ u)
Resolviendo: cy − 2 y − x = 0
2

MG NOLAN JARA JARA


EDO REDUCIBLES A HOMOGENEAS
En la EDO:
(A1x + B1y + C1)dx + (A2x + B2y + C2)dy = 0
Hacer el cambio de variables:
x = u + r1;y = v + r2
(r1,r2):intersección de las rectas no
paralelas
A1x + B1y + C1=0; A2x + B2y + C2=0
En la EDO:
(Ax + By + C1)dx + (Ax + By + C2)dy = 0
Hacer el cambio de variables: ax+by = u
MG NOLAN JARA JARA
RESOLVER :
1)(2x − 5y + 3)dx − (5x − 12y + 8)dy = 0
2) ( x − 2 y + 1)dx + ( 2x − 4 y + 3)dy = 0
3x − y + 2
3) y  =
6x − 2 y

MG NOLAN JARA JARA


ECUACION DIFERNCIAL LINEAL
ORDINARIA
La forma general es

n−1
n
d y d y dy
a n (x) n + an-1 (x) n−1 + … + a1 (x) + a0 (x) y = f (x)
dx dx dx

MG NOLAN JARA JARA


ECUACIÓN DIFERENCIAL LINEAL DE PRIMER
ORDEN
dy
Una EDO de la forma: a1 (x) + a0 (x) y = g(x)
dx
se llama ecuación diferencial ordinaria lineal de
primer orden (EDOL).
Realizando las operaciones adecuadas se
escribe en la forma estándar:

y'+ p(x) y = q(x)


MG NOLAN JARA JARA
FACTOR INTEGRANTE DE UNA EDOL

Un factor integrante
para una EDOL de
primer orden es:

u(x) = e p( x)dx

Resolver un ejemplo de la pág 55 del Zill.

MG NOLAN JARA JARA


DEMOSTRACION DE LA EDOL DE
PRIMER ORDEN DE LA FORMA
y  + p(x) y = q(x)...(*)
Para resolver esta ecuación, escribirlo como:
( p(x) y − q(x) ) dx + dy = 0
tal que : P(x, y) = p(x) y − q(x); Q(x, y) = 1
El Factor Integrante sera:
 P −  Q 
 
y x
 
 Q
dx

u(x) = e



 =e  p ( x)dx

MG NOLAN JARA JARA


en (*) Multiplicando por el factor Integrante
e  p( x)dx
y + e  p( x)dx
p(x) y = e  p( x)dx
q(x)

 e  p( x)dx y  = e  p( x)dx q(x)
 
 
e  p( x)dx
y = e  p( x)dx
q(x)dx + c

 − p( x)dx  e  p( x)dx q(x)dx+ c 


y=e  
 

MG NOLAN JARA JARA


RESOLVER LAS EDOL
1) y  − 3 y = 6
2)xy − 4 y = x 6ex
( )
3) x 2 − 9 y  + x y = 0
4) y + y = x, y(0) = 4
5) y  + 2 x y = 2 x e − x
2

(
6 ) 2 x d y = 2x 3− y d x
2y
)

7) y + = x3
x

MG NOLAN JARA JARA


Pr oblema : Re solver
1) y - 2 x y = 2 x e − x2

2 ) yx  ( y ) + (1 + y ) x = e y

3) 2 x d y = (2 x +y )d x
3

2 y
4) y - = x 3
x

MG NOLAN JARA JARA


MG NOLAN JARA JARA
y´+2y=4x;y(0)=3…EDOL; u(x)=𝑒 ‫ ׬‬2𝑑𝑥 = 𝑒 2𝑥 ; 𝑝𝑜𝑟 𝐸𝐷𝑂𝐿
𝑒 𝑒 2𝑥 2𝑥
𝑦. 𝑒 2𝑥 ´= 4𝑥𝑒 2𝑥 → 𝑦. 𝑒 2𝑥 = 2𝑥
‫ ׬‬4𝑥. 𝑒 dx=4(𝑥 2 − ‫ ׬‬2 𝑑𝑥) →
2𝑥 2𝑥𝑒 2𝑥 −2𝑥
𝑦. 𝑒 =2 𝑥𝑒 − 2 +c→ 𝑦 = 2𝑥 − 1 + 𝑐 𝑒 ;3=y(0)=-1+c
𝒚= 𝟐𝒙 − 𝟏 + 𝟒 𝒆−𝟐𝒙

2 2
‫ ׬‬− 𝑥𝑑𝑥 −2𝑙𝑛𝑥 1 1 ´ 4
y´− 𝑥y=4x; u(x)=𝑒 =𝑒 =𝑥 2 → 𝑦. 𝑥 2 = 𝑥
1
𝑦. 𝑥 2 =4lnx+c→ 𝑦 = 4𝑥²lnx+cx²

MG NOLAN JARA JARA


P1) Considere que un tanque contiene 1000 litros de agua pura y esta
conectada a líneas de abasto y de descarga (ver figura).en t=0, tanto la
línea de abasto como la de descarga están abiertas y la salmuera que
contiene 0.1 kg de sal por litro de agua entra y sale del tanque a razón
de 50 lt/min después de mezclarse perfectamente con el agua del
tanque. Suponga que la sal disuelta no cambia el volumen del agua, el
contenido de sal en el tanque aumenta con el tiempo, aun cuando el
volumen de agua permanezca constante. Obtenga la relación para la
cantidad de sal en el tanque en cualquier tiempo t, y determine la
cantidad máxima de sal que el tanque contendrá finalmente.
50 lt/min
0.1 kg sal/lt

50 lt/min
Agua

MG NOLAN JARA JARA


Sea x la cantidad de sal en el tanque en un determinado tiempo t. la
variación de la cantidad de sal en el tanque puede expresarse como:
𝒅𝒙 𝒙
= 𝟎. 𝟏 𝒌𝒈/𝒍𝒕 𝟓𝟎 𝒍𝒕/𝒎𝒊𝒏 − 𝟏𝟎𝟎𝟎 𝒌𝒈/𝒍𝒕 𝟓𝟎 𝒍𝒕/𝒎𝒊𝒏
𝒅𝒕
𝒅𝒙 𝒙 𝒅𝒙 𝟏
𝒅𝒕
+ 50 𝟏𝟎𝟎𝟎
=50(0.1) → + 𝒙=5→ 𝒙 𝒕 = 𝟏𝟎𝟎
𝒅𝒕 𝟐𝟎
𝟏 − 𝒆−𝟎.𝟎𝟓𝒕

MG NOLAN JARA JARA


P2) Suponga que en el circuito simple
la resistencia R es de 12 ohms Ω y que
la inductancia L es de 4 henrys H.
Si una batería proporciona un voltaje
Constante E de 60 volts V y el
Interruptor está cerrado cuando t = 0,
de modo que la corriente comience
con I(0) = 0, determine: I(t), la corriente
después de 1 segundo, El valor límite de
la corriente
dI
Solución Si L = 4, R = 12 y E(t) = 60, en la EDO Lineal L dt + RI = E
obtenemos el problema con valor inicial
dI
4 dt + 12I = 60, I(0) =0→ I(t) = 5 1 − 𝑒 −3𝑡 ; I(1) = 5 1 − 𝑒 −3 ≈ 4,7
Amperes A
lim 𝐼 𝑡 = lim 5 1 − 𝑒 −3𝑡 =5
𝑡→∞ 𝑡→∞

MG NOLAN JARA JARA


APLICACIONES DE ECUACIONES DIFERENCIALES ORDINARIAS
DE PRIMER ORDEN A PROBLEMAS DE CIRCUITOS

MG NOLAN JARA JARA


MG NOLAN JARA JARA
EJERCICIOS DE APLICACIONES DE ECUACIONES DIFERENCIALES
ORDINARIAS DE PRIMER ORDEN A PROBLEMAS DE CIRCUITOS
1) Un generador con una fem de 100 voltios se conecta en serie con
una resistencia de 10 ohmios y un inductor de 2 henrios. Si el
interruptor K se cierra en tiempo t = 0, establezca una ecuación
diferencial para la corriente y determine la corriente en tiempo t.

E=

R=

L=

Establezca y resuelva una ecuación diferencial para este circuito


eléctrico si el generador de 100 voltios se remplaza por otro con una
fem de 20cos5t volt. di
L + Ri(t) = E(t) ...EDO LINEAL
dt
MG NOLAN JARA JARA
E=
1
Rq´(t) + q(t) = E(t) c=
R=
c

MG NOLAN JARA JARA


1) Un inductor de 0,5 henrios es conectado en serie con una resistencia
de 6 ohmios, un condensador de 0,02 faradios, un generador con un
voltaje alterno dado por 24 sen 10t, t 0, y un interruptor K (Figura).

Establezca una ecuación diferencial para la carga instantánea en el


condensador. Encuentre la carga y la corriente al tiempo t si la carga
en el condensador es cero cuando el interruptor K se cierra en t = 0.

MG NOLAN JARA JARA


Formulación matemática. La caída de voltaje a través de la resistencia es
6I. La caída de voltaje a través del inductor es 0.5 I´(t). La caída de
𝑄
voltaje a través del condensador es 0.02 = 50Q. por la ley de Kirchhoff:
dI
6 I + 0.5 + 50Q = 24sen10t ; I=Q´(t)
dt
2
d Q dQ
0.5 2 + 6 + 50Q = 24sen10t
dt dt
2
d Q dQ
2
+ 12 + 100Q = 48sen10t
dt dt

Las condiciones son Q = 0 e I = Q´(t) = 0 en t = 0.

MG NOLAN JARA JARA


La solución complementaria es.
Q(t)=𝒆−𝟔𝒕 𝒄𝟏 𝐜𝐨𝐬 𝟖𝒕 + 𝒄𝟐 𝒔𝒆𝒏(𝟖𝒕)
Asumiendo la solución particular
2
Q(t)= 𝐴cos 10𝑡 + Bsen(10t); A= - ,B= 0
5
la solución general es:

𝟐
Q(t)=𝒆−𝟔𝒕 𝒄𝟏 𝐜𝐨𝐬 𝟖𝒕 + 𝒄𝟐 𝒔𝒆𝒏 𝟖𝒕 − 𝒄𝒐𝒔(𝟏𝟎𝒕)
𝟓

𝟒 𝟑
De las condiciones iniciales encontramos. 𝒄𝟏 = ; 𝒄𝟐 = 𝟏𝟎
𝟏𝟎
Donde, la solución requerida es:
𝒆−𝟔𝒕 𝟐
Q(t)= 𝟒 𝐜𝐨𝐬 𝟖𝒕 + 𝟑𝒔𝒆𝒏 𝟖𝒕 − 𝒄𝒐𝒔(𝟏𝟎𝒕)
𝟏𝟎 𝟓

MG NOLAN JARA JARA


La corriente I está definida por:

𝑰 𝒕 = −𝟓𝒆−𝟔𝒕 𝒔𝒆𝒏 𝟖𝒕 + 𝟒𝒔𝒆𝒏(𝟏𝟎𝒕)


Gráfica de la Carga

MG NOLAN JARA JARA


PROBLEMAS DE MEZCLAS
Se tiene un tanque con volumen inicial V0 de solución (mezcla de soluto y
solvente). Hay un flujo de entrada y de salida, se quiere calcular la cantidad de
soluto x(t) que hay en el tanque en cualquier instante de tiempo t, en función
de la cantidad inicial de soluto x0 al tiempo de iniciar el proceso de mezclado.

Q2
MG NOLAN JARA JARA
Si la solución que se inyecta al tanque tiene una concentración de C1
gramos de soluto por litro, y fluye al mismo con una tasa de Q1 litros por segundo,
en tanto que la sustancia contenida en el tanque se mantiene bien mezclada por
agitación y fluye hacia fuera de este a una tasa de Q2 litros por segundo.
Sea x(t) la cantidad de soluto en el tanque en un instante de tiempo t. La cantidad de
soluto que fluye hacia el tanque durante ∆t segundos es (Q1C1∆t) gramos. La
cantidad de soluto que fluye hacia fuera del tanque durante el mismo intervalo de
tiempo, depende de la concentración de soluto C2(t) en el tanque al instante t.

𝑪𝟐 (t)

El volumen de líquido en el tanque, en cualquier instante de tiempo


t, viene dado por la ecuación V(t) = V0 + (Q1 – Q2 ) t

MG NOLAN JARA JARA


la variación de la cantidad de soluto en un instante t, es igual a la
diferencia entre la cantidad de líquido que fluye hacia el tanque (Q1C1 ∆t) y la
cantidad de líquido que fluye fuera del tanque (Q2C2∆t):
∆x = ( gramos que ingresan) - (gramos que salen)
∆x = (Q1C1∆t) - (Q2C2∆t) = (Q1C1 -Q2C2) ∆t
ya que ∆t ≠ 0, dividiendo entre ∆t→0; x´(t)= Q1C1 - Q2C2 …3
donde Q1, C1 y Q2 son constantes

MG NOLAN JARA JARA


P1

MG NOLAN JARA JARA


Evaluando la ecuación (2) para t = 30 min
V (30) = V0 + ( Q1 – Q2) 30 sustituyendo V(30), V0 y Q1

120 = 90 + ( 4 - Q2 ) 30 resolviendo = − 30 120 90 4 – Q2


despejando Q2 = 3 gal/min

b) Para determinar la cantidad de sal cuando el tanque se


llena se debe determinar la ley de variación de la cantidad
de sal en el tanque, en cualquier instante t. Para ello, se
sustituyen todos los datos en la ecuación (1)

MG NOLAN JARA JARA


𝑑𝑥 3 𝑑𝑥 3
+ 𝑥 = 4(2) → + 𝑥 =8
𝑑𝑡 90+ 4−3 𝑡 𝑑𝑡 90+𝑡

MG NOLAN JARA JARA


EJERCICIOS SOBRE APLICACIONES DE LAS ECUACIONES DIFERENCIALES
ORDINARIAS DE PRIMER ORDEN A PROBLEMAS DE MEZCLAS
1.Un tanque contiene originalmente 100 galones de agua fresca. Se vierte dentro
del tanque, agua que contiene ½ libra de sal por galón a una rapidez de 2 gal/min y
se permite que salga la mezcla con la misma rapidez. Después de 10 min se para el
proceso y se vierte agua fresca dentro del tanque a la rapidez de 2 gal/min, dejando
salir la mezcla a la misma rapidez. Encontrar la cantidad de sal en el tanque al final
de los 20 min.
2.Un tanque con capacidad de 500 galones contiene inicialmente 200 galones de
agua con 100 lb de sal en solución. Se inyecta al tanque agua cuya concentración de
sal es de 1 lb/gal, a razón de 3 gal/min. La mezcla debidamente agitada y
homogeneizada sale del tanque a razón de 2 gal/min.
a) Encuentre la cantidad de sal y la concentración de sal en el tanque para
cualquier tiempo
b) Determine la concentración de sal en el instante justo en que la solución alcanza
el volumen total del tanque

MG NOLAN JARA JARA


3.Un tanque contiene 450 lt de líquido en el que se disuelven 30 gr de sal: Una
salmuera que contiene 3 gr/lt se bombea al tanque con una intensidad de 6 lt/min,
la solución adecuadamente mezclada se bombea hacia fuera con una intensidad de 8
lt/min. Encuentre el número de gramos de sal y la concentración de sal, que hay en
el tanque en un instante cualquiera
4.Un gran depósito está lleno de 500 gal de agua pura. Una salmuera que contiene 2
lb/gal se bombea al tanque a razón de 5 gl/min. La salmuera, adecuadamente
mezclada, se bombea hacia fuera con la misma rapidez.
a) Halle el número de libras de sal y la concentración de sal en el tanque en un
instante t cualquiera
b) Determine la cantidad de sal y la concentración al cabo de hora y media de
iniciado el proceso de mezclado
c)¿Cuánto tiempo debe transcurrir para que la cantidad de sal en el tanque sea de
632,12 lb?

MG NOLAN JARA JARA


ECUACION DE BERNOULLI
La ecuación diferencial de Bernoulli, es una generalización
de la ecuación diferencial lineal de primer orden :
y + p(x) y = q(x) y n
La sustitución z = y1−n , la ecuaciónse transforma en una ecuación
lineal con z como variable dependiente, resolviéndose de manera
análoga.
PROBLEMA : 1)y - 3xy = y 2 x
2) y - e x y = e 3x y 2
3)y - e x y = e2 x y 3
4
4)xy + 6y = 3xy 3
MG NOLAN JARA JARA
ECUACION DIFER ENCIAL DE RICATTI
Una ecuacióndiferencial tiene la forma de la introducida
por Jacobo Francesco Riccati cuando presenta la estructura
y = P( x) y 2 + Q( x) y + R( x)
1
Para resolverla, se debe hacer la sustitución , y = y p +
z
donde y p es cualquier solución particular.
1)y = y 2 − 2xy + x 2 + 1, y = px
y 4 2
2) y = y − − 2 , y p =
2
x x x

MG NOLAN JARA JARA


ECUACION DIFE RENCIAL DE LAGRANGE
Una ecuación diferencial de Lagrange presenta la forma :
y = x f ( y ) + g ( y )
Resolviéndose conla sustitución y' = r  dy = rdx
obteniéndose una EDO LINEAL
x (r) + p(r)x = q(r)
1)2y= xy + ylny
2)y = 2xy + lny
3)y = x(1 + y) + y 2
4) y = 2xy  + seny 
2 − 1
5) y = xy 
y
MG NOLAN JARA JARA
ECUACION DIFERENCIAL DE CLAIRAUT
Una ecuacióndiferencial de Clairaut, llamada así en honor a
Alexis - Claude Clairaut, tiene la forma :
y = x y + g ( y )
Comose puede apreciar, esta ecuaciónes una forma particular
de la ecuacióndiferencial de Lagrange, con f (y' ) = y' , por lo
cual, su resoluciónes análoga a la anterior.
1)y = xy + y 2
2)xy  2 − yy − y  + 1 = 0
ay 
3) y = xy +
1 + y 2

MG NOLAN JARA JARA


USO DE DERIVE

3
dy tan x
Resuelva la EDO: =
dx y

MG NOLAN JARA JARA


USO DE CLASSPAD
La EDO es de variable separable. Se
separan las variables y luego se integra :

y dy = tan 3 x d x

MG NOLAN JARA JARA


MODELOS LINEALES

Plantear los ejemplos del 1 al 5 de las


páginas 83 a 87 del texto.

Comentar los resultados (soluciones)


que da el texto.

MG NOLAN JARA JARA


di R E(t)
L + Ri = E  i (t) + i(t) =
dt L L
Si suponemos que el circuito fue exitado en el instante t = 0. El flujo de
corriente se obtiene Como solución del problema de valor inicial
R E(t)
i (t) + i(t) = ;i (0) = 0… (*)
L L
E  − t 
R
1)Si E(t) = E 0 ⎯⎯
→i (t) = 0  1 − e L 
R  
𝐸0
…termino independiente del tiempo de estado estacionario
𝑅

𝑅
𝐸0 −𝐿𝑡
− 𝑒 … termino transitorio cuyo efecto disminuye con el tiempo
𝑅
t L − RLt  E 0 L − RL t
2)Si E(t) = E 0 t ⎯⎯→i (t) = 0 

− 2e − 2 e
E R R  R
MG NOLAN JARA JARA
−at  𝑅
3)Si E(t) = E 0 e →i (t)
𝐸0
= 𝑅 −𝑎𝐿 𝑒 −𝑎𝑡 −
𝐸0
𝑒 − 𝑡
𝐿
𝑅 −𝑎𝐿

4)Si E(t) = E 0 sen t  𝐸0 𝐿 𝑅


𝐸0
→i (t) = 𝑒 −𝐿𝑡
+ 𝑠𝑒 𝑛 ( t-𝛼)
𝑧² 𝑧

MG NOLAN JARA JARA


MODELADO CON ECUACIONES DIFERENCIALES DE
ORDEN SUPERIOR
considerar sistemas dinámicos lineales en los que cada modelo
matemático es una Ecuación Diferencial de Segundo Orden con
Coeficientes Constantes junto con condiciones iniciales especificadas
en un tiempo que tomaremos como t=0

ax´´(t) + bx(t) + cx(t) =f(t); x( 0 ) = x0 , x( 0 ) = x1


f: Es la entrada, función de conducción o función forzada del sistema.
Una solución x(t) de la EDO en un intervalo I que contiene a t=0 que
satisface las condiciones iniciales se llama Salida o Respuesta del
Sistema.

23/06/2021 Mg NOLAN JARA 1


JARA
ECUACIONES DIFERENC IALES LIN EALES DE ORDEN n
𝑎0 (𝑥)𝑦 (𝑛) + 𝑎1 (𝑥)𝑦 (𝑛−1) +𝑎2 (𝑥)𝑦 𝑛−2 +
 ⋯ + 𝑎𝑛−1 (𝑥)𝑦´ + 𝑎𝑛 (𝑥)𝑦 = 𝑔 𝑥
a0(x),a1(x),a2(x), … ,an−1(x),an(x),y,g(x) son funciones que
dependen de x.

𝑑 𝑛𝑦 𝑑 𝑛−1 𝑦 𝑑 𝑛−2 𝑦 𝑑𝑦
(𝑛) (𝑛−1) 𝑛−2
𝑦 = 𝑛 ,𝑦 = 𝑛−1 , 𝑦 = 𝑛−2 , 𝑦´ = ,𝑦 𝑥
𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑑𝑥
Si n = 2  a0(x)y´´ + a1(x)y + a2(x)y = g(x)

Si a0 ,a1 ,a 2 , … ,a n−1 ,an son constan tes, la ecuación

𝑎0 𝑦 (𝑛) + 𝑎1 𝑦 (𝑛−1) +𝑎2 𝑦 (𝑛−2) + ⋯ + 𝑎𝑛−1 𝑦´ + 𝑎𝑛 𝑦= 𝑔 𝑥 … (∗)


se llama ECUACION D IFERENCIAL ORDINARIA LINEAL
CON COEFICIENTES CONSTANTES
23/06/2021 Mg NOLAN JARA 2
JARA
Si en (*) g(x) = 0, entoncesla Ecuacion

𝑎0 𝑦 (𝑛) + 𝑎1 𝑦 (𝑛−1) +𝑎2 𝑦 (𝑛−2) + 


⋯ + 𝑎𝑛−1 𝑦´ + 𝑎𝑛 𝑦 = 0 
… (∗ 1)

se llama ECUACION D IFERENCIAL LINEAL


HOMOGENEA DE COEF ICIENT ES
CONSTANTES

23/06/2021 Mg NOLAN JARA 3


JARA
POLINOMIO CARACTERISTICO Y ECUACION
CARACTERISTICA
Dada la ecuación diferencial lineal homogénea de
coeficientes constantes
𝑎0 𝑦 (𝑛) + 𝑎1 𝑦 (𝑛−1) +𝑎2 𝑦 (𝑛−2) + ⋯ + 𝑎𝑛−1 𝑦´ + 𝑎𝑛 𝑦 = 0
𝑦 = 𝑒 𝑟𝑥 → 𝑒 𝑟𝑥 𝑎0 𝑟 𝑛 + 𝑎1 𝑟 𝑛−1 + ⋯ + 𝑎𝑛−1 𝑟 + 𝑎𝑛 = 0

1º )El Polinomio
𝑃𝑛 𝑥 = 𝑎0 𝑟 𝑛 + 𝑎1 𝑟 𝑛−1 + 𝑎2 𝑟 𝑛−2 + ⋯ + 𝑎𝑛−1 𝑟 + 𝑎𝑛

es llamado Polinomio Caracteristico

23/06/2021 Mg NOLAN JARA 4


JARA
2º )A la Ecuacion
 + 𝑎𝑛−1 𝑟 + 𝑎𝑛 = 0 … ∗ 2
𝑎0 𝑟 𝑛 + 𝑎1 𝑟 𝑛−1 +𝑎2 𝑟 𝑛−2 + ⋯

Se le llama Ecuacion Caracteristica

23/06/2021 Mg NOLAN JARA 5


JARA
CLASES DE RAICES DE LA ECUACION CARACTERISTICA
1º)Si las n raíces r1 , r2 , r3 
, … , rn de la ecuación caracteristica
son reales y distintas, las soluciones basicas de la ecuacion
𝑎0 𝑦 (𝑛) + 𝑎1 𝑦 (𝑛−1) +𝑎2 𝑦 (𝑛−2) + ⋯ + 𝑎𝑛−1 𝑦´ + 𝑎𝑛 𝑦 = 0

son : e r1 x , e r 2 x , e r 3,
x
,… , e r n x , la soluciongeneral es :
y = c1e r1 x + c 2 e r2 x + c 3 e r3x + … + c n e rn x
2º )Si las raices de la ecuacion caracteristica son reales pero
algunas de ellas se repiten.
Supongamos que r1 = r2 = r3 = … = rk =  y las demas n - k
raíces sean reales y distintas
en (*2): (r -  ) (r - r k+1 )(r- r k +2 ) … (r - rn ) = 0
k

las soluciones Basicas de la ecuación


23/06/2021 Mg NOLAN JARA 6
JARA
𝑎0 𝑦 (𝑛) + 𝑎1 𝑦 (𝑛−1) +𝑎2 𝑦 (𝑛−2) + ⋯ + 𝑎𝑛−1 𝑦´ + 𝑎𝑛 𝑦 = 0


son : e  x , xe x , x 2 e  x ,…
,, x k −1 e  x , e r k +1 x
,er
k+2x
,… , e r x
n

y la solución general es :
𝛼𝑥 𝛼𝑥 2𝑒 𝛼𝑥
𝑦 = 𝑐1 𝑒 + 𝑐2 𝑥𝑒 + 𝑐3 𝑥 + ⋯ + 𝑐𝑘 𝑥 𝑘−1 𝑒 𝛼𝑥 +
𝑐𝑘+1 𝑒 𝑟𝑘+1

𝑥 +𝑐
𝑘+2 𝑒 𝑟𝑘+2 𝑥 + ⋯ + 𝑐 𝑒 𝑟𝑛 𝑥
𝑛
3º )Si alg unas de las raíces de la ecuación caracteristica
son complejas. Supongamos que :
r1 = a + ib, r2 = a − ib, r3 = c + id, r4 = c − id y las demas
(n − 4) raices r5 , r6 , r7 , … , rn sean reales y distintas  las
soluciones basicas de laecuacion

𝑎0 𝑦 (𝑛) + 𝑎1 𝑦 (𝑛−1) +𝑎2 𝑦 (𝑛−2) + ⋯ + 𝑎𝑛−1 𝑦´ + 𝑎𝑛 𝑦 = 0



son :
23/06/2021 Mg NOLAN JARA 7
JARA
ea x cosbx, ea x senbx
ec x cosdx, ec x sendx
e r5 x , e r6 x , 
… , e rn x

y la solución general es :
y = c ea x cosbx + c ea x senbx+ c ec x cosdx + c ec x sendx+
1 2 3 4

c5e r5 x + c6e r6 x + 
… +c e n
n
r x

y= 𝑐1 𝑒 𝑎𝑥 𝑐𝑜𝑠𝑏𝑥 + 𝑐2 𝑒 𝑎𝑥 𝑠𝑒𝑛𝑏𝑥 + 𝑐3 𝑒 𝑐𝑥 𝑐𝑜𝑠𝑑𝑥 + 𝑐4 𝑒 𝑐𝑥 𝑠𝑒𝑛𝑑𝑥 +


𝑐5 𝑒 𝑟5 𝑥 + 𝑐6 𝑒 𝑟6 𝑥 + ⋯ + 𝑐𝑛 𝑒 𝑟𝑛 𝑥

23/06/2021 Mg NOLAN JARA 8


JARA
4º ) Si r1 = a + ib es una raiz que se repite k veces 
n
r2 = a − ib tambien se repite k veces, k  
2
las soluciones basicas son :

e a x cosbx, e a xsenbx, xe a x cosbx,xe a x senbx, … , x k −1e a x cosbx,


k −1 a x 𝑟 𝑥 𝑟 𝑥 
x e senbx,𝑒 2𝑘+1 ,…,𝑒 𝑛

y la solución general es :

y = 𝑐1 e a x cosbx + 𝑐2 e a x senbx+ 𝑐3 xe a x cosbx + 𝑐4 xe a x senbx+...+


𝑐2𝑘−1 x k −1ea x cosbx + 𝑐2𝑘 x k −1ea x senbx + 𝑐2𝑘+1 𝑒𝑟2𝑘+1𝑥 +...+ 𝑐𝑛 𝑒 𝑟𝑛 𝑥

23/06/2021 Mg NOLAN JARA 9


JARA
1) y´´´´´´ - y´´ = 0
𝑟6 − 𝑟2= 𝑟2 𝑟 + 1 𝑟 − 1 𝑟 + 𝑖 𝑟 − 𝑖 = 0

2) x´´´´+ 2x´´+ x=0

23/06/2021 Mg NOLAN JARA 10


JARA
ECUACIONES DIFERENCIALES LINEALES NO
HOMOGENEAS DE COEFICIENTES CONSTANTES
𝑎0 𝑦 (𝑛) + 𝑎1 𝑦 (𝑛−1) +𝑎2 𝑦 (𝑛−2) + ⋯ + 𝑎𝑛−1 𝑦´ + 𝑎𝑛 𝑦 = 𝑔 𝑥

1º Se halla la solución complementaria (solución general de la Ec. Homogenea)
2º Se halla la solución particular de la Ec. no Homogenea
3º Se halla la solución general sumando la solución complementaría
mas la solución particular. yg = yc + yp
DISCUSION PARA HALLAR LA SOLUCION PARTICULAR: METODO DE SELECCIÓN
O COEFICIENTES INDETERMINADOS

1.1) g(x)=𝑃𝑚 𝑥 , 𝑟 ≠ 0; 𝑦𝑝 = 𝐴0 𝑥 𝑚 + 𝐴1 𝑥 𝑚−1 + 𝐴2 𝑥 𝑚−2 + ⋯ + 𝐴𝑚−1 𝑥 + 𝐴𝑚

g(x)=𝑃𝑚 𝑥 , 𝑟 = 0; 𝑦𝑝 = 𝑥 𝑠 (𝐴0 𝑥 𝑚 + 𝐴1 𝑥 𝑚−1 + 𝐴2 𝑥 𝑚−2 + ⋯ + 𝐴𝑚−1 𝑥 + 𝐴𝑚 )



s: grado de la raiz r = 0 de la ecuación caracteristica
: orden

23/06/2021 Mg NOLAN JARA 11


JARA
1.2)

1.3)

1.4)

23/06/2021 Mg NOLAN JARA 12


JARA
+
2) g(x)=𝑃𝑚 𝑥 𝑒 𝑎𝑥 , 𝑎 ∈ 𝑅, 𝑟 ≠ 𝑎; 𝑦𝑝 = 𝑒 𝑎𝑥 (𝐴0 𝑥 𝑚 + 𝐴1 𝑥 𝑚−1 + 𝐴2 𝑥 𝑚−2 + ⋯ + 𝐴𝑚−1 𝑥 + 𝐴𝑚 )

g(x)=𝑃𝑚 𝑥 𝑒 𝑎𝑥 , 𝑎 ∈ 𝑅, 𝑟 = 𝑎; 𝑦𝑝 = 𝑥 𝑠 𝑒 𝑎𝑥 (𝐴0 𝑥 𝑚 + 𝐴1 𝑥 𝑚−1 + 𝐴2 𝑥 𝑚−2 + ⋯ + 𝐴𝑚−1 𝑥 + 𝐴𝑚 )

s: grado de la raiz r = a de la ecuacion caracteristica

ത ത
orden𝑛 𝑥 cosbx+ 𝑄𝑚 𝑥 senbx, 𝑟 ≠ ±𝑏𝑖;𝑦𝑝 = 𝑃𝑘 𝑥 𝑐𝑜𝑠𝑏𝑥 + 𝑄𝑘 𝑥 𝑠𝑒𝑛𝑏𝑥
3) g(x)=𝑃
g(x)=𝑃𝑛 𝑥 cosbx+ 𝑄𝑚 𝑥 senbx, 𝑟 = ±𝑏𝑖;𝑦𝑝 = 𝑥 𝑠 ( 𝑃ത𝑘 𝑥 𝑐𝑜𝑠𝑏𝑥 + 𝑄ത𝑘 𝑥 𝑠𝑒𝑛𝑏𝑥 )
s:grado de la raiz r = ib de la ecuacion característica;k = maxm,n

4) g(x)=𝑒 𝑎𝑥 (𝑃𝑛 𝑥 cosbx+ 𝑄𝑚 𝑥 senbx), 𝑟 ≠ (𝑎 ± 𝑏𝑖);𝑦𝑝 = 𝑒 𝑎𝑥 (𝑃ത𝑘 𝑥 𝑐𝑜𝑠𝑏𝑥 + 𝑄ത𝑘 𝑥 𝑠𝑒𝑛𝑏𝑥)


orden
g(x)=𝑒 𝑎𝑥 (𝑃𝑛 𝑥 cosbx+ 𝑄𝑚 𝑥 senbx), 𝑟 = (𝑎 ± 𝑏𝑖);𝑦𝑝 = 𝑥 𝑠 𝑒 𝑎𝑥 (𝑃ത𝑘 𝑥 𝑐𝑜𝑠𝑏𝑥 + 𝑄ത𝑘 𝑥 𝑠𝑒𝑛𝑏𝑥)

s : grado de la raiz r = (a  ib) de la ecuacion caracteristica ; k = maxm, n

23/06/2021 Mg NOLAN JARA 13


JARA
vi
vi

vi

vi

23/06/2021 Mg NOLAN JARA 14


JARA
vi

vi

vi

vi

23/06/2021 Mg NOLAN JARA 15


JARA
vi
vi

vi

23/06/2021 Mg NOLAN JARA 16


JARA
DETERMINAR LA SOLUCIÓN GENERAL DE LAS
SIGUIENTES ECUACIONES

x + x − 2x = e2t sen2t x + 2x − x = e−t cos2t


xiv + 3x−+4x
4x = − 4 e t (3cos2t + 5sen2t) x iv + 5x + 4x = 2et (5t 2 + 14t + 11)

x iv + 2x + x = cost y´´ + 2y  + 4 y = 2senhx

u´´ + 𝑤02 u = cosw 0 t y ´´ + 9 y = x2 e3x + 6


y ´ ´ + iy  + 2 y = 0 y ´´ + 2 y + iy = 0

y´ ´ − 2 y = ex senx y ´´ − 3y  + 2 y = 14sen2x − 18cos2x

y´´+ 4 y = x2 + 3ex ; y(0) = 0; y(0) =2 y´´ + 4 y = 2 cos x cos3x


− 9 y  + 18y = e e
−3 x
y ´´
23/06/2021 Mg NOLAN JARA 17
JARA
−7 1
y´´ − y − 2 y = cos x − sen2x; y(0) = , y(0) =
20 5

y´´ + y  − 12y = e x+e 2 x − 1; y(0) = 1; y(0) = 3

y´´ − y = senx − e 2 x ; y(0) = 1; y(0) = −1

y´´ − 7 y  + 10y = x 2 − 4 + e x ; y(0) = 3; y(0) = −3

23/06/2021 Mg NOLAN JARA 18


JARA
WRONSKIANO

Supongamos que cada una de las funciones f1(x), f2(x), …, fn(x)


posee al menos n – 1 derivadas. El determinante

f1 f2 
… fn
f1 ' f2 ' 
… fn '
W ( f1,..., f n ) =

… 
… …

f (n−1)
1
f (n−1)
2

… f n(n−1)

se llama el Wronskiano de las funciones.

23/06/2021 Mg NOLAN JARA 19


JARA
CRITERIO PARA SOLUCIONES LINEALMENTE INDEPENDIENTES
Sean y1(x), y2(x), …, yn(x) soluciones de una ED homogénea de n-ésimo orden en
un intervalo I. Este conjunto de soluciones es linealmente independiente si y sólo
si W(y1, y2, …, yn)  0 para todo x en el intervalo.
CONJUNTO FUNDAMENTAL DE SOLUCIONES
Cualquier conjunto y1(x), y2(x), …, yn(x) de n soluciones linealmente
independientes de una ED homogénea de n-ésimo orden se llama conjunto
fundamental de soluciones.
EXISTENCIA DE UN CONJUNTO FUNDAMENTAL
Existe un conjunto fundamental de soluciones para una ED lineal homogénea de
orden n en un intervalo I.
SOLUCIÓN GENERAL (ECUACIONES HOMOGÉNEAS)
Sea y1(x), y2(x), …, yn(x) un conjunto fundamental de soluciones de nuestra ED
lineal homogénea en un intervalo I. Entonces la solución general es
y = c1y1(x) + c2y2(x) + … + cnyn(x) donde
ci son constantes arbitrarias.

23/06/2021 Mg NOLAN JARA 20


JARA
SOLUCIÓN GENERAL ECUACIONES EDOL NO HOMOGENEA

Sea yp cualquier solución particular de una EDO no homogénea en un


intervalo I.
Sean y1(x), y2(x), …, yk(x)
un conjunto fundamental de soluciones de su EDO homogénea asociada,
entonces la solución general de la ecuación en el intervalo es

y= c1y1 + c2y2 +… + ckyk + yp


donde las ci , i= 1,2,….,n son constantes arbitrarias

y = c1y1 + c2y2 +… + ckyk + yp = yc + yp


= solución complementaria + una solución particular

23/06/2021 Mg NOLAN JARA 21


JARA
METODO VARIACION DE PARAMETROS
Consideremos la EDO lineal completa
a 0 y (n) + a1 y (n−1) + a 2 y (n−2) + … + an−1 y  + a n y = g(x)
De la solución complementaria
yc = c1 y1 (x) + c 2 y 2 (x)+ … + c n y n (x)
El metodo de variacion de parámetros consiste en el reemplazo de las
constantes 𝒄 𝒌 por funciones 𝒖𝒌(𝒙). k=1,2,3,…,n
y = u (x) y (x) + u (x) y (x) + + u (x) y (x)
p 1 1 2 2 … n n

Buscamos las funciones u(x) resolviendo el sistema


u1 (x) y1 (x) + u 2  (x) y 2 (x) + …+ u n (x) y n (x) = 0
u  (x) y  (x) + u  (x) y  (x) +  … + u  (x) y  (x) = 0
1 1 2 2 n n

...........................................................................
 ( n −1)  +  (x) y (n−1) (x) = g(x)
(x) + u 2 (x) y 2 (x) ++… u
(n−1)
u1 (x) y 1 n n

23/06/2021 Mg NOLAN JARA 22
JARA
La soluciongeneral de la E D O lineal
a 0 y ( n ) + a 1 y ( n−1) + a 2 (n − 2 ) +
+… + a n −1 y  + a n y = g ( x )
sera : y = yc + y p
Ejemplo : Re solver y + y = sec x
Solucion
yc = A cos x + Bsenx ; yp = u(x) cos x + v(x)senx
 u(x) cos x + v(x)senx = 0
− u(x)senx + v(x) cos x = sec x

0 sen x cos x 0
sec x cos x − tgx  − senx sec x 1
u(x) = = ; v (x) = =
cos x sen x 1 cos x sen x 1
− sen x cos x − sen x cos x
cos x sen x
w ( y1 , y 2 ) = Wr o n s k i a n o d e y 1 y y 2
− sen x cos x
u (x) = ln cos x ; v(x) = x
y p = ln cos x cos x + xsenx
y = A cos x + Bsen x + ln cos x cos x + xsenx
23/06/2021 Mg NOLAN JARA 23
JARA
1
Ejemplo2.Re solver: y  + 4 y  + 4 y = 2 2x
; x 0
x e
y  + 9 y = 9 sec 2 3x;0  x 

6
y  + y = tgx
y  + 4y = 3cosec2x
y  − 4 y  + 4 y = (x + 1)e2 x

23/06/2021 Mg NOLAN JARA 24


JARA
METODO DE LOS OPERADRES DIFERENCIALAES

Consideremos la EDO lineal completa


a 0 y (n) + a1 y (n−1) + a 2 y (n−2) +... + a n −1 y  + a n y = g(x)… (a)
Medianteun conveniente cambio de notacion, escribiendo
y(x) = Dy, y  (x) = D.Dy = D 2 y, ..., y (n) (x) = D n y
(a) toma la forma
a D n

+ a1 D n −1 + a 2 D n − 2 + + a n −1 D + a n y = g(x)… (b)
+a D
0

P(D) = a 0 D n 1
n −1
+ a 2 D n − 2 +... + a n −1 D + a n
Si; P(D) = (D − m1 )(D − m 2 )… (D − m n )
en (b) (D − m1 )(D − m 2 )… (D − m n ) y = g(x)
1
 P(D)y = g(x)  y = g (x)
P(D)

23/06/2021 Mg NOLAN JARA 25


JARA
1 1 1
 y=  .... g(x)… (c)
(D − m1 ) (D − m 2 ) (D −m n )
(c) se resuelve como una sucesión de ecuaciones diferenciales lineales
de la forma
1 1 1
u1 = g(x); u 2 = u1;… ; y = u n−1
(D − m n ) (D − m n −1 ) (D − m1 )
tambien, si :
 A1 A2 An 
y= + +...+  g(x)… (c)
 (D − m1 ) + (D − m 2 ) (D − m n ) 
(c) se resuelve como una suma de ecuaciones diferenciales lineales
(c) Solución particular de la EDO (a)
PROBLEMAS
1) y + 2y - y - 2y = e 2x 2) y - 3y + 2y = e x
3) y + 3y - 4y = xe -2x 4) y + 5y + 4y = 3 − 2x
5) y - 5y + 8y - 4y = e 2x 6) y - 3y + 2y = sene -x
23/06/2021 Mg NOLAN JARA 26
JARA
ECUACIONES DIFERENCIALES LINEALES CON
COEFICIENTES VARIABLES

a 0 (x) y (n) + a1 (x) y (n−1) + a 2 (x) y (n−2) + + a n−1 (x) y + a n (x) y = g(x)
a0 (x), a1 (x), a2 (x), , an−1 (x), an (x),g(x) son funciones que
dependen de x.
n−1
d n
y d y dy
y = n ,y
(n) (n−1)
= n−1 , , y  = ; y = y(x)
dx dx dx
Si n = 2  a0 (x) y´´ + a1 (x) y + a2 (x) y = g(x)

23/06/2021 Mg NOLAN JARA 27


JARA
ECUACIÓN DE CAUCHY-EULER
FORMA DE LA ECUACIÓN DE CAUCHY-EULER
n n−1
d y d y dy
+a
n−1
ax
n
n
x
n−1 n−1
+ + a x + a0 y = g(x)
1
dx n dx dx
MÉTODO DE SOLUCIÓN 1: la sustitución x = et
Convierte la ecuación de Cauchy-Euler en una ecuación diferencial lineal
de coeficientes constantes
d ny d n−1 y dy
n + n−1
bn b n−1
+  + b 1 + b 0 y = G(t)
dt dt dt
1
y( x) = y(t)
x
1
y( x) = 2
( y(t ) − y(t))
x
1
y    (x) = 3
( y(t ) − 3 y(t ) + 2 y(t ))
x
23/06/2021 Mg NOLAN JARA 28
JARA
RESOLOVER : x 2 y ´´ − xy + y = ln x
Solucion
x = e t  t = ln x
dy dy dt 1 dy
= =
dx dt dx x dt
d 2 y d  1 dy  1 dy 1  d  dy   1 dy 1  d  dy  
=  =− 2 +    = − 2 +    =
dx 2
dx  x dt  x dt x  dx  dt   x dt x  dt  dx  
1 dy 1  d  1 dy   1  d 2 y dy 
− 2 +    = 2 − 
x dt x  dt  x dt   x  dt 2
dt 

23/06/2021 Mg NOLAN JARA 29


JARA
x 2 y´´− xy + y = ln x
d2y dy
−2 + y =t
dt 2 dt
y = c1et + c2tet + 2 + t
y = c1x + c2 x ln x + 2 + ln x
RESOLVER : x y ´´ − xy + 2 y = x ln x
2

23/06/2021 Mg NOLAN JARA 30


JARA
ECUACIÓN DE CAUCHY-EULER
OTRA FORMA DE LA ECUACIÓN DECAUCHY-EULER

n−1
d ny n−1 d y dy
an (ax +b) n
+ an−1 (ax + b) n−1
+ + a1 (ax + b) + a0 y = g(x)
dx n dx dx

la sustitución ax + b = e t
También convierte la ecuación de Cauchy-Euler en una ecuación diferencial lineal de
coeficientes constantes

RESOLOVER : (x + 2)2 y´´ + 3(x + 2) y − 3y = 0


(2x +1)2 y´´ − 2(2x +1)y + 4y = 0

23/06/2021 Mg NOLAN JARA 31


JARA
ECUACIÓN DE CAUCHY-EULER
FORMA DE LA ECUACIÓN DE CAUCHY-EULER
n n−1
d y d y dy
+ a n−1 x + + a1 x + a0 y = 0
n n−1
an x n−1
dx n
dx dx
MÉTODO DE SOLUCIÓN 2: Se buscan soluciones de la forma y = x
m

Obteniéndose para m una ecuación que coincide con La ecuación característica de


la ecuación n n−1
d y d y dy
bn
dt n
+b n−1
dt n−1
+  + b1 + a 0 y = 0
dt
Observación.

d ky
ak x k
k
= ak x k
m(m −1)(m − 2) (m − k +1)x m−k

dx
= akm(m −1)(m − 2) (m − k +1)x m

23/06/2021 Mg NOLAN JARA 32


JARA
Re solver : x 2 y + 2xy (x) − 6 y = 0 (E)
solucion
Metodo1 : x= e t ⎯(⎯E)→ y(t) + y(t) − 6 y = 0
A
 y = 3 + Bx 2
x
Metodo2 : y =x m ⎯(⎯
E)
→m(m −1) + 2m − 6x m = 0
A
 y = 3 + Bx 2
x

23/06/2021 Mg NOLAN JARA 33


JARA
Re s o l v e r : 2 x 2 y  + 3 xy  − y = 0 … ( E )
so lu cion
y = x m ⎯ (⎯
E)

→ 2m 2 + m − 1 x m = 0 
B
 y = A x + ;x  0
x
Re s o l v e r : x 2 y  + 5 x y  + 4 y = 0  ( E )
Re s o l v e r : x 2 y  + xy  + y = 0  ( E )

23/06/2021 Mg NOLAN JARA 34


JARA
MODELADO CON ECUACIONES DIFERENCIALES DE
ORDEN SUPERIOR
considerar sistemas dinámicos lineales en los que cada modelo
matemático es una Ecuación Diferencial de Segundo Orden con
Coeficientes Constantes junto con condiciones iniciales especificadas
en un tiempo que tomaremos como t=0

Ax´´(t) + bx(t) + cx(t) =f(t); x( 0 ) = x0 , x( 0 ) = x1


f: Es la entrada, función de conducción o función forzada del sistema.
Una solución x(t) de la EDO en un intervalo I que contiene a t=0 que
satisface las condiciones iniciales se llama Salida o Respuesta del
Sistema.

23/06/2021 Mg NOLAN JARA 35


JARA
SISTEMA RESORTE-MASA: MOVIMIENTO LIBRE NO AMORTIGUADO

x(0): cantidad de
desplazamiento inicial
x´(0): rapidez inicial de la
masa
x(0)>0: la masa parte de un
X=0 Punto debajo de la posición
de equilibrio
x(0)<0: la masa parte de un
punto arriba de la posición de
equilibrio

x´(0)>0: rapidez inicial de la masa hacia abajo


x´(0)<0: rapidez inicial de la masa hacia arriba
x´(0)=0: la masa se suelta partiendo del reposo
23/06/2021 Mg NOLAN JARA 36
JARA
SISTEMA RESORTE-MASA:
MOVIMIENTO LIBRE NO
AMORTIGUADO
LEY DE HOOKE (LH): A un resorte suspendido verticalmente
de un cuerpo rígido se le fija una masa m en su extremo libre.
El alargamiento del resorte depende de la masa m. Por la LH
el resorte ejerce una fuerza restauradora F opuesta a la
dirección del alargamiento y proporcional a la cantidad
alargada S. F = kS; K constante de proporcionalidad
( Constante de Resorte )

23/06/2021 Mg NOLAN JARA 37


JARA
Movimiento libre no amortiguado

Ley de Hooke:
F=ks
F: fuerza restauradora
k= Constante de resorte
2da ley de Newton:
d 2x
m 2 = −k(s + x) + mg =−kx
dt
PVI del movimiento:
x ´´ (t) + kx = 0; x(0) = x 0 , x  (0) = x1
23/06/2021 Mg NOLAN JARA 38
JARA
SISTEMA RESORTE-MASA: MOVIMIENTO
LIBRE AMORTIGUADO
El concepto de movimiento armónico libre no es realista por que
el movimiento que describe la ecuación (I) supone que no hay
fuerzas retardadoras actuando sobre la masa m en movimiento.
A menos que la masa se suspenda en un vacío perfecto, habrá
por lo menos una fuerza de resistencia debida al medio
circundante. La masa podría estar suspendida en un medio
viscoso o conectada a un dispositivo amortiguador.

23/06/2021 Mg NOLAN JARA 39


JARA
SISTEMA RESORTE-MASA: MOVIMIENTO LIBRE
AMORTIGUADO
ECUACIÓN DIFERENCIAL DE UN MOVIMIENTO LIBRE
AMORTIGUADO
En mecánica, las fuerzas de amortiguamiento que actúan sobre
un cuerpo se consideran proporcionales a una potencia de la
velocidad instantánea. En particular, en el análisis posterior se
supone que esta fuerza está expresada por un múltiplo constante
de x’(t).
Cuando no hay otras fuerzas externas aplicadas al sistema, se
sigue por la segunda ley de Newton:

mx ´´ (t) = −kx −  x  (t)


donde  es una constante de amortiguamiento positiva y el
signo negativo es consecuencia del hecho de que la fuerza
amortiguadora actúa en dirección opuesta a la del movimiento.
23/06/2021 Mg NOLAN JARA 40
JARA
Movimiento libre amortiguado
•En un resorte, se observa que la amplitud de la vibración
disminuye con el tiempo, ya que hay una pérdida de energía.
Se dice que la oscilación está amortiguada.
•Para el análisis del oscilador dinámico, se puede suponer
que además de la fuerza elástica, también actúa una
fuerza disipativa que se opone a la velocidad, de la forma

Fd =  v
2da Ley de Newton:
dx d2 x
− kx − v = ma, entonces − kx−  =m 2
dt dt
PVI deI mov. libre amortiguado:
mx ´´ (t) +  x  (t) + kx = 0; x(0) = x0 , x(0) = x1
23/06/2021 Mg NOLAN JARA 41
JARA
SISTEMA RESORTE-MASA: MOVIMIENTO LIBRE
AMORTIGUADO
ECUACIÓN DIFERENCIAL DE UN MOVIMIENTO LIBRE
AMORTIGUADO

mx ´´ (t) +  x  (t) + kx = 0
 2 − 4mk  0 El sistema está sobre
amortiguado

 2 − 4mk = 0 El sistema está críticamente amortiguado

 2 − 4mk  0 El sistema está subamortiguado

23/06/2021
06/07/2020 Mg NOLAN
NOLAN JARAJARA
JARA 42
39
JARA
EJEMPLO 1

Un objeto de 5 kg de masa estira un resorte


0,2 m sobre su longitud natural.
Si el resorte se estira hasta medir 0,6 m más
que su longitud en el equilibrio y luego el
objeto se suelta con velocidad inicial de 0,7
m/s hacia arriba, encuentre la posición del
objeto en el tiempo t.
Considere que no hay fuerzas amortigüadoras.
23/06/2021 Mg NOLAN JARA 43
JARA
Solución.

mg = ks  5(9.8) = k (0.2)  k = 245


mx = −  x  − kx  5x  = −0x  − 245x
x  + 49x = 0; x(0) = 0.6, x  (0) = −0.7
m 2 + 49 = 0  m = 7i
x(t) = A cos(7t) + Bsen(7t) (i);0.6 = x(0) = A  A = 0.6
x  (t) = −7 Asen(7t) + 7B cos(7t);−0.7 = x  (0) = 7B  B = −0.1
6 1
en i : x(t) = cos(7t) − sen(7t)
10 10

23/06/2021 Mg NOLAN JARA 44


JARA
EJEMPLO 2

Suponga que el sistema del Ejemplo 1


se sumerge en un fluido con constante
de amortiguamiento  = 70. Halle la
posición del objeto en el instante t, si
este parte de la posición de equilibrio y
recibe un empujón hacia abajo que le
imparte una velocidad inicial de 0,6 m/s.

23/06/2021 Mg NOLAN JARA 45


JARA
Solución.

mg = ks  5(9.8) = k (0.2)  k = 245


mx = −x − kx  5x = −70x − 245x
x + 14x + 49x = 0; x(0) = 0, x(0) = 0.6
m +
2 14m + 49 = 0  (m + 7) 2
= 0  m = −7
2

x(t) = Ae −7t + Bte−7t  (i);0 = x(0) = A  A = 0


x(t) = −7 Ae −7t +Be −7t − 7Bte −7t
0.6 = x(0) = B  B = 0.6
6 −7t
en i : x(t) = te
23/06/2021
10 Mg NOLAN JARA 46
JARA
EJEMPLO 3

Un objeto de 10 kg de masa se une a un


resorte de 1,50 m de longitud natural, que en
el equilibrio mide 1,99 m. Si el resorte se
comprime hasta medir 1,70 m y luego el
objeto se suelta a una velocidad de 1 m/s
hacia abajo, encuentre la posición del objeto
en cualquier tiempo t. Considere que el
resorte se encuentra en un medio que ofrece
una resistencia numéricamente igual a 80
veces su velocidad instantánea en m/s.
23/06/2021 Mg NOLAN JARA 47
JARA
Mg. Nolan Jara Jara 1
Mg. Nolan Jara Jara 2
Mg. Nolan Jara Jara 3
Mg. Nolan Jara Jara 4
Mg. Nolan Jara Jara 5
Mg. Nolan Jara Jara 6
LA TRANSFORMADA DE LAPLACE
Sea f(t) una función definida para t ≥ 0, la transformada de Laplace
de f(t) se denota y define como:

∞ 𝑏
𝐿 𝑓 𝑡 = න 𝑒 −𝑠𝑡 𝑓 𝑡 𝑑𝑡 = lim න 𝑒 −𝑠𝑡 𝑓 𝑡 𝑑𝑡
0 𝑏→∞ 𝑡=0

Donde sR
Se dice que la transformada de Laplace de f(t) existe si la integral
converge para algunos valores de s.

Mg. Nolan Jara Jara 7


LA TRANSFORMADA DE LAPLACE ES UNA INTEGRAL
IMPROPIA, UNO DE SUS LÍMITES ES INFINITO:

∞ 𝑏
න 𝑒 −𝑠𝑡 𝑓 𝑡 𝑑𝑡 = lim න 𝑒 −𝑠𝑡 𝑓 𝑡 𝑑𝑡
0 𝑏→∞ 𝑡=0

Notación: L  f (t ) = F ( s ),
L  y (t ) = Y ( s ),
L  x(t ) = X ( s ), etc.
L f (t ) = F ( s)  L {F ( s)} = f (t )
−1

Mg. Nolan Jara Jara 8


CONDICIONES SUFICIENTES DE EXISTENCIA
DE LA TRANSFORMADA DE LAPLACE

F ( s ) = L{ f (t )} = 
t =0
f (t )e − st dt

Si f(t) es continua a trozos en [0, ∞) y

e −t f (t )  M o | f (t ) | Met , t  [0, )


Es decir, f(t) es de orden exponencial  en el infinito:

  R tal que lim | f (t )e − t


|= 0
t →
Entonces: L{f(t)} = F(s) existe s > 

lim F(s)=0
𝑏→∞
Mg. Nolan Jara Jara 9
Hallar la transformada de Laplace de: f(t) = 1
∞ b
−1 −𝑠𝑡 1
𝐿 1 =න 1𝑒 −𝑠𝑡 𝑑𝑡 = lim 𝑒 ቤ = = 𝐹(𝑠); 𝑠 > 0
0 𝑏→∞ 𝑠 𝑠
t= 0

1
∴ 𝑆𝑖 𝑓 𝑡 = 1 ⇒ 𝐿 1 = = 𝐹 𝑠 ; 𝑠 > 0.
𝑠

𝐴
𝑂𝑏𝑠𝑒𝑟𝑣𝑎𝑐𝑖𝑜𝑛: 𝐿 𝐴 = , 𝑠 > 0; 𝐿 0 = 0
𝑠

1
𝐿 𝑡 = 2 = 𝐹 𝑠 ;𝑠 > 0
𝑠
Mg. Nolan Jara Jara 10
Hallar la transformada de Laplace de: f(t) = t n
− st 
e − st
 
L t n = F ( s) = 

0
n − st
t e dt = t n e

−s

−  nt n −1
0 −s
dt
0

n  n −1 − st n
=  t e dt = L t n −1
s 0 s
 

 
L t n n
= L t
s
 
n −1 n!
 L t = n +1
n

s
 

𝒏 + 𝒏
𝒏!
𝑺𝒊 𝒇 𝒕 = 𝒕 ; 𝒏 ∈ 𝒁𝟎 ⇒𝑳 𝒕 = = 𝑭 𝒔 ;𝒔 > 𝟎
𝒔𝒏+𝟏

Mg. Nolan Jara Jara 11


Hallar la transformada de Laplace de: f(t) = e- t

F ( s) = L e = 
−t

0
e e dt =  e −( s +1)t dt
−t − st

0

− 1 −( s +1)t 1
e =
s +1 0 s +1

1
∴ 𝑆𝑖 𝑓 𝑡 = 𝑒 −𝑡 ⇒𝐿 𝑒 −𝑡 = = 𝐹(𝑠); 𝑠 > −1
𝑠+1

Si f (t ) = e  F ( s ) = L e
t
 
t
=
1
s −1
;s 1

Mg. Nolan Jara Jara 12


Hallar la transformada de Laplace de: f(t) = eat

F ( s) = L e  = at

0
e e dt =  e
at − st

0
− ( s − a )t
dt
− ( s − a )t 
−e 1
= , sa
( s − a) 0
s−a

𝒂𝒕 𝒂𝒕
𝟏
∴ 𝑺𝒊 𝒇 𝒕 = 𝒆 ⇒𝑳 𝒆 = = 𝑭(𝒔); 𝒔 > 𝒂
𝒔−𝒂

Mg. Nolan Jara Jara 13


Hallar la transformada de Laplace de: f(t) = sen(at)
 e − st  e − st
F ( s ) = L sen(at) =  − st
sen(at)e dt = sen(at) − a cos(at) dt
0 −s 0
0 −s

a e − st 

 e − st   a 
cos(at) − (− a ) sen(at) a
dt  = 2 − 2  sen(at)e − st dt =
a
s −s 0 −s  s s 0 s2 + a2
 0 
𝒂
𝑺𝒊 𝒇 𝒕 = 𝒔𝒆𝒏 𝒂𝒕 ⇒ 𝑳 𝒔𝒆𝒏(𝒂𝒕) = 𝟐 𝟐
= 𝑭(𝒔); 𝒔 > 𝟎
𝒔 +𝒂
Si f (t ) = sen(t )  F ( s ) = Lsent =
1
;s  0
s +1
2

f (t ) = cos(t )  F ( s ) = Lcos t  =
s
Si ;s  0
s +1
2

𝒔
𝑺𝒊 𝒇 𝒕 = 𝒄𝒐𝒔 𝒂𝒕 ⇒ 𝑳 𝒄𝒐𝒔(𝒂𝒕) = 𝟐 𝟐
= 𝑭(𝒔); 𝒔 > 𝟎
𝒔 +𝒂
Mg. Nolan Jara Jara 14
LA FUNCIÓN HEAVISIDE O ESCALÓN UNITARIO

La función escalón unitario u(t – a) se define Como

0 , 0  t  a
u(t − a ) = 
1 , ta

Mg. Nolan Jara Jara 15


TRANSFORMADA DE LAPLACE DE LA FUNCIÓN
HEAVISIDE O ESCALÓN UNITARIO

0; si t  a
u (t − a) = 
1
1; si t  a
 
Lu (t − a ) =  e − st
u (tc − a )dt =  e − st dt
t
0 a

=
s
e
− 1 − st

t =a =
e − as
s
;s  0

𝒆−𝒂𝒔
∴ 𝑺𝒊 𝒇 𝒕 = 𝒖 𝒕 − 𝒂 ⇒ 𝑳 𝒖 𝒕 − 𝒂 = = 𝑭(𝒔); 𝒔 > 𝟎
𝒔
Mg. Nolan Jara Jara 16
L{u(t-2)-u(t-5)}
0, 𝑠𝑖 𝑡 ∈ [0,2 >∪ [5, ∞ >
Sea: f(t)= u(t-2)-u(t-5)=ቊ
1, 𝑠𝑖 𝑡 ∈ [2,5 >

L {f(t)}= ‫׬‬0 𝑓 𝑡 𝑒 −𝑠𝑡 𝑑𝑡 =
2 5 ∞
න 0. 𝑒 −𝑠𝑡 𝑑𝑡 + න 1. 𝑒 −𝑠𝑡 𝑑𝑡 + න 0. 𝑒 −𝑠𝑡 𝑑𝑡 =
0 2 5
5 −𝑠𝑡 −1 𝑒 −2𝑠 −𝑒 −5𝑠
‫=𝑡׬‬2 𝑒 𝑑𝑡= 𝑠 𝑒 −5𝑠 − 𝑒 −2𝑠
= 𝑠
=L{u(t−2)−u(t−5)}
f(t)

0 t
2 5
Mg. Nolan Jara Jara 17
𝟏
¿𝒈(𝒕) = 𝒖(𝒕 − 𝒂) − 𝒖 𝒕 − (𝒂 + 𝜺) ?
𝜺 Mg. Nolan Jara Jara 18
Función impulso unitario-Función delta de Dirac
0, 𝑠𝑖 𝑡 ∈ [0, 𝑎 >∪ [𝑎 + 𝜀, ∞ >
1
𝑔(𝑡) = 𝑢(𝑡 − 𝑎) − 𝑢 𝑡 − (𝑎 + 𝜀) = ቐ1
𝜀 , 𝑠𝑖 𝑡 ∈ [𝑎, 𝑎 + 𝜀 >
g(t) 𝜀
1
−− −
𝜀

𝑨 𝑹 →𝟏

𝜀→0
0 t
a a+𝜀
Mg. Nolan Jara Jara 19
FUNCION DE IMPULSO UNITARIO - DELTA DE DIRAC
𝟏
𝜹 𝒕 − 𝒂 = 𝒍𝒊𝒎+ 𝒈 𝒕 = 𝐥𝐢𝐦+ 𝒖(𝒕 − 𝒂) − 𝒖 𝒕 − (𝒂 + 𝜺)
𝜺→𝟎 𝜺→𝟎 𝜺
0 ; si t  0,a )  a + ε, )

1
u(t − a) − u(t − (a +  )) =  1 ; si t  a,a + ε )
 
ε

0 ; si t  0,a )  a + ε, )

 (t − a) =  1 ; ε → 0+
; si t  a,a + ε )

ε

0 ; si t  ε, )

 (t ) =  1 ; ε → 0+
; si t  0,ε )

ε
Mg. Nolan Jara Jara 20
FUNCION DE IMPULSO UNITARIO - DELTA DE DIRAC
𝟏
𝜹 𝒕 − 𝒂 = 𝒍𝒊𝒎+𝒈 𝒕 = 𝐥𝐢𝐦+ 𝒖(𝒕 − 𝒂) − 𝒖 𝒕 − (𝒂 + 𝜺)
𝜺→𝟎 𝜺→𝟎 𝜺
𝟏
𝑳 𝜹 𝒕 − 𝒂 = 𝑳 𝒍𝒊𝒎+ 𝒈 𝒕 = 𝑳 𝐥𝐢𝐦+ 𝒖(𝒕 − 𝒂) − 𝒖 𝒕 − (𝒂 + 𝜺)
𝜺→𝟎 𝜺→𝟎 𝜺
𝟏
= 𝐥𝐢𝐦+𝑳 𝒖(𝒕 − 𝒂) − 𝒖 𝒕 − (𝒂 + 𝜺)
𝜺→𝟎 𝜺
𝟏
= 𝐥𝐢𝐦+ 𝜺 𝑳 𝒖(𝒕 − 𝒂) − 𝒖 𝒕 − (𝒂 + 𝜺)
𝜺→𝟎
𝟏
= 𝐥𝐢𝐦+ 𝜺 𝑳 𝒖(𝒕 − 𝒂) − 𝑳 𝒖 𝒕 − (𝒂 + 𝜺)
𝜺→𝟎
𝟏 𝒆−𝒂𝒔 𝒆− 𝒂+𝜺 𝒔 𝟏 𝒆−𝒂𝒔 −𝒆− 𝒂+𝜺 𝒔
= 𝐥𝐢𝐦+ − 𝒔 = 𝐥𝐢𝐦+ 𝜺
𝜺→𝟎 𝜺 𝒔 𝜺→𝟎 𝒔
𝟏−𝒆−𝜺𝒔 𝒆 −𝜺𝒔 −𝟏
= 𝒆−𝒂𝒔 𝐥𝐢𝐦+ 𝜺𝒔 = 𝒆−𝒂𝒔 𝐥𝐢𝐦+ −𝜺𝒔 = 𝒆−𝒂𝒔
𝜺→𝟎 𝜺→𝟎

∴ 𝑳 𝜹 𝒕 − 𝒂 = 𝒆−𝒂𝒔
Mg. Nolan Jara Jara 21
TRANSFORMADA DE LAPLACE DE LA FUNCIÓN
DELTA DE DIRAC

𝑳 𝜹(𝒕 − 𝒂) = 𝒆−𝒂𝒔

𝑳 𝜹(𝒕) = 𝟏  (t )  (t − a )

t
a

Mg. Nolan Jara Jara 22


FUNCION DE IMPULSO UNITARIO -
f  (t )
DELTA DE DIRAC area = 1
Sea la función :
1/ 
f  (t ) = u (t − a ) − u (t − (a +  ) )
1
 t
a a +

lim →0 f  (t ) =  (t − a)
1  e − as e − ( a + ) s  − as 1 − e −s 
L  f  (t ) =  − =e  
 s s   s 
− s − s
 −   
lim →0 L f  (t ) = e lim →0 
− as 1 e − as se − as
 = e lim →0   = e
 s   s 
Mg. Nolan Jara Jara 23
FUNCION DE IMPULSO UNITARIO O DELTA DE DIRAC


න 𝛿(𝑡 − 𝑎)𝑔(𝑡) 𝑑𝑡 = 𝑔 𝑎
0


න 𝛿(𝑡)𝑔(𝑡) 𝑑𝑡 = 𝑔(0)
0


0
 (t − a ) dt = 1
Mg. Nolan Jara Jara 24
Mg. Nolan Jara Jara 25
Mg. Nolan Jara Jara 26
FUNCIONES PERIÓDICAS
Supongamos que f (t) es una función periódica de periodo T.
Entonces:

F ( s ) = L  f (t ) =
1
F (s)
− sT 1
1− e
donde F1(s) es la transformada de Laplace de la función f(t)
sobre el primer periodo y cero fuera.
T
F1 ( s ) =  e − st f (t )dt
0

t Mg. Nolan Jara Jara t 27


¿L 𝒕 − 𝒕 ?
Mg. Nolan Jara Jara 28
PROPIEDADES

1. Linealidad: Si c1 y c2 son constantes, f1(t) y f2(t) son


funciones cuyas transformadas de Laplace son F1(s) y F2(s),
respectivamente; entonces:

L{c1 f1 (t ) + c2 f 2 (t )} = c1L f1 (t )+ c2 L f 2 (t ) = c1F1 (s) + c2 F2 (s)

La transformada de Laplace es un operador lineal.

Mg. Nolan Jara Jara 29


2. Desplazamiento temporal
Si : L{ f (t )} = F ( s ) y
 0, 0  t  a
g (t ) = f (t − a)u (t − a) =  
 f (t − a), t  a
 
G ( S ) =  f (t − a)u (t − a)e dt =  f (t − a)e dt; v = t − a
− st − st
0 a


− as
G(s) = e f (v)e − sv dv = e −as F ( s )
0

𝑳{𝒈(𝒕)} = 𝑳{𝒖(𝒕 − 𝒂)𝒇(𝒕 − 𝒂)} = 𝒆−𝒂𝒔 𝑭(𝒔)

𝑳{𝒖(𝒕 − 𝒂)𝒇(𝒕)} = 𝒆−𝒂𝒔 𝑳{𝒇(𝒕 + 𝒂)}

Mg. Nolan Jara Jara 30


Mg. Nolan Jara Jara 31
Mg. Nolan Jara Jara 32
3. Desplazamiento en frecuencias

Si : L f (t ) =  e − st f (t )dt =F ( s ) 
0
 
  ( )
L e at f (t ) =  e at f (t ) e − st dt =  e −( s − a ) t f (t )dt
0 0

= F ( s − a); s  a

∴ 𝑳{𝒆𝒂𝒕 𝒇(𝒕)} = 𝑭(𝒔 − 𝒂); 𝒔 > 𝒂

Ejemplo:
L  t = 2  L e t =
s
1 at
 
1
( s − a) 2

Mg. Nolan Jara Jara 33


4. Cambio de escala en tiempo

Si : L f (t ) =  e − st f (t )dt = F ( s) 
0
  s
− v
L f (at) =  f (at)e dt =  e f (v)dv; (v = at )
− st1 a

0
a0
1 s
= F( )
a a
𝟏 𝒔
𝑳{𝒇(𝒕)} = 𝑭(𝒔) ⇒ 𝑳{𝒇(𝒂𝒕)} = 𝑭
𝒂 𝒂

Mg. Nolan Jara Jara 34


5. Derivada de la transformada de Laplace

𝑆𝑖: 𝐿 𝑓(𝑡) = න 𝑒 −𝑠𝑡 𝑓(𝑡)𝑑𝑡 = 𝐹(𝑠)


0

𝑑 𝑑 −𝑠𝑡
𝐹(𝑠) = න 𝑒 𝑓(𝑡)𝑑𝑡
𝑑𝑠 𝑑𝑠
0

= − න 𝑒 −𝑠𝑡 𝑡𝑓(𝑡) 𝑑𝑡
0
= −𝐿 𝑡𝑓(𝑡)

𝑺𝒊: 𝑭 𝒔 = 𝑳 𝒇 𝒕 ⇒ 𝑭′ 𝒔 = −𝑳 𝒕𝒇 𝒕 ; 𝑳 𝒕𝒇 𝒕 = −𝑭′ 𝒔

Mg. Nolan Jara Jara 35


Si : F ( s) = L{ f (t )}  L tf (t ) = − F ( s)

𝑺𝒊: 𝑭(𝒔) = 𝑳{𝒇(𝒕)} ⇒ 𝑳 𝒕𝟐 𝒇 𝒕 = (−𝟏)𝟐 𝑭´´(𝒔)

 
Si : F ( s ) = L{ f (t )}  L t f (t ) = (−1) F ( s )
3 3

𝑺𝒊: 𝑭(𝒔) = 𝑳{𝒇(𝒕)} ⇒ 𝑳 𝒕𝒏 𝒇(𝒕) = (−𝟏)𝒏 𝑭(𝒏) (𝒔)

Mg. Nolan Jara Jara 36


Mg. Nolan Jara Jara 37
6. Transformada de Laplace de las derivadas de
una función
La transformada de Laplace de la derivada de una función está dada
por:
𝑳{𝒇′(𝒕)} = 𝒔𝑭(𝒔) − 𝒇(𝟎)

La transformada de Laplace de la segunda derivada de una función


está dada por:

𝑳{𝒇′′(𝒕)} = 𝒔𝟐 𝑭(𝒔) − 𝒔𝒇(𝟎) − 𝒇′(𝟎)

Mg. Nolan Jara Jara 38


En general:
𝑳{𝒇(𝒏) (𝒕)} = 𝒔𝒏 𝑭(𝒔) − 𝒔𝒏−𝟏 𝒇(𝟎) − 𝒔𝒏−𝟐 𝒇′(𝟎) − ⋯ − 𝒇(𝒏−𝟏) (𝟎)

x´´(t)-2x´(t)-3x(t)=1;x(0)=0,x´(0)=1. Solución
1 1
s²X(s)-sx(0)-x´(0)-2(sX(s)-x(0))-3X(s)= → s²X(s)-1-2(sX(s))-3X(s)=
𝑠 𝑠
1 𝑠+1 1 1 1 1
(s²-2s-3)X(s)= +1→(s-3)(s+1)X(s)= →X(s)= 𝑠(𝑠−3) = ( - )
𝑠 𝑠 3 𝑠−3 𝑠
1
x(t)= 3

Mg. Nolan Jara Jara 39


7. Transformada de Laplace de la integral de una función

𝑆𝑖 𝐿{𝑓(𝑡)} = න 𝑒 −𝑠𝑡 𝑓(𝑡)𝑑𝑡 = 𝐹(𝑠) ⇒


0

𝑡 𝑡
𝐿 න 𝑓(𝜏)𝑑𝜏 = න 𝑒 −𝑠𝑡 න 𝑓(𝜏)𝑑𝜏 𝑑𝑡
0 0
0
∞ ∞
𝑡
1 −𝑠𝑡 1 −𝑠𝑡
1
= න 𝑓(𝜏)𝑑𝜏 − 𝑒 อ + න 𝑒 𝑓(𝑡)𝑑𝑡 = 𝐹(𝑠)
0 𝑠 𝑠 𝑠
0 0
𝒕
𝟏 𝑭(𝒔)
𝑳 න 𝒇(𝝉)𝒅𝝉 = 𝑳{𝒇(𝒕)} =
𝟎 𝒔 𝒔
Mg. Nolan Jara Jara 40
𝑓(𝑡) ∞
,𝐿 = ‫𝑢𝑑 𝑢 𝐹 𝑠׬‬, 𝐹(𝑠) = 𝐿 𝑓(𝑡)
𝑡

Mg. Nolan Jara Jara 41


8. Transformada de Laplace de f(t)/t


𝒇(𝒕)
𝑳 = න 𝑭(𝒖)𝒅𝒖
𝒕 𝒔

con 𝑭(𝒔) = 𝑳 𝒇(𝒕)

Mg. Nolan Jara Jara 42


9. Teorema del valor final

Si lim f (t ) existe, entonces:


t→

lim t → f (t ) = lim s→0 sF ( s)

10. Teorema del valor inicial


El valor inicial f(0) de la función f(t) cuya transformada de Laplace
es F(s), entonces:

f (0) = lim t →0+ f (t ) = lim s→ sF ( s )

Mg. Nolan Jara Jara 43


𝑡 ∞ 𝑡
𝐿 ‫׬‬0 𝑓 𝜏 𝑔 𝑡 − 𝜏 𝑑𝜏 =‫=𝑡׬‬0 𝑒 −𝑠𝑡 ‫=𝜏׬‬0 𝑓 𝜏 𝑔 𝑡 − 𝜏 𝑑𝜏 dt=
∞ 𝑡 ∞ ∞
‫=𝑡׬‬0 ‫=𝜏׬‬0 𝑒 −𝑠𝑡 𝑓 𝜏 𝑔 𝑡 − 𝜏 𝑑𝜏 𝑑𝑡 = ‫=𝝉׬‬0 ‫ 𝑒 𝜏=𝑡׬‬−𝑠𝑡 𝑓 𝜏 𝑔 𝑡 − 𝜏 𝑑𝑡𝑑𝝉
𝜏
R:0≤ 𝒕 < ∞ R:0≤ 𝝉 < ∞
0≤ 𝝉 ≤ 𝒕 𝝉≤𝒕≤∞

𝜏=t

t→ ∞
t=0
R
t
𝜏=0

Mg. Nolan Jara Jara 44


𝑡 ∞ ∞
𝐿 ‫׬‬0 𝑓 𝜏 𝑔 𝑡 − 𝜏 𝑑𝜏 = ‫=𝝉׬‬0 ‫ 𝑒 𝜏=𝑡׬‬−𝑠𝑡 𝑓 𝜏 𝑔 𝑡 − 𝜏 𝑑𝑡𝑑𝝉; 𝑣 =𝑡−𝜏
∞ ∞ ∞ ∞
= ‫=𝜏׬‬0 ‫=𝑣׬‬0 𝑒 −𝑠(𝜏+𝑣) 𝑓 𝜏 𝑔 𝑣 𝑑𝑣𝑑𝜏 = ‫=𝜏׬‬0 ‫=𝑣׬‬0 𝑒 −𝑠𝜏 𝑒 −𝑠𝑣 𝑓 𝜏 𝑔 𝑣 𝑑𝑣𝑑𝜏
∞ ∞
= ‫=𝜏׬‬0 𝑒 𝑓 𝜏 𝑑𝜏 ‫=𝑣׬‬0 𝑒 −𝑠𝑣 𝑔 𝑣 𝑑𝑣 = 𝐹 𝑠 𝐺(𝑠)
−𝑠𝜏

𝒕
𝑳 ‫𝒇 𝟎׬‬ 𝝉 𝒈 𝒕 − 𝝉 𝒅𝝉 =𝑭 𝒔 . 𝑮 𝒔 = 𝑳{𝒇(𝒕) ∗ 𝒈(𝒕)} →
𝒕
‫𝒇 𝟎׬‬ 𝝉 𝒈 𝒕 − 𝝉 𝒅𝝉 =𝑳−𝟏 𝑭 𝒔 . 𝑮(𝒔) =f(t)*g(t)

Mg. Nolan Jara Jara 45


11. Integral de convolución


la integral න 𝑓(𝜏)𝑔(𝑡 − 𝜏)𝑑𝜏
−∞

se conoce como la convolución de 𝑓 𝑡 𝑦 𝑔(𝑡)

y se denota como


න 𝑓(𝜏)𝑔(𝑡 − 𝜏)𝑑𝜏 = 𝑓(𝑡) ∗ 𝑔(𝑡)
−∞

Mg. Nolan Jara Jara 46


Si trabajamos con funciones que son cero para para t < 0,
entonces la convolución queda:

 t f ( ) g (t −  )d , t  0
f (t ) * g (t ) = 0
 0, t0
Así que para estas funciones podemos definir la convolución
como:
𝒕
𝒇(𝒕) ∗ 𝒈(𝒕) = න 𝒇(𝝉)𝒈(𝒕 − 𝝉)𝒅𝝉, (𝒕 ≥ 𝟎)
𝟎

𝒕
𝑳 ‫𝒇 𝟎׬‬ 𝝉 𝒈 𝒕 − 𝝉 𝒅𝝉 =𝑭 𝒔 . 𝑮 𝒔 = 𝑳{𝒇(𝒕) ∗ 𝒈(𝒕)} →
𝒕
‫𝒇 𝟎׬‬ 𝝉 𝒈 𝒕 − 𝝉 𝒅𝝉 =𝑳−𝟏 𝑭 𝒔 . 𝑮(𝒔) =f(t)*g(t)

Mg. Nolan Jara Jara 47


Ejemplo: Verificar que funciona para f(t) = t y g(t) = e-2t
con valores 0 para t < 0.
𝑡 𝑡
𝑓(𝑡) ∗ 𝑔(𝑡) = න 𝑓(𝜏)𝑔(𝑡 − 𝜏)𝑑𝜏 = න 𝜏𝑒 −2(𝑡−𝜏) 𝑑𝜏 =
0 0
𝑡 −2𝑡
𝑡 1 𝑒
𝑒 −2𝑡 න 𝜏𝑒 2𝜏 𝑑𝜏 = − + ⇒
0 2 4 4

 t 1 e −2t  1 1 1 − 2t
L − +  = L{t} − L{1} + L{e } =
2 4 4  2 4 4
1 1 11 1 1 1 1 1
− + = 2 = 2
2s 2
4 s 4 ( s + 2) s ( s + 2) s ( s + 2)
1 − 2t 1
L{t} = 2 ; L{e } =
s ( s + 2)
Mg. Nolan Jara Jara 48
TABLA DE TRANSFORMADAS DE LAPLACE

 (t) 1 sen t

s2 +  2
1
1 cos t
s
s s2 +  2
1 
t e − at sen t
s2 ( ) + 2
+
2
s a
n! s+a
tn + e − at
cos t
sn 1 (s + a) + 2
2

−at 1
e t n e − at
n!
s+a
(s + a)
n +1

Mg. Nolan Jara Jara 49


TABLA DE TRANSFORMADAS DE LAPLACE

1
L{e } =
at
;s  a
s−a
𝑎
L{senℎ(𝑎𝑡)} = 2 2
;𝑠 > 𝑎
𝑠 −𝑎
s
L{cosh(at )} = 2 ;s  a
s −a 2

L{e a t f (t )} = F ( s − a )

1 s
L{ f ( at )} = F ( )
a a

Mg. Nolan Jara Jara 50


TRANSFORMADA INVERSA DE LAPLACE

Al proceso inverso de encontrar f(t) a partir de F(s) se le conoce


como transformada inversa de Laplace y se obtiene mediante:

𝐿 𝑓(𝑡) = 𝐹 𝑠 ↔ 𝑓 𝑡 = 𝐿−1 {𝐹(𝑠)}

conocida también como integral de Bromwich o integral de


Fourier-Mellin.

Mg. Nolan Jara Jara 51


TRANSFORMADAS INVERSAS DE LAPLACE

Algunas transformadas inversas


−11 
(a) 1= L  
s 
−1
1 
(b) −1 n! 
t = L  n +1  , n  Z +
n (c) e =L 
at

s  s − a 
𝑘 −1 s 
(d) sen(𝑘𝑡) = L−1 (e) cos kt = L  2 2
𝑠2 + 𝑘2 s + k 

𝑘 −1s 
(f) senh(kt) = L−1 (g) cosh kt = L  2 2
𝑠2 − 𝑘2 s − k 

Mg. Nolan Jara Jara 52


Hallar las transformadas inversas de

1  −1
1 −1
(a) L  5  (b) L  2 
s  s + 7
Solución
(a) L −1 1  = 1 L −1 4!  = 1 t 4
s 
5
4! s 
5
24

(b) L−1
1
=
1 −1
L
7
=
1
sen( 7 𝑡)
2
𝑠 +7 2
𝑠 +7
7 7

Mg. Nolan Jara Jara 53


Hallar
−1 − 2 s
+ 6
L  2
 s +4 
Solución

−1 − 2 s
+ 6 −2𝑠 6
 2  = L 𝑠2 + 4 + 𝑠2 + 4
−1
L
 s +4  𝑠 2
−1 −1
= −2L + 3L
𝑠2 +4 𝑠2 + 4

= −2 cos( 2 𝑡) + 3 sen( 2 𝑡)

Mg. Nolan Jara Jara 54


Hallar
−1 s 2
+ 6s + 9 
L  
( s − 1)( s − 2)( s + 4) 
Solución
Usando fracciones parciales
s 2 + 6s + 9 A B C
= + +
( s − 1)(s − 2)(s + 4) s −1 s − 2 s + 4

Luego

s + 6s + 9 = A(s − 2)(s + 4) + B( s − 1)(s + 4) + C ( s − 1)(s − 2)


2

Cuando s = 1, 2, −4, entonces

Mg. Nolan Jara Jara 55


16 25 1
A= − , B = ,c =
5 6 30

−1 s 2
+ 6s + 9 
L  
( s − 1)( s − 2)( s + 4) 

16 −1 1  25 −1 1  1 −1 1 


=− L  + L  + L  
5  s − 1 6  s − 2  30 s + 4
16 t 25 2t 1 −4t
=− e + e + e
5 6 30

Mg. Nolan Jara Jara 56


utilizar la convolución para encontrar
transformadas inversas de Laplace:

 1  −1  1 1 
−1
=L  2  = t e
t
L  2
 s ( s − 1)   s s − 1
t
=  e d = e − t − 1
t − t

Mg. Nolan Jara Jara 57


SOLUCIÓN DE EDO LINEALES CON TRANSFORMADA DE
LAPLACE
n n−1
d y d y
an n + an−1 n−1 +  + a0 y = g (t )
dt dt
y(0) = y0 , y(0) = y1 , , y ( n−1) (0) = yn−1
Aplicando la transformada de Laplace

d n y   d n −1 y 
anL  n  + an −1L  n −1  +  + a0L{ y} = L{g (t )} 
 dt   dt 
an [ s nY ( s) − s n−1 y(0) −  − y ( n−1) (0)]
+ an −1[ s n −1Y ( s ) − s n − 2 y (0) −  − y ( n − 2 ) (0)]
+  + a0Y ( s ) = G ( s )
Mg. Nolan Jara Jara 58
SOLUCIÓN DE EDO LINEALES CON TRANSFORMADA DE
LAPLACE

Entonces
P ( s )Y ( s ) = Q ( s ) + G ( s )

Q( s ) G ( s )
Y ( s) = +
P( s) P( s)

𝑷𝟏 ) 𝒚´´ + 𝟐𝒚´ + 𝟐𝒚 = 𝜹 𝒕 − 𝝅 ; 𝒚 𝟎 = 𝟎; 𝒚´ 𝟎 = 𝟎

Solución Mg. Nolan Jara Jara 59


𝑳 𝒚´´ + 𝟐𝒚´ + 𝟐𝒚 = 𝑳 𝜹 𝒕 − 𝝅 ; 𝒚 𝟎 = 𝟎; 𝒚´ 𝟎 = 𝟎
𝑳 𝒚´´ + 2𝑳 𝒚´ + 𝟐𝑳 𝒚 = 𝑳 𝜹 𝒕 − 𝝅
(s²Y(s)-sy(0)-y´(0))+2(SY(s)-y(0))+2Y(s)=𝒆−𝝅𝒔
s²Y(s)+2SY(s) +2Y(s)= 𝒆−𝝅𝒔
𝒆−𝝅𝒔 𝒆−𝝅𝒔 𝒆 −𝝅𝒔
Y(s)=𝒔2+𝟐𝒔+𝟐= (𝒔+𝟏)2 +𝟏 → 𝑳−𝟏 Y(s) = 𝑳−𝟏 (𝒔+𝟏)2+𝟏
−𝟏 −𝝅𝒔 𝟏 −(𝒕−𝝅) sen(t-𝝅).u(t-𝝅)
y(t)= 𝑳 𝒆 = 𝒆
(𝒔+𝟏)2 +𝟏
𝟎 ; 𝒔𝒊 𝟎 < 𝒕 < 𝝅
y(t)=𝒆−(𝒕−𝝅) sen(t-𝝅).u(t-𝝅)=ቊ
𝒆−(𝒕−𝝅) sen(t−𝝅) ; 𝒔𝒊 𝒕 > 𝝅

𝟏
𝑳−𝟏 (𝒔+𝟏)2 +𝟏
=𝒆−𝒕 sent=g(t)
𝐿 𝑢 𝑡 − 𝑎 𝑔 𝑡 − 𝑎 = 𝑒 −𝑎𝑠 𝐺 𝑠 ,
𝑢 𝑡 − 𝑎 𝑔 𝑡 − 𝑎 = 𝐿−1 𝑒 −𝑎𝑠 𝐺 𝑠
Mg. Nolan Jara Jara 60
Mg. Nolan Jara Jara 61
″ sen 𝑡 0≤𝑡<𝜋
𝑃2 ) 𝑦 + 𝑦 = ቊ , 𝑦(0) = 𝑦 ′ (0) = 0
0 𝑡>𝜋

𝑷𝟑 ) 𝒚´´ − 𝟒𝒚´ + 𝟒𝒚 = 𝒆𝟐𝒕 𝒕𝟑 ; 𝒚 𝟎 = 𝟎; 𝒚´ 𝟎 = 𝟎

−𝒕 𝜋
𝑷𝟒 ) 𝒚´´ + 𝟐𝒚´ + 𝟓𝒚 = 𝒆 sent ; 𝒚 𝜋 = 𝟏; 𝒚 =𝟎
𝟒
Solución

sen 𝑡 0≤𝑡<𝜋
𝑃2 ) 𝑦 ″ + 𝑦 = ቊ , 𝑦(0) = 𝑦 ′ (0) = 0
0 𝑡>𝜋

sen 𝑡 , 0≤𝑡<𝜋
f(t)=ቊ
0, 𝑡>𝜋

𝐿 𝑦´´ + 𝑦 =L 𝑓(𝑡) → 𝑠 2 𝑌 𝑠 + 𝑌 𝑠 = L 𝑓(𝑡)


Mg. Nolan Jara Jara 62
∞ −st π
L{f(t)}=‫׬‬0 e f t dt = ‫׬‬0 sente−st dt=uv− ‫ ׬‬vdu
π π
e−st (cost) se−st (sent) e−st (cost+ssent)
L{f(t)}=- s2+1 + s2+1 =-
t=0 s2 +1 t=0
1 e−πs +1 1 e−πs
=- 2 e−πs −1 − 1 = 2 = 2 + 2
s +1 s +1 s +1 s +1
1 e−πs
L{f(t)} = s2+1+ s2+1

1 e−πs 1 e−πs
𝑠2𝑌 𝑠 +𝑌 𝑠 = 2 + 2 →𝑌 𝑠 = 2 +
s +1 s +1 (s +1)² (s2 +1)²
−1 −1 1 −1 e−πs
𝐿 𝑌 𝑠 =𝐿 (s2 +1)²
+𝐿 (s2 +1)²
1 1 e −πs e −πs
y(t)= 𝐿−1 s2+1 . s2 +1 + 𝐿−1 (s2+1)² =sent*sent+ 𝐿−1 (s2 +1)²

Mg. Nolan Jara Jara 63


𝑡
sent*sent=‫׬‬0 𝑠𝑒𝑛𝑢. 𝑠𝑒𝑛 𝑡 − 𝑢 𝑑𝑢

1 𝑡
= න (cos 2𝑢 − 𝑡 − cos 𝑡)𝑑𝑢
2 0
1 𝑠𝑒𝑛 2𝑢−𝑡 𝑡 1 𝑠𝑒𝑛𝑡 −𝑠𝑒𝑛𝑡
=2 2
− 𝑢𝑐𝑜𝑠𝑡 =2 2
− 𝑡𝑐𝑜𝑠𝑡 −( 2 )
𝑢=0
1
= 2 𝑠𝑒𝑛𝑡 − 𝑡𝑐𝑜𝑠𝑡
−1 1 1
𝐿 2 = 𝑠𝑒𝑛𝑡 − 𝑡𝑐𝑜𝑠𝑡
(s +1)² 2
−1 e−πs
𝐿 =
(s2 +1)²
−1 −πs 1 1
𝐿 e . 2 = 𝑠𝑒𝑛(𝑡 − π) − (𝑡 − π)cos(𝑡 − π )u(t- π)
(s +1)² 2

Mg. Nolan Jara Jara 64


e−πs 1
−1
𝐿 2 = 𝑠𝑒𝑛(𝑡 − π) − (𝑡 − π)cos(𝑡 − π) u(t− π)
(s +1)² 2
−1 e−πs 1
𝐿 2 = −𝑠𝑒𝑛𝑡 + (𝑡 − π)cost u(t− π)
(s +1)² 2
1 1
y(t)= 𝑠𝑒𝑛𝑡 − 𝑡𝑐𝑜𝑠𝑡 + −𝑠𝑒𝑛𝑡 + (𝑡 − π)cost u(t− π)
2 2
1
𝑠𝑒𝑛𝑡 − 𝑡𝑐𝑜𝑠𝑡 ; 0 ≤ 𝑡 < 𝜋
2
y(t)=൞1 1
2
𝑠𝑒𝑛𝑡 − 𝑡𝑐𝑜𝑠𝑡 + 2
−𝑠𝑒𝑛𝑡 + (𝑡 − π)cost ; 𝑡 ≥ 𝜋

1
2
𝑠𝑒𝑛𝑡 − 𝑡𝑐𝑜𝑠𝑡 ; 0 ≤ 𝑡 < 𝜋
y(t)=ቐ π
− cost; 𝑡 ≥𝜋
2

Mg. Nolan Jara Jara 65


SOLUCIÓN DE EDO LINEALES CON TRANSFORMADA
DE LAPLACE
1)Resolver 𝑑𝑦
+ 3𝑦 = 13 sen 2 𝑡 , 𝑦(0) = 6
𝑑𝑡
Solución
 dy 
L   + 3L { y} = 13L {sin 2t}
 dt 
26
sY ( s) − 6 + 3Y ( s) = 2
s +4
26
( s + 3)Y ( s) = 6 + 2
s +4
6 26 6s 2 + 50
Y ( s) = + =
s + 3 ( s + 3)( s + 4) ( s + 3)( s 2 + 4)
2

Mg. Nolan Jara Jara 66


6s 2 + 50 A Bs + C
Y (s) = = + 2
( s + 3)(s + 4) s + 3 s + 4
2

6s + 50 = A( s + 4) + ( Bs + C ) 
2 2

A = 8, B = −2, C = 6. Entonces
6s 2 + 50 8 − 2s + 6
Y ( s) = = + 2
( s + 3)( s + 4) s + 3 s + 4
2

−1 1  −1 s  −1 2 


y (t ) = 8L   − 2L  2  + 3L  2 
 s + 3 s + 4 s + 4
𝑦(𝑡) = 8𝑒 −3𝑡 − 2 cos 2𝑡 + 3 sen 2𝑡

Mg. Nolan Jara Jara 67


2) Resolver: −4t
y"−3 y '+2 y = e , y (0) = 1, y ' (0) = 5
Solución
d 2 y   dy  + 2L { y} = L {e −4t }
L  2  − 3L  
 dt   dt 
1
s Y ( s ) − sy (0) − y(0) − 3[ sY ( s ) − y (0)] + 2Y ( s ) =
2
s+4
1
( s − 3s + 2)Y ( s ) = s + 2 +
2
s+4
s+2 1 s 2 + 6s + 9
Y ( s) = 2 + 2 = 
s − 3s + 2 ( s − 3s + 2)(s + 4) ( s − 1)(s − 2)(s + 4)
−1 16 t 25 2t 1 − 4t
y (t ) = L {Y ( s )} = − e + e + e
5 6 30
Mg. Nolan Jara Jara 68
3) RESOLVER : f ' (t ) + 2 f (t ) = e −3t ; (t  0 y f (0) = 4)

f ' (t ) + 2 f (t ) − e −3t = 0 ; L{ f ' (t ) + 2 f (t ) − e −3t } = 0


− 3t
L{ f ' (t )} + 2 L{ f (t )} − L{e } = 0
1
( sF ( s ) − f (0)) + 2 F ( s ) − =0
s+3
1
sF ( s ) − 4 + 2 F ( s ) − =0
s+3
5 1
F (s) = − → f (t ) = 5e − 2t − e −3t
s+2 s+3
Mg. Nolan Jara Jara 69
4) Obtener, mediante la Transformada de Laplace, la
solución del problema de Cauchy

𝑦 ″ + 𝑦 = sen 𝑡 …E
𝑦(0) = 1
𝑦 ′ (0) = 0
Solución

L y ″ = 𝑠 2 𝐿 𝑦 − 𝑠𝑦(0) − 𝑦 ′ (0) = 𝑠 2 𝑌(𝑠) − 𝑠; 𝐿 𝑦 = 𝑌 𝑠 ;


1
𝐿 𝑠𝑒𝑛𝑡 = 1+𝑠2 𝑒𝑛 𝐸:
1
S²Y(S)-S+Y(s)= 2
𝑆 +1

Mg. Nolan Jara Jara 70


Transformada de la ecuación:

″ 2
1
𝐿 𝑦 + 𝑦 = 𝐿 sin 𝑡 ⇒ (𝑠 + 1)𝑌(𝑆) − 𝑠 = 2 ⇒
1+𝑠
𝑠 1
⇒𝑌 𝑠 = 2 + 2 2
⇒𝑦 𝑡
𝑠 +1 𝑠 +1
−1
𝑠 −1
1
=𝐿 2
+𝐿
𝑠 +1 𝑠2 + 1 2

−1
𝑠
𝐿 2
= cos 𝑡
𝑠 +1
1 1 1
𝐿−1 2 2
=𝐿 −1
2 2
= sen 𝑡 ∗ sen 𝑡
𝑠 +1 𝑠 +1𝑠 + 1

Mg. Nolan Jara Jara 71


𝑡
sen 𝑡 ∗ sen 𝑡 = න sen 𝑢 sen( 𝑡 − 𝑢)𝑑𝑢 =
0
1 𝑡 1 𝑡
= න (cos( 2𝑢 − 𝑡) − cos 𝑡)𝑑𝑢 = sen 𝑡 − cos 𝑡
2 0 2 2

1 𝑡
𝑦(𝑡) = cos 𝑡 + sen 𝑡 − cos 𝑡
2 2

Mg. Nolan Jara Jara 72


5) En la ED di
L + Ri = E (t )
dt
Hallar i(t) donde i(0) = 0, E(t) es Como ilustra la Figura

Mg. Nolan Jara Jara 73


Solución

1
LsI ( s ) + RI ( s ) =
s (1 + e − s )
1/ L 1
I (s) = 
s ( s + R / L) 1 + e − s

1 −s −2 s −3 s
= 1 − e + e − e +
1+ e -s

1 L/ R L/ R
= −
s ( s + R / L) s s + R/ L

Mg. Nolan Jara Jara 74


1 1
I ( s) =  −
1 
(
 1 − e−s + e−2 s − ...
R s s + R L
)

Luego i(t) se escribe de la siguiente manera y se ilustra en la Fig *1

1 − 𝑒 −𝑡 , 0≤𝑡<1
−𝑒 −𝑡 + 𝑒 −(𝑡−1) , 1≤𝑡<2
𝑖(𝑡) =
1 − 𝑒 −𝑡 + 𝑒 −(𝑡−1) − 𝑒 −(𝑡−2) , 2≤𝑡<3
−𝑒 −𝑡 + 𝑒 −(𝑡−1) − 𝑒 −(𝑡−2) + 𝑒 −(𝑡−3) , 3≤𝑡<4

Mg. Nolan Jara Jara Ch4_75


Fig *1

Mg. Nolan Jara Jara 76


sin t 0  t 
6)Resolver y + y =  , y(0) = y(0) = 0
 0 t 
s 2Y ( s ) + Y ( s) = L sin t − u (t −  ) sin t
= L sin t + u (t −  ) sin(t −  )
−s
1 e
= 2 + 2
s +1 s +1
−s
1 e
Y (s) = 2 + 2
( s + 1) ( s + 1)
2 2

y(t ) = 12 sin t − t cost  + u (t −  ) 12 sin(t −  ) − (t −  ) cos(t −  )


 12 sin t − t cost  0  t  
=
− 2  cos t t 
1

Mg. Nolan Jara Jara 77
7) Resolver
y  + 3 y  + 2 y =  (t − 1), y (0) = y (0) = 0

s 2Y ( s) + 3sY ( s) + 2Y ( s) = e − s
−s 1 −s  1 1 
Y ( s) = e 2 =e  − 
s + 3s + 2  s + 1 s + 2 


y (t ) = u (t − 1) e − ( t −1) − e −2( t −1) 

Mg. Nolan Jara Jara 78


ty + y + 4ty = 0, y(0) = 0, y (0) = 3.Y = cj0 (2t )
sent
ty + 2 y + ty = 0, y (0) = 0, y ( ) = 0.Y =
t
y + ty − y = 0, y(0) = 1, y (0) = 0.Y = t
ty + (1 − 2t ) y − 2 y = 0, y(0) = 2, y (0) = 1.Y = e 2t
ty + (t − 1) y − y = 0, y () = 0, y (0) = 5.Y = 5e −t

Mg. Nolan Jara Jara 79


Función Gamma

 :  0,  ) → R

t
x −1 − t
( x ) = e dt
t =0
Esta integral impropia es convergente

L f (t ) =  f (t )e − st dt = F ( s )
0

Mg. Nolan Jara Jara 80


Propiedades de la Función Gamma

 (1) = 1
( x + 1) = x ( x )
 ( n + 1) = n !
1
( ) = 
2

1 − 1 −t
( ) =  t 2 e dt = 
2 0

Mg. Nolan Jara Jara 81


t
−t
 ( p + 1) = p
e dt
0

 ( p + 1)
 
L t p
=
s p +1
; p  −1, s  0

 = s
L t
−1
2

Lt =
1 2
3
2
2s

Mg. Nolan Jara Jara 82


SERIES

Si trata de sumar los términos de una sucesión infinita a 


n n =1
a1 + a2 + a3 + ... + an + ... ,
tal expresión se denomina serie infinita, o sólo serie, y se denota
con el símbolo

a
n =1
n
o a n

Mg Nolan Jara Jara 1


DEFINICIÓN


Dada la serie: n =1
an = a1 + a2 + a3 + ...

Denotemos con Sn la n-ésima suma parcial:


n
S n =  ak = a1 + a2 + a3 + ... + an
k =1

Si la sucesión S n  es convergente y nlím


→
Sn =S
existe, la serie a n se denomina convergente y se
escribe a n = S.

El número S se llama suma de la serie.


Si no fuese así, la serie se llama divergente.
Mg Nolan Jara Jara 2
EJEMPLOS
Diga si las series convergen o divergen.
Si alguna converge, calcule su suma.

1.  n
n =1


2.  (−1) n
n =1


1
3. 
n =1 n( n + 1)

Mg Nolan Jara Jara 3


SERIE GEOMÉTRICA
Un ejemplo importante de una serie infinita es la serie
geométrica:

Mg Nolan Jara Jara 4


SUMA DE UNA SERIE
GEOMÉTRICA
Suma de la serie geométrica:

 ar
n −1 2
= a + ar + ar + ...
n =1
• Si | r |  1, la serie geométrica converge y la suma es
igual a:

S =  ar
n −1
=
a
n =1 1− r
• Si | r |  1 y a  0, la serie diverge.

Mg Nolan Jara Jara 5


EJEMPLOS

1. Calcule la suma se la serie geométrica:


10 20 40
5− + − + ...
3 9 27

2. ¿Es convergente o divergente la serie dada?


2
n =1
2n
3 1− n

Mg Nolan Jara Jara 6


TEOREMA

Si la serie a
n =1
n es convergente, entonces:

lím an = 0
n→

¡Pero el enunciado inverso no es cierto!

Mg Nolan Jara Jara 7


PRUEBA DE LA DIVERGENCIA
Si lím an no existe o lím an  0 , entonces
n → n→

a
n =1
n es divergente.

Mg Nolan Jara Jara 8


EJEMPLO

Diga si la serie es convergente o divergente.



n2

n =1 5n 2
+4


n2

n =1 5n 2
+4
divergente

Mg Nolan Jara Jara 9


TEOREMA

Si a n y b n son series convergentes, entonces


también lo son las series c a n (donde c es constante),
y  (.a n + bn ).

    
(i)  ca
n =1
n = c  an
n =1
(ii)  (a
n =1
n + bn ) =  an +  bn
n =1 n =1

Mg Nolan Jara Jara 10


EJEMPLO

Determine la suma de la serie:



 3 1 
  + n 
n =1  n( n + 1) 2 

Mg Nolan Jara Jara 11


SERIES ALTERNANTES

Una seriede la forma :


 n − 1b = b − b + b − b + ...
 ( −1) n 1 2 3 4
n =1
o también
 n
 ( −1) bn = −b1 + b2 − b3 + b4 − ...
n =1
se llama seriealternante.

Mg Nolan Jara Jara 12


PRUEBA DE LA SERIE
ALTERNANTE
Si la serie alternante

 (− 1)n −1
bn = b1 − b2 + b3 − b4 + b5 − b6 + ... (bn  0)
n =1

cumple con:

bn +1
 bn

(ii) lím b n
=0
n →

entonces la serie es convergente

Mg Nolan Jara Jara 13


EJEMPLOS
Determine si las series convergen o divergen:
1 1 1 1 1
1. 1 − + − + − + ... = 

(− 1)
n −1

Convergente
2 3 4 5 6 n =1 n


(− 1)n 3n
2. 
n =1 4n − 1
Divergente


n2
 (− 1)
n +1
3. Convergente
n =1 n3 + 1

Mg Nolan Jara Jara 14


CONVERGENCIA ABSOLUTA

Una serie a
n =1
n es denominada absolutamente convergente

si la serie de los valores absolutos a
n =1
n
es convergente.

Una serie a
n =1
n se llama condicionalmente convergente

si es convergente, pero no absolutamente convergente.

Mg Nolan Jara Jara 15


EJEMPLOS
Determine si las series son condicionalmente
convergentes, absolutamente convergentes
o divergentes.
 
1
Nota: 
1
divergey 
n =1 n
2
converge
n =1 n


(− 1)n −1 1 1 1 1
1. 
n =1 n2
= 1 − 2 + 2 − 2 + 2 − ... Absolutamente
2 3 4 5 convergente


(− 1)n −1 1 1 1 1
2. 
n =1 n
= 1 − + − + − ...Condicionalmente
2 3 4 5 convergente
Mg Nolan Jara Jara 16
TEOREMA

Si una serie es absolutamente


convergente, entonces es convergente.

Mg Nolan Jara Jara 17


CRITERIO DE LA RAZÓN

an +1
Sea lím =L
n → a
n


- Si L  1, entonces la serie  an es convergente.
n =1

- Si L  1, o L = , entonces la serie  an es divergente.
n =1

- Si L = 1, el criterio no da información.

Mg Nolan Jara Jara 18


EJEMPLOS

n3
1. Diga si la serie  (− 1)
n =1
n

3n
es
absolutamente convergente.


nn
2. Diga si la serie 
n = 1 n!
es convergente.

Mg Nolan Jara Jara 19


SERIES DE POTENCIAS
Una serie de potencias tiene la forma:

 c (x − a ) = c0 + c1 (x − a ) + c2 (x − a ) + ...
k 2
k
k =0

Potencias en: x−a


con centro en: a
Nota: Se toma la convención de que (x – a)0 = 1
incluso si x = a.
Mg Nolan Jara Jara 20
RADIO DE CONVERGENCIA
Teorema 

c ( x − a)
k
Para una serie de potencias k
k =0

Se cumple una y sola una de las


siguiente posibilidades
i)La serie converge sólo cuando x=a
ii)La serie converge para todo x
iii) Existe un número R>0 tal que la serie
converge si x − a  R y diverge si x − a  R
Mg Nolan Jara Jara 21
RADIO DE CONVERGENCIA
(Continua)
convergencia para x−a  R
R R

a−R a a+R

divergencia para
x−a  R
Mg Nolan Jara Jara 22
RADIO DE CONVERGENCIA
(Continua)
El número R>0 en el caso (iii) es llamado
radio de convergencia de la serie de
potencias

c ( x − a)
k
k
k =0

Por convención, el radio de convergencia

es R=0 en el caso (i) y R= ∞ en el caso (ii)


Mg Nolan Jara Jara 23
Convergencia absoluta
Una serie de potencias converge absolutamente
Si, la siguiente serie converge:

n=0 | cn ( x − a)
 n
|
Prueba de convergencia (criterio del cociente)
Si cn  0 para todo n, y
cn +1 ( x − a) n +1 cn +1
lim = | x − a | lim =L
n → cn ( x − a ) n n → c
n

Si L < 1, la serie converge absolutamente


si L > 1, la serie diverge
y si L = 1, el criterio no es concluyente.
•Ch5_24
an
Si  R
r = lim  rRes el radio de convergencia de la serie
n → a
n +1

Si el limite es infinito  el radio de convergencia es infinito

1
Si rR== lim n  rRes el radio de convergencia de la serie
n→ an
si el limite es infinito  el radio de convergencia es infinito
 n
x
La serie de potencias  tiene radio de convergencia infinito
n = 0 n!

( x − 1)
 n

La serie de potencias  (− 1)
n +1
tiene radio de convergencia
n =0 n
Rr = 1,0  x  2
 n

La serie de potencias  cn ( x − a ) es absolutamente convergente


n =0

cn +1 ( x − a )
n +1

si L  1, L= n lim
cn ( x − a )
n

Mg Nolan Jara Jara 25


RADIO DE CONVERGENCIA

Ejemplo
Determine el radio de convergencia
de las siguientes series de potencias
 
1 k
1)  k ! x k
2)  x
k =0 k =0 k!
( −3)
n
 n
x
3) 
n =0 n +1
Mg Nolan Jara Jara 26
DERIVACIÓN E INTEGRACIÓN
DE SERIES

Si una serie de potencias f ( x) =  c (x − a )
k
k
k =0
es convergente en un intervalo abierto I,
entonces es derivable en I y se cumple:

f ( x) =  kck (x − a ) ;
k −1
1. x  I.
k =1

ck ( x − a )
 k +1

2.  f ( x)dx = C +  ; x  I.
k =0 k +1
Mg Nolan Jara Jara 27
EJEMPLOS

Se sabe que:

1
= 1 + x + x + x + ... =  x , − 1  x  1
2 3 k

1− x k =0

En base a lo anterior, deduzca una


representación en serie de potencias de
las siguientes funciones:
1 f ( x) = arctan x
f ( x) =
(1 − x )2
Mg Nolan Jara Jara 28
SERIE DE TAYLOR

Si f tiene derivadas de todo orden en


a  I, entonces:

f ( x) = 
f (k )
(a ) (x − a )k
k =0 k!
se llama la serie de Taylor de f
alrededor de a.
Mg Nolan Jara Jara 29
OBSERVACIONES

1. La serie de Taylor de alrededor de a.


converge a f en un punto x  I

 lím Rn ( x) = 0
n →
Rn (x) es el residuo de orden n.

2. Si en la Serie de Taylor a = 0,

Mg Nolan Jara Jara 30


ALGUNAS SERIES DE TAYLOR
 k
x
e =
x
; x  
k = 0 k!

 2 k −1
sen x =  (− 1)
x
 x 
k
;
k =0 (2k − 1)!
¿Qué serie le corresponderá a cos x = ?

Mg Nolan Jara Jara 31


EJEMPLOS
1) Utilice una serie conocida y
represente f como serie de f ( x) = ln(1 − x)
potencias:

2) Desarrolle f en serie de
Taylor alrededor de  / 3 :
f ( x) = cos x

Mg Nolan Jara Jara 32


ANALITICIDAD DE FUNCIONES

Definición:

Una función f es analítica en un punto a


si se puede expresar mediante una
serie de potencias alrededor de a.

Mg Nolan Jara Jara 33


PUNTOS ORDINARIOS Y SINGULARES
DE UNA EDOL

Sea la EDOL: a2 ( x) y ' '+ a1 ( x) y '+ a0 ( x) y = 0


Se puede escribir
en la forma: y ' '+ P ( x) y '+Q( x) y = 0

Definición:
Se dice que x0 es un punto ordinario de la
ecuación dada si P(x) y Q(x) son analíticas
en x0. Un punto que no es ordinario se llama
punto singular de la ecuación.
Mg Nolan Jara Jara 34
TEOREMA

Sea la EDOL: a2 ( x) y ' '+ a1 ( x) y '+ a0 ( x) y = 0


Si x0 es un punto ordinario de la EDOL,
existen dos soluciones L.I. en forma de series
de potencias centradas en x0, es decir, dos
soluciones de la forma:

y =  cn (x − x0 )
n

n =0

Mg Nolan Jara Jara 35


CORRIMIENTO DE ÍNDICES

Escriba la expresión siguiente


como una sola serie de la forma  an x :
n

 

 n(n − 1)c x
n=2
n
n−2
+  cn x
n =0
n +1

Mg Nolan Jara Jara 36


CORRIMIENTO DE ÍNDICES

Escriba la expresión siguiente


como una sola serie de la forma  an x :
n

(x + 1) n(n − 1)cn x


  
2 n−2
+ x ncn x n −1
− cn x n

n=2 n =1 n =0

Mg Nolan Jara Jara 37


EJEMPLO

Encuentre dos soluciones L.I. en forma


de series de potencia alrededor de 0 de
la EDOL:

y ' '+ xy = 0

Mg Nolan Jara Jara 38


EJEMPLO

Encuentre dos soluciones L.I. en forma


de series de potencia alrededor de 0 de
la EDOL:
( x + 1) y ' '+ xy'− y = 0
2

Mg Nolan Jara Jara 39


SOLUCION DE ECUACIONES
DIFERENCIALES MEDIANTE
SERIES DE POTENCIAS

8/10/2022 Mg. Nolan Jara Jara 1


ANALITICIDAD DE FUNCIONES

Definición: Una función f es analítica en un


punto a si se puede expresar mediante una serie
de potencias alrededor de a.

8/10/2022 Mg. Nolan Jara Jara 2


PUNTOS ORDINARIOS Y SINGULARES
DE UNA EDOL
Dada la EDOL: a2 ( x) y ' '+ a1 ( x) y '+ a0 ( x) y = 0
Que se puede escribir en la forma:
y ' '+ P ( x) y '+Q( x) y = 0
Definición:
Se dice que x0 es un punto ordinario de la
ecuación dada si P(x) y Q(x) son analíticas en x0.
Un punto que no es ordinario se llama punto
singular de la ecuación.
8/10/2022 Mg. Nolan Jara Jara 3
TEOREMA
Sea la EDOL: a2 ( x) y ' '+ a1 ( x) y '+ a0 ( x) y = 0

Si x0 es un punto ordinario de la EDOL, existen dos


soluciones L.I. en forma de series de potencias centradas
en x0, es decir, dos soluciones de la forma:


y =  cn (x − x0 )
n

n =0

8/10/2022 Mg. Nolan Jara Jara 4


CORRIMIENTO DE ÍNDICES
Escriba la expresión siguiente
 

 n(n − 1)c
n=2
n x n−2
+  cn x
n =0
n +1

como una sola serie de la forma

a x n
n
:

8/10/2022 Mg. Nolan Jara Jara 5


 

 n(n − 1)c
n=2
n x n−2
+  cn x
n =0
n +1

∞ ∞

෍ 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛−2 + ෍ 𝑐𝑛 𝑥 𝑛+1 =
𝑛=2 ∞ 𝑛=0 ∞

෍ (𝑁 + 2) 𝑁 + 1 𝑐𝑁+2 𝑥 𝑁 + ෍ 𝑐𝑀−1 𝑥 𝑀 =
𝑁=0 𝑀=1
∞ ∞

2𝑐2 + ෍ (𝑁 + 2) 𝑁 + 1 𝑐𝑁+2 𝑥 𝑁 + ෍ 𝑐𝑀−1 𝑥 𝑀 =


𝑁=1 𝑀=1
∞ ∞

2𝑐2 + ෍ (𝑛 + 2) 𝑛 + 1 𝑐𝑛+2 𝑥 𝑛 + ෍ 𝑐𝑛−1 𝑥 𝑛 =


𝑛=1 𝑛=1

2𝑐2 + ෍ (𝑛 + 2) 𝑛 + 1 𝑐𝑛+2 +𝑐𝑛−1 𝑥 𝑛


𝑛=1
8/10/2022 Mg. Nolan Jara Jara 6
CORRIMIENTO DE ÍNDICES
Escriba la expresión siguiente
∞ ∞ ∞

𝑥 2 + 1 ෍ 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛−2 + 𝑥 ෍ 𝑛𝑐𝑛 𝑥 𝑛−1 − ෍ 𝑐𝑛 𝑥 𝑛


𝑛=2 𝑛=1 𝑛=0

𝑥 2 + 1 𝑦´´ + 𝑥𝑦´ − 𝑦 ∞

෍ 𝑎𝑛 𝑥 𝑛 . 𝑆𝑜𝑙𝑢𝑐𝑖𝑜𝑛
como una sola serie de la forma 𝑛=0

(x²+1)σ∞
𝑛=2 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛−2 + 𝑥 σ∞ 𝑛𝑐 𝑥 𝑛−1 − σ∞ 𝑐 𝑥 𝑛 =
𝑛=1 𝑛 𝑛=0 𝑛

x²σ∞
𝑛=2 𝑛 𝑛 − 1 𝑐𝑛 𝑥
𝑛−2
+ σ∞
𝑛=2 𝑛 𝑛 − 1 𝑐𝑛 𝑥
𝑛−2
+ 𝑥 σ∞
𝑛=1 𝑛𝑐𝑛 𝑥
𝑛−1
− σ∞ 𝑛
𝑛=0 𝑐𝑛 𝑥 =

∞ ∞ ∞ ∞

෍ 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛 + ෍ 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛−2 + ෍ 𝑛𝑐𝑛 𝑥 𝑛 − ෍ 𝑐𝑛 𝑥 𝑛 =
𝑛=2 𝑛=2 𝑛=1 𝑛=0
8/10/2022 Mg. Nolan Jara Jara 7
∞ ∞ ∞ ∞

෍ 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛 + ෍ (𝑁 + 2) 𝑁 + 1 𝑐𝑁+2 𝑥 𝑁 + ෍ 𝑛𝑐𝑛 𝑥 𝑛 − ෍ 𝑐𝑛 𝑥 𝑛 =
𝑛=2 𝑁=0 𝑛=1 𝑛=0

(2c2 + 6c3 x)+c1 x − c0 + c1 x +


∞ ∞ ∞ ∞

෍ 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛 + ෍ (𝑛 + 2) 𝑛 + 1 𝑐𝑛+2 𝑥 𝑛 + ෍ 𝑛𝑐𝑛 𝑥 𝑛 − ෍ 𝑐𝑛 𝑥 𝑛 =
𝑛=2 𝑛=2 𝑛=2 𝑛=2

2c2 + 6c3 x − c0 + σ∞ 𝑛
𝑛=2 𝑛 𝑛 − 1 𝑐𝑛 + 𝑛 + 2 𝑛 + 1 𝑐𝑛+2 + 𝑛𝑐𝑛 − 𝑐𝑛 𝑥 =

2c2 + 6c3 x − c0 + ෍ 𝑛² − 1 𝑐𝑛 + 𝑛 + 2 𝑛 + 1 𝑐𝑛+2 𝑥 𝑛


𝑛=2

8/10/2022 Mg. Nolan Jara Jara 8


EJEMPLO
1) Encuentre la solución en forma de series de potencia alrededor de
0 de la EDOL: 𝑦´ − 2𝑥𝑦 = 0

2) Encuentre dos soluciones L.I. en forma de series de


potencia alrededor de 0 de la EDOL:
𝑦′′ + 𝑥𝑦 = 0
3) Encuentre dos soluciones L.I. en forma de series de
potencia alrededor de 0 de la EDOL:
(𝑥 2 + 1)𝑦′′ + 𝑥𝑦′ − 𝑦 = 0. Solución

𝑦 = σ∞
𝑛=0 𝑐𝑛 𝑥 𝑛
; 𝑦´ = σ∞
𝑛=1 𝑛𝑐𝑛 𝑥 𝑛−1
; 𝑦´´ = ෍ 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛−2

8/10/2022 Mg. Nolan Jara Jara 𝑛=2 9


(x²+1)σ∞
𝑛=2 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛−2
+ 𝑥 σ ∞
𝑛=1 𝑛𝑐𝑛 𝑥 𝑛−1
− σ ∞
𝑛=0 𝑐𝑛 𝑥 𝑛
=0

x²σ∞
𝑛=2 𝑛 𝑛 − 1 𝑐𝑛 𝑥
𝑛−2 + σ∞ 𝑛 𝑛 − 1 𝑐 𝑥 𝑛−2 + 𝑥 σ∞ 𝑛𝑐 𝑥 𝑛−1 − σ∞ 𝑐 𝑥 𝑛 = 0
𝑛=2 𝑛 𝑛=1 𝑛 𝑛=0 𝑛

∞ ∞ ∞ ∞

෍ 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛 + ෍ 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛−2 + ෍ 𝑛𝑐𝑛 𝑥 𝑛 − ෍ 𝑐𝑛 𝑥 𝑛 = 0
𝑛=2 𝑛=2 𝑛=1 𝑛=0

∞ ∞ ∞ ∞

෍ 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛 + ෍ (𝑁 + 2) 𝑁 + 1 𝑐𝑁+2 𝑥 𝑁 + ෍ 𝑛𝑐𝑛 𝑥 𝑛 − ෍ 𝑐𝑛 𝑥 𝑛 = 0
𝑛=2 𝑁=0 𝑛=1 𝑛=0

(2c2 + 6c3 x)+c1 x − c0 + c1 x +


∞ ∞ ∞ ∞

෍ 𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛 + ෍ (𝑛 + 2) 𝑛 + 1 𝑐𝑛+2 𝑥 𝑛 + ෍ 𝑛𝑐𝑛 𝑥 𝑛 − ෍ 𝑐𝑛 𝑥 𝑛 = 0
𝑛=2 𝑛=2 𝑛=2 𝑛=2

𝑐0 = 2𝑐2 , 𝑐3 = 0; ෍ 𝑛 𝑛 − 1 𝑐𝑛 + 𝑛 + 2 𝑛 + 1 𝑐𝑛+2 + 𝑛𝑐𝑛 − 𝑐𝑛 𝑥 𝑛 = 0


𝑛=2
8/10/2022 Mg. Nolan Jara Jara 10

𝑐0 = 2𝑐2 , 𝑐3 = 0; ෍ 𝑛² − 1 𝑐𝑛 + 𝑛 + 2 𝑛 + 1 𝑐𝑛+2 𝑥 𝑛 = 0
𝑛=2

1 1−𝑛
𝑐2 = 2 𝑐0 , 𝑐3 = 0; 𝑐𝑛+2 =𝑛+2 𝑐𝑛 ; 𝑛 = 2,3, …

−1 −1 1 −1 −1 3
𝑛= 2; 𝑐4 = 4 𝑐2 = 4 2 𝑐0 = 8 𝑐0 = 2 𝑐0
−2
𝑛= 3; 𝑐5 = 5 𝑐3 = 0
−3 −1 −1 1 −1 4
𝑛= 4; 𝑐6 = 6 𝑐4 = 2 ( 8 𝑐0 )= 16 𝑐0 = 2 𝑐0
−4
𝑛= 5; 𝑐7 = 7 𝑐5 = 0; 𝑛 = 7; 𝑐9 = 0; …
−5 −5 −1 4 −1 7
𝑛= 6; 𝑐8 = 𝑐6 = 𝑐0 = 5 𝑐0
8 8 2 2
−7 −7 −1 7 −1 8
𝑛= 8; 𝑐10 = 𝑐8 = 5 𝑐0 = 7 𝑐0
10 10 2 2

8/10/2022 Mg. Nolan Jara Jara 11


𝑦 = σ∞ 𝑐
𝑛=0 𝑛 𝑥 𝑛 = 𝑐 + 𝑐 𝑥 + 𝑐 𝑥 2 + 𝑐 𝑥 3 + 𝑐 𝑥 4+ 𝑐 𝑥 5 + 𝑐 𝑥 6 + ⋯
0 1 2 3 4 5 6
3 4 7
1 2
−1 −1 −1
𝑦 = 𝑐0 + 𝑐1 𝑥 + 𝑐0 𝑥 + 𝑐0 𝑥 4+ 𝑐0 𝑥6 +5 𝑐0 𝑥 8 + ⋯
2 2 2 2

3 4 7 8
1 −1 −1 −1 −1
𝑦 = 𝑐0 1 + 𝑥2 + 𝑥4 + 𝑥6 + 5 𝑥8 + 7 𝑥 10 + ⋯ + 𝑐1 𝑥
2 2 2 2 2

3 4 7 8
1 2 −1 4
−1 6
−1 8
−1
𝑦 = 𝑐0 1+ 𝑥 + 𝑥 + 𝑥 +5 𝑥 +7 𝑥 10 + ⋯ + 𝑐1 𝑥
2 2 2 2 2

2 ∞
𝑥 𝑛−1
1.3.5. ⋯ 2𝑛 − 3 2𝑛
𝑦 = 𝑐0 1 + + ෍ −1 𝑛 𝑥 + 𝑐1 𝑥; ; 𝑥 < 1
2 2 𝑛!
𝑛=2

8/10/2022 Mg. Nolan Jara Jara 12


SOLUCIONES RESPECTO A PUNTOS SINGULARES
Un punto singular x0 de una EDO lineal
a2 ( x) y + a1 ( x) y + a0 ( x) y = 0
puede ser regular o irregular. La clasificación
depende de y + P( x) y + Q( x) y = 0

PUNTOS SINGULARES REGULARES E IRREGULARES


Se dice que un punto singular x0 es un punto singular
regular si p(x) = (x – x0)P(x) y q(x) = (x – x0)2Q(x)
son analíticas en x0 ,es decir admiten desarrollos en series
de potencias centradas en x0. Un punto singular que no
es regular es un punto singular irregular de la ecuación.

8/10/2022 Mg. Nolan Jara Jara 13


COEFICIENTES POLINOMIALES
• Si (x – x0) aparece a lo sumo a la primera potencia en el
denominador de P(x) y a lo sumo a la segunda potencia en el
denominador de Q(x), entonces x0 es un punto singular
regular.
• Para resolver la EDO la multiplicaremos en forma estándar
por (x – x0)2:

( x − x0 ) 2 y + ( x − x0 ) 2 P( x) y + ( x − x0 ) 2 Q( x) y = 0
( x − x0 ) y + ( x − x0 ) p( x) y + q( x) y = 0
2

donde p(x) = (x – x0)P(x) y q(x) = (x – x0)2Q(x) son analíticas en x = x0. Observe


que de esta manera hemos "matado" las singularidades.
8/10/2022 Mg. Nolan Jara Jara 14
• Ejemplo: observemos que x = 2, x = – 2 son puntos
singulares de:
(x2 – 4)2y” + 3(x – 2)y’ + 5y = 0
Entonces:
3 5
P( x) = Q( x) =
( x − 2)( x + 2) 2
( x − 2) 2 ( x + 2) 2

Para x = 2, la potencia de (x – 2) en el denominador de P(x)


es 1, y la potencia de (x – 2) en el denominador de Q(x) es 2.
Así que x = 2 es un punto singular regular.

Para x = −2, la potencia de (x + 2) en el denominador de P(x)


y Q(x) es 2. Así que x = − 2 es un punto singular irregular.

8/10/2022 Mg. Nolan Jara Jara 15


TEOREMA DE FROBENIUS

Si x = x0 es un punto singular regular, entonces


existe al menos una solución de la forma
 
y = ( x − x0 ) r  cn ( x − x0 ) n =  cn ( x − x0 ) n+ r
n =0 n =0

donde el número r es una constante por determinar.


La serie converge al menos en algún intervalo
0 < x – x0 < R.

8/10/2022 Mg. Nolan Jara Jara 16


Encuentre la solución de las siguientes ecuaciones
diferenciales mediante la utilización de series de potencias:

y´ - y = 0
Solución
+
Sea y =  n la solución general de la ED.
C
n=0
x n

+

Derivándola: y =  nCn x
 n −1

n = 01
Sustituyendo en la ecuación diferencial: y´− y = 0
+ +

 n
nC x
n = 01
n −1

n=0
 n =0
C x n

Para poder sumar las series, requerimos efectuar cambio de variable:


8/10/2022 Mg. Nolan Jara Jara 17
+ +
 +
n −1   + n
 nC x
n = 01
n
n −1
−  Cn x
n =0
n
= 0 +  nCn x  −  Cn x 
 n =1   n =0 
+ +
=  (k + 1)Ck +1 x −  Cn x n
k

k =0 n =0
+ +
=  (k + 1)Ck +1 x k −  Cn x n
k =0 n =0

  
+
=  (k + 1)Ck +1 x k − C K x k
k =0

Como:

Ck
x  0 entonces ( k + 1)Ck +1 − Ck = 0
k
y Ck +1 =
k +1
8/10/2022 Mg. Nolan Jara Jara 18
Para C0
k = 0 → C1 = = C0
1
C1 C0
k = 1 → C2 = =
2 2
C2 C1 C0
k = 2 → C3 = = =
3 3( 2) 3( 2)
C3 C2 C1 C0
k = 3 → C4 = = = =
4 4(3) 4(3)( 2) 4(3)( 2)
Con estos coeficientes tenemos:

8/10/2022 Mg. Nolan Jara Jara 19


+
y =  Cn x n = C0 + C1 x + C2 x 2 + C3 x 3 + C4 x 4 + 
n =0

C0 2 C0 3 C0 4
= C0 + C0 x + x + x + x +
2 6 24
x 2 x3 x 4
= C0 (1 + x + + + + )
2! 3! 4!
= C0 (e x )

Que representa la solución de la ecuación diferencial:


y´ - y = 0

8/10/2022 Mg. Nolan Jara Jara 20


xy = y + 1
( x − 1) y + ( 2 x + 1) y = 0
xy − y = 0
y + y = 0
x( x 3 − 1) y − ( 4 x 3 − 1) y = 0
y(1 + x ) = 1

8/10/2022 Mg. Nolan Jara Jara 21


xy  = y + 1
+ +
x  nC n x n −1
=  Cn x n + 1
n =1 n =0
+ +

 nC
n =1
n x =  Cn x n + 1
n

n =0
+ +

 nC
n =1
n x = C0 +  Cn x n + 1
n

n =1
+

 (nC
n =1
n − Cn ) x n = C0 + 1
+

 (n. − 1)C
n =1
n x n = C0 + 1

C 0 + 1 = 0 → C 0 = −1
( n − 1)Cn = 0 pero C n  0 ya que pertenece a la solución → n − 1 = 0  n = 1
Por ello, sólo queda C1 x
La solición es : y = C 0 + C1 x
y = −1 + C1 x
8/10/2022 Mg. Nolan Jara Jara 22
( x − 1) y  + ( 2 x + 1) y = 0
+ + + +
x  nC n x n −1
−  nC n x n −1
+ 2 x C n x +  C n x n = 0
n

n =1 n =1 n =0 n =0
+ + + +

 nC
n =1
n x −  ( k + 1)C k +1 x + 2 C n x
n

k =0
k

n =0
n +1
+  Cn x n = 0
n =0
+ + + +

 nC
n =1
n x −  ( k + 1)C k +1 x + 2 C k −1 x +  C n x n = 0
n

k =0
k

k =1
k

n =0
+
 +
k 
+
 +
n

n =1
nC n x −  1 

Cn
+
k =1
( k + 1)C k +1 x 

+ 2 
n =1
C k −1 x k
+  0  n =0

C +
n =1
C x

+ +

 ( kC
k =1
k − ( k + 1)C k +1 + 2C k −1 + C n ) x = (C1 − C 0 ) +  0 x k
n

k =1

C1 − C 0
2C k −1 + C k + kC k
C k +1 = para k = 1,2,3, 
´k + 1
8 11 71
C 2 = 2C 0 , C 3 = C 0 , C 4 = C0 , C5 = C 0 , etc.
3 3 15
8 11 4 71 5
La solución es : y = C 0 (1 + x + 2 x 2 + x 3 + x + x + )
3 3 15
8/10/2022 Mg. Nolan Jara Jara 23
Re solver y − xy − y = 0 en un entorno de x = 0

y(0) = 1; y (0) = 0.sol. y = e


x2
2

Resolver: (1+x²)y´´+2xy´-2y=0 en un entorno de x=0;


𝑦1 (0)=1, 𝑦1 ´ 0 = 0; 𝑦2 (0)=0, 𝑦2 ´ 0 = 1

8/10/2022 Mg. Nolan Jara Jara 24


Resolver las siguientes ecuaciones diferenciales en un entorno
del punto x=0, mediante el uso de series de potencias.

1)(1 + x ) y + 2 xy − 2 y = 0; y1 (0) = 1, y1 (0) = 0


2

y2 (0) = 0, y2 (0) = 1


Sol. y ( x) = Ax + B(1 + xarctgx )
2) y + xy = 0
3)(1 + x ) y + xy − y = 0; y2 = x
2

n −1 1.3.5. (2n − 3) 2 n
2 
+  (− 1)
x
y1 = 1 + n
x ; x 1
2 n=2 2 n!
4) y − (1 + x) y = 0
8/10/2022 Mg. Nolan Jara Jara 25
5) y  − xy = 0
6) y  + x y = 0
2

7) y  − 2 xy  + y = 0
8) y  − xy  + 2 y = 0
9) y  + x y  + xy = 0
2

10) y  + 2 xy  + 2 y = 0
8/10/2022 Mg. Nolan Jara Jara 26
EJEMPLO 1 DEL MÉTODO DE FROBENIUS
1
𝑥 2 𝑦´´+𝑥𝑦´+ 𝑥2 −4 𝑦 = 0 … 𝐸𝑐𝑢𝑎𝑐𝑖𝑜𝑛 𝑑𝑒 𝐵𝑒𝑠𝑠𝑒𝑙
x = 0 un punto singular. Hallaremos una solución en Serie de la forma:

y = n=0 cn x n+r


y =  ( n + r ) cn x n + r −1

n =0


y =  (n + r )(n + r − 1)cn x n+ r −2
n =0

8/10/2022 Mg. Nolan Jara Jara 27


∞ ∞ ∞ ∞
1
෍ 𝑛 + 𝑟 𝑛 + 𝑟 − 1 𝑐𝑛 𝑥 𝑛+𝑟 + ෍ 𝑛 + 𝑟 𝑐𝑛 𝑥 𝑛+𝑟 + ෍ 𝑐𝑛 𝑥 (𝑛+2)+𝑟 − ෍ 𝑐𝑛 𝑥 𝑛+𝑟 = 0
4
𝑛=0 𝑛=0 𝑛=0 𝑛=0
∞ ∞ ∞ ∞
1
෍ 𝑛 + 𝑟 𝑛 + 𝑟 − 1 𝑐𝑛 𝑥 𝑛+𝑟 + ෍ 𝑛 + 𝑟 𝑐𝑛 𝑥 𝑛+𝑟 + ෍ 𝑐𝑛−2 𝑥 𝑛+𝑟 − ෍ 𝑐𝑛 𝑥 𝑛+𝑟 = 0
4
𝑛=0 𝑛=0 𝑛=2 𝑛=0

1 1
𝑐0 𝑟 𝑟 − 1 + 𝑟 − 4 𝑥 𝑟 + 𝑐1 1 + 𝑟 𝑟 + 1 + 𝑟 − 4 𝑥 1+𝑟 +
∞ ∞ ∞ ∞
1
෍ 𝑛 + 𝑟 𝑛 + 𝑟 − 1 𝑐𝑛 𝑥 𝑛+𝑟 + ෍ 𝑛 + 𝑟 𝑐𝑛 𝑥 𝑛+𝑟 + ෍ 𝑐𝑛−2 𝑥 𝑛+𝑟 − ෍ 𝑐𝑛 𝑥 𝑛+𝑟 = 0
4
𝑛=2 𝑛=2 𝑛=2 𝑛=2

1 1
𝑟² − 4 𝑐0 𝑥 𝑟 + 1 + 𝑟 ² − 4 𝑐1 𝑥 1+𝑟 +
∞ ∞ ∞ ∞
1
෍ 𝑛 + 𝑟 𝑛 + 𝑟 − 1 𝑐𝑛 𝑥 𝑛+𝑟 + ෍ 𝑛 + 𝑟 𝑐𝑛 𝑥 𝑛+𝑟 + ෍ 𝑐𝑛−2 𝑥 𝑛+𝑟 − ෍ 𝑐𝑛 𝑥 𝑛+𝑟 = 0
4
𝑛=2 𝑛=2 𝑛=2 𝑛=2
1 1
𝑟 2 − 4 = 0; 1 + 𝑟 ² − 4=0
∞ ∞ ∞ ∞
1
෍ 𝑛 + 𝑟 𝑛 + 𝑟 − 1 𝑐𝑛 𝑥 𝑛+𝑟 + ෍ 𝑛 + 𝑟 𝑐𝑛 𝑥 𝑛+𝑟 + ෍ 𝑐𝑛−2 𝑥 𝑛+𝑟 − ෍ 𝑐𝑛 𝑥 𝑛+𝑟 = 0
4
𝑛=2 𝑛=2 𝑛=2 𝑛=2

8/10/2022 Mg. Nolan Jara Jara 28


1 −1 −3 −1
𝑟 = ± ;𝑟 = , ;𝑟 = ; 𝑐 ∈ 𝑅, 𝑐1 ∈ 𝑅
2 2 2 2 0
∞ ∞ ∞ ∞
1 3 1
𝑛−2 1 1
𝑛−2
1
𝑛−2 1 1
𝑛−2
෍ 𝑛− 𝑛− 𝑐 𝑥 +෍ 𝑛− 𝑐 𝑥 + ෍ 𝑐𝑛−2 𝑥 − ෍ 𝑐𝑛 𝑥 =0
2 2 𝑛 2 𝑛 4
𝑛=2 𝑛=2 𝑛=2 𝑛=2


1 3 1 1 1
𝑛−2
෍ 𝑛− 𝑛− + 𝑛− − 𝑐𝑛 + 𝑐𝑛−2 𝑥 =0↔
2 2 2 4
𝑛=2
1 3 1 1 1 1 1
𝑛− 𝑛 −2 + 𝑛 − 2 − 4 𝑐𝑛 + 𝑐𝑛−2 =0↔ 𝑛−2 𝑛 − 2 − 4 𝑐𝑛 + 𝑐𝑛−2 =0
2
𝒄𝒏−𝟐
↔ 𝑛(𝑛 − 1) 𝑐𝑛 + 𝑐𝑛−2 =0↔ 𝒄𝒏 = − 𝒏 𝒏−𝟏 ; 𝒏 = 𝟐, 𝟑, 𝟒, …
𝑐0 𝑐0 𝑐1 𝑐1 𝑐2 𝑐0
𝑛 = 2 → c2 =− =− ; 𝑛=3 → c3 =− =− ; 𝑛=4 → c4 =− =
2 2! 3(2) 3! (4)3 4!

𝑐3 𝑐1 𝑐4 𝑐 𝑐 𝑐0
𝑛=5 → c5 =− = ; 𝑛=6 → c6 =− = − 6!0 ; 𝑛 = 7 → c7 = − 7!1 ; 𝑛 = 8 → c8 =
5(4) 5! 6(5) 8!
𝑐1
n=9→ c9 = …
9!

8/10/2022 Mg. Nolan Jara Jara 29


1
∞ 𝑛− 2 1 1
𝑦 = σ𝑛=0 𝑐𝑛 𝑥 = 𝑥 σ∞
𝑛=0 𝑐𝑛 𝑥
𝑛
= [𝑐 + 𝑐1 𝑥 + 𝑐2 𝑥 2 + 𝑐3 𝑥 3 +
𝑥 0
𝑐4 𝑥 4 + 𝑐5 𝑥 5 + ⋯]
1 𝑐0 2 𝑐1 3 𝑐0 4 𝑐1 5 𝑐0 6 𝑐1 7
𝑦 = ቀ𝑐0 + 𝑐1 𝑥 − 𝑥 − 𝑥 + 𝑥 + 𝑥 − 𝑥 − 𝑥 +
𝑥 2! 3! 4! 5! 6! 7!
𝑐0 8 𝑐1 9
𝑥 + 𝑥 +⋯ቁ
8! 9!
𝑐0 𝑥2 𝑥4 𝑥6 𝑥8 𝑐1 𝑥3 𝑥5 𝑥7 𝑥9
𝑦= 𝑥
1− 2!
+ 4! − 6! + 8! +⋯ + 𝑥
𝑥− 3!
+ 5! − 7! + 9! +⋯

∞ 2𝑛 ∞
𝑐0 𝑛
𝑥 𝑐1 𝑛
𝑥 2𝑛+1
𝑦= ෍ −1 + ෍ −1
𝑥 2𝑛 ! 𝑥 2𝑛 + 1 !
𝑛=0 𝑛=0

𝑐𝑜𝑠𝑥 𝑠𝑒𝑛𝑥
𝑦= 𝑐0 + 𝑐1
𝑥 𝑥
8/10/2022 Mg. Nolan Jara Jara 30
∞ ∞ ∞ ∞
1
෍ 𝑛 + 𝑟 𝑛 + 𝑟 − 1 𝑐𝑛 𝑥 𝑛+𝑟 + ෍ 𝑛 + 𝑟 𝑐𝑛 𝑥 𝑛+𝑟 + ෍ 𝑐𝑛 𝑥 𝑛+𝑟+2 − ෍ 𝑐𝑛 𝑥 𝑛+𝑟 = 0
4
𝑛=0 𝑛=0 𝑛=0 𝑛=0
∞ ∞

෍𝑓 𝑛 = ෍𝑓 𝑛+𝑘
𝑛=𝑘 𝑛=0
∞ ∞

෍ 𝑛 + 2 + 𝑟 𝑛 + 1 + 𝑟 𝑐𝑛+2 𝑥 𝑛+𝑟+2 + ෍ 𝑛 + 2 + 𝑟 𝑐𝑛+2 𝑥 𝑛+2+𝑟


𝑛=−2 𝑛=−2
∞ ∞
𝑛+𝑟+2
1
+ ෍ 𝑐𝑛 𝑥 − ෍ 𝑐𝑛+2 𝑥 𝑛+2+𝑟 = 0
4
𝑛=0 𝑛=−2
1 1
𝑐0 𝑟 𝑟 − 1 + 𝑟 − 4 𝑥 𝑟 + 𝑐1 𝑟 1 + 𝑟 + 1 + 𝑟 − 4 𝑥1+𝑟 +
∞ ∞

෍ 𝑛 + 𝑟 + 2 𝑛 + 𝑟 + 1 𝑐𝑛+2 𝑥 𝑛+𝑟+2 + ෍ 𝑛 + 𝑟 + 2 𝑐𝑛+2 𝑥 𝑛+𝑟+2 +


𝑛=0 𝑛=0
∞ ∞
1
෍ 𝑐𝑛 𝑥 𝑛+𝑟+2 − ෍ 𝑐𝑛+2 𝑥 𝑛+𝑟+2 = 0
4
𝑛=0 𝑛=0
1 1
𝑟 𝑟 − 1 𝑐0 𝑥𝑟 + 𝑟𝑐0 𝑥𝑟 − 𝑐0 𝑥 𝑟 + 1 + 𝑟 𝑟𝑐1 𝑥 1+𝑟 + (1 + 𝑟)𝑐1 𝑥 1+𝑟 − 𝑐1 𝑥 1+𝑟 +
4 4
8/10/2022 Mg. Nolan Jara Jara 31
1 −1 −3 −1
𝑟 = ± , 𝑐0 ∈ 𝑅; 𝑟 = , ;𝑟 = ; 𝑐 ∈ 𝑅, 𝑐1 ∈ 𝑅
2 2 2 2 0
3 1 3
σ∞
n=0 n + 2 n + 2 cn+2 x n+r+2 + σ∞
n=0 n + 2 cn+2 x
n+r+2
+ σ∞
𝑛=0 𝑐𝑛 𝑥
𝑛+𝑟+2

1
− ෍ 𝑐𝑛+2 𝑥 𝑛+𝑟+2 = 0
4
𝑛=0

3 1 3 1
σ∞
𝑛=0 n+ n+ cn+2 + n + cn+2 − 𝑐𝑛+2 + 𝑐𝑛 𝑥 𝑛+𝑟+2 =0
2 2 2 4

3 2 1
σ∞
𝑛=0 n+2 − 4 cn+2 + 𝑐𝑛 𝑥 𝑛+𝑟+2 =0
𝑐𝑛
σ∞
𝑛=0 𝑛 + 2 𝑛 + 1 c n+2 + 𝑐𝑛 𝑥 𝑛+𝑟+2 =0 ↔ c
n+2 =- 𝑛+2 𝑛+1

𝑐0 𝑐0 𝑐1 𝑐1 𝑐2 𝑐0
𝑛 = 0, c2 =− =− ; 𝑛=1, c 3 =− =− ; 𝑛=2, c 4 =− =
2 2! 3(2) 3! (4)3 4!
𝑐3 𝑐1 𝑐4 𝑐0 𝑐1 𝑐0 𝑐1
𝑛=3, c5 =− = 5! ; 𝑛=4, c6 =− 6(5) = − 6! ; c7 = − 7! ; c8 = 8! ; c9 = 9! …
5(4)
8/10/2022 Mg. Nolan Jara Jara 32
1
𝑛− 2 1 1
𝑦 = σ∞ 𝑐
𝑛=0 𝑛 𝑥 𝑛+𝑟
= σ∞
𝑛=0 𝑛𝑐 𝑥 = σ∞ 𝑐𝑛 𝑥 𝑛
= [𝑐 + 𝑐1 𝑥 +
𝑥 𝑛=0 𝑥 0
𝑐2 𝑥 2 + 𝑐3 𝑥 3 + 𝑐4 𝑥 4 + 𝑐5 𝑥 5 + ⋯]

1 𝑐 𝑐 𝑐0 4 𝑐1 5 𝑐0 6 𝑐1 7
𝑦= ቀ𝑐0 + 𝑐1 𝑥 − 0 𝑥 2 − 1 𝑥 3 + 𝑥 + 𝑥 − 𝑥 − 𝑥 +
𝑥 2! 3! 4! 5! 6! 7!
𝑐0 8 𝑐1 9
8!
𝑥 + 9! 𝑥 + ⋯ ቁ

𝑐0 𝑥2 𝑥4 𝑥6 𝑥8 𝑐1 𝑥3 𝑥5 𝑥7 𝑥9
𝑦= 1 − 2! + 4! − 6! + 8! +⋯ + 𝑥− + − 7! + +⋯
𝑥 𝑥 3! 5! 9!

∞ ∞
𝑐0 𝑛
𝑥 2𝑛 𝑐1 𝑛
𝑥 2𝑛+1
𝑦= ෍ −1 + ෍ −1
𝑥 2𝑛 ! 𝑥 2𝑛 + 1 !
𝑛=0 𝑛=0

𝑐𝑜𝑠𝑥 𝑠𝑒𝑛𝑥
𝑦= 𝑐0 + 𝑐1
𝑥 𝑥
8/10/2022 Mg. Nolan Jara Jara 33
EJEMPLO 2 DEL MÉTODO DE FROBENIUS
3𝑥𝑦 ″ + 𝑦 ′ + 𝑦 = 0

x = 0 es un punto singular
hallaremos una solución en Serie de la forma:

𝑦 = ෍ 𝑐𝑛 𝑥 𝑛+𝑟
𝑛=0


y =  (n + r )cn x n+r −1
n =0

y =  (n + r )(n + r − 1)cn x n+r −2
n =0

8/10/2022 Mg. Nolan Jara Jara 34


3 xy + y − y
  
= 3 (n + r )(n + r − 1)cn x n+ r −1 +  (n + r )cn x n+ r −1 −  cn x n+ r
n =0 n =0 n =0
 
=  (n + r )(3n + 3r − 2)cn x n + r −1
−  cn x n+ r

n =0 n =0

 
  
n
= x r (3r − 2)c0 x +  (n + r )(3n + 3r − 2)cn x −  cn x 
r −1 n −1

 n =1
     n =0
 
 k = n −1 k =n 
 
k
= x r (3r − 2)c0 x +  [( k + r + 1)(3k + 3r + 1)ck +1 − ck ]x  = 0
r −1

 k =0 
8/10/2022 Mg. Nolan Jara Jara 35
Que implica:
r(3r – 2)c0 = 0
(k + r + 1)(3k + 3r + 1)ck+1 – ck = 0, k = 0, 1, 2, …
Haciendo c0 = 1, r(3r – 2) = 0. Entonces: r = 0, 2/3.
Cada valor de r nos transforma la ecuación de
recurrencia
ck
ck +1 = , k = 0, 1, 2, 
(k + r + 1)(3k + 3r + 1)
en:
ck
r1 = 2/3, ck +1 = , k = 0, 1, 2,…
(3k + 5)(k + 1)
ck
r2 = 0, ck +1 = , k = 0, 1, 2,…
8/10/2022
(k + 1)(3kMg.+Nolan
1) Jara Jara 36
r1 = 2/3, r2 = 0,
ck ck
ck +1 = , ck +1 = ,
(3k + 5)(k + 1) (k + 1)(3k + 1)
k = 0, 1, 2,… k = 0, 1, 2,…

c1 c0 c1 c0
c2 = = c2 = =
8.2 2! 5.8 2.4 2!1.4
c2 c0 c2 c0
c3 = = c3 = =
11.3 3! 5.8.11 3.7 3!1.4.7
c3 c0 c3 c0
c4 = = c4 = =
14.4 4! 5.8.11.14 4.10 4!1.4.7.10
 
c0 (−1) n c0
cn = cn =
n! 5.8.11  (3n + 2) n!1.4.7  (3n − 2)
8/10/2022 37
los dos conjuntos contienen el mismo múltiplo
c0. Si se omite este término, haciéndolo igual
a ,1 tenemos:
 
1 n
y0 ( x) = x 1 + 
2/3
x 
 n =1 n! 5 .
8.11 (3n + 2) 

 
1 n
y1 ( x) = x 1 + 
0
x 
 n =1 n!1 .
4.7  (3n − 2) 

que convergen para |x| <  y son linealmente


independientes. Así que la solución es:
y(x) = C0y0(x) + C1y1(x), 0 < x < 
8/10/2022 Mg. Nolan Jara Jara 38
EJEMPLO 3 DEL MÉTODO DE FROBENIUS
𝑥𝑦 ″ + 2𝑦 ′ + 𝑥𝑦 = 0

x = 0 es un punto singular
hallaremos una solución en Serie de la forma:

𝑦 = ෍ 𝑐𝑛 𝑥 𝑛+𝑟
𝑛=0

8/10/2022 Mg. Nolan Jara Jara 39

También podría gustarte