一种基于可靠度下多测点实测值的噪声源声功率反推方法技术

本发明专利技术涉及一种基于可靠度下多测点实测值的噪声源声功率反推方法，其包括的步骤为：1)测量接受点的声压级数据；2)构建声源声功率反演数学模型；3)可靠度分析；4)使用Matlab优化工具箱解算反演模型获得声源的声功率级。上述方法根据国际标准ISO9613

【技术实现步骤摘要】
一种基于可靠度下多测点实测值的噪声源声功率反推方法


[0001]本专利技术涉及噪声预测领域，具体是一种通过可靠度将用于反推声源声功率实测点中的坏点剔除以此提高工业企业噪声预测精度的方法。

技术介绍

[0002]随着城市和工业化进程的推进，我国环境问题日益凸显，工业噪声已经成为一大焦点。在工业企业噪声预测及控制中，方便、准确地获得工业设备的声功率级，提高噪声预测精度尤为重要。目前，在预测模型实际搭建中，声源声功率级大多根据声源附近测点声压级与实测值反复对比来调整估算，不能确保所有测点预测值与实测值之间的误差在一定范围内，均是人为确定其中一些测点满足误差范围，而一些误差较大的点部分不考虑，该方法无科学依据，随意性较大，噪声预测精度难以得到保证。获得声源声功率级的理论方法主要有声压法、声强法、振速法：其中声强法、振速法对操作要求高，不适用于户外实际工程应用；而声压法应用最为普通，但在获得应用于声压法所需的测点声压级时，难免有的测点数据会由于环境等因素存在偏差，导致最后结果不正确，需要将这些存在误差的测点剔除。因此提供一种基于可靠度下多测点实测值的噪声源声功率反推方法显得尤为重要。

技术实现思路

