CN104204850A - 信号搜索方法、信号搜索程序、信号搜索装置、gnss信号接收装置及信息设备终端 - Google Patents

信号搜索方法、信号搜索程序、信号搜索装置、gnss信号接收装置及信息设备终端 Download PDF

Info

Publication number
CN104204850A
CN104204850A CN201380015191.2A CN201380015191A CN104204850A CN 104204850 A CN104204850 A CN 104204850A CN 201380015191 A CN201380015191 A CN 201380015191A CN 104204850 A CN104204850 A CN 104204850A
Authority
CN
China
Prior art keywords
frequency
signal
mentioned
correlation
peak
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.)
Granted
Application number
CN201380015191.2A
Other languages
English (en)
Other versions
CN104204850B (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.)
Furuno Electric Co Ltd
Original Assignee
Furuno Electric Co Ltd
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 Furuno Electric Co Ltd filed Critical Furuno Electric Co Ltd
Publication of CN104204850A publication Critical patent/CN104204850A/zh
Application granted granted Critical
Publication of CN104204850B publication Critical patent/CN104204850B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details 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/69—Spread spectrum techniques
    • H04B1/713—Spread spectrum techniques using frequency hopping
    • H04B1/715—Interference-related aspects
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/30—Acquisition or tracking or demodulation of signals transmitted by the system code related
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details 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/69—Spread spectrum techniques
    • H04B1/707—Spread spectrum techniques using direct sequence modulation
    • H04B1/7073—Synchronisation aspects
    • H04B1/7075—Synchronisation aspects with code phase acquisition
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details 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/69—Spread spectrum techniques
    • H04B1/713—Spread spectrum techniques using frequency hopping
    • H04B1/7136—Arrangements for generation of hop frequencies, e.g. using a bank of frequency sources, using continuous tuning or using a transform
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details 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/69—Spread spectrum techniques
    • H04B1/713—Spread spectrum techniques using frequency hopping
    • H04B1/7143—Arrangements for generation of hop patterns

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)
  • Monitoring And Testing Of Transmission In General (AREA)

Abstract

判别是否能够比以往方法更正确地捕捉目的的GPS信号。将由GPS信号的码周期决定的1000Hz间隔的多个频率作为1个组，以全部的码相位取得累积相关值(S101)。取得构成组的各频率i下的最大的累积相关值即峰相关值Peak(i)(S102)。将各频率i的峰相关值Peak(i)比较，将它们的最大值设定为最大峰值Tpeak，将第2高的值的峰相关值Peak(i)设定为最大噪声值Np。如果最大峰值Tpeak与最大噪声值Np的比Ratio比阈值Th低，则判定为交叉相关(S104：是)，如果较高，则判定最大峰值Tpeak是由目的的GPS信号带来的(S104：否)，捕捉GPS信号(S105)。

Description

信号搜索方法、信号搜索程序、信号搜索装置、GNSS信号接收装置及信息设备终端
技术领域
本发明涉及从接收信号中搜索希望的信号的信号搜索方法，特别涉及搜索GNSS中的GPS信号的方法。
背景技术
目前，作为GNSS(Global Navigation Satellite System：全球导航卫星系统)之一，有GPS(Global Positioning System)。
在GPS中，接收从多个GPS卫星发送的GPS信号，使用接收到的GPS信号的码相位及载波相位进行接收机的测位。在GPS中，按照GPS卫星设定不同的扩散码，将各GPS信号用不同的扩散码进行码调制。
在这样的GPS中，有将与作为目的的GPS卫星不同的GPS卫星的GPS信号误捕捉为作为目的的GPS信号而进行跟踪处理的情况。将这样的现象称作交叉相关。
作为防止交叉相关的方法，在专利文献1中，通过比较单一的搜索频率的码相位上的各码相位点的累积相关值比较，判定所取得的累积相关值是由作为目的的GPS信号带来的、还是由交叉相关带来的。图1是表示专利文献1所记载的信号搜索方法的概念的图。
在专利文献1的方法中，取得单一的搜索频率的各码相位下的累积相关值，将这些值比较。在专利文献1的方法中，检测全部的码相位下的累积相关值中的作为最大值的累积相关值Vpk1。接着，在专利文献1的方法中，检测全部的码相位下的累积相关值中的值第2高的累积相关值Vpk2。在专利文献1的方法中，计算作为最大值的累积相关值Vpk1与值第2高的累积相关值Vpk2的比。在专利文献1的方法中，如果该比是规定阈值以上，则认为与作为最大值的累积相关值Vpk1对应的码相位Cpk1是由目的信号带来的，捕捉该码相位Cpk1的信号。另一方面，如果比是阈值以下，则继续信号搜索。
此外，在专利文献2中，使用以下的方法。首先，根据以规定频率间隔且一定的累积时间计算累积相关值的结果检测强信号。利用检测出的强信号的频率与捕捉对象的频率的差与信号水平有规定的关系，进行交叉相关的检测。并且，如果将检测出的强信号判定为交叉相关，则继续信号搜索。
如上述那样，在以往的信号搜索方法中，在能够将搜索范围即执行信号搜索的码相位范围及频率范围的整体搜索之前，各码相位、频率下的累积时间是一定的。并且，在使这样的整体搜索至少结束1次(根据情况是多次)后，如果能够检测出作为目的的信号，则使累积时间变长，执行信号搜索。
此外，在检测弱信号的情况下，与强信号检测用的累积时间相比，将累积时间设定得较长。
现有技术文献
专利技术
专利文献1：美国专利7161977号说明书
专利文献2：美国专利7623070号说明书
发明内容
发明概要
发明要解决的课题
但是，在专利文献1的方法中，必须在检测作为最大的累积相关值Vpk1后，再次在码相位轴上扫描，检测值第2高的累积相关值Vpk2。因而，与单单检测作为最大的累积相关值Vpk1的码相位Cpk1的情况相比，需要2倍的搜索时间。
此外，在专利文献1的方法中，有第2高的累积相关值为朝向包括作为最大的累积相关值Vpk1的极大值的中途的码相位(如果是图1的例子，则为码相位Cpk2w)的情况。因而，必须以不包含包括累积相关值Vpk1的码相位Cpk1的规定码相位范围的方式进行第2高的累积相关值的检测，处理复杂化。
进而，如果将码相位轴上的第2高的累积相关值的检测范围不是设为码相位全域，而是如图1所示那样限制为部分性的码相位范围Cov，则有检测出值比真的第2高的累积相关值Vpk2低的累积相关值Vpk2'的情况。在使用该累积相关值Vpk2'的情况下，与作为最大值的累积相关值Vpk1的比变得比真的比大。因此，即便是交叉相关，也容易发生判定为目的的GPS信号的误判定。
此外，在专利文献2的信号搜索方法中，如上述那样，遍及搜索范围全域以相同的累积时间计算累积相关值，所以即使例如将累积时间设定得较短，也不能检测弱信号等，所以必须将搜索范围全域反复搜索，结果累积时间变长。另一方面，如果将累积时间设定得较长，则码相位、频率的二维范围中的各相关累积值计算点处的累积时间变长，对于搜索范围全域的搜索时间变长。
此外，在遍及搜索范围全域以相同的累积时间计算累积相关值的情况下，容易受到交叉相关的影响，不能使发生上述误检测的概率降低。
这样，在以往的专利文献1、专利文献2所示的方法中，有关交叉相关的判定及真的信号的捕捉的时间变长，并且有可能误判定交叉相关。
因而，本发明的目的是提供一种能够通过比以往方法简单的处理、比以往更正确地捕捉目的信号的信号搜索方法。
解决课题的手段
本发明是用来捕捉目的信号的信号搜索方法，具有相关值取得工序和判定工序。在相关值取得工序中，将以根据将目的信号调制的扩散码的码周期决定的频率间隔设定的多个频率作为1个组，进行接收信号与以多个频率生成的目的信号的副本(replica复制)信号的相关处理而取得相关值。在判定工序中，根据该组内的多个相关值进行是否捕捉到目的信号的判定。
在该方法中，基于以下的两个相关值的频率特性，判定是否捕捉到目的信号。
(i)交叉相关以与扩散码的码周期(频率)对应的间隔发生相关值的峰。即，以特定的频率间隔检测峰相关值。
(ii)在真的信号(目的信号)中，利用仅峰相关值突出而较高的特性。即，即使是以上述特定的频率间隔设定的多个频率，也只有这些多个频率中的特定的1个频率的相关值比其他频率的相关值大幅地变高。
如果利用这样的相关值的频率特性，则仅通过将由同时取得的多个频率下的最大相关值构成的峰相关值比较，就能够进行是否捕捉到目的信号的判定。由此，能够以比较简单的处理且正确地进行是否捕捉到目的信号的判定。
此外，在该发明的信号搜索方法中，目的信号是GPS信号的L1波；频率间隔是1000Hz间隔。这表示作为目的信号而具体地捕捉GPS信号的L1波的情况。作为将GPS信号的L1波调制的扩散码的C/A码由1023码片(chip)构成，比特率是1.023Mbps，码周期是1msec.。通过使用这样的C/A码的L1波发生的交叉相关以约1000Hz间隔具有相关值的极大值。因而，如果以1000Hz间隔将相关值比较，则如上述那样，交叉相关与目的信号的相关特性的差变得明确。
此外，在该发明的信号搜索方法中，在判定工序中，检测组内的各频率的最大相关值，将该最大值设定为该频率的峰相关值。判定工序基于各频率的峰相关值的大小关系，判定是否捕捉到目的信号。
在该方法中，表示是否捕捉到目的信号的判定的具体的方法例。如上述那样，交叉相关中，基于码周期的频率间隔下的相关值大致相同。目的信号仅在特定的一频率下相关值变高。因而，如果将各频率的峰相关值的大小关系比较，则能够正确地进行是否捕捉到目的信号的判定。
此外，在该发明的信号搜索方法中，判定工序具有：检测作为各频率的峰相关值的最大值的最大峰值的工序；根据除了最大峰值以外的各频率的峰相关值决定最大噪声值的工序；计算最大峰值与最大噪声值的比的工序。并且，在判定工序中，基于这样得到的比，判定是否捕捉到上述目的信号。
在该方法中，表示更具体的判定方法的形态。
此外，在该发明的信号搜索方法中，当通过判定工序判定为没有捕捉到目的信号时，使多个频率以比频率间隔短的一定的频率宽度进行频率变动。
如该方法所示，通过一边使频率以一定的频率变动，一边重复是否捕捉到目的信号的判定，能够执行是否捕捉到对于要进行捕捉用的扫描的全部频率带域的目的信号的判定，能够正确且可靠地捕捉目的信号。
此外，在该发明的信号搜索方法中，相关值取得工序具有搜索基准频率设定工序、累积时间决定工序和相关值计算工序。在搜索基准频率设定工序中，根据跟踪中的信号的频率和频率间隔设定搜索基准的频率。在累积时间决定工序中，根据搜索基准频率与搜索对象的频率的频率差，设定该搜索对象的频率下的累积时间。在相关值计算工序中，通过以设定的累积时间将相关值累积，计算累积相关值。
在该方法中，利用跟踪中的信号为作为目前搜索对象而捕捉目的信号时的交叉相关的原因的特性。因而，通过根据跟踪中的信号频率和如上述那样设定的频率间隔得到的容易检测交叉相关的峰相关值的搜索基准频率与作为目前搜索对象的频率的差，决定作为目前搜索对象的频率下的相关值的累积时间。即，根据搜索基准频率与作为目前搜索对象的频率的差，使累积时间变化。由此，对于目的信号的捕捉，能够减小由跟踪中的信号带来的交叉相关的影响。进而，通过适当设定累积时间，与单单将累积时间设为一定的情况相比，能够缩短遍及频率带域全域的总的累积时间。
此外，在该发明的信号搜索方法中，累积时间决定工序具有频率区间决定工序、区分判定工序和决定工序。在频率区间决定工序中，根据搜索基准频率与搜索对象的频率的频率差，决定搜索对象的频率能够有选择地对应的多个频率区间。在区分判定工序中，判定搜索对象的频率属于多个频率区间的哪个。在决定工序中，按每个频率区间设定累积时间。
在该方法中，将搜索对象的频率分为多个频率区间，按每个频率区间决定累积时间。由此，相对于按照单独的搜索频率详细地决定累积时间的情况，能够不怎么使信号搜索性能降低而使累积时间的决定工序简单化、缩短。
发明的效果
根据该发明，是比以往方法简单的方法，并且能够比以往更正确地捕捉目的信号。
附图说明
图1是表示专利文献1所记载的信号搜索方法的概念的图。
图2是本发明的实施方式的信号搜索方法的流程图。
图3是用来说明相关值计算的码相位、频率设定的概念的图。
图4表示捕捉到目的的GPS信号的情况下的峰相关值Peak(i)的频率特性及交叉相关的情况下的峰相关值Peak(i)的频率特性。
图5是表示目的的GPS信号的相关特性及交叉相关的相关特性的图。
图6是表示目的的GPS信号及交叉相关的判定流程的流程图。
图7是表示使用平均值及标准偏差的目的的GPS信号及交叉相关的判定流程的流程图。
图8是表示切换使用两种最大噪声值的情况下的目的的GPS信号及交叉相关的判定流程的流程图。
图9是本实施方式的累积相关值计算方法的流程图。
图10是用来说明频率区间的设定概念的图。
图11是表示累积时间设定用映射表的一例的图。
图12是层叠时间变更用映射表的形成工序的流程图。
图13是表示频率区间的判定处理及累积时间的设定处理的流程图。
图14是用来说明频率区间的判定概念及累积时间的设定概念的图。
图15是用来说明信号强度区分的设定概念的图。
图16是表示也包括信号强度区分的累积时间设定用映射表的一例的图。
图17是表示有关本发明的实施方式的GPS信号接收装置1的结构的块图。
图18是表示具备GPS信号接收装置1的信息设备终端100的结构的块图。
具体实施方式
参照附图有关本发明的实施方式的信号搜索方法进行说明。另外，在本实施方式中，对搜索从GPS卫星发送的GPS信号的L1波(以下，单称作GPS信号)的方法进行说明。但是，并不限定于GPS信号的L1波，在捕捉以一定的码周期进行了码调制的信号的情况下，可以采用以下的方法。
图2是本实施方式的信号搜索方法的流程图。图3是用来说明相关值计算的码相位、频率设定的概念的图。
如图2所示，在本实施方式的信号搜索方法中，首先，对以GPS信号的频率1575.42MHz为扫描频带的大致中心的频率的8kHz的频率区间，进行接收信号和副本信号的相关处理。所谓副本信号，是将在GPS信号中使用的C/A码再现的信号。并且，在本实施方式的信号搜索方法中，通过按照各码相位及频率的组合、遍及规定时间长将相关值累积，来计算累积相关值。
进而，在本实施方式的信号搜索方法中，如图3所示，对于码相位轴方向，对于1个频率的1023码片的码相位，以0.5码片的分辨率，将全部的码相位点同时相关处理，按照码相位点计算累积相关值。将这样以1个频率得到的累积相关值群称作频率上的累积相关值。
此外，在本实施方式的信号搜索方法中，对于频率轴方向，将以1000Hz间隔设定的8个频率作为1组，同时并行执行相关处理而计算累积相关值。作为具体的例子，如图3所示，在由1000Hz间隔构成的频率Fsig011、Fsig012、Fsig013、Fsig014、Fsig015、Fsig016、Fsig017、Fsig018的组Gr1中，同时计算累积相关值(S101)。
接着，在相对于组Gr1使频率变动了50Hz的组Gr2(包含频率Fsig021的组)中，在8点同时执行相关处理。
接着，在相对于组Gr2使频率变动了50Hz的组Gr3(包含频率Fsig031的组)中，在8点同时执行相关处理。
通过将这样的相关处理从组Gr1到组Gr19依次执行，得到包罗扫描频带的全域(8kHz)的1个频道的累积相关值。并且，如果1个频道的相关处理结束，则再次将从组Gr1起的相关处理依次反复执行。另外，组Gr间的转移也可以与后述的交叉相关的判定同时并行进行，以下，表示在交叉相关判定中的组Gr中进行交叉相关判定的情况下、换言之在没有捕捉到目的的信号的情况下向下个组Gr转移的处理。
接着，按照构成1个组Gr的频率i(i＝1～8)，将在该频率i上存在的全部的累积相关值比较。并且，取得作为最大值的累积相关值作为该频率i的峰相关值Peak(i)(S102)。
这样取得的各个频率i的峰相关值Peak(i)成为图4所示那样的特性。图4(A)表示捕捉到目的的GPS信号的情况下的峰相关值Peak(i)的频率特性，图4(B)表示交叉相关的情况下的峰相关值Peak(i)的频率特性。
在捕捉到目的的GPS信号的情况下，特定的频率的峰相关值相比其他频率的峰相关值大幅变高。具体而言，如果是图4(A)的例子，则频率Fsig012的频率i＝2的峰相关值Peak(2)相比组Gr1的其他频率Fsig011，Fsig013～Fsig018的频率i＝1，3～8的峰相关值Peak(1)，Peak(3)～Peak(8)大幅变大。
在交叉相关的情况下，所取得的全部的频率i＝1～8的峰相关值Peak(1)～Peak(8)大致相同。
这是因为，对于目的的GPS信号的相关特性和交叉相关的相关特性具有图5所示那样的频率特性。图5(A)表示目的的GPS信号的相关特性，图5(B)是表示交叉相关的相关特性的图。
如图5(A)所示，在目的的GPS信号的情况下，累积相关值仅在目的的GPS信号的接收时的频率下累积相关值变得非常高。另一方面，在其他频率下累积相关值大幅变低。
如图5(B)所示，在交叉相关的情况下，累积相关值以1000Hz间隔为极大值，成为大致相同的累积相关值。这是因为，C/A码由1023码片构成，比特率是1.023Mbps。即，因为C/A码的码周期是1msec.(毫秒)，码周期是1000Hz。
利用这样的相关特性，在本实施方式的信号搜索方法中，接着以以下所示的流程进行目的的GPS信号的捕捉及交叉相关的判定。
如上述那样，将所取得的各频率i的峰相关值Peak(i)比较，进行交叉相关的判定(S103)。图6是表示目的的GPS信号及交叉相关的判定流程的流程图。
首先，将各频率i的峰相关值Peak(i)比较，检测它们的最大值作为最大峰值Tpeak(S301)。
接着，将各频率i的峰相关值Peak(i)比较，将其中第2高的值的峰相关值Peak(i)设定为最大噪声值Np(S302)。换言之，在各频率i的峰相关值Peak(i)中，将除了最大峰值Tpeak以外为最大值的峰相关值Peak(i)设定为最大噪声值Np。
接着，通过Tpeak/Np的运算式计算最大峰值Tpeak与最大噪声值Np的比Ratio(S303)。
接着，将预先设定的阈值Th与比Ratio比较，如果比Ratio比阈值Th高(S304：是)，则判定最大峰值Tpeak是由目的的GPS信号带来的(S306)。另一方面，如果比Ratio是阈值Th以下(S304：否)，则判定这些峰相关值Peak(i)群是由交叉相关带来的(S305)。
具体而言，在图4(A)的例子中，将频率Fsig012的峰相关值Peak(2)设定为最大峰值Tpeak。将频率Fsig014的峰相关值Peak(4)设定为最大噪声值Np。比Ratio(Tpeak/Np)为Peak(2)/Peak(4)。这里，如上述那样，根据GPS信号的相关特性，峰相关值Peak(2)相比峰相关值Peak(4)大幅变高。因而，比Ratio为非常高的正值。
在图4(B)的例子中，也将频率Fsig012的峰相关值Peak(2)设定为最大峰值Tpeak。将频率Fsig014的峰相关值Peak(4)设定为最大噪声值Np。比Ratio(Tpeak/Np)为Peak(2)/Peak(4)。但是，如上述那样，在交叉相关的情况下，峰相关值Peak(2)和包括峰相关值Peak(4)的其他峰相关值Peak为大致相同的值。因而，比Ratio为约1左右。
因此，如果将阈值Th设定为例如2～3左右的规定的正值，则在GPS信号的情况下，比Ratio变得比阈值Th高，能够正确地判定最大峰值Tpeak是由目的的GPS信号带来的。另一方面，在交叉相关的情况下，比Ratio变得比阈值Th低，能够正确地判定为交叉相关。
这样，判定是目的的GPS信号还是交叉相关，如果能够判定为从目的的GPS信号得到了累积相关值(S104：否)，则通过取得得到了作为最大峰值Tpeak的累积相关值的频率和码相位，进行目的的GPS信号的捕捉(S105)，转移到跟踪处理。
这里，如果判定为交叉相关(S104：是)，则将构成组的各频率分别各变动50Hz，即，将组移动，同样进行以上述累积相关值的取得为开始的交叉相关的判定处理。将该处理例如持续地反复进行，直到捕捉到目的的GPS信号。但是，在即使将扫描频带的全域扫描规定次数(例如2、3次)也没有得到目的的GPS信号的情况下，也可以进行例如进行表示没有捕捉到目的的GPS信号的显示等的处理。另外，在本实施方式中表示了进行50Hz的频率变动的例子，但也可以根据执行运算的资源量及需要的捕捉速度等而进行其他频率(例如100Hz)的频率变动。
如以上这样，如果使用本实施方式的信号搜索方法，则能够正确地判定交叉相关，能够正确且可靠地捕捉目的的GPS信号。
进而，通过使用本实施方式的方法，仅通过对离散设定的各频率计算最大的累积相关值(分频率最大相关值)的工序、和从多个分频率最大相关值再检测交叉相关判定用的最大累积相关值和第2高的累积相关值，就能够进行交叉相关的判定。即，不需要如以往技术所示那样对码相位轴的整体进行两次最大值检测，使处理简单化。例如在将1023码片的码相位以0.5码片的分辨率扫描的情况下，如果是以往例，则必须进行1023×2×2＝4092次的累积相关值的取得，在将其以8个不同的频率取得的情况下，需要4092×8＝32736次的累积相关值的取得。另一方面，在本实施方式的方法中，对于一个频率是1023×2＝2046次，在8个频率下为2046×8＝16368次。并且，为了从8个频率的累积相关值决定最大峰相关值，为了决定8次第2高的峰相关值而为7次的扫描。因而，为16368+8+7＝16383次，为以往的32736次的一半左右就足够。而且，在本实施方式中，由于在独立的8个频率下平行地执行对累积相关值的各处理，所以能够实现更高速的处理。
此外，在以往技术中，为了找到在码相位轴上存在的第2高的相关峰值，必须进行码相位轴上的范围限制，而在本实施方式中，由于根据相互独立的多个频率的峰相关值间决定第2高的峰相关值，所以不需要设置限制，能够进行比以往技术简单的处理。
此外，在已经取得导航电文、根据该导航电文得到了目的的GPS信号的信息的情况下，不论是以往技术还是本实施方式的方法都能够缩小码相位轴上的扫描范围，但在以往技术的情况下，如上述那样考虑到真的第2高的峰相关值不在扫描范围内的情况，有交叉相关的判定准确度下降的情况。但是，在本实施方式的方法中，由于在各频率的码相位轴上取得最大的累积相关值，根据这些相互独立的多个频率中的最大的累积相关值得到最大峰值和第2高的峰相关值(最大噪声值)，所以能够得到为了进行交叉相关判定而需要的真的第2高的峰相关值(累积相关值)，不会使交叉相关判定准确度下降。
如以上这样，如果使用本实施方式的方法，是比以往技术简单的处理，并且能够正确地判定交叉相关，能够正确地捕捉目的的GPS信号。
另外，在上述说明中，表示了使最大噪声值Np为第2个的峰相关值的例子，但也可以如以下所示那样，根据各频率i的峰相关值Peak(i)的平均值及标准偏差设定最大噪声值Np'。图7是表示使用平均值及标准偏差的目的的GPS信号及交叉相关的判定流程的流程图。
首先，将各频率i的峰相关值Peak(i)比较，检测这些最大值作为最大峰值Tpeak(S311)。
接着，计算同时取得的全部频率i(i＝1～8)的累积相关值的平均值E(S312)。例如，在通过上述8频率同时处理以对于1023码片的0.5码片分辨率取得了累积相关值的情况下，通过将1023×2×8＝16368个累积相关值相加、将该加算值用模数16368去除，计算平均值E。
接着，计算同时取得的全部频率i(i＝1～8)的累积相关值的标准偏差σ(S313)。例如，在通过上述8频率同时处理以对于1023码片的0.5码片分辨率取得了累积相关值的情况下，计算1023×2×8＝16368个累积相关值的平方和，将该平方和用模数16368去除，减去上述平均值E的平方。并且，通过取该减法结果的平方根，计算标准偏差σ。
接着，根据平均值E和标准偏差σ计算最大噪声值Np'(S314)。最大噪声值Np'是对标准偏差σ乘以常数N、通过乘法结果与平均值E的相加来计算。即，根据Np'＝E+Nσ得到。这里，N是依据沿着码相位轴方向的扫描范围设定的常数。
例如，扫描范围越是扩大则N越小，扫描范围越窄则N设定得越大。作为具体例，在将码相位的全码范围(1023码片)作为扫描范围的情况下设定N＝4，在将码相位的前码范围的1/4(256码片)作为扫描范围的情况下设定N＝5。
通过根据扫描范围这样设定常数N，能够抑制由样本数(扫描范围的大小)带来的对于噪声计算误差的影响，即使使扫描范围变化也能够设定可靠性较高的最大噪声值Np'。
接着，通过Tpeak/Np'的运算式计算最大峰值Tpeak与最大噪声值Np'的比Ratio'(S315)。
接着，将预先设定的阈值Th'与比Ratio'比较，如果比Ratio'比阈值Th'高(S316：是)，则判定最大峰值Tpeak是由目的的GPS信号带来的(S318)。另一方面，如果比Ratio'是阈值Th'以下(S316：否)，则判定这些峰相关值Peak(i)群是由交叉相关带来的(S317)。
这样，如果使用设定使用平均值E和标准偏差σ的最大噪声值Np'的方法，则能够使最大噪声值的统计学上的可靠性。因而，能够更正确地判定交叉相关。
另外，也可以将上述两种最大噪声值的设定方法组合使用。图8是表示将两种最大噪声值切换使用的情况下的目的的GPS信号及交叉相关的判定流程的流程图。
首先，将各频率i的峰相关值Peak(i)比较，检测这些最大值作为最大峰值Tpeak(S321)。
取得码相位的扫描范围，在码相位的扫描范围较大(1023码片的全码等)情况下(S322：是)，将各频率i的峰相关值Peak(i)比较，将它们之中第2高的值的峰相关值Peak(i)设定为最大噪声值Np(S323)。
接着，通过Tpeak/Np的运算式计算最大峰值Tpeak与最大噪声值Np的比Ratio(S324)。
接着，将预先设定的阈值Th与比Ratio比较，如果比Ratio比阈值Th高(S325：是)，则判定最大峰值Tpeak是由目的的GPS信号带来的(S326)。另一方面，如果比Ratio是阈值Th以下(S325：否)，则判定这些峰相关值Peak(i)群是由交叉相关带来的(S327)。
取得码相位的扫描范围，在码相位的扫描范围较窄(1023码片的全码的1/4(256码片)等)情况下(S322：否)，计算同时取得的全部的频率i(i＝1～8)的累积相关值的平均值E和标准偏差σ(S328)。这些累积相关值的平均值E及标准偏差σ的计算方法使用与上述方法相同的方法。
接着，根据平均值E和标准偏差σ计算最大噪声值Np'(S329)。最大噪声值Np'是对标准偏差σ乘以常数N，通过乘法结果与平均值E的相加来计算。即，根据Np'＝E+Nσ得到。这里，N是依据沿着码相位轴方向的扫描范围设定的常数。
接着，通过Tpeak/Np'的运算式计算最大峰值Tpeak与最大噪声值Np'的比Ratio'(S330)。
接着，将预先设定的阈值Th'与比Ratio'比较，如果比Ratio'比阈值Th'高(S331：是)，则判定最大峰值Tpeak是由目的的GPS信号带来的(S326)。另一方面，如果比Ratio'是阈值Th'以下(S331：否)，则判定这些峰相关值Peak(i)群是由交叉相关带来的(S327)。
在使用这样的两个最大噪声值的设定方法的组合的情况下，在模数较多的情况下，由于使用第2高的值的峰相关值，所以对于最大噪声值的设定不易受到因模数的较多带来的影响。进而，由于模数较多，所以即使用第2高的值的峰相关值设定最大噪声值，该值也具有充分的可靠性。另一方面，在模数较少的情况下，通过使用平均值E和标准偏差σ的统计学的方法，由于设定最大噪声值，所以能够抑制因模数的较少带来的最大噪声值的可靠性下降。进而，由于模数较少，所以几乎不给运算处理负荷的增加带来影响。因而，通过做成使用两个最大噪声值的设定方法的组合的方法，能够维持最大噪声值的可靠性，并且几乎不增加运算处理负荷而进行交叉相关的判定。
在上述信号搜索方法中，没有特别详细地定义用来计算相关值的累积时间，不取决于频率而设为一定。但是，在已经存在跟踪中的GPS信号的情况下，也可以如以下所示那样按照组单独设定累积时间。
参照附图对有关本发明的实施方式的信号搜索方法进行说明。另外，在本实施方式中，对从GPS卫星发送的GPS信号的搜索方法进行说明。此外，在本实施方式中，为了使说明简略化，对在本实施方式所示的信号搜索的开始时点、跟踪中的GPS信号是1个的情况进行说明。
图9是本实施方式的累积相关值计算方法的流程图。如图9所示，在本实施方式的信号搜索方法中，首先，取得跟踪中的GPS信号的频率FB，设定为搜索基准频率FSB(S601)。搜索基准频率FSB由所取得的跟踪中的GPS信号的频率FB、和以该频率FB为基准设定为存在于上述8kHz的频率范围内的1000Hz间隔的频率构成。该搜索基准频率FSB相当于在本实施方式中表示的信号搜索时的发生因交叉相关带来的累积相关值的极大值的频率。
接着，利用这样的交叉相关特有的累积相关值特性，基于搜索基准频率FSB，如图10所示那样决定频率区间(S602)。图10是用来说明频率区间的设定概念的图。
如图10所示，在本实施方式的信号搜索方法中，以搜索基准频率FSB为基准的频率(0Hz)，将从该基准的频率离开500Hz的频率范围划分为第1频率区间ARc、第2频率区间ARn及第3频率区间ARf而设定。
第1频率区间ARc在概念上，是最容易受交叉相关的影响、相对于交叉相关的频率极接近的频率区间。第1频率区间ARc以从基准的频率到第1阈值频率Fc的频率带域设定。
第2频率区间ARn在概念上，是在某种程度上有受到交叉相关的影响的可能性、相对于交叉相关的频率比极接近区间远离但某种程度上接近的频率区间。第2频率区间ARn以从第1阈值频率Fc到第2阈值频率Fn(<Fc)的频率带域设定。
第3频率区间ARf在概念上，是几乎不受交叉相关的影响、相对于交叉相关的频率离开的频率区间。第3频率区间ARf以从第2阈值频率Fn到500Hz的频率带域设定。
另外，第1阈值频率Fc、第2阈值频率Fn基于累积时间适当设定。例如，设定第1阈值频率Fc，以使与交叉相关的主波束大致对应的频率带域为第1频率区间ARc。此外，设定第2阈值频率Fn，以使与相邻于该主波束的侧波束大致对应的频率带域为第2频率区间ARn。
并且，如果累积时间变化，则构成交叉相关的相关特性的各波束的宽度变化，所以只要对应于该相关特性的变化而使第1阈值频率Fc、第2阈值频率Fn变化就可以。例如，如果累积时间变长，则主波束及侧波束的宽度变窄，所以只要将第1阈值频率Fc、第2阈值频率Fn设定为更小的值就可以。
接着，在本实施方式的信号搜索方法中，取得通过预备搜索检测出的预备检测信号的频率，进行图11所示那样的累积时间变更用的映射表形成(S603)。图11是表示累积时间设定用映射表的一例的图。此外，图12是层叠时间变更用映射表的形成工序的流程图。
如果取得通过预备搜索检测出的预备检测信号的频率F，则计算该预备检测信号的频率F与搜索基准频率FSB的差分值(频率差分值Df(F))。此时，如果有多个预备检测信号，则按照各个预备检测信号计算频率差分值Df(F)(S701)。
预备检测信号的频率差分值Df(F)可以根据下式得到。Df(F)＝(ABS(F－FSB))/1000[Hz]这里，ABS()是表示绝对值运算的符号。
接着，进行标准化，以使预备检测信号的频率差分值为从0Hz到500Hz之间的值，计算标准化频率差分值Δf(F)。预备检测信号的标准化频率差分值Δf(F)根据下式得到。另外，Mod()是表示余数的符号。如果Mod(Df(F))≥500[Hz]则Δf(F)＝ABS(Df(F)－1000)此外Mod(Df(F))<500[Hz]则Δf(F)＝Mod(Df(F))将这样计算出的各预备检测信号的标准化频率差分值Δf(F)与第1阈值频率Fc、第2阈值频率Fn比较(S702)。
如果预备检测信号的标准化频率差分值Δf(F)比第1阈值频率Fc低(S703：是)，则判定该预备检测信号为第1频率区间ARc内(S705)，将第1频率区间ARc的预备检测信号数计数+1(S706)。
如果预备检测信号的标准化频率差分值Δf(F)是第1阈值频率Fc以上(S703：否)，比第2阈值频率Fn低(S704：是)，则将该预备检测信号判定为第2频率区间ARn内(S707)，将第2频率区间ARn的预备检测信号数计数+1(S708)。
如果预备检测信号的标准化频率差分值Δf(F)是第1阈值频率Fc以上(S703：否)、是第2阈值频率Fn以上(S704：否)，则将该预备检测信号判定为第3频率区间ARf内(S709)，将第3频率区间ARf的预备检测信号数计数+1(S710)。
将这样的预备检测信号对应于哪个频率区间的处理对全部的预备检测信号进行(S711：是)，在有没有被进行该判断处理的预备检测信号的情况下(S711：否)，重复从上述步骤S702起的处理，直到对全部的预备检测信号的判断处理结束。
通过进行这样的处理，按照第1频率区间ARc、第2频率区间ARn、第3频率区间ARf形成存在的预备检测信号数的映射表(参照图11)。
接着，基于映射表，设定第1频率区间ARc的累积时间Tc、第2频率区间ARn的累积时间Tn及第3频率区间ARf的累积时间Tf(S604)，通过以设定的累积时间将相关值累积，计算累积相关值(S605)。
这里，如图11所示，在第1频率区间ARc、第2频率区间ARn的两者中有预备检测信号的情况下，设定第1频率区间ARc的累积时间Tc、第2频率区间ARn的累积时间Tn及第3频率区间ARf的累积时间Tf，以成为Tc>Tn>Tf的关系。
通过进行这样的设定，设定为，越是容易受到交叉相关的影响的频率区间则累积时间越长。这里，根据累积时间越长则由交叉相关带来的相关特性的峰越陡峭也可知，使累积时间越长，则频率轴上的信号检测的分辨率越提高。
因而，通过如本实施方式那样设定累积时间，即使是与交叉相关的峰频率极接近的频率带、或与交叉相关的峰频率接近的频率带，也对应于各自的频率区间而适当地降低交叉相关的影响，能够更正确地检测、捕捉作为目的的GPS信号。
这样，如果设定频率区间和累积时间，则对于各搜索频率Fsig判定与上述哪个频率区间对应，决定搜索频率Fsig的累积时间。此时，如上述那样构成组Gr的搜索频率Fsig处于1000Hz间隔，所以包含在1个组Gr中的全部的搜索频率Fsig与搜索基准频率FSB的差相同。
因而，与1个组Gr中包含的全部的搜索频率Fsig对应的累积时间的决定只要由1个搜索频率Fsig进行、将所决定的累积时间也设定为属于相同的组Gr的其他的搜索频率Fsig就可以。图13是表示频率区间的判定处理及累积时间的设定处理的流程图。
计算搜索基准频率FSB与代表组Gr的1个搜索频率Fsig的频率差分值Df，计算标准化频率差分值Δfsig(S801)。具体而言，进行与上述预备检测信号的标准化频率差分值Δf(F)同样的处理。
计算搜索基准频率FSB与代表组Gr的1个搜索频率Fsig的频率差分值Df。
频率差分值Df根据下式得到。Df＝(ABS(Fsig－FSB))/1000[Hz]接着，进行标准化以使频率差分值成为0Hz到500Hz之间的值，计算标准化频率差分值Δfsig。标准化频率差分值Δfsig根据下式得到。如果Mod(Df)≥500[Hz]则Δfsig＝ABS(Df－1000)此外Mod(Df)<500[Hz]则Δfsig＝Mod(Df)将这样计算出的相对于搜索频率Fsig的标准化频率差分值Δfsig与上述第1阈值频率Fc比较(S802)。
如果标准化频率差分值Δfsig比第1阈值频率Fc低(S802：是)，则判定搜索频率Fsig为第1频率区间ARc内，对于搜索频率Fsig采用累积时间Tc(S804)。
如果标准化频率差分值Δfsig是第1阈值频率Fc以上(S802：否)，则将相对于搜索频率Fsig的标准化频率差分值Δfsig与上述第2阈值频率Fn比较(S803)。
如果标准化频率差分值Δfsig比第2阈值频率Fn低(S803：是)，则判定搜索频率Fsig为第2频率区间ARn内，对于搜索频率Fsig采用累积时间Tn(S805)。
如果标准化频率差分值Δfsig是第2阈值频率Fn以上(S803：否)，则判定搜索频率Fsig为第3频率区间ARf内，对于搜索频率Fsig采用累积时间Tf(S806)。
这样，设定与搜索频率Fsig相对于交叉相关的峰频率的离开状况对应的累积值。
具体而言，如上述那样，在使搜索频率按照组以50Hz间隔变动的情况下，如图14那样设定代表各组Gr1～Gr18的搜索频率Fsig011，Fsig021，Fsig031，Fsig041，Fsig051，Fsig061，Fsig071，Fsig081，Fsig091，Fsig101，Fsig111，Fsig121，Fsig131，Fsig141，Fsig151，Fsig161，Fsig171，Fsig181，Fsig191。图14是用来说明频率区间的判定概念及累积时间的设定概念的图。另外，在图14中，表示各组Gr1的搜索频率Fsig011距交叉相关的峰频率最远的例子。
在图14的例子中，搜索频率Fsig011，Fsig021，Fsig031，Fsig041，Fsig051，Fsig061对应于第3频率区间ARf。搜索频率Fsig071，Fsig081对应于第2频率区间ARn。搜索频率Fsig091，Fsig101，Fsig111对应于第1频率区间ARc。搜索频率Fsig121，Fsig131对应于第2频率区间ARn。搜索频率Fsig141，Fsig151，Fsig161，Fsig171，Fsig181，Fsig191对应于第3频率区间ARf。
因而，在搜索频率Fsig011，Fsig021，Fsig031，Fsig041，Fsig051，Fsig061及搜索频率Fsig141，Fsig151，Fsig161，Fsig171，Fsig181，Fsig191中，采用累积时间Tf。在搜索频率Fsig071，Fsig081及搜索频率Fsig121，Fsig131中，采用累积时间Tn。在搜索频率Fsig091，Fsig101，Fsig111中采用累积时间Tc。
以该结果为基础，将对于属于分别包括搜索频率Fsig011，Fsig021，Fsig031，Fsig041，Fsig051，Fsig061及搜索频率Fsig141，Fsig151，Fsig161，Fsig171，Fsig181，Fsig191的组Gr1，Gr2，Gr3，Gr4，Gr5，Gr6，Gr14，Gr15，Gr16，Gr17，Gr18，Gr19的全部搜索频率Fsig的累积时间设定为Tf。
将对于属于分别包括搜索频率Fsig071，Fsig081及搜索频率Fsig121，Fsig131的组Gr7，Gr8，Gr12，Gr13的全部搜索频率Fsig的累积时间设定为Tn。
将对于属于分别包括搜索频率Fsig091，Fsig101，Fsig111的组Gr9，Gr10，Gr11的全部搜索频率Fsig的累积时间设定为Tc。
根据这样设定的累积时间，以各搜索频率进行相关处理，计算与累积时间对应的累积相关值。
通过进行这样的处理，如上述那样，能够根据频率区间适当地减小交叉相关的影响，使信号搜索性能比以往方法提高。
此外，通过使不易受到交叉相关的影响的第2、第3频率区间ARn，ARf的累积时间Tn，Tf比容易受到交叉相关的影响的第1频率区间的累积时间Tc短，与在全频率带域中采用容易受到交叉相关的影响的第1频率区间的累积时间Tc相比，能够缩短1个频道的信号搜索时间。
如以上这样，如果使用本实施方式的累积相关值计算方法，则还能够不产生由交叉相关带来的误捕捉、而以比以往方法短时间进行GPS信号的捕捉。
另外，在上述说明中，表示了根据计数数设定累积时间Tc，Tn，Tf的例子，但也可以不考虑计数数，而按照频率区间，以距交叉相关的峰频率较近的顺序使累积时间变长。如果进行这样的累积时间的设定，则不需要预备搜索。但是，在根据计数数设定累积时间的情况下，在没有预备检测信号的情况下能够进行使第2频率区间ARn的累积时间Tn及第3频率区间ARf的累积时间Tf同等地变短、或使各累积时间均匀地变短等的设定，能够根据状况进行更适当的累积时间的设定。
此外，在上述说明中，表示了根据从交叉相关的峰频率(搜索基准频率FSB)的离开量设定累积时间的例子，但也可以也使用预备检测信号的信号强度设定累积时间。图15是用来说明信号强度区分的设定概念的图。图16是表示也包括信号强度区分的累积时间设定用映射表的一例的图。
如图15所示，信号强度区分根据预备检测信号的C/No(载噪比)而设定为3个级别。具体而言，将不到第1阈值C/N0n设定为第1信号强度区分ZONEw，将第1阈值C/N0n以上不到第2阈值C/N0s设定为第2信号强度区分ZONEn，将第2阈值C/N0s以上设定为第3信号强度区分ZONEs。
在取得了预备检测信号的情况下，按照预备检测信号检测C/N0，判定对应于哪个信号强度区分。
通过进行该处理，能够与上述频率区间的判定结果一起形成图16所示那样的累积时间设定用映射表。
并且，在各频率区间中，参照与各信号强度区分对应的预备检测信号数的分布，设定各频率区间的累积时间。例如，如图16所示，如果判明了在第1频率区间ARc中存在C/N0的较高的预备检测信号，则该预备检测信号是产生交叉相关的峰频率的跟踪中的GPS信号的可能性较高，所以只要将累积时间修正得较长、以便不易受到该跟踪中的GPS信号的影响就可以。
由此，还能够根据接收状况设定更适当的累积时间，能够更加使得不会发生因交叉相关造成的误捕捉。
另外，在上述说明中，以GPS信号的交叉相关的情况为例进行了说明，但对于其他GNSS信号的捕捉也同样能够采用。进而，对于以根据将目的信号调制的扩散码的周期决定的频率间隔在相关值中出现峰值那样的无线通信信号的捕捉也同样能够采用。
此外，在上述频率区间的设定处理中，表示了设定3个频率区间的例子，但也可以设定为两个以上的频率区间。同样，可以将信号强度区分也设定为两个以上的信号强度区分。
此外，在上述说明中没有详细表示，但累积时间只要通过将相干累积时间、非相干累积时间中的某一方或两者设定变更，来实现累积时间的修正就可以。
此外，在上述说明中，表示了将搜索对象频率Fsig分配给多个频率区间、按照频率区间设定累积时间的例子，但也可以根据搜索对象频率Fsig与搜索基准频率FSB的频率差分值，按照搜索对象频率Fsig设定累积时间。在此情况下，只要例如进行设定，以使搜索对象频率Fsig与搜索基准频率FSB的频率差分值越大则使累积时间越短就可以。
以上那样的处理可以通过以下所示的结构的GPS信号接收装置实现。图17是表示有关本发明的实施方式的GPS信号接收装置1的结构的块图。
GPS信号接收装置1具备GPS接收天线10、RF处理部20、基带处理部30及测位运算部40。
GPS接收天线10接收从各GPS卫星广播(发送)的GPS信号，向RF处理部20输出。RF处理部20将接收到的GPS信号降频变换，生成中间频率信号(IF信号)，向基带处理部30输出。
基带处理部30与本发明的具备“相关值取得部”及“判定部”的“信号搜索装置”对应。此外，基带处理部30也与本发明的“捕捉跟踪部”对应。另外，基带处理部30也可以将相当于“判定部”的硬件和相当于“相关值取得部”及“捕捉跟踪部”的硬件，也可以由一体的硬件实现。基带处理部30通过对IF信号乘以载波频率信号而生成基带信号，进行该基带信号的GPS信号的捕捉处理及跟踪处理。此时，在捕捉处理中使用上述信号搜索方法。由此，能够抑制交叉相关的误捕捉，可靠地捕捉作为目的的GPS信号。
这样的对于捕捉到的GPS信号的捕捉处理向跟踪处理转移。将该跟踪得到的码相关结果及载波相关结果、还有根据码相关结果得到的模拟距离向测位运算部40输出。
测位运算部40基于码相关结果将导航电文解调，并根据码相关结果、载波相位结果、模拟距离进行GPS信号接收装置1的测位。
使用这样的结构，通过使用上述信号搜索方法，能够抑制误捕捉且进行高速的捕捉，所以GPS信号的跟踪精度提高，结果能够使测位结果的精度提高。
另外，执行上述信号搜索方法的基带处理部30既可以通过执行各处理的硬件群实现，也可以通过以将上述信号搜索方法的各处理程序化的状态存储在存储媒体中、用计算机将该程序读出而执行的形态实现。
此外，这样的GPS信号接收装置1及GPS信号接收功能在图18所示那样的信息设备终端100中使用。图18是表示具备本实施方式的GPS信号接收装置1的信息设备终端100的主要结构的块图。
图18所示那样的信息设备终端100具备例如便携电话机、汽车导航装置、PND、照相机、时钟等，具备天线10、RF处理部20、基带处理部30、测位运算部40、应用处理部130。天线10、RF处理部20、基带处理部30、测位运算部40是上述结构，由此，如上述那样构成GPS信号接收装置1。
应用处理部130基于从GPS信号接收装置1输出的测位结果显示自装置位置及自装置速度、或执行用于导航等的处理。
在这样的结构中，通过如上述那样得到高精度的测位结果，能够实现高精度的位置显示及导航等。
标号说明
1：GPS信号接收装置；10：GPS接收天线；20：RF处理部；30：基带处理部；40：测位运算部；100：信息设备终端；130：应用处理部

