一种导航卫星星间Ka波段双向测量伪距的预处理方法与流程

本发明属于航天测量与控制技术领域，涉及一种导航卫星星间ka波段双向测量伪距的预处理方法。

背景技术：

导航星座采用时分多址星间测量和通信技术，把时间分割成互不重叠的时段(帧)，再将帧分割成互不重叠的时隙，与连接对象一一对应，依据时隙区分来自不同地址的信号，从而完成多址的连接。假设导航星座由24颗卫星组成，每颗导航卫星分配一个时隙，每个时隙维持1.5s，在轨工作的24颗导航卫星均分配一个不同的时隙，36s时间可以对星座所有卫星轮询一遍，该时间定义为一帧。在测距帧期间，任意一颗导航卫星在分配的时隙里发射测距信号，其他所有可见的导航卫星进行接收，完成伪距的测量。

导航卫星利用星载ka波段无线电天线可以建立星间链路(intersatellitelink，isl)，实现双向距离测量，用于支持卫星轨道确定。相比传统的地面测距技术，该技术具有以下优势：一是可实现卫星全弧段连续跟踪测量，不受地面站可见弧段的约束；二是不依赖地面测控站支持，节省了地面测控资源，提高了卫星自主生存的能力。

另外，导航卫星采用时分多址(tdma)体制与不同的卫星建立星间链路，同一条链路的双向测量实际上包含同一个测距帧内完成的两次单向测量。具体的测量过程描述如下：

如图2所示，卫星a在其卫星时ta,1时刻(对应系统时间t1)发出信号，经发射系统的电路延迟ca,1、无线电的空间延迟τab、卫星b的电路处理延迟cb,2，在卫星b的卫星时tb,2时刻(对应系统时间t2)被检测到。于是信号从发出到接收的实际总延时为t2-t1，而卫星b测量到的时间差tab为tb,2-ta,1。空间延迟τab与两卫星相对距离相关，是时间的函数，考虑到电路延迟都是小量，因此τab可记为t1、t2时刻的函数τab(t1,t2)，该延迟包括了路径延迟和电离层延迟等所有空间路径上的延迟。

另一方面，卫星b也在钟面时tb,3(系统时t3)时刻发出测距信号，经发射系统的电路延迟cb,3、无线电的空间延迟τba、卫星a的电路处理延迟ca,4，在卫星a的钟面时ta,4(系统时t4)时刻被检测到。同理，卫星a测到的时间差tba为ta,4-tb,3，τba则记为τba(t3,t4)。

上述过程的基本测量量为

tab＝tb,2-ta,1＝(tb,2-tb,1)+(tb,1-ta,1)(10)

tba＝ta,4-tb,3＝(ta,4-ta,3)+(ta,3-tb,3)(11)

注意到卫星时等于系统时加上钟差：ta＝t+δta(t)，tb＝t+δtb(t)，则有

tab＝[δtb(t2)-δta(t1)]+ca,1+τab(t1,t2)+cb,2(12)

tba＝[δta(t4)-δtb(t3)]+cb,3+τba(t3,t4)+ca,4(13)

上述ka星间链路测量时延tab、tba乘以光速包含两颗卫星的位置信息和钟差信息，得到星间双向测量伪距

ρab＝c·{[δtb(t2)-δta(t1)]+ca,1+τab(t1,t2)+cb,2}(14)

ρba＝c·{[δta(t4)-δtb(t3)]+cb,3+τba(t3,t4)+ca,4}(15)

式中，c表示真空中的光速。

在实际测控应用过程中，上述双向测量伪距主要用于计算卫星轨道，传统的处理方法有两种：一是建立卫星轨道和钟差的联合估计方程，同时估计轨道和钟差参数，这种处理过程较为复杂；二是引入独立的钟差信息，对观测数据中的卫星钟差进行消除，这需要引入其它精确的钟差数据源。

技术实现要素：

本发明的目的是提供一种导航卫星星间ka波段双向测量伪距的预处理方法，能够生成对应于两颗卫星质心位置的测距观测数据，用于支持后续轨道计算。

本发明所采用的技术方案是，一种导航卫星星间ka波段双向测量伪距预处理方法，具体按照下述步骤进行：

步骤1，在遥测数据中提取卫星的星间测量数据及接收端时标并解码，读取星间测量数据及接收端时标，并对星间测量数据中每个单向测量数据进行时间排序、重点剔重后，将单向测量数据进行双向数据配对，得到多对星间双向测量值；

步骤2，将每对星间双向测量值进行时标统一改正；

步骤3，将每对时标统一改正后的星间双向测量值进行消钟差组合，得到双向距离观测值；

步骤4，将双向距离观测值进行天线安装偏差改正，得到对应于卫星质心的双向距离观测值；

步骤5，利用地面标校值对对应于卫星质心的双向距离观测值进行星载测量天线的通道延时改正，得到卫星的星间距。

本发明的特点还在于：

根据双单向测量原理，所述步骤1中的星间双向测量值包括：

ρab＝c·{[δtb(t2)-δta(t1)]+ca,1+τab(t1,t2)+cb,2}(1)

其中，ρab为卫星a发b收的单向伪距，c为真空中的光速，t2为卫星b接收到测距信号时对应的系统时间，δtb(t2)为系统时t2时刻卫星b的钟差，t1为卫星a发出测距信号时对应的系统时间，δta(t1)为系统时t1时刻卫星a的钟差，ca,1为卫星a发射通道时延，τab(t1,t2)为卫星a、b空间路径传播时延，cb,2为卫星b接收通道时延；

ρba＝c·{[δta(t4)-δtb(t3)]+cb,3+τba(t3,t4)+ca,4}(2)

其中，ρba为卫星b发a收的单向伪距，t4为卫星a接收到测距信号时对应的系统时间，δta(t4)为t4时刻卫星b的钟差，t3为卫星b发出测距信号时对应的系统时间，δtb(t3)为系统时t3时刻卫星b的钟差，cb,3为卫星b发射通道时延，τba(t3,t4)为卫星a、b空间路径传播时延，ca,4为卫星a接收通道时延。

步骤2中按照下述步骤将每对星间双向测量值进行时标统一改正：

步骤2.1，将星间双向测量值归算至同一测距历元t0：

其中，a1,a为卫星a的卫星钟频率偏差，a1,b为卫星b的卫星钟频率偏差，xb(t2)为系统时t2时刻卫星b的位置，xa(t2)为系统时t2时刻卫星a的位置，xb(t4)为系统时t4时刻卫星b的位置，xa(t4)为系统时t4时刻卫星a的位置，xb(t0)为系统时t0时刻卫星b的位置，xa(t0)为系统时t0时刻卫星a的位置；

步骤2.2，归算至同一测距历元t0的星间双向测量值为：

其中，表示考虑光行时的发射端时刻，δtb(t0)为系统时t0时刻卫星b的钟差，δta(t0)为系统时t0时刻卫星a的钟差，为系统时t0时刻卫星a、b之间的空间路径传播时延。

步骤3中具体按照下述方法将时标统一改正后的星间双向测量值进行消钟差组合：

将时标统一改正后的星间双向测量值相加求平均得到双向距离观测值：

步骤4中按照下述步骤将所述双向距离观测值进行天线安装偏差改正，得到对应于卫星质心的双向距离观测值：

步骤4.1，计算将双向距离观测值改正到卫星质心的改正量：

其中，为卫星a本体坐标系到惯性坐标系的转换矩阵，δra,pco为卫星a天线安装位置相对于质心偏差，为卫星b本体坐标系到惯性坐标系的转换矩阵，δrb,pco为卫星b天线安装位置相对于质心偏差，e表示惯性系中卫星a到卫星b的方向矢量；

步骤4.2，根据所述改进量计算对应于卫星质心的双向距离观测值：

步骤5中通过利用卫星方提供的通道时延参数作为初始值的方法对对应于卫星质心的双向距离观测值进行通道延时改正，得到卫星的星间距。

本发明的有益效果是

本发明一种导航卫星星间ka波段双向测量伪距的预处理方法，操作方法简单，能够更精确地测量两颗卫星质心位置的测距观测数据，测量得到的颗卫星质心位置的距离观测数据文件能够直接用于完成卫星轨道的确定。

附图说明

图1是本发明一种导航卫星星间ka波段双向测量伪距的预处理方法的流程图；

图2是本发明一种导航卫星星间ka波段双向测量伪距的预处理方法背景技术中计算双向测量伪距的示意图。

具体实施方式

下面结合附图和具体实施方式对本发明进行详细说明。

一种导航卫星星间ka波段双向测量伪距的预处理方法，如图1所示，具体按照下述步骤进行：

步骤1，在遥测数据中提取卫星的星间测量数据及接收端时标并解码，读取星间测量数据及接收端时标，并对星间测量数据中每个单向测量数据进行时间排序，然后将距离不在[10000km，70000km]区间内的数据进行重点剔冲后，将±3s时隙内的单向测量数据进行双向数据配对，得到多对星间双向测量值：

根据双单向测量原理，步骤1中的星间双向测量值包括：

ρab＝c·{[δtb(t2)-δta(t1)]+ca,1+τab(t1,t2)+cb,2}(1)

其中，ρab为卫星a发b收的单向伪距，c为真空中的光速，t2为卫星b接收到测距信号时对应的系统时间，δtb(t2)为系统时t2时刻卫星b的钟差，t1为卫星a发出测距信号时对应的系统时间，δta(t1)为系统时t1时刻卫星a的钟差，ca,1为卫星a发射通道时延，τab(t1,t2)为卫星a、b空间路径传播时延，cb,2为卫星b接收通道时延；

ρba＝c·{[δta(t4)-δtb(t3)]+cb,3+τba(t3,t4)+ca,4}(2)

其中，ρba为卫星b发a收的单向伪距，t4为卫星a接收到测距信号时对应的系统时间，δta(t4)为t4时刻卫星b的钟差，t3为卫星b发出测距信号时对应的系统时间，δtb(t3)为系统时t3时刻卫星b的钟差，cb,3为卫星b发射通道时延，τba(t3,t4)为卫星a、b空间路径传播时延，ca,4为卫星a接收通道时延。

步骤2，将每对星间双向测量值进行时标统一改正，聚天下好下述步骤进行：

将双单向测量值ρab(t2)和ρba(t4)(时标t4)的时标算到同一测距帧历元t0，其中，t0与t2、t4归间隔不大于3s：

其中，a1,a为卫星a的卫星钟频率偏差，a1,b为卫星b的卫星钟频率偏差，xb(t2)为系统时t2时刻卫星b的位置，xa(t2)为系统时t2时刻卫星a的位置，xb(t4)为系统时t4时刻卫星b的位置，xa(t4)为系统时t4时刻卫星a的位置，xb(t0)为系统时t0时刻卫星b的位置，xa(t0)为系统时t0时刻卫星a的位置；

步骤2.2，归算至同一测距历元t0的星间双向测量值为：

其中，表示考虑光行时的发射端时刻，δtb(t0)为系统时t0时刻卫星b的钟差，δta(t0)为系统时t0时刻卫星a的钟差，为系统时t0时刻卫星a、b之间的空间路径传播时延。

步骤3，将每对时标统一改正后的星间双向测量值进行消钟差组合，得到双向距离观测值，具体按照下述方法进行：

将时标统一改正后的星间双向测量值相加求平均得到双向距离观测值：

步骤4，将所述双向距离观测值进行天线安装偏差改正，得到对应于卫星质心的双向距离观测值，具体按照下述步骤进行：

步骤4.1，计算将双向距离观测值改正到卫星质心的改正量：

其中，为卫星a本体坐标系到惯性坐标系的转换矩阵，δra,pco为卫星a天线安装位置相对于质心偏差，为卫星b本体坐标系到惯性坐标系的转换矩阵，δrb,pco为卫星b天线安装位置相对于质心偏差，e表示惯性系中卫星a到卫星b的方向矢量；

步骤4.2，根据所述改进量计算对应于卫星质心的双向距离观测值：

步骤5，通过利用卫星方提供的通道时延参数作为初始值的方法对对应于卫星质心的双向距离观测值进行通道延时改正，得到卫星的星间距：

对应于卫星质心的双向距离观测值中还包含了a、b两颗卫星天线的收发通道时延，需要进一步修正，对于星载天线设备来说，它的通道时延具备一定的稳定性，因此在一次轨道计算过程中，可以将收发通道时延作为常量进行改正。对于每颗卫星，建立一个收发通道时延和参数估计模型，在轨道计算过程中与轨道参数一起估计。