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

基于MATLAB的实时声信号与谱分析设计 本文介绍了基于MATLAB的声信号采集与谱分析的设计过程，在阐述频率分析及功率谱估计原理的基础上，分析了经典功率谱估计和现代功率谱估计的两大类算法，主要采用了经典功率谱估计的周期法，间接法和改进的周期法及现代功率谱估计的AR模型参数法对声信号进行谱分析，并通过大量试验对各种算法进行验证对比。该设计利用串口将采集到的声信号上传，并在MATLAB环境下进行谱分析。 声 传 感 器 A/D SPI DSP 串 口 通 信 谱 分 析 谱分析 本设计的主要任务是通过串口把实时采集的数据传至计算机后，在MATLAB环境下将数据读回并实现信号的频谱分析。本设计主要研究了各种频谱分析算法实现后的效果图，直观比较不同的算法的不同效果。其中包括经典功率谱分析和现代功率谱分析。经典谱估计主要方法有直接法，间接法和改进的直接法。现代谱估计主要是针对经典谱估计的分辨率低和方差性能不好等问题提出的，主要的参数模式是自回归（AR）模型、移动平均（MA）模型和自回归/移动平均（ARMA）模型，其中AR模型由线性方程描述，而MA和ARMA模型则由非线性方程描述。由于MA和ARMA模型均可用高阶的AR模型来近似，本文使用的是AR参数模型。 设计的主要任务 串口传输模块 （1）本设计通过TMS320F2812内部的A/D转换器实现信号的采集，理解串口通信的原理并编写串口读数的程序，将采样结果通过串口上传到计算机。 （2）MS320F2812内部具有两个相同的SCI模块，SCIA和SCIB，每一个SCI模块都各有一个接收器和发送器。SCI的接收器和发送器各具有一个16级深度的FIFO（First in fist out先入先出）队列，它们还都有自己独立的使能位和中断位。 （3）SCI发送数据的速度是由波特率来决定。所谓的波特率就是指每秒所能发送的位数。 （4）在进行通信的时候，双方都必须以相同的数据格式和波特率进行通信，否则通信会失败。 本次设计采用的波特率是115200bps，无校验，8 个 数据位，1个停止位。下图是显示串口调试的界面。 在本次设计中要求对从串口接收的数据进行谱分析，由于MATLAB软件相对其他软件具有许多优点，例如：语言简洁紧凑，使用方便灵活，库函数极其丰富，拥有功能强劲的工具箱等，所以本设计选择MATLAB软件进行谱分析。 谱分析在MATLAB中的实现 经典功率谱估计 经典谱估计具有物理概念明确、算法简单的特点，是目前经常使用的谱估计方法。在经典谱估计中，Bartlett算法和Welch算法都是周期图法的改进算法。 周期图法是把随机序列y（n）的N个观测数据视为一个能量有限的序列，直接计算y（n）的离散傅立叶变换得Y（k），然后再取其幅值的平方，并除以N，作为序列y（n）真实功率谱的估计。 间接法先由序列x（n）估计出自相关函数R（n），然后R （n）进行傅立叶变换，便得到x（n）的功率谱估计。 Bartlett平均周期图法将具有N个观测点的可用样本分成L=N/M个子样本，每个子样本有M个观测点，然后在每个w值上对所有子样本的周期图进行平均，一次来减小周期图中较大的波动 Welch平滑平均周期图法方法综合了Bartlett改进周期图法的优点，是通过先对数据分段加窗处理，然后再求平均的方法进行的。 几种经典功率谱估计方法的实验仿真 原始信号 直接法 间接法 Bartlett法 Welch法 现代功率谱估计即参数谱估计方法是通过观测数据估计参数模型再按照求参数模型输出功率的方法估计信号功率谱。主要是针对经典谱估计的分辨率低和方差性能不好等问题提出的。主要方法有最大的谱分析法(AR模型法)、Pisarenko谐波分解法、Prony提取极点法、Prony谱线分解法以及Capon最大似然法等。其中AR模型由线性方程描述，而MA和ARMA模型则由非线性方程描述。由于MA和ARMA模型均可用高阶的AR模型来近似，本文使用的是AR参数模型。 AR模型的物理意义是认为序列y(n)是白噪声作用于线性时不变系统时的系统响应。 本设计中采用AR模型的协方差算法，自相关算法，Burg算法。 现代估率谱估计 几种现代功率谱估计的实验仿真 仿真结果直观地说明了自相关算法、Burg算法和改进协方差算法各自的优缺点。对相隔仅有5HZ的2个混合信号进行谱估计时，自相关法不容易看出其频率成份，而Burg算法和改进协方差算法提高了参数估计的精度和频率分辨率。而Burg方法，用其进行功率谱估计时令前后向预测误差功率之和最小，即对，前后都不加窗，Burg算法是