[发明专利]一种基于粒子滤波的相干分布源波达方向跟踪方法

说明书

技术领域

本发明属于阵列天线信号处理技术领域，特别是涉及一种基于粒子滤波的相干分布源波达方向跟踪方法。

背景技术

波达方向估计是阵列天线信号处理领域的重要研究内容。针对波达估计方法，相关研究人员已经提出了Capon波束形成、最大似然估计、MUSIC谱以及ESPRIT等方法。但是这些方法主要针对点目标源，即信号由一个远场质点目标发出。这种点目标源模型在雷达、声呐以及无线通信等领域并不适用，这些场景中的目标源通常在角度空间上存在一定的扩展，因此利用空间分布源模型作为这类场景中目标的波达方向估计模型更为适合。根据分布源的散射特性，可以将其分为相干分布源和非相干分布源两类。

针对相干分布源的波达方向估计方法，研究人员提出了广义Capon波束形成方法、DPSE方法、广义ESPRIT方法等。但是由于相干分布源参数包括中心波达方向及角分布参数等至少两个参数，所以上述方向均需要进行多维搜索处理，这使得其计算量大大增加。

上述针对点目标源及分布源的波达估计方法只适用于目标位置保持不变的情况，对于目标运动的非平稳信号的波达方向估计则需要对于每一个快拍的阵列天线接收信号进行一次波达方向估计，并对目标源的波达方向进行跟踪。目前这类方法主要针对点目标源，包括投影近似子空间跟踪算法(PAST)、正交投影近似子空间跟踪算法(OPAST)等。

对于分布源的波达方向跟踪方法仅有FAPI-TLS-ESPRIT方法，该方法利用FAPI算法实现噪声子空间在每个快拍的更新，再利用TLS-ESPRIT算法实现每个快拍的分布源的波达方向估计。该方法仅能实现中心波达方向的跟踪，不能进行角分布参数的跟踪。

粒子滤波算法是一种贝叶斯递推估计方法，能够实现较强噪声条件下的不确定参数准确估计。但到目前为止尚未发现利用粒子滤波算法实现相干分布源波达方向跟踪方面的报道。

发明内容

为了解决上述问题，本发明的目的在于提供一种基于粒子滤波的相干分布源波达方向跟踪方法。

为了达到上述目的，本发明提供的基于粒子滤波的相干分布源波达方向跟踪方法包括按顺序进行的下列步骤：

2)建立相干分布源的广义阵列流型及阵列天线接收信号模型的S1阶段；

2)基于上述相干分布源的广义阵列流型及阵列天线接收信号模型，建立基于相干分布源参数谱的观测似然函数的S2阶段；

3)构造相干分布源参数状态方程的S3阶段；

4)基于步骤2)获得的观测似然函数和步骤3)获得的相干分布源参数的状态方程，利用粒子滤波算法实现相干分布源波达方向跟踪的S4阶段。

在步骤1)中，所述的建立相干分布源的广义阵列流型及阵列天线接收信号模型的方法是：在等距线阵情况下，首先将分布源的阵列天线接收信号模型描述为角信号密度函数在分布空间积分的形式；然后对于相干分布源模型，则角信号密度函数可以表示为随机信号幅度与分布源的空间分布函数乘积的形式；假设相干分布源的空间分布函数为高斯分布，则可以建立起相干分布源的广义阵列流型及阵列天线接收信号模型。

在步骤2)中，所述的基于上述相干分布源的广义阵列流型及阵列天线接收信号模型，建立基于相干分布源参数谱的观测似然函数的方法是：在步骤1)建立的相干分布源的广义阵列流型及阵列天线接收信号模型的基础上，利用DSPE算法估计出接收天线信号的协方差矩阵，并对接收天线信号的协方差矩阵进行特征值分解，选取小特征值对应的特征向量构建噪声子空间，利用阵列天线导向矢量与噪声子空间估计相干分布源参数谱，并以该相干分布源参数谱作为观测似然函数。

在步骤3)中，所述的构造相干分布源参数状态方程的方法是：以相干分布源的波达方向、波达方向的速度、角分布参数、角分布速度构成相干分布源的状态向量，然后利用匀速运动模型和上述相干分布源的状态向量建立相干分布源参数的状态方程。

在步骤4)中，所述的基于步骤2)获得的观测似然函数和步骤3)获得的相干分布源参数的状态方程，利用粒子滤波算法实现相干分布源波达方向跟踪的方法是：使用粒子滤波算法，利用波达方向变化的时间相关性，结合每个快拍的阵列天线接收信号，实现波达方向跟踪。

本发明提供的基于粒子滤波的相干分布源波达方向跟踪方法首先建立相干分布源的阵列天线接收信号模型，然后在此模型基础上利用DSPE算法建立相干分布源参数谱的观测似然函数，并构造相干分布源参数状态方程，进而使用粒子滤波算法实现相干分布源的波达方向跟踪。本发明方法具有不需要搜索处理，同时估计中心波达方向及角分布参数，估计精度高、估计性能好等优点。