关于1/3倍频程的计算方法

        对于人耳听音频率的分辨能力可以用频带来形容，而1/3倍频程被认为是比较符合人耳特性的频带划分方法。因此在噪声测量中，采用1/3倍频程频谱分析方法能够更加详细地反映出噪声源的频谱特性，较为全面地了解声源产生机理和提出最佳的降噪对策。比如针对一段信号，要将其用1/3倍频程频谱展示出来，首先要知道1/3倍频程的划分方法；然后，使用傅里叶变换（FFT）将时间域的声信号变换到频率域，计算所定义的每个1/3倍频带的声压等效连续声压级；最后，得到1/3倍频程声压级。

        任意一段声信号的1/3倍频程谱图的大致计算步骤如下：

        Step1: 读取声信号，进行FFT变换；

        Step2: 按照上述表格划分频带：若信号频率处于20~16000Hz范围，则可以得到30个频带；若信号频率范围处于1000~8000Hz范围，则可以得到10个频带；

        Step3: 对于第i个频带，首先根据该频带截止频率fl、fu及采样频率fs计算其在频率坐标轴上对应的点数，由此从FFT得到的阵列中截取出[fl,fu]区间内的频谱数据（其他区间幅值均为0）并还原为双边谱；然后进行IFFT变换，得到该频带对应的时域信号；最后，计算该频带内的声压有效值yc(i)；

        Step4: 根据声压有效值yc(i)，计算第i个频带的未计权声压级LPl(i)；

        Step5: 根据A计权修正值cf(i)，计算第i个频带的A计权声压级Lap (i)；

        Step6: 当所有频带内的声压级都依次计算完成后，即可得到1/3倍频程谱图。 

图1  计算流程图

 

图2  1/3倍频程的中心频率和带宽

 

图3  时域波形图

 

图4  频谱图

 

图5  未计权1/3倍频程声压级图谱

 

图6  A计权1/3倍频程声压级图谱