[0003]本专利技术的目的在于：针对上述噪声源声功率级难以准确获取问题，提供一种基于可靠度下多测点实测值的噪声源声功率反推方法，从而提高噪声预测精度。
[0004]本专利技术实现专利技术目的采用如下技术方案：
[0005]一种基于可靠度下噪声源声功率反推方法，包括以下实施步骤：
[0006]步骤1：测量接受点的声压级数据；
[0007]步骤2：构建声源声功率反演数学模型；
[0008]步骤3：可靠度分析；
[0009]步骤4：使用Matlab优化工具箱解算反演模型获得声源的声功率级。
[0010]作为优选，合理布置测点，在设备运行稳定时，使用噪声采集设备声级计采集测点的声压级，注意记录采集数据时温度、风速、测点高度。
[0011]作为优选，声功率级的反推原理根据ISO 9613
‑
2户外声传播方法，具体表示为：
[0012]L
p
＝L
W
+D
c
‑
A
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式(1)
[0013]其中，L
w
为点声源声功率，单位：dB(A)；
[0014]D
c
为声源指向性校正，单位：dB(A)；
[0015]A为各类衰减，其中包括大气衰减，几何衰减，地面效应衰减，屏蔽衰减等，单位：dB(A)。
[0016]作为优选，计算室外多个声源，测点的声压级为各声源在此点的能量叠加。具体表示为：
[0017][0018]其中，L
pj
为测点j的声压级，在本文中即表示测点实测值数据，单位：dB(A)；
[0019]L
Wi
为声源i的声功率级，单位：dB(A)；
[0020]A
ij
表示为声源i到测点j的声衰减，单位：dB(A)。
[0021]作为优选，在实际工程应用中，由于噪声预测软件Cadna/A的计算准则参考9613
‑
2标准，而衰减量和几何模型相关，因此通过Cadna/A软件直接获取，可大大减少计算工作量。
[0022]作为优选，根据噪声源附近测点声压级数据，判断出声源声功率级范围，利用计算机在此范围内随机取点，结合matlab优化工具箱中fmincon()函数，匹配出目标函数最小时的各声源声功率级。
[0023]作为优选，测量设备精度、声衰减计算、噪声预测模型中的房屋以及地形设置等导致结果存在偏差，但这些因素导致的误差在工程计算中允许存在，一般为2
‑
3dB。
[0024]作为优选，可靠度为实测值与预测值误差在一定范围的测点个数占总的测点个数比例，表达式为：
[0025]R＝m'/m
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式(3)
[0026]其中，m'为实测值和预测值误差在给定范围内的测点个数，m为总的测点个数。
[0027]作为优选，以测点声压级实测值和根据求解出的声源声功率级反推出测点声压级预测值的均差为目标函数，表达式为：
[0028][0029]作为优选，将所求的点声源作为设计变量。
[0030]作为优选，约束条件为可靠度约束，表达式为：
[0031][0032]s.t.
[0033]g1(X)＝[R]‑
R≤0
[0034]g2(X)＝L
W1min
≤x1≤L
W1max
[0035]g3(X)＝L
W2min
≤x2≤L
W2max
[0036]…
[0037]g
n+1
(X)＝L
Wnmin
≤x
n
≤L
Wnmax
[0038]其中，A
ij
为声源i到测点j之间总的声衰减量，单位：dB(A)；
[0039][R]为给定的许用可靠度；
[0040]R为实际计算出的可靠度；
[0041]m为测点个数；
[0042]n为声源个数；
[0043]L
Wimax
为声源i声功率级给定的最大限值，L
Wimin
为声源i声功率级给定的最小限值，单位：dB(A)，为常量。
[0044]作为优选，测点数据若是不存在坏点，则可靠度可选择100％，反之存在坏点，当可靠度选择100％时，方程无解。根据上节中提到的m个测点,n个声源，表1列出存在坏点的个数与可设置的最高可靠度。
[0045]表1：坏点的个数与可设置的最高可靠度
[0046]坏点个数123...s可靠度(％)(m
‑
1)/m(m
‑
2)/m(m
‑
3)/m...(m
‑
s)/m
[0047]本专利技术与现有技术相比，其有益效果体现在：
[0048]本专利技术基于ISO 9613
‑
2户外声传播计算方法，提出可靠度计算方法以及最大可靠度设置，构建反推声功率的数学模型。通过不同可靠度设置，本专利技术可自动识别并剔除坏点，最大的减少有效测点与实测值的综合误差，从而得到较为准确的声源声功率级。
附图说明
[0049]图1是本专利技术实施例中提供的一种基于可靠度下多测点实测值的噪声源声功率反推方法流程图；
[0050]图2是本专利技术实施例中测点布置示意图；
[0051]图3是本专利技术声源与测点数目不同的声学模型示意图。
具体实施方式
[0052]以下结合具体实施例和说明书附图对本专利技术具体、详细地解释说明。
[005本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于可靠度下多测点实测值的噪声源声功率反推方法，包括以下实施步骤：步骤1：测量接受点的声压级数据；步骤2：构建声源声功率反演数学模型；步骤3：可靠度分析；步骤4：使用Matlab优化工具箱解算反演模型获得声源的声功率级。2.根据权利要求1所述的一种基于可靠度下多测点实测值的噪声源声功率反推方法，其特征在于：所述步骤1中的测点选在敏感点附近，且距离声源要足够远，方便将声源等效为点声源。3.根据权利要求1所述的一种基于可靠度下多测点实测值的噪声源声功率反推方法，其特征在于：所述步骤2中的声功率级和声压级的关系式为：L
p
＝L
W
+D
c
‑
A
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式(1)其中，L
w
为点声源声功率,单位：dB(A)；D
c
为声源指向性校正,单位：dB(A)；A为各类衰减，其中包括大气衰减，几何衰减，地面效应衰减，屏蔽衰减等，单位：dB(A)。4.根据权利要求1所述的一种基于可靠度下多测点实测值的噪声源声功率反推方法，其特征在于：所述步骤2中的声衰减量根据噪声预测软件Cadna/A直接获得。5.根据权利要求1所述的一种基于可靠度下多测点实测值的噪声源声功率反推方法，其特征在于：计算室外多个声源，测点的声压级为各声源在此点的能量叠加。具体表示为：其中，L
pj
为测点j的声压级，在本文中即表示测点实测值数据，单位：dB(A)；L
Wi
为声源i的声功率级，单位：dB(A)；A
ij
表示为声源i到测点j的声衰减，单位：dB(A)。6.根据权利要求1所述的一种基于可靠度下多测点实测值的噪声源声...

【专利技术属性】
技术研发人员：阮学云，邵良友，张馨月，贾世林，李佳静，傅军杰，龚莹，王旭东，
申请(专利权)人：安徽理工大学，
类型：发明
国别省市：