PRÁCTICA 3
Aplicación de la Transformada Z a señales y sistemas.
Tratamiento y Tranmisión de señales
Rubén Esteban Sánchez
14/04/2025
1
Índice
1. Introducción 3
2. Ejercicio 1 3
2.1. Calcule y represente gráficamente la respuesta impulsiva h(n) para n= -20, . . . , 100. Uti-
lice la función impz o filter de Matlab. . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.2. Obtenga H(z) y su diagrama de ceros y polos. ¿Es estable el sistema? Utilice el comando
zplane y roots de Matlab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.3. Represente la magnitud y fase de H(z) como función de la frecuencia. Utilice comando
freqz de Matlab. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.4. Calcule y represente gráficamente la respuesta del sistema al escalón unidad para n= -20,
. . . , 100. Utilice la función stepz o filter de Matlab.. . . . . . . . . . . . . . . . . . . . . 5
3. Ejercicio 2 6
3.1. La respuesta al impulso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3.2. El diagrama de ceros y polos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3.3. La salida y(n) cuando la entrada es de la forma x(n) = 20,9nu(n). . . . . . . . . . . . . 8
4. Ejercicio 3 8
4.1. Utilizando el comando freqz, represente gráficamente la magnitud y fase de la respuesta
en frecuencia del filtro. Obtenga el valor numérico de la magnitud y fase en ω = π/3 y
en ω = π . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
4.2. Genere 200 muestras de la señal x(n) = sin(πn/3) + 5cos(πn) y fı́ltrela mediante el
sistema anteriormente descrito. ¿Cómo se modifican las amplitudes y fases de las dos
sinusoides? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
5. Efecto Wah Wah 10
6. Conclusiones 11
2
1. Introducción
El objetivo de esta práctica consiste en familializarse con las herrmientas disponibles en MATLAB
para el procesado básico de señales, además de la implementacion de un efecto electrónico muy utilizado
en difrentes instrumentos electrónicos como guitarra y bajo, conocido como efecto Wah-Wah.
2. Ejercicio 1
Sea el sistema LTI:
y(n) − y(n − 1) + 0,9y(n − 2) = x(n)
2.1. Calcule y represente gráficamente la respuesta impulsiva h(n) para n= -20, . . . , 100.
Utilice la función impz o filter de Matlab.
Se calcula la respuesta al impulso del sistema, que es su salida cuando la entrada es una señal
delta ([n], 1 en n=0 y 0 en el resto). Usamos filter para controlar los ı́ndices negativos (n = -20 a 100).
La entrada x(n) es un impulso centrado en n=0, desplazado adecuadamente para representar el rango
deseado. Por último, se representa la respuesta al impulso. Como vemos, la señal se mantiene en cero
hasta que entra el impulso, entonces la salida oscila y va convergiendo poco a poco, a partir de unas 80
muestras podemos ver como se va estabilizando, lo que en principio puede ser una respuesta algo lenta
ya que estamos ante un impulso unitario.
Figura 1: Respuesta al Impulso
2.2. Obtenga H(z) y su diagrama de ceros y polos. ¿Es estable el sistema? Utilice el co-
mando zplane y roots de Matlab .
Para realizar este ejercicio podemos representar los polos y ceros con zplane (Figura 3), donde
podemos estudiar la estabilidad visualmente comprobando que todos los polos estén o no dentro de la
3
circunferencia unidad. Otra forma de hacerlo es usando la función roots (Figura 4) para el numerador y
el denominador, la cual nos devuelve el valor numérico de los polos y los ceros. En este caso, se comple-
mentan ambas funciones, además de añadir una tercera que se encarga de comprobar automáticamente
si el sistema es estable y nos lo indica en la ventana de comandos (Figura 5). De la forma que sea, la
respuesta es que el sistema es estable.
Figura 2: Estabilidad: Zplane
Figura 3: Estabilidad: Roots
Figura 4: Estabilidad: Comprobación Automática
4
2.3. Represente la magnitud y fase de H(z) como función de la frecuencia. Utilice coman-
do freqz de Matlab.
En este apartado se evalúa la respuesta en frecuencia del sistema mediante freqz, que calcula H(e jω )
para un rango de frecuencias ω ∈ [0, π]. La respuesta en frecuencia indica cómo el sistema modifica la
amplitud y la fase de los componentes en frecuencia de la señal. En este caso, apreciamos que se trata de
un filtro pasa banda (Figura ), el cual permite al paso de frecuencias en un rango especifico y atenua las
demás.
Figura 5: Respuesta e frecuencia
2.4. Calcule y represente gráficamente la respuesta del sistema al escalón unidad para
n= -20, . . . , 100. Utilice la función stepz o filter de Matlab..
En la figura 6, podemos apreciar la respuesta al escalón unitario y ver como el sistema aunque en el
primer instante alcanza el doble de la amplitud de entrada y tiene bastantes oscilaciones, al final converge
al valor 1. Al igual que en la respuesta al impulso, se va estabilizando en torno a las 80 muestras.
5
Figura 6: Respuesta Escalón Unitario
3. Ejercicio 2
Para el siguiente sistema LTI:
y(n) = x(n) + 0,5x(n − 1) − 0,5y(n − 1) + 0,25y(n − 2)
3.1. La respuesta al impulso
Procediendo como en el ejercicio anterior, la respuesta al impulso queda reflejada en la figura 7. En
este sistema vemos que tiene menos oscilaciones que el anterior y que además se estabiliza mucho más
rápido, en torno a las 20 muestras.
6
Figura 7: Respuesta al Impulso 2
3.2. El diagrama de ceros y polos
En este caso representando el diagrama de polos y ceros, obtenemos que tambı́en el sistema es
estable.
Figura 8: Diagrama de Polos y Ceros
7
3.3. La salida y(n) cuando la entrada es de la forma x(n) = 20,9nu(n).
Como la señal de entrada no es fácil de imaginar de primeras, en la figura 9, se han representado
tanta la entrada como la salida, para que se pueda entender con más fácilidad. En este caso, vemos como
la señal de salida tras unas pequeñas oscilaciones, se va acercando a la de entrada y al cabo de unas 20
muestras, toma totalmente comportamiento exponencial sin oscilaciones.
Figura 9: Respuesta Exponencial
4. Ejercicio 3
Un filtro digital se encuentra descrito a partir de la siguiente ecuación en diferencias
y(n) = x(n) + x(n − 1) + 0,9y(n − 1) − 0,81y(n − 2)
4.1. Utilizando el comando freqz, represente gráficamente la magnitud y fase de la res-
puesta en frecuencia del filtro. Obtenga el valor numérico de la magnitud y fase en
ω = π/3 y en ω = π
. El codigo y los valores obtenidos se reflejan en la figura 10.
8
Figura 10: Código y valores obtenidos
4.2. Genere 200 muestras de la señal x(n) = sin(πn/3) + 5cos(πn) y fı́ltrela mediante el
sistema anteriormente descrito. ¿Cómo se modifican las amplitudes y fases de las
dos sinusoides?
En la figura 11, se muestra la señal de entrada y la de salida, podemos ver que las amplitudes de la
señal Sen(π/3) son amplificadas, mientras que las de la señal Cos(π) quedan atenuadas. Esto es porque
tenemos un filtro Pasa Banda, que amplifica las frecuencias que están muy proxima al ω = π/3 y atenúa
tanto las que están por encima como por abajo de dicha frecuencia.
9
Figura 11: Señal de entrada y Señal filtrada
5. Efecto Wah Wah
El efecto Wah-Wah consiste en un filtro pasa banda cuya frecuencia central es variable. Cuando esa
frecuencia central del filtro varı́a de un valor bajo a uno alto, produce un sonido similar a una voz humana
pronunciando la sı́laba ”wah”, simulando los formantes de la voz humana.
La implementación de este efecto en MATLAB se lleva a cabo con el siguiente código
10
Figura 12: Efecto Wah Wah
6. Conclusiones
Esta práctica ha permitido aplicar conceptos de la transformada Z para el análisis de sistemas LTI
discretos. Se comprobó cómo una ecuación en diferencias puede representarse mediante una función
de transferencia H(z), y cómo esta describe completamente el comportamiento del sistema. Además ha
permitido, conocer uno de los diferentes efectos que se utilizan en el mundo de la música, algo que me
ha llamado bastante la atención. Serı́a muy interesante conocer otros tipos de filtros similares
11
