数字信号处理： FIR数字滤波器的设计实验

FIR数字滤波器的设计实验

1、实验目的
（1）学习采用Matlab程序应用窗函数法设计FIR数字低通滤波器；
（2）学习使用Matlab工具fdatool设计数字滤波器；
2、实验原理
窗函数法设计FIR数字低通滤波器的步骤：
（1）根据要求的滤波器阻带衰减值，确定窗函数的类型；
（2）根据要求的滤波器的通带截止角频率和阻带截止角频率，确定过渡带和理想低通滤波器的截止频率；
（3）由窗函数的类型和过渡带，确定滤波器单位冲激响应的长度；
（4）计算得到理想低通滤波器的单位冲激响应；
（5）对理想低通滤波器加窗处理；
（6）得到加窗后的数字低通滤波器的频率响应
3、实验内容及要求
（1）利用窗函数法设计一个偶对称的线性相位FIR数字低通滤波器，给定的技术指标为：通带截止频率，最大通带衰减，阻带截止频率，阻带最小衰减
（2）请选定窗形状，求出及，并检验设计结果。
相关Matlab函数：hamming,hanning, chebwin, sin, freqz,abs,log,min,max等。
3、实验结果

根据实验要求给出相关图表及数据
实验源程序

wp=0.3*pi;
ws=0.45*pi;
deltaw =ws-wp;
NO=ceil(8*pi/deltaw);
N=NO+mod(NO+1,2);
n=0:N-1;
wd=hamming(N)';
wc=(wp+ws)/2;

tao=(N-1)/2;
n=[0:(N-1)];
m=n-tao+eps;
hd=sin(wc*m) ./(pi*m);

h=hd.*wd;
L=1000;
[H,w]=freqz(h,1,L,'whole');
mag=abs(H);
db=20*log10((mag+eps)/max(mag));
dw=2*pi/L;
Rp=-min(db(1:wp/dw+1));
As=-max(db(ws/dw+1:L/2+1));
ang=unwrap(angle(H)/pi*180,180);
figure,
subplot(131);plot(n,wd);grid on;
title('海明窗');xlabel('n');ylabel('w(n)');
subplot(132);plot(w/pi,db);grid on;
title('幅度相应(dB)');xlabel('\omega/\pi');
ylabel('20log|H(e^j^\omega)|(dB)');
axis([0,1,-80,5]);
subplot(133);plot(w/pi,ang);grid on;
title('相位相应(degree)');xlabel('\omega/\pi');
[N,Rp,As]