实验三：FIR数字滤波器设计与软件实现.ppt

;1．实验目的 （1）掌握用窗函数法设计FIR数字滤波器的原理和方法。 （2）掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。 （3）掌握FIR滤波器的快速卷积实现原理。 （4）学会调用MATLAB函数设计与实现FIR滤波器。;FS=1000 T=1/FS;实验内容及步骤 （1）认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理； （2）调用信号产生函数xtg产生具有加性噪声的信号xt，并自动显示xt及其频谱，如图10.5.1所示； 图10.5.1 具有加性噪声的信号x(t)及其频谱如图 （3）请设计低通滤波器，从高频噪声中提取xt中的单频调幅信号，要求信号幅频失真小于0.1dB，将噪声频谱衰减60dB。先观察xt的频谱，确定滤波器指标参数。;（4）根据滤波器指标选择合适的窗函数，计算窗函数的长度N，调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序，调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。 （5）重复（3），滤波器指标不变，但改用等波纹最佳逼近法，调用MATLAB函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。; 窗函数法设计FIR滤波器 窗函数法是通过对理想滤波器的单位取样响应加窗来逼近理想滤波器的。 函数fir1用于设计标准的低通、带通、高通和带阻滤波器。 函数fir1的调用格式为 b= fir1(n, Wc, ftype, Windows) 其中，n为滤波器阶数，Wc为截止频率 ftype决定滤波器类型，ftype= high，设计高通FIR滤波器，ftype= stop，设计带阻FIR滤波器。 Windows指定窗函数类型，默认为Hamming窗；可选Hanning、Hamming、Blackman、triangle、bartlett和boxcar窗，每种窗都可以由Matlab的相应函数生成。;前面介绍了FIR滤波器的优化设计方法，通过迭代的方法求解FIR滤波器，过程十分复杂。 在Matlab中，可以调用函数remez实现滤波器的设计。 b= remez(n, f, m) 函数remez采用Parks-McClellan算法设计线性相位FIR滤波器，n为滤波器阶数，其幅频特性由f和m指定。;实验结果与分析 具有加性噪声的信号x(t)及其频谱如图;分析： 用窗函数法设计滤波器，滤波器长度 Nb=184。滤波器损耗函数和滤波器输出yw(nT)分别如图所示。 用等波纹最佳逼近法设计滤波器，滤波器长度 Ne=91。滤波器损耗函数和滤波器输出ye(nT)分别如图所示。 两种方法设计的滤波器都能有效地从噪声中提取信号，但等波纹最佳逼近法设计的滤波器阶数低得多，当然滤波实现的运算量以及时延也小得多，从图10.5.3(b)和（d）可以直观地看出时延差别。;1、如果给定通带截止频率和阻带截止频率以及阻带最小衰减，如何用窗函数设计线形相位低通滤波器？请写出设计步骤。 （1）根据对阻带衰减及过度带的指标要求，选择窗函数的类型，并估计窗口长度N; （2）构造希望逼近的频率响应函数Hd(ejw); （3）计算hd(n)。;2、如果要求用窗函数法设计带通滤波器，且给定通带上下截止频率为wp1和wpu,阻带上下截止频率为ws1和wsu，试求理想带通滤波器的截止频率 希望逼近的理想带通滤波器的截止频率 分别为： Wcl=（Wsl+Wpl)/2 Wcu=(Wsu+Wpu)/2 ;3、解释为什么对同样的技术指标，用等波纹最佳逼近法设计的滤波器阶数低 ①用窗函数法设计的滤波器，如果在阻带截止频率附近刚好满足，则离开阻带截止频率越远，阻带衰减富裕量越大，即存在资源浪费； ② 几种常用的典型窗函数的通带最大衰减和阻带最小衰减固定，且差别较大???又不能分别控制。所以设计的滤波器的通带最大衰减和阻带最小衰减通常都存在较大富裕。如本实验所选的blackman窗函数，其阻带最小衰减为74dB,而指标仅为60dB。 ③ 用等波纹最佳逼近法设计的滤波器，其通带和阻带均为等波纹特性，且通带最大衰减和阻带最小衰减可以分别控制，所以其指标均匀分布，没有资源浪费，所以期阶数低得多。