Universidad Nacional de Ingenierı́a

Facultad de Ingenierı́a Mecánica

PRÁCTICA CALIFICADA N°3

Curso: Procesamiento Digital de Señales – MT417

Villafani Benites, Diego Alonso Código: 20192063K

Jorge Edinson Reátegui Llanos Código: 20127028J
Alex Enrique Saldivar Loayza Código: 20200006G
Valderrama Damian Rezo Lee Código: 20182527D
Lima, Perú
2025-I

Docente: JOO AGUAYO, ALCIDES GUILLERMO

Sección: A
Problema 1
Se nos da una función de transferencia discreta en el dominio Z:
z − 0.4
H(z) =
z(z 2 − 1.42z + 0.85)
Se pide analizar las caracterı́sticas del sistema y su comportamiento en
frecuencia.

Análisis de Polos y Ceros

Factorizamos el denominador para encontrar los polos:

p
1.42 ± 1.422 − 4(0.85)
z 2 − 1.42z + 0.85 = 0 ⇒ z = ⇒ z = 0.7 ± 0.6i
2
• Cero: z = 0.4

• Polos: z = 0, z = 0.7 + 0.6i, z = 0.7 − 0.6i

a) Propiedades del Sistema

• Estable: Todos los polos están dentro del cı́rculo unitario, por lo tanto
el sistema es estable.

• Causal: El número de polos es mayor al de ceros, y la función está

expresada como cociente de polinomios en z, por lo tanto es causal.

• No es FIR: El sistema tiene polos, por lo tanto no es de respuesta

finita al impulso (FIR).

• No tiene fase lineal: No hay simetrı́a entre polos y ceros respecto al

eje real.

• Coeficientes reales: Sı́, ya que los polos complejos están en pares

conjugados.

• No es todo polos .

2
Plano-Z: Diagrama de Polos y Ceros
Im(z)

Cı́rculo unitario
0.7 + 0.6i

Re(z)
0 0.4 (cero)

0.7 − 0.6i

b) Evaluación en z = ejπ/6

Figure 1: Distancias a Z.

3
Calculamos módulo y ángulos entre el punto de evaluación y los polos/ceros:

|z − 0.4| ≈ 0.6835, ∠(z − 0.4) ≈ 0.8206 rad

|z − (0.7 + 0.6i)| ≈ 0.1938, ∠ ≈ 0.5421 rad
|z − (0.7 − 0.6i)| ≈ 1.125, ∠ ≈ 1.4210 rad
La magnitud total:
0.6835
|H(ejπ/6 )| = ≈ 3.1701
1 · 0.1938 · 1.125
La fase total:

∠H(ejπ/6 ) = 0.8206 − 0.5236 − 0.5421 − 1.4210 ≈ −1.6661 rad

c) Evaluación en z = ejπ/4

Figure 2: Distancias a Z.

Análogamente:

4
 |z − 0.4| ≈ 0.8025, ∠(z − 0.4) ≈ 0.8961 rad
|z − (0.7 + 0.6i)| ≈ 0.1977, ∠ ≈ 1.1005 rad
|z − (0.7 − 0.6i)| ≈ 1.0726, ∠ ≈ 2.4897 rad
Entonces:
0.8025
|H(ejπ/4 )| = ≈ 3.7455
1 · 0.1977 · 1.0726
Y la fase:

∠H(ejπ/4 ) = 0.8961 − 0.7854 − 1.1005 − 2.4897 ≈ −3.4795 rad

Problema 2
a) Armando la ecuacion de diferencias a partir de los
datos brindados :
Sea:

H(z) = h(0) 1 + z̄ 3 + h(1) z̄ + z̄ 2

Evaluamos la respuesta en frecuencia:

H(ejω ) = h(0) 1 + e−j3ω + h(1) e−jω + e−j2ω

Agrupando términos y usando propiedades de exponenciales:

= e−j1.5ω h(0)(ej1.5ω + e−j1.5ω ) + h(1)(ej0.5ω + e−j0.5ω )

 

= e−j1.5ω [2h(0) cos(1.5ω) + 2h(1) cos(0.5ω)]

5
 Figure 3: Diagrama de polos.

b) Análisis de polos y ceros:

