基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法与流程

本发明属于目标辐射源定位技术领域，尤其涉及一种基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法。

背景技术：

众所周知，目标辐射源定位技术已广泛应用于通信、雷达、目标监测、导航遥测、地震勘测、射电天文、紧急救助、安全管理等领域，其在工业生产与军事应用中都发挥着重要作用。目标辐射源定位技术是指在观测站(也称传感器)不主动发射电磁信号的情况下，通过接收目标辐射的无线信号来确定目标位置参数(有时也包括速度参数)。该类技术属于无源定位范畴，由于其系统并不主动发射电磁信号，因此具有生存能力强、侦察作用距离远等优势。依据观测站的个数进行划分可以将辐射源定位系统分成单站定位系统和多站定位系统两大类，其中单站定位系统具有灵活度高、机动性强、系统简洁以及无需站间通信与同步等优点，本专利主要涉及单站无源定位体制。

在现有的单站无源定位体制中，短波单站定位是应用较为广泛的一类定位方法，该方法主要是针对超视距短波目标源进行定位。其基本原理是利用单个观测站测得的信号方位角和仰角以及电离层虚高参数对短波辐射源进行定位。然而，在实际应用中电离层虚高参数是通过有源探测所获得的，因此难免会存在一定的偏差，除此以外，电离层倾斜角也会对方位角估计产生一定偏差。不难想象，无论是电离层虚高误差还是电离层倾斜角偏差都会对短波单站定位产生较大的影响。

另一方面，传统的无源定位技术大多采用两步估计方式，即首先从接收信号中提取出用于定位的相关参数(主要包括空域、时域、频域以及能量域等参量)，然后利用这些中间参数确定目标位置参数或者速度参数。虽然这种两步定位模式在现代无源定位系统中被广泛使用，但以色列学者a.j.weiss和a.amar却指出了其中所存在的诸多缺点，并提出了直接定位的思想(amara,weissaj.localizationofnarrowbandradioemittersbasedondopplerfrequencyshifts[j].ieeetransactionsonsignalprocessing,2008,56(11):5500-5508.)(weissaj.directgeolocationofwidebandemittersbasedondelayanddoppler[j].ieeetransactionsonsignalprocessing,2011,59(6):2513-5520.)，其基本理念是从信号采集的数据域中直接估计目标的位置参数，而无需估计其它中间定位参数。显然，这种直接定位体制同样适用于短波单站定位场景中。不幸的是，短波单站直接定位方法仍然会受到电离层虚高误差和电离层倾斜角偏差的影响，从而会产生较大的定位偏差。针对此问题，本专利通过利用短波目标源附近区域的短波校正源信息，提出了一种基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法，可以大幅提升短波单站直接定位精度。

技术实现要素：

本发明针对现有单站无源定位方法中短波单站直接定位方法受到电离层虚高误差和电离层倾斜角偏差的影响，从而会产生较大的定位偏差的问题，提出一种基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法。

为了实现上述目的，本发明采用以下技术方案：

一种基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法，包括：

步骤1：在短波目标源所在空域位置所在区域分时放置n个地理坐标精确已知的短波校正源；

步骤2：利用观测站的地理坐标和已知的电离层虚高信息，建立第n个短波校正源的地理坐标与其辐射信号到达观测站二维波达方向之间的代数关系式，其中1≤n≤n；

步骤3：针对第n个短波校正源，利用k元均匀圆阵对其辐射信号进行接收和采集，并且利用子空间方法获得关于短波校正源地理坐标空域位置谱；

步骤4：针对第n个短波校正源，按照一定步长将其地理坐标空域位置谱主峰所在区域网格化，并且利用每个网格的谱值构造关于短波校正源地理坐标空域位置谱矩阵

步骤5：利用n个短波校正源的地理坐标空域位置谱矩阵以及对应的真实地理坐标训练径向基神经网络；

步骤6：利用观测站的地理坐标和电离层虚高信息，建立短波目标源的地理坐标与其辐射信号到达观测站二维波达方向之间的代数关系式；

步骤7：针对短波目标源，利用k元均匀圆阵对其辐射信号进行接收和采集，并且利用子空间方法获得关于短波目标源地理坐标空域位置谱；

步骤8：针对短波目标源，按照步骤4中步长将其地理坐标空域位置谱主峰附近区域网格化，并且利用每个网格的谱值构造关于短波目标源地理坐标空域位置谱矩阵p^(e)；

步骤9：利用向量化算子将步骤8获得的地理坐标空域位置谱矩阵转化成向量，并输入到步骤5所训练好的径向基神经网络中，网络的输出值即为短波目标源地理坐标的最终估计值。

进一步地，所述步骤2中第n个短波校正源的地理坐标与其辐射信号到达观测站二维波达方向之间的代数关系式为：

式中

其中θ^(s)和β^(s)分别为观测站的经、纬度，分别为第n个短波校正源的经、纬度，h为电离层虚高，r为地球半径，分别为第n个短波校正源辐射信号到达观测站的方位角和仰角，t1、t2为坐标转换向量，为短波单站与短波目标源之间的地心角的1/2。

进一步地，所述步骤3包括：

步骤3.1：利用k元均匀圆阵对第n个短波校正源的辐射信号进行接收和采集，采集m个信号样本点并且构造协方差矩阵

步骤3.2：对协方差矩阵进行奇异值分解，将奇异值由大到小进行排列，利用其后面k-1个小奇异值对应的左奇异向量构造矩阵

步骤3.3：通过构造关于短波校正源地理坐标空域位置谱函数从而得出关于短波校正源地理坐标空域位置谱，其中b(θ,β)表示以短波校正源地理坐标为函数的阵列流形向量。

进一步地，所述步骤5包括：

利用向量化算子vec(·)将第n个短波校正源对应的谱图矩阵转化成向量并将作为径向基神经网络的输入值，然后将第n个短波校正源的真实经度和纬度作为径向基神经网络的输出值，采用自组织选取中心法对径向基神经网络进行训练。

进一步地，所述步骤6中短波目标源的地理坐标与其辐射信号到达观测站二维波达方向之间的代数关系式为：

式中

其中θ^(e)和β^(e)分别为短波目标源的经、纬度，α^(e)、γ^(e)分别为短波目标源辐射信号到达观测站的方位角和仰角。

进一步地，所述步骤7包括：

步骤7.1：利用k元均匀圆阵对第n个短波校正源的辐射信号进行接收和采集，采集m个信号样本点{x^(e)(tm)}1≤m≤m，并且构造协方差矩阵

步骤7.2：对协方差矩阵进行奇异值分解，将奇异值由大到小进行排列，利用其后面k-1个小奇异值对应的左奇异向量构造矩阵

步骤7.3：通过构造关于短波校正源地理坐标空域位置谱函数从而得出关于短波校正源地理坐标空域位置谱。

进一步地，所述步骤9包括：

利用向量化算子vec(·)将步骤8获得的p^(e)转化成向量p^(e)，并输入到步骤5中训练好的径向基神经网络中，网络的输出值即为短波目标源地理坐标的最终估计值，通过该估计值补偿由电离层虚高误差和电离层倾斜角偏差所引起的定位偏差。

与现有技术相比，本发明具有的有益效果：

本发明提出了基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法，利用短波目标源附近的短波校正源地理坐标空域位置谱矩阵训练径向基神经网络，基于该神经网络可以有效消除由电离层虚高误差和电离层倾斜角偏差所引起的定位偏差，从而大幅度提高短波单站直接定位的精度。

附图说明

图1为本发明实施例一种基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法的基本流程图；

图2为本发明实施例一种基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法的地理坐标空域位置谱；

图3为本发明实施例一种基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法的直接定位结果散布对比图；

图4为本发明实施例一种基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法的短波目标源定位均方根误差随着信噪比的变化曲线对比图；

图5为本发明实施例一种基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法的短波目标源定位均方根误差随着信号样本点数的变化曲线对比图；

图6为本发明实施例一种基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法的短波目标源定位均方根误差随着电离层虚高误差的变化曲线对比图；

图7为本发明实施例一种基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法的短波目标源定位均方根误差随着电离层倾斜角偏差的变化曲线对比图。

具体实施方式

下面结合附图和具体的实施例对本发明做进一步的解释说明：

如图1所示，一种基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法，包括：

