全球定位系统定位结果的自适应数据窗长度平滑方法
技术领域
本发明涉及全球卫星定位与导航系统领域,尤其涉及针对全球定位系统定位结果的自适应地确定选取数据的窗口长度,进而实现对定位结果平滑方法的领域。
背景技术
全球卫星定位与导航系统,例如美国的全球定位系统(GPS),包括一组发送GPS信号的一个卫星星座(又被称为Navstar卫星),该GPS信号能被接收机用来确定该接收机的位置。卫星轨道被安排在多个平面内,以便在地球上任何位置都能从至少四颗卫星接收该种信号。更典型的情况是,在地球上绝大多数地方都能从六颗以上卫星接收该种信号。
每一颗GPS卫星所传送的GPS信号都是直接序列扩展频信号。商业上使用的信号与标准定位服务(SPS)有关,而且被称之为粗码(C/A码)的直接序列二相扩展信号,在1575.42MHz的载波下,具有每秒1.023兆码片的速率。伪随机噪声(PN)序列长度是1023个码片,对应于1毫秒的时间周期。每一颗卫星发射不同的PN码(Gold码),使得信号能够从几颗卫星同时发送,并由一接收机同时接收,相互间几乎无干扰。术语“卫星星号”和这个PN码相关,可以用以标示不同的GPS卫星。
GPS的调制信号是导航电文(又被称为D码)和PN码的组合码。导航电文的速率为每秒50比特。D码的基本单位是一个1500比特的主帧,主帧又分为5个300比特的子帧。其中子帧一包含了标识码,星种数据龄期,卫星时钟修正参数信息。子帧二和子帧三包含了实时的GPS卫星星历(ephemeris),星历是当前导航定位信息的最主要内容。子帧四和子帧五包含了1-32颗卫星的健康状况,UTC校准信息和电离层修正参数及1-32颗卫星的历书(alamanc)。历书是卫星星历参数的简化子集。其每12.5分钟广播一次,寿命为一周,可延长至2个月。
由于各种系统误差的综合作用,接收机获得的原始定位结果存在一定相对误差。这些误差包括卫星时钟误差、星历预测误差、相对论效应、电离层效应、对流层效应、接收机噪声和多径效应。由于这些误差的存在,往往需要进一步对原始数据进行处理以提高精度。一个常用的方法是对原始定位数据流进行平滑,也即低通滤波。常规的方式使用固定长度的数据进行平滑,虽然简单易行,但存在一些显而易见的缺点。平滑所用的数据长度太短,则平滑后结果方差依然很大,精度较低;平滑用的数据长度太长,平滑后结果延迟较大,降低了导航的实时性。对于高速运动的接收机,这个问题尤为明显,限制了很多导航应用。如何高效的实现高精度、低方差、小延迟的平滑是GPS接收机一个重要的研究方向。
发明内容
本发明目的是提供一种全球定位系统接收机定位结果的自适应数据窗长度的平滑方法。
本发明提出的方法利用计算出的接收机速度,自适应地计算出用于对定位结果进行平滑的数据窗长度,根据该窗口获取初始数据进行平滑。包所述括以下步骤:
获得接收机的原始定位结果;
写入平滑缓存;
计算接收机的速度;
根据接收机速度自适应地确定数据窗的长度;
利用该长度的数据窗获得用以平滑的数据进行平滑。
进一步,方法中所述数据窗的长度与接收机的速度成反比。
进一步,是利用载波相位方法或者卡尔曼滤波方法计算接收机的速度。
进一步,所述接收机原始定位结果写入平滑缓存时,更新的原始定位结果和计算出的平滑值差距大于阈值时,则舍弃该值,不写入平滑缓存。
利用该方法可以实现低速情况下的高精度、低方差和高速情况下平滑结果的实时性,降低数据输出延迟。这种自适应数据窗长度的平滑适用于任意的多通道全球定位系统接收机,相对其他方法简单易用,复杂度大大降低,同时方便和其他各种平滑方法组合。
附图说明
图1是典型的全球定位系统接收机实现定位的示意图;
图2是典型的GPS接收机定位流程;
图3是针对全球定位系统定位结果的自适应数据窗长度平滑方法流程框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
图1描述了一个典型的全球定位系统接收机实现定位的示意图。已知四颗卫星的位置以及这四颗卫星发射信号到达接收机的时间,即可以计算出接收机的位置和接收机时间相对卫星星载时钟时间的误差值。应当理解,由于各种系统误差的综合作用,接收机获得的原始定位结果存在一定相对误差,这些误差包括卫星时钟误差、星历预测误差、相对论效应、电离层效应、对流层效应、接收机噪声和多径效应。由于这些误差的存在,往往需要进一步对原始数据进行处理以提高精度。一个常用的方法是对原始定位数据流进行平滑,也即低通滤波。常规的方式使用固定长度的数据进行平滑,虽然简单易行,但存在一些显而易见的缺点:平滑所用的数据长度太短,则平滑后结果方差依然很大,精度较低;平滑用的数据长度太长,平滑后结果延迟较大,降低了导航的实时性。对于高速运动的接收机,这个问题尤为明显,限制了很多应用。本发明提出了一种自适应数据窗长度的平滑方法,能够高效地解决上述问题。
图2描述了一个典型的GPS接收机定位流程。图2所示过程是从接收机上电初始化201开始,直至解算出用户位置106结束。通常接收机上电后进入捕获状态202,搜索视线内的卫星、该卫星的载波频率和PN码的码相位。这个状态下对载波频率的搜索是粗糙的,通常在数百赫兹的量级。之后进入频率牵引状态203,把本地频率牵引到和卫星载波频率相差几个赫兹的量级,同时进行比特同步。完成频率牵引后,接收机进入跟踪状态204,完成帧同步,即可以进入解调电文状态205,将载波上调制的电文解调出来用于在随后的解算状态206下计算接收机位置。
图3描述了针对全球定位系统定位结果的自适应数据窗长度平滑方法的细节。术语“数据窗长度”是指对原始数据流进行加窗操作获取一定长度数据过程中,所用窗的长度。数据窗的类型不是限定性的,可以是最简单的矩形窗,也可以是三角窗,汉明窗等等,往往由所采用的平滑算法决定。接收机首先通过图2所示的流程解算出接收机位置信息301。进而将该结果写入平滑缓存302,用于随后的平滑操作。该操作还包含了排除异常点的步骤。当更新的原始定位结果和计算出的平滑值差距太大时,则认为是错误定位点,放弃写入平滑缓冲区的操作,将该值排除出平滑的计算过程。该平滑缓存的长度与数据窗的最大长度相同。进一步地,计算接收机的速度303。计算接收机的方法有很多种,大致可以分为载波相位法,卡尔曼滤波法和定位结果差分法。定位结果差分法由于误差较大且和定位结果的误差相关,因此不推荐使用。而载波相位法和卡尔曼滤波法都能很好的满足本方法的要求。本发明的一个较优实施例使用载波相位法计算速度,本发明的另一个较优实施例使用卡尔曼方法同时计算接收机位置和速度。计算接收机位置301和速度303的相对顺序不是限定性的,可以调整。事实上,如上文所诉,本发明的一个较优实施例利用卡尔曼方法同时获得接收机位置和速度信息。对于两种计算接收机速度的方法的具体细节不再赘述,本领域熟练人员都应熟悉。
进一步地,利用获得的接收机速度自适应地调整数据窗的长度304,数据窗的长度由接收机的速度决定,与接收机的速度成反比。以实现在速度较低的时候数据窗较长,获得较小的方差的定位结果,达到较高的精度;在速度较高时数据窗较短,获得较低的延迟。本发明的一个较优实施例使用一个一阶方程获得数据窗长度。本发明的另一个较优实施例使用一个二阶方程获得数据窗长度。本发明对计算数据窗长度的方法不是限定性的,但对接收机速度和数据窗的相对关系是限定性的。
进一步地,利用自适应长度的数据窗选取该长度的数据305,进而进行平滑306。本发明对平滑的方法不是限定性的。本专利的一个较优实施例使用了简单的平均实现平滑。本专利的另一个较优实施例使用三次拟合实现平滑。获得平滑的结果后即可以输出平滑定位结果307。
尽管本发明的方法和装置是参照GPS卫星来描述的,但应当理解,这些原理同样适用于采用假卫星(pseudolites)或卫星与假卫星的组合的定位系统。假卫星是一种基于地面的发射机,它传播调制在L频段在波信号上PN码(与GPS信号相似),并且通常是与GPS时间同步的。每一发射机可以被赋予一个独特的PN码,从而允许由远端接收机进行识别。假卫星用在这样的情况下,即,来自轨道卫星的GPS信号缺失,如隧道、矿山、建筑物或者其他的封闭区及明显遮挡。这里所使用的术语“卫星”包括假卫星或假卫星的等效,而这里所使用的术语GPS信号包括来自假卫星或者假卫星等效的类似GPS的信号。
在前面的讨论中,本发明是参照美国全球定位系统(GPS)来描述的。然而,应当理解,这些方法同样适用于类似的卫星定位系统,如俄罗斯的格洛纳斯(Glonass)系统,欧洲的伽利略(Galileo)系统和中国的北斗1及北斗2系统。所使用的术语“GPS”还包括这样一些卫星定位系统,如俄罗斯的格洛纳斯(Glonass)系统,欧洲的伽利略(Galileo)系统和中国的北斗1及北斗2系统。术语“GPS信号”包括来自另一些卫星定位系统的信号。
上文中,已经描述了通过快速启动与定位方法实现GPS接收机的系统。尽管本发明是参照特定实施例来描述的,但很明显,本领域的普通技术人员,在不偏移权利要求书所限定的发明范围和精神的情况下,还可以对这些实施例作各种修改和变更。因此,说明书和附图是描述性的,而不是限定性的。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可轻易想到的变换或替换,都应涵盖在本发明的包含范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。