Notas 2 Ecuaciones Diferenciales
Notas 2 Ecuaciones Diferenciales
Notas 2 Ecuaciones Diferenciales
ECUACIONES DIFERENCIALES
ORDINARIAS DE PRIMERORDEN
dy
x + 2y = 0 xdx + (2xy −1)dy = 0
dx
MG NOLAN JARA JARA
ECUACIÓN DIFERENCIAL
Definición
EDO EDP
w w
dy
= x 2 y − cosy − = x2 + y
dx x y
y''' = y − y'
3
( , ( n)
F x, y, y y´´ ,..., y = 0 )
MG NOLAN JARA JARA
ORDEN
𝒚´´
xexy``` − + sen(xy)= 0 es de 3er orden.
𝒙
(𝑛)
an (x)𝑦 + 𝑎𝑛−1(𝑥)𝑦 (𝑛−1) +…+𝑎0 𝑥 𝑦 = g(𝑥)
( )
F x, y, y, y ,..., y(n ) = 0
1
función : y = EDO : xy'+ y =0
x
1 xy + y = 0
y = f (x) = 𝑥
es solución de:
1
𝑦= 2 es solución general de: y + 2xy 2
=0
𝑥 +𝑐
Si y(0)=-1
1
𝑦 = 𝑓(𝑥) = 2
𝑥 −1
x 2 + y 2 = 25 x
relación: EDO : y' =−
y
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é?
y´= 2x + 3y
y(1) = 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
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
𝒅𝒚
7. 𝒚 𝒚
𝒆 𝒄𝒐𝒔𝒙 + 𝒄𝒐𝒔𝒙 + 𝒆 𝒔𝒆𝒏𝒙 + 𝒆 𝒚
= 𝟎; 𝒚 𝝅 = 𝟎
𝒅𝒙
𝟏−𝒔𝒆𝒏𝒙
Solución. 𝒚 = 𝒍𝒏 𝟏+𝒔𝒆𝒏𝒙
𝑵𝟎 𝟏 𝟏 𝟏
= N(1500)= 𝑵𝟎 𝒆𝒌.𝟏𝟓𝟎𝟎 →𝒆 𝒌.𝟏𝟓𝟎𝟎 = →k= ln
𝟐 𝟐 𝟏𝟓𝟎𝟎 𝟐
−𝒍𝒏𝟐
−𝒍𝒏𝟐
k=𝟏𝟓𝟎𝟎 → N(t)=𝑵𝟎 𝒆𝟏𝟓𝟎𝟎𝒕
−𝒍𝒏𝟐
(𝟒𝟓𝟎𝟎) 𝑵 𝑵𝟎
N(4500)= 𝑵𝟎 𝒆 𝟏𝟓𝟎𝟎 = 𝑵𝟎 𝒆−𝟑𝒍𝒏𝟐 = 𝒆𝒍𝒏𝟖
𝟎
= =12.5% 𝑵𝟎
𝟖
1)y’= ax+by+c
2) y − e x e y +1 = 0
𝝏𝒇 𝝏𝒇
=𝑷; =𝑸
𝝏𝒙 𝝏𝒚
P Q
=
y x
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
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
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
u(x) = e g ( x )dx
Py − Q x
dx
u(x) = e
Q
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
f ( x, y) = f ( x, y)
n
Ejemplo:
f ( x, y) = x − y + xy
2 2
f( x, y) =
n f(x, y)
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
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
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
Un factor integrante
para una EDOL de
primer orden es:
u(x) = e p( x)dx
u(x) = e
=e p ( x)dx
(
6 ) 2 x d y = 2x 3− y d x
2y
)
7) y + = x3
x
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
2 2
− 𝑥𝑑𝑥 −2𝑙𝑛𝑥 1 1 ´ 4
y´− 𝑥y=4x; u(x)=𝑒 =𝑒 =𝑥 2 → 𝑦. 𝑥 2 = 𝑥
1
𝑦. 𝑥 2 =4lnx+c→ 𝑦 = 4𝑥²lnx+cx²
50 lt/min
Agua
E=
R=
L=
𝟐
Q(t)=𝒆−𝟔𝒕 𝒄𝟏 𝐜𝐨𝐬 𝟖𝒕 + 𝒄𝟐 𝒔𝒆𝒏 𝟖𝒕 − 𝒄𝒐𝒔(𝟏𝟎𝒕)
𝟓
𝟒 𝟑
De las condiciones iniciales encontramos. 𝒄𝟏 = ; 𝒄𝟐 = 𝟏𝟎
𝟏𝟎
Donde, la solución requerida es:
𝒆−𝟔𝒕 𝟐
Q(t)= 𝟒 𝐜𝐨𝐬 𝟖𝒕 + 𝟑𝒔𝒆𝒏 𝟖𝒕 − 𝒄𝒐𝒔(𝟏𝟎𝒕)
𝟏𝟎 𝟓
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)
3
dy tan x
Resuelva la EDO: =
dx y
y dy = tan 3 x d x
𝑅
𝐸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
𝑒 − 𝑡
𝐿
𝑅 −𝑎𝐿
𝑑 𝑛𝑦 𝑑 𝑛−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
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 laecuacion
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 la solución general es :
1.3)
1.4)
ത ത
orden𝑛 𝑥 cosbx+ 𝑄𝑚 𝑥 senbx, 𝑟 ≠ ±𝑏𝑖;𝑦𝑝 = 𝑃𝑘 𝑥 𝑐𝑜𝑠𝑏𝑥 + 𝑄𝑘 𝑥 𝑠𝑒𝑛𝑏𝑥
3) g(x)=𝑃
g(x)=𝑃𝑛 𝑥 cosbx+ 𝑄𝑚 𝑥 senbx, 𝑟 = ±𝑏𝑖;𝑦𝑝 = 𝑥 𝑠 ( 𝑃ത𝑘 𝑥 𝑐𝑜𝑠𝑏𝑥 + 𝑄ത𝑘 𝑥 𝑠𝑒𝑛𝑏𝑥 )
s:grado de la raiz r = ib de la ecuacion característica;k = maxm,n
vi
vi
vi
vi
vi
vi
f1 f2
… fn
f1 ' f2 '
… fn '
W ( f1,..., f n ) =
…
… …
f (n−1)
1
f (n−1)
2
… f n(n−1)
...........................................................................
( 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
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)
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)
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
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
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
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.
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
23/06/2021
06/07/2020 Mg NOLAN
NOLAN JARAJARA
JARA 42
39
JARA
EJEMPLO 1
∞ 𝑏
𝐿 𝑓 𝑡 = න 𝑒 −𝑠𝑡 𝑓 𝑡 𝑑𝑡 = lim න 𝑒 −𝑠𝑡 𝑓 𝑡 𝑑𝑡
0 𝑏→∞ 𝑡=0
Donde sR
Se dice que la transformada de Laplace de f(t) existe si la integral
converge para algunos valores de s.
∞ 𝑏
න 𝑒 −𝑠𝑡 𝑓 𝑡 𝑑𝑡 = 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
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
𝒏 + 𝒏
𝒏!
𝑺𝒊 𝒇 𝒕 = 𝒕 ; 𝒏 ∈ 𝒁𝟎 ⇒𝑳 𝒕 = = 𝑭 𝒔 ;𝒔 > 𝟎
𝒔𝒏+𝟏
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
F ( s) = L e = at
0
e e dt = e
at − st
0
− ( s − a )t
dt
− ( s − a )t
−e 1
= , sa
( s − a) 0
s−a
𝒂𝒕 𝒂𝒕
𝟏
∴ 𝑺𝒊 𝒇 𝒕 = 𝒆 ⇒𝑳 𝒆 = = 𝑭(𝒔); 𝒔 > 𝒂
𝒔−𝒂
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 ) = Lsent =
1
;s 0
s +1
2
f (t ) = cos(t ) F ( s ) = Lcos t =
s
Si ;s 0
s +1
2
𝒔
𝑺𝒊 𝒇 𝒕 = 𝒄𝒐𝒔 𝒂𝒕 ⇒ 𝑳 𝒄𝒐𝒔(𝒂𝒕) = 𝟐 𝟐
= 𝑭(𝒔); 𝒔 > 𝟎
𝒔 +𝒂
Mg. Nolan Jara Jara 14
LA FUNCIÓN HEAVISIDE O ESCALÓN UNITARIO
0 , 0 t a
u(t − a ) =
1 , ta
0; si t a
u (t − a) =
1
1; si t a
Lu (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
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
= F ( s − a); s a
Ejemplo:
L t = 2 L e t =
s
1 at
1
( s − a) 2
0
a0
1 s
= F( )
a a
𝟏 𝒔
𝑳{𝒇(𝒕)} = 𝑭(𝒔) ⇒ 𝑳{𝒇(𝒂𝒕)} = 𝑭
𝒂 𝒂
= − න 𝑒 −𝑠𝑡 𝑡𝑓(𝑡) 𝑑𝑡
0
= −𝐿 𝑡𝑓(𝑡)
𝑺𝒊: 𝑭 𝒔 = 𝑳 𝒇 𝒕 ⇒ 𝑭′ 𝒔 = −𝑳 𝒕𝒇 𝒕 ; 𝑳 𝒕𝒇 𝒕 = −𝑭′ 𝒔
Si : F ( s ) = L{ f (t )} L t f (t ) = (−1) F ( s )
3 3
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
∞
𝒇(𝒕)
𝑳 = න 𝑭(𝒖)𝒅𝒖
𝒕 𝒔
𝜏=t
t→ ∞
t=0
R
t
𝜏=0
𝒕
𝑳 𝒇 𝟎 𝝉 𝒈 𝒕 − 𝝉 𝒅𝝉 =𝑭 𝒔 . 𝑮 𝒔 = 𝑳{𝒇(𝒕) ∗ 𝒈(𝒕)} →
𝒕
𝒇 𝟎 𝝉 𝒈 𝒕 − 𝝉 𝒅𝝉 =𝑳−𝟏 𝑭 𝒔 . 𝑮(𝒔) =f(t)*g(t)
∞
la integral න 𝑓(𝜏)𝑔(𝑡 − 𝜏)𝑑𝜏
−∞
y se denota como
∞
න 𝑓(𝜏)𝑔(𝑡 − 𝜏)𝑑𝜏 = 𝑓(𝑡) ∗ 𝑔(𝑡)
−∞
t f ( ) g (t − )d , t 0
f (t ) * g (t ) = 0
0, t0
Así que para estas funciones podemos definir la convolución
como:
𝒕
𝒇(𝒕) ∗ 𝒈(𝒕) = න 𝒇(𝝉)𝒈(𝒕 − 𝝉)𝒅𝝉, (𝒕 ≥ 𝟎)
𝟎
𝒕
𝑳 𝒇 𝟎 𝝉 𝒈 𝒕 − 𝝉 𝒅𝝉 =𝑭 𝒔 . 𝑮 𝒔 = 𝑳{𝒇(𝒕) ∗ 𝒈(𝒕)} →
𝒕
𝒇 𝟎 𝝉 𝒈 𝒕 − 𝝉 𝒅𝝉 =𝑳−𝟏 𝑭 𝒔 . 𝑮(𝒔) =f(t)*g(t)
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
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
𝑘 −1s
(f) senh(kt) = L−1 (g) cosh kt = L 2 2
𝑠2 − 𝑘2 s − k
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
−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 𝑡)
Luego
−1 s 2
+ 6s + 9
L
( s − 1)( s − 2)( s + 4)
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
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)
𝑷𝟏 ) 𝒚´´ + 𝟐𝒚´ + 𝟐𝒚 = 𝜹 𝒕 − 𝝅 ; 𝒚 𝟎 = 𝟎; 𝒚´ 𝟎 = 𝟎
𝟏
𝑳−𝟏 (𝒔+𝟏)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, 𝑡>𝜋
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)²
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
1
2
𝑠𝑒𝑛𝑡 − 𝑡𝑐𝑜𝑠𝑡 ; 0 ≤ 𝑡 < 𝜋
y(t)=ቐ π
− cost; 𝑡 ≥𝜋
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
𝑦 ″ + 𝑦 = sen 𝑡 …E
𝑦(0) = 1
𝑦 ′ (0) = 0
Solució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
1 𝑡
𝑦(𝑡) = cos 𝑡 + sen 𝑡 − cos 𝑡
2 2
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
1 − 𝑒 −𝑡 , 0≤𝑡<1
−𝑒 −𝑡 + 𝑒 −(𝑡−1) , 1≤𝑡<2
𝑖(𝑡) =
1 − 𝑒 −𝑡 + 𝑒 −(𝑡−1) − 𝑒 −(𝑡−2) , 2≤𝑡<3
−𝑒 −𝑡 + 𝑒 −(𝑡−1) − 𝑒 −(𝑡−2) + 𝑒 −(𝑡−3) , 3≤𝑡<4
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)
: 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
(1) = 1
( x + 1) = x ( x )
( n + 1) = n !
1
( ) =
2
1 − 1 −t
( ) = t 2 e dt =
2 0
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
Lt =
1 2
3
2
2s
a
n =1
n
o a n
2. (−1) n
n =1
1
3.
n =1 n( n + 1)
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.
2
n =1
2n
3 1− n
Si la serie a
n =1
n es convergente, entonces:
lím an = 0
n→
a
n =1
n es divergente.
n2
n =1 5n 2
+4
divergente
(i) ca
n =1
n = c an
n =1
(ii) (a
n =1
n + bn ) = an + bn
n =1 n =1
(− 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 →
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
Una serie a
n =1
n se llama condicionalmente convergente
(− 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
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.
3n
es
absolutamente convergente.
nn
2. Diga si la serie
n = 1 n!
es convergente.
c (x − a ) = c0 + c1 (x − a ) + c2 (x − a ) + ...
k 2
k
k =0
c ( x − a)
k
Para una serie de potencias k
k =0
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
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
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
cn +1 ( x − a )
n +1
si L 1, L= n lim
cn ( x − a )
n
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
lím Rn ( x) = 0
n →
Rn (x) es el residuo de orden n.
2. Si en la Serie de Taylor a = 0,
2 k −1
sen x = (− 1)
x
x
k
;
k =0 (2k − 1)!
¿Qué serie le corresponderá a cos x = ?
2) Desarrolle f en serie de
Taylor alrededor de / 3 :
f ( x) = cos x
Definición:
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
n =0
n(n − 1)c x
n=2
n
n−2
+ cn x
n =0
n +1
n=2 n =1 n =0
y ' '+ xy = 0
y = cn (x − x0 )
n
n =0
n(n − 1)c
n=2
n x n−2
+ cn x
n =0
n +1
a x n
n
:
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 + 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
𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛 + (𝑛 + 2) 𝑛 + 1 𝑐𝑛+2 𝑥 𝑛 + 𝑛𝑐𝑛 𝑥 𝑛 − 𝑐𝑛 𝑥 𝑛 =
𝑛=2 𝑛=2 𝑛=2 𝑛=2
2c2 + 6c3 x − c0 + σ∞ 𝑛
𝑛=2 𝑛 𝑛 − 1 𝑐𝑛 + 𝑛 + 2 𝑛 + 1 𝑐𝑛+2 + 𝑛𝑐𝑛 − 𝑐𝑛 𝑥 =
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
𝑛 𝑛 − 1 𝑐𝑛 𝑥 𝑛 + (𝑛 + 2) 𝑛 + 1 𝑐𝑛+2 𝑥 𝑛 + 𝑛𝑐𝑛 𝑥 𝑛 − 𝑐𝑛 𝑥 𝑛 = 0
𝑛=2 𝑛=2 𝑛=2 𝑛=2
∞
𝑐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
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
( 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
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
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:
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 )
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
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
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
∞
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!
∞ 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
∞ ∞
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
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)
x = 0 es un punto singular
hallaremos una solución en Serie de la forma:
∞
𝑦 = 𝑐𝑛 𝑥 𝑛+𝑟
𝑛=0