CN114063125A - 一种基于分时双向测距的精密轨道确定和时间同步方法 - Google Patents
一种基于分时双向测距的精密轨道确定和时间同步方法 Download PDFInfo
- Publication number
- CN114063125A CN114063125A CN202111349471.XA CN202111349471A CN114063125A CN 114063125 A CN114063125 A CN 114063125A CN 202111349471 A CN202111349471 A CN 202111349471A CN 114063125 A CN114063125 A CN 114063125A
- Authority
- CN
- China
- Prior art keywords
- time
- spacecraft
- equation
- satellite
- navigation satellite
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000005259 measurement Methods 0.000 claims abstract description 28
- 230000002457 bidirectional effect Effects 0.000 claims abstract description 19
- 238000004891 communication Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 36
- 230000014509 gene expression Effects 0.000 claims description 23
- 238000012937 correction Methods 0.000 claims description 20
- 230000005540 biological transmission Effects 0.000 claims description 5
- 238000009825 accumulation Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 3
- 230000000712 assembly Effects 0.000 description 2
- 238000000429 assembly Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
- G01S19/256—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to timing, e.g. time of week, code phase, timing offset
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
- G01S19/258—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to the satellite constellation, e.g. almanac, ephemeris data, lists of satellites in view
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种基于分时双向测距的精密轨道确定和时间同步方法,包括下列步骤:利用误差公里级的历书,按照收发时隙来执行待确定轨道和时间同步的航天器与相应导航卫星之间的双向测距,以获得星间伪距;所述导航卫星播发星历;所述航天器通过Ka星间链路的通信功能接收所述星历;使用未知系数的多项式表达所述航天器的位置、速度和钟差;归算每一组双向测距的星间伪距并建立双向测量方程;以及使用最小二乘法多次迭代来确定所述航天器的位置多项式、速度多项式和钟差多项式的系数。
Description
技术领域
本发明总的来说涉及航天测量与控制技术领域,具体而言涉及一种基于分时双向测距的精密轨道确定和时间同步方法。
背景技术
目前传统全球定位系统(GNSS)的定位方式一般是由GNSS卫星持续播发无线信号,GNSS接收机同时测得可见卫星的伪距,利用相同时标的伪距、误差和计算得到的发射时刻GNSS卫星瞬时位置构建方程,解出自身的位置和钟差。这种方法依赖于GNSS接收机同时获得多个卫星伪距,从而建立方程。
北斗系统配置Ka波段星间链路载荷。Ka波段的频率更高,波束角小,采用相控阵发射接收。这种体制采用分时双向建链,完成通信和测距。对于某颗星或地面站来说,只需知道一个误差在10公里以内的简约历书或站址,即可完成与对方导航卫星的双向建链。在一个收发时隙内,分接收和发射两个等时段。若某颗星或地面站自身先发射,则对方卫星先接收,后一时段自身接收,对方卫星发射,反之亦然。因此,某颗星或地面自身获得的与一个对方卫星的两个伪距值是在一个收发时隙内两个不同时段完成的。一个收发时隙与一颗卫星建链,之后的时隙才能与其他卫星建链,自身获得的与不同卫星的伪距是在不同收发时隙完成的。
随着低轨卫星或其他配置有Ka星间链路载荷并可接入导航卫星星间网络的远行航天器的发展,如地月空间等GNSS信号很弱但可保证Ka信号强度的情况下,利用Ka波段测得的伪距值完成轨道确定和时间同步成为一项最优选择。
发明内容
从现有技术出发,本发明的任务是提供一种基于分时双向测距的精密轨道确定和时间同步方法,通过该方法,可以在GNSS信号很弱但Ka信号强度稳定的情况下,利用Ka波段分时双向测距的星间伪距,完成航天器的精密卫星轨道确定和时间同步。
根据本发明,前述任务通过一种基于分时双向测距的精密轨道确定和时间同步方法来解决,该方法包括下列步骤:
利用误差公里级的历书,按照收发时隙来执行待确定轨道和时间同步的航天器与相应导航卫星之间的双向测距,以获得星间伪距;
所述导航卫星播发星历;
所述航天器通过Ka星间链路的通信功能接收所述星历;
使用未知系数的多项式表达所述航天器的位置、速度和钟差;
归算每一组双向测距的星间伪距并建立双向测量方程;以及
使用最小二乘法多次迭代来确定所述航天器位置多项式、速度多项式和钟差多项式的系数。
在本发明的一个优选方案中规定,使用未知系数的多项式表达所述航天器的位置、速度和钟差包括:
设参考时刻为tref,将所述航天器的在t时刻的位置表达为表达式F1:
其中bi、ci、di为位置的待确定系数,x、y、z表示所述航天器的位置,参考时刻tref为设定值;
将表达式F1对t求导,得到所述航天器的速度表达式F2:
其中vx、vy、vz为所述航天器的速度;
将所述航天器的钟差表达为表达式F3:
其中,δt(t)为钟差,ai为钟差的待确定系数,tref为参考时刻;所述航天器可根据自身的位置和时钟特性,选择位置和钟差多项式的阶数;
阶数越高,所述多项式的拟合精度越高,所需的数据积累越多,计算复杂度也越高;以及
对于每一组测距对,时隙0时刻t0已知,可用表达式F1、F2和F3分别表示所述航天器在所述时隙0时刻的位置、速度和钟差。
在本发明的另一优选方案中规定,归算每一组双向测距的星间伪距并建立双向测量方程包括:
计算Ka信号的传播时延;
将双向测距的星间伪距归算到所述时隙0时刻;
建立归算后的双向测量方程F6,并改写为方程F7;
将方程F7中的两式相加、相减得到方程F8;
根据所述导航卫星的位置、钟差和所述航天器的位置、钟差改写方程F8得到方程F9;
在轨道待确定系数的近似值处和钟差待确定系数的近似值处,分别将方程F9中两式的右侧线性化展开得到方程F10;
舍去高阶项,将方程F10改写为矩阵形式得到方程F11;以及
对于M组建链测距对,联立M组方程F11得到方程F12。
在本发明的又一优选方案中规定,计算Ka信号的传播时延包括:
根据方程F4计算Ka信号从所述导航卫星到所述航天器的传播时延τ1:
其中,Tr1为航天器测量时刻,Te1为所述导航卫星发送Ka信号的时刻,R(Tr1)代表所述航天器在Tr1时刻的位置,Rj(Te1)代表所述导航卫星在Te1时刻的位置;
根据方程F5计算的Ka信号从所述航天器到所述导航卫星的传播时延τ2:
其中,Tr2为导航卫星测量时刻,Te2为所述航天器发送Ka信号的时刻,Rj(Tr2)代表所述导航卫星在Tr2时刻的位置,R(Te1)代表所述航天器在Te1时刻的位置;
通过迭代方程F4和方程F5计算传播时延τ1和τ2。
在本发明的另一个优选方案中规定,将双向测距的星间伪距归算到所述时隙0时刻包括:
根据下列公式计算归算时间Δt:
其中Δt1为所述航天器的归算时间,Δt2为所述导航卫星的归算时间;
所述航天器和所述导航卫星基于钟面时测量伪距,以δt(t)表示所述航天器在t时刻的钟差,δtj(t)表示所述导航卫星在t时刻的钟差;
其中Δρ1为所述航天器的星间伪距归算修正量,Δρ2为所述导航卫星的星间伪距归算修正量,R(t0)为所述航天器在所述时隙0时刻的位置,Rj(t0)为所述导航卫星在所述时隙0时刻的位置;
其中ρji为所述航天器测量所述导航卫星的星间伪距归算到所述时隙0时刻的值,ρij为所述导航卫星测量所述航天器的星间伪距归算到时所述隙0时刻的值。
在本发明的又一个优选方案中规定,建立归算后的双向测量方程F6,并改写为方程F7包括:
利用归算后的伪距建立双向测量方程F6:
δt(t0)为航天器在时隙0时刻的钟差,δtj(t0)为导航卫星j在时隙0时刻的钟差,c为光速,分别为航天器和导航卫星j的收发通道时延,lrelji、lrelij为相对论效应改正,ltroji和ltroij为大气对流层时延,为相位中心改正,εji和εij为测量噪声;以及
将方程F6改写为方程F7:
在本发明的另一个优选方案中规定,将方程F7中的两式相加、相减得到方程F8包括:
在本发明的又一个优选方案中规定,根据所述导航卫星的位置、钟差和所述航天器的位置、钟差改写方程F8得到方程F9包括:
利用所述星历计算所述导航卫星在所述时隙0时刻的位置、速度和钟差:
根据所述导航卫星的位置、钟差和航天器的位置、钟差改写方程F8得到方程F9:
其中方程F9右侧的系数ai,bi,ci,di未知。
在本发明的另一个优选方案中规定,在轨道待确定系数的近似值处和钟差待确定系数的近似值处,分别将方程F9中两式的右侧线性化展开得到方程F10包括:
其中Aorb为轨道待确定系数,Aclk为钟差待确定系数;
以所述轨道待确定系数和所述钟差待确定系数为变量,将方程F9中两式的右侧分别在所述轨道待确定系数的近似值处和所述钟差待确定系数的近似值处线性化展开得到方程F10:
其中Aorb*为所述轨道待确定系数的近似值,Aclk*为所述钟差待确定系数的近似值,ρ0 ij(Aorb*)为以所述轨道待确定系数的近似值Aorb*计算的t0时刻航天器与导航卫星的距离。
在本发明的另一个优选方案中规定,舍去高阶项,将将方程F10改写为矩阵形式得到方程F11包括:
在本发明的又一个优选方案中规定,对于M组建链测距对,联立M组方程F11得到方程F12包括:
对于M条建链测距对,联立M个方程式F11,并以X1=Aorb-Aorb*和X2=Aclk-Aclk*作为变量,得到方程F12:
其中Borb为M×1矩阵,Horb为M×(3n+3)矩阵,X1为(3n+3)×1矩阵,Bclk为M×1矩阵,Hclk为M×(p+1)矩阵,X2为(p+1)×1矩阵。
在本发明的一个优选方案中规定,使用最小二乘法多次迭代确定所述航天器的位置多项式、速度多项式和钟差多项式的系数包括:
设置一个合理的初值X1 1,X2 1,根据最小二乘法原理迭代公式F12,第k步最优解为:
其中Aorb+1为第K+1步的轨道的多项式系数,Aclkk+1为第K+1步的钟差的多项式系数;以及
设置迭代门限,多次迭代后得到一组位置、速度和钟差的多项式系数解。
本发明至少具有下列有益效果:在GNSS信号很弱但Ka信号强度稳定的情况下,可以通过采用本发明的方法,利用Ka波段分时双向测距的星间伪距,完成航天器的轨道确定和时间同步,而不需要知道航天器的位置、速度和钟差信息,方法简单,精度高。
附图说明
下面结合附图参考具体实施例来进一步阐述本发明。
图1示出了根据本发明的双向伪距归算原理图;
图2示出了根据本发明的一种基于分时双向测距的精密轨道确定和时间同步方法的流程图;
图3示出了根据本发明的多项式阶数与位置的最大误差;
图4示出了根据本发明的5分钟内四阶多项式表达的位置的误差;以及
图5示出了根据本发明的某北斗卫星在5分钟内的定位误差及时间同步误差。
具体实施方式
应当指出,各附图中的各组件可能为了图解说明而被夸大地示出,而不一定是比例正确的。在各附图中,给相同或功能相同的组件配备了相同的附图标记。
在本发明中,各实施例仅仅旨在说明本发明的方案,而不应被理解为限制性的。
在本发明中,除非特别指出,量词“一个”、“一”并未排除多个元素的场景。
在此还应当指出,在本发明的实施例中,为清楚、简单起见,可能示出了仅仅一部分部件或组件,但是本领域的普通技术人员能够理解,在本发明的教导下,可根据具体场景需要添加所需的部件或组件。
在此还应当指出,在本发明的范围内,“相同”、“相等”、“等于”等措辞并不意味着二者数值绝对相等,而是允许一定的合理误差,也就是说,所述措辞也涵盖了“基本上相同”、“基本上相等”、“基本上等于”。
另外,本发明的各方法的步骤的编号并未限定所述方法步骤的执行顺序。除非特别指出,所述方法步骤可以以不同顺序执行。
本发明的方案基于发明人的如下洞察:对于某个航天器,只需知道一个误差公里级历书,即可完成与相应导航卫星的双向建链。航天器获得的与一个导航卫星的两个星间伪距是在一个收发时隙内两个不同时段完成的。航天器获得的与不同导航卫星的星间伪距是在不同收发时隙完成的。若要使用这些分时双向测距的星间伪距进行轨道确定和时间同步,则必须已知航天器的位置、速度和钟差,将不同时段测得的星间伪距归算到同一时刻,然后再根据双向测量方程的计算结果对航天器的位置、速度、钟差进行修正。但在航天器的精确位置、速度和钟差未知的情况下,无法对不同收发时隙获得的星间伪距进行精确归算,也就无法利用这些分时双向测距的星间伪距确定航天器的位置、速度和钟差。为了实现在航天器的位置、速度和钟差信息未知的情况下,利用分时双向测距的星间伪距完成航天器的精密轨道确定和时间同步,发明人提出了一种基于分时双向测距的精密轨道确定和时间同步方法,下面结合具体实施例参考附图进一步阐述本发明。
实施例一
图1示出了根据本发明的双向伪距归算原理图。
航天器获得的与导航卫星j的两个星间伪距是在一个收发时隙内两个不同时段完成的。需要将两个星间伪距归算到同一时刻得到归算后的伪距才能进行计算。如图1所示,Tr1为航天器测量时刻,Tr2为导航卫星j测量时刻,将Tr1时刻获得的航天器测量导航卫星j的星间伪距和Tr2时刻获得的导航卫星j测量航天器的星间伪距都归算到时隙0时刻。
第一步,计算传播时延,τ1为Ka信号导航卫星j到航天器的传播时延,τ2为Ka信号从航天器到导航卫星j的传播时延。第二步,计算归算时间,归算时间是测量时刻与时隙0时刻的差值。第三步,计算伪距归算修正量,利用航天器和导航卫星j的归算时间、航天器和导航卫星在时隙0时刻的钟速、传播时延τ1、传播时延τ2、航天器和导航卫星j在时隙0时刻的位置、航天器在Tr1时刻的位置、导航卫星j在Te1时刻的位置、导航卫星j在Tr2时刻的位置、航天器在时刻Te2的位置计算伪距归算修正量。第四步,利用航天器测量导航卫星j的星间伪距、导航卫星j测量航天器的星间伪距和伪距归算修正量计算出归算后的伪距。
实施例二
图2示出了根据本发明的一种基于分时双向测距的精密轨道确定和时间同步方法的流程图。
一种基于分时双向测距的精密轨道确定和时间同步方法,具体按照下述步骤进行:
步骤1,利用误差公里级的历书,按照收发时隙来执行待确定轨道和时间同步的航天器与相应导航卫星之间的双向测距,以获得星间伪距。
步骤2,导航卫星播发星历,待确定轨道和时间同步的飞行器通过Ka星间链路的通信功能接相应导航卫星星历。
步骤3,使用未知系数的多项式表示待确定轨道和时间同步的航天器的位置、速度和钟差,具体按照下述步骤进行:
步骤3.1,设参考时刻为tref,将待确定轨道和时间同步的航天器在t时刻的位置表达为表达式F1:
其中bi、ci、di为位置的待确定系数,x、y、z表示航天器的位置,参考时刻tref为设定值;以及
航天器可根据自身位置,选择位置多项式的阶数,多项式阶数越高,多项式的拟合精度越高,所需的数据积累越多,计算复杂度也越高。
步骤3.2,将表达式F1对t求导,得到待确定轨道和时间同步的航天器的速度表达式F2:
其中vx、vy、vz为航天器的速度;
步骤3.3,将航天器的钟差表达为表达式F3:
其中,δt(t)为钟差,ai为钟差的待确定系数,tref为参考时刻;
航天器可根据自身时钟特性,选择钟差多项式阶数,阶数越高,钟差多项式的拟合精度越高,所需的数据积累越多,计算复杂度也越高;以及
对于每一组测距对,时隙0时刻t0已知,可用表达式F1、F2和F3分别表示航天器在时隙0时刻的位置、速度和钟差。
步骤4,归算每一组双向测距的星间伪距并建立双向测量方程,具体按照下述步骤进行:
步骤4.1,计算Ka信号的传播时延:
根据方程F4计算Ka信号从导航卫星到航天器j的传播时延τ1:
其中,Tr1为航天器测量时刻,Te1为导航卫星j发送Ka信号的时刻,R(Tr1)代表航天器在Tr1时刻的位置,Rj(Te1)代表导航卫星j在Te1时刻的位置;
根据方程F5计算的Ka信号从航天器到导航卫星j的传播时延τ2:
其中,Tr2为导航卫星j测量时刻,Te2为航天器发送Ka信号的时刻,Rj(Tr2)代表导航卫星j在时刻Tr2的位置,R(Te2)代表航天器在时刻Te2的位置;
通过迭代方程F4和方程F5计算传播时延τ1和τ2;
步骤4.2,将双向测距的星间伪距归算到时隙0时刻,具体按照下述步骤进行:
步骤4.2.1,根据下列公式计算归算时间Δt:
其中Δt1为航天器的归算时间,Δt2为导航卫星j的归算时间;
步骤4.2.2,计算星间伪距归算修正量:
星间伪距是航天器和导航卫星基于钟面时测量的,以δt(t)表示航天器在t时刻的钟差,δtj(t)表示导航卫星在t时刻的钟差;
其中Δρ1为航天器的星间伪距归算修正量,Δρ2为导航卫星j的星间伪距归算修正量,R(t0)为航天器在时隙0时刻的位置,Rj(t0)为导航卫星j在时隙0时刻的位置;
步骤4.2.3,根据星间伪距归算修正量和星间伪距计算归算后的伪距:
其中ρji为航天器测量导航卫星j的星间伪距归算到时隙0时刻的值,ρij为导航卫星j测量航天器的星间伪距归算到时隙0时刻的值;
步骤4.3,建立归算后的双向测量方程F6,并改写为方程F7:
利用归算后的伪距建立双向测量方程F6:
δt(t0)为航天器在时隙0时刻的钟差,δtj(t0)为导航卫星j在时隙0时刻的钟差,c为光速,分别为航天器和导航卫星j的收发通道时延,lrelji、lrelij为相对论效应改正,ltroji和ltroij为大气对流层时延,为相位中心改正,εji和εij为测量噪声;
将方程F6改写为方程F7:
步骤4.4,将方程F7中的两式相加、相减得到方程F8:
步骤4.5,利用星历计算导航卫星j在时隙0时刻的位置、速度和钟差:
步骤4.6,根据导航卫星j的位置、钟差和航天器的位置、钟差改写方程F8得到方程F9:
其中方程F9右侧的系数ai,bi,ci,di未知;
步骤4.7,在轨道待确定系数的近似值处和钟差待确定系数的近似值处,分别将方程F9中两式的右侧线性化展开得到方程F10:
其中Aorb为轨道待确定系数,Aclk为钟差待确定系数;
以轨道待确定系数和钟差待确定系数为变量,将方程F9中两式的右侧分别在轨道待确定系数的近似值处和钟差待确定系数的近似值处线性化展开得到方程F10:
其中Aorb*为轨道待确定系数的近似值,Aclk*为钟差待确定系数的近似值,ρ0 ij(Aorb*)为以轨道待确定系数的近似值Aorb*计算的t0时刻航天器与导航卫星j的距离;
步骤4.8,舍去高阶项,将方程F10改写为矩阵形式得到方程F11:
步骤4.9,对于M组建链测距对,联立M组方程F11得到方程F12包括:
对于M条建链测距对,联立M个方程式F11,并以X1=Aorb-Aorb*和X2=Aclk-Aclk*作为变量,得到方程F12:
其中Borb为M×1矩阵,Horb为M×(3n+3)矩阵,X1为(3n+3)×1矩阵,Bclk为M×1矩阵,Hclk为M×(p+1)矩阵,X2为(p+1)×1矩阵。
步骤5,使用最小二乘法多次迭代来确定航天器的位置多项式、速度多项式和钟差多项式的系数:
设置一个合理的初值X1 1,X2 1,根据最小二乘法原理迭代公式F12,第k步最优解为:
其中Aorbk+1为第K+1步的轨道的多项式系数,Aclkk+1为第K+1步的钟差的多项式系数;以及
设置迭代门限,多次迭代后得到一组轨道和钟差的多项式系数解。
实施例三
图3示出了根据本发明的多项式阶数与位置的最大误差,并且图4示出了根据本发明的5分钟内四阶多项式表达的位置的误差。
随机选择某北斗卫星系统中的一颗导航MEO卫星作为待确定轨道和钟差的航天器,其按照时隙表与其他导航卫星在规定的时间内完成测距和通信。要达到一定的精度,位置、速度和钟差的拟合多项式需满足一定的阶数。对于导航MEO卫星来说,在长时段的情况下,位置的多项式一般要达到8-10阶。数据时段越长,数据积累越多,计算复杂度越高,计算所需的时间也越长。在本发明设定的场景中,一项隐含的需求是尽量快速的实现定位和时间同步,所以数据时段不能过长。因此,多项式阶数可以大大的降低,在计算时,系数矩阵的维度M×(3n+3)也会降低。
在本实施例中,数据时段设置为5分钟。数据时段设置为5分钟时,位置多项式阶数可以是3阶、4阶、5阶,钟差多项式一般设置为2阶。如图3所示,当数据时段为5分钟时,阶数越高,位置多项式的拟合精度越高,多项式表达的位置的误差越小。当多项式阶数为3阶,x的误差为0.0475米,y的误差为0.0013米,z的误差为0.0273米;当多项式阶数为4阶,x的误差为5.9330e-06米,y的误差为1.6153e-04米,z的误差为4.1530e-05米;当多项式阶数为5阶,x的误差为7.1450e-07米,y的误差为5.5134e-07米,z的误差为5.2899e-07米。如图4所示,四阶多项式表达的位置的误差在300秒(5分钟)内的变化,其中位置误差xerr、yerr和zerr在300秒处最小。如图3和图4所示,在5分钟内,如果多项式阶数在4阶及4阶以上,多项式的拟合误差能够降低到毫米以下。
M≥3n+3=15;
当阶数n=5时,建链测距对数应满足:
M≥3n+3=18。
在实际的航天器运行中,时隙表的设计按照天线可见性设置,可能无法满足上述需求。因此选择四阶多项式来表达航天器的位置,此时,建链测距对应满足至少15个的需求。
实施例四
本实施例中针对具体型号卫星,数据时段选择5分钟数据时段,位置多项式的阶数选择4阶,钟差多项式的阶数选择2阶,步骤1-2、步骤4中的步骤4.1-4.5与实施例一中的步骤1-2、步骤4中的步骤4.1-4.5一致。
步骤3,使用未知系数的多项式表示待确定轨道和时间同步的航天器的位置、速度和钟差,具体按照下述步骤进行:
步骤3.1,将待确定轨道和时间同步的航天器的位置表达为表达式F1:
其中bi、ci、di为位置待确定系数,x、y、z表示航天器的位置,n=4,t-tref为5分钟时段,tref为参考时刻;
步骤3.2,将表达式F1对t求导,得到待确定轨道和时间同步的航天器的速度表达式F2:
其中vx、vy、vz为航天器的速度,n=4,t-tref为5分钟时段;
步骤3.3,将航天器的钟差表达为表达式F3:
其中,δt(t)为钟差,ai为钟差的待确定系数,p=2,tref为参考时刻,t-tref为5分钟时段;
对于每一组测距对,时隙0时刻t0已知,可用表达式F1、F2和F3分别表示航天器在所述时隙0时刻的位置、速度和钟差。
步骤4.6,根据导航卫星的位置、钟差和航天器的位置、钟差改写方程F8得到方程F9:
其中方程F9右侧的系数ai,bi,ci,di未知,t0-tref=5分钟,n=4,
p=2。
步骤4.7,在轨道待确定系数的近似值处和钟差待确定系数的近似值处,分别将方程F9中两式的右侧线性化展开得到方程F10:
其中Aorb为轨道待确定系数,Aclk为钟差待确定系数;
以轨道待确定系数和钟差待确定系数为变量,将方程F9中两式的右侧分别在轨道待确定系数的近似值处和钟差待确定系数的近似值处线性化展开得到方程F10:
其中Aorb*为轨道待确定系数的近似值,Aclk*为钟差待确定系数的近似值,ρ0 ij(Aorb*)为以轨道待确定系数的近似值Aorb*计算的t0时刻航天器与导航卫星的距离,t-tref为5分钟时段,n=4,p=2。
步骤4.8,舍去高阶项,将方程F10改写为矩阵形式得到方程F11:
其中t-tref为5分钟时段。
步骤4.9,对于M组建链测距对,联立M组方程F11得到方程F12包括:
对于M条建链测距对,联立M个方程式F11,并以X1=Aorb-Aorb*和X2=Aclk-Aclk*作为变量,得到方程F12:
其中Borb为M×1矩阵,Horb为M×15矩阵,X1为15×1矩阵,Bclk为M×1矩阵,Hclk为M×3矩阵,X2为3×1矩阵。
步骤5,使用最小二乘法多次迭代来确定航天器的位置多项式、速度多项式和钟差多项式的系数:
设置迭代初值:
根据最小二乘法原理迭代公式F12,第k步最优解为:
其中Aorbk+1为第K+1步的轨道的多项式系数,Aclkk+1为第K+1步的钟差的多项式系数;以及
设置迭代门限,多次迭代后得到一组位置、速度和钟差的多项式系数解。
实施例五
图5示出了根据本发明的某北斗卫星在5分钟内的定位误差及时间同步误差。
根据本发明的方法,以某北斗卫星系统内的某导航MEO卫星在某5分钟内得到的双向测距值进行计算,得出一组轨道和钟差的多项式系数解,最终得出某导航MEO卫星的位置和钟差多项式后与标准轨道作差得到误差曲线。如图5所示,其中xerr曲线、yerr曲线和zerr曲线代表某导航MEO卫星的位置误差,clkerr曲线代表某导航MEO卫星的钟差误差。clkerr和xerr的值在0至-0.1米范围内,yerr和zerr的值在0至0.12米范围内,可见无论位置还是钟差,5分钟内的误差都在0.12米范围以内。上述数据表明,采用根据本发明的基于分时双向测距的定位方法,可以较好地完成卫星的轨道确定和时间同步。
本发明至少具有下列有益效果:在GNSS信号很弱但Ka信号强度稳定的情况下,可以通过采用本发明的方法,利用Ka波段分时双向测距的星间伪距,完成航天器的轨道确定和时间同步,而不需要知道航天器的位置、速度和钟差信息,方法简单,精度高。
虽然本发明的一些实施方式已经在本申请文件中予以了描述,但是本领域技术人员能够理解,这些实施方式仅仅是作为示例示出的。本领域技术人员在本发明的教导下可以想到众多的变型方案、替代方案和改进方案而不超出本发明的范围。所附权利要求书旨在限定本发明的范围,并藉此涵盖这些权利要求本身及其等同变换的范围内的方法和结构。
Claims (12)
1.一种基于分时双向测距的精密轨道确定和时间同步方法,该方法包括下列步骤:
利用误差公里级的历书,按照收发时隙来执行待确定轨道和时间同步的航天器与相应导航卫星之间的双向测距,以获得星间伪距;
所述导航卫星播发星历;
所述航天器通过Ka星间链路的通信功能接收所述星历;
使用未知系数的多项式表达所述航天器的位置、速度和钟差;
归算每一组双向测距的星间伪距并建立双向测量方程;以及
使用最小二乘法多次迭代来确定所述航天器的位置多项式、速度多项式和钟差多项式的系数。
2.根据权利要求1所述的方法,其中使用未知系数的多项式表达所述航天器的位置、速度和钟差包括:
设参考时刻为tref,将所述航天器在t时刻的位置表达为表达式F1:
其中bi、ci、di为位置的待确定系数,x、y、z表示所述航天器的位置,参考时刻tref为设定值;
将表达式F1对t求导,得到所述航天器的速度表达式F2:
其中vx、vy、vz为所述航天器的速度;
将所述航天器的钟差表达为表达式F3:
其中,δt(t)为钟差,ai为钟差的待确定系数,tref为参考时刻;
所述航天器可根据自身的位置和时钟特性,选择位置和钟差多项式的阶数;
阶数越高,所述多项式的拟合精度越高,所需的数据积累越多,计算复杂度也越高;以及
对于每一组测距对,时隙0时刻t0已知,可用表达式F1、F2和F3分别表示所述航天器在所述时隙0时刻的位置、速度和钟差。
3.根据权利要求1所述的方法,其中归算每一组双向测距的星间伪距并建立双向测量方程包括:
计算Ka信号的传播时延;
将双向测距的星间伪距归算到所述时隙0时刻;
建立归算后的双向测量方程F6,并改写为方程F7;
将方程F7中的两式相加、相减得到方程F8;
根据所述导航卫星的位置、钟差和所述航天器的位置、钟差改写方程F8得到方程F9;
在轨道待确定系数的近似值处和钟差待确定系数的近似值处,分别将方程F9中两式的右侧线性化展开得到方程F10;
舍去高阶项,将方程F10改写为矩阵形式得到方程F11;以及
对于M组建链测距对,联立M组方程F11得到方程F12。
5.根据权利要求3所述的方法,其中将双向测距的星间伪距归算到所述时隙0时刻包括:
根据下列公式计算归算时间Δt:
其中Δt1为所述航天器的归算时间,Δt2为所述导航卫星的归算时间;
所述航天器和所述导航卫星基于钟面时测量伪距,以δt(t)表示所述航天器在t时刻的钟差,δtj(t)表示所述导航卫星在t时刻的钟差;
其中Δρ1为所述航天器的星间伪距归算修正量,Δρ2为所述导航卫星的星间伪距归算修正量,R(t0)为所述航天器在所述时隙0时刻的位置,Rj(t0)为所述导航卫星在所述时隙0时刻的位置;
其中ρji为所述航天器测量所述导航卫星的星间伪距归算到所述时隙0时刻的值,ρij为所述导航卫星测量所述航天器的星间伪距归算到时所述隙0时刻的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111349471.XA CN114063125A (zh) | 2021-11-15 | 2021-11-15 | 一种基于分时双向测距的精密轨道确定和时间同步方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111349471.XA CN114063125A (zh) | 2021-11-15 | 2021-11-15 | 一种基于分时双向测距的精密轨道确定和时间同步方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114063125A true CN114063125A (zh) | 2022-02-18 |
Family
ID=80272381
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111349471.XA Pending CN114063125A (zh) | 2021-11-15 | 2021-11-15 | 一种基于分时双向测距的精密轨道确定和时间同步方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114063125A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114894199A (zh) * | 2022-06-16 | 2022-08-12 | 中国科学院空间应用工程与技术中心 | 一种地月空间航天器的天基测定轨方法 |
CN115598676A (zh) * | 2022-10-17 | 2023-01-13 | 北京航天飞行控制中心(Cn) | 星载多模gnss融合精密定轨方法和装置 |
-
2021
- 2021-11-15 CN CN202111349471.XA patent/CN114063125A/zh active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114894199A (zh) * | 2022-06-16 | 2022-08-12 | 中国科学院空间应用工程与技术中心 | 一种地月空间航天器的天基测定轨方法 |
CN114894199B (zh) * | 2022-06-16 | 2023-04-11 | 中国科学院空间应用工程与技术中心 | 一种地月空间航天器的天基测定轨方法 |
CN115598676A (zh) * | 2022-10-17 | 2023-01-13 | 北京航天飞行控制中心(Cn) | 星载多模gnss融合精密定轨方法和装置 |
CN115598676B (zh) * | 2022-10-17 | 2023-05-05 | 北京航天飞行控制中心 | 星载多模gnss融合精密定轨方法和装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110031881B (zh) | 高精度星间激光测距辅助精密单点定位的方法 | |
US9365303B2 (en) | Position and elevation acquisition for orbit determination | |
CN114063125A (zh) | 一种基于分时双向测距的精密轨道确定和时间同步方法 | |
CN113671543B (zh) | 一种基于三频模式的高精度星地时间比对方法及系统 | |
CN107664947B (zh) | 一种卫星双向时间比对系统校准的装置及方法 | |
CN101082663A (zh) | 转发式卫星导航中的虚拟原子钟方法 | |
EP2699934B1 (en) | Process and system to determine temporal changes in retransmission and propagation of signals used to measure distances, synchronize actuators and georeference applications | |
CN111308524B (zh) | 一种面向微小卫星多星编队的测量方法 | |
US7920650B2 (en) | Method and apparatus for evaluating a clock in a satellite | |
CN109752737B (zh) | 一种导航卫星星间Ka波段双向测量伪距的预处理方法 | |
WO2009130305A1 (en) | Method of positioning using satellites | |
CN109738914B (zh) | 应用于隧道内导航仿真系统的本地时钟频率偏差修正方法 | |
CN114935767A (zh) | 基于干涉时差测量的卫星无源定位授时方法及系统 | |
CN115877415B (zh) | 一种超高精度空间航天器钟差估计方法及系统 | |
KR101608809B1 (ko) | Gbas 운용범위 확장을 위한 벡터형 오차 보정 장치 및 방법 | |
CN113382466B (zh) | 一种基于时分双工体制下的新型时延监测及校准方法 | |
CN112782734B (zh) | 一种导航星座自主导航电文系统差检测与处理方法 | |
CN111077759B (zh) | 星载原子钟性能评估方法 | |
CN114280650A (zh) | 基于leo获取gnss卫星频间偏差的方法和系统 | |
WO2012093249A1 (en) | A system and method for determining the position of an object | |
RU2402786C1 (ru) | Способ определения местоположения объектов-потребителей навигационной информации и устройство для его реализации | |
Kaufmann et al. | Signal transceiver transit times and propagation delay corrections for ranging and georeferencing applications | |
CN115390104B (zh) | 一种导航卫星时延偏差建模方法 | |
CN117471506B (zh) | 一种构建导航虚拟时空基准的方法 | |
Ko | Analysis of the characteristics of GRACE dual one-way ranging system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |