Señales y Sistemas en Tiempo Discreto
Señales y Sistemas en Tiempo Discreto
Trabajo Práctico N° 3
Módulo C: Señales y Sistemas en Tiempo Discreto
Integrantes:
Fiorentino, Ignacio Legajo: 96367
Ressia, Lorenzo Legajo: 94609
Badano Brito, Valentín Legajo: 96804
Pappalardo, Santiago David Legajo: 94606
Sisterna, José Ignacio Legajo: 97491
Martínez, Mariano Farid Legajo: 81690
Curso: 2R4
Docentes: Modesti, Mario Benasulin, Dimas
Índice general
1.2. b) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.3. c) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
3.2. b) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3.3. c) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3.4. d) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
4.2. b) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
5. Espectro de Frecuencias 17
5.1. a) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
5.2. b) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
6. Sistema LTI 21
6.1. a) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
6.2. b) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
6.3. c) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
6.4. d) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
6.5. e) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
6.6. f) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
6.7. g) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
6.8. h) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2
ÍNDICE GENERAL 3
Rúbrica:
Ejercicio 1
1.1. a)
Una señal analógica, en tiempo continuo x(t) = e−2t µ(t) se muestrea para generar la se-
cuencia de tiempo discreto x[n], considerar n = nTs , con Ts = 0,05, Ts = 0,01, Ts = 0,1.
4
1.2. B) 5
1.2. b)
Una señal analógica, en tiempo continuo x(t) = 10cos(2t) µ(t) se muestrea para generar la
secuencia de tiempo discreto x[n], con: Ts = 0,1 y Ts = 0,01.
De misma forma que para 1.a) :
1.3. c)
En este caso, podemos aplicar una técnica de Regresión Polinómica computarizada, ya que
podemos prescindir de las herramientas de análisis y álgebra. Para esto, introducimos la se-
cuencia en un array que almacena los valores de imagen y denimos uno correspondiente a los
valores del dominio discreto n, con lo cual nos generará un polinomio de grado 10, igual que el
tamaño de la muestra.
6 EJERCICIO 1. MUESTREO DE SEÑALES Y REGRESIÓN POLINÓMICA
Ya que la secuencia de datos es lo único que sabemos sobre la señal original, este polinomio
es la mejor aproximación con la que podemos contar. Si modicáramos el grado del polinomio,
nos veremos obligados a reducir la precisión con la que representa los datos adquiridos, o
el polinomio perderá capacidad de aproximación a futuro. La Regresión Polinómica permite
ajustarse también de acuerdo al proceso de adquisición de datos, lo cual, no es sensato aplicar
o asumir en este escenario.
Ejercicio 2
Planteo inicial
El sistema h(n) se obtiene de dos sistemas en cascada:
Donde:
Entonces:
h(n) = ((3− nu(n) ∗ 3− nu(n)) − (4n u(n)) ∗ 2n u(n))
Transformada Z
Aplicando la transformada Z:
z z z z
h1 (z) = , h2 (z) = , h3 (z) = , h4 (z) =
z − 1/3 z − 1/3 z − 1/2 z−1
Por tanto:
z2 z z
h(z) = ( − )·
(z − 1/3)(z − 1/3) z − 1/2 z − 1
Simplicando:
z2
h(z) =
(z − 1)(z − 1/3)2
h(z) z z
= 2
− 1
z (z − 1/2)(z − 1/3) (z − 2 )(z − 1)
7
8 EJERCICIO 2. CONVOLUCIÓN TIEMPO DISCRETO
Queremos expresar:
h(z) z z
= 1 2 1 − 1
z (z − 3 ) (z − 2 )(z − 1) (z − 2 )(z − 1)
z2 z
(a) (b) −
(z − 31 )2 (z − 1) (z − 1
2
)(z − 1)
Parte (a)
Queremos expresar:
z2 A B C
1 2 = 1 2 + 1 +
(z − 3 ) (z − 1) (z − 3 ) z− 3
z−1
2
z=1: C = (1 − 31 )2 = 23 = 94
1 (1/3)2
z= 3
: A= = (se obtiene por derivación, omitir cálculo directo)
(1/3 − 1)(1/3 − 1/3)2
z = otros valores: B = por comparación de coecientes.
(El desarrollo completo se deja indicado en la hoja: los valores nales aproximados eran
A = − 61 , B = − 54 , C = 49 ).
h(z) 1 5 9
=− 1 2 − 1 +
2 6(z − 3 ) 4(z − 3 ) 4(z − 1)
Parte (b)
z A B
− 1 = 1 +
(z − 2
)(z − 1) z− 2
z−1
Multiplicando por el denominador común:
−z = A(z − 1) + B(z − 12 )
Evaluando:
z = 1 : −1 = A(1 − 1) + B(1 − 12 ) ⇒ B = −2
z = 12 : − 12 = A( 21 − 1) ⇒ A = −1
z −1 −2
− 1 = 1 +
(z − 2
)(z − 1) z−2 z−1
9
Resultado nal
Sumando las dos partes:
h(z) 1 5 9 1 2
=− 1 2 − 1 + − 1 +
z 6(z − 3 ) 4(z − 3 ) 4(z − 1) z − 2
z−1
n n n
1 1 5 1 9 1
H[n] = − n µ[n] − µ[n] + µ[n] − µ[n] + (2µ[n])
2 3 4 3 4 2
n k k k !
X 1 1 5 1 7 1
Y [n] = − k − + − +2
k=0
4 3 4 3 4 2
n k n k n k X n
1X 1 5X 1 7X 1
Y [n] = − k − + − + 2
4 k=0 3 4 k=0 3 4 k=0 2 k=0
" n 1 n+1
# " n+1 # "
1 n+1
#
1 1 − (n + 1) 31 + n 5 1 − 13
3 7 1 − − 2
Y [n] = − − + + 2(n + 1)
4 (1 − 13 )2 4 1− 1
3
4 1 − (− 12 )
z z z
H(z) = 1 − 1 · 1
z− 2
z− 3
z− 3
z z − 31 − z z − 21
z
H(z) = 1 1 · 1
(z − 2 )(z − 3 ) z− 3
z z − 31 z z − 21
H(z) = −
(z − 12 )(z − 13 )2 (z − 12 )(z − 13 )2
z z
H(z) = 1 1 −
(z − 2
)(z − 3
) (z − 31 )2
z z
H(z) = 1 1 −
(z − 2
)(z − 3
) (z − 31 )2
(a)
z A B
1 1 = 1 + 1
(z − 2
)(z − 3
) z− 2
z− 3
z = A(z − 13 ) + B(z − 12 )
1
Para z= 2
:
1
1 1 1 2
2
=A 2
− 3
⇒ A= 1 =3
6
1
Para z= 3
:
1
1 1 1 3
=B − ⇒ B= = −2
3 3 2
− 16
z 3 2
⇒ 1 1 = 1 − 1
(z − 2
)(z − 3
) z− 2
z− 3
(b)
z A B
1 2 = 1 2 + 1
(z − 3 ) (z − 3 ) z− 3
z = A + B(z − 13 )
1
Para z= 3
:
1 1
3
=A ⇒ A= 3
Para z = 1:
1
1= 3
+ B(1 − 13 ) ⇒ B=1
1
z 3 1
⇒ 1 2 = 1 2 +
(z − 3 ) (z − 3 ) z − 31
3z 2z 1 z z
H(z) = − − · +
z − 12 z − 13 3 (z − 13 )2 z − 13
z z 1 z z
H(z) = 3 · 1 −2· 1 − · 1 2 − 1
z− 2
z− 3
3 (z − 3 ) z− 3
n
X
y[n] = h[k]
k=0
n k n k n k X n k
X 1 X 1 X 1 1
y[n] = 3 −2 − k −
k=0
2 k=0
3 k=0
3 k=0
3
n
! n
!
X 1 − rn+1 X r − (n + 1)rn+1 + nrn+2
rk = krk =
k=0
1−r k=0
(1 − r)2
n+1 n+1
1 − 21 1 − 13
y[n] = 3 · 1 −2· 1
1− 2
1− 3
n+1
h
1 n+1
i
1 n+2 1 1 − 13
− (n + 1) −n −
3 3
(1 − 31 )2 1− 1
3
1 n+1
1 n+1 2 1− 1
n+1 1 n+2
n+1
− (n + 1) 31 1 − 13
3 3 2 3 3
+n 3
y[n] = 1 − 1 − 1 −2· 1 2 + 1
2 2
1− 3
(1 − 3 ) 1− 3
1 n+1
1 n+1
3 1 n+1
3 2 n+1
3 2 1 n+1
y[n] = 6 − 6 2
−3+3 3
+ (n + 1) 3
− 3
+ + 3
2 2 2 3
3 1 n+1
1 n+1
3 1 n+1
3 2 n+1
y[n] = −6 2
−3 3
+ (n + 1) 3
− 3
2 2 4
3 7 1 n+1
3 1 n+1
3 2 n+1
1 n+1
y[n] = − 3
+ (n + 1) 3
− 3
−6 2
2 3 2 4
Ejercicio 3
3.1. a)
Identicar en el problema las variables de entrada y salida del mismo. Identicación de
variables.
3.2. b)
Plantear un modelo en ecuaciones de diferencias capaz de caracterizar el sistema propuesto.
En cada periodo, el saldo anterior y[n − 1] acumula interés y luego se suma el depósito x[n].
Por tanto la ecuación en diferencias es
12
3.2. B) 13
yh [n] − (1 + r)yh [n − 1] = 0,
C − (1 + r)C = 10 =⇒ C 1 − (1 + r) = 10.
10
−rC = 10 =⇒ C=− .
r
Para r = 0,02:
10
C=− = −500.
0,02
Por lo tanto la solución general es
10
y[n] = A(1 + r)n −
r
10 10
y[0] = A(1 + r)0 − =⇒ 10 = A − .
r r
Despejando A:
10
A = 10 +
r
Con r = 0,02:
10
A = 10 + = 10 + 500 = 510
0,02
Finalmente, la solución explícita es
y sustituyendo r = 0,02:
y[n] = 510 (1,02)n − 500
14 EJERCICIO 3. ECUACIONES EN DIFERENCIAS - SISTEMAS
3.3. c)
Representar la EdeD por medio de un diagrama de bloques con retardos unitarios.
3.4. d)
Evaluar el saldo para un periodo de capitalización de 48 meses.
Usando la fórmula obtenida,
Calculamos numéricamente:
(1,02)48 ≈ 2,5870703855
510 · (1,02)48 ≈ 1319,4059
y[48] ≈ 1319,4059 − 500 ≈ 820
Por tanto el saldo aproximado al cabo de 48 meses es
y[48] ≈ $ 820 .
Ejercicio 4
El objetivo de este Trabajo Práctico es poner en valor los modos de expresar la ecuación de
diferencias que caracteriza un sistema en tiempo discreto y expresar la síntesis de la expresión
matemática en lenguaje de diagrama de bloques. Considerar este último como la expresión
mínima que permita desarrollar un algoritmo de cálculo numérico con la menor cantidad de
operaciones posibles.
4.1. a)
Diagramar en diagrama de bloques normalizado representativos del sistema para la siguiente
ecuación en diferencias:
Para obtener el diagrama de bloques normalizado, primero debemos obtener la salida actual
del sistema, puesto que en la ecuación dada la salida está adelantada en [n + 1].
Utilizando el operador de retardo D tal que Dy[n] = y[n − 1]:
2 1 1
y[n] = Dx[n] + D 2
x[n] − 2
y[n]
h i
1
y[n] = D x[n] + D 2
x[n] − 12 y[n]
15
16 EJERCICIO 4. ECUACIONES EN DIFERENCIAS - DIAGRAMAS EN BLOQUE
4.2. b)
Dado el sistema simulado por el diagrama de bloques a continuación, determinar la ecuación
en diferencias que lo describe.
Podemos ver que el sistema tiene 2 retardos, por lo que la ecuación de diferencias deberá
ser de segundo orden.
Leyendo el diagrama de izquierda a derecha obtenemos:
1
y[n] = −5x[n] + D 3y[n] + D 2
x[n] − 3y[n]
Distribuyendo el operador D:
Espectro de Frecuencias
5.1. a)
2
2π
X
x(n) = ak ej 3
kn
k=0
2
1 X 2π
2πn
e−j kn
ak = 1 + cos 3
3
3 n=0
17
18 EJERCICIO 5. ESPECTRO DE FRECUENCIAS
Empezando por a0 :
2
1 X 2π
2πn
e−j (0)n
a0 = 1 + cos 3
3
3 n=0
2
1X 2πn
= 1 + cos 3
3 n=0
1
= 2 + 12 + 21
3
= a0 = 1
Ahora para a1 :
2
1 X 2π
2πn
e−j n
a1 = 1 + cos 3
3
3 n=0
Para este cálculo y otros similares se utilizará la relación de Euler ejθ = cos(θ) + j sin(θ).
1h √ √ i
a1 = 2 + − 4 − j 4 + − 4 + j 43
1 3 1
3
1
a1 = 2 − 12
3
a1 = 21
2π 4π
x(n) = 1 + 12 ej 3
n
+ 21 ej 3
n
5.2. b)
Obtener el espectro de frecuencias para N = 6
x[n] = 1 + cos πn
3
5
π
X
x(n) = ak ej 3 kn
k=0
5
1 X π
x(n)e−j 3 kn
6 n=0
5.2. B) 19
Empezando por a0 :
5
1X πn
a0 = 1 + cos
6 n=0 3
1 3 1 1 3
a0 = [2 + + + 0 + + ]
6 2 2 2 2
1
a0 = [6]
6
a0 = 1
Cálculo de a1
5
1X πn π
a1 = [1 + cos( ]e−j 3 n
6 n=0 3
√ √ √ √
1 3 3 3 1 3 1 3 3 3 3
a1 = [2 + (− − j )+( −j )+( +j ) + (− + j )]
6 4 4 4 4 4 4 4 4
1
a1 = [1]
6
1
a1 =
6
Cálculo de a2
5
1X πn 2π
a2 = [1 + cos( ]e−j 3 n
6 n=0 3
√ √ √ √
1 3 3 3 1 3 1 3 3 3 3
a2 = [2 + ( − j )+( +j )+( −j )+( +j )]
6 4 4 4 4 4 4 4 4
1
a2 = [4]
6
2
a2 =
3
Cálculo de a3
5
1X πn
a3 = [1 + cos( ]e−jπn
6 n=0 3
1 3 1 1 3
a3 = [2 − + + 0 + − ]
6 2 2 2 2
a3 = 0
Cálculo de a4
5
1X πn 4π
a4 = [1 + cos( ]e−j 3 n
6 n=0 3
√ √ √ √
1 3 3 3 1 3 1 3 3 3 3
a4 = [2 + ( − j )+( −j )+( +j )+( +j )]
6 4 4 4 4 4 4 4 4
2
a4 =
3
20 EJERCICIO 5. ESPECTRO DE FRECUENCIAS
Cálculo de a5
5
1X πn 5π
a5 = [1 + cos( ]e−j 3 n
6 n=0 3
√ √ √ √
1 3 3 3 1 3 1 3 3 3 3
a5 = [2 + (− + j )+( +j )+( −j ) + (− − j )]
6 4 4 4 4 4 4 4 4
1
a5 =
6
5
X 1 π 5π 2 2π 4π
x(n) = 1 + (ej 3 n + ej 3 n ) + (ej 3 n + ej 3 n )
k=0
6 3
Ejercicio 6
Sistema LTI
6.1. a)
Determinar Y (z) = Z{y[n]}, la solución completa, con las siguientes CI:
6.2. b)
Obtener la FdT H(z), respuesta al impulso con CI nulas, por medio de residuos.
Y (z)
Ya que H(z) = X(z)
1 z3
H(z) = = →
1 − 1,2z −1 − 0,13z −2 + 0,36z −3 z 3 (1 − 1,2z −1 − 0,13z −2 + 0,36z −3 )
z3 z3
H(z) = =
z 3 − 1,2z 2 − 0,13z + 0,36 (z − 0,5)(z + 0,9)(z + 0,8)
6.3. c)
Obtener h[n] = Z −1 {H(z)}. A partir de
I I
1 n−1 1 X
x[n] = X(z)z dz → h[n] = H(z)z n−1 dz = Res(H(z)z n−1 , zm )
2πj 2πj
z 3 z n−1 z2zn
= ≡ S(z)
(z − 0,5)(z + 0,9)(z + 0,8) (z − 0,5)(z + 0,9)(z + 0,8)
21
22 EJERCICIO 6. SISTEMA LTI
Res(S(z),-0.5)
0,25(−0,5)n 25
lı́m S(z)(z + 0,5) = = (−0,5)n
z→−0,5 (−1,4)(−1,3) 182
Res(S(z),0.9)
0,81(0,9)n 81
lı́m S(z)(z + 0,9) = = (0,9)n
z→0,9 (1,4)(0,1) 14
Res(S(z),0.8)
0,64(0,8)n −64
lı́m S(z)(z + 0,8) = = (0,8)n
z→0,8 (1,3)(−0,1) 13
6.4. d)
Obtener la respuesta al escalón unitario con CI nulas, por medio de fracciones parciales:
y[n] = Z −1 {Y (z)}
Queremos desarrollar en serie el siguiente cociente en la transformada Z:
z4
Y (z) = 4
z − 3,2z 3 + 3,77z 2 − 1,93z + 0,36
Comenzamos la división larga:
z4 ÷ z4 = 1
Multiplicamos: (z 4 − 3,2z 3 + 3,77z 2 − 1,93z + 0,36)(1) ⇒ R1 = 3,2z 3 − 3,77z 2 + 1,93z − 0,36
R1 /z 3 = 3,2
Multiplicamos y restamos: 3,2(z 4 − 3,2z 3 + 3,77z 2 − 1,93z + 0,36)
= 3,2z 4 − 10,24z 3 + 12,064z 2 − 6,176z + 1,152
R2 /z 2 = 6,47
⇒ R3 = 7,7319z − 9,24537
6.5. E) 23
R3 /z = 7,7319
⇒ R4 = 9,24537
6.5. e)
Obtener y[n] = µ[n] ∗ h[n], la respuesta al escalón unitario sin CI por medio de convolución
temporal.
9,24537 7,7319
,
z4 z3
Estos valores se emplean en los pasos sucesivos de la división continua para obtener los coe-
cientes de la expansión en potencias de 1/z .
6.6. f)
Vericar y[n] para la respuesta al escalón unitario por división directa hasta la 5° muestra
y reconstruir Y (z) = Z{y[n]}.
Dado:
z4
X(z) =
z 4 − 3,2z 3 + 3,77z 2 − 1,93z + 0,36
Evaluamos el límite:
z4
lı́m X(z) = lı́m
z→∞ z→∞ z 4 − 3,2z 3 + 3,77z 2 − 1,93z + 0,36
1
lı́m 3,2 3,77 1,93 0,36 =1
z→∞ 1− z
+ z2
− z3
+ z4
Por lo tanto:
x(0) = 1
24 EJERCICIO 6. SISTEMA LTI
6.7. g)
Vericar el TVI (régimen transitorio) y TVF (régimen de estado permanente) en ambos
dominios {z, n} para la respuesta al escalón unitario con CI nulas.
Aplicando a X(z):
z4
lı́m (1 − z −1 )
z→1 z 4 − 3,2z 3 + 3,77z 2 − 1,93z + 0,36
Multiplicamos por z/z para simplicar:
z 3 (z − 1)
lı́m
z→1 z 4 − 3,2z 3 + 3,77z 2 − 1,93z + 0,36
Numerador: (1)3 (1 − 1) = 0
Denominador: 1 − 3,2 + 3,77 − 1,93 + 0,36 = 0
Aplicamos la regla de L'Hôpital:
4z 3 − 3z 2
TVF = lı́m
z→1 4z 3 − 9,6z 2 + 7,54z − 1,93
Sustituyendo z = 1:
4−3 1
TVF = = ≈ 100
4 − 9,6 + 7,54 − 1,93 0,01
Por tanto:
x(∞) = 100
6.8. h)
Elaborar la síntesis del sistema por medio de diagrama de bloques con retardos unitarios.