(z + 1)(z 2 + 1)
H(z) = 1 + z + z −1 + z −2 + z −3 =
z3
• Ceros: z1 = −1, z2 , z3 = ±j
• Los ceros están en el Cı́rculo Unitario (C.R.U.).
• Conclusión: Es un sistema de fase no lineal.

c) Condición para fase lineal:

H(z) = (1 + z −1 )(1 − z −1 + z −2 )
Desarrollamos:

H(z) = h(0)(1 + z −1 ) + h(1)(−z −1 + z −2 ) = α(1 + z −1 ) + βz −1 + αz −2

Donde:
α = h(0), β = h(1) − h(0)

6
 ω 
jω −j1.5ω
H(e ) = 2e cos (β + 2 cos(ω))
2
La magnitud es real solo si:
ω 
cos (β + 2 cos(ω)) > 0
2
Esto se cumple si:

2 cos(ω) + β > 0 ⇒ h(1) − h(0) > −2h(0) ⇒ h(1) < 3h(0)

Condición: h(1) < 3h(0) ⇒ fase lineal

d) Análisis de fase lineal mediante raı́ces

H(z) = (α + βz −1 + αz −2 )(1 + z −1 )

β
z1 = −1, z2 + z+1=0
α
p
−β ± (β)2 − 4α2
z2,3 =
2α
Nota: Para fase lineal, las raı́ces deben estar en el Cı́rculo Unitario
(C.R.U.).

β h(0) − h(1)
z + z −1 = − ⇒ Usamos:
α h(0)
Si z0 = 0, entonces z10 → ∞, no se cumple.
Conclusión: Bajo esta condición, el sistema sı́ cumple fase lineal.

PROBLEMA 3 (5P)
Diseñe un filtro digital de fase lineal tipo rechaza-banda de 5 coeficientes
utilizando una ventana de Blackman. El filtro debe tener:

• Una frecuencia de corte inferior de 40 Hz

• Una frecuencia de corte superior de 80 Hz
• Una frecuencia de muestreo de 360 muestras por segundo

7
Se pide determinar:
a) Los coeficientes del filtro h[n] (2p)

b) La respuesta en frecuencia del sistema H(ejω ) (1p)

c) La ecuación en diferencias del sistema (1p)

Datos proporcionados:
La ventana de Blackman está definida como:

   
0.42 − 0.5 cos 2πn + 0.08 cos 4πn , 0 ≤ n ≤ N


w[n] = N N
0, en otro caso


donde N es el orden del filtro.

SOLUCIÓN:
Datos del Problema
• Tipo de filtro: Rechaza-banda (band-stop)

• Longitud del filtro: 5 coeficientes (N = 4, ya que MATLAB cuenta

desde 0)

• Frecuencia de corte inferior: 40 Hz

• Frecuencia de corte superior: 80 Hz

• Frecuencia de muestreo: 360 Hz

• Ventana utilizada: Blackman

8
 Determinamos las frecuencias digitales nor-
malizadas
Primero convertimos las frecuencias analógicas a frecuencias digitales nor-
malizadas (ω):
2πf
ω=
fs
Donde fs es la frecuencia de muestreo (360 Hz).
• Frecuencia de corte inferior (ω1 ):
2π × 40 2π
ω1 = = ≈ 0.6981 rad/muestra
360 9
• Frecuencia de corte superior (ω2 ):
2π × 80 4π
ω2 = = ≈ 1.3963 rad/muestra
360 9

Obtenemos la respuesta al impulso ideal del

filtro rechaza-banda
La respuesta al impulso ideal para un filtro rechaza-banda es:
ω2 − ω1

1 − , n=0


π
hideal [n] =
 sin(ω2 n) sin(ω1 n)
δ[n] −
 + , n ̸= 0
πn πn

Cálculo para n = 0:
4π 2π
9
− 9
hideal [0] = 1 −
π
2π
9
=1−
π
2
=1−
9
≈ 0.7778

9
Cálculo para n ̸= 0:
4πn 2πn



sin 9
sin 9
hideal [n] = − +
πn πn

Cálculos detallados:

Table 1: Cálculo de coeficientes ideales

n Expresión Valor
−0.9848 −0.6428
≈− +
− 8π − 4π



sin 9
sin 9
−2π −2π
−2 − + ≈ −0.1567 + 0.1023
−2π −2π
≈ −0.0544
−0.6428 −0.3420
≈− +
sin − 4π 2π



9
sin − 9
−π −π
−1 − + ≈ −0.2046 + 0.1089
−π −π
≈ −0.0957
0.6428 0.3420
4π

2π

≈− +
sin 9
sin 9 π π
1 − + ≈ −0.2046 + 0.1089
π π
≈ −0.0957
0.9848 0.6428
≈− +
sin 8π sin 4π



9 9 2π 2π
2 − + ≈ −0.1567 + 0.1023
2π 2π
≈ −0.0544

Resumen de coeficientes:
 
hideal = −0.0544 −0.0957 0.7778 −0.0957 −0.0544 para n = −2 a 2

Aplicamos la ventana de Blackman

La ventana de Blackman para N = 4 (5 puntos) es: La ventana de Blackman
para N = 4 (5 puntos) es:

10
    
2πn 4πn
w[n] = 0.42 − 0.5 cos + 0.08 cos para n = 0 a 4
N N

Cálculo para cada valor de n:

Table 2: Cálculo de los valores de la ventana de Blackman

n Expresión Valor
0 0.42 − 0.5 cos(0) + 0.08 cos(0) = 0.42 − 0.5 + 0.08 0
1 0.42 − 0.5 cos(π/2) + 0.08 cos(π) = 0.42 − 0 + .08(−1) 0.34
2 0.42 − 0.5 cos(π) + 0.08 cos(2π) = 0.42 − (−0.5) + 0.08 1.0
3 0.42 − 0.5 cos(3π/2) + 0.08 cos(3π) = 0.42 − 0 + 0.08(−1) 0.34
4 0.42 − 0.5 cos(2π) + 0.08 cos(4π) = 0.42 − 0.5 + 0.08 0

Aplicación centrada de la ventana:

Como nuestro filtro es simétrico y va de n = −2 a n = 2, aplicamos la
ventana centrada:

wwindow = [w[4], w[3], w[2], w[1], w[0]] = [0, 0.34, 1.0, 0.34, 0]

Multiplicamos la respuesta ideal por la ven-

tana
h[n] = hideal [n] × wwindow [n]

h = [−0.0544 × 0, −0.0957 × 0.34, 0.7778 × 1.0, −0.0957 × 0.34, −0.0544 × 0]

= [0, −0.037, 0.7778, −0.037, 0]

11
Respuestas
a) Coeficientes del filtro h[n] (2p)
Los coeficientes del filtro son:

h[0] = −0.037
h[1] = 0.7778
h[2] = −0.037
h[3] = 0
h[4] = 0

O en orden causal (n = 0 a 4):

h[n] = [0, −0.037, 0.7778, −0.037, 0]

b) Respuesta en frecuencia H(ejω ) (1p)

La respuesta en frecuencia se obtiene mediante la DTFT de h[n]:
4
X
H(e ) =jω
h[n]e−jωn
n=0
−jω
jω
H(e ) = 0 − 0.037e + 0.7778e−j2ω − 0.037e−j3ω + 0

c) Ecuación en diferencias del sistema (1p)

La ecuación en diferencias se obtiene directamente de los coeficientes:

y[n] = 0 · x[n] − 0.037x[n − 1] + 0.7778x[n − 2] − 0.037x[n − 3] + 0 · x[n − 4]

= −0.037x[n − 1] + 0.7778x[n − 2] − 0.037x[n − 3]

PROBLEMA 4
Diseñar un filtro Butterworth (IIR), digital que cumpla con las siguientes
especificaciones:

12
• Tipo de filtro: Rechaza-banda

• Frecuencia de muestreo: fs = 1000 Hz

• Frecuencias de paso:

– fp1 = 100 Hz (inicio de la banda de paso)

– fp2 = 300 Hz (fin de la banda de paso)

• Frecuencias de rechazo:

– fr1 = 150 Hz (inicio de la banda rechazada)

– fr2 = 250 Hz (fin de la banda rechazada)

• Atenuación máxima en banda de paso: Ap = 1 dB

• Atenuación mı́nima en banda de rechazo: As = 40 dB

Se pide:
Dar su respuesta en la representación de forma directa II o en grafo de
flujo de señales. (4P)

