GPS卫星信号的跟踪.pptVIP

6.3.3　BASS方法中的载波频率更新 为简化起见，采用1 ms的电文，利用离散傅立叶变换(DFT)，如果输入信号是x(n), DFT输出就是： （6－40） 其中，k表示某一特定的频率分量，N是数据点总数。如果x(n)是从数字化的正弦波中得到的，那么最大幅值|X(ki)|就代表了输入信号的频率，由|X(ki)|的实部Re和虚部Im得到相位角为： (6－41) 　　我们在5.4.5节讨论过，利用相位角θ来推算输入信号的精频。图6－5所示为跨越两个不同时域的数据点。在每个时域， 计算同一个X(ki)，对应的相位角为θn和θn+m，它们在时间上相隔m，得到的精频值为： （6－42） 这个关系得到了比DFT结果还要精细得多的精频，得到的频率精确度取决于测量的角精度。 图6－5　从连续两组数据中得到的相位角 6.3.4　Kernel函数的不连续性 　在传统DFT过程中，方程(6－40)中的k值是一个整数。但是，在实际应用方程(6－40)于跟踪过程时，k值往往是非整数，因为在使用整数k值时，从Kernel函数e-j2πkn/N中产生的频率远离输入频率。如果k值远离输入信号，方程(6－40)给出的值会很小，这也就意味着跟踪过程的灵敏度将大大降低。为了避免这个问题，k值应尽可能地保持与输入频率相近。与输入频率非常接近的k值同时还可减小频率不确定度。这样，k值就往往不是一个整数了。 　　当k值是整数时，Kernel函数e-j2πkn/N的初相是零，且从两组连续计算中，得到的初相角是连续的。第一组的起始点是n=0，第二组的起始点是n=N，如果k是整数，则很容易得到： （6－43） 如果k不是整数，上式就不成立了。下面举例说明这一点。 　　假设N=256，取n=0～255，对于任何一个整数k值，由e-j2πkn/N以产生256个数据点。同样的设置放在另一组数据中，总共产生512个数据点，这样使得从第一组数据的最后一个数据点， 即数据点256，到第二组的第一个数据点，即数据点257，是连续的，因为由e-j2πkn/N产生的值是复数，连续性只能以图表的形式通过e-j2πkn/N的实部和虚部显现。图6－6(a)给出了k=20时的实部和虚部的结果，图中仅画出了从数据点240到270的图形，它是连续的。图6－6(b)给出了k=20.5时的结果，在点256到257之间e-j2πkn/N的实部和虚部均是不连续的。不连续性将影响方程(6－42)的应用。 图6－6　e-j2πkn/N的实部与虚部图（N=256，n=0～255） (a)k=20; (b) k=20.5 　　图6－7给出了一个连续波形信号和e-j2πkn/N的实部、虚部两部分。如果e-j2πkn/N是连续的，两组DFT信号可以看做是输入信号与一连续复信号的相关，这样，就可以用式(6－41)与式(6－42)来计算精频。如果Kernel函数不连续，则两组DFT数据可看做是输入信号与两组Kernel函数的相关。为了确保无误地应用式(6－42)， 必须考虑到Kernel函数的不连续性。 图6－7　 输入信号和e-j2πkn/N的两组实部 （a）输入信号; （b）两组连续Kernel函数e-j2πkn/N的实部，无不连续; （c） 两组连续Kernel函数e-j2πkn/N的实部， 有不连续 　　不连续性出现在计算n=N时的相位角时，在前述例子中， 用N=256和n=0～255生成Kernel函数值。为了产生连续Kernel函数，e-j2πkn/N=e-j2πk(n=256)必须等于e0（或0°）。 如果k是整数，这个关系就成立；如果k不是整数，这个关系就不成立，且e-j2πk和e0（或0°）之间的相位差是不连续的，这个相位必须在应用式(6－42)之前去掉。 　　如果在每毫秒结束时，减去Kernel函数的相位差，在相邻两个毫秒之间会出现两种情况： 一是它们之间无相位改变， 即是连续的；另一种是导航电文引起两者之间相位相差π。 由于输入信号中含有噪声，因此相位差不会精确等于0或　π。 如果相位转移接近0或2π，则认为无相位偏移；如果相位偏移接近π，则认为相位偏移为π。通常，在±π/2处设置一个门限。 如图6－8所示，门限设置在π/2处和3π/2处，如果相位差的绝对值在π/2～3π/2之间，就可以归类为相位变化π，否则，认为无相位变化。相位相差π的情况不会发生在20 ms内，而只会发生在20 ms的倍数时。 图6－8　相位偏移的门限示意图 6.3.5　C/A码起始点的精确性 输入信号以5 MHz的频率数字化，或者说，每个数据点相隔200 ns。这样的时间分辨率，与之相应的距离分辨率大约是60 m（3×108×200×10-9），这个