第5章GS基本原理及其Matlab仿真ppt课件

1、第5章 GPS卫星信号的捕获 第5章 GPS卫星信号的捕获 5.1 5.1 概述概述 5.2 GPS5.2 GPS卫星信号的多普勒效应卫星信号的多普勒效应5.3 GPS5.3 GPS卫星信号捕获的思索卫星信号捕获的思索5.4 GPS5.4 GPS卫星信号的捕获方法卫星信号的捕获方法 5.5 Matlab5.5 Matlab与与SimulinkSimulink仿真简介仿真简介 5.6 GPS5.6 GPS卫星信号捕获的例子卫星信号捕获的例子 5.7 5.7 关于捕获的一些子程序关于捕获的一些子程序 第5章 GPS卫星信号的捕获 5.1 概概 述述 为了跟踪和解码GPS信号, 首先要捕获到GPS信

2、号。将捕获到的GPS信号的必要参数立刻传送给跟踪过程，再经过跟踪过程便可得到卫星的导航电文。GPS卫星处于高速运动中，因此，其频率会产生多普勒频移。载波频率与C/A码的多普勒频移将在下面详细讲述。为覆盖高速卫星预期中的一切多普勒频率范围，捕获方法覆盖的频率范围必需在10 kHz之内。一旦捕获到GPS信号，立刻去丈量两个重要参数： C/A码的起始点和载波频率由于多普勒频移而变化。接纳机接纳到的一系列数据往往包含多个卫星信号，每个信号具有不同的C/A码的不同起始点和不同的多普勒频率。 第5章 GPS卫星信号的捕获 针对某个特定的卫星信号， 捕获过程就是要找到C/A码的起始点，并利用找到的起始点展开

3、C/A码频谱，一旦复现了C/A码的频谱， 输出信号将变成延续波Continuous Wave， CW，于是便得到其载波频率。也就是说，捕获过程就是要获得输入信号的C/A码的起始点和载波频率，然后传送给跟踪过程。 捕获与跟踪过程所用到的数据都是从原始的卫星信号经过下变频器即与中频混频之后搜集到的，其中频IF为21.25 MHz, 采样频率为5 MHz， 信号的中心频率为1.25 MHz。这组数据来源于我们设计的卫星信号接纳安装，但是用模拟仿真算法也可以产生类似的数据，以验证本章引见的卫星信号的捕获算法。 第5章 GPS卫星信号的捕获 5.2 GPS卫星信号的多普勒效应卫星信号的多普勒效应GPS卫

4、星绕地球旋转一周的时间是11 h 58 min 2.05 s， 从卫星轨道的近似半径可以求得卫星的角速度d/dt和运动速度vs： (5-1) 4265601.458 103874/ssr dvkmm sdt421.458 10/11 360058 602.05drad sdt第5章 GPS卫星信号的捕获 这里rs表卫星轨道的平均半径。我们知道，一个太阳日和一个恒星日之间相差3 min55.91 s，在这段时间里，卫星大约运转了914 km(3874 m/s235.91 s)，对应地球外表与卫星的最高点，相应的角度近似为0.045 rad(914/20 192)或2.6。假设卫星接近地平线，相应

5、的角度为0.035 rad或者2。 因此我们可以看出，对于地球外表的固定一点，在每天的同一时间里， 卫星位置大约改动22.6。 第5章 GPS卫星信号的捕获 如图5-1所示，卫星在位置S处，用户在位置A处，相对用户的卫星角速度vd呵斥了多普勒频移， 其中vd值为： (5-2) sindsvv根据卫星轨道速度，取程度方向的最大值， 得到多普勒角速度的最大值为： m/h2078s/m9295602663683874sedmaxrrvvs(5-3) 第5章 GPS卫星信号的捕获 图图5-1 5-1 卫星运动引起的多普勒频移卫星运动引起的多普勒频移 第5章 GPS卫星信号的捕获 这个值与高速航空卫星的

6、速度相等。由地面设备引起的多普勒频移通常非常小，即使它直接相对卫星运动引起了最大多普勒效应。对由C/A码调制过的频率L1(f=1575.42 MHz，最大的频移为： 81575.42 9294.93 10rdmdrf vfkHzc(5-4) 这里c是光速。 因此， 对于一固定观测器来说，最大的多普勒频移是5 kHz。 第5章 GPS卫星信号的捕获 假设地面GPS接纳设备高速挪动，就要思索进多普勒效应了。要单独产生一个5 kHz的频率范围，接纳设备需以2078 m/h的速度相对卫星运动，这个速度将包含最高速度的航空器。 因此，在设计GPS接纳机时，假设接纳机用在低速设备，那么以为载波频率的多普勒

7、频移范围在5 kHz；相反地，假设用在高速设备，就要合理假定其多普勒频移的范围在10 kHz。 这些值对于确定捕获过程的搜索频率范围是至关重要的。因此，就有了引言中我们提到的为覆盖高速卫星预期中的一切多普勒频率范围，捕获方法覆盖的频率范围必需在10 kHz之内之说。 第5章 GPS卫星信号的捕获 由于C/A码频率非常低，所以C/A码上的多普勒频移也是非常小的。C/A码的频率是1.023 MHz，比载波频率低1540(1575.42/1.023)倍。 C/A码的多普勒频率为： 681.023 109293.23 10chdcf vfHzc(5-5) 假设接纳机高速挪动，这个值就得取双倍，为6.4

8、 Hz，得到的这个值对跟踪过程BASS方法来说是至关重要的，这点我们将在下一章讨论。在BASS跟踪过程中，输入信号与本地产生数据必需求严密对齐，而C/A码的多普勒频移能够呵斥输入码与本地生成码的不对齐。 第5章 GPS卫星信号的捕获 在数字化卫星信号中，假设数据用5 MHz采样(称为采样频率)， 那么每个采样之间相隔200 ns称为采样时间。 在跟踪过程中，我们期望本地生成信号与输入信号未对齐的长度在半个采样时间或近似100 ns之内，假设两个信号之间相差超越这个长度，跟踪将失锁，即失去跟踪灵敏性。C/A码的基波时间为977.5 ns或1/1.023106 s，它挪动一个周期需破费156.3

9、ms1/6.4。对这句话的了解是：多普勒频移使C/A码的频率由1.023106变为JP21.023106+6.4，即每秒内多变化了6.4个周期，那么多变化一个周期的时间为(1/6.4) ms。 第5章 GPS卫星信号的捕获 因此，它挪动100ns的数据长度，近似破费16 ms(100156.3/977.5)。在高速导航器中，需每16 ms选择一批数据，以保证输入信号与本地生成码更好地匹配。 但是由于输入信号中有噪声，用1 ms的数据能够不能非常准确地匹配， 因此可以拓宽输入信号的调整时间，每20 ms选择一批输入信号。 对一个低速导航器来说，这个时间可拓宽到40 ms。 第5章 GPS卫星信号

10、的捕获 从以上讨论我们知道，输入信号的调整取决于采样频率。 采样频率高将缩短调整时间，由于采样时间变短，在半个周期内匹配输入码与本地生成码的时间变短。假设输入信号的强弱和跟踪灵敏度不成问题，那么输入信号的调整时间可以拓宽。然而，输入信号与本地生成信号相离必需在半个基波或488.75 ns(977.5/2)之内，这个时间看做是输入信号与本地生成信号相离允许的最大值。C/A码的多普勒频率是6.4 Hz， 调整时间可以拓宽到78.15 ms1/26.4。 第5章 GPS卫星信号的捕获 5.3 GPS卫星信号捕获的思索卫星信号捕获的思索 5.3.1 捕获时的最大电文长度捕获时的最大电文长度 C/A码长

11、码长1 ms，那么至少要用，那么至少要用1ms的电文来捕获，甚至只用的电文来捕获，甚至只用1 ms的电文来捕获时，都能够发生导航电文相位偏移。假设电文的电文来捕获时，都能够发生导航电文相位偏移。假设电文中有数据偏移，电文的下一个中有数据偏移，电文的下一个1ms将不会含有偏移。因此，为将不会含有偏移。因此，为了保证捕获电文中不含有数据偏移，需求用两组延续的电文来了保证捕获电文中不含有数据偏移，需求用两组延续的电文来捕获，捕获， 这个电文最大长度是这个电文最大长度是10 ms。假设运用两组延续的。假设运用两组延续的10 ms电文来捕获，电文来捕获， 就保证了在某一组电文中不含相位偏移。就保证了在某

12、一组电文中不含相位偏移。 第5章 GPS卫星信号的捕获 限制电文长度的第二个要素是C/A码的多普勒效应。假设理想相关峰值是1，当C/A码超前或滞后1/2码时，相关的峰值将降为0.5，对应地，幅值下降6 dB。如上节所讨论的，假定未对准的C/A码允许在半个基波0.489 s内, 基波频率是1.023 MHz，C/A码上预期的最大多普勒频移是6.4 Hz， 那么两个相差6.4 Hz的频率交换半个基波需求78 ms1/26.4， 这个数据远大于10 ms。 第5章 GPS卫星信号的捕获 5.3.2 5.3.2 捕获中的频率步长捕获中的频率步长 捕获时的另一个思索要素是捕获中剥离载波所需的频率。捕获时

13、的另一个思索要素是捕获中剥离载波所需的频率。 在前文曾提过，需被覆盖的多普勒频率范围是在前文曾提过，需被覆盖的多普勒频率范围是10 kHz10 kHz， 决决议覆盖这议覆盖这20 kHz20 kHz范围的步进频率非常重要。步进频率的大小与范围的步进频率非常重要。步进频率的大小与捕获中的电文长度严密相关。假设输入信号与本地混频信号相捕获中的电文长度严密相关。假设输入信号与本地混频信号相距距1 1个周期，它们两者没有相关性；假设两者小于个周期，它们两者没有相关性；假设两者小于1 1个周期，个周期， 那么它们有部分相关性。随意选取两个信号允许的最大频率相那么它们有部分相关性。随意选取两个信号允许的最

14、大频率相离是离是0.50.5个周期，个周期， 假设电文记录长度为假设电文记录长度为1 ms1 ms，一个，一个1 kHz1 kHz的信的信号将在号将在1 ms1 ms内变换一个周期，为了坚持内变换一个周期，为了坚持1 ms1 ms内的最大频率相离内的最大频率相离在在0.50.5个周期，步进频率须为个周期，步进频率须为1 kHz1 kHz。在这种情况下，在输入信。在这种情况下，在输入信号与卷积信号之间最大限制的频率相离是号与卷积信号之间最大限制的频率相离是500 Hz/s500 Hz/s或或0.5 0.5 Hz/msHz/ms， 且输入信号刚好在两个频率点之间。假设电文记录长且输入信号刚好在两个

15、频率点之间。假设电文记录长度为度为10 ms10 ms，100 Hz100 Hz的搜索步进频率就可以满足要求。对待这的搜索步进频率就可以满足要求。对待这个问题的一个简一方法是频率相离与电文长度成反向关系。个问题的一个简一方法是频率相离与电文长度成反向关系。 第5章 GPS卫星信号的捕获 由以上的讨论可以得出如下结论：当捕获输入电文为1 ms长时，步进频率是1kHz；当电文为10 ms时，步进频率是100 Hz。从这个简单的讨论可以明显看出，捕获中的操作执行次数与总的数据点并不是成线性比例关系，当电文长度从1 ms上升到10 ms时，电文长度上升10倍，频率点数目也上升了10倍， 捕获所需的操作

16、次数上升了不止10倍。因此，实践捕获时， 假设强调捕获的速度的话，电文长度需坚持在最小值。 执行次数的添加，取决于实践捕获所用的方法，下面将讨论详细的捕获方法。 第5章 GPS卫星信号的捕获 5.4 GPS卫星信号的捕获方法卫星信号的捕获方法5.4.1 5.4.1 传统捕获方法传统捕获方法 捕获的根本思想就是展开输入信号，找到载波频率。假设捕获的根本思想就是展开输入信号，找到载波频率。假设相位正确的相位正确的C/AC/A码与输入信号相乘，输入信号将如图码与输入信号相乘，输入信号将如图5-25-2中所示中所示变成一个延续信号。图变成一个延续信号。图5-25-2中最上端是输入信号，这个信号是中最上

17、端是输入信号，这个信号是C/AC/A码编码的射频信号相位，留意到射频信号与码编码的射频信号相位，留意到射频信号与C/AC/A码是恣意选码是恣意选择的，它们并不代表卫星传输过来的真正的信号。中间的图表择的，它们并不代表卫星传输过来的真正的信号。中间的图表示示C/AC/A码，其值在码，其值在1 1之间。之间。 最下端的图是延续信号，代表输最下端的图是延续信号，代表输入信号与入信号与C/AC/A码的乘积，其对应的频谱不再是扩频，码的乘积，其对应的频谱不再是扩频， 而是一个而是一个延续信号，这个过程有时称做输入信号的延续信号，这个过程有时称做输入信号的C/AC/A码剥离。码剥离。 第5章 GPS卫星信

18、号的捕获 图5-2 C/A编码的输入信号与本地C/A码相乘 第5章 GPS卫星信号的捕获 一旦输入变成延续信号，经过傅立叶变换可以知晓其频率。 假设输入电文长是1 ms，那么傅立叶变换的频率分辨率为1 kHz。 可以设置一个阈值，即一个门限来判别一个频率分量能否足够大。高于这个门限的最高频率分量就是我们要找的频率。 假设用5 MHz使输入信号数字化，1 ms的电文将含有5000个数据点， 一个5000点的FFT(快速傅立叶变换)将产生5000个频率分量。 然而， 5000个频率分量中只需前2500个频率分量包含有用信息， 后2500个是前2500个频率分量的复共轭。频率分辨率是1 kHz， 因

19、此， FFT覆盖的总频率范围是2.5 MHz，是采样频率的一半。 但我们感兴趣的频率范围只需20 kHz，而非2.5 MHz，因此， 为了节省计算时间， 可以用DFT(离散傅立叶)方法， 只计算相隔1 kHz的21个频率分量。 第5章 GPS卫星信号的捕获 由于输入信号中C/A码的起始点还未知，因此，必需先找到这个点。为了找到这个点，我们产生本地的一个C/A码，将其数字化为5000个点，且与输入信号点对点相乘，对此乘积的结果进展FFT或DFT，可以找到其频率。为了搜索1 ms的电文， 输入电文与本地产生数字化的C/A码必需彼此相对挪动5000次。 假设运用FFT，需操作5000次，每次操作包含

20、5000个点对点的乘积和一个5000点的FFT。由于只需2500个频率分量提供信息， 其他2500个分量提供的是多余的信息，因此，输出的5000点数据 每 个 含 有 2 5 0 0 个 频 率 分 量 ， 这 样 ， 在 频 域 共 有1.25107(50002500)个输出。 在这1.25107个输出中最高幅值假设大于预先设定的门限值，那么该分量就是期望值。在如此多的数据中寻觅最高频率分量也是非常费时的。由于在期望的20 kHz的FFT输出中，只需21个频率分量是我们感兴趣的， 总输出可以降低到105 000(500021)个。经过这种方法，可以用时间分辨率为200 ns(1/5 MHz)

21、和1 kHz的频率分辨率找出C/A码的起始点。 第5章 GPS卫星信号的捕获 假设运用10 ms的电文，由于信号只需求卷积1 ms，因此需求5000次操作。每次操作包含50 000个点对点乘积和50 000次FFT，总共有1.25108(500025 000)个输出结果。假设在预期的20 kHz范围内只思索201个频率分量，需从1 005 000(5000201)个结果中挑选，因此可以看出，从1 ms到10 ms操作的增长是非常显著的。 寻觅C/A码起始点的时间分辨率仍是200 ns， 但频率分辨率提高到了100 Hz。 第5章 GPS卫星信号的捕获 传统捕获方法是由硬件实现的，硬件主要用以执

22、行以上所讲的处置。假定输入信号用5 MHz数字化，一个可行方法是产生5000点的数字C/A码，并与输入信号点对点相乘，每200 ns执行一次5000点相乘，每200 ns对一个5000点相乘的结果进展FFT频率分析，图5-3表示了该过程。假设C/A码与输入电文匹配， FFT结果中将含有一个强分量，如前所述，这个方法将产生 1.25107(50002500)个输出结果，然而，只需在10 kHz频率范围内的输出结果中挑选，这个约束使得排序过程变简单了。 第5章 GPS卫星信号的捕获 图5-3 用C/A码和频率分析进展捕获 第5章 GPS卫星信号的捕获 完成这个操作的另一个方法是经过离散傅立叶变换(

23、DFT), 将本地产生码变换成含有C/A码和射频信号的电文。射频是复数，可用ejwt表示。本地码信号由射频与C/A码乘积获得，因此， 本地码也是复数。假定L1频率(1575.42 MHz)被转换到21.25 MHz，并用5 MHz数字化，输出频率将会是1.25 MHz; 同时假定捕获过程以1 kHz的步长搜索125010 kHz 的频率范围，且总共有21个频率分量，本地码lsi可用下式表示： 2ijf tsislC e(5-6) 第5章 GPS卫星信号的捕获 其中，s表示卫星的编号；i=1, 2, 3, , 21；Cs表示编号为s的卫星的C/A码；fi=1250-10,1250-9,1250-

24、8,1250+10 kHz。 本地信号须用5 MHz数字化，产生5000个数据点， 这21组数据代表相距1 kHz的21个频率。这些数据与输入信号有关，假设本地产生信号包含正确的C/A码和正确的频率分量，当遇到正确的C/A码相位时，输出结果将会特别高。 第5章 GPS卫星信号的捕获 图5-4 经过本地产生C/A码和射频码捕获 第5章 GPS卫星信号的捕获 5.4.2 循环相关捕获方法循环相关捕获方法 1. 循环卷积与循环相关循环卷积与循环相关 在讲述循环相关捕获方法之前，在讲述循环相关捕获方法之前， 首先简单引见执行相关操首先简单引见执行相关操作的主要数学运算。作的主要数学运算。 假设一个输入

25、信号经过一个线性时不变系统，可由卷积运假设一个输入信号经过一个线性时不变系统，可由卷积运算得到时域结果，也可由傅立叶变换得到频域结果。假设系统算得到时域结果，也可由傅立叶变换得到频域结果。假设系统脉冲呼应是脉冲呼应是h(t)，输入信号，输入信号x(t)可经过卷积产生如下的输出可经过卷积产生如下的输出y(t)： ( )() ( )( ) ()y tx thdxh td (5-7) 第5章 GPS卫星信号的捕获 经过傅立叶变换，y(t)的频域呼应为： 22( )( ) ()( )()jftjftY fxh td edtxh td edt d (5-8) 令变量t-=u，那么 222( )( )(

26、)( )( )( )( )jfufjfY fxh u edu edH fxedH f X f (5-9) 第5章 GPS卫星信号的捕获 为得到时域输出，需对Y(f)进展反傅立叶变换， 结果可写成： 1( )( )( )( )( )y tx th tFX fH f(5-10) 此处，*表示卷积，F-1表示傅立叶反变换。 可以发现一个简单的关系，时域的卷积等效于频域的傅立叶变换，或者说频域卷积运算等效于时域傅立叶变换， 这两种关系可以写成： (5-11) ( )( )( )( )x th tX f H f( )( )( ) ( )X fH fx t h t通常称其为傅立叶变换中卷积的二元性。 第5

27、章 GPS卫星信号的捕获 这个思想还可用于离散系统。与延续时域表达式不同，呼应y(n)可以表达成： 10( )( ) ()Nmy nx m h nm (5-12) 此处，x(m)表示输入信号，h(n-m)是离散时域系统呼应。留意到，该式子中，h(n-m)中的时移是循环的，由于离散操作是周期性的。对上式进展离散傅立叶变换(DFT)， 结果为： 11(2)/0011(2 () )/(2)/001(2)/0( )( ) ()( )()( )( )( )( )NNjknNnmNNjn m kNjmkNmnNjmkNmY kx m h nm ex mh nm eeH kx m eX k H k(5-13

28、) 第5章 GPS卫星信号的捕获 然而，捕获时用到的并不是卷积运算，而是相关运算， 它与卷积运算不同。x(n)与h(n)之间的相关表示为： (5-14) 10( )( ) ()Nmz nx m h nm上式与式(5-12)的独一不同点就在于h(n+m)中的m前面的符号不同。h(n)是另一信号， 而非线性系统的脉冲呼应。 z(n)的DFT变换为： 11(2)/0011(2 () )/( 2)/001( 2)/01( )( ) ()( )()( )( )( )( )NNjknNnmNNjn m kNjmkNmnNjmkNmZ kx m h nm ex mh nm eeH kx m eH k Xk(

29、5-15) 第5章 GPS卫星信号的捕获 这里X-1(k)表示离散傅立叶反变换。上式还可写成： 11(2)/100( )() ( )( )( )NNjknNnmZ kx nm h m eHk X k(5-16) 假设x(n)是实数，x(n)*=x(n)，其中*表示复共轭。 有了这个关系， Z(k)的幅值为： *( )( )( )( )( )Z kHk X kH k Xk(5-17) 利用上式的这个关系，可得到输入信号与本地产生信号的相关性。如前所述，这个式子提供了循环相关性。 第5章 GPS卫星信号的捕获 2. 循环相关捕获的方法 上面讨论了循环相关的运算法那么，利用这些运算法那么就可以执行捕

30、获了。然而实践捕获中， 输入电文不是延续到达接纳机的，因此循环相关操作适宜一组或一批电文。输入电文经过5 MHz的ADC(模拟/数字转换器)采样存储在存储器中，只需1 ms的输入电文用来寻觅C/A码的起始点，其搜索频率分辨率即步进频率是1 MHz。 第5章 GPS卫星信号的捕获 捕获输入电文， 需经过以下几个步骤： (1) 对1ms的输入电文x(n)进展快速傅立叶变换FFT， 将输入转换到频域，值为X(k)，这里， n=k=0, 1, 2, , 4999。 (2) 取X(k)的复共轭，值为X(k)*。 (3) 利用式(5-6)，产生21个本地码lsi(n)(i=1, 2, , 21)，本地码包

31、含了卫星s的C/A码和一个复射频信号的乘积， 然后用5 MHz采样本地码，本地码的频率fi相距1MHz。 (4) 对lsi(n)进展FFT，变换到频域中，值为Lsi(k)。 (5) 将X(k)*与Lsi(k)点对点相乘，结果为Rsi(k)。 第5章 GPS卫星信号的捕获 (6) 对Rsi(k)进展反傅立叶变换，变换到时域中的值为rsi(n)，得到其绝对值|rsi(n)|，总共有105 000(500021)个|rsi(n)|。 (7) 在输入电文200ns的时间分辨率和载波频率为1 kHz分辨率的条件下，|rsi(n)|最大值中的第n位和第i个载波频率芯给出了C/A码的初始点。 第5章 GPS

32、卫星信号的捕获 以上的操作可由图5-5阐明。 图5-5所示结果是时域中的， 且只给出21个本地码的其中一个。可以为输入电文与本地电文位于两个圆柱体外表，为了去匹配输入电文， 本地码要旋转5000次。换句话说，一个圆柱体相对于另一个圆柱体旋转5000次。在每一步，5000个输入电文与5000个本地电文点对点相乘，相乘结果加到一同。 包含本地码与输入码一切能够的乘积需5000步，乘积中最高幅值将被记录下。最高幅值假好像时大于预置门限值的话， 就是我们的期望值。 第5章 GPS卫星信号的捕获 图5-5 循环相关捕获表示图 第5章 GPS卫星信号的捕获 3. 循环相关捕获的改良方法 循环相关捕获的改良

33、方法与上面引见的捕获方法的独一区别就是傅立叶变换FFT的长度减小了一半。上一节的循环相关方法的第(3)步产生本地码lsi(n)，由于lsi(n)是复量，故其频谱是不对称的， 如图5-6所示。 第5章 GPS卫星信号的捕获 图5-6 本地产生信号频谱图 第5章 GPS卫星信号的捕获 从图5-6中显而易见，其信息包含在前半频谱图中，后半部分只含有非常少量的信息，从而，循环相关的捕获方法可以做如下修正： (1) 对1ms的输入电文x(n)进展快速傅立叶变换FFT， 将输入转换到频域，值为X(k)，这里，n=k=0, 1, 2, , 4999。 (2) 取X(k)的前2500(k=0, 1, 2, ,

34、 4999)个值的复共轭， 值为X(k)*。 (3) 利用式(5-6)，产生21个本地码lsi(n)(i=1, 2, , 21)， 每个lsi(n)有5000个数据点。 第5章 GPS卫星信号的捕获 (4) 对lsi(n)进展FFT，变换到频域中，值为Lsi(k)。 (5) Lsi(k)的后半部只含有少量信息， 因此，取前半部分的Lsi(k) ，与X(k)*点对点相乘， 结果为Rsi(k)(k=0, 1, 2, , 4999)。 (6) 对Rsi(k)进展反傅立叶变换，变换到时域中的值为rsi(n)， 并求出其绝对值| rsi(n) |，总共有52 500(250021)个| rsi(n) |

35、。 第5章 GPS卫星信号的捕获 (7) |rsi(n) |的最大值假设大于预置门限值的话就是我们的期望值。第i个频率就是分辨率1kHz的载波频率，第n位就是分辨率400ns的C/A码的初始点。 (8) 由于这种方法寻觅C/A码起始点的时间分辨率是400 ns，可以经过比较第n位的幅值与(n-1)和(n+1)位幅值， 将分辨率提高到200ns。 该方法中，从第(5)步到第(7)步，只需2500个数据点的操作，而非5000个。第(7)步挑选最大幅值的过程由于只用到一半的输出而简化了。因此这种方法节省了操作时间。仿真结果阐明，这种方法的信噪比略低于常规循环相关捕获方法1.1dB，这能够是由于丢弃频

36、域内后半部分信号引起的。 第5章 GPS卫星信号的捕获 5.4.3 5.4.3 延迟与累积捕获方法延迟与累积捕获方法 基于延迟与累积的捕获方法是为了去除输入信号中的频基于延迟与累积的捕获方法是为了去除输入信号中的频率信息。去除了频率信息，只需用本地率信息。去除了频率信息，只需用本地C/AC/A产生码去寻觅输入产生码去寻觅输入信号信号C/AC/A码的起始点，一旦找到其起始点，就可由码的起始点，一旦找到其起始点，就可由FFTFFT或或DFTDFT得到其频率。从实际角度看，这是个非常吸引人的方法，然得到其频率。从实际角度看，这是个非常吸引人的方法，然而，运用到实践中处置而，运用到实践中处置GPSGP

37、S信号，还需进一步研讨。这种方法信号，还需进一步研讨。这种方法讨论如下。讨论如下。 首先，我们假定输入信号首先，我们假定输入信号s(t)s(t)是复量，因此表示为：是复量，因此表示为： 2( )( )jftss tC t e(5-18) 第5章 GPS卫星信号的捕获 这里，Cs(t)表示卫星s的C/A码。 输入信号的延迟表达式为： 2()()( )jf tss tC t e (5-19) *22()2( ) ()( )()( )jftjf tjfssns t s tC t C teeC t e ( )( )()nssC tC t C t表示延迟时间，s(t)与延迟表达式的复共轭的乘积为： (5

38、-20) 这里的 (5-21) 第5章 GPS卫星信号的捕获 戈尔德码(Gold Code)与其延迟的乘积可以看做是一个“新码。这个“新码也属于戈尔德码家族。仿真结果阐明，利用其自相关和交叉相关可找到“新码的起始点。“新码的起始点与C/A码的起始点一样。式5-20有意思的地方在于它的频率自变性，即频率变化不受其他变量支配，由于f和t都是常量， 因此ej2f项是常量，因此，只需搜索“新码的起始点。虽然这种方法非常吸引人，但它要求输入信号必需是复量，由于搜集到的输入电文都是实数， 因此必需将它们转化成复数。 第5章 GPS卫星信号的捕获 为防止将实数信号转化成复数信号的额外计算， 我们可以采用如下

39、方法。 输入信号为 ( )( )sin(2)ss tC tft这里，Cs(t)表示卫星s的C/A码。输入信号的延迟表达式为： ()()sin 2()ss tC tf t(5-23) s(t)与延迟s(t-)的乘积是： ( ) ()( )()sin(2)sin 2()( )cos(2)cos2(2)2ssns t s tC t C tftf tC tfft (5-24) 第5章 GPS卫星信号的捕获 Cn(t)是式(5-21)的定义。式(5-24)含有两项：一个直流项和一个高频项。 通常，可以滤去高频。为使式(5-24)可用， |cos(2f)|必需接近单位1。实际上，很难做到这点，由于频率f未

40、知，然而，由于频率在125010 kHz之内，选择一个延迟时间来满足这个条件是完全能够的。例如，可选择21250103=，如此=0.410-6s=400 ns。输入电文用5MHz数字化，采样时间是200 ns(1/5 MHz)。 在这个条件下，|cos(2f)|=|cos()|=1。假设频率偏离10 kHz，其相关值|cos(2f)|=|cos(2126010340010-6)|=0.9997， 非常接近单位1。 第5章 GPS卫星信号的捕获 因此这个方法可以运用到实践电文中，独一的约束是不可以像式(5-20)中那样恣意选择延迟时间，但可以选择延迟不是很长的0.4s的倍数，例如=1.6 s，当

41、频率偏移10kHz时，|cos(2f)|=0.995。可以看出，假设对中心频率1250kHz用一个长延迟的话，|cos(2f)|会下降很快。 假设延迟过长，|cos(2f)|能够不再接近1。 这个方法的问题是当两个带噪声的信号相乘时，噪底会上升，搜索1ms的电文将捕获不到某个卫星。 第5章 GPS卫星信号的捕获 图 5-7 相位偏移对延迟和累积捕获方法的影响(a) 无相位偏移； (b) 相位偏移两位 第5章 GPS卫星信号的捕获 5.4.4 长记录电文的相关处置长记录电文的相关处置寻觅弱信号的方法通常是添加捕获的电文长度，这样做寻觅弱信号的方法通常是添加捕获的电文长度，这样做的益处是提高了信噪

42、比，的益处是提高了信噪比， 对此的一个简单解释是，长对此的一个简单解释是，长2 ms的的电文的电文的FFT产生产生500 Hz的频率分辨率，而的频率分辨率，而1 ms的电文的的电文的FFT产生产生1 kHz的频率分辨率。由于信号在频谱展开后带宽窄，因的频率分辨率。由于信号在频谱展开后带宽窄，因此信号强度并不会由于相对窄的频率分辨率而降低，降低频此信号强度并不会由于相对窄的频率分辨率而降低，降低频率分辨率带宽反而将噪声降低了一半，因此，信噪比上升了率分辨率带宽反而将噪声降低了一半，因此，信噪比上升了3 dB。 第5章 GPS卫星信号的捕获 长记录电文相关处置的思想是利用很少的数据点执行快速傅立叶

43、变换(FFT)。利用10ms的电文(或50 000个数据点)来阐明这个思想。电文的中心频率是1.25 MHz，输入信号将被转换成一个基带信号和一个2.5 MHz的高频带。假设将高频带用滤波器滤去，将只需处置基带信号。假定此时曾经滤去高频信号， 基带信号是带C/A码的输入的下变频信号，将这个信号与C/A码点对点相乘。假设到达了C/A码的正确相位，输出将变成延续信号，由多普勒效应引起的最大频率范围在10 kHz。既然这个信号的带宽是20 kHz，我们就可以用其2.5倍的频率来采样， 即用50 kHz采样该信号。10 ms的电文在此采样频率下只产生500个数据点。然而，我们是用5 MHz采样的， 产

44、生了50 000个数据点。我们可以每100个数据点取一次平均值作为一个新的数据点，这个取均值的过程等价于一个低通滤波器，因此， 它去除了乘法运算后的1.25 MHz延续信号中的高频分量，同样， 也去除了接纳机搜集到的信号中含有的噪声信号。 第5章 GPS卫星信号的捕获 这个方法可以用如下10ms的电文的例子来阐明： (1) 将10ms长的输入信号与本地生成1.25 MHz的复合延续信号相乘，并用5 MHz采样。由于输出最大频率是10 kHz，我们可以将此时的输出称做低频输出。频率近似为2.5 MHz的高频分量将被滤波除去，因此可以忽略它。这时，得到了50 000个数据点。 (2) 将这些输出数

45、据与10个C/A码点对点相乘，总共得到50 000个数据点。 (3) 将每100个相邻数据点取其均值，变成一个新的数据点，这个过程滤去了频率近似为2.5 MHz的高频分量。 第5章 GPS卫星信号的捕获 (4) 对新产生的500个数据点进展FFT变换，找到频率内的最高输出，这个操作产生的频率输出中，只需250个输出有用。 (5) 将本地生长码相对低频输出挪动一个数据点，然后反复步骤(3)和(4)。由于每毫秒产生一次C/A码，我们需求操作5000次，而非50 000次。 (6) 在频域中，总共得到1.250106(2505000)个输出， 大于预置门限值的最高幅值就是期望值，从这个值可以得到C/

46、A码的起始点和多普勒频移。 由此得到的频率分辨率是100 Hz。 虽然上述方法简单易懂，但循环相关的捕获方法也可以用很少的操作来到达同样的目的。 第5章 GPS卫星信号的捕获 5.4.5 5.4.5 精频估计的根本概念精频估计的根本概念 从从1 ms1 ms电文中得到的频率分辨率大约是电文中得到的频率分辨率大约是1 kHz1 kHz，对跟踪环，对跟踪环来说，这个值太粗糙了。适宜跟踪过程的频率必需在几十来说，这个值太粗糙了。适宜跟踪过程的频率必需在几十HzHz之之内。通常，跟踪环的频宽是几内。通常，跟踪环的频宽是几HzHz，用，用DFT(DFT(或或FFT)FFT)计算精细频率计算精细频率Fin

47、e FrequencyFine Frequency是不可取的。由于为了得到是不可取的。由于为了得到10 Hz10 Hz的分辨的分辨率，需率，需100 ms100 ms的电文记录，假设每毫秒有的电文记录，假设每毫秒有50005000个数据点，个数据点，100 100 msms就有就有500 000500 000个数据点，其个数据点，其FFTFFT操作将是非常费时的。另外，操作将是非常费时的。另外， 100 ms100 ms电文的相位偏移的概率也相对高。电文的相位偏移的概率也相对高。 第5章 GPS卫星信号的捕获 得到适宜跟踪环的精细频率分辨率(Fine Frequency Resolution)

48、以后简称精频的方法是经过相位关系。 一旦输入信号中剥离了C/A码，输入就将变成延续波形。假设在m时辰，1ms电文中最高频率分量是Xm(k)(k表示输入信号的频率分量), 那么可由DFT输出得到输入信号的初始相位m(k)： )(Re()(Im(tan)(1kXkXkmmm (5-25) 第5章 GPS卫星信号的捕获 Im和Re分别表示虚部和实部。让我们假定在m时辰之后很短时间的n时辰，1ms电文的DFT分量Xn(k)也是最强分量，由于输入分量在很短时间内不会迅速变化。n时辰输入信号的初始相位角和频率分量k为： )(Re()(Im(tan)(1kXkXknnn(5-26) 这两个相位角可用来计算精

49、频： )-(2)()(mnkkfmn(5-27) 第5章 GPS卫星信号的捕获 这个方程式给出了一个比从DFT得到的要准确得多的频率分辨率。为了坚持其值的独一性，n(k)-m(k)的相位差必需小于2。假设相位差是最大值2，带宽就是1/(n-m)。这里的(n-m)是两组延续电文之间的延时。 第5章 GPS卫星信号的捕获 5.4.6 5.4.6 消除精频丈量中的不确定性消除精频丈量中的不确定性虽然求解精频的根本方法都是基于式虽然求解精频的根本方法都是基于式(5-27)(5-27)，但是详细实，但是详细实现它，那么有几种略微相异的方法。假设每毫秒取现它，那么有几种略微相异的方法。假设每毫秒取DFTD

50、FT的第的第k k个个分量，频率分辨率是分量，频率分辨率是1 kHz1 kHz，精频丈量的不确定宽度也是，精频丈量的不确定宽度也是1 kHz1 kHz。 在图在图5-8(a)5-8(a)中，给出了中，给出了5 5个频率分量，它们之间相隔个频率分量，它们之间相隔1 kHz1 kHz。 假设输入信号如图假设输入信号如图5-8(b)5-8(b)所示落在了两个频率分量之间，由于所示落在了两个频率分量之间，由于系统中噪声的影响，系统中噪声的影响， 将不能确定其相位。将不能确定其相位。 第5章 GPS卫星信号的捕获 图5-8 频域内的不确定区域a 频率相离; b 输入信号 第5章 GPS卫星信号的捕获 为

51、了去除这种不确定性，有两个处理方法。第一个方法是加速DFT操作过程。假设每0.5 ms执行一次DFT操作，不确定带宽是2 kHz。在频率分辨率为1 kHz和不确定带宽为2 kHz的条件下，确定一个精细频率就不存在不确定的问题了。 但是这种方法使DFT操作添加了一倍。处理这个问题的第二个方法是， 假设输入是延续信号的话，采用一种幅值比较方案，而非加快DFT操作速度。 如图5-8(b)所示，输入信号落在两个频率槽之间。假定X(k)的幅值略高于X(k-1)，将X(k)代入式(5-26)和式(5-27)，得到精细分辨率。输入频率与X(k)的值相差接近500 Hz。 第5章 GPS卫星信号的捕获 正确的

52、结果是输入频率比X(k)低500Hz，由于噪声的存在， 我们能够错误地认定是输入频率比X(k)高了500 Hz。然而， 输入频率X(k)和X(k-1)是非常接近的，且它们比X(k-1)要强得多。因此，假设输入最高频率槽是X(k)，且相位计算在不确定区域范围内，这个区域接近于X(k)与X(k-1)的中心或者X(k)与X(k+1)的中心，这时就比较X(k-1)与X(k+1)两者之间的相位。假设X(k-1)大于X(k+1)，输入频率就是低于X(k)的;否那么，就以为输入频率高于X(k)。如此，精频的准确性由相位决议，但是频率差值的符号由最高频率分量临近的两个频率分量的幅值决议。 第5章 GPS卫星信

53、号的捕获 然而，问题比这要略为复杂一些，由于两组延续数据间能够由于导航电文引起了180的相位偏移。假设发生了这种情况，输入信号就不能被看做延续信号对待了。这种能够性使不确定宽度延迟了1ms，且限制在250Hz。频率相距250Hz对应的相位角是/2。假设频率相距+250 Hz，相位角是+/2，然而导航电文引起的相位偏移将使相位角变成+/2(+/2-)，相应地引起-250 kHz的频率变换。假设在寻觅精频时不思索相位偏移，这个结果将变成500 kHz。 第5章 GPS卫星信号的捕获 为了防止这个问题，频率不确定的最大值必需小于250Hz。 假设频率最大相差是200 Hz(这个值是实验选取的)，那么

54、相应的相位角的差是2/5，如图5-9所示。假设存在相位偏移，相位差了3/5(|(2/5)|)，这个值要比2/5大得多。 这样来看，相位差可以用来确定精频而不产生错误的频率偏移。 假设相位差大于2/5，就从结果中减去以保证频率差小于2/5。为了坚持频率在200 kHz之内，X(k)中的k值之间的最大间隔将是400 kHz。假设输入是在两个临近的k值当中， 那么输入信号跟两边的k值都相距200 kHz。第5章 GPS卫星信号的捕获 图 5-9 相位差小于2/5将不会呵斥频率误差 第5章 GPS卫星信号的捕获 图 5-10 从两个相位角得到的角度差 第5章 GPS卫星信号的捕获 讨论的最后一点是将X

55、(k)的实部和虚部转化成相位角。 通常， 相位角的范围在之间。式(5-26)和式(5-27)中的两个相位角都将用这种方式得到，两个相位角的角度差可以是02间的任何值，如图5-10所示。由于200 kHz允许的最大相位角差是2/5，这两个相位角的角度差必需等于或小于2/5。假设结果大于2/5，可以从结果中加上或减去2，其绝对值就一定小于2/5了。假设将噪声思索进去，2/5门限值可以适当放宽，如用2.3/5，意味着其差必需等于或小于2.3/5这个值。假设这样处置后的值依然大于这个门限值，就意味着两个延续毫秒数据间有相位偏移了，必需从结果中减去。当然，得到的角度值也还要经过加减2，以保证最后得到的值

56、小于门限值。 第5章 GPS卫星信号的捕获 由以上讨论可知，在某一确定卫星中需求经过下面几个步骤来找出它的C/A码起始点和载波频率： (1) 对1 ms的输入执行循环相关操作，某个确定的C/A码起始点可以从这些循环相关里找到，载波频率可以以1 kHz的分辨率得到。 (2) 找到最高频率分量X(k)，在同一个毫秒数据内的两个分量，一个比X(k)中的k值低400 kHz，一个比k值高400 kHz， 对这两个分量执行DFT操作。三个输出X(k-1), X(k), X(k+1)中的最高输出将被指定为新的X(k)，并用这个新的X(k)作为DFT分量来求解精频。 第5章 GPS卫星信号的捕获 (3) 从

57、C/A码起始点处开场选择延续几个毫秒的数据，我们随意选择为5ms，将这些数据与5组延续C/A码相乘，结果一定是一个5 ms长的延续信号。但是在恣意1 ms数据中，都有能够含有一个相位偏移。 (4) 从一切输入数据中找出Xn(k)，其中n=1, 2, 3, 4, 5， 然后由式(5-25)找出相位角。 角度差值定位为： nn1(5-28) 第5章 GPS卫星信号的捕获 (5) 角度差的绝对值必需小于门限值(2.3/5)，假设不能实现，就要从上加上或减去2。假设结果依然大于门限值，就要调整相位偏移，在上加上或减去，得到的值还要与门限值2.3/5进展比较，来决议能否要再次加减2。经过这些调整后，最终

58、的角度值就是期望值。 (6) 式(5-27)用来计算精频。由于有5ms的数据，将得到4个精频值。为了提高准确度，我们将这4个精频的平均值作为要求解的值。 第5章 GPS卫星信号的捕获 5.5 Matlab与与Simulink仿真简介仿真简介 5.5.1 5.5.1 仿真工具简介仿真工具简介 MatlabMatlab是由美国是由美国MathsworksMathsworks公司开发的大型软件。在公司开发的大型软件。在MatlabMatlab软件中，包括了两大部分：数学计算和工程仿真。其数软件中，包括了两大部分：数学计算和工程仿真。其数学计算部分提供了强大的矩阵处置和绘图功能。在工程仿真方学计算部分

59、提供了强大的矩阵处置和绘图功能。在工程仿真方面，面，MatlabMatlab提供的软件几乎遍及各个工程领域，并且不断加以提供的软件几乎遍及各个工程领域，并且不断加以完善，现已成为国际上最为流行的科学计算与工程计算软件之完善，现已成为国际上最为流行的科学计算与工程计算软件之一。一。第5章 GPS卫星信号的捕获 5.5.2 Simulink5.5.2 Simulink的运用的运用本小节将利用本小节将利用BPSKBPSK传输作为例子对传输作为例子对SimulinkSimulink进展简单引进展简单引见。但是，我们只局限在见。但是，我们只局限在SimulinkSimulink的原理及在这个仿真辅导的原

60、理及在这个仿真辅导中所要求的特点上，至于更为详尽的内容，建议读者参考中所要求的特点上，至于更为详尽的内容，建议读者参考SimulinkSimulink手册。这里给出的手册。这里给出的WindowsWindows界面的图坚持了原样。界面的图坚持了原样。 SimulinkSimulink是是MatlabMatlab软件的扩展，它是实现动态建模和仿软件的扩展，它是实现动态建模和仿真的一个软件包。它与真的一个软件包。它与Matlab Matlab 言语的主要区别在于，它与用言语的主要区别在于，它与用户接口是基于户接口是基于WindowsWindows的模型化图形输入的，从而使得用户可的模型化图形输入的