人体感知的声音和振动的频率计权

TestGuru

人体感知的声音和振动的频率计权

一、简介

人总是暴露于各种声音和振动中，其来源包括家用电器、交通系统、机械设备、手持电动工具，以及打桩和爆破等工业活动。对噪声和振动的过度暴露可能引起身体不适或健康问题，例如听力丧失或脊柱疾病。因此，必须按照相关国际或国家标准，来正确地测量声音和振动，并评估其风险。

人对声音和振动的感知并不随频率呈线性变化。人耳对从500Hz到8kHz之间的声音频率最敏感，而对更低和更高的声音频率的敏感度较弱。为了将测得的声音与人耳对声音的响应关联起来，可采用频率计权。IEC 61672定义了三种常用的计权函数：A、C和Z。相似地，人体也是对一些振动频率更为敏感，而对另一些振动频率则不那么敏感。因此，应对测得的数据进行“计权”，以给人更敏感的频率更高的权重。与声音不同，人体对振动的敏感度还取决于振动的方向和身体的姿势。因此，ISO 2631和ISO 5349规定了多种频率计权函数，以将这些因素考虑在内。

二、声音的频率计权

人对声压的主观感受可用响度来描述。人耳的敏感度不仅取决于声音频率，还取决于声级大小，如下图中的等响度曲线所示。总体上，随着声级的提高，等响度曲线变得更加平坦。

等响度曲线

A计权大致是在1 kHz穿过40 dB（即40 phon的响度级）的等响度曲线的倒置，通常用于测量环境噪音并评估其对人感知的影响；而C计权则基于在1 kHz穿过100 dB（即100 phon的响度级）的等响度曲线，通常用于测量高声级的工业噪音。与A和C计权不同，Z计权在10Hz至20 kHz之间具有平坦的频率响应，它用于未计权的测量。还有一些其它的频率计权方案，例如B、D、G和ITU-R468计权。B计权基于70 phon的等响度曲线。D计权专用于采用非旁路型喷气发动机的军用飞机。G计权用于8Hz到约40 Hz的次声波测量。ITU-R 468噪声计权是为了能更准确地反映人耳对各种噪声而非纯音的主观感受而设计的，它在欧洲得到广泛采用，特别是在广播行业中。下图显示的是最常用的A、C和Z计权曲线。

A、C 和 Z 计权曲线

TestGuru

IEC 61672和ANSI S1.43均提供了C和A频率计权的解析表达式，单位为dB，如下所示。C计权由频率为f1的两个低频极点、频率为f4的两个高频极点、和0 Hz的两个零点实现。通过这些极点和零点，相对于参考频率fr = 1 kHz，C计权的功率响应将在fL = 10^1.5 Hz和fH = 10^3.9 Hz处下降1/2（约为-3dB）。A计权通过在C计权的基础上添加两个耦合的一阶高通滤波器来实现。这相当于在0 Hz处添加两个零点，并在频率f2和f3处各添加一个极点。对于每个高通滤波器，截止频率由fA = 10^2.45 Hz给出。

其中：

fr= 1000 Hz，fA = 10^2.45 Hz，fL = 10^1.5 Hz，fH = 10^3.9 Hz。A1000和C1000 为归一化常数，分别用于将A和C计权在1000Hz处的增益置于0dB。由上面的公式可计算出f1 = 20.60 Hz， f2 = 107.7 Hz，f3 = 737.9 Hz，f4 = 12194 Hz，C1000 = -0.062 dB，A1000 = -2.000 dB。

ANSI S1.42提供了s域传递函数，以用模拟方法实现C和A计权的上述幅频响应。

其中：ω1= 2πf1，ω2 = 2πf2，ω3= 2πf3，ω4 = 2πf4，Gc =100^0.062/20，GA = 10^2/20。

下表列出了按前面公式计算出来的在10Hz到20kHz之间的1/3倍频带中心频率处（f = fr ×10^(n-30)/10), n = 10,11, …, 43）的A、C、Z计权值。

1/3倍频带的A、C和Z计权值

需要注意的是，上面提及的每种频率计权方案只是频率的函数，并未考虑到听觉灵敏度随声级变化的特性以及听觉掩蔽效应。这些频率计权函数相对较容易实现，在各种标准中有明确定义，因此得到了广泛应用。然而，长期以来人们知道，这些传统的声学度量指标并不能很好地反映人耳感觉到的声音品质。这里说的声音品质指的是人们听到声音后舒服或不舒服的整体感觉，或者对发出这种声音的设备的声学品质以外的质量的感觉，例如对引擎的功率、结构的稳定性的感觉等等。对于一些要求严格的应用场景，需要更高级的声音品质度量指标，例如响度、响度级、尖锐度和语言清晰度指数等。这些指标基于更复杂的心理声学模型，将另行讨论。

