一种基于观测站数目最小化的TDOA/FDOA无源定位算法

在三维（3D）运动目标无源定位系统中，无模糊定位最少需要4个观测站。而传统的两步加权最小二乘（TSWLS）及其改进的闭式算法至少需要5个观测站进行求解，当减少一个观测站时，这些闭式算法往往无法提供可靠解。针对这一问题，提出一种最小化观测站数目的到达时间差（TDOA）与到达频率差（FDOA）定位算法。该算法是一种闭式解法并且能够在三维场景下仅使用4个观测站进行定位。该算法分为两步：第1步分离传统的TSWLS算法中未知参数空间，建立了一组新的等式，并且利用加权最小二乘（WLS）算法得到目标位置与速度的初始值；第2步利用泰勒级数展开算法将中间变量线性化，对目标位置和速度初始值进一步校正。理论分析证明了在适当的噪声水平下该算法能够达到克拉美罗界（CRLB）。此外，计算机仿真表明仅使用4个观测站时，该算法对于近场以及远场目标参数的估计精度在测量噪声较小时可以实现CRLB；并且还表明使用5个观测站估计时，该算法比TSWLS及其改进算法能更好地适应大的测量噪声。     

传感器参数误差下的运动目标TDOA/FDOA无源定位算法

在运动目标无源定位系统中，许多算法的前提是精确已知传感器的位置以及速度，但实际情况下可利用的传感器的参数均会存在一些噪声扰动。针对这一问题，提出一种改进的两步加权最小二乘（TSWLS）时差（TDOA）与频差（FDOA）定位算法。该算法是一种闭式算法并且分为2步。第1步与经典的两步加权最小二乘算法相同，第2步进一步研究了额外变量与目标参数之间的关系并且建立了新的矩阵方程。随后，利用加权最小二乘技术给出了最终解。理论分析证明了在测量噪声较小时该算法能够达到克拉美罗界（CRLB）。所提算法具有计算复杂度低，实时性高的优点；另外，经过适当的维度调整，该算法同样适用于对多非相交源进行定位求解。计算机仿真进一步证明了理论分析的正确性。     

一种针对恒模信号的运动单站直接定位算法

相比于常规的"测向+位置估计"两步定位模式,以Weiss等提出的目标直接位置确定(DPD)算法具有估计精度高、分辨能力强和无需数据关联等诸多优点。基于该类定位算法的基本理念,提出了一种利用单个运动天线阵列对恒模(即相位调制)信号的DPD算法。首先,依据最大似然(ML)准则以及恒模信号的恒包络特征,建立了相应的直接定位优化模型;接着,根据优化函数的代数特征提出了一种有效的多参量交替迭代算法,用以获得ML估计器的最优数值解;此外,推导了针对恒模信源的位置直接估计方差的克拉美罗界(CRB),从而为新算法的定位精度提供定量的理论下界。仿真实验表明:相比于已有的基于单个运动天线阵列的直接定位算法以及传统的两步定位算法,通过利用恒模信号的恒包络特征可以明显提高目标直接定位的估计精度。     

一种运动双站多脉冲无源定位算法

针对运动双站对已知高度的地面目标辐射源高精度定位问题,借鉴SAR成像原理和直接定位法思想,提出了一种利用空中两运动侦察站接收到目标多脉冲信号间的到达时间差(TDOA)和到达频率差(FDOA)信息的无源定位算法。采样时将观测时间分为快时间和慢时间,在快时间域估计时差,慢时间域估计频差,具体方法为:首先将两侦察站的采样信号在频域互相关,然后利用参考函数和广义Keystone变换消除距离单元徙动,再通过二维傅里叶变换得到目标位置点的TDOA和FDOA联合估计值,最后将时差频差联合估计值通过几何关系映射到目标的空间位置。仿真结果表明信噪比较高时逼近直接定位法的克拉美罗下限(CRLB),信噪比较低时仍然可以定位。此外,本文算法不仅计算量小,且适用于不可区分多目标辐射源定位问题,具有高精度和超分辨特性。     

一种适用于机场场面MLAT监视系统的定位算法及其性能分析

与通常的三维目标不同,机场场面目标的高度是已知固定值,这点在设计应用于民用机场场面监视的无源多点时差定位(MLAT)系统时应予以考虑。基于机场场面目标及MLAT系统的特点,提出了一种改进的时差定位算法,该算法特别适合应用于场面监视的MLAT系统中。相比于传统的三维定位算法,该算法不但可以提高定位精度,还可增强系统的鲁棒性。给出了该改进时差定位算法的克拉美罗界(CRLB),并进一步推导出了该算法与三维定位算法精度之间的理论关系,从理论上证明了该算法的定位精度高于三维定位算法。与此同时,基于条件数理论还证明了该算法的鲁棒性优于三维定位算法。仿真实验及实测数据证明了所得结论的正确性和有效性。     

一种4站情况下基于TDOA/FDOA的无源定位方法

在基于时频差的三维运动目标无源定位系统中,针对在4个接收站的情况下搜索法实时性低的问题,提出了一种基于改进的加权最小二乘法（MWLS）与萤火虫算法（FA）相结合的无源定位方法（MWLS-FA）。该方法的第1步通过构造一组新的方程来对加权最小二乘（WLS）方法进行改进,使得改进后的WLS方法在4站情况下也能得到目标位置和速度的初始值,第2步利用这个初始值为FA方法提供一个动态的搜索区域,同时在约束条件的添加和参数选择两个方面针对性地对FA方法做出了调整和改进。仿真结果表明,该方法在4站情况下对目标的定位精度可以达到克拉美罗下限（CRLB）,而且在实时性和抗噪性方面优于传统的搜索法,同时该方法在5站情况下的抗噪性能优于两步加权最小二乘法（TSWLS）和约束加权最小二乘（CWLS）法。     

基于定位误差修正的运动目标TDOA/FDOA无源定位方法

针对时差(TDOA)、频差(FDOA)无源定位的两步加权最小二乘(TSWLS)方法定位均方根误差(RMSE)和定位偏差适应测量噪声能力差的问题,在分析了影响两步法定位性能的因素基础上提出一种基于一阶泰勒级数展开的定位误差修正方法。该方法的第1步和两步法相同;其第2步避免了两步法第2步中引入估计偏差的平方运算,利用一阶泰勒级数展开得到第1步定位误差的线性最小均方估计,修正第1步定位结果得到目标位置和速度的最终估计,从理论上证明了该方法可以达到定位的克拉美罗下限(CRLB)。计算机仿真对比了新方法和TSWLS方法、基于泰勒级数(TS)展开的迭代极大似然(ML)方法以及约束总体最小二乘(CTLS)方法的定位性能,新算法复杂度和两步法相当,且均方误差和定位偏差低于两步法、泰勒级数法和CTLS方法。     

一种基于多普勒频率的恒模信号直接定位方法

相比于传统的差分多普勒(DD)两步定位方法,以Amar和Weiss提出的基于多普勒频率的单步直接定位方法在低信噪比和小样本条件下具有更高的定位精度。在该类新型定位体制的基础上,提出了一种基于多普勒频率的恒模信号直接定位方法。首先,依据最大似然(ML)准则以及恒模信号的恒包络特征,建立相应的直接定位优化模型。然后,根据目标函数的代数特征将全部未知参量分成两组,并提出一种有效的多参量交替迭代算法,用以获得该优化问题的最优数值解。新算法包含了针对这两组未知参量的Newton型迭代公式,用以避免网格搜索,并能实现多维参数的"解耦合"估计。最后,推导出针对恒模信号的目标位置直接估计方差的克拉美罗界(CRB)。数值实验验证了新方法的优越性。     

一种考虑传感器位置误差的改进源定位算法

在无线传感器网络的定位问题中,通过移动目标与传感器之间的到达时差(TDOA)与到达频差(FDOA)测量量可以估计目标的位置和速度。但是,当传感器自身信息存在误差时,传统的定位方法将失去精确的定位效果。针对带有传感器误差的源定位问题,基于极大似然(ML)法获取一个封闭的近似解。提出了一种改进的近似极大似然(AML)算法,更新带有传感器误差的代价函数,不仅能实现实时定位,而且保证全局收敛。仿真结果表明,本文算法在存在传感器误差场景下依然可以达到克拉美-罗下限(CRLB),比已有的改进两步加权最小二乘(2-step WLS)法更有效。     

基于MCIS技术的相位差变化率单站无源定位

针对最大似然估计（ML）方法求解测相位差变化率单站无源定位问题计算量大、定位慢的问题,本文提出一种利用蒙特卡洛重要性抽样技术（MCIS）高精度、低复杂度的估计方法。根据Pincus定理推导出ML问题的近似全局解,利用重要性抽样（IS）技术构建符合高斯分布概率密度（PDF）的重要性函数,作为样本选取的依据,通过逆变换采样获得样本集,统计样本均值直接得到辐射源位置估计结果。MCIS方法简单易实现且运算量低,能够克服传统ML估计多维网格搜索耗时较长的缺陷,而且对目标位置初始估计误差有较低的敏感性。实验结果表明,MCIS算法在相同测量噪声水平下,定位精度优于EKF、NLS算法,有效减小了初始化估计误差对算法定位精度的影响,也进一步讨论分析了算法参数和不同观测条件对定位性能的影响。     

DOA estimation under unknown mutual coupling and multipath

We propose a new method for direction of arrival (DOA) estimation in the presence of multipath propagation and mutual coupling for a frequency hopping (FH) system. With the use of pilot symbols and assuming perfect time-frequency synchronization for a linear array, we take mutual coupling and multipath propagation into account, and derive a maximum likelihood (ML) estimator for both the mutual coupling matrix and DOA estimation. We then formulate an iterative alternating minimization (AM) algorithm for finding the mutual coupling and DOA parameters in an alternate manner. Simulation results illustrating the performance of the algorithm and comparison with the Cramer-Rao bound (CRB) are also presented.     

Direct position determination using single moving rotating linear array: Noncoherent and coherent processing

To improve the resolution and accuracy of Direct Position Determination (DPD), this paper investigates the problem of positioning multiple emitters directly with a single moving Rotating Linear Array (RLA). Firstly, the geometry of the RLA is formulated and analysed. According to its geometry, the intercepted noncoherent signals in multiple interception intervals are modeled. Correspondingly, the MUltiple SIgnal Classification (MUSIC) based noncoherent DPD approach is proposed. Secondly, the synchronous coherent pulse signals are individually considered and formulated. And the coherent DPD approach which aims for localizing this special type of signal is presented by stacking all array responses at different interception intervals. Besides, we also derive the constrained Cramér-Rao Lower Bound (CRLB) expression for both noncoherent and coherent DPD with RLA under the constraint that the altitudes of the emitters are known. At last, computer simulations are included to examine the performance of the proposed approach. The results demonstrate that the localization accuracy and resolution of DPD with single moving linear array can be significantly improved by the array rotation. In addition, coherent DPD with RLA further improves the resolution and increases the maximum emitter number that can be localized compared with the noncoherent DPD with RLA.     

Beam resource scheduling for a VHF-MIMO radar network in low-angle tracking

The low-angle tracking in multipath interference is a challenging problem for the Very High Frequency (VHF) radar. The colocated Multi-Input Multi-Output (MIMO) technique can remedy such a defect. In this paper, a Joint Beam-Target Assignment and Power Allocation (JBTAPA) strategy is proposed for the VHF-MIMO radar network tracking low-angle targets. The core of the JBTAPA strategy is to improve the worst tracking accuracy among multiple targets by assigning appropriate beams to targets and allocating the power resource in each beam using the feedback information in the tracking cycle. Taking into account the transmit multipath and receive multipath, we derive the Cramér-Rao Lower Bound (CRLB) on angle estimate, which is then incorporated in the Predicted Conditional CRLB (PC-CRLB). A more accurate and consistent lower bound is provided as the optimization metric since the PC-CRLB is based on the most recently realized measurements. A two-stage-based technique is proposed to solve the JBTAPA problem, which is originally NP-hard. Simulation results verify the effectiveness and efficiency of the proposed method. The results also imply that the target reflectivity plays one of the important roles in resource allocation.     

A direct position determination method with combined TDOA and FDOA based on particle filter

The localization of a stationary transmitter using moving receivers is considered. The original Direct Position Determination (DPD) methods, with combined Time Difference of Arrival (TDOA) and Frequency Difference of Arrival (FDOA), do not perform well under low Signal-to-Noise Ratio (SNR), and worse still, the computation cost is difficult to accept when the computational capabilities are limited. To get better positioning performance, we present a new DPD algorithm that proves to be more computationally efficient and more precise for weak signals than the conventional approach. The algorithm partitions the signal received with the same receiver into multiple non-overlapping short-time signal segments, and then uses the TDOA, the FDOA and the coherency among the short-time signals to locate the target. The fast maximum likelihood estimation, one iterative method based on particle filter, is designed to solve the problem of high computation load. A secondary but important result is a derivation of closed-form expressions of the Cramer-Rao Lower Bound (CRLB). The simulation results show that the algorithm proposed in this paper outperforms the traditional DPD algorithms with more accurate results and higher computational efficiency, and especially at low SNR, it is more close to the CRLB.     

Analysis of an adaptive algorithm for unbiased multipath time delay estimation

The multipath equalization time delay estimator (METDE) provides an adaptive approach for estimating the difference in arrival times of a signal received at spatially separated sensors as well as the multipath channel characteristics. However, the parameter estimates of the METDE are biased in the presence of noise. The METDE algorithm is improved for unbiased parameter estimation via minimizing a modified cost function. Convergence behaviors and variances of the system variables of the amended method are derived. Computer simulations are included to corroborate the theoretical analysis and to evaluate the adaptive multipath delay estimation performance of the proposed algorithm.     

Clutter and jammer multipath cancellation in airborne adaptiveradar

Airborne surveillance radars must detect and localize targets in diverse interference environments consisting of ground clutter, conventional jamming, and terrain scattered jammer multipath. Multidimensional adaptive filtering techniques have been proposed to adaptively cancel this interference. However, a detailed analysis that includes the effects of multipath nonstationarity has been elusive. This work addresses the nonstationary nature of the jammer multipath and its impact on clutter cancellation and target localization. It is shown that the weight updating needed to track this interference will also modulate sidelobe signals. At the very least, this complicates the localization of targets. At the worst, it also greatly complicates the rejection of clutter. Several techniques for improving cancellation of jammer multipath and clutter are proposed, including 1) weight vector interpolation, extrapolation, and updating; 2) filter architecture, constraint, and beamspace selection; 3) prefilters; 4) 3-D STAP architectures; and 5) multidimensional sidelobe target editing     

Radar measurement extraction in the presence of sea-surface multipath

In the presence of sea-surface multipath monopulse radar signals from a low elevation target have three alternative paths in addition to the direct (radar-to-target) path due to reflections from the sea surface. The specular reflection causes significant signal fading. The diffuse reflection causes an approximately constant bias to the in-phase component of the monopulse ratio, which is the standard extractor of the direction of arrival (DOA) in the monopulse processing. The diffuse reflection also causes higher standard deviation to the in-phase component of the monopulse ratio. We propose a maximum likelihood (ML) angle extraction technique for low elevation targets of known average signal strength having a Rayleigh fluctuation. The results show that this method reduces the error of the estimated angle compared with the conventional monopulse ratio estimator. Subsequently, the ML angle extractor is modified for the unknown average signal strength case. This modified angle extractor has only a small performance degradation compared with the known average signal strength case, but it performs much better than the monopulse ratio based estimator. An algorithm to calculate the accuracy of the estimated angle (or height) is also presented. This angle extractor reduces the root-mean-square error (RMSE) by more than 50% in the signal processing stage when used in a low flying target tracking scenario. The same algorithm can be used to track sea skimmers.     

密集杂波背景下的水下多目标跟踪方法

水下多目标跟踪是水声信号处理领域研究的热点和难点问题。高斯混合概率假设密度（Gaussian mixture probability hypothesis density, GM-PHD）滤波器以其高效的计算效率为解决水下多目标跟踪问题提供了保证。然而，GM-PHD滤波器在跟踪目标时需要先验已知新生目标的强度，否则其性能会出现严重退化。针对该问题，提出一种滑动窗两步初始化高斯混合概率假设密度（sliding window two step initialization GM-PHD, SWTSI-GMPHD）滤波器。将提出的滑动窗两步初始化方法嵌入GM-PHD滤波器，利用滑动窗两步初始化方法估计新生目标强度，减少杂波干扰导致跟踪结果中出现的虚假目标。仿真实验表明，在杂波密集环境下，相较于其他跟踪方法，提出方法将跟踪精度提高69.84%，52.62%和41.05%。