lOMoARcPSD|17850602

Tarea 3 Juan Contreras (1)

Catedra unadista (Universidad Nacional Abierta y a Distancia)

Studocu is not sponsored or endorsed by any college or university

Downloaded by ronald garzon ([email protected])
 lOMoARcPSD|17850602

Tarea 3 – Planificar métodos y herramientas para el diseño de filtros digitales

Contreras, Juan Erasmo

Universidad Nacional Abierta y a Distancia
Escuela de Ciencias Básicas, Tecnología e Ingeniería
Curso de Procesamiento Digital de Señales
Bucaramanga - Colombia
[email protected]

Abstract: In this document the different points proposed for step

3 are developed, where the methods and tools for the design of FIR
and IIR digital filters are planned, mainly a difference equation
that is stipulated in the work is chosen, then a block diagram of
the difference equation in draw.io, the Z transform of the
difference equation is developed, the transfer function of the
system H (Z) is found, the frequency response of the system is
performed, and the magnitude of the frequency response, the
Euler identity is applied, the represented function of the phase
response of the system is found. The simulation is also developed
in MATLAB (Simulink), for the impulse response diagrams of the
system, diagram of poles and zeros, Bode diagram and finally the
theoretical contributions are provided.

Keywords: FIR, IIR, Z transform, frequency response, impulse

response, Bode diagram.

INTRODUCCIÓN Figura 1. Diagrama de bloques IIR

En este documento se encuentra los principios de la • Cada estudiante realizará la transformada Z de la
transformada Z como elemento principal para evaluar la ecuación de diferencias. Esta debe realizarse en el editor de
estabilidad de los sistemas frente a diferentes estímulos y ecuaciones de Word. No se aceptan pantallazos.
determinaciones, se verán las respuestas ante entrada impulso,
su diagrama de Bode y su respuesta en fase. 𝒀(𝒁) = 𝒃𝟎 𝑿(𝒁) + 𝒃𝟏 𝑿(𝒁)𝒛−𝟏 − 𝒂𝟏 𝒀(𝒁)𝒛−𝟏 − 𝒂𝟐 𝒀(𝒁)𝒛−𝟐

Actividades a desarrollar
• Una vez se tenga la transformada Z de la ecuación de
Cada estudiante escogerá una (1) ecuación de diferencia de diferencia, cada estudiante hallará la función de transferencia
las expuestas a continuación, luego reportará en el foro su del sistema H(Z). Esto también se realizará con el editor de
decisión esto con el fin de que cada estudiante tenga una ecuaciones de Word. Recordar que la función de transferencia
ecuación diferente es:
1- Ecuaciones diferencia: 𝒚(𝒛)
𝒉(𝒛) =
𝒙(𝒛)
Ejercicio 3
Separando variables se tiene;
𝒚[𝒏] = 𝒃𝟎 𝒙[𝒏] + 𝒃𝟏 𝒙[𝒏 − 𝟏] − 𝒂𝟏 𝒚[𝒏 − 𝟏] − 𝒂𝟐 𝒚[𝒏 − 𝟐]

• Cada estudiante realizará el diagrama de bloques de su 𝑌(𝑍) + 𝑎1 𝑌(𝑍)𝑧 −1 + 𝑎2 𝑌(𝑍)𝑧 −2 = 𝑏0 𝑋(𝑍) + 𝑏1 𝑋(𝑍)𝑧 −1
ecuación de diferencia en la página de internet:
Aplicando factor común de Y(Z) y X(Z), se obtiene;

𝑌(𝑍)(1 + 𝑎1 𝑧 −1 + 𝑎2 𝑧 −1 ) = 𝑋(𝑍)(𝑏0 + 𝑏1 𝑧 −1 )

Por lo tanto, la función de transferencia H(Z), será;

Downloaded by ronald garzon ([email protected])

 lOMoARcPSD|17850602

𝑌(𝑍) 𝒃𝟎 + 𝒃𝟏 𝒛−𝟏 La magnitud de la respuesta en frecuencia es;

𝐻(𝑍) = =
𝑋(𝑍) 𝟏 + 𝒂𝟏 𝒛−𝟏 + 𝒂𝟐 𝒛−𝟐
𝐴𝐶 + 𝐵𝐷 2 𝐴𝐷 − 𝐵𝐶 2
|𝐻(𝜔)| = √( 2 ) + ( )
• Una vez se tenga la función de transferencia, se hallará 𝐶 + 𝐷2 𝐶 2 + 𝐷2
la respuesta en frecuencia del sistema, remplazando: 𝐴2 𝐶 2 + 2𝐴𝐵𝐶𝐷 + 𝐵2 𝐷2 + 𝐴2 𝐷2 − 2𝐴𝐵𝐶𝐷 + 𝐵2 𝐶 2
= √
𝒛 = 𝒆𝒋𝒘 (𝐶 2 + 𝐷2 )2

La respuesta en frecuencia de la función de transferencia es;

𝐴2 𝐶 2 + 𝐵2 𝐷2 + 𝐴2 𝐷2 + 𝐵2 𝐶 2
|𝐻(𝜔)| = √
𝒃𝟎 + 𝒃𝟏 𝒆−𝒋𝝎 (𝐶 2 + 𝐷2 )2
𝑯(𝝎) =
𝟏 + 𝒂𝟏 𝒆−𝒋𝝎 + 𝒂𝟏 𝒆−𝟐𝒋𝝎
𝐴2 (𝐶 2 + 𝐷2 ) + 𝐵2 (𝐷2 + 𝐶 2 )
|𝐻(𝜔)| = √
• Una vez se cuente con la respuesta en frecuencia del (𝐶 2 + 𝐷2 )2
sistema, se hallará la magnitud de la respuesta en frecuencia,
para ello se aplicará la identidad de Euler, que según el caso se
podría utilizar cualquiera de las siguientes ecuaciones: (𝐴2 + 𝐵2 )(𝐷2 + 𝐶 2 )
|𝐻(𝜔)| = √
(𝐶 2 + 𝐷2 )2
𝒆𝒋𝒘 = 𝐜𝐨𝐬(𝒘) + 𝒋𝒔𝒊𝒏(𝒘)

𝒆−𝒋𝒘 = 𝐜𝐨𝐬(𝒘) − 𝒋𝒔𝒊𝒏(𝒘) 𝐴2 + 𝐵2

|𝐻(𝜔)| = √
𝐶 2 + 𝐷2
Para hallar la función de magnitud, recordar utilizar la siguiente
ecuación:

|𝒂 + 𝒃𝒋| = √𝒂𝟐 + 𝒃𝟐
• Se hallará la función que represente la respuesta en
Donde a y b son los coeficientes del número imaginario (𝑎 +
Fase del sistema, recordar utilizar la siguiente ecuación:
𝑏𝑗)
𝒃
La respuesta en frecuencia aplicando la identidad de Euler es; 𝜽(𝒂 + 𝒃𝒋) = 𝐚𝐫𝐜𝐭𝐚𝐧 ( )
𝑏0 + 𝑏1 (cos 𝜔 − 𝑗 sin 𝜔) 𝒂
𝐻(𝜔) =
1 + 𝑎1 (cos 𝜔 − 𝑗 sin 𝜔) + 𝑎2 (cos 2𝜔 − 𝑗 sin 2𝜔) Donde b son los coeficientes de numero (𝒂 + 𝒃𝒋)
𝑏0 + 𝑏1 cos 𝜔 − 𝑗𝑏1 sin 𝜔
=
1 + 𝑎1 cos 𝜔 − 𝑗𝑎1 sin 𝜔 + 𝑎2 cos 2𝜔 − 𝑗𝑎2 sin 2𝜔 Para la fase se tiene;
Para simplificar un poco más la ecuación se emplearán las
siguientes sustituciones; 𝐴𝐷 − 𝐵𝐶
2 2 𝐴𝐷 − 𝐵𝐶
𝜃 = tan −1
( 𝐶 + 𝐷 ) = 𝐭𝐚𝐧−𝟏 ( )
𝐴 = 𝑏0 , , 𝐵 = 𝑏1 sin 𝜔 , 𝐶 = 1 + 𝑎1 cos 𝜔 + 𝑎2 cos 2𝜔 𝑦 𝐷 𝐴𝐶 + 𝐵𝐷 𝐴𝐶 + 𝐵𝐷
= 𝑎1 sin 𝜔 + 𝑎2 sin 2𝜔 𝐶 2 + 𝐷2
Por lo tanto, se tiene;
• Grupal
𝐴 − 𝑗𝐵
𝐻(𝜔) =
𝐶 − 𝑗𝐷 𝒚[𝒏] = 𝒃𝟎 𝒙[𝒏] + 𝒃𝟏 𝒙[𝒏 − 𝟏] + 𝒃𝟐 𝒙[𝒏 − 𝟐] + 𝒃𝟑 𝒙[𝒏 − 𝟑]
+ 𝒃𝟒 𝒙[𝒏 − 𝟒] + 𝒃𝟓 𝒙[𝒏 − 𝟓]
Racionalizando el denominador se obtiene;
Diagrama de bloques
𝐴 − 𝑗𝐵 𝐶 + 𝑗𝐷 (𝐴 − 𝑗𝐵)(𝐶 + 𝑗𝐷)
𝐻(𝜔) = ∗ =
𝐶 − 𝑗𝐷 𝐶 + 𝑗𝐷 𝐶 2 + 𝐷2
𝐴𝐶 + 𝑗𝐴𝐷 − 𝑗𝐵𝐶 + 𝐵𝐷
=
𝐶 2 + 𝐷2
𝐴𝐶 + 𝐵𝐷 𝐴𝐷 − 𝐵𝐶
𝐻(𝜔 ) = +𝑗 2
𝐶 2 + 𝐷2 𝐶 + 𝐷2

