CN110907975A - 一种基于序贯最小二乘的模糊度固定方法 - Google Patents
一种基于序贯最小二乘的模糊度固定方法 Download PDFInfo
- Publication number
- CN110907975A CN110907975A CN201911285957.4A CN201911285957A CN110907975A CN 110907975 A CN110907975 A CN 110907975A CN 201911285957 A CN201911285957 A CN 201911285957A CN 110907975 A CN110907975 A CN 110907975A
- Authority
- CN
- China
- Prior art keywords
- matrix
- antenna
- ambiguity
- satellite
- measurement
- 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
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
- 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
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
本发明一种基于序贯最小二乘的模糊度固定方法,能够提高在动态环境下使用单频或双频导航卫星信号实时测出两个天线间模糊度的成功率,为后续的高精度基线测量、姿态确定等提供数据支持。在模糊度实数解估计过程中引入先验信息作为状态新息传递给滤波器,在减小设计矩阵病态性的同时加速收敛。本发明在序贯条件下应用基于测量误差建模理论和矩阵三角化分解技术的思想,实现了累加测量新息的同时不增加额外的计算负担,通过构建最小二乘形式的模糊度滤波器实现对模糊度实数解和协因数矩阵信息的计算,最后利用LAMBDA法对模糊度实数解信息进行处理完成模糊度整数解的搜索,进而解决了传统方法中导航卫星信号在动态环境下求解模糊度成功率不高的问题。
Description
技术领域
本发明一种基于序贯最小二乘的模糊度固定方法,涉及基于卫星导航系统的载波差分定位技术。
背景技术
当前基于北斗/GPS/Galileo卫星导航系统进行高精度相对定位主要通过载波相位差分技术来实现,而差分技术的关键在于能否快速可靠的在航解算(On The Fly,OTF)载波相位的双差整周模糊度(以下简称为模糊度)。对于双差方程中的模糊度而言,当不考虑其整数性质时,求得的解常称为实数解(或浮点解);当顾及其整数性质时,求得的解称为整数解(或固定解),经典的最小二乘方法只能求出模糊度的实数解,而模糊度的整数解的求取方法大致可分为如下四类:一般的最小二乘搜索法、模糊度函数法、模糊度协方差法以及直接取整法。目前,国际公认最小二乘降相关平差法(LAMBDA)是理论最严密、模糊度搜索成功率最高整周模糊度解算方法。近年来随着小卫星、无人机等相关应用市场的兴起,低成本、低功耗、小体积的单频接收机越来越多的用于相对定位(包括测姿、测向)过程中,其中单频模糊度的求解是该应用过程中的核心话题。
采用载波相位观测量过程中,由于其波长和码精度的限制,使双差观测方程的法方程(以下简称为法方程)的病态性很严重,系数矩阵求逆会不稳定,又由于观测噪声不可避免,采用现有的在航解算模糊度方法会导致模糊度实数解与准确值之间的偏差较大,即使使用LAMBDA方法也很难搜索出正确的模糊度值。
发明内容
本发明的目的在于克服现有技术的不足,提供了一种基于序贯最小二乘的模糊度固定方法,该方法能够确保模糊度实数解及其协方差阵的精度和收敛速度,而且可以减小模糊度搜索空间,从而减少了LAMBDA方法搜索的难度,从而达到提高模糊度固定速度和可靠性的目的。
本发明的上述目的通过以下技术方案实现:
一种基于序贯最小二乘的模糊度固定方法,包括以下步骤:
1)地面接收机利用位于地面不同位置的天线A和主天线B捕获跟踪卫星信号;
其中,为第p颗卫星相对于主天线B的方向向量,p∈[1,P];为接收机测量得到的第j颗卫星和第k颗卫星到天线A和主天线B的伪距双差测量量,所述伪距双差测量量根据伪距测量量确定;为测量得到的第j颗卫星和第k颗卫星到天线A和主天线B的载波双差测量量,所述载波双差测量量根据载波测量量确定,j∈[1,P],k∈[1,P],且j≠k;为双差测量的线性化初值; 为天线A与卫星j的空间距离,为天线A与卫星k的空间距离,为天线A与卫星k的空间距离,为天线A与卫星k的空间距离,均由系统外部直接输入获得;
Z=A·X;
获得测量方程的协因数矩阵∑Z,具体为:
构建测量方程的误差模型:
3)对矩阵∑Z使用乔里斯基分解方法,求解设计矩阵R,使矩阵R和矩阵∑Z两者之间的数学关系如下:
∑Z=(R-1)(R-1)T
4)对步骤2)中测量方程Z=A·X的等号两边同时左乘矩阵R,整理可得满足最小二乘形式的测量方程:
本发明与现有技术相比有益效果为:
1)本发明在模糊度滤波计算过程中,将模糊度随时间不变的特性作为滤波器的先验信息,降低了滤波方程的奇异性,提高了载波相位模糊度实数解以及实数解协因数阵的精度和收敛速度,能够更快地完成载波相位模糊度的分解、搜索以及整数估计,实现快速的高精度相对定位。
2)本发明在模糊度实数解及其协因数阵求解过程中,应用基于测量误差建模理论和矩阵三角化分解技术的思想,实现了累加测量新息的同时不增加额外的计算负担,避免了传统方法中高阶矩阵求逆运算带来的庞大计算量。
附图说明
图1为GNSS接收机测量示意图;
图2为本发明中应用的双差观测值的示意图。
具体实施方式
本发明提供的一种基于序贯最小二乘的模糊度固定方法,可以应用于基于码分多址的卫星导航系统,如北斗B1、B2、B3频点信号,GPS L1、L2频点信号,以及Galileo E1、E5、E6频点信号,进行单频、双频、三频形式的模糊度在航解算。
其中,利用载波相位测量值进行相对定位时,如果采用LAMBDA方法进行整周模糊度解算,则需要满足以下几点要求:一、需要确保模糊度实数解的稳定性和计算效率,也就是要尽可能减小系数矩阵发生奇异的可能性,并降低计算矩阵的维度;二、需要确保模糊度实数解及其协因数阵的精度和收敛速度。本发明采用如下措施达到以上两点要求:首先在模糊度实数解求解过程中,引入了双差伪距测量信息、双差载波相位测量信息,在增加观测方程个数的基础上可以有效降低模糊度滤波器中系数矩阵发生奇异性的可能,有效增强滤波器抗发散的能力;另外本发明在每个历元根据模糊度的不变特性通过引入前一历元的模糊度实数解信息作为先验信息,使之作为一个“伪测量”辅助解算当前模糊度的实数解,可以有效提高模糊度实数解及其协因数阵的精度和收敛速度。针对高阶矩阵的数学计算方面引入基于测量误差建模理论和矩阵三角化分解技术的思想构件模糊度分解技术,规避了传统方法中模糊度实数解计算过程中涉及的高阶矩阵求逆运算,提高了滤波器的计算效率。
采用以上技术措施后,相对于现有技术来说,本发明计算得到的模糊度搜索空间和产生的计算负担会明显减小,进而减小了LAMBDA方法搜索的难度。
本发明的基于序贯最小二乘的模糊度搜索方法,其原理推导过程如下:
如图1所示,GNSS接收机测量的是导航卫星与接收天线在空间通视方向上的信号传输时延,本质上体现的是两者间相对运动状态在通视方向的投影。由于GNSS接收机导航解算结果存在一定程度的误差,其原始观测量的测量值与理论真值之间存在以下关系:
其中,理论真值代表卫星Sj与天线A的真实几何距离;卫星j与天线A之间测量的几何距离Sj为卫星j的位置坐标,一般由广播星历计算得到;为天线A的初始位置坐标,一般通过导航解算得到;卫星j相对于天线A的方向向量 为天线A初始位置的一阶修正值。
同理可得,主天线B的测量值与理论真值的关系如下所示:
如图2所示,卫星导航系统进行载波相位差分定位时,接收机利用设置在不同位置的两个天线接收卫星信号,分别为天线A和主天线B,对应的测量接收机分别为接收机A和接收机B。采用站际-星间形式的双差模型的双差真值具体表示为
GNSS接收机的测量分为伪距和载波相位两种,两者的主要区别在于载波相位测量存在模糊度。
假设,在某一历元时刻,GNSS接收机观测P+1颗卫星时,选取其中一颗作为主星,可建立P个载波相位双差方程和P个伪距双差方程,即:
其中:伪距双差测量量 和分别为接收机A测量的卫星j和卫星k伪距测量值,和分别为接收机B测量的卫星j和卫星k伪距测量值;载波相位双差测量量 为,和表示接收机A测量的卫星j和卫星k的载波相位测量值,和表示接收机B测量的卫星j和卫星k的载波相位测量值;和分别表示伪距和载波相位双差测量量的噪声;λ为载波相位的波长;为双差模糊度,具有整数特性。
将式(5)代入式(6),可得:
为了简化式(8),设定如下:
由于在没有周跳的情况下,整周模糊度具有稳定不变的特性,因此滤波器的状态空间的历元更新可以引入整周模糊度的均方根信息,即:
Z=A·X (11)
由于双差观测方程具有高度耦合特性,式(10)的设计矩阵A属于稀疏病态阵,对矩阵A的求逆操作需要额外增加大量的计算负担,且计算精度有限。本文采用基于测量误差建模理论和矩阵三角化分解技术避免上述计算负担的增加,提高计算精度。下文对测量误差的建模与矩阵三角化分解技术进行详细描述。
普遍认为,站际-星间双差形式的伪距、载波观测方程的协因数矩阵形式∑P和∑Φ如下所示:
结合公式(11)、(13),针对本发明的测量方程可以完整描述为:
Z=A·X,∑Z (14)
由于∑P、∑Φ、均为非对角矩阵,因此测量方程(14)的待估参数使用传统的最小二乘算法进行计算是较为困难且计算负担较大的。最小二乘算法只有在针对零均值高斯白噪声的测量问题才能取得最优解,显而易见公式(14)不属于这种情况。通过公式(13)的推导可以得知,测量方程的协因数阵为正定对称矩阵,因此使用乔里斯基分解的方法对∑Z进行处理,即
∑Z=(R-1)(R-1)T (15)
然后公式(14)的等号两边同时左乘矩阵R,公式(14)变化为:
R·Z=R·A·X,∑R·(R-1)(R-1)T·RT (16)
其中∑R·(R-1)(R-1)T·RT=E,新的观测方程(16)的协因数矩阵为单位阵,满足对其采用最小二乘方法进行参数估计。
为了进一步减小计算量,对公式(16)的做进一步的整理可得:
根据以上的理论推导过程可知,本发明的提高单频模糊度在航计算成功率的方法,包括如下步骤:
1)地面接收机利用位于地面不同位置的天线A和主天线B捕获跟踪卫星信号;
其中,为第p颗卫星相对于主天线B的方向向量,p∈[1,P];为接收机测量得到的第j颗卫星和第k颗卫星到天线A和主天线B的伪距双差测量量,所述伪距双差测量量根据伪距测量量确定;为测量得到的第j颗卫星和第k颗卫星到天线A和主天线B的载波双差测量量,所述载波双差测量量根据载波测量量确定,j∈[1,P],k∈[1,P],且j≠k;为双差测量的线性化初值; 为天线A与卫星j的空间距离,为天线A与卫星k的空间距离,为天线A与卫星k的空间距离,为天线A与卫星k的空间距离,均由系统外部直接输入获得;
Z=A·X;
获得测量方程的协因数矩阵∑Z,具体为:
构建测量方程的误差模型:
4)对矩阵∑Z使用乔里斯基分解方法,求解设计矩阵R,使矩阵R和矩阵∑Z两者之间的数学关系如下:
∑Z=(R-1)(R-1)T
5)对步骤4)中测量方程Z=A·X的等号两边同时左乘矩阵R,整理可得满足最小二乘形式的测量方程:
以上所述,仅为本发明一个具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (1)
1.一种基于序贯最小二乘的模糊度固定方法,其特征在于,包括如下步骤:
1)地面接收机利用位于地面不同位置的天线A和主天线B捕获跟踪卫星信号;
其中,为第p颗卫星相对于主天线B的方向向量,p∈[1,P];为接收机测量得到的第j颗卫星和第k颗卫星到天线A和主天线B的伪距双差测量量,所述伪距双差测量量根据伪距测量量确定;为测量得到的第j颗卫星和第k颗卫星到天线A和主天线B的载波双差测量量,所述载波双差测量量根据载波测量量确定,j∈[1,P],k∈[1,P],且j≠k;为双差测量的线性化初值; 为天线A与卫星j的空间距离,为天线A与卫星k的空间距离,为天线A与卫星k的空间距离,为天线A与卫星k的空间距离,均由系统外部直接输入获得;
Z=A·X;
获得测量方程的协因数矩阵∑Z,具体为:
构建测量方程的误差模型:
4)对矩阵∑Z使用乔里斯基分解方法,求解设计矩阵R,使矩阵R和矩阵∑Z两者之间的数学关系如下:
∑Z=(R-1)(R-1)T
5)对步骤4)中测量方程Z=A·X的等号两边同时左乘矩阵R,整理可得满足最小二乘形式的测量方程:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911285957.4A CN110907975B (zh) | 2019-12-13 | 2019-12-13 | 一种基于序贯最小二乘的模糊度固定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911285957.4A CN110907975B (zh) | 2019-12-13 | 2019-12-13 | 一种基于序贯最小二乘的模糊度固定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110907975A true CN110907975A (zh) | 2020-03-24 |
CN110907975B CN110907975B (zh) | 2021-10-01 |
Family
ID=69824445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911285957.4A Active CN110907975B (zh) | 2019-12-13 | 2019-12-13 | 一种基于序贯最小二乘的模糊度固定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110907975B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111381266A (zh) * | 2020-05-12 | 2020-07-07 | 泉州信息工程学院 | 整周模糊度确定方法、装置、电子设备及计算机可读介质 |
CN112462396A (zh) * | 2020-11-20 | 2021-03-09 | 中国人民解放军战略支援部队信息工程大学 | 一种高采样率导航卫星钟差的实时并行确定方法 |
CN113296138A (zh) * | 2021-04-25 | 2021-08-24 | 北京遥测技术研究所 | 一种陆基导航高精度定位方法 |
CN115235513A (zh) * | 2022-09-15 | 2022-10-25 | 中国船舶重工集团公司第七0七研究所 | 一种基于卫导伪距和伪距率的惯导校正方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103487819A (zh) * | 2013-09-06 | 2014-01-01 | 哈尔滨工程大学 | 一种基于升序排列优化的整周模糊度去相关方法 |
US20140240172A1 (en) * | 2009-11-17 | 2014-08-28 | Topcon Positioning Systems, Inc. | Detection and Correction of Anomalous Measurements and Ambiguity Resolution in a Global Navigation Satellite System Receiver |
CN105182378A (zh) * | 2015-07-20 | 2015-12-23 | 武汉大学 | 一种lll模糊度降相关算法 |
CN106289154A (zh) * | 2016-07-19 | 2017-01-04 | 中国科学院重庆绿色智能技术研究院 | 一种通信基站天线方位角监测方法 |
CN108254774A (zh) * | 2018-03-29 | 2018-07-06 | 千寻位置网络有限公司 | 基于gnss多频信号的单基站长距离实时定位方法 |
CN108508467A (zh) * | 2018-03-07 | 2018-09-07 | 上海海积信息科技股份有限公司 | 卫星导航系统静态定位方法 |
CN108802781A (zh) * | 2018-05-03 | 2018-11-13 | 广州市中海达测绘仪器有限公司 | 整周模糊度的获取方法 |
CN109085628A (zh) * | 2018-08-27 | 2018-12-25 | 桂林电子科技大学 | 一种整周模糊度的固定方法及系统 |
-
2019
- 2019-12-13 CN CN201911285957.4A patent/CN110907975B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140240172A1 (en) * | 2009-11-17 | 2014-08-28 | Topcon Positioning Systems, Inc. | Detection and Correction of Anomalous Measurements and Ambiguity Resolution in a Global Navigation Satellite System Receiver |
CN103487819A (zh) * | 2013-09-06 | 2014-01-01 | 哈尔滨工程大学 | 一种基于升序排列优化的整周模糊度去相关方法 |
CN105182378A (zh) * | 2015-07-20 | 2015-12-23 | 武汉大学 | 一种lll模糊度降相关算法 |
CN106289154A (zh) * | 2016-07-19 | 2017-01-04 | 中国科学院重庆绿色智能技术研究院 | 一种通信基站天线方位角监测方法 |
CN108508467A (zh) * | 2018-03-07 | 2018-09-07 | 上海海积信息科技股份有限公司 | 卫星导航系统静态定位方法 |
CN108254774A (zh) * | 2018-03-29 | 2018-07-06 | 千寻位置网络有限公司 | 基于gnss多频信号的单基站长距离实时定位方法 |
CN108802781A (zh) * | 2018-05-03 | 2018-11-13 | 广州市中海达测绘仪器有限公司 | 整周模糊度的获取方法 |
CN109085628A (zh) * | 2018-08-27 | 2018-12-25 | 桂林电子科技大学 | 一种整周模糊度的固定方法及系统 |
Non-Patent Citations (4)
Title |
---|
刘立龙 等: "单频GPS整周模糊度动态快速求解的研究", 《武汉大学学报·信息科学版》 * |
刘立龙: "动态对动态GPS高精度定位理论及其应用研究", 《中国优秀博硕士学位论文全文数据库 (博士) 基础科学辑》 * |
沈雪峰 等: "一种序贯平差中长基线模糊度的解算方法", 《测绘科学》 * |
石小丽: "载波相位测量的定位解算算法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111381266A (zh) * | 2020-05-12 | 2020-07-07 | 泉州信息工程学院 | 整周模糊度确定方法、装置、电子设备及计算机可读介质 |
CN112462396A (zh) * | 2020-11-20 | 2021-03-09 | 中国人民解放军战略支援部队信息工程大学 | 一种高采样率导航卫星钟差的实时并行确定方法 |
CN112462396B (zh) * | 2020-11-20 | 2023-06-16 | 中国人民解放军战略支援部队信息工程大学 | 一种高采样率导航卫星钟差的实时并行确定方法 |
CN113296138A (zh) * | 2021-04-25 | 2021-08-24 | 北京遥测技术研究所 | 一种陆基导航高精度定位方法 |
CN113296138B (zh) * | 2021-04-25 | 2022-12-13 | 北京遥测技术研究所 | 一种陆基导航高精度定位方法 |
CN115235513A (zh) * | 2022-09-15 | 2022-10-25 | 中国船舶重工集团公司第七0七研究所 | 一种基于卫导伪距和伪距率的惯导校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110907975B (zh) | 2021-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110907975B (zh) | 一种基于序贯最小二乘的模糊度固定方法 | |
CN111239787B (zh) | 一种集群自主协同中的gnss动态卡尔曼滤波方法 | |
Giorgi et al. | Instantaneous ambiguity resolution in Global-Navigation-Satellite-System-based attitude determination applications: A multivariate constrained approach | |
CN109541663B (zh) | 一种gnss定位测姿多路径误差的纠正方法 | |
CN113466903B (zh) | 一种顾及观测值系统误差的部分模糊度固定算法 | |
CN105445772A (zh) | 多gnss天线组合平台位姿一体化的确定装置及方法 | |
CN111399020A (zh) | 一种定向测姿系统及方法 | |
US20140197988A1 (en) | Method of estimating a quantity associated with a receiver system | |
CN110954932B (zh) | 一种三频rtk单历元模糊度快速搜索方法 | |
CN110133702B (zh) | 一种基于正交变换的姿态测量方法和设备 | |
CN110646822A (zh) | 一种基于惯导辅助的整周模糊度Kalman滤波算法 | |
US11821999B2 (en) | Attitude determination based on global navigation satellite system information | |
CN115755115A (zh) | 基于gnss对流层层析技术的ppp改善方法 | |
CN109752747B (zh) | 差分数据确定方法、装置、服务器及存储介质 | |
CN106199670B (zh) | 一种基于蒙特卡洛采样的gnss单频单历元姿态确定方法 | |
US20110181462A1 (en) | System and Method for Positioning with GNSS Using Multiple Integer Candidates | |
CN116106953B (zh) | 一种gnss相对定位半参数估计的时间差分载波相位增强方法 | |
CN113064195A (zh) | 一种利用多天线几何特征的高精度低计算载体测姿方法 | |
CN112230254A (zh) | 一种gps载波相位多径误差的校正方法及装置 | |
Wu et al. | Particle filter-based real-time phase line bias estimation for GNSS-based attitude determination with common-clock receivers | |
Ning et al. | A modified geometry-and ionospheric-free combination for static three-carrier ambiguity resolution | |
Rapoport | General purpose kinematic/static GPS/GLONASS postprocessing engine | |
Jia et al. | Mitigating multipath errors using semi-parametric models for high precision static positioning | |
Teunissen et al. | Low-complexity instantaneous GNSS attitude determination with multiple low-cost antennas | |
CN111158031A (zh) | 电离层延迟的提取方法及装置、用户及服务端、定位系统 |
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 |