信号处理实验七音频频谱分析报告仪设计与实现

《信号处理实验七音频频谱分析报告仪设计与实现》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实用标准哈尔滨工程大学实验报告实验名称：离散时间滤波器设计班级：电子信息工程4班学号：姓名：实验时间：2016年10月31日18：30成绩：________________________________指导教师：栾晓明实验室名称：数字信号处理实验室哈尔滨工程大学实验室与资产管理处制文案大全实用标准实验七音频频谱分析仪设计与实现一、实验原理MATLAB是一个数据分析和处理功能十分强大的工程实用软件，其数据采集工具箱为实现数据的输入和输出提供了十分方便的函数命令。本实验要求基于声卡和MTLAB实现音频信号频谱分析仪的设计原理与实现，功能包括：(1)音频信号输入，从

2、声卡输入、从WAV文件输入、从标准信号发生器输入；(2)信号波形分析，包括幅值、频率、周期、相位的估计、以及统计量峰值、均值、均方值和方差的计算。(3)信号频谱分析，频率、周期的统计，同行显示幅值谱、相位谱、实频谱、虚频谱和功率谱的曲线。1、频率(周期)检测对周期信号来说，可以用时域波形分析来确定信号的周期，也就是计算相邻的两个信号波峰的时间差、或过零点的时间差。这里采用过零点(ti)的时间差T(周期)。频率即为f=1/T，由于能够求得多个T值(ti有多个)，故采用它们的平均值作为周期的估计值。2、幅值检测在一个周期内，求出信号最大值ymax与最小值ymin的差的一

3、半，即A=(ymax-ymin)/2，同样，也会求出多个A值，但第1个A值对应的ymax和ymin不是在一个周期内搜索得到的，故以除第1个以外的A值的平均作为幅值的估计值。3、相位检测采用过零法，即通过判断与同频零相位信号过零点时刻，计算其时间差，然后换成相应的相位差。φ=2π(1-ti/T)，{x}表示x的小数部分，同样，以φ的平均值作为相位的估计值。频率、幅值和相位估计的流程如图1所示。4、数字信号统计量估计(1)峰值P的估计在样本数据x中找出最大值与最小值，其差值为双峰值，双峰值的一半即为峰值。P=0.5[max(yi)-min(yi)](2)均值估计式中，N

4、为样本容量，下同。(3)均方值估计(4)方差估计文案大全实用标准开始n=1,ymax=ymin=0i=3Yyi=1<0且yi=2<0且yi>=0且yi+1>0NNymaxN-1?i++NYT为（ti-ti-1)的平均值，i=1.2,…f=1/Tn++ymin=yiA为Ai的平均值，i=2,3,…ymin>yiφ为2π（1-{ti/Ti})的平均值结束图1　频率,幅值和相位估计的流程图其中tin表示第n个过零点

5、，yi为第i个采样点的值，Fs为采样频率。文案大全实用标准5、频谱分析原理时域分析只能反映信号的幅值随时间的变化情况，除单频率分量的简单波形外，很难明确提示信号的频率组成和各频率分量大小，而频谱分析能很好的解决此问题。由于从频域能获得的主要是频率信息，所以本节主要介绍频率(周期)的估计与频谱图的生成。(1)DFT与FFT对于给定的时域信号y，可以通过Fourier变换得到频域信息Y，Y可按下式计算式中，N为样本容量，Δt=1/Fs为采样间隔。采样信号的频谱是一个连续的频谱，不可能计算出所有的点的值，故采用离散Fourier变换(DFT)，即式中，Δf=Fs/N。但上

6、式的计算效率很低，因为有大量的指数(等价于三角函数)运算，故实际中多采用快速Fourier变换(FFT)。其原理即是将重复的三角函数算计的中间结果保存起来，以减少重复三角函数计算带来的时间浪费。由于三角函数计算的重复量相当大，故FFT能极大地提高运算效率。(2)频率、周期的估计对于Y(kΔf)，如果当kΔf=f时，Y(kΔf)取最大值，则f为频率的估计值，由于采样间隔的误差，f也存在误差，其误差最大为Δf/2，周期T=1/f。从原理上可以看出，如果在标准信号中混有噪声，用上述方法仍能够精确地估计出原标准信号的频率和周期。(3)频谱图为了直观地表示信号的频率特性，工程

7、上常常将Fourier变换的结果用图形的方式表示，即频谱图。以频率f为横坐标，

8、Y(f)

9、为纵坐标，可以得到幅值谱；以频率f为横坐标，arg Y(f)为纵坐标，可以得到相位谱；以频率f为横坐标，Re Y(f)为纵坐标，可以得到实频谱；以频率f为横坐标，Im Y(f)为纵坐标，可以得到虚频谱。根据采样定理，只有频率不超过Fs/2的信号才能被正确采集，即Fourier变换的结果中频率大于Fs/2的部分是不正确的部分，故不在频谱图中显示。即横坐标f∈[0,Fs/2](4)频谱图文案大全实用标准模块化就是把程序划分成独立命名且可独立访问的模块，每个模块完成一个子功能，把