CN109061701A - 一种适用于星载多天线融合载波相位差分测量方法 - Google Patents
一种适用于星载多天线融合载波相位差分测量方法 Download PDFInfo
- Publication number
- CN109061701A CN109061701A CN201810815780.3A CN201810815780A CN109061701A CN 109061701 A CN109061701 A CN 109061701A CN 201810815780 A CN201810815780 A CN 201810815780A CN 109061701 A CN109061701 A CN 109061701A
- Authority
- CN
- China
- Prior art keywords
- antenna
- aircraft
- satellite
- main
- flying
- 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.)
- Granted
Links
- 238000000691 measurement method Methods 0.000 title abstract description 5
- 239000013598 vector Substances 0.000 claims abstract description 83
- 230000004927 fusion Effects 0.000 claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 16
- 238000000034 method Methods 0.000 claims description 35
- 239000011159 matrix material Substances 0.000 claims description 29
- 238000006243 chemical reaction Methods 0.000 claims description 17
- 230000009466 transformation Effects 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000007667 floating Methods 0.000 claims description 7
- 230000000007 visual effect Effects 0.000 abstract 1
- 230000003287 optical effect Effects 0.000 description 3
- 238000004804 winding Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000003032 molecular docking Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
- G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
-
- 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/40—Correcting position, velocity or attitude
- G01S19/41—Differential correction, e.g. DGPS [differential GPS]
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)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Details Of Aerials (AREA)
Abstract
本发明提供了一种适用于星载多天线融合载波相位差分测量方法:解算绕飞飞行器主天线与主飞行器天线之间的基线矢量根据绕飞飞行器主天线与主飞行器天线之间的基线矢量ECEF坐标系下绕飞飞行器主天线与从天线之间基线矢量解算绕飞飞行器从天线与主飞行器天线之间的基线矢量估计值当绕飞飞行器主天线视场变差,绕飞飞行器主天线和从天线测量得到的观测量均无法独立完成RTK运算时,根据绕飞飞行器从天线与主飞行器天线之间的基线矢量估计值辅助解算绕飞飞行器从天线与主飞行器天线的载波相位模糊度整数解,求解绕飞飞行器从天线与主飞行器天线的基线矢量bB解决主飞行器与绕飞飞行器的精密相对测量问题。
Description
技术领域
本发明涉及一种适用星载多天线融合载波相位差分测量方法,属于空间飞行器技术领域。
背景技术
RTK(Real Time Kinematic)是载波相位差分测量技术,可快速测量高精度相对位置,因而可应用于星间测量系统或者飞行器交会对接控制系统。 RTK技术广泛应用于航天航空等高科技领域中。
随着轨道飞行器的技术发展,编队飞行出现了伴飞,绕飞等复杂飞行模式。在伴飞,绕飞等飞行过程中,控制失当将引起飞行器相互碰撞,因此精确测定飞行器间相对位置是保证飞行器安全的重要技术保障。与交会对接不同的是,在绕飞过程中,绕飞飞行器往往承担着光学拍照摄像的任务,在光学拍照摄像任务中,绕飞飞行器光学镜头需一直指向主飞行器,这意味着绕飞飞行器会随着绕飞过程而旋转,因而姿态不再呈三轴稳定对地状态而呈现大幅度姿态机动。姿态机动过程会造成安装在绕飞飞行器上的卫星导航天线视场不断变化,从而造成在某些角度导航天线只能接收到少量导航卫星信号或者根本接收不到卫星信号。在卫星星体上布置多个天线可以解决卫星星体转动时连续观测的问题,但是又带来了绕飞飞行器切换选取参加RTK运算天线的问题,在切换过程中,新参与RTK运算的天线需要重新搜索载波相位模糊度,这会造成切换过程中RTK中断,从而造成主飞行器与绕飞飞行器的观测过程中断,给飞行器的安全运行带来一定的威胁。
发明内容
本发明所解决的技术问题是:克服现有技术不足,提供一种适用于星载多天线融合载波相位差分测量方法,根据卫星的在轨姿态任务,在卫星星体上安装两个或者三个天线,保证在任何时间任何姿态下,至少有一个天线可与主飞行器完成RTK运算。通过硬件设计和算法设计,保证姿态转动过程中主飞行器与绕飞飞行器的精密相对测量。
本发明的技术解决方案是:一种适用于星载多天线融合载波相位差分测量方法,适用于星间测量系统,所述星间测量系统包括主飞行器和绕飞飞行器,绕飞飞行器上布置多个天线,选取视场内收星状态良好,天线为绕飞飞行器的主天线,其余天线为从天线,当绕飞飞行器主天线能够独立完成RTK 运算时,执行如下步骤:
(1.1)、采用绕飞飞行器主天线与主飞行器天线测量得到的观测量完成RTK运算,得到绕飞飞行器主天线与主飞行器天线间载波相位模糊度整数解,解算绕飞飞行器主天线与主飞行器天线之间的基线矢量
(1.2)、根据绕飞飞行器主天线与主飞行器天线之间的基线矢量 ECEF坐标系下绕飞飞行器主天线与从天线之间基线矢量解算绕飞飞行器从天线与主飞行器天线之间的基线矢量估计值;
当绕飞飞行器从天线视场较差时,绕飞飞行器从天线测量得到的观测量均无法独立完成RTK运算,执行如下步骤,直至绕飞飞行器任一天线收星状态良好,能够独立完成RTK运算时,将该天线切换为主天线:
(1.3)、根据绕飞飞行器从天线与主飞行器天线之间的基线矢量估计值辅助解算绕飞飞行器从天线与主飞行器天线的载波相位模糊度整数解,进而求解绕飞飞行器从天线与主飞行器天线的基线矢量bB。
所述绕飞飞行器从天线与主飞行器天线间的基线矢量估计值的计算方法为:
所述步骤(1.3)的具体计算过程为:
(1.3.1)、选取绕飞飞行器从天线与主飞行器天线的N颗共视卫星作为观测卫星,将其中一颗卫星记为观测主星,其余各颗卫星记为观测从星,N≥4;
(1.3.2)、采用主飞行器天线位置观测主星的位置计算主飞行器到观测主星间单位方向矢量采用主飞行器天线位置各观测从星的位置计算主飞行器天线到各观测从星间的单位方向矢量k=2~N;
(1.3.3)、根据主飞行器天线和绕飞飞行器从天线测量得到的观测卫星的载波相位,计算观测主星与观测从星之间的载波相位双差k=2~N;
(1.3.4)、采用绕飞飞行器从天线与主飞行器天线基线矢量估计值主飞行器天线到观测主星间单位方向矢量主飞行器天线到观测从星间单位方向矢量以及观测主星与各观测从星之间的载波相位双差得到观测主星与各观测从星之间的双差载波相位模糊度浮点解k=2~N;
(1.3.5)、采用步骤(1.3.4)中的载波相位模糊度浮点解计算载波相位模糊度整数解k=2~N;
(1.3.6)、根据主飞行器天线到观测主星间单位方向矢量主飞行器天线到观测从星间单位方向矢量观测主星与各观测从星之间的载波相位双差和载波相位模糊度整数解求解绕飞飞行器从天线与主飞行器天线的基线矢量bB。
第k颗主飞行器天线到观测从星间单位方向矢量k=2~N的计算公式为:
所述步骤(1.3.4)中观测主星与各观测从星之间的双差载波相位模糊度浮点解的具体计算公式为:
所述绕飞飞行器从天线与主飞行器天线的基线矢量bB的具体计算公式为:
所述ECEF坐标系下绕飞飞行器主天线与从天线之间基线矢量根据星体坐标系到ECEF坐标系的转换矩阵以及星体坐标系下其他天线到主天线基线矢量计算得到。
所述星体坐标系到ECEF坐标系的转换矩阵根据下列方法获得:
(a)、根据绕飞飞行器所在的纬度、经度、高度(Lb,λb,hb),计算卫星导航坐标系到ECEF坐标系的转换矩阵
(b)、获取绕飞飞行器卫星姿态角计算卫星本体坐标系到卫星导航坐标系的转化矩阵
(c)、根据卫星导航坐标系到ECEF坐标系的转换矩阵和卫星本体坐标系到卫星导航坐标系的转化矩阵得到星体坐标系到ECEF坐标系的转换矩阵
所述卫星导航坐标系到ECEF坐标系的转换矩阵为:
卫星本体坐标系到卫星导航坐标系的转化矩阵为:
本发明与现有技术相比的优点在于:
(1)、本发明利用绕飞飞行器主天线与主飞行器天线的载波相位模糊度固定信息来计算绕飞飞行器从天线与主飞行器天线的载波相位模糊度,避免了切换过程中,从天线与主飞行器天线重新搜搜载波相位模糊度。
(2)、本发明绕飞飞行器主天线与主飞行器天线RTK运算中,载波相位模糊度已固定时,绕飞飞行器从天线与主飞行器天线可在1个历元时间内固定载波相位模糊度。
(3)、本发明绕飞飞行器主天线与主飞行器天线RTK运算中,载波相位模糊度已固定时,绕飞飞行器从天线与主飞行器天线只需有4颗共视卫星,即可单独解算基线。
(4)、本发明绕飞飞行器接收卫星数大于等于4的天线均可与主飞行器天线保持RTK载波相位模糊度固定的状态,在天线切换过程中,可保证无缝基线测量。
附图说明
图1本发明实施例绕飞飞行器星体表面安装双天线示意图;
图2本发明实施例多天线融合RTK示意图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细描述。
本发明提供了一种适用于星载多天线融合RTK方法,适用于星间测量系统,所述星间测量系统包括主飞行器和绕飞飞行器,绕飞飞行器上布置多个天线,单独接收导航卫星信号,计算每个卫星的伪距和载波相位,选取视场内收星状态良好,天线为绕飞飞行器的主天线,其余天线为从天线。
本发明通过信息融合方法,减少绕飞飞行器切换参与RTK运算天线时新天线完成载波相位模糊度固定的时间,实现了绕飞飞行器姿态机动过程中切换天线时RTK运算不中断,保证在任何姿态下至少有一个天线可始终与主飞行器的接收机数据完成RTK运算,从而完成在各种复杂姿态下主飞行器与绕飞飞行器的精密相对测量问题。该方法具体步骤为:
当绕飞飞行器主天线能够独立完成RTK运算时,执行如下步骤:
(1.1)、采用绕飞飞行器主天线与主飞行器天线测量得到的观测量完成RTK运算,得到绕飞飞行器主天线与主飞行器天线间载波相位模糊度整数解,解算绕飞飞行器主天线与主飞行器天线之间的基线矢量
(1.2)、根据绕飞飞行器主天线与主飞行器天线之间的基线矢量 ECEF坐标系下绕飞飞行器主天线与从天线之间基线矢量解算绕飞飞行器从天线与主飞行器天线之间的基线矢量估计值
所述ECEF坐标系下绕飞飞行器主天线与从天线之间基线矢量根据星体坐标系到ECEF坐标系的转换矩阵以及星体坐标系下其他天线到主天线基线矢量计算得到:
所述星体坐标系到ECEF坐标系的转换矩阵根据下列方法获得:
(a)、根据绕飞飞行器所在的纬度、经度、高度(Lb,λb,hb),计算卫星导航坐标系到ECEF坐标系的转换矩阵
(b)、获取绕飞飞行器卫星姿态角计算卫星本体坐标系到卫星导航坐标系的转化矩阵
(c)、根据卫星导航坐标系到ECEF坐标系的转换矩阵和卫星本体坐标系到卫星导航坐标系的转化矩阵得到星体坐标系到ECEF坐标系的转换矩阵
当绕飞飞行器从天线视场较差,绕飞飞行器从天线测量得到的观测量均无法独立完成RTK运算时,执行如下步骤,直至绕飞飞行器任一天线收星状态良好,能够独立完成RTK运算时,将该天线切换为主天线:
(1.3)、根据绕飞飞行器从天线与主飞行器天线之间的基线矢量估计值辅助解算绕飞飞行器从天线与主飞行器天线的载波相位模糊度整数解,进而求解绕飞飞行器从天线与主飞行器天线的基线矢量bB,具体计算过程为:
(1.3.1)、选取绕飞飞行器从天线与主飞行器天线的N颗共视卫星作为观测卫星,将其中一颗卫星记为观测主星,其余各颗卫星记为观测从星,N≥4;
(1.3.2)、采用主飞行器天线位置观测主星的位置计算主飞行器到观测主星间单位方向矢量采用主飞行器天线位置各观测从星的位置计算主飞行器天线到各观测从星间的单位方向矢量k=2~N;具体计算方法为:
(1.3.3)、根据主飞行器天线和绕飞飞行器从天线测量得到的观测卫星的载波相位,计算观测主星与观测从星之间的载波相位双差k=2~N;所述载波相位双差为主飞行器天线测量得到的观测主星载波相位与绕飞飞行器从天线测量得到的观测主星载波相位之差减去主飞行器天线测量得到的观测从星载波相位与绕飞飞行器从天线测量得到的观测从星载波相位之差;
(1.3.4)、采用绕飞飞行器从天线与主飞行器天线基线矢量估计值主飞行器天线到观测主星间单位方向矢量主飞行器天线到观测从星间单位方向矢量以及观测主星与各观测从星之间的载波相位双差得到观测主星与各观测从星之间的双差载波相位模糊度浮点解k=2~N:
(1.3.5)、采用步骤(1.3.4)中的载波相位模糊度浮点解计算载波相位模糊度整数解k=2~N;
具体计算公式为:
其中floor函数为向下取整函数,为载波相位模糊度整数解;
(1.3.6)、根据主飞行器天线到观测主星间单位方向矢量主飞行器天线到观测从星间单位方向矢量观测主星与各观测从星之间的载波相位双差和载波相位模糊度整数解求解绕飞飞行器从天线与主飞行器天线的基线矢量bB:
实施例:
以下以主飞行器和双天线绕飞航天器为例进行阐述,所述方法同样适用于三天线或者多天线。如图1所示。
如图2所示,在绕飞飞行器卫星星体的上表面布置两个天线,分别为天线A与天线B,卫星的星体坐标系为OXYZ,在卫星的星体坐标系中,天线 A与天线B间基线矢量为在绕飞飞行器天线A视场处于最优状态时,选择绕飞飞行器天线A与主飞行器天线完成RTK运算,同时解算绕飞飞行器天线B与主飞行器天线的基线,当绕飞飞行器天线A视场比绕飞飞行器天线B差时,选取天线B为主天线,完成RTK运算,同时,辅助解算绕飞飞行器天线A与主飞行器天线的基线。保证在任一时刻,有不少于1个天线可与主飞行器完成RTK运算。
在主飞行器上表面安装导航天线,装配卫星导航接收机,利用传统方法求取主飞行器在ECEF坐标系中的位置坐标得到卫星星历文件 Nb,获取伪距以及载波相位其中为主飞行器收到的第i颗卫星的伪距,为主飞行器收到的第i颗卫星的载波相位。
天线A和天线B分别进行收星处理,分别获取伪距与载波相位等测量量,天线A的伪距为载波相位为其中为绕飞飞行器天线A收到的第j颗卫星的伪距,为绕飞飞行器天线A收到的第j颗卫星的载波相位。天线B的伪距为载波相位为其中为绕飞飞行器天线B收到的第k颗卫星的伪距,为绕飞飞行器天线B收到的第k颗卫星的载波相位。
比较天线A和天线B的接收卫星数目,若j≥k,则选取天线A与主飞行天线开展载波相位模糊度搜索,若j,则选取天线B与主飞行天线开展载波相位模糊度搜索。为便于论述,假设选取天线A与主飞行器天线进行 RTK运算。
(1.1)、采用传统载波相位差分方法,搜索天线A与主飞行器天线间载波相位模糊度,固定后获取天线A与主飞行器天线间基线矢量bA。
(1.2)、解算绕飞飞行器从天线与主飞行器天线之间的基线矢量估计值
(1.2.1)、将主飞行器在ECEF坐标系中的位置坐标转化为纬度、经度、高度(Lb,λb,hb)。;
(1.2.2)、计算由主飞行器所在位置导航坐标系到ECEF坐标系的转换矩阵,具体由公式如下:
给出。
(1.2.3)、收绕飞飞行器星务计算机给出的卫星姿态角计算由卫星本体坐标系到卫星导航坐标系的转化矩阵,具有由公式如下:
给出。
(1.2.4)采用步骤(1.2.2)以及步骤(1.2.3)得到的计算由星体坐标系到ECEF坐标系的转换矩阵,具体公式如下:
(1.2.5)采用步骤(1.2.4)得到的以及基线矢量计算在ECEF 坐标系下,天线A到天线B之间的基线矢量具体公式如下:
其中,天线A到天线B之间的基线矢量为天线A相位中心与天线B 相位中心在ECEF坐标系下的向量。
(1.2.6)如图2所示,依据空间矢量关系可以得到绕飞飞行器天线B 与主飞行器天线间的基线矢量,具体公式如下:
(1.3)、根据绕飞飞行器天线B与主飞行器天线之间的基线矢量估计值辅助解算绕飞飞行器天线B与主飞行器天线的载波相位模糊度整数解,进而求解绕飞飞行器天线B与主飞行器天线的基线矢量bB。
(1.3.1)、选取绕飞飞行器天线B与主飞行器天线的共视卫星,为便于阐述,以共视卫星为1号卫星、2号卫星、3号卫星、4号卫星为例,其中1号卫星为仰角最高卫星。
(1.3.2)选取最高仰角卫星为主星,选取1号卫星为主星。利用步骤 (1)中得到的星历文件Nb,获取1号卫星的位置
(1.3.3)采用主飞行器位置以及步骤(13)中的1号导航卫星的位置计算主飞行器到1号导航卫星间的单位方向矢量,具体公式如下:
其中为主飞行器到1号导航卫星的单位方向矢量。
(1.3.4)利用星历文件Nb,获取k号卫星的位置采用主飞行器位置以及第k号导航卫星的位置计算主飞行器到k号导航卫星间的单位方向矢量,具体公式如下:
其中为主飞行器到k号导航卫星的单位方向矢量。
主飞行器到1号卫星、2号卫星、3号卫星、4号卫星的单位方向矢量分别记为:
(1.3.5)采用主飞行器天线接收到的载波相位与绕飞飞行器天线B接收到的载波相位,计算1号星与第k号星的载波相位双差,具体公式如下:
其中,为主飞行器天线接收到的1号卫星载波相位,为主飞行器天线接收到的k号卫星载波相位,为绕飞飞行器天线B接收到的1号卫星载波相位,为绕飞飞行器天线B接收到的k号卫星载波相位。
2号卫星、3号卫星、4号卫星与1号卫星的载波相位双差值分别为:
(1.3.6)采用天线B与主飞行器天线基线矢量以及步骤(1.3.3)得到的求取载波相位模糊度浮点解,具体公式如下:
(1.3.7)采用载波相位模糊度浮点解计算载波相位模糊度整数解,具体由公式如下:
给出,其中floor函数为向下取整函数,为载波相位模糊度整数解。
2号卫星、3号卫星、4号卫星与1号卫星的载波相位模糊度整数解分别为:
(1.3.8)利用单位方向矢量、载波相位双差值、载波相位模糊度整数解求解绕飞飞行器天线B与主飞行器天线的基线矢量,具体由公式如下:
给出。其中,bB为绕飞飞行器天线B与主飞行器RTK运算得到的基线矢量。
以上所述,仅为本发明一个具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (10)
1.一种适用于星载多天线融合载波相位差分测量方法,适用于星间测量系统,所述星间测量系统包括主飞行器和绕飞飞行器,绕飞飞行器上布置多个天线,选取视场内收星状态良好,天线为绕飞飞行器的主天线,其余天线为从天线,其特征在于:
当绕飞飞行器主天线能够独立完成RTK运算时,执行如下步骤:
(1.1)、采用绕飞飞行器主天线与主飞行器天线测量得到的观测量完成RTK运算,得到绕飞飞行器主天线与主飞行器天线间载波相位模糊度整数解,解算绕飞飞行器主天线与主飞行器天线之间的基线矢量
(1.2)、根据绕飞飞行器主天线与主飞行器天线之间的基线矢量ECEF坐标系下绕飞飞行器主天线与从天线之间基线矢量解算绕飞飞行器从天线与主飞行器天线之间的基线矢量估计值
当绕飞飞行器从天线视场较差时,绕飞飞行器从天线测量得到的观测量均无法独立完成RTK运算,执行如下步骤,直至绕飞飞行器任一天线收星状态良好,能够独立完成RTK运算时,将该天线切换为主天线:
(1.3)、根据绕飞飞行器从天线与主飞行器天线之间的基线矢量估计值辅助解算绕飞飞行器从天线与主飞行器天线的载波相位模糊度整数解,进而求解绕飞飞行器从天线与主飞行器天线的基线矢量bB。
2.根据权利要求1所述的一种适用于星载多天线融合载波相位差分测量方法,其特征在于:所述绕飞飞行器从天线与主飞行器天线间的基线矢量估计值的计算方法为:
3.根据权利要求1所述的一种适用于星载多天线融合载波相位差分测量方法,其特征在于:所述步骤(1.3)的具体计算过程为:
(1.3.1)、选取绕飞飞行器从天线与主飞行器天线的N颗共视卫星作为观测卫星,将其中一颗卫星记为观测主星,其余各颗卫星记为观测从星,N≥4;
(1.3.2)、采用主飞行器天线位置观测主星的位置计算主飞行器到观测主星间单位方向矢量采用主飞行器天线位置各观测从星的位置计算主飞行器天线到各观测从星间的单位方向矢量k=2~N;
(1.3.3)、根据主飞行器天线和绕飞飞行器从天线测量得到的观测卫星的载波相位,计算观测主星与观测从星之间的载波相位双差k=2~N;
(1.3.4)、采用绕飞飞行器从天线与主飞行器天线基线矢量估计值主飞行器天线到观测主星间单位方向矢量主飞行器天线到观测从星间单位方向矢量以及观测主星与各观测从星之间的载波相位双差得到观测主星与各观测从星之间的双差载波相位模糊度浮点解k=2~N;
(1.3.5)、采用步骤(1.3.4)中的载波相位模糊度浮点解计算载波相位模糊度整数解k=2~N;
(1.3.6)、根据主飞行器天线到观测主星间单位方向矢量主飞行器天线到观测从星间单位方向矢量观测主星与各观测从星之间的载波相位双差和载波相位模糊度整数解求解绕飞飞行器从天线与主飞行器天线的基线矢量bB。
4.根据权利要求3所述的一种适用于星载多天线融合载波相位差分测量方法,其特征在于:第k颗主飞行器天线到观测从星间单位方向矢量k=2~N的计算公式为:
5.根据权利要求3所述的一种适用于星载多天线融合载波相位差分测量方法,其特征在于:所述步骤(1.3.4)中观测主星与各观测从星之间的双差载波相位模糊度浮点解的具体计算公式为:
6.根据权利要求3所述的一种适用于星载多天线融合载波相位差分测量方法,其特征在于:所述绕飞飞行器从天线与主飞行器天线的基线矢量bB的具体计算公式为:
7.根据权利要求1所述的一种适用于星载多天线融合载波相位差分测量方法,其特征在于:所述ECEF坐标系下绕飞飞行器主天线与从天线之间基线矢量根据星体坐标系到ECEF坐标系的转换矩阵以及星体坐标系下其他天线到主天线基线矢量计算得到。
8.根据权利要求7所述的一种适用于星载多天线融合载波相位差分测量方法,其特征在于:所述星体坐标系到ECEF坐标系的转换矩阵限据下列方法获得:
(a)、根据绕飞飞行器所在的纬度、经度、高度(Lb,λb,hb),计算卫星导航坐标系到ECEF坐标系的转换矩阵
(b)、获取绕飞飞行器卫星姿态角(θ,ψ),计算卫星本体坐标系到卫星导航坐标系的转化矩阵
(c)、根据卫星导航坐标系到ECEF坐标系的转换矩阵和卫星本体坐标系到卫星导航坐标系的转化矩阵得到星体坐标系到ECEF坐标系的转换矩阵
9.根据权利要求8所述的一种适用于星载多天线融合载波相位差分测量方法,其特征在于:所述卫星导航坐标系到ECEF坐标系的转换矩阵为:
10.根据权利要求8所述的一种适用于星载多天线融合载波相位差分测量方法,其特征在于卫星本体坐标系到卫星导航坐标系的转化矩阵为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810815780.3A CN109061701B (zh) | 2018-07-24 | 2018-07-24 | 一种适用于星载多天线融合载波相位差分测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810815780.3A CN109061701B (zh) | 2018-07-24 | 2018-07-24 | 一种适用于星载多天线融合载波相位差分测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109061701A true CN109061701A (zh) | 2018-12-21 |
CN109061701B CN109061701B (zh) | 2021-02-09 |
Family
ID=64836167
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810815780.3A Active CN109061701B (zh) | 2018-07-24 | 2018-07-24 | 一种适用于星载多天线融合载波相位差分测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109061701B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111624584A (zh) * | 2020-03-20 | 2020-09-04 | 中国人民解放军火箭军工程大学 | 一种非协作目标激光诱偏距离测量系统及方法 |
CN112068166A (zh) * | 2020-09-14 | 2020-12-11 | 天津云遥宇航科技有限公司 | 一种用于双天线卫星的卫星定位系统 |
CN113064194A (zh) * | 2020-08-12 | 2021-07-02 | 中国科学院微小卫星创新研究院 | 基于互联网星座的皮纳卫星通信、导航、定姿方法 |
CN113094371A (zh) * | 2021-04-14 | 2021-07-09 | 嘉兴毕格智能科技有限公司 | 一种用户自定义坐标系的实现方法 |
CN113126022A (zh) * | 2021-04-14 | 2021-07-16 | 成都金诺信高科技有限公司 | 一种双天线定位测向方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6313788B1 (en) * | 1998-08-14 | 2001-11-06 | Seagull Technology, Inc. | Method and apparatus for reliable inter-antenna baseline determination |
US20020029110A1 (en) * | 2000-06-01 | 2002-03-07 | Masaru Fukuda | System for determining the heading and/or attitude of a body |
CN101614802A (zh) * | 2009-07-28 | 2009-12-30 | 中国电子科技集团公司第二十八研究所 | 一种导航卫星姿态测量方法 |
CN105807293A (zh) * | 2016-05-27 | 2016-07-27 | 重庆卓观科技有限公司 | 基于soc的单板多天线定姿接收机 |
CN107102346A (zh) * | 2017-06-08 | 2017-08-29 | 中国电子科技集团公司第五十四研究所 | 一种基于北斗系统的多天线测姿方法 |
CN108205150A (zh) * | 2016-12-19 | 2018-06-26 | 千寻位置网络有限公司 | 差分定位方法及系统 |
-
2018
- 2018-07-24 CN CN201810815780.3A patent/CN109061701B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6313788B1 (en) * | 1998-08-14 | 2001-11-06 | Seagull Technology, Inc. | Method and apparatus for reliable inter-antenna baseline determination |
US20020029110A1 (en) * | 2000-06-01 | 2002-03-07 | Masaru Fukuda | System for determining the heading and/or attitude of a body |
CN101614802A (zh) * | 2009-07-28 | 2009-12-30 | 中国电子科技集团公司第二十八研究所 | 一种导航卫星姿态测量方法 |
CN105807293A (zh) * | 2016-05-27 | 2016-07-27 | 重庆卓观科技有限公司 | 基于soc的单板多天线定姿接收机 |
CN108205150A (zh) * | 2016-12-19 | 2018-06-26 | 千寻位置网络有限公司 | 差分定位方法及系统 |
CN107102346A (zh) * | 2017-06-08 | 2017-08-29 | 中国电子科技集团公司第五十四研究所 | 一种基于北斗系统的多天线测姿方法 |
Non-Patent Citations (2)
Title |
---|
唐兴: "GPS在姿态测量中的应用", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
高源骏: "GPS测姿算法与天线布局研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111624584A (zh) * | 2020-03-20 | 2020-09-04 | 中国人民解放军火箭军工程大学 | 一种非协作目标激光诱偏距离测量系统及方法 |
CN111624584B (zh) * | 2020-03-20 | 2023-06-20 | 中国人民解放军火箭军工程大学 | 一种非协作目标激光诱偏距离测量系统及方法 |
CN113064194A (zh) * | 2020-08-12 | 2021-07-02 | 中国科学院微小卫星创新研究院 | 基于互联网星座的皮纳卫星通信、导航、定姿方法 |
CN112068166A (zh) * | 2020-09-14 | 2020-12-11 | 天津云遥宇航科技有限公司 | 一种用于双天线卫星的卫星定位系统 |
CN113094371A (zh) * | 2021-04-14 | 2021-07-09 | 嘉兴毕格智能科技有限公司 | 一种用户自定义坐标系的实现方法 |
CN113126022A (zh) * | 2021-04-14 | 2021-07-16 | 成都金诺信高科技有限公司 | 一种双天线定位测向方法 |
CN113126022B (zh) * | 2021-04-14 | 2023-11-03 | 成都金诺信高科技有限公司 | 一种双天线定位测向方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109061701B (zh) | 2021-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109061701B (zh) | 一种适用于星载多天线融合载波相位差分测量方法 | |
Lau et al. | An innovative deep space application of GPS technology for formation flying spacecraft | |
CN106842268B (zh) | 双gnss接收机载波相位双差整周模糊度浮点解向量估计方法 | |
US6072433A (en) | Autonomous formation flying sensor | |
US7817087B2 (en) | Method and apparatus for relative navigation using reflected GPS signals | |
US11609346B2 (en) | GNSS-based attitude determination algorithm and triple-antenna GNSS receiver for its implementation | |
CN107003386A (zh) | 一种卫星导航测姿方法和装置及无人机 | |
Kahr et al. | GPS relative navigation for the CanX-4 and CanX-5 formation-flying nanosatellites | |
Huang et al. | Improvement of orbit determination for geostationary satellites with VLBI tracking | |
CN103675861A (zh) | 一种基于星载gnss多天线的卫星自主定轨方法 | |
CN104459751B (zh) | 基于gnss反射信号的双站雷达空间目标相对导航方法 | |
CN101614802A (zh) | 一种导航卫星姿态测量方法 | |
US6594582B1 (en) | GPS compound eye attitude and navigation sensor and method | |
Medina et al. | On the Kalman filtering formulation for RTK joint positioning and attitude quaternion determination | |
Zhou et al. | Relative position determination between Chang'E-3 lander and rover using in-beam phase referencing. | |
CN106199670B (zh) | 一种基于蒙特卡洛采样的gnss单频单历元姿态确定方法 | |
CN113064195B (zh) | 一种利用多天线几何特征的高精度低计算载体测姿方法 | |
CN108205151B (zh) | 一种低成本gps单天线姿态测量方法 | |
WO2017066750A1 (en) | Triple difference formulation for formation flight | |
Sollie et al. | Pose estimation of UAVs based on INS aided by two independent low-cost GNSS receivers | |
Yang et al. | Availability analysis of GNSS signals above GNSSs constellation | |
KR102233123B1 (ko) | 소형 위성군을 위한 네비게이션 모듈 및 방법 | |
Zheng et al. | Interferometry imaging technique for accurate deep-space probe positioning | |
Fu et al. | Performance Analysis of Kinematic-to-Kinematic Relative Positioning Based on GPS/BDS-3 | |
Liu et al. | GNSS ambiguity resolution and attitude determination by leveraging relative baseline and frequency information |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |