CN103926597B - 一种基于北斗rdss双向通信功能的多路径检测方法 - Google Patents
一种基于北斗rdss双向通信功能的多路径检测方法 Download PDFInfo
- Publication number
- CN103926597B CN103926597B CN201410172299.9A CN201410172299A CN103926597B CN 103926597 B CN103926597 B CN 103926597B CN 201410172299 A CN201410172299 A CN 201410172299A CN 103926597 B CN103926597 B CN 103926597B
- Authority
- CN
- China
- Prior art keywords
- satellite
- gdop
- signal
- carrier
- noise ratio
- 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.)
- Expired - Fee Related
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/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/22—Multipath-related issues
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
本发明涉及的是一种基于北斗RDSS双向通信功能的多路径检测方法。本发明在半径为30km的目标区域内授权合作基准站,收集各颗卫星信号的原始观测量,通过北斗RDSS上行信道回传给北斗地面控制站,在控制站对原始观测量做统计分析处理,获得原始观测量载噪比的阈值。本发明充分抑制受多路径影响比较严重的卫星信号,从而提高用户定位的精度。
Description
技术领域
本发明涉及的是一种基于北斗RDSS双向通信功能的多路径检测方法。
背景技术
基于卫星导航定位是通过接收机天线接收卫星发射的电磁波信号并跟踪卫星而实现的。很多情况下,到达接收机天线的信号为直射信号和反射信号的叠加信号,甚至只有反射信号,这就是多路径效应产生的根源。
目前在数据处理中还难以准确估计多路径效应并消除其影响,为减小或消除多路径误差的影响,目前广泛采取的方法包括:(1)天线设置在远离反射体的地点,但是这种方法限制了卫星导航系统的使用环境,适用性无法得到满足;(2)改进天线与底座组合的性能,这种方法需要极大提高接收机硬件成本,普通用户难以承受;(3)设计更为先进的码和载波跟踪环路,这是以复杂的实现算法为代价来实现的,而且所设计方法的鲁棒性还有待进一步验证。综上所述,有必要研究一种新方法,对卫星信号进行多路径检测与排除,以充分抑制受多路径影响比较严重的卫星信号。
卫星信号发生多路径效应时,单路反射波信号的功率小于直射波信号的功率,直接表现之一就是接收机接收到的卫星信号的载噪比值降低。多路反射信号的组合虽然不能确定使得载噪比值降低,但是也将服从某种统计分布模型[1]。具体而言,当某颗卫星信号载噪比值高于正常值时,意味着该卫星信号受多路径影响较小或不受多路径影响;而当该卫星信号载噪比值明显低于正常值时,意味着该卫星信号受多路径影响的可能性较大。因此,卫星信号的载噪比值可以作为衡量卫星信号受多路径影响的指标之一。
北斗二代卫星导航系统集卫星无线电测定业务(Radio DeterminationSatellite System,RDSS)和无线电导航卫星业务(Radio Navigation Satellite System,RNSS)两种体制于一身,因此具有RDSS双向通信和RNSS导航定位两种功能。通过RDSS双向通信功能,可以实现用户与地面控制站之间的双向信息传递,即用户可以通过上行信道向地面控制站传送信息,地面控制站可以通过下行信道向用户播发信息。利用北斗RDSS双向通信功能,收集目标区域内的观测信息,通过在北斗地面控制站的数据统计分析结果,为用户接收机提供参考信息,用户接收机根据该信息做出判断,对在当前解算历元中对受多路径影响的卫星信号进行相应处理,以充分抑制多路径误差的影响。
发明内容
本发明的目的在于提供一种提高用户定位精度的基于北斗RDSS双向通信功能的多路径检测方法。
本发明的目的是这样实现的:
(1)在半径为30km的目标区域内授权合作基准站,收集各颗卫星信号的原始观测量,通过北斗RDSS上行信道回传给北斗地面控制站,在控制站对原始观测量做统计分析处理,获得原始观测量载噪比的阈值;Si表示i号卫星的载噪比值,则n个基准站观测到的i号卫星的载噪比值S1i、S2i、S3i、…、Sni是Si的观测值,Si服从正态分布,即Si~N(μi,σi 2),其中μi为载噪比均值,为载噪比方差,
(2)在获得i号卫星载噪比值分布函数后,为均值μi设定置信水平α,对应分位数为Thi,从而获得置信区间[μi-Thi,μi+Thi],将置信区间的下限μi-Thi设定为载噪比阈值,在有效历元周期内,若用户接收机观测到i号卫星,将卫星信号载噪比与阈值作比较,当载噪比值大于或等于阈值时,卫星信号基本不受多路径影响,为合格信号;否则,卫星信号受多路径影响比较严重,为不合格信号;
(3)获得各颗卫星信号的载噪比阈值后,通过北斗RDSS下行信道向目标区域内的用户接收机播发对应每颗卫星的载噪比阈值信息;
(4)用户接收机接收到载噪比阈值信息后,根据定位精度要求,与对应卫星的载噪比阈值比较:
(4.1)判断历元接收到的所有共视卫星信号数目是否少于4颗,若少于4颗,则无法对该历元实现定位解算;若不少于4颗,则执行下一步骤;
(4.2)若所有共视卫星数目n只有4颗,执行步骤(4.4);否则,执行步骤(4.3);
(4.3)若所有共视卫星数目n大于4颗,使用所有的卫星信号进行加权最小二乘法定位解算,以获得几何精度因子GDOPo值,根据GDOPo值设定GDOP值损失门限值,记为ThGDOP,GDOP值损失门限值ThGDOP由所需导航性能定义的完好性参数推导获取,GDOPo值为:
对于观测到的n颗卫星,可得观测量的线性模型为:
式中,ρ表示伪距观测量,α、β、γ为对应卫星和用户之间的方向余弦矢量,xu、yu和zu为规定坐标系下的用户坐标,c为光速时间常量,δt代表用户钟差项,ε为观测量噪声,线型模型的等效矩阵形式为:
Y=GX+ε;
得到权系数矩阵为:
H=(GTG)-1;
从而获得定位几何精度因子为:
式中,trace表示取矩阵的迹,将接收到的卫星信号的载噪比与相应卫星信号的载噪比阈值作比较,检测出合格信号卫星集Qg和不合格信号卫星集,将所有不合格卫星集卫星所引入的GDOP损失与GDOP门限值ThGDOP进行比较,第i颗卫星引入的GDOP损失GDOPloss为:
GDOPloss=GDOPi-GDOPo;
式中,GDOPi为排除第i颗卫星后的GDOP值,比较GDOPloss与ThGDOP,若GDOPloss≤ThGDOP,即排除掉第i颗卫星后,GDOP损失GDOPloss在可接受的范围之内,纳入卫星集Qb,设有l颗;否则,将其纳入卫星集Qd,设有m颗;若l+m≥4,则直接将Qg和Qd中的卫星加入到最终参与定位的卫星集,对Qd中的卫星进行降权处理,执行步骤(4.4);若l+m<4,需将Qb中的卫星进行排序处理,并从中选取GDOP损失最小的4-(l+m)颗卫星,对Qd中的卫星和从Qb中挑选出的卫星进行降权处理,并加入到最终参与定位的卫星集;
(4.4)利用所挑选出的卫星信号进行加权最小二乘法定位解算,获得用户的位置用户位置的加权最小二乘解为:
X=(GTWG)-1GTWY;
式中,W为加权对角矩阵:
其中
式中,分别是卫星集Qg,l颗卫星、卫星集Qd,m颗卫星和从卫星集Qb中挑选出的4-(l+m)颗卫星的加权系数;卫星对GDOP值的贡献越大,对应的加权系数值也越大,应有加权系数值;最终检测并抑制用户端受多路径影响的卫星信号,提高用户定位的精度。
本发明的有益效果在于:
本发明利用接收到的阈值信息,用户接收机通过比较所接收卫星信号的原始观测信息进行检测处理,对当前历元中对载噪比值小于阈值的卫星信号做剔除或降权处理,以充分抑制受多路径影响比较严重的卫星信号。利用接收到的阈值信息,用户接收机通过比较所接收卫星信号的原始观测信息进行检测处理,对当前历元中对载噪比值小于阈值的卫星信号做剔除或降权处理,以充分抑制受多路径影响比较严重的卫星信号,从而提高用户定位的精度。
附图说明
图1为本设计所提出的检测与抑制多路径效应定位方法执行流程图。
具体实施方式
下面结合附图对本发明做进一步描述。
本发明利用北斗RDSS通信功能实现收集目标区域内各授权合作基准站原始观测量的方法;基于以上内容在北斗地面控制站实现目标区域内的原始观测量的质量统计分析,获得针对各颗卫星信号在某些置信水平下的原始观测量载噪比阈值信息,并通过北斗RDSS通信功能向目标区域内播发;用户根据接收到的载噪比阈值信息,并通过几何精度因子损失准则对卫星信号做剔除或降权处理。
在半径为30km的目标区域内授权合作多个基准站,利用这些基准站收集各颗卫星信号的原始观测量,并将这些原始观测量通过北斗RDSS上行信道回传给北斗地面控制站,在控制站对这些原始观测量做统计分析处理,获得原始观测量载噪比的阈值信息。
假设某时刻地面控制站接收到的原始观测量中,在某一历元,这n个基准站都观测到了i号卫星,令Si表示i号卫星的载噪比值,那么可以认为这n个基准站观测到的i号卫星的载噪比值S1i、S2i、S3i、…、Sni是Si的观测值。根据经验,可以认为Si近似服从正态分布,即Si~N(μi,σi 2),其中μi为载噪比均值,为载噪比方差,两者可按以下两式计算获得。
在获得i号卫星载噪比值分布函数后,为均值μi设定置信水平α,对应分位数为Thi,从而获得置信区间[μi-Thi,μi+Thi],可将置信区间的下限μi-Thi设定为载噪比阈值。该阈值只适用于与该历元相隔较近的历元,历元周期的设定与卫星几何变化特性相关,可设定为30s。在有效历元周期内,若用户接收机观测到i号卫星,将该卫星信号载噪比与其阈值作比较,当载噪比值大于或等于阈值时,认为该卫星信号基本不受多路径影响,称之为合格信号;反之,则认为该卫星信号受多路径影响比较严重,称之为不合格信号。
在获得各颗卫星信号的载噪比阈值后,将这些信息通过北斗RDSS下行信道向目标区域播发,这样,目标区域内的用户接收机就可以接收到对应每颗卫星的载噪比阈值信息。
用户接收机接收到这些信息后,根据自己的定位精度要求,与对应卫星的载噪比阈值比较,在单次解算周期内,按照如下步骤执行定位解算。
步骤一:判断该历元接收到的所有共视卫星信号数目是否少于4颗,若少于4颗,则无法对该历元实现定位解算;若不少于4颗,则执行以下步骤。
步骤二:若所有共视卫星数目n只有4颗,跳至步骤四;否则,跳至步骤三。
步骤三:所有共视卫星数目n大于4颗。首先,使用所有的卫星信号进行加权最小二乘法定位解算,以获得几何精度因子GDOPo值,根据GDOPo人为设定一个GDOP值损失门限值,记为ThGDOP。GDOP值损失门限值ThGDOP也可以由所需导航性能定义的完好性参数推导获取,具体可参考文献[2]。GDOPo值计算方法如下:
对于观测到的n颗卫星,可得观测量的线性模型为:
式中,ρ表示伪距观测量,α、β、γ为对应卫星和用户之间的方向余弦矢量,xu、yu和zu为规定坐标系下的用户坐标,c为光速时间常量,δt代表用户钟差项,ε为观测量噪声。该线型模型的等效矩阵形式为:
Y=GX+ε (4)
根据式(4),得到权系数矩阵为:
H=(GTG)-1 (5)
从而获得定位几何精度因子为:
式中,trace表示取矩阵的迹。将接收到的卫星信号的载噪比与相应卫星信号的载噪比阈值作比较,检测出合格信号卫星集Qg(假设有l颗)和不合格信号卫星集。将所有不合格卫星集卫星所引入的GDOP损失与GDOP门限值ThGDOP进行比较,第i颗卫星引入的GDOP损失GDOPloss可表示为:
GDOPloss=GDOPi-GDOPo (7)
式中,GDOPi为排除第i颗卫星后的GDOP值。比较GDOPloss与ThGDOP,若GDOPloss≤ThGDOP,即排除掉第i颗卫星后,GDOP损失(GDOPloss)在可接受的范围之内,那么将其纳入卫星集Qb;否则,将其纳入卫星集Qd(假设有m颗)。若l+m≥4,则直接将Qg和Qd中的卫星加入到最终参与定位的卫星集,但是对Qd中的卫星进行降权处理,进行步骤四;若l+m<4,需将Qb中的卫星进行排序处理,并从中选取GDOP损失最小的4-(l+m)颗卫星,对Qd中的卫星和从Qb中挑选出的卫星进行降权处理,并加入到最终参与定位的卫星集。
步骤四:利用所挑选出的卫星信号进行加权最小二乘法定位解算,获得用户的位置。根据式(4)可得用户位置的加权最小二乘解为:
X=(GTWG)-1GTWY (8)
式中,W为加权对角矩阵,其形式如下:
其中
式中,分别是卫星集Qg(l颗卫星)、卫星集Qd(m颗卫星)和从卫星集Qb中挑选出的4-(l+m)颗卫星的加权系数。卫星对GDOP值的贡献越大,其对应的加权系数值也越大,根据上文可知,应有
通过以上处理,就可以检测并充分抑制用户端受多路径影响比较严重的卫星信号,从而提高用户定位的精度。
本设计所提出的检测与抑制多路径效应方法实现流程图如图1所示。
Claims (1)
1.一种基于北斗RDSS双向通信功能的多路径检测方法,其特征在于:
(1)在半径为30km的目标区域内授权合作基准站,收集各颗卫星信号的原始观测量,通过北斗RDSS上行信道回传给北斗地面控制站,在控制站对原始观测量做统计分析处理,获得原始观测量载噪比的阈值;Si表示i号卫星的载噪比值,则n个基准站观测到的i号卫星的载噪比值S1i、S2i、S3i、…、Sni是Si的观测值,Si服从正态分布,即其中μi为载噪比均值,为载噪比方差,
(2)在获得i号卫星载噪比值分布函数后,为均值μi设定置信水平α,对应分位数为Thi,从而获得置信区间[μi-Thi,μi+Thi],将置信区间的下限μi-Thi设定为载噪比阈值,在有效历元周期内,若用户接收机观测到i号卫星,将卫星信号载噪比与阈值作比较,当载噪比值大于或等于阈值时,卫星信号基本不受多路径影响,为合格信号;否则,卫星信号受多路径影响比较严重,为不合格信号;
(3)获得各颗卫星信号的载噪比阈值后,通过北斗RDSS下行信道向目标区域内的用户接收机播发对应每颗卫星的载噪比阈值信息;
(4)用户接收机接收到载噪比阈值信息后,根据定位精度要求,与对应卫星的载噪比阈值比较:
(4.1)判断历元接收到的所有共视卫星信号数目是否少于4颗,若少于4颗,则无法对该历元实现定位解算;若不少于4颗,则执行下一步骤;
(4.2)若所有共视卫星数目n只有4颗,执行步骤(4.4);否则,执行步骤(4.3);
(4.3)若所有共视卫星数目n大于4颗,使用所有的卫星信号进行加权最小二乘法定位解算,以获得几何精度因子GDOPo值,根据GDOPo值设定GDOP值损失门限值,记为ThGDOP,GDOP值损失门限值ThGDOP由所需导航性能定义的完好性参数推导获取,GDOPo值为:
对于观测到的n颗卫星,可得观测量的线性模型为:
式中,ρ表示伪距观测量,α、β、γ为对应卫星和用户之间的方向余弦矢量,xu、yu和zu为规定坐标系下的用户坐标,c为光速时间常量,δt代表用户钟差项,ε为观测量噪声,线型模型的等效矩阵形式为:
Y=GX+ε;
得到权系数矩阵为:
H=(GTG)-1;
从而获得定位几何精度因子为:
式中,trace表示取矩阵的迹,将接收到的卫星信号的载噪比与相应卫星信号的载噪比阈值作比较,检测出合格信号卫星集Qg和不合格信号卫星集,将所有不合格卫星集卫星所引入的GDOP损失与GDOP门限值ThGDOP进行比较,第i颗卫星引入的GDOP损失GDOPloss为:
GDOPloss=GDOPi-GDOPo;
式中,GDOPi为排除第i颗卫星后的GDOP值,比较GDOPloss与ThGDOP,若GDOPloss≤ThGDOP,即排除掉第i颗卫星后,GDOP损失GDOPloss在可接受的范围之内,纳入卫星集Qb;否则,将其纳入卫星集Qd,设有m颗;若l+m≥4,则直接将Qg和Qd中的卫星加入到最终参与定位的卫星集,对Qd中的卫星进行降权处理,执行步骤(4.4);若l+m<4,需将Qb中的卫星进行排序处理,并从中选取GDOP损失最小的4-(l+m)颗卫星,对Qd中的卫星和从Qb中挑选出的卫星进行降权处理,并加入到最终参与定位的卫星集;
(4.4)利用所挑选出的卫星信号进行加权最小二乘法定位解算,获得用户的位置用户位置的加权最小二乘解为:
X=(GTWG)-1GTWY;
式中,W为加权对角矩阵:
其中
式中,分别是卫星集Qg,l颗卫星、卫星集Qd,m颗卫星和从卫星集Qb中挑选出的4-(l+m)颗卫星的加权系数;卫星对GDOP值的贡献越大,对应的加权系数值也越大,应有加权系数值;最终检测并抑制用户端受多路径影响的卫星信号,提高用户定位的精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410172299.9A CN103926597B (zh) | 2014-04-25 | 2014-04-25 | 一种基于北斗rdss双向通信功能的多路径检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410172299.9A CN103926597B (zh) | 2014-04-25 | 2014-04-25 | 一种基于北斗rdss双向通信功能的多路径检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103926597A CN103926597A (zh) | 2014-07-16 |
CN103926597B true CN103926597B (zh) | 2016-08-17 |
Family
ID=51144885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410172299.9A Expired - Fee Related CN103926597B (zh) | 2014-04-25 | 2014-04-25 | 一种基于北斗rdss双向通信功能的多路径检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103926597B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106487547B (zh) * | 2015-08-26 | 2019-06-21 | 中国移动通信集团公司 | 一种评估网络性能指标变化趋势的方法及装置 |
CN108627857B (zh) * | 2017-03-17 | 2021-10-26 | 展讯通信(上海)有限公司 | 多径检测方法、装置及gnss接收机 |
CN108279425A (zh) * | 2018-01-29 | 2018-07-13 | 鄢名扬 | 一种多频观测过程中多径误差的修正方法 |
CN111798648B (zh) * | 2019-04-08 | 2022-04-08 | 千寻位置网络有限公司 | 智能报警的方法及装置、报警平台及终端 |
CN112904390B (zh) * | 2021-01-14 | 2022-03-25 | 腾讯科技(深圳)有限公司 | 定位方法、装置、计算机设备和存储介质 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7095370B1 (en) * | 2005-10-27 | 2006-08-22 | Global Locate, Inc. | Method and apparatus for mitigating multipath effects at a satellite signal receiver using a sequential estimation filter |
CN103064092B (zh) * | 2012-12-28 | 2014-05-21 | 中国科学院光电研究院 | 一种导航卫星的选择方法 |
CN103499822B (zh) * | 2013-10-14 | 2016-04-13 | 北京大学 | 一种基于最优gdop和牛顿恒等式的双星座快速选星方法 |
-
2014
- 2014-04-25 CN CN201410172299.9A patent/CN103926597B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN103926597A (zh) | 2014-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103926597B (zh) | 一种基于北斗rdss双向通信功能的多路径检测方法 | |
CN104714244B (zh) | 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法 | |
CN102273081B (zh) | 使用广义误差分布来进行位置估计的方法 | |
US20140306842A1 (en) | Accuracy and performance of the hybrid positioning system | |
CN102149192B (zh) | 一种基于移动台合作的蜂窝网络无线定位方法 | |
CN102638888B (zh) | 基于信号统计的室内定位指纹分组方法 | |
CN103616700B (zh) | 接收机和接收机评估所处环境的卫星信号遮挡状况的方法 | |
CN109672973B (zh) | 一种基于最强ap的室内定位融合方法 | |
CN104902562B (zh) | 一种基于多层指纹匹配的室内定位方法 | |
CN103168251B (zh) | 增强系统和增强方法 | |
CN103209474B (zh) | 一种移动终端定位的方法、定位服务器及服务基站 | |
CN106793082A (zh) | 一种在wlan/蓝牙异构网络环境中的移动设备定位方法 | |
EP3115807A2 (en) | Gnss receiver with an on-board capability to implement an optimal error correction mode | |
CN106054216A (zh) | 基于gdop和uere的多模gnss选星方法 | |
CN103428850A (zh) | 基于压缩感知的分布式多区域定位方法 | |
CN103618997B (zh) | 一种基于信号强度概率的室内定位方法和装置 | |
CN113031031B (zh) | 一种城市峡谷内基于gnss信号精确分类的加权定位方法 | |
TW201445168A (zh) | 一種接收機和衛星定位及測速方法 | |
CN106842250B (zh) | 一种基于遗传粒子滤波的多径抑制码环的设计方法 | |
CN102111873B (zh) | 选取可见基站的方法及装置、终端定位方法及装置 | |
CN106304330A (zh) | 一种减轻背景电磁波影响的射频层析成像定位方法 | |
CN105472733A (zh) | 一种应用于室内定位中的基于ap选择定位方法 | |
CN108196273A (zh) | 一种组合导航系统的定位卫星选择方法 | |
CN107153208A (zh) | 一种gps载波相位周跳探测与修复的方法 | |
CN111220065A (zh) | 一种多模式联合的北斗变形监测数据处理系统及方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 |
|
CF01 | Termination of patent right due to non-payment of annual fee |