基于MATLAB的语音信号录制采集和窗函数法滤波器设计(优秀篇)

PAGEword文档可自由复制编辑PAGE1课程设计说明书题目：基于MATLAB的语音信号录制采集和窗函数法滤波器设计姓名：院（系）：电气信息工程学院专业班级：电子信息工程11-1班学号：指导教师：成绩：摘要随着计算机和信息科学的飞速发展，数字信号处理逐渐发展成为一门独立的学科，成为信息科学的重要组成部分，在语音处理、雷达、图像处理、通信、生物医学工程等众多领域中得到广泛应用。Matlab语言是一种广泛应用于工程计算及数值分析领域的新型高级语言，Matlab功能强大、简单易学、编程效率高，深受广大科技工作者的喜爱。特别是Matlab还具有信号分析工具箱，不需具备很强的编程能力，就可以很方便地进行语音信号分析、处理和设计。语音信号的抽样、频谱分析，滤波讨论分析，通过数字信号处理课程的理论知识的综合运用。从实践上初步实现对数字信号的处理。关键词：语音处理抽样频谱分析滤波器word文档可自由复制编辑目录TOC\o"1-3"\h\u17091绪论 122875第一章课程设计方案提出 219610第一节设计要求 27516第二节设计方案 223875第二章语音信号的采集及加噪 312789第一节语音信号的录制 313791第二节语音信号的时域分析 49879第三节语音信号加噪和频谱分析 66074第三章滤波器的选择 74187第一节滤波器的基本概念 716523一滤波器的功能 710035二四种基本的滤波器 732068第二节数字滤波器的基本知识 819946一FIR数字滤波器 820633二IIR数字滤波器 925362第三节数字滤波器的选取 912999第四章心得体会 1027197参考文献 11绪论MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室），是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。在开发环境中，使用户更方便地控制多个文件和图形窗口；在编程方面支持了函数嵌套，有条件中断等；在图形化方面，有了更强大的图形标注和处理功能，包括对性对起连接注释等；在输入方面，可以直接向Excel和HDF5进行连接。第一章课程设计方案提出第一节设计要求要求录制一段自己的语音信号后，在MATLAB软件中采集语音信号、回放语音信号并画出语音信号的时域波形和频谱图；对所采集的信号加入干扰噪声，对加入噪声进行播放，并进行时域和频谱分析；对比及早前后的时域图和频谱图，分析讨论采用什么的滤波器进行滤除噪声。第二节设计方案Windows自带的录音机Windows自带的录音机声音麦克风声卡滤波采样A/D转换Wav首先，录制一段语音信号（如图1.1）；然后，利用MATLAB中的wavread命令来读入（采集）语音信号，将它赋值给某一向量。再将该向量看作一个普通的信号，对其进行FFT变换实现频谱分析；随后，利用MATLAB中的随机函数(rand或randn)产生噪声加入到语音信号中，模仿语音信号被污染，并对其频谱分析，在讨论用什么滤波器去除噪声。

第二章语音信号的采集及加噪第一节语音信号的录制利用PC机上的声卡和WINDOWS操作系统可以进行数字信号的采集。将话筒输入计算机的语音输入插口上,启动录音机。按下录音按钮,接着对话筒说话“语音信号处理”,说完后停止录音,屏幕左侧将显示所录声音的长度，点击放音按钮,可以实现所录音的重现。以文件名“shijiebei”保存入D:\MATLAB\work中。文件存储器的后缀默认为.wav,这是WINDOWS操作系统规定的声音文件存的标准。具体操作如下：a.打开开始/程序/附件/娱乐/录音机（如图2.1）;图2.1b.用windowsmediaplayer播放一首音乐并用MIC对着耳机录音或自已说话录音（按键），到10秒时停止（按键）（如图2.2）；图2.2c.将录制文件加存为D:/MATLAB7/work中，文件名为shijiebei.wav（如图2.2、图2.3）；图2.3第二节语音信号的时域分析在MATLAB软件平台下，利用MATLAB中的“wavread”命令来读入（采集）语音信号，将它赋值给某一向量。再对其进行采样，记住采样频率和采样点数。通过使用wavread函数，可以深入理解采样频率、采样位数等概念，下面介绍Wavread函数几种调用格式。（1）y=wavread（file）；读取file所规定的wav文件，返回采样值放在向量y中。（2）[y,fs,nbits]=wavread(file)；采样值放在向量y中，fs表示采样频率（hz），nbits表示采样位数。（3）y=wavread（file，N）；读取钱N点的采样值放在向量y中。（4）y=wavread（file，[N1,N2]）；读取从N1到N2点的采样值放在向量y中。接下来，对语音信号shijiebei.wav进行采样。其采样程序如下：[y,fs,nbits]=wavered(‘E:\sheji\shijiebei’);即把语音信号加载入Matlab仿真软件平台中。然后，画出语音信号的时域波形，再对语音信号进行频谱分析。MATLAB提供了快速傅里叶变换算法FFT计算DFT的函数fft,其调用格式如下：Xk=fft(xn,N)参数xn为被变换的时域序列向量，N是DFT变换区间长度，当N大于xn的长度时，fft函数自动在xn后面补零。，当N小于xn的长度时，fft函数计算xn的前N个元素，忽略其后面的元素。在本次设计中，我们利用fft对语音信号进行快速傅里叶变换，就可以得到信号的频谱特性。其程序如下：[y,fs,nbits]=wavread('shijiebei');sound(y,fs,nbits);%回放语音信号N=length(y);%求出语音信号的长度Y=fft(y,N);%傅里叶变换subplot(2,1,1);plot(y);title('原始信号波形');subplot(2,1,2);plot(abs(Y));title('原始信号频谱')程序运行结果如下图：图2.4原始信号波形及频谱第三节语音信号加噪和频谱分析MATLAB中产生高斯白噪声非常方便，有两个产生高斯白噪声的两个函数。我们可以直接应用两个函数：一个是WGN，另一个是AWGN。WGN用于产生高斯白噪声，AWGN则用于在某一信号中加入高斯白噪声。也可直接用randn函数产生高斯分布序列。在本次课程设计中，我们是利用MATLAB中的随机函数(rand或randn)产生噪声加入到语音信号中，模仿语音信号被污染，并对其频谱分析。Randn函数有两种基本调用格式：Randn(n)和Randn（m,n）,前者产生n×n服从标准高斯分布的随机数矩阵，后者产生m×n的随机数矩阵。在这里，我们选用Randn（m,n）函数。语音信号添加噪声及其频谱分析的主要程序如下：[y,fs,nbits]=wavread('shijiebei’);N=length(y);%求出语音信号的长度Noise=0.01*randn(n,2); %随机函数产生噪声Si=y+Noise;%语音信号加入噪声sound(Si);subplot(2,1,1);plot(Si);title('加噪语音信号的时域波形');S=fft(Si);%傅里叶变换subplot(2,1,2);plot(abs(S));title('加噪语音信号的频域波形');程序运行结果如下图：图2.5加噪后的波形及频谱分析第三章滤波器的选择第一节滤波器的基本概念1、滤波器的功能滤波器的功能是对输入信号进行滤波以增强所需信号部分，抑制不要的部分。2、四种基本的滤波器四种基本滤波器为低通（LP）、高通(HP)、带通(BP)和带阻滤波器(BRP)，这四种滤波器特性如下图3-1所示：图3.1第二节数字滤波器的基本知识数字滤波器分为FIR数字滤波器和IIR数字滤波器两种，即有限冲激响应滤波器（FIR，FiniteImpulseResponse）滤波器和无线冲激响应（IIR，InfiniteImpulseResponse）滤波器。1、FIR数字滤波器（1）FIR的特点：（3-1）不存在极点（z=0除外），系统函数|z>0|处收敛。系统单位冲激响应在有限个n值处不为零。结构上主要是非递归结构，没有输出到输入的反馈。故只能用较高的阶数达到高的选择性。FIR数字滤波器的幅频特性精度较之于IIR数字滤波器低，但是线性相位即不同频率分量的信号经过FIR数字滤波器后他们的时间差不变。FIR数字滤波器具有系统稳定，易实现相位控制，允许设计多带通滤波器等优点。（2）FIR结构FIR数字滤波器的结构有四种：横截性；级联型结构；频率采样型结构；线性相位FIR滤波器的结构。二IIR数字滤波器1、IIR数字滤波器的特点电位冲激响应h(n)是无限长的；系统函数H(z)在有限z平面上（0<|z|<∞）有极点存在；结构上存在着输出到输入的反馈，也就是结构上是递归型的。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成。2、有限阶IIR的表达式：（3-2）（3-3）3、IIR滤波器的四种结构IIR数字滤波器的结构也有四种：直接I型；直接II型—典范型；级联型；并联型。第三节数字滤波器的选取由图2.4和图2.5可以看出，语音信号的能量集中在低频部分，所以根据滤波器的特性和功能，所选择的滤波器应该具有滤除无用的高频成分的作用，所以可以选用低通滤波器。又由于IIR数字滤波器的输出不但取决于过去和现在的输入，还取决于过去的输出，所以，选择IIR数字滤波器比选择FIR数字滤波器要好。通过徐靖涛.基于MATLAB的语音信号分析与处理[J].重庆科技学院学报2008,1中的2.4的阅读，可知，本语音信号的滤波选取椭圆型数字低通滤波器效果最好。第四章心得体会这次课程设计历时一个星期，可以说是苦多于甜，但是可以学的到很多很多的东西，同时不仅可以巩固以前所学过的知识，而且学到了很多在书