CN111475961B - 一种传声器阵列的自适应阵型优化设计方法 - Google Patents

一种传声器阵列的自适应阵型优化设计方法 Download PDF

Info

Publication number
CN111475961B
CN111475961B CN202010316228.7A CN202010316228A CN111475961B CN 111475961 B CN111475961 B CN 111475961B CN 202010316228 A CN202010316228 A CN 202010316228A CN 111475961 B CN111475961 B CN 111475961B
Authority
CN
China
Prior art keywords
array
microphone
annular
weight factor
type
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
CN202010316228.7A
Other languages
English (en)
Other versions
CN111475961A (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.)
Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center
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 Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202010316228.7A priority Critical patent/CN111475961B/zh
Publication of CN111475961A publication Critical patent/CN111475961A/zh
Application granted granted Critical
Publication of CN111475961B publication Critical patent/CN111475961B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
    • G01S5/22—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Circuit For Audible Band Transducer (AREA)

Abstract

本发明公开了一种传声器阵列的自适应阵型优化设计方法，该方法推导了环形传声器阵列的响应解析表达式和分辨率约束函数，构建了由多圈环形阵列组成的传声器阵列性能评价函数，建立了基于环形阵列的阵型优化算法，实现了噪声源的高分辨率、高动态范围识别和定位测量。该方法的优势在于传声器阵列的阵型优化在全频段进行，并以空间分辨率为约束条件；通过对阵列中的传声器进行归类，对每一类传声器设置一个权重因子，大幅减小算法模型中自变量数目，并能根据测量频率范围，智能选取权重因子值大的阵元组成新的阵列进行噪声源测量和定位，提高阵列优化处理速度和传声器利用率。

Description

