HUAMANI MESAJEL Jhon Manuel
Codigo: 1111039
UNIVERSIDAD TECNOLÓGICA
DEL PERÚ
TRABAJO AUTÓNOMO
FILTROS DIGITALES
ESTUDIANTE:
HUAMANI MESAJEL, Jhon
CODIGO:
1111039
CURSO: PROCESAMIENTO DIGITAL DE SEÑALES
PROFESOR: GOMEZ VERNAL, Francisco Armando
AÑO:
2017
1
HUAMANI MESAJEL Jhon Manuel
Codigo: 1111039
I. INTRODUCCIÓN
En este informe, se trabajará con conceptos básicos de filtros, y
herramientas de diseño y análisis de MATLAB. Así mismo mencionaremos de
manera resumida la definición de un filtro, los tipos de filtros, sus aplicaciones,
las diferencias entre un filtro digital y un filtro analógico y sus caracterizaciones.
Dentro del campo del procesamiento digital de señales, el concepto de filtro
digital es todo aquel procesamiento que permite la transformación de los datos
digitalizados en otros datos empleando un determinado algoritmo. En esta
ocasión, nos vamos a referirnos al filtro digital como aquel sistema lineal e
invariante en el tiempo: que permite el paso de los componentes de la señal
existentes en un determinado intervalo de frecuencia, y elimina las demás.
Nos abocaremos al análisis de los filtros no recursivos tipo FIR. A
continuación se presenta la formula del filtro digital definido por la siguiente
ecuación.
La transformada Z permite analizar la respuesta de frecuencia de los filtros.
La magnitud de esta respuesta determina cómo se altera la amplitud en función
de la frecuencia, y su respuesta de fase define los retardos en el tiempo
producto del filtro en función de la frecuencia.
2
HUAMANI MESAJEL Jhon Manuel
Codigo: 1111039
II. DEFINICIÓN
Un filtro digital es un sistema que, dependiendo de las variaciones de las
señales de entrada en el tiempo y amplitud, se realiza un procesamiento
matemático sobre dicha señal; generalmente mediante el uso de la
Transformada rápida de Fourier; obteniéndose en la salida el resultado del
procesamiento matemático o la señal de salida.
Los filtros digitales tienen como entrada una señal analógica o digital y en su
salida tienen otra señal analógica o digital, pudiendo haber cambiado en
amplitud, frecuencia o fase dependiendo de las características del filtro digital.
El filtrado digital es parte del procesado de señal digital. Se le da la
denominación de digital más por su funcionamiento interno que por su
dependencia del tipo de señal a filtrar, así podríamos llamar filtro digital tanto a
un filtro que realiza el procesado de señales digitales como a otro que lo haga
de señales analógicas.
III. TIPOS DE FILTROS
Hay varios tipos de filtros así como distintas clasificaciones para estos filtros:
De acuerdo con la parte del espectro que dejan pasar y que atenúan
hay:
Filtros pasó alto.
Filtros paso bajo.
Filtros pasó banda.
Banda eliminada.
Multibanda.
Pasa todo.
Resonador.
Oscilador.
Filtro peine (Comb filter).
Filtro ranura o filtro rechaza banda (Notch filter).
De acuerdo con su orden:
primer orden
segundo orden
De acuerdo con el tipo de respuesta ante entrada unitaria:
FIR (Finite Impulse Response)
IIR (Infinite Impulse Response)
TIIR (Truncated Infinite Impulse Response)
De acuerdo con la estructura con que se implementa:
Laticce
Varios en cascada
Varios en paralelo
3
HUAMANI MESAJEL Jhon Manuel
Codigo: 1111039
IV. APLICACIONES
Separación de señales que fueron combinadas desafortunadamente (ruido,
Interferencias provenientes de otros sistemas)
Recuperación de señales distorsionadas de alguna forma (por ejemplo, al
Ser trasmitidas)
Síntesis de sonido: creación o modificación de señales para moldear
Espectros o formas de onda y lograr el efecto auditivo buscado.
Efectos de audio: chorus, flanger, phaser, reverb
Separación de señales que fueron combinadas desafortunadamente (ruido,
Interferencias provenientes de otro sistema)
4
HUAMANI MESAJEL Jhon Manuel
Codigo: 1111039
Filtros digitales vs. Filtros analógicos
El desempeño de los filtros digitales es ampliamente superior a los filtros
analógicos. En muchas ocasiones, la motivación para muestrear una señal es
emplear un filtro digital.
La gran ventaja de los filtros digitales sobre los analógicos es que presentan
una gran estabilidad de funcionamiento en el tiempo.
El filtrado digital consiste en la realización interna de un procesado de datos de
entrada.
5
HUAMANI MESAJEL Jhon Manuel
Codigo: 1111039
En general el proceso de filtrado consiste en el muestreo digital de la señal de
entrada, el procesamiento considerando el valor actual de entrada y
considerando las entradas anteriores. El último paso es la reconstrucción de la
señal de salida.
En general la mecánica del procesamiento es:
1. Tomar las muestras actuales y algunas muestras anteriores (que
previamente habían sido almacenadas) para multiplicadas por unos
coeficientes definidos.
2. También se podría tomar valores de la salida en instantes pasados y
multiplicarlos por otros coeficientes.
3. Finalmente todos los resultados de todas estas multiplicaciones son
sumados, dando una salida para el instante actual.
El procesamiento interno y la entrada del filtro serán digitales, por lo que puede
ser necesaria una conversión analógica-digital o digital-analógica para uso de
filtros digitales con señales analógicas.
Un tema muy importante es considerar las limitaciones del filtro de entrada
debido a Teorema de muestreo de Nyquist-Shannon que en pocas palabras; si
quiero procesar hasta una frecuencia de 10KHz, debo muestrear a por lo
menos 20 KHz.
Los filtros digitales se usan frecuentemente para tratamiento digital de la
imagen o para tratamiento del sonido digital.
Caracterización de un filtro
Hay tres formas equivalentes de caracterizar un filtro:
Respuesta al impulso
Respuesta en frecuencia
Respuesta al escalón
6
HUAMANI MESAJEL Jhon Manuel
Codigo: 1111039
V. CÓDIGO DE FUENTE PARA UN FILTRO
%% Procesamiento de una señal de audio usando MATLAB
%% Selección del tipo de filtrado
% 1 -> Pasa bajo
% 2 -> Pasa alto
% 3 -> Pasa banda
tipo=3;
%% Leer la señal de audio
[y,fs]=wavread('Yabu_mono');
% y -> muestras de la señal
% fs-> frecuencia de muestreo
% Graficar la señal original
subplot(411)
T=length(y)/fs;
t=linspace(0,T,T*fs);
plot(t,y)
title('SEÑAL ORIGINAL')% Título
xlabel('Tiempo (s)') % Etiqueta del eje X
ylabel('Amplitud (V)') % Etiqueta del eje Y
% xlim([0 20/]) % Límite de la señal
%% Reproducir la señal de audio
% wavplay(0.1*y,fs)
%% FFT de la señal
subplot(412)
% Llamado a la función que calcula la FFT
fft_signal(y,fs);title('ESPECTRO DE LA SEÑAL ORIGINAL')
xlim([0 4e3])
%% Filtrado de la señal
switch tipo
case 1
% Cálculo de los coeficientes del filtro (filtro pasa bajas)
% Este filtrado deja solo la señal por debajo de 500 Hz
titulo='FILTRO PASA BAJAS';
% Frecuencia normalizada
fNorm = 1000 / (fs/2);
[b,a] = butter(10, fNorm, 'low');
case 2
% Cálculo de los coeficientes del filtro (filtro pasa altas)
% Este filtrado deja solo la señal por encima de 500 Hz
titulo='FILTRO PASA ALTAS';
% Frecuencia normalizada
fNorm = 1000 / (fs/2);
[b,a] = butter(10, fNorm, 'high');
otherwise
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% Cálculo de los coeficientes del filtro (filtro pasa banda)
% Este filtrado deja solo la señal de 2KHz a 3KHz
% Frecuencias normalizadas
titulo='FILTRO PASA BANDA';
Wp = [2e3 3e3]/(fs/2); Ws = [1.5e3 3.5e3]/(fs/2);
Rp = 3; Rs = 40; % Rizado de la banda de paso y de parada (s)
[n,Wn] = buttord(Wp,Ws,Rp,Rs);% Orden del filtro y frecuencia de corte óptima
7
HUAMANI MESAJEL Jhon Manuel
Codigo: 1111039
[b,a] = butter(n,Wn); % Coeficientes del filtro
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
end
% Filtrado de la señal
y_Low = filtfilt(b, a, y);
% Graficación de la señal en el tiempo
subplot(413)
plot(t,y_Low)
title('SEÑAL FILTRADA')
xlabel('Tiempo (s)')
ylabel('Amplitud (V)')
% xlim([0 20/f0])
% Graficación de la señal en frecuencia
subplot(414)
% Llamado a la función que calcula la FFT
fft_signal(y_Low,fs);title('ESPECTRO DE LA SEÑAL FILTRADA')
xlim([0 4e3])
%% Gráficas del filtro
% Respuesta en frecuencia del filtro
[H,w]=freqz(b,a,512,fs);
figure(2)
%Trazado de la respuesta en Magnitud
subplot(221)
plot(w,20*log10(abs(H)));
grid on;
title ([titulo, ' Respuesta en magnitud']);
xlabel('Frecuencia (Hz)');
ylabel('H(f) db')
xlim([0 5e3])
% Respuesta en fase
subplot(222)
plot(w,angle(H));
grid on;
title ([titulo,' Respuesta en fase']);
xlabel('Frecuencia (Hz)')
ylabel('ángulo de H rad')
xlim([0 5e3])
%Respuesta al impulso
subplot(223)
[y_eje,tt]= impz(b,a,60);
stem(tt,y_eje);
title ([titulo,' Respuesta al impulso']);
xlabel('n')
ylabel('h[n]')
%Ploteo de los polos y ceros
z= roots(b); % Ceros
p = roots(a); % Polos
subplot(224)
zplane(z,p)
title('Polos y ceros')
legend('Ceros','Polos')
%% Reproducción de audio de entrada y salida
pause(2)
disp('Audio de entrada')
% Se multiplica por 0.2 para atenuar la salida del tono por la bocina
% wavplay(0.2*y,fs)
disp('Audio de salida (señal filtrada)')
% Se multiplica por 0.2 para atenuar la salida del tono por la bocina
% wavplay(0.2*y_Low,fs)
8
HUAMANI MESAJEL Jhon Manuel
Codigo: 1111039
Gráficos