1、第 卷 第 期兵 器 装 备 工 程 学 报 年 月 收稿日期：；修回日期：基金项目：国防科技 计划技术领域基金项目（）作者简介：李徽（），男，硕士，高级工程师，：。：一种基于距离参数化 的单站无源定位方法李 徽（江苏自动化研究所，江苏 连云港）摘要：针对单站无源定位系统中，观测站机动导致的滤波精度降低问题，提出了一种基于距离参数化 的单站无源定位方法（，）。引入距离参数化的思想，结合观测站观测范围作为先验值，将观测范围划分为若干个区间，并赋予初始权重，在各个区间引入后向平滑容积卡尔曼滤波（，），利用各时刻的预测值与观测值的比值来更新区间权重，最后对各区间状态信息加权融合来实现目标状态的获取。
2、仿真结果表明：该方法能够有效降低滤波全局对观测站机动的敏感性，提高滤波稳定性与定位精度。关键词：单站无源定位；容积卡尔曼滤波；距离参数化；后向平滑本文引用格式：李徽 一种基于距离参数化 的单站无源定位方法 兵器装备工程学报，（）：，（）：中图分类号：；文献标识码：文章编号：（）（，）：，（），（），：；引言无源定位技术主要指观测系统在不发出电磁波信号，通过被动接收来自目标的辐射信号，并根据这些信号的参数来解算目标位置信息的技术。具有功耗低，隐身性强的优点，。单站无源定位系统由于其设备简单、机动性强的特点，在无源定位领域应用广泛。对单站无源目标定位问题，目前应用较为广泛的主要有基于扩展卡尔曼滤波
3、的算法（，）、基于无迹卡尔曼滤波的算法（，）、基于容积卡尔曼滤波的算法（，）等。算法具有原理简单、易于实现的特点，但是该算法在滤波过程中通过计算非线性函数的雅克比矩阵来进行协方差矩阵更新，当系统非线性程度较高时，该算法的精度和稳定性无法保证。算法通过设置一组加权采样点，来逼近估计的均值、协方差矩阵或更高阶矩阵，相对于，该算法具有更优异的跟踪性能，但是该算法在迭代的过程中需要对状态估计协方差矩阵进行分解运算，而状态估计协方差矩阵是难以保持正定的，限制了该算法的应用。是研究者们近年来重点关注的对象，该算法通过容积点规则实现泰勒公式二阶近似，相对于 与，具有更高好的数值稳定性与滤波精度。由于只有一个
4、观测站，单站系统无法通过测向结果及几何关系获取目标初始信息，只能通过观测站的机动，积累多个时刻的测量参数，然后滤波来获取目标状态信息，这种定位方式也导致了滤波易受观测站机动的影响而精度降低的先天缺陷。鉴于 算法优异的鲁棒性，霍光结合 与平滑思想提出了一种后向平滑 算法，提高了目标跟踪的精度，然而该方法并没有摆脱由于观测站机动导致的滤波精度降低问题。樊垚等面向水面目标典型机动形式，提出了一种距离参数化 滤波器的定位方法，采用多个区间分段滤波的方式，降低了滤波对观测站机动的敏感性，提高了定位精度。但是该算法在野战环境下观测站机动较多时存在鲁棒性差，滤波容易发散的问题。结合 算法与距离参数化的优点，
5、提出了一种基于距离参数化 的单站无源定位算法。引入观测范围这一先验值，将目标距离范围划分为若干个独立区间，并赋予初始权重。然后，在各区间引入多组 滤波器并对各区间权重进行迭代，最后，通过各区间结果的加权融合来实现目标状态信息的获取。系统模型 几何关系目标与观测站之间的几何关系图如图 所示，和 分别表示观测站获得的方位角与俯仰角，是观测站的位置，观测站坐标为（，），目标位置信息为（，）。图 目标和观测站的几何关系图 状态模型假设目标在 时刻的状态记为，（），?（），（），?（），（），?（），分别表示目标在，三个方向的位置和速度。在匀速运动模型（）下状态方程为 （）在这里，表示目标在 时刻的状态
6、转移矩阵，通过目标的运动学特性获取，在此处为一个 的矩阵：|（）表示过程噪声，服从高斯分布。观测模型系统的观测方程为（）（）为 时刻的量测信息，（）为一个非线性的量测方程，在此处量测信息分别为观测站获得的方位角与俯仰角，观测方程为（）（）（）（）（）（）（）（）（）其中，（）和（）分别表示观测站在 时刻获取的方位角与俯仰角，表示测量噪声，方差为。基于 的单站无源定位方法将探测范围划分为 个区间，并认为目标在区间内服从均匀分布，由该信息可获取目标在该区间的均值与协方差信息。对每个区间赋予初始权重，利用各时刻的预测值与观测值的比值来更新区间权重，最后通过各区间的加权融合来获取目标状态信息。步骤 假
7、设起始目标点位与观测站在（，）上分布，将该区间分为 小区间，其中第 个小区间为（，李 徽：一种基于距离参数化 的单站无源定位方法），其中 的表达式为()（）在该区间内目标初始距离的均值与协方差分别为?（）（）步骤 基于各区间设置目标距离统计信息，并基于量测信息与目标速度统计信息，实现目标初始信息确认。目标位置信息可由测量方程确认。（）?（）?（）?（）目标速度信息可由相邻时刻角度测量值代入测量方程进行估算：（）?（）（）?（）（）?（）（）其中，、和、分别是相邻时刻的目标方位角与俯仰角。步骤 利用容积卡尔曼滤波算法，时刻状态 进行第一次前向滤波。容积卡尔曼滤波算法的核心思想是利用 准则，采用
8、个等权值的容积点，来实现非线性估计，算法具体为）计算容积点：，（）其中，为 时刻目标状态的协方差矩阵，表示对矩阵进行 分解，即 ，表示第 个容积点，表示等值权值点个数，表示矩阵 ，第 列元素，表示 维单位阵，。计算通过状态方程传播的容积点，（）状态预测 ，（）方差预测 ，（）量测更新。计算容积点：，（）其中，。计算通过测量方程传播的容积点，（，）（）量测预测 ，（）估计新息协方差，“”（）估计互协方差，（）计算卡尔曼滤波增益，（）状态更新 （）（）估计协方差 ，（）步骤 对滤波结果进行后向平滑处理，后向平滑过程如下：（）（）（）（）（）（）（）步骤 第二次向前滤波，将平滑后得到的 和 代替第一
9、次前向滤波中的 和 ，在式（）式（）中，进行第二次前向滤波，得到滤波结果和。步骤 距离参数化容积卡尔曼滤波需给每个区间分配初始权重，由于假设目标距离在区间内服从均匀分布，因此第 个小区间的初始权重为（）（）为了获取下一时刻的权重可根据量测似然函数对区间权重进行更新。，（）（，）（），（）（，）（）（）其中，为 时刻第 个区间的权重，为 时刻第 个区间获取的量测协方差，为 时刻第 个区间获取的预测量测，为 时刻的测量值。利用新的权重与各区间状态信息加权求和实现 时刻目标状态的获取（）其中，为第 个区间获取的 时刻目标状态。兵 器 装 备 工 程 学 报：仿真与分析为了验证算法的有效性，参照霍光所
10、提出的 算法的仿真条件设计一组对照试验。假设目标初始位置为（，），速度为（，），初始速度误差协方差是，目标做匀速直线运动（）。观测站探测范围是（，），最大移动速度为 ，从（，）开始运动，在 做加速度为（.，）的匀加速运动，在 、做转弯速度为 的匀转弯运动，其余时间做匀速运动。采样周期为 ，总时间为，蒙特卡洛次数为。考虑目标最大速度小于 ，目标初始速度误差协方差的公式如下：|（）是调节参数，假设 分别取，。仿真结果如图 所示。采用位置和速度的均方根误差（，）及时间平均均方根误差（，）来评估仿真性能，与 定义如下：（）（）（?）（）（）（）（）（）其中，、和?分别表示在第 次蒙特卡洛试验中第 时刻
11、目标的位置或速度在，上的估计结果，是蒙特卡洛总次数，是一次蒙特卡洛试验采样总次数。图 仿真结果 由图 可以看出，在 、区间附近，的 值产生了较大波动，与 相对 的波动幅度降低，这说明采用了分区间滤波的方式能够有效降低观测站对滤波精度的影响，提高算法的稳定性；同时，当滤波收敛后，的位置精度与速度精度都明显优于 算法与 算法。表 和 表 的结果表明：的跟踪精度高于 和。表 位置的 李 徽：一种基于距离参数化 的单站无源定位方法表 速度的 结论针对单站无源定位系统中，由于观测站机动导致的滤波精度降低问题，提出了一种基于距离参数化 的单站无源定位方法（）。通过引入观测范围进行分段滤波，降低了算法对观测
12、站机动的敏感性，相对与传统算法，本文方法具有以下优势：采用区间分段滤波的方式，降低了观测站机动对滤波全局的影响；相对于 算法，采用 具有更好的稳定性与鲁棒性。本文在分段滤波的过程中，引入了后向平滑，目的是进一步提高滤波的精度，然而这是以牺牲更多的时间实现的，下一步将考虑 的改进或引入新的算法，以期在应用中能取得更好的效果。参考文献：，（）：，（）：，（）：，（）：，（）：，（）：，（）：，（，），（）：，（）：，（）：，（）：，（）：，（）：，：，（）：霍光，李冬海 基于后向平滑容积卡尔曼滤波的单站无源定位算法 信号处理，（）：，（）：樊垚，费玮玮，杨洪康，等 基于距离参数化 滤波器的纯方位目标运动分析 舰船科学技术，（）：，（）：杨明亮，龙华，邵玉斌，杜庆治 基于 全局参数联合的说话人识别 重庆邮电大学学报（自然科学版），（）：，（），（）：科学编辑 杨继森 博士（重庆理工大学教授）责任编辑 唐定国兵 器 装 备 工 程 学 报：
