CN100516927C - 用于卫星定位系统的校准和校正系统 - Google Patents
用于卫星定位系统的校准和校正系统 Download PDFInfo
- Publication number
- CN100516927C CN100516927C CNB2003801058048A CN200380105804A CN100516927C CN 100516927 C CN100516927 C CN 100516927C CN B2003801058048 A CNB2003801058048 A CN B2003801058048A CN 200380105804 A CN200380105804 A CN 200380105804A CN 100516927 C CN100516927 C CN 100516927C
- Authority
- CN
- China
- Prior art keywords
- frequency
- error
- gps
- oscillator
- signal
- 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.)
- Expired - Lifetime
Links
- 238000012937 correction Methods 0.000 title claims description 49
- 230000001413 cellular effect Effects 0.000 claims abstract description 102
- 238000005259 measurement Methods 0.000 claims abstract description 29
- 230000010267 cellular communication Effects 0.000 claims abstract description 12
- 238000000034 method Methods 0.000 claims description 87
- 238000004891 communication Methods 0.000 claims description 27
- 230000004044 response Effects 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000007620 mathematical function Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000000052 comparative effect Effects 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims 1
- 230000008859 change Effects 0.000 description 17
- 238000010586 diagram Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 17
- 238000012545 processing Methods 0.000 description 13
- 230000015572 biosynthetic process Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 7
- 238000003786 synthesis reaction Methods 0.000 description 7
- 238000005070 sampling Methods 0.000 description 6
- 230000000875 corresponding effect Effects 0.000 description 4
- 239000000654 additive Substances 0.000 description 3
- 230000000996 additive effect Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000012417 linear regression Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000013139 quantization Methods 0.000 description 3
- YBJHBAHKTGYVGT-ZKWXMUAHSA-N (+)-Biotin Chemical compound N1C(=O)N[C@@H]2[C@H](CCCCC(=O)O)SC[C@@H]21 YBJHBAHKTGYVGT-ZKWXMUAHSA-N 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- YGANSGVIUGARFR-UHFFFAOYSA-N dipotassium dioxosilane oxo(oxoalumanyloxy)alumane oxygen(2-) Chemical compound [O--].[K+].[K+].O=[Si]=O.O=[Al]O[Al]=O YGANSGVIUGARFR-UHFFFAOYSA-N 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 229910052627 muscovite Inorganic materials 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 238000005057 refrigeration Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- FEPMHVLSLDOMQC-UHFFFAOYSA-N virginiamycin-S1 Natural products CC1OC(=O)C(C=2C=CC=CC=2)NC(=O)C2CC(=O)CCN2C(=O)C(CC=2C=CC=CC=2)N(C)C(=O)C2CCCN2C(=O)C(CC)NC(=O)C1NC(=O)C1=NC=CC=C1O FEPMHVLSLDOMQC-UHFFFAOYSA-N 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 240000004859 Gamochaeta purpurea Species 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 229910052792 caesium Inorganic materials 0.000 description 1
- TVFDJXOCXUVLDH-UHFFFAOYSA-N caesium atom Chemical compound [Cs] TVFDJXOCXUVLDH-UHFFFAOYSA-N 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000003292 diminished effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000001915 proofreading effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 125000006850 spacer group Chemical group 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Radio Relay Systems (AREA)
Abstract
为了在实际的移动定位系统中提供准确和快速的位置测量,对GPS接收机进行校准,使用锁定到外部发射信号的第一频率预测下一个时间周期中的频率误差,然后,利用GPS振荡器产生第二频率。为了预测下一个时间周期中的误差,进行多次相对时间的测量,进行误差估计,响应于所述误差估计的集合,近似出一个误差函数。然后,在下一个时间周期内使用该预测误差,校正GPS接收机。在一种实现中,公开了一种多功能便携设备,其使用多个工作在预定的理想蜂窝频率上的蜂窝基站来提供蜂窝通信,并且还使用在预定GPS频率上发射GPS信号的GPS卫星来提供定位。
Description
相关申请
本申请要求2002年12月13日提交的美国临时专利申请No.60/433,138的优先权。
发明领域
本发明涉及一种通过使用无线信号来计算移动设备的位置的装置和方法,如GPS系统。
技术背景
基于GPS的现有定位技术使用多个卫星的网络,这些卫星在已知时间在预定的GPS频率上发射信号。地面上的GPS接收机测量来自可被“看到”的天空中的每个卫星的信号的到达时间。信号的到达时间以及卫星的确切位置和从各卫星发射信号的确切时间被用于对该GPS接收机的位置进行三角测量。标准的GPS接收机包括本地振荡器,其用于接收来自卫星的GPS信号。
对来自GPS卫星的信号的获取对GPS接收机的本地振荡器中的频率变化非常敏感。很多因素会造成使GPS信号的获取变得困难。以相对低的功率发射GPS信号,并且GPS卫星处于地球轨道中。在GPS信号到达地面之前,其起始的低功率已经被大大降低了,这使得信号非常弱。结果是,如果GPS接收机的本振频率即使稍微偏离GPS频率,那么有效地接收GPS信号也将变得困难且耗时。
在很多包括GPS接收系统的通信系统中,存在有一个主本地振荡器,被称为“参考本地振荡器”。该参考本地振荡器的输出信号被依次馈送到一个或多个频率合成电路中,所述频率合成电路依次产生处于其他频率上的另外的信号,该另外的信号被提供给该系统中的各电路。作为一个示例,对GPS接收系统而言,使用其额定输出频率为16.368MHz的参考本地振荡器是常见的。该振荡器的输出通常被馈送给频率合成电路,该频率合成电路使用该参考频率来产生处于1575.42MHz附近的本振频率,随后在用于将输入GPS信号转换为接近基带的下变频电路中使用该本振频率。同样,该参考振荡器的输出可用于合成一个采样时钟,其频率通常为1.023MHz的倍数,其中,该采样时钟用作对接近基带的转换后的GPS信号进行采样的数字化电路的一部分。同样,在蜂窝电话接收机中,参考振荡器的振荡频率经常在10MHz至20MHz范围内(取决于设计),该振荡频率用于产生另外的频率,以对信号进行下变频和采样。
在一些情况下,对GPS系统而言,采用如美国专利No.5,841,396、6,421,002等中披露的频率校准方法是有利的。在一种方法中,测量蜂窝电话的本地振荡器(VCTCXO)的平均频率,然后,使用该平均频率来校准GPS接收机的振荡器的频率误差。该VCTCXO通常被频率锁定到非常稳定的所接收的蜂窝信号。另一种方法是将该GPS振荡器频率锁定到蜂窝电话的本地振荡器。如果该蜂窝电话振荡器的短期稳定性差,那么后一种方法的性能会受到影响。此种蜂窝电话振荡器的稳定性取决于许多因素,如蜂窝电话设计和平台移动。
所述校准方法有时会存在性能问题,包括:
(A)由于与蜂窝电话发射机开/关相关的加热/致冷而引起的频率“竞争(racing)”问题;
(B)还是由于与蜂窝电话发射机开/关相关的电压波动而引起的频率变化;以及
(C)由于环境影响引起的其他加热/致冷情况。
虽然在静止或以相对低的速度(如步行速度)移动的平台中,上述校准测量方法可以充分地工作,但是,在车用中所观察到的蜂窝振荡器的快速频率波动会导致校准性能的恶化。也就是说,由于蜂窝振荡器被频率锁定到网络,所以所接收的蜂窝信号的频率中的快速波动会导致很大的校准误差。例如,当车辆接近服务小区基站然后经过该基站时,这样的快速波动就会发生。在这个短时间内,所观察的来自该小区基站的直接导频信号的频率会在大约一秒或几秒内从大的正多普勒频率(如+100Hz)快速地变化到大的负多普勒频率(如-100Hz)。所接收的蜂窝载波频率中的快速变化,结合蜂窝电话发射机的开或关,会导致频率锁定校准方法的效率恶化。
发明内容
为了在实际的移动定位设备中提供准确和快速的位置测量,这里描述的一种系统通过响应于锁定到外部发射信号的第一频率和GPS振荡器产生的第二频率预测下一个时间周期(period)中的频率误差,来校准GPS接收机。具体地讲,该校准系统进行数次相对时间的测量,估计每次测量中的误差,近似出一个误差函数,并预测下一个时间周期的误差。然后,在下一个时间周期中,使用该预测误差来校正GPS接收机。
公开了一种使用外部发射信号(不是GPS信号)来校准和校正移动设备中的GPS接收机的方法和装置,其中该外部发射信号具有预定的精确载波频率。为了接收在预定GPS频率上发射的GPS信号,所述方法包括:响应于该精确载波频率,产生第一频率信号;以及在GPS接收机中产生被应用来处理GPS信号的第二频率信号。例如,第一频率可以是蜂窝接收机所接收的频率的分谐波,而第二频率可以由GPS接收机中的参考GPS振荡器直接产生,或者,可以由参考GPS振荡器得到。所述校准方法包括:估计第一时间周期中的第一和第二频率之间的误差;对于至少一个另外的时间周期,重复该误差估计步骤,以提供一个误差估计集合;响应于该误差估计集合,近似出第二频率的误差函数;以及使用该近似出的函数,预测下一个时间周期中的误差。产生一个校正信号,并且响应于该预测误差,在下一个时间周期中对该GPS接收机中的振荡器进行校正,以处理GPS信号。
通常,所述校准方法包括:在多个时间周期内,重复测量第二和第一频率的比值,并且,对于每个时间周期,将该比值与一个预定数进行比较,以分别对每个时间周期估计误差。所述误差预测步骤通常包括:通过响应于所述误差估计集合,拟合一个数学函数,如对该误差估计集合求均值、使用该误差估计集合执行数学回归方法以产生该误差估计集合的最小均方拟合,从而确定第二频率的频率误差与时间的关系。
在公开的一些方法中,所述校正步骤包括:校正GPS本地振荡器,该GPS本地振荡器提供第二频率以将在该GPS频率上的GPS信号转换为预定的中频。而在公开的其他方法中,该GPS接收机包括数字处理系统,该数字处理系统包括数字本地振荡器,并且所述校正步骤包括校正该数字本地振荡器。
在一种实现中,公开了一种多功能便携设备,其使用工作在预定理想蜂窝频率上的多个蜂窝基站来提供蜂窝通信,并且,还使用在预定GPS频率上发射GPS信号的GPS卫星来提供定位。该多功能便携设备包括:蜂窝通信系统,其包括第一本地振荡器,该第一本地振荡器产生用于对与多个蜂窝基站之间的通信进行解调的第一频率;GPS系统,其包括第二本地振荡器,该第二本地振荡器产生第二频率以处理所述GPS信号;误差预测与校正系统,用于响应于在多个先前时间周期内的第一和第二频率的测量值,预测下一个时间周期中的误差,该系统还响应于该预测误差,产生下一个时间周期的校正信号;以及该GPS系统中的一个本地振荡器,用于接收校正信号,并与此响应,在下一个时间周期内处理所述GPS信号。
附图简述
为了更充分地理解本发明,下面参考附图,对实施例进行详细描述,其中:
图1是多个蜂窝基站、GPS卫星和一个手持诸如蜂窝电话之类的移动设备的用户的透视图;
图2是一个包含有通信和定位系统的实施例中的移动设备的方框图;
图3是通信系统、定位系统以及误差预测与校准系统的一个实施例的方框图;
图4是通信系统、定位系统以及误差预测与校准系统的另一个实施例的方框图;
图5是示出在一个实施例中被执行来预测误差并产生校正信号的一系列步骤的流程图;
图6A是一个示例中蜂窝振荡器频率和GPS振荡器频率的典型频率变化的图示;
图6B是使用图6A的频率观察结果的GPS振荡器频率与蜂窝振荡器频率的比值的图示;
图6C是使用图6B的比值的误差估计的图示;
图6D是示出根据图6C中的误差估计而近似出的频率校正的图示;
图7是由接近即时的+100Hz到-100Hz频率跳变得到的线性斜率估计的图示;
图8是示出一个其中小区基站塔距公路大约为15米、高度大约为15米的示例的线性斜率估计的图示;和
图9是正态分布的随机频率数据的线性斜率估计的示意图。
发明详述
下面结合附图对本发明进行描述,在附图中,相同的数字表示相同或相似的部件。
术语和字母缩写表
在详细描述中,使用以下术语和字母缩写:
基站:发射精确载波频率的站,比如蜂窝电话网络中的基站;
校准:关于测量当前时间的误差和预测未来时间的误差的处理过程;
校正:关于对在校准处理过程中预测的误差进行补偿的方法;
CDMA:码分多址;
FSK:频移键控;
GPS:全球定位系统,虽然术语GPS经常被采用来表示美国全球定位系统,该术语在这里的含义和这里描述的所有装置的方法同样很好地适用于其他全球定位系统,如俄罗斯的Glonass系统和规划的欧洲Galileo系统。在一些系统中,用术语卫星定位系统(SPS)替换术语GPS;
GSM:全球移动通信系统;
移动设备:通常由其位置待确定的用户携带的如蜂窝电话之类的便携设备;
多功能设备:具有两种或更多操作功能(如GPS系统和蜂窝通信系统)的设备;
导频信号:一种信号,通常利用伪随机序列对其进行调制,该信号由蜂窝基站发射,用于与远程设备建立通信。虽然术语“导频”经常用于CDMA蜂窝系统的环境中,但该术语在这里同样可广泛地适用于所有其他蜂窝通信系统;
精确载波频率:一种从外部信源发射的可高度预测、高度精确的频率,精确载波频率可以直接由外部信源发射,或者从外部信源发射的信号中获取,如美国专利No.6,421,002中所披露的那样。
变量表
下表列出了这里讨论的一些变量:
变量 描述
δf: (蜂窝或GPS频率的)理想频率和实际频率之间的差值;
fc: 实际的蜂窝振荡器频率(例如第一振荡器频率);
fc0: 理想的蜂窝振荡器频率(例如精确载波频率);
fg: GPS振荡器频率(例如第二振荡器频率);
fg *: 一个实施例中由频率合成器提供的GPS混合频率;
fg0: 理想GPS振荡器频率
f1: 时间t=0时的初始GPS振荡器频率(例如16.368MHz);
f2: 时间t=0时的初始蜂窝振荡器频率(例如19.2MHz);
R: 第一和第二振荡器频率之间的测量比值;
R0: 蜂窝振荡器频率和GPS振荡器频率之间的理想比值;
S: 一个实施例中用作校正信号的缩放(scale)因子。
如同背景部分所述,来自参考本地振荡器的参考频率可用于合成(或更通常地,“获取”)在该系统中使用的其他频率。为便于描述这里的一些实施例,使用术语“本地振荡器”或“振荡器”,在一个实施例中,该术语与参考本地振荡器是同义的。但是,其他实施例可以将术语“本地振荡器”或“振荡器”关联到例如经由频率合成过程从此种参考本地振荡器获取的任何振荡信号。由于此种所获取的信号的平均频率相对于参考振荡器的频率具有一种已知的关系,所以,将本发明的方法和装置应用于这些振荡器输出的任何一个中是可以理解的。因此,在本发明描述中,术语“本地振荡器”或“振荡器”被广泛地定义成包括参考振荡器或其频率是从参考振荡器获取的任何振荡器。同样,当使用术语“振荡器频率”时,所指的是参考振荡器的频率或其频率是从参考振荡器获取的任何振荡器的频率。在一种诸如蜂窝电话收发器和GPS接收机的组合设备之类的多功能通信设备中,可能存在有两个或多个参考振荡器,例如,一个用于蜂窝电话电路,一个用于GPS电路。使用术语GPS本地振荡器和蜂窝电话本地振荡器来在不同的振荡器链之间进行区分。
还应当注意的是,参考振荡器的频率可以稍微变化,其中,为了补偿相对于理想频率的各种误差,可以按照一些方式对该频率进行控制。例如,在蜂窝电话中,将压控温度补偿晶体振荡器(VCTCXO)用作频率参考是很常见的。在这里,使用补偿过程,在补偿过程中,所接收的非常稳定的蜂窝信号用于测量和校正VCTCXO频率中的误差。此外,温度测量电路被用于补偿由于环境稳定波动而导致的频率误差。
概述
如同背景技术中所述,所接收的蜂窝振荡器频率中的快速变化会导致校准方法的效率恶化,其中所接收的蜂窝振荡器频率中的快速变化与蜂窝电话发射机的开或关相一致,在所述校准方法中,使用由蜂窝系统的本地振荡器观察的外部精确载波频率来校准GPS本地振荡器。使用这里描述的系统,通过在足够长时间周期内观察蜂窝电话和/或GPS系统中的本振频率的频率变化,可以改善对GPS振荡器的校准。如同这里所述,可以定期地执行频率比较过程,其间隔可能是0.5至1.0秒,从而产生一个观察结果集合。然后,使用如线性(或更高阶)回归等曲线拟合技术,将这些观察结果(例如10至20个最新的观察结果)用来确定频率的变化率与时间的关系。该过程可以产生一个误差预测,然后将该误差预测用于校准GPS系统中的本地振荡器,这将消除由快速频率跳变而导致的误差,比如,当车辆中的移动设备经过小区基站时会发生的+100至-100Hz的频率跳变。
这里披露的方法和装置可用于GPS振荡器频率相对于时间漂移的情形。因此,这里讨论的各种误差和校正是时间的函数,例如,作为环境温度快速改变的结果。如同这里所述,随后将多个此种测量值(至少两个)用于预测未来的GPS振荡器频率与时间的关系,因此可以用于对该振荡器的漂移进行校正。通常,在相对短的时间周期(如30秒)内,将这些漂移适当地模型化为时间的线性函数,而在更长的时间周期内,二阶曲线拟合(或更高阶的拟合)可能更合适。这里描述的对测量误差的曲线拟合可以用于对这种GPS振荡器漂移进行补偿。
描述
图1示出了一种可实现本发明的环境。在所示出的环境中,GPS接收机和蜂窝电话在一个移动单元中一起实现,但是,显而易见的是,本发明也可用于其他系统中,在所述其他系统中,被锁定到一个精确载波频率的第一振荡器的频率可用于校准第二振荡器。
图1是所示统称为10的多个蜂窝基站、统称为11的GPS卫星和一个手持移动设备14的用户12的透视图。如结合图2进行的更详细的描述,移动设备14可以包括诸如蜂窝电话之类的通信系统和诸如GPS系统之类的定位系统。用户12例如可以如图所示的步行,或者也可以在汽车中或公共交通系统中移动。
蜂窝基站10包括任意数量的蜂窝基站,该蜂窝基站用作与该移动设备进行连接的通信网络的一部分。所述蜂窝基站连接到蜂窝基础结构网络(infrastructure network)15,蜂窝基础网络15提供与多个其他通信网络之间的通信服务,所述其他通信网络比如是公共电话系统16、如因特网之类的计算机网络17以及如18所示的其他通信系统。蜂窝基础结构网络通常提供通信服务,该通信服务使蜂窝电话用户通过使用电话系统16与另一电话连接,但是,也可以使用这些蜂窝基站与其他设备进行通信和/或用于其他通信目的,如与手持个人数字助理(PDA)进行因特网连接。
在一个实施例中,蜂窝基站10是CDMA通信网络的一部分,但在其他实施例中,可以使用其他类型的通信网络,如GSM网络。在CDMA系统中,每个蜂窝基站定期地发射用于惟一标识该蜂窝基站的伪随机序列。该伪随机序列是一系列用于锁定接收机的比特。在CDMA术语中,该伪随机序列被称为“导频信号”,如同这里所使用的那样,术语“导频信号”可用于任何蜂窝系统以及CDMA系统。例如,由GSM蜂窝基站提供的各种同步信号。在预定的、非常精确的频率上,从蜂窝基站发射CDMA和GSM信号,其中所述预定的、非常精确的频率提供在这里描述的一个实施例中的精确载波频率。
GPS卫星11包括任意数量的卫星,它们用于对GPS接收机进行定位。这些卫星与基于地面的GPS通信系统19进行通信,该GPS通信系统19为GPS系统提供控制功能。这些卫星被同步来同时定期地发送无线信号。当GPS接收机检测到这些信号时,GPS系统计算所检测到的GPS信号从各自的卫星传输到该接收机所需的时间量。这假设GPS接收机具有准确的本地时钟。但是,只要同时接收到四个或更多个此种信号,该时钟也可由所接收到的GPS信号自身设定。
图2是包含通信系统和定位系统的移动设备14的一个实施例的方框图。图2示出了与天线21连接的蜂窝通信系统22。蜂窝通信系统22包括用于与蜂窝基站进行通信和/或检测来自蜂窝基站的信号的合适设备、硬件和软件,其包括本地振荡器23和将本地振荡器23锁定到蜂窝基站频率的自动频率控制(AFC)系统24。在理想情况下,该频率为fco。如上所述,虽然在AFC的控制下,本地振荡器可以被锁定到所接收的蜂窝基站频率,但由于比如移动设备的移动的原因,本地振荡器产生的实际蜂窝频率fc可能与理想的基站频率fco不同。换言之,所发射的基站频率和实际的本振频率可能没有很好地同步。蜂窝通信系统22被耦合到误差预测与校准系统30,由此允许传送诸如蜂窝本振频率fc之类的数据,该本振频率fc用于误差预测,如别处详细描述。
在一个实施例中,蜂窝通信系统22包括适合于与基站的CDMA网络进行通信的CDMA通信系统,但是,在其他实施例中,该蜂窝通信系统可以包括诸如GSM之类的另一类型的网络。
移动设备控制系统25与通信系统22和定位系统27连接。该移动设备控制系统25包括任何合适的结构,如微处理器、存储器、其他硬件、固件和软件,以为与该控制系统25连接的系统提供合适的控制功能。控制系统25连接到用户接口26,该用户接口26包括任何适合于与用户交互的组件,如键盘、用于语音通信服务的麦克风/扬声器以及诸如背亮式(backlit)LCD显示器之类的显示器。连接到定位系统27的移动设备控制系统25和用户接口26为GPS接收机提供合适的功能,如控制用户输入和显示结果。
定位系统27,在该实施例中为GPS系统,其与GPS天线28连接以接收在理想GPS频率或附近发射的GPS信号。GPS系统27包括用于接收和处理GPS信号的任何合适的硬件和软件,并执行使用合适的定位算法来确定位置所必需的任何计算。GPS系统的示例如图3和4所示,并参照图3和图4进行描述。在Norman F.Krasner的美国专利No.5,841,396、6,002,363、6,421,002中公开了GPS系统的其他示例。GPS系统27包括一个或多个本地振荡器29,所述本地振荡器29用于对来自GPS卫星的GPS信号进行下变频。
GPS系统27被耦合到误差预测与校准系统30,由此可以用于传送诸如GPS振荡器频率fg之类的数据。如同本说明书中其他地方详细描述的那样,将GPS振荡器频率fg与蜂窝振荡器频率fc一起使用,进行误差预测并产生一个校正信号,该校正信号被提供给GPS系统27。
现在参考图3,在图3中示出了定位系统27a、通信系统22a和误差预测与校准系统30a的一个实施例。通过向前面讨论的参考标记附加一个小写字母(如向“21”附加“a”构成“21a”)来指定前述部件的实现。
通信系统22a包括一个具有一个AFC电路的无线调制解调器32,该AFC电路频率锁定到所接收的通信信号,从而提供线路33上的蜂窝振荡器频率输出fc,该蜂窝振荡器频率输出fc被提供给误差预测与校准系统30a。
在该实施例中的定位系统27a包括GPS RF-IF下变频器34,该GPS RF-IF下变频器34接收GPS信号并将其转换为线路35上的处于中频fIF上的信号。将中频(IF)信号提供给模/数转换器36,然后,将该数字输出提供给GPS数字处理系统37,GPS数字处理系统37用于接收该数字信号并以一种方式对其进行处理,从而提供定位信息。
RF-IF下变频器34将所接收的GPS信号与频率合成器38提供的合成频率fg *组合。下变频器中的组合处理的结果是一个处于中频fIF上的信号。例如,如果期望的是IF为0,那么该频率合成器必须提供一个与所接收的GPS信号的频率相等的频率。但是,如果频率合成器提供的合成频率fg *不等于所接收的GPS信号的频率,那么中频fIF就不为0,而是有一定量的差别。如果中频fIF与预期的GPS频率差别足够大,那么获取GPS信号可能非常困难或者甚至不可能。
GPS本地振荡器39被耦合来向频率合成器38提供GPS本振频率fg,频率合成器38合成本振频率fg以提供合成频率fg *,随后该合成频率fg *被下变频器34利用来处理GPS信号。由于合成频率fg *确定中频fIF,并且该合成频率是由GPS振荡器频率fg合成的,所以,GPS振荡器的本振频率fg中的任何变化都会影响中频fIF。由于乘法器特征,该效果可能非常显著:典型的频率合成器可能将本振频率乘以100倍或更多。因此,本振频率fg中的任何误差将被放大来提供fg *,并且,该放大的误差会使中频fIF与其预期值明显不同。
为了观察和预测本地振荡器中的误差并相对于GPS信号来对GPS接收机进行校准,误差预测与校准系统30a通过线路40从GPS本地振荡器39接收GPS本振频率fg,并且还接收线路33上的蜂窝振荡器频率。频率比较电路41监视该两个频率,并定期将其测量值提供给后续电路,例如,该频率比较电路可以在固定时间间隔内计数这两个频率中每个频率的周期数目,并将这些数目直接提供给作为时间函数的误差的监视与预测系统42。可选地,频率比较电路41可以计数零交叉点(zero crossing)的数目,它等于周期数目的2倍,并且由于比较次数翻倍,所以准确度将改善。可以以任何合适的形式来周期性地提供所计数的数量,所述合适的形式例如是原始计数(rawcount)、该两个频率之间的比值或者频率差δf。也可以使用各种频率比较方法,比如使用频谱分析方法的哪些方法,比如快速傅立叶变换方法等。
作为时间函数的误差的监视和预测系统42被耦合到频率比较电路41,以接收这两个频率的测量值。通常,定期进行测量和提供测量值,如每半秒或每1秒。作为时间函数的误差的监视和预测系统42包括任何合适的硬件或软件,并且可以包括模拟或数字设备。例如,可以通过包括微处理器、存储器、固件和/或软件来数字地实现系统42,以提供其预期的功能。如参照图5至9进行的更详细讨论,例如,误差监视系统接收频率测量值,估计并存储每个测量的误差,然后响应于一个误差估计集合来确定一个数学函数。根据该函数,以任何合适的形式计算下一个时间间隔的误差预测值。
将该误差预测值提供给频率校正产生电路43,频率校正产生电路43产生在线路44上的校正信号。在一个实施例中,校正信号的形式为缩放因子S。将校正信号提供给GPS本地振荡器39,GPS本地振荡器39响应于此来调整输出的GPS本振频率fg,从而校正预测误差。
现在参考图4,图4例示了定位系统27和误差预测与校准系统30的另一个实施例。在图4中,很多组件与图3中的组件相似,但是,替代如图3所示将校正信号提供给GPS本地振荡器39,将线路45上的校正信号提供给数字处理系统37a中的数字振荡器46。在图4中,数字处理系统37a包括与数字振荡器46耦合的数字下变频器47,以协助GPS下变频器34解调GPS信号。数字下变频器47从模/数转换器35接收数字信号,利用数字振荡器46的频率对其进行数字化组合,并产生处于数字中频fIFD上的转换后的数字化信号。然后,GPS数字处理系统48接收转换后的数字化信号,试图获取其中所包含的GPS信号,并对其进行适当的处理。
应当注意的是,在图4的实施例中,GPS本地振荡器39没有被耦合到误差校正电路30b,并且因此没有锁定到蜂窝振荡器频率。相反,GPS本地振荡器被允许有正常的频率变化,该频率变化具有改变来自GPS下变频器的信号的中频fIF的效果。使用数字本地振荡器46和数字下变频器47,将中频fIF转换为一个新的频率fIFD,其响应该校正信号来校正中频fIF中的误差。因此,在图4的实施例中,校正信号在形式上可以不同于图3的实施例中的校正信号,因为它们用来控制不同的元件。此外,除误差校正调整之外,任一实施例中的校正信号可以包括其他合适的调整。
现在参照图5,图5例示了本发明一个实施例中的被执行来预测误差和产生校正信号的一系列步骤的流程图。参照图5的下面讨论用来提供误差预测和校正方法的概述。在其他地方还详细描述了该方法的各个方面,如参考图6A至6D、7、8和9。
在51中,开始一个观察时间间隔(周期)。
在52中,观察GPS振荡器和蜂窝本地振荡器的频率。在一个实施例中,通过比较这两个频率之间的差别,来观察GPS振荡器和本地振荡器的频率。在另一个实施例中,比较这两个频率的比值。
在53中,在观察间隔结束时,进行频率测量。例如,如果该观察涉及观察GPS振荡器信号和本地振荡器信号之间的频率比值,那么,将该间隔内观察的平均频率比值(比如通过计算)作为测量值。然后,观察过程返回到51,以开始下一个间隔并重复每个后续周期的观察和测量。通常定期进行测量,例如每0.5秒到2秒的间隔。
在54中,频率测量被用来估计其各个时间周期的频率误差。54中的频率误差估计过程包括将该周期内的频率比值与一个预定常量进行比较,例如,由于GPS和理想蜂窝频率的比值是一个预定常量R0,所以,通过考虑预定常量和这两个频率之间的观察比值之间的差值,可以估计测量误差。
在55中,在已经估计误差之后,将误差存储到一个合适的位置中,比如以数字方式实现的存储器。例如,误差监视系统42(图3和4)可以包括一个数据库,用于存储多个时间周期的误差估计的相关信息。
在56中,选择一个误差估计集合,以便为下一步骤提供数据,该下一步骤是函数近似步骤。可以预先选择该集合,使其包括多个最新的误差估计,如最新的10到20个误差估计。
在57中,使用该误差估计集合来近似一个函数。例如,可以使用以下方法来近似出一个数学函数:比如对所述集合中的所述误差估计求均值;使用所述集合中的误差估计来执行线形回归法;或者,使用所述集合中的误差估计来计算最小均方曲线。
在58中,使用该近似函数,确定下一个时间间隔的预测频率误差。例如,可以直接使用该近似函数来预测下一个时间间隔中的频率误差。
在59中,使用所预测的频率误差,产生一个在下一个时间间隔被施加到合适的本地振荡器上的校正信号。例如,在图3的实现中,校正信号被提供给GPS本地振荡器39,在图4的实现中,校正信号被提供给数字本地振荡器46。
GPS振荡器和蜂窝电话主振荡器中的频率误差
现在参考图6A、6B、6C和6D。下面对两种情况下的频率误差进行分析:对于情况1,平台移动的影响被假设为最小,并且基站发射机载波频率的精确度被假设为近乎理想,正好被蜂窝电话振荡器跟踪;情况2的不同之处在于,假设蜂窝电话振荡器中存在变化。
为便于参考,被频率锁定到所接收的蜂窝电话信号(或信号)的振荡器被称为蜂窝电话主振荡器,其频率被称为“蜂窝参考”或“蜂窝振荡器”频率,用于驱动GPS电路的振荡器被称为“GPS参考”或“GPS振荡器”。
图6A是在一个示例中在紧张情况(stressed condition)下蜂窝振荡器频率和GPS振荡器频率的典型频率变化的图示。如图6A所示,蜂窝振荡器频率fc从其初始值f2开始增加,与此同时,GPS振荡器频率fg从其初始值f1开始增加,但速率更快。这是很典型的,因为主要的频率误差经常存在于GPS振荡器而非蜂窝振荡器中。
图6B是使用图6A的频率观察结果的GPS振荡器频率与蜂窝振荡器频率的比值与时间的关系的图示。从中可以看出,该比值是这两个频率的函数,但通常由GPS振荡器频率起主导作用。尽管将该比值表示为GPS振荡器频率除以蜂窝振荡器频率,但也可以使用其倒数(即蜂窝振荡器频率除以GPS振荡器频率)。在其他地方还提供了有关频率测量和比值计算的详细讨论。
图6C是使用图6B中的比值的误差估计的图示。在图6C中,在每1秒间隔的中点执行误差估计。在其他实施例中,其他合适的点可用于误差估计。此外,在一些实施例中,可以使用频率测量和/或误差校正来计算校准校正的集合,然后在曲线拟合过程中用它来近似出一个函数。
图6D是根据图6C中的误差估计近似出的频率校正曲线的图示。下面将更详细地描述一个用于产生该频率校正曲线的过程。应当注意的是,y轴的单位是无量纲的“GPS振荡器频率校正因子”,这是提供给GPS振荡器频率的一个缩放因子,用于补偿其误差。
情况1:理想的蜂窝振荡器、非理想的GPS振荡器
假设锁定到蜂窝基站的移动振荡器(蜂窝振荡器)具有第一频率fc并且用于驱动GPS电路的移动振荡器(GPS振荡器)具有第二频率fg。假设写成:
fg=fg0+δf (1)
其中,fg0为理想的GPS频率,这是期望的GPS振荡器频率,δf被定义为理想GPS频率和实际GPS振荡器频率之间的差值。图6A示出了蜂窝振荡器频率和GPS振荡器频率的一个示例。
理想GPS频率相对于理想的蜂窝频率具有特定的比值:
R0=fg0/fc0 (2A)
如果假设蜂窝振荡器频率fc是理想的(即假设fc=fc0),则:
R0=fg0/fc (2B)
图6B示出了使用图6A中的观察结果计算的比值的一个示例。在其他实施例中,也可以采样类似的处理,其中,测量蜂窝振荡器频率与GPS参考振荡器频率的比值,即,计算值fc/fg,而非R(参见(3))。
在一个示例中,所期望的GPS振荡器频率fg0可能是16.368MHz,而fc(理想的)可能是19.2MHz,在这种情况下,R0=0.8525。图6C是通过产生比值来示出连续间隔的误差测量的图示:
R=fg/fc=(fg0+δf)/fc=R0+δf/fc (3)
由于已知R0,所以可以将其从该比值中减去(如下所述,例如通过使用频率计数型操作来测量该比值),剩下误差比值δf/fc。当然,在知道fc后,就可以确定δf。
在一个实施例中,为了校正GPS误差,可以以一个能产生频率fg0的量对GPS振荡器频率fg进行缩放,也就是说,希望找到一个缩放因子S,其为S=fg0/fg。然后,该缩放因子被计算为:
S=fg0/fg=fg0/(fg0+δf)=1/(1+δf/fg0)=(1-δf/fg0-(δf/fg0)2-...) (4)
≈(1-δf/fg0)
其中,后面的近似值遵循这样的事实:该振荡器频率中的微小误差通常非常小(在10和100PPM之间)。在一个示例中对其进行计算,如图6D所示。我们发现,出于实际应用的目的,GPS频率的修改是因子(1-δf/fg0)。应当注意的是,缩放因子校正的符号与误差比值校正相反(参见等式(3))。
情况2:非理想的蜂窝振荡器和GPS振荡器
这里,对于GPS振荡器,我们写成:
fg=fg0+δfg
对于蜂窝振荡器,我们写成:
fc=fc0+δfc
其中,fg0和fc0是理想的GPS振荡器频率和蜂窝振荡器频率,δfg和δfc是各振荡器中存在的误差。
通常,这些误差是作为时间的函数而变化。在进行如上述情况1的分析之后,得到这两个频率的比值为:
R=fg/fc=(fg0+δfg)/(fc0+δfc)=fg0/fc0(1+δfg/fg0)/(1+δfc/fc0)
≈R0(1+δfg/fg0)(1-δfc/fc0)
≈R0(1+δfg/fg0-δfc/fc0) (5)
=R0+δfg/fc0-δfcfg0/fc0 2
其中,我们忽略了更高阶项,因为δfc/fc0和δfg/fc0很小,所以这种忽略是有效的。查看(5),发现与(2)基本相同,区别在于最后一项-δfcfg0/fc0 2,因为我们假设δfc不是0,所以最后一项不为0。我们还看到,蜂窝振荡器误差相对于GPS振荡器误差的影响的比值与这两个误差的比值成正比,即与δfc/δfg成正比。因此,例如,如果蜂窝振荡器误差(用Hertz表示)为GPS振荡器误差的10%并且这两个标称振荡器频率相同,那么,蜂窝振荡器误差的影响将是由于GPS振荡器自身导致的影响的1/10。
当通过使用一个仅基于(5)的右侧的前两个量的公式来补偿GPS振荡器误差时,将会出现误差。也就是说,由于通常无法测量δfc,因此假设它是0(尽管在某些情况下,长期平均会导致该量中的短期变化非常小)。如果δfc为0,与(4)类似,对于这种情况的校正缩放因子的计算将是(1-δfg/fg0)。但是,从(5)可以看出,残留误差(减去R0之后)是+δfg/fc0-δfcfg0/fc0 2,而不是量+δfg/fc0。因此,我们应用的是不正确的缩放因子(1-δfg/fg0δfc/fc0),而不是正确的缩放因子(1-δfg/fg0)。因此,最后一项是误差项。该量的该误差项与第二项的比值为-(δfc/δfg)(fg0/fc0)。通常,后一项大约为1,因此,误差比值-(δfc/δfg)确定补偿过程之后得到的残留误差。
频率测量
在一个实施例中,使用频率计数型处理来计算频率比值,其中,将蜂窝振荡器用作“时间基准(timebase)”,将GPS振荡器用作待计数的频率。通常情况下,计数该时间基准的固定数量的周期,以产生一个时间间隔。
例如,如果蜂窝振荡器频率标称为19.2MHz,那么,计数19.2百万个周期将产生1秒的时间间隔。在此时间间隔内,计数周期的数目或GPS振荡器的零交叉点的数目。在这1秒内计数的GPS振荡器周期的数目与GPS振荡器频率成正比。实际上,对于该具体示例,在1秒间隔内计数的周期数目是以HZ计量的GPS振荡器频率的估计。因为在任何间隔内零交叉点大约为完整周期的两倍,所以对该振荡器的零交叉点进行计数能够提供改进的精确度。
上述的一种替代方法是将GPS振荡器用作时间基准,即产生计数间隔,然后计数该蜂窝振荡器的零交叉点或周期的数目。显然,该比值与GPS振荡器频率的频率倒数成正比,即与GPS振荡器频率的周期成正比。该方法也被称为倒数计数法,或倒数法,而前一种方法则被称为直接法。
当GPS振荡器频率超过蜂窝振荡器频率时,直接计数法通常优于倒数计数法,相反,当蜂窝振荡器频率超过GPS振荡器频率时,倒数计数法通常优于直接计数法。
在另一个实施例中,可以得到计算频率比值的更精确的方法,代价是更加复杂。在上述计算过程中误差的一个主要来源是由于以下事实:一个时间间隔内振荡器的计数数目是以一个周期或半个周期为单位来进行量化的,这取决于使用的方法。可以采用几种方法来降低这种时间量化的影响。在一种方法中,取代在一个时间周期T内执行单频率计数处理,在每个相邻的时段T/M内执行一次,共执行M次计数处理而不停止计数器。这M个间隔中每一个间隔的计数数目将被获取,并且可用一个线形回归过程(假设频率漂移是线形的)来对该数目集合进行拟合,以确定该频率的最小均方拟合。相对于使用长度为T的时间间隔的处理,该方法将由于量化(甚至其他噪声源)而引起的RMS误差降低一个大约等于M1/2的量。
在另一个实施例中,可以将来自待测量的频率源的时钟作为模拟信号处理,并在由该时间基准确定的时间对其进行采样(例如,在等于该时间基准的N个时钟周期的时间)。对如此获得的数据进行傅立叶分析,在与该参考时钟的频率处相关的频率上产生一条谱线。可以使用各种插入技术对该频率进行精确的估计。该方法还提供高于单间隔计数法的改善的性能。
也可以使用其他频率估计方法,但最终结果是对蜂窝振荡器和GPS振荡器的频率比值进行估计。
用于近似函数的线形回归示例
下面是近似出频率误差与时间的函数的方法的一个示例。例如,可以定期地执行频率校准处理,可能以1秒的间隔。使用线形(或甚至更高阶)回归,将误差估计集合用于确定频率变化率与时间的关系曲线,这将趋向于消除快速的频率跳变,如驾车经过蜂窝基站时发生的+100-至-100Hz的频率跳变。
如果最近的m个误差估计包括用于确定频率漂移分布图(profile)的集合,那么,在一个示例中,如下给出对线形漂移的最小均方误差拟合:
其中,yn是M个频率测量值的集合,每个测量值的方差为σ2,n为测量次数,是相对于时间特性的线形频率的估计,是与该m个频率测量值中的第一个测量值对应的频率的估计。上划线表示对应量的均值。如果存在随机RMS误差/幅度测量值σ,则(6)的测量均方误差为:
在一个多模式设备的实施例中,测量的量化步长是所提供的校正信号(例如,来自蜂窝振荡器)的1/2周期。如果该信号为19MHz,那么,1秒测量间隔的量化步长为1/38e6=2.6×10-8。如果在该步长上误差均匀分布,那么,1秒测量间隔的频率中的RMS误差只有 这对应于GPS L1频率上的12.0Hz。从(7B)可以看出,如果处理12个这样的测量值,那么RMS频率误差将大约是3.45Hz。此外,从(7A)得到频率斜率的误差为RMS误差/测量值的0.084倍,或0.084×7.6×10-9=6.35×10-10Hz/秒。这对应于GPS L1频率上的大约1.0Hz/秒。显然,使用12个测量值的校准方法可以产生较好的频率斜率测量。我们发现,该校准方法本身足以提供确定频率偏移和GPS振荡器改变率所需的频率测量精确度。在线形斜率模型不能提供精确拟合的情况下,二次斜率拟合通常是合适的。具体而言,在那些情况下,替换(6)的二次回归法会工作得更好。
当然,在一些实施例中,上述误差可能主要由蜂窝振荡器而非GPS振荡器中的误差决定。但是,对于CDMA系统,蜂窝振荡器跟踪一个蜂窝信号,该蜂窝信号基本上实现与铯标准相称的长期稳定性。在GSM系统中,蜂窝信号也非常稳定,同步于一个高质量的TCXO。通常,蜂窝振荡器中的误差的主要来源是与平台移动相关联的。
这里描述的方法的一个优点在于,正常的频率变化通常由与该车辆移动相关联的多普勒频率限制来界定。因此,下面的讨论将分析平台移动对估计的不利影响。
在一些情况下,当使用如上面讨论的线形回归方法时会出现问题。假设未校正的GPS本地振荡器表现出频率误差与时间的关系为线性。如果某人在M个时间周期的时间内观察该误差,那么,得到的校正将是一个常量加上一个时间的线形函数。因此,在时间周期M+1,校正后的GPS本地振荡器将表现出非常小的误差。所以,如果某人然后使用前面M个时段的误差(包括在时段M+1的微小误差)来执行线形回归,则会得到错误的结果。该问题的一个解决办法是,对于每个时间周期(M+1、M+2等),计算在不进行频率校正时存在的GPS本振频率。由于只需要累加GPS本振频率校正的集合,然后将它们从当前振荡器设定中减去,以得到未校正的频率,所以这是非常容易理解的。然后,将这些未校正频率应用于线形(或更高阶)的回归过程中。避免上述影响的其他方法为基于闭环校正过程的方法,其中,当前测量误差的一部分被馈送到一个滤波器,其中该滤波器的输出表示对该振荡器控制的当前校正。
由于移动引起的跟踪误差
现在结合图7、8和9,说明对在使用这里描述的误差校准方法的环境中的车辆移动进行的分析。运动应当只是车辆速度的结果;在步行速度下,如1.0米/秒,L1频率上的最大多普勒频率只有5.3Hz,并且,明显的加速是不常见的。这样小的多普勒频率对跟踪性能的影响可以忽略不计。
由于车辆靠近蜂窝塔引起的快速频率跳变
图7是示出由近似即时+100Hz到-100Hz频率跳变所导致的线形斜率估计的图示。具体而言,在图7中,从蜂窝塔观察到的多普勒频率,如同针对GPS频率所计算的那样,在一秒钟内从+100Hz变化到-100Hz。例如,如果蜂窝基站塔紧靠着一条其上汽车正以大约42英里/小时通过的公路,那么将发生这种现象。当然,对于更高速车辆,多普勒频率可能更高。但是,蜂窝塔通常离公路远得多,因此,并不经常观察到此种频率的突然改变,相反,更可能的是观察到可能在数秒内平滑改变的多普勒频率。对这种更现实的情况进行简略地分析。
假设使用前面部分中描述的处理,并且使用12个校准测量值以及线形回归拟合。根据上面的假设,如果不使用线形回归拟合,可能观察到最大值为200Hz/秒的快速频率改变。然而,在图7中示出了线形回归的平滑效果。
为了获得图7的曲线,在不存在附加的噪声或任何其他线形斜率的情况下,使用12个测量值的集合,执行频率斜率与时间的关系的连续(running)估计。从图中可以看出,斜率估计中的最大误差大约是-25.1Hz。
如果存在+100到-100Hz的跳变以及附加的线形斜率(例如由于振荡器漂移),那么,估计误差仍与图7所示相同。在一个实验中,分析25Hz频率斜率误差的性能损失,并且输出的SNR损失被发现仅仅大约为用于高灵敏度跟踪期间使用的典型设定的0.6dB。前面部分的讨论还表明,与该误差相比,斜率估计中的量化相关误差是微不足道的。
图8是示出更现实情况的线形斜率估计的图示,其中,蜂窝基站塔距离公路大约15米,其高度大约为15米,车辆行驶速度为60英里/小时(27米/秒)。这是靠近公路的较低蜂窝塔。这导致在数秒时间内从+142Hz到-142的多普勒漂移。所得到的线形斜率估计如图8所示。还假设只有来自蜂窝塔的直接导频信号用于AFC测量(即没有利用其他反射的导频信号加权)。如果有其他的导频信号,通常会改善这种情形。在单个导频信号的情况下,可以看出,误差的最差情况估计大约为34.2Hz。分析表明,这种情况的输出SNR的恶化大约为1.13dB。
由于导频音调功率变化引起的随机波动
上面结合图8分析的蜂窝塔情形被期望为多数情况下的最糟情形。其他关心的情形包括这样一些情形,在这些情形中,接收几个导频信号并将其组合,以控制AFC环路。一个或其他这样的导频信号主要体现为时间函数,导致AFC漂移。但是,可以预期的是,这种快速的漂移在严重的市区环境中是很常见的,其中车辆速度通常被期望为远低于上面所分析的。此外,作为时间函数的AFC随机波动趋向于导致能够降低频率斜率漂移的平均影响。
图9是描述正态分布的随机频率数据的线形斜率估计的示意图。假设车辆移动的影响导致AFC频率在[-50Hz,50Hz]范围内游动,每秒一个随机值。那么,对该情形进行的1000秒内的模拟产生如图9所示的频率斜率估计。使用前面12个测量值的线形回归对该模拟进行的分析,表现出4.2Hz的RMS斜率估计误差,这与(7A)对于这种情况的预测是一致的,(7A)的预测是4.18Hz。这么小的误差导致可以忽略性能损失。
3、进一步的评论和改进
可以使用上述线形回归的各种修改来进一步改进上述结果。例如,可以对校准采样进行特殊滤波。在上述一种方案中,频率与时间的关系的函数基本为一个阶跃函数,因此,如果对该阶跃函数进行求微分,结果将是一个很大的脉冲,可以经由合适的门限技术将该脉冲滤去。进一步的集成操作产生校准采样,其频率阶跃函数被移除。可以基于识别校准值与时间关系的特定“模式”来研究其他方法,其中所述特定“模式”很可能与期望的情况相关联。此外,应当明白的是,这里描述的系统通常适用于任何定位系统,如俄罗斯的Glosnass系统、规划的欧洲Galileo系统、美国的全球定位系统以及任何其他卫星定位系统。
鉴于这些教导,本领域技术人员应当认识到,可以在不偏离本发明精神或保护范围的前提下,实现其他实施例。本发明仅由下面的权利要求书进行限定,当结合上述说明书和附图考虑时,该权利要求书包括所有此种实施例和修改。
Claims (24)
1、一种用于校正GPS接收机的方法,包括:
产生锁定到蜂窝网络信号的第一频率;
在第一时间周期内,在GPS振荡器中产生第二频率;
在所述第一时间周期内,基于所述第一和第二频率的比较结果来计算第一误差;
在第二时间周期内,基于所述第一和第二频率的比较结果来计算第二误差;
根据所述第一误差和第二误差来提供误差函数;
使用所述误差函数预测下一个时间周期的GPS振荡器误差;以及
响应于所述预测的GPS振荡器误差,在所述下一个时间周期中校正所述GPS振荡器。
2、如权利要求1所述的方法,还包括在操作期间重复所述方法。
3、如权利要求1所述的方法,其中,所述提供误差函数步骤包括对所述第一和第二误差计算拟合一个数学函数。
4、如权利要求1所述的方法,其中,所述校正步骤包括:校正GPS本地振荡器,该GPS本地振荡器提供所述第二频率以将GPS信号转换为预定的中频。
5、如权利要求1所述的方法,其中,所述GPS振荡器是数字本地振荡器,并且,所述校正步骤包括校正所述数字本地振荡器。
6、一种使用具有预定精确载波频率的蜂窝网络信号来校正移动设备中的GPS接收机,以接收预定GPS频率上发射的GPS信号的方法,所述方法包括以下步骤:
响应于所述精确载波频率,产生第一频率信号;
在所述GPS接收机中产生第二频率信号,以处理所述GPS信号;
在第一时间周期内计算所述第一和第二频率之间的误差;
在至少一个不同的第二时间周期内,至少重复所述误差计算步骤一次,以提供一个误差估计集合;
直接根据所述误差估计集合,提供所述第二频率的误差函数;
使用所述误差函数,遵循所述重复步骤,预测下一个时间周期中的误差;以及
校正所述第二频率,以在所述下一个时间周期中处理所述GPS信号。
7、如权利要求6所述的方法,其中,所述计算误差步骤包括:测量每个所述时间周期内的所述第二和所述第一频率的比值,并且,对于每个时间周期,将该比值与一个预定数值进行比较。
8、如权利要求6所述的方法,其中,所述误差预测步骤包括:通过对所述误差估计集合拟合数学函数,来确定所述第二频率的频率误差与时间的关系。
9、如权利要求6所述的方法,其中,所述误差预测步骤包括以下之一:
对所述误差估计集合求均值;以及
使用所述误差估计集合执行数学回归方法,以产生所述误差估计集合的最小均方拟合。
10、如权利要求6所述的方法,其中,所述蜂窝网络信号由蜂窝基站提供,所述移动设备包括蜂窝收发器,并且还包括:
频率锁定所述精确载波频率,以提供所述第一频率;以及
在所述移动设备和所述蜂窝基站之间进行通信。
11、如权利要求6所述的方法,还包括:从所述蜂窝网络信号中提取出精确载波频率,以提供所述第一频率。
12、如权利要求6所述的方法,其中,所述校正步骤包括:校正GPS本地振荡器,该GPS本地振荡器提供所述第二频率以将所述GPS频率上的所述GPS信号转换为预定的中频。
13、如权利要求6所述的方法,其中,所述GPS接收机包括数字本地振荡器,并且所述校正步骤包括校正所述数字本地振荡器。
14、一种使用从蜂窝基站发射的通信信号来校正移动设备中的GPS接收机以接收在预定GPS频率上发射的GPS信号的方法,其中所述通信信号包括预定精确载波频率,所述方法包括以下步骤:
将第一本地振荡器同步到所述精确载波频率,以产生第一频率;
在所述GPS接收机中产生第二频率信号,第二频率信号具有标定值,所述第二频率信号用于处理所述GPS信号;
在多个时间周期内监视所述第一和第二频率,以提供所述多个时间周期中每一个时间周期的频率测量结果;
对每个监视时间周期计算所述第一和第二频率中的误差,以提供一个误差估计集合;
通过响应于所述误差估计集合拟合出一个数学函数,来提供所述第二频率的频率误差与时间关系的函数;
使用所述函数,预测下一个时间周期中的误差;
响应于所述预测误差,产生所述下一个时间周期的校正信号;以及
响应于所述校正信号,处理所述GPS信号。
15、如权利要求14所述的方法,其中,所述监视步骤包括:在多个时间周期内重复测量所述第二和所述第一频率的比值,并且,对于每个时间周期将该比值与一个预定数值进行比较,以计算所述误差。
16、如权利要求14所述的方法,其中,所述误差预测步骤包括以下之一:
对所述集合中的所述误差计算求均值;以及
使用所述集合中的所述误差计算执行数学回归方法以计算出所述误差估计计算的最小均方曲线。
17、如权利要求14所述的方法,还包括步骤:响应于所述校正信号,校正GPS本地振荡器的频率,该GPS本地振荡器提供所述第二频率以将所述GPS频率上的所述GPS信号转换为预定的中频。
18、如权利要求14所述的方法,其中,所述GPS接收机包括数字本地振荡器,并且,还包括步骤:响应于所述校正信号,校正从所述数字本地振荡器输出的频率。
19、一种多功能便携设备,用于使用在预定蜂窝频率上工作的蜂窝基站网络来提供蜂窝通信,并且,还使用在预定GPS频率上发射GPS信号的GPS卫星来提供定位,包括:
蜂窝通信系统,包括第一本地振荡器,用于产生用来对与所述蜂窝基站网络之间的通信进行解调的第一频率;
GPS系统,包括第二本地振荡器,该第二本地振荡器产生被耦合来处理所述GPS信号的第二频率;
误差预测与校正系统,用于根据表示多个先前时间周期内的所述第一和第二频率的测量结果的误差函数,预测下一个时间周期中的误差,所述系统还响应于所述预测误差,产生所述下一个时间周期的校正信号;以及
所述第二本地振荡器,用于接收所述校正信号并响应于该校正信号在所述下一个时间周期内处理所述GPS信号。
20、如权利要求19所述的多功能便携设备,其中,所述误差预测与校正系统包括:
频率比较电路;
误差监视和预测电路;以及
频率校正产生电路。
21、如权利要求20所述的多功能便携设备,其中,所述频率比较电路包括:
用于监视所述第一和第二频率并且定期地提供表示观察时期内所述第一和第二频率的测量结果的装置,以及
所述误差监视和预测电路定期地将每个所述测量结果与一个预定数值进行比较,以提供每个时间周期的误差估计。
22、如权利要求20所述的多功能便携设备,其中,所述误差函数是数学函数,所述误差监视和预测电路包括用于根据所述多个误差估计来近似出所述数学函数的系统。
23、如权利要求22所述的多功能便携设备,其中,所述近似系统包括以下之一:
用于在预定数量的周期内对所述误差估计求均值的装置;以及
用于使用预定数量时间周期的所述误差估计执行数学回归方法的装置,以计算所述多个时间周期内所述误差估计的最小均方。
24、如权利要求19所述的多功能便携设备,其中,所述第二本地振荡器包括在GPS信号的数字频率转换中使用的数字本地振荡器。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US43313802P | 2002-12-13 | 2002-12-13 | |
US60/433,138 | 2002-12-13 | ||
US10/336,125 | 2003-01-03 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009101413956A Division CN101604010B (zh) | 2002-12-13 | 2003-12-10 | 用于卫星定位系统的校准和校正系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1726405A CN1726405A (zh) | 2006-01-25 |
CN100516927C true CN100516927C (zh) | 2009-07-22 |
Family
ID=35151621
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2003801058048A Expired - Lifetime CN100516927C (zh) | 2002-12-13 | 2003-12-10 | 用于卫星定位系统的校准和校正系统 |
CN2009101413956A Expired - Lifetime CN101604010B (zh) | 2002-12-13 | 2003-12-10 | 用于卫星定位系统的校准和校正系统 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009101413956A Expired - Lifetime CN101604010B (zh) | 2002-12-13 | 2003-12-10 | 用于卫星定位系统的校准和校正系统 |
Country Status (2)
Country | Link |
---|---|
CN (2) | CN100516927C (zh) |
BR (1) | BRPI0317181B1 (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7613468B2 (en) * | 2006-02-28 | 2009-11-03 | Trimble Navigation Limited | Method and a system for communicating information to a land surveying rover located in an area without cellular coverage |
CN101059560B (zh) * | 2006-04-17 | 2011-04-20 | 中国科学院空间科学与应用研究中心 | 一种检测掩星双频gps接收机观测量测量误差的方法 |
CN102135623B (zh) * | 2011-03-23 | 2014-11-05 | 中兴通讯股份有限公司 | 移动终端和全球定位系统参数校准方法 |
US10107917B2 (en) * | 2011-12-05 | 2018-10-23 | Mediatek Inc. | Method of inter-channel bias calibration in a GNSS receiver and related device |
US8923778B2 (en) * | 2012-08-20 | 2014-12-30 | Google Technology Holdings LLC | Method for automatic frequency correction in a multi-carrier communications device |
US20150118980A1 (en) * | 2013-10-29 | 2015-04-30 | Qualcomm Incorporated | Transmitter (tx) residual sideband (rsb) and local oscillator (lo) leakage calibration using a reconfigurable tone generator (tg) and lo paths |
US9184786B2 (en) * | 2013-11-19 | 2015-11-10 | Marvell World Trade Ltd. | Systems and methods for clock calibration for satellite navigation |
CN105353602B (zh) * | 2015-11-21 | 2018-01-12 | 广西南宁至简至凡科技咨询有限公司 | 一种基于gps的校时系统 |
US10856203B2 (en) * | 2017-01-19 | 2020-12-01 | Qualcomm Incorporated | Signaling for link aggregation setup and reconfiguration |
GB2564406B (en) * | 2017-07-06 | 2022-09-07 | Focal Point Positioning Ltd | Method and system for correcting the frequency or phase of a local signal generated using a local oscillator |
CN109743710B (zh) * | 2019-01-08 | 2022-02-01 | 合肥移顺信息技术有限公司 | 提高gps授时精度的方法 |
CN112230256B (zh) * | 2019-07-15 | 2024-04-09 | 苏州宝时得电动工具有限公司 | 自主机器人及其定位校准方法、装置和存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5623414A (en) * | 1995-01-24 | 1997-04-22 | Massachusetts Inst Technology | Clock-aided satellite navigation receiver system for enhanced position estimation and integrity monitoring |
US6208290B1 (en) * | 1996-03-08 | 2001-03-27 | Snaptrack, Inc. | GPS receiver utilizing a communication link |
-
2003
- 2003-12-10 CN CNB2003801058048A patent/CN100516927C/zh not_active Expired - Lifetime
- 2003-12-10 CN CN2009101413956A patent/CN101604010B/zh not_active Expired - Lifetime
- 2003-12-10 BR BRPI0317181A patent/BRPI0317181B1/pt active IP Right Grant
Also Published As
Publication number | Publication date |
---|---|
CN1726405A (zh) | 2006-01-25 |
BR0317181A (pt) | 2005-10-25 |
CN101604010B (zh) | 2012-08-22 |
BRPI0317181B1 (pt) | 2016-10-18 |
CN101604010A (zh) | 2009-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101201974B1 (ko) | 위성 위치 확인 시스템들을 위한 교정 및 정정 시스템 | |
EP3458876B1 (en) | Position estimation in a low earth orbit satellite communications system | |
US6525688B2 (en) | Location-determination method and apparatus | |
AU665571B2 (en) | Navigation and positioning system and method using uncoordinated beacon signals | |
US6525687B2 (en) | Location-determination method and apparatus | |
CN102272625B (zh) | 使用近地轨道(leo)卫星的本地时钟频率校准 | |
CN101299886B (zh) | 使用蜂窝通信系统的位置确定系统 | |
CN100384104C (zh) | 在gps信号捕获中的振荡器频率校正 | |
CN101099089B (zh) | 用于在接收定位信号时增加相干积分长度的方法和设备 | |
CN100516927C (zh) | 用于卫星定位系统的校准和校正系统 | |
US20110109506A1 (en) | Simulating Phase-Coherent GNSS Signals | |
CN1202968A (zh) | 定位系统 | |
JP2010528495A (ja) | 移動端末において時間変化を決定する方法 | |
KR100713456B1 (ko) | 어시스트 전세계 위치 확인 시스템 수신기에서 타임 동기유지 장치 및 방법 | |
CN109782314B (zh) | Gnss卫星信号接收分级处理仿真实验平台 | |
US7068710B2 (en) | Receiver for a spread spectrum system | |
US20080123718A1 (en) | Positioning apparatus and control method thereof | |
CN106932795B (zh) | 一种gnss信号的矢量和标量混合跟踪方法及跟踪环路 | |
US6714159B1 (en) | Method for performing positioning and an electronic device | |
Kast et al. | Timing offset and timing stability for a dual-clock channel sounder | |
RU2282881C1 (ru) | Способ измерения абсолютного времени распространения коротких радиоволн в ионосфере с помощью радиосигналов с линейно-частотной модуляцией | |
Xu et al. | Design of a high fidelity signal simulator aided with Precise Point Positioning | |
CN117130025A (zh) | 一种星间链路系统在轨基准频率稳定度估计方法 | |
Kou et al. | A software-based receiver sampling frequency calibration technique and its application in GPS signal quality monitoring | |
JPH09127230A (ja) | Gps受信機の衛星信号追尾方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 1084733 Country of ref document: HK |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: GR Ref document number: 1084733 Country of ref document: HK |
|
CX01 | Expiry of patent term | ||
CX01 | Expiry of patent term |
Granted publication date: 20090722 |