GPS软件接收机基础精.docx

1、GPS软件接收机基础精第 2章 坐标系统和时间系统GPS 卫星环绕地球运行, 卫星的观测量是卫星位置 /速度和接收机位置 /速度 函数,为了能够描述观测量,有必要定义适宜的坐标系统和时间系统。2.1 参考坐标系统2.1.1 地心惯性坐标系统 (ECI为了测量和确定 GPS 卫星的轨道,使用地心惯性坐标系统是很方便的 6。 在 ECI 系统中,原点位于地球的质心, GPS 卫星遵循牛顿的运动和引力定律。 在一个典型的 ECI 坐标系中, xy 平面选为与地球的赤道面重合, x 轴相对于天 球永久固定地指向某一方向, z 轴垂直于 xy 平面指向北极点, y 轴的选择使得 坐标系构成一右手坐标系。

2、地面站确定与预测 GPS 卫星的轨道就是在 ECI 坐 标系中进行的。地球的形状是一个扁状的椭球,由于日月对于地球赤道的突出部分的引力 作用,使得地球的赤道面相对于天球运动。因为 x 轴的定义是相对于天球的, 而 z 轴的定义是相对于赤道面的,这样地球运动的不规则性将使得如上定义的 ECI 系统并非真正的惯性系统。此问题的解决办法为:定义某一瞬间的坐标轴 指向。 GPS 的 ECI 坐标系统以 UTC (USNO 时间 2000年 1月 1日凌晨的赤 道面方向作为其基准, x 轴选为由地球质心指向春分点方向, y 轴和 z 轴的定 义和上面一样。由于坐标轴的指向是固定的,因此以此种方式定义的

3、ECI 坐标 系统对于 GPS 来说可以认为是惯性的。2.1.2 地心地固坐标系统 (ECEFGPS 接收机的位置计算常在一个随地球旋转的坐标系统中进行,此种坐标 系统称为地心地固坐标系统。在这样的坐标系统中,可方便地计算接收机的纬 度,经度和高度。如同 ECI 坐标系统一样, GPS 使用的 ECEF 坐标系统的 xy 平面也与地球的赤道面重合。不同的是,在 ECEF 坐标系统中, x 轴指向 0经 度, y 轴指向东经 90度。因而, x 轴和 y 轴将随地球自转而旋转。在 ECEF 系统中, z 轴垂直赤道面指向地理北极点。在计算 GPS 接收机位置之前, 有必要把根据星历信息计算出来的

4、 ECI 坐标 系中的卫星位置和速度转换到 ECEF 中,这种转换可以通过旋转矩阵来实现, 具体细节请参考 2.1.4小节。 一般情况下, 计算出来的接收机位置都是用 ECEF 坐标表示,对于某些应用需要把这种坐标转换成接收机的经度、纬度和高度。 为实现这种转换, 需要定义一个描述地球的标准物理模型, 在 GPS 中选用的是 世界大地 84坐标系统(WGS-84 。下面给出的是该标准下的地球物理模型的 参数:地球长半轴长:6378.137a km=地球扁率 :1/298.257 223 563f =各类参数间变换关系为:( 11bf b a faee=其中, e 为地球的偏心率, e 为地球的

5、二阶偏心率。2.1.3 站心坐标系站心坐标系又称地方坐标系 (Local Coordinate System,属于左手笛卡尔坐 标系统。如 图 2-1所示, (, , x y z 表示站心坐标系,它的原点位于地方点(1111, ,P x y z , z 轴垂直于该点地球椭球的切平面指向上方, x 轴在该切平面内指向正北, y 轴在该切平面内指向东。对于该坐标系中的任意一点2P , A 表示方位角, Z 是天顶距, E 是仰角, d 是2P 点在站心坐标系中的向径 d 的长度。其 中, A 从正北方向按顺时针方向开始计起的; Z 是 z 轴和向径 d 的夹角; E 的取值为:当2P 点在切平面以

6、上时为正值,在切平面以下时为负值。 图 2-1站心坐标系统 2.1.4 不同坐标系间坐标变换在 GPS 数据处理中, 常常需要把一个向量的坐标在不同坐标系统间进行转 换。下面首先介绍坐标变换的基础,然后给出几种常见的坐标变换公式。2.1.4.1 坐标变换与旋转矩阵如果两个笛卡尔坐标系统共原点且都是左手或都是右手坐标系,那么二者 之中的任何一个坐标系都可通过三次连续的旋转转换到另一坐标系统中。这三 个旋转矩阵为:(123100R 0cos sin , 0sin cos cos 0sin R 010, sin 0cos cos sin 0R sin cos 0001=(2.1 其中, 代表旋转角,

7、 从旋转轴的正方向看如果逆时针旋转则该角取正值。 1R ,2R 和 3R 分 别 称 作 x , y 和 z 轴 的 旋 转 矩 阵 。 对 于 任 何 旋 转 矩 阵 R , 有 (1T R R R =;也就是说,旋转矩阵是正交矩阵,其中 1R 和 T R 分 别表示矩阵的 R 的逆和转置。对于两个具有不同原点和不同单位长度的笛卡尔坐标系统,一般的转换公 式表示为0R n old =+x x x (2.2或者 000R n old n old n old x x x y y y z z z =+其中, 是比例因子, R 是转换矩阵。 n x 和 old x 分别表示同一向量新的和旧的坐标;

8、0x 表示平移向量,它是旧的坐标系的原点在新的坐标系中的坐标。2.1.4.2 ECEF与大地坐标坐标系间的坐标变换ECEF 坐标与大地坐标间转换有不止一种实现方法,比如文献 4、 6和 7中分别介绍了三种不同的实现,其中文献 7中的推导最为详细。为方便计,这 里仅给出文献 6中的变换公式。如 图 2-2所示。 x图 2-2 ECEF与大地坐标系(, , u u u U x y z 表示接收机在 ECEF 坐标系中的位置, 表示接收机所处位置的经度, 表示接收机所处位置的纬度, h 表示接收机所处位置相对于参考地球球面的高度。已知接收机在 ECEF 中的坐标, 求接收机在大地坐标系中的坐标按 表

9、 2-1进 行计算 6。已知接收机在大地坐标系中的坐标参数 , 和 h ,按照下式可以计算得 接收机在 ECEF 坐标系中的坐标参数。 2cos cos sin cos 1sin sin u u u h x y h z a e h +=+u (2.3 其中, , a e 分别表示地球的长半轴长和偏心率。表 2-1 根据 ECEF 参数计算大地坐标系下的经度,纬度和高度 2.1.4.3 ECEF与站心坐标系的变换为了表示 ECEF 与站心坐标系的变换关系, 两个坐标系一起绘制于 图 2-3中。 其中, x, y,z 对应 ECEF 坐标系, x, y,z 对应站心坐标系,接收机的位置位于(u u

10、 u U x , y ,z ,其大地坐标为 (, ,h 。 S 表示卫星的位置,它在 ECEF 中的坐标用 (x, y,z 表示,在站心坐标系中的坐标用 (x, y,z 表示。坐标变换过程如下:首先计算向量 US在 ECEF 坐标系中的坐标,然后把 ECEF 坐标系绕 z 轴逆时针旋转 ,接下来再绕 y 轴逆时针旋转 (o 90, 由于 站心坐标系是左手坐标系,因此最后还要把 x 轴的正轴反向,上述过程用矩阵 表示为 图 2-3 ECEF与站心坐标系(900900010000190090o ou u o o u cos sin x cos sin x x y sin cos y y z z z

11、 sin cos =-100010001(2.4 化简为u u u x sin cos sin sin cos x x y sin cos 0y y z cos cos cos sin sin z z =(2.52.2 时间系统GPS 分发一种形式的协调世界时 (UTC时间,能够使世界范围内时间得以 同步。下面分别介绍 UTC 时间和 GPS 系统时的概念。2.2.1 UTC时间协调世界时(UTC 是一种复合的时间系统,它的秒长严格等于原子时的 秒长。原子时是一种以物质内部原子运动的特征为基础的时间系统。协调世界 时采用闰秒的办法,使之与世界时的时刻接近。2.2.2 GPS系统时GPS 系统参

12、考美国海军天文台(USNO 维护的 UTC 时间。 GPS 系统时是 一种书面定义的时间标尺 (time scale,它是对卫星内部的原子钟时间和各不同 地面控制段的时间分量统计处理的结果。 GPS 系统时是一个连续的时间标尺, 它不存在跳秒。 GPS 系统时和 UTC (USNO 在 1980年 1月 6日零时重合。 GPS 系统时的一个历元用 GPS 星期数和周内秒计数表示。 GPS 星期数从 1980年 1月 6日零时开始计起,初值为 0;周内秒计数从每个星期的星期日零时计 起,范围是 0604800秒。第 3章 观测量GPS 接收机要实现定位功能, 除了必备处理 GPS 扩频信号的能力

13、外, 还应 能提供观测量给导航解算模块。常用的观测量包括伪距、伪距变换量、载波相 位和多普勒频率。 在这一章里首先给出一般 GPS 接收机的结构, 然后给出不同 类型观测量的定义及常用的获取方法。现代 GPS 接收机设计大部分都是数字接收机。一直以来,这类接收机的设 计向着组件集成度越来越高的方向发展,而且这种趋势仍在继续。基于这一原 因,这里用一个现代数字 GPS 接收机的高层结构图来代表一般的 GPS 接收机结 构,如 图 3-1所示。 图 3-1 GPS接收机结构图 【6】接收机的信号处理是 GPS 接收机的核心,它执行如下的功能:(1把接收信号分配到多个信号处理通道以进行多颗卫星的信号

14、处理;(2产生 PRN 码;(3捕获卫星信号;(4进行码跟踪和载波跟踪;(5从卫星信号中解调出导航电文;(6从卫星信号的 PRN 码中提取码相位(伪距观测量;(7从卫星信号的载波中提取载波频率(伪距变化率 、载波相位和伪距 变化量等观测量;(8估计 GPS 系统时间。信号处理模块的输出是伪距、 伪距变化率、 伪距变化量和卫星的导航电文, 所有的数据都为导航应用处理模块使用。导航应用处理模块控制信号处理过程,并对信号处理模块输出的观测量进 行数据处理以满足不同应用的需求。3.1 伪距及其测量方法 3.1.1 伪距定义假设一颗 GPS 卫星在系统时刻 s T 发射某一码相位, GPS 接收机在系统

15、时刻 u T 精确地接收到该码相位。如果在卫星和接收机天线相位中心间没有传播延 迟 , 则 理 想 的 真 实 距 离 将 为 (u s r c T T =米 , 其 中 c 表 示 光 速(829979245810. m /s , s T 和 uT 以秒为单位。然而, GPS 信号传播过程中有几种延迟, 如果不消去它们的影响, 那么将会使上面的测量距离比真实的距离长。 这些延迟包含与频率无关的对流层延迟 (trop t 和与频率有关的电离层延迟(ion t ,它们共同合成了大气层延迟 (atm t 效应。同时还有天线相位中心与接 收机内部码相关点的接收机延迟 (dr t 。把这些延迟考虑进测

16、距方程中,方程 变为 (u s atm dr r c T T t t n =米,其中 r 表示真实距离, atm trop ion t t t =+以秒为单位, (dr t n 的索引 n 表示接收机的第 n 个通道 8。图 3-2说明了 GPS 测量过程中的时间关系。关于大气层等误差源的影响会 第 4章在做进一步说明。一般地, GPS 卫星时钟并非完美地与 GPS 系统时对齐, 它们的时钟相对于 真实的 GPS 系统时有微小的偏移。 不过, 卫星的这种时间偏移和漂移特性精确 地为地面控制段测量,并通过导航电文把时钟的修正以多项式系数 0f a , 1f a 和2f a 的形式发送给用户。这样

17、,测距方程需进一步修改,从而包含进卫星发射 时间偏移项,即 (u s s atm dr r c T t t t t n =+米,其中 s s s T t t =,且 s t 为 卫星信号依据自己内部时钟得出的发射时间,当用 s t 对它修正即等于信号真实的 GPS 发射时间 s T 。 sReceive 符号说明:s u u atm dr T T t t t =s s u t t t 真实的信号的发射时间 卫星时间相对于系统时的钟差 真实的信号的接收时间 接收机时间相对于系统时的钟差卫星钟面信号的发射时间 卫星与接收机间的所有大气延迟 接收机钟面的信号接收时间 接收机天线与相关器间的所有延迟图

18、 3-2 GPS发射时间和接收时间关系用户接收机时钟相对于真实的 GPS 系统时间也是没有对齐而且是缓慢漂移 的。因此,用户必须使用至少 4颗卫星来计算除了接收机位置 (u u u x , y ,z 以外 的接收机钟差。最后,修改测距方程使它反应用户时间相对于系统时的偏差, 即(u s s atm dr u r c t t t t t t =+ (3.1方程 (3.1反映了用户接收机用以确定与卫星间真实距离的测量量。伪距是 指进行延迟改正及决定接收机钟差之前的原始测量量,它的定义如下:(u s c t t =(3.2它可以看成是接收机与卫星间的视在距离 (apparent range, 使用发

19、射时间 st 和接收时间 u t 的差乘以光速 c 得到。两个时间偏差修正必须考虑,同时也必须 修正路径延迟才能近似真实的距离,方程 (3.1以伪距的形式可以表示如下:(s atm dr u r c t t t t =+ (3.33.1.3 伪距测量方法对于观测量的测量,不同的接收机设计有不同的实现。这里为了说明伪距 的测量方法, 首先给出一个 GPS 接收机的基带信号处理框图, 该图是在文献 10中的一个基带信号处理框图的基础上修改而成。 图 3-3 GPS接收机基带信号处理框图GPS 数字接收机并不是直接地测量伪距。它直接测量的是码跟踪环的再生 伪码, 因为再生伪码所表示就是卫星信号的发射

20、时刻 s t 。 对于 P 码来说, 捕获成 功后,我们获得了接收信号的 P 码的码相位估计,而 P 码的码相位在一星期内与 时间是一一对应的,因而我们可以把估计的 P 码相位所对应的时间转换为码累 积器的形式,码累积器在此基础上进行累积表示的就是信号的发射时刻。码累 积器的实现表示为 图 3-4。 码发生器时钟32321212c , , f 秒 121c c, f ms f 秒 2992, ,量化单位:毫秒 1007996, ,量化单位:秒 CLOCK图 3-4码累积器当伪码为 P 码时, 6102310c f . Hz =,发射时刻可以表示为3211200026s c c t CODE _

21、NCO CHIP _cnt f f BIT _cnt . SUBF _cnt , =+单位:秒(3.4当伪码为 C/A码时,由于 C 码的码周期为 1毫秒,所以仅由码相位并不能 获得绝对的时间。 解决的办法是, 码环在首次找到帧头后, 就自动把 CHIP_cnt,BIT_cnt,和 SUBF_cnt三个计数器设置为零,然后在此基础上计数,此后码累 积器的计数值所对应的时间实际上就是相对于该帧头发射时刻的相对时间,而 此帧头的发射时刻在该帧解调后即可从导航电文中解析出, 若该时间用符号 0s t 表示,则 C/A码的发射时刻可以表示为3211200026604800604800s c cs s

22、st CODE _NCO CHIP _cnt f f t BIT _cnt . SUBF _cnt modt , if t+=+=(3.5 其中,对于 C/A码, 6102310c f . Hz =。对于接收时刻 u t ,由接收机的本地时钟提供,回顾式(2-8 ,伪距表示为(u s u dr atm s c t t r c t t t t =+ (3.6假设通道延迟 dr t 在形成伪距观测量时,我们已经提前消除掉了,则上式 变为(u s u s ion trop c t t r c t t c t t =+ (3.7这里定义的伪距还有一些其它误差因素没有考虑,这些误差将在第 4章予 以说明

23、。3.2 载波相位及其测量方法 3.2.1 载波相位的定义载波相位,又称载波差拍相位,是接收机接收的带有多普勒的卫星载波信 号与接收机产生的恒定频率信号的差拍信号的相位。这种观测量作为相关通道 的副产品而得到或产生自平方通道。平方通道用接收到信号去乘以自己从而得 到载波的二次谐波,并且平方后的信号不再含有码的调制 11。由于载波的波长比 P 码和 C/A码的波长都短,因而载波相位的测量量的精 度要比用伪码测到的伪距的精度高。对于 GPS 的 L1载波,波长大约是 20厘 米。根据经验法则,相位测量可精确到 1%的波长,这意味着可获得 2毫米的 精度。载波相位测量的主要缺点与整周模糊度有关。获得

24、初始的卫星与接收机间 的整周模糊度是比较困难的。 一个高质量的 GPS 接收机大多数时候, 在卫星与 接收机间相对位置变化时仍能维护整周计数。然而,由于各种原因,像信号噪 声、天线阻挡,都可以导致整周跳变。在多数时候,艰苦的后处理允许检测周 跳和对周跳进行修正。可是,整周跳变限制了载波相位测量量的实时应用。根据定义,载波相位的表达式如下(j j u u u s t t = (3.8其中, (j s t 表示卫星在卫星时间 s t 发射的载波信号的相位, (u u t 表示接收 机在接收机时间 u t 的载波信号的相位,二者的单位都是周 (cycles。在有些的文 献中,载波相位定义成 (j s

25、 t 减去 (u u t ,不过这并不重要。高稳振荡器在短时间间隔内,相位和频率满足如下关系(0t t t f t +=+ (3.9令 u s t t t =,有(0j u u s u s t t f t t =+ (3.10根据上式我们可以得到(0j u u s u s t t f t t = (3.11根据伪距测量一节的说明,我们知道接收时刻和发射时刻满足如下关系s s ion trop u u rt t t t t t += (3.12 需要注意是,载波相位测量中的电离层对于相位传播的影响与对伪码传播 的影响符号相反,相对于光在真空中的传播速度一个超前一个滞后,关于这一 点的详细推导可以

26、参考文献 4、 5和 6等。根据上式,有u s u s ion trop rt t t t t t c=+ (3.13这样,带有误差的载波相位的数学模型为(000u s ion trop fr f t t f t t c =+ (3.14 在实际应用中,载波相位在某一历元的测量基于接收机再生载波与卫星信 号载波的对齐,而并不知道哪一周期表示理想的周期同步 12。因此,总的相位total 包含测量得到的小数相位分量 (Fr 和从初始锁定历元 0t 到时刻 t 的相位的整数周期计数 (Int ,以及在初始历元 0t 未知的整数周期 N(00total Fr Int ;t ,t N t =+ (3.

27、15未知的周期计数 N 通常被称作整周模糊度 (integer cycle ambiguity。只要接 收机在观测阶段保持对载波的连续跟踪,那么在接收机与每颗卫星间只有一个 整周模糊度。然而,如果有失锁,将会有周跳引入。真 正 为 接 收 机 观 测 的 是 (0measured Fr Int ;t ,t =+。 因 此 ,(0total measured N t =+。在某一瞬间,某颗卫星与接收机间的观测方程可以写作(0000total u s ion trop measured fr f t t f t t N t c=+=+进而(0000measured u s ion trop fr

28、f t t f t t N t c =+ (3.16等式两边乘以载波波长 0L c f =,并定义L measured = (3.17则可得到以米为单位的相应的载波相位方程(u s ion trop L r c t t c t t N =+ (3.18与式 (3.16相比,方程增加了一个整周模糊度,而且电离层的误差项符号相反。3.2.2 载波相位的测量方法这一小节在上面的基础上给出一种载波相位的测量方法。假设接收到的第 j 颗卫星信号表示为(12j Rj s t A (t cos f t = (3.19 其中, Rj Tj dj f f f =+表示接收到的信号频率, 其中 T f 是发射频率

29、, 对于 L1载波, 不考虑钟差,该值为 0f =1575.42MHz,对于 L2载波,该值为 1227600f =. MHz , d f 是多普勒频率, j 表示第 j 颗卫星, 1A (t 是信号幅度、伪码和数据位合成的 信号的简写,而且这里只代表信号的某一支路。假设接收机产生的本地载波信号为(202r g s t A cos f t =+ (3.20 其中 0是初始相位, g f 表示本地载波的频率。若 0T g f f f ,测量的载波相位根据前面的定义表示(0000t t t measured 00dj dj t t t f dt f f dt f dt =+=+ (3.21 其中, 0 对应 0,单位为 cycle 。根据多普勒的定义,参见 3.3小节。 measured 实际上表示的是接收机与卫星 j 间的距离自 0t 开始的变化量,正值表示距离在增加,负值表示距离在减小。从 这里可以知道,接收机只要对载波的多普勒连续积分就得到了载波相位观测 量。参照图 3-3,如果输入的信号是零中频信号,即 (r s t 与 (j s t 混频后的差频 信号。实际上