• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 17
    • 下载费用:30 金币  

    重庆时时彩有121期吗: 基于封闭空间几何信息建模的单传声器声源定位方法.pdf

    关 键 词:
    基于 封闭 空间 几何 信息 建模 传声器 声源 定位 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201611230305.7

    申请日:

    2016.12.28

    公开号:

    CN106772220A

    公开日:

    2017.05.31

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G01S 1/72申请日:20161228|||公开
    IPC分类号: G01S1/72 主分类号: G01S1/72
    申请人: 西北工业大学
    发明人: 王海涛; 曾向阳; 杜博凯; 刘延善; 王璐; 陈克安
    地址: 710072 陕西省西安市友谊西路127号
    优先权:
    专利代理机构: 西北工业大学专利中心 61204 代理人: 王鲜凯
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201611230305.7

    授权公告号:

    |||

    法律状态公告日:

    2017.06.23|||2017.05.31

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明公开了一种基于封闭空间几何信息建模的单传声器声源定位方法,用于解决现有封闭空间声源定位方法实用性差的技术问题。技术方案是在封闭空间边界所包围的流体区域内布置n个节点,将此n个节点从1到n进行编号。利用此n个节点构建出封闭空间及内部物体的几何形状,并根据此n个节点坐标信息建立描述封闭空间数学性质的系统矩阵;在实际声源定位时,利用单传声器测得多个频率下的声信号,将此声信号与系统矩阵共同运算,得到声信号的位置信息,从而实现声源定位,实用性好。

    权利要求书

    1.一种基于封闭空间几何信息建模的单传声器声源定位方法,其特征在于包括以下步
    骤:
    步骤一、在封闭空间边界所包围的流体区域内布置n个节点,将此n个节点从1到n进行
    编号;n个节点一方面作为封闭空间的几何信息在后期进行建模计算,一方面用来表示未来
    声源定位时的位置;
    n个节点近似均匀分布于整个流体区域内,封闭空间的边界处需布有节点;当封闭空间
    内部存在物体时,物体内部不布置节点,且物体边界上应布置有节点,即节点应勾勒出物体
    形状;节点布置完毕后,封闭空间及其内部物体的边界用Γ表示,边界所包围的流体区域用
    Ω表示;
    步骤二、假定在封闭空间内声源的位置为r,它在单位时间内向单位体积内的空间提供
    了ρ0q(r,t)的媒质质量;根据质量守恒定律,媒质中声波的连续方程写为:
    <mrow> <mfrac> <mrow> <mo>&part;</mo> <msup> <mi>&rho;</mi> <mo>&prime;</mo> </msup> </mrow> <mrow> <mo>&part;</mo> <mi>t</mi> </mrow> </mfrac> <mo>+</mo> <mi>d</mi> <mi>i</mi> <mi>v</mi> <mrow> <mo>(</mo> <msub> <mi>&rho;</mi> <mn>0</mn> </msub> <mi>v</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>&rho;</mi> <mn>0</mn> </msub> <mi>q</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
    式中,ρ'为媒质密度增量,ρ0表示媒质静态密度,q为q(r,t)的简写,v为媒质质点速度,t
    表示时间,div为散度算子,在三维空间笛卡尔坐标系中,
    除了连续性方程之外,用来描述媒质声波的基本方程还有两个,它们不受声源的影响,
    分别为运动方程:
    <mrow> <msub> <mi>&rho;</mi> <mn>0</mn> </msub> <mfrac> <mrow> <mo>&part;</mo> <mi>v</mi> </mrow> <mrow> <mo>&part;</mo> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mi>g</mi> <mi>r</mi> <mi>a</mi> <mi>d</mi> <mi> </mi> <mi>p</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
    和物态方程:
    <mrow> <mi>p</mi> <mo>=</mo> <msup> <msub> <mi>c</mi> <mn>0</mn> </msub> <mn>2</mn> </msup> <msup> <mi>&rho;</mi> <mo>&prime;</mo> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
    以上两式中,p代表声压,c0代表声速,grad为梯度算子,在三维空间笛卡尔坐标系中,

    由媒质中声波的三个基本方程得到有源情况下封闭空间中有关声压p的波动方程:
    <mrow> <msup> <mo>&dtri;</mo> <mn>2</mn> </msup> <mi>p</mi> <mo>-</mo> <mfrac> <mn>1</mn> <msubsup> <mi>c</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mfrac> <mfrac> <mrow> <msup> <mo>&part;</mo> <mn>2</mn> </msup> <mi>p</mi> </mrow> <mrow> <mo>&part;</mo> <msup> <mi>t</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <msub> <mi>&rho;</mi> <mn>0</mn> </msub> <mfrac> <mrow> <mo>&part;</mo> <mi>q</mi> </mrow> <mrow> <mo>&part;</mo> <mi>t</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
    式中,为拉普拉斯算子,在三维空间笛卡尔坐标系中,
    在频域,声源强度q(r,t)表示为
    q(r,t)=qω(r)ejωt (5)
    式中,ω为谐振频率,qω(r)为在位置r处频域内的声源强度;
    封闭空间声场各点的声压的频率与声源相同,声压表示为:
    p(r,t)=pω(r)ejωt (6)
    式中,pω(r)为在位置r处频域内的声压;
    将式(6)及式(7)代入式(5)中,得到简谐声源激励下的声波波动方程为:
    <mrow> <msup> <mo>&dtri;</mo> <mn>2</mn> </msup> <msub> <mi>p</mi> <mi>&omega;</mi> </msub> <mrow> <mo>(</mo> <mi>r</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mi>&omega;</mi> <mi>t</mi> </mrow> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mi>&omega;</mi> <msub> <mi>c</mi> <mn>0</mn> </msub> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msub> <mi>p</mi> <mi>&omega;</mi> </msub> <mrow> <mo>(</mo> <mi>r</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mi>&omega;</mi> <mi>t</mi> </mrow> </msup> <mo>+</mo> <msub> <mi>j&rho;</mi> <mn>0</mn> </msub> <msub> <mi>&omega;q</mi> <mi>&omega;</mi> </msub> <mrow> <mo>(</mo> <mi>r</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mi>&omega;</mi> <mi>t</mi> </mrow> </msup> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
    令式中k=ω/c0,称其为波数,并且消去ejωt,得到只依赖于空间坐标的那部分方程,即
    室内有源Helmholtz方程:
    <mrow> <msup> <mo>&dtri;</mo> <mn>2</mn> </msup> <msub> <mi>p</mi> <mi>&omega;</mi> </msub> <mrow> <mo>(</mo> <mi>r</mi> <mo>)</mo> </mrow> <mo>+</mo> <msup> <mi>k</mi> <mn>2</mn> </msup> <msub> <mi>p</mi> <mi>&omega;</mi> </msub> <mrow> <mo>(</mo> <mi>r</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>j</mi> <msub> <mi>&rho;</mi> <mn>0</mn> </msub> <mi>&omega;</mi> <msub> <mi>q</mi> <mi>&omega;</mi> </msub> <mrow> <mo>(</mo> <mi>r</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
    这样就将声压的时域问题转换为频域问题,式(9)即为封闭空间声场的控制方程;
    在封闭空间中,边界具有吸声能力,其声压梯度表示为:
    <mrow> <mfrac> <mrow> <mo>&part;</mo> <mi>p</mi> </mrow> <mrow> <mo>&part;</mo> <mi>n</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mi>j</mi> <mi>k</mi> <mi>p</mi> </mrow> <mi>&zeta;</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
    式中,n为封闭空间壁面外法线方向,ζ称为比声阻抗,满足下式:
    <mrow> <mi>&zeta;</mi> <mo>=</mo> <mfrac> <mi>Z</mi> <mrow> <msub> <mi>&rho;</mi> <mn>0</mn> </msub> <msub> <mi>c</mi> <mn>0</mn> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
    式中,Z为界面声阻抗;
    根据Galerkin型加权残量法,为了求解式(8),首先设一试函数为代入有源
    Helmholtz方程及其边界条件,试函数产生残量R和
    <mrow> <mi>R</mi> <mo>=</mo> <msup> <mo>&dtri;</mo> <mn>2</mn> </msup> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>+</mo> <msup> <mi>k</mi> <mn>2</mn> </msup> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>+</mo> <msub> <mi>j&rho;</mi> <mn>0</mn> </msub> <msub> <mi>&omega;q</mi> <mi>&omega;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <mover> <mi>R</mi> <mo>&OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <mo>&part;</mo> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> </mrow> <mrow> <mo>&part;</mo> <mi>n</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <mi>j</mi> <mi>k</mi> </mrow> <mi>&zeta;</mi> </mfrac> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
    根据伽辽金法确定权函数,有
    <mrow> <msub> <mo>&Integral;</mo> <mi>&Omega;</mi> </msub> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>&CenterDot;</mo> <mrow> <mo>(</mo> <msup> <mo>&dtri;</mo> <mn>2</mn> </msup> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>+</mo> <msup> <mi>k</mi> <mn>2</mn> </msup> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>+</mo> <msub> <mi>j&rho;</mi> <mn>0</mn> </msub> <msub> <mi>&omega;q</mi> <mi>&omega;</mi> </msub> <mo>)</mo> </mrow> <mi>d</mi> <mi>v</mi> <mo>-</mo> <msub> <mo>&Integral;</mo> <mi>&Gamma;</mi> </msub> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>&CenterDot;</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&part;</mo> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> </mrow> <mrow> <mo>&part;</mo> <mi>n</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <mi>j</mi> <mi>k</mi> </mrow> <mi>&zeta;</mi> </mfrac> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mi>d</mi> <mi>s</mi> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
    由格林第一公式

    式(13)简化为
    <mrow> <msub> <mo>&Integral;</mo> <mi>&Omega;</mi> </msub> <mrow> <mo>(</mo> <mo>&dtri;</mo> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>&CenterDot;</mo> <mo>&dtri;</mo> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>-</mo> <msup> <mi>k</mi> <mn>2</mn> </msup> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>&CenterDot;</mo> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>-</mo> <msub> <mi>j&rho;</mi> <mn>0</mn> </msub> <mi>&omega;</mi> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <msub> <mi>q</mi> <mi>&omega;</mi> </msub> <mo>)</mo> </mrow> <mi>d</mi> <mi>v</mi> <mo>+</mo> <msub> <mo>&Integral;</mo> <mi>&Gamma;</mi> </msub> <mfrac> <mrow> <mi>j</mi> <mi>k</mi> </mrow> <mi>&zeta;</mi> </mfrac> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>&CenterDot;</mo> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mi>d</mi> <mi>s</mi> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
    在声场中任意一点的声压用各节点声压来表示,即
    <mrow> <mover> <mi>p</mi> <mo>&OverBar;</mo> </mover> <mo>=</mo> <mi>N</mi> <mi>p</mi> <mo>=</mo> <mo>&lsqb;</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>N</mi> <mn>2</mn> </msub> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>N</mi> <mi>n</mi> </msub> <mo>&rsqb;</mo> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <msub> <mi>p</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>p</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> </mtr> <mtr> <mtd> <msub> <mi>p</mi> <mi>n</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
    式中,Ni为节点i处的形函数,pi为节点i处的声压;
    将式(16)代入式(15),得
    <mrow> <msub> <mo>&Integral;</mo> <mi>&Omega;</mi> </msub> <mo>&lsqb;</mo> <msup> <mi>p</mi> <mi>T</mi> </msup> <msup> <mrow> <mo>(</mo> <mo>&dtri;</mo> <mi>N</mi> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mrow> <mo>(</mo> <mo>&dtri;</mo> <mi>N</mi> <mo>)</mo> </mrow> <mi>p</mi> <mo>-</mo> <msup> <mi>k</mi> <mn>2</mn> </msup> <msup> <mi>p</mi> <mi>T</mi> </msup> <msup> <mi>N</mi> <mi>T</mi> </msup> <mi>N</mi> <mi>p</mi> <mo>-</mo> <msub> <mi>j&rho;</mi> <mn>0</mn> </msub> <msup> <mi>&omega;p</mi> <mi>T</mi> </msup> <msup> <mi>N</mi> <mi>T</mi> </msup> <msub> <mi>q</mi> <mi>&omega;</mi> </msub> <mo>&rsqb;</mo> <mi>d</mi> <mi>&Omega;</mi> <mo>+</mo> <msub> <mo>&Integral;</mo> <mi>&Gamma;</mi> </msub> <mfrac> <mrow> <mi>j</mi> <mi>k</mi> </mrow> <mi>&zeta;</mi> </mfrac> <mrow> <mo>(</mo> <msup> <mi>p</mi> <mi>T</mi> </msup> <msup> <mi>N</mi> <mi>T</mi> </msup> <mi>N</mi> <mi>p</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&Gamma;</mi> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>
    式中,为形函数的导数矩阵,其表达式为:
    <mrow> <mo>&dtri;</mo> <mi>N</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>z</mi> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>N</mi> <mn>2</mn> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>N</mi> <mn>2</mn> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>N</mi> <mn>2</mn> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>z</mi> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>N</mi> <mi>n</mi> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>N</mi> <mi>n</mi> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>N</mi> <mi>n</mi> </msub> </mrow> <mrow> <mo>&part;</mo> <mi>z</mi> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <msub> <mi>N</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>N</mi> <mn>2</mn> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>N</mi> <mi>n</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>z</mi> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
    整理式(17),得到
    <mrow> <msub> <mo>&Integral;</mo> <mi>&Omega;</mi> </msub> <mrow> <mo>&lsqb;</mo> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <mo>&dtri;</mo> <mi>N</mi> </mrow> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mrow> <mo>(</mo> <mrow> <mo>&dtri;</mo> <mi>N</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>&rsqb;</mo> </mrow> <mi>d</mi> <mi>v</mi> <mo>&CenterDot;</mo> <mi>p</mi> <mo>-</mo> <msub> <mo>&Integral;</mo> <mi>&Omega;</mi> </msub> <mrow> <mo>(</mo> <mrow> <msup> <mi>k</mi> <mn>2</mn> </msup> <msup> <mi>N</mi> <mi>T</mi> </msup> <mi>N</mi> </mrow> <mo>)</mo> </mrow> <mi>d</mi> <mi>v</mi> <mo>&CenterDot;</mo> <mi>p</mi> <mo>-</mo> <msub> <mo>&Integral;</mo> <mi>&Omega;</mi> </msub> <mrow> <mo>(</mo> <mrow> <msub> <mi>j&rho;</mi> <mn>0</mn> </msub> <msup> <mi>&omega;N</mi> <mi>T</mi> </msup> <msub> <mi>q</mi> <mi>&omega;</mi> </msub> </mrow> <mo>)</mo> </mrow> <mi>d</mi> <mi>v</mi> <mo>+</mo> <msub> <mo>&Integral;</mo> <mi>&Gamma;</mi> </msub> <mfrac> <mrow> <mi>j</mi> <mi>k</mi> </mrow> <mi>&zeta;</mi> </mfrac> <mrow> <mo>(</mo> <mrow> <msup> <mi>N</mi> <mi>T</mi> </msup> <mi>N</mi> </mrow> <mo>)</mo> </mrow> <mi>d</mi> <mi>s</mi> <mo>&CenterDot;</mo> <mi>p</mi> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>

    <mrow> <msub> <mo>&Integral;</mo> <mi>&Omega;</mi> </msub> <msup> <mrow> <mo>(</mo> <mo>&dtri;</mo> <mi>N</mi> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mrow> <mo>(</mo> <mo>&dtri;</mo> <mi>N</mi> <mo>)</mo> </mrow> <mi>d&Omega;</mi> <mo>=</mo> <mi>K</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <mfrac> <mn>1</mn> <mrow> <msup> <msub> <mi>c</mi> <mn>0</mn> </msub> <mn>2</mn> </msup> </mrow> </mfrac> <msub> <mo>&Integral;</mo> <mi>&Omega;</mi> </msub> <msup> <mi>N</mi> <mi>T</mi> </msup> <mi>N</mi> <mi>d</mi> <mi>&Omega;</mi> <mo>=</mo> <mi>M</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <mfrac> <mn>1</mn> <mrow> <msub> <mi>c</mi> <mn>0</mn> </msub> <mi>&zeta;</mi> </mrow> </mfrac> <msub> <mo>&Integral;</mo> <mi>&Gamma;</mi> </msub> <msup> <mi>N</mi> <mi>T</mi> </msup> <mi>N</mi> <mi>d</mi> <mi>&Gamma;</mi> <mo>=</mo> <mi>C</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>22</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <msub> <mo>&Integral;</mo> <mi>&Omega;</mi> </msub> <msub> <mi>j&rho;</mi> <mn>0</mn> </msub> <msup> <mi>&omega;N</mi> <mi>T</mi> </msup> <msub> <mi>q</mi> <mi>&omega;</mi> </msub> <mi>d</mi> <mi>&Omega;</mi> <mo>=</mo> <mi>G</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>23</mn> <mo>)</mo> </mrow> </mrow>
    其中,K称为刚度矩阵,M称为质量矩阵,C称为阻尼矩阵,G称为载荷矩阵;当声源位于
    位置r0(x0,y0,z0)处时,频域内的声源强度表示为:
    qω(r)=qωδ(r-r0) (24)
    其中
    <mrow> <mi>&delta;</mi> <mrow> <mo>(</mo> <mi>r</mi> <mo>-</mo> <msub> <mi>r</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>r</mi> <mo>&NotEqual;</mo> <msub> <mi>r</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mi>r</mi> <mo>=</mo> <msub> <mi>r</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>25</mn> <mo>)</mo> </mrow> </mrow>
    将式(24)代入式(19)中,得
    G=∫Ω-jρ0ωqωδ(r-r0)NTdv=-jρ0ωqωNT (26)
    最后,将式(20)、(21)、(22)、(23)代入式(19)并整理得到
    <mrow> <mfrac> <mrow> <mo>(</mo> <mi>K</mi> <mo>+</mo> <mi>j</mi> <mi>&omega;</mi> <mi>C</mi> <mo>-</mo> <msup> <mi>&omega;</mi> <mn>2</mn> </msup> <mi>M</mi> <mo>)</mo> <msup> <mi>N</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow> <mi>&omega;</mi> </mfrac> <mi>p</mi> <mo>=</mo> <mi>F</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>27</mn> <mo>)</mo> </mrow> </mrow>
    式中,F=jρ0NTqω;K、C、M均为n×n阶的系数矩阵,各自的表达式分别为:
    M=∫ΩNTNdΩ/c02,C=∫ΓNTNdΓ/c0ζ,各式中的N为形函数,在实际
    求解中,K、C、M中的积分运算用求和运算代替
    其中m为积分点的数量,mb为边界上积分点的数量,ξi为积分系数;ω为
    圆频率;c0为空气中的声速;ζ称为比声阻抗,满足ζ=Z/ρ0c0,ρ0为空气密度,Z为界面材料的
    声阻抗;p为封闭空间内任意位置处的声压,实际定位时为单传声器所测得的声信号;F为n
    ×1阶的列向量,表示声源相关信息,其表达式为F=jρ0NTqω,其中qω表示声源强度;
    步骤三、根据节点坐标,利用移动最小二乘法,获得步骤二中所涉及的形函数N;利用
    移动最小二乘法构建形函数;一个场函数u(x)在一点的近似值表示为:
    <mrow> <msup> <mi>u</mi> <mi>h</mi> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>p</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <msub> <mi>a</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>p</mi> <mi>T</mi> </msup> <mrow> <mo>(</mo> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mi>a</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>28</mn> <mo>)</mo> </mrow> </mrow>
    其中是计算点x的邻域范围内各节点的坐标,为
    基函数向量,m为基函数的个数,a(x)=[a1(x),a2(x),…am(x)]为待定系数向量;使用单项
    式基函数做运算,在三维空间中常用的线性及二次单项式基函数分别为:
    <mrow> <mtable> <mtr> <mtd> <mrow> <mi>p</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mrow> <mo>&lsqb;</mo> <mn>1</mn> <mo>,</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>z</mi> <mo>&rsqb;</mo> </mrow> <mi>T</mi> </msup> <mo>,</mo> <mi>m</mi> <mo>=</mo> <mn>4</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>p</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mrow> <mo>&lsqb;</mo> <mrow> <mn>1</mn> <mo>,</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>z</mi> <mo>,</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> <mo>,</mo> <mi>x</mi> <mi>y</mi> <mo>,</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>,</mo> <mi>y</mi> <mi>z</mi> <mo>,</mo> <msup> <mi>z</mi> <mn>2</mn> </msup> <mo>,</mo> <mi>x</mi> <mi>z</mi> </mrow> <mo>&rsqb;</mo> </mrow> <mi>T</mi> </msup> <mo>,</mo> <mi>m</mi> <mo>=</mo> <mn>10</mn> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>29</mn> <mo>)</mo> </mrow> </mrow>
    将求解域用节点离散后,在每个节点处定义一个权函数该函数只在支
    撑域内不为零,在支撑域之外为零,在三维情况下,权函数的支撑域为球形;选定权函数后,
    就求得近似函数在节点处的误差加权平方和:
    <mrow> <mi>J</mi> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>I</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>w</mi> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <msup> <mrow> <mo>&lsqb;</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>p</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mover> <msub> <mi>x</mi> <mi>I</mi> </msub> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <msub> <mi>a</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>u</mi> <mi>I</mi> </msub> <mo>&rsqb;</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>30</mn> <mo>)</mo> </mrow> </mrow>
    令J取最小值,即
    <mrow> <mtable> <mtr> <mtd> <mrow> <mfrac> <mrow> <mo>&part;</mo> <mi>J</mi> </mrow> <mrow> <mo>&part;</mo> <msub> <mi>a</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>=</mo> <mn>2</mn> <munderover> <mi>&Sigma;</mi> <mrow> <mi>I</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>w</mi> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mrow> <mo>&lsqb;</mo> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>p</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mover> <msub> <mi>x</mi> <mi>I</mi> </msub> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <msub> <mi>a</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>u</mi> <mi>I</mi> </msub> </mrow> <mo>&rsqb;</mo> </mrow> <msub> <mi>p</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mover> <msub> <mi>x</mi> <mi>I</mi> </msub> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mi>j</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>m</mi> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>31</mn> <mo>)</mo> </mrow> </mrow>
    经过整理后,得到下式:
    A(x)a(x)=B(x)u (32)
    式中,A(x),B(x)的含义为:
    <mrow> <mtable> <mtr> <mtd> <mrow> <mi>A</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>I</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>w</mi> <mi>I</mi> </msub> <mo>(</mo> <mi>x</mi> <mo>)</mo> <mi>p</mi> <mo>(</mo> <mover> <msub> <mi>x</mi> <mi>I</mi> </msub> <mo>&OverBar;</mo> </mover> <mo>)</mo> <msup> <mi>p</mi> <mi>T</mi> </msup> <mo>(</mo> <mover> <msub> <mi>x</mi> <mi>I</mi> </msub> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>B</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>&lsqb;</mo> <msub> <mi>w</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mi>p</mi> <mo>(</mo> <mover> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>&OverBar;</mo> </mover> <mo>)</mo> <mo>,</mo> <msub> <mi>w</mi> <mn>2</mn> </msub> <mo>(</mo> <mi>x</mi> <mo>)</mo> <mi>p</mi> <mo>(</mo> <mover> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>&OverBar;</mo> </mover> <mo>)</mo> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>w</mi> <mi>N</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mi>p</mi> <mrow> <mo>(</mo> <mover> <msub> <mi>x</mi> <mi>N</mi> </msub> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mo>&rsqb;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>33</mn> <mo>)</mo> </mrow> </mrow>
    由公式(32)得到a(x),将其代入式(28)得:
    <mrow> <msup> <mi>u</mi> <mi>h</mi> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>p</mi> <mi>T</mi> </msup> <mrow> <mo>(</mo> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <msup> <mi>A</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mi>B</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mi>u</mi> <mo>=</mo> <mi>N</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mover> <mi>x</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mi>u</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>34</mn> <mo>)</mo> </mrow> </mrow>
    步骤四、在封闭空间内部任意位置处设置一传声器;当声源有声音发出时,传声器拾
    取到一段音频信号f(t),将其进行傅里叶变换,得到此音频信号的频域信号F(ω);
    步骤五、将步骤四中得到的频域信号F(ω)作为式(29)中的p值代入式(29),求解后得
    到n×1阶的列向量F,由于F中的j、ρ0、qω均为常数,因此,求得的列向量中值最大的元素所
    代表的节点的位置即为声源的位置。

    说明书

    基于封闭空间几何信息建模的单传声器声源定位方法

    技术领域

    本发明涉及一种封闭空间声源定位方法,特别涉及一种基于封闭空间几何信息建
    模的单传声器声源定位方法。

    背景技术

    封闭空间内的声源定位问题是声学研究中的一项重要内容,在人们生活中的很多
    方面都有较广的应用。例如飞机舱室的噪声源定位问题,会议室的智能话筒传音控制系统,
    以及地下车库的车位识别系统等。这类空间的主要特点是环境复杂且混响程度较高,某些
    小尺度的空间还有严重的声波干涉、衍射效应,对声源定位技术的精度、稳定性要求较高。

    文献“基于传声器阵列的声源定位,信息技术,2016(6):136-138”采用了基于声达
    时间差的定位技术,这种声源定位方法一般分为2个步骤进行:先进行声达时间差估计,并
    从中获取传声器阵列中阵元间的声延迟;再利用获取的声达时间差,结合已知的传声器阵
    列的空间位置进一步定出声源的方向。这种方法相对于可控波束形成技术和基于高分辨率
    的谱估计技术的精度较好,计算量较小。

    从算法和实际应用来看,这种声源定位技术还存在着一些缺陷。首先,要获得较高
    的定位精度,需要大规模的传声器阵列,这在经济上会造成很大的负担;对于基于声达时间
    差的定位技术,声源定位实际上只是定位出了声源的方向,具体的距离则无法给出,另外,
    这种方法受环境混响影响很大,当混响程度较高时,精度很低。同时该方法在定位阶段用的
    参数已经是对过去时间的估计,因此不是最优估计。

    发明内容

    为了克服现有封闭空间声源定位方法实用性差的不足,本发明提供一种基于封闭
    空间几何信息建模的单传声器声源定位方法。该方法在封闭空间边界所包围的流体区域内
    布置n个节点,将此n个节点从1到n进行编号。利用此n个节点构建出封闭空间及内部物体的
    几何形状,并根据此n个节点坐标信息建立描述封闭空间数学性质的系统矩阵;在实际声源
    定位时,利用单传声器测得多个频率下的声信号,将此声信号与系统矩阵共同运算,得到声
    信号的位置信息,从而实现声源定位,实用性好。

    本发明解决其技术问题所采用的技术方案:一种基于封闭空间几何信息建模的单
    传声器声源定位方法,其特点是包括以下步骤:

    步骤一、在封闭空间边界所包围的流体区域内布置n个节点,将此n个节点从1到n
    进行编号。n个节点一方面作为封闭空间的几何信息在后期进行建模计算,一方面用来表示
    未来声源定位时的位置。

    n个节点近似均匀分布于整个流体区域内,封闭空间的边界处需布有节点。当封闭
    空间内部存在物体时,物体内部不布置节点,且物体边界上应布置有节点,即节点应勾勒出
    物体形状。节点布置完毕后,封闭空间及其内部物体的边界用Γ表示,边界所包围的流体区
    域用Ω表示。

    步骤二、假定在封闭空间内声源的位置为r,它在单位时间内向单位体积内的空间
    提供了ρ0q(r,t)的媒质质量。根据质量守恒定律,媒质中声波的连续方程写为:


    式中,ρ'为媒质密度增量,ρ0表示媒质静态密度,q为q(r,t)的简写,v为媒质质点
    速度,t表示时间,div为散度算子,在三维空间笛卡尔坐标系中,

    除了连续性方程之外,用来描述媒质声波的基本方程还有两个,它们不受声源的
    影响,分别为运动方程:


    和物态方程:

    p=c02ρ' (3)以上两式中,p代表声压,c0代表声速,grad为
    梯度算子,在三维空间笛卡尔坐标系中,

    由媒质中声波的三个基本方程得到有源情况下封闭空间中有关声压p的波动方
    程:


    式中,为拉普拉斯算子,在三维空间笛卡尔坐标系中,

    在频域,声源强度q(r,t)表示为

    q(r,t)=qω(r)ejωt (5)

    式中,ω为谐振频率,qω(r)为在位置r处频域内的声源强度。

    封闭空间声场各点的声压的频率与声源相同,声压表示为:

    p(r,t)=pω(r)ejωt (6)

    式中,pω(r)为在位置r处频域内的声压。

    将式(6)及式(7)代入式(5)中,得到简谐声源激励下的声波波动方程为:


    令式中k=ω/c0,称其为波数,并且消去ejωt,得到只依赖于空间坐标的那部分方
    程,即室内有源Helmholtz方程:


    这样就将声压的时域问题转换为频域问题,式(9)即为封闭空间声场的控制方程。

    在封闭空间中,边界具有吸声能力,其声压梯度表示为:


    式中,n为封闭空间壁面外法线方向,ζ称为比声阻抗,满足下式:


    式中,Z为界面声阻抗。

    根据Galerkin型加权残量法,为了求解式(8),首先设一试函数为代入有源
    Helmholtz方程及其边界条件,试函数产生残量R和



    根据伽辽金法确定权函数,有


    由格林第一公式


    式(13)简化为


    在声场中任意一点的声压用各节点声压来表示,即


    式中,Ni为节点i处的形函数,pi为节点i处的声压。

    将式(16)代入式(15),得



    式中,为形函数的导数矩阵,其表达式为:


    整理式(17),得到








    其中,K称为刚度矩阵,M称为质量矩阵,C称为阻尼矩阵,G称为载荷矩阵。当声源位
    于位置r0(x0,y0,z0)处时,频域内的声源强度表示为:

    qω(r)=qωδ(r-r0) (24)

    其中


    将式(24)代入式(19)中,得


    最后,将式(20)、(21)、(22)、(23)代入式(19)并整理得到


    式中,F=jρ0NTqω。K、C、M均为n×n阶的系数矩阵,各自的表达式分别为:
    M=∫ΩNTNdΩ/c02,C=∫ΓNTNdΓ/c0ζ,各式中的N为形函数,在实
    际求解中,K、C、M中的积分运算用求和运算代替
    其中m为积分点的数量,mb为边界上积分点的数量,ξi为积分系数;ω为
    圆频率;c0为空气中的声速;ζ称为比声阻抗,满足ζ=Z/ρ0c0,ρ0为空气密度,Z为界面材料的
    声阻抗;p为封闭空间内任意位置处的声压,实际定位时为单传声器所测得的声信号;F为n
    ×1阶的列向量,表示声源相关信息,其表达式为F=jρ0NTqω,其中qω表示声源强度。

    步骤三、根据节点坐标,利用移动最小二乘法,获得步骤二中所涉及的形函数N。利
    用移动最小二乘法构建形函数。一个场函数u(x)在一点的近似值表示为:


    其中是计算点x的邻域范围内各节点的坐标,
    为基函数向量,m为基函数的个数,a(x)=[a1(x),a2(x),…am(x)]为待定系数向量。使用单
    项式基函数做运算,在三维空间中常用的线性及二次单项式基函数分别为:


    将求解域用节点离散后,在每个节点处定义一个权函数该函数只在
    支撑域内不为零,在支撑域之外为零,在三维情况下,权函数的支撑域为球形。选定权函数
    后,就求得近似函数在节点处的误差加权平方和:


    令J取最小值,即


    经过整理后,得到下式:

    A(x)a(x)=B(x)u (32)

    式中,A(x),B(x)的含义为:


    由公式(32)得到a(x),将其代入式(28)得:


    步骤四、在封闭空间内部任意位置处设置一传声器。当声源有声音发出时,传声器
    拾取到一段音频信号f(t),将其进行傅里叶变换,得到此音频信号的频域信号F(ω)。

    步骤五、将步骤四中得到的频域信号F(ω)作为式(29)中的p值代入式(29),求解
    后得到n×1阶的列向量F,由于F中的j、ρ0、qω均为常数,因此,求得的列向量中值最大的元
    素所代表的节点的位置即为声源的位置。

    本发明的有益效果是:该方法在封闭空间边界所包围的流体区域内布置n个节点,
    将此n个节点从1到n进行编号。利用此n个节点构建出封闭空间及内部物体的几何形状,并
    根据此n个节点坐标信息建立描述封闭空间数学性质的系统矩阵;在实际声源定位时,利用
    单传声器测得多个频率下的声信号,将此声信号与系统矩阵共同运算,得到声信号的位置
    信息,从而实现声源定位,实用性好。本发明将几何信息建模引入到封闭空间的声源定位之
    中,仅需单传声器即可完成精度较高的定位。解决了背景技术声源定位方法需要使用传声
    器阵列进行声源定位的技术问题,仅仅利用单传声器测得多个频率下的声信号,将此声信
    号与系统矩阵共同运算,得到声信号的位置信息,从而实现快速实时声源定位。

    下面结合附图和具体实施方式对本发明作详细说明。

    附图说明

    图1是本发明基于封闭空间几何信息建模的单传声器声源定位方法的流程图。

    图2是本发明方法所涉及的封闭空间及传声器、声源示意图。

    图3是本发明方法所涉及的声源点及接收点设置示意图。

    具体实施方式

    参照图1-3。本发明基于封闭空间几何信息建模的单传声器声源定位方法具体步
    骤如下:

    实际封闭空间为一矩形空间,其长宽高分别为0.9m、1.0m、1.2m,其边界为亚克力
    玻璃结构。坐标原点设置在封闭空间的底角处,传声器吊装在封闭空间顶部,其坐标为
    (0.5m,0.5m,1m),声源位置坐标为(0.7m,0.8m,0.3m)。

    步骤1、建立包含封闭空间几何信息的节点模型,由于此空间形状规则,因此节点
    为均匀分布,每个坐标轴向上有9个节点,共有729个节点。

    步骤2、获取计算声源位置的系统方程。

    步骤3、计算形函数N。

    步骤4、在(0.5m,0.5m,1m)处吊装单传声器。当声源有声音发出时,传声器拾取到
    一段音频信号f(t),将其进行傅里叶变换,得到此音频信号的频域转换信号,为F(ω)。

    步骤5、将步骤4中得到的频域信号F(ω)作为式(29)中的p值代入式(29),求解后
    得到n×1阶的列向量F,由于F中的j、ρ0、qω均为常数,列向量F实际由N控制。因此,求得的列
    向量F中值最大的元素所代表的节点的位置即为声源的位置。

    在本实施例中,对10个声源位置分别利用不同的声音作为测试信号进行了定位测
    试,每个声源位置测试10次。最终测试正确率为96%,有效地证明了本发明的有效性。

    更具体的方法步骤如下:

    步骤一、设有一封闭空间,在其内部要实现声源定位。在封闭空间边界所包围的流
    体区域内布置n个节点,将此n个节点从1到n进行编号。n个节点一方面作为封闭空间的几何
    信息在后期进行建模计算,一方面用来表示未来声源定位时的位置。

    节点布置方式没有严格规定,但是应近似均匀分布于整个区域内,封闭空间的边
    界处需布有节点。当封闭空间内部存在物体时,物体内部不布置节点,且物体边界上应布置
    有节点,即节点应勾勒出物体形状。节点布置完毕后,封闭空间及其内部物体的边界用Γ表
    示,边界所包围的流体区域用Ω表示。

    步骤二、获取计算声源位置的系统方程。假定在封闭空间内声源的位置为r,它在
    单位时间内向单位体积内的空间提供了ρ0q(r,t)的媒质质量。根据质量守恒定律,媒质中
    声波的连续方程写为:


    式中ρ'为媒质密度增量,ρ0表示媒质静态密度,q为q(r,t)的简写,v为媒质质点速
    度,t表示时间,div为散度算子,在三维空间笛卡尔坐标系中,

    除了连续性方程之外,用来描述媒质声波的基本方程还有两个,它们不受声源的
    影响,分别为运动方程:


    和物态方程:

    p=c02ρ' (3)

    上两式中,p代表声压,c0代表声速,grad为梯度算子,在三维空间笛卡尔坐标系
    中,

    利用与推导无源波动方程相类似的方法,由媒质中声波的三个基本方程得到有源
    情况下封闭空间中有关声压p的波动方程:


    式中为拉普拉斯算子,在三维空间笛卡尔坐标系中,

    在频域,声源强度q(r,t)表示为

    q(r,t)=qω(r)ejωt (5)

    式中ω为谐振频率,qω(r)为在位置r处频域内的声源强度。

    由于通常将封闭空间声场作为线性系统考虑,因此空间内部各点的声压的频率与
    声源相同,声压表示为:

    p(r,t)=pω(r)ejωt (6)

    式中pω(r)为在位置r处频域内的声压。

    将式(6)及式(7)代入式(5)中,得到简谐声源激励下的声波波动方程为:


    令式中k=ω/c0,称其为波数,并且消去ejωt,得到只依赖于空间坐标的那部分方
    程,即室内有源Helmholtz方程:


    这样就将声压的时域问题转换为频域问题,式(9)即为封闭空间声场的控制方程。

    在封闭空间中,边界通常具有一定的吸声能力,其声压梯度表示为:


    式中n为封闭空间壁面外法线方向,ζ称为比声阻抗,满足下式:


    式中Z为界面声阻抗,查表得到。

    根据Galerkin型加权残量法,为了求解式(8),首先设一试函数为代入有源
    Helmholtz方程及其边界条件,由于试函数通常不是精确解,因此将产生残量R和



    根据伽辽金法确定权函数,有


    由格林第一公式


    式(13)简化为


    在声场中任意一点的声压用各节点声压来表示,即


    式中Ni为节点i处的形函数,pi为节点i处的声压。

    将式(16)代入式(15),得


    式中为形函数的导数矩阵,其表达式为:


    整理式(17),得到






    ∫Ωjρ0ωNTqωdΩ=G (23)

    其中K称为刚度矩阵,M称为质量矩阵,C称为阻尼矩阵,G称为载荷矩阵。当声源位
    于某一特定位置r0(x0,y0,z0)处时,频域内的声源强度表示为:

    qω(r)=qωδ(r-r0) (24)

    其中


    将式(24)代入式(19)中,得

    G=∫Ω-jρ0ωqωδ(r-r0)NTdv=-jρ0ωqωNT (26)

    最后,将式(20)、(21)、(22)、(23)代入式(19)并整理得到


    式中F=jρ0NTqω。

    式中K、C、M均为n×n阶的系数矩阵,各自的表达式分别为:M
    =∫ΩNTNdΩ/c02,C=∫ΓNTNdΓ/c0ζ,各式中的N为形函数,在实际求解中,K、C、M中的积分运算
    用求和运算代替其中m为
    积分点的数量,mb为边界上积分点的数量,ξi为积分系数;ω为圆频率;c0为空气中的声速;ζ
    称为比声阻抗,满足ζ=Z/ρ0c0,ρ0为空气密度,Z为界面材料的声阻抗,查表得到;p为封闭空
    间内任意位置处的声压,实际定位时为单传声器所测得的声信号;F为n×1阶的列向量,表
    示声源相关信息,其表达式为F=jρ0NTqω,式中qω表示声源强度。

    步骤三、计算形函数。根据节点坐标,利用移动最小二乘法,获得步骤2中所涉及的
    形函数N。利用移动最小二乘法构建形函数。一个场函数u(x)在一点的近似值表示为:


    其中是计算点x的邻域范围内各节点的坐标,
    为基函数向量,m为基函数的个数,a(x)=[a1(x),a2(x),…am(x)]为待定系数向量。通常使
    用单项式基函数做运算,在三维空间中常用的线性及二次单项式基函数分别为:


    将求解域用节点离散后,在每个节点处定义一个权函数该函数只在
    一个有限区域(支撑域)内不为零,在区域之外为零,在三维情况下,权函数的支撑域通常为
    球形。常用的权函数有高斯函数、样条函数等。选定权函数后,就求得近似函数在节点处的
    误差加权平方和:


    令J取最小值,即


    经过整理后,得到下式:

    A(x)a(x)=B(x)u (32)

    式中A(x),B(x)的含义为:


    由公式(32)得到a(x),将其代入式(28)得:


    步骤四、在封闭空间内部任意位置处设置一传声器。当声源有声音发出时,传声器
    拾取到一段音频信号f(t),将其进行傅里叶变换,得到此音频信号的频域转换信号,为F
    (ω)。

    步骤五、将步骤四中得到的频域信号F(ω)作为式(29)中的p值代入式(29),求解
    后得到n×1阶的列向量F,由于F中的j、ρ0、qω均为常数,因此,求得的列向量中值最大的元
    素所代表的节点的位置即为声源的位置。

    关于本文
    本文标题:基于封闭空间几何信息建模的单传声器声源定位方法.pdf
    链接地址://www.4mum.com.cn/p-6001404.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 摔了没游戏规则 压单双大小还有时时彩 赛车不定位八码技巧 中国女篮对日女篮比赛 盛兴秒速时时开奖号码 奇妙趋势判断三期内中 11选5利用和值投注技巧 七星彩包头尾奖技巧 实战后二直选45至50注 幸运飞艇不倍投 幸运飞艇单双怎么玩 6码复式三中三论坛资料 快三单双计算方法 黑龙江时时平台网址 赛车pk10定位技巧 后三复试每天稳赚技巧