Downloaded by ronald garzon ([email protected])

 lOMoARcPSD|17850602

𝐻(𝜔) = 𝑏0 + 𝑏1 (𝐜𝐨𝐬(𝝎) − 𝒋𝒔𝒊𝒏(𝝎))

+ 𝑏2 (𝐜𝐨𝐬(𝟐𝝎) − 𝒋𝒔𝒊𝒏(𝟐𝝎))
+ 𝑏𝟑 (𝐜𝐨𝐬(𝟑𝝎) − 𝒋𝒔𝒊𝒏(𝟑𝝎))
+ 𝑏𝟒 (𝐜𝐨𝐬(𝟒𝝎) − 𝒋𝒔𝒊𝒏(𝟒𝝎))
+ 𝑏𝟓 (𝐜𝐨𝐬(𝟓𝝎) − 𝒋𝒔𝒊𝒏(𝟓𝝎))

𝐻(𝜔) = 𝑏0 + 𝑏1 𝐜𝐨𝐬(𝝎) + 𝑏2 𝐜𝐨𝐬(𝟐𝝎) + 𝑏𝟑 𝐜𝐨𝐬(𝟑𝝎)

+ 𝑏𝟒 𝐜𝐨𝐬(𝟒𝝎) + 𝑏𝟓 𝐜𝐨𝐬(𝟓𝝎) − 𝑗𝑏1 𝒔𝒊𝒏(𝝎)
− 𝑗𝑏2 𝒔𝒊𝒏(𝟐𝝎) − 𝑗𝑏𝟑 𝒔𝒊𝒏(𝟑𝝎)
− 𝑗𝑏𝟒 𝒔𝒊𝒏(𝟒𝝎) − 𝑗𝑏𝟓 𝒔𝒊𝒏(𝟓𝝎)

𝐻(𝜔) = (𝑏0 + 𝑏1 𝐜𝐨𝐬(𝝎) + 𝑏2 𝐜𝐨𝐬(𝟐𝝎) + 𝑏𝟑 𝐜𝐨𝐬(𝟑𝝎)

+ 𝑏𝟒 𝐜𝐨𝐬(𝟒𝝎) + 𝑏𝟓 𝐜𝐨𝐬(𝟓𝝎))
− 𝑗(𝑏1 𝒔𝒊𝒏(𝝎) + 𝑏2 𝒔𝒊𝒏(𝟐𝝎) + 𝑏𝟑 𝒔𝒊𝒏(𝟑𝝎)
+ 𝑏𝟒 𝒔𝒊𝒏(𝟒𝝎) + 𝑏𝟓 𝒔𝒊𝒏(𝟓𝝎))

Para simplificar un poco más la ecuación se emplearán las

siguientes sustituciones;

𝐴 = 𝒃𝟎 + 𝒃𝟏 cos(𝜔) + 𝒃𝟐 cos(2𝜔) + 𝒃𝟑 cos(3𝜔)

+ 𝒃𝟒 cos(4𝜔) + 𝒃𝟓 cos(5𝜔) , 𝐵
= 𝒃𝟏 𝑠𝑖𝑛(𝜔) + 𝒃𝟐 𝑠𝑖𝑛(2𝜔) + 𝒃𝟑 𝑠𝑖𝑛(3𝜔)
+ 𝒃𝟒 𝑠𝑖𝑛(4𝜔) + 𝒃𝟓 𝑠𝑖𝑛(5𝜔)

Por lo tanto, se tiene;

𝐻(𝜔) = 𝐴 − 𝑗𝐵
La magnitud de la respuesta en frecuencia es;
Figura 2. Diagrama de bloques FIR
|𝑯(𝝎)| = √𝑨𝟐 + 𝑩𝟐

