哈工大数字信号处理上机实验..docx

Harbin Institute of Technology实验报告课程名称：数字信号处理实验题目：用窗函数法设计FIR数字滤波器院系：电子与信息工程学院班级：姓 名：学号：指导教师：冀振元、宿富林 实验时间：2015 年11月哈尔滨工业大学实验原理及目的实验原理：如果所希望的滤波器的理想频率响应函数为，则其对应的单位脉冲响应为用窗函数w(n)将hd(n)截断，并进行加权处理，得到h(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列， 其频率响应函数为如果要求线性相位特性， 则h(n)还必须满足：根据上式中的正、负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。要根据所设计的滤波特性正确选择其中一类。例如，要设计线性相位低通特性，可选择h(n)=h(N-1-n)一类，而不能选h(n)=-h(N-1-n)一类。实验目的：(1)熟悉矩形窗、汉宁窗、海明窗和布莱克曼窗。(2) 掌握用上述窗函数法设计FIR数字滤波器的原理和方法。(3) 熟悉线性相位FIR数字滤波器特性。(4) 了解各种窗函数对滤波特性的影响。二、实验步骤及内容实验步骤：(1) 复习用窗函数法设计FIR数字滤波器一节内容，阅读本实验原理，掌握设计步骤。(2) 编写程序。①编写能产生四种窗函数的子程序。②编写主程序。主程序框图如图 10.6.1 所示，仅供参考。其中幅度特性要求用dB表示。图10.6.1 用窗函数法设计滤波器主程序框图实验内容：设计低通FIR数字滤波器时，一般以理想低通滤波特性为逼近函数，即用四种窗函数设计线形相位低通滤波器。要求在两种窗口长度下，绘制相应的幅频和相频特性曲线，观察3dB和20dB带宽以及阻带最小衰减，比较四种窗函数对滤波器特性的影响。三、实验结果及分析1、矩形窗窗口长度N=15时幅频和相频特性曲线窗口长度N=33时幅频和相频特性曲线2、汉宁窗窗口长度N=15时幅频和相频特性曲线窗口长度N=33时幅频和相频特性曲线3、海明窗窗口长度N=15时幅频和相频特性曲线窗口长度N=33时幅频和相频特性曲线4、布莱克曼窗窗口长度N=15时幅频和相频特性曲线窗口长度N=33时幅频和相频特性曲线分析：N=33的幅频特性优于N=15的幅频特性，更为接近理想滤波器。N的值增大会使过渡带宽减小但是对于阻带最小衰减影响并不大（吉布斯效应）。增大N值，则通带、阻带波纹变化越快，抽样点更密，则平坦区域逼近误差越小。只有选择合适的窗函数才可以满足阻带衰减指标，然后再选择N满足过渡带宽指标。四、实验结论1.利用窗函数设计FIR滤波器，原理简单，计算简单，只要选择合理的窗函数就能接近理想滤波器，但是，其很难控制滤波器的通带边缘。即通带阻带的截止频率不易控制。2. 对于低通滤波器来说：过渡带宽大小依次为布莱克曼窗海明窗汉宁窗矩形窗。阻带衰减的绝对值大小依次为布莱克曼窗海明窗汉宁窗矩形窗。3. 验证了吉布斯效应。4. 设计关键在于窗函数及N的选择。一般根据要求查阅现有的窗函数指标表，可以找到需要的窗函数，及各项指标以及窗函数表达式。窗函数的选择原则是：①具有较低的旁瓣幅度尤其是第一旁瓣幅度。②旁瓣幅度下降速度要大以利增加阻带衰减。③主瓣的宽度要窄以获得较陡的过渡带。通常上述三点很难同时满足。当选用主瓣宽度较窄时，虽然得到较陡的过渡带，但通带和阻带的波动明显增加；当选用最小的旁瓣幅度时，虽能得到匀滑的幅度响应和较小的阻带波动，但过渡带加宽。因此实际选用的窗函数往往是它们的折衷。在保证主瓣宽度达到一定要求的条件下适当牺牲主瓣宽度来换取旁瓣波动的减少。五、思考题1.如果给定通带截止频率和阻带截止频率以及阻带最小衰减，如何用窗函数法设计线性相位低通滤波器? 写出设计步骤答：(1).选择理想低通滤波器作逼近函数，为保证线性相位，取a=(N-1)/2。(2).求单位冲激响应。(3).根据阻带衰减指标，选择窗函数形状，根据允许的过渡带宽Δω选定N值（查表）(4).将与窗函数相乘得FIR数字滤波器的单位冲激响应h(n)。(5).计算FIR数字滤波器的频率响应，并验证是否达到所要求的指标。(6).由计算幅度函数和相位函数θ(ω)。若不满足指标要求，重复步骤(3)-(5)，直到满足要求。2. 如果要求用窗函数法设计带通滤波器，且给定上、下边带截止频率为和，试求理想带通的单位脉冲响应。答：已知、，求(1).选择窗函数的类型和长度点数N。(2).求(3).检验是否满足要求，若不满足要求，则重新选用窗形式或N，再进行计算。注意。六、程序说明：为方便报告分析，使更直观利于观察，报告所用图为单独处理截取图片。%程序说明：用窗函数法设计FIR数字滤波器clear all;close all;clc;N = input(请输入窗口长度N(33或15):);n = 0:N-1; %长度a =