合肥工业大学宣城校区信号与处理实验解剖.docVIP

合肥工业大学 实验报告 自动化 专业 2 班级 学号201221**** 姓名 大方 日期 指导教师 共 页 第 页 成绩 实验一 含噪声语音信号频谱分析及滤波 一、实验原理 1.1数字滤波器 在工程实际中遇到的信号经常伴有噪声。为了消除或减弱噪声，提取有用的信号，必须进行滤波。数字滤波器，是数字信号处理中极其重要的一部分。随着信息时代和数字技术的发展，收到人们越来越多的重视。数字滤波器可以通过数值运算实现滤波，所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。数字滤波器的种类很多，根据其实现的网络结构或者其冲激响应函数的时域特性，可分为两种，即有限冲激响应（FIR）滤波器和无限冲激响应（IIR）滤波器。 1.2 采样定理 采样定理又称奈奎斯特定理，在进行模拟/数字信号的转换过程中，当采样频率fs不小于信号中最高频率fm的2倍时，采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍。 1.3 FIR数字滤波器设计原理 FIR滤波器结构上主要是非递归结构，没有输出到输入的反馈，系统函数H(z)处于收敛，极点全在z=0处（因果系统），因而只能用较高的阶数达到高的选择性。FIR数字滤波器的幅频特性精度较之于IIR数字滤波器低，但是线性相位，就是不同频率分量的信号经过FIR滤波器后他们的时间差不变，这是很好的性质。FIR滤波器是有限的单位响应也有利于对数字信号的处理，便于编程，用于计算的时延也小，这对实时的信号处理很重要。FIR滤波器因具有系统稳定，以实现相位控制，允许设计多通带滤波器等优点收到人们的青睐。基于窗函数的FIR 数字滤波器的设计方法通常也称之为傅立叶级数法，是用一定宽度窗函数截取无限脉冲响应序列，获得有限长的脉冲响应序列，从而得到FIR 滤波器。它是在时域进行的，由理想滤波器的频率响应推导出其单位冲激响应hd (n)，再设计一个FIR数字滤波器的单位冲激响应h(n)去逼近hd (n)表示，由此得到的离散滤波器的系统传递函数Hd (z)，该hd (n) 为无限长序列，因此Hd (z)是物理不可实现的。为了使系统变为物理可实现的，且使实际的FIR 滤波器频率响应尽可能逼近理想滤波器的频率响应，采用窗函数将无限脉冲响应hd (n)截取一段h(n)来近似表示hd (n)，可得：h (n) = hd (n)w(n) ,从而有：式中N 表示窗口长度，这样H(z)就是物理可实现的系统。并且从线性相位FIR滤波器的充要条件可知，为了获得线性相位FIR 数字滤波器的冲激响应h(n) ，那么序列h(n) 应有τ = (N ?1) / 2的延迟。由于窗函数的选择对结果起着重要的作用，针对不同的信号和不同的处理目的来确定窗函数的选择才能收到良好的效果。 1.4 IIR数字滤波器设计原理 IIR滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由时延、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。同时，IIR滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、切比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，在设计一个IIR滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。利用双线性变换设计IIR滤波器（巴特沃斯数字低通滤波器的设计），首先要设计出满足指标要求的模拟滤波器的传递函数Ha(s),然后由Ha(s)通过双线性变换可得所要设计的IIR滤波器的系统函数H(z)。如果给定的指标为数字滤波器的指标，则首先要转换成模拟滤波器的技术指标，这里主要是边界频率Wp和Ws的转换，对ap和as指标不作变化。边界频率的转换关系为 ∩=2/T tan(w/2)。接着，按照模拟低通滤波器的技术指标根据相应设计公式求出滤波器的阶数N和3dB截止频率 ∩c ；根据阶数N查巴特沃斯归一化低通滤波器参数表，得到归一化传输函数Ha(p)；最后，将p=s/ ∩c 代入Ha(p)去归一，得到实际的模拟滤波器传输函数Ha(s)。之后，通过双线性变换法转换公式s=2/T((1-1/z)/(1+1/z))得到所要设计的IIR滤波器的系统函数H(z)。 实验内容 录制一段个人自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；并在语音信号中自行加入高斯白噪声信号，对加入噪声信号后的语音信号进行频谱分析；给定滤波器的性能