CN114858271B - 一种声音探测用的阵列放大方法 - Google Patents

一种声音探测用的阵列放大方法 Download PDF

Info

Publication number
CN114858271B
CN114858271B CN202210781867.XA CN202210781867A CN114858271B CN 114858271 B CN114858271 B CN 114858271B CN 202210781867 A CN202210781867 A CN 202210781867A CN 114858271 B CN114858271 B CN 114858271B
Authority
CN
China
Prior art keywords
array
sound
sound pressure
pressure value
received
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210781867.XA
Other languages
English (en)
Other versions
CN114858271A (zh
Inventor
曹祖杨
包君康
张鑫
侯佩佩
闫昱甫
黄铖栋
方吉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Crysound Electronics Co Ltd
Original Assignee
Hangzhou Crysound Electronics Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hangzhou Crysound Electronics Co Ltd filed Critical Hangzhou Crysound Electronics Co Ltd
Priority to CN202210781867.XA priority Critical patent/CN114858271B/zh
Publication of CN114858271A publication Critical patent/CN114858271A/zh
Application granted granted Critical
Publication of CN114858271B publication Critical patent/CN114858271B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

一种声音探测用的阵列放大方法，属于信号处理技术领域。方法包括步骤S10，获取对声源测量后的实际阵列接收到的声压值
Figure 724301DEST_PATH_IMAGE001
；步骤S20，基于噪声因素，计算扩展阵列接收到的第一声压值
Figure DEST_PATH_IMAGE002
，以及基于声源因素，计算扩展阵列接收的第二声压值
Figure 438179DEST_PATH_IMAGE003
；步骤S30，进行多次声源测量，重复步骤S10‑步骤S20，获得多个第一声压值和第二声压值；步骤S40，基于公式
Figure 483495DEST_PATH_IMAGE004
估算出最优的扩展阵列接收到的声压值。本发明方法简单，能兼顾噪声和声源双方面因素来扩展获得虚拟阵列，解决了现有阵列扩展方法运算量庞大的问题；还能适用规则或不规则阵列的放大，适用性广。

Description

一种声音探测用的阵列放大方法
技术领域
本发明涉及信号处理技术领域，尤其涉及一种声音探测用的阵列放大方法。
背景技术
在现有技术中，声源定位方法可以通过利用声音传感器阵列来采集声音信号，进行定位。这种利用声音传感器阵列的声学检测手段具有测量范围广、非接触测量、故障预警能力强等又是，其可以实现对电力放电、气体泄漏、车辆噪声中异常声源的主动定位，确定故障区段，防止故障进一步蔓延。
对于给定的阵列配置，物理阵元数和阵列孔径是确定的，测量声音的准确度受阵列面积影响。如何在不改变阵列配置下提升检测性能是目前重点研究方向，为此现有技术提出了阵列扩展技术，通过虚拟增加阵元数目，实现阵列孔径扩展。现有的阵列扩展技术主要有：基于高阶累积量法、基于内插变换法、基于信源和接收信号特性的重构数据法的阵列扩展方法等。现有方法均存在计算量庞大的问题，且适用于均匀阵列，如线性阵列的扩展。
发明专利申请CN202010708433.8公开了一种线列阵阵列扩展方法，并具体公开了方法包括：获取线列阵的接收信号；将线列阵按奇、偶阵元分成两个子阵；计算奇、偶阵元构成的两个子阵的接收信号协方差矩阵；根据接收信号协方差矩阵构造扩展接收阵列信号；利用常规波束形成或最小方差无畸变响应波束形成方法进行波束形成和目标检测。该发明只适用于线性阵列，无法适用不规则阵列。
发明专利申请CN202110308374.X公开了一种基于虚拟阵列扩展的变电站声源定位方法及系统，并具体公开了方法包括步骤：(1)设置声音传感器均匀平面阵列，以采集K个声音信号，所述声音传感器均匀平面阵列具有M×M个阵元；(2)采用四阶累积量阵列扩展方法对声音传感器均匀平面阵列进行扩展；(3)基于扩展阵列构建空间谱函数；(4)对空间谱函数进行极大值搜索，并基于该极大值获得与其对应的入射方位角和入射俯仰角θ的估计值，以对变电站设备声源进行定位。该发明只能扩展均匀阵列。
发明内容
本发明针对现有技术存在的问题，提出了一种简单且适用性广的声音探测用的阵列放大方法，不仅能兼顾噪声和声源双方面因素来扩展获得虚拟阵列，解决了现有阵列扩展方法运算量庞大的问题，还能适用规则或不规则阵列。
本发明提供一种声音探测用的阵列放大方法，包括：
步骤S10，获取对声源测量后的实际阵列接收到的声压值
Figure 100002_DEST_PATH_IMAGE001
；
Figure 897157DEST_PATH_IMAGE002
其中，S p （t）为P个声源发出的信号，n N （t）为N个阵元接收到的噪声，
Figure 100002_DEST_PATH_IMAGE003
，d表示阵元之间的间距，c表示声速，
Figure 969018DEST_PATH_IMAGE004
表示声源的方位角，
Figure 100002_DEST_PATH_IMAGE005
表示声源的俯仰角；
步骤S20，利用公式（1）获得扩展阵列接收到的第一声压值
Figure 100002_DEST_PATH_IMAGE007
，以及利用公式（2）获得扩展阵列接收的第二声压值
Figure 409227DEST_PATH_IMAGE008
；
Figure 100002_DEST_PATH_IMAGE009
公式（1）
其中，P为转换矩阵，
Figure 541131DEST_PATH_IMAGE010
，
Figure 100002_DEST_PATH_IMAGE011
，
Figure 801211DEST_PATH_IMAGE012
，
Figure 100002_DEST_PATH_IMAGE013
为AN的转置矩阵，AN为实际阵列中P个声源所在空间位置参量矩阵；
Figure 676763DEST_PATH_IMAGE014
，AM为扩展阵列中M个声源所在空间位置参量矩阵；
Figure 100002_DEST_PATH_IMAGE015
公式（2）
其中，
Figure 705899DEST_PATH_IMAGE016
，
Figure 100002_DEST_PATH_IMAGE017
为向量系数，
Figure 274283DEST_PATH_IMAGE018
表示空气密度，
Figure 490501DEST_PATH_IMAGE020
通过
Figure 100002_DEST_PATH_IMAGE021
计算获得，
Figure 169744DEST_PATH_IMAGE022
表示声源的方位角，
Figure 100002_DEST_PATH_IMAGE023
表示声源的俯仰角，r表示声源到阵列中心的距离；
Figure 318966DEST_PATH_IMAGE024
表示第b阶第一类球Hankel函数，
Figure 100002_DEST_PATH_IMAGE025
表示伴随拉格朗日函数，k表示波数，
Figure 100002_DEST_PATH_IMAGE027
表示相数；
步骤S30，进行多次声源测量，重复步骤S10-步骤S20，获得多个第一声压值
Figure 58252DEST_PATH_IMAGE028
和第二声压值
Figure 100002_DEST_PATH_IMAGE029
；
步骤S40，基于公式（3）估算出最优的扩展阵列接收到的声压值
Figure 292924DEST_PATH_IMAGE030
；
Figure 100002_DEST_PATH_IMAGE031
公式（3）
其中，
Figure 244699DEST_PATH_IMAGE032
为步骤S30获得的多个第一声压值
Figure 100002_DEST_PATH_IMAGE033
的平均值，
Figure 100002_DEST_PATH_IMAGE035
为步骤S30获得的多个第二声压值
Figure 779586DEST_PATH_IMAGE036
的平均值，
Figure 100002_DEST_PATH_IMAGE037
，
Figure 695632DEST_PATH_IMAGE038
为通过步骤S30获得的多个第一声压值计算得到的方差，
Figure 100002_DEST_PATH_IMAGE039
为通过步骤S30获得的多个第二声压值计算得到的方差。
本发明利用两种方式计算扩展阵列的声压值，如基于噪声的转换关系与声源的转换关系一致的原则获得转换矩阵，继而确定第一声压值，又模拟声源的构成来计算第二声压值。这样，本发明在考虑噪声和声源双层因素的前提下，最后进行最优估计获得融合后的新的扩展阵列。本发明方法简单，能适用多种形状的阵列；本发明为声音探测提供了一个准确度高、计算量无需庞大复杂的阵列扩展方法。
作为优选，所述转换矩阵在计算第一声压值
Figure 152021DEST_PATH_IMAGE040
前预先计算获得。
作为优选，所述转换矩阵预先计算获得过程如下：
基于实际阵列形状构建虚拟阵列；
利用波束形成算法获得实际阵列中实际阵元的空间位置参量，并根据实际阵列和扩展阵列的物理位置，确定虚拟阵列中虚拟阵元的空间位置参量；
根据
Figure 100002_DEST_PATH_IMAGE041
计算获得转换矩阵。
作为优选，转换矩阵计算公式确定过程如下：
定义N个阵元的实际阵列接收到的声压值为：
Figure 173067DEST_PATH_IMAGE042
；
定义扩展M个阵元的扩展阵列接收到的声压值为：
Figure 100002_DEST_PATH_IMAGE043
；
设定扩展阵列与实际阵列的关系为：
Figure 765722DEST_PATH_IMAGE044
；
将上述三个公式转换为
Figure 100002_DEST_PATH_IMAGE045
；
假定噪声的转换关系与声源的转换关系一致，则获得转换矩阵
Figure 846811DEST_PATH_IMAGE046
。
作为优选，所述步骤S20中第二声压值
Figure DEST_PATH_IMAGE047
的计算过程具体为：
计算向量系数
Figure 790496DEST_PATH_IMAGE048
；
计算第二声压值
Figure 67763DEST_PATH_IMAGE049
，并将扩展阵列的实际阵元的位置接收到的数据，用实际阵列进行代替；
重复上述过程，直至计算获得的第二声压值
Figure 514924DEST_PATH_IMAGE050
趋于稳定。
作为优选，在计算第二声压值
Figure DEST_PATH_IMAGE051
的过程中，第一次计算的向量系数
Figure DEST_PATH_IMAGE053
是基于实际阵列接收到的声压值
Figure DEST_PATH_IMAGE055
计算获得，
Figure 32493DEST_PATH_IMAGE056
；随后重复过程中计算的向量系数
Figure DEST_PATH_IMAGE057
时候基于前次计算获得的第二声压值
Figure 729054DEST_PATH_IMAGE058
计算获得，
Figure 560744DEST_PATH_IMAGE059
。
作为优选，方法还包括：
步骤S50，将步骤S40获得的最优的扩展阵列接收到的信号送入声音检测设备中。
作为优选，所述步骤S50具体为：将步骤S40获得的最优的扩展阵列接收到的信号作为波束形成的初始数据，送入声音检测设备中，并得到优化后的成像云图。
作为优选，该方法用于手持式声像仪。
作为优选，所述步骤S20中第一声压值在第二声压值计算前计算获得；或者，所述步骤S20中第一声压值在第二声压值计算后计算获得；或者，所述步骤S20中第一声压值和第二声压值同时进行计算。
本发明具有以下有益效果：
本发明一种声音探测用的阵列放大方法，能适用各种规则的阵列放大，能在不改变阵列实际面积的情况下，扩大阵列的有效计算面积来提高声音的准确度；该方法主要用于声音探测领域，例如手持式声像仪，采用此方法的手持式声像仪不仅能满足探测过程中体积小方便携带的需求外，还能高准确度探测声音。
附图说明
图1为本发明一种声音探测用的阵列放大方法的流程图。
具体实施方式
以下是本发明的具体实施例并结合附图，对本发明的技术方案作进一步的描述，但本发明并不限于这些实施例。
如图1，本发明一种声音探测用的阵列放大方法，包括：
步骤S10，获取对声源测量后的实际阵列接收到的声压值
Figure DEST_PATH_IMAGE060
；
Figure 659150DEST_PATH_IMAGE061
其中，S p （t）为P个声源发出的信号，n N （t）为N个阵元接收到的噪声，
Figure 285303DEST_PATH_IMAGE063
，d表示阵元之间的间距，c表示声速，
Figure 938001DEST_PATH_IMAGE004
表示声源的方位角，
Figure DEST_PATH_IMAGE064
表示声源的俯仰角；
步骤S20，利用公式（1）获得扩展阵列接收到的第一声压值
Figure 104541DEST_PATH_IMAGE033
，以及利用公式（2）获得扩展阵列接收的第二声压值
Figure 260715DEST_PATH_IMAGE029
；
Figure 323349DEST_PATH_IMAGE065
公式（1）
其中，P为转换矩阵，
Figure DEST_PATH_IMAGE066
，
Figure 728923DEST_PATH_IMAGE067
，
Figure DEST_PATH_IMAGE068
，
Figure 433574DEST_PATH_IMAGE013
为AN的转置矩阵，AN为实际阵列中P个声源所在空间位置参量矩阵；
Figure 975414DEST_PATH_IMAGE069
，AM为扩展阵列中M个声源所在空间位置参量矩阵；
Figure DEST_PATH_IMAGE070
公式（2）
其中，
Figure 474528DEST_PATH_IMAGE016
，
Figure 836239DEST_PATH_IMAGE071
为向量系数，
Figure DEST_PATH_IMAGE072
表示空气密度，
Figure 344581DEST_PATH_IMAGE073
通过
Figure DEST_PATH_IMAGE074
计算获得，
Figure 735068DEST_PATH_IMAGE022
表示声源的方位角，
Figure DEST_PATH_IMAGE076
表示声源的俯仰角，r表示声源到阵列中心的距离；
Figure 670663DEST_PATH_IMAGE077
表示第b阶第一类球Hankel函数，
Figure DEST_PATH_IMAGE078
表示伴随拉格朗日函数，k表示波数，
Figure 785249DEST_PATH_IMAGE079
表示相数；
步骤S30，进行多次声源测量，重复步骤S10-步骤S20，获得多个第一声压值
Figure DEST_PATH_IMAGE080
和第二声压值
Figure 831703DEST_PATH_IMAGE081
；
步骤S40，基于公式（3）估算出最优的扩展阵列接收到的声压值
Figure 82555DEST_PATH_IMAGE082
；
Figure DEST_PATH_IMAGE083
公式（3）
其中，
Figure 189052DEST_PATH_IMAGE084
为步骤S30获得的多个第一声压值
Figure 525355DEST_PATH_IMAGE085
的平均值，
Figure 844341DEST_PATH_IMAGE087
为步骤S30获得的多个第二声压值
Figure 949700DEST_PATH_IMAGE029
的平均值，
Figure DEST_PATH_IMAGE088
，
Figure 227098DEST_PATH_IMAGE089
为通过步骤S30获得的多个第一声压值计算得到的方差，
Figure DEST_PATH_IMAGE090
为通过步骤S30获得的多个第二声压值计算得到的方差。
在步骤S10中，利用声音检测设备对声源进行探测，探测后获得N个阵元接收到的p个声源的信号（包括空间位置参数，如声源的方位角和俯仰角）。所述声音检测设备可以是现有任意检测设备，如手持式声像仪。根据获取的参数，公式基于
Figure 50697DEST_PATH_IMAGE091
确定实际阵列接收到的声压值
Figure 438953DEST_PATH_IMAGE093
。之后进行步骤S20的扩展阵列的声压值计算。
在计算第一声压值前，预先计算转换矩阵。所述转换矩阵预先计算获得过程如下：基于实际阵列形状构建虚拟阵列；利用波束形成算法获得实际阵列中实际阵元的空间位置参量，并根据实际阵列和扩展阵列的物理位置，确定虚拟阵列中虚拟阵元的空间位置参量；根据
Figure DEST_PATH_IMAGE094
计算获得转换矩阵。例如，有多个实际阵元排布成圆圈状，由里到外分布成层层圈状；基于实际阵列形状，在最外层按实际阵元排布方式扩展虚拟阵元，形成直径更大的圆圈，可根据需要扩展多圈虚拟阵元，以确保整体阵列的均匀性。之后利用现有的波束形成算法获得实际阵列中实际阵元的空间位置参量，并根据实际阵列和扩展阵列的物理位置，确定虚拟阵列中虚拟阵元的空间位置参量。最后利用转换矩阵公式来计算获得转换矩阵。所述转换矩阵根据扩展的数量、位置不同，根据扩展的实际阵列形状不同而发生变化。为此，需要在计算第一声压值前进行预计算。
上述转换矩阵计算公式也是基于一定考虑因素确定的。具体确定过程如下：
定义N个阵元的实际阵列接收到的声压值为：
Figure 664398DEST_PATH_IMAGE095
；
定义扩展M个阵元的扩展阵列接收到的声压值为：
Figure 581539DEST_PATH_IMAGE096
；
设定扩展阵列与实际阵列的关系为：
Figure 158014DEST_PATH_IMAGE097
；
将上述三个公式转换为
Figure 818802DEST_PATH_IMAGE098
；
当假定噪声的转换关系与声源的转换关系一致时，则获得转换矩阵
Figure 633174DEST_PATH_IMAGE099
。
上述假定主要考虑了噪声因素，通过以此方式获得的转换矩阵来计算扩展阵列的第一声压值。
在计算第二声压值时，主要模拟声源的构成，为此通过公式（2）来计算第二声压值。由于不清楚阵列的形状，需要获得稳定而准确的第二声压值。所述第二声压值的计算过程具体为：首先，计算向量系数
Figure 986795DEST_PATH_IMAGE101
；其次，计算第二声压值
Figure 519408DEST_PATH_IMAGE049
，并将扩展阵列的实际阵元的位置接收到的数据，用实际阵列进行代替；重复上述过程，直至计算获得的第二声压值
Figure 249466DEST_PATH_IMAGE049
趋于稳定。
在计算第二声压值
Figure 918345DEST_PATH_IMAGE102
的过程中，第一次计算的向量系数
Figure 442867DEST_PATH_IMAGE103
是基于实际阵列接收到的声压值
Figure 462776DEST_PATH_IMAGE104
计算获得，
Figure 730946DEST_PATH_IMAGE105
；随后重复过程中计算的向量系数
Figure 519911DEST_PATH_IMAGE106
时候基于前次计算获得的第二声压值计算获得，
Figure 215334DEST_PATH_IMAGE107
。
为了减少误差，当实际阵列和扩展阵列的位置存在重叠部分时，假定认为实际阵列收集到的声压是最真实的，为此，将计算出来的扩展阵列中实际阵元位置的声压值替换为实际阵列接收到的声音值。
所述步骤S20中对于第一声压值和第二声压值的计算可不分先后进行，也可在第二声压值计算前先计算第一声压值，也可在第一声压值计算前先计算第二声压值。
在所述步骤S30后，且在所述步骤S40 前，本发明方法还包括计算多个第一声压值的平均值、方差，以及多个第二声压值的平均值、方差。之后，步骤S40计算公式(3)中的最优g值。
计算最优的g值如下表示：
Figure 456960DEST_PATH_IMAGE108
；
则最优的g表示为
Figure 528821DEST_PATH_IMAGE109
。
之后，基于公式（3）估算出最优的扩展阵列接收到的声压值
Figure 437871DEST_PATH_IMAGE110
；
Figure 38617DEST_PATH_IMAGE111
公式（3）；
Figure 298697DEST_PATH_IMAGE112
为步骤S30获得的多个第一声压值
Figure 643090DEST_PATH_IMAGE033
的平均值，
Figure 406647DEST_PATH_IMAGE034
为步骤S30获得的多个第二声压值
Figure 178294DEST_PATH_IMAGE113
的平均值，
Figure 660091DEST_PATH_IMAGE089
为通过步骤S30获得的多个第一声压值计算得到的方差，
Figure 73755DEST_PATH_IMAGE114
为通过步骤S30获得的多个第二声压值计算得到的方差。
本发明方法还包括：步骤S50，将步骤S40获得的最优的扩展阵列接收到的信号送入声音检测设备中。具体地，将步骤S40获得的最优的扩展阵列接收到的信号作为波束形成的初始数据，送入声音检测设备中，并得到优化后的成像云图。
本领域的技术人员应理解，上述描述及附图中所示的本发明的实施例只作为举例而并不限制本发明。本发明的目的已经完整有效地实现。本发明的功能及结构原理已在实施例中展示和说明，在没有背离所述原理下，本发明的实施方式可以有任何变形或修改。

Claims (10)

1.一种声音探测用的阵列放大方法，其特征在于，包括：
步骤S10，获取对声源测量后的实际阵列接收到的声压值
Figure DEST_PATH_IMAGE001
；
Figure 163384DEST_PATH_IMAGE002
其中，S p （t）为P个声源发出的信号，n N （t）为N个阵元接收到的噪声，
Figure DEST_PATH_IMAGE003
，d表示阵元之间的间距，c表示声速，
Figure 817219DEST_PATH_IMAGE004
表示声源的方位角，
Figure DEST_PATH_IMAGE005
表示声源的俯仰角；
步骤S20，利用公式（1）获得扩展阵列接收到的第一声压值
Figure DEST_PATH_IMAGE007
，以及利用公式（2）获得扩展阵列接收的第二声压值
Figure 573822DEST_PATH_IMAGE008
；
Figure DEST_PATH_IMAGE009
公式（1）
其中，P为转换矩阵，
Figure 756542DEST_PATH_IMAGE010
，
Figure DEST_PATH_IMAGE011
，
Figure 333017DEST_PATH_IMAGE012
，
Figure DEST_PATH_IMAGE013
为A N 的转置矩阵，A N 为实际阵列中P个声源所在空间位置参量矩阵；
Figure 259384DEST_PATH_IMAGE014
，A M 为扩展阵列中M个声源所在空间位置参量矩阵；
Figure DEST_PATH_IMAGE015
公式（2）
其中，
Figure 843730DEST_PATH_IMAGE016
，
Figure DEST_PATH_IMAGE017
为向量系数，
Figure 462931DEST_PATH_IMAGE018
表示空气密度，
Figure DEST_PATH_IMAGE019
通过
Figure 261122DEST_PATH_IMAGE020
计算获得，
Figure DEST_PATH_IMAGE021
表示声源的方位角，
Figure 991181DEST_PATH_IMAGE022
表示声源的俯仰角，r表示声源到阵列中心的距离；
Figure DEST_PATH_IMAGE023
表示第b阶第一类球Hankel函数，
Figure 925639DEST_PATH_IMAGE024
表示伴随拉格朗日函数，k表示波数，
Figure DEST_PATH_IMAGE025
表示相数；
步骤S30，进行多次声源测量，重复步骤S10-步骤S20，获得多个第一声压值
Figure 715740DEST_PATH_IMAGE007
和第二声压值
Figure 1228DEST_PATH_IMAGE008
；
步骤S40，基于公式（3）估算出最优的扩展阵列接收到的声压值
Figure 3819DEST_PATH_IMAGE026
；
Figure DEST_PATH_IMAGE027
公式（3）
其中，
Figure 58363DEST_PATH_IMAGE028
为步骤S30获得的多个第一声压值
Figure 255251DEST_PATH_IMAGE007
的平均值，
Figure DEST_PATH_IMAGE029
为步骤S30获得的多个第二声压值
Figure 762456DEST_PATH_IMAGE008
的平均值，
Figure 568738DEST_PATH_IMAGE030
，
Figure DEST_PATH_IMAGE031
为通过步骤S30获得的多个第一声压值计算得到的方差，
Figure 743367DEST_PATH_IMAGE032
为通过步骤S30获得的多个第二声压值计算得到的方差。
2.根据权利要求1所述的一种声音探测用的阵列放大方法，其特征在于，所述转换矩阵在计算第一声压值
Figure 344113DEST_PATH_IMAGE007
前预先计算获得。
3.根据权利要求2所述的一种声音探测用的阵列放大方法，其特征在于，所述转换矩阵预先计算获得过程如下：
基于实际阵列形状构建虚拟阵列；
利用波束形成算法获得实际阵列中实际阵元的空间位置参量，并根据实际阵列和扩展阵列的物理位置，确定虚拟阵列中虚拟阵元的空间位置参量；
根据
Figure DEST_PATH_IMAGE033
计算获得转换矩阵。
4. 根据权利要求3所述的一种声音探测用的阵列放大方法，其特征在于， 转换矩阵计算公式确定过程如下：
定义N个阵元的实际阵列接收到的声压值为：
Figure 869772DEST_PATH_IMAGE034
；
定义扩展M个阵元的扩展阵列接收到的声压值为：
Figure DEST_PATH_IMAGE035
；
设定扩展阵列与实际阵列的关系为：
Figure 479745DEST_PATH_IMAGE036
；
将上述N个阵元的实际阵列接收到的声压值公式、M个阵元的扩展阵列接收的声压值公式以及扩展阵列与实际阵列的关系公式转换为
Figure DEST_PATH_IMAGE037
；
假定噪声的转换关系与声源的转换关系一致，则获得转换矩阵
Figure 508881DEST_PATH_IMAGE038
。
5.根据权利要求1所述的一种声音探测用的阵列放大方法，其特征在于，所述步骤S20中第二声压值
Figure DEST_PATH_IMAGE039
的计算过程具体为：
计算向量系数
Figure 546107DEST_PATH_IMAGE040
；
计算第二声压值
Figure 526439DEST_PATH_IMAGE039
，并将扩展阵列的实际阵元的位置接收到的数据，用实际阵列进行代替；
重复上述过程，直至计算获得的第二声压值
Figure DEST_PATH_IMAGE041
趋于稳定。
6. 根据权利要求5所述的一种声音探测用的阵列放大方法，其特征在于， 在计算第二声压值
Figure 940103DEST_PATH_IMAGE041
的过程中，第一次计算的向量系数
Figure 558166DEST_PATH_IMAGE042
是基于实际阵列接收到的声压值
Figure DEST_PATH_IMAGE043
计算获得，
Figure 31873DEST_PATH_IMAGE044
；随后重复过程中计算的向量系数
Figure 735387DEST_PATH_IMAGE042
时候基于前次计算获得的第二声压值
Figure 421583DEST_PATH_IMAGE041
计算获得，
Figure DEST_PATH_IMAGE045
。
7.根据权利要求1所述的一种声音探测用的阵列放大方法，其特征在于，还包括：
步骤S50，将步骤S40获得的最优的扩展阵列接收到的信号送入声音检测设备中。
8.根据权利要求7所述的一种声音探测用的阵列放大方法，其特征在于，所述步骤S50具体为：将步骤S40获得的最优的扩展阵列接收到的信号作为波束形成的初始数据，送入声音检测设备中，并得到优化后的成像云图。
9.根据权利要求1所述的一种声音探测用的阵列放大方法，其特征在于，该方法用于手持式声像仪。
10.根据权利要求1所述的一种声音探测用的阵列放大方法，其特征在于，所述步骤S20中第一声压值在第二声压值计算前计算获得；或者，所述步骤S20中第一声压值在第二声压值计算后计算获得；或者，所述步骤S20中第一声压值和第二声压值同时进行计算。
CN202210781867.XA 2022-07-05 2022-07-05 一种声音探测用的阵列放大方法 Active CN114858271B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210781867.XA CN114858271B (zh) 2022-07-05 2022-07-05 一种声音探测用的阵列放大方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210781867.XA CN114858271B (zh) 2022-07-05 2022-07-05 一种声音探测用的阵列放大方法

Publications (2)

Publication Number Publication Date
CN114858271A CN114858271A (zh) 2022-08-05
CN114858271B true CN114858271B (zh) 2022-09-23

Family

ID=82626224

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210781867.XA Active CN114858271B (zh) 2022-07-05 2022-07-05 一种声音探测用的阵列放大方法

Country Status (1)

Country Link
CN (1) CN114858271B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6914854B1 (en) * 2002-10-29 2005-07-05 The United States Of America As Represented By The Secretary Of The Army Method for detecting extended range motion and counting moving objects using an acoustics microphone array
CN111243615A (zh) * 2020-01-08 2020-06-05 环鸿电子（昆山）有限公司 麦克风阵列信号处理方法及手持式装置
CN111736118A (zh) * 2020-07-22 2020-10-02 东南大学 一种线列阵阵列扩展方法
KR102197812B1 (ko) * 2019-09-30 2021-01-04 재단법인 다차원 스마트 아이티 융합시스템 연구단 위치 지정을 통한 소음 모니터링 시스템 및 그 방법
CN112739996A (zh) * 2018-07-24 2021-04-30 弗兰克公司 用于分析和显示声学数据的系统和方法
CN113866718A (zh) * 2021-10-09 2021-12-31 西北工业大学 一种基于互质阵的匹配场被动定位方法
CN113884986A (zh) * 2021-12-03 2022-01-04 杭州兆华电子有限公司 波束聚焦增强的强冲击信号空时域联合检测方法及系统
CN114325584A (zh) * 2022-03-14 2022-04-12 杭州兆华电子股份有限公司 基于合成孔径的多阵元超声波声源三维成像方法及系统
CN114527427A (zh) * 2022-01-27 2022-05-24 华南理工大学 一种基于球形麦克风阵列的低频波束形成声源定位方法
CN114624689A (zh) * 2022-05-12 2022-06-14 杭州兆华电子股份有限公司 一种基于声像仪的近场聚焦声源距离计算方法及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8861756B2 (en) * 2010-09-24 2014-10-14 LI Creative Technologies, Inc. Microphone array system
KR101282673B1 (ko) * 2011-12-09 2013-07-05 현대자동차주식회사 음원 위치 추정 방법

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6914854B1 (en) * 2002-10-29 2005-07-05 The United States Of America As Represented By The Secretary Of The Army Method for detecting extended range motion and counting moving objects using an acoustics microphone array
CN112739996A (zh) * 2018-07-24 2021-04-30 弗兰克公司 用于分析和显示声学数据的系统和方法
KR102197812B1 (ko) * 2019-09-30 2021-01-04 재단법인 다차원 스마트 아이티 융합시스템 연구단 위치 지정을 통한 소음 모니터링 시스템 및 그 방법
CN111243615A (zh) * 2020-01-08 2020-06-05 环鸿电子（昆山）有限公司 麦克风阵列信号处理方法及手持式装置
CN111736118A (zh) * 2020-07-22 2020-10-02 东南大学 一种线列阵阵列扩展方法
CN113866718A (zh) * 2021-10-09 2021-12-31 西北工业大学 一种基于互质阵的匹配场被动定位方法
CN113884986A (zh) * 2021-12-03 2022-01-04 杭州兆华电子有限公司 波束聚焦增强的强冲击信号空时域联合检测方法及系统
CN114527427A (zh) * 2022-01-27 2022-05-24 华南理工大学 一种基于球形麦克风阵列的低频波束形成声源定位方法
CN114325584A (zh) * 2022-03-14 2022-04-12 杭州兆华电子股份有限公司 基于合成孔径的多阵元超声波声源三维成像方法及系统
CN114624689A (zh) * 2022-05-12 2022-06-14 杭州兆华电子股份有限公司 一种基于声像仪的近场聚焦声源距离计算方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Plane-wave expansion for arrays of arbitrarily oriented piecewise linear elements and its application in determining the impedance of a single linear antenna in a lossy half-space;B. Munk;《IEEE Transactions on Antennas and Propagation》;19790531;第27卷(第3期);331-343 *
基于声阵列技术的低频噪声源识别定位研究;朱亚辉;《硕士电子期刊出版工程科技Ⅱ辑》;20210815(第8期);1-82 *
汽车通过噪声声源识别方法研究进展;张鑫等;《汽车实用技术》;20220430;第47卷(第8期);205-209 *

Also Published As

Publication number Publication date
CN114858271A (zh) 2022-08-05

Similar Documents

Publication Publication Date Title
JP5496192B2 (ja) 音響場を再構成するための方法
CN107247251B (zh) 基于压缩感知的三维声源定位方法
CN108375763B (zh) 一种应用于多声源环境的分频定位方法
CN101995574B (zh) 一种近场聚焦波束形成定位法
Kook et al. An efficient procedure for visualizing the sound field radiated by vehicles during standardized passby tests
CN109683134B (zh) 一种面向旋转声源的高分辨率定位方法
CN113868583B (zh) 一种子阵波束聚焦的声源距离计算方法及系统
Yang et al. Functional delay and sum beamforming for three-dimensional acoustic source identification with solid spherical arrays
CN111798869B (zh) 一种基于双麦克风阵列的声源定位方法
CN102147458A (zh) 一种针对宽带声源的波达方向估计方法及其装置
CN110927669A (zh) 一种用于无线声传感器网络的cs多声源定位方法及系统
CN111812581B (zh) 基于原子范数的球面阵列声源波达方向估计方法
CN106997037A (zh) 声矢量传感器阵列空间旋转解相干到达角估计方法
CN108981903A (zh) 一种变压器振动噪声源全频域识别方法
CN109343003A (zh) 一种快速迭代收缩波束形成声源识别方法
CN114858271B (zh) 一种声音探测用的阵列放大方法
CN109375197B (zh) 一种小尺寸矢量阵低频散射校正方法
CN110907892A (zh) 一种球麦克风阵列语音信号到达角估计方法
CN106124044A (zh) 实心球声源识别低旁瓣超高分辨率声学图像快速获取方法
CN113075645B (zh) 一种基于主成分分析-密度聚类的畸变阵形线谱增强方法
Li et al. Wind parameters measurement method based on co-prime array signal processing
CN115061089A (zh) 一种声源定位方法、系统、介质、设备及装置
CN114609580A (zh) 一种基于非圆信号的无孔互质阵列设计方法
CN112558009A (zh) 基于相位模态理论的高阶声场传感器阵列的方位估计方法
CN113296056A (zh) 一种声阵列构型及声源定位方法、系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant