基于MATLAB的声信号采集与频谱谱分析设计论文答辩

1、基于MATLAB的实时声信号与谱分析设计 本文介绍了基于MATLAB的声信号采集与谱分析的设计过程，在阐述频率分析及功率谱估计原理的基础上，分析了经典功率谱估计和现代功率谱估计的两大类算法，主要采用了经典功率谱估计的周期法，间接法和改进的周期法及现代功率谱估计的AR模型参数法对声信号进行谱分析，并通过大量试验对各种算法进行验证对比。该设计利用串口将采集到的声信号上传，并在MATLAB环境下进行谱分析。声传感器A/DSPIDSP串口通信谱分析谱分析 本设计的主要任务是通过串口把实时采集的数据传至计算机后，在MATLAB环境下将数据读回并实现信号的频谱分析。本设计主要研究了各种频谱分析算法实现后的

2、效果图，直观比较不同的算法的不同效果。其中包括经典功率谱分析和现代功率谱分析。经典谱估计主要方法有直接法，间接法和改进的直接法。现代谱估计主要是针对经典谱估计的分辨率低和方差性能不好等问题提出的，主要的参数模式是自回归（AR）模型、移动平均（MA）模型和自回归/移动平均（ARMA）模型，其中AR模型由线性方程描述，而MA和ARMA模型则由非线性方程描述。由于MA和ARMA模型均可用高阶的AR模型来近似，本文使用的是AR参数模型。 设计的主要任务串口传输模块 （1）本设计通过TMS320F2812内部的A/D转换器实现信号的采集，理解串口通信的原理并编写串口读数的程序，将采样结果通过串口上传到计

3、算机。 （2）MS320F2812内部具有两个相同的SCI模块，SCIA和SCIB，每一个SCI模块都各有一个接收器和发送器。SCI的接收器和发送器各具有一个16级深度的FIFO（First in fist out先入先出）队列，它们还都有自己独立的使能位和中断位。 （3）SCI发送数据的速度是由波特率来决定。所谓的波特率就是指每秒所能发送的位数。 （4）在进行通信的时候，双方都必须以相同的数据格式和波特率进行通信，否则通信会失败。 本次设计采用的波特率是115200bps，无校验，8 个数据位，1个停止位。下图是显示串口调试的界面。 在本次设计中要求对从串口接收的数据进行谱分析，由于MATL

4、AB软件相对其他软件具有许多优点，例如：语言简洁紧凑，使用方便灵活，库函数极其丰富，拥有功能强劲的工具箱等，所以本设计选择MATLAB软件进行谱分析。谱分析在谱分析在MATLAB中的实现中的实现 经典谱估计具有物理概念明确、算法简单的特点，是目前经常使用的谱估计方法。在经典谱估计中，Bartlett算法和Welch算法都是周期图法的改进算法。 周期图法是把随机序列y（n）的N个观测数据视为一个能量有限的序列，直接计算y（n）的离散傅立叶变换得Y（k），然后再取其幅值的平方，并除以N，作为序列y（n）真实功率谱的估计。 间接法先由序列x（n）估计出自相关函数R（n），然后R （n）进行傅立叶变换

5、，便得到x（n）的功率谱估计。 Bartlett平均周期图法将具有N个观测点的可用样本分成L=N/M个子样本，每个子样本有M个观测点，然后在每个w值上对所有子样本的周期图进行平均，一次来减小周期图中较大的波动 Welch平滑平均周期图法方法综合了Bartlett改进周期图法的优点，是通过先对数据分段加窗处理，然后再求平均的方法进行的。 020040060080010001200-8-6-4-20246050100150200250300350400450500-50-40-30-20-1001005010015020025030035040045050051015202530350501001

6、50200250300350400450500-30-20-10010203040 几种经典功率谱估计方法的实验仿真原始信号 直接法间接法BartlettBartlett法法050100150200250300350400450500-40-200矩 形 窗050100150200250300350400450500-40-200海 明 窗050100150200250300350400450500-40-200blackman 窗Welch法 现代功率谱估计即参数谱估计方法是通过观测数据估计参数模型再按照求参数模型输出功率的方法估计信号功率谱。主要是针对经典谱估计的分辨率低和方差性能不好等问

