CN105204047B - 一种卫星导航系统中观测量单个粗差的探测与修复方法 - Google Patents
一种卫星导航系统中观测量单个粗差的探测与修复方法 Download PDFInfo
- Publication number
- CN105204047B CN105204047B CN201510659018.7A CN201510659018A CN105204047B CN 105204047 B CN105204047 B CN 105204047B CN 201510659018 A CN201510659018 A CN 201510659018A CN 105204047 B CN105204047 B CN 105204047B
- Authority
- CN
- China
- Prior art keywords
- observed quantity
- error
- value
- detection
- component
- 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.)
- Active
Links
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
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
本发明公开了一种卫星导航系统中观测量单个粗差的探测与修复方法,包括以下步骤:根据随机模型和当前历元的卫星观测量,建立并求解伪距观测方程,获取观测量对应的残差;根据步骤一中的观测量和观测量对应的残差获取观测量各分量的粗差初值;根据步骤二中的粗差初值,对所述粗差初值进行迭代搜索,获取观测量各分量的粗差探测值;对步骤三中的观测量各分量的粗差探测值进行判断,获取粗差探测位置,并对相应观测量分量进行粗差修复。本发明的有益效果是,在可能的粗差范围内可以成功搜索到粗差估值,通过位置判断可以准确获得粗差的实际位置,该方法与现有方法相比,探测出的粗差探测值更加接近粗差实际值,且探测位置更加准确,可靠性高。
Description
技术领域
本发明涉及卫星导航技术领域,尤其是一种卫星导航系统中观测量单个粗差的探测与修复方法。
背景技术
目前,GNSS可为陆海空提供全天候和全球性的导航、定位和测量服务,目前已被广泛应用于交通、测绘等许多行业。由于其高精度和自动测量的特点,作为先进的测量手段和新的生产力,已经融入了国民经济建设、国防建设和社会发展的各个应用领域。
GNSS虽然具有上述优点,但由于其观测量受卫星数据质量和外界环境的影响较大,且随着GNSS系统的增多,观测数据也变得越来越多,观测数据出现粗差的风险越来越大。这些粗差如果直接参与数据处理,将会大大降低GNSS导航和定位结果应有的精度,甚至会带来破坏性影响。为了保证观测数据的质量,经常要对数据进行处理,探测并修复或剔除里面的粗差,这是GNSS导航与定位中非常重要、必须解决且最具挑战性的问题之一。
但是在粗差探测中,不论是独立观测数据,还是相关观测数据,都将面临着粗差不能被正确识别的问题。GNSS数据中包括存在多个粗差和单个粗差的情况。针对存在粗差的情况,目前常用的单个粗差的探测与修复方法主要是根据残差平方和或个别残差绝对值的大小来探测粗差,安全可靠性较差。
发明内容
本发明的目的是为克服上述现有技术的不足,提供一种观测量中的单个粗差探测可靠性较高且适合于多个卫星导航系统中的粗差探测与修复方法。
为实现上述目的,本发明采用下述技术方案:
一种卫星导航系统中观测量单个粗差的探测与修复方法,包括以下步骤:
步骤一:根据随机模型和当前历元的卫星观测量,建立并求解伪距观测方程,获取观测量对应的残差;
步骤二:根据步骤一中的观测量和观测量对应的残差获取观测量各分量的粗差初值;
步骤三:根据步骤二中的粗差初值,对所述粗差初值进行迭代搜索,获取观测量各分量的粗差探测值;
步骤四,对步骤三中的观测量各分量的粗差探测值进行判断,获取粗差探测位置,并对相应观测量分量进行粗差修复。
优选的,所述步骤一中,根据最小二乘估计法求解伪距观测方程。
优选的,所述步骤二包括以下子步骤:
S201:根据伪距观测方程,计算单位向量对应的残差;
S202:取单位向量对应的残差各分量符号,与观测量权阵P对角线对应分量平方根相乘,构成加权向量;
S203:将加权向量各分量分别与观测量对应的残差和单位向量对应的残差相乘,分别得到实际偏差和单位偏差,将实际偏差和单位偏差各分量分别相加分别得到单位残差对应的统计量和实际残差对应的统计量;
S204:将单位残差对应的统计量和实际残差对应的统计量作比值,得到粗差初值。
优选的,所述步骤三包括以下子步骤:
S301:构造搜索区间(Biask0-R,Biask0+R),其中R>0,Biask0表示第k个观测量上的粗差初值;
S302:计算搜索区间内粗差探测可能值以及粗差探测可能值对应的观测量残差平方和,并根据观测量残差平方和从粗差探测可能值中获取粗差探测值;
在搜索区间(Biask0-R,Biask0+R)内,R>0,计算观测量各分量对应的粗差探测可能值,针对观测量各分量分别进行迭代搜索,将搜索区间内残差平方和最小值对应的粗差探测可能值作为该观测量分量的粗差探测值。
选择残差平方和最小值的原因是:
观测量L的第k个分量lk上修复的粗差越接近观测量真实值,修复后的观测量所计算出的残差越小,因此观测量分量中对应的粗差探测可能值中,最接近观测量真实值的粗差对应的残差平方和必然是该观测量分量中所有的粗差探测可能值对应的残差平方和最小值,该粗差探测可能值即为粗差探测值。
所述步骤四中,对观测量各分量的粗差探测值进行判断是指根据步骤三中得到的n个粗差探测位置通过判断准则唯一得确定粗差位置。
优选的,所述步骤四中,根据粗差探测值获取探测位置的判断标准为:
第k个观测量lk含有粗差,且其大小为判断标准为同时满足绝对值准则和差值准则,即:
1.绝对值准则:绝对值大于阈值;
2.差值准则:任意一个观测量分量修复粗差后的残差平方和与第k个观测量lk修复粗差后的残差平方和的差值大于差值阈值。
优选的,所述步骤四中,根据粗差探测值和探测位置,对观测量各分量进行粗差修复具体为:在探测位置上,将对应的观测量分量减去步骤三中粗差探测探测值,实现粗差修复。
本发明的有益效果是,在可能的粗差范围内可以成功搜索到粗差估值,通过位置判断可以准确获得粗差的实际位置,该方法与现有方法相比,探测出的粗差探测值更加接近粗差实际值,且探测位置更加准确,可靠性高。
附图说明
图1是本发明流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
如图1所示,一种卫星导航系统中观测量单个粗差的探测与修复方法,包括以下步骤:
步骤一:根据随机模型和当前历元的卫星观测量,建立并求解伪距观测方程,获取观测量对应的残差;
步骤二:根据步骤一中的观测量和观测量对应的残差获取观测量各分量的粗差初值;
步骤三:根据步骤二中的粗差初值,对所述粗差初值进行迭代搜索,获取观测量各分量的粗差探测值;
步骤四:对步骤三中的观测量各分量的粗差探测值进行判断,获取粗差探测位置,并对相应观测量分量进行粗差修复。
实施例:步骤一:全球卫星导航系统GNSS的伪距观测方程为AX=L,随机模型中,权阵为P,表示为
式中,X为位置参数向量,A为系数矩阵,观测量为L,表示为L=[l1l2…lk…ln]T。根据最小二乘估计方法,得到位置参数向量X的估值,观测量L对应的残差表达为:
V=(I-A(ATPA)-1ATP)L=L-AX=e=[v1v2…vk…vn]T,式中,I为单位阵。
所述步骤二包括步骤:
S201:根据伪距观测方程,计算单位向量对应的残差;
单位向量组为bk,k=1,2,…,n,单位向量组包括n个向量,单位向量bk为bk=[00…1…0]T,表示第k个分量为1,其余分量为0;
根据单位向量组建立伪距观测方程AX=bk,k=1,2,…,n,依次计算单位向量bk,k=1,2,…,n对应的残差为:
求解单位向量bk对应的残差的目的是计算出第k个观测量lk上为一个单位粗差对应的残差。
S202:根据单位向量对应的残差偏差和权阵对角线各分量获得加权向量;
取单位向量各分量符号,与观测量权阵P对角线对应分量平方根相乘,构成加权向量Sk,表示为:
S203:根据加权向量与观测量对应的残差得到单位残差对应的统计量根据加权向量和单位向量对应的残差得到实际残差对应的统计量Qk;
将加权向量Sk各分量分别与观测量对应的残差和单位向量对应的残差相乘,分别得到实际偏差Bk和单位偏差
将实际偏差Bk和单位偏差各分量分别相加分别得到单位残差对应的统计量和实际残差对应的统计量Qk,即
S204:将单位残差对应的统计量和实际残差对应的统计量Qk作比值,得到粗差初值。
所述步骤三包括以下步骤:
S301:构造搜索区间(Biask0-R,Biask0+R)其中R>0,Biask0表示第k个观测量上的粗差初值;
S302:计算搜索区间内粗差探测可能值以及粗差探测可能值对应的观测量残差平方和,并根据观测量残差平方和从粗差探测可能值中获取粗差探测值;
进一步,在各观测量分量对应的搜索区间(Biask0-R,Biask0+R)内搜索粗差探测可能值Biass,并计算每一个粗差探测可能值修复后的观测量残差,将残差平方和最小值对应的粗差探测估计值作为粗差探测值。
以观测量分量l2为例:
设观测量为L=[l1,l2,…ln]T,l2分量上的粗差初值为Bias20=5,R=3,设定步长值为1,则搜索区间为(2,8);
该搜索区间(2,8)内粗差探测可能值为Biass,Biass∈{3,4,5,6,7},此时,根据粗差探测可能值得到修复后的观测量L=[l1,l2-Bias20,l3,…ln]T,分别计算Biass∈{3,4,5,6,7}对应的观测量的残差平方和若Bias20=5时,第2个位置的残差平方和最小,则观测量分量l2上粗差探测值为Bias2=5。
采用上述算法原因是:
以观测量分量l2为例,若观测量分量ls中粗差探测可能值为Biass,那么该粗差探测可能值Biass对应的残差平方和必然是残差平方和最小值,该残差平方和最小值对应的粗差探测可能值为观测量分量ls的粗差探测值,并用表示,对应的残差和残差平方和分别记和
步骤四,对步骤三中的观测量各分量的粗差探测值进行判断,获取粗差探测位置,并对相应观测量分量进行粗差修复。
进一步,所述步骤四中,对观测量各分量的粗差探测值进行判断是指根据步骤三中得到的n个粗差探测位置通过判断准则唯一得确定粗差位置。
设第k个观测量lk含有粗差,且其大小为判断标准为同时满足绝对值准则和差值准则,即:
1.绝对值准则:绝对值大于阈值;
2.差值准则:任意一个观测量分量修复粗差后的残差平方和与第k个观测量lk修复粗差后的残差平方和的差值大于差值阈值。
所述步骤四中,根据粗差探测值和探测位置,对观测量分量进行粗差修复指在探测位置上,将观测量分量减去粗差探测值,实现观测量分量的粗差修复。
实施例采用澳大利亚科廷大学内某一空旷区域,选择相隔4米左右为一条基线,选取2014年1月1日数据中的第2个历元作为试验数据,定位方式采用GPS卫星和北斗卫星双差联合定位,伪距类型选用“C1”,试验数据总共包含12颗卫星,其中前6颗是GPS卫星,后6颗是北斗卫星。依次在12颗卫星上人为加入粗差,分别根据本发明的方法与巴尔达法获取粗差探测值,并将本发明的方法与巴尔达法进行比较。
表1是本发明的粗差探测试验结果,表2是巴尔达法的探测试验结果。
以表1中第1列数据为例,卫星号1表示所加粗差的卫星号,4.8表示人为加入的粗差值,即粗差实际值,探测位置1表示新方法探测出的粗差位置,探测值4.74表示本发明的方法获取的粗差探测值。表2同理,探测位置为0表示巴尔达法未能探测到粗差。
表1本发明的粗差探测试验结果
表2巴尔达法粗差探测试验结果
从表1、表2可以明显看出:
一方面,与巴尔达法相比,本发明的粗差探测与修复方法可以准确探测出所有粗差的位置,并且所探测出粗差最终估计值更加接近粗差实际值,同时,对于较小的粗差,如1.0、2.3,本发明的粗差探测与修复方法也能成功探测。
另一方面,本发明的粗差探测与修复方法针对GPS卫星探测出来得粗差比针对北斗卫星探测出的粗差最终估计值更加接近实际值,原因是北斗卫星的伪距噪声水平比GPS卫星高,因此对粗差探测影响更大,所以导致探测出来的粗差最终估计值偏离粗差实际值略大。
综上所述,本发明的粗差探测与修复方法能有效探测粗差的位置及大小,验证了本发明的可靠性,所探测出粗差的大小与实际值较接近,并且对较小粗差也能成功探测。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (5)
1.一种卫星导航系统中观测量单个粗差的探测与修复方法,其特征是,包括以下步骤:
步骤一:根据随机模型和当前历元的卫星观测量,建立并求解伪距观测方程,获取观测量对应的残差;
步骤二:根据步骤一中的观测量和观测量对应的残差获取观测量各分量的粗差初值;
步骤三:根据步骤二中的粗差初值,对所述粗差初值进行迭代搜索,获取观测量各分量的粗差探测值;
步骤四,对步骤三中的观测量各分量的粗差探测值进行判断,获取粗差探测位置,并对相应观测量分量进行粗差修复;
所述步骤二包括以下子步骤:
S201:根据伪距观测方程,计算单位向量对应的残差;
S202:取单位向量对应的残差各分量符号,与观测量权阵P对角线对应分量平方根相乘,构成加权向量;
S203:将加权向量各分量分别与观测量对应的残差和单位向量对应的残差相乘,分别得到实际偏差和单位偏差,将实际偏差和单位偏差各分量分别相加分别得到单位残差对应的统计量和实际残差对应的统计量;
S204:将单位残差对应的统计量和实际残差对应的统计量作比值,得到粗差初值。
2.如权利要求1所述的卫星导航系统中观测量单个粗差的探测与修复方法,其特征是,所述步骤一中,根据最小二乘估计法求解伪距观测方程。
3.如权利要求1所述的卫星导航系统中观测量单个粗差的探测与修复方法,其特征是,所述步骤三包括以下子步骤:
S301:构造搜索区间(Biask0-R,Biask0+R),其中R>0,Biask0表示第k个观测量上的粗差初值;
S302:计算搜索区间内粗差探测可能值以及粗差探测可能值对应的观测量残差平方和,并根据观测量残差平方和从粗差探测可能值中获取粗差探测值。
4.如权利要求3所述的卫星导航系统中观测量单个粗差的探测与修复方法,其特征是,所述步骤S302中,在搜索区间(Biask0-R,Biask0+R)内,针对观测量各分量分别进行迭代搜索,将搜索区间内残差平方和最小值对应的粗差探测可能值作为该观测量分量的粗差探测值。
5.如权利要求1所述的卫星导航系统中观测量单个粗差的探测与修复方法,其特征是,所述步骤四中,第k个观测量lk含有粗差,且其大小为判断标准为同时满足绝对值准则和差值准则,即:
1.绝对值准则:绝对值大于阈值;
2.差值准则:任意一个观测量分量修复粗差后的残差平方和与第k个观测量lk修复粗差后的残差平方和的差值大于差值阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510659018.7A CN105204047B (zh) | 2015-10-13 | 2015-10-13 | 一种卫星导航系统中观测量单个粗差的探测与修复方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510659018.7A CN105204047B (zh) | 2015-10-13 | 2015-10-13 | 一种卫星导航系统中观测量单个粗差的探测与修复方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105204047A CN105204047A (zh) | 2015-12-30 |
CN105204047B true CN105204047B (zh) | 2016-04-06 |
Family
ID=54951824
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510659018.7A Active CN105204047B (zh) | 2015-10-13 | 2015-10-13 | 一种卫星导航系统中观测量单个粗差的探测与修复方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105204047B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111736188B (zh) * | 2020-03-09 | 2022-05-17 | 腾讯科技(深圳)有限公司 | 一种卫星定位方法、装置、电子设备和存储介质 |
CN111505687B (zh) * | 2020-04-17 | 2021-12-21 | 中国科学院国家授时中心 | 一种基于gps卫星导航系统的原始观测值粗差剔除方法 |
CN112799106B (zh) * | 2021-03-30 | 2021-06-29 | 腾讯科技(深圳)有限公司 | 卫星观测信息处理方法、移动终端及存储介质 |
CN113281796B (zh) * | 2021-07-23 | 2021-10-15 | 腾讯科技(深圳)有限公司 | 位置确定方法、速度确定方法、装置、设备和存储介质 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5424338B2 (ja) * | 2010-03-18 | 2014-02-26 | 日本電気株式会社 | 衛星測位システムの異常値検出装置、異常値検出方法及び異常値検出プログラム |
CN101833101B (zh) * | 2010-05-05 | 2013-06-12 | 北京航空航天大学 | 基于局域增强系统的完好性监测方法及装置 |
CN103163542B (zh) * | 2013-03-26 | 2014-12-10 | 东南大学 | 一种gnss基线解算中基于观测常量的粗差探测方法 |
-
2015
- 2015-10-13 CN CN201510659018.7A patent/CN105204047B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105204047A (zh) | 2015-12-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106168672B (zh) | 一种gnss多模单频rtk周跳探测方法及装置 | |
CN108169774B (zh) | 支持rtppp和rtk的多模gnss单频周跳探测与修复方法 | |
CN104181562B (zh) | 一种全球导航卫星系统的卫星优选和定位方法 | |
CN110007317B (zh) | 一种选星优化的高级接收机自主完好性监测方法 | |
CN105242293B (zh) | 一种全球导航卫星系统的高精度厘米级定位方法 | |
CN105204047B (zh) | 一种卫星导航系统中观测量单个粗差的探测与修复方法 | |
CN106885576B (zh) | 一种基于多点地形匹配定位的auv航迹偏差估计方法 | |
CN109521443A (zh) | 一种探测星历异常的方法 | |
CN106842236B (zh) | Gnss接收机周跳探测与修复处理装置 | |
CN110018507B (zh) | 一种基于星座间作差的组合精密单点定位方法及系统 | |
CN105158778B (zh) | 多系统联合实施载波相位差分故障卫星剔除方法及其系统 | |
CN104680008A (zh) | 一种基于多参考站的网络rtk区域大气误差建模方法 | |
CN105738927B (zh) | 顾及卫星导航系统间偏差的融合定位导航授时计算模型构建方法 | |
CN105301617A (zh) | 一种卫星导航系统中的整周模糊度有效性检验方法 | |
CN104459722B (zh) | 一种基于多余观测分量的整周模糊度可靠性检验方法 | |
KR101433908B1 (ko) | Gnss 데이터 품질 평가 방법 및 시스템 | |
CN104215977A (zh) | 一种基于卫星导航系统的精度评估方法及系统 | |
CN110531383A (zh) | 一种卫星导航定位中异常卫星剔除方法 | |
CN106066485A (zh) | 一种周跳探测方法及装置 | |
CN106547002A (zh) | 一种提高rtk定位精度的选星方法 | |
CN110879407A (zh) | 一种基于完好性风险模型的卫星导航观测量新息检测方法 | |
CN104504247B (zh) | 一种gps双卫星故障的raim方法 | |
CN105549046A (zh) | Gnss接收机周跳探测与修复处理方法 | |
CN105511481A (zh) | 一种星载定轨优化方法 | |
KR20140026902A (ko) | Gps 데이터 분석에 의한 플래그 데이터 제공 방법 및 시스템 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |