Ecuaciones en Diferencias

Dr. Hipólito Aguilar Sierra

Señales y Sistemas
November 19, 2020

Ejemplo 1: Encuentre la solución a la siguiente ecuación en diferencias

f [n + 2] 4f [n + 1] + 3f [n] = 0

con

f [0] = 0
f [1] = 1

Solución: Buscamos la ecuación característica

rn+2 4rn+1 + 3rn = 0

rn r2 4r + 3 = 0
r2 4r + 3 = 0

Se buscan las raíces del polinomio anterior, así encontramos que

r1 = 3
r2 = 1

Planteamos la solución:

f [n] = C1 r1n + C2 r2n

f [n] = C1 3n + C2 1n
f [n] = C1 3n + C2

Evaluamos las condiciones de frontera para determinar las constantes C1 y C2

f [0] = C1 3(0) + C2 = 0
= C1 + C2 = 0
! C1 = C2
f [1] = C1 3(1) + C2 = 1
= 3C1 + C2 = 1

remplazando C1

3C1 + C2 = 1
3 ( C2 ) + C2 = 1
2C2 = 1
1
! C2 =
2
1
! C1 =
2

1
 Finalmente, sustituimos los valoes de C1 y C2 en la solución planteada

f [n] = C1 3n + C2
1 n 1
f [n] = 3
2 2
3n 1
f [n] =
2
Ejemplo 2: Encuentre la solución a la siguiente ecuación en diferencias

y[n + 2] y[n + 1] + y[n] = 0

con

y[0] = 0
y[1] = 1

Solución: Buscamos la ecuación característica

rn+2 rn+1 + rn = 0
rn r2 r + 1 = 0
r2 r+1 = 0

Se buscan las raíces del polinomio anterior, así encontramos que

p
1 3
r1 = +j
2 p2
1 3
r2 = j
2 2
Utilizando la identidad de Euler tenemos que

jr1 j = jr2 j = jrj

v
u
u 1 2 p !2
3
jrj = t +
2 2
jrj = 1

por otro lado

0p 1
3
B 2 C
' = B
arctan @ C
1 A
2
=
3
Entonces, planteamos la solución
n n
y[n] = C1 jrj cos (n') + C2 jrj sin (n')
n n
y[n] = C1 1n cos + C2 1n sin
3 3
n n
y[n] = C1 cos + C2 sin
3 3

2
 Evaluamos las condiciones de frontera para determinar las constantes C1 y C2

(0) (0)
y[0] = C1 cos + C2 sin =0
3 3
= C1 (1) + C2 (0) = 0
! C1 = 0
(1) (1)
y[1] = C1 cos + C2 sin =1
3 3
= (0) cos + C2 sin =1
p 3 3
3
= C2 =1
2
2
! C2 = p
3
Finalmente, sustituimos los valoes de C1 y C2 en la solución planteada
n n
y[n] = C1 cos+ C2 sin
3 3
n 2 n
y[n] = (0) cos + p sin
3 3 3
2 n
y[n] = p sin
3 3

Ejemplo 3: Encuentre la solución a la siguiente ecuación en diferencias

f [n + 3] 7f [n + 2] + 15f [n + 1] 9f [n] = 0

con

f [0] = 1
f [1] = 2
f [2] = 17

Solución: Buscamos la ecuación característica

rn+3 7rn+2 + 15rn+1 9rn = 0

rn r3 7r2 + 15r 9 = 0
r3 7r2 + 15r 9 = 0

Se buscan las raíces del polinomio anterior, así encontramos que

r1 = 1
r2 = 3
r3 = 3

Entonces, planteamos la solución

f [n] = C1 r1n + C2 r2n + C3 nrnn

f [n] = C1 1n + C2 3n + C3 n3n
f [n] = C1 + C2 3n + C3 n3n

3
 Evaluamos las condiciones de frontera para determinar las constantes C1 , C2 y C3
