用于DOA估计的阵列信号处理方法和系统与流程

文档序号:11457834阅读:447来源:国知局
用于DOA估计的阵列信号处理方法和系统与流程

本发明涉及阵列信号处理领域,具体涉及一种用于doa估计的阵列信号处理方法和系统。



背景技术:

波达方向(doa,directionofarrival)估计在过去几十年一直是阵列信号处理领域的重要研究方向,已广泛应用于雷达、通信、侦查等多个领域。最经典的doa估计方法是schmidt提出的music算法,众多研究人员对该算法进行了改进与提高,也提出了一些新的算法。在模型准确的前提下,music算法能高精度的估计出窄带独立信号的波达方向,空间平滑music算法能精确估计出窄带相干信号的波达方向。

但是music算法在计算过程中,需要对空间向量矩阵进行特征值分解,这就导致由c语言实现计算复杂度高,针对这种情况,在专利号为cn200510021874.6的文献中,公开了利用阵列天线并行采集数据,由实数接收矩阵计算空间向量矩阵,通过改进的低复杂度坐标旋转算法实现jacobi特征值分解方法,求解噪声子空间,进而由谱峰搜索估计信源角度。这样的方法在向量空间的特征值分解部分会有较大的运算,其运算量与阵元(传感器)数有关,当阵元数为n时,其运算量约为50n3,在现有的设计中,阵元数大多为8个(通常称为八元阵列)。

在现有的波达方向估计的阵列中,根据实际需要(如待测波源数,精度),会设计出不同数目的阵元的阵列。现有的用于doa估计的阵列信号处理方法或系统,并未对阵元数提出要求或者判断,这就出现了无论阵元数目是多少,均采用同样方法,在一些情况下,就提升了运算量。

因此,开发一种能够根据实际情况,减少运算量的用于doa估计的阵列信号处理方法和系统,在阵列信号处理领域还是非常有市场的。



技术实现要素:

本发明针对现有的music算法,在向量空间的特征分解部分运算量大的技术问题,提供了一种用于doa估计的阵列信号处理方法和系统。

本发明提供的基础方案为:用于doa估计的阵列信号处理方法,包括如下步骤:

接收步骤:利用传感器组成阵列接收目标信号;

处理步骤:对接收到的目标信号进行滤波、放大,然后构造空间向量矩阵;

计数步骤:对用于接收目标信号的传感器进行计数;

求解步骤:根据计数步骤所得的值,选取qr算法或者jacobi算法,求解空间向量矩阵特征值和特征向量;

搜索步骤:谱峰搜索获得定位参数;

显示步骤:将定位结果由显示电路显示。

本发明的工作原理及优点在于:在接收步骤中,传感器组成一个接收传感器阵列,接收目标信号。处理步骤中,能够对接收传感器阵列(或者说传感器)接收到的信号进行滤波、放大处理,然后构造接收矩阵(空间向量矩阵)。在计数步骤中,能够得到用于接收目标信号的传感器的数量,求解步骤中,根据传感器的数量,选取运算量小的算法,达到空间向量矩阵特征值和特征向量,然后通过谱峰搜索获得目标信号的定位参数,最后将定位结果用显示电路显示出来。

本发明用于doa估计的阵列信号处理方法,增加了计数步骤,让后面的求解步骤可以选取运算量小的方法获得空间向量矩阵特征值和特征向量。

进一步,求解步骤中,根据计数步骤所得的值具体为:当计数步骤所获取的值为0-16时,采用qr算法,当计数步骤所获取的值大于17时,采用jacobi算法。

在music类算法中,空间向量矩阵是为hermite矩阵,有ah=a的性质,且其主对角线上元素为实数,特征值均为实数

jacobi算法的思想是通过一系列正交变换不断矫正,减小非对角元素的范数,使其更加对角化。jacobi算法是旋转变换,当阵元数(用于接收目标源信号的传感器数量)为n时,其n阶旋转矩阵形式如下:

即旋转矩阵的对角元素除位置为(p,p)与(q,q)之外其余元素都为1,非对角元素除位置为(p,q)与(q,p)的元素外都为0。c与s为旋转角的余弦与正弦,因而存在关系c2+s2=1。平面旋转则按照下面方式:

显然,旋转后只改变原矩阵的p和q的行与列的值。

下面先介绍c、s、θ值的选择。先考虑2阶方阵,由hermite矩阵性质,对角线元素为实数,故θ2=0,有:

其中带入,则有:

为了使非对角元素零化,令a'p,q=0,设s=tc,由于旋转角的限制,t选择为两个根里绝对值最小的,则得出:

θ1=θ3

这样就完成了c、s、θ值的选择。

qr算法是将矩阵分解为正交矩阵q与上三角矩阵r,并应用位移思想,多次迭代约化hermite矩阵,使其以立方速度收敛到对角型,即:

t=rq=qh(qr)q=qhtq

本申请中,先用householder变换的qr分解将hermite矩阵约化为三对角hermite矩阵,再用隐式位移的givens变换的qr分解进行约化,约化后其仍为三对角矩阵,但几近对角阵。householder变换

形如:

的n阶方阵p记为householder矩阵。矩阵p是正交的,且有p=ph。对于n阶非零列向量x,使p满足:

px=αe0

其中e0=[e,0,...,0]t。代入上式,则有:

因此,有v∈span{x,e0}。令v=x+αe0,带入,为使向量x的系数为0,得α=±||x||2,且e=x(1)/|x(1)|,其中,x(1)为列向量x的第一个元素。这时有:

这样则可以通过连续的householder变换使任意矩阵化为上三角矩阵。

householder变换可使hermite矩阵引入大量零元素,但不能选择性的消去指定元素,而givens变换可以。givens变换与jacobi变换相类似,givens变换的旋转矩阵形式也为:

不同的是,jacobi变换的旋转矩阵jp,q是作用在hermite矩阵的两个位置对称、数值共轭的元素上,故jp,q的对角元素为实数。而givens变换的旋转矩阵gp,q则是作用在hermite矩阵的任意两个元素上,故gp,q的全部元素都为复数。

givens变换是将向量某一指定分量化0。考虑其二阶情况,有:

则有:

θ1=θ3

θ2=θ4

这样则可使二维向量的某一分量化为零。

进一步,处理步骤中,在对接收到的目标信号进行滤波、放大前,先对信号进行幅值调整。幅值调整能够将各路信号调整至相同幅值,能够更好适应后续电路的动态范围,并且调整至相同幅值之后,能够更好的进行滤波处理。

进一步,用于doa估计的阵列信号处理系统,包括:

接收传感器阵列,接收传感器阵列包括计数器和多个传感器,传感器用于接收目标信号,计数器记录传感器的数量;

信号预处理模块,信号处理模块包括增益控制电路和低噪滤波电路,用于对接收信号进行幅值调整和低噪滤波;

算法执行模块,算法执行模块包括多路ad采样电路、处理器及外围电路和存储器,用于将模拟接收信号转换为数字量,处理器用于将数字量转换成空间向量矩阵,处理器用于接收计数器记录的传感器的数量信息,处理器根据传感器的数量对该空间向量矩阵以不同的方法进行特征分解,预估信号源数,构成噪声子空间,然后进行谱峰搜索获得定位参数;

定位结果显示模块,用于显示波达方向信息。

本系统的优势在于:在接收传感器阵列中,加入了计数器以记录传感器的数量。当接收传感器阵列需要检修时,某一传感器被拆卸下后,计数器能够将这个情况反馈到算法执行模块中,这样在执行运算时,能够忽略被拆卸下来的传感器对结果的影响。即,原有的传感器为n时,算法执行模块执行n个传感器的算法,当某一传感器被拆卸下来后,计数器将这一信息反馈到算法执行模块,算法执行模块执行(n-1)个传感器的算法。这样就保证了在接收传感器阵列中,某些传感器失效后doa估计的准确性。

算法执行模块在计数器的反馈下,根据传感器的数量,执行不同的算法,在保证精度的情况下,极大降低了运输量。

进一步,传感器数量与处理器对空间向量矩阵进行特征分解方法的关系为:当传感器数量为0-16时,采用qr算法,当传感器数量大于17时,采用jacobi算法。在运算量上:jacobi算法的运算量约为50n3,qr分解法的运算量则约为3n4,n表示阵元数,因此,算法执行模块可根据传感器的数量,执行不同的算法,在保证精度的情况下,极大了降低了运输量。

进一步,接收传感器阵列中用于接收目标信号的传感器的数目可调。在接收传感器阵列中某一传感器出现问题需要检修时,接收传感器阵列依然可以正常工作。

附图说明

图1是jacobi算法与qr算法误差仿真图;

图2是本发明用于doa估计的阵列信号处理方法的误差仿真图;

图3是本发明用于doa估计的阵列信号处理方法的doa估计的仿真图;

图4是本发明用于doa估计的阵列信号处理方法的流程图。

具体实施方式

下面通过具体实施方式对本发明作进一步详细的说明:

实施例1

基本如附图4所示:用于doa估计的阵列信号处理方法,包括如下步骤:

接收步骤:利用传感器组成阵列接收目标信号;

处理步骤:对接收到的目标信号进行先进行幅值调整,再进行滤波、放大,然后构造空间向量矩阵;

计数步骤:对用于接收目标信号的传感器进行计数;

求解步骤:当计数步骤所获取的值为0-16时,采用qr算法,求解空间向量矩阵特征值和特征向量;当计数步骤所获取的值大于17时,采用jacobi算法,求解空间向量矩阵特征值和特征向量;

搜索步骤:谱峰搜索获得定位参数;

显示步骤:将定位结果由显示电路显示。

对比例1

与实施例1相比,不同之处仅在于,不含计数步骤,在求解步骤中,采用qr算法求解空间向量矩阵特征值和特征向量。

对比例2

与实施例1相比,不同之处仅在于,不含计数步骤,在求解步骤中,采用jacobi算法求解空间向量矩阵特征值和特征向量。

实验过程,对对比例1和对比例2提供的方法进行仿真实验,实验中阵元数范围为[6,25],波长λ=1m,粗搜步长q=1.2°,细搜步长p=0.048°,快拍数n=300,均匀线阵,阵元间距d=0.5m,信噪比为10db,信号源数为三个且信号之间相互独立,信号间角度差为15°,分别为5°、20°、35°,分别实验500次。

实施例1在实验过程中,阵元数为3,即采用的是qr算法。

实验结果如图1所示,从实验结果上来看,jacobi算法与qr算法造成的角度估计的误差较小,误差的数量级仅为10-3,但明显qr算法的误差更小,精度更高。本实施例提供的方法,等同于以阵元数16为分割线,在阵元数为0-16时,采用的是qr算法;当阵元数大于17时,采用的是jacobi算法。实施例1的仿真结果如图2所示。

图3为用本方法作为doa估计的仿真结果。从图3中可以明显的看到三个波峰非常尖锐、明显,可以准确的估计波达方向。

在运算量上:jacobi算法的运算量约为50n3,qr分解法的运算量则约为3n4,n表示阵元数。在不同的阵元数的情况下,达到了既保证精度,也减小运算量的效果。

用于doa估计的阵列信号处理系统,包括:

接收传感器阵列,接收传感器阵列包括计数器和多个可拆卸的传感器,传感器用于接收目标信号,计数器记录传感器的数量;

信号预处理模块,信号处理模块包括增益控制电路和低噪滤波电路,用于对接收信号进行幅值调整和低噪滤波;

算法执行模块,算法执行模块包括多路ad采样电路、处理器及外围电路和存储器,用于将模拟接收信号转换为数字量,处理器用于将数字量转换成空间向量矩阵,处理器用于接收计数器记录的传感器的数量信息,当传感器数量为0-16时,处理器该空间向量矩阵采用qr算法进行特征分解;当传感器数量大于17时,处理器该空间向量矩阵采用jacobi算法进行特征分解,然后预估信号源数,构成噪声子空间,然后进行谱峰搜索获得定位参数;

定位结果显示模块,用于显示波达方向信息。

以上所述的仅是本发明的实施例,方案中公知的具体结构及特性等常识在此未作过多描述。应当指出,对于本领域的技术人员来说,在不脱离本发明结构的前提下,还可以作出若干变形和改进,这些也应该视为本发明的保护范围,这些都不会影响本发明实施的效果和专利的实用性。本申请要求的保护范围应当以其权利要求的内容为准,说明书中的具体实施方式等记载可以用于解释权利要求的内容。

当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1