数字信号处理FIR数字滤波器的矩形窗函数法设计.docVIP

设计题目 语音信号的数字滤波—FIR数字滤波器的矩形窗函数法设计 一、 课程设计的目的 通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；掌握利用数字滤波器对信号进行滤波的方法。。[s,fs] = wavread(lishuang.wav); 2）下面代码用于产生信号的时域波形，波形如图2里的1 n = 0:length(s)-1; subplot(411); plot(n,s); title(声音信号的时域波形); 3）由上图可看出，此声音信号在0到1s内基本为空白，即没录上声音，故将其截去，截短后的声音信号的代码如下，声音信号的截短如图2里的2，代码如下： L = 2^floor(log2(length(s))); s1=s(10000:10000+L-1); wavwrite(s1,fs,s1.wav); n = 0:length(s1)-1; subplot(412); plot(n,s1); title(截短的声音信号的时域波形); 图2 原始语音信号的时域图及截短后的时域及频域图 4）下面代码用于产生声音信号的频谱，从图中可以看出信号的主要能量分布，如图2的3 s1 = fft(s1); L1 = 2^floor(log2(length(s1))); k = 0:L1-1; subplot(413); plot(k,abs(s1)); 5）下面代码用于产生信号的单边谱，即取上图中的一半，此图是确定滤波频率范围的参考图，如图2的4 s2 = awgn(s1,27); wavwrite(s2,fs,s2.wav); detf = fs/L1; subplot(414); plot(k(1:L1/2)*detf,abs(S1(1:L1/2))); title(单边谱); （4） 由图2的2可以看出，所截取的语音信号的长度为N = 35000，由图2的4可以看出语音信号的能量主要频段为fL = 200Hz,fH = 1100Hz； 2.3 滤波器的设计 2.3.1确定滤波器的参数 （1）确定参数： 根据2.2中的（4）可知语音信号的主要能量分布的频段为，,,通带截止频率为，阻带截止频率为; 模拟角频率： 数字角频率： （2）数字滤波器类型的选择： 数字滤波器从实现方法上可以分为有限长冲激响应FIR数字滤波器和无限长冲激响应IIR数字滤波器 IIR滤波器的单位冲激响应h(n)是无限长的，即n趋于无穷大。其系统函数H(Z)在有限长Z平面有极点存在，其结构是递归的，即存在输出到输入的反馈。由于IIR滤波器能够保留一些模拟滤波器的优良特性，因此应用很广。但是这些特性是以牺牲线性相位频率特性为代价的，即用巴特沃斯，切比雪夫，椭圆设计的数字滤波器逼近理想的滤波器的幅度频率特性，得到的滤波器往往是非线性的。在许多电子系统中，对幅度频率特性和线性相位特性都有较高的要求，所以IIR在这些滤波器中往往难以胜任。有限长单位冲激响应（FIR）数字滤波器具有以下特点： 1）可在设计任意幅度频率特性滤波器的同时，保证精确严格的线性相位特性 2）FIR数字滤波器的单位冲激响应h(n)是有限长的，可以用一个固定的系统来实现，因而FIR数字滤波器可以做成因果稳定系统 3）允许做成多通带系统； 有以上原因选择FIR数字滤波器进行设计； （3）窗函数的选择 设计FIR数字滤波器最简单的方法是用窗函数法，设定通带允许的最大衰减为-2dB,阻带必须达到的最小衰减为-21dB,由通带允许的最大衰减并结合窗函数的阻带最小衰减的特性表，可选择矩形窗，其过渡带为; （4）为了改善滤波器的性能， 通常要求窗函数具有一些好的特性， 即 　　（1） 主瓣宽度窄， 以获得较陡的过渡带。  　　（2） 最大旁瓣相对主瓣值（简称为旁瓣电平）尽可能小， 以改善通带的平稳度和增大阻带中的衰减。  　　但是，这两者之间常常有矛盾，难以同时满足。常用的窗函数在这两个因素之间取得适度的折衷时，往往需要增加主瓣宽度来换取旁瓣的抑制。如果选用一个窗函数的主要目的是为了得到较窄的过渡带， 就选用主瓣较窄的窗函数， 但通带和阻带内的振荡将会较大。相反，如果主要目的是为了得到平坦的幅度响应和较小的阻带滤纹，这时选用的窗函数的旁瓣电平就要小，但所设计的FIR滤波器的过渡带就可能较宽。 （5）选择低通滤波器进行设计，因为低通滤波器对语音信号无失真恢复； 2.3.2 设计滤波器 （1）根据上一步对滤波器进行设计； 1）根据2.3.1的（1）以及矩形窗的过渡带公式可得：滤波器的阶数，故取N=12； 2）滤波器的频率响应 数字滤波器的频率响应为Hd(ejω)，hd(n)是与其对应的单位采样响应。但在一般情况下，Hd(ejω)通常是逐段恒定的，在边界频率处有不连