matlab声音信号的采集与滤波处理.doc

MATLAB实验报告(题目二)(题目二)声音信号的采集与滤波处理采用IIR滤波器或FIR滤波器参考资料信号的采集、数字信号处理及滤波实例要求(1)采集声音信号或打开已录好的声音文件，并显示其信号图与频域图。(2)根据信号的特点，选用合适的滤波器，给定滤波器的规一化性能指标(参考指标,实际中依据每个同学所叠加噪声情况而定)例如通带截止频率WP025PI,阻通带截止频率WS03PI通带最大衰减RP1DB阻带最小衰减RS15DB，对信号进行滤波。在MATLAB中，可以利用函数FIR1设计FIR滤波器，可以利用函数BUTTE，CHEBY1和ELLIP设计IIR滤波器；利用MATLAB中的函数FREQZ画出各滤波器的频率响应，滤波器设计完后，用FILTER函数用这些数字滤波器对含噪语音信号分别进行滤波处理。(3)还原音乐信号，并画出其时域图与频域图，并与原始信号比较，且回放音乐信号。(1)打开一个自己录制的音乐文件进行实验，这是实验程序FS22050语音信号采样频率为22050X1WAVREAD E\威尼斯的泪WAV 读取语音信号的数据，赋给变量X1SOUNDX1,22050播放语音信号Y1FFTX1,1024对信号做1024点FFT变换FFS0511/1024FIGURE1PLOTX1做原始语音信号的时域图形TITLE 原始语音信号 XLABEL TIMEN YLABEL FUZHIN FIGURE2FREQZX1绘制原始语音信号的频率响应图TITLE 频率响应图 FIGURE3SUBPLOT2,1,1PLOTABSY11512做原始语音信号的FFT频谱图TITLE 原始语音信号FFT频谱 SUBPLOT2,1,2PLOTF,ABSY11512TITLE 原始语音信号频谱 XLABEL HZ YLABEL FUZHI 实验效果(2)实验程序CLEARFS22050X1WAVREAD E\威尼斯的泪WAV FFS0511/1024T01/22050LENGTHX11/22050将所加噪声信号的点数调整到与原始信号相同AU1D05COS2PI1000T 噪声为1KHZ的余弦信号X2X1DSOUNDX1,8000PAUSE50SOUNDX2,22050播放加噪声后的语音信号Y2FFTX2,1024FIGURE1PLOTT,X2TITLE 加噪后的信号 XLABEL TIMEN YLABEL FUZHIN FIGURE2SUBPLOT2,1,1PLOTF,ABSX11512TITLE 原始语音信号频谱 XLABEL HZ YLABEL FUZHI SUBPLOT2,1,2PLOTF,ABSX21512TITLE 加噪后的信号频谱 XLABEL HZ YLABEL FUZHI 实验效果(3)实验程序FS22050X1WAVREAD E\威尼斯的泪WAV T01/22050LENGTHX11/22050AU05DAUCOS2PI8000T X2X1DWP025PIWS03PIRP1RS15FS22050TS1/FSWP12/TSTANWP/2将模拟指标转换成数字指标WS12/TSTANWS/2N,WNBUTTORDWP1,WS1,RP,RS, S 选择滤波器的最小阶数Z,P,KBUTTAPN创建BUTTERWORTH模拟滤波器BAP,AAPZP2TFZ,P,KB,ALP2LPBAP,AAP,WNBZ,AZBILINEARB,A,FS用双线性变换法实现模拟滤波器到数字滤波器的转换H,WFREQZBZ,AZ绘制频率响应曲线FIGURE1PLOTWFS/2PI,ABSHGRIDXLABEL 频率／HZ YLABEL 频率响应幅度 TITLE BUTTERWORTH F1FILTERBZ,AZ,X2FIGURE2SUBPLOT2,1,1PLOTT,X2画出滤波前的时域图TITLE 滤波前的时域波形 SUBPLOT2,1,2PLOTT,F1画出滤波后的时域图TITLE 滤波后的时域波形 SOUNDF1,22050播放滤波后的信号F0FFTF1,1024FFS0511/1024FIGURE3Y2FFTX2,1024SUBPLOT2,1,1PLOTF,ABSY21512画出滤波前的频谱图TITLE 滤波前的频谱 XLABEL HZ YLABEL FUZHI SUBPLOT2,1,2F1PLOTF,ABSF01512画出滤波后的频谱图TITLE 滤波后的频谱 XLABEL HZ YLABEL FUZHI 实验结果