基于音调能量和人耳频率选择性的减变速机啸叫评判方法与流程

本发明涉及减变速机啸叫异响评价技术领域，特别涉及一种基于音调能量和人耳频率选择性的减变速机啸叫评判方法。

背景技术：

减变速机是用来改变来自发动机转速及转矩的机构，能够实现不同转速的输出，具有效率高、构造简单、使用方便的优点，广泛应用于船舶、车辆等领域中。减变速机的齿轮系统一般是由常啮合齿轮对、挡位啮合齿轮对及惰轮啮合齿轮对组成，齿轮在生产过程中，由于加工设备精度不足、装配操作不合理，齿轮齿面出现微观形貌不平整及局部缺陷等问题，进而导致齿轮对在啮合过程中出现刺耳的啸叫异响，其信号频谱通常表现为某一主要频率阶次能量突出。准确识别引起啸叫问题的阶次特征对指导齿轮微观修形，降低啸叫异音具有重要意义。

啸叫异响是一种典型的主观评价指标，通过人耳评定，人耳是一个强非线性系统，对不同的频率敏感度不同，即使相同能量，不同频率的声音听觉感受不同。直接通过1/3倍频程分析中心频率能量的大小来评判啸叫仅考虑噪声信号能量，忽略了人耳对频率的选择性，无法根据人耳主观感觉给出引起啸叫的频率阶次的主次关系。

技术实现要素：

为了克服上述现有技术的缺点，本发明的目的在于提供了一种基于音调能量和人耳频率选择性的减变速机啸叫评判方法，使低频区成为非敏感区，使高频区凸显其尖锐属性，从而有效模拟人耳对啸叫的主观感觉。

为了达到上述目的，本发明采取的技术方案为：

一种基于音调能量和人耳频率选择性的减变速机啸叫评判方法，包括以下步骤：

步骤1：使用噪声传感器高频采样减变速机噪声信号；

步骤2：对噪声信号进行a计权处理，并将处理后信号转换到对数域，然后进行傅里叶变换，得到其能量幅频曲线a(f)；a计权频响函数为：

其中：

c1＝3.5041384×10¹⁶

c2＝20.598997²

c3＝107.65265²

c4＝737.86223²

c5＝12194.217²

f—频率(hz)

步骤3：根据iso532-1计算噪声信号特征响度与临界频带率，然后计算每个特征响度的频带尖锐度s(z)，并对频带尖锐度s(z)进行归一化处理，频带尖锐度s(z)计算公式为：

其中，z代表临界频带率，n'(z)代表随临界频带率变化的特征响度，g(z)代表随临界频带率变化的尖锐度权值系数：

归一化处理为：

smin为s(z)中最小值，smax为s(z)中最大值；

步骤4：对步骤2得到的能量幅频曲线a(f)进行1/24倍频程高斯平滑处理，将处理后的数据作为背景噪声b(f)，选取{a(fi)＝a(fi)|a(f)-b(f)≥6}作为候选音调，其中a(fi)为满足上式的第i个候选音调；

步骤5：根据临界频带率z与频率f的对应关系映射a(fi)→a(zi)，计算步骤3中所得频带尖锐度s(z)加权因子:

其中，s(zi)是根据步骤3中计算出的频带尖锐度，n为选取的临界频带个数；令a(fi)＝a(fi)＝s·a(zi)，将加权后的a(f)作为噪声信号的能量幅频曲线，进行后续计算；

步骤6：根据a(fi)在频谱上对应的双边谷底位置确定音调带宽上下限[fi1,fi2]，计算音调带宽内的音调能量eti；

其中，ak代表能量幅频曲线的幅值；

步骤7：根据a(fi)在频谱上对应的中心频率fi，计算其对应的临界频带zi＝[fli,fhi]，计算临界频带内噪声总能量etotali；

其中，ak代表能量幅频曲线的幅值；

步骤8：计算背景噪声能量eni：

步骤9：计算音调能量与背景噪声能量的差eri＝eti-eni，并根据eri的大小来进行影响啸叫异响频率的主次排序。

所述的步骤1-步骤5的a计权与频带尖锐度加权，体现了人耳频率选择性。

所述的步骤6-步骤9的能量计算，用于量化不同频率音调对啸叫的贡献程度。

本发明的有益效果为：

本发明中a计权特性为：1000hz以内对噪声幅值具有较强的衰减作用，1000hz-6000hz对噪声幅值具有微弱的增益作用，6000hz-20000hz对噪声幅值具有较强的衰减作用；由于现实中减变速机的齿轮啮合频率一般不会高于6000hz，所以本发明方法主要应用了a计权在1000hz以内的低频衰减特性(低频选择性)。本发明中频带尖锐度随频率的增大而增大，将其加权到噪声信号中，对小于1000hz的低频成分有轻微的衰减作用，对大于1000hz的高频成分具有明显的增益效应，且随着频率的增加，增益效应越强，所以本发明方法主要应用了频带尖锐度的高频增益效应(高频选择性)。与传统1/3倍频程分析方法评判啸叫相比，本发明方法可有效模拟人耳对整个频带的选择性，并指出影响啸叫的主次频率，所得到的客观评价数值更加符合人耳主观感受。

附图说明

图1为本发明方法的流程图。

图2为本发明实施例实验变速器结构简图。

图3为本发明实施例噪声信号时域图、频谱图(对数域)和经过a计权后的频谱图(对数域)。

图4为本发明实施例噪声信号计算所得特征响度图。

图5为本发明实施例a计权后频谱图经过1/24倍频程后所得背景噪声曲线图(对数域)。

图6为本发明实施例所选音调最终评判结果图。

具体实施方式

下面结合附图和实施例对本发明进一步详细说明。

如图1所示，一种基于音调能量和人耳频率选择性的减变速机啸叫评判方法，包括以下步骤：

步骤1：通过噪声传感器高频采样减变速机噪声信号，其采样频率为40960hz，减变速机每个挡位在加速过程中均存在较为不同程度的啸叫异音；该减变速机结构简图与齿轮参数如图2所示，图2中标出了每个齿轮的齿数、输入轴、输出轴和中间轴，并在下方标出了各个齿轮副对应的挡位；选择啸叫较为严重的3挡，转速n＝2337rpm，其各转轴转频和挡位啮合频率如表1，

表1减变速机特征频率

步骤2：对噪声信号进行a计权处理，将处理后信号转换到对数域，然后进行傅里叶变换，得到能量幅频特性曲线a(f)；a计权频响函数为：

其中：

c1＝3.5041384×10¹⁶

c2＝20.598997²

c3＝107.65265²

c4＝737.86223²

c5＝12194.217²

f—频率(hz)

本实施例噪声信号、未进行a计权幅频曲线(对数域)和a计权后幅频曲线(对数域)如图3所示,从图3中可以看出a计权对1000hz以内的低频区能量幅值衰减较大，对1000hz以上高频区能量幅值影响微小，体现了a计权在低频区的主要选择特性；

步骤3：根据iso532-1计算噪声信号的特征响度与临界频带率，在此基础上，然后计算每个特征响度的频带尖锐度s(z)：

其中，z代表临界频带率，n'(z)代表随临界频带率变化的特征响度，g(z)代表随临界频带率变化的尖锐度权值系数：

根据公式将频带尖锐度s(z)归一化到[0,1]之间，作为后续音调能量加权因子，本实施例特征响度图如图4所示，图4中横坐标为临界频带率，纵坐标为特征响度，该图表示响度在单位临界频带率上的分布，是计算频带尖锐度的基础；

步骤4：对步骤2中能量幅频函数a(f)进行1/24倍频程高斯平滑处理，将处理后的数据作为背景噪声b(f)，选取作)为候选音调，其中a(fi)为满足上式的第i个候选音调；将b(f)、a(f)、候选音调绘制在同一张图上，如图5所示，图5中黑色为能量幅频曲线，灰色为背景噪声，黑点为选出的候选音调；

步骤5：依据临界频带率z与频率f的对应关系映射a(fi)→a(zi)，计算步骤3中所得频带尖锐度s(z)加权因子:

其中，s(zi)时根据步骤3中计算出的频带尖锐度，n为选取的临界频带个数，此处取n＝1，可根据不同效果进行更改；令a(fi)＝a(fi)＝s·a(zi)，将加权后的a(f)作为噪声信号的能量幅频曲线，进行后续计算；考虑频带尖锐度加权后能量幅频曲线如图6上图所示，图6中上图为频带尖锐度加权后的能量幅频曲线，对比图5与图6中能量幅频曲线可以看出频带尖锐度加权对1170hz和5834hz的音调能量幅值均有增益效应，但是对5834hz音调能量幅值增益效应更强，使原本幅值较大的1170hz音调变为幅值第二大，充分表现了频带尖锐度的高频选择特性；

步骤1-步骤5的a计权与频带尖锐度加权，体现了人耳频率选择性；

步骤6：根据a(fi)在频谱上对应的双边谷底位置确定音调带宽上下限[fi1,fi2]，计算音调带宽内的音调能量eti：

其中，ak代表能量幅频曲线的幅值；

步骤7：根据a(fi)在频谱上对应的中心频率fi，计算其对应的临界频带zi＝[fli,fhi]，计算临界频带内噪声总能量etotali：

其中，ak代表能量幅频曲线的幅值；

步骤8：计算背景噪声能量eni：

步骤9：计算音调能量与背景噪声能量的差eri＝eti-eni，并根据eri的大小来进行影响啸叫异响频率的主次排序，其结果如图6所示，由图6可以看出5834hz处音调的er值最大，1170hz处音调的er值第二大，所以该实例中5834hz处音调对啸叫的贡献最大，1170hz处音调对啸叫的贡献第二大。

步骤6-步骤9的能量计算，用于量化不同频率音调对啸叫的贡献程度。

本实施例经5人专家评审团的人耳主观评价检验，检验方式为：每个噪声样本中依次滤除候选音调成分，然后与没有滤除任何音调成分的噪声样本进行比较，主观判断出对啸叫贡献大小的音调成分的主次顺序。其检验结果如表2所示，从表2可以得出，在22个噪声样本中，有21个噪声样本的人耳评价结果与本发明方法判断结果一致，正确率达到95.45％，说明本发明方法具有较高的正确率。

表2

以上所述仅是本发明的优选实施例，应当指出，对于本技术领域的普通人员来说，在不脱离本发明原理的前提下，还可以做出若干改进，这些改进也应视为本发明的保护范围。