基于拉普拉斯范数的快速迭代收缩阈值的声源识别方法.pdf

《基于拉普拉斯范数的快速迭代收缩阈值的声源识别方法.pdf》由会员分享，可在线阅读，更多相关《基于拉普拉斯范数的快速迭代收缩阈值的声源识别方法.pdf（14页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201911348842.5 (22)申请日 2019.12.24 (71)申请人 重庆大学 地址 400044 重庆市沙坪坝区沙正街174号 (72)发明人 黄琳森徐中明李怡张志飞 贺岩松 (74)专利代理机构 重庆缙云专利代理事务所 (特殊普通合伙) 50237 代理人 王翔 (51)Int.Cl. G01S 5/22(2006.01) (54)发明名称 一种基于拉普拉斯范数的快速迭代收缩阈 值的声源识别方法 (57)摘要 本发明公开了一种基于拉普拉斯范数的快 速迭代收缩阈值。

2、的声源识别方法， 主要步骤如 下： 1)搭建基于拉普拉斯范数的快速迭代收缩阈 值的声源识别系统； 2)每个麦克风分别监测K个 等效声源的时域模拟声压信号， 并发送至多通道 的信号采集器； 3)多通道的信号采集器将接收到 的时域模拟声压信号转换为数字声压信号p， 并 发送至数据处理器； 4)基于接收到的数字声压信 号p， 数据处理器建立声源和麦克风阵列之间的 传递矩阵A； 5)数据处理器建立声场声源分布模 型； 6)数据处理器利用快速迭代收缩阈值算法对 声场声源分布模型进行迭代解算， 得到声源识别 结果x。 该算法有效的克服了低频范围( 0 500Hz)声源识别的不稳定性， 扩大了声源识别的 频。

3、率范围。 权利要求书2页 说明书7页 附图4页 CN 111257833 A 2020.06.09 CN 111257833 A 1.一种基于拉普拉斯范数的快速迭代收缩阈值的声源识别方法， 其特征在于， 主要包 括以下步骤： 1)搭建基于拉普拉斯范数的快速迭代收缩阈值的声源识别系统， 主要包括麦克风阵 列、 多通道的信号采集器和数据处理器； 所述麦克风阵列包括分布在声源检测空间内的M个 麦克风。 2)每个麦克风分别监测K个等效声源的时域模拟声压信号， 并发送至多通道的信号采 集器； 所述K个等效声源为随机分布于声源平面周围的等效声源； 3)所述多通道的信号采集器将接收到的时域模拟声压信号转换为。

4、数字声压信号p， 并 发送至数据处理器； 4)基于接收到的数字声压信号p， 所述数据处理器建立声源和麦克风阵列之间的传递 矩阵A。 5)所述数据处理器建立声场声源分布模型； 6)所述数据处理器利用快速迭代收缩阈值算法对声场声源分布模型进行迭代解算， 得 到声源识别结果x。 2.根据权利要求1所述的一种基于拉普拉斯范数的快速迭代收缩阈值的声源识别方 法， 其特征在于， 建立声源和麦克风阵列之间的传递矩阵A的主要步骤如下： 1)确定第m个麦克风测量到的声压信号p(m)如下所示： 式中， m1,2,M；为自由场格林函数； k为波数,为声源到全息 面的距离； qk为假想的单极子声源强度； 2)将公式(。

5、1)转化成向量-矩阵形式， 即： pAq； (2) 式中， A为MK声场传递矩阵； p为M维测量声压向量， 其元素为对应传感器在单一频率f 对应的频域响应； q为K维声源强度列向量， 其元素成分代表声源平面中对应点处的强度分 布。 3.根据权利要求1或2所述的一种基于拉普拉斯范数的快速迭代收缩阈值的声源识别 方法， 其特征在于， 声场声源分布模型的目标函数如下所示： 式中， |2指的是l2范数； |1指的是l1范数； 为Tikhonov正则化参数； 其中， 正则化参数 如下所示： 式中， SNR为信噪比； smax为声场传递矩阵A奇异值SVD分解后的最大值。 4.根据权利要求3所述的一种基于拉。

6、普拉斯范数的快速迭代收缩阈值的声源识别方 法， 其特征在于， 声场传递矩阵A奇异值SVD分解如下： 权利要求书 1/2 页 2 CN 111257833 A 2 sdiag(S)； (5) 式中， diag为将向量转化成对角矩阵形式。 5.根据权利要求1所述的一种基于拉普拉斯范数的快速迭代收缩阈值的声源识别方 法， 其特征在于， 利用快速迭代收缩阈值算法对声场声源分布模型进行迭代解算的主要步 骤如下： 1)计算第n次最速下降值y(n)， 即： 式中， L为Lipschitz常数； L2smax2； 迭代初始值y(1)0； t(1)1； y(n)为x(n-1)和x(n-2) 的线性连接方式； n。

7、为迭代次数； 梯度函数T表示矩阵的转 置； 其中， 线性连接方式y(n+1)满足下式： 其中， 迭代步长t满足下式： 2)设定拉普拉斯范数阈值 ； 当第n次迭代的声源识别结果|y(n)| /(2*L)时， 输出第n次迭代的声源识别结果x(n)， 反之进入步骤3； 3)利用拉普拉斯范数对第n次迭代的声源识别结果x(n)进行更新， 并返回步骤1； 第n+1 次迭代的声源识别结果x(n+1)如下所示： 式中， 为拉普拉斯范数的阶数。 权利要求书 2/2 页 3 CN 111257833 A 3 一种基于拉普拉斯范数的快速迭代收缩阈值的声源识别方法 技术领域 0001 本发明涉及声源识别领域， 具体是。

8、一种基于拉普拉斯范数的快速迭代收缩阈值的 声源识别方法。 背景技术 0002 当前基于麦克风阵列的两种主流声源识别方法分别为波束形成与声全息。 其中波 束形成的中远距离声源识别能力强于声全息， 而声全息在低频范围内的识别性能则更优于 波束形成， 二者各有千秋， 相辅相成。 声全息技术是近年来快速发展的一种声源识别方法， 基本原理是在紧靠被测声源物体表面的测量面上记录噪声数据， 然后通过空间声场变换算 法重构空间声场， 它通过在非常靠近声源的二维或三维面上测量声压， 重建出声音在三维 空间传播的声学量， 如声压、 声强、 声功率。 声全息在振动、 声学辐射和周边媒介中的声学能 量流之间建立起直接。

9、的关联。 由于其具有灵活的动态显示范围、 分辨率高等优点， 因此近年 受到研究人员以及从业人员的广泛关注， 其应用范围也从航空航海领域逐渐发展过渡到汽 车领域， 与此同时衍生出了各种不同的近场声全息算法， 如宽带声全息、 快速宽带声全息、 两步迭代收缩阈值等算法， 不断地拓展着其应用范围和领域。 0003 现有技术中存在一种基于麦克风阵列的声源识别方法， 公开了一种两步迭代收缩 阈值算法(MTwIST)。 其核心思想是在靠近声源面较近的等效源平面构造一系列假想的等效 源， 由这些等效源来替代实际声源分布。 通过声源识别理论， 建立一个等效源面到麦克风阵 列的声学传递方程， 并通过传统的梯度下降。

10、算法求解该传递方程， 从而迭代求解等效源强 度， 然后通过该等效源强度， 基于声传播过程重构在等效源面与阵列面之间的声源分布， 而 且该方法在等效源求解过程中， 引入一个迭代滤波过程， 从而保证迭代过程的收敛效率与 精准度。 该算法由于引入梯度下降算法和滤波过程， 能够提高传统等效源方法在中低频的 分辨率， 同时该算法也能够实现与其他声全息算法之间的结合或者切换， 从而实现在中低 频声源的重构。 然而， 该算法缺陷在于在低频范围(0500Hz)重建性能不稳定， 存在较大误 差。 发明内容 0004 本发明的目的是解决现有技术中存在的问题。 0005 为实现本发明目的而采用的技术方案是这样的， 。

11、一种基于拉普拉斯范数的快速迭 代收缩阈值的声源识别方法， 主要包括以下步骤： 0006 1)搭建基于拉普拉斯范数的快速迭代收缩阈值的声源识别系统， 主要包括麦克风 阵列、 多通道的信号采集器和数据处理器。 所述麦克风阵列包括分布在声源检测空间内的M 个麦克风。 0007 2)每个麦克风分别监测K个等效声源的时域模拟声压信号， 并发送至多通道的信 号采集器。 所述K个等效声源为随机分布于声源平面周围的等效声源。 0008 3)所述多通道的信号采集器将接收到的时域模拟声压信号转换为数字声压信号 说明书 1/7 页 4 CN 111257833 A 4 p， 并发送至数据处理器。 0009 4)基于。

12、接收到的数字声压信号p， 所述数据处理器建立声源和麦克风阵列之间的 传递矩阵A。 0010 进一步， 建立声源和麦克风阵列之间的传递矩阵A的主要步骤如下： 0011 4.1)确定第m个麦克风测量到的声压信号p(m)如下所示： 0012 0013式中， m1,2,M。为自由场格林函数。 k为波数,为声源 到全息面的距离。 qk为假想的单极子声源强度。 0014 4.2)将公式(1)转化成向量-矩阵形式， 即： 0015 pAq。 (2) 0016 式中， A为MK声场传递矩阵。 p为M维测量声压向量， 其元素为对应传感器在单一 频率f对应的频域响应。 q为K维声源强度列向量， 其元素成分代表声源。

13、平面中对应点处的强 度分布。 0017 5)所述数据处理器建立声场声源分布模型。 0018 进一步， 声场声源分布模型的目标函数如下所示： 0019 0020 式中， |2指的是l2范数。 |1指的是l1范数。 为Tikhonov正则化参数。 0021 其中， 正则化参数 如下所示： 0022 0023 式中， SNR为信噪比。 smax为声场传递矩阵A奇异值SVD分解后的最大值。 0024 声场传递矩阵A奇异值SVD分解如下： 0025 sdiag(S)。 (5) 0026 式中， diag为将向量转化成对角矩阵形式。 0027 6)所述数据处理器利用快速迭代收缩阈值算法对声场声源分布模型进。

14、行迭代解 算， 得到声源识别结果x。 0028 进一步， 利用快速迭代收缩阈值算法对声场声源分布模型进行迭代解算的主要步 骤如下： 0029 6.1)计算第n次最速下降y(n)， 即： 0030 0031 式中， L为Lipschitz常数。 L2smax2。 迭代初始值y(1)0； t(1)1。 y(n)为x(n-1)和x (n-2)的线性连接方式。 n为迭代次数。 梯度函数 T表示矩阵的 转置。 说明书 2/7 页 5 CN 111257833 A 5 0032 其中， 线性连接方式y(n+1)满足下式： 0033 0034 其中， 迭代步长t满足下式： 0035 0036 2)设定拉普拉。

15、斯范数阈值 。 当第n次迭代的声源识别结果|y(n)| /(2* L)时， 输出第n次迭代的声源识别结果x(n)， 反之进入步骤3。 0037 3)利用拉普拉斯范数对第n次迭代的声源识别结果x(n)进行更新， 并返回步骤1。 第 n+1次迭代的声源识别结果x(n+1)如下所示： 0038 0039 式中， 为拉普拉斯范数的阶数。 0040 本发明的技术效果是毋庸置疑的。 本发明提供了一种声源重建精度高、 稳定性较 好、 误差较小的声源识别算法LFISTA。 本发明首先基于声源和麦克风阵列构造传递矩阵， 再 将其转化为1-范数的最小化问题， 最后通过快速迭代收缩阈值， 结合拉普拉斯范数， 构造新。

16、 的声源识别算法。 该算法在运用快速迭代收缩阈值的高分辨率优势的同时， 能够充分利用 拉普拉斯范数的性质， 有效的克服了两步迭代收缩阈值算法(MTwIST)在低频范围(0 500Hz)的不稳定性， 以及扩大了声源识别的频率范围。 附图说明 0041 图1为声源重构示意图； 0042 图2为理论重建图； 0043 图3为传统Tikhonov方法重建图； 0044 图4为两步迭代收缩阈值算法MTwIST方法重建图； 0045 图5为基于拉普拉斯范数的快速迭代收缩阈值算法的声源识别方法LFISTA方法重 建图； 0046 图6为重建误差曲线图。 具体实施方式 0047 下面结合实施例对本发明作进一步。

17、说明， 但不应该理解为本发明上述主题范围仅 限于下述实施例。 在不脱离本发明上述技术思想的情况下， 根据本领域普通技术知识和惯 用手段， 做出各种替换和变更， 均应包括在本发明的保护范围内。 0048 实施例1： 0049 参见图1至图2， 一种基于拉普拉斯范数的快速迭代收缩阈值的声源识别方法， 主 说明书 3/7 页 6 CN 111257833 A 6 要包括以下步骤： 0050 1)搭建基于拉普拉斯范数的快速迭代收缩阈值的声源识别系统， 主要包括麦克风 阵列、 多通道的信号采集器和数据处理器。 所述麦克风阵列包括分布在声源检测空间内的M 个麦克风。 0051 数据处理器为电脑。 0052。

