语音信号滤波去噪――使用双线性变换法设计的切比雪夫II型滤波器

目录一、引言.....................................................31.1课程设计目的..........................................31.2课程设计要求..........................................3二、设计原理.................................................42.1IIR滤波器.............................................42.2切比雪夫I型滤波器....................................42.3双线性变换法..........................................7三、设计步骤.................................................83.1设计流程图............................................83.2语音信号的采集........................................93.3语音信号的频谱分析...................................103.4滤波器设计...........................................123.5完整的滤波程序及滤波效果图...........................153.6结果分析.............................................18四、出现的问题及解决方法....................................20五、课程设计心得体会........................................21六、参考文献................................................22成绩2语音信号滤波去噪——使用双线性变换法设计的切比雪夫I型滤波器摘要随着信息技术的发展，现代信号处理正向着数字化、软件化方向发展。滤波器设计是信号处理的重要组成部分，而研究语音信号的滤波设计是现代信息处理的基本内容。本设计利用计算机WINDOWS下的录音机录入一句语音信号，用MATLAB软件对其进行频谱分析，然后加入一干扰信号，利用设计好的滤波器将干扰信号去除，最后对各部分的频谱进行分析比较。关键词滤波设计；MATLAB；一、引言用麦克风采集一段8000Hz，8k的单声道语音信号，绘制波形并观察其频谱，给定通带截止频率为2000Hz，阻带截止频率为2100Hz，通带波纹为1dB，阻带波纹为60dB，用双线性变换法设计的一个满足上述指标的切比雪夫II型IIR滤波器，对该语音信号进行滤波去噪处理。1.1课程设计目的《数字信号处理》课程设计是在学生完成数字信号处理和MATLAB的结合后的基本实验以后开设的。本课程设计的目的是为了让学生综合数字信号处理和MATLAB并实现一个较为完整的小型滤波系统。这一点与验证性的基本实验有本质性的区别。开设课程设计环节的主要目的是通过系统设计、软件仿真、程序安排与调试、写实习报告等步骤，使学生初步掌握工程设计的具体步骤和方法，提高分析问题和解决问题的能力，提高实际应用水平。1.2课程设计的要求（1）学会MATLAB的使用，掌握MATLAB的程序设计方法；3（2）滤波器指标必须符合工程实际，根据模拟滤波器的性能指标，确定数字滤波器指标；（3）采用双线性变换法，设计满足上述性能指标要求的ChebyshevII型数字低通滤波器；（4）设计完后应检查其频率响应曲线是否满足指标；（5）处理结果和分析结论应该一致，而且应符合理论；（6）独立完成课程设计并按要求编写课程设计报告书；二、设计原理用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析。2.1IIR滤波器从离散时间来看，若系统的单位抽样(冲激)响应延伸到无穷长，称之为“无限长单位冲激响应系统”，简称为IIR系统。无限长单位冲激响应（IIR）滤波器有以下几个特点:(1)系统的单位冲激响应h(n)是无限长；(2)系统函数H(z)在有限z平面（0z∞）；(3)结构上存在着输出到输入的反馈，也就是结构上是递归型的。IIR滤波器采用递归型结构，即结构上带有反馈环路。同一种系统函数H（z）可以有多种不同的结构，基本网络结构有直接Ⅰ型、直接Ⅱ型、级联型、并联型四种，都具有反馈回路。同时，IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等，这些典型的滤波器各有特点。有现成的设计数据或图表可查，在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。42.2切比雪夫I型滤器2.2.1切比雪夫滤波器简介切比雪夫滤波器（又译车比雪夫滤波器）是在通带或阻带上频率响应幅度等波纹波动的滤波器。在通带波动的为“I型切比雪夫滤波器”，在阻带波动的为“II型切比雪夫滤波器”。切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快，但频率响应的幅频特性不如后者平坦。切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动。这种滤波器来自切比雪夫多项式，因此得名，用以纪念俄罗斯数学家巴夫尼提·列波维其·切比雪夫2.2.2切比雪夫滤波器原理巴特沃兹滤波器在通带内幅度特性是单调下降的，如果阶次一定，则在靠近截止处，幅度下降很多，或者说，为了使通带内的衰减足够小，需要的阶次N很高，为了克服这一缺点，采用切比雪夫多项式来逼近所希望的。切比雪夫滤波器的在通带范围内是等幅起伏的，所以在同样的通常内衰减要求下，其阶数较巴特沃兹滤波器要小。切比雪夫滤波器的振幅平方函数为（1）式中Ωc—有效通带截止频率—与通带波纹有关的参量，大,波纹大01VN（x）—N阶切比雪夫多项式5（2）|x|≤1时,|VN(x)|≤1|x|1时,|x|↗,VN(x)↗切比雪夫滤波器的振幅平方特性如图所示,通带内，的变化范围为1(max)→211(min)时，|x|1,随↗，→0(迅速趋于零)当=0时，（3）N为偶数，cos2()=1，得到min，，（4）N为奇数，cos2(，得到max，（5）切比雪夫滤波器的振幅平方特性如图1所示。6图1切比雪夫滤波器的振幅平方特性2.3双线性变换法双线性变换法是使数字信号滤波器的频率响应与模拟滤波器的频率响应相似的一种变换方法。为了客服多值映射这一缺点，我们首先把整个s平面压缩变换到某一中介的s1平面的一条横带里（宽度为T2,即从T到T),其次再通过上面讨论过的标准变换关系Tsez1将此横带变换到这个z平面上去，这样就使s平面与z平面式一一对应的关系，消除了多值变换性，也就消除了频谱混叠现象。将s平面整个j轴压缩变换到s1平面j轴上的T到T一段，可以采用以下变换关系：)2tan(1T（6）这样，变换到T1，0变到01可将（6）式写成22221111TjTjTjTjeeeej（7）解析延拓到整个s平面和s1平面，令sj，11sj，则得TsTsTjTjTjTjeeTstheeeej11111111212222（8）再将s1平面通过以下标准变换关系映射到z平面：Tsez1（9）从而得到s平面和z平面的单值映射的关系为1111zzs（10）ssz11（11）一般来说，为了使模拟滤波器的某一频率与数字滤波器的任一频率有对应的关系，可以引入待定常数c，使（6）式和（7）式变换成7)2tan(1Tc（12）TsTseecTscths111121（13）仍将Tsez1代入（13）式，可得1111zzcs（14）scscz（15）（14）式和（15）式是s平面与z平面之间的单值映射关系，这种变换就称为双线性变换。三、设计步骤3.1设计流程图语音信号滤波去噪——使用双线性变换法设计的切比雪夫I型滤波器的设计流程如图3.1.1所示：语音信号的采集（wavread函数），画时域图快速傅里叶变换，并且画频谱图设定滤波器性能指标，通带截止频率fb=2000，阻带截止频率fc=2100，通带波纹Ap=1，阻带波纹As=60双线性变换法设计切比雪夫I型滤波器验证并进行频谱分析设计好的滤波器进行滤波处理比较滤波前后语音信号的波形及频谱回放语音信号结束开始8双线性变换法切比雪夫II型滤波器对语音信号去噪流程图3.2语言信号的采集点击windows系统桌面的“开始”按钮，点击开始菜单栏里的“附件”，选择“录音机”选项，点击录音机“文件”选项，进入“声音选定”设置，把属性一栏设置成“8000Hz，8位，单声道，7KB/秒”(见图3.2.1)。点击确定，然后开始语言信号的采集，采集时间为1秒左右为最佳。采集的声音文件以“.wav”格式存储（见图3.2.2）。图3.2.1采集声音的参数设置图3.2.2采集声音93.3语音信号的频谱分析在MATLAB中编辑m函数，使用wavread函数读取采集的声音文件（.wav)将它赋值给某一向量，再对其进行采样,然后使用plot语句画出相关的频谱图形在figure(1）上。（1）Wavread函数调用格式：[y,Fs,nbits]=wavread(file)功能说明：采样值放在向量y中，Fs表示采样频率（Hz），nbits表示采样位数。（2）快速傅里叶变换算法FFT计算DFT的函数fft,其调用格式如下：Xk=fft(x,n)参数x为被变换的时域序列向量，N是DFT变换区间长度，当n大于x的长度时，fft函数自动在x后面补零。，当n小于xn的长度时，fft函数计算x的前n个元素，忽略其后面的元素。在本次课程设计中，我们利用fft函数对语音信号进行快速傅里叶变换，就可以得到信号的频谱特性。（3）声音采样文件读取的程序（文件名：123.wav)%用麦克风采集一段8000Hz，8k的单声道语音信号，绘制波形并观察其频谱Fs=8000;%语音信号采样频率为8000x1=wavread('123.wav');%读取语音信号的数据，赋给变量x1sound(x1,8000);%播放语音信号y1=fft(x1,8000);%对信号做FFT变换subplot(1,2,1);plot(x1)%做原始语音信号的时域图形title('原始语音信号时域图');xlabel('时间n');ylabel('音量n');subplot(1,2,2);10plot(abs(y1(1:8000)))%做原始语音信号的FFT频谱图axis([0,4000,0,10])title('原始语音信号FFT频谱')figurefreqz(x1)%绘制原始语音信号的频率响应图title('原始语音信号频率响应图')所得图形图3.3.1原始语音信号的时域图和FFT频谱11图3.3.2原始语音信号频率响应图3.4滤波器设计设计指标：通带截止频率为1100Hz，阻带截止频率为1200Hz，通带波纹为1dB，阻带波纹为20dB，用脉冲（冲激）响应不变法设计的一个满足上述指标的切比雪夫I型滤波器双线性变换法设计切比雪夫I型滤波器Fs=8000;%采样频率wp=1100*2/Fs;ws=1200*2/Fs;%根据采样频率将滤波器边界进行转化Rp=1;%通带波纹Rs=20;%阻带波纹滤波程序如下：fp=2000;fs=2100;%给