[发明专利]汽车车内双耳异响时变噪声响度评价方法

权利要求书

1.汽车车内双耳异响时变噪声响度评价方法，其特征在于：包括如下步骤，

步骤一、利用声音信号采集设备采集非稳态工况下车内目标点处的双耳噪声信号；

步骤二、在数值计算软件中对左右耳噪声信号分别进行瞬时特征响度计算；

步骤三、由步骤二得到的左右耳瞬时特征响度分别计算左右耳瞬时特征响度；在每个erb_m子频带内，比较第二帧的瞬时特征响度与第一帧的瞬时特征响度的大小，如果第二帧大于第一帧，按照公式(1)计算，否则按照公式(2)计算；从第三帧开始依次比较每帧的瞬时特征响度与前一帧的短期特征响度的大小，如果大于前一帧，按照公式(3)计算，否则按照公式(4)计算，即可分别得到左右耳短期特征响度；

S′₂＝α_αS₂+(1-α_α)S₁……(1)

S′₂＝α_rS₂+(1-α_r)S₁……(2)

S′_n＝α_αS_n+(1-α_α)S′_n-1……(3)

S′_n＝α_rS_n+(1-α_r)S′_n-1……(4)

其中，S′₂表示第二帧的短期特征响度，S₂表示第二帧的瞬时特征响度，S₁表示第一帧的瞬时特征响度，S′_n表示正在计算的这一帧信号的短期特征响度，S_n表示正在计算的这一帧信号的瞬时特征响度，S′_n-1表示前一帧信号的短期特征响度，α_α与每帧的时长和上升时间有关，α_r与每帧的时长和下降时间有关；

步骤四、将步骤三得到的左右耳短期特征响度分别与相应的高斯平滑函数进行卷积计算，分别得到左右耳平滑的短期特征响度，左耳平滑的短期特征响度由式(5)计算得到，右耳平滑的短期特征响度由式(6)计算得到；

其中N′_L(i)_smoothed为平滑后的左耳短期特征响度，N′_R(i)_smoothed为平滑后的右耳短期特征响度，N′_L为每帧内的左耳短期特征响度，N′_R为每帧内的右耳短期特征响度；i为正在进行平滑处理的子频带Cam值，D_i为正在计算的子频带Cam值与i之间的差值，当i+D_i小于人耳听阈下限或大于人耳听阈上限时，N′_L、N′_R为0，B为影响耳间抑制作用曲线形状的参数；

步骤五、根据公式(7)和(8)，由步骤四得到的左右耳平滑的短期特征响度，分别计算左右耳抑制因子；左耳短期特征响度除以左耳抑制因子，得到受抑制的左耳短期特征响度，右耳短期特征响度除以右耳抑制因子，得到受抑制的右耳短期特征响度，分别在人耳听阈下限至人耳听阈上限的范围内积分得到左右耳短期响度；

INH_L(i)＝2/[1+{sech(N′_R(i)_smoothed/N′_L(i)_smoothed)}^θ]……(7)

INH_R(i)＝2/[1+{sech(N′_L(i)_smoothed/N′_R(i)_smoothed)}^θ]……(8)

其中INH_L(i)为左耳抑制因子，INH_R(i)为右耳抑制因子，θ为幂指数；为了防止N′_L(i)_smoothed或N′_R(i)_smoothed为0时出现除以0的情况，将所述两个参数分别增加一个预设的很小的数a_m；

步骤六、在步骤五得到的左右耳短期响度的基础上，由式(9)～(12)分别计算得到左右耳长期响度；

沿着时间轴，首先比较第二帧的短期响度与第一帧的瞬时特征响度的大小，如果第二帧大于第一帧，按照公式(9)计算，否则按照公式(10)计算；从第三帧开始依次比较每帧的短期响度与前一帧的长期响度的大小，如果大于前一帧，按照公式(11)计算，否则按照公式(12)计算，即可分别得到左右耳长期响度；

S″_2l＝α_αlS′_2s+(1-α_αl)S_1s……(9)

S″_2l＝α_rlS′_2s+(1-α_rl)S_1s……(10)

S″_nl＝α_αlS′_ns+(1-α_αl)S″_(n-1)l……(11)

S″_nl＝α_rlS′_ns+(1-α_rl)S″_(n-1)l……(12)

其中，S″_2l表示第二帧的长期响度，S′_2s表示第二帧的短期响度，S_1s表示第一帧的瞬时特征响度，S″_nl表示正在计算的这一帧信号的长期响度，S′_ns表示正在计算的这一帧信号的短期响度，S″_(n-1)l表示前一帧信号的长期响度，α_αl与每帧的时长和上升时间有关，α_rl与每帧的时长和下降时间有关；

步骤七、对步骤五得到的左右耳短期响度相加得到双耳短期响度，即可利用得到的双耳短期响度实现对车内双耳异响时变噪声响度详细变化过程的评价；对步骤六得到的左右耳长期响度相加得到双耳长期响度，双耳长期响度的最大值即为此噪声样本在采样时间段内的双耳总响度，即能够利用得到的双耳长期响度实现对车内双耳异响时变噪声响度整体感知的评价。