语音信号数字滤波基于MATLAB的语音信号去噪及实现.docVIP

《数字信号处理》 课程设计报告 基于MATLAB的语音信号去噪及实现 专 业： 通信工程 班 级： 通信08-1BF 语音信号的数字滤波-基于MATLAB的语音信号去噪及实现 摘要 本课程设计主要运用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用一个满足指标的，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析，根据结果和学过的理论得出合理的结论。MATLAB 1 设计目的 本课程设计主要利用麦克风采集一段语音信号， 图1 脉冲响应不变法巴特沃斯滤波器对语音信号去噪流程图 3 技术实现 3.1语言信号的采集 图3是基于PC机的语音信号采集过程:通过计算机录音系统录制1秒： 图2 语音采集过程 图3 录音过程截图 3.2语音信号的频谱分析 3.3滤波器设计 将数字滤波器的设计指标设为通带截止频率fb=1100HZ,阻带频率fc=1200HZ,通带波纹Ap=1dB，阻带波纹As=20dB,要求确定H(z)。设计程序如下： Fp=1200;Ft=8000;As=20;Ap=1; wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft; fp=2*Ft*tan(wp/2);fs=2*Ft*tan(ws/2); [n,wn]=buttord(wp,ws,Ap,As,s); [b,a]=butter(n,wn,s); [num,den]=bilinear(b,a,1); [h,w]=freqz(num,den); Buttor波形如下图5： 图5 IIR低通滤波器 3.4 信号滤波过程 [y,fs,nbits]=wavread(3333.wav);%载入波形 N=length(y)%计算输入波形长度 Y=fft(y,N);%进行快速傅里叶变换 sound(y);%播放声音 figure（4）；%时域图 figure（5）；plot（abs（Y））；%滤波前频谱图 Fp=1200;%阻带截止频率 Fs=1100;%通带截止频率 Ft=8000;%采集频率 As=20;%通带波纹Ap=1 Ap=1;%阻带波纹As=20 wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft; fp=2*Ft*tan(wp/2);fs=2*Ft*tan(ws/2); [n,wn]=buttord(wp,ws,Ap,As,s); %求低通滤波器的阶数和截止频率 [b,a]=butter(n,wn,s); %求S域的频率响应的参数 [num,den]=bilinear(b,a,1); %双线性变换实现S域到Z域的变换 [h,w]=freqz(num,den); %根据参数求出频率响应 figure(3) plot(w*8000*0.5/pi,abs(h));z=filter(num,den,y); sound(z); m=z; %求滤波后的信号 figure(1) subplot(2,2,3);plot(abs(m),r);title(滤波后信号的频谱); grid; subplot(2,2,4);plot(z,b);title(滤波后的信号波形); grid;subplot(2,2,2); plot(y,b);title(滤波前信号的波形); grid;subplot(2,2,1); plot(abs(Y),r);title(滤波前信号的频谱); grid; figure(2);p=angle(m);q=angle(Y); subplot(2,1,1);plot(q,b);title(滤波前相位); grid;subplot(2,1,2);plot(p,b);title(滤波后相位);grid; 程序运行结果如下所示： 图6.1 滤波前后的时域对比波形和滤波前后频域频谱图 图6.2滤波前后相位比较图 3.5结果分析 由图6.1中滤波前后波形比较可看出，经过滤波后的波形比原波形的振幅有所减小，去除了很多由于噪声所产生的干扰；从滤波前后的频谱比较可以看出经过滤波后除了原本的声音外，中间由于噪声产生的频谱波形已经滤除；由图6.2滤波前后相位比较图可看出由于经过滤波，相位变得稀疏；经过MATLAB仿真，听滤波前后的声音，可以听出有明显的滤波效果。因此利用脉冲响应不变法设计的巴特沃斯滤波器已经达到了设计的要求。