18、 2)每个麦克风分别监测K个等效声源的时域模拟声压信号， 并发送至多通道的信 号采集器。 所述K个等效声源为随机分布于声源平面周围的等效声源。 0053 3)所述多通道的信号采集器将接收到的时域模拟声压信号转换为数字声压信号 p， 并发送至数据处理器。 0054 4)基于接收到的数字声压信号p， 所述数据处理器建立声源和麦克风阵列之间的 传递矩阵A。 0055 进一步， 建立声源和麦克风阵列之间的传递矩阵A的主要步骤如下： 0056 4.1)确定第m个麦克风测量到的声压信号p(m)如下所示： 0057 0058式中， m1,2,M。为自由场格林函数。 k为波数,为声源 到全息面的距离。 qk为。

19、假想的单极子声源强度。 K为波的总数， 也即空间维度。 j为虚数。 0059 4.2)将公式(1)转化成向量-矩阵形式， 即： 0060 pAq。 (2) 0061 式中， A为MK声场传递矩阵。 p为M维测量声压向量， 其元素为对应传感器在单一 频率f对应的频域响应。 q为K维声源强度列向量， 其元素成分代表声源平面中对应点处的强 度分布。 0062 5)所述数据处理器建立声场声源分布模型。 0063 进一步， 声场声源分布模型的目标函数如下所示： 0064 0065 式中， |2指的是l2范数。 |1指的是l1范数。 为Tikhonov正则化参数。 0066 其中， 正则化参数 如下所示：。

20、 0067 0068 式中， SNR为信噪比。 smax为声场传递矩阵A奇异值SVD分解后的最大值。 0069 声场传递矩阵A奇异值SVD分解如下： 0070 sdiag(S)。 (5) 0071 式中， diag为将向量转化成对角矩阵形式。 S表示声场传递矩阵A奇异值SVD， s为分 解后的声场传递矩阵A奇异值。 0072 6)所述数据处理器利用快速迭代收缩阈值算法对声场声源分布模型进行迭代解 算， 得到声源识别结果x。 说明书 4/7 页 7 CN 111257833 A 7 0073 进一步， 利用快速迭代收缩阈值算法对声场声源分布模型进行迭代解算的主要步 骤如下： 0074 6.1)计。

21、算第n次最速下降y(n)， 即： 0075 0076 式中， L为Lipschitz常数。 L2smax2。 迭代初始值y(1)0； t(1)1。 y(n)为x(n-1)和x (n-2)特定的线性连接方式。 n为迭代次数。 梯度函数 T表示矩 阵的转置。 x(n)和x(n-1)表示第n次和第n-1次迭代解算出的声源识别结果。 0077 其中， x(n)和x(n-1)的线性连接方式y(n+1)满足下式： 0078 0079 其中， 迭代步长t满足下式： 0080 0081 式中， t(n+1)和t(n)表示第n+1次和第n次迭代的步长。 0082 2)设定拉普拉斯范数阈值 。 当第n次迭代的声源。

22、识别结果|y(n)| /(2* L)时， 输出第n次迭代的声源识别结果x(n)， 反之进入步骤3。 0083 3)利用拉普拉斯范数对第n次迭代的声源识别结果x(n)进行更新， 并返回步骤1。 第 n+1次迭代的声源识别结果x(n+1)如下所示： 0084 0085 式中， 为拉普拉斯范数的阶数。 0086 实施例2： 0087 一种基于拉普拉斯范数的快速迭代收缩阈值的声源识别方法， 主要包括以下步 骤： 0088 1)搭建基于拉普拉斯范数的快速迭代收缩阈值的声源识别系统， 主要包括麦克风 阵列、 多通道的信号采集器和数据处理器。 所述麦克风阵列包括分布在声源检测空间内的M 个麦克风。 0089。

23、 2)每个麦克风分别监测K个等效声源的时域模拟声压信号， 并发送至多通道的信 号采集器。 所述K个等效声源为随机分布于声源平面周围的等效声源。 0090 3)所述多通道的信号采集器将接收到的时域模拟声压信号转换为数字声压信号 p， 并发送至数据处理器。 0091 4)基于接收到的数字声压信号p， 所述数据处理器建立声源和麦克风阵列之间的 传递矩阵A。 说明书 5/7 页 8 CN 111257833 A 8 0092 5)所述数据处理器建立声场声源分布模型。 0093 6)所述数据处理器利用快速迭代收缩阈值算法对声场声源分布模型进行迭代解 算， 得到声源识别结果x。 0094 实施例3： 00。

24、95 一种基于拉普拉斯范数的快速迭代收缩阈值的声源识别方法， 主要步骤见实施例 2， 其中， 声场声源分布模型的目标函数如下所示： 0096 0097 式中， |2指的是l2范数。 |1指的是l1范数。 为Tikhonov正则化参数。 0098 其中， 正则化参数 如下所示： 0099 0100 式中， SNR为信噪比。 smax为声场传递矩阵A奇异值SVD分解后的最大值。 0101 声场传递矩阵A奇异值SVD分解如下： 0102 sdiag(S)。 (3) 0103 式中， diag为将向量转化成对角矩阵形式。 0104 实施例4： 0105 一种基于拉普拉斯范数的快速迭代收缩阈值的声源识别。

25、方法， 主要步骤见实施例 2， 其中， 利用快速迭代收缩阈值算法对声场声源分布模型进行迭代解算的主要步骤如下： 0106 6.1)计算第n次最速下降y(n)， 即： 0107 0108 式中， L为Lipschitz常数。 L2smax2。 迭代初始值y(1)0； t(1)1。 y(n)为x(n-1)和x (n-2)的线性连接方式。 n为迭代次数。 梯度函数 T表示矩阵的 转置。 0109 其中， 线性连接方式y(n+1)满足下式： 0110 0111 其中， 迭代步长t满足下式： 0112 0113 2)设定拉普拉斯范数阈值 。 当第n次迭代的声源识别结果|y(n)| /(2* L)时， 输。

26、出第n次迭代的声源识别结果x(n)， 反之进入步骤3。 0114 3)利用拉普拉斯范数对第n次迭代的声源识别结果x(n)进行更新， 并返回步骤1。 第 n+1次迭代的声源识别结果x(n+1)如下所示： 说明书 6/7 页 9 CN 111257833 A 9 0115 0116 式中， 为拉普拉斯范数的阶数。 0117 实施例5： 0118 一种基于拉普拉斯范数的快速迭代收缩阈值的声源识别系统， 主要包括麦克风阵 列、 多通道的信号采集器和数据处理器。 0119 所述麦克风阵列包括声压或者声速传感器、 麦克风支架以及相关线束。 在麦克风 阵列中， 麦克风的分布可以是规则化的或者是非规则化的一、。

27、 二、 三维空间分布， 根据实际 测量需求而选择。 传感器数目的选择和传声器间距的选择则根据需要测量的频率范围而决 定。 0120 所述多通道的信号采集器用于测量时域模拟信号， 并将其转化成数字信号传输给 计算机进行后处理。 其应该包括多通道的前置放大器以及多通道的具有高精度的模拟/数 字转换器， 甚至应该包括相关的滤波器， 从而保证输出到计算机的数字信号足够精确。 0121 所述计算机的主要功能主要是存储以及后处理采集到的声场信号， 并能高效显示 重构声场分布情况。 0122 实施例6： 0123 一种验证基于拉普拉斯范数的快速迭代收缩阈值的声源识别方法的实验， 主要包 括以下步骤： 012。

28、4 1)参见图3至图5， 设置对比方法： 0125 M1： 基于拉普拉斯范数的快速迭代收缩阈值的声源识别方法(即实施例1公开方 法) 0126 M2： 两步迭代收缩阈值算法方法 0127 M3： Tikhonov方法(基于波叠加和正则化的声源识别方法) 0128 2)分别利用M1、 M2、 M3方法对声源进行识别。 0129 参见图6， 当声源频段为0500Hz时， M2方法误差超过100， 而M1方法误差低于 40。 故， 证明M1方法比M2方法在低频效果下识别效果好。 0130 当声源频率为500Hz时， 理论值显示(-0.2,0)和(0.2,0)分别有两个声源， M3方法 和M2方法均无法有效识别出两个声源， 而M1方法则有效的识别到两个声源。 故， 证明M1方法 比M2方法、 M3方法声源识别效果好。 说明书 7/7 页 10 CN 111257833 A 10 图1 图2 说明书附图 1/4 页 11 CN 111257833 A 11 图3 说明书附图 2/4 页 12 CN 111257833 A 12 图4 说明书附图 3/4 页 13 CN 111257833 A 13 图5 图6 说明书附图 4/4 页 14 CN 111257833 A 14 。