CN106199670A - 一种基于蒙特卡洛采样的gnss单频单历元姿态确定方法 - Google Patents
一种基于蒙特卡洛采样的gnss单频单历元姿态确定方法 Download PDFInfo
- Publication number
- CN106199670A CN106199670A CN201610487545.9A CN201610487545A CN106199670A CN 106199670 A CN106199670 A CN 106199670A CN 201610487545 A CN201610487545 A CN 201610487545A CN 106199670 A CN106199670 A CN 106199670A
- Authority
- CN
- China
- Prior art keywords
- fuzziness
- matrix
- attitude
- gnss
- lambda
- 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
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/53—Determining attitude
- G01S19/54—Determining attitude using carrier phase measurements; using long or short baseline interferometry
-
- 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/53—Determining attitude
- G01S19/54—Determining attitude using carrier phase measurements; using long or short baseline interferometry
- G01S19/55—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种基于蒙特卡洛采样的GNSS单频单历元姿态确定方法,其步骤如下:一:准备工作:主要是给出线性化的站际‑星际双差模型;二:建立多变量的GNSS姿态模型;三:用蒙特卡洛取样方法构造模糊度的概率分布函数;四:用LAMBDA算法搜索出模糊度的候选值;五:计算模糊度最优整数解并确定姿态;通过以上步骤,用蒙特卡洛方法求得了模糊度的期望和协方差,然后将其用于求解模糊度的LAMBDA方法中进行模糊度的解算,减少了计算模糊度的计算量,摆脱了对伪距测量量的依赖,从而可以简便快捷的确定姿态,减弱了对环境条件的要求。本发明在确定姿态的过程中,可以快速固定模糊度,从而可以快速确定姿态。
Description
技术领域
本发明提供一种基于蒙特卡洛采样的GNSS的单频单历元姿态确定方法,它涉及一种基于蒙特卡洛采样的GNSS单频单历元模糊度固定和姿态估计方法,即一种对静止或者运动载体利用GNSS进行姿态确定的单频单历元载波模糊度固定和姿态估计方法,属于导航技术领域。
背景技术
全球导航卫星系统(GNSS,Global Navigation Satellite System)目前主要包括美国的全球定位系统(GPS)、俄罗斯的格洛纳斯系统(GLONASS)、中国的北斗导航系统(BeiDou)以及欧洲的伽利略系统(Galileo)。GNSS系统具有实时、高精度和全球连续覆盖等优势,可以为地面车辆、舰船、航空飞行器以及低轨卫星提供精密授时、定位和定姿等导航服务。近年来,基于GNSS的姿态确定技术在导航领域得到了广泛的关注。通过在载体平台上安装至少三个GNSS接收机,构建单差或双差观测方程,可以获得载体的完整姿态信息。GNSS姿态确定具有零漂移、低成本、功耗低等优点,对现有的定姿系统可实现很好的替代、完善或补充。
GNSS系统的基本观测量包含伪距和载波两种,二者的观测噪声分别在分米和毫米水平。高精度GNSS姿态确定主要依赖载波信号,尤其是对于中小尺寸(米级或分米级)的运动载体,比如中小型车辆、无人机、探空火箭、卫星等,必须借助载波信号。但是,载波信号存在整周模糊度问题。因此,如何快速准确的固定载波相位整数模糊度,是实现高精度GNSS姿态确定的核心和关键。根据模糊度解算过程中是否利用载体和GNSS卫星之间的相对运动信息,可以将姿态确定模糊度解算方法分为多历元法和单历元法两种。多历元法又称运动法,其基本原理是利用观测信息的时间累积获取模糊度浮点解和协方差矩阵,从而得到固定整数解。多历元法容易受到载波信号失锁和模糊度周跳的影响,不利于在复杂的导航环境中(信号遮蔽、高动态飞行等)使用。相比之下,单历元法可以不受信号失锁和周跳的影响。而且,单历元法本身的固定速度要快于多历元法。实现单历元模糊度固定,可以采取若干种途径,比如遍历搜索、多频信号组合和多接收机长-短基线配置等。对于低成本单频GNSS定姿,搜索法是最佳选择。其中,LAMBDA方法(最小二乘模糊度去相关调整,Least squaresAMBiguity Decorrelation Adjustment)是目前效率最高、应用最广的一种方法。LAMBDA方法在本质上是一种整数最小二乘搜索法,通过Z变换使不同通道的模糊度的相关程度降到最低,同时仍保留其整数特性,因此具有更高的搜索效率和准确度。然而,应用LAMBDA方法固定模糊度需要同时利用伪距和载波测量。伪距的主要作用是缩小模糊度的搜索范围,而不是直接提供姿态解算。当伪距噪声较小时(5~30cm),LAMBDA方法可以实现100%的模糊度固定成功率;当伪距噪声过大(>50cm)或者多径误差存在时,性能会随之下降。
本发明在LAMBDA方法的基础上提出一种基于蒙特卡洛采样的GNSS单频单历元模糊度固定和姿态确定的方法。该方法从模糊度概率分布的角度出发,利用事先已知的多天线几何分布信息对模糊度的搜索空间进行限定,从而摆脱对伪距信号的依赖,因此可以用于高伪距噪声(低端接收机、弱信号环境等)和高多径环境(城市街道、载体形状不规则等)的姿态确定。该方法适合于利用短基线和超短基线定姿的中小型(静止或运动)载体。
发明内容
(一)发明目的:本发明以LAMBDA固定模糊度方法为基础,针对短基线和超短基线姿态确定问题,提出了一种基于蒙特卡洛采样的GNSS单频单历元模糊度固定和姿态估计方法。
(二)技术方案
本发明一种基于蒙特卡洛采样的GNSS单频单历元姿态确定方法,其步骤如下:
步骤一:准备工作
首先,给出m+1个接收机对n+1个GNSS卫星的载波观测量给出站际-星际双差观测线性化方程如下所示:
式(1)中,Φ为m条基线向量的载波双差观测量组成的n×m阶矩阵,每一列向量表示一条基线对应的n个双差观测量;B为m条未知基线向量在GNSS参考坐标系下的坐标值组成的3×m阶矩阵;G是从卫星到接收机方向的双差单位向量在GNSS参考坐标系下的表示组成的n×3阶矩阵(由于各个副接收机之间距离很小,而卫星距离各副接收机很远,所以这里认为同一个卫星到任何一个副接收机方向的单位向量相同);λ是载波波长;Z是双差整周模糊度组成的n×m阶矩阵,每一列向量表示一条基线对应的n个双差整周模糊度;V观测噪声组成的n×m阶矩阵;符号vec(·)表示把矩阵的列向量按列序号从小到大的顺序依次从上到下重新排列成一列,组成新的列向量;Q是vec(V)的协方差矩阵;取Q为:
其中σ是载波噪声的标准差;符号表示克罗内克积;
蒙特卡洛取样方法:表示概率分布函数px(x)的随机测量,是取样点,是各点对应的权重,并且Ns是取样次数,那么概率分布函数px(x)可近似表示为
其中,δ(·)表示狄拉克δ函数;
步骤二:建立多变量的GNSS姿态模型
定义本地坐标系:原点在主接收机位置,x轴沿着第一条基线方向,y轴在第一条基线和第二条基线决定的平面内并且垂直于x轴,z轴方向由右手定则决定;m条基线向量在本地坐标系中的表示已知,则由m条基线向量组成的3×m坐标矩阵F如下所示:
其中,F中每一列向量表示一条基线向量在本地坐标系下的表示,即:在本地坐标系下,副接收机1的坐标为(f11,0,0),副接收机2的坐标为(f21,f22,0),副接收机3的坐标为(f31,f32,f33),副接收机m的坐标为(fm1,fm2,fm3),其中各个坐标值的含义见说明书附图中的图2;
令R表示从本地坐标系到GNSS参考坐标系的坐标转移矩阵,则
B=RF (5)
将其代入式(1)即可得到多变量的GNSS姿态模型:
其中,3阶标准正交方阵R和双差模糊度矩阵为未知量;将R中的参数用四元数表示,则(6)式可写为:
其中,q=(q0,q4)T,q0=(q1,q2,q3)T,并且满足qTq=1,
q为未知量;
綜上所述,所述的“建立多变量的GNSS姿态模型”,其建立的过程是将(4)式和(5)式代入(1)式得到(6)式,再将(8)式代入(6)式得到(7)式,(7)式即为建立的多变量GNSS姿态模型;
步骤三:蒙特卡洛取样方法构造模糊度的概率分布函数
根据(7)式,双差模糊度可以表示为:
观测噪声遵从均值为0协方差为Q的正态分布;vec(V)的概率分布函数为:
如果已知四元数q有上界qu和下界ql,根据先验姿态信息,q遵从均匀分布:
其中,
根据qTq=1,可得
并且,
其中,是ql的前三个分量,是qu的前三个分量。在缺少先验姿态信息的情况下,可取ql=[1,1,1]T,qu=[1,1,1]T;根据概率分布函数pq(q)和pvec(V)(v)对q和v进行Ns次取样
每次取样点对应的权重为
根据(9)式可以得出模糊度的样本点则
模糊度的期望和协方差为
这里得出的期望和方差将用于步骤四的LAMBDA算法中;
綜上所述,所述的“构造模糊度的概率分布函数”,是指先根据姿态转换矩阵R(q)中的四元数元素q和噪声的概率分布函数对q和噪声分别进行取样,由取样结果得到模糊度的样本,这样可以根据蒙特卡洛方法得到模糊度的概率分布函数;
步骤四:LAMBDA算法搜索出模糊度的候选值
由于LAMBDA算法是已有的现成方法,这里不做详细说明;通过LAMBDA算法可以搜索出Nc个模糊度的候选值
步骤五:计算模糊度最优整数解并确定姿态
由模糊度的候选值根据(9)式可以算得各候选值对应的姿态矩阵
X+表示矩阵X的伪逆。其中四元数的候选值可由下式得到
其中{Rij,i,j=1,2,3}是矩阵R的元素;则模糊度的最终整数解和姿态矩阵为
其中,在步骤四中所述的“LAMBDA算法搜索出模糊度的候选值”,其搜索的过程如下:
先不考虑模糊度的整数限制,直接用最小二乘法解出模糊度的浮点解,把上文得出的协方差矩阵构造转换矩阵Z,采用Z整数转换矩阵将模糊度和协方差矩阵进行转换,进行整周模糊度搜索,求得模糊度的整数解及其协方差矩阵,再对求得的模糊度和协方差矩阵进行Z逆变换,得到模糊度的整数解及其协方差矩阵;
通过以上步骤,用蒙特卡洛方法求得了模糊度的期望和协方差,然后将其用于求解模糊度的LAMBDA方法中进行模糊度的解算,减少了计算模糊度的计算量,摆脱了对伪距测量量的依赖,从而可以简便快捷的确定姿态,减弱了对环境条件的要求。
(三)优点
本发明提供的一种基于蒙特卡洛采样的GNSS的单频单历元姿态确定方法的优点在于:
①本发明对于基线的姿态确定和模糊度的解算不需要伪距测量量的辅助,很好地适用于在高伪距测量噪声,多路径条件下。
②本发明在确定姿态的过程中,计算模糊度方差的计算量小,从而减小总体计算量。
③本发明在确定姿态的过程中,可以快速固定模糊度,从而可以快速确定姿态。
附图说明
图1是本发明所述方法流程图。
图2是本地坐标系下各基线向量的示意图。
具体实施方式
本发明一种基于蒙特卡洛采样的GNSS的单频单历元姿态确定方法,见图1所示,其具体实施步骤如下:
步骤一:给出线性化的星际-站际双差模型
m+1个接收机对n+1个GNSS卫星的载波观测量的站际-星际双差观测线性化方程如下所示:
式(23)中,Φ为m条基线向量的载波双差观测量组成的n×m阶矩阵,每一列向量表示一条基线对应的n个双差观测量;B为m条未知基线向量在GNSS参考坐标系下的坐标值组成的3×m阶矩阵;G是从卫星到接收机方向的双差单位向量在GNSS参考坐标系下的表示组成的n×3阶矩阵(由于各个副接收机之间距离很小,而卫星距离各副接收机很远,所以这里认为同一个卫星到任何一个副接收机方向的单位向量相同);λ是载波波长;Z是双差整周模糊度组成的n×m阶矩阵,每一列向量表示一条基线对应的n个双差整周模糊度;V观测噪声组成的n×m阶矩阵;符号vec(·)表示把矩阵的列向量按列序号从小到大的顺序依次从上到下重新排列成一列,组成新的列向量;Q是vec(V)的协方差矩阵。取Q为:
其中σ是载波噪声的标准差;符号表示克罗内克积。
步骤二:建立多变量的GNSS姿态模型
定义本地坐标系:原点在主接收机位置,x轴沿着第一条基线方向,y轴在第一条基线和第二条基线决定的平面内并且垂直于x轴,z轴方向由右手定则决定。m条基线向量在本地坐标系中的表示已知,则由m条基线向量组成的3×m坐标矩阵F如下所示:
其中,F中每一列向量表示一条基线向量在本地坐标系下的表示,其中各个坐标值的含义见附图2。
令R表示从本地坐标系到GNSS参考坐标系的坐标转移矩阵,则B=RF。那么,多变量的GNSS姿态模型可由式(23)得到:
其中,3阶标准正交方阵R和双差模糊度矩阵Z为未知量。将R中的参数用四元数表示,则(26)式可写为:
其中,q=(q0,q4)T,q0=(q1,q2,q3)T,并且满足qTq=1,
q为未知量。
根据(27)式,双差模糊度可以表示为:
步骤三:蒙特卡洛取样方法构造模糊度的概率分布函数
观测噪声遵从均值为0协方差为Q的正态分布,则vec(V)的概率分布函数为:
如果已知四元数q有上界qu和下界ql,根据先验姿态信息,q遵从均匀分布:
其中,
根据qTq=1,可得
并且,
其中,是ql的前三个分量,是qu的前三个分量。在缺少先验姿态信息的情况下,可取ql=[1,1,1]T,qu=[1,1,1]T。
根据概率分布函数pq(q)和pvec(V)(v)对q和v进行Ns次取样
每次取样点对应的权重为
根据(29)式可以得出模糊度的样本点则
模糊度的期望和协方差为
这里得出的期望和方差将用于步骤四的LAMBDA算法中。
步骤四:LAMBDA算法搜索出模糊度的候选值
由于LAMBDA算法是已有的现成方法,这里不做详细说明。通过LAMBDA算法可以搜索出Nc个模糊度的候选值
步骤五:计算各模糊度候选值对应的姿态矩阵并筛选出最优解
由模糊度的候选值根据(27)式可以算得各候选值对应的姿态矩阵
符号X+表示矩阵X的伪逆。其四元数的候选值可由下式得到
其中{Rij,i,j=1,2,3}是矩阵R的元素。
模糊度的最终整数解和姿态矩阵为
Claims (2)
1.一种基于蒙特卡洛采样的GNSS单频单历元姿态确定方法,其特征在于:其步骤如下:
步骤一:准备工作
首先,给出m+1个接收机对n+1个GNSS卫星的载波观测量给出站际-星际双差观测线性化方程如下所示:
式(1)中,Φ为m条基线向量的载波双差观测量组成的n×m阶矩阵,每一列向量表示一条基线对应的n个双差观测量;B为m条未知基线向量在GNSS参考坐标系下的坐标值组成的3×m阶矩阵;G是从卫星到接收机方向的双差单位向量在GNSS参考坐标系下的表示组成的n×3阶矩阵;λ是载波波长;Z是双差整周模糊度组成的n×m阶矩阵,每一列向量表示一条基线对应的n个双差整周模糊度;V观测噪声组成的n×m阶矩阵;符号vec(·)表示把矩阵的列向量按列序号从小到大的顺序依次从上到下重新排列成一列,组成新的列向量;Q是vec(V)的协方差矩阵;取Q为:
其中σ是载波噪声的标准差;符号表示克罗内克积;
蒙特卡洛取样方法:表示概率分布函数px(x)的随机测量,是取样点,是各点对应的权重,并且Ns是取样次数,那么概率分布函数px(x)近似表示为
其中,δ(·)表示狄拉克δ函数;
步骤二:建立多变量的GNSS姿态模型
定义本地坐标系:原点在主接收机位置,x轴沿着第一条基线方向,y轴在第一条基线和第二条基线决定的平面内并且垂直于x轴,z轴方向由右手定则决定;m条基线向量在本地坐标系中的表示已知,则由m条基线向量组成的3×m坐标矩阵F如下所示:
其中,F中每一列向量表示一条基线向量在本地坐标系下的表示,即:在本地坐标系下,副接收机1的坐标为(f11,0,0),副接收机2的坐标为(f21,f22,0),副接收机3的坐标为(f31,f32,f33),副接收机m的坐标为(fm1,fm2,fm3);
令R表示从本地坐标系到GNSS参考坐标系的坐标转移矩阵,则
B=RF···············(5)
将其代入式(1)即得到多变量的GNSS姿态模型:
其中,3阶标准正交方阵R和双差模糊度矩阵为未知量;将R中的参数用四元数表示,则(6)式写为:
其中,q=(q0,q4)T,q0=(q1,q2,q3)T,并且满足qTq=1,
q为未知量;
步骤三:蒙特卡洛取样方法构造模糊度的概率分布函数
根据(7)式,双差模糊度表示为:
观测噪声遵从均值为0协方差为Q的正态分布;vec(V)的概率分布函数为:
如果已知四元数q有上界qu和下界ql,根据先验姿态信息,q遵从均匀分布:
其中,
根据qTq=1,得
并且,
其中,是ql的前三个分量,是qu的前三个分量;在缺少先验姿态信息的情况下,取ql=[1,1,1]T,qu=[1,1,1]T;根据概率分布函数pq(q)和pvec(V)(v)对q和v进行Ns次取样
每次取样点对应的权重为
根据(9)式得出模糊度的样本点则
模糊度的期望和协方差为
这里得出的期望和方差将用于步骤四的LAMBDA算法中;
步骤四:LAMBDA算法搜索出模糊度的候选值
由于LAMBDA算法是已有的现成方法,这里不做详细说明;通过LAMBDA算法能搜索出Nc个模糊度的候选值
步骤五:计算模糊度最优整数解并确定姿态
由模糊度的候选值根据(9)式算得各候选值对应的姿态矩阵
X+表示矩阵X的伪逆;其中四元数的候选值由下式得到
其中{Rij,i,j=1,2,3}是矩阵R的元素;则模糊度的最终整数解和姿态矩阵为
通过以上步骤,用蒙特卡洛方法求得了模糊度的期望和协方差,然后将其用于求解模糊度的LAMBDA方法中进行模糊度的解算,减少了计算模糊度的计算量,摆脱了对伪距测量量的依赖,从而能简便快捷的确定姿态,减弱了对环境条件的要求。
2.根据权利要求1所述的一种基于蒙特卡洛采样的GNSS单频单历元姿态确定方法,其特征在于:在步骤四中所述的“LAMBDA算法搜索出模糊度的候选值”,其搜索的过程如下:
先不考虑模糊度的整数限制,直接用最小二乘法解出模糊度的浮点解,把得出的协方差矩阵构造转换矩阵Z,采用Z整数转换矩阵将模糊度和协方差矩阵进行转换,进行整周模糊度搜索,求得模糊度的整数解及其协方差矩阵,再对求得的模糊度和协方差矩阵进行Z逆变换,得到模糊度的整数解及其协方差矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610487545.9A CN106199670B (zh) | 2016-06-28 | 2016-06-28 | 一种基于蒙特卡洛采样的gnss单频单历元姿态确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610487545.9A CN106199670B (zh) | 2016-06-28 | 2016-06-28 | 一种基于蒙特卡洛采样的gnss单频单历元姿态确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106199670A true CN106199670A (zh) | 2016-12-07 |
CN106199670B CN106199670B (zh) | 2018-12-18 |
Family
ID=57461627
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610487545.9A Active CN106199670B (zh) | 2016-06-28 | 2016-06-28 | 一种基于蒙特卡洛采样的gnss单频单历元姿态确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106199670B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109932735A (zh) * | 2019-03-25 | 2019-06-25 | 中国铁路设计集团有限公司 | 北斗短基线单频单历元解算的定位方法 |
CN110133702A (zh) * | 2019-05-13 | 2019-08-16 | 桂林电子科技大学 | 一种基于正交变换的姿态测量方法和设备 |
CN115184895A (zh) * | 2022-09-07 | 2022-10-14 | 北京海兰信数据科技股份有限公司 | 一种基于导航雷达的静态目标识别方法及系统 |
CN115877431A (zh) * | 2023-01-04 | 2023-03-31 | 中国民航大学 | 基于阵列天线无整周模糊策略的低运算量测向装置及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4954837A (en) * | 1989-07-20 | 1990-09-04 | Harris Corporation | Terrain aided passive range estimation |
CN101614802A (zh) * | 2009-07-28 | 2009-12-30 | 中国电子科技集团公司第二十八研究所 | 一种导航卫星姿态测量方法 |
CN102156478A (zh) * | 2010-12-28 | 2011-08-17 | 北京航空航天大学 | 一种基于蚁群Unscented粒子滤波算法的组合定姿方法 |
CN103149936A (zh) * | 2013-03-01 | 2013-06-12 | 国家测绘地理信息局卫星测绘应用中心 | 一种基于由dna算法优化的upf算法的组合定姿方法 |
-
2016
- 2016-06-28 CN CN201610487545.9A patent/CN106199670B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4954837A (en) * | 1989-07-20 | 1990-09-04 | Harris Corporation | Terrain aided passive range estimation |
CN101614802A (zh) * | 2009-07-28 | 2009-12-30 | 中国电子科技集团公司第二十八研究所 | 一种导航卫星姿态测量方法 |
CN102156478A (zh) * | 2010-12-28 | 2011-08-17 | 北京航空航天大学 | 一种基于蚁群Unscented粒子滤波算法的组合定姿方法 |
CN103149936A (zh) * | 2013-03-01 | 2013-06-12 | 国家测绘地理信息局卫星测绘应用中心 | 一种基于由dna算法优化的upf算法的组合定姿方法 |
Non-Patent Citations (2)
Title |
---|
FRANCESCO CAPPELLO ET AL.: "Particle filter based multi-sensor data fusion techniques for RPAS navigation and guidance", 《METROLOGY FOR AEROSPACE (METROAEROSPACE)》 * |
逢淑涛 DENG: "基于改进LAMBDA算法的GNSS载波相位姿态测量方法研究", 《第二届中国卫星导航学术年会 CSNC2011》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109932735A (zh) * | 2019-03-25 | 2019-06-25 | 中国铁路设计集团有限公司 | 北斗短基线单频单历元解算的定位方法 |
CN110133702A (zh) * | 2019-05-13 | 2019-08-16 | 桂林电子科技大学 | 一种基于正交变换的姿态测量方法和设备 |
CN110133702B (zh) * | 2019-05-13 | 2022-12-27 | 桂林电子科技大学 | 一种基于正交变换的姿态测量方法和设备 |
CN115184895A (zh) * | 2022-09-07 | 2022-10-14 | 北京海兰信数据科技股份有限公司 | 一种基于导航雷达的静态目标识别方法及系统 |
CN115877431A (zh) * | 2023-01-04 | 2023-03-31 | 中国民航大学 | 基于阵列天线无整周模糊策略的低运算量测向装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106199670B (zh) | 2018-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106842268B (zh) | 双gnss接收机载波相位双差整周模糊度浮点解向量估计方法 | |
CN106569241B (zh) | 一种基于gnss的单频高精度定位方法 | |
CN101403790B (zh) | 单频gps接收机的精密单点定位方法 | |
CN103746757B (zh) | 一种基于星载多波束天线的单星干扰源定位方法 | |
CN110531392A (zh) | 一种基于ppp算法的高精度定位方法和系统 | |
CN101825717B (zh) | 一种基于载波平滑码伪距技术的动态定姿方法 | |
CN105223598B (zh) | 一种gnss载波相位整周模糊度单历元解算方法 | |
CN101539627B (zh) | 一种电离层立体探测星载sar成像处理平台的构建方法 | |
CN102928858B (zh) | 基于改进扩展卡尔曼滤波的gnss单点动态定位方法 | |
CN101833080A (zh) | 一种利用gps系统附加约束条件的载体姿态测量方法 | |
CN107193028A (zh) | 基于GNSS的Kalman相对定位方法 | |
CN107193023B (zh) | 一种具有闭式解的高精度北斗卫星系统单点定位方法 | |
CN106405589A (zh) | 一种全球电离层格网模型的确定方法及装置 | |
CN106199670A (zh) | 一种基于蒙特卡洛采样的gnss单频单历元姿态确定方法 | |
CN103197335B (zh) | 采用改进正则化方法抑制dgps整周模糊度病态性的方法 | |
CN101614802A (zh) | 一种导航卫星姿态测量方法 | |
CN105372691A (zh) | 一种模糊度固定的长基线卫星编队gnss相对定位方法 | |
CN101710179A (zh) | 一种gnss三频动对动定位方法 | |
CN104597471A (zh) | 面向时钟同步多天线gnss接收机的定向测姿方法 | |
CN108181630B (zh) | 一种北斗双天线旋转快速定向方法 | |
CN102998690A (zh) | 一种基于gps载波双差方程的姿态角直接求解方法 | |
CN102749637A (zh) | 一种车载gps精确定位的实现方法 | |
CN103033822B (zh) | 移动信息确定装置、方法以及接收机 | |
CN104537202B (zh) | 基于卫星编队协作的空间天线阵列合成方法 | |
CN102435194A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |