一、实验目的:
研究数字滤波器设计思想,理解数字频率与模拟频率的关系,掌握数字系统处理模拟信号的方法。掌握窗函数设计FIR数字滤波器的方法,理解FIR滤波器的线性相位重要意义。
二、实验要求:
设计FIR数字低通滤波器,截止频率,在不同窗口长度(N=15,33)下,分别求出h(n),通过幅频特性和相频特性,观察3dB带宽和20dB带宽,总结窗口长度N对滤波特性的影响;
三、技术方案
设计滤波器并滤波的主要步骤有八步,首先写出,注意需要将模拟频率对采样率Fs归一化,然后计算出过度带宽度,根据阻带最小衰减选择窗函数,再求窗长度,求6dB截止频率,计算FIR滤波器系数向量,绘制滤波器频率特性波特图,完成滤波。
在滤波器设计实验中,第一部分是设计FIR数字低通滤波器,截止频率,在不同窗口长度(N=15,33)下,分别求出h(n),通过幅频特性和相频特性,观察3dB带宽和20dB带宽,总结窗口长度N对滤波特性的影响,由于已知,窗长度已知,我们可以直接列出FIR滤波器系数向量,然后绘制滤波器频率特性波特图,比较不同窗口长度下的波特图。
四、实验程序
clc;
clear;
close all;
N1=15;
wc1=pi/4;
hn1=fir1(N1-1,wc1/pi,'low',bartlett(N1));
figure(1),freqz(hn1,1);
N2=33;
wc2=pi/4;
hn2=fir1(N2-1,wc2/pi,'low',bartlett(N2));
figure(2),freqz(hn2,1);
该代码以三角窗为例,读者可以按需自行修改。
五、实验结果
该图为矩形窗窗口长度15时的幅频特性和相频特性曲线
该图为矩形窗窗口长度33时的幅频特性和相频特性曲线
该图为汉宁窗窗口长度15时的幅频特性和相频特性曲线
该图为汉宁窗窗口长度33时的幅频特性和相频特性曲线
该图为三角窗窗口长度15时的幅频特性和相频特性曲线
该图为三角窗窗口长度33时的幅频特性和相频特性曲线
通过对比可以看出:窗口的长度越长,截取的信号就越长,而信号越长,频率分辨率就越高。而不同的窗函数过渡带宽度,旁瓣峰值,阻带最小衰减不同,导致同样的窗口长度下幅频特性和相频特性曲线仍具有较大差异。