基于MATLAB的语音信号去噪及实现.docVIP

《数字信号处理》 课程设计报告 基于MATLAB的语音信号去噪及实现 专 业： 通信工程 班 级： 通信08-1BF 组 次： 第6组 姓 名： 王 丹 学 号： 14082302499 语音信号的数字滤波-基于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仿真，听滤波前后的