Transformada Z.
Fase;
𝒀(𝒁) = 𝒃𝟎 𝑿(𝒁) + 𝒃𝟏 𝑿(𝒁)𝒛−𝟏 + 𝒃𝟐 𝑿(𝒁)𝒛−𝟐 + 𝒃𝟑 𝑿(𝒁)𝒛−𝟑
+ 𝒃𝟒 𝑿(𝒁)𝒛−𝟒 + 𝒃𝟓 𝑿(𝒁)𝒛−𝟓 𝑩
𝜽 = 𝐭𝐚𝐧−𝟏 ( )
𝑨
Función de transferencia
• Función de transferencia en el Command Windows
Factorizando se obtiene;
IIR
−1 −2 −3 −4
𝑌(𝑍) = 𝑋(𝑍)(𝑏0 + 𝑏1 𝑧 + 𝑏2 𝑧 + 𝑏3 𝑧 + 𝑏4 𝑧
+ 𝑏5 𝑧 −5 ) 𝑌(𝑍) 𝑏0 𝑧 2 + 𝑏1 𝑧
𝐻(𝑍) = = 2
𝑋(𝑍) 𝑧 + 𝑎1 𝑧 + 𝑎2
Por lo tanto;
FIR
𝒀(𝒁)
𝑯(𝒁) = = 𝒃𝟎 + 𝒃𝟏 𝒛−𝟏 + 𝒃𝟐 𝒛−𝟐 + 𝒃𝟑 𝒛−𝟑 + 𝒃𝟒 𝒛−𝟒 𝒀(𝒁) 𝑏0 𝑧 5 + 𝑏1 𝑧 4 + 𝑏2 𝑧 3 + 𝑏3 𝑧 2 + 𝑏4 𝑧 + 𝑏5
𝑿(𝒁)
𝑯(𝒁) = =
+ 𝒃𝟓 𝒛−𝟓 𝑿(𝒁) 𝑧5

Respuesta en frecuencia.

𝑯(𝝎) = 𝒃𝟎 + 𝒃𝟏 𝒆−𝒋𝝎 + 𝒃𝟐 𝒆−𝟐𝒋𝝎 + 𝒃𝟑 𝒆−𝟑𝒋𝝎 + 𝒃𝟒 𝒆−𝟒𝒋𝝎

+ 𝒃𝟓 𝒆−𝟓𝒋𝝎

Magnitud de la respuesta en frecuencia;

Downloaded by ronald garzon ([email protected])

 lOMoARcPSD|17850602

Figura 3. Respuesta al impulso sistema IIR

Figura 5. Diagrama de polos y ceros IIR

En la figura 5 se puede apreciar que esta función de trasferencia

posee 2 polos y 2 ceros, los polos se encuentran dentro del
semiplano izquierdo por lo que se puede decir que es un sistema
estable.

Figura 4. Respuesta al impulso sistema FIR

• Diagrama de polos y ceros

Figura 6. Diagrama de polos y ceros FIR

En la figura 6 se puede ver que esta función de transferencia

posee 5 ceros y se encuentran dentro del círculo unitario

• Diagrama de bode

Downloaded by ronald garzon ([email protected])

 lOMoARcPSD|17850602

Figura 7. Diagrama de Bode IIR

En la figura 7 el filtro IIR tiene un compartimiento de filtro pasa

bajas en donde la frecuencia de corte es aproximadamente 2 Figura 9. Código empleado en el Command Windows de
rad/s, Matlab.

Marco Teórico

1-) ¿Qué es la transformada de Z?

La transformada es una herramienta matemática con la cual

podemos trabajar de una forma más sencilla en los sistemas
discretos básicamente se asigna una frecuencia discreta {𝒇(𝒌)}
[1].
∞

𝑭(𝒛) = ∑ 𝒇(𝒌) ∗ 𝒛−𝒌

𝒌=𝟎

2-) ¿Qué representa Z en una función?

Se puede identificar como un punto de relación en la función

que nos permite tener un punto de referencia con relación a
dicha ecuación [1].

Figura 8. Diagrama de Bode FIR 𝒛 = 𝒆𝒋𝒘

El filtro FIR tiene un comportamiento de rechaza banda. 3-) ¿cuál es la diferencia entre la transformada Z bilateral y la
unilateral?
- Código utilizado en el Command Window
La TZ bilateral de una señal definida en el dominio del tiempo
discreto x[n] es una función X(z) que se define;