TestGuru

三、振动的频率计权

人体振动关注的是机械振动对人体的影响。人体振动主要分为两种类型：全身振动和手传振动。全身振动以整个人体为传导对象，通常通过支撑表面传输：站者的脚部、坐者的臀部、背部和脚部，或躺者的支撑区域。这种振动常见于车辆、机械设备、建筑物内部以及工作机械周边环境。最突出的例子包括越野车、喷气艇、直升机和喷气式飞机。手传振动传输到手和手臂，常见于手持电动工具的操作者。全身系统和手臂系统 “在机械特性上有所不同”，因此需分别研究。人体可视为由一系列的弹簧-质量-阻尼系统构成，每个系统具有自身的固有频率。下图是一个人体的简化机械模型。

人体的机械模型及各部位的共振频率范围

该模型显示人体的不同部位具有不同的共振频率。某些振动频率可能对某些身体部位造成的损害更大。例如，20Hz到30Hz之间的振动可能会引起头部和肩部之间的共振，而振动可被放大至350%！人体是一个有强大阻尼的系统，因此，当其某部分被激发到其固有频率时，它将在一个频率范围内共振，而不是在单一频率上。

人体对机械振动的敏感度取决于振动频率、人体暴露位置、平移或直线振动的方向、旋转振动的轴以及人体的姿势。人们设计了各种频率计权函数来将这些因素考虑在内。下图显示了人体和手部的基本中心坐标系。

人体的基本中心轴

手的基本中心轴

在评估振动对健康、舒适度、感知和运动病的影响时，根据不同位置、姿势和方向，应采用不同的计权函数，见下表。

全身振动有三种主要的频率计权：Wk、Wd和Wf。Wk适用于Z方向和躺姿的竖直向（头部除外），而Wd适用于X和Y方向以及躺姿的水平向。这两种计权都用于评估与健康、舒适度和感知相关的全身振动。然而，对于与运动病相关的评估，应该采用Wf，它侧重于低频全身振动。

Wc、We和Wj是针对全身振动的特殊情况的附加计权，分别适用于座椅靠背振动、旋转振动和躺姿的头下振动。

Wm是针对建筑物中全身振动和冲击的频率计权，涉及到居住者的舒适感和烦扰感，无需定义居住者的姿势。

Wb是用于评估具有固定导轨的运输系统中乘客和机组人员舒适度的全身振动频率计权。

Wh是用于评估振动工具、机械或工件的操作者的舒适感、熟练度和健康状况的手传振动频率计权。

总体来说，上面提到的所有频率计权都可以通过级联四个模拟滤波器来实现：一个高通滤波器、一个低通滤波器、一个加速度-速度过渡滤波器和一个向上阶跃滤波器。这些在相应的ISO标准以及ISO8041中有定义。

带限是通过高通滤波器和低通滤波器的组合来实现的。它们都是二阶Butterworth滤波器，其传递函数如下。

高通滤波器：

低通滤波器：

加速度-速度过渡滤波器在低频时与加速度成正比，在高频时与速度成正比。其传递函数如下。

加速度-速度过渡滤波器：

注意：当f3和f4 （ω3和ω4）为无穷大时，Ht(s) = K。

向上阶跃滤波器的斜率大约为每倍频程6dB，与加速度变化率成正比。其传递函数如下。

向上阶跃滤波器：

注意：当f5和f6 （ω5和ω6）为无穷大时，Hs(s) = 1。

总的频率计权传递函数是以上四个传递函数的乘积，即：

H(s) = Hh(s).Hl(s).Ht(s).Hs(s)

其中：s为拉普拉斯变换的变量。Hh(s).Hl(s)代表带限传递函数，除Wf、Wh和Wm之外，其余的计权皆采用相同的带限设置。Ht(s).Hs(s)代表某应用的实际的频率计权传递函数。Wj 的Ht(s) = 1，而Wc、Wd、We、Wm、Wh 的Hs(s) = 1。令s = jw = j2πf,可得总的频率响应如下：

H(jw) = Hh(jw).Hl(jw).Ht(jw).Hs(jw)

则幅频响应为：

|H(jw)| = |Hh(jw)|.|Hl(jw)|.|Ht(jw)|.|Hs(jw)|

或

|H(f)| = |Hh(f)|.|Hl(f)|.|Ht(f)|.|Hs(f)|

可推得：

上述传递函数的参数列于下表。

振动频率计权的传递函数的参数