f [0] = C1 + C2 3(0) + C3 (0) 3(0) = 1
= C1 + 3C2 + 0 = 1
! C1 = 1 3C2
f [1] = C1 + C2 3(1) + C3 (1) 3(1) = 2
= C1 + 3C2 + 3C3 = 2
= ( 1 3C2 ) + 3C2 + 3C3 = 2
= 1 + 3C3 = 2
! C3 = 1
f [2] = C1 + C2 3(2) + C3 (2) 3(2) = 17
= C1 + 9C2 + 18C3 = 17
= ( 1 3C2 ) + 9C2 + 18 (1) = 17
= 6C2 + 17 = 17
! C2 = 0
! C1 = 1
Finalmente, sustituimos los valoes de C1 , C2 y C3 en la solución planteada
f [n] = C1 + C2 3n + C3 n3n
f [n] = ( 1) + (0) 3n + (1) n3n
f [n] = 1 + n3n
f [n] = n3n 1
Ejemplo 4: Encuentre la solución a la siguiente ecuación en diferencias
y[n + 2] 2y[n + 1] + y[n] = n + 1

Solución: Aplicando transformada Z a ambos lados y con condiciones de frontera iguales a cero, tenemos
z 2 Y (z) 2zY (z) + Y (z) = Zfn + 1g
2
z 2z + 1 Y (z) = Zfn + 1g
Usando las tablas de transformada encontramos que el miembro del lado derecho es:
z2
Zfn + 1g = 2
(z 1)
z
= z 2
(z 1)
Remplazando
z
z2 2z + 1 Y (z) = z 2
(z 1)
Y (z) z
= 2
z (z 1) (z 2 2z + 1)
Y (z) z
= 4
z (z 1)
Utilizando la expansión en fracciones parciales
z A1 A2 A3 A4
4 = + 2 + 3 + 4
(z 1) (z 1) (z 1) (z 1) (z 1)
Y (z) A1 A2 A3 A4
= + 2 + 3 + 4
z (z 1) (z 1) (z 1) (z 1)

4
 4
Multiplicando ambas partes por (z 1)
3 2
z = A1 (z 1) + A2 (z 1) + A3 (z 1) + A4

Cuando z = 1
3 2
1 = A1 ((1) 1) + A2 ((1) 1) + A3 ((1) 1) + A4
! A4 = 1

Remplazando A4 = 1
3 2
z = A1 (z 1) + A2 (z 1) + A3 (z 1) + 1
z = A1 z 3 3A1 z 2 + 3A1 z A1 + A2 z 2 2A2 z + A2 + A3 z A3 + 1
z = A1 z 3 (3A1 A2 ) z 2 + (3A1 2A2 + A3 ) z A1 + A2 A3 + 1

Igualando los coe…cientes de los terminos del lado izquierdo con los del lado derecho, tenemos

A1 = 0
(3A1 A2 ) = 0
3A1 2A2 + A3 = 1
A1 + A2 A3 + 1 = 0

Remplazando A1 = 0

(3A1 A2 ) = 0
(3 (0) A2 ) = 0
! A2 = 0

Reamplazando A1 = A2 = 0

3A1 2A2 + A3 = 1
3 (0) 2 (0) + A3 = 1
! A3 = 1

Remplazando los valores de Ai en la expasión en fracciones parciales

Y (z) A1 A2 A3 A4
= + 2 + 3 + 4
z (z 1) (z 1) (z 1) (z 1)
Y (z) (0) (0) (1) (1)
= + + +
z (z 1) (z 1)2 (z 1)
3
(z 1)
4

Y (z) 1 1
= 3 + 4
z (z 1) (z 1)
z z
Y (z) = 3 + 4
(z 1) (z 1)

De las tablas de transformadas Z comunes, sabemos que

m+1 m
Y2
1 z kn
Z m = (n h)
(z k) (m 1)!
h=0

5
 Para nuestro caso, k = 1 y m = 3 y m = 4. Entonces para m = 3, tenemos
