搜索
课程论文基于MATLAB的语音信号的分析和处理研究课程名称语音信号处理学院信息工程学院专业电子信息工程班级电子信息工程信息处理姓名杨歌论文分数评阅教师中国传媒大学2011年7月24日一、研究综述本次研究题目主要是利用MATLAB对语音信号进行分析和处理。具体来说，是在采集语音信号后，在MATLAB软件平台进行频谱分析，并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器滤除噪声，恢复原信号。语音信号处理是语音学与数字信号处理技术相结合的交叉学科，因此在这里将语音当做一种特殊的信号，即一种“复杂向量”来看待。也就是说，题目更多的还是体现了数字信号处理技术。相关的资料检索如下：以上是国内基于MATLAB对语音信号进行处理的一些专利，具体包括电子耳蜗听觉仿真平台、仿真试验方法及算法比较方法、客观音质评价归一化主客观统计相关模型的建立方法、一种基于矢量量化的声纹识别方法，从专利的题目和内容来看，可看到语音处理包含在方方面面，具有很大的潜力，而利用MATLAB对语音信号进行识别、处理等有这广泛的发展前景。基于这些我将我的的研究题目放在很小很基本但必不可少的利用MATLAB对语音信号进行处理识别滤波。二、研究目的本次研究题目的目的主要是利用MATLAB对语音信号进行分析和处理。MATLAB是一个很好的信号处理平台，通过在采集语音信号后，在MATLAB软件平台进行频谱分析，并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，从而设计合适的滤波器滤除噪声，恢复原信号。其中，待处理语音信号是一个在20Hz～20kHz频段的低频信号。三、设计原理1．采样原理在进行模拟/数字信号的转换过程中，当采样频率fs.max大于信号中最高频率fmax的2倍时，即：fs.max=2fmax,则采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定理又称奈奎斯特定理。1924年奈奎斯特(Nyquist)就推导出在理想低通信道的最高大码元传输速率的公式:理想低通信道的最高大码元传输速率=2W*log2N(其中W是理想低通信道的带宽,N是电平强度)2.采样频率采样频率是指计算机每秒钟采集多少个声音样本，是描述声音文件的音质、音调，衡量声卡、声音文件的质量标准。采样频率越高，即采样的间隔时间越短，则在单位时间内计算机得到的声音样本数据就越多，对声音波形的表示也越精确。采样频率与声音频率之间有一定的关系，根据奎斯特理论，只有采样频率高于声音信号最高频率的两倍时，才能把数字信号表示的声音还原成为原来的声音。这就是说采样频率是衡量声卡采集、记录和还原声音文件的质量标准。3.采样位数与采样频率采样位数即采样值或取样值，用来衡量声音波动变化的参数，是指声卡在采集和播放声音文件时所使用数字声音信号的二进制位数。采样频率是指录音设备在一秒钟内对声音信号的采样次数，采样频率越高声音的还原就越真实越自然。采样位数和采样率对于音频接口来说是最为重要的两个指标，也是选择音频接口的两个重要标准。无论采样频率如何，理论上来说采样的位数决定了音频数据最大的力度范围。每增加一个采样位数相当于力度范围增加了6dB。采样位数越多则捕捉到的信号越精确。对于采样率来说你可以想象它类似于一个照相机，44.1kHz意味着音频流进入计算机时计算机每秒会对其拍照达441000次。显然采样率越高，计算机摄取的图片越多，对于原始音频的还原也越加精确。三、研究内容1.程序源代码利用MATLAB对语音信号进行处理的代码及注释如下：%读取语音信号fs=22050;%语音信号采样频率为22050x1=wavread('D:\课程设计\shuzi.wav');%读取语音信号的数据，赋给变量x1sound(x1,22050);%播放语音信号y1=fft(x1,4096);%对信号做4096点FFT变换figure(1)subplot(321);plot(x1);title('原始信号波形');xlabel('timen');ylabel('fuzhin');subplot(322);plot(y1);title('原始信号频谱');xlabel('Hz');ylabel('fuzhi');pause(5);%延时%加入噪声n=length(x1);Y=0.03*randn(n,1);y=real(fft(Y,4096));Y=Y+x1;y=y+y1;sound(Y,22050);figure(2)subplot(323);plot(Y);title('加噪后信号波形');xlabel('timen');ylabel('fuzhin');subplot(324);plot(y);title('加噪后信号频谱');xlabel('Hz');ylabel('fuzhi');pause(5);%延时%滤波器设计wp=0.25*pi;ws=0.3*pi;Rp=1;Rs=15;Fs=22050;Ts=1/Fs;wp1=2/Ts*tan(wp/2);%将模拟指标转换成数字指标ws1=2/Ts*tan(ws/2);[N,Wn]=buttord(wp1,ws1,Rp,Rs,'s');%选择滤波器的最小阶数[Z,P,K]=buttap(N);%创建butterworth模拟滤波器[Bap,Aap]=zp2tf(Z,P,K);[b,a]=lp2lp(Bap,Aap,Wn);[bz,az]=bilinear(b,a,Fs);%用双线性变换法实现模拟滤波器到数字滤波器的转换[H,W]=freqz(bz,az);%绘制频率响应曲线figure(3)plot(W*Fs/(2*pi),abs(H))gridxlabel('频率／Hz')ylabel('频率响应幅度')title('Butterworth')f1=filter(bz,az,Y);figure(4)subplot(211)plot(Y)%画出滤波前的时域图title('滤波前的时域波形');subplot(212)plot(f1);%画出滤波后的时域图title('滤波后的时域波形');sound(f1,22050);%播放滤波后的信号F0=fft(f1,4096);f=fs*(0:511)/4096;figure(5)subplot(211);plot(y);title('加噪后滤波前信号频谱');xlabel('Hz');ylabel('fuzhi');subplot(212)F1=plot(F0);%画出滤波后的频谱图title('滤波后的频谱')xlabel('Hz');ylabel('fuzhi');2.调试分析过程描述I.语音信号的分析及处理方法a.语音的录入与打开在MATLAB中，[y,fs,bits]=wavread('Blip',[N1N2]);用于读取语音，采样值放在向量y中，fs表示采样频率(Hz)，bits表示采样位数。[N1N2]表示读取从N1点到N2点的值（若只有一个N的点则表示读取前N点的采样值）。sound(x,fs,bits);用于对声音的回放。向量y则就代表了一个信号（也即一个复杂的“函数表达式”）也就是说可以像处理一个信号表达式一样处理这个声音信号。b.时域信号的FFT分析FFT即为快速傅氏变换，是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。在MATLAB的信号处理工具箱中函数FFT和IFFT用于快速傅立叶变换和逆变换。函数FFT用于序列快速傅立叶变换，其调用格式为y=fft(x)，其中，x是序列，y是序列的FFT，x可以为一向量或矩阵，若x为一向量，y是x的FFT且和x相同长度；若x为一矩阵，则y是对矩阵的每一列向量进行FFT。如果x长度是2的幂次方，函数fft执行高速基－2FFT算法，否则fft执行一种混合基的离散傅立叶变换算法，计算速度较慢。函数FFT的另一种调用格式为y=fft(x,N)，式中，x，y意义同前，N为正整数。函数执行N点的FFT，若x为向量且长度小于N，则函数将x补零至长度N；若向量x的长度大于N，则函数截短x使之长度为N；若x为矩阵，按相同方法对x进行处理。II.数字滤波器设计原理数字滤波是数字信号分析中最重要的组成部分之一，与模拟滤波相比，它具有精度和稳定性高、系统函数容易改变、灵活性强、便于大规模集成和可实现多维滤波等优点。在信号的过滤、检测和参数的估计等方面，经典数字滤波器是使用最广泛的一种线性系统。数字滤波器的作用是利用离散时间系统的特性对输入信号波形(或频谱)进行加工处理，或者说利用数字方法按预定的要求对信号进行变换。数字滤波器的设计步骤不论是IIR滤波器还是FIR滤波器的设计都包括三个步骤：a.按照实际任务的要求，确定滤波器的性能指标。b.用一个因果、稳定的离散线性时不变系统的系统函数去逼近这一性能指标。根据不同的要求可以用IIR系统函数，也可以用FIR系统函数去逼近。c.利用有限精度算法实现系统函数，包括结构选择、字长选择等。III.IIR滤波器与FIR滤波器的性能比较FIR:FiniteImpulseresponse，有限冲击响应IIR:InfiniteImpulseresponse，无限冲击响应从性能上来说，IIR滤波器传输函数的极点可位于单位圆内的任何地方，因此可用较低的阶数获得高的选择性，所用的存贮单元少，所以经济而效率高。但是这个高效率是以相位的非线性为代价的。选择性越好，则相位非线性越严重。相反，FIR滤波器却可以得到严格的线性相位，然而由于FIR滤波器传输函数的极点固定在原点，所以只能用较高的阶数达到高的选择性；对于同样的滤波器设计指标，FIR滤波器所要求的阶数可以比IIR滤波器高5~10倍，结果，成本较高，信号延时也较大；如果按相同的选择性和相同的线性要求来说，则IIR滤波器就必须加全通网络进行相位较正，同样要大增加滤波器的节数和复杂性。整体来看，IIR滤波器达到同样效果阶数少，延迟小，但是有稳定性问题，非线性相位；FIR滤波器没有稳定性问题，线性相位，但阶数多，延迟大。3.数据、结果和分析a.数据和结果b.结果分析:本设计完成了对语音信号的读取与打开，与课题的要求相符；本设计也较好的完成了对语音信号的频谱分析，通过fft变换，得出了语音信号的频谱图；在滤波这一块，课题主要是从巴特沃斯滤波器入手来设计滤波器，也从一方面基本实现了滤波。四、参考文献[1]高西全、丁玉美编著.数字信号处理.西安:西安电子科技大学出版社，2008.[2]丁玉美、高西全编著.数字信号处理学习指导.西安：西安电子科技大学出版社，2001.[3]郑君里等编.信号与系统.北京:高等教育出版社，2000.[4]刘树棠译.数字信号处理——使用MATLAB.西安:西安交通大学出版社，2002.[5]导向科技编著.MATLAB程序设计与实例应用.北京：中国铁道出版社，2001.[6]罗军辉等编著.MATLAB7.0在数字信号处理中的应用.北京：机械工业出版社，2005.[7]陈怀琛等编著.MATLAB及在电子信息课中的应用.北京：电子工业出版社，2002.[8]胡广书编著.数字信号处理――理论、算法与实现.北京：清华大学出版社，2002.[9]梁虹等编.信号与线性系统分析――基于MATLAB的方法与实现.北京：高等教育出版社，2006.[10]刘卫国主编.MATLAB程序设计与应用（第二版）.北京：高等教育出版社，2006.基于MATLAB的语音信号的分析和处理研究