Solución:
Diseño de filtro Butterworth IIR rechaza-banda (Notch)
Especificaciones del filtro
• Frecuencia de muestreo: fs = 1000 Hz

• Frecuencias de paso: fp1 = 100 Hz, fp2 = 300 Hz

• Frecuencias de rechazo: fr1 = 150 Hz, fr2 = 250 Hz

• Atenuación máxima en banda de paso: Ap = 1 dB

• Atenuación mı́nima en banda de rechazo: As = 40 dB

13
1. Conversión de frecuencias a frecuencias normalizadas
fs
La frecuencia de Nyquist es fN = = 500 Hz.
2
Se normalizan todas las frecuencias respecto a fN :
100
fp1 = 100 Hz =⇒ ωp1 = = 0.20
500
150
fr1 = 150 Hz =⇒ ωr1 = = 0.30
500
250
fr2 = 250 Hz =⇒ ωr2 = = 0.50
500
300
fp2 = 300 Hz =⇒ ωp2 = = 0.60
500
Por lo tanto:

Banda de paso: [0, 0.20] ∪ [0.60, 1.00] Banda de rechazo: [0.30, 0.50]

2. Cálculo de parámetros de atenuación

• Rizado en banda de paso:
p
ε = 10Ap /10 − 1
p √ √
ε = 101/10 − 1 = 1.2589 − 1 = 0.2589 ≈ 0.508

• Parámetro de banda de rechazo:

δ = 10−As /20 = 10−40/20 = 10−2 = 0.01

3. Cálculo del orden mı́nimo del filtro (N)

La fórmula para el orden de un filtro Butterworth rechaza-banda es:
 As /10 
10 −1
log10
10Ap /10 − 1
N≥
2 log10 (Ωs /Ωp )

Donde:

14
 - Usamos los dos pares de transiciones más estrechas:

Transición izquierda: ωp1 = 0.20 → ωr1 = 0.30

Transición derecha: ωr2 = 0.50 → ωp2 = 0.60

Calculamos ambos órdenes:

a) Orden en la transición izquierda

 40/10 
10 −1
log10
101/10 − 1
N1 =
2 log10 (0.30/0.20)

1040/10 − 1 = 10000 − 1 = 9999

101/10 − 1 = 1.2589 − 1 = 0.2589
9999
= 38616.5
0.2589
log10 (38616.5) ≈ 4.587
2 log10 (0.30/0.20) = 2 log10 (1.5) = 2 × 0.1761 = 0.3522
4.587
N1 = ≈ 13.02
0.3522

b) Orden en la transición derecha

 
9999
log10
0.2589
N2 =
2 log10 (0.50/0.40)

0.50/0.40 = 1.25 =⇒ log10 (1.25) = 0.0969

2 × 0.0969 = 0.1938
4.587
N2 = ≈ 23.67
0.1938

c) Orden final

N = ⌈max(N1 , N2 )⌉ = ⌈23.67⌉ = 24

15
4. Resumen del diseño
• Orden mı́nimo del filtro: N = 24

• Frecuencias normalizadas:

– Banda de paso: [0, 0.20] y [0.60, 1.00]

– Banda de rechazo: [0.30, 0.50]

• Atenuación máxima en banda de paso: Ap = 1 dB

• Atenuación mı́nima en banda de rechazo: As = 40 dB

Por lo tanto, el filtro requerido es un Butterworth digital rechaza-banda

de orden 24 con las frecuencias y atenuaciones especificadas.

En software, estos parámetros se utilizan directamente en la función de diseño

(butter, scipy.signal.butter, MATLAB butter), por ejemplo:

[b, a] = butter(24, [0.30, 0.50], ’stop’)

PROBLEMA 5 (4P)
Se tiene un filtro analógico Butterworth de segundo orden, tipo pasa bajas,
con una frecuencia de media potencia (half-power ) de 1 rad/s. Este filtro se
convierte en un filtro digital H(z), mediante la transformación bilineal, con
una frecuencia de muestreo de S = 1 Hz.

Se pide:
a) ¿Cuál es la función de transferencia H(s) del filtro analógico? (1p)

b) ¿Cuál es la función de transferencia H(z) del filtro digital resultante,

usando la transformación bilineal? (1p)

