[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN102081161B - 卫星信号捕获方法及卫星信号接收装置 - Google Patents

卫星信号捕获方法及卫星信号接收装置 Download PDF

Info

Publication number
CN102081161B
CN102081161B CN201010568318.1A CN201010568318A CN102081161B CN 102081161 B CN102081161 B CN 102081161B CN 201010568318 A CN201010568318 A CN 201010568318A CN 102081161 B CN102081161 B CN 102081161B
Authority
CN
China
Prior art keywords
code
differential
signal
satellite
replica
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 - Fee Related
Application number
CN201010568318.1A
Other languages
English (en)
Other versions
CN102081161A (zh
Inventor
山形整功
照内则生
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Seiko Epson Corp
Original Assignee
Seiko Epson Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Seiko Epson Corp filed Critical Seiko Epson Corp
Publication of CN102081161A publication Critical patent/CN102081161A/zh
Application granted granted Critical
Publication of CN102081161B publication Critical patent/CN102081161B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/709Correlator structure
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/30Acquisition or tracking or demodulation of signals transmitted by the system code related

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种用于缩短捕获卫星信号所需时间的卫星信号捕获方法及卫星信号接收装置。在基带处理电路部(20)的卫星信号捕获部(21)中,接收了从GPS卫星发送的GPS卫星信号的接收信号r(t)在接收信号差动电路部(211)中被差动,生成CA码的差动码信号w(t-T)。由复制码信号发生部(213)生成的复制码信号CAR(t)在复制码信号差动电路部(215)中被差动,生成差动复制码信号wR(t-T)。在相关处理部(217)中进行差动码信号w(t-T)和差动复制码信号wR(t-T)的相关处理,相关值P(Δt)输出到CPU 25。利用CPU 25进行基于相关值P(Δt)的GPS卫星信号捕获。

Description

卫星信号捕获方法及卫星信号接收装置
技术领域
本发明涉及卫星信号捕获方法及卫星信号接收装置。
背景技术
GPS(Global Positioning System,全球定位系统)作为利用定位用信号的定位系统广为人知,其被用于便携式电话机及汽车导航系统等中内置的GPS接收装置。GPS接收装置进行如下的位置计算处理:根据多个GPS卫星的位置及从各GPS卫星到接收装置之间的伪距等信息,计算表示接收装置的位置的三维坐标值和时钟误差。
GPS卫星信号是使用作为扩频调制方式早就已知的CDMA(CodeDivision Multiple Access,码分多址)方式进行扩展调制的通信信号的一种。GPS卫星发送GPS卫星信号时的频率(载波频率)是1.57542[GHz],然而,由于GPS卫星及接收装置的移动而产生的多普勒影响等,GPS接收装置实际接收到的频率并不一定与载波频率一致。因此,GPS接收装置为了从接收信号中捕获GPS卫星信号而进行频率搜索,确定捕获频率。
一般来说,捕获频率是这样确定的:使频率及码相位发生变化(频率方向及相位方向的相关运算。也称为频率搜索及相位搜索),同时进行接收了GPS卫星信号的接收信号与作为GPS卫星信号的扩展码的复制码的复制码信号的相关运算(例如专利文献1)。
现有技术文献
专利文献
专利文献1:日本特开2007-256111号公报
在频率方向的相关运算中,因GPS卫星以高速移动,所以不能忽视伴随GPS卫星移动产生的多普勒效应的影响,在GPS接收装置实际接收的频率中可能包括由多普勒频率引起的较大的误差。也会产生由GPS接收装置内的本地时钟误差(时钟误差)引起的频率误差。因此,GPS接收装置要在假定的频率误差范围内使频率发生细微变化,同时进行接收信号及复制码信号的相关运算。这样,存在如下的问题:在频率搜索上要花费很多时间,且使捕获GPS卫星信号的所需时间增大。
发明内容
本发明的目的在于提供一种用于缩短捕获卫星信号所需的时间的新的卫星信号捕获方法及卫星信号接收装置。
为了解决上述问题,本发明第一方式的卫星信号捕获方法包括:使用接收来自定位用卫星的卫星信号所得的第一接收信号至第三接收信号,生成对上述卫星信号的扩展码实施差动后的差动码;进行使用上述差动码和作为上述扩展码的复制码的扩展码复制码的规定的相关处理,取得相关值;以及根据上述相关值捕获上述卫星信号。
作为本发明其他方式,可以构建卫星信号接收装置,所述卫星信号接收装置包括:接收部,接收来自定位用卫星的卫星信号;差动部,使用通过上述接收部接收到的接收信号和使该接收信号延迟规定的延迟时间的延迟信号,生成对上述卫星信号的扩展码实施差动后的差动码;相关处理部,进行使用上述差动码和作为上述扩展码的复制码的扩展码复制码的规定的相关处理,取得相关值;以及捕获部,根据上述相关值捕获上述卫星信号。
根据本发明的第一方式等,使用接收来自定位用卫星的卫星信号所得的第一接收信号至第三接收信号,生成对卫星信号的扩展码实施差动后的差动码。进行使用生成的差动码及扩展码复制码的规定的相关处理,取得相关值,并根据该相关值捕获卫星信号。
使卫星信号的扩展码发生差动的差动码与接收信号的频率无关,可利用第一至第三这三种接收信号来产生。因此,通过进行利用差动码及扩展码复制码的相关处理,并根据该相关值捕获卫星信号,从而可以不进行频率方向的相关运算就能捕获到卫星信号,并可缩短捕获卫星信号所需要的时间。
根据本发明第一方式的卫星信号捕获方法得到本发明的第二方式,在本发明的第二方式中,也可以构成为如下的卫星信号捕获方法:所述第二接收信号是与所述第一接收信号相差时间T的信号,所述第三接收信号是与所述第一接收信号相差时间2T的信号。
根据本发明的第二方式,利用第一接收信号、与第一接收信号相差时间T的第二接收信号及与第一接收信号相差时间2T的第三接收信号来生成差动码。例如,通过使用接收信号(第一接收信号)及使第一接收信号延迟而得到的两种延迟信号(第二接收信号、第三接收信号),从而可以简单地生成使卫星信号的扩展码发生差动的差动码。
根据本发明第二方式的卫星信号捕获方法得到本发明的第三方式,在本发明的第三方式中,也可以构成为如下的卫星信号捕获方法:生成所述差动码包括:通过自乘所述第二接收信号来计算复共轭;以及把所述第一接收信号、所述复共轭的信号、所述第三接收信号相乘来计算所述差动码。
根据本发明的第三方式,通过自乘第二接收信号来计算复共轭。通过把第一接收信号、计算出的复共轭的信号、第三接收信号相乘来计算差动码。
根据本发明第二方式或第三方式的卫星信号捕获方法得到本发明的第四方式,在本发明的第四方式中,也可以构成为如下的卫星信号捕获方法:所述卫星信号捕获方法还包括:使用所述第一接收信号、与所述第一接收信号相差时间T′的第四接收信号以及与所述第一接收信号相差时间2T′的第五接收信号,生成第二差动码,其中,所述时间T′与所述时间T不同;进行使用所述第二差动码和所述扩展码复制码的规定的相关处理,取得第二相关值;以及根据所述相关值及所述第二相关值捕获所述卫星信号。
根据本发明的第四方面,使用第一接收信号、与第一接收信号相差时间T′的第四接收信号以及与第一接收信号相差时间2T′的第五接收信号,生成第二差动码。进行使用第二差动码和扩展码复制码的规定的相关处理,取得第二相关值,并根据相关值及第二相关值捕获卫星信号。使用第四接收信号、第五接收信号生成第二种差动码,并使用通过使用第二接收信号、第三接收信号而生成的差动码、及第二差动码来进行相关处理,从而可以得到互相独立的多个相关值。通过根据得到的多个相关值来捕获卫星信号,从而可更高精度地进行卫星信号的捕获。
根据本发明第一方式至第四方式中任一方式的卫星信号捕获方法得到本发明的第五方式,在本发明的第五方式中,也可以构成为如下的卫星信号捕获方法:取得所述相关值包括:对所述扩展码复制码实施差动,生成差动复制码;以及把所述差动码和所述差动复制码进行相关运算。
根据本发明的第五方式,对扩展码复制码实施差动,生成差动复制码,并把差动码和差动复制信号进行相关运算。把通过使卫星信号的扩展码发生差动而生成的差动码、与使扩展码复制码发生差动而生成的差动复制码相乘,从而可以适当得到用于捕获卫星信号的相关值。
根据本发明第二方式至第四方式中任一方面的卫星信号捕获方法得到本发明的第六方式,在本发明的第六方式中,也可以构成为如下的卫星信号捕获方法:所述扩展码复制码包括第一扩展码复制码、与所述第一扩展码复制码相差所述时间2T的第二扩展码复制码,取得所述相关值包括:通过把所述第一扩展码复制码和所述第二扩展码复制码相乘,从而生成对所述扩展码复制码实施差动后的差动复制码;以及把所述差动码和所述差动复制码进行相关运算。
根据本发明的第六方式,通过把第一扩展码复制码和与第一扩展码复制码相差时间2T的第二扩展码复制码相乘,从而生成对扩展码复制码实施差动后的差动复制码。把差动码和差动复制码进行相关运算。
根据本发明第一方式至第四方式中任一方式的卫星信号捕获方法得到本发明的第七方式,在本发明的第七方式中,也可以构成为如下的卫星信号捕获方法:取得所述相关值包括:根据所述差动码预测所述扩展码;以及把作为预测到的扩展码的预测码与所述扩展码复制码进行相关运算。
根据本发明的第七方式,根据差动码预测扩展码。把作为预测到的扩展码的预测码与扩展码复制码进行相关运算。通过根据差动码适当地预测扩展码,从而可提高对预测码及扩展码复制码的相关运算的正确性。
根据本发明第一方式至第七方式中任一方式的卫星信号捕获方法得到本发明的第八方式,在本发明的第八方式中,也可以构成为如下的卫星信号捕获方法:所述卫星信号捕获方法还包括:通过错开所述接收信号的相位并进行所述规定的相关处理来搜索相关峰值,捕获所述卫星信号是根据所述相关峰值来捕获所述卫星信号,而不需要错开频率地搜索相关峰值来捕获卫星信号。
根据本发明的第八方式,通过错开接收信号的相位并进行规定的相关处理来搜索相关峰值。根据相关峰值来捕获卫星信号,而不需要错开频率地搜索相关峰值来捕获卫星信号。通过进行使用差动码和扩展码复制码的相关处理而得到的相关值是与对接收信号进行频率方向的搜索而得到的相关值相等的相关值。因此,无需频率方向的搜索,就可以缩短现有技术中频率方向搜索中所消耗的时间。
附图说明
图1是GPS卫星信号捕获原理的说明图;
图2是GPS卫星信号捕获原理的说明图;
图3是GPS卫星信号捕获原理的说明图;
图4是GPS卫星信号捕获原理的说明图;
图5是GPS卫星信号捕获原理的说明图;
图6是GPS卫星信号捕获原理的说明图;
图7是GPS卫星信号捕获原理的说明图;
图8是GPS卫星信号捕获原理的说明图;
图9是示出便携式电话机的功能结构的框图;
图10是示出基带处理电路部的电路结构的一例的图;
图11是示出接收信号差动电路部的电路结构的一例的图;
图12是示出复制码信号差动电路部的电路结构的一例的图;
图13是示出基带处理的流程的流程图;
图14是示出码相位检测处理的流程的流程图;
图15是示出变形例中的卫星信号捕获部的电路结构的一例的图;
图16是示出变形例中的卫星信号捕获部的电路结构的一例的图;以及
图17是示出变形例中的卫星信号捕获部的电路结构的一例的图。
具体实施方式
下面,参照附图来说明本发明的优选实施方式。下面对将本发明适用于接收、捕获从GPS(Global Positioning System,全球定位系统)卫星发送来的GPS卫星信号的GPS接收装置(卫星信号接收装置)进行说明。可适用本发明的实施方式不仅限于下面的实施方式。
1.原理
首先,对GPS接收装置执行的GPS卫星信号捕获处理(卫星信号捕获方法)的原理进行说明。
GPS卫星是定位用卫星的一种,在地球外围的六个轨道面上分别配置有四颗以上的卫星,原则上其使用方式是,根据几何学的配置,在地球上的任何地方总是可以观测到四颗以上的卫星。
GPS卫星将历书或星历表等导航信息包含在GPS卫星信号中进行发送。GPS卫星信号是一种使用作为扩展码之一的CA(Course andAcquisition,粗获)码,根据作为扩频方式公知的CDMA(码分多址)方式调制的1.57542[GHz]的通信信号。CA码是将1023码片的码长作为1PN帧的、重复周期为1ms的伪随机噪声码,对每一个卫星都不同。
虽然将在GPS卫星发送GPS卫星信号时的频率(载波频率)预定为1.57542[GHz],但是受GPS卫星和GPS接收装置的移动而产生的多普勒效应的影响等,GPS接收装置在接收GPS卫星信号时的频率不一定与载波频率一致。因此,现有的GPS接收装置要进行频率搜索(使频率反复变化而进行相关运算的搜索)及相位搜索(使相位反复变化而进行相关运算的搜索),频率搜索是用于从接收信号中捕获GPS卫星信号的频率方向的相关运算,相位搜索是相位方向相关运算。由于与本实施方式特别相关的是频率搜索,因此,以下对频率搜索进行说明。
例如,受多普勒效应影响,假定接收频率以载波频率为中心在±5[kHz]范围内变化的情况下,在以载波频率即1.57542[GHz]为中心的±5[kHz]的范围内,例如每次以100[Hz]这样微小量使频率发生变化来进行频率搜索。这时,由于必须在频率的假定误差范围内对频率进行拉网式搜索,因此,产生了计算量增大且捕获GPS卫星信号所需的时间也增加的问题。
鉴于上述问题,本发明人提出了如下全新的方法:不进行上述拉网式的频率搜索,在多普勒频率未知的情况下,在短时间内捕获GPS卫星信号。在以下的说明中,使用了多种变量来说明,某变量“A”在时间“t”的值用“A(t)”表示。需要按通常方法进行相位方向的搜索。
可用下式(1)来表示GPS接收装置在时间“t”接收的接收信号“r(t)”。
【式1】
r(t)=I(t)+iQ(t)
=CA(t)·eiωt…(1)
在式(1)中,“I(t)”和“Q(t)”分别表示接收信号“r(t)”的IQ成分。即,“I(t)”表示接收信号“r(t)”的同相成分(实部),“Q(t)”表示“r(t)”的正交成分(虚部)。“CA(t)”表示GPS卫星信号的CA码,是“+1”或“-1”的任一个值。“eiωt”是表示输送GPS卫星信号的载波的项。
在式(1)中,“ω”是接收信号的频率,通过式(2)表示。
【式2】
ω=ωcd…(2)
其中,“ωc”是GPS卫星信号的载波频率,“ωd”是多普勒频率。
将式(1)的两边二乘后,可以得到式(3)。
【式3】
r(t)2=CA(t)2·e2iωt
=e2iωt…(3)
“CA(t)”是“+1”或“-1”,因此利用了“CA(t)2=1”的性质。
接着,如果计算式(3)的复共轭(使复数的虚部符号反转而得出),从而导出式(4)。
【式4】
{r(t)2}*=e-2iωt…(4)
复共轭用上标“*”表示。
另一方面,考虑到比接收信号“r(t)”滞后相当于规定时间T的信号“r(t-T)”和比接收信号“r(t)”快了相当于规定时间T的信号“r(t+T)”,如果把这些信号相乘,则导出式(5)。
【式5】
r(t-T)·r(t+T)
=CA(t-T)·eiω(t-T)·CA(t+T)·eiω(t+T)
=CA(t-T)·CA(t+T)·e2iωt…(5)
在本实施方式中,把对由1023码片构成的CA码进行采样时的采样时间间隔称为“时钟”,一时钟的经过时间用“T”表示。采样的时间间隔(即时间T)既可以是相当于一码片的时间间隔,也可以是把一码片更细化后所得的时间间隔。
用式(4)和式(5)可导出式(6)。
【式6】
w(t)=CA(t-T)·CA(t+T)
=r(t-T)·r(t+T)·e-2iωt
=r(t-T)·r(t+T)·{r(t)2}*…(6)
在式(6)中,“w(t)”是用比CA码“CA(t)”滞后一时钟的CA码“CA(t-T)”和比CA码“CA(t)”快了一时钟的CA码“CA(t+T)”的乘积表示的码。在本实施方式中,将用该“w(t)”表示的码称为CA码的“差动码”。
从式(6)可看出,CA码的差动码“w(t)”是用时间不同的多个信号“r(t-T)”、“r(t+T)”及“{r(t)2}*”的乘积来表示的。它意味着仅通过接收信号(式(6)中为r(t+T))、将接收信号延迟一时钟(时间T)而得的第一延迟信号(式(6)中为{r(t)2}*)和将接收信号延迟两时钟(时间2T)而得的第二延迟信号(式(6)中为r(t-T)),就能生成CA码的差动码的信号。
特征点在于,包括接收信号的频率“ω”的项“e-2iωt”是由接收信号的自乘的复共轭的信号“{r(t)2}*”来替代的。由此,通过进行利用差动码“w(t)”的相关处理,从而就能得到如同进行频率方向的相关运算这样的效果,且也可以不依赖接收信号的频率“ω”而捕获到GPS卫星信号。换言之,不论多普勒频率“ωd”为何值,都不用频率方向的相关运算,而仅进行相位方向的相关运算即可。
虽然以“ωd”是多普勒频率的情况进行说明,但也可包括多普勒频率以外的所有频率误差。比如,也可以包括GPS接收装置内的本地时钟误差(时钟误差)引起的频率误差。这是因为不论“ωd”为何值不进行频率搜索都能捕获到GPS卫星信号。
本发明人提出了执行利用了上述差动码“w(t)”的相关处理来捕获GPS卫星信号的两种方法。下面,对这两种方法进行详细说明。
1.1使CA码复制码发生差动的方法
第一种方法是如下的方法:使CA码复制码发生差动,生成差动复制码(differential replica code)“wR(t)”,从而进行差动码“w(t)”与差动复制码“wR(t)”的相关运算,计算相关值。
图1是该第一种方法的原理的说明图。图1示出了表示各种码的时序变化的表。最上面表示CA码“CA(t)”,第二行表示CA码的差动码“w(t)”,第三行表示CA码复制码“CAR(t)”,第四行表示差动复制码“wR(t)”,最下面表示CA码的差动码“w(t)”与CA码复制码的差动码“wR(t)”的乘积值“g(t)”。
首先,利用时刻“t-T”的CA码复制码“CAR(t-T)”和时刻“t+T”的CA码复制码“CAR(t+T)”,根据式(7)算出时间“t”的差动复制码“wR(t)”。
【式7】
wR(t)=CAR(t-T)·CAR(t+T)…(7)
错开差动复制码“wR(t)”的相位“Δt”,同时进行差动码“w(t)”和差动复制码“wR(t+Δt)”的相关运算。在本实施方式中,CA码复制码的相位用“Δt”表示。准确地说,CA码复制码的相位换算成时间的值是“Δt”。
具体地说,根据式(8)算出相关值“P(Δt)”。
【式8】
P ( Δt ) = ( Σ i g ( t ) ) 2 = ( Σ i w ( t ) · w R ( t + Δt ) ) 2
= ( Σ i CA ( t - T ) · CA ( t + T ) · CA R ( t - T + Δt ) · CA R ( t + T + Δt ) ) 2 . . . ( 8 )
“g(t)”是时刻“t”的差动码“w(t)”与差动复制码“wR(t+Δt)”的乘积值。算出的相位“Δt”的相关值用“P(Δt)”表示。
在使差动复制码“wR(t)”的相位“Δt”变化的同时,算出相关值“P(Δt)”。检测相关值“P(Δt)”变为最大的相位“Δt”。把检测出的相位“Δt”确定为该GPS卫星信号的码相位。
参照图2和图3,对第一种方法的具体实例进行说明。在此,以将接收信号的CA码的时序变化表示为“CA(t)=1,-1,-1,1,-1,1,…”的情况为例进行说明。作为捕获对象的GPS卫星信号的CA码是已知的,CA码复制码的时序变化也用“CAR(t)=1,-1,-1,1,-1,1,…”来表示,进行说明。
(i)CA码复制码的相位为“Δt=0”的情况
图2示出以CA码复制码的相位设定为“Δt=0”时的计算结果的一个例子。如果用“CA(t-1)”及“CA(t+1)”,根据式(6)来算出差动码“w(t)”,则可得到“w(t)=(N/A),-1,-1,1,1,(N/A),…”。“(N/A)”部分因为根据图来看前后的CA码的值不明,所以其是未记载的部分。
如果用“CAR(t-1)”及“CAR(t+1)”,根据式(7)算出差动复制码“wR(t)”,则可得到“wR(t)=(N/A),-1,-1,1,1,(N/A),…”。
接着,通过把CA码的差动码“w(t)”与差动复制码“wR(t)”相乘,从而得到乘积值“g(t)=(N/A),1,1,1,1,(N/A),…”。这时,比如只取出“g(2)”~“g(5)”这部分来计算相关值“P(Δt=0)”,则相关值变为“P(Δt=0)=(1+1+1+1)2=16”,可得到大的值。这是因为,从图2可知,CA码“CA(t)”的相位和CA码复制码“CAR(t-1)”的相位正好完全一致,乘积值“g(t)”全部都是“1”。
(ii)CA码复制码相位为“Δt=1”时
图3示出了CA码复制码的相位设定为“Δt=1”时的计算结果。与(i)的情况相同,如果对差动复制码“wR(t)”进行计算,则可得到“wR(t)=(N/A),1,-1,-1,1,(N/A),…”。这时,比如如果只着眼于“g(2)”~“g(5)”部分,并计算相关值“P(Δt=1)”,则相关值变为“P(Δt=1)=(-1+1-1+1)2=0”,则可得到小的值。如从图3可知的那样,这是因为由于CA码“CA(t)”的相位与CA码复制码“CAR(t-1)”的相位偏离,所以乘积值“g(t)”反复“-1”和“1”,计算和时互相抵消了。
这里,虽然只图示了在(i)Δt=0及(ii)Δt=1这两个情况下的结果,但无论使相位“Δt”怎样变化地进行计算,都不会比(i)Δt=0时的相关值“P(Δt=0)”大。如果相位不是正好一致,则这是因为乘积值“g(t)”中混有“-1”和“1”,互相抵消了的缘故。由此,在使相位“Δt”变化的同时,计算相关值“P(Δt)”,通过检测相关值“P(Δt)”变为最大的相位“Δt”,从而能正确地检测码相位。
1-2.根据差动码“w(t)”预测原CA码的预测方法
第二种方法是如下的方法:根据差动码“w(t)”来预测(推测)接收信号的CA码“CA(t)”,并进行使用预测到的CA码“CAP(t)”(预测CA码)的相关运算,从而计算相关值。
图4是该第二种方法的说明图。在图4的表中,最上面表示接收信号的CA码“CA(t)”,第二行表示CA码的差动码“w(t)”。第三行表示预测CA码“CAP(t)”,第四行表示CA码的复制码即CA码复制码“CAR(t)”。最下面表示预测CA码与CA码复制码的乘积值“f(t)=CAP(t)··CAR(t)”。
首先,进行利用了差动码“w(t)”的时序变化的预测处理来预测接收信号的CA码。具体地说,设定{(CAP(1),CAP(2))=(1,1),(1,-1),(-1,1),(-1,-1)}这四个组合作为时刻“1”的预测CA码“CAP(1)”与时刻“2”的预测CA码“CAP(2)”组合。
对于设定的四个组合中的每个组合,分别从时刻“1”开始依次进行把预测CA码“CAP(t)”与差动码“w(t+1)=CA(t)×CA(t+2)”相乘而获得的值作为预测CA码“CAP(t+2)”的运算,对各时刻“t”的预测CA码“CAP(t)”进行运算。具体地说,按“CAP(3)=CAP(1)··w(2)”、“CAP(4)=CAP(2)··w(3)”、“CAP(5)=CAP(3)··w(4)”、…这样的顺序依次对预测CA码“CAP(t)”进行运算。
根据运算式“CAP(t+2)=CAP(t)··w(t+1)=CAP(t)··CA(t)··CA(t+2)”来进行CA码的预测。它的根据如下:CA码的预测正确时,“CAP(t)=CA(t)”、“CAP(t)×CA(t)=1”成立,因此上述的运算式变为“CAP(t+2)=CA(t+2)”,CA码与预测CA码相等。
这时,对四个组合分别运算预测CA码“CAP(t)”,则对应四个组合中的任一个组合而运算出的预测CA码“CAP(t)”与接收信号的CA码“CA(t)”正好一致。然而,因为CA码“CA(t)”不明,所以无法得知哪个组合是正确的。
因此,对对应四个各组合而运算出的预测CA码“CAP(t)”进行与CA码复制码“CAR(t)”之间的相关运算。即,对预测CA码“CAP(t)”,在使CA码复制码“CAR(t)”的相位“Δt”发生变化的同时,进行相关运算,从而计算出相关值。
具体而言,根据下式(9)计算相关值“P(Δt)”。
【式9】
P ( Δt ) = ( Σ i f ( t ) ) 2 = ( Σ i CA P ( t ) · CA R ( t + Δt ) ) 2 . . . ( 9 )
在式(9)中,“f(t)”是把预测CA码“CA(t)”与CA码复制码“CAR(t+Δt)”相乘得到的乘积值。
如果在使CA码复制码“CAR(t)”的相位“Δt”发生变化的同时,计算出相关值“P(Δt)”,则可检测出相关值“P(Δt)”变为最大时的相位“Δt”。把检测出的相位“Δt”确定为该GPS卫星信号的码相位。
参照图5~图8来说明该第二种方法的具体实例。这里,以通过“CA(t)=1,-1,-1,1,-1,1,…”表示接收信号的CA码的时序变化的情况为例来说明。为使说明简单,以CA码复制码“CAR(t)”的相位“Δt=0”,通过“CAR(t)=1,-1,-1,1,-1,1,…”表示CA码复制码的时序变化的情况为例进行说明。
(A)(CAP(1),CAP(2))=(1,1)的组合
图5示出了用(CAP(1),CAP(2))=(1,1)的组合的计算结果。首先,根据接收信号“r(t)”,按照式(6)计算出差动码“w(t)”。结果,可以得到“w(t)=(N/A),-1,-1,1,1,(N/A),…”。
其次,用“(CAP(1),CAP(2))=(1,1)”计算出预测CA码“CAP(t)”。具体地说,如“CAP(3)=CAP(1)×w(2)=-1”、“CAP(4)=CAP(2)×w(3)=-1”、“CAP(5)=CAP(3)×w(4)=-1”、“CAP(6)=CAP(4)×w(5)=-1”这样,依次计算预测CA码“CAP(t)”。这样,可得到“CAP(t)=1,1,-1,-1,-1,-1,…”。
接着,进行预测CA码“CAP(t)”和“CAR(t)”的相关运算。具体地说,如“f(1)=CAP(1)×CAR(1)=1”、“f(2)=CAP(2)×CAR(2)=-1”、“f(3)=CAP(3)×CAR(3)=1”、“f(4)=CAP(4)×CAR(4)=-1”、“f(5)=CAP(5)×CAR(5)=1”、“f(6)=CAP(6)×CAR(6)=-1”这样,依次计算乘积值“f(t)”。由此,可得到“f(t)=1,-1,1,-1,1,-1,…”。
利用这个结果,根据式(9)计算出相关值“P(Δt=0)”。这时,例如只取出“f(1)”~“f(6)”部分,试算一下其相关值“P(Δt=0)”,则可得到相关值“P(Δt=0)=(1-1+1-1+1-1)2=0”。由于乘积值“f(t)”相互反复“1”和“-1”,所以值可以相抵,相关值变成小值。
相关值“P(Δt=0)”变成小值表示未正确进行CA码“CAP(t)”的预测。实际上,在图5中,若把预测CA码“CAP(t)”和实际CA码“CA(t)”比较,则可知没有正确预测CA码。
(B)(CAP(1),CAP(2))=(1,-1)的组合
图6示出了(CAP(1),CAP(2))=(1,-1)的组合的计算结果。用“(CAP(1),CAP(2))=(1,-1)”的组合,和(A)情况同样地,计算出预测CA码“CAP(t)”。由此,可得到“CAP(t)=1,-1,-1,1,-1,1,…”。
接着,如果用预测CA码“CAP(t)”和CA码复制码“CAR(t)”计算出乘积值“f(t)”,则可以得到“f(t)=1,1,1,1,1,1,…”。这时,例如也只取出“f(1)”~“f(6)”部分,试计算一下相关值P(Δt=0),则相关值为“P(Δt=0)=(1+1+1+1+1+1)2=36”。这时,因乘积值“f(t)”全部都为“1”,所以相关值“P(Δt=0)”变为大值。
相关值“P(Δt=0)”变为大值表示正确进行了CA码“CAP(t)”的预测。实际上,在图6中,若把预测CA码“CAP(t)”和实际CA码“CA(t)”比较,则可知正确预测了CA码。
(C)(CAP(1),CAP(2))=(-1,1)的组合
图7示出了(CAP(1),CAP(2))=(-1,1)的组合的计算结果。如果用(CAP(1),CAP(2))=(-1,1)的组合计算出预测CA码“CAP(t)”,则可以得到“CAP(t)=-1,1,1,-1,1,-1,…”。
接着,如果用CA码“CAR(t)”和CA码复制码“CAR(t)”计算乘积值“f(t)”,则可以得到“f(t)=-1,-1,-1,-1,-1,-1,…”。这时,例如也只取出“f(1)”~“f(6)”部分,试计算一下相关值,则变为“P(Δt=0)=(-1-1-1-1-1-1)2=36”。这时,因乘积值“f(t)”全部都为“-1”,所以相关值“P(Δt=0)”变为大值。
相关值变为大值表示正确进行了CA码“CAP(t)”的预测。在图7中,若比较CA码“CA(t)”和预测CA码“CAP(t)”,可知符号正好变为反转的关系。可以获知与GPS卫星信号重叠的导航信息每隔20毫秒极性可反转一次。因此,接收信号的CA码可能因为导航信息的极性反转而每隔20毫秒发生一次符号反转。因此,如果CA码“CA(t)”和预测CA码“CAP(t)”之间的差异是符号的反转的差异,则可以适当地认为正确地进行了CA码的预测。
(D)(CAP(1),CAP(2))=(-1,-1)的组合
图8示出了(CAP(1),CAP(2))=(-1,-1)的组合的计算结果。若用(CAP(1),CAP(2))=(-1,-1)的组合计算出预测CA码“CAP(t)”,可得到“CAP(t)=-1,-1,1,1,1,1,…”。
接着,如果用预测CA码“CAP(t)”和CA码复制码“CAR(t)”计算出乘积值“f(t)”,则得到“f(t)=-1,1,-1,1,-1,1,…”。这时,例如也只取出“f(1)”~“f(6)”部分来计算出相关值,则变为“P(Δt=0)=(-1+1-1+1-1+1)2=0”。在这种情况下,由于乘积值“f(t)”相互反复“-1”和“1”,所以计算出相关值“P(Δt=0)”变为小值。
这表示没有正确地进行使用差动码“w(t)”的CA码“CAP(t)”的预测。实际上,在图8中,如果比较预测CA码“CAP(t)”和实际CA码“CA(t)”,则可知没有正确地进行CA码的预测。
根据以上的结果,在(B)的情况或(C)的情况中,可以通过检测出相关值变为最大的相位来正确地检测出码相位。因此,对于(A)~(D)的各组合,在使CA码复制码的相位“Δt”变化的同时,计算相关值“P(Δt)”,进行相关值“P(Δt)”变为最大的相位“Δt”的检测处理。
根据上述结果可知,正是由于有符号的差异,(A)组合和(D)组合才赋予相同的结果、(B)组合和(C)组合才赋予相同的结果。因此,无需对四个组合全部都计算。实际上,只计算(A)与(B)这两个组合和(C)与(D)这两个组合中的任一个即可。
2.实施例
下面,就适用以上所述原理的GPS接收装置的实施例进行说明。这里,以安装了GPS接收装置的一种电子设备、即便携式电话机1作为具体例进行说明。
2-1.第一实施例
第一实施例是适用了以上所述原理的“1-1.使CA码复制码发生差动的方法”的实施例。
(1)结构
图9是示出在第一实施例中的便携式电话机1的功能结构的框图。便携式电话机1构成为包括:GPS天线9、GPS接收部10、主CPU(CentralProcessing Unit,中央处理器)30、操作部40、显示部50、便携式电话用天线60、便携式电话用无线通信电路部70和存储部80。
GPS天线9是用于接收包括从GPS卫星发送的GPS卫星信号的RF(Radio Frequency,射频)信号的天线,并向GPS接收部10输出所接收的信号。
GPS接收部10是根据从GPS天线9输出的信号计测便携式电话机1的位置的位置计算电路,其是相当于GPS接收装置的功能模块。GPS接收部10构成为包括RF(Radio Frequency)接收电路部11和基带处理电路部20。RF接收电路部11和基带处理电路部20可分别作为LSI(LargeScale Integration,大规模集成电路)来制造,也可作为一个码片来制造。
射频接收电路部11(相当于接收部)是RF信号的处理电路模块,通过把规定的振荡信号分频或者递增,从而生成RF信号乘算用的振荡信号。通过把生成的振荡信号乘以从GPS天线9输出的RF信号,从而RF信号降频变频为中频信号(IF(Intermediate Frequency)信号),再把IF信号进行放大等之后,用A/D转换器转换成数字信号,向基带处理电路20输出。下面说的所谓“接收GPS卫星信号”是指从用GPS天线9接收射频信号的处理开始直至生成已被降频变频的IF信号的处理为止。所谓“接收信号”表示IF信号。
基带处理电路部20是对从RF接收电路部11输出的IF信号进行相关运算处理等,并捕获、提取GPS卫星信号,对数据进行译码,从而取出导航信息、时间信息等的电路部。
图10示出了基带处理电路部20的电路结构的一例。基带处理电路部20由卫星信号捕获部21、CPU 25和存储部27构成。
卫星信号捕获部21是从RF接收电路部11输出的IF信号即接收信号中捕获GPS卫星信号的电路部,其构成为包括:接收信号差动电路部211(相当于差动部)、复制码信号发生部213、复制码信号差动电路部215和相关处理部217。
图11示出了接收信号差动电路部211的电路结构的一例。接收信号差动电路部211是使从射频接收电路部11输入的接收信号“r(t)”发生差动以生成差动码信号“w(t-T)”的电路部,是实现式(6)的电路部。但是,所谓差动码信号,是相当于以原理说明的差动码的信号。
在接收信号差动电路部211中,以图中的“T”表示的模块是延迟电路部DC,它使输入的信号仅延迟时间“T”来输出延迟信号。复共轭运算部PS是运算输入的信号的复共轭的电路部,例如包括移相电路。把接收信号“r(t)”(第一接收信号)、使接收信号延迟相当于时间T的第一延迟信号(第二接收信号)的自乘的复共轭的信号“{r(t-T)2}”和使接收信号延迟相当于时间2T的第二延迟信号(第三接收信号)“r(t-2T)”在乘算部M中进行乘算,从而生成差动码信号“w(t-T)=r(t-2T)··r(t)··{r(t-T)2}”,将其向相关处理部217输出。
复制码信号发生部213是生成GPS卫星信号的CA码的扩展码复制码的发生信号、即复制码信号的电路部。复制码信号发生部213根据从CPU 25输出的CA码指示信号(捕获对象卫星的指示信号)生成复制码信号“CAR(t)”,将其向复制码信号差动电路部215输出。
图12示出了复制码信号差动电路部215的电路结构的一例。复制码信号差动电路部215是使从复制码信号发生部213输入的复制码信号“CAR(t)”发生差动,从而生成差动复制码信号“wR(t-T)”的电路部,其是实现式(7)的电路部。但是,所谓差动复制码信号是相当于以原理说明的差动复制码的信号。
在复制码信号差动电路部215中,以图中的“T”表示的模块是延迟电路部DC。与图11的延迟电路部DC一样,使信号延迟相当于时间T。在乘法计算部M中,把输入的复制码信号“CAR(t)”与使复制码信号延迟相当于时间2T的信号“CAR(t-2T)”相乘,从而生成差动复制码信号“wR(t-T)=CAR(t-2T)··CAR(t)”,将其向相关处理部217输出。
回到图10的说明。相关处理部217是把从接收信号差动电路部211输入的差动码信号“w(t)”和从复制码信号差动电路部215输入的差动复制码信号“wR(t)”进行相关处理的电路部。相关处理部217根据从CPU 25输入的相位指示信号,在使差动复制码信号的相位“Δt”发生变化的同时,计算“w(t)”和“wR(t+Δt)”的相关,将相关值“P(Δt)”向CPU 25输出。
CPU 25是根据存储在存储部27中的系统程序等各种程序对基带处理电路部20的各部分进行统一控制的处理器。对于各捕获对象卫星,CPU25(相当于捕获部)进行如下的处理:对于各捕获对象卫星,根据从相关处理部217输出的相关值“P(Δt)”,检测码相位。用检测出的码相位计算出该捕获对象卫星与便携式电话机1之间的伪距,利用算出的伪距进行位置算出计算,从而算出便携式电话机1的位置。
为了向复制码信号发生部213输出用于指示捕获对象卫星的CA码(捕获对象卫星的PRN号)的CA码指示信号,而使复制码信号发生部213生成捕获对象卫星的复制码信号“CAR(t)”。向相关处理部217输出用于指示差动复制码信号“wR(t)”的相位“Δt”的相位指示信号,在相关处理部217中,在使差动复制码信号“wR(t)”的相位“Δt”发生变化的同时,执行相关处理。
存储部27由ROM(Read Only Memory,只读存储器)、闪存ROM和RAM(Random Access Memory,随机存取存储器)等存储装置构成,CPU 25存储用于控制基带处理电路部20的系统程序、用于实现位置计算功能的各种程序、数据等。形成暂时存储由CPU 25执行的系统程序、各种处理程序、各种处理的处理中数据、处理结果等的工作区。
回到图9的说明。主CPU 30是根据存储部80存储的系统程序等各种程序,统一控制便携式电话机1的各部分的处理器。主CPU 30进行把从基带处理电路部20输入的位置信息显示在显示部50上的处理、或进行利用该位置信息的各种应用处理。
操作部40是由如触摸面板和按钮开关构成的输入装置,其把被按下的键或按钮的信号向主CPU 30输出。通过该操作部40的操作,进行通话请求、电子邮件收发请求以及位置算出请求等的指令输入。
显示部50由LCD(Liquid Crystal Display,液晶显示器)等构成,显示部50是根据从主CPU30输入的显示信号进行显示的显示装置。在显示部50中显示有位置显示画面和时间信息等。
便携式电话用天线60是在与便携式电话机1的通信服务企业设置的无线基站之间发送和接收便携式电话用无线信号的天线。
便携式电话用无线通信电路部70是由RF转换电路、基带处理电路等构成的便携式电话的通信电路部,通过调制/解调便携式电话用无线信号,从而实现通话和电子邮件的收发等。
存储部80是存储主CPU 30用于控制便携式电话机1的系统程序、用于实现位置计算功能的程序及数据等的存储装置。
(2)处理流程
图13是表示基带处理电路部20的CPU 25执行的基带处理流程的流程图。虽然没有明确提出,但在以下基带处理的执行过程中,进行基于GPS天线9的RF信号的接收、基于RF接收电路部11的RF信号向IF信号的降频变频,设定为使IF信号处于随时向基带处理电路部20输出的状态。
首先,CPU 25进行捕获对象卫星的判断处理(步骤A1)。具体地说,在由图中未显示的钟表部所计时的当前时刻,用历书、星历表等卫星轨道数据判断位于所给基准位置的天空中的GPS卫星,并将其作为捕获对象卫星。基准位置是用如下方法来确定的:例如,在对接通电源后的初始位置进行计算的情况下,把基准位置设定为通过辅助服务器从便携式电话机1的基站取得的位置,在计算第二次以后的位置的情况下,把基准位置设定为根据上一次位置计算求得的最新GPS计算位置。
接着,CPU 25对在步骤A1中判断出的各捕获对象卫星执行循环A的处理(步骤A3~A11)。在循环A的处理中,CPU 25根据该卫星轨道数据,算出该捕获对象卫星的卫星位置、卫星移动速度及卫星移动方向等卫星信息(步骤A5)。CPU 25进行码相位检测处理(步骤A7)。
图14是表示码相位检测处理的流程的流程图。
首先,CPU 25把该捕获对象卫星的CA码指示信号向复制码信号发生部213输出(步骤B1)。CPU 25设定相位的搜索范围及相位的搜索幅宽(步骤B3,B5)。CPU 25根据设定的相位的搜索范围及搜索幅宽,确定作为搜索对象的相位(搜索相位)(步骤B7)。
接着,CPU 25对在步骤B7设定的各搜索相位执行循环B的处理(步骤B9~B15)。在循环B的处理中,CPU 25将该搜索相位“Δt”的指示信号向相关处理部217输出(步骤B11)。若执行步骤B11,则相关处理部217根据在图1~图3中说明的原理进行相关处理,并把相关值“P(Δt)”向CPU 25输出。CPU 25当从相关处理部217输入相关值“P(Δt)”时,把该相关值“P(Δt)”存储到存储部27中(步骤B13)。CPU 25将处理移至下一个搜索相位。
在对所有的搜索相位进行步骤B11及步骤B13的处理之后,CPU 25结束循环B的处理(步骤B15)。CPU 25把存储在存储部27中的相关值“P(Δt)”变为最大的搜索相位“Δt”确定为码相位(步骤B17)。CPU 25结束码相位的检测处理。
回到图13的基带处理,在完成了码相位检测处理之后,CPU 25用在步骤A5中算出的卫星信息和在步骤A7中检测出的码相位,算出该捕获对象卫星与便携式电话机1之间的伪距(步骤A9)。可用例如最新的GPS算出位置和卫星位置来算出伪距的整数部分,可用码相位来算出伪距的小数部分(fractional part)。在算出伪距后,CPU 25进入对下一个捕获对象卫星的处理。
在对所有的捕获对象卫星进行步骤A5~A9的处理之后,CPU 25结束循环A的处理(步骤A11)。然后,CPU 25利用在步骤A9中对各捕获对象卫星算出伪距进行GPS位置计算处理,并算出便携式电话机1的位置(步骤A13)。由于关于利用伪距的位置算出计算是已知的,所以省略对其的详细说明。
接着,CPU 25把通过GPS位置计算处理算出的位置向主CPU30输出(步骤A15)。CPU 25判断位置计算是否结束(步骤A17),在判断出还未结束的情况下(步骤A17中为否),回到步骤A1。在判断出结束了位置计算的情况下(步骤A17中为是),则结束基带处理。
(3)作用效果
在基带处理电路部20的卫星信号捕获部21中,在接收信号差动电路部211中使接收了从GPS卫星发送的GPS卫星信号而获得的接收信号r(t)差动,生成CA码的差动码信号w(t-T)。由复制码信号发生部213生成的复制码信号CAR(t)在复制码信号差动电路部215中差动,生成差动复制码信号wR(t-T)。在相关处理部217中进行差动码信号w(t-T)和差动复制码信号wR(t-T)之间的相关处理,相关值P(Δt)被输出到CPU 25。由CPU 25进行基于相关值P(Δt)的GPS卫星信号的捕获。
把GPS卫星信号的扩展码、即CA码差动后的差动码信号是不依赖于接收信号的频率,而用接收信号和使该接收信号延迟规定的时间的延迟信号就可以生成的信号。因此,通过进行使用差动码信号和差动复制码信号的相关处理,并依照该相关值捕获GPS卫星信号,从而无需进行频率方向的相关运算,并可在短时间内捕获GPS卫星信号。
也就是说,进行使用差动码信号和差动复制码信号的相关处理而得到的相关值是相当于对接收信号进行频率方向的搜索得到的相关值的相关值。因此,在捕获卫星信号时,不用频率方向的搜索,就可削减现有技术中在频率方向搜索上花费的时间,所以只需在相位方向的搜索所需的时间即可。
本发明人在通过现有的方法对20个频率进行频率搜索(频率方向的相关运算)来捕获GPS卫星信号的情况下、和在通过本发明的方法来捕获GPS卫星信号的情况下,进行分别调查捕获GPS卫星信号所需的时间的实验。其结果是,相对于通过现有方法需要“2.8秒”,用本发明的方法“0.12秒”即可。该时间是现有方法时间的约“1/20”,所以无需对20个频率进行频率搜索就能捕获GPS卫星信号。也对通过本发明方法检测出的码相位的正确性进行了调查,确认了可以比现有方法高的正确性检测码相位。
2-2.第二实施例
第二实施例是适用了上述原理“1-2.用差动码预测CA码的预测方法”的实施例。第一实施例的便携式电话机1以及卫星信号捕获部之外的功能模块及处理流程是共通的,因此,对不同部分即卫星信号捕获部进行说明。
(1)结构
图15示出了第二实施例中的卫星信号捕获部22的电路结构的一例。对与图10的卫星信号捕获部21相同的结构元件标注了相同的标号并省略对其的说明。
卫星信号捕获部22构成为包括:接收信号差动电路部211、CA码预测处理部221、复制码信号发生部213、相关处理部217和延迟电路部DC。
CA码预测处理部221用从接收信号差动电路部211输出的差动码信号“w(t-T)”,根据图4~图8说明的原理,进行预测CA码的处理。生成预测到的CA码信号、即预测CA码信号(相当于预测码),向相关处理部217输出。
由复制码信号发生部213生成的复制码信号“CAR(t)”被延迟电路部DC延迟一时钟后,向相关处理部217输出。相关处理部217进行从CA码预测处理部221输入的预测CA码信号“CAP(t-T)”与从复制码信号发生部213输入的复制码信号“CAR(t-T)”的相关处理。这时,根据CPU 25的相位指示信号,在使复制码信号“CAR(t-T)”的相位“Δt”发生变化的同时,进行“CAP(t-T)”与“CAR(t-T+Δt)”之间的相关运算。把得到的相关值“P(Δt)”向CPU 25输出。
(2)作用效果
在第二实施例中,在CA码预测处理部221中进行预测接收信号的CA码“CA(t)”的处理,在相关处理部217中,进行预测CA码信号“CAP(t-T)”与复制码信号“CAR(t-T)”之间的相关处理。通过根据上述原理的预测处理,可正确预测CA码“CA(t)”。因此,通过相关处理得到的相关值可信度变高,可以提高GPS卫星信号捕获的可靠性。
3.变形例
3-1.适用系统
在上述实施方式中,举例说明了GPS卫星信号的捕获,然而本发明可同样适用于接收GPS卫星信号之外的信号的接收装置。即,只要是接收通过规定的扩展码被扩展调制的信号的接收装置,则本发明都能适用。只要是在信号捕获上需要频率搜索的接收装置,则应用本发明的效果更好。
3-2.电子设备
在上述实施方式中,以在作为电子设备的一种的便携式电话机中适用本发明的情况为例进行了说明,然而,可适用本发明的电子设备并不仅限于便携式电话机。例如,对于汽车导航装置或便携式导航设备、个人计算机、PDA(Personal Digital AssistanT,个人数字助理)、手表等其他的电子设备也同样能适用。
3-3.卫星定位系统
在上述实施方式中,虽然已经将GPS作为卫星定位系统的情况为例进行了说明,但是也可以是WAAS(Wide Area Augmentation System,广域增强系统)、QZSS(Quasi Zenith Satellite System,准天顶卫星系统)、GLONASS(GLObal NAvigation Satellite System,全球导航卫星系统)、GALILEO(伽俐略系统)等其他卫星定位系统。
3-4.相关部的并行
在第一实施例中,已经对如下的情况进行了说明:相关处理部217根据来自CPU 25的相位指示信号,在使差动复制码信号“wR(t)”的相位“Δt”发生变化的同时,串行执行相关处理。但是,也可以是并列进行差动码信号“w(t)”和差动复制码信号“wR(t)”之间的相关运算的电路结构。
图16示出了该变形例中的卫星信号捕获部23的电路构成的一例。对于与图10的卫星信号捕获部21相同的构成元件标注了相同的符号,并省略对其的说明,将以与卫星信号捕获部21不同的部分为中心进行说明。
在卫星信号捕获部23中构成为:通过延迟电路部DC将由复制码信号发生部213生成的复制码信号“CAR(t)”逐个相位“Δt”地多级延迟。构成第一卫星信号处理电路部231-1~第N卫星信号处理电路部231-N这N个卫星信号处理电路部231,以便将通过延迟电路部DC延迟为“n级(n=1,2,…N)”的复制码信号“CAR(t-N··Δt)”输入给第n卫星信号处理电路部231-n。
各卫星信号处理电路部231构成为包括:复制码信号差动电路部215和相关部233。向第n卫星信号处理电路部231-n输入的复制码信号“CAR(t-N··Δt)”通过复制码信号差动电路部215转换为差动复制码信号“wR(t-T-N··Δt)”。利用相关部233,进行差动码信号“w(t-T)”与差动复制码信号“wR(t-T-N··Δt)”的相关运算,相关值“P(N··Δt)”被输出到CPU 25。
3-5.差动码信号的生成
在上述实施例中,虽然已经对以时间“t”为基准,用接收信号“r(t)”、使接收信号延迟一时钟得到的第一延迟信号“r(t-T)”和使接收信号延迟两时钟得到的第二延迟信号“r(t-2T)”生成的差动码信号“w(t)”的情况为例进行了说明,但是延迟的时钟“T”是可变的,也能生成多个差动码信号“w(t)”。
原理为,在“m”是一个整数的情况下,用接收信号“r(t)”、第一信号“r(t-m··T)”和第二信号“r(t+m··T)”,根据下式(10)和(11),计算出差动码“w(t,m)”。
【式10】
r(t-m·T)·r(t+m·T)
=CA(t-m·T)·eiω(t-m·T)·CA(t+m·T)·eiω(t+m·T)
=CA(t-m·T)·CA(t+m·T)·e2iωt…(10)
【式11】
w(t,m)=CA(t-m·T)·CA(t+m·T)
=r(t-m·T)·r(t+m·T)·e-2iωt
=r(t-m·T)·r(t+m·T)·{r(t)2}*…(11)
对于CA码复制码,也同样根据下式(12)计算CA码复制码的差动码“wR(t,m)”。
【式12】
wR(t,m)=CAR(t-m·T)·CAR(t+m·T)…(12)
这时,用CA码的差动码“w(t,m)”和CA码复制码的差动码“wR(t+Δt,m)”,根据下式(13),算出相关值“P(Δt,m)”。
【式13】
P ( Δt , m ) = ( Σ i g ( t , m ) ) 2 = ( Σ i w ( t , m ) · w R ( t + Δt , m ) ) 2
= ( Σ i CA ( t - m · T ) · CA ( t + m · T ) · CA R ( t - m · T + Δt ) · CA R ( t + m · T + Δt ) ) 2 . . . ( 13 )
如果对于多个“m”而计算相关值“P(Δt,m)”,则根据下式(14),算出合计相关值“PA(Δt)”。
【式14】
P A ( Δt ) = ( Σ m P ( Δt , m ) ) 2 . . . ( 14 )
在使CA码复制码的差动码“wR(t,m)”的相位“Δt”发生变化的同时,根据式(14)算出合计相关值“PA(Δt)”,通过检测出合计相关值“PA(Δt)”变为最大的相位“Δt”,从而检测出码相位。
图17示出了该变形例中的卫星信号捕获部24的电路结构的一例。对于与图10的卫星信号捕获部21相同的构成元件标注了相同的标号,并省略对其的说明。为了方便观看附图,省略了从CPU 25向复制码信号发生部213发出的CA码指示信号及向相关处理部217发出的相位指示信号的图示。
在卫星信号捕获部24中,设置了N个接收信号差动电路部211-1~211-N。接收信号差动电路部211-1~211-N是分别与上述整数“m=1~N”对应的差动电路部,延迟电路部DC的延迟时间分别是“T~NT”。对于复制码信号差动电路部,也同样地设置有N个复制码信号差动电路部215-1~215-N。设置有与各个差动电路部相对应的N个相关处理部217-1~217-N。
通过将接收信号“r(t)”向接收信号差动电路部211-1~211-N的输入,从而可以生成N个差动码信号“w(t-T,1)”~“w(t-T,N)”,并输出到对应的相关处理部217-1~217-N。通过将复制码信号“CAR(t)”输入到复制码信号差动电路部215-1~215-N,从而可以生成N个差动复制码信号“wR(t-T,1)”~“wR(t-T,N)”,并被输出到对应的相关处理部217-1~217-N。
在各相关处理部217-1~217-N中,算出相关值“P(Δt,1)”~“P(Δt,N)”,并输出给合计部241。在合计部241中,计算对相关值“P(Δt,1)”~“P(Δt,N)”进行合计后得到的合计相关值“PA(Δt)”,输出到CPU 25。
如上所述,利用以不同的时间延迟的多种延迟信号生成多个差动码信号“w(t-T,m)”,再分别用这些多个的差动码信号进行相关处理,就能得到互相独立的多个相关值“P(Δt,m)”。把上述相关值合计,得到合计相关值“PA(Δt)”,根据合计相关值“PA(Δt)”检测出码相位,这样可提高GPS卫星信号捕获的准确性。
虽然已经对使接收信号逐个延迟把时钟“T”扩大“m=1~N”倍后而得到时间“m··T”,生成多个差动码信号“w(t,m)”的情况进行了说明,但也可以任意地设定各差动电路部的延迟时间。即,也可以在接收信号差动电路部211-1~211-N中,对各自的延迟电路DC设定任意的延迟时间“T1”~“TN”,生成差动码信号“w(t,Tm)”。对于复制码信号差动电路部215-1~215-N也一样。
3-6.信号的差动(1)
虽然在上述实施例中已经对通过把接收信号输入接收信号差动电路部211,从而用电路进行接收信号的差动的情况进行了说明,但也可以使用软件进行数字信号处理。对于复制码信号的差动也一样地可以用软件加以进行。
3-7.信号的差动(2)
虽然在上述实施例中已经对通过用接收信号差动电路部211使接收信号延迟而发生差动的情况进行了说明。但是,也可以把从RF接收电路部11输入的接收信号暂时保存到采样存储器中,并从采样存储器中读出“r(t)”(第一接收信号)、与第一接收信号相差时间T的“r(t-T)”(第二接收信号)、与第一接收信号相差时间2T的“r(t-2T)”(第三接收信号),利用该读出的各接收信号来产生差动。除此之外,可能适用可以获得“r(t)”、与该接收信号相差时间T的“r(t-T)”及与该接收信号相差时间2T的“r(t-2T)”的各种方法。
符号的说明
1   便携式电话机                10  GPS接收部
11  RF接收电路部                20  基带处理电路部
21  卫星信号捕获部              25  CPU
27  存储部                        30  主CPU
40  操作部                        50  显示部
60  便携式电话用天线
70  便携式电话用无线通信电路部
80  存储部                        211 接收信号差动电路部
213 复制码信号发生部              215 复制码信号差动电路部
217 相关处理部                    DC  延迟电路部
PS  复共轭运算部                  M   乘算部

Claims (10)

1.一种卫星信号捕获方法,其特征在于,包括:
使用接收来自定位用卫星的卫星信号所得的第一接收信号至第三接收信号,生成对所述卫星信号的扩展码实施差动后的差动码;
进行使用所述差动码和作为所述扩展码的复制码的扩展码复制码的规定的相关处理,取得相关值;以及
根据所述相关值捕获所述卫星信号,
其中,取得所述相关值包括:
根据所述差动码预测所述扩展码;以及
把作为预测到的所述扩展码的预测码与所述扩展码复制码进行相关运算。
2.根据权利要求1所述的卫星信号捕获方法,其特征在于,
所述第二接收信号是与所述第一接收信号相差时间T的信号,
所述第三接收信号是与所述第一接收信号相差时间2T的信号。
3.根据权利要求2所述的卫星信号捕获方法,其特征在于,
生成所述差动码包括:
通过自乘所述第二接收信号来计算复共轭;以及
把所述第一接收信号、所述复共轭的信号、所述第三接收信号相乘来计算所述差动码。
4.根据权利要求2或3所述的卫星信号捕获方法,其特征在于,
所述卫星信号捕获方法还包括:
使用所述第一接收信号、与所述第一接收信号相差时间T'的第四接收信号以及与所述第一接收信号相差时间2T'的第五接收信号,生成第二差动码,其中,所述时间T'与所述时间T不同;
进行使用所述第二差动码和所述扩展码复制码的规定的相关处理,取得第二相关值;以及
根据所述相关值及所述第二相关值捕获所述卫星信号。
5.根据权利要求1至3中任一项所述的卫星信号捕获方法,其特征在于,
取得所述相关值包括:
对所述扩展码复制码实施差动,生成差动复制码;以及
把所述差动码和所述差动复制码进行相关运算。
6.根据权利要求4所述的卫星信号捕获方法,其特征在于,
取得所述相关值包括:
对所述扩展码复制码实施差动,生成差动复制码;以及
把所述差动码和所述差动复制码进行相关运算。
7.根据权利要求2或3所述的卫星信号捕获方法,其特征在于,
所述扩展码复制码包括第一扩展码复制码、与所述第一扩展码复制码相差所述时间2T的第二扩展码复制码,
取得所述相关值包括:
通过把所述第一扩展码复制码和所述第二扩展码复制码相乘,从而生成对所述扩展码复制码实施差动后的差动复制码;以及
把所述差动码和所述差动复制码进行相关运算。
8.根据权利要求4中所述的卫星信号捕获方法,其特征在于,
所述扩展码复制码包括第一扩展码复制码、与所述第一扩展码复制码相差所述时间2T的第二扩展码复制码,
取得所述相关值包括:
通过把所述第一扩展码复制码和所述第二扩展码复制码相乘,从而生成对所述扩展码复制码实施差动后的差动复制码;以及
把所述差动码和所述差动复制码进行相关运算。
9.根据权利要求1至3中任一项所述的卫星信号捕获方法,其特征在于,
所述卫星信号捕获方法还包括:通过错开所述接收信号的相位并进行使用所述差动码和所述扩展码复制码的所述规定的相关处理来搜索相关峰值,
捕获所述卫星信号是根据所述相关峰值来捕获所述卫星信号,而不需要错开频率地搜索相关峰值来捕获卫星信号。
10.根据权利要求4所述的卫星信号捕获方法,其特征在于,
所述卫星信号捕获方法还包括:通过错开所述接收信号的相位并进行使用所述差动码和所述扩展码复制码的所述规定的相关处理来搜索相关峰值,
捕获所述卫星信号是根据所述相关峰值来捕获所述卫星信号,而不需要错开频率地搜索相关峰值来捕获卫星信号。
CN201010568318.1A 2009-11-30 2010-11-30 卫星信号捕获方法及卫星信号接收装置 Expired - Fee Related CN102081161B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2009-271897 2009-11-30
JP2009271897A JP4911220B2 (ja) 2009-11-30 2009-11-30 衛星信号捕捉方法及び衛星信号受信装置

Publications (2)

Publication Number Publication Date
CN102081161A CN102081161A (zh) 2011-06-01
CN102081161B true CN102081161B (zh) 2015-01-14

Family

ID=44068893

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010568318.1A Expired - Fee Related CN102081161B (zh) 2009-11-30 2010-11-30 卫星信号捕获方法及卫星信号接收装置

Country Status (3)

Country Link
US (2) US8462833B2 (zh)
JP (1) JP4911220B2 (zh)
CN (1) CN102081161B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102419444B (zh) * 2011-09-06 2013-09-04 厦门大学 基于统计的捕获判别方法
FR3003357B1 (fr) * 2013-03-15 2015-03-20 Thales Sa Procede et dispositif de compression d'un signal de radio-navigation large bande, procede et dispositif associes de calcul de la fonction de correlation du code d'etalement dudit signal compresse
GB201401237D0 (en) * 2014-01-24 2014-03-12 Qinetiq Ltd Improvements to satellite positioning methods
FR3060245B1 (fr) * 2016-12-12 2019-05-17 Uwinloc Procede et dispositif pour la detection d’une impulsion d’un signal

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1303483A (zh) * 1998-04-14 2001-07-11 施耐普特拉克股份有限公司 快速捕获的高灵敏度gps接收机
CN1615596A (zh) * 2002-08-30 2005-05-11 英特尔公司 用于接收差分超宽带信号的方法和装置
CN101042430A (zh) * 2006-03-23 2007-09-26 精工爱普生株式会社 定位装置、定位控制方法、定位控制程序及存储介质
CN101088025A (zh) * 2004-12-21 2007-12-12 诺基亚公司 代码调制信号的获取
CN101356449A (zh) * 2006-05-26 2009-01-28 全球定位有限公司 为接收自多种卫星系统的卫星信号执行信号相关的方法和装置
CN101568851A (zh) * 2006-12-27 2009-10-28 英特尔公司 用于捕获gps信号的方法及具有采样时间误差和频率偏移补偿的gps接收机

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2721474B2 (ja) * 1992-06-29 1998-03-04 三菱電機株式会社 スペクトル拡散通信用受信装置
JP3298597B2 (ja) * 1993-08-03 2002-07-02 日本無線株式会社 Gps受信装置
JP2005011977A (ja) * 2003-06-18 2005-01-13 Ebara Corp 基板研磨装置および基板研磨方法
EP1545019A1 (en) * 2003-12-19 2005-06-22 Telefonaktiebolaget LM Ericsson (publ) GPS receiver using differential correlation
JP4595855B2 (ja) 2006-03-23 2010-12-08 セイコーエプソン株式会社 測位装置、測位装置の制御方法、測位装置の制御プログラム、測位装置の制御プログラムを記録したコンピュータ読み取り可能な記録媒体
US7570206B2 (en) * 2006-03-23 2009-08-04 Seiko Epson Corporation Positioning device, positioning control method, positioning control program, and computer-readable recording medium having positioning control program recorded thereon
CA2669065A1 (en) * 2006-11-07 2008-05-15 Meadwestvaco Packaging Systems, Llc Integrated secondary and tertiary packaging machine
JP4359638B2 (ja) * 2007-08-24 2009-11-04 Okiセミコンダクタ株式会社 相関演算器及び相関演算装置
JP5526492B2 (ja) * 2008-04-22 2014-06-18 セイコーエプソン株式会社 擬似距離算出方法、測位方法、プログラム及び測位装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1303483A (zh) * 1998-04-14 2001-07-11 施耐普特拉克股份有限公司 快速捕获的高灵敏度gps接收机
CN1615596A (zh) * 2002-08-30 2005-05-11 英特尔公司 用于接收差分超宽带信号的方法和装置
CN101088025A (zh) * 2004-12-21 2007-12-12 诺基亚公司 代码调制信号的获取
CN101042430A (zh) * 2006-03-23 2007-09-26 精工爱普生株式会社 定位装置、定位控制方法、定位控制程序及存储介质
CN101356449A (zh) * 2006-05-26 2009-01-28 全球定位有限公司 为接收自多种卫星系统的卫星信号执行信号相关的方法和装置
CN101568851A (zh) * 2006-12-27 2009-10-28 英特尔公司 用于捕获gps信号的方法及具有采样时间误差和频率偏移补偿的gps接收机

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JP特开2007-256111A 2007.10.04 *
高动态单频干扰环境下的FFT直扩捕获方案;张一等;《哈尔滨工业大学学报》;20070930;第39卷(第9期);1431-1435 *

Also Published As

Publication number Publication date
US20110129003A1 (en) 2011-06-02
US8462833B2 (en) 2013-06-11
US20130251005A1 (en) 2013-09-26
US8750348B2 (en) 2014-06-10
CN102081161A (zh) 2011-06-01
JP2011112616A (ja) 2011-06-09
JP4911220B2 (ja) 2012-04-04

Similar Documents

Publication Publication Date Title
US8634504B2 (en) Correlation calculating method, satellite signal capturing method, and correlation calculating circuit
US8855173B2 (en) Method of capturing satellite signal and device for capturing satellite signal
US8362953B2 (en) Sequential chip correlation array
CN102081161B (zh) 卫星信号捕获方法及卫星信号接收装置
US8842718B2 (en) Method of capturing satellite signal and device for capturing satellite signal
CN101329396A (zh) 相干累积增进方法及电路、定位方法及电路、电子设备
US8174440B2 (en) Positioning method, program thereof, and positioning device
CN101487890B (zh) 定位方法、程序以及定位装置
JP2008281552A (ja) 初回測位出力位置演算決定方法、プログラム、記憶媒体、測位装置及び電子機器
US20110235687A1 (en) Signal acquisition method, signal acquisition apparatus and electronic device
US7630449B2 (en) Method and receiver capable of increasing acquisition efficiency in a code division multiple access communication system
US8767891B2 (en) Reception signal integrating method and receiver
JP2011203234A (ja) 信号捕捉方法、信号捕捉装置及び電子機器
US8761232B2 (en) Bit change determination method and bit change determination device
JP5488380B2 (ja) 衛星信号捕捉方法及び衛星信号受信装置
JP2011137802A (ja) 受信信号積算方法及び受信装置
JP2008281551A (ja) 現在位置測位方法、プログラム、記憶媒体、測位装置及び電子機器
JP6011042B2 (ja) 受信部駆動制御方法及び受信装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150114

Termination date: 20161130

CF01 Termination of patent right due to non-payment of annual fee