低轨卫星TDMA动中通系统的返向突发发送时间估计方法与流程

专利检索2022-05-10  37


低轨卫星tdma动中通系统的返向突发发送时间估计方法
技术领域
1.本发明属于数字通信技术领域,涉及tdma卫星通信系统的时间同步技术,尤其涉及一种透明转发低轨卫星tdma动中通系统的返向突发发送时间估计方法。


背景技术:

2.tdma系统是一种基于时分多址技术的多用户通信系统。在tdma卫星通信系统中,主站对返向链路的各种时隙进行规划,并为每个入网的端站分配专用时隙进行数据传输。不同端站则根据时隙规划分时共享同一个物理信道。因此,所有端站都需要建立与主站一致的计时系统,以便于准确识别主站规划的时隙及其范围,以及恰当选择返向突发的发送时间,以使返向突发在期望的时间达到期望的时隙。这就是tdma系统的全网时间同步,其包含两个部分:前向链路时间同步和返向链路时间同步。前向链路时间同步的作用是实现网络时钟参考(ncr)同步,即建立端站计时系统与主站计时系统之间的映射关系,进而达到端站通过本地计时系统识别主站时隙位置的目的。而返向链路时间同步的作用是实现返向突发的到达时间(toa)同步,即端站准确估计出返向突发的发送时间,并在此时间发送返向突发,以使其准确地落入目标时隙。
3.如图1所示,估算返向突发发送时间的一般方法是:首先,估算出传输返向突发bf[n]的返向链路的传输时延τ
rl
[n]。然后,用返向突发bf[n]到达主站目标时隙的期望时间t
n
(已知量)减去返向链路的传输时延τ
rl
[n],得到返向突发的发送时间t

n

[0004]
在低轨卫星通信系统中,由于卫星的位置是随时间变化的,所以返向链路不是一条固定不变的传输链路,它是随时间变化而变化的。准确估计返向链路传输时延的首要前提是正确找到待估计的返向链路是哪一条。对于动中通场景,主站的位置固定不变,而卫星和端站的位置随时间变化,所以返向链路由卫星和端站的位置唯一确定。如图2所示,返向突发bf[n]在t

n
时刻从端站发出,在t

n
时刻达到卫星,随后被卫星转发给主站,并在t
n
时刻到达主站。在t

n
时刻,端站的位置是r(t

n
);在t

n
时刻,卫星的位置是p(t

n
)。因此,传输bf[n]的返向链路就是r(t

n
)

p(t

n
)

h。由于动中通端站处于运动状态,所以,在t

n
之前,不可能准确知道端站的位置r(t

n
),进而,无法准确估计出返向链路的传输时延τ
rl
[n]。因此,采用上述方法难以准确估计出返向突发的发送时间。


技术实现要素:

[0005]
为了解决上述相关现有技术问题,本发明提供一种透明转发低轨卫星tdma动中通系统的返向突发发送时间估计方法,无需估算返向链路的传输时延便可直接估计出返向突发的发送时间,准确度高且计算简便。
[0006]
为了实现本发明的目的,本发明拟通过以下技术方案实现:
[0007]
一种透明转发低轨卫星tdma动中通系统的返向突发发送时间估计方法,包括以下步骤:
[0008]
s1、建立直角坐标系以星历时间t为横坐标,并以返向突发bf
[n]到达主站目标时隙的期望星历时间t
n
为横坐标原点;以卫星/主站间的传输时延τ

sat/hub
(t)为纵坐标,并以0为纵坐标原点。
[0009]
s2、在直角坐标系中,建立两个线性方程。一个是反恒等直线方程l1:τ

sat/hub
(t)=

t,另一个是卫星/主站传输时延与星历时间的函数τ
sat/hub
(t)=g1(t)的坐标平移版本τ

sat/hub
(t)=τ
sat/hub
(t t
n
)=g1(t t
n
)在区间[

t,0]上的局部近似直线段方程l2:后者的建立方法是:首先,在曲线τ

sat/hub
(t)=g1(t t
n
),t∈[

t,0]上选取两个相邻的点(0,g1(t
n
))和(

t,g1(t
μ
)),其中,t
μ
=t
n

t,t是一个时间增量,且个时间增量,且表示返向下行链路自由空间传输时延的最大值。然后,根据上述两点的坐标建立直线段方程,得
[0010][0011]
s3、求解直线l1和直线段l2的交点坐标,得
[0012][0013]
s4、估算返向突发bf[n]到达卫星的星历时间,得
[0014]
t

n
=t
n
t
intersect

[0015]
s5、建立直角坐标系以星历时间t为横坐标,并以返向突发bf[n]到达卫星的星历时间t

n
为横坐标原点;以卫星/端站间的传输时延τ

sat/rcst
(t|t

n
)为纵坐标,并以0为纵坐标原点。
[0016]
s6、在直角坐标系中,建立两个线性方程。一个是反恒等直线方程l3:τ

sat/rcst
(t|t

n
)=

t,另一个是卫星/端站传输时延与星历时间的函数τ
sat/rcst
(t|t

n
)=g2(t)的坐标平移版本τ

sat/rcst
(t|t

n
)=τ
sat/rcst
(t t

n
|t

n
)=g2(t t

n
)在区间[

t
α
,0]上的局部近似直线段方程l4:后者的建立方法是:首先,在曲线τ

sat/rcst
(t|t

n
)=g2(t t

n
),t∈[

t
α
,0]上选取两个点(

t
α
,g2(t
α
))和(

t
β
,g2(t
β
)),其中,t
α
=t

n

t
α
,t
β
=t

n

t
β
,t
α
和t
β
是两个时间增量,且是两个时间增量,且表示返向上行链路自由空间传输时延的最大值。然后,根据上述两点的坐标建立直线段方程,得
[0017][0018]
s7、求解直线l3和直线段l4的交点坐标,得
[0019][0020]
s8、估算发送返向突发bf[n]的星历时间,得
[0021]
t

n
=t

n
t

intersect

[0022]
s9、根据端站侧星历时间与ncr时间的映射关系,将发送返向突发bf[n]的星历时间t

n
映射成端站本地的ncr时间。
[0023]
本发明有益效果在于:
[0024]
无需估算返向链路的传输时延便可直接估计出返向突发的发送时间,准确度高且计算简便。
附图说明
[0025]
图1是透明转发低轨卫星tdma通信系统中,返向突发发送时间、返向突发到达时间和返向链路传输时延三者关系的示意图。
[0026]
图2是透明转发低轨卫星tdma通信系统的返向链路示意图。
[0027]
图3是本发明中估算返向突发bf[n]到达卫星的星历时间t

n
的方法的原理示意图。
[0028]
图4是本发明的方法中步骤s1~s3的几何示意图。
[0029]
图5是本发明中估算发送返向突发bf[n]的星历时间t

n
的方法的原理示意图。
[0030]
图6是本发明的方法中步骤s5~s7的几何示意图。
具体实施方式
[0031]
为了使本技术的目的、技术方案和具体实施方法更为清楚,结合附图实例对本技术进行进一步详细说明。
[0032]
如图1所示为透明转发低轨卫星tdma通信系统中,返向突发发送时间、返向突发到达时间和返向链路传输时延三者关系的示意图。返向突发的发送时间等于返向突发到达主站目标时隙的期望时间减去返向链路的传输时延。
[0033]
如图2所示为透明转发低轨卫星tdma通信系统的返向链路示意图。返向链路是一条从端站到卫星,再到主站的传输链路,由返向上行链路和返向下行链路两部分组成。其中,返向上行链路是指从端站到卫星的传输链路,返向下行链路是指从卫星到主站的传输链路。由于低轨卫星的位置是随时间变化的,所以返向链路不是一条固定不变的传输链路,它是随时间变化而变化的。对于动中通场景,主站的位置固定不变,而卫星和端站的位置随时间变化,所以返向链路由卫星和端站的位置唯一确定。
[0034]
本技术实施例提供一种透明转发低轨卫星tdma动中通系统的返向突发发送时间估计方法,该方法的设计思路如下:
[0035]
首先,根据卫星的星历信息和主站的gnss位置信息预先计算出卫星可视窗口内,卫星/主站传输时延与星历时间的函数τ
sat/hub
(t)=g1(t)。
[0036]
然后,根据返向突发bf[n]到达主站目标时隙的期望星历时间t
n
(已知量)和函数τ
sat/hub
(t)=g1(t),估算出返向突发bf[n]到达卫星的星历时间t

n

[0037]
然后,估算出端站在时间段t∈[t

n

t
α
,t

n
]内的移动过程中,返向上行链路r(t)

p(t

n
)的传输时延与星历时间的函数τ
sat/rcst
(t|t

n
)=g2(t)的线性近似函数其中t
α
>0是一个时间增量,形如f(t|t

n
)的函数是关于t

n
的条件函数,表示在参数t

n
确定的条件下,建立的关于变量t的函数。
[0038]
然后,根据返向突发bf[n]到达卫星的星历时间t

n
和函数
估算出发送返向突发bf[n]的星历时间t

n

