MATLAB课程设计

郑州轻工业学院MATLAB课程设计题目：基于MATLAB的语音信号滤波处理和程序设计姓名：XXX院（系）：电气信息工程专业班级：电信11-01学号：11111111111111指导教师：XXX成绩：时间：2014年6月09日至2014年6月13日目录摘要............................................................................................................................................11、绪论......................................................................................................................................21.1语音信号分析技术.......................................................................................................21.2数字滤波技术...............................................................................................................21.3主要内容和基本要求...................................................................................................32、实验内容..............................................................................................................................42.1语音信号的采集...........................................................................................................42.2语音信号的频谱分析...................................................................................................42.3构建干扰信号...............................................................................................................52.4信号合成.......................................................................................................................72.5用数字滤波器对加噪的语音信号进行滤波...............................................................8心得体会..................................................................................................................................11参考文献..................................................................................................................................121基于MATLAB的语音信号滤波处理和程序设计摘要语音信号处理属于信息科学的一个重要分支，大规模集成技术的高度发展和计算机技术的飞速前进，推动了这一技术的发展。Matlab是一种数据分析和处理功能十分强大的计算机应用软件,信号处理是Matlab重要应用的领域之一。本设计是基于Matlab对语音信号进行处理，先完成语音信号的采集，利用MATLAB对所录制的信号进行时域和频域分析；构建一个干扰信号，把语音信号和干扰信号叠加，利用MATLAB对干扰后的语音信号进行时域和频域分析，对干扰后的信号进行滤波，画出滤波后的信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化。关键词MATLAB语音信号干扰信号滤波信号叠加21、绪论1.1语音信号分析技术语音信号分析是语音信号处理的前提和基础，只有分析出可表示语音信号本质特征的参数，才有可能利用这些参数进行高效的语音通信、语音合成和语音识别等处理。而且，语音合成的音质好坏，语音识别率的高低，也都取决于对语音信号分桥的准确性和精确性。因此语音信号分析在语音信号处理应用中具有举足轻重的地位。根据所分析出的参数的性质的不同，可将语音信号分析分为时域分析、频域分析、倒领域分析等；时域分析方法具有简单、计算量小、物理意义明确等优点，但由于语音信号最重要的感知特性反映在功率谱中，而相位变化只起着很小的作用，所以相对于时域分析来说频域分析更为重要。通过对大量语音信号的观察和分析发现，语音信号主要有下面两个特点：①在频域内，语音信号的频谱分量主要集中在300～3400Hz的范围内。利用这个特点，可以用一个防混迭的带通滤波器将此范围内的语音信号频率分量取出，然后按8kHz的采样率对语音信号进行采样，就可以得到离散的语音信号。②在时域内，语音信号具有“短时性”的特点，即在总体上，语音信号的特征是随着时间而变化的，但在一段较短的时间间隔内，语音信号保持平稳。在浊音段表现出周期信号的特征，在清音段表现出随机噪声的特征。1.2数字滤波技术数字滤波是数字信号处理的重要基础，数字信号处理主要是研究用数字或符号的序列来表示信号波形，并用数字的方式去处理这些序列，把它们改变成在某种意义上更希望的形式，以便估计信号的特征参量，或削弱信号中的多余分量和增强信号中的有用分量。数字滤波器在对信号的过滤、检测与参数估计等处理过程中，是使用最为广泛的一种线性系统。数字滤波器（DigitalFilter，简称为DF）是指完成信号滤波处理功能的、用有限精度算法实现的离散时间线性非时变系统。数字滤波器的数学运算通常有两种实现方式。一种是频域法，即利用FFT快速运算办法对输入信号进行离散傅立叶变换，分析其频谱，3然后根据希望的频率特性进行滤波，再利用傅立叶反变换恢复出时域信号。这种方法具有较好的频域选择特性和灵活性，并且由于信号频率与所希望的频谱特性是简单的相乘关系，所以它比计算等价的时域卷积要快得多。另一种方法是时域法，这种方法是通过对离散抽样数据作差分数学运算来达到滤波目的的。数字滤波器的输入是一组（由模拟信号取样和量化的）数字量，其输出是经过数字变换的另一组数字量。数字滤波器具有稳定性高、精度高、灵活性大等突出优点。随着数字技术的发展，用数字技术实现滤波器的功能愈来愈受到人们的重视，并得到了广泛的应用。数字信号处理学科的一项重大进展是关于数字滤波器设计方法的研究。60年代中期，开始形成关于数字滤波器的一套完整的正规理论。这一时期，提出了各种各样的数字滤波器结构，有的以运算误差最小为特点，有的则以运算速度高见长；出现了数字滤波器的各种逼近方法和实现方法，对递归和非递归两类滤波器作了全面的比较；统一了数字滤波器的基本概念和理论，对有限冲激响应（IIR）和无限冲激响应（FIR）的认识有了完整理论。70年代后，科学技术蓬勃发展，数字信号处理开始与大规模和超大规模集成电路技术、微处理技术等新工艺新技术结合起来，并引进计算机辅助设计方法，大大丰富了数字滤波器的分析与设计，各种新的数字信号处理系统，也都能用专用数字硬件实时加以实现。相信在未来，随着电子仪器与电子技术应用系统朝着数字化、小型化、自动化以及多功能化等方向发展，包括数字滤波器在内的数字信号处理技术会有以惊人的速度进行飞跃式发展。1.3主要内容和基本要求主要内容：要求用Windows下的录音机，录制一段自己的话语“大家好”。利用MATLAB对所录制的信号进行时域和频域分析；构建一个干扰信号，把语音信号和干扰信号叠加，利用MATLAB对干扰后的语音信号进行时域和频域分析，对干扰后的信号进行滤波，画出滤波后的信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化。基本要求：1、掌握在MATLAB环境下语音信号的采集方法。2、构建的干扰信号为f(t)=10sin(4000𝜋t)。43、掌握数字信号处理的基本概念、基本理论、基本方法。4、利用MATLAB设计滤波器对干扰的语音信号进行滤波。5、掌握MATLAB的程序设计方法。2、实验内容2.1语音信号的采集利用Windows下的录音机,录制一句自己的话，保存为“大家好.wav”文件。然后在Matlab软件平台下,利用函数wavread对语音信号进行采样,记住采样频率和采样点数。使用[y,fs,bits]=wavread();读取音频文件。录制的音频文件为wam格式，需转换为wav格式的文件。2.2语音信号的频谱分析首先画出语音信号的时域波形，然后对语音信号进行频谱分析。在Matlab中,可以利用函数fft对信号进行快速傅里叶变换,得到信号的频谱特性。程序：[y,fs,bits]=wavread('C:\Users\Administrator\Desktop\matlab资料\大家好.wav');n=length(y)y_p=fft(y,n);%对n点进行傅里叶变换到频域f=fs*(0:n/2-1)/n;%对应点的频率figure(1)subplot(2,1,1);plot(y);%语音信号的时域波形图title('原始语音信号时域波形');xlabel('时间')ylabel('幅度')5subplot(2,1,2);plot(f,abs(y_p(1:n/2)));%语音信号的频谱图axis([0200000600]);title('原始语音信号频谱图');xlabel('频率Hz');ylabel('幅度');图1原始语音信号波形图2.3构建干扰信号构建的干扰信号为f(t)=10sin(4000𝜋t)，此信号为正弦信号。程序：[y,fs,bits]=wavread('C:\Users\Administrator\Desktop\matlab资料\大家好.wav');n=1000t=(0:n-1)/100000;010002000300040005000600070008000900010000050010001500原始语音信号频谱图频率Hz幅度0123456789x104-1-0.500.51原始语音信号时域波形时间幅度6f=fs*(0:n/2-1)/n;y1=10*sin(2*pi*2000*t);y2=fft(y1,10*n);subplot(2,1,1);plot(y1)%干扰信号的时域图title('干扰信号时域图');xlabel('时间')ylabel('幅度')subplot(2,1,2);plot(f,abs(y2(1:n/2)));%干扰信号的频谱图title('干扰信号频谱图');xlabel('频率');ylabel('幅度');图2干扰信号波形图01002003004005006007008009001000-10-50510干扰信号时域图时间幅度00.511.522.5x1040200040006000干扰信号频谱图频率幅度72.4信号合成用MATLAB将干扰信号和原始信号进行合成。程序：L=lengt