Claims (17)

1.一种信号搜索方法，用来捕捉目的信号，具有：
相关值取得工序，将以根据将上述目的信号调制而得到的扩散码的码周期决定的频率间隔设定的多个频率作为1个组，进行接收信号与以上述多个频率生成的上述目的信号的副本信号的相关处理而取得相关值；以及
判定工序，根据该组内的多个相关值进行是否捕捉到上述目的信号的判定。
2.如权利要求1所述的信号搜索方法，其特征在于，
上述目的信号是GPS信号的L1波；
上述频率间隔是1000Hz间隔。
3.如权利要求1或2所述的信号搜索方法，其特征在于，
上述判定工序中，检测上述组内的各频率的最大相关值，将该最大相关值设定为该频率的峰相关值，基于上述各频率的峰相关值的大小关系，判定是否捕捉到上述目的信号。
4.如权利要求3所述的信号搜索方法，其特征在于，
上述判定工序具有：
检测作为上述各频率的峰相关值的最大值的最大峰值的工序；
根据除了上述最大峰值以外的上述各频率的峰相关值决定最大噪声值的工序；以及
计算上述最大峰值与上述最大噪声值的比的工序；
基于上述比，判定是否捕捉到上述目的信号。
5.如权利要求1～4中任一项所述的信号搜索方法，其特征在于，
当通过上述判定工序判定为没有捕捉到上述目的信号时，使上述多个频率以比上述频率间隔短的一定的频率宽度进行频率变动。
6.如权利要求1～5中任一项所述的信号搜索方法，其特征在于，
上述相关值取得工序具有：
搜索基准频率设定工序，根据跟踪中的信号的频率和上述频率间隔设定搜索基准的频率；
累积时间决定工序，根据上述搜索基准频率与搜索对象的频率的频率差，设定该搜索对象的频率下的累积时间；以及
相关值计算工序，通过以设定的累积时间将上述相关值累积，计算累积相关值。
7.如权利要求6所述的信号搜索方法，其特征在于，
上述累积时间决定工序具有：
频率区间决定工序，根据上述搜索基准频率与搜索对象的频率的频率差，决定上述搜索对象的频率能够有选择地对应的多个频率区间；
区分判定工序，判定上述搜索对象的频率属于上述多个频率区间的哪个；以及
决定工序，按每个上述频率区间设定累积时间。
8.一种信号搜索程序，用来使计算机执行捕捉目的信号的处理，具有：
相关值取得处理，将以根据将上述目的信号调制而得到的扩散码的码周期决定的频率间隔设定的多个频率作为1个组，进行以该多个频率生成的上述目的信号的副本信号与接收信号的相关处理而取得相关值；以及
判定处理，根据该组内的多个相关值判定是否捕捉到上述目的信号。
9.如权利要求8所述的信号搜索程序，其特征在于，
上述目的信号是GPS信号的L1波；
上述频率间隔是1000Hz间隔。
10.如权利要求8或9所述的信号搜索程序，其特征在于，
在上述判定处理中，检测上述组内的各频率的最大相关值，使该最大相关值设定为该频率的峰相关值，基于上述各频率的峰相关值的大小关系，判定是否捕捉到上述目的信号。
11.如权利要求8～10中任一项所述的信号搜索程序，其特征在于，
在上述相关值取得处理中，具有：
搜索基准频率设定处理，根据跟踪中的信号的频率和上述频率间隔设定搜索基准的频率；
累积时间决定处理，根据上述搜索基准频率与搜索对象的频率的频率差，设定该搜索对象的频率下的累积时间；
相关值计算处理，通过以设定的累积时间将上述相关值累积，计算累积相关值；
在上述判定处理中，使用上述层叠相关值作为判定是否捕捉到上述目的信号的相关值。
12.一种信号搜索装置，捕捉目的信号，具备：
相关值取得部，将以根据将上述目的信号调制而得到的扩散码的码周期决定的频率间隔设定的多个频率作为1个组，进行以该多个频率生成的上述目的信号的副本信号与接收信号的相关处理而取得相关值；以及
判定部，根据该组内的多个相关值，判定是否捕捉到上述目的信号。
13.如权利要求12所述的信号搜索装置，其特征在于，
上述目的信号是GPS信号的L1波；
上述频率间隔是1000Hz间隔。
14.如权利要求12或13所述的信号搜索装置，其特征在于，
上述判定部检测上述组内的各频率的最大相关值，将该最大相关值设定为该频率的峰相关值，基于上述各频率的峰相关值的大小关系，判定是否捕捉到上述目的信号。
15.如权利要求12～14中任一项所述的信号搜索装置，其特征在于，
上述相关值取得部具备：
搜索基准频率设定部，根据跟踪中的信号的频率和上述频率间隔设定搜索基准的频率；
累积时间决定部，根据上述搜索基准频率与搜索对象的频率的频率差，设定该搜索对象的频率下的累积时间；
相关值计算部，通过以设定的累积时间将上述相关值累积，计算累积相关值；
上述判定部使用上述层叠相关值作为判定是否捕捉到上述目的信号的相关值。
16.一种GNSS信号接收装置，具备：
权利要求12～15中任一项所述的信号搜索装置；
捕捉跟踪部，根据上述相关值捕捉、跟踪上述目的的信号；以及
测位运算部，基于跟踪结果进行测位。
17.一种信息设备终端，具备：
权利要求16所述的GNSS信号接收装置；以及
应用处理部，使用上述测位运算部的测位运算结果执行规定的应用。
CN201380015191.2A 2012-03-22 2013-02-18 信号搜索方法、信号搜索装置、gnss信号接收装置及信息设备终端 Expired - Fee Related CN104204850B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2012064634 2012-03-22
JP2012-064634 2012-03-22
PCT/JP2013/053890 WO2013140911A1 (ja) 2012-03-22 2013-02-18 信号サーチ方法、信号サーチプログラム、信号サーチ装置、ｇｎｓｓ信号受信装置、および情報機器端末

