一种基于粒子滤波的相干分布源波达方向跟踪方法技术

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

A particle filter based coherent source direction of arrival tracking method

The invention discloses a method for tracking the direction of arrival of coherent distributed sources based on particle filtering. This method first establishes the antenna array of coherently distributed source model of received signals, and then based on this model to establish the likelihood function observation of coherently distributed source spectrum parameters using DSPE algorithm, and construct coherent distributed source parameters of state equation, and then use the DOA tracking particle filter algorithm to achieve the phase distribution of dry source. The method of the invention has the advantages of no search processing, simultaneous estimation of the central direction of arrival and angular distribution parameters, high estimation accuracy and good estimation performance, etc..

【技术实现步骤摘要】
一种基于粒子滤波的相干分布源波达方向跟踪方法
本专利技术属于阵列天线信号处理
，特别是涉及一种基于粒子滤波的相干分布源波达方向跟踪方法。
技术介绍
波达方向估计是阵列天线信号处理领域的重要研究内容。针对波达估计方法，相关研究人员已经提出了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算法建立相干分布源参数谱的观测似然函数，并构造相干分布源参数状态方程，进而使用粒子滤波算法实现相干分布源的波达方向跟踪。本专利技术方法具有不需要搜索处理，同时估计中心波达方向及角分布参数，估计精度高、估计性能好等优点。附图说明图1为本专利技术提供的基于粒子滤波的相干分布源波达方向跟踪方法流程图。图2为SNR＝3dB时本专利技术方法获得的相干分布源中心波达方向估计值。图3为SNR＝3dB时本专利技术方法获得的相干分布源中心波达方向估计值的RMSE。图4SNR＝3dB时本专利技术方法获得的相干分布源角分布参数估计值。图5为相干分布源中心波达方向估计值随信噪比变化的平均RMSE。具体实施方式下面结合附图和具体实施例对本专利技术提供的基于粒子滤波的相干分布源波达方向跟踪方法进行详细说明。图1为本专利技术提供的基于粒子滤波的相干分布源波达方向跟踪方法流程图。其中的全部操作都是在计算机系统中完成的，操作的主体均为计算机系统。如图1所示，本专利技术提供的基于粒子滤波的相干分布源波达方向跟踪方法包括按顺序进行的下列步骤：3)建立相干分布源的广义阵列流型及阵列天线接收信号模型的S1阶段：本阶段是在等距线阵情况下，利用分布源阵列天线接收信号模型，结合高斯角分布函数，建立相干分布源的广义阵列流型及阵列天线接收信号模型，然后进入下一步S2阶段。在此阶段中，首先将如式(1)所示的分布源的阵列天线接收信号模型描述为角信号密度函数在分布空间积分的形式，即：其中，y(t)为分布源的阵列天线接收信号，且y(t)＝[y1(t),y2(t),…,yL(t)]T；a(θ)为阵列天线导向矢量；si(θ-θi,t)为第i个分布源的角信号密度函数，nt为观测噪声。对于相干分布源模型，则角信号密度函数可以表示为随机信号幅度与分布源的空间分布函数乘积的形式，即：si(θ-θi,t)＝si(t)gi(θ-θi)(2)其中，si(t)为随机信号幅度；gi(θ-θi)为分布源的空间分布函数。将公式(2)带入公式(1)，并假设相干分布源的空间分布函数为高斯分布，即：其中，Δi为未知角分布参数，则有：其中，将公式(3)带入公式(5)则可以建立相干分布源的广义阵列流型及阵列天线接收信号模型，即：yt＝Btst+nt(6)其中，Bt＝[b1(θ1),…,bq(θq)]，且2)基于上述相干分布源的广义阵列流型及阵列天线接收信号模型，建立基于相干分布源参数谱的观测似然函数的S2阶段：本阶段是利用DSPE(DistributedSourceParameterEstimation，分布源参数估计)算法，基于上述相干分布源的广义阵列流型及阵列天线接收信号模型，先估计阵列天线接收信号的协方差矩阵，并对接收信号的协方差矩阵进行特征值分解，选取小特征值对应的特征向量构建噪声子空间，利用阵列天线导向矢量与噪声子空间估计相干分布源参数谱，并以此作为观测似然函数，然后进入下一步S3阶段。在此阶段中，阵列天线接收信号的协方差矩阵可以表示为：对上述协方差矩阵Rt进行特征值分解，选取L-q个较本文档来自技高网...

【技术保护点】
一种基于粒子滤波的相干分布源波达方向跟踪方法，其特征在于：所述的基于粒子滤波的相干分布源波达方向跟踪方法包括按顺序进行的下列步骤：1)建立相干分布源的广义阵列流型及阵列天线接收信号模型的S1阶段；2)基于上述相干分布源的广义阵列流型及阵列天线接收信号模型，建立基于相干分布源参数谱的观测似然函数的S2阶段；3)构造相干分布源参数状态方程的S3阶段；4)基于步骤2)获得的观测似然函数和步骤3)获得的相干分布源参数的状态方程，利用粒子滤波算法实现相干分布源波达方向跟踪的S4阶段。

【技术特征摘要】
1.一种基于粒子滤波的相干分布源波达方向跟踪方法，其特征在于：所述的基于粒子滤波的相干分布源波达方向跟踪方法包括按顺序进行的下列步骤：1)建立相干分布源的广义阵列流型及阵列天线接收信号模型的S1阶段；2)基于上述相干分布源的广义阵列流型及阵列天线接收信号模型，建立基于相干分布源参数谱的观测似然函数的S2阶段；3)构造相干分布源参数状态方程的S3阶段；4)基于步骤2)获得的观测似然函数和步骤3)获得的相干分布源参数的状态方程，利用粒子滤波算法实现相干分布源波达方向跟踪的S4阶段。2.根据权利要求1所述的基于粒子滤波的相干分布源波达方向跟踪方法，其特征在于：在步骤1)中，所述的建立相干分布源的广义阵列流型及阵列天线接收信号模型的方法是：在等距线阵情况下，首先将分布源的阵列天线接收信号模型描述为角信号密度函数在分布空间积分的形式；然后对于相干分布源模型，则角信号密度函数可以表示为随机信号幅度与分布源的空间分布函数乘积的形式；假设相干分布源的空间分布函数为高斯分布，则可以建立起相干分布源的广义阵列流型及阵列天线接收信号模型。3.根据权利要求1所述的基于粒子滤波的相干分布源波达方向跟踪方法，其特征在于：在步骤2)中，所...

【专利技术属性】
技术研发人员：章涛，李海，吴仁彪，
申请(专利权)人：中国民航大学，
类型：发明
国别省市：天津,12