一种双五元十字阵被动声定位融合方法技术

一种双五元十字阵被动声定位融合方法，属于声音信号处理及测量技术领域。本发明专利技术解决了现有的四元十字阵的定位方法定向和定距精度和目标的方位有关，无法判断目标在哪个空域；现有五元十字阵测距和侧向误差依然存在，尤其是当目标与平面阵的夹角较小时，测距和侧向误差变大，以及在推导公式时忽略了一些微小距离量也引入了一定的误差，估计的精度下降，不利于实际的工程应用的问题。本发明专利技术用七个传声器构造了双五元十字阵传声器阵列，利用两个五元十字阵分别推导出声源点准确坐标计算公式，然后利用两个五元十字阵推出的计算公式引入加权算数平均法来构造声源点坐标值的计算公式。本发明专利技术适用于基于时延估计的被动声声源定位领域。

【技术实现步骤摘要】
一种双五元十字阵被动声定位融合方法
本专利技术涉及五元十字阵被动声定位融合方法，尤其涉及一种双五元十字阵被动声定位融合方法，属于声音信号处理及测量

技术介绍
在被动声定位方法中，波束形成方法或时延估计方法对声目标定位性能的影响都和传声器阵列形式相关。阵列设计方式的不同，得到的点声源目标定位性能也不一样。传声器阵列设计可分为线型阵列，平面阵列和立体阵列等。线型阵列可确定目标的二维参量，平面阵列能确定目标的三维参量，而且在时延估计精度相同的情况下，对于不同的传声器阵型，所得到的点声源目标定位性能是不同的。四元十字阵的定位方法的缺点是：定向和定距精度和目标的方位有关，四元十字阵定位所解得的目标坐标是以阵平面为分界的上、下半球内的两个解,无法判断目标在哪个空域；现有五元十字阵虽然具备全空域声源定位能力,但与平面四元阵相比,定位精度并没有得到实际的提高，测距和侧向误差依然存在，尤其是当目标与平面阵的夹角较小时,测距和侧向误差变大，以及大家在推导公式时忽略了一些微小距离量也引入了一定的误差，估计的精度下降，不利于实际的工程应用。现有各种时延估计方法包括GCC方法、LMS方法、倒谱预滤波技术、基于EVD以及基于ATF比的时延估计方法等，不管哪一种方法在估计时延值时都会产生时延估计误差，这就使得利用平面十字阵推导的被动声源点计算公式计算出的坐标值不够准确。
技术实现思路
本专利技术的目的是提出一种双五元十字阵被动声定位融合方法，以解决现有的四元十字阵的定位方法定向和定距精度和目标的方位有关，无法判断目标在哪个空域；现有五元十字阵测距和侧向误差依然存在，尤其是当目标与平面阵的夹角较小时,测距和侧向误差变大，以及在推导公式时忽略了一些微小距离量也引入了一定的误差，估计的精度下降，不利于实际的工程应用的问题。本专利技术所述的一种双五元十字阵被动声定位融合方法，所述双五元十字阵被动声的传声器阵列模型由7个微传声器组成，S0、S1、S2、S3、S4、S5、S6分别表示所述的7个传声器，位于笛卡尔坐标系XYZ坐标轴上；设S0(0,0,0)为基准传声器，其余六个传声器位于XOY和XOZ两个平面上，设6个传声器与坐标原点O的距离都是D；声源P在球坐标系下的坐标为与坐标原点之间的距离PO为r，俯仰角POZ为θ，方位角P0OX为c为声波传播速度；双五元十字阵都以S0为基准传声器，获得6个相互独立的时间延迟值为τ10、τ20、τ30、τ40、τ50、τ60；所述定位融合方法包括以下步骤：步骤一、利用位于笛卡尔坐标系XOY平面内第一个五元十字阵，推导出该五元十字阵准确的声源点坐标计算公式；步骤二、利用位于笛卡尔坐标系XOZ平面内第二个五元十字阵，推导出该五元十字阵准确的声源点坐标计算公式；步骤三、以步骤一和步骤二推导出的准确声源点坐标计算公式，对上述两个步骤得到的声源点坐标计算公式用加权算术平均法来加权处理，构造出双五元十字阵的声源点坐标计算公式，声源P在第一个五元十字阵中以P0表示，在第二个五元十字阵中以P1表示。本专利技术的有益效果是：1、本专利技术能够克服单一平面四元十字阵和平面五元十字阵在定位过程中测距和侧向的精度都与目标的方位有关，尤其是当目标与平面阵的夹角较小时,测距和侧向误差变大的缺点。(当目标与平面阵XOY的夹角较小时,目标与垂直阵XOZ的夹角较大，这样通过调整加权系数就可以准确测量到声源点的位置信息)2、本专利技术克服了现有方法中，在推导公式时忽略了一些微小距离量引入的误差，估计的精度下降以及不利于实际的工程应用的缺点。3、本专利技术专利设计的双五元十字阵传声器阵列具有全空域声源定位的能力,在时延估计误差一定的情况下能够获得最准确的声源点的位置信息。4、本专利技术方法不仅可以消除目标方位变化时对测距测向精度的影响，还可以抑制时延估计误差对测距测向性能的影响，能够更加准确的确定出声源点的位置信息。5、本专利技术性能优良、结构简单、仿真过程中理想情况下计算得到结果没有误差；当时延估计引入0.1％的误差时，若用本专利技术仿真计算结果，声源点与坐标原点的距离误差为0.1％，方位角误差为0.41％，俯仰角误差为0.0255％；若用单一五元十字阵仿真计算结果，声源点与坐标原点的距离误差为0.1％，方位角误差为0.6684％，俯仰角误差为-0.1273％，可见本专利技术的计算结果更加精确。附图说明图1为本专利技术的双五元十字被动声定位融合方法阵型，其中S0、S1、S2、S3、S4、S5、S6分别为传声器，位于笛卡尔坐标系XYZ坐标轴上，P为声源点，P’为声源点P投影到XOY平面的投影点；图2为本专利技术的双五元十字被动声定位融合方法第一个五元十字阵型，其中S0、S1、S2、S3、S4、S5、S6分别为传声器，位于笛卡尔坐标系XYZ坐标轴上，P为声源点，P’为声源点P投影到XOY平面的投影点；图3为本专利技术的双五元十字被动声定位融合方法第二个五元十字阵型，其中S0、S1、S2、S3、S4、S5、S6分别为传声器，位于笛卡尔坐标系XYZ坐标轴上，P为声源点，P’为声源点P投影到XOY平面的投影点；图4为本专利技术的双五元十字被动声定位融合方法流程图。具体实施方式结合图1、图2、图3和图4说明本专利技术的具体实施方式。具体实施方式一：本专利技术所述的一种双五元十字阵被动声定位融合方法，所述双五元十字阵被动声的传声器阵列模型由7个微传声器组成，S0、S1、S2、S3、S4、S5、S6分别表示所述的7个传声器，位于笛卡尔坐标系XYZ坐标轴上，用来接收声源发出的声信号；设S0(0,0,0)为基准传声器，其余六个传声器位于XOY和XOZ两个平面上，设6个传声器与坐标原点O的距离都是D，如图1，则各个传声器的坐标可表示为：S1(D,0,0)、S2(0,-D,0)、S3(0,0,D)、S4(-D,0,0)、S5(0,D,0)、S6(0,0,-D)；设声源P，笛卡尔坐标系下的坐标为(x,y,z)，声源P在球坐标系下的坐标为与坐标原点之间的距离PO为r，俯仰角POZ为θ，方位角P0OX为声源P产生的声波以球面波形式传播，c为声波传播速度；双五元十字阵都以S0为基准传声器，获得6个相互独立的时间延迟值为τ10、τ20、τ30、τ40、τ50、τ60；6个时延值之间是相互独立的。所述定位融合方法包括以下步骤：步骤一、利用位于笛卡尔坐标系XOY平面内第一个五元十字阵，推导出该五元十字阵准确的声源点坐标计算公式；步骤二、利用位于笛卡尔坐标系XOZ平面内第二个五元十字阵，推导出该五元十字阵准确的声源点坐标计算公式；步骤三、以步骤一和步骤二推导出的准确声源点坐标计算公式，对上述两个步骤得到的声源点坐标计算公式用加权算术平均法来加权处理，构造出双五元十字阵的声源点坐标计算公式，声源P在第一个五元十字阵中以P0表示，在第二个五元十字阵中以P1表示。所述的声源点P的坐标在笛卡尔坐标系中P(x，y，z)与球坐标系中关系为：其中：0°≤θ≤180°，具体实施方式二：所述的双五元十字阵分别位于笛卡尔坐标系XOY和XOZ两个平面上，其中在笛卡尔坐标系中以XOY平面内的四个传声器S1、S2、S4、S5与S0构成第一个五元十字阵，在笛卡尔坐标系中以XOZ平面内的四个传声器S1、S3、S4、S6与S0构成第二个五元十字阵。具体实施方式本文档来自技高网...

【技术保护点】
一种双五元十字阵被动声定位融合方法，其特征在于所述双五元十字阵被动声的传声器阵列模型由7个微传声器组成，S0、S1、S2、S3、S4、S5、S6分别表示所述的7个传声器，位于笛卡尔坐标系XYZ坐标轴上；设S0(0,0,0)为基准传声器，其余六个传声器位于XOY和XOZ两个平面上，设6个传声器与坐标原点O的距离都是D；声源P在球坐标系下的坐标为与坐标原点之间的距离PO为r，俯仰角POZ为θ，方位角P0OX为c为声波传播速度；双五元十字阵都以S0为基准传声器，获得6个相互独立的时间延迟值为τ10、τ20、τ30、τ40、τ50、τ60；所述定位融合方法包括以下步骤：步骤一、利用位于笛卡尔坐标系XOY平面内第一个五元十字阵，推导出该五元十字阵准确的声源点坐标计算公式；步骤二、利用位于笛卡尔坐标系XOZ平面内第二个五元十字阵，推导出该五元十字阵准确的声源点坐标计算公式；步骤三、以步骤一和步骤二推导出的准确声源点坐标计算公式，对上述两个步骤得到的声源点坐标计算公式用加权算术平均法来加权处理，构造出双五元十字阵的声源点坐标计算公式，声源P在第一个五元十字阵中以P0表示，在第二个五元十字阵中以P1表示。...

