CN102749633B - 一种卫星导航接收机的动态定位解算方法 - Google Patents
一种卫星导航接收机的动态定位解算方法 Download PDFInfo
- Publication number
- CN102749633B CN102749633B CN2012102256259A CN201210225625A CN102749633B CN 102749633 B CN102749633 B CN 102749633B CN 2012102256259 A CN2012102256259 A CN 2012102256259A CN 201210225625 A CN201210225625 A CN 201210225625A CN 102749633 B CN102749633 B CN 102749633B
- Authority
- CN
- China
- Prior art keywords
- error
- matrix
- model
- pseudorange
- vector
- 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
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提出一种卫星导航接收机的动态定位解算方法,属于卫星导航技术领域,该动态定位解算方法包括以下几个步骤:建立以机动加速度作为未知输入的载体运动模型作为状态方程;建立导航接收机伪距和伪距率的测量方程;基于预测滤波估计[tk,tk+1]区间的未知机动加速度;实现载体导航定位参数的预测;利用测量残差对预测状态进行修正,实现载体导航定位参数的估计。本发明采用的改进的预测滤波可以在线实时估计机动加速度并修正载体运动模型,克服了传统动态定位解算过程中将系统误差处理为零均值白噪声的缺点;不仅可以在线估计载体机动加速度,而且可以获得较好的定位效果,特别适合于载体机动未知或较大机动的情况。
Description
技术领域
本发明涉及一种卫星导航接收机的动态定位解算方法,属于卫星导航技术领域。
背景技术
导航定位解算是卫星导航接收机的重要组成部分,利用伪距和伪距率测量来计算载体的位置、速度、加速度等导航定位参数。常用的定位方法有最小二乘法(LS)和扩展卡尔曼滤波法(EKF)。其中最小二乘法(LS)简单不考虑载体的运动状态,定位精度较低;而扩展卡尔曼滤波法(EKF)需要对载体运动进行建模,定位精度优于最小二乘法(LS)。测量模型的非线性和系统模型的不精确是制约卫星导航接收机动态导航定位精度的两个主要因素。由于伪距和伪距率测量模型的非线性强,无迹卡尔曼滤波(UKF)对非线性具有良好的逼近效果,因此基于无迹卡尔曼滤波(UKF)的定位解算精度优于扩展卡尔曼滤波法(EKF)。但是载体在实际运动中存在多种形式的机动,很难建立精确的模型,特别是对于机动性强的载体,即使采用Singer模型、“当前”统计模型也难以准确描述载体运动的动态过程。如果动态定位所用运动模型不准确,尽管可以通过调整系统噪声的方差阵来尽可能的减少模型误差的影响,但是由于机动带来的模型误差并不满足白噪声特性,必将导致基于EKF和UKF的动态定位精度下降,因此该方法难以满足任意机动形式下的导航解算的需要。非线性预测滤波(NPF)是可以在线估计和修正模型误差的非线性递推状态估计方法。但NPF收敛速度较慢,且对初值误差较敏感,若NPF对模型误差估计不准确,通过状态传播还会导致较大的状态误差,甚至造成滤波发散。
发明内容
本发明的目的是为了解决载体机动情况未知或机动较大时扩展卡尔曼滤波法(EKF)的卫星导航接收机动态定位解算精度降低的问题,提出一种基于改进的预测滤波(IPF)的卫星导航接收机动态定位解算方法。
本发明是一种卫星导航接收机的动态定位解算方法,包括以下几个步骤:
步骤一:建立以机动加速度作为输入的载体运动方程,用该载体运动方程构成系统状态方程:
采用以机动加速度作为输入的载体运动模型,同时利用导航接收机的时钟误差来建立动态定位解算的状态方程,如下:
其中,P=[x y z]T和V=[vx vy vz]T分别为接收机在导航参考坐标系下位置向量和速度向量;x、y、z和vx、vy、vz分别表示接收机在导航参考坐标系中的沿x、y、z轴的位置分量和速度分量,ain=[ax ay az]T为由载体机动引起的未建模加速度向量,其中ax、ay、az表示在导航参考坐标系的沿x、y、z轴的加速度分量;tb为接收机的等效时钟误差;td为接收机的等效时钟频率误差;系统噪声向量由速度噪声wv、加速度噪声wa、时钟误差噪声wtb和时钟频率误差wtd构成,均为零均值的高斯白噪声,方差分别为σv、σa、σtb、σtd。
步骤二:建立卫星导航接收机伪距和伪距率的观测方程:
利用接收机测量得到的伪距和伪距率构建测量向量y如下:
设第i颗卫星与接收机的几何距离Ri为:
其中,xsi、ysi、zsi为在导航参考坐标系下第i颗卫星沿x、y、z轴的位置分量;
ρi=Ri+tb+vi (4)
步骤三:利用tk+1时刻的伪距、伪距率观测量,基于预测滤波获得[tk,tk+1]区间的未知机动加速度:
将动态导航定位的系统模型表示为如下的非线性模型:
其中,为模型向量,表示n维实空间;t为连续系统的时间变量;对于动态导航定位,f=[vx vy vz 000td 0]T,n=8;为状态向量;为模型误差向量,表示p维实空间,d(t)=ain,p=3;为模型误差扰动矩阵,表示n×p维的实空间, I3×3表示3×3的单位阵;为测量向量,表示m维的实空间;为测量输出向量,m=2×num;系统噪声w(t)是均值为零,方差为Q的高斯白噪声;测量噪声v(t)是均值为零、方差为R的高斯白噪声。
其中为tk时刻的系统状态估计;为由根据测量方程得到的tk时刻的测量估计;R为观测噪声方差阵;yk+1为tk+1时刻的实际测量向量;T为滤波周期;为灵敏度矩阵;Λ(T)为对角矩阵;为列向量;W为加权矩阵。
灵敏度矩阵 其中分别与第i颗卫星的伪距、伪距率测量对应的行向量为:
对角阵Λ(T)为:
列向量 其中分别与第i颗卫星的伪距、伪距率测量对应的元素为:
最终根据公式(6)得到[tk,tk+1]区间的机动加速度估计。
步骤四:利用[tk,tk+1]区间的机动加速度估计修正系统模型,并进行时间更新,预测载体导航定位参数:
将得到的机动加速度作为系统输入代入系统模型公式(1)中,并对状态进行时间更新,更新过程如公式(10)、(11)所示,状态及其误差协方差阵的一步预测为:
Pk+1/k=(Ak-BkMkHk)Pk(Ak-BkMkHk)T+Q′k (11)
其中为tk+1时刻的状态一步预测;Pk+1/k为tk+1时刻的状态误差协方差阵一步预测;为tk时刻模型向量的估计值;f(x,tk)为tk时刻的模型向量;为[tk,tk+1]区间的模型误差估计;gd为误差扰动矩阵;Ak、Bk和Mk为实现表达而出现的变量,它们是
步骤五:进行测量更新,利用测量残差修正预测状态,实现载体导航定位参数的估计:
如式(12)、(13)、(14)所示,tk+1时刻的增益矩阵Kk+1为:
本发明的优点在于:
(1)本发明提出一种卫星导航接收机的动态定位解算方法,采用的改进的预测滤波(IPF)可以在线实时估计机动加速度并修正载体运动模型,克服了传统动态定位解算过程中将系统误差处理为零均值白噪声的缺点;
(2)本发明提出一种卫星导航接收机的动态定位解算方法,采用了改进的预测滤波(IPF)不仅可以在线估计载体机动加速度,而且可以获得较好的定位效果,特别适合于载体机动未知或较大机动的情况,并可以推广用于以载波相位和载波相位变化率作为观测量的接收机动态定位的场合。
附图说明
图1是本发明提出的一种卫星导航接收机的动态定位解算方法的流程图;
图2a是无初始误差时IPF(CV)与EKF(CV)和EKF(Singer)方法仿真结果的位置误差曲线对比图;
图2b是无初始误差时IPF(CV)与EKF(CV)和EKF(Singer)方法仿真结果的速度误差曲线对比图;
图2c是无初始误差时IPF(CV)与EKF(CV)和EKF(Singer)方法仿真结果的加速度误差曲线对比图;
图3a是有初始误差时IPF(CV)与EKF(CV)和EKF(Singer)方法仿真结果的位置误差曲线对比图;
图3b是有初始误差时IPF(CV)与EKF(CV)和EKF(Singer)方法仿真结果的速度误差曲线对比图;
图3c是有初始误差时IPF(CV)与EKF(CV)和EKF(Singer)方法仿真结果的加速度误差曲线对比图;
图4是仿真时的载体飞行三维轨迹图;
图5是系统加速度输入(ain)的三轴变化曲线图;
具体实施方式
下面将结合附图对本发明作进一步的详细说明。
本发明提出一种卫星导航接收机的动态定位解算方法,利用预测滤波来实时在线估计并补偿所建载体运动模型误差—机动加速度,再利用EKF进行时间更新和测量更新,得到卫星接收机的导航参数,如图1所示,包括以下几个步骤:
步骤一:建立以机动加速度作为输入的载体运动方程,用该载体运动方程构成系统状态方程。
卫星导航接收机导航方法中常用的动态定位模型有常速度模型(CV模型)和Singer模型。其中,CV模型假设载体速度变化满足随机游走模型,当载体的速度发生较大变化时,将会造成较大的模型误差。Singer模型假设加速度变化满足一阶马尔科夫随机过程模型,用机动频率α表示载体的机动程度。当载体以复杂的机动形式运动时,将导致所用载体运动模型存在模型误差。
本发明中采用以机动加速度作为输入的载体运动模型,同时考虑导航接收机的时钟误差来建立动态定位解算的状态方程,如下:
其中,P=[x y z]T和V=[vx vy vz]T分别为接收机在导航参考坐标系下位置向量和速度向量;x、y、z和vx、vy、vz分别表示接收机在导航参考坐标系中的沿x、y、z轴的位置分量和速度分量,ain=[ax ay az]T为由载体机动引起的未建模加速度向量,其中ax、ay、az表示在导航参考坐标系的沿x、y、z轴的加速度分量;tb为接收机的等效时钟误差;td为接收机的等效时钟频率误差;系统噪声向量由速度噪声wv、加速度噪声wa、时钟误差噪声wtb和时钟频率误差wtd构成,均为零均值的高斯白噪声,方差分别为σv、σa、σtb、σtd。
步骤二:建立卫星导航接收机伪距和伪距率的观测方程。
利用接收机测量得到的伪距和伪距率构建测量向量y如下:
其中,ρi、分别为第i颗卫星的伪距、伪距率测量值(i=1,…num);m为用于动态导航解算的可见卫星颗数,如果采用几何精度因子(GDOP)最小的4颗卫星的测量用于动态定位,则num=4。
设第i颗卫星与接收机的几何距离Ri为:
其中,xsi、ysi、zsi为在导航参考坐标系下第i颗卫星沿x、y、z轴的位置分量。
利用接收机测量得到的伪距和伪距率作为观测量,由此建立由第i颗卫星数据得到的伪距ρi和伪距率测量方程分别如下:
ρi=Ri+tb+vi (4)
步骤三:利用tk+1时刻的伪距、伪距率观测量,基于预测滤波获得[tk,tk+1]区间的未知机动加速度估计。
将动态导航定位的系统模型表示为如下的非线性模型:
其中,为模型向量,表示n维实空间;t为连续系统的时间变量;对于动态导航定位,f=[vx vy vz 000td 0]T,n=8;为状态向量;为模型误差向量,表示p维实空间,d(t)=ain,p=3;为模型误差扰动矩阵,表示n×p维的实空间, I3×3表示3×3的单位阵;为测量向量,表示m维的实空间;为测量输出向量,m=2×num;系统噪声w(t)是均值为零,方差为Q的高斯白噪声;测量噪声v(t)是均值为零、方差为R的高斯白噪声。
以下标k来表示tk时刻,根据预测滤波方法可得[tk,tk+1]区间的模型误差为:
其中,为tk时刻的系统状态估计;为由根据测量方程得到的tk时刻的测量估计;R为观测噪声方差阵;yk+1为tk+1时刻的实际测量向量;T为滤波周期;为灵敏度矩阵;Λ(T)为对角矩阵;为列向量;W为加权矩阵。
加权矩阵W可预先设定数值,其选取原则是首先确定滤波器稳定的W的取值范围,然后在这个范围内平衡滤波估计的噪声影响和估计偏差的大小,来选择W矩阵。本发明中考虑这些因素后设定
灵敏度矩阵 其中分别与第i颗卫星的伪距、伪距率测量对应的行向量为:
对角阵Λ(T)为:
列向量 其中分别与第i颗卫星的伪距、伪距率测量对应的元素为:
最终根据公式(6)可以得到[tk,tk+1]区间的机动加速度估计。
步骤四:利用[tk,tk+1]区间的机动加速度估计修正系统模型,并进行时间更新,预测载体导航定位参数。
将得到的机动加速度作为系统输入代入系统模型即公式(1)中,并对状态进行时间更新,更新过程如公式(10)、(11)所示,状态及其误差协方差阵的一步预测为:
Pk+1/k=(Ak-BkMkHk)Pk(Ak-BkMkHk)T+Q′k (11)
其中为tk+1时刻的状态一步预测;Pk+1/k为tk+1时刻的状态误差协方差阵一步预测;为tk时刻模型向量的估计值;f(x,tk)为tk时刻的模型向量;为[tk,tk+1]区间的模型误差估计;gd为误差扰动矩阵;Ak、Bk和Mk为实现表达而出现的变量,它们是 I为单位阵,d为模型误差;Hk为基于状态得到的测量模型的雅可比矩阵,Q′k为对称正定矩阵;Pk为tk时刻状态误差方差阵。
式(11)中的Q′k是一个特别引入的对称正定矩阵,可以看作经过估计模型误差补偿后的等效系统噪声的方差阵,由于一部分系统噪声的影响在模型误差估计中被等效估计并补偿,因此Q′k的取值略小于系统噪声方差阵Qk,Q′k=λQk,λ=0.8~1,λ为修正系数。
步骤五:进行测量更新,利用测量残差修正预测状态,实现载体导航定位参数的估计。
该过程如式(12)、(13)、(14)所示,tk+1时刻的增益矩阵Kk+1为:
本发明提出的一种卫星导航接收机的动态定位解算方法,具有对载体机动加速度实时估计并利用估计的模型误差修正状态估计的能力,因此可以得到更高精度的导航参数估计。同时,本方法还可以适用于其他类似带有模型误差的系统。
在存在模型误差时,本发明提出的一种卫星导航接收机的动态定位解算方法为CV模型下基于IPF的定位方法IPF(CV),与CV模型下基于EKF的定位方法EKF(CV)、Singer模型下基于EKF的定位方法EKF(Singer)在系统无初始误差、有初始误差条件下的仿真结果误差曲线对比分别如图2和图3所示。图4和图5分别为仿真时载体轨迹的三维图和系统加速度输入(ain)的三轴变化曲线图。
仿真参数设置如下:
(1)滤波周期T=1s。
(2)
(3)σv=10m,σa=6m/s,σtb=0.8m,σtd=0.05m/s。
(4)仿真次数en=200。
当存在模型误差而无初始状态误差时,在同样的仿真条件下,本发明提出的CV模型下基于IPF的定位方法IPF(CV)、CV模型下基于EKF的定位方法EKF(CV),Singer模型下基于EKF的定位方法EKF(Singer)的位置、速度、加速度的误差曲线分别如图2a、图2b、图2c所示。当同时存在模型误差和初始误差时,在同样地仿真条件下,本发明提出的CV模型下基于IPF的定位方法IPF(CV)、CV模型下基于EKF的定位方法EKF(CV)Singer模型下基于EKF的定位方法EKF(Singer)的位置、速度、加速度的误差曲线分别如图3a、图3b、图3c所示。
图2a和图3a中横坐标均是仿真时间,纵坐标是位置估计误差曲线(真实值与估计值之差),从图中可以看出,EKF(CV)和Singer模型下基于EKF的定位方法EKF(Singer)的位置估计误差发散,而本发明提出的CV模型下基于IPF的定位方法IPF(CV)的位置估计误差在10m以内。
图2b和图3b中横坐标均是仿真时间,纵坐标是速度的估计误差曲线;图2c和图3c的横坐标均是仿真时间,纵坐标是加速度的估计误差曲线;从图中可以看出,CV模型下基于EKF的定位方法EKF(CV)和Singer模型下基于EKF的定位方法EKF(Singer)的速度估计误差震荡比较剧烈,而本发明CV模型下基于IPF的定位方法IPF(CV)估计误差因为能够实时估计模型误差并补偿,所以速度估计误差仍然在1m/s以内,且加速度估计误差在0.7m/s2以内;因此本发明提出CV模型下基于IPF的定位方法IPF(CV)精度有明显提高。三种方法的估计效果如表1、表2所示:
表1:无初始状态误差的滤波精度
表2:有初始状态误差的滤波精度
从表1和表2中可以看出,无论是否有初始误差,本发明提出的CV模型下基于IPF的定位方法IPF(CV)的定位精度明显优于CV模型下基于EKF的定位方法EKF(CV)和Singer模型下基于EKF的定位方法EKF(Singer)。
Claims (4)
1.一种卫星导航接收机的动态定位解算方法,其特征在于:包括以下几个步骤:
步骤一:建立以机动加速度作为输入的载体运动方程,用该载体运动方程构成系统状态方程:
采用以机动加速度作为输入的载体运动模型,同时利用导航接收机的时钟误差来建立动态定位解算的状态方程,如下:
其中,P=[x y z]T和V=[vx vy vz]T分别为接收机在导航参考坐标系下位置向量和速度向量;x、y、z和vx、vy、vz分别表示接收机在导航参考坐标系中的沿x、y、z轴的位置分量和速度分量,ain=[ax ay az]T为由载体机动引起的未建模加速度向量,其中ax、ay、az表示在导航参考坐标系的沿x、y、z轴的加速度分量;tb为接收机的等效时钟误差;td为接收机的等效时钟频率误差;系统噪声向量由速度噪声wv、加速度噪声wa、时钟误差噪声wtb和时钟频率误差wtd构成,均为零均值的高斯白噪声,方差分别为σv、σa、σtb、σtd;
步骤二:建立卫星导航接收机伪距和伪距率的观测方程:
利用接收机测量得到的伪距和伪距率构建测量向量y如下:
设第i颗卫星与接收机的几何距离Ri为:
其中,xsi、ysi、zsi为在导航参考坐标系下第i颗卫星沿x、y、z轴的位置分量;
ρi=Ri+tb+vi (4)
步骤三:利用tk+1时刻的伪距、伪距率观测量,基于预测滤波获得[tk,tk+1]区间的未知机动加速度:
将动态导航定位的系统模型表示为如下的非线性模型:
其中, 为模型向量, 表示n维实空间;t为连续系统的时间变量;对于动态导航定位,f=[vx vy vz 000td 0]T,n=8; 为状态向量; 为模型误差向量, 表示p维实空间,d(t)=ain,p=3; 为模型误差扰动矩阵, 表示n×p维的实空间, I3×3表示3×3的单位阵; 为测量向量, 表示m维的实空间; 为测量输出向量, 系统噪声w(t)是均值为零,方差为Q的高斯白噪声;测量噪声v(t)是均值为零、方差为R的高斯白噪声;
其中, 为tk时刻的系统状态估计; 为由 根据测量方程得到的tk时刻的测量估计;R为观测噪声方差阵;yk+1为tk+1时刻的实际测量向量;T为滤波周期; 为灵敏度矩阵;Λ(T)为对角矩阵; 为列向量;W为加权矩阵;
对角阵Λ(T)为:
列向量其中分别与第i颗卫星的伪距、伪距率测量对应的元素为:
最终根据公式(6)得到[tk,tk+1]区间的机动加速度估计;
步骤四:利用[tk,tk+1]区间的机动加速度估计修正系统模型,并进行时间更新,预测载体导航定位参数;
将得到的机动加速度作为系统输入代入系统模型公式(1)中,并对状态进行时间更新,更新过程如公式(10)、(11)所示,状态及其误差协方差阵的一步预测为:
Pk+1/k=(Ak-BkMkHk)Pk(Ak-BkMkHk)T+Q′k (11)
其中 为tk+1时刻的状态一步预测;Pk+1/k为tk+1时刻的状态误差协方差阵一步预测; 为tk时刻模型向量的估计值;f(x,tk)为tk时刻的模型向量; 为[tk,tk+1]区间的模型误差估计;gd为误差扰动矩阵;Ak、Bk和Mk为实现表达而出现的变量,它们是 I为单位阵,d为模型误差;Hk为基于状态 得到的测量模型的雅可比矩阵, Q′k为对称正定矩阵;Pk为tk时刻状态误差方差阵;
步骤五:进行测量更新,利用测量残差修正预测状态,实现载体导航定位参数的估计;
如式(12)、(13)、(14)所示,tk+1时刻的增益矩阵Kk+1为:
Hk+1/k为基于预测状态 得到的测量模型的雅可比矩阵, Rk+1为tk+1时刻观测的噪声方差阵;基于测量残差 对tk+1时刻的预测状态 进行修正得到tk+1时刻的状态估计 及其误差协方差阵Pk+1为:
2.根据权利要求1所述的一种卫星导航接收机的动态定位解算方法,其特征在于:当采用几何精度因子最小的4颗卫星的测量用于动态定位时,可见卫星颗数num=4。
4.根据权利要求1所述的一种卫星导航接收机的动态定位解算方法,其特征在于:Qk′=λQk,对称正定矩阵Qk为系统噪声方差阵,λ=0.8~1 ,λ为修正系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012102256259A CN102749633B (zh) | 2012-06-29 | 2012-06-29 | 一种卫星导航接收机的动态定位解算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012102256259A CN102749633B (zh) | 2012-06-29 | 2012-06-29 | 一种卫星导航接收机的动态定位解算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102749633A CN102749633A (zh) | 2012-10-24 |
CN102749633B true CN102749633B (zh) | 2013-11-27 |
Family
ID=47029977
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012102256259A Expired - Fee Related CN102749633B (zh) | 2012-06-29 | 2012-06-29 | 一种卫星导航接收机的动态定位解算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102749633B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103389502B (zh) * | 2013-07-29 | 2015-05-27 | 中国人民解放军国防科学技术大学 | 基于多个地面基站高精度确定载体加速度的方法 |
CN103529482B (zh) * | 2013-10-25 | 2016-05-11 | 中国人民解放军国防科学技术大学 | 一种高精度确定载体动态加速度的方法 |
CN103630912B (zh) * | 2013-11-26 | 2016-04-13 | 中国科学院嘉兴微电子与系统工程中心 | 一种卫星接收机静止的检测方法 |
CN104678382B (zh) * | 2013-11-29 | 2017-07-14 | 中国航天科工集团第三研究院第八三五七研究所 | 一种适用于扩频体制通信测控系统下的远程高精度测距方法 |
CN104035110B (zh) * | 2014-06-30 | 2016-08-17 | 北京理工大学 | 应用于多模卫星导航系统中的快速卡尔曼滤波定位方法 |
GB2566748B (en) * | 2017-09-26 | 2022-08-17 | Focal Point Positioning Ltd | A method and system for calibrating a system parameter |
CN106338753B (zh) * | 2016-09-22 | 2019-03-12 | 北京航空航天大学 | 一种基于地面站/星间链路/gnss联合测量的地球同步轨道星座定轨方法 |
CN109903545A (zh) * | 2019-03-27 | 2019-06-18 | 湖北医药学院 | 一种车联网数据传输的方法及系统 |
CN111274715A (zh) * | 2020-03-16 | 2020-06-12 | 深圳见炬科技有限公司 | 一种基于多维时空数据的高精度动态轨迹约束方法 |
CN114397681B (zh) * | 2021-11-26 | 2023-02-28 | 中国人民解放军军事科学院国防科技创新研究院 | 基于鲁棒预测变结构滤波的gnss接收机载波跟踪方法 |
CN115494527B (zh) * | 2022-04-13 | 2023-10-31 | 无锡奇芯科技有限公司 | 一种基于相关系数的卫星系统故障排除的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1840592A1 (en) * | 1998-05-05 | 2007-10-03 | Snaptrack, Inc. | Method and system for using altitude information in a satellite positioning system |
CN101609140A (zh) * | 2009-07-09 | 2009-12-23 | 北京航空航天大学 | 一种兼容导航接收机定位系统及其定位方法 |
CN101819278A (zh) * | 2010-03-25 | 2010-09-01 | 北京航空航天大学 | 一种gps接收机的信号捕获方法 |
-
2012
- 2012-06-29 CN CN2012102256259A patent/CN102749633B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1840592A1 (en) * | 1998-05-05 | 2007-10-03 | Snaptrack, Inc. | Method and system for using altitude information in a satellite positioning system |
CN101609140A (zh) * | 2009-07-09 | 2009-12-23 | 北京航空航天大学 | 一种兼容导航接收机定位系统及其定位方法 |
CN101819278A (zh) * | 2010-03-25 | 2010-09-01 | 北京航空航天大学 | 一种gps接收机的信号捕获方法 |
Non-Patent Citations (2)
Title |
---|
EKF与UKF在紧耦合组合导航系统中的应用;赵思浩等;《系统工程与电子技术》;20091031;第31卷(第10期);全文 * |
赵思浩等.EKF与UKF在紧耦合组合导航系统中的应用.《系统工程与电子技术》.2009,第31卷(第10期), |
Also Published As
Publication number | Publication date |
---|---|
CN102749633A (zh) | 2012-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102749633B (zh) | 一种卫星导航接收机的动态定位解算方法 | |
Shen et al. | Dual-optimization for a MEMS-INS/GPS system during GPS outages based on the cubature Kalman filter and neural networks | |
Hu et al. | Robust unscented Kalman filtering with measurement error detection for tightly coupled INS/GNSS integration in hypersonic vehicle navigation | |
Georgy et al. | Modeling the stochastic drift of a MEMS-based gyroscope in gyro/odometer/GPS integrated navigation | |
Chang et al. | Initial alignment by attitude estimation for strapdown inertial navigation systems | |
Shen et al. | Multi-rate strong tracking square-root cubature Kalman filter for MEMS-INS/GPS/polarization compass integrated navigation system | |
CN102928858B (zh) | 基于改进扩展卡尔曼滤波的gnss单点动态定位方法 | |
CN101979277B (zh) | 卫星磁测磁控系统的全实物验证平台与工作方法 | |
CN102620748B (zh) | 捷联惯导系统晃动基座条件下杆臂效应的估计和补偿方法 | |
CN103033186A (zh) | 一种用于水下滑翔器的高精度组合导航定位方法 | |
CN106291645A (zh) | 适于高维gnss/ins深耦合的容积卡尔曼滤波方法 | |
Williamson et al. | An instrumentation system applied to formation flight | |
CN104019828A (zh) | 高动态环境下惯性导航系统杆臂效应误差在线标定方法 | |
CN103439731A (zh) | 基于无迹卡尔曼滤波的gps/ins组合导航方法 | |
Zhang et al. | A hybrid intelligent algorithm DGP-MLP for GNSS/INS integration during GNSS outages | |
CN101082493A (zh) | 一种农业机械导航的组合定位方法 | |
CN108931799A (zh) | 基于递归快速正交搜索的列车组合定位方法及系统 | |
CN104061932A (zh) | 一种利用引力矢量和梯度张量进行导航定位的方法 | |
CN102353378A (zh) | 一种矢量形式信息分配系数的自适应联邦滤波方法 | |
CN115267863A (zh) | 一种精密单点定位逐级模糊度固定方法 | |
CN103017787A (zh) | 适用于摇摆晃动基座的初始对准方法 | |
CN104374405A (zh) | 一种基于自适应中心差分卡尔曼滤波的mems捷联惯导初始对准方法 | |
Gao et al. | Modeling of multi-sensor tightly aided BDS triple-frequency precise point positioning and initial assessments | |
CN102508217B (zh) | 建立雷达测量误差标定模型的方法 | |
Liu et al. | Robust train localisation method based on advanced map matching measurement-augmented tightly-coupled GNSS/INS with error-state UKF |
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: 20131127 Termination date: 20150629 |
|
EXPY | Termination of patent right or utility model |