双星对地面未知辐射源直接定位方法研究

罗 迪,尹灿斌,李 智

(航天工程大学,北京 101400)

目前，基于卫星平台的对地定位方法一般分为有源定位和无源定位[1]。其中，有源定位通常运用GNSS系统。该系统由卫星发射定位信号,接收机接收信号并根据信号解算自身位置,现行系统有美国GPS、俄罗斯GLONASS、欧洲Galileo和中国北斗导航定位系统[2]。这种定位方法通常只能对合作目标进行定位。为解决非合作目标的定位问题,学者们提出无源定位方法,通过接收目标发射的电磁波对目标进行定位[3]。这种方法一定程度上解决了非合作目标的定位问题。

无源定位一般采用基于参数测量的两步定位法,即在正确分选信号的前提下,首先估计信号参数,然后根据信号参数估计辐射源位置[4-6]。这种方法受到诸多制约,主要包括:特征参数估计误差,比如天线阵误差等造成的角度测量偏差[7],时间基准不一致等导致的时差测量偏差[8],频率基准不一致等导致的频差测量偏差[9];卫星导航、定轨误差,如卫星轨道位置误差、速度误差、时延误差[10-11]。此外,卫星与目标的几何构型变化、电磁波在大气中的折射、系统本身存在的热噪声、空间环境造成的随机误差、观测方程的线性化近似处理等也是产生误差的原因。而且,当存在多组定位观测量时,各观测量通常是分别独立进行估计的,忽略了所有参数必须对应于相同位置的约束,因此，基于参数测量的两步定位法是次优的[12]。

直接定位方法突破了两步法定位的局限性[13]。该方法免去测量参数提取过程,利用接收的信号数据,对目标位置进行直接估计。相关研究表明，直接定位技术相比于传统两步定位方式，具有精度优势,特别是在低信噪比下，具有更好的精度,并具备多目标的定位能力[14-15]。但因为需要对多路观测信号统一进行处理,该方法需要更大的站间和下行数据传输量和计算复杂度[16]。受限于算法自由度难以满足当前复杂电磁环境、大通信压力(多站数据互传、地面数据接收)以及高计算复杂度等实际问题,算法还难以在工程中应用[17]。

当前，直接定位方法主要应用于二维视距场景[12，16]。本文受相关研究启发,提出一种针对未知信号的脉冲相参累积直接定位算法,并将其应用拓展至卫星平台。首先根据轨道动力学模型建立双星星座,构造了发射线性调频脉冲信号的地面辐射源,利用信号中时差与相位差信息,基于最大似然估计原理构造代价函数,建立针对未知辐射源的脉冲相参累积直接定位算法,运用坐标变换在地球表面划分网格并进行搜索定位。最后通过仿真验证了脉冲相参累积直接定位算法的有效性,并分析了在不同信噪比条件下,脉冲累积量的变化对地面辐射源定位精度的影响。研究表明，在低信噪比条件下，通过脉冲相参累积可以获得更好的定位精度。

1 信号模型

为简化分析,假设辐射源发射第一个脉冲的起始时刻为0时刻,令各脉冲的初始相位均为零,且脉冲之间是相参的。对于单个辐射源l个接收机的情况,辐射源位于p0,接收机保持时频同步,并于pl,n处接收到第n个脉冲信号。对于短脉冲辐射源而言,多普勒频移可以忽略,因此仅考虑接收信号的时延和相位。接收脉冲信号模型可表示为

rl,n(t)=bl,nsn(t-τl,n)e-j2πf0τl,n+wl,n(t),0

(1)

其中

τl,n=‖pl,n-p0‖/c

(2)

假设信号的中心频率f0已通过载频估计得到,则经过接收机下变频采样处理后,上式可化为

(3)

(4)

以上完成了信号时域模型的建立。对式(3)进行离散傅里叶变换,可以得到信号频域模型

(5)

写成矩阵形式为

rl,n=bl,nal,nsn+wl,n

(6)

其中

(7)

2 相参累积定位算法

2.1 算法流程

定位算法基于极大似然估计原理,首先根据信号的时差和相位差构建代价函数,再在目标区域划分网格点,计算每两个网格点之间的时差、相位差的估计值,用估计值补偿实际接收到的信号。将补偿后的信号代入代价函数中计算,补偿后的信号相似程度越高,代价函数的输出值越大,因此会在辐射源所在网格点处出现峰值。最后通过对网格点的遍历寻找到峰值点,实现对辐射源的定位。步骤如图1所示。

图1 算法流程图

2.2 代价函数

下面对定位算法中的代价函数进行推导。式(6)可以写为

wl,n=rl,n-bl,nal,nsn

(8)

(9)

上式两边取自然对数得

(10)

(11)

代入式(10)中可得代价函数C(p)

(12)

其中，b为信道衰减,对sn进行归一化,再求式(12)的最小值,则b的极大似然估计值为

(13)

再代入式(12)得

(14)

其中，rl,n是独立参数,因此可以将最小化式(14)转化为最大化式(15)

(15)

其中，定义N×N的Hermitian矩阵Qn为：

(16)

(17)

对于未知信号的情况,最大化代价函数式(15)可以通过最大化每一个关于sn的二次型来实现。因此要选择一组向量sn使得代价函数最大,也就是需要sn作为Qn的最大特征值λmax{Qn}对应的特征向量。因此代价函数可以转化为

