UNIVERSIDAD CATÓLICA DE SANTA MARÍA”
ESCUELA PROFESIONAL DE INGENIERIA MECÁNICA,
MECÁNICA ELECTRICA Y MECATRÓNICA
PROCESAMIENTO DIGITAL DE SEÑALES
GRUPO : 01
TEMA : RESPUESTA DE FRECUENCIA
PROFESOR : ING. JUAN CARLOS CUADROS MACHUCA
ALUMNO : RAMIREZ ZAMORA JOSE FELIPE BLAIR
2020
1. OBJETIVO
1.1. Obtener la respuesta en frecuencia de sistemas discretos aplicando diferentes
herramientas de software y funciones
contenidas en el toolbox de procesamiento digital de señales de Matlab.
2. MARCO TEORICO
2.1. Respuesta en frecuencia
El interés y uso en el PDS de la respuesta en frecuencia de sistemas discretos es el
mismo que en sistemas continuos:
determinar la salida de un sistema (en estado estacionario; cuando t ➝∞) cuando la
entrada es una combinación de
sinusoides. En tiempo discreto será cuando n → y la entrada es una armónica
discreta.
Una de las características de un sistema de tiempo discreto lineal e invariante en el
tiempo es que la respuesta en estado
estacionario del sistema a una entrada sinusoidal es otra señal de tipo sinusoidal; la
diferencia entre estas es solamente de
magnitud y fase.
Si conociéramos la forma en que afecta el sistema a una entrada sinusoidal de
cualquier frecuencia podríamos determinar
su respuesta a cualquier señal de entrada ya que todas las señales se pueden
considerar como una combinación lineal de
señales sinusoidales. Por eso resulta conveniente caracterizar a los sistemas con su
respuesta en frecuencia, es decir con la
información sobre el cambio que produce en la magnitud y en la fase de las señales
sinusoidales de entrada antes de llevarlas
a la salida.
Alternativas para encontrar la respuesta en estado estacionario son la convolución y
usar Transformadas Z→ Y(z)=H(z)·X(z)
para, posteriormente, aplicar transformadas inversas. De esta forma se obtiene el
término transitorio y estacionario (para
n→0). Esta forma de cálculo tiene una alta complejidad en la mayoría de las
ocasiones (recordar el análisis en el domnio
temporal n).
La definición de la respuesta en frecuencia de un sistema L.T.I con una respuesta
impulsional h[n] es:
La respuesta en frecuencia es la evaluación de la Transformada Z en la
circunferencia de radio unidad. Para tener dicha
respuesta en frecuencia la Transformada Z debe converger en esa circunferencia; es
decir, nuestra R.O.C debe incluir dicha
circunferencia
De la expresión de la respuesta en frecuencia o de su relación con la Transformada
Z se comprueba que es periódica de
periodo 2π. Además, recordando nuestro rango de trabajo con las frecuencias
digitales sólo habrá que evaluarla en el rango
0≦ w ≦ π. Es decir, hay que evaluar la Transformada Z sólo en media
circunferencia.
De la relación con la Transformada Z es inmediato obtener la siguiente relación:
PROCEDIMIENTO
4.1. Con la ayuda del help de Matlab y bibliografía especializada explicar el
funcionamiento de las funciones de Matlab
mencionadas en la última parte del marco teórico, incluida la herramienta
FVTOOL.
4.2. Determinar en forma analítica la respuesta en frecuencia del sistema dado por
la EDD siguiente: 𝑦[𝑛] = 𝛼𝑦[𝑛 − 1] +
𝑥[𝑛]. Empleando Matlab, realice las gráficas del módulo y la fase en frecuencia.
4.3. A partir de la función compleja, determinada en el ítem IV.1 y con el empleo
de las funciones Matlab: abs y angle ,
realice las gráficas del módulo y la fase en frecuencia y compárelas con las
obtenidas en el ítem 4.2.
4.4. Del sistema del ítem IV.1 dado por la EDD, determine la función de
transferencia del sistema. Luego emplee la función
Matlab freqz para graficar la respuesta en frecuencia del sistema. Compare con los
resultados anteriores.
4.5. Repita los pasos 4.2 al 4.4 para los siguientes sistemas de tiempo discreto
4.6. Dibuje la respuesta del sistema dado por la siguiente FT:
Emplee por lo menos dos métodos distintos, luego compare ambas respuestas.
4.7. Obtener la respuesta en frecuencia y la respuesta al impulso de los siguientes
sistemas. Utilizar preferentemente la
función FVTOOL. En algún caso, obtener lo solicitado empelando métodos vistos
en prácticas pasadas y realizar una
comparación y explicación de los resultados obtenidos
AVANCE
clear all;close all;clc;
%EDD: y[n]-y[n-1]-y[n-2]+y[n-3]=12x[n-2]
a=[1 -1 -1 1];%Definimos el vector A
b=[0 0 12 0];%Definimos el vector B
w=-pi:pi/100:pi;%Intervalo
H=freqz(b,a,w);%Frecuencia
figure(1)
Mod=20*log10(abs(H));
plot(w,Mod);%graficamos
title('Módulo');%titulo
xlabel('Muestras');%eje x
ylabel('Decibelios');%eje y
grid on;
figure(2)
Fas=angle(H);%angulo
plot(w,Fas);%grafiquemos
title('Fase');%titulo
xlabel('Muestras');%eje x
ylabel('Radianes');%eje y
grid on;
figure(3)
impz(b,a);%graficamos el impulso unitario
grid on;
clear all;close all;clc;
%EDD: y[n]-0.5y[n-1]-0.9375y[n-2]+0.5625y[n-3]=x[n-1]
a=[1 -0.5 -0.9375 0.5625];%definimos vecto A
b=[0 1];%definimos vector B
w=-pi:pi/100:pi;%intervalo
H=freqz(b,a,w);%frecuencia
figure(1)
Mod=20*log10(abs(H));
plot(w,Mod);%graficamos
title('Respuesta de frecuencia');%titulo
xlabel('Muestras');%eje x
ylabel('Decibelios');%eje y
grid on;
figure(2)
Fas=angle(H);
plot(w,Fas);%graficamos
title('Fase');%titulo
xlabel('Muestras');%eje x
ylabel('Radianes');%eje y
grid on;
figure(3)
impz(b,a);%graficamos impulso unitario
grid on;
