低轨卫星TDMA静中通系统的前向链路时延估计方法与流程

低轨卫星tdma静中通系统的前向链路时延估计方法
技术领域
1.本发明涉及tdma卫星通信系统的时间同步技术，尤其涉及一种透明转发低轨卫星tdma静中通系统的前向链路时延估计方法。


背景技术：

2.tdma系统的全网时间同步包含两个部分：前向链路时间同步和返向链路时间同步。前向链路时间同步的作用是实现网络时钟参考(ncr)同步，即建立端站计时系统与主站计时系统之间的映射关系，进而达到端站通过本地计时系统识别主站时隙位置的目的。
3.前向链路传输时延一般包含：前向链路的发送侧处理时延、前向链路的接收侧处理时延、星上处理时延和前向链路的自由空间传输时延。对于透明转发卫星而言，星载转发器仅对信号进行滤波、放大、变频等简单处理，其时延远小于前向链路的自由空间传输时延，可以忽略；前向链路的发送侧处理时延和接收侧处理时延均可以根据具体的处理机制准确地估计。因此，本技术主要讨论前向链路自由空间传输时延的估计方法。
4.透明转发卫星tdma通信系统中的前向链路是一条从主站到卫星，再到端站的传输链路，由前向上行链路和前向下行链路两部分组成，如图1所示。其中，前向上行链路是指从主站到卫星的传输链路，前向下行链路是指从卫星到端站的传输链路。通过分别估计出前向上行链路和前向下行链路的自由空间传输时延，并将二者相加，即可得到前向链路的自由空间传输时延。
5.在低轨卫星通信系统中，由于卫星的位置是随时间变化的，所以前向链路不是一条固定不变的传输链路，它是随时间变化而变化的。准确估计前向链路自由空间传输时延的首要前提是正确找到待估计的前向链路是哪一条。
6.在静中通场景中，主站和端站的位置固定不变，前向链路由卫星的位置唯一确定；而卫星的位置是时间的函数，若是确定了前向信息到达卫星的时间，就能确定相应时刻卫星的位置，进而唯一确定传输前向信息的前向链路。在如图1所示的场景中，前向信息d[n]在t
n
时刻从主站发出，在t
′
n
时刻到达卫星，随后被卫星转发给端站。在t
′
n
时刻，卫星的位置是p(t
′
n
)，传输d[n]的前向链路就是h
→
p(t
′
n
)
→
r。
[0007]
因此，确定前向信息到达卫星的时间是准确估计前向链路自由空间传输时延的关键。


技术实现要素：

[0008]
为了解决上述相关现有技术问题，本发明提供一种透明转发低轨卫星tdma静中通系统的前向链路时延估计方法，该方法采用了基于线性近似的外推估计法，计算简便；在卫星星历信息准确的条件下，可以获得良好的估计性能。
[0009]
为了实现本发明的目的，本发明拟通过以下技术方案实现：
[0010]
一种透明转发低轨卫星tdma静中通系统的前向链路时延估计方法，包括步骤：
[0011]
s1、以星历时间t为横坐标，并以发送前向信息d[n]的星历时间t
n
为横坐标原点，
以卫星/主站间的传输时延τ
′
sat/hub
(t)为纵坐标，并以0为纵坐标原点，建立直角坐标系
[0012]
s2、在直角坐标系中，建立两个线性方程：一个是恒等直线方程l1:τ
′
sat/hub
(t)＝t，另一个是卫星/主站传输时延与星历时间的函数τ
sat/hub
(t)＝g1(t)的坐标平移版本τ
′
sat/hub
(t)＝τ
sat/hub
(t t
n
)＝g1(t t
n
)在区间[0,t]上的局部近似直线段方程t是一个时间增量，且t是一个时间增量，且表示前向上行链路自由空间传输时延的最大值；
[0013]
其中，直线段方程l2的建立方法是：
[0014]
在曲线τ
′
sat/hub
(t)＝g1(t t
n
),t∈[0,t]上选取两个相邻的点(0,g1(t
n
))和(t,g1(t
μ
))，其中，t
μ
＝t
n
t；
[0015]
根据(0,g1(t
n
))和(t,g1(t
μ
))两点的坐标建立直线段方程，得
[0016][0017]
s3、求解直线l1和直线段l2的交点坐标，得
[0018][0019]
s4、估算前向信息d[n]到达卫星的星历时间，得
[0020]
t
′
n
＝t
n
t
intersect
；
[0021]
s5、根据预先计算的两组参数：卫星/主站传输时延与星历时间的函数τ
sat/hub
(t)和卫星/端站传输时延与星历时间的函数τ
sat/rcst
(t)，以及前向信息d[n]到达卫星的星历时间t'
n
，分别估算出前向上行链路的自由空间传输时延τ
fl/ul
(t
′
n
)和前向下行链路的自由空间传输时延τ
fl/dl
(t
′
n
)：
[0022][0023]
s6、将前向上行链路和前向下行链路的自由空间传输时延相加，得到前向链路的自由空间传输时延τ
fl/fs
(t
′
n
)，即
[0024]
τ
fl/fs
(t
′
n
)＝τ
fl/ul
(t
′
n
) τ
fl/dl
(t
′
n
)。
[0025]
本发明有益效果在于：
[0026]
提供了一种透明转发低轨卫星tdma静中通系统的前向链路时延估计方法，该方法采用了基于线性近似的外推估计法，计算简便；在卫星星历信息准确的条件下，可以获得良好的估计性能。
附图说明
[0027]
图1是透明转发低轨卫星tdma通信系统的前向链路示意图。
[0028]
图2是本技术中估算前向信息d[n]到达卫星的星历时间t
′
n
的方法的原理示意图。
[0029]
图3是本技术所述方法中步骤s1～s3的几何示意图。
具体实施方式
[0030]
为了使本技术的目的、技术方案和具体实施方法更为清楚，结合附图实例对本技术进行进一步详细说明。
[0031]
本技术实施例提供一种透明转发低轨卫星tdma静中通系统的前向链路时延估计方法，该方法基于两个前提条件：
[0032]
1)在主站侧，建立了ncr时间与星历时间的映射关系。
[0033]
2)主站将前向信息d[n]及其发送时刻(此处的发送时刻特指信息d[n]离开主站发射天线的时刻)对应的星历时间t
n
一并发送给端站。
[0034]
该方法的设计思路如下：
[0035]
首先，根据卫星的星历信息以及主站和端站的gnss位置信息预先计算出卫星可视窗口内，卫星/主站传输时延与星历时间的函数τ
sat/hub
(t)＝g1(t)和卫星/端站传输时延与星历时间的函数τ
sat/rcst
(t)＝g2(t)。
[0036]
然后，根据主站发送前向信息d[n]的星历时间t
n
(已知量)和函数τ
sat/hub
(t)＝g1(t)，估算出前向信息d[n]到达卫星的星历时间t
′
n
。
[0037]
然后，通过函数τ
sat/hub
(t)＝g1(t)和τ
sat/rcst
(t)＝g2(t)分别估算出前向上行链路的自由空间传输时延τ
fl/ul
(t
′
n
)＝τ
sat/hub
(t
′
n
)和前向下行链路的自由空间传输时延τ
fl/dl
(t
′
n
)＝τ
sat/rcst
(t
′
n
)。
[0038]
最后，将前向上行链路和前向下行链路的自由空间传输时延相加，得到前向链路的自由空间传输时延。
[0039]
为了简化前向链路时延估计的计算复杂度，需要预先准备两组参数：卫星/主站传输时延与星历时间的函数和卫星/端站传输时延与星历时间的函数。
[0040]
计算卫星/主站传输时延与星历时间的函数的方法如下：
[0041]
1)根据卫星的星历信息计算出卫星可视窗口内，卫星在地心地固(ecef)坐标系下的位置与星历时间的函数。
[0042]
2)根据主站的gnss位置换算得到主站在ecef坐标系下的位置。
[0043]
3)在ecef坐标系下，根据卫星和主站的位置信息，计算出卫星/主站距离与星历时间的函数d
sat/hub
(t)＝f1(t)。
[0044]
4)用d
sat/hub
(t)＝f1(t)除以电磁波的传播速度c，得到卫星/主站传输时延与星历时间的函数τ
sat/hub
(t)＝d
sat/hub
(t)/c。
[0045]
采用类似的方法，可以计算出卫星/端站传输时延与星历时间的函数τ
sat/rcst
(t)，其中τ
sat/rcst
(t)＝d
sat/rcst
(t)/c，d
sat/rcst
(t)＝f2(t)为卫星/端站距离与星历时间的函数。
[0046]
在上述条件下，进行前向链路的自由空间传输时延估计：
[0047]
s1、以星历时间t为横坐标，并以发送前向信息d[n]的星历时间t
n
为横坐标原点，以卫星/主站间的传输时延τ
′
sat/hub
(t)为纵坐标，并以0为纵坐标原点，建立直角坐标系
[0048]
s2、在直角坐标系中，建立两个线性方程。
[0049]
一个是恒等直线方程l1:τ
′
sat/hub
(t)＝t；
[0050]
另一个是卫星/主站传输时延与星历时间的函数τ
sat/hub
(t)＝g1(t)的坐标平移版
本τ
′
sat/hub
(t)＝τ
sat/hub
(t t
n
)＝g1(t t
n
)在区间[0,t]上的局部近似直线段方程
[0051]
后者的建立方法如下：首先，在曲线τ
′
sat/hub
(t)＝g1(t t
n
),t∈[0,t]上选取两个相邻的点(0,g1(t
n
))和(t,g1(t
μ
))，其中，t
μ
＝t
n
t，t是一个时间增量，且 t，t是一个时间增量，且表示前向上行链路自由空间传输时延的最大值。然后，根据上述两点的坐标建立直线段方程，得
[0052][0053]
s3、求解直线l1和直线段l2的交点坐标，得
[0054][0055]
如图3所示，是本方法中步骤s1～s3的几何示意图，其代表的场景是卫星与主站的距离逐渐减小，因此前向上行链路的传输时延也逐渐减小。
[0056]
s4、估算前向信息d[n]到达卫星的星历时间，得
[0057]
t
′
n
＝t
n
t
intersect
。
[0058]
估算前向信息d[n]到达卫星的星历时间t
′
n
的方法如下：如图1所示，在t
n
时刻，主站发出前向信息d[n]，此时，卫星的位置是p(t
n
)。之后，d[n]在前向上行链路中传输，而卫星继续沿着箭头指向的方向在轨道上运动。在t
′
n
时刻，卫星运动到新的位置p(t
′
n
)，而此时，d[n]也刚好到达卫星，即经过时延t
′
n
‑
t
n
后d[n]与卫星相遇。在上述相遇问题中，卫星在轨道p(t
n
)
→
p(t
′
n
)上的运动时间与d[n]在前向上行链路h
→
p(t
′
n
)上的传输时延恰好相等，即直线τ
sat/hub
(t)＝t
‑
t
n
与曲线τ
sat/hub
(t)＝g1(t)在t
n
之后必然相交，且交点就是t
′
n
(如图2所示)。因此，求解非线性方程组
[0059][0060]
就可以得到前向信息d[n]到达卫星的星历时间t
′
n
。
[0061]
通常，在低轨卫星通信系统中，前向上行链路的自由空间传输时延很小，卫星在这段时间的移动距离很短，其运动轨迹可以近似为直线段。因此，曲线τ
sat/hub
(t)＝g1(t)在t
n
附近的区域也可以近似为直线段。令此直线段的方程为其中，t是一个时间增量。进而，上述问题可以简化为直线τ
sat/hub
(t)＝t
‑
t
n
和直线段的相交问题。综述所述，求解线性方程组
[0062][0063]
就可以得到前向信息d[n]到达卫星的星历时间t
′
n
。
[0064]
s5、根据下式
[0065][0066]
分别估算出前向上行链路的自由空间传输时延τ
fl/ul
(t
′
n
)和前向下行链路的自由空间传输时延τ
fl/dl
(t
′
n
)。
[0067]
s6、将前向上行链路和前向下行链路的自由空间传输时延相加，得到前向链路的自由空间传输时延τ
fl/fs
(t
′
n
)，即
[0068]
τ
fl/fs
(t
′
n
)＝τ
fl/ul
(t
′
n
) τ
fl/dl
(t
′
n
)。