CN105510935A - 一种用于导航的装置 - Google Patents
一种用于导航的装置 Download PDFInfo
- Publication number
- CN105510935A CN105510935A CN201510937376.XA CN201510937376A CN105510935A CN 105510935 A CN105510935 A CN 105510935A CN 201510937376 A CN201510937376 A CN 201510937376A CN 105510935 A CN105510935 A CN 105510935A
- Authority
- CN
- China
- Prior art keywords
- satellite
- module
- navigation
- calculate
- observed quantity
- 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.)
- Pending
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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
-
- 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/42—Determining position
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种用于导航的装置,包括:捕获跟踪处理模块,报文提取模块,观测量提取模块,定位解算模块,输出模块。通过捕获跟踪处理模块对卫星信号进行持续跟踪得到卫星的导航电文比特数据;通过报文提取模块,观测量提取模块,获得卫星导航电文和观测量,然后通过定位解算模块进行定位解算,并通过输出模块根据海用电子设备制定的标准格式NMEA格式进行输出,可以自动持续的获得卫星的导航电文比特数据进行定位,进而避免出现导航不及时的情况发生,具有及时准确的技术效果。
Description
技术领域
本申请涉及卫星通信技术领域,尤其涉及一种用于导航的装置。
背景技术
随着科学技术的不断发展,卫星通信技术也得到了飞速的发展,电子产品的种类也越来越多,人们也享受到了科技发展带来的各种便利。现在人们可以通过各种类型的电子设备,享受随着科技发展带来的舒适生活。
例如导航仪便是现代科技的产物,导航仪多用于汽车上,用于定位、导航和娱乐,随着汽车的普及和道路的建设,城际间的经济往来更加频繁,车载导航仪显得很重要,准确定位、导航、娱乐功能集于一身的导航更能满足车主的需求,成为车上的基本装备。
但是,目前的导航仪的导航并不精确,经常会出现导航路线错误或者导航不及时的状况,例如在导航的过程中出现导航延迟,不能及时获得后续行进路线,使车主不知所措;或者将车主错误引导至海里,增加车主危险。
发明内容
为解决上述技术问题,本发明提供了一种用于导航的装置,包括:
捕获跟踪处理模块,用于对卫星信号进行持续跟踪,并得到卫星的导航电文比特数据;
报文提取模块,用于接收所述导航电文比特数据,并从所述导航电文比特数据中提取出卫星导航电文;
观测量提取模块,用于从所述导航电文比特数据中提取出观测量;
定位解算模块,用于根据所述观测量和所述卫星导航电文进行联合定位,获得导航结果;
输出模块,用于根据海用电子设备制定的标准格式NMEA格式输出所述导航结果。
优选的,所述卫星导航电文包括卫星的星历、卫星工作状态、时间***、卫星历书。
优选的,所述观测量提取模块,具体用于根据基带处理数据得到被跟踪卫星的伪距观测量、或多普勒观测量、或载波相位观测量。
优选的,所述伪距观测量是利用卫星信号的传输时间乘以光速得到信号的传输距离;所述卫星信号的传输时间指卫星信号发送时间tsv减去接收机本地时间tu。
优选的,所述多普勒观测量是在得到载波NCO的频率值后,将所述载波NCO的频率值与射频前端的标称中频值fc相减获得。
优选的,所述定位解算模块具体包括:
计算模块,用于根据所述卫星的星历计算出卫星的位置;
定位解算模块,用于利用所述卫星的位置和所述观测量对用户位置进行定位解算。
优选的,所述计算模块具体用于:
第一步:计算卫星运动的平均角速度n;
第二步:计算观测瞬间卫星的平近点角Ms;
Ms=M0+n(t-toe);式中,M0为参考时刻toe时的平近点角,由广播星历给出;
第三步:计算偏近点角:Es=Ms+essinEs;
第四步:计算近点角fs:
第五步:计算升交距角式中,ω为近地点角距,由广播星历给出;
第六步:计算摄动改正项;
第七步:对升交距角卫星矢径r、轨道倾角i进行摄动改正;
第八步:计算卫星在轨道面坐标系中的位置;
第九步:计算观测瞬间升交点的经度L;
第十步:计算卫星在瞬时地球坐标系中的位置;
第十一步:计算卫星在协议地球坐标系中的位置。
优选的,所述平均角速度n通过以下步骤获得:
首先根据广播星历中给出的参数计算出参考时刻toe的平均角速度n0:式中,GM为万有引力常数G与地球总质量M的乘积,GM=3.9860047×1014m3/s2,然后根据广播星历中给定的摄动参数△n计算观测时刻卫星的平均角速度n:n=n0+△n。
优选的,所述计算摄动改正项包括以下步骤:
通过广播星历中给出的Cuc、Cus、Crc、Crs、Cic、Cis6个摄动参数,求出由于地球引力场位函数的二阶带谐系数项而引起的升交距角μ的摄动改正项δμ、卫星矢径r的摄动改正项δr,和卫星轨道倾角i的摄动改正项δi,计算公式如下:
优选的,所述对升交距角卫星矢径r、轨道倾角i进行摄动改正的步骤为:
r=as(1-escosEs)+δr
式中:as为卫星轨道的长半径,i0为toe时刻的轨道倾角,由广播星历中的开普勒参数给出;为i的变化率,由广播星历中的摄动参数给出。
通过本发明的一个或者多个技术方案,本发明具有以下有益效果或者优点:
本发明提供了一种用于导航的装置,包括:捕获跟踪处理模块,报文提取模块,观测量提取模块,定位解算模块,输出模块。通过捕获跟踪处理模块对卫星信号进行持续跟踪得到卫星的导航电文比特数据;通过报文提取模块,观测量提取模块,获得卫星导航电文和观测量,然后通过定位解算模块进行定位解算,并通过输出模块根据海用电子设备制定的标准格式NMEA格式进行输出,可以自动持续的获得卫星的导航电文比特数据进行定位,进而避免出现导航不及时的情况发生,具有及时准确的技术效果。
附图说明
图1为本发明实施例中用于导航的装置的结构示意图。
具体实施方式
为了使本申请所属技术领域中的技术人员更清楚地理解本申请,下面结合附图,通过具体实施例对本申请技术方案作详细描述。
本发明描述的一种用于导航的装置,从原理上主要分为四部分:捕获跟踪处理部分、报文和观测量提取部分、定位解算及***输出。
基于上述原理,用于导航的装置的结构如图1所示。用于导航的装置主要包括:捕获跟踪处理模块11,报文提取模块12,观测量提取模块13,定位解算模块14,输出模块15。
捕获跟踪处理模块11,用于对卫星信号进行持续跟踪,并得到卫星的导航电文比特数据;
报文提取模块12,用于接收所述导航电文比特数据,并从所述导航电文比特数据中提取出卫星导航电文;
观测量提取模块13,用于从所述导航电文比特数据中提取出观测量。
定位解算模块14,用于根据所述观测量和所述卫星导航电文进行联合定位,获得导航结果;
输出模块15,用于根据海用电子设备制定的标准格式NMEA格式输出所述导航结果。
下面分别介绍各个模块的具体功能。
在捕获跟踪处理模块11中,捕获跟踪数据处理包括转换通道状态、调整跟踪载波环路和码环路状态、维持正确跟踪状态、得出正确多普勒、码相位以及载波相位等信息,还会进行比特同步和卫星的导航电文比特数据。
在测量时刻,捕获跟踪数据处理模块会把得到的观测量值进行封装,通过***接口传递给下一个单元(报文提取模块12,观测量提取模块13),进行报文解算和观测量提取。
报文提取模块12用于接收所述导航电文比特数据,并从所述导航电文比特数据中提取出卫星星历、历书、UTC参数、电离层参数等导航电文。而观测量提取模块13,则主要求得卫星的伪距、多普勒、载波相位等观测量。
下面具体介绍报文提取模块12和观测量提取模块13。
报文提取模块12。
报文提取模块12,主要根据GPS和BD2(一种定位***)接口控制文件对比特流信息进行解码,获取卫星导航电文。卫星导航电文为用户提供有关卫星的星历、卫星工作状态、时间***、卫星历书等数据,是用户进行定位和导航所必须的基础数据。针对不同卫星定位***导航电文的结构,需要采用不同的解码方式,以正确高效地获取导航电文。
GPS的L1波段上民码导航电文的传输速率为50bit/s,以“帧”为单位向外发送。每帧的长度为1500bit,播发完一个帧需30s。每帧包括5个子帧,每个子帧均包含300bit,播发时间为6s。每个子帧又可分为10个字,每个字都由30bit组成。导航电文第一、二、三子帧为卫星的星历参数,每30s重复一次,其内容每隔2h更新一次。第四、五子帧为卫星的历书数据,每30s翻转一页,其内容包含有25个页面,所以发送完一套完整的导航电文需要12.5min,然后再重复。
GPS导航电文比特流信息中包含前导字符、子帧的Z计数,以及卫星的星历、历书数据等内容。对比特流进行同步搜索,直至找到前导字符“10001011b”。然后按照校验算法对信息进行校验,若未通过校验则丢弃该数据;若通过校验,则当前数据有效,可以根据导航电文格式,解析出电离层参数、星历以及历书数据等等。
BD2导航电文的提取与GPS类似,也包含解帧结构和信息解析。
观测量提取模块13。
观测量提取模块13,根据基带处理数据得到被跟踪卫星的伪距、多普勒、载波相位等观测量。
伪距观测量是利用卫星信号的传输时间乘以光速得到信号的传输距离。这里的传输时间指卫星信号发送时间tsv减去接收机本地时间tu。本地时间tu直接由接收机本地时钟提供,而卫星信号发送时间tsv需要通过子帧或字串的计数,并由码跟踪环提供本地伪码相位以后才能准确计算得到。
多普勒观测量的提取是在得到载波NCO的频率值后,将其与射频前端的标称中频值fc相比,其差就是多普勒频移观测量。
定位解算模块14。
定位解算模块14,根据观测量信息和得到的导航电文进行联合定位,并利用滤波算法使定位结果更稳定可靠,满足高灵敏度需求。
定位解算模块14具体包括:
计算模块,用于根据所述卫星的星历计算出卫星的位置;
定位解算模块14,用于利用所述卫星的位置和所述观测量对用户位置进行定位解算。
下面具体介绍计算模块的具体实施原理,其解算卫星的距离是分步骤进行的,另外,可以解算GPS卫星的距离或者BD2(一种定位***)卫星的距离。
1)GPS卫星空间位置的计算
第一步:计算卫星运动的平均角速度n。
首先,根据广播星历中给出的参数计算出参考时刻toe的平均角速度n0:式中,GM为万有引力常数G与地球总质量M的乘积,其值为GM=3.9860047×1014m3/s2。然后根据广播星历中给定的摄动参数△n计算观测时刻卫星的平均角速度n:n=n0+△n。
第二步:计算观测瞬间卫星的平近点角Ms。
由于卫星的运行周期为12小时左右,采用卫星过近地点时刻t0来计算平近点角M时,外推间隔最大可达6小时。而广播星历每2小时更新一次,将参考时刻设在中央时刻时,外推间隔小于等于1小时。所以用toe来取代卫星过近地点时刻t0后,外推间隔将大大减小,用较简单的模型也能获得精度较高的结果。
Ms=M0+n(t-toe)
式中,M0为参考时刻toe时的平近点角,由广播星历给出。
第三步:计算偏近点角。
用弧度表示的开普勒方程为:Es=Ms+essinEs。
第四步:计算真近点角。
根据开普勒轨道方程,可得近点角fs与偏近点角Es之间的关系:
式中,es为卫星轨道的偏心率,由广播星历给出。由此可得真近点角计算常用公式:
第五步:计算升交距角
利用公式计算升交距角,式中,ω为近地点角距,由广播星历给出。
第六步:计算摄动改正项。
广播星历中给出了Cuc、Cus、Crc、Crs、Cic、Cis6个摄动参数,据此可求出由于地球引力场位函数的二阶带谐系数项而引起的升交距角的摄动改正项δμ、卫星矢径r的摄动改正项δr,和卫星轨道倾角i的摄动改正项δi。计算公式如下:
第七步:对升交距角卫星矢径r、轨道倾角i进行摄动改正。
r=as(1-escosEs)+δr,
式中:as为卫星轨道的长半径,i0为toe时刻的轨道倾角,由广播星历中的开普勒参数给出。为i的变化率,由广播星历中的摄动参数给出。
第八步:计算卫星在轨道面坐标系中的位置。
在轨道平面直角坐标系中(坐标原点位于地心,X轴指向升交点)卫星的平面直角坐标为:
xp=rcosμ
yp=rsinμ
第九步:计算观测瞬间升交点的经度L。
若参考时刻toe时升交点的赤经为升交点对时间的变化率为△Ω,那么观测瞬间t的升交点赤经Ω应为:
式中,△Ω可从广播星历的摄动参数中给出。设本周开始时刻(星期日0时)格林尼治恒星时为GASTweek,则观测瞬间的格林尼治恒星时为:
GAST=GASTweek+ωet,
式中ωe为地球自转角速度,其值为ωe=7.292115×l0-5rad/s,t为本周内的时间(s)。这样就可求得观测瞬间升交点的经度值为:
令 则有:
L=Ω0+△Ω(t-toe)-ωet=Ω0+(△Ω-ωe)t-△Ωtoe。
第十步:计算卫星在瞬时地球坐标系中的位置。
已知升交点的大地经度L以及轨道平面的倾角i后,就可通过两次旋转求得卫星在地固坐标系中的位置:
第十一步:计算卫星在协议地球坐标系中的位置。
观测瞬间卫星在协议地球坐标系中的位置
BD2卫星广播星历的格式与GPS卫星的广播星历基本类似,因此其卫星空间位置的计算也与上节类似,这里不再赘述。
定位解算模块14。
定位解算模块14,用到卫星位置和伪距观测量,通过求解方程来计算接收机当前时刻的位置。如果已知卫星i的位置为ri(xi,yi,zi),该卫星的伪距为pi,设接收机位置为ru(xu,yu,zu),接收机钟差为bu。对于双模定位***而言,***时差△T也作为一个未知量,则可得到测量方程:
对于双模定位***,已知至少五颗卫星的位置和伪距,代入上式并做线性化,即可联合求解,求出ru(xu,yu,zu)和接收机钟差的值。
基于上述定位基本原理,在软件程序中使用最小二乘和卡尔曼滤波方法进行接收机位置和速度的求解。然后,对定位结果进行检测。如果通过该检测,则输出此次的定位结果,定位标志置1,否则不输出本次定位结果,定位标志置0。
最小二乘法,是定位解算的基本方法,它能在含有误差与噪声的各个测量值之间求得最优解,使所有测量值的偏差平方和最小。
在满足单***定位模式下至少四颗星,双***定位模式下至少五颗星的情况下,根据GPS和BD2接口控制文件提供的方法,利用星历参数计算出参与定位卫星的位置和速度。由于伪距中包含卫星钟差和电离层、对流层延迟误差,定位解算前先利用星历参数计算出卫星钟差;对于GPS和BD2卫星,均可根据电离层延迟参数计算出电离层延迟;并根据相关模型对流层延迟进行估算。之后利用这些校正量校正伪距。将测量误差较大的伪距剔除后,使用牛顿迭代法进行求解,直至迭代结束。
卡尔曼滤波算法,是一种线性最小方差滤波方法,它采用状态方程和量测方程来描述随机线性***,按照估计状态误差的方差最小的准则,从被测量噪声污染了的量测值中实时地估计出***每一时刻的各个状态。它包括预测和修正两个过程。卡尔曼滤波器在计算方法上采用递推形式,易于在计算机上实现。
设离散化后的***状态方程和量测方程分别为:
式中:
Xk表示k时刻(n×1)维***状态矢量;Xk-1则是k-1时刻的***状态矢量;Φk|k-1为k-1时刻***的转移矩阵;Wk-1是k-1时刻(n×1)维***噪声矢量;Γk-1是***噪声驱动矩阵;Zk表示k时刻的(m×1)维测量矢量;Hk为k时刻的(m×n)维观测矩阵;Vk是k时刻(m×1)维量测噪声向量。
同时,Wk-1和Vk是互不相关的零均值白噪声序列,即:
E[Wk-1]=0,
E[Vk]=0,
其中,Qk为***噪声方差矩阵,是非负定的,Rk为测量噪声方差矩阵,是正定阵;δkj是Kroneker函数:
如果离散***满足上述要求,离散***的状态最优估计可以采用如下步骤求得:
状态一步预测:
状态估计:
滤波增益:
一步预测均方误差:
估计均方误差:
输出模块15,用于根据海用电子设备制定的标准格式NMEA格式输出所述导航结果。
NMEA格式:
1)GGA-PsitionData。
功能描述:输出语句。描述定位数据。本语句包含与接收机定位、测时相关的数据。如果只将BD、GPS卫星用于位置解算,传送标识符为BD、GP,如果使用了双***卫星取得位置解算,传送标识符用GN。
格式:
$--GGA,hhmmss.ss,llll.ll,a,yyyyy.yy,a,x,xx,x.x,x.x,U,x.x,U,xxxx,x.x,x.x*hh<CR><LF>。
表1
编号 | 含义 | 取值范围 | 单位 | 备注 |
1 | 定位时刻(UTC时间) | |||
2 | 纬度 | |||
3 | 纬度方向 | N/S | N-北纬,S-南纬 | |
4 | 经度 | |||
5 | 经度方向 | E/W | E-东经,W-西经 | |
6 | 状态指示 | 0-8 | 注1 | |
7 | 参与定位的卫星数 | |||
8 | HDOP值 | |||
9 | 天线大地高 | |||
10 | 天线大地高单位 | 米 | ||
11 | 高程异常(CGS2000) | |||
12 | 高程异常单位 | 米 | ||
13 | 差分数据龄期 | |||
14 | 差分站台ID号 | |||
15 | VDOP值 |
注1:状态指示(该数据字段不能为空)
(1)当该语句标识符为GP时,状态指示:0-定位模式不可用或无效;1-GPSSPS模式,定位有效;2-差分GPSSPS模式,定位有效;3-GPSPPS模式,定位有效;4-实时动态(RTK),***处于RTK模式中,有固定的整周数;5-浮动的RTK,***处于RTK模式中,整周数是浮动的;6-估算模式(航位推算);7-手动输入模式;8-模拟器模式。
(2)当该语句标识符为BD时,状态指示:0-定位不可用或无效;1-无差分定位,定位有效;2-差分定位,定位有效;3-双频定位,定位有效。
(3)当该语句标识符为GN时,状态指示:0-定位不可用或无效;1-兼容定位,定位有效。
(4)无定位结果时,定位信息字段为空。
2)RMC-MinimumSpecificData。
功能描述:最简导航传输数据。格式:
$--RMC,hhmmss.ss,A,llll.ll,a,yyyyy.yy,a,x.x,x.x,ddmmyy,x.x,a,a*hh<CR><LF>。
表2
3)GSV-VisibleSatellitesData。
功能描述:输出语句。本语句包含可视的卫星数、卫星标识号、仰角、方位角及信噪比(C/N0)值。每次传送最多为4颗卫星,传送的语句总数和传送的语句号在前两个字段中显示。如果可以看到多颗BD、GPS卫星,分别使用GSV语句,用标识符BD标识看到的BD卫星,用GP表示看到GPS卫星。GN标识符不与本语句一起使用。
格式:$--GSV,x,x,xx,xx,xx,xxx,x.x,……*hh<CR><LF>。
表3
编号 | 含义 | 取值范围 | 单位 | 备注 |
1 | GSV语句总数 | |||
2 | 当前GSV语句序号 | |||
3 | 视野内卫星数 | |||
4 | 卫星号 | |||
5 | 卫星仰角 | 度 | ||
6 | 卫星方位角 | 度 | ||
7 | 信噪比 | dB-Hz | ||
…… | 重复4~7字段 | 其它卫星信息 |
4)VTG-ForwardDirectionandGroundSpeedData。
功能描述:本语句包含相对于地面的实际航向和速度。如果只将BD、GPS卫星用于位置解算,传送标识符为BD、GP,如果使用了双***卫星取得位置解算,传送标识符用GN。
格式:$GPVTG,xxx.x,T,xxx.x,M,xxx.x,N,xxx.x,K*hh<CR><LF>。
表4
编号 | 含义 | 取值范围 | 单位 | 备注 |
1,2 | 对地航向 | 度 | 真北 | |
3,4 | 对地航向 | 度 | 磁北 | |
5,6 | 对地速度 | 节 | 10 --> | |
7,8 | 对地速度 | 千米/小时 | ||
5 | 卫星仰角 | 度 | ||
6 | 卫星方位角 | 度 | ||
7 | 信噪比 | dB-Hz | ||
…… | 重复4~7字段 | 其它卫星信息 |
5)GLL-GeographicTimeandPositionData。
功能描述:输出语句。大地坐标定位信息,载体的纬度、经度、定位时间与状态。如果只将BD、GPS卫星用于位置解算,传送标识符为BD、GP,如果使用了双***卫星取得位置解算,传送标识符用GN。
格式:$--GLL,llll.ll,a,yyyyy.yy,a,hhmmss.ss,A,x*hh<CR><LF>。
表5
注1:模式指示:0-自动模式;1-差分模式;2-估算(航位推算)模式;3-手动输入模式;4-模拟器模式。
6)GSA-DOPandVisibleSatellitesData。
功能描述:输出语句。本语句包含用户设备工作模式、GGA语句报告的导航解算中用到的卫星以及精度因子(DOP)值。当只用BD、GPS卫星***解算位置时,分别用标识符BD表示传送BD卫星精度因子和有效卫星号,用GP表示传送GPS卫星精度因子和有效卫星号。当综合运用BD、GPS以获得位置解算时,会产生多条GSA语句,每一条GSA语句分别用BD、GP作为标识符,以区分综合解算中用到的不同卫星***的卫星,且每条语句都有用于位置解算的组合卫星***的PDOP、HDOP、VDOP和TDOP。GN标识符不与本语句一起使用。
格式:$--GSA,a,x,xx,…,xx,x.x,x.x,x.x*hh<CR><LF>。
表6
注1:模式指示:
M-手动,强制用于2D或3D模式;
A-自动,允许2D/3D自动变换。
注2:选用模式:
1-定位不可用或无效;
2-2D;
3-3D。
7)ZDA-DataandTime
功能描述:输出语句。大地坐标定位信息,载体的纬度、经度、定位时间与状态。如果只将BD、GPS卫星用于位置解算,传送标识符为BD、GP,如果使用了双***卫星取得位置解算,传送标识符用GN。
格式:$--ZDA,hhmmss.ss,dd,mm,yyyy,xx,yy*CC。
表7
编号 | 含义 | 取值范围 | 单位 | 备注 |
1 | 时间(UTC) | 时分秒 | 注1 | |
2 | 时间(UTC) | 日月年 | 注2 | |
3 | 本地时区小时 | -13~13 | 小时 | |
4 | 本地时区分钟 | 0~59 | 分钟 |
通过本发明的一个或者多个实施例,本发明具有以下有益效果或者优点:
本发明提供了一种用于导航的装置,包括:捕获跟踪处理模块,报文提取模块,观测量提取模块,定位解算模块,输出模块。通过捕获跟踪处理模块对卫星信号进行持续跟踪得到卫星的导航电文比特数据;通过报文提取模块,观测量提取模块,获得卫星导航电文和观测量,然后通过定位解算模块进行定位解算,并通过输出模块根据海用电子设备制定的标准格式NMEA格式进行输出,可以自动持续的获得卫星的导航电文比特数据进行定位,进而避免出现导航不及时的情况发生,具有及时准确的技术效果。
尽管已描述了本申请的优选实施例,但本领域内的普通技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本申请范围的所有变更和修改。
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包含这些改动和变型在内。
Claims (10)
1.一种用于导航的装置,其特征在于,包括:
捕获跟踪处理模块,用于对卫星信号进行持续跟踪,并得到卫星的导航电文比特数据;
报文提取模块,用于接收所述导航电文比特数据,并从所述导航电文比特数据中提取出卫星导航电文;
观测量提取模块,用于从所述导航电文比特数据中提取出观测量;
定位解算模块,用于根据所述观测量和所述卫星导航电文进行联合定位,获得导航结果;
输出模块,用于根据海用电子设备制定的标准格式NMEA格式输出所述导航结果。
2.如权利要求1所述的装置,其特征在于,所述卫星导航电文包括卫星的星历、卫星工作状态、时间装置、卫星历书。
3.如权利要求1所述的装置,其特征在于,所述观测量提取模块,具体用于根据基带处理数据得到被跟踪卫星的伪距观测量、或多普勒观测量、或载波相位观测量。
4.如权利要求3所述的装置,其特征在于,所述伪距观测量是利用卫星信号的传输时间乘以光速得到信号的传输距离;所述卫星信号的传输时间指卫星信号发送时间tsv减去接收机本地时间tu。
5.如权利要求3所述的装置,其特征在于,所述多普勒观测量是在得到载波NCO的频率值后,将所述载波NCO的频率值与射频前端的标称中频值fc相减获得。
6.如权利要求2所述的装置,其特征在于,所述定位解算模块具体包括:
计算模块,用于根据所述卫星的星历计算出卫星的位置;
定位解算模块,用于利用所述卫星的位置和所述观测量对用户位置进行定位解算。
7.如权利要求6所述的装置,其特征在于,所述计算模块具体用于:
第一步:计算卫星运动的平均角速度n;
第二步:计算观测瞬间卫星的平近点角Ms;
Ms=M0+n(t-toe);式中,M0为参考时刻toe时的平近点角,由广播星历给出;
第三步:计算偏近点角:Es=Ms+essinEs;
第四步:计算近点角fs:
第五步:计算升交距角式中,ω为近地点角距,由广播星历给出;
第六步:计算摄动改正项;
第七步:对升交距角卫星矢径r、轨道倾角i进行摄动改正;
第八步:计算卫星在轨道面坐标系中的位置;
第九步:计算观测瞬间升交点的经度L;
第十步:计算卫星在瞬时地球坐标系中的位置;
第十一步:计算卫星在协议地球坐标系中的位置。
8.如权利要求7所述的装置,其特征在于,所述平均角速度n通过以下步骤获得:
首先根据广播星历中给出的参数计算出参考时刻toe的平均角速度n0:式中,GM为万有引力常数G与地球总质量M的乘积,GM=3.9860047×1014m3/s2,然后根据广播星历中给定的摄动参数△n计算观测时刻卫星的平均角速度n:n=n0+△n。
9.如权利要求7所述的装置,其特征在于,所述计算摄动改正项包括以下步骤:
通过广播星历中给出的Cuc、Cus、Crc、Crs、Cic、Cis6个摄动参数,求出由于地球引力场位函数的二阶带谐系数项而引起的升交距角μ的摄动改正项δμ、卫星矢径r的摄动改正项δr,和卫星轨道倾角i的摄动改正项δi,计算公式如下:
10.如权利要求7所述的装置,其特征在于,所述对升交距角卫星矢径r、轨道倾角i进行摄动改正的步骤为:
r=as(1-escosEs)+δr,式中:as为卫星轨道的长半径,i0为toe时刻的轨道倾
角,由广播星历中的开普勒参数给出;为i的变化率,由广播星历中的摄动参数给出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510937376.XA CN105510935A (zh) | 2015-12-15 | 2015-12-15 | 一种用于导航的装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510937376.XA CN105510935A (zh) | 2015-12-15 | 2015-12-15 | 一种用于导航的装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105510935A true CN105510935A (zh) | 2016-04-20 |
Family
ID=55719044
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510937376.XA Pending CN105510935A (zh) | 2015-12-15 | 2015-12-15 | 一种用于导航的装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105510935A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107643527A (zh) * | 2017-09-08 | 2018-01-30 | 桂林电子科技大学 | Gps共视信号模拟生成方法及gps共视信号模拟器 |
CN108267765A (zh) * | 2018-03-13 | 2018-07-10 | 北京沙谷科技有限责任公司 | 使用不完整gnss信息的mems组合导航建模方法 |
CN108415044A (zh) * | 2018-02-09 | 2018-08-17 | 广州吉欧电子科技有限公司 | 一种处理***、处理方法和gnss接收设备 |
CN109085627A (zh) * | 2018-06-11 | 2018-12-25 | 北京讯腾智慧科技股份有限公司 | 一种提高全球导航卫星***单频接收机定位精度的方法 |
CN112255648A (zh) * | 2020-10-20 | 2021-01-22 | 腾讯科技(深圳)有限公司 | 运动状态检测方法、装置、电子设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101132191A (zh) * | 2007-10-15 | 2008-02-27 | 北京航空航天大学 | 一种gnss接收机基带信号处理方法 |
CN101592723A (zh) * | 2008-05-30 | 2009-12-02 | 凹凸电子(武汉)有限公司 | Gps接收机及其定位方法 |
CN101609140A (zh) * | 2009-07-09 | 2009-12-23 | 北京航空航天大学 | 一种兼容导航接收机定位***及其定位方法 |
-
2015
- 2015-12-15 CN CN201510937376.XA patent/CN105510935A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101132191A (zh) * | 2007-10-15 | 2008-02-27 | 北京航空航天大学 | 一种gnss接收机基带信号处理方法 |
CN101592723A (zh) * | 2008-05-30 | 2009-12-02 | 凹凸电子(武汉)有限公司 | Gps接收机及其定位方法 |
CN101609140A (zh) * | 2009-07-09 | 2009-12-23 | 北京航空航天大学 | 一种兼容导航接收机定位***及其定位方法 |
Non-Patent Citations (2)
Title |
---|
工测顽童: "GPS卫星的坐标计算", 《百度文库》 * |
谢钢: "《GPS原理与接收机设计》", 31 July 2009 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107643527A (zh) * | 2017-09-08 | 2018-01-30 | 桂林电子科技大学 | Gps共视信号模拟生成方法及gps共视信号模拟器 |
CN108415044A (zh) * | 2018-02-09 | 2018-08-17 | 广州吉欧电子科技有限公司 | 一种处理***、处理方法和gnss接收设备 |
CN108415044B (zh) * | 2018-02-09 | 2020-09-15 | 广州吉欧电子科技有限公司 | 一种处理***、处理方法和gnss接收设备 |
CN108267765A (zh) * | 2018-03-13 | 2018-07-10 | 北京沙谷科技有限责任公司 | 使用不完整gnss信息的mems组合导航建模方法 |
CN108267765B (zh) * | 2018-03-13 | 2021-08-03 | 北京沙谷科技有限责任公司 | 使用不完整gnss信息的mems组合导航建模方法 |
CN109085627A (zh) * | 2018-06-11 | 2018-12-25 | 北京讯腾智慧科技股份有限公司 | 一种提高全球导航卫星***单频接收机定位精度的方法 |
CN109085627B (zh) * | 2018-06-11 | 2019-07-12 | 北京讯腾智慧科技股份有限公司 | 一种提高全球导航卫星***单频接收机定位精度的方法 |
CN112255648A (zh) * | 2020-10-20 | 2021-01-22 | 腾讯科技(深圳)有限公司 | 运动状态检测方法、装置、电子设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Skog et al. | In-car positioning and navigation technologies—A survey | |
US8788200B2 (en) | Method and system for a data interface for aiding a satellite positioning system receiver | |
US7940210B2 (en) | Integrity of differential GPS corrections in navigation devices using military type GPS receivers | |
CN105510935A (zh) | 一种用于导航的装置 | |
CN100437142C (zh) | 基于地基增强***的电离层延迟误差分离方法及地基增强*** | |
CN101950025B (zh) | 用于局域增强***的数据质量监测方法 | |
CN103837879B (zh) | 基于北斗***民用载波相位组合实现高精度定位的方法 | |
CN105467415B (zh) | 一种基于差分气压高度约束的小型无人机rtk相对定位方法 | |
CN101395443B (zh) | 混合定位方法和设备 | |
CN102253399B (zh) | 一种利用载波相位中心值的多普勒差分补偿测速方法 | |
CN101833080A (zh) | 一种利用gps***附加约束条件的载体姿态测量方法 | |
Johnson et al. | Where in the world are my field plots? Using GPS effectively in environmental field studies | |
Dovis et al. | Recent advancement on the use of global navigation satellite system-based positioning for intelligent transport systems [guest editorial] | |
CN106707322A (zh) | 基于rtk/sins的高动态定位定姿***及方法 | |
Bhardwaj et al. | Satellite navigation and sources of errors in positioning: a review | |
CN105510946A (zh) | 一种bds卫星载波相位整周模糊度快速解算方法 | |
Wright et al. | The effectiveness of global positioning system electronic navigation | |
CN105158781A (zh) | 一种地球gnss卫星实现外层空间用户定位的方法 | |
CN105572705A (zh) | 一种导航装置 | |
CN108351423A (zh) | 铁路定位方法和*** | |
Wells et al. | NAVSTAR performance analysis | |
CN105510936A (zh) | 星载gnss联合定轨方法 | |
CN105572706A (zh) | 一种自动导航驾驶*** | |
CN103645488B (zh) | 用于确定gps接收器位置的***和方法 | |
Yu | Improved positioning of land vehicle in ITS using digital map and other accessory 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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160420 |