步骤s101：在短波目标源所在空域位置附近分时放置n个地理坐标精确已知的短波校正源。

步骤s102：利用观测站的地理坐标和已知的电离层虚高信息(含有误差)，依次建立第n(1≤n≤n)个短波校正源的地理坐标与其辐射信号到达观测站二维波达方向之间的代数关系式。

步骤s103：依次针对第n(1≤n≤n)个短波校正源，利用k元均匀圆阵对其辐射信号进行接收和采集，并且利用子空间方法获得关于短波校正源地理坐标空域位置谱。

步骤s104：依次针对第n(1≤n≤n)个短波校正源，按照一定步长将其地理坐标空域位置谱主峰附近区域网格化，并且利用每个网格的谱值构造关于短波校正源地理坐标空域位置谱矩阵。

步骤s105：利用n个短波校正源的地理坐标空域位置谱矩阵以及它们的真实地理坐标训练径向基神经网络。

步骤s106：利用观测站的地理坐标和电离层虚高信息，建立短波目标源的地理坐标与其辐射信号到达观测站二维波达方向之间的代数关系式。

步骤s107：针对短波目标源，利用k元均匀圆阵对其辐射信号进行接收和采集，并且利用子空间方法获得关于短波目标源地理坐标空域位置谱。

步骤s108：针对短波目标源，按照步骤s104中步长将其地理坐标空域位置谱主峰附近区域网格化，并且利用每个网格的谱值构造关于短波目标源地理坐标空域位置谱矩阵。

步骤s109：利用向量化算子将步骤s108获得的地理坐标空域位置谱矩阵转化成向量，并输入到步骤s105所训练好的径向基神经网络中，网络的输出值即为短波目标源地理坐标的最终估计值。

具体地，在所述步骤s101中，在短波目标源所在空域位置附近分时放置n个地理坐标精确已知的短波校正源，其中第n(1≤n≤n)个短波校正源的经度为纬度为

具体地，在所述步骤s102中，假设观测站的经纬度分别为θ^(s)和β^(s)，电离层虚高为h，地球半径为r，第n(1≤n≤n)个短波校正源辐射信号到达观测站的方位角和仰角分别为和根据信号传播的几何关系可以建立如下代数关系式：

式中

其中，t1、t2为坐标转换向量，为短波单站与短波目标源之间的地心角的1/2，g()为中间参量。

具体地，在所述步骤s103中，依次针对第n(1≤n≤n)个短波校正源，利用k元均匀圆阵对其辐射信号进行接收和采集，其阵列接收信号模型为：

式中表示均匀圆阵针对第n个短波校正源的接收信号；表示第n个短波校正源辐射信号的复包络；表示阵列加性噪声；表示以短波校正源辐射信号二维波达方向为函数的阵列流形向量；表示以短波校正源地理坐标为函数的阵列流形向量，其满足

接着利用子空间方法获得关于短波校正源地理坐标空域位置谱，其计算过程如下：

(1)采集m个信号样本点并且构造协方差矩阵

(2)对协方差矩阵进行奇异值分解，将奇异值由大到小进行排列，利用其后面k-1个小奇异值对应的左奇异向量构造矩阵

(3)构造关于短波校正源地理坐标空域位置谱函数从而得出关于短波校正源地理坐标空域位置谱。

具体地，在所述步骤s104中，依次针对第n(1≤n≤n)个短波校正源，按照一定步长将其地理坐标空域位置谱主峰附近区域网格化，如图2所示，并且利用每个网格的谱值构造关于短波校正源地理坐标空域位置谱矩阵为了减少运算复杂度，步长可以适当放宽，无需过于精细。

具体地，在所述步骤s105中，利用向量化算子vec(·)将第n(1≤n≤n)个短波校正源对应的谱图矩阵转化成向量(即有)，并将其作为径向基神经网络的输入值，然后将第n个短波校正源的真实经纬度和作为径向基神经网络的输出值，共有n组这样的输入-输出对，利用它们对径向基神经网络进行训练，学习算法采用自组织选取中心法。利用径向基神经网络对学习样本进行训练以后，该网络就具备了对短波单站直接定位偏差的自动补偿功能，可以有效弥补由电离层虚高误差和电离层倾斜角偏差所带来的影响。

