CN107003386A - 一种卫星导航测姿方法和装置及无人机 - Google Patents
一种卫星导航测姿方法和装置及无人机 Download PDFInfo
- Publication number
- CN107003386A CN107003386A CN201580060906.5A CN201580060906A CN107003386A CN 107003386 A CN107003386 A CN 107003386A CN 201580060906 A CN201580060906 A CN 201580060906A CN 107003386 A CN107003386 A CN 107003386A
- Authority
- CN
- China
- Prior art keywords
- antenna
- satellite
- double
- difference
- antennas
- 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 51
- 239000013598 vector Substances 0.000 claims abstract description 158
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 claims description 25
- 238000004364 calculation method Methods 0.000 claims description 25
- 238000005259 measurement Methods 0.000 claims description 21
- 238000000691 measurement method Methods 0.000 claims description 11
- 238000012360 testing method Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 5
- 230000004927 fusion Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 7
- 239000011159 matrix material Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000009434 installation Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000005433 ionosphere Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000005436 troposphere Substances 0.000 description 1
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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/14—Determining absolute distances from a plurality of spaced points of known location
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
一种卫星导航测姿方法和装置及无人机,以提高得到的姿态数据的精度,并简化测姿方法,降低成本。在本发明一些可行的实施方式中,方法包括:获得卫星信号的原始观测量、星历以及天线的位置,所述原始观测量包括载波观测值和伪距观测值;根据星历和天线的位置得出天线至卫星之间的单位矢量;基于载波观测值得出载波双差,基于星历和天线的位置得出天线到卫星的单位矢量双差;基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度;基于所述整周模糊度、天线到卫星的单位矢量双差和载波双差得出三个天线之间的基线矢量;基于天线的位置和三个天线之间的基线矢量得出姿态角。
Description
本发明涉及卫星导航技术领域,具体涉及一种卫星导航测姿方法和装置及无人机。
小型无人机等飞行器基于体积、成本和重量等考虑,普遍使用指南针获取航向信息,使用INS(Inertial Navigation System,即惯性导航系统)获取飞行器姿态信息。其中,指南针通过地球磁场的作用,达到指向的效果,但指南针测得的航向包含磁偏角,且在有磁场干扰时无法使用。而INS的精度不高,获取的姿态数据难以直接使用。
发明内容
本发明实施例提供一种卫星导航测姿方法和装置及无人机,以提高得到的姿态数据的精度,并简化测姿方法,降低成本。
本发明第一方面提供一种卫星导航测姿方法,用于测试无人机的姿态,所述无人机上设置有不在同一条直线上的三个天线,所述方法包括:
获得卫星信号的原始观测量、星历以及天线的位置,所述原始观测量包括载波观测值和伪距观测值;
根据星历和天线的位置得出天线至卫星之间的单位矢量;
基于载波观测值得出载波双差,基于星历和天线的位置得出天线到卫星的单位矢量双差;
基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度;
基于所述整周模糊度、天线到卫星的单位矢量双差和载波双差得出三个天线之间的基线矢量;
基于天线的位置和三个天线之间的基线矢量得出姿态角。
本发明第二方面提供一种卫星导航装置,用于测试无人机的姿态,所述无人机上设置有不在同一条直线上的三个天线,所述装置包括:
获取模块,用于获得卫星信号的原始观测量、星历以及天线的位置,所述原始观测量包括载波观测值和伪距观测值;
计算模块,用于根据星历和天线的位置得出天线至卫星之间的单位矢量;基于载波观测值得出载波双差,基于星历和天线的位置得出天线到卫星的单位矢量双差;基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度;基于所述整周模糊度、天线到卫星的单位矢量双差和载波双差得出三个天线之间的基线矢量;基于天线的位置和三个天线之间的基线矢量得出姿态角。
本发明第三方面提供一种无人机,包括:机身,设置在所述机身且不在同一条直线上的至少三个天线,所述至少三个天线分别与至少三个接收机连接,所述至少三个接收机与一飞控系统连接,所述飞控系统用于:
通过所述至少三个接收机和至少三个天线,获得卫星信号的原始观测量、星历以及天线的位置,所述原始观测量包括载波观测值和伪距观测值;根据星历和天线的位置得出天线至卫星之间的单位矢量;基于载波观测值得出载波双差,基于星历和天线的位置得出天线到卫星的单位矢量双差;基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度;基于所述整周模糊度、天线到卫星的单位矢量双差和载波双差得出三个天线之间的基线矢量;基于天线的位置和三个天线之间的基线矢量得出姿态角。
由上可见,本发明实施例卫星导航测姿方法采用基于卫星信号得出的载波双差等数据,求解整周模糊度和天线之间的基线矢量,根据基线矢量计算姿态角的技术方案,取得了以下技术效果:
一、采用卫星信号测试无人机姿态,相对于INS定位提高了姿态数据的精度,并且避免了指南针指向存在的磁偏角和磁场干扰等问题;
二、通过构建载波相位双差方程,采用基线距离约束的最小二乘搜索方法求解整周模糊度,一方面,适用于单频接收机,相对于现有技术中采用的双频接收机,降低了成本;另一方面,简化了计算方法,降低了计算复杂度。
为了更清楚地说明本发明实施例技术方案,下面将对实施例和现有技术描
述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1是本发明实施例提供的卫星导航测姿方法的流程图;
图2是基线矢量与卫星距离的关系示意图;
图3是无人机及天线的示意图;
图4是无人机姿态解算处理过程的示意图;
图5是本发明第二实施例提供一种卫星导航装置的示意图。
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
本发明实施例技术方案适用于各种无人机,例如无人飞行器,无人汽车等,尤其是各种小型无人机的导航测姿。所述无人飞行器可以是直升机无人机或固定翼无人机等。下面通过具体实施例,进行详细的说明。
本发明第一实施例提供一种卫星导航测姿方法,用于测试无人机的姿态。本发明实施例中,使用GPS(Global Positioning System,全球定位系统)定位获取飞行器例如无人机的姿态信息。所述无人机上可设置有不在同一条直线上的至少三个天线,且每个天线与一个接收机相连。以无人机上设置有三个天线为例,三个天线分别与三个接收机相连,三个接收机则与一个卫星导航装置连接,卫星导航装置与一飞控系统连接。其中,接收机通过天线接收卫星信号,获得原始观测量、星历以及天线的位置,并传送给卫星导航装置。卫星导航装置可包括一姿态解算处理器,用于根据获取的卫星信号的原始观测量、星历以及天线的位置,解算出无人机的姿态数据,并传送给飞控系统。飞控系统可根据得到的姿态数据控制无人机飞行。其中,卫星导航装置也可以视为飞控系统的一部分。本发明实施例方法可由所述卫星导航装置执行,或者说,由所述姿
态解算处理器执行。
请参考图1,本发明实施例提供的卫星导航测姿方法的具体过程可包括:
101、获得卫星信号的原始观测量、星历以及天线的位置,所述原始观测量包括载波观测值和伪距观测值。
卫星导航测姿是基于载波相对定位技术上发展起来的,载波相对定位是用两个天线同时接收GPS信号,再通过相对定位算法计算两个天线的基线矢量(位置矢量)。安装3个GPS天线在特定位置可以计算出两个基线矢量,从而计算出载体的姿态。载波误差为毫米级别,因此基线矢量的精度也在毫米级,在保证一定基线长度(即天线距离)的前提下,GPS测姿可以达到很高的精度。
GPS接收机可以获取卫星的伪距观测值,载波观测值,以及卫星星历。卫星位置可以通过接收到的星历计算得出,只要有4颗以上卫星的伪距,就可以计算出接收机的时间、位置和速度。因此,卫星导航装置可以从GPS接收机获取卫星信号的载波观测值和伪距观测值,星历和天线的位置,以及接收机的时间和速度等。
102、根据星历和天线的位置得出天线至卫星之间的单位矢量;基于载波观测值得出载波双差,基于星历和天线的位置得出天线到卫星的单位矢量双差。
载波观测值具有较高的测量精度,接收机稳定跟踪卫星信号后,可以得到载波周期的小数部分以及后续的载波变化值,但存在一个固定不变的载波整周数无法确定,无法直接用于单点定位。该无法确定的载波整周数称为整周模糊度。由于伪距和载波测量值包含有误差,要想达到高精度的基线矢量测量,必须用载波差分定位。载波差分定位的关键是确定整周模糊度。
伪距和载波观测值包含三类的误差,列举如下:
一、与卫星相关的误差,包括卫星星历误差、卫星钟差误差、卫星天线相位中心偏差以及相对论效应误差;
二、与信号传播有关的误差,包括电离层延迟误差、对流层延迟误差以及信号多路径误差;
三、与接收设备有关的误差,包括接收机钟差、天线相位中心偏差以及接
收机测量噪声。
单点定位时,上述第一类误差只能修正相对论效应误差,第二类误差能修复部分电离层延迟误差和部分对流层延迟误差,第三类误差靠接收机的硬件和软件设计来降低。单点定位的误差只能保证在10m以内的精度,无人机安装两个天线的距离不会超过10m,无法通过单点定位的方式确定两个天线坐标,再计算方向。
对于相距不远的两个天线而言,第一类和第二类误差有较强的空间相关性,同一时间的误差影响基本一致。测量值差分是通过测量值之间作差的方式来获取新的观测值,达到消除误差的目的。差分相对定位是利用差分观测值来计算天线之间的矢量,包括伪距差分定位和载波差分定位两种。由于伪距的测量精度不高,为米级,因此伪距差分定位的误差也是米级。要想达到高精度的基线矢量测量,必须用载波差分定位。
载波差分观测值是将相同频率的载波相位观测值依据某种方式求差所获得的新的组合观测值。一般使用的是载波相位双差观测值,即天线之间做一次差,可以消除卫星相关误差和信号传播误差,卫星之间再做一次差,可以消除接收机钟差。
103、基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度。
本发明实施例中,可基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度。优选的,采用基线距离约束的最小二乘搜索方法求解所述整周模糊度。该方法适用于短距离固定长度基线测姿,且可使用单频接收机。相对于使用双频接收机,使用单频接收机可有效降低成本。
104、基于所述整周模糊度、天线到卫星的单位矢量双差和载波双差得出三个天线之间的基线矢量。
将所述整周模糊度、天线到卫星的单位矢量双差和载波双差代入构建的载波相位双差方程进行求解,即可得出三个天线之间的基线矢量。
105、基于天线的位置和三个天线之间的基线矢量得出姿态角。
结合天线的位置和三个天线之间的基线矢量即可得出无人机的姿态角,所
说的姿态角包括航向角和俯仰角以及滚转角。
进一步的,本发明一些实施例中,所述的采用基线距离约束的最小二乘搜索方法求解所述整周模糊度可以包括:将跟踪到的卫星分为主星和冗余星,所述主星包括四颗,且四颗主星中的一颗作为参考星,所述参考星的仰角最高;构造所述参考星与其他主星的载波双差方程,确定主星载波双差方程的双差模糊度搜索范围,计算基线矢量可能解;构造所述参考星与冗余星的载波双差方程,将所述基线矢量可能解代入冗余星载波双差方程,确定冗余星的双差模糊度;根据所有卫星的双差模糊度,确定整周模糊度。
本发明一些实施例中,所述根据所有卫星的双差模糊度,确定整周模糊度可以包括:将所有卫星的双差模糊度代入中载波双差方程中,计算方程的残差平方和,在所有可能解中将最小残差平方和对应的整周模糊度取为待定解,当多个历元取得相同的待定解,确认此待定解即为固定解。
本发明一些实施例中,所述基线矢量包括:从所述三个天线中的第一天线至第二天线的第一基线矢量,以及第一天线至第三天线的第二基线矢量;所述基于天线的位置和三个天线之间的基线矢量得出姿态角包括:基于所述第一天线和第二天线的位置以及所述第一基线矢量计算得到所述无人机的航向角和俯仰角;基于所述第一天线和第三天线的位置以及所述第二基线矢量计算得到所述无人机的滚转角。
本发明一些实施例中,所述获得接收卫星信号的原始观测量之前,还包括:利用单频接收机通过所述天线接收卫星信号。
本发明一些实施例中,上述方法还包括:将得出的姿态角和惯导数据融合,得到姿态角的优化估值,用于进行无人机姿态控制。
在本发明的一些实施例中,所述无人机还包括一惯性测量装置,所述惯导数据是通过所述惯性测量装置测量。
由上可见,本发明实施例公开一种卫星导航测姿方法,采用基于卫星信号得出的载波双差等数据,求解整周模糊度和天线之间的基线矢量,根据基线矢量计算姿态角的技术方案,取得了以下技术效果:
一、采用卫星信号测试无人机姿态,相对于INS定位提高了姿态数据的精
度,并且避免了指南针指向存在的磁偏角和磁场干扰等问题;
二、通过构建载波相位双差方程,采用基线距离约束的最小二乘搜索方法求解整周模糊度,一方面,适用于单频接收机,相对于现有技术中采用的双频接收机,降低了成本;另一方面,简化了计算方法,降低了计算复杂度。
下面,对本发明实施例技术方案做进一步详细说明。
一、阐述载波相位差分方程的列取。
1、载波相位单差方程
如图2所示,设天线1至卫星k之间的单位矢量为天线2至卫星k之间的单位矢量为为天线1到卫星k的距离减去天线2到卫星k的距离,b为两天线之间的矢量,当两天线之间的距离小于40Km时,远小于天线到卫星的距离(20000Km左右),可以认为则
设天线1和天线2之间的钟差为Δt,则单差观测量模型可以表示为
(式1)
其中λk为卫星k的载波波长,fk为载波频率,λk·fk=c,c为光速。为单差模糊度。考虑b=(Δx12,Δy12,,Δz12,),式1可以改写为
(式2)
2、载波相位双差方程
基准星设为j,对于GPS系统,相同频点的频率一致,由多普勒影响造成的频率变化可以忽略,则双差观测模型可以表示为
(式3)
其中为双差模糊度。
若有多颗星则可组成矩阵形式的载波相位双差方程
y=Bb+Aa+e 式(4)
其中y为双差观测值,可基于原始观测量得到,a为整周模糊度,e为噪声,A和B可理解为系数矩阵,具体的,A为模糊度权阵。
噪声e的主要部分为观测误差,跟踪环路稳定的情况下可认为小于1mm,且为无偏白噪声,在计算时可取0。接收机对卫星载波保持连续观测时,矢量a中的各元素为整数,且不随时间变化,保持不变。解算式(4)的关键是求解整周模糊度a,a求解后即可计算基线矢量b,误差为毫米级。
二、求解所述整周模糊度。
求解载波相位双差方程的关键是确定整周模糊度。本发明根据无人机天线安装位置固定,且基线长度较短的特点,采用基线距离约束的最小二乘搜索方法求解整周模糊度,适用于短距离固定长度基线测姿,且可使用单频接收机。
载波相位双差观测方程中,仅有三个双差模糊度参数独立。求解出这三个双差模糊度即可以确定基线矢量,再反求其他未知模糊度。最小二乘搜索方法将观测卫星分为四个主星和冗余星,四个主星包括一颗参考星,其他卫星跟参考星构造双差方程。确定主星载波相位双差方程的模糊度搜索范围后,计算基
线可能解,再代入冗余星观测方程,确定冗余卫星的双差模糊度。
四颗主星的选取非常关键,参考星一般选择仰角最高的卫星,其余三颗选择符合一定仰角和载噪比要求,且GDOP(Geometric Dilution Precision,几何精度因子)值较小的卫星。
假设基线矢量从天线1指向天线2,基线长度范围为[dmin,dmax],天线2被约束在一个空心球上,球的内壳半径为dmin,外壳半径为dmax。下面分析主星双差模糊度搜索范围。四颗主星构造的载波相位双差观测方程如下:
式(5)
式中,λ为载波波长,基线矢量b=(Δx,Δy,Δz)T ,(xi,yi,zi)T为天线到卫星的单位矢量双差,Ni为双差模糊度,为以周期为单位的载波相位双差观测值,其中i=1,2,3。基线矢量和双差模糊度为未知量,其余为已知量。
式(5)的矩阵方式表示如下:
Bb+W=0 式(6)
其中,W=(w1,w2,w3)T,i=1,2,3。
由式(6)可得
b=-B-1W 式(7)
假设基线长度为d,据式(7)则有
d2=bTb=WT(BBT)-1W=WTQW 式(8)
其中,Q为3阶对称正定矩阵。对矩阵Q做乔里斯基分解,即Q=LTL,代入式(8)可得
d2=WTLTLW=(LW)T(LW) 式(9)
L是下三角矩阵,可以表示为
式(10)
定义
l1=l11w1,l2=l21w1+l22w2,l3=l31w1+l32w2+l33w3 式(11)
代入式(9)有
式(12)
由基线长度约束可得
式(13)
由上式可以推得
-dmax<l1<dmax 式(14)
式(15)
由式(14)可得
式(16)
式(16)为双差模糊度N1的搜索范围。N1选择一个具体的值后,w1和l1就确定,根据式(12)就得到N2的搜索范围为
式(17)
N1和N2都确定后,w1和w2也就确定,从而计算得到l1和l2的值。根据式(13)可推得
式(18)
若则N3的搜索范围包括两个区间
式(19)
若则N3的搜索范围为
式(20)
于是得到整周模糊度N3。
该方法在短基线时较为有效,将整周模糊度约束在基线限制的范围,适用于无人机测姿。
进一步的,得到整周模糊度N3后,计算出基线矢量b,剔除不符合基线长度要求的解,将符合长度要求的可能基线向量解代入冗余星的观测方程,代入冗余星的载波相位双差方程,即刻计算冗余星的整周模糊度;然后,根据所有卫星的整周模糊度,计算载波相位双差方程的残差平方和,在所有可能解中将最小残差平方和对应的整周模糊度取为待定解,当多个历元取得相同的待定解,可以认为此待定解即为固定解。
进一步的,所述基线矢量包括:从所述三个天线中的第一天线至第二天线的第一基线矢量,以及第一天线至第三天线的第二基线矢量;所述根据得到的基线矢量计算所述无人机的姿态数据包括:根据所述第一基线矢量计算得到所述无人机的航向角和俯仰角;根据所述第二基线矢量计算得到所述无人机的滚转角。下面举例说明:
如图3和图4所示,本方案使用3个GPS天线,安装在机头、机尾和机翼(或称为机臂),其中无人机中心在天线1和天线2的中点。天线1和天线2的方位角即无人机的偏航角,倾斜角即俯仰角。天线3至中心的倾斜角即滚转角。若两点坐标已知,可以根据公式求解方位角和倾斜角。因此,只需有天线1的位置,同时计算出天线2、天线3以及中心相对天线1的精确矢量,即可计算出精确的
姿态信息。
如图4所示,接收机接收天线信号,计算天线位置并将位置信息和原始观测量发送至姿态解算处理器。姿态解算处理器以天线1坐标为基准,计算天线1至天线2矢量,天线1至天线3矢量,最终通过公式计算出无人机的航向角、俯仰角和滚转角。该姿态信息可作为卡尔曼(Kalman)滤波等最优估计算法的量测信息,与惯导等其他姿态数据作融合,以得到姿态的最优估计值,进一步提高姿态估计精度。计算得到的飞控数据被传送至飞控系统。
需要说明的是,目前全球卫星导航系统除GPS外,还有俄罗斯的GLONASS,中国的北斗,以及欧洲的GALILEO等。除了使用GPS实施本方案外,还可以用其他卫星导航系统,或者用多个卫星导航系统实施本方案。
由上可见,本发明实施例提供了一种无人机上卫星导航测姿的实施方案,该方法提高了飞行姿态估计的准确性,解决了指南针的磁场干扰问题,可为无人机的飞行控制提供更可靠的姿态数据;本发明实施例采用基线距离约束的最小二乘搜索方法求解整周模糊度和载波相位双差方程,适用于单频和双频接收机,在使用单频接收机时,可有效降低产品成本,同时简化了计算方法,降低了计算复杂度。
为了更好的实施本发明实施例的上述方案,下面还提供用于配合实施上述方案的相关装置。
请参考图5,本发明第二实施例提供一种卫星导航装置,用于测试无人机的姿态,所述无人机上设置有不在同一条直线上的三个天线。需要说明的是,无人机上设置的天线个数也可以多于三个,超过三个的天线可作为冗余,也可以参与测姿计算处理。所述装置可包括:
获取模块510,用于获得卫星信号的原始观测量、星历以及天线的位置,所述原始观测量包括载波观测值和伪距观测值;
计算模块520,用于根据星历和天线的位置得出天线至卫星之间的单位矢量;基于载波观测值得出载波双差,基于星历和天线的位置得出天线到卫星的单位矢量双差;基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度;基于所述整周模糊度、天线到卫星的单位矢量双差和
载波双差得出三个天线之间的基线矢量;基于天线的位置和三个天线之间的基线矢量得出姿态角。
在本发明的一些实施例中,所述计算模块520包括:
模糊度计算单元,用于采用基线距离约束的最小二乘搜索方法求解所述整周模糊度。
在本发明的一些实施例中,所述模糊度计算单元具体用于:
将跟踪到的卫星分为主星和冗余星,所述主星包括四颗,且四颗主星中的一颗作为参考星,所述参考星的仰角最高;
构造所述参考星与其他主星的载波双差方程,确定主星载波双差方程的双差模糊度搜索范围,计算基线矢量可能解;
构造所述参考星与冗余星的载波双差方程,将所述基线矢量可能解代入冗余星载波双差方程,确定冗余星的双差模糊度;
根据所有卫星的双差模糊度,确定整周模糊度。
在本发明的一些实施例中,所述模糊度计算单元还用于:
将所有卫星的双差模糊度代入中载波双差方程中,计算方程的残差平方和,在所有可能解中将最小残差平方和对应的整周模糊度取为待定解,当多个历元取得相同的待定解,确认此待定解即为固定解。
在本发明的一些实施例中,所述基线矢量包括:从所述三个天线中的第一天线至第二天线的第一基线矢量,以及第一天线至第三天线的第二基线矢量;所述计算模块520包括:姿态数据计算单元,用于基于所述第一天线和第二天线的位置以及所述第一基线矢量计算得到所述无人机的航向角和俯仰角;基于所述第一天线和第三天线的位置以及所述第二基线矢量计算得到所述无人机的滚转角。
在本发明的一些实施例中,所述装置还包括:
融合模块530,用于将得出的姿态角和惯导数据融合,得到姿态角的优化估值,用于进行无人机姿态控制。
在本发明的一些实施例中,所述获取模块510,用于利用单频接收机通过所述天线接收卫星信号。
可以理解,本发明实施例的卫星导航装置的各个功能模块的功能可根据上述方法实施例中的方法具体实现,其具体实现过程可参照上述方法实施例中的相关描述,此处不再赘述。
由上可见,在本发明实施例的一些可行的实施方式中,提供了一种卫星导航装置,该装置基于卫星信号得出的载波双差等数据,求解整周模糊度和天线之间的基线矢量,根据基线矢量计算姿态角的技术方案,取得了以下技术效果:
一、采用卫星信号测试无人机姿态,相对于INS定位提高了姿态数据的精度,并且避免了指南针指向存在的磁偏角和磁场干扰等问题;
二、通过构建载波相位双差方程定位,采用基线距离约束的最小二乘搜索方法求解整周模糊度,一方面,适用于单频接收机,相对于现有技术中采用的双频接收机,降低了成本;另一方面,简化了计算方法,降低了计算复杂度。
本发明第四实施例还提供一种无人机,该无人机包括:
机身,设置在所述机身且不在同一条直线上的至少三个天线,所述至少三个天线分别与至少三个接收机连接,所述至少三个接收机与一飞控系统连接,所述飞控系统用于:
通过所述至少三个接收机和至少三个天线,获得卫星信号的原始观测量、星历以及天线的位置,所述原始观测量包括载波观测值和伪距观测值;根据星历和天线的位置得出天线至卫星之间的单位矢量;基于载波观测值得出载波双差,基于星历和天线的位置得出天线到卫星的单位矢量双差;基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度;基于所述整周模糊度、天线到卫星的单位矢量双差和载波双差得出三个天线之间的基线矢量;基于天线的位置和三个天线之间的基线矢量得出姿态角。
容易理解的是,该无人机还可以包括动力系统和电气系统等。
在本发明的一些实施例中,所述无人机为固定翼无人机,直升机无人机,或无人汽车。
在本发明的一些实施例中,如图3所示,所述机身包括中心部、在中心部两侧的机臂、机头和机尾;三个天线分别设置在所述机头和机尾及一侧机臂。
在本发明的一些实施例中,所述飞控系统还用于采用基线距离约束的最小
二乘搜索方法求解所述整周模糊度。
在本发明的一些实施例中,所述飞控系统还具体用于:
将跟踪到的卫星分为主星和冗余星,所述主星包括四颗,且四颗主星中的一颗作为参考星,所述参考星的仰角最高;
构造所述参考星与其他主星的载波双差方程,确定主星载波双差方程的双差模糊度搜索范围,计算基线矢量可能解;
构造所述参考星与冗余星的载波双差方程,将所述基线矢量可能解代入冗余星载波双差方程,确定冗余星的双差模糊度;
根据所有卫星的双差模糊度,确定整周模糊度。
在本发明的一些实施例中,所述飞控系统还用于:将所有卫星的双差模糊度代入中载波双差方程中,计算方程的残差平方和,在所有可能解中将最小残差平方和对应的整周模糊度取为待定解,当多个历元取得相同的待定解,确认此待定解即为固定解。
在本发明的一些实施例中,所述基线矢量包括:从所述三个天线中的第一天线至第二天线的第一基线矢量,以及第一天线至第三天线的第二基线矢量;所述飞控系统还用于基于所述第一天线和第二天线的位置以及所述第一基线矢量计算得到所述无人机的航向角和俯仰角;基于所述第一天线和第三天线的位置以及所述第二基线矢量计算得到所述无人机的滚转角。
在本发明的一些实施例中,所述飞控系统还用于利用单频接收机通过所述天线接收卫星信号。
在本发明的一些实施例中,所述飞控系统用于将得出的姿态角和惯导数据融合,得到姿态角的优化估值,用于进行无人机姿态控制。
在本发明的一些实施例中,所述无人机还包括一惯性测量装置,所述惯导数据是通过所述惯性测量装置测量。
由上可见,在本发明实施例的一些可行的实施方式中,提供了一种无人机,该无人机取得了以下技术效果:
一、采用卫星信号测试无人机姿态,相对于INS定位提高了姿态数据的精度,并且避免了指南针指向存在的磁偏角和磁场干扰等问题;
二、通过构建载波相位双差方程定位,采用基线距离约束的最小二乘搜索方法求解整周模糊度,一方面,适用于单频接收机,相对于现有技术中采用的双频接收机,降低了成本;另一方面,简化了计算方法,降低了计算复杂度。
本发明第五实施例还提供一种处理器,用于测试无人机的姿态,所述无人机上设置有不在同一条直线上的三个天线,该处理器可执行以下操作步骤:
获得卫星信号的原始观测量、星历以及天线的位置,所述原始观测量包括载波观测值和伪距观测值;
根据星历和天线的位置得出天线至卫星之间的单位矢量;
基于载波观测值得出载波双差,基于星历和天线的位置得出天线到卫星的单位矢量双差;
基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度;
基于所述整周模糊度、天线到卫星的单位矢量双差和载波双差得出三个天线之间的基线矢量;
基于天线的位置和三个天线之间的基线矢量得出姿态角。
在本发明的一些实施例中,所述处理器还用于:采用基线距离约束的最小二乘搜索方法求解所述整周模糊度。
在本发明的一些实施例中,所述处理器还用于:
将跟踪到的卫星分为主星和冗余星,所述主星包括四颗,且四颗主星中的一颗作为参考星,所述参考星的仰角最高;
构造所述参考星与其他主星的载波双差方程,确定主星载波双差方程的双差模糊度搜索范围,计算基线矢量可能解;
构造所述参考星与冗余星的载波双差方程,将所述基线矢量可能解代入冗余星载波双差方程,确定冗余星的双差模糊度;
根据所有卫星的双差模糊度,确定整周模糊度。
在本发明的一些实施例中,所述处理器还用于:
将所有卫星的双差模糊度代入中载波双差方程中,计算方程的残差平方和,在所有可能解中将最小残差平方和对应的整周模糊度取为待定解,当多个
历元取得相同的待定解,确认此待定解即为固定解。
在本发明的一些实施例中,所述基线矢量包括:从所述三个天线中的第一天线至第二天线的第一基线矢量,以及第一天线至第三天线的第二基线矢量;所述处理器还用于:基于所述第一天线和第二天线的位置以及所述第一基线矢量计算得到所述无人机的航向角和俯仰角;基于所述第一天线和第三天线的位置以及所述第二基线矢量计算得到所述无人机的滚转角。
所述处理器还用于:利用单频接收机通过所述天线接收卫星信号。
所述处理器还用于:将得出的姿态角和惯导数据融合,得到姿态角的优化估值,用于进行无人机姿态控制。
本发明第六实施例还提供一种计算机存储介质,该计算机存储介质可存储有程序,该程序执行时包括上述方法实施例中记载的卫星导航测姿方法的部分或全部步骤。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详细描述的部分,可以参见其它实施例的相关描述。
需要说明的是,对于前述的各方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本发明并不受所描述动作顺序的限制,因为依据本发明,某些步骤可以采用其它顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作和模块并不一定是本发明所必须的。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储介质中,存储介质可以包括:ROM、RAM、磁盘或光盘等。
以上对本发明实施例所提供的卫星导航测姿方法和装置及无人机进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (24)
- 一种卫星导航测姿方法,其特征在于,用于测试无人机的姿态,所述无人机上设置有不在同一条直线上的三个天线,所述方法包括:获得卫星信号的原始观测量、星历以及天线的位置,所述原始观测量包括载波观测值和伪距观测值;根据星历和天线的位置得出天线至卫星之间的单位矢量;基于载波观测值得出载波双差,基于星历和天线的位置得出天线到卫星的单位矢量双差;基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度;基于所述整周模糊度、天线到卫星的单位矢量双差和载波双差得出三个天线之间的基线矢量;基于天线的位置和三个天线之间的基线矢量得出姿态角。
- 根据权利要求1所述的方法,其特征在于,所述基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度包括:采用基线距离约束的最小二乘搜索方法求解所述整周模糊度。
- 根据权利要求2所述的方法,其特征在于,所述采用基线距离约束的最小二乘搜索方法求解所述整周模糊度包括:将跟踪到的卫星分为主星和冗余星,所述主星包括四颗,且四颗主星中的一颗作为参考星,所述参考星的仰角最高;构造所述参考星与其他主星的载波双差方程,确定主星载波双差方程的双差模糊度搜索范围,计算基线矢量可能解;构造所述参考星与冗余星的载波双差方程,将所述基线矢量可能解代入冗余星载波双差方程,确定冗余星的双差模糊度;根据所有卫星的双差模糊度,确定整周模糊度。
- 根据权利要求3所述的方法,其特征在于,所述根据所有卫星的双差模糊度,确定整周模糊度包括:将所有卫星的双差模糊度代入中载波双差方程中,计算方程的残差平方 和,在所有可能解中将最小残差平方和对应的整周模糊度取为待定解,当多个历元取得相同的待定解,确认此待定解即为固定解。
- 根据权利要求1至4任一所述的方法,其特征在于,所述基线矢量包括:从所述三个天线中的第一天线至第二天线的第一基线矢量,以及第一天线至第三天线的第二基线矢量;所述基于天线的位置和三个天线之间的基线矢量得出姿态角包括:基于所述第一天线和第二天线的位置以及所述第一基线矢量计算得到所述无人机的航向角和俯仰角;基于所述第一天线和第三天线的位置以及所述第二基线矢量计算得到所述无人机的滚转角。
- 根据权利要求1至4任一所述的方法,其特征在于,所述获得接收卫星信号的原始观测量之前,还包括:利用单频接收机通过所述天线接收卫星信号。
- 根据权利要求1至4任一所述的方法,其特征在于,还包括:将得出的姿态角和惯导数据融合,得到姿态角的优化估值,用于进行无人机姿态控制。
- 一种卫星导航装置,其特征在于,用于测试无人机的姿态,所述无人机上设置有不在同一条直线上的三个天线,所述装置包括:获取模块,用于获得卫星信号的原始观测量、星历以及天线的位置,所述原始观测量包括载波观测值和伪距观测值;计算模块,用于根据星历和天线的位置得出天线至卫星之间的单位矢量;基于载波观测值得出载波双差,基于星历和天线的位置得出天线到卫星的单位矢量双差;基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度;基于所述整周模糊度、天线到卫星的单位矢量双差和载波双差得出三个天线之间的基线矢量;基于天线的位置和三个天线之间的基线矢量得出姿态角。
- 根据权利要求8所述的装置,其特征在于,所述计算模块包括:模糊度计算单元,用于采用基线距离约束的最小二乘搜索方法求解所述整周模糊度。
- 根据权利要求9所述的装置,其特征在于,所述模糊度计算单元具体用于:将跟踪到的卫星分为主星和冗余星,所述主星包括四颗,且四颗主星中的一颗作为参考星,所述参考星的仰角最高;构造所述参考星与其他主星的载波双差方程,确定主星载波双差方程的双差模糊度搜索范围,计算基线矢量可能解;构造所述参考星与冗余星的载波双差方程,将所述基线矢量可能解代入冗余星载波双差方程,确定冗余星的双差模糊度;根据所有卫星的双差模糊度,确定整周模糊度。
- 根据权利要求10所述的装置,其特征在于,所述模糊度计算单元还用于:将所有卫星的双差模糊度代入中载波双差方程中,计算方程的残差平方和,在所有可能解中将最小残差平方和对应的整周模糊度取为待定解,当多个历元取得相同的待定解,确认此待定解即为固定解。
- 根据权利要求8至11任一所述的装置,其特征在于,所述基线矢量包括:从所述三个天线中的第一天线至第二天线的第一基线矢量,以及第一天线至第三天线的第二基线矢量;所述计算模块包括:姿态数据计算单元,用于基于所述第一天线和第二天线的位置以及所述第一基线矢量计算得到所述无人机的航向角和俯仰角;基于所述第一天线和第三天线的位置以及所述第二基线矢量计算得到所述无人机的滚转角。
- 根据权利要求8至11任一所述的装置,其特征在于,所述获取模块,用于利用单频接收机通过所述天线接收卫星信号。
- 根据权利要求8至11任一所述的装置,其特征在于,还包括:融合模块,用于将得出的姿态角和惯导数据融合,得到姿态角的优化估值,用于进行无人机姿态控制。
- 一种无人机,其特征在于,包括:机身,设置在所述机身且不在同一条直线上的至少三个天线,所述至少三个天线分别与至少三个接收机连接,所述至少三个接收机与一飞控系统连接,所述飞控系统用于:通过所述至少三个接收机和至少三个天线,获得卫星信号的原始观测量、星历以及天线的位置,所述原始观测量包括载波观测值和伪距观测值;根据星历和天线的位置得出天线至卫星之间的单位矢量;基于载波观测值得出载波双差,基于星历和天线的位置得出天线到卫星的单位矢量双差;基于载波双差、天线到卫星的单位矢量双差和两个天线之间的基线长度求整周模糊度;基于所述整周模糊度、天线到卫星的单位矢量双差和载波双差得出三个天线之间的基线矢量;基于天线的位置和三个天线之间的基线矢量得出姿态角。
- 根据权利要求15所述的无人机,其特征在于,所述机身包括中心部、在中心部两侧的机臂、机头和机尾;三个天线分别设置在所述机头和机尾及一侧机臂。
- 根据权利要求15所述的无人机,其特征在于,所述无人机为固定翼无人机,直升机无人机,或无人汽车。
- 根据权利要求15所述的无人机,其特征在于,所述飞控系统还用于采用基线距离约束的最小二乘搜索方法求解所述整周模糊度。
- 根据权利要求18所述的无人机,其特征在于,所述飞控系统还具体用于:将跟踪到的卫星分为主星和冗余星,所述主星包括四颗,且四颗主星中的一颗作为参考星,所述参考星的仰角最高;构造所述参考星与其他主星的载波双差方程,确定主星载波双差方程的双差模糊度搜索范围,计算基线矢量可能解;构造所述参考星与冗余星的载波双差方程,将所述基线矢量可能解代入冗余星载波双差方程,确定冗余星的双差模糊度;根据所有卫星的双差模糊度,确定整周模糊度。
- 根据权利要求19所述的无人机,其特征在于,所述飞控系统还用于:将所有卫星的双差模糊度代入中载波双差方程中,计算方程的残差平方和,在所有可能解中将最小残差平方和对应的整周模糊度取为待定解,当多个历元取得相同的待定解,确认此待定解即为固定解。
- 根据权利要求15至20任一所述的无人机,其特征在于,所述基线矢量 包括:从所述三个天线中的第一天线至第二天线的第一基线矢量,以及第一天线至第三天线的第二基线矢量;所述飞控系统还用于基于所述第一天线和第二天线的位置以及所述第一基线矢量计算得到所述无人机的航向角和俯仰角;基于所述第一天线和第三天线的位置以及所述第二基线矢量计算得到所述无人机的滚转角。
- 根据权利要求15至20任一所述的无人机,其特征在于,所述飞控系统还用于利用单频接收机通过所述天线接收卫星信号。
- 根据权利要求15至20任一所述的无人机,其特征在于,所述飞控系统用于将得出的姿态角和惯导数据融合,得到姿态角的优化估值,用于进行无人机姿态控制。
- 根据权利要求23所述的无人机,其特征在于,所述无人机还包括一惯性测量装置,所述惯导数据是通过所述惯性测量装置测量。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2015/092243 WO2017066915A1 (zh) | 2015-10-20 | 2015-10-20 | 一种卫星导航测姿方法和装置及无人机 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107003386A true CN107003386A (zh) | 2017-08-01 |
CN107003386B CN107003386B (zh) | 2019-06-28 |
Family
ID=58556617
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580060906.5A Expired - Fee Related CN107003386B (zh) | 2015-10-20 | 2015-10-20 | 一种卫星导航测姿方法和装置及无人机 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN107003386B (zh) |
WO (1) | WO2017066915A1 (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107632316A (zh) * | 2017-09-11 | 2018-01-26 | 北京数码汇博科技有限公司 | 一种利用卫星定位测量车身姿态的方法及结构 |
CN109154815A (zh) * | 2017-11-30 | 2019-01-04 | 深圳市大疆创新科技有限公司 | 最高温度点跟踪方法、装置和无人机 |
CN110412638A (zh) * | 2019-08-16 | 2019-11-05 | 中国科学院重庆绿色智能技术研究院 | 一种低成本三天线gnss rtk定位及测姿方法 |
CN110673172A (zh) * | 2019-09-17 | 2020-01-10 | 闽江学院 | 一种接收机静态相对定位精度的测试方法及终端 |
CN111295567A (zh) * | 2018-12-03 | 2020-06-16 | 深圳市大疆创新科技有限公司 | 航向的确定方法、设备、存储介质和可移动平台 |
CN112068167A (zh) * | 2020-09-14 | 2020-12-11 | 天津云遥宇航科技有限公司 | 一种星载双天线融合定位方法 |
CN112325842A (zh) * | 2020-10-30 | 2021-02-05 | 中国电子科技集团公司第五十四研究所 | 一种多天线平面投影加权测姿方法 |
CN112703425A (zh) * | 2018-08-20 | 2021-04-23 | 半球全球卫星导航系统公司 | 用于检测伪全球导航卫星系统卫星信号的系统和方法 |
CN112771411A (zh) * | 2020-04-24 | 2021-05-07 | 深圳市大疆创新科技有限公司 | 定位方法、系统及存储介质 |
CN112924999A (zh) * | 2021-01-14 | 2021-06-08 | 华南理工大学 | 一种无人机的定位方法、系统、装置及介质 |
CN114355422A (zh) * | 2021-12-14 | 2022-04-15 | 中船航海科技有限责任公司 | 一种基于姿态辅助的多天线定位系统及定位方法 |
CN115061166A (zh) * | 2022-06-16 | 2022-09-16 | 湘潭大学 | 一种载波相位重构方法、装置、电子设备及介质 |
CN116321418A (zh) * | 2023-03-02 | 2023-06-23 | 中国人民解放军国防科技大学 | 一种基于节点构型优选地集群无人机融合估计定位方法 |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107248891B (zh) * | 2017-06-13 | 2023-08-04 | 施浒立 | 一种用于移动通信天线指向监测的测向测姿装置 |
CN108230722A (zh) * | 2017-12-02 | 2018-06-29 | 山东大学 | 北斗精准时空公交站台实时信息服务融合处理工作方法及系统装置 |
CN108363084B (zh) * | 2018-01-18 | 2022-04-08 | 和芯星通科技(北京)有限公司 | 利用卫星定位的方法和装置、卫星导航接收机、存储介质 |
WO2019203677A1 (en) * | 2018-04-18 | 2019-10-24 | Limited Liability Company "Topcon Positioning Systems" | Attitude determination with multiple antennas using semi-definite programming relaxation |
CN109444924B (zh) * | 2018-12-24 | 2023-02-14 | 哈尔滨工程大学 | 一种基于短基线多参考接收机的卫星星历故障监测方法 |
CN110554373B (zh) * | 2019-08-25 | 2022-11-15 | 中国科学院国家授时中心 | 干涉时间测量与测距方法 |
CN111221013B (zh) * | 2019-11-29 | 2023-06-09 | 中国兵器装备集团自动化研究所 | 一种可变基线双天线定向系统及其使用方法 |
CN110986937B (zh) * | 2019-12-19 | 2022-05-17 | 北京三快在线科技有限公司 | 用于无人设备的导航装置、方法及无人设备 |
CN111308533B (zh) * | 2020-03-04 | 2022-01-18 | 北京建筑大学 | 建筑塔机塔顶三维动态检测与分级预警装置 |
CN111381256B (zh) * | 2020-03-10 | 2022-07-26 | 上海卫星工程研究所 | 主动遥感卫星天线相位中心偏移误差计算的方法和系统 |
CN111624584B (zh) * | 2020-03-20 | 2023-06-20 | 中国人民解放军火箭军工程大学 | 一种非协作目标激光诱偏距离测量系统及方法 |
CN111323804B (zh) * | 2020-04-22 | 2023-08-29 | 北京国泰星云科技有限公司 | 一种基于北斗系统的船舶姿态测量设备及测量方法 |
CN113064194A (zh) * | 2020-08-12 | 2021-07-02 | 中国科学院微小卫星创新研究院 | 基于互联网星座的皮纳卫星通信、导航、定姿方法 |
CN112327342A (zh) * | 2020-10-16 | 2021-02-05 | 山东省科学院海洋仪器仪表研究所 | 一种基于平滑器滤波模型的mimu辅助gnss测姿方法 |
CN112382844B (zh) * | 2020-11-13 | 2023-04-25 | 成都戎星科技有限公司 | 一种低轨卫星通信系统的天线伺服电机控制方法和系统 |
CN113064195B (zh) * | 2021-03-16 | 2023-09-01 | 西南交通大学 | 一种利用多天线几何特征的高精度低计算载体测姿方法 |
CN113253308B (zh) * | 2021-05-12 | 2022-03-11 | 中国人民解放军陆军工程大学 | 一种无人机卫星导航终端定位性能确定方法及系统 |
CN115480274A (zh) * | 2021-05-31 | 2022-12-16 | 千寻位置网络有限公司 | 一种参考星选择方法、装置及接收机 |
CN115267837B (zh) * | 2022-06-22 | 2024-05-28 | 中国科学院国家空间科学中心 | 一种星间通信测距时间同步的地面验证方法 |
CN116232439B (zh) * | 2023-03-21 | 2024-04-26 | 中国民航大学 | 基于北斗的5g atg地面基站对空覆盖姿态测量与调控方法 |
CN116030646B (zh) * | 2023-03-30 | 2023-06-16 | 南昌航天广信科技有限责任公司 | 一种路况信息播报方法、系统、计算机及存储介质 |
CN116256782B (zh) * | 2023-05-15 | 2023-09-01 | 南京信息工程大学 | 一种基于双天线gnss单差算法的多路径误差消除方法 |
CN116879927B (zh) * | 2023-09-06 | 2023-11-21 | 智慧司南(天津)科技发展有限公司 | 基于三天线共线共钟架构的船用卫星罗经艏向测定方法 |
CN117405108B (zh) * | 2023-11-08 | 2024-05-07 | 中国人民解放军63620部队 | 目标姿态测量方法、系统、电子设备和介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7411545B2 (en) * | 2006-09-29 | 2008-08-12 | Honeywell International Inc. | Carrier phase interger ambiguity resolution with multiple reference receivers |
CN101825717A (zh) * | 2010-04-16 | 2010-09-08 | 北京航空航天大学 | 一种基于载波平滑码伪距技术的动态定姿方法 |
CN101833080A (zh) * | 2009-03-12 | 2010-09-15 | 周迅 | 一种利用gps系统附加约束条件的载体姿态测量方法 |
CN102736094A (zh) * | 2012-06-19 | 2012-10-17 | 哈尔滨工程大学 | 一种基于自适应遗传算法的单频gnss整周模糊度获取方法 |
CN104407354A (zh) * | 2014-11-03 | 2015-03-11 | 中国人民解放军63961部队 | 一种北斗卫星导航系统双频定位定向仪 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5561432A (en) * | 1995-05-12 | 1996-10-01 | Trimble Navigation | Out of plane antenna vector system and method |
US8120527B2 (en) * | 2008-01-30 | 2012-02-21 | Javad Gnss, Inc. | Satellite differential positioning receiver using multiple base-rover antennas |
CN102230971B (zh) * | 2011-03-29 | 2012-10-31 | 哈尔滨工程大学 | Gps多天线测姿方法 |
CN102998690B (zh) * | 2012-11-26 | 2014-04-16 | 哈尔滨工程大学 | 一种基于gps载波双差方程的姿态角直接求解方法 |
-
2015
- 2015-10-20 WO PCT/CN2015/092243 patent/WO2017066915A1/zh active Application Filing
- 2015-10-20 CN CN201580060906.5A patent/CN107003386B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7411545B2 (en) * | 2006-09-29 | 2008-08-12 | Honeywell International Inc. | Carrier phase interger ambiguity resolution with multiple reference receivers |
CN101833080A (zh) * | 2009-03-12 | 2010-09-15 | 周迅 | 一种利用gps系统附加约束条件的载体姿态测量方法 |
CN101825717A (zh) * | 2010-04-16 | 2010-09-08 | 北京航空航天大学 | 一种基于载波平滑码伪距技术的动态定姿方法 |
CN102736094A (zh) * | 2012-06-19 | 2012-10-17 | 哈尔滨工程大学 | 一种基于自适应遗传算法的单频gnss整周模糊度获取方法 |
CN104407354A (zh) * | 2014-11-03 | 2015-03-11 | 中国人民解放军63961部队 | 一种北斗卫星导航系统双频定位定向仪 |
Non-Patent Citations (2)
Title |
---|
史琳: "GPS整周模糊度及其在姿态测量中的应用研究", 《中国优秀硕士学位论文全文数据库 基础科技辑》 * |
秦红磊 等: "GPS单频单历元定姿算法性能分析与改进方法", 《解放军理工大学学报(自然科学版)》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107632316A (zh) * | 2017-09-11 | 2018-01-26 | 北京数码汇博科技有限公司 | 一种利用卫星定位测量车身姿态的方法及结构 |
US11153494B2 (en) | 2017-11-30 | 2021-10-19 | SZ DJI Technology Co., Ltd. | Maximum temperature point tracking method, device and unmanned aerial vehicle |
CN109154815A (zh) * | 2017-11-30 | 2019-01-04 | 深圳市大疆创新科技有限公司 | 最高温度点跟踪方法、装置和无人机 |
US11798172B2 (en) | 2017-11-30 | 2023-10-24 | SZ DJI Technology Co., Ltd. | Maximum temperature point tracking method, device and unmanned aerial vehicle |
CN112703425A (zh) * | 2018-08-20 | 2021-04-23 | 半球全球卫星导航系统公司 | 用于检测伪全球导航卫星系统卫星信号的系统和方法 |
CN111295567A (zh) * | 2018-12-03 | 2020-06-16 | 深圳市大疆创新科技有限公司 | 航向的确定方法、设备、存储介质和可移动平台 |
CN110412638A (zh) * | 2019-08-16 | 2019-11-05 | 中国科学院重庆绿色智能技术研究院 | 一种低成本三天线gnss rtk定位及测姿方法 |
CN110412638B (zh) * | 2019-08-16 | 2023-03-28 | 中国科学院重庆绿色智能技术研究院 | 一种低成本三天线gnss rtk定位及测姿方法 |
CN110673172A (zh) * | 2019-09-17 | 2020-01-10 | 闽江学院 | 一种接收机静态相对定位精度的测试方法及终端 |
CN112771411A (zh) * | 2020-04-24 | 2021-05-07 | 深圳市大疆创新科技有限公司 | 定位方法、系统及存储介质 |
CN112068167A (zh) * | 2020-09-14 | 2020-12-11 | 天津云遥宇航科技有限公司 | 一种星载双天线融合定位方法 |
CN112325842A (zh) * | 2020-10-30 | 2021-02-05 | 中国电子科技集团公司第五十四研究所 | 一种多天线平面投影加权测姿方法 |
CN112924999A (zh) * | 2021-01-14 | 2021-06-08 | 华南理工大学 | 一种无人机的定位方法、系统、装置及介质 |
CN112924999B (zh) * | 2021-01-14 | 2023-08-22 | 华南理工大学 | 一种无人机的定位方法、系统、装置及介质 |
CN114355422A (zh) * | 2021-12-14 | 2022-04-15 | 中船航海科技有限责任公司 | 一种基于姿态辅助的多天线定位系统及定位方法 |
CN115061166A (zh) * | 2022-06-16 | 2022-09-16 | 湘潭大学 | 一种载波相位重构方法、装置、电子设备及介质 |
CN115061166B (zh) * | 2022-06-16 | 2024-04-05 | 湘潭大学 | 一种载波相位重构方法、装置、电子设备及介质 |
CN116321418A (zh) * | 2023-03-02 | 2023-06-23 | 中国人民解放军国防科技大学 | 一种基于节点构型优选地集群无人机融合估计定位方法 |
CN116321418B (zh) * | 2023-03-02 | 2024-01-02 | 中国人民解放军国防科技大学 | 一种基于节点构型优选的集群无人机融合估计定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107003386B (zh) | 2019-06-28 |
WO2017066915A1 (zh) | 2017-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107003386A (zh) | 一种卫星导航测姿方法和装置及无人机 | |
US11585951B1 (en) | Heading or pitch determination systems and methods with high confidence error bounds | |
US9035826B2 (en) | Satellite differential positioning receiver using multiple base-rover antennas | |
CN109477900B (zh) | 全球导航卫星系统接收器中用于模糊度解算的频率间偏差的估算 | |
CN111751853B (zh) | 一种gnss双频载波相位整周模糊度解算方法 | |
EP3109672B1 (en) | Gnss receiver with a capability to resolve ambiguities using an uncombined formulation | |
CN110531392A (zh) | 一种基于ppp算法的高精度定位方法和系统 | |
CN109613585A (zh) | 一种对基站天线超短基线gnss双天线实时测向的方法 | |
US11609346B2 (en) | GNSS-based attitude determination algorithm and triple-antenna GNSS receiver for its implementation | |
CN108051840B (zh) | 一种基于gnss的含约束ekf相对定位方法 | |
EP2044457A2 (en) | A method for increasing the reliability of position information when transitioning from a regional, wide-area, or global carrier-phase differential navigation (wadgps) to a local real-time kinematic (rtk) navigation system | |
EP1906201A1 (en) | Carrier phase interger ambiguity resolution with multiple reference receivers | |
CN105425261A (zh) | 基于GPS/Beidou2/INS的组合导航与定位方法 | |
US20210072408A1 (en) | System and method for satellite positioning | |
US8922426B1 (en) | System for geo-location | |
Willi et al. | GNSS attitude determination with non-synchronized receivers and short baselines onboard a spacecraft | |
CN110186463A (zh) | 一种基于一致性滤波的航天器集群仅测距相对导航方法 | |
Tolman et al. | Absolute precise kinematic positioning with GPS and GLONASS | |
CN113064195B (zh) | 一种利用多天线几何特征的高精度低计算载体测姿方法 | |
CN108205151B (zh) | 一种低成本gps单天线姿态测量方法 | |
WO2017066750A1 (en) | Triple difference formulation for formation flight | |
US8963764B1 (en) | Ship heading and pitch using satellite ephemerides and radar range measurement of satellite | |
CN112415540A (zh) | 多源定位数据的无人机自主飞行系统 | |
CN115184977A (zh) | 一体化组合导航装置及导航系统 | |
CN113759404B (zh) | 一种差分定位方法及设备 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190628 |
|
CF01 | Termination of patent right due to non-payment of annual fee |