c) ¿Son idénticas las ganancias en DC (frecuencia cero) de H(s) y H(z)?

¿Deberı́an serlo? Justifique su respuesta. (1p)

d) ¿Son idénticas las ganancias en las frecuencias de media potencia (half-

power ) de H(s) y H(z)? ¿Deberı́an serlo? Justifique. (1p)

16
Solución:
Datos del Problema
Filtro analógico Butterworth de segundo orden, tipo pasa-bajas
Frecuencia de media potencia (half-power): 1 rad/s
Transformación bilineal con frecuencia de muestreo S = 1 Hz

a) Función de transferencia H(s) del filtro analógico

(1p)
Para un filtro Butterworth pasa-bajas de segundo orden con frecuencia de
corte ωc = 1 rad/s, la función de transferencia normalizada es:
1
H(s) = √
s2 + 2s + 1
Esta es la forma estándar de un filtro Butterworth de segundo orden con
ωc = 1 rad/s.

b) Función de transferencia H(z) usando transformación

bilineal (1p)
La transformación bilineal está dada por:
2 z−1
s= ·
T z+1
Donde T = 1/S = 1/1 = 1 (ya que S = 1 Hz)
Por lo tanto:
z−1
s=2·
z+1
Sustituyendo en H(s):
1
H(z) = √
z−1 2 z−1



2· z+1
+ 2· 2· z+1
+1
Desarrollando:
Término cuadrático:
2
4(z − 1)2

z−1
2· =
z+1 (z + 1)2

17
 Término lineal: √
√ z−1 2 2(z − 1)
2·2· =
z+1 z+1
Reescribimos H(z):
1
H(z) = √
4(z−1)2 2 2(z−1)
(z+1)2
+ z+1
+1

Multiplicamos numerador y denominador por (z + 1)2 para eliminar de-

nominadores:
(z + 1)2
H(z) = √
4(z − 1)2 + 2 2(z − 1)(z + 1) + (z + 1)2

Desarrollamos cada término:

4(z − 1)2 = 4(z 2 − 2z + 1) = 4z 2 − 8z + 4

√ √ √ √
2 2(z − 1)(z + 1) = 2 2(z 2 − 1) = 2 2z 2 − 2 2
(z + 1)2 = z 2 + 2z + 1

Sumando todos los términos del denominador:

√ √
(4z 2 − 8z + 4) + (2 2z 2 − 2 2) + (z 2 + 2z + 1)
√ √
= (4 + 2 2 + 1)z 2 + (−8 + 2)z + (4 − 2 2 + 1)
√ √
= (5 + 2 2)z 2 − 6z + (5 − 2 2)

Por lo tanto:
z 2 + 2z + 1
H(z) = √ √
(5 + 2 2)z 2 − 6z + (5 − 2 2)

(z + 1)2
H(z) ≈
7.8284z 2 − 6z + 2.1716

18
c) Comparación de ganancias en DC (frecuencia cero)
(1p)
Para el filtro analógico H(s):
En DC (s = 0):
1
H(0) = =1
0+0+1
Para el filtro digital H(z):
En DC (z = 1):
1+2+1 4 4
H(1) = √ √ = = =1
(5 + 2 2) − 6 + (5 − 2 2) 10 − 6 4
Respuesta:
Las ganancias en DC son idénticas (ambas igual a 1). Esto deberı́a ser ası́
porque la transformación bilineal preserva la respuesta en frecuencia en DC.

d) Comparación de ganancias en las frecuencias de me-

dia potencia (1p)
Para el filtro analógico H(s):
En la frecuencia de media potencia (s = jωc = j1):
1 1 1
|H(j1)| = √ =q √ = √ ≈ 0.7071 (-3 dB)
| − 1 + j 2 + 1| 02 + ( 2)2 2

Para el filtro digital H(z):

La frecuencia de media potencia digital θ = 1 ⇒ z = ej
Evaluamos H(z) en z = ej :
ej + 2ej + 1
H(ej ) = √ √
(5 + 2 2)ej − 6ej + (5 − 2 2)
Calculando la magnitud:

|H(ej )| ≈ 0.6421

Respuesta:
La transformación bilineal sin ajuste distorsiona la respuesta en frecuencia,
lo que explica la discrepancia inicial (0.707 vs 0.0421).

19