有限冲激响应数字滤波器设计

1、数字信号处理课程设计 有限冲激响应数字滤波器设计一、实验目的：加深对数字滤波器的常用指标和设计过程的理解。二、实验原理：低通滤波器的常用指标：数字滤波器技术指标：通带边缘频率，阻带边缘频率 ，通带最大衰减，最小阻带衰减通带峰值起伏，最小阻带衰减。数字滤波器有IIR和FIR两种类型，它们的特点和设计方法不同。 三、设计内容： 设计FIR低通滤波器，通带边界频率 ，阻带边界频率 ，阻带衰减不小于50dB。选择一个合适的窗函数，确定脉冲响应，并给出所设计的滤波器的频率响应图。（要求设计完后，再选择其他的窗函数，并比较结果，分析原因）设计步骤：¨ 给定理想的频率响应函数及技术指 通带截止频率

2、为 通带起始频率为 阻带衰减相当于 =50dB低通滤波器的过度带要求理想低通滤波器的截止频率 ¨ 根据过渡带宽度确定N值 由于=50dB，查表可选海明窗，其阻带最小衰减-53dB满足要求。 所要求的过度带宽 由于海明窗过度带宽满足 所以¨ 求所设计的FIR滤波器的单位抽样响应海明窗 所以窗系数需要实现用窗函数blackman（N）, hamming(N),hanning(N) wd=hanning(N);%数组wd返回N点汉宁窗函数wd=hamming(N);%数组wd返回N点汉明窗函数wd=blackman(N);%数组wd返回N点布莱克曼窗函数程序如下wn=hammin

3、g(66);%产生窗函数，N=66nn=0:1:65;alfa=(66-1)/2;hd=sin(0.3*pi*(nn-alfa+eps)./(pi*(nn-alfa+eps);%eps为一个非常小的数，防止出%现零为除数h=hd.*wn' %hd为一个行向量，wn为一个列向量，需要转置h1,w1=freqz(h,1);plot(w1/pi,20*log10(abs(h1);%通过求对数得到以db为单位的幅度特性axis(0,1,-100,10); grid;xlabel('归一化频率/pi');ylabel('幅度/dB');再分别用布拉克曼窗和汉宁窗得

4、到图形如下布拉克曼窗wn=blackman(66);%产生窗函数，N=66nn=0:1:65;alfa=(33-1)/2;hd=sin(0.3*pi*(nn-alfa+eps)./(pi*(nn-alfa+eps);%eps为一个非常小的数，防止出%现零为除数h=hd.*wn' %hd为一个行向量，wn为一个列向量，需要转置h1,w1=freqz(h,1);plot(w1/pi,20*log10(abs(h1);%通过求对数得到以db为单位的幅度特性axis(0,1,-100,10); grid;xlabel('归一化频率/pi');ylabel('幅度/dB');汉宁窗wn=hanning(33);%产生窗函数，N=33nn=0:1:32;alfa=(33-1)/2;hd=sin(0.3*pi*(nn-alfa+eps)./(pi*(nn-alfa+eps);%eps为一个非常小的数，防止出%现零为除数h=hd.*wn' %hd为一个行向量，wn为一个列向量，需要转置h1,w1=freqz(h,1);plot(w1/pi,20*log10(abs(h1);%通过求对数得到以db为单位的幅度特性axis(0,1,-100,10); grid;xlabel('归一化频率/pi