[0039]
最后,根据端站侧星历时间与ncr时间的映射关系,将发送返向突发bf[n]的星历时间映射成端站本地的ncr时间。
[0040]
首先,为使用本发明所述方法准备三个前提条件:
[0041]
1、在主站侧,建立ncr时间与星历时间的映射关系。
[0042]
2、端站完成前向链路时间同步,即ncr同步。
[0043]
3、端站能够实时获取到自身的gnss位置信息。
[0044]
在上述前提条件下,为了简化返向突发发送时间估计的计算复杂度,预先准备一组参数,即卫星/主站传输时延与星历时间的函数,其具体计算方法如下:
[0045]
1)根据卫星的星历信息计算出卫星可视窗口内,卫星在地心地固(ecef)坐标系下的位置与星历时间的函数。
[0046]
2)根据主站的gnss位置换算得到主站在ecef坐标系下的位置。
[0047]
3)在ecef坐标系下,根据卫星和主站的位置信息,计算出卫星/主站距离与星历时间的函数d
sat/hub
(t)=f1(t)。
[0048]
4)用d
sat/hub
(t)=f1(t)除以电磁波的传播速度c,得到卫星/主站传输时延与星历时间的函数τ
sat/hub
(t)=d
sat/hub
(t)/c。
[0049]
在上述条件下,进行返向突发发送时间估计:
[0050]
s1、建立直角坐标系以星历时间t为横坐标,并以返向突发bf[n]到达主站目标时隙的期望星历时间t
n
为横坐标原点;以卫星/主站间的传输时延τ

sat/hub
(t)为纵坐标,并以0为纵坐标原点。
[0051]
s2、在直角坐标系中,建立两个线性方程。一个是反恒等直线方程l1:τ

sat/hub
(t)=

t,另一个是卫星/主站传输时延与星历时间的函数τ
sat/hub
(t)=g1(t)的坐标平移版本τ

sat/hub
(t)=τ
sat/hub
(t t
n
)=g1(t t
n
)在区间[

t,0]上的局部近似直线段方程l2:后者的建立方法是:首先,在曲线τ

sat/hub
(t)=g1(t t
n
),t∈[

t,0]上选取两个相邻的点(0,g1(t
n
))和(

t,g1(t
μ
)),其中,t
μ
=t
n

t,t是一个时间增量,且时间增量,且表示返向下行链路自由空间传输时延的最大值。然后,根据上述两点的坐标建立直线段方程,得
[0052][0053]
s3、求解直线l1和直线段l2的交点坐标,得
[0054][0055]
如图4所示为s1~s3的几何示意图,其代表的场景是卫星与主站的距离逐渐减小,因此返向下行链路的传输时延也逐渐减小。
[0056]
s4、估算返向突发bf[n]到达卫星的星历时间,得
[0057]
t

n
=t
n
t
intersect

[0058]
具体的,估算返向突发bf[n]到达卫星的星历时间t

n
的方法如下:如图2所示,在t

n
时刻,返向突发bf[n]从端站发出,并在t

n
时刻到达卫星。此时,卫星的位置是p(t

n
)。随后,bf[n]被卫星转发给主站,并在t
n
时刻到达主站。而在t
n

t

n
的时间内,卫星也从位置p(t

n
)移动到了新的位置p(t
n
)。卫星在轨道p(t

n
)

p(t
n
)上的运动时间与bf[n]在返向下行链路p(t

n
)

h上的传输时延恰好相等,即直线τ
sat/hub
(t)=

t t
n
与曲线τ
sat/hub
(t)=g1(t)在t
n
之前必然相交,且交点就是t

n
(如图3所示,其代表的场景是卫星与主站的距离逐渐减小,因此返向下行链路的传输时延也逐渐减小。在本图中,g1(t

n
)=t
n

t

n
)。因此,求解非线性方程组
[0059][0060]
就可以得到返向突发bf[n]到达卫星的星历时间t

n

[0061]
通常,在低轨卫星通信系统中,返向下行链路的自由空间传输时延很小,卫星在这段时间的移动距离很短,其运动轨迹可以近似为直线段。因此,曲线τ
sat/hub
(t)=g1(t)在t
n
附近的区域也可以近似为直线段。令此直线段的方程为
[0062][0063]
其中,t是一个时间增量。进而,上述问题可以简化为:
[0064]
直线τ
sat/hub
(t)=

t t
n
和直线段的相交问题。
[0065]
综上所述,求解线性方程组
[0066][0067]
就可以得到返向突发bf[n]到达卫星的星历时间t

n

[0068]
s5、建立直角坐标系以星历时间t为横坐标,并以返向突发bf[n]到达卫星的星历时间t

