CN102183298B - 不规则单全息声压测量面分离非自由声场的方法 - Google Patents
不规则单全息声压测量面分离非自由声场的方法 Download PDFInfo
- Publication number
- CN102183298B CN102183298B CN 201110049761 CN201110049761A CN102183298B CN 102183298 B CN102183298 B CN 102183298B CN 201110049761 CN201110049761 CN 201110049761 CN 201110049761 A CN201110049761 A CN 201110049761A CN 102183298 B CN102183298 B CN 102183298B
- Authority
- CN
- China
- Prior art keywords
- sound
- point
- holographic
- sound source
- item number
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000001788 irregular Effects 0.000 title claims abstract description 14
- 238000009530 blood pressure measurement Methods 0.000 title claims abstract description 8
- 238000005259 measurement Methods 0.000 claims abstract description 40
- 239000013598 vector Substances 0.000 claims abstract description 17
- 238000000926 separation method Methods 0.000 claims description 18
- 230000005855 radiation Effects 0.000 claims description 9
- 238000013178 mathematical model Methods 0.000 abstract 2
- 239000002131 composite material Substances 0.000 abstract 1
- 206010038743 Restlessness Diseases 0.000 description 6
- 239000002245 particle Substances 0.000 description 6
- 230000010349 pulsation Effects 0.000 description 5
- 239000011796 hollow space material Substances 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000005428 wave function Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001093 holography Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
不规则单全息声压测量面分离非自由声场的方法,按如下步骤进行:a、在由目标声源和背景噪声源构成的非自由声场中,布置不规则全息测量面;b、在不规则全息测量面上放置传声器,测量全息面上的声压值;c、选取不规则全息测量面部分测点的声压值重构其余测点的声压值,确定由来波和去波构成的非自由声场复声压的数学模型,和数学模型的系数矢量;d、重构任意指定空间坐标处的声场;分离重构的由目标声源和背景噪声源组成的声场。
Description
技术领域
本发明涉及一种不规则单全息声压测量面分离非自由声场的方法。
背景技术
一般情况下,测试一个声源的声辐射强度,或评估某个噪声源对周围环境的干扰和影响程度,需要把该声源结构移至具有自由声场的消声室内进行测量。但很多情况下这种做法是不现实的,无法实现。如流水线上产品噪声的实时监测和评估。为了能够在声源工作的状态下,对其有效测量和评估,声场分离方法是最终的解决方案。
1956年,Pachner首次提出使用两个不同半径的同心球形全息测量面,分离声场中的行波和驻波。1980年,Weinreich和Arnold同样使用了两个半径分别为58.1cm和72.6cm的同心球形全息测量面,分离和重构在小全息面内的声源和在大全息面外的声源。2007年,国内学者李卫兵等用同Weinreich和Arnold相似的研究方法、分离由刚性散射球源和点声源组成的声场,得到一些仿真计算结果。2008年C.X.Bi、X.Z.Chen和J.Chen使用系列等效声源法表达声场的声压,在近场布置两个间隔7.5cm的平行全息平面,分离和重构位于全息面两侧两个喇叭辐射的声场。
以上提到的各位学者研究的声场分离方法,都是使用两个间隙很近的、共形的全息测量面。如果声源的结构非常大,精确布置两层平行、共形、规则几何形状的传感器阵列,难度大、成本高、时间长。2007年,FinnJacobsen和Virginie Jaud提出了只用一个全息测量面的声场分离方法。但是,为了识别和分离全息面两侧的声源,在全息面的每个测量点上必须同时测量声压和空气粒子速度。空气粒子速度是矢量,可以用来指出声源在全息面的哪一侧。使用平面波函数作为基函数表达声场的声压,即基于统计最优的近场声全息方法,并在实验室对一个点声源和一个直升机的齿轮箱形成的声场进行了声场分离实验。
Finn Jacobsen等提出的声场分离方法虽然只使用一个全息测量面,但需要在所有测点上逐点扫描,测量其声压和空气粒子速度。空气粒子速度是矢量,测量难度比测量声压大,耗时且成本高。由于声压和空气粒子速度不是同时同步测量获得的,因此,该方法只适用于稳态声场分离情况,不适用于瞬态声场分离情况。
发明内容
本发明所解决的现有声场分离方法存在测量难度大、耗时且成本高、不适用于瞬态声场分离情况的问题,提出一种测量方便、成本低、适用于瞬态声场分离情况的不规则单全息声压测量面分离非自由声场的方法。
本发明解决技术问题所采用的技术方案,
不规则单全息声压测量面分离非自由声场的方法,其特征在于:按如下步骤进行:
a、测量全息测量面上的声压值,在由目标声源和背景噪声源构成的声场中布置全息测量面,在全息测量面节点上放置传声器;
b、选取全息测量面部分测点的声压值重构其余测点的声压值,确定最优展开项数;
把波动方程在球坐标下的完全解为:
在全息测量面上选出部分测点,选取方法是每间隔一点选择一点,共选择N点,把选出点对应的声压按球坐标下声波方程行波解从低阶到高阶分别展开得式(2):
简记为:
其中:p1 p2...pN代表选取点上测得的声压值;代表当目标声源产生去波时,第N点坐标值按球坐标下声波方程行波解展开M项的表达;代表当背景噪声源产生来波时,第N点坐标值按球坐标下声波方程行波解展开M项的表达;c1 c2...cM代表选取的N个点上对应的声压去波展开系数;d1 d2...dM代表选取的N个点上对应的声压来波展开系数;
系数向量:
其中:右上角“+”代表求伪逆运算;
把每次的重构结果与真实测量值做误差运算,求出重构误差Δ:
获得展开项数-误差值(M-Δ)曲线,由此可以找出最小误差对应的展开项数,即为最优展开项数Mopti;
c、分离由目标声源和背景噪声源组成的声场;
把最优展开项数带入式(2)可以确定系数向量,目标声源产生的声波是去波,由此可以获得目标声源在空间任意点的声压值,故分离出的目标声源声压辐射为:
本发明的技术构思,在全息测量面上合理选取部分测点重构其余测点的声压,通过重构值与实际测量值的误差运算获得扩展项数-误差值(M-Δ)曲线,由M-Δ曲线中最小误差对应的展开项数确定最优展开项数Mopti,并由最优展开项数确定系数向量,再分离出目标声源声压辐射。本发明在基于球坐标下声波方程行波解从低阶到高阶的扩展项预先把声场定义为去波和来波区别运算。
通过本发明所述的方法,能够将混合声场中的目标声源所辐射的声压分离开,大幅减弱了背景噪声对目标声源的干扰,能够对混合声场进行有效的分离。由于声场分离基于正交的球波函数扩展,所以全息测量面可以是不规则的,因而该方法即可应用于平面、柱面和球面等规则声源面和背景噪声产生的声场分离,也可以应用于其它非规则几何面的声源和背景噪声产生的声场。
本发明的有益效果:
1、本发明的方法只需要全息面的声压值,不需要介质粒子速度矢量,进行非自由声场分离的计算,因而本发明可应用于稳态、非稳态和瞬态情况,并且计算容易、速度快。
2、本发明只需要单层全息测量面,并且全息面可以是不规则的,即全息面可以不是平面、柱面或球面。因此,本发明的方法易于实施,具有广泛用于环境噪声中目标声源的测量和评估的前景。
附图说明
图1为全息测量面测点分布示意图;
图2为声场分离方法示意图;
图3为实例中声场分离示意图。
图4为任意一空间平面上-声场分离面的总声压幅值分布图;
图5为目标声源单独辐射到声场分离面上的声压幅值分布图);
图6为采用本发明分离出来的目标声源辐射到声场分离面上的声压幅值分布图。
具体实施方式
参照图1、图2,不规则单全息声压测量面分离非自由声场的方法,按如下步骤进行:
a、测量全息测量面1上的声压值,在由目标声源3和背景噪声源4构成的声场中布置全息测量面1,在全息测量面1节点上放置传声器2;
b、选取全息测量面1部分测点5的声压值重构其余测点的声压值,确定最优展开项数;
把波动方程在球坐标下的完全解为:
在全息测量面1上选出部分测点,选取方法是每间隔一点选择一点,共选择N点,把选出点对应的声压按球坐标下声波方程行波解从低阶到高阶分别展开得式(2):
简记为:
其中:p1 p2...pN代表选取点上测得的声压值;代表当目标声源产生去波时,第N点坐标值按球坐标下声波方程行波解展开M项的表达;代表当背景噪声源产生来波时,第N点坐标值按球坐标下声波方程行波解展开M项的表达;c1 c2...cM代表选取的N个点上对应的声压去波展开系数;d1 d2...dM代表选取的N个点上对应的声压来波展开系数;
系数向量:
其中:右上角“+”代表求伪逆运算;
把每次的重构结果与真实测量值做误差运算,求出重构误差Δ:
获得展开项数-误差值(M-Δ)曲线,由此可以找出最小误差对应的展开项数,即为最优展开项数Mopti;
c、分离由目标声源和背景噪声源组成的声场;
把最优展开项数带入式(2)可以确定系数向量,目标声源产生的声波是去波,由此可以获得目标声源在空间任意点的声压值,故分离出的目标声源声压辐射为:
本发明的技术构思,在全息测量面上合理选取部分测点重构其余测点的声压,通过重构值与实际测量值的误差运算获得扩展项数-误差值(M-Δ)曲线,由M-Δ曲线中最小误差对应的展开项数确定最优展开项数Mopti,并由最优展开项数确定系数向量,再分离出目标声源声压辐射。本发明在基于球坐标下声波方程行波解从低阶到高阶的扩展项预先把声场定义为去波和来波区别运算。
通过本发明所述的方法,能够将混合声场中的目标声源所辐射的声压分离开,大幅减弱了背景噪声对目标声源的干扰,能够对混合声场进行有效的分离。由于声场分离基于正交的球波函数扩展,所以全息测量面可以是不规则的,因而该方法即可应用于平面、柱面和球面等规则声源面和背景噪声产生的声场分离,也可以应用于其它非规则几何面的声源和背景噪声产生的声场。
应用例如下:
在空间布置两个同心、半径不等的脉动膨胀球和脉动球形空腔。由于目标声源3(脉动膨胀球)和背景噪声源4(脉动球形空腔)都是简单声源,由本发明算法可知,只需两个到球心距离不等的测点,即可获得球形来波和去波扩展函数的系数。按图3布置测点A、B。采用本发明的声场分离方法实施声场分离,并与其理论解析解进行比较。
对于半径为a的脉动膨胀球,其在半径r处声压的解析解为:
式中:空气密度ρ=1.29kg/m3,脉动膨胀球均布径向振速v=0.5m/s、脉动球形空腔均布径向振速v=0.05m/s,两声源振动频率均为500Hz。
为便于表达,在空腔中取一空间平面作为声场分离面,来检验声场分离的效果。声场分离面大小为0.52m×0.52m,到坐标原点的垂直距离r0=0.2m。目标声源和背景噪声源的球心均在坐标原点,脉动膨胀球的半径r1=0.05m,脉动球形空腔的半径r2=0.5m。测点A到坐标原点的距离rA=0.3m,测点B到坐标原点的距离rB=0.37m。
图4是目标声源(脉动膨胀球)和背景噪声源(脉动球形空腔)共同辐射到声场分离面上的声压幅值分布;图5是当参数保持不变,目标声源单独辐射到声场分离面上的声压幅值分布;采用本发明对声场实施分离,得出图6,即采用本发明分离出来的目标声源,辐射到声场分离面上声压幅值分布。由图5和图6可以看出,分离后的结果与理论精确值完全相同。故本发明能够把目标声源和背景噪声源完全分离开。
对于简单的声源,本发明能够以很少的传声器分离声场;对于较复杂的声场环境,通过适量增加阵列的传声器数量并计算出最优展开项数,仍然能够有效地将声场分离。
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举,本发明的保护范围的不应当被视为仅限于实施例所陈述的具体形式,本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。
Claims (1)
1.不规则单全息声压测量面分离非自由声场的方法,其特征在于:按如下步骤进行:
a、测量全息测量面上的声压值,在由目标声源和背景噪声源构成的声场中布置全息测量面,在全息测量面节点上放置传声器;
b、选取全息测量面部分测点的声压值重构其余测点的声压值,确定最优展开项数;
把波动方程在球坐标下的完全解为:
在全息测量面上选出部分测点,选取方法是每间隔一点选择一点,共选择N点,把选出点对应的声压按球坐标下声波方程行波解从低阶到高阶分别展开得式(2):
简记为:
其中:p1 p2...pN代表选取点上测得的声压值;代表当目标声源产生去波时,第N点坐标值按球坐标下声波方程行波解展开M项的表达;代表当背景噪声源产生来波时,第N点坐标值按球坐标下声波方程行波解展开M项的表达;c1 c2...cM代表选取的N个点上对应的声压去波展开系数;d1 d2...dM代表选取的N个点上对应的声压来波展开系数;
系数向量:
其中:右上角“+”代表求伪逆运算;
把每次的重构结果与真实测量值做误差运算,求出重构误差Δ:
获得展开项数-误差值(M-Δ)曲线,由此可以找出最小误差对应的展开项数,即为最优展开项数Mopti;
c、分离由目标声源和背景噪声源组成的声场;
把最优展开项数带入式(2)可以确定系数向量,目标声源产生的声波是去波,由此可以获得目标声源在空间任意点的声压值,故分离出的目标声源声压辐射为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110049761 CN102183298B (zh) | 2011-03-02 | 2011-03-02 | 不规则单全息声压测量面分离非自由声场的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110049761 CN102183298B (zh) | 2011-03-02 | 2011-03-02 | 不规则单全息声压测量面分离非自由声场的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102183298A CN102183298A (zh) | 2011-09-14 |
CN102183298B true CN102183298B (zh) | 2012-12-12 |
Family
ID=44569530
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110049761 Active CN102183298B (zh) | 2011-03-02 | 2011-03-02 | 不规则单全息声压测量面分离非自由声场的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102183298B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217214B (zh) * | 2013-03-18 | 2014-09-17 | 哈尔滨工程大学 | 一种采用双面振速测量和局部声全息法的分离声场方法 |
CN103389155B (zh) * | 2013-06-26 | 2015-05-27 | 浙江工业大学 | 声品质客观参量三维空间分布数字图像生成方法 |
CN103674223B (zh) * | 2013-12-13 | 2017-02-01 | 中联重科股份有限公司渭南分公司 | 噪声源测试装置、***及方法 |
CN103674232B (zh) * | 2013-12-13 | 2017-04-12 | 中联重科股份有限公司渭南分公司 | 用于噪声源测试的声强获取装置 |
CN105136280B (zh) * | 2015-09-18 | 2018-04-20 | 华晨汽车集团控股有限公司 | 一种多源噪声环境下测试单一噪声品质的***及方法 |
CN107478325B (zh) * | 2017-09-07 | 2019-07-26 | 安徽大学 | 一种非稳态平面声源的自由场还原方法 |
CN108051800B (zh) * | 2017-12-13 | 2022-01-28 | 贵州航天计量测试技术研究所 | 基于球面近场声全息重构无功声强的室内噪声源定位方法 |
CN110579275B (zh) * | 2019-10-21 | 2022-03-11 | 南京南大电子智慧型服务机器人研究院有限公司 | 一种基于球形矢量传声器阵列实现声场分离的方法 |
CN110824427A (zh) * | 2019-10-21 | 2020-02-21 | 中国人民解放军陆军勤务学院 | 一种斜三棱锥声压测试装置及其空间次声源定向定位方法 |
CN110763328B (zh) * | 2019-11-18 | 2021-10-22 | 湖北文理学院 | 一种半空间声场重建方法和装置 |
CN113591303A (zh) * | 2021-07-29 | 2021-11-02 | 广东电网有限责任公司 | 一种噪声预测方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7068801B1 (en) * | 1998-12-18 | 2006-06-27 | National Research Council Of Canada | Microphone array diffracting structure |
CN1952627A (zh) * | 2006-10-27 | 2007-04-25 | 合肥工业大学 | 一种声场分离方法 |
CN101403634A (zh) * | 2008-10-20 | 2009-04-08 | 合肥工业大学 | 压力速度法声场分离方法 |
CN101539455A (zh) * | 2009-04-22 | 2009-09-23 | 合肥工业大学 | 采用移动等效源法重建移动声源的方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7330396B2 (en) * | 2004-10-13 | 2008-02-12 | Wayne State University | Farfield analysis of noise sources |
WO2007045016A1 (en) * | 2005-10-20 | 2007-04-26 | Personal Audio Pty Ltd | Spatial audio simulation |
-
2011
- 2011-03-02 CN CN 201110049761 patent/CN102183298B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7068801B1 (en) * | 1998-12-18 | 2006-06-27 | National Research Council Of Canada | Microphone array diffracting structure |
CN1952627A (zh) * | 2006-10-27 | 2007-04-25 | 合肥工业大学 | 一种声场分离方法 |
CN101403634A (zh) * | 2008-10-20 | 2009-04-08 | 合肥工业大学 | 压力速度法声场分离方法 |
CN101539455A (zh) * | 2009-04-22 | 2009-09-23 | 合肥工业大学 | 采用移动等效源法重建移动声源的方法 |
Non-Patent Citations (4)
Title |
---|
JP特表2009-512364A 2009.03.19 |
张永斌等.基于声压-振速测量的单面声场分离技术.《物理学报》.2009,第58卷(第12期),全文. * |
李卫兵等.球坐标系下入射与散射声场分离理论.《中国科学 E辑:技术科学》.2007,第37卷(第1期),全文. * |
毕传兴等.半自由声场的全息重建和预测技术研究.《声学学报》.2004,第29卷(第4期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN102183298A (zh) | 2011-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102183298B (zh) | 不规则单全息声压测量面分离非自由声场的方法 | |
CN107576388B (zh) | 一种浅海信道下三维结构声源辐射声场预报方法 | |
CN103592022B (zh) | 采用声压和质点加速度测量的实时声场分离方法 | |
CN101403634B (zh) | 压力速度法声场分离方法 | |
CN103267571B (zh) | 采用单面声压和振速测量的实时声场分离方法 | |
CN100442030C (zh) | 一种声场分离方法 | |
CN101566496B (zh) | 采用双面振速测量和等效源法分离声场的方法 | |
CN102121847B (zh) | 一种瞬态声场重建方法 | |
CN103743469B (zh) | 压电水声换能器声辐射模态测量方法及*** | |
CN102494755B (zh) | 一种非稳态声场分离方法 | |
CN107478325B (zh) | 一种非稳态平面声源的自由场还原方法 | |
CN104596636B (zh) | 声场分离方法 | |
CN109764956B (zh) | 一种基于组合优化正则化方法的近场声全息方法 | |
Vold et al. | Axisymmetry and azimuthal modes in jet noise | |
CN109211399A (zh) | 一种锥面声源统计近场声全息方法 | |
CN112926231B (zh) | 一种基于等效源法的有限空间中近场声全息测量方法 | |
CN108154557B (zh) | 一种基于家居环境的静音区域的三角面片化方法 | |
CN114252148B (zh) | 一种基于长椭球波叠加的声场重建方法 | |
Kletschkowski et al. | Identification of noise sources in an aircraft fuselage using an inverse method based on a finite element model | |
Kasatkin et al. | Spectral-Correlation Signal Processing in the Infrasonic Frequency Range | |
Yang et al. | A hybrid wave superposition method based on particle filter | |
Comesaña et al. | Designing an array for performing Near-field Acoustic Holography with a small number of pu probes | |
Ostashev et al. | Overview of acoustic tomography of the atmosphere | |
CN116559850A (zh) | 入射声场与散射声场分离方法、***、装置及其存储介质 | |
Sinha et al. | Study of acoustic vector sensor based direction of arrival estimation of in-air maneuvering tonal source |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |