炮位侦校雷达定位精度影响因素分析

2022-05-23秦鹏程

雷达科学与技术 2022年2期

秦鹏程，王锐，姜洋

(1. 陆军炮兵防空兵学院，安徽合肥 230031; 2. 中国兵器工业第206研究所，陕西西安 710000)

0 引言

炮位侦校雷达是保障炮兵全天候、全天时获取战场情报信息不可或缺的重要侦察装备。弹道外推算法是炮位侦校雷达实现炮位侦察校射功能的关键，算法的适用性和优劣直接影响弹道外推的精度。在弹道外推过程中，任何数据都要受到设备或仪器、方法、环节和人员等因素的影响，因此所得到的外推结果存在误差。

现将主要误差因素描述如下：

1) 地球曲率：在建立炮位侦校雷达数据处理系统时，通常将地球表面看作水平面，而实际上由于地球表面是球面，外推炮位与真实炮位存在高程差；2) 时间不同步误差：测量数据因采样时间不同步间隔(不一致)引起的误差称为时间不同步误差。当前弹道测量系统的各设备和计算机均应工作在统一时间基准上，时间不同步会产生测量误差和数据处理误差；3) 弹道模型误差：建立弹道方程最初的目的主要是为了编制射表，因此弹道方程只是较好地拟合了火炮射程，但对整条弹道曲线的空间位置拟合效果并不理想。即弹道模型和真实弹道并不完全重合，存在模型误差；4) 信噪比变化：雷达的检测能力测量误差实质上受信噪比影响，雷达量测信噪比受雷达与目标距离的影响，并随着距离的增加而减小。

本文结合精密雷达测量数据和七态滤波外推算法对以上误差因素进行详细建模和量化分析，为算法优化奠定了基础。

1 建模分析1.1 地球曲率对定位精度的影响

炮位侦察过程中，考虑地球曲率所引起的外推落点高程差可概略用下式计算。

(1)

式中，为地球半径，为雷达到炮位的距离。

如图1所示，雷达在利用…采样点进行外推时，将地球表面当做一个水平面。外推炮位在，而实际上由于地球表面是球面，炮位并未在点而是继续下落到点，在原来假设落点处实际还存在一个弹道高度'。从平面几何可知：

≈=

(2)

(3)

解算弹道方程时，根据落点高程来确定落点。因此，高程差的存在对定位精度存在较大影响。考虑到地球曲率的影响，在外推到炮位高程时，再继续向下递推。

图1 地球表面曲率对坐标转换的影响

1.2 时间不同步误差

理想条件下，某炮位侦校雷达采样点间隔取固定值0.3 s，但由于发射机电路及波导系统对发射脉冲的延时作用，造成扫描起始时刻超前于天线口辐射的时刻，可能造成实际采样间隔大于0.3 s，以某炮位侦校雷达为例，其实际采样步长接近0.4 s。靶场试验证明：采样点间隔取0.4 s，外推结果更接近真实炮位，由于雷达系统不同，采样步长可能会有所不同。

1.3 真实弹道与模型匹配情况

射表用的弹道模型存在模型误差，并且外推距离越远，外推误差越大，如图2所示。

图2 模型误差影响外推精度示意图

外推距离受遮蔽角影响，遮蔽角越大，雷达起始波束俯仰角增大，开始采样时间距发射时间越长，外推距离增加，模型误差累积越多，雷达的侦察定位精度误差就会越明显。

图 3为遮蔽角影响外推精度示意图。正常情况下，起始波束到截止波束的角度为固定值。和分别为不同遮蔽角下起始波束俯仰角低和高时的外推起点，外推距离分别为点、到外推炮位的弹道曲线长度。从图中可以看出，在作用距离和火炮射角相同的情况下，遮蔽角越大，截止波束俯仰角就越大，外推距离也就越远，导致外推误差越大。其弹丸位置对应的俯仰角可用下式表示：

(4)

式中，为弹丸北向坐标，为弹丸高度，为炮位距雷达距离，为弹丸位置对应的俯仰角。

图3 遮蔽角影响外推精度示意图

图4为雷达作用距离影响外推精度示意图。和分别为近作用距离和远作用距离时的外推起点，外推距离分别为点、到外推炮位的弹道曲线长度。从图中可以看出，在火炮射角和雷达起始波束俯仰角相同的条件下，雷达的作用距离越远，外推距离就会越远，导致外推误差也就越大。

图4 雷达作用距离影响外推精度示意图

以上情况均以侦察作业为例，其分析结果同样适用于校射作业。

1.4 弹丸与雷达距离决定信噪比

设雷达的发射功率为，天线增益为，则在自由空间工作时，距离天线处目标的功率密度为

(5)

假设目标可将接收到的功率无损耗、均匀地辐射出来，用目标的散射截面积来表征其散射特性，接收天线的有效接收面积为，则在雷达接收处接收回波功率为

(6)

由式(6)可以看出，接收的回波功率反比于目标和雷达站间的距离的四次方。