n
为横坐标原点;以卫星/端站间的传输时延τ

sat/rcst
(t|t

n
)为纵坐标,并以0为纵坐标原点。
[0069]
s6、在直角坐标系中,建立两个线性方程。一个是反恒等直线方程l3:τ

sat/rcst
(t|t

n
)=

t,另一个是卫星/端站传输时延与星历时间的函数τ
sat/rcst
(t|t

n
)=g2(t)的坐标平移版本τ

sat/rcst
(t|t

n
)=τ
sat/rcst
(t t

n
|t

n
)=g2(t t

n
)在区间[

t
α
,0]上的局部近似直线段方程l4:后者的建立方法是:首先,在曲线τ

sat/rcst
(t|t

n
)=g2(t t

n
),t∈[

t
α
,0]上选取两个点(

t
α
,g2(t
α
))和(

t
β
,g2(t
β
)),其中,t
α
=t

n

t
α
,t
β
=t

n

t
β
,t
α
和t
β
是两个时间增量,且是两个时间增量,且表示返向上行链路自由空间传输时延的最大值。然后,根据上述两点的坐标建立直线段方程,得
[0070][0071]
其中,g2(t
α
)和g2(t
β
)的计算方法如下:
[0072]
1)根据卫星的星历信息和返向突发bf[n]到达卫星的星历时间t

n
,估算出此时卫星在ecef坐标系下的位置p(t

n
)。
[0073]
2)在时间段[t

n

t
α
,t

n
]内,分别选取两个时刻t
α
=t

n

t
α
和t
β
=t

n

t
β
获取端站在这两个时刻的gnss位置r
gnss
(t
α
)和r
gnss
(t
β
),并换算出它们对应的ecef坐标r(t
α
)和r(t
β
)。
[0074]
3)在ecef坐标系下,分别计算出卫星位置p(t

n
)到端站位置r(t
α
)和r(t
β
)间的传输距离。
[0075]
4)用传输距离除以电磁波的传播速度,得到卫星位置p(t

n
)到端站位置r(t
α
)和r(t
β
)间的传输时延g2(t
α
)和g2(t
β
)。
[0076]
s7、求解直线l3和直线段l4的交点坐标,得
[0077][0078]
如图6所示,为s5~s7的几何示意图,其代表的场景是动中通端站与卫星的距离逐渐减小,因此返向上行链路的传输时延也逐渐减小。
[0079]
s8、估算发送返向突发bf[n]的星历时间,得
[0080]
t

n
=t

n
t

intersect

[0081]
具体的,估算发送返向突发bf[n]的星历时间t

n
的方法如下:如图2所示,在t

n
时刻,端站发出返向突发bf[n],此时,卫星位于p(t

n
)的位置。之后,bf[n]在返向上行链路中传输,而卫星继续沿着箭头指向的方向在轨道上运动。在t

n
时刻,卫星运动到新的位置p(t

n
),而此时,bf[n]也刚好到达卫星,即经过时延t

n

t

n
后bf[n]与卫星相遇。在上述相遇问题中,卫星在轨道p(t

n
)

p(t

n
)上的运动时间与bf[n]在返向上行链路r(t

n
)

p(t

n
)上的传输时延恰好相等,即直线τ
sat/rcst
(t|t

n
)=

t t

n
与曲线τ
sat/rcst
(t|t

n
)=g2(t)在t

n
之前必然相交,且交点就是t

n
(如图5所示,其代表的场景是动中通端站与卫星的距离逐渐减小,因此返向上行链路的传输时延也逐渐减小。在本图中,g2(t

n
)=t

n

t

n
)。因此,求解非线性方程组
[0082][0083]
就可以得到发送返向突发bf[n]的星历时间t

n

[0084]
通常,在t

n
附近一个很小的时间段内,端站的移动距离很短,其运动轨迹可以近似为直线段。因此,曲线τ
sat/rcst
(t|t

n
)=g2(t)在t

n
附近的区域也可以近似为直线段。令此直线段的方程为其中,t
α
是一个时间增量。进而,上述问题可以简化为:
[0085]
直线τ
sat/rcst
(t|t

n
)=

t t

n
和直线段的相交问题;
[0086]
综上所述,求解线性方程组
[0087][0088]
就可以得到发送返向突发bf[n]的星历时间t

n

[0089]
s9、根据端站侧星历时间与ncr时间的映射关系,将发送返向突发bf[n]的星历时间t

n
映射成端站本地的ncr时间。
转载请注明原文地址:https://win.8miu.com/read-250093.html

最新回复(0)