CN104392108B

CN104392108B - 一种采用迭代差分算法的远程定位***及方法

Info

Publication number: CN104392108B
Application number: CN201410609871.3A
Authority: CN
Inventors: 陆岩; 汤锦辉; 程俊斌; 方雷; 兰洪亮; 邵欣; 刘晓丽; 刘兆云; 周臣
Original assignee: INSTITUTE OF RADAR AND ELECTRONIC COUNTERMEASURE OF CHINESE PLA AIR FORCE EQUIPM
Current assignee: INSTITUTE OF RADAR AND ELECTRONIC COUNTERMEASURE OF CHINESE PLA AIR FORCE EQUIPM
Priority date: 2014-11-03
Filing date: 2014-11-03
Publication date: 2018-05-01
Anticipated expiration: 2034-11-03
Also published as: CN104392108A

Abstract

本发明涉及远程定位领域，具体涉及一种采用迭代差分算法的远程定位***及方法，包括地面控制中心、远程网络和用户终端；远程网络由两个空间目标及一个空间备用目标构成，其中，空间目标上载荷有变频转发器、两个波束的S频段天线和L频段天线，两颗卫星的4个S频段波束分区覆盖全服务区，空间目标的两个L频段波束分区覆盖全服务区，本发明的有益效果：可实现授时、定位、GPS授时、GPS定位的功能，以摆脱对GPS的以来，提高我国民用基础行业的安全性和稳定性。

Description

一种采用迭代差分算法的远程定位***及方法

技术领域

本发明属于远程定位领域，具体涉及一种采用迭代差分算法的远程定位***及方法。

背景技术

目前全球的导航卫星主要有美国的全球卫星导航***GPS(Global PositioningSystem)、俄罗斯的全球导航卫星***GLONASS(Global Navigation Satellite System)、中国的北斗卫星导航定位***(Compass Navigation Satellite System)和欧盟的伽利略全球导航定位***(European Satellite Navigation System)。现阶段，GPS定位以及授权功能在我国的交通、电力、通信、金融、水利、石油等领域有着广泛的应用；但因GPS受美国限制，采用单一的GPS存在风险，其安全可靠性低，自主性差。采用迭代差分算法的远程定位***是我国正在实施的自主发展、独立运行的全球卫星导航***。

发明内容

本发明是所要解决的技术问题是针对现有技术的不足提供一种采用迭代差分算法的远程定位***及方法。

一种采用迭代差分算法的远程定位***，所述***包括：

控制中心：连接远程网络，同时管理整个***及采用迭代差分算法进行导航定位解算；

远程网络：覆盖全服务区，用于转发控制中心与用户终端之间的双向信号传输；

用户终端：连接远程网络，向远程网络接收及发送信息。

进一步地，所述远程网络包括两个工作单元及一个备用单元。

进一步地，所述两个空间目标及一个空间备用目标均载荷有变频转发设备、两个波束的S频段天线及L频段天线。

进一步地，所述控制中心包括信号收发子***、信息处理分***、时间分***、监控子***及轨道监控分***。

一种采用迭代差分算法的远程定位方法，应用上述的***，包括以下步骤：

a)建立一个位置已知的空间参考系，测量待定位点到各参考点的相对位置，其中控制中心及远程网络的空间目标为参考点，用户终端为待定位点；

b)创建待定位点观测量，定位***中的控制中心采集获得用户终端的信号强度、载波相位、信号到达角度和时间观测量，创建待定位点观测量，表示为：

s⁽ⁱ⁾为待定位点的观测量；为控制中心到空间目标i的几何距离；待定位点到空间目标i的几何距离；空间目标i的转发延迟；t接收机钟差；c光速信号；的空间传输延迟；u⁽ⁱ⁾测量误差及噪声；

c)获得空间目标到待定位点的几何距离为：

(x_i，y_i，z_i)为空间目标i的坐标；

d)将步骤b)中的公式变换得到：

e)通过气压计及高度计测量得到高度信息后，d)步骤公式可解，获得定位结果。

进一步地，所述空间参考系包括地心惯性坐标系和地心地球固连坐标系。

进一步地，当e)步骤中无法通过气压计及高度计测量得到高度信息时，可构造方程，包括椭球观测方程及虚拟卫星观测方程，并通过线性迭代的方法求解，获得定位结果。

进一步地，所述椭球观测方程为：

其中高度为h；a，b分别为地球椭球形的长轴和短轴。

进一步地，所述虚拟卫星观测方程为：

本发明的有益效果：可实现北斗卫星授时、北斗定位、GPS授时、GPS定位的功能，以摆脱对GPS的以来，提高我国民用基础行业的安全性和稳定性。

附图说明

图1为本发明***的结构示意图；

图2为本发明的授时示意图；

图3为本发明的BD1+GPS联合授时策略框图。

具体实施方式

以下结合附图和具体实施例，对本发明进行详细说明。

本部分的发明内容是对空射击手持终端中的一种采用迭代差分算法的远程定位***及方法。

定位***是由地面控制中心(地面段)、远程网络(空间段)和大量用户终端(用户段)三部分构成。空间段由两个空间目标及一个空间备用目标组成，其空间目标及备用目标可为三颗地球赤道同步轨道卫星，其中，两颗工作星，一颗备用星，但不限于为卫星，两颗工作星分别定点于东经80°和120°，轨道高度为35，000km。卫星的主要功能是转发地面控制指挥中心和用户机之间的双向信号传输。因而卫星上的主要载荷为变频转发器、两个波束的S频段天线和L频段天线。两颗卫星的4个S频段波束分区覆盖全服务区，每颗卫星的两个L频段波束分区覆盖全服务区。

地面控制中心是整个***的控制、管理及导航定位解算的核心，负责整个***的信息生成、收集、处理及***状况的检测，不仅要完成对卫星及用户的管理，还要完成定位解算及授时等业务。主要包括信号收发子***，信息处理分***，时间分***，监控子***及轨道监控分***等。

根据导航定位的基本原理，建立一个位置已知的空间参考系，测量待定位点到各参考点的相对位置，从而确定待测点在空间参考系中的位置。因而首先需要一个空间参考基准，在远程定位***中即为卫星网络。一般要把空间参考系级待测点放到一定的空间坐标系中，以便于描述它们的空间位置。采用的空间坐标系包括地心惯性坐标系(ECI)和地心地球固连坐标系(ECEF)。因为无论是卫星的位置还是用户的位置都可能随着时间发生变化，因而精确描述卫星位置或用户位置均需要一个精确的时间基准。可提供的时间参考基准包括世界时***、使用原子钟提供的地方原子时及国际原子时、以及协调世界时(UTC)等，各个不同的卫星定位***也有自己的时间参考***，作为整个***的时间同步标准。

定位***可以利用控制中心获得的信号强度(SS:Signal Strength)、载波相位(CP:Carrier Phase)、信号到达角度(AOA:Angle Of Arrival；DOA:Direction OfArrival)和时间测量值(TOA:Time Of Arrival；TDOA:Time Difference Of Arrival)以及它们的组合进行定位估算。在卫星定位***中，多利用时间测量值(可换算为待测点与参考点间的相对距离)实现定位，由于观测点时间无法与***时间完全同步，因此测量值中包含观测点钟差信息，观测量一般称为伪距。通常伪距可以表示为：

r⁽ⁱ⁾待定位点到卫星i的伪距离；待定位点到卫星i的几何距离，T接收机钟差；T⁽ⁱ⁾卫星i的钟差；c光速；信号的空间传输延迟；u⁽ⁱ⁾测量误差及噪声。

由公式可见，如果卫星位置已知的情况下，观测方程中包含观测点位置和钟差在内的四个未知量，卫星钟差传输延迟等都当作测量误差，其中有一些量是可以进行校正的。对于四个未知量的方程组，需要四个观测量才可得解，即需要至少四颗卫星，也就是通常卫星定位***的工作模式。如果不仅高程已知，测量点时间也与***钟同步，知需两个观测量即可，此为北斗一代***中有源定位模式。如果高程信息已知，则只需要三颗卫星，为本文中面临的模式。

如图2所示，北斗***中，信号在地面站生成，上行至空间段，卫星进行变频转发，用户机的测距观测量实际是地面站到卫星再到用户的传输延迟。

用户机观测量可以表示为：

s⁽ⁱ⁾为待定位点的观测量；为地面站到卫星i的几何距离，***通过广播信号上行延迟来描述此项；观测点到卫星i的几何距离；卫星i的转发延迟；t接收机钟差；c光速；信号的空间传输延迟；u⁽ⁱ⁾测量误差及噪声。由于所有卫星信号都来自地面站信号，因此***时钟是同步的。卫星到用户的几何距离可以表示为：

(x_i，y_i，z_i)为卫星i的坐标。在(2.2)和(2.3)中(x_i，y_i，z_i)均可通过***广播信息和本地估计得到。公式(2.2)变为：

(2.4)

当高程已知的情况下，(2.4)可解。

用户位置信息中，高度是最有可能通过高度计、气压计等方式获取的。并且在有些应用中，例如海上用户，高度信息基本不变。高度是经纬高坐标系的描述，而式(2.4)中未知量(x,y,z,t)中的位置信息是ECEF坐标系的描述，两种坐标系不正交，因此高度信息不能直接在公式(2.4)中作为已知量使用。利用高度，可以构造出一个观测方程，与公式(2.4)联立即可得解。

观测方程构造方法有两种：椭球方程、虚拟卫星。

椭球方程。地球可以模拟化成一个椭球，高度如果为零，则用户应该位于地球面上。高度为h时，观测方程可以近似为公式(2.5)。a，b分别为地球椭球模型的长轴和短轴。

构造虚拟卫星。与建立椭球方程在本质上是一样的。这种方法的基本思想为：过观测点作地球表面的垂线，与赤道面产生交点，交点即虚拟卫星位置。对于椭球模型，交点一般并不位于椭球中心。根据虚拟卫星位置和用户位置估计值可得几何距离。有了虚拟卫星的位置及“伪距”值，则可以构造与式(2.4)类似的观测方程，如下：

线性化迭代的方法求解方程(2.5)(2.6)，即可获得定位结果。

用户机将在位置已知和未知两种模式下工作，通过用户界面实现切换。

用户机接收天线相位中心的位置是精确已知的，信号从地面站发出时其含有时间标志与***时是一致的。经过地面站到卫星、卫星到用户机的两段空间延迟及卫星的转发延迟后，到达用户机接收天线，此时信号包含的时间信息相对***时是延后的。地面站到卫星的空间延迟、卫星转发延迟均包含在导航电文中。而卫星到用户机的空间延迟需要利用卫星位置和用户机位置计算得到。因此授时的关键都是要获知用户机的位置。位置已知的模式下，用户机位置通过用户界面置入；位置未知模式下，用户机位置用定位解算得到。

另外，用户机输出的授时结果包含由线延迟、处理延迟等组成的设备零值，在天线馈线、用户机硬件结构不改变的情况下，设备零值为常数，可通过用户界面进行调整。

用户机配备较高精度频标，如恒温晶振或高精度温补晶振，分频产生本地1pps信号。将其作为授时的参考时标，伪距及载波相位等的测量均基于参考时标进行。本地1pps信号与地面控制中心1pps存在钟差t，t校正本地1pps信号，产生授时结果。未知位置授时模式，通过PVT解算，得到本地钟差项，当做校正项t。已知位置授时模式下，可以得到用户机与卫星的几何距离，校正项t包括：地面站到卫星的上行延迟，卫星转发延迟，用户机到卫星的几何距离造成的延迟，下行大气延迟，用户机零值延迟。利用时间延迟电路将本地1pps延迟t，使得通过延迟输出产生的1pps信号与***1pps信号同步对齐。年、月、日、时、分、秒可以从导航电文中得到。未知及已知位置授时原理如图2所示。

(3)BD1+GPS联合授时调度策略

单独使用BD1或GPS进行授时均存在的一定的局限性，信号容易被遮挡、***不可用等因素均可能导致授时中断。对于需要连续授时服务的用户来说，授时中断是不允许的。如果采用BD1+GPS联合的方式，授时可靠性将得到极大的增强。

BD1与GPS既可以在信号处理阶段进行融合，也可以在各自产生最终结果后再进行融合。融合节点越靠前称融合愈紧，节点越靠后称融合愈松。更紧的融合对于软硬件的改动都很大，要在设计之初进行统一考虑；如果两种***下的设备都已经成型，则适合采用松融合的方式。

具体见图3所示，本方法采用的是一种板级联合的松融合方式，已有的BD1无源终端设备与第三方提供的GPS授时设备作为两块背板，统一受控制板的管理。

BD1无源授时终端、GPS授时终端、BD1天线、GPS天线均为已有模块，需要重新设计控制板。控制板完成的功能分为三个部分：

(1)为BD1无源授时终端、GPS授时终端供电。

(2)BD1无源授时终端、GPS授时终端可以通过控制板与外部实现数据交互。

(3)选择BD1无源授时终端、GPS授时终端之中的一路1pps进行输出，选择规则如下：

选取规则1：根据外部设定的对象输出。

选取规则2：当发生异常时，选择非异常的1pps输出。

选取规则3：当两路1pps均正常，且外部没有设定输出对象时，按照默认选择输出。

Claims

1.一种采用迭代差分算法的远程定位***，其特征在于，所述***包括：

控制中心：连接远程网络，同时管理整个***、采用迭代差分算法进行导航定位解算；

用户终端：连接远程网络，向远程网络接收及发送信息；

所述控制中心采集获得用户终端的信号强度、载波相位、信号到达角度和时间观测量，创建待定位点观测量，表示为：

s⁽ⁱ⁾为待定位点的观测量；为控制中心到空间目标i的几何距离；待定位点到空间目标i的几何距离；空间目标i的转发延迟；t接收机钟差；c光速；信号的空间传输延迟；u⁽ⁱ⁾测量误差及噪声。

2.根据权利要求1所述的远程定位***，其特征在于，所述远程网络包括两个空间目标及一个空间备用目标。

3.根据权利要求2所述的远程定位***，其特征在于，所述两个空间目标及一个空间备用目标均载荷有变频转发设备、两个波束的S频段天线及L频段天线。

4.根据权利要求1所述的远程定位***，其特征在于，所述控制中心包括信号收发子***、信息处理分***、时间分***、监控子***及轨道监控分***。

5.一种采用迭代差分算法的远程定位方法，应用权利要求1-4之一的***，其特征在于，包括以下步骤：

a)建立一个位置已知的空间参考系，测量待定位点到各参考点的相对位置，其中控制中心及远程网络的工作单元为参考点，用户终端为待定位点；

s⁽ⁱ⁾为待定位点的观测量；为控制中心到空间目标i的几何距离；待定位点到空间目标i的几何距离；空间目标i的转发延迟；t接收机钟差；c光速；信号的空间传输延迟；u⁽ⁱ⁾测量误差及噪声；

c)获得空间目标到待定位点的几何距离为：

(x_i，y_i，z_i)为空间目标i的坐标；

d)由待定位点观测量s⁽ⁱ⁾和通过***广播信息和本地估计得到的(x_i，y_i，z_i)，将步骤b)中的公式变换得到：

；r⁽ⁱ⁾为待定位点到空间目标i的伪距离；

6.根据权利要求5所述的定位方法，其特征在于，所述空间参考系包括地心惯性坐标系和地心地球固连坐标系。

7.根据权利要求5所述的方法，其特征在于，当e)步骤中无法通过气压计及高度计测量得到高度信息时，可构造方程，包括椭球观测方程及虚拟卫星观测方程，并通过线性迭代的方法求解，获得定位结果。

8.根据权利要求7所述的方法，其特征在于，所述椭球观测方程为：

其中高度为h；a，b分别为地球椭球形的长轴和短轴。

9.根据权利要求7所述的方法，其特征在于，所述虚拟卫星观测方程为：

<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msup> <mi>r</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>+</mo> <mi>t</mi> <mo>*</mo> <mi>c</mi> <mo>+</mo> <msup> <mi>u</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mi>r</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>+</mo> <mi>t</mi> <mo>*</mo> <mi>c</mi> <mo>+</mo> <msup> <mi>u</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mi>r</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </msup> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>+</mo> <mi>t</mi> <mo>*</mo> <mi>c</mi> <mo>+</mo> <msup> <mi>u</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mrow> <mi>v</mi> <mi>s</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>x</mi> <mrow> <mi>v</mi> <mi>s</mi> </mrow> </msub> <mo>-</mo> <mi>x</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>y</mi> <mrow> <mi>v</mi> <mi>s</mi> </mrow> </msub> <mo>-</mo> <mi>y</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>z</mi> <mrow> <mi>v</mi> <mi>s</mi> </mrow> </msub> <mo>-</mo> <mi>z</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>+</mo> <mn>0</mn> <mo>*</mo> <mi>c</mi> <mo>+</mo> <msup> <mi>u</mi> <mo>&prime;</mo> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow>