基于LabVIEW的多采样率数字滤波器1／3倍频程计算方法-Open.PDF

No．20 第15卷第20期2015年7月科学技术与工程 V01．15 Jul．2015 Science and 1671—1815(2015)20—0095—07 TechnologyEngineering ⑥2015Sci．Tech．Engrg． 基于LabVIEW的多采样率数字 滤波器1／3倍频程计算方法 蔡兴唐 求‘ (湖南大学电气与信息工程学院，长沙410082) 摘要 利用多采样率数字滤波器法计算噪声1／3倍频程谱在硬件上设计困难。因此为降低设计的难度，提出基于Lab— VIEW实现多采样率数字滤波器法计算1／3倍频程。信号首先通过抗混叠滤波器降低信号频谱混叠，然后经过降采样减少数 据计算量，最后利用数字带通滤波器实现通带内信号的能量衰减并进行均方根计算得到1／3倍频程谱值。通过LabVIEW设 计出FIR抗混叠滤波器和8阶巴特沃斯IIR数字带通滤波器，并在此基础上结合声卡搭建噪声分析系统完成1／3倍频程计 算。结果表明，所设计的系统能快速准确地计算1／3倍频程谱值，得到理想的低频处频谱分析。 关键词1／3倍频程 多采样率数字滤波器LabVIEW噪声分析系统 中图法分类号TN713．7； 文献标志码A 噪声测量中，1／3倍频程占有重要地位。相比 问题，减少计算量，得到较为理想的低频处的频谱 倍频程而言，1／3倍频程对频带的划分更细，详细的 分析[3】。 反映噪声信号频谱特性和能量分布情况，便于观察 针对频域算法和多采样率数字滤波器算法的优 噪声信号宏观上的能量分布，消除信号频率或相位 缺点，采用LabVIEW设计多采样率数字滤波器组实 信息微小变化对观察结果的影响‘1J。 现1／3倍频程分析，降低多采样率数字滤波器算法 1／3倍频程谱值的计算多采用频域算法和多采 的设计难度，同时提高声音信号低频段1／3倍频程 样率数字滤波器法。频域算法首先对噪声采样信号 频分析的准确度。 进行快速傅里叶变换，计算出功率谱，然后计算每个 多采样率数字滤波器原理 中心频率的带宽内功率谱总和，得到1／3倍频程频 谱值…。为对噪声信号进行准确分析，通常采用很 多采样率数字滤波器法计算噪声1／3倍频程谱 大的FFT长度(大于或等于8192点)来提高分析的值时，将20Hz到20kHz的声频范围划分为30个 准确性旧J。而多采样率数字滤波器算法对信号进 频带j，其原理框图如图1所示。图1中，h(n)为 行分析前会先进行信号抽取，减少需要处理的数据 低通滤波器即抗混叠滤波器，J2表示对信号进行 和滤波运算量，然后输入带通滤波器进行计算求出 抽取，IIR为数字带通滤波器。由图1可知，中心频 率为25Hz和31．5Hz两频带数据经过四次抗混叠 1／3倍频程频谱旧’41。1／3倍频程分析中，有30个 频带谱值，因此频域算法在硬件上比数字滤波器算 滤波和信号抽取，40一80Hz四频带数据经过三次 法实现更简单，计算量相同的情况下运算速度比数 抗混叠滤波和信号抽取，100～160Hz三频带数据 字滤波器要快∞J。但是1／3倍频程分析中心频率为 经过两次抗混叠滤波和信号抽取，200Hz频带的数 Hz一20 20Hz的带宽只有4．6Hz，在低频率段，频域算法对 据经过一次抗混叠滤波和信号抽取，250 噪声信号的频率分辨率无法满足测量要求。只有提 kHz二十频带数据均不经过抗混叠滤波和信号抽 高采样点数才能做到，而通常噪声测量中每次测量