Publications (2)

Publication Number Publication Date
CN104204850A true CN104204850A (zh) 2014-12-10
CN104204850B CN104204850B (zh) 2016-08-17

Family

ID=49222376

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201380015191.2A Expired - Fee Related CN104204850B (zh) 2012-03-22 2013-02-18 信号搜索方法、信号搜索装置、gnss信号接收装置及信息设备终端

Country Status (5)

Country Link
US (1) US9231652B2 (zh)
EP (1) EP2829895B1 (zh)
JP (1) JP5918351B2 (zh)
CN (1) CN104204850B (zh)
WO (1) WO2013140911A1 (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104155662B (zh) * 2014-08-05 2016-08-24 中国空间技术研究院 基于gnss相关峰值探测器的自适应互干扰抑制方法
CN107422341B (zh) * 2016-05-23 2020-02-07 大唐半导体设计有限公司 一种实现信号捕获的方法和装置
CN110068841A (zh) * 2019-05-06 2019-07-30 西安开阳微电子有限公司 一种卫星信号高灵敏度快速捕获方法
CN110501729B (zh) * 2019-06-18 2023-03-31 山东大学 一种基于fpga分步码相位细化的gnss信号的捕获方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6047016A (en) * 1997-06-23 2000-04-04 Cellnet Data Systems, Inc. Processing a spread spectrum signal in a frequency adjustable system
US6282231B1 (en) * 1999-12-14 2001-08-28 Sirf Technology, Inc. Strong signal cancellation to enhance processing of weak spread spectrum signal
US6810072B1 (en) * 2000-05-30 2004-10-26 Nokia Corporation System for acquiring spread spectrum signals
US7110435B1 (en) * 1999-03-15 2006-09-19 Parkervision, Inc. Spread spectrum applications of universal frequency translation
US20090207076A1 (en) * 2008-02-19 2009-08-20 Honeywell International Inc. System and method for gnss position aided signal acquisition
US20100117895A1 (en) * 2008-11-11 2010-05-13 Seiko Epson Corporation Position calculating method and position calculating device
CN102193097A (zh) * 2010-02-23 2011-09-21 精工爱普生株式会社 信号捕捉方法及信号捕捉装置

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003098244A (ja) * 2001-09-26 2003-04-03 Japan Radio Co Ltd 偽信号相互相関検出方法、送信源選択制限方法及び衛星選択制限方法
JP3792145B2 (ja) * 2001-09-28 2006-07-05 日本無線株式会社 相関値識別方法
US7209076B2 (en) 2002-07-10 2007-04-24 Qualcomm Incorporated Cross-correlation mitigation method and apparatus for use in a global positioning system receiver
MXPA05000474A (es) * 2002-07-10 2005-04-19 Qualcomm Inc Metodo y aparato de mitigacion de correlacion cruzada para utilizarse en un receptor del sistema de posicionamiento global.
US7161977B1 (en) 2003-01-28 2007-01-09 Trimble Navigation Limited Receiver having a ratio-based signal acquisition method
US7443339B2 (en) * 2006-01-26 2008-10-28 Qualcomm Incorporated Cross-correlation suppression technique for position location receivers
EP1916535B1 (en) * 2006-10-26 2015-11-18 Qualcomm Incorporated GNSS receiver with cross-correlation rejection
JP4888110B2 (ja) * 2006-12-26 2012-02-29 セイコーエプソン株式会社 相関演算制御回路及び相関演算制御方法
JP5483750B2 (ja) * 2009-11-30 2014-05-07 古野電気株式会社 不要信号判別装置、不要信号判別方法、不要信号判別プログラム、ｇｎｓｓ受信装置および移動端末

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6047016A (en) * 1997-06-23 2000-04-04 Cellnet Data Systems, Inc. Processing a spread spectrum signal in a frequency adjustable system
US7110435B1 (en) * 1999-03-15 2006-09-19 Parkervision, Inc. Spread spectrum applications of universal frequency translation
US6282231B1 (en) * 1999-12-14 2001-08-28 Sirf Technology, Inc. Strong signal cancellation to enhance processing of weak spread spectrum signal
US6810072B1 (en) * 2000-05-30 2004-10-26 Nokia Corporation System for acquiring spread spectrum signals
US20090207076A1 (en) * 2008-02-19 2009-08-20 Honeywell International Inc. System and method for gnss position aided signal acquisition
US20100117895A1 (en) * 2008-11-11 2010-05-13 Seiko Epson Corporation Position calculating method and position calculating device
CN102193097A (zh) * 2010-02-23 2011-09-21 精工爱普生株式会社 信号捕捉方法及信号捕捉装置

Also Published As

Publication number Publication date
EP2829895A4 (en) 2015-12-02
EP2829895A1 (en) 2015-01-28
JPWO2013140911A1 (ja) 2015-08-03
US9231652B2 (en) 2016-01-05
EP2829895B1 (en) 2019-11-27
WO2013140911A1 (ja) 2013-09-26
CN104204850B (zh) 2016-08-17
US20150049789A1 (en) 2015-02-19
JP5918351B2 (ja) 2016-05-18

Similar Documents

Publication Publication Date Title
CN102608626A (zh) 一种高灵敏度卫星导航信号捕获方法及系统
US9829581B2 (en) Enhancing search capacity of global navigation satellite system (GNSS) receivers
US20090153397A1 (en) Gnss satellite signal interference handling method and correlator implementing the same
CN104204850A (zh) 信号搜索方法、信号搜索程序、信号搜索装置、gnss信号接收装置及信息设备终端
CN103954977A (zh) 一种gnss欺骗干扰感知方法和系统
JP2014228537A (ja) マルチモード測位のための受信機と方法
KR100953820B1 (ko) 지피에스 신호의 록 상태 검출 방법 및 장치
CN102037658A (zh) 码空间搜索中的多相关处理
CN106526632A (zh) 一种gnss捕获相干器和信号捕获方法
CN103837877B (zh) 一种卫星识别的方法和装置
US11585943B2 (en) Detection and elimination of GNSS spoofing signals with PVT solution estimation
CN107820212B (zh) 一种基于移动多媒体广播系统实现的定位方法及定位终端
CN105699992B (zh) 高动态gnss载波信号跟踪方法及系统
CN105301611A (zh) 用于将卫星信号分八频点处理的基带部分结构
JP4859790B2 (ja) Ｇｐｓ受信機
JP2008209287A (ja) 衛星航法受信機
EP1808967A1 (en) A method for code alignment for DSSS signal processing
CN101262246A (zh) 用于获取和追踪展频信号的方法及装置
Li et al. Highly sensitive weak signal acquisition method for GPS/compass
JP2013228380A (ja) ナビゲーションのための方法およびシステム
WO2013140910A1 (ja) 信号サーチ方法、信号サーチプログラム、信号サーチ装置、ｇｎｓｓ信号受信装置、および情報機器端末
JP6118535B2 (ja) 測位信号サーチ装置、測位装置、測位信号サーチ方法、および、測位信号サーチプログラム
RU2310212C1 (ru) Цифровой коррелятор приемника сигналов спутниковых радионавигационных систем
CN116774256B (zh) 一种信号频点捕获方法、系统、电子装置及存储介质
CN116953748A (zh) 多路径环境中的直接及反射导航卫星信号参数的检测及估计

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
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: 20160817

Termination date: 20220218

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