( ) 3 2
1 z 1n 3+1 Y
Z 3 = (n h)
(z 1) (3 1)!
h=0
1
1n+2 Y
= (n h)
(2)!
h=0
1
= (n (n 1))
2
1 2
= n n
2
Para m = 4
( )
4+1 4Y2
1 z 1n
Z 4 = (n h)
(z 1) (4 1)!
h=0
2
1n+3 Y
= (n h)
(3)!
h=0
1
= (n (n 1) (n 2))
6
1 3
= n 3n2 + 2n
6
Sustituyendo ambas transformadas Z inversa
1 2 1 3
y[n] = n n + n 3n2 + 2n
2 6
1 2 1 1 3 2 2
y[n] = n n + n3 n + n
2 2 6 6 6
1 3 1 3 2 2 1
y[n] = n + n + n
6 2 6 6 2
1 3 1
y[n] = n n
6 6
1 3
y[n] = n n
6
Ejemplo 5: Encuentre la solución a la siguiente ecuación en diferencias
y[n + 2] y[n + 1] + y[n] = 0

con
y[0] = 0
y[1] = 1
Solución: Aplicando transformada Z a ambos lados, tenemos
z 2 Y (z) z 2 y[0] zy[1] zY (z) + zy[0] + Y (z) = 0
Evaluando las condiciones de frontera y agrupando los términos de Y (z)
z 2 Y (z) z 2 (0) z (1) zY (z) + z (0) + Y (z) = 0
z 2 Y (z) z zY (z) + Y (z) = 0
z 2 z + 1 Y (z) z = 0
z2 z + 1 Y (z) = z
z
Y (z) =
z2 z+1
Y (z) 1
=
z z2 z+1

6
Como sabemos las raíces son:
p
1 3
z1 = +i
2 p2
1 3
z2 = i
2 2
Aplicando expansión en fracciones parciales:
Y (z) 1
=
z (z z1 ) (z
z2 )
A1 A2
= +
(z z1 ) (z z2 )
Determinamos los coe…cientes A1 y A2
Y (z)
A1 = (z z1 )
z z=z1
1
= (z z1 ) p
(z z1 ) (z z2 ) z= 12 +i 3
2

1
= p p
1 3 1 3
2 +i 2 2 i 2
1 p
= i 3
3
Y (z)
A2 = (z z2 )
z z=z2
1
= (z z2 ) p
(z z1 ) (z z2 ) z= 12 i 3
2

1
= p p
1 3 1 3
2 i 2 2 +i 2
1 p
= i 3
3
Escribiendo las raices en su forma polar
z1 = ei 3
z2 = e i3
Remplazando los coe…cientes y las raices en su forma polar
1
p 1
p
Y (z) 3i 3 3i 3
= +
z z ei 3 z e i3

Y (z) 1 p 1 1
= i 3 +
z 3 z ei 3 z e i3

1 p z z
Y (z) = i 3 +
3 z ei 3 z e i3

1
z
Mutiplicando al lado derecho por 1
z
1 p z z 1
z z 1
Y (z) = i 3 +
3 z ei 3 z 1 z e i3 z 1

1 p 1 1
= i 3 +
3 1 ei 3 z 1 1 e i3 z 1

7
De las tablas de transformadas Z conocidas tenemos que
1
Zfan u[n]g = 1
1 az

Para nuestro caso si a1 = ei 3 y a2 = e i3

la transformada Z inversa de Y (z), resulta en

1 1 p n i3 n
Z fY (z)g = i 3 ei 3 u[n] + e u[n]
3
1 p i n3 n
y[n] = i 3 e ei 3 u[n]
3
Utilizando la identidad de Euler
1 p n n n n
y[n] = i 3 cos i sin cos i sin u[n]
3 3 3 3 3
1 p n
y[n] = i 3 2i sin u[n]
3p 3
2 3 n
y[n] = sin u[n]
3 3