法律状态公告日
法律状态信息
法律状态
2023-08-22
专利权的转移 IPC(主分类):G01S 3/80 专利号:ZL201710171734X 登记生效日:20230807 变更事项:专利权人 变更前权利人:陕西理工大学 变更后权利人:深圳万知达科技有限公司 变更事项:地址 变更前权利人:723001 陕西省汉中市汉台区朝阳路东关正街505号 变更后权利人:518000 广东省深圳市龙华区大浪街道横朗社区华兴路13号智云产业园A栋1002
专利申请权、专利权的转移
2020-06-19
著录事项变更 IPC(主分类):G01S3/80 变更前: 变更后: 申请日:20170321
著录事项变更
2019-10-18
授权
授权
2017-08-25
实质审查的生效 IPC(主分类):G01S3/80 申请日:20170321
实质审查的生效
2017-08-01
公开
公开
技术领域
本发明属于信号处理技术领域,尤其涉及一种声矢量传感器阵列的近场源频率、二维到达角和距离的估计方法。
背景技术
确定声波的到达方向是声信号处理的一个重要应用,早期的声波测向采用的是仅能测量声压强度信息的声压传感器,随着矢量传感器技术的发展,利用声矢量传感器进行波达方向估计受到国内外诸多学者的关注。不同于无指向性的标量声压传感器,声矢量传感器由一个声压传感器和可选的三个空间轴向垂直的质点振速传感器构成,正是由于这种特殊结构,声矢量传感器具备空间同步、共点测量声场中的标量声压信息和矢量振速信息的能力,因而能获得比传统声压传感器更多的信息量,从而使得声矢量传感器阵列能够获得比声压传感器阵列更好的定位精度和空间分辨率。
传统的声矢量传感器信号处理采用的是长矢量信号模型,该模型下的每一个声矢量传感器输出数据是用一个复数矢量来描述,声矢量传感器阵列中的多个输出数据顺序排列构成复数数据矢量。基于长矢量模型的数据处理并没有利用声矢量传感器各个分量间的正交关系,本发明利用声矢量传感器自身具有的正交矢量特性,构建了由x轴、y轴和z轴三个振速传感器分量和一个声压传感器分量组成的声矢量传感器的四元数新模型。因为近场的波面是球面波。阵元间的相位不仅与阵元间距有关而且与声源到阵元的距离有关,所以远场条件下具有的平移不变结构对于近场不再适用,标量传感器阵列无法利用旋转不变技术估计信号参数(ESPRIT)方法估计近场声源信号的参数,现有的文献亦未涉及声矢量传感器阵列的近场源ESPRIT参数估计方法,本发明利用声矢量传感器自身具有的旋转不变结构进行参数估计,进而给出了一种近场的声矢量传感器阵列参数估计ESPRIT方法。
发明内容
本发明的目的是提供一种近场源任意空间结构的声矢量传感器阵列ESPRIT多参数估计方法。
为了实现上述目的,本发明采取如下的技术解决方案:
任意声矢量传感器阵列近场源ESPRIT参数估计方法,K个不同频率、互不相关窄带、随机平稳近场声源信号从不同的方向和距离(θk,φk,rk)入射到该接收阵列,所述阵列由M个空间任意分布的阵元构成,所述阵元是具有同步共点测量声压以及x轴、y轴和z轴方向振速分量的声矢量传感器,所有传感器的对应通道相互平行:所有的声压传感器相互平行,所有的x轴方向振速传感器相互平行,所有的y轴方向振速传感器相互平行,以及所有的z轴方向振速传感器相互平行;相邻阵元间隔与入射声波信号的波长和声源的距离之间满足近场条件;
任意声矢量传感器阵列近场源ESPRIT参数估计方法的步骤如下:
步骤一、M个空间任意分布的阵元构成空间任意阵接收近场信号,利用排布好的接收阵列获取数据,所有阵元同步采样一次的输出数据称为一次快拍数据,N次快拍数据构成阵列接收数据Z;
步骤二、计算数据相关矩阵,由数据相关矩阵获取信号子空间Us;
计算数据相关矩阵
步骤三、通过信号子空间Us的分块运算获取重排后子阵间的旋转不变关系矩阵
根据数据排列方式将信号子空间Us分成M×K的U1,U2,U3和U4四个子空间,即Us=[U1,U2,U3,U4]T,其中Us=AT,U1=A1T,U2=A2T,U3=A3T,U4=A4T,A就是步骤二中的阵列导向矢量矩阵,A1是x轴方向振速分量采样数据阵列导向矢量矩阵,A2是y轴方向振速分量采样数据阵列导向矢量矩阵,A3是z轴方向振速分量采样数据阵列导向矢量矩阵,A4是声压分量采样数据阵列导向矢量矩阵,T是阵列导向矢量和信号子空间之间的K×K的非奇异变换矩阵,可以得到Ψ1T1=Ω1T1,Ψ2T2=Ω2T2,Ψ3T3=Ω3T3,其中
其中,
步骤四、利用旋转不变关系矩阵
其中,
前述步骤中的k=1,...,K,l=1,...,K,j表示虚数单位。
本发明采用的任意阵型的阵列,阵列的阵元为由声压传感器和x轴、y轴及z轴方向的振速传感器构成的声矢量传感器,并且所有的声压传感器相互平行,所有的x轴方向振速传感器相互平行,所有的y轴方向振速传感器相互平行,所有的z轴方向振速传感器相互平行。
本发明利用声矢量传感器自身具有的正交矢量特性,给出了声矢量传感器二维到达角和距离的近场ESPRIT估计方法,突破了现有线阵近场源参数估计方法只能估计一维到达角的局限性,解决了近场条件下声矢量传感器阵列不具有的平移不变结构,无法利用ESPRIT方法的难题,该发明方法不限定阵列阵型,根据同一个列矢量与自身转置复共轭相乘是最大的原则进行参数配对。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中需要使用的附图做简单介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例声矢量传感器阵列的示意图;
图2为本发明方法的流程图;
图3为本发明方法的到达角估计散布图;
图4为本发明方法的俯仰角估计均方根误差随信噪比的变化曲线图;
图5为本发明方法的方位角估计均方根误差随信噪比的变化曲线图;
图6为本发明方法的到达角估计均方根误差随信噪比的变化曲线图;
图7为本发明方法的距离估计均方根误差随信噪比的变化曲线图;
图8为本发明方法的到达角估计成功概率随信噪比的变化曲线图。
具体实施方式
为了让本发明的上述和其它目的、特征及优点能更明显,下文特举本发明实施例,并配合所附图示,做详细说明如下。
图1所示为本发明实施例的声矢量传感器阵列的示意图。本发明的声矢量传感器阵列由M个空间任意分布的阵元构成,所述阵元是具有同步共点测量声压以及x轴、y轴和z轴方向振速分量的声矢量传感器,阵元间隔与入射声波信号的波长和声源的距离之间满足近场条件;
参照图2,本发明的近场源ESPRIT参数估计方法的步骤如下:任意阵型声矢量传感器阵列接收K个不同频率、互不相关窄带、随机平稳近场声源信号,K为入射声源信号的数量,
步骤一、M个空间任意分布的阵元构成空间任意阵接收近场信号,利用排布好的接收阵列获取数据,所有阵元同步采样一次的输出数据称为一次快拍数据,N次快拍数据构成阵列接收数据Z;
步骤二、计算数据相关矩阵,由数据相关矩阵获取信号子空间Us;
计算数据相关矩阵
步骤三、通过信号子空间Us的分块运算获取重排后子阵间的旋转不变关系矩阵
根据数据排列方式将信号子空间Us分成M×K的U1,U2,U3和U4四个子空间,即Us=[U1,U2,U3,U4]T,其中Us=AT,U1=A1T,U2=A2T,U3=A3T,U4=A4T,A就是步骤二中的阵列导向矢量矩阵,A1是x轴方向振速分量采样数据阵列导向矢量矩阵,A2是y轴方向振速分量采样数据阵列导向矢量矩阵,A3是z轴方向振速分量采样数据阵列导向矢量矩阵,A4是声压分量采样数据阵列导向矢量矩阵,T是阵列导向矢量和信号子空间之间的K×K的非奇异变换矩阵,可以得到Ψ1T1=Ω1T1,Ψ2T2=Ω2T2,Ψ3T3=Ω3T3,其中
其中,
步骤四、利用旋转不变关系矩阵
其中,
前述步骤中的k=1,...,K,l=1,...,K,j表示虚数单位;
本发明给出了基于任意阵型的近场源参数估计的ESPRIT方法,根据空间上共点的声矢量传感器的各个子阵对应的空域导向矢量相同的特点,利用同一个列矢量与自身转置复共轭相乘是最大的原则,对特征矢量和特征值重排,从而获取子阵间的旋转不变关系矩阵,利用旋转不变关系矩阵得到信号的到达角和声源到坐标原点的距离估计值,该方法突破了现有线阵参数估计方法的局限性以及多重信号分类(MUSIC)方法的全域三维搜索的超大计算量难题,参数配对运算简单;
本发明的效果可以通过以下的仿真结果进一步说明:
仿真实验条件如下:
不失一般性,假设两个不同频率、互不相关窄带、随机平稳近场声源信号入射到由9个阵元构成的空间任意排布的声矢量传感器阵列,各阵元坐标分别为(0.2λmin,0.2λmin,0),(0.3λmin,0.2λmin,0),(0.4λmin,0.3λmin,0),(0.2λmin,0,0.2λmin),(0.2λmin,0,0.4λmin),(0.3λmin,0,0.5λmin),(0,0.4λmin,0.4λmin),(0,0.5λmin,0.5λmin),(0,0.7λmin,0.2λmin),λmin为入射声波信号的最小波长,如图1所示,入射信号的参数为:(θ1,φ1)=(20°,24°),(θ2,φ2)=(50°,27°),其归一化频率为(f1,f2)=(0.3,0.4),快拍数为512次,200次独立实验。
仿真实验结果如图3至图8所示,图3为信噪比是15dB时,本发明方法到达角估计的散布图,从图3可以看出本发明方法可以估计出到达角参数,且本发明方法有较高的到达角参数估计精度;从图4和图7可以看出本发明方法,俯仰角、方位角、到达角和距离估计的均方根误差较小,也就是估计值在真值附近的较小范围内扰动;到达角估计成功概率是指在200次独立试验中俯仰角和方位角估计值满足关系式
以上所述,仅是本发明的较佳实施例而已,并非对本发明做任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
机译: 声矢量传感器阵列的定位和校准
机译: 声矢量传感器阵列的定位和校准
机译: 声矢量传感器阵列的定位和校准