7、题提出的。主要方法有最大的谱分析法(AR模型法)、Pisarenko谐波分解法、Prony提取极点法、Prony谱线分解法以及Capon最大似然法等。其中AR模型由线性方程描述，而MA和ARMA模型则由非线性方程描述。由于MA和ARMA模型均可用高阶的AR模型来近似，本文使用的是AR参数模型。 AR模型的物理意义是认为序列y(n)是白噪声作用于线性时不变系统时的系统响应。 本设计中采用AR模型的协方差算法，自相关算法，Burg算法。现代估率谱估计 几种现代功率谱估计的实验仿真 仿真结果直观地说明了自相关算法、Burg算法和改进协方差算法各自的优缺点。对相隔仅有5HZ的2个混合信号进行谱估计时，

8、自相关法不容易看出其频率成份，而Burg算法和改进协方差算法提高了参数估计的精度和频率分辨率。而Burg方法，用其进行功率谱估计时令前后向预测误差功率之和最小，即对，前后都不加窗，Burg算法是建立在数据基础之上的，避免了先计算自相关函数从而提高计算速度，计算不太复杂，分辨率优于自相关法。 MATLAB提供了一套可视化的创建图形窗口的工具，使用用户界面开发环境可方便的创建GUI应用程序，它可以根据用户设计的GUI布局，自动生成M文件的框架，用户使用这一框架编制自己的应用程序。图形用户界面就是通过窗口、选单、按钮、文字说明等对象构成一个美观的界面，用户利用鼠标或键盘可以方便地实现操作。 基于MA

9、TLAB用户界面的谱分析 本设计中用户界面，大体分为信号输入模块，原始信号显示模块，频谱分析模块，经典谱分析模块和现代谱分析模块四部分。经典谱分析模块包含直接法、间接法、Welch法和Bartlett法四个小模块。现代谱分析模块则包括Burg法和自相关法两个小模块。 根据不同的信号来源，界面主要采用两种不同输入方式。一种信号来源是由串口实时采集到的信号，把该信号作为一种输入方式。然后点击打开串口，其次点击频谱分析按钮，坐标轴会显示信号的幅度谱。 用户界面的实现过程 最后就是功率谱分析部分， 串口开始读取数据同时把数据保存在D盘下的一个文件中。接着关闭串口，点击打开按钮，即可读回由串口实时采集到

10、的声信号数据。打开文件的同时首先在坐标轴中会显示原始信号波形。其次点击频谱分析按钮，坐标轴会显示信号的幅度谱。最后就是功率谱分析部分。 另一种信号来源是打开mat文件中已存的两个不同频率叠加的声音信号，把该信号作为一种输入方式。然后点打开按钮，在打开文件名的同时首先会在坐标轴中显示原始信号波形。其他的操作如上第一种方式。 输入信号为串口接收的数据界面 输入信号为mat文件已存的数据 通过实验仿真可以直观地看出以下特性： （1）经典功率谱估计中的周期图法所得到的结果特点是离散性大，曲线粗糙，方差较大，但是分辨率较高。 （2）Bartlett平均周期图法和Welch平滑平均周期图法的收敛性较好，曲线平滑，估计的结果方差较小，但是功率谱主瓣较宽，分辨率低。这是由于对随机序列的分段处理引起了长度有限所带来的Gibbs现象而造成的。 （3）与Bartlett法相比，Welch法的估计曲线比较粗糙，但是分辨率较好，原因是Welch法中对数据进行截断时加的是Hanning窗，而在Bartlett法中使用的是矩形窗，相对于矩形窗，窗的主瓣包含更多的能量，因而使功率谱的主瓣较窄，分辨率较高。 总结 本设计主要是对于信号的频率分析算法进行研究，所谓的频率分析即