【技术特征摘要】
1.一种双五元十字阵被动声定位融合方法，其特征在于所述双五元十字阵被动声的传声器阵列模型由7个微传声器组成，S0、S1、S2、S3、S4、S5、S6分别表示所述的7个传声器，位于笛卡尔坐标系XYZ坐标轴上；设S0(0,0,0)为基准传声器，其余六个传声器位于XOY和XOZ两个平面上，设6个传声器与坐标原点O的距离都是D；声源P在球坐标系下的坐标为与坐标原点之间的距离PO为r，俯仰角POZ为θ，方位角PaOX为c为声波传播速度；双五元十字阵都以S0为基准传声器，获得6个相互独立的时间延迟值为τ10、τ20、τ30、τ40、τ50、τ60；所述定位融合方法包括以下步骤：步骤一、利用位于笛卡尔坐标系XOY平面内第一个五元十字阵，推导出该五元十字阵准确的声源点坐标计算公式；步骤二、利用位于笛卡尔坐标系XOZ平面内第二个五元十字阵，推导出该五元十字阵准确的声源点坐标计算公式；步骤三、以步骤一和步骤二推导出的准确声源点坐标计算公式，对上述两个步骤得到的声源点坐标计算公式用加权算术平均法来加权处理，构造出双五元十字阵的声源点坐标计算公式，声源P在第一个五元十字阵中以P0表示，在第二个五元十字阵中以P1表示；所述的双五元十字阵分别位于笛卡尔坐标系XOY和XOZ两个平面上，其中在笛卡尔坐标系中以XOY平面内的四个传声器S1、S2、S4、S5与S0构成第一个五元十字阵，在笛卡尔坐标系中以XOZ平面内的四个传声器S1、S3、S4、S6与S0构成第二个五元十字阵；步骤一所述的位于XOY平面的第一个五元十字阵推导准确的声源点坐标计算公式的具体过程为：假设在XOY平面的第一个五元...

【专利技术属性】
技术研发人员：叶东，路敬祎，陈刚，郭玉波，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江;23