CN101609140A - 一种兼容导航接收机定位系统及其定位方法 - Google Patents
一种兼容导航接收机定位系统及其定位方法 Download PDFInfo
- Publication number
- CN101609140A CN101609140A CNA200910089221XA CN200910089221A CN101609140A CN 101609140 A CN101609140 A CN 101609140A CN A200910089221X A CNA200910089221X A CN A200910089221XA CN 200910089221 A CN200910089221 A CN 200910089221A CN 101609140 A CN101609140 A CN 101609140A
- Authority
- CN
- China
- Prior art keywords
- mtd
- mrow
- mtr
- msub
- msubsup
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000004807 localization Effects 0.000 title abstract 2
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 239000000284 extract Substances 0.000 claims abstract description 11
- 238000012937 correction Methods 0.000 claims description 77
- 238000001914 filtration Methods 0.000 claims description 32
- 238000006243 chemical reaction Methods 0.000 claims description 24
- 239000005436 troposphere Substances 0.000 claims description 24
- 239000005433 ionosphere Substances 0.000 claims description 23
- 230000001133 acceleration Effects 0.000 claims description 15
- 238000005259 measurement Methods 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000013179 statistical model Methods 0.000 claims description 5
- 230000003044 adaptive effect Effects 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical compound C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 5
- 230000009977 dual effect Effects 0.000 description 3
- 238000013213 extrapolation Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种兼容导航接收机定位系统及其定位方法。本系统包括参数配置模块、数据预处理器、卫星参数解算模块、伪距修正模块、坐标转换模块、定位解算模块、卡尔曼滤波器和报文形成模块。首先读取接收机的配置参数,然后提取星历数据、历书数据和观测数据,接着获取有效卫星的参数信息,再获取接收机位置和接收机速度,最后将导航信息生成报文。本发明实现了卫星导航系统的选择使用,北斗二代导航系统的独立导航摆脱了非常时期对其他导航系统的依赖,多卫星导航系统的组合导航增加了有效卫星的数量,保证了卫星定位服务的连续性、可用性和安全性,增强了观测卫星的几何图形强度,提高了整个卫星定位系统的可靠性。
Description
技术领域
本发明属于卫星导航领域,具体涉及一种兼容导航接收机定位系统及其定位方法。
背景技术
卫星导航能够向各类用户和运动平台实时提供准确、连续的位置、速度和时间信息。全球卫星导航定位系统(Global Navigation Satellite System,简称为GNSS)是第二代卫星导航定位系统,具有陆海空全能性、全球性、全天候、连续性和实时性的特点。
目前世界上已经存在的两大全球卫星导航定位系统是美国的全球定位系统(GlobalPositioning System,简称为GPS)系统和俄罗斯的全球导航卫星系统(Global NavigationSatellite System,简称为GLONASS),正在设计建设阶段的有欧洲的“伽利略”卫星定位系统(Galileo GNSS,简称为Galileo)和中国的北斗二代导航系统(BeiDou-2,简称为BD2)。
GPS是目前发展最为成熟的导航系统,但始终受美军方控制,在发生危机时,随时可能被关闭;BD2目前为区域导航系统,存在12颗定位卫星,单BD2在高遮蔽角的情况下定位精度较差,难以满足高精度定位要求,但其作为我国具有独立知识产权的全球卫星导航定位系统,它的应用在非常时期是必要的。
各系统的卫星星座规模有限,独立应用时还存在一定问题。首先,定位精度低,并且位置精度随用户接收机所处的位置和时间产生较大变化。其次,可靠性差,特别是在干扰较强的环境中,可能出现短暂信号中断或捕获卫星数少的情况,定位精度将受到很大影响。组合导航的应用得到广泛的研究,但迄今为止,尚未有包含BD2的多卫星导航系统兼容导航接收机定位系统,尚未提出包含BD2的多卫星导航系统兼容导航接收机定位方法。
GPS、GLONASS、Galileo、BD2系统定义的参考坐标系基本是一致的,但由于不可避免的存在测轨跟踪站站址坐标误差和测量误差,定义的坐标系与实际使用的坐标系存在一定的差异。因此,实际上各系统所采用的坐标系两两之间都存在差异。因此,在不同卫星导航系统组合应用时,需要对坐标系进行转换。确定两个坐标系转换参数的数学模型有很多。具有代表性的模型为布尔莎模型(Bursa Model,简称Bursa模型)和莫洛金斯基模型(Molodensky Model,简称Molodensky模型),Bursa模型较为常用。
7参数Bursa模型:
其中,ΔX为地球中心偏移量在B系统参考坐标系X轴方向上的分量,ΔY为地球中心偏移量在B系统参考坐标系Y轴方向上的分量,ΔZ为地球中心偏移量在B系统参考坐标系Z轴方向上的分量,θx为B系统参考坐标系X轴到A系统参考坐标系X轴的旋转角,θy为B系统参考坐标系Y轴到A系统参考坐标系Y轴的旋转角,θz为B系统参考坐标系Z轴到A系统参考坐标系Z轴的旋转角, 为比例因子,la和lb分别是同一条基线在A系统和B系统两坐标系中的基线长度。
GPS、GLONASS、Galileo、BD2系统皆为测时-测距导航系统,以测码伪距观测方式为例,用户接收机通过观测由卫星发射的测距码到用户接收机的传播时间,乘以光速c后获得接收机与卫星间通视距离的观测量。基于这一测量原理,四系统对时间测量以及对时间同步的要求都比较高,使得不同导航系统为了确保各自系统运行的独立性、可靠性以及稳定性,都独立维持其系统时间参考系统的运行。即使采用相同精度的设备,采用同样的技术也很难实现严格同步,导致多星座组合导航应用中需要进行时间系统统一,通常的方法是增加观测量,将时间同步参数作为未知数求解,便可以实现不同系统间的时间系统统一。
目前,普通接收机通用的单机定位方法是基于线性化的迭代求解算法,因其在大部分情况下都能取得良好的定位效果,并且算法简单,计算量相对较小,便于计算机实现。
基于线性化的迭代求解算法基于接收机与卫星Sj间的伪距公式构造导航定位方法:
其中,j为观测卫星Sj的编号,为t时刻卫星Sj的观测伪距,Rj(t)为用户接收机与卫星Sj间的几何距离,δ(t)为t时刻用户接收机时钟与卫星导航系统的系统时之间的钟差,δtj(t)为t时刻卫星Sj的星钟与卫星导航系统的系统时之间的钟差,cδtj(t)称为星钟修正误差,ΔI j(t)为t时刻电离层折射对测码伪距的影响,称为电离层修正误差,ΔT j(t)为t时刻大气对流层折射对测码伪距的影响,称为对流层修正误差,(xj(t),yj(t),zj(t))T为卫星Sj的位置坐标,(x(t),y(t),z(t))T为用户接收机的位置坐标,c为光速;其中x(t)、y(t)、z(t)、δ(t)为未知量。
设进行星钟误差、电离层误差、对流层误差修正后的观测伪距为ρj(t),则有:
即:
由于上述方程是非线性的,需运用泰勒一阶展开来线性化。
设用户接收机的近似位置的位置坐标为真实位置和近似位置之间的偏离为(δx(t),δy(t),δz(t))T,将Rj(t)按泰勒级数绕近似位置展开,将位置偏离(δx(t),δy(t),δz(t))T表示为已知卫星坐标和伪距测量值的线性函数:
即:
其中, 为接收机近似位置到卫星Sj的几何距离;
当有效卫星数目满足系统定位要求时,可以联立线性方程生成矩阵,然后用最小二乘法迭代求解。
但是,用一阶泰勒展开不能完全表征信号的真实特性,在理论上导致明显的信号误差,以至于导致使用该方法进行定位时存在定位精度不高、恶劣定位条件下鲁棒性不强的不足。因此,为满足对定位精度要求较高的用户,考虑对线性迭代求解方法的结果进行滤波。
滤波的目的是估计当前载体的运动状态,对于导航接收机来说就是估计天线相位中心的运动状态,包括其在导航解算坐标系内的位置、速度等。目前,最基本的数据滤波方法包括线性回归滤波,两点外推滤波、为纳滤波、加权最小二乘滤波、α-β与α-β-γ滤波、卡尔曼滤波和简化的卡尔曼滤波等。从跟踪精度来看,以卡尔曼滤波为最好,两点外推滤波为最差;从计算量来看,以卡尔曼滤波为最大,两点外推滤波最小。虽然,随着现代微处理技术突飞猛进的发展,计算要求已不再成为应用卡尔曼滤波的主要障碍,但在多星座组合定位应用中,随着观测向量或状态向量维数的扩充,卡尔曼滤波的计算量将会有大量增长。因此,卡尔曼滤波的使用,需根据用户对跟踪精度和计算量的要求进行折衷考虑。
常规卡尔曼滤波需要准确地知道系统噪声和观测噪声的统计特性,而实际应用中观测噪声的随机特性取决于测量设备以及系统的工作环境,很难准确描述,通常采用自适应卡尔曼滤波解决这一问题。为合理的应用卡尔曼滤波技术进行动态定位,最重要的是要建立载体的运动模型。常速度模型和常加速度模型构造简单易于实现,但适用匀速或匀加速等无明显机动的场合;一阶时间相关模型对一般机动的载体具有很好地描述效果,但如突发剧烈机动效果迅速恶化;机动目标“当前”统计模型本质上是非零均值时间相关模型,其机动加速度的“当前”概率密度用修正的瑞利分布描述,均值为“当前”加速度预测值,相比之下能更为真实地反映目标机动范围和强度的变化。
常速度模型:
常加速度模型:
一阶时间相关模型:
机动目标“当前”统计模型:
发明内容
本发明的目的是提供一种包含BD2的多卫星导航系统兼容导航接收机定位系统及其定位方法,可以选择用于进行定位的卫星导航系统,单BD2定位可以在非常时期摆脱对其他卫星导航系统的依赖,多卫星导航系统的组合定位可以增加有效卫星的数量、保证卫星定位服务的连续性和可用性、增强观测卫星的几何图形强度、提高整个卫星定位系统的可靠性。
本发明提供的一种兼容导航接收机定位系统,包括参数配置模块、数据预处理器、卫星参数解算模块、伪距修正模块、坐标转换模块、定位解算模块、卡尔曼滤波器和报文形成模块。其中,参数配置模块的第一个输出端与数据预处理器的第一个输入端相连,参数配置模块的第二个输出端与卫星参数解算模块的第一个输入端相连,参数配置模块的第三个输出端与伪距修正模块的第一个输入端相连,参数配置模块的第四个输出端与坐标转换模块的第一个输入端相连,参数配置模块的第五个输出端与定位解算模块的第一个输入端相连,参数配置模块的第六个输出端与卡尔曼滤波器的第一个输入端相连,参数配置模块的第七个输出端与报文形成模块的第一个输入端相连;导航电文从数据预处理器的第二个输入端输入,观测数据从数据预处理器的第三个输入端输入,数据预处理器的第一个输出端与卫星参数解算模块的第二个输入端相连,数据预处理器的第二个输出端与伪距修正模块的第二个输入端相连;卫星参数解算模块的输出端与坐标转换模块的第二个输入端相连,坐标转换模块的输出端与定位解算模块的第二个输入端相连;伪距修正模块的输出端与定位解算模块的第三个输入端相连;定位解算模块的输出端与卡尔曼滤波器的第二个输入端相连;卡尔曼滤波器的第一个输出端与伪距修正模块的第三个输入端相连,卡尔曼滤波器的第二个输出端与报文形成模块的第二个输入端相连。其中,参数配置模块读取兼容导航接收机的配置参数,并将配置参数输出给定位系统的各模块以控制其工作状态;数据预处理器采集有效卫星的导航电文和观测数据帧,并从中提取星历数据、历书数据和观测信息;卫星参数解算模块获取有效卫星的卫星位置和卫星速度;伪距修正模块获取有效卫星的观测伪距和伪距误差修正信息,并对观测伪距进行修正;坐标转换模块将来自不同卫星导航系统的有效卫星的坐标统一到同一参考坐标系中;定位解算模块获取接收机位置和接收机速度信息;卡尔曼滤波器对接收机位置和接收机速度信息进行滤波;报文形成模块形成导航报文。
本发明提出的一种兼容导航接收机定位方法,具体包含以下步骤:
步骤一:读取接收机的配置参数;
参数配置模块读取接收机的配置参数,
首先,根据用户需求,确定接收机的工作状态和使用的卫星导航系统;工作状态包括:单卫星导航系统定位、双卫星导航系统组合定位、三卫星导航系统组合定位和四卫星导航系统组合定位,使用的卫星导航系统包括:全球定位系统GPS、全球导航卫星系统GLONASS、北斗二代导航系统BD2和“伽利略”卫星定位系统Galileo;
然后,根据用户需求,确定卡尔曼滤波器的工作状态,确定导航信息的数据更新频率,确定导航报文的格式;
最后,将接收机的配置参数输出给数据预处理器、卫星参数解算模块、伪距修正模块、坐标转换模块、定位解算模块、卡尔曼滤波器和报文形成模块,以控制各模块的工作状态;
其中,接收机的配置参数包括:接收机的工作状态、使用的卫星导航系统、卡尔曼滤波器的工作状态、导航信息的数据更新频率以及导航报文的格式;
步骤二:提取星历数据、历书数据和观测数据;
根据参数配置模块输出的接收机的工作状态和使用的卫星导航系统,数据预处理器依据有效卫星所属的卫星导航系统的说明书对采集的有效卫星的导航电文进行拼接、校验、纠错,从中提取有效卫星的星历数据和历书数据;
根据参数解算模块输出的接收机工作状态和使用的卫星导航系统,数据预处理器依据捕获跟踪系统提供的观测数据帧格式从采集的有效卫星观测数据帧中提取毫秒数、码片偏移、相位偏移和多普勒频率;
步骤三:获取有效卫星的参数信息;
首先,根据参数配置模块输出的接收机的工作状态和使用的卫星导航系统,卫星参数解算模块依据有效卫星所属的卫星导航系统的说明书利用有效卫星的星历数据获取其位置坐标和速度坐标;
然后,根据参数配置模块输出的接收机的工作状态和使用的卫星导航系统,伪距修正模块依据有效卫星所属卫星导航系统的说明书利用有效卫星的星历数据和历书数据获取其星钟修正误差、电离层修正误差和对流层修正误差;
之后,伪距修正模块通过将有效卫星的周内秒计数、比特数、毫秒数、码片偏移、相位偏移进行累加获取其观测伪距;
最后,伪距修正模块利用获取的有效卫星的星钟修正误差、电离层修正误差和对流层修正误差对该卫星的观测伪距进行误差修正,获取修正后的观测伪距;
步骤四:获取接收机位置和接收机速度;
首先,根据参数配置模块输出的接收机的工作状态和使用的卫星导航系统,在单卫星导航系统定位时坐标转换模块不参与系统工作,在多卫星导航系统组合定位时坐标转换模块依据各卫星导航系统参考坐标系的关系利用7参数Bursa模型将来自不同卫星导航系统的有效卫星的坐标转换到同一参考坐标系,获取各有效卫星在同一参考坐标系中的三维坐标;
然后,根据参数配置模块输出的接收机的工作状态和使用的卫星导航系统,定位解算模块利用有效卫星的参数信息获取接收机位置和接收机速度;
利用有效卫星的位置坐标和修正后的观测伪距通过线性迭代方法获取接收机的位置,利用有效卫星的速度坐标和多普勒频率通过线性迭代方法获取接收机的速度;
最后,根据参数配置模块输出的接收机的工作状态、使用的卫星导航系统和卡尔曼滤波器的工作状态,确定卡尔曼滤波器是否对定位解算模块输出的接收机位置和接收机速度信息进行滤波处理;
如果卡尔曼滤波器的工作状态是工作,基于机动目标“当前”统计模型建立系统模型和测量模型,并通过自适应卡尔曼滤波获取接收机位置和接收机速度的滤波结果;
步骤五:将导航信息生成报文;
根据参数配置模块输出的用户对导航报文的要求,报文形成模块将导航信息以用户需要的格式形成报文。
本发明的优点在于:
(1)本发明采用了多卫星导航系统兼容定位技术,单BD2独立定位可以在非常时期摆脱对其他导航系统的依赖;
(2)本发明采用了多卫星导航系统兼容定位技术,多卫星导航系统组合定位可以大大增加有效卫星的数量,提高整个卫星定位系统的完善性监测性能,保证卫星定位服务的连续性、可用性,提高卫星导航定位应用的安全性。
(3)本发明采用了多卫星导航系统兼容定位技术,多卫星导航系统组合定位可以增强观测卫星的几何图形强度,提高整个卫星定位系统的可靠性。
附图说明
图1是本发明的系统框图;
图2是本发明的方法流程图;
图3是本发明组合定位与单GPS定位仰角限为10度时的有效卫星数目曲线图;
图4是本发明组合定位与单GPS定位仰角限为30度时的有效卫星数目曲线图;
图5是本发明组合定位与单GPS定位仰角限为10度时的DOP曲线图;
图6是本发明组合定位与单GPS定位仰角限为30度时的DOP曲线图;
图7是本发明组合定位与单GPS定位时的水平定位误差曲线图;
图8是本发明组合定位与单GPS定位时的垂直定位误差曲线图。
图中:1.参数配置模块 2.数据预处理器 3.卫星参数解算模块 4.伪距修正模块5.坐标转换模块 6.定位解算模块 7.卡尔曼滤波器 8.报文形成模块
具体实施方式
下面将以选择BD2和GPS进行双卫星导航系统组合定位为例,结合附图和实施例对本发明作进一步的详细说明。
本发明是一种兼容导航接收机定位系统及其方法,系统框图如图1所示,它包括参数配置模块1、数据预处理器2、卫星参数解算模块3、伪距修正模块4、坐标转换模块5、定位解算模块6、卡尔曼滤波器7和报文形成模块8。其中,参数配置模块1的第一个输出端与数据预处理器2的第一个输入端相连,参数配置模块1的第二个输出端与卫星参数解算模块3的第一个输入端相连,参数配置模块1的第三个输出端与伪距修正模块4的第一个输入端相连,参数配置模块1的第四个输出端与坐标转换模块5的第一个输入端相连,参数配置模块1的第五个输出端与定位解算模块6的第一个输入端相连,参数配置模块1的第六个输出端与卡尔曼滤波器7的第一个输入端相连,参数配置模块1的第七个输出端与报文形成模块8的第一个输入端相连;导航电文从数据预处理器2的第二个输入端输入,观测数据从数据预处理器2的第三个输入端输入,数据预处理器2的第一个输出端与卫星参数解算模块3的第二个输入端相连,数据预处理器2的第二个输出端与伪距修正模块4的第二个输入端相连;卫星参数解算模块3的输出端与坐标转换模块5的第二个输入端相连,坐标转换模块5的输出端与定位解算模块6的第二个输入端相连;伪距修正模块4的输出端与定位解算模块6的第三个输入端相连;定位解算模块6的输出端与卡尔曼滤波器7的第二个输入端相连;卡尔曼滤波器7的第一个输出端与伪距修正模块4的第三个输入端相连,卡尔曼滤波器7的第二个输出端与报文形成模块8的第二个输入端相连;
其中,参数配置模块1读取兼容导航接收机的配置参数,并将配置参数输出给定位系统的各模块以控制其工作状态;
数据预处理器2采集有效卫星的导航电文和观测数据帧,并从中提取星历数据、历书数据和观测信息;
卫星参数解算模块3获取有效卫星的卫星位置和卫星速度;
伪距修正模块4获取有效卫星的观测伪距和伪距误差修正信息,并对观测伪距进行修正;
坐标转换模块5将来自不同卫星导航系统的有效卫星的坐标统一到同一参考坐标系中;
定位解算模块6获取接收机位置和接收机速度信息;
卡尔曼滤波7对接收机的位置和速度信息进行滤波;
报文形成模块8形成导航报文。
本发明所述兼容导航接收机定位方法流程图如图2所示。
接收机完成了捕获跟踪,获得了一个完整的导航电文子帧并且观测数据已经更新后,需要根据用户的需要输出导航信息。本发明所述兼容导航接收机定位方法的具体实施方式如下:
步骤一:读取接收机的配置参数;
参数配置模块1读取接收机的配置参数,
首先,根据用户需求,确定接收机的工作状态和使用的卫星导航系统;工作状态为:双卫星导航系统组合定位,使用的卫星导航系统包括:北斗二代导航系统BD2和全球定位系统GPS;
然后,根据用户需求,确定卡尔曼滤波器7的工作状态为工作,确定导航信息的数据更新频率为10Hz,确定导航报文的格式为RINEX格式;
最后,将接收机的配置参数输出给数据预处理器2、卫星参数解算模块3、伪距修正模块4、坐标转换模块5、定位解算模块6、卡尔曼滤波器7和报文形成模块8,以控制各模块的工作状态;
其中,接收机的配置参数包括:接收机的工作状态、使用的卫星导航系统、卡尔曼滤波器7的工作状态、导航信息的数据更新频率以及导航报文的格式;
步骤二:提取星历数据、历书数据和观测数据;
根据参数配置模块1输出的接收机的工作状态和使用的卫星导航系统,接收机使用BD2和GPS进行双卫星导航系统组合定位,数据预处理器1依据BD2说明书,对采集的BD2有效卫星的导航电文子帧进行电文拼接和纠错,并从纠错后的电文中提取星历数据和历书数据;数据预处理器1依据GPS说明书,对采集的GPS有效卫星的导航电文子帧进行电文拼接和校验,并从校验后的电文中提取星历数据和历书数据;
其中,星历数据包括轨道长半轴方根轨道偏心率e,参考时刻轨道倾角i0,近地点角距ω,参考时刻升交点赤经Ω0,参考时刻的平近点角M0,平均运行速度差Δn,轨道倾角变率升交点赤经变率升交距角的调和改正项中的正弦项振幅Cus,升交距角的调和改正项中的余弦项振幅Cuc,卫星地心距的调和改正项中的正弦项振幅Crs,卫星地心距的调和改正项中的余弦项振幅Crc,轨道倾角的调和改正项中的正弦项振幅Cis,轨道倾角的调和改正项中的余弦项振幅Cic,星历参数的参考历元toe,星历数据的龄期AODE;
其总,历书数据包括时钟数据基准时间toc,卫星时间修正参数a0,卫星时间修正参数a1、卫星时间修正参数a2,卫星钟数据有效龄期AODC,延迟改正参数(GPS单频接收机延迟改正参数TGD;BD2星上设备时延差TGD1,BD2星上设备时延差TGD2,TGD1为BD2卫星发射的B1频点与B3频点信号之间存在的设备时延之差,TGD2为BD2卫星发射的B2频点与B3频点信号之间存在的设备时延之差;电离层修正参数GPS:αn、βn,其中n=1,2,3;BD2D1码:αn、βn,其中n=1,2,3,D2码:A1、B、αn、βn、γn,其中n=1,2,3;
根据参数配置模块1输出的接收机的工作状态和使用的卫星导航系统,接收机使用BD2和GPS进行双卫星导航系统组合定位,数据预处理器2依据捕获跟踪系统提供的观测数据帧格式,从采集的观测数据帧中提取BD2和GPS双系统有效卫星的毫秒数、码片偏移和相位偏移以及观测信息多普勒频率;
步骤三:获取有效卫星的参数信息;
首先,根据参数配置模块1输出的接收机的工作状态和使用的卫星导航系统,接收机使用BD2和GPS进行双卫星导航系统组合定位,卫星参数解算模块3分别依据BD2说明书和GPS说明书,利用双系统有效卫星的星历数据获取其位置坐标和速度坐标;
然后,根据参数配置模块1输出的接收机的工作状态和使用的卫星导航系统,接收机使用BD2和GPS进行双卫星导航系统组合定位,伪距修正模块4分别依据BD2说明书和GPS说明书,利用双系统有效卫星的星历数据和历书数据获取其星钟修正误差、电离层修正误差和对流层修正误差;
之后,伪距修正模块4通过将双系统有效卫星的周内秒计数、比特数、毫秒数、码片偏移、相位偏移进行累加获取其观测伪距;
最后,伪距修正模块4利用获取的双系统有效卫星的星钟修正误差、电离层修正误差和对流层修正误差对其观测伪距进行误差修正,获取修正后的观测伪距;
步骤四:获取接收机位置和接收机速度;
首先,根据参数配置模块1输出的接收机的工作状态和使用的卫星导航系统,接收机使用BD2和GPS进行双卫星导航系统组合定位,BD2参考坐标系CGS2000与GPS参考坐标系WGS-84间的差异很小,在10cm以内,认为是等同的。故使用GPS与BD2双系统进行定位时,坐标转换模块5不需要参与工作,直接将卫星参数解算模块3输出的有效卫星坐标输入给定位解算模块6即可;
然后,根据参数配置模块1输出的接收机的工作状态和使用的卫星导航系统,接收机使用BD2和GPS进行双卫星导航系统组合定位,定位解算模块6利用BD2和GPS双系统的有效卫星的参数信息获取接收机位置和接收机速度;
A、获取接收机位置
大气折射对伪距观测量的影响已改正,且卫星钟差可由导航电文给出的参数加以修正,则站星伪距观测为:
BD2和GPS双系统工作模式下,
其中,j=1,2,…,n,j为观测卫星Sj的编号,n为观测卫星总数;ρu j为卫星Sj的观测伪距进行星历误差、电离层误差和对流层误差修正后的伪距,Ru j为接收机与卫星Sj间的几何距离,Δt1为接收机时钟与BD2系统时钟之间的钟差,Δt2为接收机时钟与GPS系统时钟之间的钟差,c为光速。
在一段时间内,接收机连续同步观测不少于5颗卫星时,利用有效卫星的位置坐标和修正后的观测伪距获取接收机的位置坐标偏移和钟差模型:
其中,上述位置坐标偏移和钟差模型等号右侧部分的
为一个n×5的矩阵模块, 为一个n×1的矩阵模块,
其中,j为观测卫星Sj的编号,j=1,2,…,n,n为观测卫星总数,S1~Sq为BD2的卫星,Sq+1~Sn为GPS的卫星;
(Δxu,Δyu,Δzu)T为用户接收机的位置坐标偏移;
Δt1为用户接收机时钟与BD2系统时钟之间的钟差,Δt2为用户接收机时钟与GPS系统时钟之间的钟差;
(xj,yj,zj)T为卫星Sj的位置坐标,j=1,…q,q+1,…,n;
为用户接收机的近似位置与卫星Sj间的几何距离,j=1,…q,q+1,…,n,为用户接收机的近似位置与BD2的卫星S1间的几何距离,为用户接收机的近似位置与BD2的卫星Sq间的几何距离,为用户接收机的近似位置与GPS的卫星Sq+1间的几何距离,为用户接收机的近似位置与GPS的卫星Sn间的几何距离;
ρu j为对卫星Sj的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距,j=1,…q,q+1,…,n,ρu 1为对BD2的卫星S1的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距,ρu q为对BD2的卫星Sq的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距,ρu q+1为对GPS的卫星Sq+1的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距,ρu n为对GPS的卫星Sn的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距;
c为光速;
用户接收机的位置坐标为:
B、获取接收机速度
由于用户接收机载体和卫星之间的相对运动,接收机接收到的载波信号与卫星发射的载波信号频率不同,产生多普勒频移。频移的大小与接收机与卫星之间距离的变率有关。改正大气折射对伪距观测量的影响,并由导航电文给出的参数修正卫星钟差,则可得站星观测伪距率为:
BD2和GPS双系统工作模式下,
其中,j为观测卫星Sj的编号,j=1,2,…,n,n为观测卫星总数,du j为卫星Sj的观测伪距率,ax j,ay j,az j各项分别表示从用户接收机位置指向卫星Sj的单位向量的方向余弦,(vx j,vy j,vz j)T为卫星Sj的速度坐标,(vux,vuy,vuz)T为用户接收机的速度坐标,为接收机时钟与BD2系统时钟之间的钟漂,为接收机时钟与GPS系统时钟之间的钟漂,c为光速。
在一段时间内,接收机连续同步观测不少于5颗卫星时,利用接收机的位置坐标和有效卫星的多普勒频率获取接收机的速度坐标和钟漂模型:
其中,上述速度坐标和钟漂模型等号右侧部分的 为一个n×5的矩阵模块,d=[d1…dq dq+1…dn]T为一个n×1的矩阵模块,
其中
j为观测卫星Sj的编号,j=1,…q,q+1,…,n,n为观测卫星总数,S1~Sq为BD2的卫星,Sq+1~Sn为GPS的卫星;
(vux,vuy,vuz)T为用户接收机的速度坐标;
(xu,yu,zu)T为用户接收机的位置坐标;
(xj,yj,zj)T为卫星Sj的位置坐标,j=1,…q,q+1,…,n;
Ru j为用户接收机与卫星Sj间的几何距离,j=1,…q,q+1,…,n;
du j为卫星Sj的观测伪距率,j=1,…q,q+1,…,n,du 1为BD2的卫星S1的观测伪距率,du q为BD2的卫星Sq的观测伪距率,j=q+1时du q+1为GPS的卫星Sq+1的观测伪距率,j=n时du n为GPS的卫星Sn的观测伪距率;
(vx j,vy j,vz j)T为卫星Sj的速度坐标,j=1,…q,q+1,…,n;
c为光速;
(3)根据参数配置模块1输出的接收机工作状态和卡尔曼滤波器7的工作状态,卡尔曼滤波器7基于机动目标“当前”统计模型建立系统模型和测量模型,并通过自适应卡尔曼滤波获取接收机位置和接收机速度的滤波结果;
采用递推方法给出滤波结果:
估计均方误差模型P(k)=[I-K(k)H(k)]P(k|k-1);
其中,状态一步预测模型
一步预测均方误差模型P(k|k-1)=Φ(k|k-1)P(k-1)ΦT(k|k-1)+Q(k-1),
滤波增益模型K(k)=P(k|k-1)HT(k)[H(k)P(k|k-1)HT(k)+R(k)]-1;
其中,系统状态模型 分别为接收机在参考坐标系x轴方向的位置、速度和加速度,分别为接收机在参考坐标系y轴方向的位置、速度和加速度,分别为接收机在参考坐标系z轴方向的位置、速度和加速度,Δt1,分别为接收机与BD2系统时钟之间的钟差和钟漂,Δt2,分别为接收机与GPS系统时钟之间的钟差和钟漂,其初值本发明取为X0=[xu,vux,0,yu,vuy,0,zu,vuz,0,0,0,0,0]T,(xu,yu,zu)T,(vux,vuy,vuz)T分别为定位解算模块6输出的接收机位置坐标和速度坐标;
其中,第一系统状态转移模型
T为采样周期,本发明取为T=0.1s,
其中,第二系统状态转移模型
αx,αy,αz为机动时间常数的倒数,本发明取为αx=αy=αz=0.2Hz;
Qy(k),Qz(k)中的元素模型与Qx(k)中的元素模型相同,将Qx(k)中各元素模型中的αx分别用αy,αz代替,即可得到Qy(k),Qz(k),
分别为与接收机在参考坐标系x,y,z坐标方向上的加速度分量对应的系统噪声方差,本发明取为 αx,αy,αz为机动时间常数的倒数,本发明取为αx=αy=αz=0.2Hz,σΔt 2,分别为与接收机钟差和钟漂对应的系统噪声方差,本发明取为
其中,第一系统量测模型H(k)=[H1…Hj…Hn]T,
当卫星Sj属于BD2,
当卫星Sj属于GPS,
j为接收机定位时使用的有效卫星Sj的编号,j=1,2,…,n,n为接收机定位使用的有效卫星的数目,
其中,量测噪声模型 分别为卫星Sj的伪距观测噪声方差和伪距率观测噪声方差,分别为卫星S1的伪距观测噪声方差和伪距率观测噪声方差,分别为卫星Sn的伪距观测噪声方差和伪距率观测噪声方差,本发明取为 j=1,2,…,n;
其中,观测模型 ρu j,du j分别为对卫星Sj的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距和观测伪距率,j=1,2,…,n,ρu 1,du 1分别为对卫星S1的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距和观测伪距率,ρu n,du n分别为对卫星Sn的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距和观测伪距率;
其中,第二系统量测模型h(X)=[h1…hj…hn]T,
当卫星Sj属于BD2时,
当卫星Sj属于GPS时,
步骤五:将导航信息生成报文;
根据参数配置模块1输出的用户对导航报文的要求,报文形成模块8将数据更新频率为10Hz的导航信息以RINEX格式形成报文。
本发明的兼容导航接收机定位方法,单BD2系统独立工作可以在非常时期摆脱对其他卫星导航系统的依赖。多卫星导航系统组合定位可以大大增加有效卫星的数量,保证卫星定位服务的连续性、可用性。由于可见卫星数量大幅度增加,多卫星导航系统组合定位大大提高整个卫星定位系统的完善性监测性能,提高卫星导航定位应用的安全性。同时,多卫星导航系统组合定位可以增强观测卫星的几何图形强度,提高整个卫星定位系统的可靠性。
图3示出了BD2-GPS组合导航与单GPS定位仰角限为10度时的有效卫星数目曲线。图3中,实线表示应用本发明提供的兼容导航接收机定位方法进行组合定位仰角限为10度时的有效卫星数目,虚线表示应用单GPS进行定位仰角限为10度时的有效卫星数目。由图3可见,仰角限为10度时,单GPS全天有效卫星数目达5颗及以上,BD2-GPS组合系统全天有效卫星数目达13颗及以上,BD2-GPS组合星座为GPS星座有效卫星数目的两倍左右。
图4示出了BD2-GPS组合导航与单GPS定位仰角限为30度时的有效卫星数目曲线。图4中,实线表示应用本发明提供的兼容导航接收机定位方法进行组合定位仰角限为30度时的有效卫星数目,虚线表示应用单GPS进行定位仰角限为30度时的有效卫星数目。由图4可见,仰角限为30度时,单GPS全天有效卫星达4颗及以上的时刻占全天24h的90%,BD2-GPS组合系统全天有效卫星达7颗及以上,与单GPS相比,BD2-GPS组合系统把高度角的限制条件扩展了,即增强了卫星导航系统的可用性。
图5示出了BD2-GPS组合定位与单GPS定位仰角限为10度时的DOP曲线。图5中,实线表示应用本发明提供的兼容导航接收机定位方法进行组合定位仰角限为10度时的DOP,虚线表示应用单GPS进行定位仰角限为10度时的DOP。由图5可见,仰角限为10度时,单GPS的DOP值分布在1~6之间,波动较大,BD2-GPS组合系统DOP值明显减小,波动更加平缓。
图6示出了BD2-GPS组合定位与单GPS定位仰角限为30度时的DOP曲线。图6中,实线表示应用本发明提供的兼容导航接收机定位方法进行组合定位仰角限为30度时的DOP,虚线表示应用单GPS进行定位仰角限为30度时的DOP。由图6可见,仰角限为30度时,单GPS系统全天有10%时段观测不到4颗卫星,其他时段DOP较大,定位精度不能达到要求,BD2-GPS组合系统除部分超过4.0外,多介于2.0~4.0之间,波动平缓,精度明显优于单GPS,适合高精度的全天候导航定位。
图7示出了BD2-GPS组合定位与单GPS定位时的水平定位误差曲线。图中,实线表示应用本发明提供的兼容导航接收机定位方法进行组合定位时的水平定位误差,虚线表示应用单GPS进行定位时的水平定位误差。由图7可见,单GPS的水平定位精度低于BD2-GPS组合系统。
图8示出了BD2-GPS组合定位与单GPS定位时的垂直定位误差曲线。图中,实线表示应用本发明提供的兼容导航接收机定位方法进行组合定位时的垂直定位误差,虚线表示应用单GPS进行定位时的垂直定位误差。由图8可见,单GPS的垂直定位精度低于BD2-GPS组合系统。
Claims (3)
1、一种兼容导航接收机定位系统,其特征在于,该系统包括:参数配置模块(1)、数据预处理器(2)、卫星参数解算模块(3)、伪距修正模块(4)、坐标转换模块(5)、定位解算模块(6)、卡尔曼滤波器(7)和报文形成模块(8);
其中,参数配置模块(1)的第一个输出端与数据预处理器(2)的第一个输入端相连,参数配置模块(1)的第二个输出端与卫星参数解算模块(3)的第一个输入端相连,参数配置模块(1)的第三个输出端与伪距修正模块(4)的第一个输入端相连,参数配置模块(1)的第四个输出端与坐标转换模块(5)的第一个输入端相连,参数配置模块(1)的第五个输出端与定位解算模块(6)的第一个输入端相连,参数配置模块(1)的第六个输出端与卡尔曼滤波器(7)的第一个输入端相连,参数配置模块(1)的第七个输出端与报文形成模块(8)的第一个输入端相连;导航电文从数据预处理器(2)的第二个输入端输入,观测数据从数据预处理器(2)的第三个输入端输入,数据预处理器(2)的第一个输出端与卫星参数解算模块(3)的第二个输入端相连,数据预处理器(2)的第二个输出端与伪距修正模块(4)的第二个输入端相连;卫星参数解算模块(3)的输出端与坐标转换模块(5)的第二个输入端相连,坐标转换模块(5)的输出端与定位解算模块(6)的第二个输入端相连;伪距修正模块(4)的输出端与定位解算模块(6)的第三个输入端相连;定位解算模块(6)的输出端与卡尔曼滤波器(7)的第二个输入端相连;卡尔曼滤波器(7)的第一个输出端与伪距修正模块(4)的第三个输入端相连,卡尔曼滤波器(7)的第二个输出端与报文形成模块(8)的第二个输入端相连;
其中,参数配置模块(1)读取兼容导航接收机的配置参数,并将配置参数输出给定位系统的各模块以控制其工作状态;
数据预处理器(2)采集有效卫星的导航电文和观测数据帧,并从中提取星历数据、历书数据和观测信息;
卫星参数解算模块(3)获取有效卫星的卫星位置和卫星速度;
伪距修正模块(4)获取有效卫星的观测伪距和伪距误差修正信息,并对观测伪距进行修正;
坐标转换模块(5)将来自不同卫星导航系统的有效卫星的坐标统一到同一参考坐标系中;
定位解算模块(6)获取接收机位置和接收机速度信息;
卡尔曼滤波器(7)对接收机的位置和速度信息进行滤波;
报文形成模块(8)形成导航报文。
2、一种利用权利要求1所述的兼容导航接收机定位系统的定位方法,其特征在于,该方法包括如下步骤:
步骤一:读取接收机的配置参数;
参数配置模块(1)读取接收机的配置参数,
首先,根据用户需求,确定接收机的工作状态和使用的卫星导航系统;工作状态包括:单卫星导航系统定位、双卫星导航系统组合定位、三卫星导航系统组合定位和四卫星导航系统组合定位,使用的卫星导航系统包括:全球定位系统GPS、全球导航卫星系统GLONASS、北斗二代导航系统BD2和“伽利略”卫星定位系统Galileo;
然后,根据用户需求,确定卡尔曼滤波器(7)的工作状态,确定导航信息的数据更新频率,确定导航报文的格式;
最后,将接收机的配置参数输出给数据预处理器(2)、卫星参数解算模块(3)、伪距修正模块(4)、坐标转换模块(5)、定位解算模块(6)、卡尔曼滤波器(7)和报文形成模块(8),以控制各模块的工作状态;
其中,接收机的配置参数包括:接收机的工作状态、使用的卫星导航系统、卡尔曼滤波器(7)的工作状态、导航信息的数据更新频率以及导航报文的格式;
步骤二:提取星历数据、历书数据和观测数据;
根据参数配置模块(1)输出的接收机的工作状态和使用的卫星导航系统,数据预处理器(2)依据有效卫星所属的卫星导航系统的说明书对采集的有效卫星的导航电文进行拼接、校验、纠错,从中提取有效卫星的星历数据和历书数据;
根据参数配置模块(1)输出的接收机的工作状态和使用的卫星导航系统,数据预处理器(2)依据捕获跟踪系统提供的观测数据帧格式从采集的有效卫星观测数据帧中提取毫秒数、码片偏移、相位偏移和多普勒频率;
步骤三:获取有效卫星的参数信息;
首先,根据参数配置模块(1)输出的接收机的工作状态和使用的卫星导航系统,卫星参数解算模块(3)依据有效卫星所属的卫星导航系统的说明书利用有效卫星的星历数据获取其位置坐标和速度坐标;
然后,根据参数配置模块(1)输出的接收机的工作状态和使用的卫星导航系统,伪距修正模块(4)依据有效卫星所属卫星导航系统的说明书利用有效卫星的星历数据和历书数据获取其星钟修正误差、电离层修正误差和对流层修正误差;
之后,伪距修正模块(4)通过将有效卫星的周内秒计数、比特数、毫秒数、码片偏移、相位偏移进行累加获取其观测伪距;
最后,伪距修正模块(4)利用获取的有效卫星的星历修正误差、电离层修正误差和对流层修正误差对该卫星的观测伪距进行误差修正,获取修正后的观测伪距;
步骤四:获取接收机位置和接收机速度;
首先,根据参数配置模块(1)输出的接收机的工作状态和使用的卫星导航系统,在单卫星导航系统定位时坐标转换模块(5)不参与系统工作,在多卫星导航系统组合定位时坐标转换模块(5)依据各卫星导航系统参考坐标系的关系利用7参数Bursa模型将来自不同卫星导航系统的有效卫星的坐标转换到同一参考坐标系,获取各有效卫星在同一参考坐标系中的三维坐标;
然后,根据参数配置模块(1)输出的接收机的工作状态和使用的卫星导航系统,定位解算模块(6)利用有效卫星的参数信息获取接收机位置和接收机速度;
A.获取接收机位置:
当有效卫星的数目不小于N+3时,利用有效卫星的位置坐标和修正后的观测伪距获取接收机的位置坐标偏移和钟差模型:
其中,上述位置坐标偏移和钟差模型等号左侧部分为一个(N+3)×1的矩阵模块;
上述位置坐标偏移和钟差模型等号右侧部分的
为一个n×(N+3)的矩阵模块,
其中,i为卫星Sj所属卫星导航系统的编号,i=1,2,…,N,N为用户接收机进行定位时使用的卫星导航系统的数目;j为观测卫星Sj的编号,j=1,…o,o+1,…,p,…,q+1,…,n,n为观测卫星总数,S1~So为第一卫星导航系统的卫星,So+1~Sp为第二卫星导航系统的卫星,Sq+1~Sn为第N卫星导航系统的卫星;
(Δxu,Δyu,Δzu)T为用户接收机的位置坐标偏移;
Δtt为用户接收机时钟与第i卫星导航系统系统时钟之间的钟差,i=1,2,…,N,Δt1为用户接收机时钟与第一卫星导航系统系统时钟之间的钟差,……,ΔtN为用户接收机时钟与第N卫星导航系统系统时钟之间的钟差;
(xj,yj,zj)T为卫星Sj的位置坐标,j=1,…o,o+1,…,p,…,q+1,…,n;
ρu j为对卫星Sj的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距,j=1,…o,o+1,…,p,…,q+1,…,n,ρu 1为对第一卫星导航系统中的卫星S1的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距,ρu o为对第一卫星导航系统中的卫星So的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距,ρu o+1为对第二卫星导航系统中的卫星So+1的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距,ρu p为对第二卫星导航系统中的卫星Sp的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距,ρu q+1为对第N卫星导航系统中的卫星Sq+1的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距,ρu n为对第N卫星导航系统中的卫星Sn的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距;
为用户接收机的近似位置与卫星Sj间的几何距离,j=1,…o,o+1,…,p,…,q+1,…,n,为用户接收机的近似位置与第一卫星导航系统中的卫星S1间的几何距离,为用户接收机的近似位置与第一卫星导航系统中的卫星So间的几何距离,为用户接收机的近似位置与第二卫星导航系统中的卫星So+1间的几何距离,为用户接收机的近似位置与第二卫星导航系统中的卫星Sp间的几何距离,为用户接收机的近似位置与第N卫星导航系统中的卫星Sq+1间的几何距离,为用户接收机的近似位置与第N卫星导航系统中的卫星Sn间的几何距离;
c为光速;
用户接收机的位置坐标为:
B.获取接收机速度:
当有效卫星的数目不小于N+3时,利用接收机的位置坐标和有效卫星的多普勒频率获取接收机的速度坐标和钟漂模型:
其中,上述速度坐标和钟漂模型等号左侧部分为一个(N+3)×1的矩阵模块;
上述速度坐标和钟漂模型等号右侧部分的 为
一个n×(N+3)的矩阵模块,
d=[d1 … do do+1 … dp … dq+1 … dn]T,为一个n×1的矩阵模块,
其中
j为观测卫星Sj的编号,j=1,…o,o+1,…,p,…,q+1,…,n,其中n为观测卫星总数,S1~So为第一卫星导航系统的卫星,So+1~Sp为第二卫星导航系统的卫星,Sq+1~Sn为第N卫星导航系统的卫星;
i为卫星Sj所属卫星导航系统的编号,i=1,2,…,N,其中N为用户接收机进行定位时使用的卫星导航系统的数目;
(xu,yu,zu)T为用户接收机的位置坐标;
(xj,yj,zj)T为卫星Sj的位置坐标,j=1,…o,o+1,…,p,…,q+1,…,n;
Ru j为用户接收机与卫星Sj间的几何距离,j=1,…o,o+1,…,p,…,q+1,…,n;
(vux,vuy,vuz)T为用户接收机的速度坐标;
为用户接收机时钟与第i卫星导航系统系统时钟之间的钟差漂移,简称钟漂,i=1,2,…,N,为用户接收机时钟与第一卫星导航系统系统时钟之间的钟漂,……,为用户接收机时钟与第N卫星导航系统系统时钟之间的钟漂;
du j为卫星Sj的观测伪距率,j=1,…o,o+1,…,p,…,q+1,…,n,du 1为第一卫星导航系统中的卫星S1的观测伪距率,du o为第一卫星导航系统中的卫星So的观测伪距率,du o+1为第二卫星导航系统中的卫星So+1的观测伪距率,du p为第二卫星导航系统中的卫星Sp的观测伪距率,du q+1为卫第N星导航系统中的卫星Sq+1的观测伪距率,du n为第N卫星导航系统中的卫星Sn的观测伪距率;
(vx j,vy j,vz j)T为卫星Sj的速度坐标,j=1,…o,o+1,…,p,…,q+1,…,n;
c为光速;
最后,根据参数配置模块(1)输出的接收机的工作状态、使用的卫星导航系统和卡尔曼滤波器(7)的工作状态,确定卡尔曼滤波器(7)是否对定位解算模块(6)输出的接收机位置和接收机速度信息进行滤波处理;如果卡尔曼滤波器(7)的工作状态是工作,基于机动目标“当前”统计模型建立系统模型和测量模型,并通过自适应卡尔曼滤波获取接收机位置和接收机速度的滤波结果;
步骤五:将导航信息生成报文;
根据参数配置模块(1)输出的用户对导航报文的要求,报文形成模块(8)将导航信息以用户需要的格式形成报文。
3、根据权利要求2所述的一种兼容导航接收机定位方法,其特征在于,步骤四中所述的滤波结果采用递推方法给出:
估计均方误差模型P(k)=[I-K(k)H(k)]P(k|k-1);
其中,状态一步预测模型 一步预测均方误差模型P(k|k-1)=Φ(k|k-1)P(k-1)ΦT(k|k-1)+Q(k-1),滤波增益模型K(k)=P(k|k-1)HT(k)[H(k)P(k|k-1)HT(k)+R(k)]-1;
其中,系统状态模型 xu,分别为接收机在参考坐标系x轴方向的位置、速度和加速度,yu,分别为接收机在参考坐标系y轴方向的位置、速度和加速度,zu,分别为接收机在参考坐标系z轴方向的位置、速度和加速度,Δti,分别为接收机时钟与第i卫星导航系统系统时钟之间的钟差和钟漂,i为接收机定位使用的卫星导航系统的编号,i=1,2,…,N,N为接收机定位时使用卫星导航系统的数目,Δt1,分别为接收机时钟与第一卫星导航系统系统时钟之间的钟差和钟漂,……,ΔtN,分别为接收机时钟与第N卫星导航系统系统时钟之间的钟差和钟漂,初值为X0=[xu,vux,0,yu,vuy,0,zu,vuz,0,0,0,…,0,0]T,(xu,yu,zu)T和(vux,vuy,vuz)T分别为定位解算模块(6)输出的接收机位置坐标和速度坐标;
其中,第一系统状态转移模型
T为采样周期;
P(k-1)为状态协方差,初值P0取为对角阵;
其中,第二系统状态转移模型
αx,αy,αz为机动时间常数的倒数;
其中,系统噪声方差模型
Qy(k),Qz(k)中的元素模型与Qx(k)中的元素模型相同,将Qx(k)中各元素模型中的αx分别用αy,αz代替,即得到Qy(k),Qz(k);
分别为与接收机在参考坐标系x,y,z坐标方向上的加速度分量对应的系统噪声方差,αx,αy,αz为机动时间常数的倒数,σΔt 2,σt 2分别为与接收机钟差和钟漂对应的系统噪声方差;
其中,第一系统量测模型H(k)=[H1…Hj…Hn]T,
当卫星Sj属于第i卫星导航系统时,
当卫星Sj不属于第i卫星导航系统时,
j为接收机定位时使用的有效卫星Sj的编号,j=1,2,…,n,n为接收机定位使用的有效卫星的数目,
i为接收机定位时使用的卫星导航系统的编号,i=1,2,…,N,N为接收机定位使用的卫星导航系统的数目,
其中,量测噪声模型 分别为卫星Sj的伪距观测噪声方差和伪距率观测噪声方差,j=1,2,…,n,分别为卫星S1的伪距观测噪声方差和伪距率观测噪声方差,分别为卫星Sn的伪距观测噪声方差和伪距率观测噪声方差;
其中,观测模型 ρu j,du j分别为对卫星Sj的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距和观测伪距率,j=1,2,…,n,ρu 1,du 1分别为对卫星S1的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距和观测伪距率,ρu n,du n分别为对卫星Sn的观测伪距进行星钟误差、电离层误差和对流层误差修正后的伪距和观测伪距率;
其中,第二系统量测模型h(X)=[h1…hj…hn]T
当卫星Sj属于第i卫星导航系统时,
当卫星Sj不属于第i卫星导航系统时,
j为接收机定位时使用的有效卫星Sj的编号,j=1,2,…,n,其中n为接收机定位使用的有效卫星的数目,
i为接收机定位时使用的卫星导航系统的编号,i=1,2,…,N,N为接收机定位使用的卫星导航系统的数目。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910089221XA CN101609140B (zh) | 2009-07-09 | 2009-07-09 | 一种兼容导航接收机定位系统及其定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910089221XA CN101609140B (zh) | 2009-07-09 | 2009-07-09 | 一种兼容导航接收机定位系统及其定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101609140A true CN101609140A (zh) | 2009-12-23 |
CN101609140B CN101609140B (zh) | 2011-07-27 |
Family
ID=41482965
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910089221XA Expired - Fee Related CN101609140B (zh) | 2009-07-09 | 2009-07-09 | 一种兼容导航接收机定位系统及其定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101609140B (zh) |
Cited By (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101975962A (zh) * | 2010-09-26 | 2011-02-16 | 东莞市泰斗微电子科技有限公司 | 一种卫星导航定位方法及相应的装置 |
CN102540220A (zh) * | 2011-12-30 | 2012-07-04 | 合众思壮北斗导航有限公司 | 卫星定位方法及设备 |
CN102540228A (zh) * | 2012-03-02 | 2012-07-04 | 重庆九洲星熠导航设备有限公司 | 一种单频gps高精度单点定位系统及方法 |
CN102749633A (zh) * | 2012-06-29 | 2012-10-24 | 北京航空航天大学 | 一种卫星导航接收机的动态定位解算方法 |
CN103037500A (zh) * | 2011-09-29 | 2013-04-10 | 北京信威通信技术股份有限公司 | 一种卫星通信中的同步方法 |
CN103487821A (zh) * | 2012-06-12 | 2014-01-01 | 广州中海达卫星导航技术股份有限公司 | 一种基线向量解算方法 |
CN104020478A (zh) * | 2014-05-13 | 2014-09-03 | 沈阳理工大学 | 基于闭合回路控制的新一代gnss同步接收模型 |
CN104412066A (zh) * | 2012-06-27 | 2015-03-11 | 三菱电机株式会社 | 定位装置 |
CN104849697A (zh) * | 2015-05-15 | 2015-08-19 | 重庆大学 | 一种基于去偏坐标转换的α-β滤波方法 |
CN104931932A (zh) * | 2015-05-28 | 2015-09-23 | 重庆大学 | 一种改进的去偏坐标转换卡尔曼滤波方法 |
CN105510935A (zh) * | 2015-12-15 | 2016-04-20 | 北京光华纺织集团有限公司 | 一种用于导航的装置 |
CN105572705A (zh) * | 2015-12-15 | 2016-05-11 | 北京光华纺织集团有限公司 | 一种导航装置 |
CN105572706A (zh) * | 2015-12-15 | 2016-05-11 | 北京光华纺织集团有限公司 | 一种自动导航驾驶系统 |
CN105842720A (zh) * | 2016-03-23 | 2016-08-10 | 哈尔滨工程大学 | 一种基于载波相位的广域精密实时定位方法 |
CN104412065B (zh) * | 2012-06-27 | 2016-10-05 | 三菱电机株式会社 | 定位装置 |
CN103037500B (zh) * | 2011-09-29 | 2016-12-14 | 北京信威通信技术股份有限公司 | 一种卫星通信中的同步方法 |
CN107238846A (zh) * | 2017-04-25 | 2017-10-10 | 清华大学 | 一种基于glonass历书参数的卫星位置与速度预报方法 |
CN107783937A (zh) * | 2017-10-19 | 2018-03-09 | 西安科技大学 | 一种求解任意旋转角三维坐标转换参数的方法 |
CN108152838A (zh) * | 2017-12-15 | 2018-06-12 | 中国航空工业集团公司洛阳电光设备研究所 | 一种基于观瞄测定目标位置的装置和方法 |
CN108279007A (zh) * | 2018-01-23 | 2018-07-13 | 哈尔滨工业大学 | 一种基于随机信号的定位方法及装置 |
CN108803374A (zh) * | 2018-06-07 | 2018-11-13 | 中国人民解放军海军工程大学 | 一种无人艇环境数据仿真方法 |
CN109257085A (zh) * | 2018-03-30 | 2019-01-22 | 北京润科通用技术有限公司 | 人造卫星与飞行设备间多普勒频移的获得方法及装置 |
CN111123320A (zh) * | 2019-12-31 | 2020-05-08 | 泰斗微电子科技有限公司 | 一种卫星定位装置、卫星信号接收机及终端设备 |
CN111694033A (zh) * | 2019-12-31 | 2020-09-22 | 泰斗微电子科技有限公司 | 移动装置的位置确定方法、装置及移动装置 |
CN111965685A (zh) * | 2020-07-07 | 2020-11-20 | 北京自动化控制设备研究所 | 一种基于多普勒信息的低轨卫星/惯性组合导航定位方法 |
CN112285749A (zh) * | 2020-10-21 | 2021-01-29 | 腾讯科技(深圳)有限公司 | 全球导航卫星系统原始观测数据处理方法、装置及存储介质 |
CN112731504A (zh) * | 2020-11-17 | 2021-04-30 | 中国人民解放军63921部队 | 对月球探测器自主定轨的方法及装置 |
CN113758373A (zh) * | 2021-09-07 | 2021-12-07 | 重庆天箭惯性科技股份有限公司 | 一种提高弹载接收机定位、测速精度的方法、装置及设备 |
CN113866807A (zh) * | 2021-12-03 | 2021-12-31 | 深圳市麦微智能电子有限公司 | 一种基于bds定位的高精度gnss导航系统 |
CN115015979A (zh) * | 2022-07-01 | 2022-09-06 | 山东大学 | 一种基于自适应动态切换的导航信号接收处理系统及方法 |
US20230003906A1 (en) * | 2021-06-30 | 2023-01-05 | Beken Corporation | Apparatus and method of calculating position-velocity-time results of receiver |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110071784B (zh) * | 2019-04-25 | 2020-08-18 | 电子科技大学 | 一种ads-b航空监视系统中的航空器定位方法 |
-
2009
- 2009-07-09 CN CN200910089221XA patent/CN101609140B/zh not_active Expired - Fee Related
Cited By (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101975962A (zh) * | 2010-09-26 | 2011-02-16 | 东莞市泰斗微电子科技有限公司 | 一种卫星导航定位方法及相应的装置 |
CN101975962B (zh) * | 2010-09-26 | 2012-11-14 | 东莞市泰斗微电子科技有限公司 | 一种卫星导航定位方法及相应的装置 |
CN103037500B (zh) * | 2011-09-29 | 2016-12-14 | 北京信威通信技术股份有限公司 | 一种卫星通信中的同步方法 |
CN103037500A (zh) * | 2011-09-29 | 2013-04-10 | 北京信威通信技术股份有限公司 | 一种卫星通信中的同步方法 |
CN102540220A (zh) * | 2011-12-30 | 2012-07-04 | 合众思壮北斗导航有限公司 | 卫星定位方法及设备 |
CN102540228A (zh) * | 2012-03-02 | 2012-07-04 | 重庆九洲星熠导航设备有限公司 | 一种单频gps高精度单点定位系统及方法 |
CN102540228B (zh) * | 2012-03-02 | 2015-09-16 | 重庆九洲星熠导航设备有限公司 | 一种单频gps高精度单点定位系统及方法 |
CN103487821B (zh) * | 2012-06-12 | 2015-10-28 | 广州中海达卫星导航技术股份有限公司 | 一种基线向量解算方法 |
CN103487821A (zh) * | 2012-06-12 | 2014-01-01 | 广州中海达卫星导航技术股份有限公司 | 一种基线向量解算方法 |
CN104412066B (zh) * | 2012-06-27 | 2016-03-16 | 三菱电机株式会社 | 定位装置 |
CN104412066A (zh) * | 2012-06-27 | 2015-03-11 | 三菱电机株式会社 | 定位装置 |
US9864064B2 (en) | 2012-06-27 | 2018-01-09 | Mitsubishi Electric Corporation | Positioning device |
CN104412065B (zh) * | 2012-06-27 | 2016-10-05 | 三菱电机株式会社 | 定位装置 |
CN102749633B (zh) * | 2012-06-29 | 2013-11-27 | 北京航空航天大学 | 一种卫星导航接收机的动态定位解算方法 |
CN102749633A (zh) * | 2012-06-29 | 2012-10-24 | 北京航空航天大学 | 一种卫星导航接收机的动态定位解算方法 |
CN104020478A (zh) * | 2014-05-13 | 2014-09-03 | 沈阳理工大学 | 基于闭合回路控制的新一代gnss同步接收模型 |
CN104849697B (zh) * | 2015-05-15 | 2017-11-14 | 上海前隆信息科技有限公司 | 一种基于去偏坐标转换的α‑β滤波方法 |
CN104849697A (zh) * | 2015-05-15 | 2015-08-19 | 重庆大学 | 一种基于去偏坐标转换的α-β滤波方法 |
CN104931932A (zh) * | 2015-05-28 | 2015-09-23 | 重庆大学 | 一种改进的去偏坐标转换卡尔曼滤波方法 |
CN105510935A (zh) * | 2015-12-15 | 2016-04-20 | 北京光华纺织集团有限公司 | 一种用于导航的装置 |
CN105572705A (zh) * | 2015-12-15 | 2016-05-11 | 北京光华纺织集团有限公司 | 一种导航装置 |
CN105572706A (zh) * | 2015-12-15 | 2016-05-11 | 北京光华纺织集团有限公司 | 一种自动导航驾驶系统 |
CN105842720A (zh) * | 2016-03-23 | 2016-08-10 | 哈尔滨工程大学 | 一种基于载波相位的广域精密实时定位方法 |
CN107238846A (zh) * | 2017-04-25 | 2017-10-10 | 清华大学 | 一种基于glonass历书参数的卫星位置与速度预报方法 |
CN107783937A (zh) * | 2017-10-19 | 2018-03-09 | 西安科技大学 | 一种求解任意旋转角三维坐标转换参数的方法 |
CN108152838A (zh) * | 2017-12-15 | 2018-06-12 | 中国航空工业集团公司洛阳电光设备研究所 | 一种基于观瞄测定目标位置的装置和方法 |
CN108152838B (zh) * | 2017-12-15 | 2024-04-16 | 中国航空工业集团公司洛阳电光设备研究所 | 一种基于观瞄测定目标位置的装置和方法 |
CN108279007A (zh) * | 2018-01-23 | 2018-07-13 | 哈尔滨工业大学 | 一种基于随机信号的定位方法及装置 |
CN108279007B (zh) * | 2018-01-23 | 2020-03-17 | 哈尔滨工业大学 | 一种基于随机信号的定位方法及装置 |
CN109257085A (zh) * | 2018-03-30 | 2019-01-22 | 北京润科通用技术有限公司 | 人造卫星与飞行设备间多普勒频移的获得方法及装置 |
CN108803374B (zh) * | 2018-06-07 | 2021-09-21 | 中国人民解放军海军工程大学 | 一种无人艇环境数据仿真方法 |
CN108803374A (zh) * | 2018-06-07 | 2018-11-13 | 中国人民解放军海军工程大学 | 一种无人艇环境数据仿真方法 |
CN111694033A (zh) * | 2019-12-31 | 2020-09-22 | 泰斗微电子科技有限公司 | 移动装置的位置确定方法、装置及移动装置 |
CN111123320A (zh) * | 2019-12-31 | 2020-05-08 | 泰斗微电子科技有限公司 | 一种卫星定位装置、卫星信号接收机及终端设备 |
CN111123320B (zh) * | 2019-12-31 | 2022-05-27 | 泰斗微电子科技有限公司 | 一种卫星定位装置、卫星信号接收机及终端设备 |
CN111965685B (zh) * | 2020-07-07 | 2023-01-13 | 北京自动化控制设备研究所 | 一种基于多普勒信息的低轨卫星/惯性组合导航定位方法 |
CN111965685A (zh) * | 2020-07-07 | 2020-11-20 | 北京自动化控制设备研究所 | 一种基于多普勒信息的低轨卫星/惯性组合导航定位方法 |
CN112285749A (zh) * | 2020-10-21 | 2021-01-29 | 腾讯科技(深圳)有限公司 | 全球导航卫星系统原始观测数据处理方法、装置及存储介质 |
CN112731504A (zh) * | 2020-11-17 | 2021-04-30 | 中国人民解放军63921部队 | 对月球探测器自主定轨的方法及装置 |
US11662476B2 (en) * | 2021-06-30 | 2023-05-30 | Beken Corporation | Apparatus and method of calculating position-velocity-time results of receiver |
US20230003906A1 (en) * | 2021-06-30 | 2023-01-05 | Beken Corporation | Apparatus and method of calculating position-velocity-time results of receiver |
CN113758373A (zh) * | 2021-09-07 | 2021-12-07 | 重庆天箭惯性科技股份有限公司 | 一种提高弹载接收机定位、测速精度的方法、装置及设备 |
CN113758373B (zh) * | 2021-09-07 | 2023-01-10 | 重庆天箭惯性科技股份有限公司 | 一种提高弹载接收机定位、测速精度的方法、装置及设备 |
CN113866807A (zh) * | 2021-12-03 | 2021-12-31 | 深圳市麦微智能电子有限公司 | 一种基于bds定位的高精度gnss导航系统 |
CN115015979A (zh) * | 2022-07-01 | 2022-09-06 | 山东大学 | 一种基于自适应动态切换的导航信号接收处理系统及方法 |
CN115015979B (zh) * | 2022-07-01 | 2024-09-27 | 山东大学 | 一种基于自适应动态切换的导航信号接收处理系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101609140B (zh) | 2011-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101609140B (zh) | 一种兼容导航接收机定位系统及其定位方法 | |
Li et al. | Review of PPP–RTK: Achievements, challenges, and opportunities | |
CN109059751B (zh) | 一种形变数据监测方法及系统 | |
CN101395443B (zh) | 混合定位方法和设备 | |
Konrad et al. | Advanced state estimation for navigation of automated vehicles | |
CN101403790A (zh) | 单频gps接收机的精密单点定位方法 | |
CN106291639A (zh) | 一种gnss接收机实现定位的方法及装置 | |
CN110133700B (zh) | 一种船载综合导航定位方法 | |
CN104316943B (zh) | 一种伪距离和多普勒组合差分定位系统及方法 | |
CN106405592B (zh) | 车载北斗载波相位周跳检测与修复方法及系统 | |
EP2193336B1 (en) | A positioning system and method | |
CN103543454B (zh) | 一种嵌入在移动通讯网中的卫星定轨系统 | |
CN112129300A (zh) | 位置间动力学约束的低轨卫星星载gnss精密定轨方法及系统 | |
CN110007326A (zh) | 一种用于星基增强系统的双频测距误差参数生成方法 | |
Yoder et al. | Low-Cost Inertial Aiding for Deep-Urban Tightly Coupled Multi-Antenna Precise GNSS | |
WO2019144480A1 (zh) | 基于速度约束的低成本接收机平滑rtd算法 | |
JP2008039691A (ja) | 搬送波位相式測位装置 | |
Wang et al. | Velocity estimations by combining time-differenced GPS and Doppler observations | |
US20240077620A1 (en) | System and method for determining gnss corrections | |
Tang et al. | GNSS localization propagation error estimation considering environmental conditions | |
Wang et al. | A study of GPS PPP velocity estimation with the constraints of Doppler observation | |
CN106054227A (zh) | 惯导辅助下的伪距差值单星高动态定位方法 | |
CN103389502B (zh) | 基于多个地面基站高精度确定载体加速度的方法 | |
Beran et al. | High-accuracy point positioning with low-cost GPS receivers: How good can It get? | |
Zhang et al. | Research on virtual reference station (VRS) technology and differential correction information |
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: 20110727 Termination date: 20140709 |
|
EXPY | Termination of patent right or utility model |