专利名称：压力速度法声场分离方法
技术领域：
本发明涉及物理专业中噪声类领域声场分离方法，具体的说是一种压力速度 法声场分离方法。
背景技术：
在实际测量时，通常会遇到测量面两侧都有声源，或是测量面的一侧存在反 射或散射。而实际工程中，为了更加准确地研究实际声源的声辐射特性或反射面 的反射特性，需要将来自测量面两侧的辐射声分离开来。G. Weinreich等在1980 年提出采用两个靠得很近测量面测量的方法来实现入射波和辐射波的分离。G. V. Frisk在E. G. Williams等提出的近场声全息技术和G. Weinreich等提出的双 面测量方法的基础上，建立了基于空间FFT法的双测量面声场分离方法。该方法 在随后的一段时间里得到进一步应用和发展。M. Tamura详细建立了通过双面测 量，再采用二维空间FFT法实现的声场分离公式，并通过数值仿真和实验成功求 得反射界面的反射系数。Z. Hu和J. S. Bolton也对采用该方法测量平面波反 射系数进行了进一步验证。M. T. Cheng等建立了迪卡尔坐标和柱面坐标下的双 测量面声场分离公式，并用于实现散射声场的分离，分析了该方法分离散射场的 敏感性。F. Yu等成功采用该方法分离近场声全息测量过程中全息面上来自背向 的噪声。基于空间FFT法的双测量面声场分离方法有其固有的缺陷 一方面对测 量面的形状有限制，即只能是平面、柱面或球面等规则形状；另一方面受到FFT 算法的影响，分离误差较大，尤其是在来自测量面两侧声压相差较大时，其误差 尤为明显。F. Jacobsen等在J. Hald提出的方法的基础上，提出了基于声压和速度测 量的统计最优声场分离方法，该方法采用Microflow公司的p-u声强探头同时测 量声场全息面上的声压和质点振速，再采用建立的联合求解公式实现来自全息面 两侧的辐射声场分离。该方法的缺陷对测量面的形状有限制，即只能是平面、 柱面或球面等规则形状。C. Langrenne等提出一种基于边界元法的双面声场分离方法。该方法首先 测量两个包络声源的平行等间距测量面上的声压；再采用Helmholtz积分法分离入射和辐射声压场。该方法的缺陷存在奇异积分、解的非唯一性等处理，计算效率低。
发明内容
本发明所解决的技术问题是避免上述现有技术所存在的不足之处，提供一种 实现方便、适用于任意形状测量面、计算稳定性好、计算精度高的采用等效源法 实现、压力速度测量的声场分离方法。本发明解决技术问题所采用的技术方案是本发明方法按如下步骤进行a、 测量测量面S上的声压信息在由声源1和声源2构成的被测声场中，位于声源1与声源2之间有测量面 S在测量面上分别分布有测量网格点，相邻网格点之间的距离小于半个波长； 测量测量面上各网格点处的声压幅值和相位信息获得测量面上的声压；所述被测 声场为稳态声场；b、 在测量面5"与声源1之间设定虚源面《，在测量面S与声源2之间设定虚源面&',并在两虚源面上分别分布有等效源，等效源的个数不大于对应测量面网络点数；所述等效源为标准点源、面源或体源；c、 建立等效源与所述测量面5"上声压和质点速度之间的传递关系A = A、v2=v;『2， 其中p,为声源1在测量面5"上所辐射的声压、 p2为声源2在测量面5"上所 辐射的声压、v,为声源1在测量面S上所辐射的质点速度、^为声源2在测量面S上辐 射的质点速度、K为虚源面《上等效源权重矢量、 罵为虚源面《上等效源权重矢且 里、/^为虚源面《上等效源与测量面S上声压之间的传递矩阵、 v:为虚源面《上等效源与测量面S上质点速度之间的传递矩阵、p;为虚源面《上等效源与测量面5"上声压之间的传递矩阵、V;为虚源面《上等效源与测量面5"上质点速度之间的传递矩阵；d、 建立测量面5"上声压和质点速度与两侧声源辐射声压和质点速度的关系/ = A + />2V二Vj+V2， 其中/7为测量面5"上测量的声压、 V为测量面S上测量的质点速度；e、 分离测量面6"上由两侧声源辐射的声压和质点速度根据步骤c所建立的传递关系和步骤d所建立的测量面5"上声压和质点速度 与两侧声源辐射声压和质点速度的关系，联合求解获得测量面5"上由声源1和声 源2分别辐射的声压和质点速度Pl = Pi* (A* - AW )+ O - G2V) Vl =《(A - G2V1* )+ _ G2V)p2 ";o;o-《v)"2=v; —)+o -《v)，其中各网格点上的声压和质点速度幅值和相位信息的测量是采用单个或多个传 声器和质点速度传感器在测量面上分别扫描、采用声强探头阵列在测量面上快 照、或采用单个或多个声强探头阵列在测量面上扫描获得。测量面^为任意形状面。本发明方法是测量测量面s上的声压和质点速度，采用等效源法来实现测量 面上由两侧声源辐射声压和质点速度的分离。 本发明的理论模型等效源法的基本思想是采用分布在声源内部的一系列等效源加权叠加来近 似实际声场，此时只需确定这些等效源的源强即可预测整个声场。在实际求解过 程中，等效源的源强可以通过测量的声源的边界条件(声压或法向振速)来确定。 对于声场中任一个测量面,也可以通过在该面的背离分析域内虚源面上分布等效 源来近似在该面前方区域的辐射声场。参见图1，测量面5"的右侧区域场点r处声学量可以由分布在该面左侧虚源 面5>御一系列等效源近似获得。设测量面S和虚源面6V上分别分布了 #个测量 点和W个等效源，第i个等效源在场点r处的辐射声压为A'(r)和粒子速度为v,'(r)，则场点r处的实际辐射声压可表示为<formula>formula see original document page 7</formula> (1)<formula>formula see original document page 7</formula>(2),=1式中w,为/个等效源所对应的源强。各等效源的源强由测量面的边界条件确定， 由等式(1)可得测量面5"上i/个测量点的声压和质点速度可以表示为<formula>formula see original document page 7</formula>(3)式中，式(3)写成矩阵形式为<formula>formula see original document page 7</formula>(4)K化)<formula>formula see original document page 7</formula>(5)w，<formula>formula see original document page 7</formula>(6)式中，P为声场中#个测点处的声压列向量；^为对应iV个等效源所占的权重系 数列向量；/为yV个等效源与i/个测点处声压之间的MxiV阶传递矩阵。同理，测量面5"上#个测量点的质点速度可以表示为 <formula>formula see original document page 8</formula> (7)式中，v为声场中#个测点处的质点速度列向量；v'为#个等效源与#个测点处声压之间的M x W阶传递矩阵。求得权重系数矩阵r后，由式(1)和式(2)就可以计算声场中任意一点的声压 与振速，实现声场的预测。由上可知，声场中测量面一侧的辐射声场可以通过在该测量面的另一侧分布 一系列等效源来近似。如果测量面的两侧都有声源，则测量面上的声压和质点速 度为两侧声源辐射声压的组合。参见图2，测量面S上的声压为<formula>formula see original document page 8</formula> (8) 式中A为声源1在测量面5!上所辐射的声压，A为声源2在测量面5"上所辐射 的声压。测量面^上的质点速度表示为<formula>formula see original document page 8</formula> (9)式中、为声源1在测量面5"上所辐射的声压，"2为声源2在测量面5"上所辐射的 声压。由于很难直接将测量面上两侧声源辐射的声压和质点速度分开。本发明 的方法是在同一个面上同时测量声压和质点速度，然后再通过等效源法来实 现分离。由上可知，测量面5"上声源1所辐射声压A和质点速度^可以通过在测 量面S与声源1之间设置的虚源面《上分布一系列等效源来近似，测量面 上声源2所辐射声压/;2和质点速度^可以通过在测量面5"与声源2之间设 置的虚源面《上分布一系列等效源来近似。(10)<formula>formula see original document page 8</formula>(11)<formula>formula see original document page 9</formula>式中，p,为声源l在测量面5"上所辐射的声压，^为声源2在测量面6"上所辐 射的声压，^为声源1在测量面5"上所辐射的质点速度，^为声源2在测量面 5"上辐射的质点速度，K为虚源面《上等效源权重矢量，然为虚源面《上等效源 权重矢量，A'为虚源面《上等效源与测量面5"上声压之间的传递矩阵，v;为虚 源面《上等效源与测量面S上质点速度之间的传递矩阵，为虚源面《上等效 源与测量面5"上声压之间的传递矩阵，v;为虚源面《上等效源与测量面S上质 点速度之间的传递矩阵。将式(IO)、 (11)、 (12)和(13)分别代入式(8)和(9)，则可获得如下关系<formula>formula see original document page 9</formula>联合式(14)和(15)，可得<formula>formula see original document page 9</formula>W2=(p*2-GlV;)+(p-GlV)式中<formula>formula see original document page 9</formula> 式中上标"+"表示矩阵广义逆。 *将式(16)和(17)分别代入式(10)、 (11)、 (12)和(13)可以分别获得测量面S 上声源1和声源2所辐射的声压和质点速度。通过上述方法，实现了测量面上声压和质点速度的分离，可以获得来自测量 面两侧声源的辐射声压和质点速度。与已有技术相比，本发明的有益效果1、 本发明所采用的测量面，可以是任意形状测量面，解决了传统方法只能 适用于平面、柱面或球面等规则形状的缺陷。2、 本发明采用等效源法来作为声场分离算法，与传统方法相比，本发明方 法具有计算稳定性好、计算精度高等优点。3、 本发明方法测量只需在一个测量面上进行。4、 本发明方法实施简单，可以广泛用于内部声场或噪声环境下的近场声全 息测量、噪声环境下的声强法声源识别、材料反射系数的测量，散射声场的分离 等。
-图1为平面声源等效源位置分布图； 图2为压力速度法声场分离示意图； 图3 (a)测量面S上实际测量的声压幅值分布； 图3 (b)测量面S上声源1辐射声压理论幅值分布； 图3 (c)采用本发明方法分离出来的声源l辐射声压幅值分布； 图3 (d)测量面S上实际测量的声压相位分布； 图3 (e)测量面S上声源1辐射声压理论相位分布； 图3 (f)采用本发明方法分离出来的声源1辐射声压相位分布； 图4 (a)测量面S上实际测量的质点速度幅值分布； 图4 (b)测量面S上声源1辐射质点速度理论幅值分布； 图4 (c)采用本发明方法分离出来的声源l辐射质点速度幅值分布； 图4 (d)测量面S上实际测量的质点速度相位分布； 图4 (e)测量面S上声源1辐射质点速度理论相位分布； 图4 (f)采用本发明方法分离出来的声源l辐射质点速度相位分布； 以下通过具体实施方式
，并结合附图对本发明作进一步描述。
具体实施例方式参见图2，本实施例中，测量面两侧均分布有声源，其中声源l为主声源， 声源2为噪声源或反射、散射源，在由声源1和声源2构成的被测声场中，位于 声源1与声源2之间有测量面&在测量面上分别分布有测量网格点，相邻网格 点之间的距离小于半个波长。具体实施步骤为a、 采用单个或多个传声器和质点速度传感器在测量面上分别扫描、采用声强探头阵列在测量面上快照、或采用单个或多个声强探头阵列在测量面上扫描获 得测量面^上的声压和质点速度信息；b、 在测量面6"与声源1之间设定虚源面《，在测量面5"与声源2之间设定 虚源面《，并在两虚源面上分别分布有等效源，等效源的个数不大于对应测量面网络点数；所述等效源为标准点源、面源或体源；c、 建立等效源与所述两测量面之间的传递关系Pi =Vl = wv2=v;W2， 其中 ^为声源1在测量面5"上所辐射的声压、/ 2为声源2在测量面5"上所 辐射的声压、、为声源1在测量面5"上所辐射的质点速度、^为声源2在测量面5"上辐 射的质点速度、K为虚源面《上等效源权重矢量、 然为虚源面S;上等效源权重矢/^为虚源面《上等效源与测量面S上声压之间的传递矩阵、 v:为虚源面《上等效源与测量面S上质点速度之间的传递矩阵、 为虚源面&'上等效源与测量面5"上声压之间的传递矩阵、v;为虚源面《上等效源与测量面6"上质点速度之间的传递矩阵；d、 建立测量面5"上声压和质点速度与两侧声源辐射声压和质点速度的关系P = A +V = V1+V2， 其中p为测量面^上测量的声压、 V为测量面5"上测量的质点速度； e、分离测量面5"上由两侧声源辐射的声压和质点速度根据步骤c所建立的传递关系和步骤d所建立的测量面^上声压和质点速度与两侧声源辐射声压和质点速度的关系，联合求解获得测量面5"上由声源1和声源2分别辐射的声压和质点速度<formula>formula see original document page 12</formula>
方法的检验在测量面两侧各布置一个脉动球,分别采用本发明的声场分离方法实施声场 分离，并与其解析解比较。对于单个半径为"的脉动球，其在场点r处声压的解析解为W,0=-v.^^.exp[汰(r-")], (22) 由)式中，均匀径向速度v: 1 m/s，空气密度为yc^ 1.2 kg/m3,声源振动频率为 1kHz.两测量面的位置关系参见图2。测量面均为lmXlm的平面，测量面上均匀 地分布21X21个测量点。声源1为位于(-O. 3, 0, O)m处的脉动球，声源2为位 于(0.3，-0.25，0.2)m处的脉动球。测量面与虚源面之间的距离1和2均为 0. lm。此处声源1为主声源，声源2为噪声源，需要将测量面5"上声源1辐射声 压和质点速度分离出来。参见图3 (其中，图3(a)和图3(d)为测量面5"上实际测量的声压幅值和相位 分布，图3(b)和图3(e)为声源l在测量面5"上辐射的声压幅值和相位分布，图 3(c)和图3(f)为采用本发明方法分离出来的声源1在测量面5"上辐射的声压幅 值和相位分布；图4(a)和图4(d)为测量面5"上实际测量的质点速度幅值和相位分布，图4(b)和图4(e)为声源1在测量面5上辐射的质点速度幅值和相位分布,图4(c)和图4(f)为采用本发明方法分离出来的声源1在测量面5"上辐射的质点速度幅值和相位分布)，采用本发明方法实施分离后，可以准确得到声源1在测量面5"上辐射信息，分离出声压和质点速度的幅值和相位分布与其理论值非常吻 合。为了更加定量描述本发明方法的分离精度，定义分离误差百分比为"=、s(w-帥2 / J力对x100 (%)， (23)V /=1 /V i=l式中，w为所有声源的表面结点总数，p,和A分别为对应第f个测量点处分离的和理论的声压或质点速度。由式(23)计算可得，声压和质点速度分离误差百分比 分别为0. 43%和0. 35%，显然采用本发明方法能获得非常精确的结果。
权利要求
1、压力速度法声场分离方法，其特征是在于按如下步骤进行a、测量测量面S上的声压信息在由声源1和声源2构成的被测声场中，位于声源1与声源2之间有测量面S；在测量面上分别分布有测量网格点，相邻网格点之间的距离小于半个波长；测量测量面上各网格点处的声压幅值和相位信息获得测量面上的声压；所述被测声场为稳态声场；b、在测量面S与声源1之间设定虚源面S1*，在测量面S与声源2之间设定虚源面S2*，在两虚源面上分别分布有等效源，等效源的个数不大于对应测量面网络点数；所述等效源为标准点源、面源或体源；c、建立等效源与所述测量面S上声压和质点速度之间的传递关系<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><msub> <mi>p</mi> <mn>1</mn></msub><mo>=</mo><msubsup> <mi>p</mi> <mn>1</mn> <mo>*</mo></msubsup><msub> <mi>W</mi> <mn>1</mn></msub> </mrow>]]></math></maths><maths id="math0002" num="0002" ><math><![CDATA[ <mrow><msub> <mi>v</mi> <mn>1</mn></msub><mo>=</mo><msubsup> <mi>v</mi> <mn>1</mn> <mo>*</mo></msubsup><msub> <mi>W</mi> <mn>1</mn></msub> </mrow>]]></math></maths><maths id="math0003" num="0003" ><math><![CDATA[ <mrow><msub> <mi>p</mi> <mn>2</mn></msub><mo>=</mo><msubsup> <mi>p</mi> <mn>2</mn> <mo>*</mo></msubsup><msub> <mi>W</mi> <mn>2</mn></msub> </mrow>]]></math></maths><maths id="math0004" num="0004" ><math><![CDATA[ <mrow><msub> <mi>v</mi> <mn>2</mn></msub><mo>=</mo><msubsup> <mi>v</mi> <mn>2</mn> <mo>*</mo></msubsup><msub> <mi>W</mi> <mn>2</mn></msub><mo>,</mo> </mrow>]]></math> id="icf0004" file="A2008101947540002C4.tif" wi="19" he="4" top= "149" left = "49" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>其中p1为声源1在测量面S上所辐射的声压、p2为声源2在测量面S上所辐射的声压、v1为声源1在测量面S上所辐射的质点速度、v2为声源2在测量面S上辐射的质点速度、W1为虚源面S1*上等效源权重矢量、W2为虚源面S2*上等效源权重矢量、p1*为虚源面S1*上等效源与测量面S上声压之间的传递矩阵、v1*为虚源面S1*上等效源与测量面S上质点速度之间的传递矩阵、p2*为虚源面S2*上等效源与测量面S上声压之间的传递矩阵、v2*为虚源面S2*上等效源与测量面S上质点速度之间的传递矩阵；d、建立测量面S上声压和质点速度与两侧声源辐射声压和质点速度的关系p＝p1+p2v＝v1+v2，其中p为测量面S上测量的声压、v为测量面S上测量的质点速度；e、分离测量面S上由两侧声源辐射的声压和质点速度根据步骤c所建立的传递关系和步骤d所建立的测量面S上声压和质点速度与两侧声源辐射声压和质点速度的关系，联合求解获得测量面S上由声源1和声源2分别辐射的声压和质点速度<maths id="math0005" num="0005" ><math><![CDATA[ <mrow><msub> <mi>p</mi> <mn>1</mn></msub><mo>=</mo><msubsup> <mi>p</mi> <mn>1</mn> <mo>*</mo></msubsup><msup> <mrow><mo>(</mo><msubsup> <mi>p</mi> <mn>1</mn> <mo>*</mo></msubsup><mo>-</mo><msub> <mi>G</mi> <mn>2</mn></msub><msubsup> <mi>v</mi> <mn>1</mn> <mo>*</mo></msubsup><mo>)</mo> </mrow> <mo>+</mo></msup><mrow> <mo>(</mo> <mi>p</mi> <mo>-</mo> <msub><mi>G</mi><mn>2</mn> </msub> <mi>v</mi> <mo>)</mo></mrow> </mrow>]]></math></maths><maths id="math0006" num="0006" ><math><![CDATA[ <mrow><msub> <mi>v</mi> <mn>1</mn></msub><mo>=</mo><msubsup> <mi>v</mi> <mn>1</mn> <mo>*</mo></msubsup><msup> <mrow><mo>(</mo><msubsup> <mi>p</mi> <mn>1</mn> <mo>*</mo></msubsup><mo>-</mo><msub> <mi>G</mi> <mn>2</mn></msub><msubsup> <mi>v</mi> <mn>1</mn> <mo>*</mo></msubsup><mo>)</mo> </mrow> <mo>+</mo></msup><mrow> <mo>(</mo> <mi>p</mi> <mo>-</mo> <msub><mi>G</mi><mn>2</mn> </msub> <mi>v</mi> <mo>)</mo></mrow> </mrow>]]></math></maths><maths id="math0007" num="0007" ><math><![CDATA[ <mrow><msub> <mi>p</mi> <mn>2</mn></msub><mo>=</mo><msubsup> <mi>p</mi> <mn>2</mn> <mo>*</mo></msubsup><msup> <mrow><mo>(</mo><msubsup> <mi>p</mi> <mn>2</mn> <mo>*</mo></msubsup><mo>-</mo><msub> <mi>G</mi> <mn>1</mn></msub><msubsup> <mi>v</mi> <mn>2</mn> <mo>*</mo></msubsup><mo>)</mo> </mrow> <mo>+</mo></msup><mrow> <mo>(</mo> <mi>p</mi> <mo>-</mo> <msub><mi>G</mi><mn>1</mn> </msub> <mi>v</mi> <mo>)</mo></mrow> </mrow>]]></math></maths><maths id="math0008" num="0008" ><math><![CDATA[ <mrow><msub> <mi>v</mi> <mn>2</mn></msub><mo>=</mo><msubsup> <mi>v</mi> <mn>2</mn> <mo>*</mo></msubsup><msup> <mrow><mo>(</mo><msubsup> <mi>p</mi> <mn>2</mn> <mo>*</mo></msubsup><mo>-</mo><msub> <mi>G</mi> <mn>1</mn></msub><msubsup> <mi>v</mi> <mn>2</mn> <mo>*</mo></msubsup><mo>)</mo> </mrow> <mo>+</mo></msup><mrow> <mo>(</mo> <mi>p</mi> <mo>-</mo> <msub><mi>G</mi><mn>1</mn> </msub> <mi>v</mi> <mo>)</mo></mrow><mo>,</mo> </mrow>]]></math> id="icf0008" file="A2008101947540003C4.tif" wi="49" he="4" top= "140" left = "52" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>其中<maths id="math0009" num="0009" ><math><![CDATA[ <mrow><msub> <mi>G</mi> <mn>1</mn></msub><mo>=</mo><msubsup> <mi>p</mi> <mn>1</mn> <mo>*</mo></msubsup><msup> <mrow><mo>(</mo><msubsup> <mi>v</mi> <mn>1</mn> <mo>*</mo></msubsup><mo>)</mo> </mrow> <mo>+</mo></msup><mo>,</mo> </mrow>]]></math> id="icf0009" file="A2008101947540003C5.tif" wi="23" he="4" top= "151" left = "35" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths><maths id="math0010" num="0010" ><math><![CDATA[ <mrow><msub> <mi>G</mi> <mn>2</mn></msub><mo>=</mo><msubsup> <mi>p</mi> <mn>2</mn> <mo>*</mo></msubsup><msup> <mrow><mo>(</mo><msubsup> <mi>v</mi> <mn>2</mn> <mo>*</mo></msubsup><mo>)</mo> </mrow> <mo>+</mo></msup><mo>.</mo> </mrow>]]></math> id="icf0010" file="A2008101947540003C6.tif" wi="26" he="4" top= "151" left = "87" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>
2、 根据权利要求i所述的声场分离方法，其特征在于所述各网格点上的声 压和质点速度幅值和相位信息的测量是采用单个或多个传声器和质点速度传感 器在测量面上分别扫描、采用声强探头阵列在测量面上快照、或采用单个或多个 声强探头在测量面上扫描获得。
3、 根据权利要求1所述的声场分离方法，其特征在于所述测量面5"为任意形状面。
4、 根据权利要求1所述的声场分离方法，其特征在于声场分离是采用测量面6"上的声压和质点速度作为输入量。
5、 根据权利要求1所述的声场分离方法，其特征在于在测量面S两侧分布等效源，采用等效源法实现声场分离。
全文摘要
本发明公开了压力速度法声场分离方法，在被测声场中设置测量面S；采用声压传感器测量测量面S上的声压，采用质点速度传感器或声强探头测量面S上的质点速度；在测量面S两侧设置虚源面S<sub>1</sub><sup>*</sup>和S<sub>2</sub><sup>*</sup>，在虚源面上分布等效源；建立等效源与测量面S上声压和质点速度之间的传递关系；根据该传递关系实现测量面上由两侧声源辐射声压和质点速度的分离。通过测量单个面上的声压和质点速度实现声场分离；采用等效源法来作为声场分离算法，计算稳定性好、计算精度高、实施简单。可以广泛用于内部声场或噪声环境下的近场声全息测量、噪声环境下的声强法声源识别、材料反射系数的测量，散射声场的分离。
文档编号G01H5/00GK101403634SQ20081019475
公开日2009年4月8日 申请日期2008年10月20日 优先权日2008年10月20日
发明者毕传兴 申请人:合肥工业大学