Filtro Fir

Universidad Autónoma de Baja California

Facultad De Ciencias Quimicas E Ingenieria

Ing. en Electrónica

Examen PDS Unidad IV

M.C. José Jaime Esqueda Elizondo

Víctor Samuel RodríguezEscalante

Tijuana, B.C., a 1 de Diciembre de 2012
Diseñar un filtro FIR tipo pasabanda de orden 50+ las últimas dos cifras de su matrícula que deje pasar señales que estén entre5 y 10 Khz, utilizando una frecuencia de muestreo de 40 Khz.

Matricula 277455
%Filtro pasabanda
fs=40000; %Frecuencia de muestreo
Orden=50+55;
Wp = [5e3 10e3]/(fs/2);flag='scale'; % Second-order section scaling.
win = hamming(Orden+1);
filtro=fir1(Orden,Wp,'bandpass',win,flag);

1. Obtener la salida del filtro ante una señal formada por lasuma de senoides de
amplitud unitaria con frecuencias de: 2500 Hz, 7500 Hz y 12 500 Hz. Utilizar 200
muestras de la señal.

% Graficación de la señal en el tiempo
N=0:1:199;y=(sin(2*pi*N*2500/fs))+(sin(2*pi*N*7500/fs))+(sin(2*pi*N*12500/fs));
% Graficar la señal original
figure(1)
subplot(211)
plot(N,y)
title('SEÑAL ORIGINAL')
xlabel('Tiempo (s)')
ylabel('Amplitud (V)')subplot(212)
y_out=conv(filtro,y);
plot(y_out)
title('SEÑAL FILTRADA')
xlabel('Tiempo (s)')
ylabel('Amplitud (V)')

2. Obtener el espectro de amplitud de la señal de entrada al filtro, así como superiodograma.

Easin=fft(y,Orden);
figure(2)
subplot(211)
stem(abs(Easin))
title('Espectro de amplitud de la señal de entrada al filtro')
xlabel('Tiempo (s)')
ylabel('Amplitud (V)')Pin=periodogram(y);
subplot(212)
stem(Pin)
title('Periodograma de entrada');
xlabel('Tiempo (s)')
ylabel('Amplitud (V)')

3. Obtener el espectro de amplitud de la señal filtrada, así como superiodograma.

Easfil=fft(y_out,Orden);
figure(3)
subplot(211)
stem(abs(Easfil))
title('Espectro de amplitud de la señal filtrada')
xlabel('Tiempo (s)')
ylabel('Amplitud (V)')
Pfil=periodogram(y_out);...