CN109946722B - 一种多系统多频段定位方法及系统 - Google Patents
一种多系统多频段定位方法及系统 Download PDFInfo
- Publication number
- CN109946722B CN109946722B CN201910258587.9A CN201910258587A CN109946722B CN 109946722 B CN109946722 B CN 109946722B CN 201910258587 A CN201910258587 A CN 201910258587A CN 109946722 B CN109946722 B CN 109946722B
- Authority
- CN
- China
- Prior art keywords
- positioning
- satellite
- detection
- fault
- information
- 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
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种多系统多频段定位方法及系统,所述定位方法包括以下步骤:S1.接收卫星广播的星历观测信息,并进行预处理;S2.根据预处理后的卫星信息对卫星进行故障检测;S3.进行定位解算并将获得的定位信息输出至用户。本发明可根据信号监收情况自动切换定位模式,如果某个关闭系统或某个频段关闭,系统定位模式会自动切换;人为再卫星信号中添加干扰误差时,本系统可以自动识别出,并将添加了干扰的卫星信号剔除,有效完成定位。
Description
技术领域
本发明涉及GPS卫星定位领域,尤其涉及一种多系统多频段定位方法及系统。
背景技术
目前,市面上大部分电子产品只有GPS定位功能,部分产品包含GPS/北斗定位功能,但是对北斗做了权重处理,降低了北斗的权重,实质上还是以GPS定位为主。依赖GPS系统最大的问题就是风险不可控,美国曾多次区域性屏蔽GPS信号并加入随机误差,导致被屏蔽区域内依赖于GPS的接收机定位精度明显下降,甚至无法定位。
发明内容
为了解决上述问题,本发明提出一种多系统多频段定位方法及系统。
本发明采用以下技术方案:一种多系统多频段定位方法,所述定位方法包括以下步骤:S1.接收卫星广播的星历观测信息,并进行预处理;
S2.根据预处理后的卫星信息对卫星进行故障检测;
S3.进行定位解算并将获得的定位信息输出至用户。所述步骤S1包括以下检测子步骤:
S11.门限检测:对频段、星号、伪距、信噪比、多普勒和卫星仰角进行检测;
S12.通道稳定检测:对接收的卫星信息进行通道稳定性检测,当通道稳定一段时间后,即为稳定信息;
S13.卫星信息计算:根据卫星和接收机位置信息,计算出卫星的仰角fi、方位角Az、PDOP、TDOP、GDOP、VDOP和HDOP的值,并将7项参数计算结果作为步骤S2故障检测的判断依据,计算方式如下:
PDOP=sqrt(h11+h22+h33);
TDOP=sqrt(h44);
GDOP=sqrt(h11+h22+h33+h44);
VDOP=sqrt(S33);
HDOP=sqrt(S11+S22);
上式中,ΔX、ΔY和ΔZ为卫星与接收机距离,hii为地心地固直角坐标系权系数矩阵元素,Sii为站心坐标系权系数矩阵元素;
S14.定位模式检测:根据接收机的监收状态,匹配定位模式,所述定位模式为GPS定位模式、北斗定位模式或GPS/北斗组合定位模式。
所述步骤S2包括以下子步骤:
S21.判断加权后验证残差是否超出门限值,进而判断系统是否存在故障,所述加权后验证残差计算方法为:
y=R-Pr+Pion+Ptrot+Prot-Vc*Δtor-Vc*Δts+Vc*tr;
Sy=S*ΔL;
式中,w为权值系数,fi为卫星仰角,PDOP(i)为卫星空间几何结构,Pion为电离层,Ptrop为对流层,Factor为通道稳定时间,H为观测向量矩阵系数,R为站星几何距离、Pr为伪距、Prot为地球旋转修正参数、Vc为光速,Sy为残差;Δts为卫星钟差修正量,Δtor为相对论修正量,tr为接收机钟差,S为复合矩阵,ΔL为经过修正的伪距与站星几何距离之差;
S22.若加权后残差超出门限,则为故障卫星;
S23.将已找出的故障卫星进行标记,被标记的卫星不参与定位解算,重复步骤S21~S23,遍历找出所有故障卫星。
所述定位解算采用加权最小二乘融合算法进行定位解算,通过步骤S21得出的权值系数w、修正后的伪距与站星距离的残差y和经过筛选的完好卫星信息,计算出用户位置信息修正量ΔH,计算方式如下:
ΔH=(ATwA)-1ATwy;
X=X0+ΔH(1);
Y=Y0+ΔH(2);
Z=Z0+ΔH(3);
上式中,X0、Y0、Z0和X、Y、Z分别为接收机上一时刻和当前时刻位置,A为观测向量矩阵系数,ΔH(1)、ΔH(2)、ΔH(3)分别代表WGS 84坐标系三个轴向上的位移修正量。
所述步骤S3还包括后处理子步骤,实际使用过程中,接收机监收信号质量会受到多种干扰因素影响,例如多径效应影响,此时卫星并未发生故障,但是信号传输过程中受高楼、山体等影响,定位结果会出现异常,这就需要结合一段时间内定位数据,对当前时刻的解算数据进行容错处理。
一种多系统多频段定位系统,定位系统包括预处理模块、故障检测模块和定位模块;
所述预处理模块包括门限检测、通道稳定检测、计算卫星信息和定位模式检测子模块;
所述故障检模块包括故障卫星检测和剔除故障卫星子模块;
所述定位模块包括定位解算和后处理子模块。
本发明的有益效果在于:
(1)本发明可根据信号监收情况自动切换定位模式,如果某个关闭系统或某个频段关闭,系统定位模式会自动切换;
(2)人为再卫星信号中添加干扰误差时,本系统可以自动识别出,并将添加了干扰的卫星信号剔除,有效完成定位。
附图说明
图1是系统整体设计框架及流程图;
图2是加入干扰定位结果图;
图3是故障卫星检测时序图;
图4是定位解算结果图。
具体实施方式
为了使本领域的技术人员更好地理解本发明的技术方案,下面结合附图和具体实施例对本发明作进一步的详细说明。
实施例一:如图1所示,一种多系统多频段定位方法,所述定位方法包括以下步骤:
S1.接收卫星广播的星历观测信息,并进行预处理;
S2.根据预处理后的卫星信息对卫星进行故障检测;
S3.进行定位解算并将获得的定位信息输出至用户。所述步骤S1包括以下检测子步骤:
S11.门限检测:对频段、星号、伪距、信噪比、多普勒和卫星仰角进行检测,门限值如表1:
检测量 | 门限值 |
频段 | [1,6] |
星号 | [1,32] |
伪距 | [1e7,5e9] |
信噪比 | [40,100] |
多普勒 | [1e-3,5e4] |
卫星仰角 | [6,90] |
表1门限列表
S12.通道稳定检测:对接收的卫星信息进行通道稳定性检测,当通道稳定一段时间后,即为稳定信息;
S13.卫星信息计算:根据卫星和接收机位置信息,计算出卫星的仰角fi、方位角Az、PDOP、TDOP、GDOP、VDOP和HDOP的值,并将7项参数计算结果作为步骤S2故障检测的判断依据,计算方式如下:
PDOP=sqrt(h11+h22+h33);
TDOP=sqrt(h44);
GDOP=sqrt(h11+h22+h33+h44);
VDOP=sqrt(S33);
HDOP=sqrt(S11+S22);
上式中,ΔX、ΔY和ΔZ为卫星与接收机距离,hii为地心地固直角坐标系权系数矩阵元素,Sii为站心坐标系权系数矩阵元素;
S14.定位模式检测:根据接收机的监收状态,匹配定位模式,所述定位模式为GPS定位模式、北斗定位模式或GPS/北斗组合定位模式。
所述步骤S2包括以下子步骤:
S21.判断加权后验证残差是否超出门限值,进而判断系统是否存在故障,所述加权后验证残差计算方法为:
y=R-Pr+Pion+Ptrot+Prot-Vc*Δtor-Vc*Δts+Vc*tr;
Sy=S*ΔL;
式中,w为权值系数,fi为卫星仰角,PDOP(i)为卫星空间几何结构,Pion为电离层,Ptrop为对流层,Factor为通道稳定时间,H为观测向量矩阵系数,R为站星几何距离、Pr为伪距、Prot为地球旋转修正参数、Vc为光速,Sy为残差;Δts为卫星钟差修正量,Δtor为相对论修正量,tr为接收机钟差,S为复合矩阵,ΔL为经过修正的伪距与站星几何距离之差;
S22.若加权后残差超出门限,则为故障卫星;
S23.将已找出的故障卫星进行标记,被标记的卫星不参与定位解算,重复步骤S21~S23,遍历找出所有故障卫星。
所述定位解算采用加权最小二乘融合算法进行定位解算,通过步骤S21得出的权值系数w、修正后的伪距与站星距离的残差y和经过筛选的完好卫星信息,计算出用户位置信息修正量ΔH,计算方式如下:
ΔH=(ATwA)-1ATwy;
X=X0+ΔH(1);
Y=Y0+ΔH(2);
Z=Z0+ΔH(3);
上式中,X0、Y0、Z0和X、Y、Z分别为接收机上一时刻和当前时刻位置,A为观测向量矩阵系数,ΔH(1)、ΔH(2)、ΔH(3)分别代表WGS 84坐标系三个轴向上的位移修正量。
所述步骤S3还包括后处理步骤,实际使用过程中,接收机监收信号质量会受到多种干扰因素影响,例如多径效应影响,此时卫星并未发生故障,但是信号传输过程中受高楼、山体等影响,定位结果会出现异常,这就需要结合一段时间内定位数据,对当前时刻的解算数据进行容错处理。容错依据见表2:
容错量 | 门限值 |
速度 | [0,300] |
加速度 | [0,20] |
位移 | [0,354] |
位置 | [6377000,20200000] |
位移推算偏差 | [0,30] |
表2容错门限表
一种多系统多频段定位系统,定位系统包括预处理模块、故障检测模块和定位模块;
所述预处理模块包括门限检测、通道稳定检测、计算卫星信息和定位模式检测子模块;
所述故障检模块包括故障卫星检测和剔除故障卫星子模块;
所述定位模块包括定位解算和后处理子模块。
基于上述实施例,选取一组数据进行说明:本组数据时长600s,包含GPS和北斗卫星数据。
第一步:卫星数据加误差;
200s-250s,在一颗GPS卫星上加入10m误差干扰;
300s-350s,在两颗GPS卫星上加入10m误差干扰。
第二步:干扰影响分析;
利用加入了干扰的数据直接进行定位解算,定位结果如图2,在200s-250s和300s-350s两个时间段内,由于加入了卫星误差,相对定位误差跳变到5m左右;未加入卫星误差时间段,相对定位误差明显较小,仅为1m。再实际使用环境中,卫星信号传播过程中会引入各种误差,导致定位精度变差甚至无法定位。
第三步:故障卫星检测;
参与定位卫星个数见图3:在图中可以看出,未加入误差时间段,参与定位卫星个数为30颗;200s-250s参与定位卫星个数为29颗,300s-350s参与定位卫星个数为28颗,加入干扰的卫星已经被剔除。
第四步:定位解算;
剔除了加入干扰的卫星后,进行定位解算,定位结果见图4。
定位结果可以看出,经过卫星信号检测再进行定位后,定位误差恢复1m左右,说明本算法可以有效的解决GPS数据人为加入干扰的影响。
本发明的有益效果在于:
(1)本发明可根据信号监收情况自动切换定位模式,如果某个关闭系统或某个频段关闭,系统定位模式会自动切换;
(2)人为再卫星信号中添加干扰误差时,本系统可以自动识别出,并将添加了干扰的卫星信号剔除,有效完成定位。
需要说明的是,对于前述的各个方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本申请并不受所描述的动作顺序的限制,因为依据本申请,某一些步骤可以采用其他顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作和单元并不一定是本申请所必须的。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详细描述的部分,可以参见其他实施例的相关描述。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、ROM、RAM等。
以上所揭露的仅为本发明较佳实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明权利要求所作的等同变化,仍属本发明所涵盖的范围。
Claims (7)
1.一种多系统多频段定位方法,其特征在于:所述定位方法包括以下步骤:
S1.接收卫星广播的星历观测信息,并进行预处理;
S2.根据预处理后的卫星信息对卫星进行故障检测,包括以下子步骤:
S21.判断加权后验证残差是否超出门限值,进而判断系统是否存在故障,所述加权后验证残差计算方法为:
Sy=S*ΔL;
y=R-Pr+Pion+Prot+Prot-Vc*Δtor-Vc*Δts+Vc*tr;
式中,w为权值系数,fi为卫星仰角,PDOP(i)为卫星空间几何结构,Pion为电离层,Ptrop为对流层,Factor为通道稳定时间,H为观测向量矩阵系数,R为站星几何距离、Pr为伪距、Prot为地球旋转修正参数、Vc为光速,Sy为残差;Δts为卫星钟差修正量,Δtor为相对论修正量,tr为接收机钟差,S为复合矩阵,ΔL为经过修正的伪距与站星几何距离之差,y为修正后的伪距与站星距离的残差;
S22.若加权后残差超出门限,则为故障卫星;
S23.将已找出的故障卫星进行标记,被标记的卫星不参与定位解算,重复步骤S21~S23,遍历找出所有故障卫星;
S3.进行定位解算并将获得的定位信息输出至用户。
2.如权利要求1所述的一种多系统多频段定位方法,其特征在于:所述步骤S1包括以下检测子步骤:
S11.门限检测:对频段、星号、伪距、信噪比、多普勒和卫星仰角进行检测;
S12.通道稳定检测:对接收的卫星信息进行通道稳定性检测,当通道稳定一段时间后,即为稳定信息;
S13.卫星信息计算:根据卫星和接收机位置信息,计算出卫星的仰角fi、方位角Az、
PDOP、TDOP、GDOP、VDOP和HDOP参数值,并将参数值结果作为步骤S2故障检测的判断依据;
S14.定位模式检测:根据接收机的监收状态,匹配定位模式。
3.如权利要求1所述的一种多系统多频段定位方法,其特征在于:所述定位解算采用加权最小二乘融合算法进行定位解算,通过步骤S21得出的权值系数w、修正后的伪距与站星距离的残差y和经过筛选的完好卫星信息,计算出用户位置信息修正量ΔH。
4.如权利要求1所述的一种多系统多频段定位方法,其特征在于:所述步骤S3还包括后处理子步骤,考虑到接收机监收信号质量会受到多种干扰因素影响,此时卫星并未发生故障,但是信号传输过程中受高楼、山体的影响,定位结果会出现异常,需要结合一段时间内定位数据,对当前时刻的解算数据进行容错处理。
5.如权利要求2所述的一种多系统多频段定位方法,其特征在于:所述定位模式为GPS定位模式、北斗定位模式或GPS/北斗组合定位模式。
6.如权利要求3所述的一种多系统多频段定位方法,其特征在于:所述用户位置信息修正量ΔH的计算式为:
ΔH=(ATwA)-1ATwy;
X=X0+ΔH(1);
Y=Y0+ΔH(2);
Z=Z0+ΔH(3);
式中,X0、Y0、Z0和X、Y、Z分别为接收机上一时刻和当前时刻位置,A为观测向量矩阵系数,ΔH(1)、ΔH(2)和ΔH(3)分别代表坐标系三个轴向上的位移修正量。
7.一种多系统多频段定位系统,用于实现权利要求1~6任意一项所述的一种多系统多频段定位方法,其特征在于:定位系统包括预处理模块、故障检测模块和定位模块;
所述预处理模块包括门限检测、通道稳定检测、计算卫星信息和定位模式检测子模块;
所述故障检模块包括故障卫星检测和剔除故障卫星子模块;
所述定位模块包括定位解算和后处理子模块。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910258587.9A CN109946722B (zh) | 2019-04-01 | 2019-04-01 | 一种多系统多频段定位方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910258587.9A CN109946722B (zh) | 2019-04-01 | 2019-04-01 | 一种多系统多频段定位方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109946722A CN109946722A (zh) | 2019-06-28 |
CN109946722B true CN109946722B (zh) | 2023-04-11 |
Family
ID=67013323
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910258587.9A Active CN109946722B (zh) | 2019-04-01 | 2019-04-01 | 一种多系统多频段定位方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109946722B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111045056B (zh) * | 2019-11-29 | 2024-04-30 | 交通运输部长江通信管理局 | 一种接收机消除干扰卫星信号的方法及装置 |
CN114485655B (zh) * | 2022-04-06 | 2022-07-29 | 深圳华大北斗科技股份有限公司 | Gnss/ins组合导航数据质量控制方法 |
CN117949977B (zh) * | 2023-12-25 | 2024-07-26 | 北京航空航天大学 | 一种gnss转发式无人机诱骗的选星和时延控制方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106707304A (zh) * | 2016-12-23 | 2017-05-24 | 天津七六四通信导航技术有限公司 | 一种卫星导航接收机故障星检测方法 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6033746A (ja) * | 1983-08-04 | 1985-02-21 | Nec Corp | 衛星回線監視方式 |
JPH0815404A (ja) * | 1994-06-30 | 1996-01-19 | Nec Home Electron Ltd | Gps受信システム |
CN100442077C (zh) * | 2007-07-12 | 2008-12-10 | 北京航空航天大学 | 接收机自主完好性监测的可用性预测方法 |
CN104181562B (zh) * | 2013-05-25 | 2016-07-13 | 成都国星通信有限公司 | 一种全球导航卫星系统的卫星优选和定位方法 |
CN103592658A (zh) * | 2013-09-30 | 2014-02-19 | 北京大学 | 多模卫星导航系统中基于选星算法的raim新方法 |
CN103592657B (zh) * | 2013-09-30 | 2016-08-31 | 北京大学 | 一种基于钟差辅助的低可见星下单模raim实现方法 |
US9678212B2 (en) * | 2014-02-28 | 2017-06-13 | Honeywell International Inc. | Satellite measurement screening to protect the integrity of existing monitors in the presence of amplitude scintillation |
CN103954982A (zh) * | 2014-04-18 | 2014-07-30 | 中国人民解放军国防科学技术大学 | 基于多模gnss接收机的可见卫星快速选择方法 |
CN105487088B (zh) * | 2015-09-12 | 2017-12-01 | 北京大学 | 一种卫星导航系统中基于卡尔曼滤波的raim算法 |
CN105549033B (zh) * | 2015-12-07 | 2018-02-23 | 郑州威科姆科技股份有限公司 | 一种基于最小二乘残差边缘检验的完好性处理方法 |
CN109085619B (zh) * | 2017-06-14 | 2020-09-25 | 展讯通信(上海)有限公司 | 多模gnss系统的定位方法及装置、存储介质、接收机 |
CN107783157B (zh) * | 2017-09-20 | 2020-04-10 | 西安空间无线电技术研究所 | 基于导航卫星的外辐射源智能选星方法、系统及雷达平台 |
-
2019
- 2019-04-01 CN CN201910258587.9A patent/CN109946722B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106707304A (zh) * | 2016-12-23 | 2017-05-24 | 天津七六四通信导航技术有限公司 | 一种卫星导航接收机故障星检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109946722A (zh) | 2019-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3627188A1 (en) | Correction information integrity monitoring in navigation satellite system positioning methods, systems, and devices | |
CN111045034B (zh) | 基于广播星历的gnss多系统实时精密时间传递方法及系统 | |
CN109946722B (zh) | 一种多系统多频段定位方法及系统 | |
CN109709579B (zh) | 一种基于用户测距误差实时估计的gnss卫星星历故障检测方法 | |
US5917445A (en) | GPS multipath detection method and system | |
AU2019204184B2 (en) | A method for operating a plurality of GNSS receivers for detecting satellite signal deformation | |
EP3012665B1 (en) | Systems and methods for averaging satellite sigmas and readmitting excluded satellite measurements into differential corrections and integrity monitors | |
CN101419275B (zh) | 基于多接收机的局域机场监视方法和系统 | |
EP3109672A1 (en) | Gnss receiver with a capability to resolve ambiguities using an uncombined formulation | |
CN102004259A (zh) | 高灵敏度环境下基于多普勒平滑伪距的卫星导航解算定位方法 | |
CN111913201B (zh) | 一种gnss差分定位方法、装置及计算机可读存储介质 | |
CN107505636A (zh) | 海基jpals的定位域mrcc方法及装置 | |
CN108363077A (zh) | 一种精密单点定位装置中的载波相位周跳修复方法 | |
CN115826016B (zh) | 一种北斗双频星基增强改正数及完好性参数解算的方法 | |
CN112130177A (zh) | 一种基于稳定分布的地基增强系统完好性监测方法 | |
CN104280741A (zh) | 电离层异常检测方法 | |
CN113126129A (zh) | 基于空间信号质量加权估计的gbas完好性监测方法 | |
CN115390096A (zh) | 基于全视角星载gnss接收系统的低轨卫星实时相对定轨方法 | |
CN114839656A (zh) | 一种部分模糊度固定的时间传递快速重收敛方法 | |
CN111142125B (zh) | 一种卫星完好性监测方法和系统 | |
CN114690220A (zh) | 一种监测数据实时解算星历异常处理方法 | |
CN111913195B (zh) | 基于陆基无线电导航信息的gps接收机抗欺骗干扰处理方法 | |
CN114325779B (zh) | 导航增强系统定位粗差检测方法及装置 | |
CN110941002B (zh) | 一种自适应抗差的序贯最小二乘精密单点定位方法 | |
US20120274511A1 (en) | Satellite Based Augmentation 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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |