窗函数法设计FIR数字滤波器报告 史佳琦 通信

1、浙江工业大学之江学院信息工程分院课程报告课程 数字信号处理 姓名 史佳琦 班级 通信902班 学号 9 教师 李春晓 题目： 窗函数法设计FIR数字滤波器日期：11月11日一、实验原理二、内容及要求三、程序清单及说明四、运行结果五、调试问题及解决方法一、实验原理 窗函数设计法是一种通过截短方法使无限长非因果序列成长为有限长脉冲响应序列的设计方法。利用加窗函数进行阶段和平滑，实现一个物理可实现且具有线性相位的FIR滤波器。二、内容及要求1.了解一个实际滤波器的设计过程。2.掌握用窗函数法设计FIR数字滤波器的原理和窗函数对数字滤波器性能的影响。题目：1.窗函数法设计线性相位FIR低通滤波器，wp

2、=0.2pi，ws=0.3pi，Rp =3db，Rs =40db2.窗函数法设计线性相位FIR带通滤波器，阶数N=15， wp1=0.3pi，wp2=0.5pi(rad)，采用hamming设计。3.编写用改进的升余弦窗函数设计FIR线性相位数字低通滤波器的程序，wc=0.5pi，N=21，滤波器阶数为N-1。 调试运行程序，要求屏幕上显示出单位冲激相应h（n）的数值，画出其幅度相应20log10|H（ejw）|及|H（ejw）|的曲线。画出窗函数w（n）及其频谱幅度相应20log10|H（ejw）|及|H（ejw）|的曲线。 三、程序清单及说明第一题的程序：clear all;fs=100;

3、wp=2*pi*10/fs;ws=2*pi*13/fs;wc=(ws+wp)/2;tr_width=ws-wp;N=ceil(1.8*pi/tr_width)+1;hd =ideal_lp(wc,N);n=0:1:N-1;stem(n,hd);figure;freqz(hd,1,512);w_box=boxcar(N);h_box=fir1(N-1,wc/pi,w_box);figure;freqz(h_box,1,512);figure;stem(n,w_box);第二题的程序：clear all;wp=0.2*pi;ws=0.3*pi;wc=(ws+wp)/2;tr_width=ws-wp

4、;N=ceil(6.2*pi/tr_width)+1;hd =ideal_lp(wc,N);n=0:1:N-1;stem(n,hd);figure;freqz(hd,1,512);w_box=hanning(N);h_box=fir1(N-1,wc/pi,w_box);figure;freqz(h_box,1,512);figure;stem(n,w_box);第三题的程序：clear all;wc=0.5*pi;N=21;hd =ideal_lp(wc,N);n=0:1:N-1;stem(n,hd);figure;freqz(hd,1,512);w_box=hamming(N);h_box=fir1(N-1,wc/pi,w_box);figure;freqz(h_box,1,512);figure;stem(n,w_box);四、运行结果第一题的图形 第二题的图形 第三题图形 五、调试问题及解决方法1.要注意题目