一种传声器阵列的自适应阵型优化设计方法
技术领域
本发明涉及声学测量和气动声学研究领域，具体涉及到是一种传声器阵列的自适应阵型优化设计方法，可用于航空飞行器、车辆、以及声学风洞试验模型等的噪声源分布测量和定位等。
背景技术
传声器阵列技术使用大量传声器识别空间声源波阵面，用相同相位检测空间声场截面，通过对测得的声信号进行相关分析，抑制不相关噪声影响，得到噪声源分布。传声器阵列克服了单个传声器测量信噪比低，不能获得噪声源分布等缺点，具有强大的噪声源识别和定位能力，是现代声学和气动声学研究领域最常用和最重要的测试技术，并广泛应用于航空、航天、地面交通等行业噪声性能评估与降噪研究。空间分辨率(主瓣宽度)和动态范围是评价传声器阵列最重要的性能参数，它们与测量频率范围、阵列孔径、阵元数量、阵元布置方式等参数密切相关。
现有的环形阵列设计方法、螺旋型阵列设计方法、随机阵列优化设计方法等在阵列设计过程中采用的阵列性能评价函数主要以阵列动态范围为目标，以阵元位置为自变量建立，没有对阵元进行归类，优化计算维度高，计算量大；而且针对不同频率声源，阵元贡献度不同，阵元冗余度未系统考虑。同时在航空飞行器气动噪声源测量过程中，测量频率范围宽，阵列动态范围和空间分辨率要求高，阵列通道数较多，一般超过100通道。
发明内容
本发明的目的是提供一种传声器阵列的自适应阵型优化设计方法，兼顾阵列动态范围和分辨率优化，提高阵列阵元利用率和优化计算速度。
为了实现上述目的，本发明采用如下技术方案：
一种传声器阵列的自适应阵型优化设计方法，包括以下过程：
S1:基于平面波理论和贝塞尔函数性质，推导环形阵列的响应函数解析表达式，建立环形阵列性能评价函数与阵元数量之间的关系式；
S2:将每一圈环形传声器阵列归为一类，设置该类阵列的权重因子，并以环形阵列半径和权重因子为该类阵列表征参数；
S3:将多圈环形传声器阵列进行组合，构建了以环形阵列半径、权重因子为自变量的传声器阵列优化评价目标函数；
S4:根据测量频率范围、阵列孔径、阵元数量、测量视角等约束条件参数，采用多目标优化算法进行优化计算，得到传声器阵列的阵元分布和各圈环形阵列的权重因子值；
S5:根据实际测量频率、动态范围要求、分辨率要求，在传声器阵列中智能选取权重因子值大的阵元组成所需的传声器阵列，并以该传声器阵列进行噪声源测量，提高传声器阵列的阵元利用率，进而提高阵列数据处理速度。
在上述技术方案中，利用平面波理论，单极子点声源位于空间
Figure BDA0002459683610000021
处声压为：
Figure BDA0002459683610000022
其中
Figure BDA0002459683610000023
为波数向量,简称波矢量，
Figure BDA0002459683610000024
为位置矢量。
在上述技术方案中，以传声器阵列中心为参考点，环形阵列响应函数为：
Figure BDA0002459683610000025
其中：θ为相对于阵列中心的视角，φ为相对于阵列中心的阵元方位角，Mn为组成环形阵列的传声器数量，阵列孔径为2Rn，m为阵元的数量。
在上述技术方案中，传声器阵列由N圈环形阵列组成，每圈环形阵列半径为RN，每圈环形阵列的阵元数为MN，那么传声器阵列的响应函数为：
Figure BDA0002459683610000031
其中：M＝M1+M2……+MN，ωn为每圈传声器阵列设置权重因子。
在上述技术方案中，传声器阵列优化设计模型为：
Figure BDA0002459683610000032
Figure BDA0002459683610000033
其中：BW0为输入的分辨率要求，BW为环形阵列的分辨率，H为测量距离，D为阵列孔径，λ为分析频率对应的声波波长。
综上所述，由于采用了上述技术方案，本发明的有益效果是：
传声器阵列的阵型优化在全频段进行，并以空间分辨率为约束条件；通过对阵列中的传声器进行归类，对每一类传声器设置一个权重因子，大幅减小算法模型中自变量数目，并能根据测量频率范围，智能选取权重因子值大的阵元组成新的阵列进行噪声源测量和定位，提高阵列优化处理速度和传声器利用率。
附图说明
本发明将通过例子并参照附图的方式说明，其中：
图1为传声器阵列的自适应阵型优化设计方法流程图；
图2是本发明设计的135通道阵列和国外某机构140通道阵列对比图；
图3是本发明设计的135通道阵列和国外某机构140通道阵列动态范围对比图；
图4是本发明设计的135通道阵列和国外某机构140通道阵列分辨率对比图。
具体实施方式
本说明书中公开的所有特征，或公开的所有方法或过程中的步骤，除了互相排斥的特征和/或步骤以外，均可以以任何方式组合。
本说明书(包括任何附加权利要求、摘要和附图)中公开的任一特征，除非特别叙述，均可被其他等效或具有类似目的的替代特征加以替换。即，除非特别叙述，每个特征只是一系列等效或类似特征中的一个例子而已。
本实施例如图1所示，为了评估阵列的性能，首先需要生成阵列响应图谱，阵列响应图谱由阵列响应函数生成，阵列响应函数可采用波束形成算法得到。
具体如下：
设一单极子点声源位于
Figure BDA0002459683610000041
那么空间r>0处声压为：
Figure BDA0002459683610000042
其中C为常数，r是声音的传播距离，ω是声波的频率，k是波数。那么阵列的响应函数为：
Figure BDA0002459683610000043
式中
Figure BDA0002459683610000044
为声源位置，r0和r0m分别为声源至传声器阵列中心和第m个阵元的距离，wm是第m个阵元的相应权重因子，可用于调整阵列的响度。
在阵列优化设计过程中，为简化问题，通常采用平面波理论,那么空间
Figure BDA0002459683610000045
处声压表示为：
Figure BDA0002459683610000046
其中
Figure BDA0002459683610000047
为波数向量,简称波矢量，
Figure BDA0002459683610000048
为位置矢量。
设有一套Mn个传声器组成的均匀环形阵列，阵列孔径为2Rn，阵元坐标为
Figure BDA0002459683610000049
选取阵列中心原点为参考点，环形阵列响应函数为
Figure BDA0002459683610000051
其中θ为相对于阵列中心的视角，φ为相对于阵列中心的阵元方位角。
当环形阵列阵元数M足够大时,环形阵列的分辨率为
Figure BDA0002459683610000052
其中Const为具体值，H为测量距离，D为阵列孔径，λ为分析频率对应的声波波长。
假设传声器阵列由N圈环形阵列组成，每圈环形阵列对应的半径分别为R1，R2……RN，每圈阵列对应的阵元数分别为M1，M2……MN，阵列总阵元数M＝M1+M2……+MN。每圈传声器阵列设置权重因子为ωn，那么传声器阵列的响应函数为
Figure BDA0002459683610000053
那么传声器阵列优化设计模型为：
Figure BDA0002459683610000054
Figure BDA0002459683610000055
其中BW0输入的分辨率要求。上述模型中自变量只包含环形阵列圈数和权重因子，自变量大幅较小。
阵列设计过程中输入频率范围、视角、测量距离、空间分辨率、动态范围等参数边可开展计算。具体根据试验条件和分辨率要求确定阵列孔径，以动态范围作为优化对象，采用数学的多目标优化设计方法在不同频段进行优化设计得到阵列分布和权重因子。
图2至图3给出了本实施例设计的135通道传声器阵列与国外某机构140通道传声器阵列对比结果。
由图中可知，本实施例设计的135通道的传声器阵列在400Hz到10kHz范围内动态范围大于12dB，优于国外某知名气动声学研究机构140通道传声器阵列；135通道直径只有3米，分辨率却接近国外某机构直径4米的140通道传声器阵列的分辨率水平。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合，以及披露的任一新的方法或过程的步骤或任何新的组合。

Claims (3)

1.一种传声器阵列的自适应阵型优化设计方法，其特征在于：
S1:基于平面波理论和贝塞尔函数性质，推导环形阵列的响应函数解析表达式，建立环形阵列性能评价函数与阵元数量之间的关系式；
S2: 将每一圈环形传声器阵列归为一类，设置该类阵列的权重因子，并以环形阵列半径和权重因子为该类阵列表征参数；
S3: 将多圈环形传声器阵列进行组合，构建以环形阵列半径、权重因子为自变量的传声器阵列优化评价目标函数；
S4: 根据测量频率范围、阵列孔径、阵元数量、测量视角的约束条件参数，采用多目标优化算法进行优化计算，得到传声器阵列的阵元分布和各圈环形阵列的权重因子值；
S5: 根据实际测量频率、动态范围要求、分辨率要求，在传声器阵列中智能选取权重因子值大的阵元，组成所需的传声器阵列，所述传声器的模型为：
Figure 474650DEST_PATH_IMAGE001
Figure 648143DEST_PATH_IMAGE002
，
Figure 330928DEST_PATH_IMAGE003
其中：
Figure 236567DEST_PATH_IMAGE004
为输入的分辨率要求，
Figure 332699DEST_PATH_IMAGE005
为环形阵列的分辨率，
Figure 247566DEST_PATH_IMAGE006
为测量距离，
Figure 283392DEST_PATH_IMAGE007
为阵列孔径，
Figure 359933DEST_PATH_IMAGE008
为分析频率对应的声波波长，
Figure 943361DEST_PATH_IMAGE009
为每圈环形阵列的权重因子，
Figure DEST_PATH_IMAGE010
为每圈环形阵列的阵元数，
Figure 566978DEST_PATH_IMAGE011
，N为环形阵列的圈数，
Figure DEST_PATH_IMAGE012
为相对于阵列中心的视角，
Figure DEST_PATH_IMAGE014
为传声器阵列的响应函数。
2.根据权利要求1所述的一种传声器阵列的自适应阵型优化设计方法，其特征在于：利用平面波理论，单极子点声源位于空间
Figure 427618DEST_PATH_IMAGE015
处声压为：
Figure DEST_PATH_IMAGE016
其中
Figure 675059DEST_PATH_IMAGE017
为波数向量,简称波矢量，
Figure 949046DEST_PATH_IMAGE015
为位置矢量。
3.根据权利要求2所述的一种传声器阵列的自适应阵型优化设计方法，其特征在于：以阵列中心原点为参考点，环形阵列响应函数为：
Figure DEST_PATH_IMAGE018
其中：
Figure 173092DEST_PATH_IMAGE012
为相对于阵列中心的视角，
Figure 216134DEST_PATH_IMAGE019
为相对于阵列中心的阵元方位角，
Figure DEST_PATH_IMAGE020
为组成环形阵列的传声器数量，阵列孔径为
Figure 368898DEST_PATH_IMAGE021
，m为阵元的数量。
CN202010316228.7A 2020-04-21 2020-04-21 一种传声器阵列的自适应阵型优化设计方法 Active CN111475961B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010316228.7A CN111475961B (zh) 2020-04-21 2020-04-21 一种传声器阵列的自适应阵型优化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010316228.7A CN111475961B (zh) 2020-04-21 2020-04-21 一种传声器阵列的自适应阵型优化设计方法

Publications (2)

Publication Number Publication Date
CN111475961A CN111475961A (zh) 2020-07-31
CN111475961B true CN111475961B (zh) 2023-01-10

Family

ID=71756010

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010316228.7A Active CN111475961B (zh) 2020-04-21 2020-04-21 一种传声器阵列的自适应阵型优化设计方法

Country Status (1)

Country Link
CN (1) CN111475961B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116976229A (zh) * 2023-09-22 2023-10-31 中国空气动力研究与发展中心低速空气动力研究所 一种湍流边界层脉动压力二维测量阵列阵型设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106125047A (zh) * 2016-06-14 2016-11-16 北京科技大学 一种气流环境下的声源定位方法
CN106709216A (zh) * 2017-02-24 2017-05-24 中国航天空气动力技术研究院 一种考虑声传播相关性丢失的麦克风阵列优化设计方法
CN110049408A (zh) * 2019-05-10 2019-07-23 苏州静声泰科技有限公司 一种麦克风传声器阵列阵型优化方法
CN110457796A (zh) * 2019-07-29 2019-11-15 南京大学 无规平面差分传声器阵列阵元位置的优化方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7783060B2 (en) * 2005-05-10 2010-08-24 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Deconvolution methods and systems for the mapping of acoustic sources from phased microphone arrays
US9502022B2 (en) * 2010-09-02 2016-11-22 Spatial Digital Systems, Inc. Apparatus and method of generating quiet zone by cancellation-through-injection techniques
CN108562871A (zh) * 2018-04-27 2018-09-21 国网陕西省电力公司电力科学研究院 基于矢量传声器阵列的低频噪声源高精度定位方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106125047A (zh) * 2016-06-14 2016-11-16 北京科技大学 一种气流环境下的声源定位方法
CN106709216A (zh) * 2017-02-24 2017-05-24 中国航天空气动力技术研究院 一种考虑声传播相关性丢失的麦克风阵列优化设计方法
CN110049408A (zh) * 2019-05-10 2019-07-23 苏州静声泰科技有限公司 一种麦克风传声器阵列阵型优化方法
CN110457796A (zh) * 2019-07-29 2019-11-15 南京大学 无规平面差分传声器阵列阵元位置的优化方法

Also Published As

Publication number Publication date
CN111475961A (zh) 2020-07-31

Similar Documents

Publication Publication Date Title
CN109669178B (zh) 一种星载三阵元单脉冲两维测向方法
Sarradj A generic approach to synthesize optimal array microphone arrangements
CN105318963B (zh) 使用自适应孔径相控阵列的噪声源分解系统和方法
CN111475961B (zh) 一种传声器阵列的自适应阵型优化设计方法
CN108447499A (zh) 一种双层圆环麦克风阵列语音增强方法
CN108957406A (zh) 一种基于神经网络的雷达主瓣干扰抑制方法
CN112146751B (zh) 一种基于时域等效源法的实时声场分离方法
US20090326870A1 (en) System and characterization of noise sources using acoustic phased arrays and time series correlations
Sijtsma et al. In-duct and far-field mode detection techniques
EP2201564A1 (en) Enhanced sound source localization system and method by using a movable microphone array
CN109633527B (zh) 基于低秩及几何约束的嵌入式平面麦克风阵列声源测向方法
CN106910493B (zh) 一种多环同轴圆环阵稳健超指向性波束形成方法
Zawodny et al. A comparative study of a 1/4-scale Gulfstream G550 aircraft nose gear model
Breen et al. Limitations of acoustic beamforming for accurate jet noise source location
Guérin et al. Frequency-domain reconstruction of the point-spread function for moving sources
Sarradj Three-dimensional acoustic source mapping
US20220308151A1 (en) Method for determining an angle of arrival, device, computer program product and non-volatile storage medium
Bahr et al. Designing shading schemes for microphone phased arrays
CN111323750B (zh) 一种基于声矢量阵列网络的直接定位方法
Kobayashi et al. Precise sound source model for aircraft noise prediction based on noise source distribution determined by phased array beamforming
Pan et al. Multiple spherical arrays design for acoustic source localization
Brausch Flight velocity influence on jet noise of conical ejector, annular plug and segmented suppressor nozzles
Altena et al. Comparison of acoustic localisation techniques for drone position estimation using real-world experimental data
Battista et al. Tuning of the functional beamforming resolution for wind tunnel measurements
CN115453503B (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