具体地，在所述步骤s106中，假设短波目标源的经纬度为θ^(e)和β^(e)，其辐射信号到达观测站的方位角和仰角分别为α^(e)和γ^(e)，根据信号传播的几何关系可以建立如下代数关系式：

式中

具体地，在所述步骤s107中，针对短波目标源，利用k元均匀圆阵对其辐射信号进行接收和采集，其阵列接收信号模型为

x^(e)(t)＝a(α^(e),γ^(e))s^(e)(t)+ξ^(e)(t)＝b(θ^(e),β^(e))s^(e)(t)+ξ^(e)(t)

式中x^(e)(t)表示均匀圆阵针对短波目标源的接收信号；s^(e)(t)表示短波目标源辐射信号的复包络；ξ^(e)(t)表示阵列加性噪声；a(α^(e),γ^(e))表示以短波目标源辐射信号二维波达方向为函数的阵列流形向量；b(θ^(e),β^(e))表示以短波目标源地理坐标为函数的阵列流形向量，其满足b(θ^(e),β^(e))＝a(α^(e),γ^(e))。

接着利用子空间方法获得关于短波目标源地理坐标空域位置谱，其计算过程如下：

(1)采集m个信号样本点{x^(e)(tm)}1≤m≤m，并且构造协方差矩阵

(2)对协方差矩阵进行奇异值分解，将奇异值由大到小进行排列，利用其后面k-1个小奇异值对应的左奇异向量构造矩阵

(3)构造关于短波目标源地理坐标空域位置谱函数从而得出关于短波校正源地理坐标空域位置谱。

具体地，在所述步骤s108中，针对短波目标源，按照一定步长将其地理坐标空域位置谱主峰附近区域网格化，并且利用每个网格的谱值构造关于短波目标源地理坐标空域位置谱矩阵p^(e)。这里所选择的区域范围和步长与步骤s104需要完全一致。

具体地，在所述步骤s109中，利用向量化算子vec(·)将步骤s108获得的地理坐标空域位置谱矩阵p^(e)转化成向量(即有p^(e)＝vec(p^(e)))，并输入到步骤s105所训练好的径向基神经网络中，网络的输出值即为短波目标源地理坐标的最终估计值，该估计值可以补偿由电离层虚高误差和电离层倾斜角偏差所引起的定位偏差。

假设观测站的经度为东经112.73°，纬度为北纬33.25°；短波目标源的经度为东经125.52°，纬度为北纬28.12°。观测站安装10元均匀圆阵，其半径与波长比为1.5，短波目标源辐射信号到达观测站所经历的电离层虚高为350公里。

(1)短波目标源辐射信号的信噪比为10db，算法采用的信号样本点数为500，电离层虚高误差为50公里，电离层倾斜角偏差为0.5°，图3给了定位结果散布图。从图3中可以看出，本专利公开的基于地理坐标空域位置谱的短波单站直接定位偏差补偿方法可以明显消除由电离层虚高误差和电离层倾斜角偏差所带来的影响，从而显著提高了对短波辐射源的定位精度。

(2)其余实验条件不变，图4和图5分别给出了短波目标源定位均方根误差随着信噪比和信号样本点数的变化曲线，从图中可以看出本专利公开方法的定位精度得到了大幅提升，并且随着信噪比和信号样本点数的提高，其定位精度会逐渐提高，而传统直接定位方法由于未考虑电离层虚高误差和电离层倾斜角偏差的影响，因此其定位精度并不会随着信噪比和信号样本点数的提高而逐渐提升，具有固定偏差。

(3)其余实验条件不变，图6和图7分别给出了短波目标源定位均方根误差随着电离层虚高误差和电离层倾斜角偏差的变化曲线，从图中可以看出本专利公开方法受到电离层虚高误差和电离层倾斜角偏差的影响很小，这是因为该方法利用短波校正源地理坐标空域位置谱对径向基神经网络进行训练，使其可以有效补偿电离层虚高误差和电离层倾斜角偏差所产生的定位偏差，而传统直接定位方法则无法消除其影响，因此其定位误差呈现出线性增长趋势。

以上所示仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。