音频频谱分析仪设计与实现剖析.pdfVIP

实验八 音频频谱分析仪设计与实现 一、实验原理 MATLAB 是一个数据分析和处理功能十分强大的工程实用软件，其数据采集工具箱为 实现数据的输入和输出提供了十分方便的函数和指令。 本实验基于声卡与 MATLAB 实现音频信号频谱分析仪。 1、频率 (周期 )检测 对周期信号来说，可以用时域波形分析来确定信号的周期，也就是计算相邻两个信号 波峰的时间差或过零点的时间差。采用过零点（ ti ）的时间差 T （周期）。频率即为 f=1/T ， 由于能够求得多个 T 值，故采用他们的平均值作为周期的估计值。 2 、峰值检测 在一个周期内，求出信号最大值 y max 与最小值 y min 的差得一半记为 A ，同样得到多个 A 值，但第一个 A 值对应的 y max 和 y min 不是在一个周期内搜索得到的，故以除第一个以外 的 A 值的平均作为幅值的估计值。 3 、相位检测 采用过零法，即通过判断与同频零相位信号过零点时刻，计算其时间差，然后换成相 应的相位差。 2π （1- ti / T ），同样以 的平均值作为相位的估计值。 4 、数字信号统计量估计 （1）峰值 P 的估计 在样本数据 x 中找出最大值与最小值，其差为双峰值，双峰值的一半即为峰值。 （2 ）均值估计 Ｎ 1 E( y ) ｙ ，N 为样本容量。 ｉ N ｉ＝０ （3）均方值估计 Ｎ 2 1 2 E ( y ) ｙ ｉ N ｉ＝０ （4 ）方差估计 Ｎ 1 2 D (y ) (ｙ - E ( y )) ｉ N ｉ＝０ 5 、频谱分析原理 时域分析只能反映信号的幅值随时间的变化情况， 除但频率分量的简单波形外， 很难 明确提示信号的频率组成和各频率分量大小，而频谱分析能很好的解决此问题。 （1）DFT 与 FFT 对于给定的时域信号 y ，可以通过 Fourier 变换得到频域信息 Y 。Y 可按下式计算 式中， N 为样本容量， Δt = 1/Fs 为采样间隔。 采样信号的频谱是一个连续的频谱，不可能计算出所有的点的值，故采用离散 Fourier 变换 (DFT) ，即 式中， Δf = Fs/N 。但上式的计算效率很低，因为有大量的指数 (等价于三角函数 )运算，故实 际中多采用快速 Fourier 变换 (FFT) 。其原理即是将重复的三角函数算计的中间结果保存起 来， 以减少重复三角函数计算带来的时间浪费。 由于三角函数计算的重复量相当大， 故 FFT 能极大地提高运算效率。 （2 ） 频率、周期的估计 ~ ~ 对于 Y(k Δf) ，如果当 k Δf = f 时， Y(k Δf) 取最大值，则 f 为频率的估计值，由于采样 ~ 间隔的误差， f 也存在误差，其误差最大为 Δf / 2 。周期 T=1/f 。