(18)

(19)

则式(18)等价于

(20)

而辐射源位置就位于代价函数的最大值处,即

(21)

3 数值仿真

本节主要通过数值仿真，验证定位方法在双星平台的可行性,并且探究其定位性能。

3.1 场景模型

1)卫星动力学计算。人造卫星的运动主要受地球引力影响,此外，在低轨情况下,还受到大气阻力、光压等的影响。这里为简化分析,忽略其他各种摄动力的影响,将卫星绕地球的运动抽象为两个质点在万有引力作用下的运动学问题,即二体问题。

三维空间中,唯一确定物体轨迹需要6个参数,如位置矢量(x,y,z)和速度矢量(vx,vy,vz)可共同确定物体轨迹。此外,用6个轨道根数也可描述。通常的轨道六根数指的是半长轴a、偏心率e、轨道倾角i、近心点幅角ω、升交点赤经Ω和真近点角φ。双星的轨道六根数见表1所示。

表1 双星轨道六根数

根据轨道六根数,通过轨道动力学方程运算,可以求出卫星任意时刻的坐标和速度[18]。

2)网格坐标转换。为计算与表示方便,辐射源和划分网格点位置采用基于WGS-84地球模型的经纬高坐标系(LLA)。卫星位置和定位计算过程采用地心球面固连坐标系(ECEF)。因此,需通过式(22)将经纬高坐标系转换到地心球面固连坐标系下进行定位计算。

(22)

其中，f为基准椭球体的极扁率,a为基准椭球体的长半轴,具体数值根据WGS-84地球模型数据获得。lon、lat、alt分别代表辐射源的经度、纬度和高程。假设地面辐射源的经度为68°E,纬度为20°N,高程为0,坐标表示为(lon=-68,lat=20,alt=0)。通过对卫星的轨道动力学计算和辐射源坐标转换,得到两颗卫星及辐射源在ECEF和LLA坐标系下的坐标分别如表2和表3所示。

表2 卫星和辐射源的初始位置与速度(ECEF)

表3 卫星和辐射源的初始位置与速度(LLA)

图2显示了卫星与辐射源在空间中的相对位置,以及在地面划分的网格点。

图2 卫星与辐射源的相对位置

3)信号数据生成。辐射源发射信号为典型的线性调频脉冲信号,设置载波频率为4 GHz,脉冲重频为5 kHz,空占比为50%,调频带宽为20 MHz。接收机采样率为50 MHz,接收信号为一个完整脉冲。图3显示了发射信号的时域波形。

图3 信号波形

3.2 单脉冲直接定位仿真分析

理想条件(无噪声影响,信号未知)下,单脉冲定位仿真得到代价函数分布如图4所示,代价函数在接近真实发射位置处有一个明显的峰值。

图4 代价函数计算结果(脉冲数=1, 无噪声)

增加噪声使得信噪比(SNR)=-5/(-10)dB,接收脉冲数依旧为1,代价函数计算结果如图5、图6所示。

图5 代价函数计算结果(脉冲数=1, SNR=-5 dB)

图6 代价函数计算结果(脉冲数=1, SNR=-10 dB)

可见，在SNR=-10 dB时,代价函数中的峰值已经被噪声淹没,导致无法成功定位。

3.3 脉冲相参累积直接定位仿真分析

上文中,低信噪比条件下定位遇到困难,因此考虑利用脉冲相参累积提升算法在低信噪比条件下的定位能力。仿真实验中,将脉冲累积数增加到10。考虑卫星飞行速度较快,不能假设卫星始终静止,因此采取“动-停”仿真,即假设卫星在接收脉冲时静止,在其余时间段内运动。通过脉冲相参累积,得到仿真结果如图7所示。

图7 代价函数计算结果(脉冲数=10, SNR=-10 dB)

由仿真结果可以发现,通过累积10个脉冲,在同样SNR=-10 dB的条件下,相比与上一节中峰值被噪声淹没的情况,代价函数中的峰值较为明显。说明通过脉冲相参累积可以有效提升直接定位的效果。

3.4 脉冲相参累积量对定位精度影响仿真分析

在验证脉冲相参累积可以提升直接定位效果后,本小节考虑研究脉冲累积量与信噪比、定位精度之间的关系。在之前同样的仿真场景下,将接收信号的信噪比设置为从-25 dB到15 dB,脉冲累积量分别为1、10和70,蒙特卡洛仿真次数为30次,由此可以得到在不同脉冲累积量下定位误差与信噪比的关系，如图8所示。

图8 不同脉冲累积量下定位误差与信噪比关系曲线

通过统计定位误差可以发现,随着接收信号信噪比的降低,定位精度也随之降低。而在相同信噪比条件下,通过脉冲累积量的增加,可以明显提升定位精度。仿真结果有效说明了脉冲相参累积对定位性能的贡献度。

4 结束语

本文针对卫星对未知地面辐射源定位,提出脉冲相参累积直接定位算法,实现了对地面未知辐射源的直接定位,通过仿真验证了脉冲相参累积直接定位算法的有效性,并进一步分析了脉冲累积量的变化在不同信噪比下对辐射源定位精度的影响。仿真结果表明,在低信噪比条件下,通过脉冲相参累积可以减少定位误差。下一步考虑将算法由对单一辐射源的定位拓展至对多辐射源的联合定位,并提升算法对不同信号的适应性和定位精度。