课程四设计_基于IIR的有噪声语音信号处理文章.doc

《基于IR的有噪声语音信号处理 第1章 绪论 1 第2章 基于滤波器的语音信号的处理 2 2.1语音信号的采集 2 2.2设计IIR数字滤波器 2 2.3用滤波器对加噪语音信号进行滤波.............................7 第3章 仿真及其结果分析 7 3.1 语音信号的时频分析 7 3.2 加噪后的语音信号及其频谱分析 9 3.3 比较滤波前后语音信号的波形及频谱.........................10 3.4 IIR滤波器处理结果的的分析 16 第二章：基于滤波器的语音信号的处理 选择一个语音信号作为分析的对象，或录制一段语音信号，对其进行频谱分析；利用MATLAB中的随机函数产生噪声加入到语音信号中，模仿语音信号被污染，并对其进行频谱分析；设计IIR数字滤波器， 并对被噪声污染的语音信号进行滤波， 分析滤波后信号的时域和频域特征，最后回放语音信号。 2.1 语音信号的采集 利用WINDOWS 操作系统可以进行数字信号的采集。将话筒输入计算机的语音输入插口上,启动录音机。按下录音按钮,接着对话筒说话“语音信号处理”,说完后停止录音,屏幕左侧将显示所录声音的长度。点击放音按钮,可以实现所录音的重现。以文件名“speech”保存入g :\ MATLAB \ work 中。可以看到,文件存储器的后缀默认为. wav ,这是WINDOWS 操作系统规定的声音文件存的标准。 1、原始信号的采集及分析 源程序： [y1,fs,bits]=wavread(D:\MATLAB\ai2.wav); sound(y1,44000); %播放语音信号 y=fft(y1); %对信号做FFT变换 f=fs*(0:511)/1024; figure(1) subplot(2,1,1); plot(abs(y(1:512))) %做原始语音信号的FFT频谱图 title(原始语音信号FFT频谱) subplot(2,1,2);%做原始语音信号的时域图形 plot(y1) title(原始语音信号波形); xlabel(时间 n); ylabel(幅值 n); IIR数字滤波器 IIR滤波器设计方法有间接法和直接法，间接法是借助于模拟滤波器的设计方法进行的。其设计步骤是：先设计过渡模拟滤波器得到系统函数H（s），然后将H（s）按某种方法转换成数字滤波器的系统函数H(z)。FIR滤波器比鞥采用间接法，常用的方法有窗函数法、频率采样发和切比雪夫等波纹逼近法。对于线性相位滤波器，经常采用FIR滤波器。 对于数字高通、带通滤波器的设计，通用方法为双线性变换法。可以借助于模拟滤波器的频率转换设计一个所需类型的过渡模拟滤波器，再经过双线性变换将其转换策划那个所需的数字滤波器。具体设计步骤如下： （1）确定所需类型数字滤波器的技术指标。 （2）将所需类型数字滤波器的边界频率转换成相应的模拟滤波器的边界频率，转换公式为Ω=2/T tan(0.5ω) (3)将相应类型的模拟滤波器技术指标转换成模拟低通滤波器技术指标。 （4）设计模拟低通滤波器。 （5）通过频率变换将模拟低通转换成相应类型的过渡模拟滤波器。 （6）采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。 我们知道，脉冲响应不变法的主要缺点是会产生频谱混叠现象，使数字滤波器的频响偏离模拟滤波器的频响特性。为了克服之一缺点，可以采用双线性变换法。 下面我们总结一下利用模拟滤波器设计IIR数字低通滤波器的步骤： （1）确定数字低通滤波器的技术指标：通带边界频率、通带最大衰减，阻带截止频率、阻带最小衰减。 （2）将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。 （3）按照模拟低通滤波器的技术指标设计及过渡模拟低通滤波器。 （4）用双线性变换法，模拟滤波器系统函数转换成数字低通滤波器系统函数。 如前所述，IIR滤波器和FIR滤波器的设计方法有很大的区别。下面我们着重介绍用窗函数法设计FIR滤波器的步骤。如下： （1）根据对阻带衰减及过渡带的指标要求，选择串窗数类型（矩形窗、三角窗、汉宁窗、哈明窗、凯塞窗等），并估计窗口长度N。先按照阻带衰减选择窗函数类型。原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣的窗函数。 （2）构造希望逼近的频率响应函数。 （3）计算h(n).。 （4）加窗得到设计结果。 接下来，我们根据语音信号的特点给出有关滤波器的技术指标： ①低通滤波器的性能指标： fp=1000Hz，fc=1200Hz，As=100db ,Ap=1dB ②高通滤波器的性能指标： fp=3500Hz，fc=4000Hz，As=100dB，Ap=1dB；