语音信号的数字滤波处理

语音信号的数字滤波处理——巴特沃思、bartlett窗滤波器设计

常用的窗函数有以下几种:矩形窗(Rectangular window)、三角窗(Triangular window)、汉宁窗(Hanning window)、海明窗(Hamming window)、布拉克曼窗(Bartlett window)、切比雪夫窗(Chebyshev window)、巴特里特窗(Bartlett window)及凯塞窗(Kaiser window)。本实验中选用布拉克曼窗(Bartlett window)设计滤波器。

2.2.1.1 Bartlett窗低通滤波器

(1)选择来自于D盘的“ding.wav”声音作为语音信号(用如下语句调用[x,FS,bits]=wavread('D:\\ding.wav'))。

(2)给信号加一个大频率的噪声(取噪声频率远大于语音信号的最大频率),产生污染信号。

(3) 设计一个Bartlett窗低通滤波器,通带范围包括语音信号,阻带频率设定为小于噪声信号频率。

(4)将设计好的Bartlett窗低通滤波器滤除被噪声污染后的语音信号。还原语音信号。

2.2.1.2 Bartlett窗高通滤波器

(1)选择来自于D盘的“ding.wav”声音作为语音信号(用如下语句调用[x,FS,bits]=wavread('D:\\ding.wav'))。

(2)给信号加一个小频率的噪声(取噪声频率远小于语音信号的最小频率),产生污染信号。

(3)设计一个Bartlett窗高通滤波器,通带范围包括语音信号,阻带频率设定为大于噪声信号频率。

(4)将设计好的Bartlett窗高通滤波器滤除被噪声污染后的语音信号。还原语音信号。

2.2.1.3 Bartlett窗带通滤波器

(1)选择来自于D盘的“ding.wav”声音作为语音信号(用如下语句调用[x,FS,bits]=wavread('D:\\ding.wav'))。

(2)给信号加一个小频率或大频率的噪声(取噪声频率远小于语音信号的最小频率或大于语音信号的最大频率),产生污染信号。本实验取小频率的噪声信号。

第 5 页 共 56 页

语音信号的数字滤波处理——巴特沃思、bartlett窗滤波器设计

(3)设计一个Bartlett窗带通滤波器,通带范围包括语音信号,阻带频率设定为不包括噪声信号频率。

(4)将设计好的Bartlett窗带通滤波器滤除被噪声污染后的语音信号。还原语音信号。

3 设计步骤和过程

3.1 设计步骤 (1)语音信号的采集

要求利用windows下的录音机(开始—程序—附件—娱乐—录音机,文件—属性—立即转换—8000Hz,8位,单声道)录制一段自己的话音,或者采用Windows自带的声音文件(默认为22050Hz),时间控制在几秒左右。然后在MATLAB软件平台下,利用函数wavread对语音信号进行采样,记住采样频率和采样点数。通过wavread函数的使用,要求理解采样频率、采样位数等概念。 wavread函数调用格式:

y=wavread(file),读取file所规定的wav文件,返回采样值放在向量y中。 [y,fs,nbits]=wavread(file),采样值放在向量y中,fs表示采样频率(Hz),nbits表示采样位数。

y=wavread(file,N),读取前N点的采样值放在向量y中。

y=wavread(file,[N1,N2]),读取从N1点到N2点的采样值放在向量y中。 (2)语音信号的频谱分析

要求首先画出语音信号的时域波形;然后对语音信号进行频谱分析,在MATLAB中,可以利用函数fft对信号进行快速付立叶变换,得到信号的频谱特性;从而加深对频谱特性的理解。 (3)设计数字滤波器和画出频率响应

根据语音信号的特点给出有关滤波器的性能指标,例如:1)低通滤波器性能指标,fp=1000Hz,fc=1200 Hz, As=100dB,Ap=1dB;2)高通滤波器性能指标,fc=2800 Hz,fp=3000 Hz,As=100dB,Ap=1dB;3)带通滤波器性能指标,fp1=1200 Hz,fp2=3000 Hz,fc1=1000 Hz,fc2=3200 Hz,As=100dB,Ap=1dB。

第 6 页 共 56 页

语音信号的数字滤波处理——巴特沃思、bartlett窗滤波器设计

(4)回放语音信号

在MATLAB中,函数sound可以对声音进行回放。其调用格式:sound(x,fs,bits);可以感觉滤波前后的声音有变化。这里可以通过两种方式进行信号采集,一是在自己电脑C:\\WINDOWS\\Media中有很多系统声音,我选择了自己电脑里的“Windows XP 叮当声”这个声音文件(.wav)。二用电脑自带录音软件(开始\\所有程序\\附件\\娱乐\\录音机),或者使用其它专业的录音软件,录制时需要配备录音硬件(如麦克风),为便于比较,需要在安静、干扰小的环境下录音。

3.2设计详细过程

3.2.1语音信号的采集与频谱分析

将电脑D:\\ding.wav这一声音文件导入MATLAB软件中,绘出采样后的语音信号的时域波形和频谱图。

MATLAB程序如下: clf;

[x,FS,bits]=wavread('C:\\WINDOWS\\Media'); x=x(:,1); figure(1);

subplot(211);plot(x); %sound(x,FS,bits);

title('语音信号时域波形图') y=fft(x,3260);

f=(FS/1630)*[1:1630]; subplot(212);

plot(f(1:1630),abs(y(1:1630))); title('语音信号频谱图');

运行程序得到语音信号时域和频谱图形如图3.1所示。

第 7 页 共 56 页

语音信号的数字滤波处理——巴特沃思、bartlett窗滤波器设计

图3.1 语音信号时域和频谱图

3.2.2 IIR数字滤波器的设计

下面分别用巴特沃思数字滤波器(低通、带通、高通)对被噪声污染后的信号进行滤波。

3.2.2.1巴特沃思低通滤波器

(1)添加高频的噪声信号zs=0.05*cos(2*pi*10000*t/22050),MATLAB程序如下。噪声信号的时域及幅频图如图3.2所示。 %噪声信号

zt=0:length(x)-1;

zs=0.05*cos(2*pi*10000*t/22050); zs0=0.05*cos(2*pi*10000*t/22050000); figure(2); subplot(2,1,1) plot(zs0)

title('噪声信号波形'); zs1=fft(zs,1200);

%sound(zs,FS,bits); %回放噪音 subplot(2,1,2)

plot(f(1:600),abs(zs1(1:600))); title('噪声信号频谱');

第 8 页 共 56 页

联系客服:779662525#qq.com(#替换为@)