实验八音频频谱分析仪设计与实现.doc

实验八：音频频谱分析仪设计与实现 一、实验原理： MATLAB是一个数据分析和处理功能十分强大的工程实用软件，其数据采集工具箱为实现数据的输入和输出提供了十分方便的函数和命令。本实验要求基于声卡与MATLAB实现音频信号频谱分析仪的设计原理与实现，功能包括： （1）音频信号输入，从声卡输入，从WAV文件输入，从标准信号发生器输入； （2）信号波形分析，包括幅值，频率，周期，相位的估计，以及统计量峰值，均值，均方值和方差的计算； （3）信号频谱分析，频率，周期的估计，图形显示幅值谱，相位谱，实频谱，虚频谱和功率谱的曲线。 1.频率（周期）检测 对周期信号来说，可以用时域波形分析来确定信号的周期，也就是计算相邻的两个信号波峰的时间差，或过零点的时间差。这里采用过零点（t）的时间差T（周期）。频率即为频率即为f=1/T，由于能够求得多个T值，故采用他们的平均值作为周期的估计值。 2.幅值检测 在一个周期内，求出信号最大值与最小值的差得一半记为A，同样得到多个A值，但第一个A值对应的和不是在一个周期内搜索得到的，故以除第一个以外的A值的平均作为幅值的估计值。 3.相位检测 采用过零法，即通过判断与同频零相位信号过零点时刻，计算其时间差，然后换成相应的相位差。，同样以的平均值作为相位的估计值。 4.数字信号统计量估计 （1）峰值P的估计 在样本数据x中找出最大值与最小值，其差为双峰值，双峰值的一半即为峰值。 P=0.5[max(y)-min(y)] （2）均值估计 N为样本容量。 （3）均方值估计 （4）方差估计 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 = 时，Y(kΔf)取最大值，则为频率的估计值，由于采样间隔的误差，也存在误差，其误差最大为Δf / 2，周期T=1/f。 从原理上可以看出，如果在标准信号中混有噪声，用上述方法仍能够精确地估计出原标准信号的频率和周期。 （3）频谱图 为了直观地表示信号的频率特性，工程上常常将Fourier变换的结果用图形的方式表示，即频谱图。 以频率f为横坐标，|Y(f)|为纵坐标，可以得到幅值谱； 以频率f为横坐标，arg Y(f)为纵坐标，可以得到相位谱； 以频率f为横坐标，Re Y(f)为纵坐标，可以得到实频谱； 以频率f为横坐标，Im Y(f)为纵坐标，可以得到虚频谱。 根据采样定理，只有频率不超过Fs/2的信号才能被正确采集，即Fourier变换的结果中频率大于Fs/2的部分是不正确的部分，故不在频谱图中显示。即横坐标f ∈[0, Fs/2] （4）模块划分 模块化就是把程序划分成独立命名且可独立访问的模块，每个模块完成一个子功能，把这些模块集成起来构成一个整体，可以完成制定的功能，满足用户需求。 在模块分化时应遵循如下原则：改进软件结构提高模块独立性；模块规模应该适中；深度、宽度、扇出和扇入都应适当；模块的作用域应该在控制域之内；力争降低模块接口的复杂程度；设计单入口单出口的模块；模块功能应该可以预测。 二、实验内容： 参考给定的代码，构建交互界面，设计一个音频频谱分析仪，实现相应的功能。 程序文本： function varargout = GGfig(varargin) % GGFIG M-file for GGfig.fig % GGFIG, by itself, creates a new GGFIG or raises the existing % singleton*. % % H = GGFIG returns the handle to a new GGFIG or the handle to % the existing singleton*. % % G