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

浙江工业大学之江学院 信息工程分院课程报告 课程 数字信号处理 姓名 史佳琦 班级 通信902班 学号 200920110219 教师 李春晓 题目： 窗函数法设计FIR数字滤波器 日期：11月11日 一、实验原理 二、内容及要求 三、程序清单及说明 四、运行结果 五、调试问题及解决方法 一、实验原理 窗函数设计法是一种通过截短方法使无限长非因果序列成长为有限长脉冲响应序列的设计方法。 利用加窗函数进行阶段和平滑，实现一个物理可实现且具有线性相位的FIR滤波器。 二、内容及要求 1.了解一个实际滤波器的设计过程。 2.掌握用窗函数法设计FIR数字滤波器的原理和窗函数对数字滤波器性能的影响。 题目： 1.窗函数法设计线性相位FIR低通滤波器，wp=0.2pi，ws=0.3pi，Rp =3db，Rs =40db 2.窗函数法设计线性相位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; 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; 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.要注意题目中给出的指标是不是可以直接带入的，还是需要转换一下的参数。 2.还要熟悉了解各种MATLAB语言工具函数，都需要凭借这些相关函数来求解问题。