若误差的度量是测量值(估计值)和真实值之间差的均方根值(RMS)，雷达测量的理论均方根误差可表示为

(7)

式中，是大约为1的常数，是接收信号能量，是单位带宽噪声功率。

从一个简单的信号脉冲波形来看，若其宽度为，则接收信号能量=·，噪声功率和噪声功率谱密度之间的关系为=。一般情况下可认为=1，这样可得到信号功率比的表达式如下：

(8)

故正比于，可表示为

(9)

引入零均值高斯白噪声，其方差为，则量测方程为

=()+

(10)

由式(9)可知

(11)

式中，、、分别为观测量中距离、方位角、俯仰角噪声方差与以弹丸距雷达距离四次方为分母函数的系数，表示两者之间的正比关系。

2 仿真实验

2.1 仿真条件

以122榴弹为例，利用靶场弹道精密测量雷达数据进行数值仿真计算。主要仿真条件如下：

1) 弹丸发射条件：初速=700 m/s，射角=195°；2) 取雷达测量误差：=65 m，==15 mil；3) 雷达距离炮位15 km，雷达初始采样点的俯仰角为3.1°，采样间隔 300 ms，从初始采样点继续探测18个弹道点后采样中止。采用两点差分法对状态进行初始化。

建立考虑弹道系数的七维状态向量、三维量测模型，在此过程中使用经典的动力学模型。分别利用扩展卡尔曼滤波(EKF)、无迹卡尔曼滤波(UKF)进行滤波处理，得到外推起点，用龙格-库塔法外推发射点。进行10 000次蒙特卡洛仿真实验，统计结果并分析。

2.2 地球曲率对定位精度的影响

对地球曲率影响下的坐标变换进行建模，代入弹道外推算法中，通过比较考虑地球曲率的弹道外推算法在不同雷达和炮位距离条件下定位的圆中心概率误差(Ecp)，研究地球曲率对定位精度的影响，如表1所示。

由表1可以看出，雷达距离炮位越远，地球曲率对定位精度的影响越大，采用考虑地球曲率模型的算法估计精度优势越明显。以UKF算法为例，当雷达距离炮位15 km时，新算法定位精度提高83.2%；当雷达距离炮位40 km时，新算法定位精度提高168.3%。

表1 雷达与炮位在不同距离条件下受地球曲率的影响

2.3 时间不同步对定位精度的影响

仿真中取采样步长0.3 s，而由于存在采样时间不同步，实际采样步长为0.4 s。对造成时间不同步的采样时间误差进行研究，通过仿真实现了采样间隔和滤波间隔不同步条件下的滤波，求其圆中心概率误差(Ecp)，得到表2。

表2 时间不同步对定位精度的影响

可以看出：考虑时间不同步的外推算法(以UKF算法为例)定位精度提高了127.2%。为进一步分析以上结论，对滤波终点的位置、速度、弹道系数误差进行比较，如图5、图6所示。

图5 EKF中考虑时间不同步误差与否的对比

图6 UKF中考虑时间不同步误差与否的对比

图5、图6分别为用EKF、UKF算法对是否考虑采样点时间不同步进行数据处理的结果，将精密测量雷达数据作为采样点真实位置、速度，结合动力学方程得到真实弹道系数。由图5、图6可以看出：时间不同步对弹道系数的估计精度没有明显影响，但随着处理采样点个数的增加，位置误差、速度误差逐渐增大，造成外推精度的降低。

2.4 弹道模型误差对定位精度的影响

为研究模型误差，利用经典动力学模型和精简动力学模型结合四阶龙格-库塔方程仿真弹道。标准条件下一条弹道由弹道系数()或阻力系数()、初速()和仰角()唯一确定。阻力系数()利用精简动力学方程求取：

(12)

初速()、仰角()取精密弹道数据采样点初值的对应值。得到各弹道和仿真弹道对应点的位置误差如图7、图8所示。

图7 仿真弹道与真实弹道对比

图8 不同模型仿真弹道的位置误差

通过计算得到图7、图8中数据弹道系数为0.571 4，阻力系数为0.07。为方便观察，图7(b)对图7(a)部分弹道进行放大。从以上仿真结果中可以看出：1)仿真弹道在上升沿时有较高的精度，与真实弹道的位置误差较小，在上升到弹道顶点附近时误差开始变大。2)在到达顶点之前，即上升沿经典和精简两种动力学模型都有较高的拟合程度。因此，相对炮位校射仿真，在进行炮位侦察仿真时，精度较高。而在全弹道的位置误差分析中，经典动力学模型位置误差较小，拟合程度较高。

为研究遮蔽角对定位精度的影响，利用精密弹道数据，变换俯仰角大小进行仿真。在炮位距雷达=15 km时，因为雷达从第三秒开始采样，可仿真最小俯仰角为2.646 0°，俯仰角从3°开始每间隔0.5°采样，得到仿真结果如表3所示。