Downloaded by ronald garzon ([email protected])

 lOMoARcPSD|17850602

∞
• Se pudo comprender el rol que cumple la transformada
𝑋(𝑧) = 𝑍{𝑥[𝑛]} = ∑ 𝑥[𝑛]𝑧 −𝑛 Z en tiempo discreto que es similar al de la
𝑛=−∞ transformada de Laplace en sistemas en tiempo
De forma alternativa, en los casos en que x[n] está definida continuo.
únicamente para n ≥ 0, la transformada Z unilateral se define
como; • Los filtros son aplicados en la separación de señales
∞
con respecto al ruido e interferencias de otros sistemas
𝑋 + (𝑧) = 𝑍 + {𝑥[𝑛]} = ∑ 𝑥[𝑛]𝑧 −𝑛 y modificación de señales para moldear espectros y
𝑛=0 formas de ondas
[2].

4-) ¿cómo se calculan los polos y ceros de una función en

términos de Z? REFERENCIAS
Se puede identificar que los diagramas de polos y ceros la
podemos interpretar en búsqueda de obtener la información
sobre la frecuencia en un estado estable, podemos ver el [1] I. J. Sacerdoti, «Transformada Z», Dep. Matemática Fac.
comportamiento del transitorio y la estabilidad del sistema. De Ing. Univ. B. Aires 2003, p. 46, 2003.
igual manera vemos un enfoque similar en el comportamiento [2] A. Zozaya, «Procesamiento digital de la señal
sistema en tiempo discreto [3]. Transformada z», Valencia - Venezuela, feb. 2007, [En
línea]. Disponible en:
𝒚(𝒛) (𝒛 − 𝒛𝟏 ). (𝒛 − 𝒛𝟐 ) … (𝒛 − 𝒛𝒎 ) https://d1wqtxts1xzle7.cloudfront.net/43750982/Procesa
𝒉(𝒛) = = 𝒛.
𝒙(𝒛) (𝒛 − 𝒑𝟏 ). (𝒛 − 𝒑𝟐 ) … (𝒛 − 𝒑𝒏 ) miento_digital_de_la_seal_Transfo20160315-20629-
14vga41.pdf?1458055295=&response-content-
5-) ¿qué es la respuesta en frecuencia de un sistema digital? disposition=inline%3B+filename%3DProcesamiento_dig
ital_de_la_senal_Transf.pdf&Expires=1604635974&Sign
La respuesta en frecuencia o respuesta frecuencial es un ature=ZDA7CMS3m~4InoObNTagmsqt6xG34nrKAd7K
parámetro que indica la variación de la salida a un estímulo de cImuSQcP20TKXD1msxkVJd8O6tYdHzycJLZsZs-
entrada con respecto a la frecuencia en el dominio discreto [4]. Jpm5EAt2Ji2O3il2m5z6d6brsWnoH6HWX-
Nk0lJP3jSnW64uhQ-ulCqAEYZtqP-
6-) ¿qué representa la respuesta en fase de un sistema digital? QztR9DNJvmtiwPkH8pW17UxWVlYbc95e1R7E~oEtV
7j3NefSxQCsIQlkg019vxetDn1bsYVnmxRHjwEuxOf6f
La respuesta en fase de un sistema se representa por el siguiente seeowKUwTUxF1VQWdnYeWjKZQVd7syQsYC79x7L
diagrama [5]; FFRF-
k8ceKfPgHAfAVjLaT0vfqxjsRXo8yulIhxv1k0HHbRJ1
L5Vw~Jf95-~uFb~eiFkvDo1H~91dXNA__&Key-Pair-
Id=APKAJLOHF5GGSLRBV4ZA.
[3] K. Ogata, Sistemas de control en tiempo discreto. Pearson
Educación, 1996.
[4] B. C. Kuo, Sistemas de control automático. Pearson
Educación, 1996.
[5] R. Palacios y J. Ernesto, «Simulación en la computadora
de filtros recursivos digitales con respuesta de fase lineal
de aplanamiento máximo», dic. 1984, Accedido: nov. 05,
2020. [En línea]. Disponible en:
Figura 10. Respuesta en fase http://bibdigital.epn.edu.ec/handle/15000/9246.

RESULTADOS Y CONCLUSIONES

• Se reconoció los diferentes métodos de diseño de

filtros como La invariancia al impulso,
Transformación Bilineal, para filtros análogos IIR y de
ventanas, rizado constante óptimo y muestreo en
frecuencia para filtros digitales FIR

Downloaded by ronald garzon ([email protected])