《信号处理消噪与播音报告》.doc

数字信号处理应用 实验报告 噪音信号处理  HYPERLINK F:\处理中心\已整理\sound\sound.wav语音信号 使用说明  HYPERLINK F:\处理中心\已整理\sound\sound.m源代码 输出图像  HYPERLINK F:\处理中心\已整理\sound\f1.jpgFigure.1 原始语音信号，噪声信号，加噪后的语音信号  HYPERLINK F:\处理中心\已整理\sound\f2.jpgFigure.2原始语音信号频谱， 噪声信号频谱，加噪后的语音信号频谱  HYPERLINK F:\处理中心\已整理\sound\f3.jpgFigure.3 滤波器脉冲响应  HYPERLINK F:\处理中心\已整理\sound\f4.jpgFigure.4滤波前的时域波形，滤波后时域波形  HYPERLINK F:\处理中心\已整理\sound\fig.jpgFigure.5 滤波前频谱，滤波后频谱 二．电话拨号音的合成与识别  HYPERLINK F:\处理中心\已整理\telephone.m(链接到telephone.m文件) 一、噪音信号处理 1、语音信号 2、使用说明 将sound文件夹放到D:\MATLAB7\work\目录下，sound文件夹里面包含sound.m 文件和sound.wav文件。 打开MATLAB，在里面OPEN FILE 打开sound.m 文件编译即可得到输出图像。 3、源代码 %%%%%对语音信号加噪声，并滤波输出%%%%%%% %%%%% 1. 对语音信号加噪声 %%%%%%%%%%% clear,clc; [x1,fs,bits]=wavread(D:\MATLAB7\work\sound\sound.wav); n=length(x1); %取读入声音信号的序列长度 f=(0:n-1)*fs/n; n1=(0:n-1); niose=0.1*sin(2*pi*3000*n1/fs);%噪声信号 x2=niose+x1; t=1/fs:1/fs:n/fs; y1=fft(x1,n);%对信号做FFT变换 y=fft(niose,n); y2=fft(x2,n); figure(1) subplot(3,1,1); plot(t,x1); title(Figure.1 原始语音信号); xlabel(t );ylabel(A ); subplot(3,1,2); plot(t,niose); title(Figure.2 噪声信号); xlabel(t );ylabel(A ); subplot(3,1,3); plot(t,x2);title(Figure.3 加噪后的语音信号); xlabel(t );ylabel(A ); figure(2) subplot(3,1,1); plot(f,abs(y1));title(Figure.1 原始语音信号频谱); xlabel(Hz); ylabel(A); subplot(3,1,2); plot(f,abs(y));title(Figure.2 噪声信号频谱); xlabel(Hz); ylabel(A); subplot(3,1,3); plot(f,abs(y2));title(Figure.3 加噪后的语音信号频谱); xlabel(Hz); ylabel(A); %%%%%%%2. 滤波输出%%%%%%%%% [b,a] = butter(2,[2500/(fs/2) 4000/(fs/2)],stop); %得到滤波器阶数为2*2，频率范围为Wn = [W1,W2]的带阻滤波器; [H,w]=freqz(b,a); figure(3);plot(w,abs(H));axis([0 3 0 2]); ylabel(H（w）); title(滤波器脉冲响应); new=filter(b,a,x2); ynew=fft(new); y3=fftshift(ynew); figure(4); subplot(2,1,1);plot(t,x2); title( Figure.1 滤波前的时域波形); subplot(2,1,2);plot(new); title( Figure.2 滤波后时域波形); figure(5); subplot(2,1,1);plot(f,abs(y2));title( Figure.1 滤波前频谱); subplot(2,1,2);plot(f,abs(y3));title( Figure.2 滤波后频谱); %sound(x1,fs);%播放原始语音信号 %sound(x2,fs);%播放加噪声后的语音信号 %sound(