Grupo 2
Grupo 2
Grupo 2
𝑿[𝒁] = ∑ 𝒙(𝒏)𝒁−𝒏
𝒏=−∞
Diagrama de bode
− 𝑎2 ∑ 𝑦[𝑛 − 2]𝑧 −𝑛
¿Qué es la respuesta en fase de un sistema?
𝑛=−∞
El retraso de grupo y de fase de un filtro
proporciona una magnitud medible para
Tomamos el segundo y tercer termino
cuantificar el grado en el cual el filtro modifica la
𝒍= 𝒏−𝟏 → 𝒏= 𝒍+𝟏
fase de una señal x[n] . Si consideramos una señal
𝒌=𝒏−𝟐 → 𝒏= 𝒌+𝟐
formada por componentes de múltiples
Remplazamos datos correspondientes
frecuencias, el retraso de fase del filtro es la ∞
magnitud temporal que el filtro introduce como 𝑌(𝑧) = 𝑏0 ∑ 𝑥[𝑛]𝑧 −𝑛
retraso a cada componente. Matemáticamente el
𝑛=−∞
retraso de fase Tf se define como el cociente ∞
cambiado de signo de la fase q(W) respecto de la − 𝑎1 ∑ 𝑦[𝑙]𝑧 −(𝑙+1)
frecuencia W 𝑛=−∞
∅(Ω) ∞
𝑇𝑓 = -
Ω
mientras que el retraso de grupo Tg se expresa − 𝑎2 ∑ 𝑦[𝑘]𝑧 −(𝑘+2)
𝑛=−∞
como:
∅(Ω)
𝑇𝑓 = - ∞
Ω
𝑌(𝑧) = 𝑏0 ∑ 𝑥[𝑛]𝑧 −𝑛
Un filtro con una característica de fase no lineal 𝑛=−∞
∞
originará una distorsión de fase puesto las
componentes de distinta frecuencia al ser − 𝑎1 ∑ 𝑦[𝑙]𝑧 −𝑙 𝑧 −1
procesadas por el filtro tendrán un retraso que no 𝑛=−∞
∞
será proporcional a la frecuencia y por lo tanto se
alterará la relación original entre los distintos − 𝑎2 ∑ 𝑦[𝑘]𝑧 −𝑘 𝑧 −2
𝑛=−∞
armónicos que la forman. En un filtro de fase
lineal no hay distorsión de fase porque el cambio
de fase del filtro es lineal
𝒃𝟎 = 𝟏
𝒂𝟏 = 𝟎. 𝟓
𝒂𝟐 = 𝟎. 𝟗
𝑏0
𝐻(𝑧) =
1 + 𝑎1 𝑧 −1 + 𝑎2 𝑧 −2
1
𝐻(𝑧) =
1 + 0.5𝑧 −1 + 0.9𝑧 −2
Ecuación matemática 𝑏0 ∏𝑀
𝑘=1(𝑧 − 𝑐𝑘 )
∞ 𝐻(𝑧) = 𝑧 (𝑁−𝑀) ( ) 𝑵
𝑎0 ∏𝑘=1(𝑧 − 𝑝𝑘 )
𝑋 + (𝑧) = 𝑍 + {𝑥[𝑛]} = ∑ 𝑥[𝑛]𝑧 −𝑛
𝑛=0
Ecuación bilateral tiempo discreto
∞
4) ¿Cómo se calcula los polos y ceros de una Cada estudiante realizará la transformada
transformada Z? Z de la ecuación de diferencias. Esta debe
El método por el cual evaluamos los polos realizarse en el editor de ecuaciones de
en las expresiones racionales y Word. No se aceptan pantallazos.
encontramos una solución en las Para realizar la transformada de Z de la
desigualdades con el beneficio de ecuación aplicamos lo siguiente:
𝑇𝑧
𝑥 (𝑛 − 𝑘) → 𝑥(𝑧) 𝑧 −𝑘 ⃓ 𝒂 + 𝒃𝒋 ⃓ = √𝒂𝟐 + 𝒃𝟐
Según la ecuación Original: Donde a y b son los coeficientes del
𝑦[𝑛] = 𝑏0 𝑥[𝑛] + 𝑏1 𝑥[𝑛 − 1] + 𝑎1 𝑦[𝑛 − número imaginario (𝒂 + 𝒃𝒋)
1] Nos queda.
𝑏0 +𝑏1 (cos(𝑤)−𝑗𝑠𝑒𝑛(𝑤))
𝑯(𝝎) =
𝟏+ 𝑎1 (cos(𝑤)−𝑗𝑠𝑒𝑛(𝑤)
𝑏0 +𝑏1 cos(𝑤)−𝑏1 𝑗𝑠𝑒𝑛(𝑤))
𝑦[𝑧] = 𝑏0 𝑥[𝑧] + 𝑏1 𝑥[𝑧] 𝑧−1 𝑯(𝝎) =
𝟏+ 𝑎1 cos(𝑤)−𝑎1 𝑗𝑠𝑒𝑛(𝑤)
+ 𝑎1 𝑦[𝑧] 𝑧 −1 Para hallar la magnitud multiplicamos por el
complejo conjugado del
Denominador y hacemos cambio de variable.
Una vez se tenga la transformada Z de la 𝐻(𝜔) =
𝐴−𝑗𝐵
∗
𝐶+𝑗𝐷
ecuación de diferencia, cada estudiante 𝐶−𝑗𝐷 𝐶+𝑗𝐷
𝐴𝐶+𝐴𝑗𝐷−𝑗𝐵𝐶+𝐵𝐷
hallará la función de transferencia del 𝐻(𝜔) =
𝐶 2 +𝑗𝐷𝐶−𝑗𝐷𝐶+𝐷2
sistema H(Z). Esto también se realizará 𝐴𝐶+𝐵𝐷+𝑗(𝐴𝐷−𝐵𝐶)
con el editor de ecuaciones de Word. 𝐻(𝜔) =
𝐶 2 +𝐷2
𝐴𝐶+𝐵𝐷 𝑗(𝐴𝐷−𝐵𝐶)
Recuerde que la función de transferencia 𝐻(𝜔) = +
𝐶 2 +𝐷2 𝐶 2 +𝐷2
es:
𝐴𝐶+𝐵𝐷 (𝐴𝐷−𝐵𝐶) 2
|𝐻(𝜔)| = √( 2 2 )2 +( )
𝐶 +𝐷 𝐶 2 +𝐷2
𝑦[𝑧] = 𝑏0 𝑥[𝑧] + 𝑏1 𝑥[𝑧] 𝑧 −1 𝐴𝐶+𝐵𝐷2
( 2 2 )
+ 𝑎1 𝑦[𝑧] 𝑧 −1 𝜃 𝐻(𝜔) = arctan 𝐶 +𝐷
(𝐴𝐷+𝐵𝐶)2
( )
𝒀(𝒁) 𝐶2 +𝐷2
𝑯(𝒁) =
𝑿(𝒛)
𝒚 (𝒛)(𝟏 + 𝑎1 𝑧 −1 ) = 𝒙(𝒛)( 𝑏0 + 𝑏1 𝑧 −1 )
𝑏0 ∏𝑚
𝑘=1(𝑧 − 𝑐𝑘 )
𝐻 (𝑧) = 𝑧 𝑁−𝑀 ( ) 𝑁
𝑎0 ∏𝑘=1(𝑧 − 𝑝𝑘 )
Una vez se cuente con la respuesta en
frecuencia del sistema, se hallará la
magnitud de la respuesta en frecuencia, 5) ¿Qué es la transformada Z?
para ello se aplicará la identidad de Euler, R/ En procesamiento de señales, la transformada Z
que según el caso se podría utilizar convierte una señal real o compleja definida en el
cualquiera de las siguientes ecuaciones: dominio del tiempo discreto en una representación
𝒆𝒋𝒘 = 𝐜𝐨 𝐬(𝒘) + 𝒋𝒔𝒊𝒏(𝒘) en el dominio de la frecuencia compleja.
𝒆−𝒋𝒘 = 𝐜𝐨 𝐬(𝒘) − 𝒋𝒔𝒊𝒏(𝒘) 𝑥(𝑧) = 3 + 2𝑧 −1 + 𝑧 −2
Para hallar la función de magnitud,
recordar utilizar la siguiente ecuación:
∞
Pueden existir ceros o polos en z = ∞ Si H(∞) = 0
𝑋[𝑧] = ∑ 𝑥[𝑛] 𝑧 −𝑛 existe un cero
𝑛=−∞ En z = ∞ y si H(∞) = ∞ existe un polo en z = ∞ .
El número total de ceros y polos de un sistema,
considerando los que están en el infinito deben
6) ¿Qué representa Z en una función? coincidir.
R/ Se puede deducir de la ecuación funcional de la
función Z de Riemann que la función Z es real para 𝑏0 ∏𝑚
𝑘=1(𝑧 − 𝑐𝑘 )
valores reales de t. Es una función impar, y 𝐻 (𝑧) = 𝑧 𝑁−𝑀 ( ) 𝑁
𝑎0 ∏𝑘=1(𝑧 − 𝑝𝑘 )
analítica para valores del argumento reales. De
hecho, se puede observar que la función theta de
Riemann-Siegel y la función zeta de Riemann son
9) ¿Qué es la respuesta en frecuencia de
ambas holomorfas en la recta crítica, y donde la
un sistema?
parte imaginaria de t está comprendida entre -1/2 y R/ Se denomina respuesta en frecuencia a la
1/2, la función Z es holomorfa en el rango crítico respuesta en estado estable de un sistema sujeto a
también. una señal sinusoidal de amplitud fija pero a una
frecuencia variable en cierto rango. Se obtiene la
𝑧 = 𝐴𝑒 𝑗𝜔 función de transferencia para el elemento o
Donde A es el módulo de z, y 𝟂 es la frecuencia sistema.
angular. 𝑏0 + 𝑏1 𝑧 −1
(𝜔) =
1 − 𝑎1 𝑧 −1
7) ¿Cuál es la diferencia entre la
transformada Z bilateral y la 10) ¿Qué es la respuesta en fase de un
unilateral? sistema?
R/ La transformada Z de una secuencia en tiempo R/ La respuesta de fase de un filtro proporciona
discreto X[n], donde Z es una variable compleja. una magnitud medible para cuantificar el grado en
Otra notación para la sumatoria es Z( X[n] ). Si la el cual el filtro modifica la fase de una señal x[n] .
secuencia es causal, la transformada Z se convierte Si consideramos una señal formada por
en Esta transformada se llama unilateral, para componentes de múltiples frecuencias, el retraso
distinguirla de la primera definición que toma el de fase del filtro es la magnitud temporal que el
nombre de la transformada Z bilateral. La filtro introduce como retraso a cada
𝒃
transformada Z unilateral es de gran utilidad en el Ɵ(𝒂 + 𝒃𝒋) = 𝒂𝒓𝒄𝒕𝒂𝒏 ( )
análisis de sistemas causales, especificados por 𝒂
ecuaciones en diferencias, con coeficientes
constantes y con condiciones iniciales, es decir,
aquellos que en su inicio no se encuentran en
reposo.
∞
+ (𝑧) + {𝑥[𝑛]}
𝑋 = 𝑍 = ∑ 𝑥[𝑛] 𝑧 −𝑛
𝑛=0
∞
𝑏0 + 𝑏1 𝑒 −𝑗𝜔
𝐻[𝜔] =
1 + 𝑎1 𝑒 −𝑗𝜔
𝑏0 + 𝑏1 (cos 𝜔 − 𝑗 sin 𝜔)
𝐻[𝜔] =
1 + 𝑎1 (cos 𝜔 − 𝑗 sin 𝜔)
III. TRANSFROMADA Z
Siendo 𝑏0 + 𝑏1 cos 𝜔 − 𝑗𝑏1 sin 𝜔
𝑦[𝑛] = 𝑏0 𝑥[𝑛] + 𝑏1 𝑥[𝑛 − 1] − 𝑎1 𝑦[𝑛 − 1] 𝐻[𝜔] =
1 + 𝑎1 cos 𝜔 − 𝑗𝑎1 sin 𝜔
Hallar la transformada z:
𝑦[𝑧] = 𝑏0 𝑥[𝑧] + 𝑏1 𝑥[𝑧]𝑧 −1 − 𝑎1 𝑦[𝑧]𝑧 −1 𝐿𝑎 𝑚𝑎𝑔𝑛𝑖𝑡𝑢𝑑 𝑑𝑒 𝑙𝑎 𝑟𝑒𝑝𝑠𝑢𝑒𝑠𝑡𝑎 𝑒𝑛 𝑓𝑟𝑒𝑐𝑢𝑒𝑐𝑖𝑎
𝑒𝑠 𝑑𝑎𝑑𝑎 𝑝𝑜𝑟 |𝐻(𝜔)|, 𝑝𝑒𝑟𝑜 𝑝𝑎𝑟𝑎 𝑚𝑎𝑛𝑒𝑗𝑜 𝑦 𝑓𝑎𝑐𝑖𝑙𝑖𝑑𝑎𝑑
IV. FUNCIÓN DE TRANSFERENCIA 𝑑𝑒𝑙 𝑝𝑟𝑜𝑐𝑒𝑑𝑖𝑚𝑖𝑒𝑛𝑡𝑜 𝑚𝑎𝑡𝑒𝑚á𝑡𝑖𝑐𝑜 𝑠𝑒 ℎ𝑎𝑐𝑒 𝑒𝑙 𝑠𝑖𝑔𝑢𝑖𝑒𝑛𝑡𝑒
Siendo 𝑐𝑎𝑚𝑏𝑖𝑜 𝑑𝑒 𝑣𝑎𝑟𝑖𝑎𝑏𝑙𝑒𝑠
𝑦[𝑧] = 𝑏0 𝑥[𝑧] + 𝑏1 𝑥[𝑧]𝑧 −1 − 𝑎1 𝑦[𝑧]𝑧 −1
𝐴 = 𝑏0 + 𝑏1 cos 𝜔
Hallar la función de transferencia H[z]: 𝐵 = 𝑏1 sin 𝜔
𝐶 = 1 + 𝑎1 cos 𝜔
𝑦[𝑧] = 𝑏0 𝑥[𝑧] + 𝑏1 𝑥[𝑧]𝑧 −1 − 𝑎1 𝑦[𝑧]𝑧 −1 𝐷 = 𝑎1 sin 𝜔
𝐴𝐶 − 𝑗𝐵𝐶 + 𝑗𝐴𝐷 + 𝐵𝐷
𝐻[𝜔] =
𝐶 2 + 𝐷2
𝐴𝐶 + 𝐵𝐷 − 𝑗(𝐵𝐶 − 𝐴𝐷)
𝐻[𝜔] =
𝐶 2 + 𝐷2
𝐴𝐶 + 𝐵𝐷 𝐵𝐶 − 𝐴𝐷
𝐻[𝜔] = 2 2
−𝑗 2
𝐶 +𝐷 𝐶 + 𝐷2
𝐴𝐶 + 𝐵𝐷 2 𝐵𝐶 − 𝐴𝐷 2
|𝐻[𝜔]| = √( 2 ) + ( )
𝐶 + 𝐷2 𝐶 2 + 𝐷2
𝐵𝐶 − 𝐴𝐷
−
2 2
𝜃 = 𝑎𝑟𝑐𝑇𝑎𝑛 ( 𝐶 + 𝐷 ) Diagrama de bode:
𝐴𝐶 + 𝐵𝐷
𝐶 2 + 𝐷2
−𝐵𝐶 − 𝐴𝐷
𝜃 = 𝑎𝑟𝑐𝑇𝑎𝑛 ( )
𝐴𝐶 + 𝐵𝐷
VI. SIMULACION EN MATLAB
Siendo:
𝑏0 + 𝑏1 𝑧 −1
𝐻[𝑧] =
1 + 𝑎1 𝑧 −1
Simulación 1:
𝑏0 = 0,5 𝑏1 = 0.8 𝑎1 = 0,3
Polos y ceros:
0,5 + 0,8𝑧 −1
𝐻[𝑧] =
1 + 0,3𝑧 −1
0,5 + 0,8𝑧 −1
𝐻[𝑧] =
1 + 1,3𝑧 −1
VII. CONCLUSIONES
REFERENCIAS BIBLIOGRAFICAS
Retomado de:
http://ocw.uv.es/ingenieria-y-
arquitectura/1-1/tema4.pdf- Año 2009-
Diagrama de bode: 2010 Autor: Marcelino Martines Sober.
Retomado de:
https://s3.amazonaws.com/academia.edu.
documents/35148115/TZ00.pdf?response
-content-
disposition=inline%3B%20filename%3D
TZ00.pdf&X-Amz-Algorithm=AWS4-
HMAC-SHA256&X-Amz-
Credential=AKIAIWOWYYGZ2Y53UL
3A%2F20191030%2Fus-east-
1%2Fs3%2Faws4_request&X-Amz-
Date=20191030T231602Z&X-Amz-
Expires=3600&X-Amz-
SignedHeaders=host&X-Amz-
Polos y ceros: Signature=e563a789602be5d9820c26c28
cd9572c4366989f2bde45a64d3266d4155
5f4df Año: 2003 Autor: Juan Sacerdoti.