(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210839159.7 (22)申请日 2022.07.18 (71)申请人 安徽理工大 学 地址 232001 安徽省淮南市山 南新区泰丰 大街168号 (72)发明人 阮学云　邵良友　张馨月　贾世林　 李佳静　傅军杰　龚莹　王旭东　 (51)Int.Cl. G06F 30/20(2020.01) G06F 111/04(2020.01) G06F 119/10(2020.01) (54)发明名称 一种基于可靠度下多测点实测值的噪声源 声功率反推方法 (57)摘要 本发明涉及一种基于可靠度下多测点实测 值的噪声源声功率反推方法， 其包括的步骤为： 1)测量接受点的声压级数据； 2)构建声源声功率 反演数学模型； 3)可靠度分析； 4)使用Matlab优 化工具箱解算反演模型获得声源的声功率级。 上 述方法根据国际标准ISO9613 ‑2给出的户外传播 计算方法， 通过多测点噪声实测数据与理论计 算， 利用可靠度理论及优化思想， 通过不同可靠 度的设置， 识别并剔除测点数据中偏差较大的 点， 以此构建了声源声功率级反推的数学模型， 实现声功率级反演计算。 权利要求书2页 说明书6页 附图1页 CN 115081241 A 2022.09.20 CN 115081241 A 1.一种基于可靠度下多测点实测值的噪声源声功率反推方法， 包括以下实施步骤： 步骤1： 测量接受点的声压级数据； 步骤2： 构建声源声功率反演数 学模型； 步骤3： 可靠度分析； 步骤4： 使用Matlab优化工具箱解 算反演模型获得声源的声功率级。 2.根据权利要求1所述的一种基于可靠度下多测点实测值的噪声源声功率反推方法， 其特征在于： 所述步骤1中的测点选在敏感点附近， 且距离声源要足够远， 方便将声源等效 为点声源。 3.根据权利要求1所述的一种基于可靠度下多测点实测值的噪声源声功率反推方法， 其特征在于： 所述 步骤2中的声功率级和声压级的关系式为： Lp＝LW+Dc‑A                       式(1) 其中， Lw为点声源声功率,单位： dB(A)； Dc为声源指向性校正,单位： dB(A)； A为各类衰减， 其中包括大气衰减， 几何衰减， 地 面效应衰减， 屏蔽衰减等， 单位： dB(A)。 4.根据权利要求1所述的一种基于可靠度下多测点实测值的噪声源声功率反推方法， 其特征在于： 所述 步骤2中的声衰减量 根据噪声预测软件Cadna/A直接获得。 5.根据权利要求1所述的一种基于可靠度下多测点实测值的噪声源声功率反推方法， 其特征在于： 计算室外多个声源， 测点的声压级为各声源在此点的能量叠加。 具体表示 为： 其中， Lpj为测点j的声压级， 在本文中即表示测点实测值数据， 单位： dB(A)； LWi为声源i的声功率级， 单位： dB(A)； Aij表示为声源i到测点j的声衰减， 单位： dB(A)。 6.根据权利要求1所述的一种基于可靠度下多测点实测值的噪声源声功率反推方法， 其特征在于： 所述步骤3中的可靠度为实测值与预测 值误差在一定范围的测 点个数占总的 测点个数比例， 表达式为： R＝m'/m                              式(3) 其中， m'为实测值和预测值 误差在给定范围内的测点个数； m为总的测点个数。 7.根据权利要求6所述的可靠度， 其特征在于： 测点数据若是不存在坏点， 则可靠度可 选择100％， 反之存在坏点， 当可靠度选择100％时， 方程无解。 所述最高的可靠度设置与测 点中坏点个数的关系见表1。 表1： 坏点的个数与可设置的最高可靠度 坏点个数 1 2 3 ...s 可靠度(％) (m‑1)/m (m‑2)/m (m‑3)/m ...(m‑s)/m 8.根据权利要求1所述的一种基于可靠度下多测点实测值的噪声源声功率反推方法，权　利　要　求　书 1/2 页 2 CN 115081241 A 2其特征在于： 根据噪声源附近测点声压级 数据， 判断出声源声功 率级范围， 利用计算机在此 范围内随机取点， 结合matlab优化工具箱中fmincon()函数， 匹配出目标函数最小时的各 声源声功率级。 9.根据权利要求8所述的求解算法， 其特征在于： 将各测点误差的均差作为目标函数， 表达式为： 将各声源等效点声源并将其作为设计 变量； 约束条件为可靠度约束， 表达式为： 其中， Aij为声源i到测点j之间总的声衰减量， 单位： dB(A)； [R]为给定的许用可靠度； R为实际计算出的可靠度； m为测点个数； n为声源个数； LWimax为声源i声功率级给定的最大限值， LWimin为声源i声功率级给定的最小限值， 单位： dB(A)。权　利　要　求　书 2/2 页 3 CN 115081241 A 3