下表列出了按前面公式计算出来的在0.02Hz到4kHz之间的1/3倍频带中心频率处的频率计权值。

1/3倍频带的Wk、Wd、Wf、Wc、We、Wj、Wm、Wb和Wh的计权值

下图显示的是相应的计权曲线。

Wk、Wd、Wf、Wc、We、Wj、Wm、Wb和Wh计权曲线

TestGuru

四、用于评估声音和振动的经频率计权的参数


4.1 经频率计权的有效值（RMS）

经频率计权的有效值（Root-Mean-Square,或称RMS或均方根值）由线性时间平均公式定义如下：

其中：aw(t)为振动测量中经频率计权后的沿某一指定方向或绕某一指定轴的平移或旋转加速度，或为声音测量中经频率计权后的声压。它是时间t的函数，而T表示测量的持续时间。

此数值是振动的基本评估标准。但是，当经频率计权后的振动峰值因数超过9时，它可能会低估振动对人体的影响。这时应采用下述的替代指标之一作为附加评估参数 - 移动有效值或四次方振动剂量 (Vibration Dose Value，VDV)。

4.2 经频率计权的移动有效值（RMS）

经频率计权的移动有效值由线性时间平均公式定义如下：

其中：T 为移动平均的积分时间，ξ 为积分时间哑变量，t 为观测时间（瞬时时间）。上述的线性时间平均可由下面的指数时间平均（时间计权）来近似。

其中：(Tao)为时间常数。当应用于与远小于(Tao)的短时冲击时，二者结果的差异非常小；而当应用于较长时间的冲击和瞬态信号时，则差异稍大（可达 30%）。

4.3 经频率计权的分贝值

经频率计权的分贝值可计算如下：

其中：A0 为振动测量的参考加速度（ISO1683定义为10-6 m/s^2 的平移加速度），或者为声音测量的参考声压（空气中20 μPa，水中1 μPa）。

在声级测量中，当Aw为线性时间平均值时，Lw称为等效连续声级，通常在下标的权重符号后加上 “T” 或 “eq” 来表示。例如，LAT 和 LAeq 都表示 A 计权的等效连续声级。Aw 也可以是带时间常数的指数时间平均（或时间计权）值，即 Awtao。通常有三种标准的时间计权：I（冲击）、F（快速）和 S（慢速），其时间常数分别为 35 ms、125 ms 和 1 s。在这种情况下，应在 Lw 上加上时间计权符号。例如，LAF 表示A 计权“快速”时间计权声级。

4.4 最大瞬态值

最大瞬态振动值(Maximum Transient Vibration Value, MTVV)定义为振动加速度的移动有效值的最大值，即：MTVV= max[AwT(t)]或MTVV= max[Awtao(t)]，其积分时间T或Tao 为1 s （对应于声级测量中的“慢速”积分时间)。当MTVV/Aw>1.5时，MTVV应作为附加评估参数报告。

作为比较，在声音测量中，通常会测量最大的时间平均声级或时间计权声级，例如LAFmax、LASmax、LCFmax、和LCSmax。

4.5 运动病剂量（MSDV）和声暴露量

运动病剂量(Motion Sickness Dose Value, MSDV)和声暴露量都可由下式定义:

其中：T 为运动或声暴露持续的总时间。

在声音测量中，声暴露量常常通过下式转化为声暴露级：

其中：A0为参考声压（空气中为20 μPa），T0 = 1s。

4.6 振动剂量（VDV）

振动剂量(Vibration Dose Value, VDV)由下式定义：

其中：T为振动暴露的总时间。当VDV/(AwT1/4)>1.75时，VDV应作为附加评估参数报告。

振动剂量VDV可通过如下线性时间平均得到经频率计权后的均四次方根值(Root-Mean-Quad，RMQ)：

RMQ对峰值比有效值RMS更灵敏。

4.7 峰值

这里的峰值指的是经过频率计权后的加速度或声压的瞬时值（或正或负）的绝对值的最大值。

在声音测量中，峰值声压常常通过参考声压转化为峰值声级，例如LCpeak 表示C计权的峰值声级。

4.8 振动总值

振动总值汇总了平移振动的三个轴向振动，由下式定义：

其中：Awx、Awy和Awz分别为三个相互垂直的轴x、y和z上的经过频率计权的加速度有效值；kx、ky和 kz为乘数因子，其值取决于具体的测量应用。例如，在手传振动中， kx = ky= kz = 1；在评估振动对坐者的健康影响时，

X轴：Wd，kx = 1.4，

Y轴：Wd，ky = 1.4，

Z轴：Wk，kz = 1。