WO2013140910A1 - 信号サーチ方法、信号サーチプログラム、信号サーチ装置、ｇｎｓｓ信号受信装置、および情報機器端末 - Google Patents

信号サーチ方法、信号サーチプログラム、信号サーチ装置、ｇｎｓｓ信号受信装置、および情報機器端末 Download PDF

Info

Publication number
WO2013140910A1
WO2013140910A1 PCT/JP2013/053889 JP2013053889W WO2013140910A1 WO 2013140910 A1 WO2013140910 A1 WO 2013140910A1 JP 2013053889 W JP2013053889 W JP 2013053889W WO 2013140910 A1 WO2013140910 A1 WO 2013140910A1
Authority
WO
WIPO (PCT)
Prior art keywords
value
peak
signal
maximum
correlation
Prior art date
Application number
PCT/JP2013/053889
Other languages
English (en)
French (fr)
Inventor
明弘 大杉
一登 多田
佐藤　美喜
Original Assignee
古野電気株式会社
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 古野電気株式会社 filed Critical 古野電気株式会社
Publication of WO2013140910A1 publication Critical patent/WO2013140910A1/ja

Links

Images

Classifications

    • 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/29—Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
    • 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

Definitions

  • the present invention relates to a signal search method for searching for a desired signal from a received signal, and more particularly to a method for searching for a GPS signal in GNSS.
  • GPS Global Positioning System
  • GNSS Global Navigation Satellite System
  • GPS GPS signals transmitted from a plurality of GPS satellites are received, and the receiver performs positioning using the code phase and carrier phase of the received GPS signals.
  • a different spreading code is set for each GPS satellite, and each GPS signal is code-modulated with a different spreading code.
  • a GPS signal from a GPS satellite different from the target GPS satellite may be erroneously captured as the target GPS signal and tracking processing may be performed.
  • Such a phenomenon is called cross-correlation.
  • Patent Document 1 compares the accumulated correlation value of each code phase point on the code phase at a single search frequency so that the obtained accumulated correlation value depends on the target GPS signal. It is determined whether it is due to cross-correlation.
  • FIG. 1 is a diagram illustrating a concept of a signal search method described in Patent Document 1. In FIG.
  • an integrated correlation value at each code phase of a single search frequency is acquired, and these values are compared.
  • an integrated correlation value Vpk1 that is the maximum value among integrated correlation values in all code phases is detected.
  • the second highest integrated correlation value Vpk2 is detected in the integrated correlation values in all code phases.
  • a ratio between the maximum correlation value Vpk1 and the second highest correlation value Vpk2 is calculated.
  • the ratio is equal to or greater than a predetermined threshold
  • the code phase Cpk1 corresponding to the integrated correlation value Vpk1 that is the maximum value is captured by the target signal, and the signal of the code phase Cpk1 is captured.
  • the ratio is less than or equal to the threshold value, the signal search is continued.
  • the second highest integrated correlation value is the code phase on the way to the maximum including the maximum integrated correlation value Vpk1 (in the example of FIG. 1, the code phase Cpk2w). There is. Therefore, the second highest integrated correlation value must be detected so that the predetermined code phase range including the code phase Cpk1 of the integrated correlation value Vpk1 is not included, and the processing becomes complicated.
  • the true second highest accumulated correlation value Vpk2 The accumulated correlation value Vpk2 ′ having a lower value may be detected.
  • the ratio to the maximum integrated correlation value Vpk1 is larger than the true ratio. For this reason, in spite of the cross-correlation, an erroneous determination that is determined as the target GPS signal is likely to occur.
  • an object of the present invention is to provide a signal search method that can capture a target signal more accurately than in the prior art by a relatively simple process compared to the conventional method.
  • the signal search method of the present invention has the following steps.
  • a correlation value is acquired by performing correlation processing between the received signal and the replica signal of the target signal generated at a plurality of frequencies.
  • the peak correlation value acquisition step a peak correlation value that is the maximum correlation value at each frequency is acquired.
  • the determination step it is determined whether or not the target signal has been captured based on the magnitude relationship between the peak correlation values of the respective frequencies.
  • the target signal is captured, for example, whether it is cross-correlation or the target signal.
  • a correlation value peak occurs at a predetermined frequency interval (specifically, an interval according to the frequency of the spreading code).
  • the true signal uses the fact that only the peak correlation value is prominently high.
  • a maximum peak value that is a maximum value of the peak correlation value of each frequency is detected, and a maximum noise value is determined from the peak correlation value of each frequency excluding the maximum peak value.
  • a step of calculating a ratio between the maximum peak value and the maximum noise value is determined whether the target signal has been captured based on the calculated ratio.
  • This method shows the specific contents of the method for determining whether or not the target signal has been captured. Since the correlation value of the cross correlation and the correlation value of the target signal have the characteristics as described above, the maximum peak value is acquired, and the maximum peak value is obtained from the peak correlation values of other frequencies acquired at the same time. If it is compared with the maximum noise value obtained, it can be determined whether the signal is a target signal or cross-correlation.
  • the maximum peak correlation value in the peak correlation value of each frequency excluding the maximum peak value is used.
  • This method shows a specific method for determining the maximum noise value.
  • the peak correlation value having the second highest correlation value following the maximum peak value is set as the maximum noise value.
  • the maximum noise value can be set by simpler processing.
  • a step of calculating an average value and a standard deviation of correlation values of a plurality of frequencies, and a constant set based on a code phase search range are determined.
  • a step of performing In the step of determining the maximum noise value, the standard noise is multiplied by a constant, and the average value is added to the multiplication result to determine the maximum noise value.
  • the maximum noise value with higher reliability can be set statistically by using the average value and the standard deviation.
  • the step of determining the maximum noise value of the signal search method of the present invention includes a first determination step of setting the maximum peak correlation value in the peak correlation value of each frequency excluding the maximum peak value as the maximum noise value. Further, the step of determining the maximum noise value includes a step of calculating an average value and a standard deviation of correlation values of a plurality of frequencies, and a step of determining a constant set based on a code phase search range. And a second determination step of determining a maximum noise value by multiplying the standard deviation by a constant and adding an average value to the multiplication result.
  • the maximum noise value is determined in the first determination step when the code phase range for calculating the correlation value is wide, and the second is determined when the code phase range for calculating the correlation value is narrow.
  • the maximum noise value is determined in the determining step.
  • code phase scanning includes a method of setting a maximum noise value with a peak correlation value having the second highest correlation value following the above-described maximum peak value, and a method of setting a maximum noise value from an average value and a standard deviation. Depending on the range, they are switched appropriately.
  • a plurality of frequencies are frequencyd with a constant frequency width shorter than the frequency interval at which the correlation value is simultaneously acquired. Shift and repeat the determination of whether or not the target signal has been captured.
  • the target signal can be captured accurately and reliably by repeating the determination of whether or not the target signal is captured while shifting the frequency at a constant frequency.
  • the target signal is a GPS signal broadcast for each GPS satellite.
  • This method specifically indicates that the target signal (search target signal) is a GPS signal. That is, it shows that the influence of cross-correlation in capturing the GPS signal is reduced, and the GPS signal is captured more reliably.
  • the target signal can be captured more accurately than the conventional method and more accurately than the conventional method.
  • FIG. It is a flowchart of the signal search method of embodiment of this invention. It is a figure for demonstrating the concept of the code phase of a correlation value calculation, and frequency setting. The frequency characteristic of the peak correlation value Peak (i) when the target GPS signal is captured and the frequency characteristic of the peak correlation value Peak (i) in the case of cross-correlation are shown. It is a figure which shows the correlation characteristic of the target GPS signal, and the correlation characteristic of cross correlation. It is a flowchart which shows the determination flow of the target GPS signal and cross correlation. It is a flowchart which shows the determination flow of the target GPS signal and cross correlation using an average value and a standard deviation.
  • FIG. 1 It is a flowchart which shows the determination flow of the target GPS signal at the time of switching and using two types of maximum noise values, and cross correlation. It is a block diagram which shows the structure of the GPS signal receiver 1 which concerns on embodiment of this invention. It is a block diagram which shows the structure of the information equipment terminal 100 provided with the GPS signal receiver 1. FIG.
  • a signal search method will be described with reference to the drawings.
  • a search method for an L1 wave hereinafter simply referred to as a GPS signal
  • a GPS signal a search method for an L1 wave (hereinafter simply referred to as a GPS signal) of a GPS signal transmitted from a GPS satellite
  • FIG. 2 is a flowchart of the signal search method of this embodiment.
  • FIG. 3 is a diagram for explaining the concept of code phase and frequency setting for correlation value calculation.
  • the correlation between the received signal and the replica signal with respect to the frequency section of 8 kHz with the frequency 1575.42 MHz of the GPS signal as the approximate center frequency of the scanning band.
  • a replica signal is a signal that reproduces a spreading code used in a GPS signal.
  • the integrated correlation value is calculated by integrating the correlation value over a predetermined time length for each combination of code phase and frequency.
  • frequencies set at 1000 Hz intervals are grouped in the frequency axis direction, and a correlation process is executed in parallel to calculate an integrated correlation value.
  • frequencies F sig 011, F sig 012, F sig 013, F sig 014, F sig 015, F sig 016, F sig 017, F sig 018 are shown.
  • the integrated correlation value is calculated simultaneously (S101).
  • the correlation processing is simultaneously performed at 8 points.
  • the correlation processing is simultaneously performed at eight points.
  • Such correlation processing is sequentially executed from the group Gr1 to the group Gr19, thereby obtaining an integrated correlation value for one channel covering the entire scanning band (8 kHz).
  • the correlation processing from the group Gr1 is sequentially repeated again.
  • the transition between groups Gr may perform the determination of the cross correlation mentioned later in parallel, below, when the cross correlation is determined in the group Gr during the cross correlation determination, in other words, The process of moving to the next group Gr when the target signal cannot be captured is shown.
  • the thus obtained peak correlation value Peak (i) for each frequency i has characteristics as shown in FIG. 4A shows the frequency characteristic of the peak correlation value Peak (i) when the target GPS signal is captured, and FIG. 4B shows the frequency characteristic of the peak correlation value Peak (i) in the case of cross-correlation. Indicates.
  • the peak correlation value of a specific frequency is significantly higher than the peak correlation values of other frequencies.
  • the peak correlation value Peak (2) of the frequency 2 of the frequency F sig 012 is the other frequencies F sig 011, F sig 013 to F sig 018 of the group Gr1.
  • the peak correlation values Peak (1) and Peak (3) to Peak (8) at frequencies 1 and 3 to 8 are significantly larger.
  • FIG. 5A shows the correlation characteristics of the target GPS signal
  • FIG. 5B shows the correlation characteristics of cross correlation.
  • the integrated correlation value becomes very high only at the frequency at the time of reception of the target GPS signal.
  • the integrated correlation value is significantly lowered at other frequencies.
  • the accumulated correlation value becomes maximum at intervals of 1000 Hz, and becomes substantially the same accumulated correlation value. This is because the C / A code is composed of 1023 chips and is superimposed on the GPS signal at a bit rate of 1.023 Mbps.
  • the target GPS signal is captured and cross-correlation is determined according to the following flow.
  • FIG. 6 is a flowchart showing a determination flow of the target GPS signal and cross correlation.
  • the peak correlation values Peak (i) of each frequency i are compared, and these maximum values are detected as the maximum peak value Tpeak (S301).
  • the peak correlation value Peak (i) of each frequency i is compared, and the peak correlation value Peak (i) having the second highest value among them is set as the maximum noise value Np (S302).
  • the peak correlation value Peak (i) that is the maximum value excluding the maximum peak value Tpeak is set as the maximum noise value Np.
  • a ratio Ratio between the maximum peak value Tpeak and the maximum noise value Np is calculated by an arithmetic expression of Tpeak / Np (S303).
  • the preset threshold Th is compared with the ratio Ratio, and if the ratio Ratio is higher than the threshold Th (S304: Yes), it is determined that the maximum peak value Tpeak is due to the target GPS signal (S306). . On the other hand, if the ratio Ratio is equal to or less than the threshold Th (S304: No), it is determined that these peak correlation value Peak (i) groups are due to cross-correlation (S305).
  • the peak correlation value Peak (2) of the frequency F sig 012 is set to the maximum peak value Tpeak.
  • the peak correlation value Peak (4) of the frequency F sig 014 is set to the maximum noise value Np.
  • the ratio Ratio (Tpeak / Np) is Peak (2) / Peak (4).
  • the peak correlation value Peak (2) is significantly higher than the peak correlation value Peak (4) due to the correlation characteristics of the GPS signal. Therefore, the ratio Ratio is a very high positive value.
  • the peak correlation value Peak (2) of the frequency F sig 012 is set to the maximum peak value Tpeak.
  • the peak correlation value Peak (4) of the frequency F sig 014 is set to the maximum noise value Np.
  • the ratio Ratio (Tpeak / Np) is Peak (2) / Peak (4).
  • the ratio Ratio is about 1.
  • the threshold value Th is set to a predetermined positive value of, for example, about 2 to 3, in the case of a GPS signal, the ratio Ratio is higher than the threshold value Th, and the maximum peak value Tpeak is determined by the target GPS signal. Can be accurately determined.
  • the ratio Ratio becomes lower than the threshold Th, and it can be accurately determined as cross-correlation.
  • the target GPS signal is a cross correlation or not, and it can be determined that the integrated correlation value is obtained from the target GPS signal (S104: No), the maximum peak value Tpeak is obtained.
  • the target GPS signal is captured (S105), and the process proceeds to the tracking process.
  • each frequency constituting the group is shifted by 50 Hz, that is, the group is shifted to start the above-described integration correlation value acquisition.
  • the relation determination process is performed in the same manner. This process is repeatedly executed continuously until, for example, a target GPS signal is captured. However, if the target GPS signal cannot be obtained even after scanning the entire scanning band a predetermined number of times (for example, a few times), for example, a display indicating that the target GPS signal could not be captured is performed. Etc. may be performed.
  • each process for the integrated correlation value is executed in parallel at eight independent frequencies, so that a higher speed process can be realized.
  • the scanning range on the code phase axis is determined by either the conventional technique or the method of the present embodiment.
  • the true second peak correlation value may not be within the scanning range, and the cross-correlation determination accuracy may decrease.
  • the maximum integrated correlation value is acquired on the code phase axis of each frequency, and the maximum peak value and the second peak are obtained from the maximum integrated correlation values at a plurality of independent frequencies.
  • the true second peak correlation value integrated correlation value
  • the cross-correlation can be accurately determined and the target GPS signal can be accurately captured, although the process is simpler than the conventional technique.
  • FIG. 7 is a flowchart showing a determination flow of a target GPS signal and cross correlation using an average value and a standard deviation.
  • the peak correlation values Peak (i) of each frequency i are compared, and these maximum values are detected as the maximum peak value Tpeak (S311).
  • the integrated correlation value is acquired with 0.5 chip resolution for 1023 chips by the above-described 8-frequency simultaneous processing
  • the standard deviation ⁇ is calculated by taking the square root of the subtraction result.
  • the maximum noise value Np ′ is calculated from the average value E and the standard deviation ⁇ (S314).
  • N is a constant set according to the scanning range along the code phase axis direction.
  • N is set smaller as the scanning range is expanded, and N is set larger as the scanning range is narrowed.
  • Np ′ By setting the constant N in this way according to the scanning range, the influence on the noise calculation error due to the number of samples (the size of the scanning range) can be suppressed, and the maximum noise value with high reliability even if the scanning range is changed Np ′ can be set.
  • a ratio Ratio ′ between the maximum peak value Tpeak and the maximum noise value Np ′ is calculated by an arithmetic expression of Tpeak / Np ′ (S315).
  • the preset threshold value Th ′ and the ratio Ratio ′ are compared, and if the ratio Ratio ′ is higher than the threshold value Th ′ (S316: Yes), it is determined that the maximum peak value Tpeak is due to the target GPS signal. (S318). On the other hand, if the ratio Ratio ′ is equal to or less than the threshold Th ′ (S316: No), it is determined that the peak correlation value Peak (i) group is due to cross-correlation (S317).
  • the method of setting the maximum noise value Np ′ using the average value E and the standard deviation ⁇ is used, the statistical reliability of the maximum noise value can be improved. Therefore, the cross correlation can be determined more accurately.
  • FIG. 8 is a flowchart showing a determination flow of a target GPS signal and cross correlation when two kinds of maximum noise values are switched and used.
  • the peak correlation values Peak (i) of each frequency i are compared, and these maximum values are detected as the maximum peak value Tpeak (S321).
  • the peak correlation values Peak (i) of each frequency i are compared, and 2 of these are compared.
  • the highest peak correlation value Peak (i) is set to the maximum noise value Np (S323).
  • a ratio Ratio between the maximum peak value Tpeak and the maximum noise value Np is calculated by an arithmetic expression of Tpeak / Np (S324).
  • the preset threshold Th is compared with the ratio Ratio, and if the ratio Ratio is higher than the threshold Th (S325: Yes), it is determined that the maximum peak value Tpeak is due to the target GPS signal (S326). . On the other hand, if the ratio Ratio is equal to or less than the threshold Th (S325: No), it is determined that these peak correlation value Peak (i) groups are due to cross-correlation (S327).
  • the maximum noise value Np ′ is calculated from the average value E and the standard deviation ⁇ (S329).
  • N is a constant set according to the scanning range along the code phase axis direction.
  • a ratio Ratio ′ between the maximum peak value Tpeak and the maximum noise value Np ′ is calculated by an arithmetic expression of Tpeak / Np ′ (S330).
  • the preset threshold value Th ′ is compared with the ratio Ratio ′, and if the ratio Ratio ′ is higher than the threshold value Th ′ (S331: Yes), it is determined that the maximum peak value Tpeak is due to the target GPS signal. (S326). On the other hand, if the ratio Ratio ′ is equal to or less than the threshold Th ′ (S331: No), it is determined that the peak correlation value Peak (i) group is due to cross-correlation (S327).
  • the second highest peak correlation value is used. Not easily affected by the volume. Further, since the number of parameters is large, even if the maximum noise value is set with the second highest peak correlation value, the value has sufficient reliability. On the other hand, when the number of parameters is small, the maximum noise value is set by a statistical method using the average value E and the standard deviation ⁇ , so that a decrease in reliability of the maximum noise value due to the small number of parameters can be suppressed. . Furthermore, since the number of parameters is small, the increase in the arithmetic processing load is not significantly affected. Accordingly, by using a combination of the two maximum noise value setting methods, the reliability of the maximum noise value can be maintained, and the cross correlation can be determined without increasing the processing load so much. it can.
  • FIG. 8 is a block diagram showing a configuration of the GPS signal receiving apparatus 1 according to the embodiment of the present invention.
  • the GPS signal receiving device 1 includes a GPS receiving antenna 10, an RF processing unit 20, a baseband processing unit 30, and a positioning calculation unit 40.
  • the GPS receiving antenna 10 receives a GPS signal broadcast (transmitted) from each GPS satellite and outputs it to the RF processing unit 20.
  • the RF processing unit 20 down-converts the received GPS signal, generates an intermediate frequency signal (IF signal), and outputs it to the baseband processing unit 30.
  • the baseband processing unit 30 corresponds to a “signal search device” including a “peak correlation value acquisition unit” and a “determination unit” of the present invention.
  • the baseband processing unit 30 also corresponds to the “capture tracking unit” of the present invention.
  • the baseband processing unit 30 may individually implement hardware corresponding to the “determination unit” and hardware corresponding to the “peak correlation value acquisition unit” and the “capture tracking unit”. It may be realized by hardware.
  • the baseband processing unit 30 generates a baseband signal by multiplying the IF signal by the carrier frequency signal, and performs a GPS signal capturing process and a tracking process using the baseband signal. At this time, the signal search method described above is used for the acquisition process. Thereby, the erroneous capture of cross-correlation can be suppressed and the target GPS signal can be reliably captured.
  • the capturing process for the GPS signal captured in this way shifts to a tracking process.
  • the code correlation result and the carrier correlation result obtained by this tracking, and the pseudo distance obtained from the code correlation result are output to the positioning calculation unit 40.
  • the positioning calculation unit 40 demodulates the navigation message based on the code correlation result, and performs positioning of the GPS signal receiving device 1 from the code correlation result, the carrier phase result, and the pseudorange.
  • the baseband processing unit 30 that executes the above-described signal search method may be realized by a hardware group that executes each process, and stores each process of the above-described signal search method in a storage medium in a programmed state. In addition, it may be realized by a mode in which the computer reads and executes the program.
  • FIG. 9 is a block diagram illustrating a main configuration of the information equipment terminal 100 including the GPS signal receiving device 1 of the present embodiment.
  • An information equipment terminal 100 as shown in FIG. 9 is, for example, a mobile phone, a car navigation device, a PND, a camera, a clock, and the like, and includes an antenna 10, an RF processing unit 20, a baseband processing unit 30, a positioning calculation unit 40, and application processing. Part 130 is provided.
  • the antenna 10, the RF processing unit 20, the baseband processing unit 30, and the positioning calculation unit 40 have the above-described configuration, and the GPS signal receiving device 1 is configured as described above.
  • the application processing unit 130 displays the own device position and the own device speed based on the positioning result output from the GPS signal receiving device 1, and executes processing for use in navigation and the like.
  • the present invention can be similarly applied to capturing other GNSS signals. Furthermore, the present invention can be similarly applied to capturing a wireless communication signal in which a peak appears in a correlation value at a predetermined frequency interval.
  • the GPS signal (L1 wave) is captured.
  • the above signal search method and configuration can be similarly applied to other GNSS signals.
  • the frequency interval may be set according to the bit rate of the spreading code used for the target GNSS signal.
  • GPS signal receiving device 10: GPS receiving antenna
  • 20 RF processing unit
  • 30 baseband processing unit
  • 40 positioning calculation unit
  • 100 information equipment terminal
  • 130 application processing unit

Abstract

【課題】従来方法よりも正確に目的のＧＰＳ信号を捕捉する。 【解決手段】所定周波数間隔で設定された各周波数ｉのピーク相関値Ｐｅａｋ（ｉ）を比較し、これらの最大値を最大ピーク値Ｔｐｅａｋとして検出する。各周波数ｉのピーク相関値Ｐｅａｋ（ｉ）を比較し、これらの中で２番目に高い値のピーク相関値Ｐｅａｋ（ｉ）を、最大ノイズ値Ｎｐに設定する。最大ピーク値Ｔｐｅａｋと最大ノイズ値Ｎｐとの比Ｒａｔｉｏを算出する。比Ｒａｔｉｏが閾値Ｔｈよりも高ければ（図（Ａ））、最大ピーク値Ｔｐｅａｋは目的のＧＰＳ信号によるものだと判定し、ＧＰＳ信号を捕捉する。一方、比Ｒａｔｉｏが閾値Ｔｈ以下であれば（図（Ｂ））、これらピーク相関値Ｐｅａｋ（ｉ）群は、クロスコリレーションによるものだと判定する。

Description

信号サーチ方法、信号サーチプログラム、信号サーチ装置、ＧＮＳＳ信号受信装置、および情報機器端末
　この発明は、受信信号から所望信号をサーチする信号サーチ方法、特に、ＧＮＳＳにおけるＧＰＳ信号をサーチする方法に関する。
　現在、ＧＮＳＳ（Ｇｌｏｂａｌ　Ｎａｖｉｇａｔｉｏｎ　Ｓａｔｅｌｌｉｔｅ　Ｓｙｓｔｅｍ）の一つとして、ＧＰＳ（Ｇｌｏｂａｌ　Ｐｏｓｉｔｉｏｎｉｎｇ　Ｓｙｓｔｅｍ）がある。
　ＧＰＳでは、複数のＧＰＳ衛星から送信されるＧＰＳ信号を受信し、受信したＧＰＳ信号のコード位相やキャリア位相を用いて、受信機の測位を行う。ＧＰＳでは、ＧＰＳ衛星毎に異なる拡散コードが設定されており、各ＧＰＳ信号は、異なる拡散コードでコード変調されている。
　このようなＧＰＳでは、目的とするＧＰＳ衛星とは異なるＧＰＳ衛星からのＧＰＳ信号を、目的とするＧＰＳ信号と誤って捕捉し、追尾処理を行ってしまうことがある。このような現象はクロスコリレーションと呼ばれている。
　クロスコリレーションを防ぐ方法として、特許文献１では、単一のサーチ周波数におけるコード位相上の各コード位相点の積算相関値を比較することで、取得した積算相関値が、目的とするＧＰＳ信号によるものか、クロスコリレーションによるものかを判定している。図１は、特許文献１に記載された信号サーチ方法の概念を示す図である。
　特許文献１の方法では、単一のサーチ周波数の各コード位相での積算相関値を取得し、これらの値を比較する。特許文献１の方法では、全てのコード位相での積算相関値における最大値となる積算相関値Ｖｐｋ１を検出する。次に、特許文献１の方法では、全てのコード位相での積算相関値における２番目に値の高い積算相関値Ｖｐｋ２を検出する。特許文献１の方法では、最大値となる積算相関値Ｖｐｋ１と２番目に値の高い積算相関値Ｖｐｋ２との比を算出する。特許文献１の方法では、当該比が所定閾値以上であれば、最大値となる積算相関値Ｖｐｋ１に対応するコード位相Ｃｐｋ１が、目的信号によるものとして、当該コード位相Ｃｐｋ１の信号を捕捉する。一方、比が閾値以下であれば、信号サーチを継続する。
米国特許７１６１９７７号明細書
　しかしながら、特許文献１の方法では、最大となる積算相関値Ｖｐｋ１を検出後に、再度コード位相軸上を走査して、２番目に値の高い積算相関値Ｖｐｋ２を検出しなければならない。したがって、単に、最大となる積算相関値Ｖｐｋ１のコード位相Ｃｐｋ１を検出する場合と比較して、２倍のサーチ時間を必要としてしまう。
　また、特許文献１の方法では、２番目に高い積算相関値は、最大となる積算相関値Ｖｐｋ１を含む極大に向かう途中のコード位相（図１の例であれば、コード位相Ｃｐｋ２ｗ）となることがある。したがって、積算相関値Ｖｐｋ１のコード位相Ｃｐｋ１を含む所定コード位相範囲を含まないように、２番目に高い積算相関値の検出を行わなければならず、処理が複雑化してしまう。
　さらに、コード位相軸上での２番目に高い積算相関値の検出範囲を、コード位相全域でなく、図１に示すようにコード位相範囲Ｃｏｖに制限すると、真の２番目に高い積算相関値Ｖｐｋ２よりも値の低い積算相関値Ｖｐｋ２'を検出することがある。この積算相関値Ｖｐｋ２'を用いた場合、最大値となる積算相関値Ｖｐｋ１との比は、真の比よりも大きくなってしまう。このため、クロスコリレーションであるにも係わらず、目的のＧＰＳ信号と判定してしまう誤判定が生じやすい。
　したがって、本発明の目的は、従来方法よりも比較的簡素な処理で、従来よりも正確に目的信号を捕捉できる信号サーチ方法を提供することにある。
　この発明の信号サーチ方法では、次の各工程を有する。相関値取得工程では、受信信号と複数の周波数で生成された前記目的の信号のレプリカ信号との相関処理を行って相関値を取得する。ピーク相関値取得工程では、各周波数での最大相関値であるピーク相関値を取得する。判定工程では、各周波数のピーク相関値の大小関係に基づいて、目的の信号を捕捉したかどうかを判定する。
　この方法では、次の２つの相関値の周波数特性に基づいて、目的の信号を捕捉したかどうか、例えば、クロスコリレーションであるか目的の信号であるかを判定している。
　（ｉ）クロスコリレーションでは、所定周波数間隔（具体的には、拡散コードの周波数に応じた間隔）で相関値のピークが生じる。
　（ｉｉ）真の信号（目的の信号）では、ピーク相関値のみが突出して高いことを利用している。
　このような相関値の周波数特性を利用すれば、同時に取得する複数の周波数での最大相関値からなるピーク相関値を比較するだけで、クロスコリレーションであるか目的の信号であるかの判定を行うことが可能になる。これにより、比較的簡素な処理で且つ正確にクロスコリレーション判定、すなわち目的の信号の検出を行うことができる。
　また、この発明の信号サーチ方法の判定工程では、各周波数のピーク相関値の最大値となる最大ピーク値を検出する工程と、最大ピーク値を除く各周波数のピーク相関値から最大ノイズ値を決定する工程と、最大ピーク値と最大ノイズ値との比を算出する工程と、を有する。そして、判定工程では、算出した比に基づいて目的の信号を捕捉したかどうかを判定する。
　この方法では、目的の信号を捕捉したかどうかの判定方法の具体的な内容を示している。クロスコリレーションの相関値と目的の信号の相関値とは、上述のような特性を有するので、最大ピーク値を取得して、当該最大ピーク値を、同時に取得した他の周波数のピーク相関値から得られる最大ノイズ値と比較すれば、目的の信号であるかクロスコリレーションであるかを判定できる。
　また、この発明の信号サーチ方法の最大ノイズ値を決定する工程では、最大ピーク値を除く各周波数のピーク相関値における最大のピーク相関値を用いる。
　この方法では、最大ノイズ値の具体的な決定方法を示している。この方法では、最大ピーク値に続く２番目に相関値が高いピーク相関値を最大ノイズ値に設定している。これにより、より簡素な処理で最大ノイズ値を設定できる。
　また、この発明の信号サーチ方法の最大ノイズ値を決定する工程では、複数の周波数の相関値の平均値および標準偏差を算出する工程と、コード位相のサーチ範囲に基づいて設定された定数を決定する工程と、を有する。そして、最大ノイズ値を決定する工程では、標準偏差に定数を乗算し、この乗算結果に平均値を加算することで、最大ノイズ値を決定する。
　この方法では、平均値と標準偏差を用いることで、統計的に、より信頼性の高い最大ノイズ値を設定できる。
　また、この発明の信号サーチ方法の最大ノイズ値を決定する工程は、最大ピーク値を除く各周波数のピーク相関値における最大のピーク相関値を最大ノイズ値に設定する第１の決定工程を有する。さらに、最大ノイズ値を決定する工程は、複数の周波数の相関値の平均値および標準偏差を算出する工程と、コード位相のサーチ範囲に基づいて設定された定数を決定する工程と、を有し、標準偏差に定数を乗算し、この乗算結果に平均値を加算することで最大ノイズ値を決定する第２の決定工程と、を有する。そして、最大ノイズ値を決定する工程では、相関値を算出するコード位相範囲が広い場合に第１の決定工程で最大ノイズ値を決定し、相関値を算出するコード位相範囲が狭い場合に第２の決定工程で最大ノイズ値を決定する。
　この構成では、上述の最大ピーク値に続く２番目に相関値が高いピーク相関値で最大ノイズ値を設定する方法と、平均値と標準偏差から最大ノイズ値を設定する方法を、コード位相の走査範囲に応じて、適切に切り替えて用いている。
　また、この発明の信号サーチ方法では、判定工程によって目的の信号を捕捉できていないと判定された時、複数の周波数を、同時に相関値を取得する前記周波数間隔よりも短い一定の周波数幅で周波数シフトさせて、上述の目的の信号を捕捉したかどうかの判定を繰り返す。
　この方法に示すように、周波数を一定の周波数でシフトさせながら、目的の信号を捕捉したかどうかの判定を繰り返すことで、目的の信号を正確且つ確実に捕捉することができる。
　また、この発明の信号サーチ方法では、目的の信号は、ＧＰＳ衛星毎に放送されるＧＰＳ信号である。
　この方法では、具体的に目的の信号（サーチ対象の信号）がＧＰＳ信号であることを示している。すなわち、ＧＰＳ信号の捕捉におけるクロスコリレーションによる影響を低減させ、ＧＰＳ信号を、より確実に捕捉することを示している。
　この発明によれば、従来方法よりも比較的簡素な方法で、且つ従来用法よりも目的信号を正確に捕捉できる。
特許文献１に記載された信号サーチ方法の概念を示す図である。 本発明の実施形態の信号サーチ方法のフローチャートである。 相関値算出のコード位相、周波数設定の概念を説明するための図である。 目的のＧＰＳ信号を捕捉した場合のピーク相関値Ｐｅａｋ（ｉ）の周波数特性、および、クロスコリレーションの場合のピーク相関値Ｐｅａｋ（ｉ）の周波数特性を示す。 目的のＧＰＳ信号の相関特性、およびクロスコリレーションの相関特性を示す図である。 目的のＧＰＳ信号およびクロスコリレーションの判定フローを示すフローチャートである。 平均値および標準偏差を用いた目的のＧＰＳ信号およびクロスコリレーションの判定フローを示すフローチャートである。 ２種の最大ノイズ値を切り替えて用いた場合の目的のＧＰＳ信号およびクロスコリレーションの判定フローを示すフローチャートである。 本発明の実施形態に係るＧＰＳ信号受信装置１の構成を示すブロック図である。 ＧＰＳ信号受信装置１を備える情報機器端末１００の構成を示すブロック図である。
　本発明の実施形態に係る信号サーチ方法について、図を参照して説明する。なお、本実施形態では、ＧＰＳ衛星から送信されるＧＰＳ信号のＬ１波（以下、単にＧＰＳ信号と称する。）をサーチ方法について説明する。
　図２は本実施形態の信号サーチ方法のフローチャートである。図３は相関値算出のコード位相、周波数設定の概念を説明するための図である。
　図２に示すように、本実施形態の信号サーチ方法では、まず、ＧＰＳ信号の周波数１５７５．４２ＭＨｚを走査帯域の略中心の周波数として８ｋＨｚの周波数区間に対して、受信信号とレプリカ信号との相関処理を行う。レプリカ信号とは、ＧＰＳ信号に用いられている拡散コードを再現した信号である。そして、本実施形態の信号サーチ方法では、各コード位相および周波数の組合せ毎に、所定時間長に亘り相関値を積算することで、積算相関値を算出する。
　さらに、本実施形態の信号サーチ方法では、図３に示すように、コード位相軸方向に対しては、１個の周波数における１０２３ｃｈｉｐのコード位相に対して、０．５ｃｈｉｐの分解能で、全てのコード位相点を同時に相関処理し、コード位相点毎に積算相関値を算出する。
　また、本実施形態の信号サーチ方法では、周波数軸方向に対しては、１０００Ｈｚ間隔で設定された８個の周波数を１グループとし、同時並行して相関処理を実行し積算相関値を算出する。具体的な例としては、図３に示すように、１０００Ｈｚ間隔からなる周波数Ｆｓｉｇ０１１，Ｆｓｉｇ０１２，Ｆｓｉｇ０１３，Ｆｓｉｇ０１４，Ｆｓｉｇ０１５，Ｆｓｉｇ０１６，Ｆｓｉｇ０１７，Ｆｓｉｇ０１８のグループＧｒ１において、同時に積算相関値が算出される（Ｓ１０１）。
　次に、グループＧｒ１に対して、周波数を５０ＨｚシフトさせたグループＧｒ２（周波数Ｆｓｉｇ０２１を含むグループ）において、８点で同時に相関処理が実行される。
　次に、グループＧｒ２に対して、周波数を５０ＨｚシフトさせたグループＧｒ３（周波数Ｆｓｉｇ０３１を含むグループ）において、８点で同時に相関処理が実行される。
　このような相関処理を、グループＧｒ１からグループＧｒ１９まで順次実行していくことで、走査帯域の全域（８ｋＨｚ）を網羅する１チャンネル分の積算相関値を得る。そして、１チャンネル分の相関処理が終了すると、再度グループＧｒ１からの相関処理が、順次繰り返して実行される。なお、グループＧｒ間の移行は、後述するクロスコリレーションの判定を同時並行して行ってもよいが、以下では、クロスコリレーション判定中のグループＧｒにおいて、クロスコリレーション判定された場合、言い換えれば目的の信号を捕捉できなかった場合に、次のグループＧｒに移行する処理を示す。
　次に、１つのグループＧｒを構成する周波数ｉ（ｉ＝１～８）毎に、当該周波数ｉ上に存在する全ての積算相関値を比較する。そして、最大値となる積算相関値を、当該周波数ｉのピーク相関値Ｐｅａｋ（ｉ）として取得する（Ｓ１０２）。
　このように取得した、周波数ｉ毎のピーク相関値Ｐｅａｋ（ｉ）は、図４に示すような特性となる。図４（Ａ）は目的のＧＰＳ信号を捕捉した場合のピーク相関値Ｐｅａｋ（ｉ）の周波数特性を示し、図４（Ｂ）はクロスコリレーションの場合のピーク相関値Ｐｅａｋ（ｉ）の周波数特性を示す。
　目的のＧＰＳ信号を捕捉した場合、特定の周波数のピーク相関値が、他の周波数のピーク相関値よりも大幅に高くなる。具体的に、図４（Ａ）の例であれば、周波数Ｆｓｉｇ０１２の周波数２のピーク相関値Ｐｅａｋ（２）が、グループＧｒ１の他の周波数Ｆｓｉｇ０１１，Ｆｓｉｇ０１３～Ｆｓｉｇ０１８の周波数１，３～８のピーク相関値Ｐｅａｋ（１），Ｐｅａｋ（３）～Ｐｅａｋ（８）よりも大幅に大きくなる。
　クロスコリレーションの場合、取得した全ての周波数１～８のピーク相関値Ｐｅａｋ（１）～Ｐｅａｋ（８）が略同じなる。
　これは、目的のＧＰＳ信号に対する相関特性と、クロスコリレーションの相関特性とが図５に示すような周波数特性を有するからである。図５（Ａ）は目的のＧＰＳ信号の相関特性を示し、図５（Ｂ）はクロスコリレーションの相関特性を示す図である。
　図５（Ａ）に示すように、目的のＧＰＳ信号の場合、積算相関値は、目的のＧＰＳ信号の受信時の周波数のみで積算相関値が非常に高くなる。その一方で、他の周波数では積算相関値が大幅に低くなる。
　図５（Ｂ）に示すように、クロスコリレーションの場合、積算相関値が１０００Ｈｚ間隔で極大となり、略同じ積算相関値となる。これは、Ｃ／Ａコードが、１０２３ｃｈｉｐからなり、ビットレート１．０２３ＭｂｐｓでＧＰＳ信号に重畳されていることに起因している。
　このような相関特性を利用し、本実施形態の信号サーチ方法では、引き続き、次に示すフローで、目的のＧＰＳ信号の捕捉およびクロスコリレーションの判定を行う。
　上述のように、取得した各周波数ｉのピーク相関値Ｐｅａｋ（ｉ）を比較してクロスコリレーションの判定を行う（Ｓ１０３）。図６は目的のＧＰＳ信号およびクロスコリレーションの判定フローを示すフローチャートである。
　まず、各周波数ｉのピーク相関値Ｐｅａｋ（ｉ）を比較し、これらの最大値を最大ピーク値Ｔｐｅａｋとして検出する（Ｓ３０１）。
　次に、各周波数ｉのピーク相関値Ｐｅａｋ（ｉ）を比較し、これらの中で２番目に高い値のピーク相関値Ｐｅａｋ（ｉ）を、最大ノイズ値Ｎｐに設定する（Ｓ３０２）。言い換えれば、各周波数ｉのピーク相関値Ｐｅａｋ（ｉ）の中で、最大ピーク値Ｔｐｅａｋを除いて最大値となるピーク相関値Ｐｅａｋ（ｉ）を、最大ノイズ値Ｎｐに設定する。
　次に、最大ピーク値Ｔｐｅａｋと最大ノイズ値Ｎｐとの比Ｒａｔｉｏを、Ｔｐｅａｋ／Ｎｐの演算式によって算出する（Ｓ３０３）。
　次に、予め設定した閾値Ｔｈと比Ｒａｔｉｏとを比較して、比Ｒａｔｉｏが閾値Ｔｈよりも高ければ（Ｓ３０４：Ｙｅｓ）、最大ピーク値Ｔｐｅａｋは目的のＧＰＳ信号によるものだと判定する（Ｓ３０６）。一方、比Ｒａｔｉｏが閾値Ｔｈ以下であれば（Ｓ３０４：Ｎｏ）、これらピーク相関値Ｐｅａｋ（ｉ）群は、クロスコリレーションによるものだと判定する（Ｓ３０５）。
　具体的に、図４（Ａ）の例では、周波数Ｆｓｉｇ０１２のピーク相関値Ｐｅａｋ（２）が最大ピーク値Ｔｐｅａｋに設定される。周波数Ｆｓｉｇ０１４のピーク相関値Ｐｅａｋ（４）が最大ノイズ値Ｎｐに設定される。比Ｒａｔｉｏ（Ｔｐｅａｋ／Ｎｐ）は、Ｐｅａｋ（２）／Ｐｅａｋ（４）となる。ここで、上述のように、ＧＰＳ信号の相関特性から、ピーク相関値Ｐｅａｋ（２）は、ピーク相関値Ｐｅａｋ（４）よりも大幅に高い。したがって、比Ｒａｔｉｏは、非常に高い正値となる。
　図４（Ｂ）の例でも、周波数Ｆｓｉｇ０１２のピーク相関値Ｐｅａｋ（２）が最大ピーク値Ｔｐｅａｋに設定される。周波数Ｆｓｉｇ０１４のピーク相関値Ｐｅａｋ（４）が最大ノイズ値Ｎｐに設定される。比Ｒａｔｉｏ（Ｔｐｅａｋ／Ｎｐ）は、Ｐｅａｋ（２）／Ｐｅａｋ（４）となる。しかしながら、上述のように、クロスコリレーションの場合、ピーク相関値Ｐｅａｋ（２）とピーク相関値Ｐｅａｋ（４）を含む他のピーク相関値Ｐｅａｋとは、略同じ値となる。したがって、比Ｒａｔｉｏは約１程度となる。
　このため、閾値Ｔｈを、例えば２～３程度の所定の正値に設定すれば、ＧＰＳ信号の場合には、比Ｒａｔｉｏが閾値Ｔｈよりも高くなり、最大ピーク値Ｔｐｅａｋが目的のＧＰＳ信号によるものと正確に判定できる。一方、クロスコリレーションの場合には、比Ｒａｔｉｏが閾値Ｔｈよりも低くなり、クロスコリレーションと正確に判定できる。
　このように、目的のＧＰＳ信号であるかクロスコリレーションであるかが判定され、目的のＧＰＳ信号から積算相関値を得ていることが判定できると（Ｓ１０４：Ｎｏ）、最大ピーク値Ｔｐｅａｋとなる積算相関値を得た周波数とコード位相を取得することで、目的のＧＰＳ信号の捕捉し（Ｓ１０５）、追尾処理に移行する。
　ここで、クロスコリレーションと判定されれば（Ｓ１０４：Ｙｅｓ）、グループを構成する各周波数をそれぞれ５０Ｈｚずつシフトして、すなわち、グループを移行して、上述の積算相関値の取得に始まるクロスコリレーションの判定処理を同様に行う。この処理は、例えば、目的のＧＰＳ信号を捕捉するまで継続的に繰り返し実行される。ただし、走査帯域の全域を所定回数（例えば２、３回）走査しても、目的のＧＰＳ信号を得られない場合には、例えば、目的のＧＰＳ信号を捕捉できなかったことを示す表示を行う等の処理を行ってもよい。
　以上のように、本実施形態の信号サーチ方法を用いれば、クロスコリレーションを正確に判定し、目的のＧＰＳ信号を正確且つ確実に捕捉することができる。
　さらに、本実施形態の方法を用いることで、離散して設定された各周波数に対して、最大の積算相関値（周波数別最大相関値）を算出する工程と、複数の周波数別最大相関値から、さらにクロスコリレーション判定用の最大積算相関値と、２番目に高い積算相関値とを検出するだけで、クロスコリレーションの判定を行うことができる。すなわち、従来技術に示すように、コード位相軸の全体に対して２度の最大値検出を行う必要がなく、処理が簡素化される。例えば、１０２３ｃｈｉｐのコード位相を０．５ｃｈｉｐの分解能で走査する場合、従来例であれば、１０２３×２×２＝４０９２回の積算相関値の取得を行わなければならず、これを８個の異なる周波数で取得する場合、４０９２×８＝３２７３６回の積算相関値の取得が必要になる。一方で、本実施形態の方法では、一周波数に対して１０２３×２＝２０４６回であり、８周波数では、２０４６×８＝１６３６８回となる。そして、８周波数の積算相関値から最大ピーク相関値を決定するために８回、２番目のピーク相関値を決定するのに７回の走査となる。したがって、１６３６８＋８＋７＝１６３８３回となり、従来の３２７３６回の半分程度で済む。その上、本実施形態では、独立する８周波数で平行して、積算相関値に対する各処理を実行するので、さらに高速な処理が実現できる。
　また、従来技術では、コード位相軸に存在する２番目の相関ピークを見つけるために、コード位相軸上での範囲制限を行わなければならなかったが、本実施形態では、互いに独立する複数の周波数のピーク相関値間から２番目のピーク相関値を決定するため、制限を設ける必要が無く、従来技術よりも簡素な処理が可能となる。
　また、既に航法メッセージを取得済みで、当該航法メッセージから目的のＧＰＳ信号の情報が得られた場合、従来技術であっても本実施形態の方法であってもコード位相軸上での走査範囲を狭くすることができるが、従来技術の場合、上述のように、真の２番目のピーク相関値が走査範囲内にないことが考えられ、クロスコリレーションの判定確度が低下することがある。しかしながら、本実施形態の方法では、各周波数のコード位相軸上では最大の積算相関値を取得し、これら互いに独立する複数の周波数での最大の積算相関値から、最大ピーク値と２番目のピーク相関値（最大ノイズ値）を得るため、クロスコリレーション判定を行うために必要とする真の２番目のピーク相関値（積算相関値）を得ることができ、クロスコリレーション判定確度を低下させることがない。
　以上のように、本実施形態の方法を用いれば、従来技術よりも簡素な処理でありながら、クロスコリレーションを正確に判定し、目的のＧＰＳ信号を正確に捕捉することができる。
　なお、上述の説明では、最大ノイズ値Ｎｐを２番目のピーク相関値とする例を示したが、次に示すように、各周波数ｉのピーク相関値Ｐｅａｋ（ｉ）の平均値および標準偏差から最大ノイズ値Ｎｐ'を設定してもよい。図７は、平均値および標準偏差を用いた目的のＧＰＳ信号およびクロスコリレーションの判定フローを示すフローチャートである。
　まず、各周波数ｉのピーク相関値Ｐｅａｋ（ｉ）を比較し、これらの最大値を最大ピーク値Ｔｐｅａｋとして検出する（Ｓ３１１）。
　次に、同時に取得した全ての周波数ｉ（ｉ＝１～８）の積算相関値の平均値Ｅを算出する（Ｓ３１２）。例えば、上述の８周波数同時処理で、１０２３ｃｈｉｐに対する０．５ｃｈｉｐ分解能で積算相関値を取得した場合、１０２３×２×８＝１６３６８個の積算相関値を加算し、この加算値を母数１６３６８で除算することで、平均値Ｅを算出する。
　次に、同時に取得した全ての周波数ｉ（ｉ＝１～８）の積算相関値の標準偏差σを算出する（Ｓ３１３）。例えば、上述の８周波数同時処理で、１０２３ｃｈｉｐに対する０．５ｃｈｉｐ分解能で積算相関値を取得した場合、１０２３×２×８＝１６３６８個の積算相関値の自乗和を算出し、当該自乗和を母数１６３６８で除算し、上述の平均値Ｅの自乗を減算する。そして、この減算結果の平方根と取ることで、標準偏差σを算出する。
　次に、平均値Ｅと標準偏差σとから最大ノイズ値Ｎｐ'を算出する（Ｓ３１４）。最大ノイズ値Ｎｐ'は、標準偏差σに定数Ｎを乗算し、乗算結果と平均値Ｅの加算により算出される。すなわち、Ｎｐ'＝Ｅ＋Ｎσから得られる。ここで、Ｎは、コード位相軸方向に沿った走査範囲に準じて設定される定数である。
　例えば、走査範囲が拡がるほどＮが小さくなり、走査範囲が狭くなるほどＮが大きく設定される。具体例としては、コード位相の全コード範囲（１０２３ｃｈｉｐ）を走査範囲とする場合にＮ＝４が設定され、コード位相の前コード範囲における１／４（２５６ｃｈｉｐ）を走査範囲とする場合にＮ＝５が設定される。
　走査範囲に応じて、このように定数Ｎを設定することで、サンプル数（走査範囲の大きさ）によるノイズ算定誤差に対する影響を抑圧でき、走査範囲を変化させても信頼性の高い最大ノイズ値Ｎｐ'を設定できる。
　次に、最大ピーク値Ｔｐｅａｋと最大ノイズ値Ｎｐ'との比Ｒａｔｉｏ'を、Ｔｐｅａｋ／Ｎｐ'の演算式によって算出する（Ｓ３１５）。
　次に、予め設定した閾値Ｔｈ'と比Ｒａｔｉｏ'とを比較して、比Ｒａｔｉｏ'が閾値Ｔｈ'よりも高ければ（Ｓ３１６：Ｙｅｓ）、最大ピーク値Ｔｐｅａｋは目的のＧＰＳ信号によるものだと判定する（Ｓ３１８）。一方、比Ｒａｔｉｏ'が閾値Ｔｈ'以下であれば（Ｓ３１６：Ｎｏ）、これらピーク相関値Ｐｅａｋ（ｉ）群は、クロスコリレーションによるものだと判定する（Ｓ３１７）。
　このように、平均値Ｅと標準偏差σを用いた最大ノイズ値Ｎｐ'を設定する方法を用いれば、最大ノイズ値の統計学上の信頼性を向上させることができる。したがって、より正確にクロスコリレーションを判定することができる。
　なお、上述の二種の最大ノイズ値の設定方法を、組み合わせて利用することもできる。図８は２種の最大ノイズ値を切り替えて用いた場合の目的のＧＰＳ信号およびクロスコリレーションの判定フローを示すフローチャートである。
　まず、各周波数ｉのピーク相関値Ｐｅａｋ（ｉ）を比較し、これらの最大値を最大ピーク値Ｔｐｅａｋとして検出する（Ｓ３２１）。
　コード位相の走査範囲を取得し、コード位相の走査範囲が広い（１０２３ｃｈｉｐのフルコード等）場合（Ｓ３２２：Ｙｅｓ）、各周波数ｉのピーク相関値Ｐｅａｋ（ｉ）を比較し、これらの中で２番目に高い値のピーク相関値Ｐｅａｋ（ｉ）を、最大ノイズ値Ｎｐに設定する（Ｓ３２３）。
　次に、最大ピーク値Ｔｐｅａｋと最大ノイズ値Ｎｐとの比Ｒａｔｉｏを、Ｔｐｅａｋ／Ｎｐの演算式によって算出する（Ｓ３２４）。
　次に、予め設定した閾値Ｔｈと比Ｒａｔｉｏとを比較して、比Ｒａｔｉｏが閾値Ｔｈよりも高ければ（Ｓ３２５：Ｙｅｓ）、最大ピーク値Ｔｐｅａｋは目的のＧＰＳ信号によるものだと判定する（Ｓ３２６）。一方、比Ｒａｔｉｏが閾値Ｔｈ以下であれば（Ｓ３２５：Ｎｏ）、これらピーク相関値Ｐｅａｋ（ｉ）群は、クロスコリレーションによるものだと判定する（Ｓ３２７）。
　コード位相の走査範囲を取得し、コード位相の走査範囲が狭い（１０２３ｃｈｉｐのフルコードの１／４（２５６ｃｈｉｐ）等）場合（Ｓ３２２：Ｎｏ）、同時に取得した全ての周波数ｉ（ｉ＝１～８）の積算相関値の平均値Ｅと標準偏差σを算出する（Ｓ３２８）。これら積算相関値の平均値Ｅおよび標準偏差σの算出方法は、上述の方法と同じ方法を用いる。
　次に、平均値Ｅと標準偏差σとから最大ノイズ値Ｎｐ'を算出する（Ｓ３２９）。最大ノイズ値Ｎｐ'は、標準偏差σに定数Ｎを乗算し、乗算結果と平均値Ｅの加算により算出される。すなわち、Ｎｐ'＝Ｅ＋Ｎσから得られる。ここで、Ｎは、コード位相軸方向に沿った走査範囲に準じて設定される定数である。
　次に、最大ピーク値Ｔｐｅａｋと最大ノイズ値Ｎｐ'との比Ｒａｔｉｏ'を、Ｔｐｅａｋ／Ｎｐ'の演算式によって算出する（Ｓ３３０）。
　次に、予め設定した閾値Ｔｈ'と比Ｒａｔｉｏ'とを比較して、比Ｒａｔｉｏ'が閾値Ｔｈ'よりも高ければ（Ｓ３３１：Ｙｅｓ）、最大ピーク値Ｔｐｅａｋは目的のＧＰＳ信号によるものだと判定する（Ｓ３２６）。一方、比Ｒａｔｉｏ'が閾値Ｔｈ'以下であれば（Ｓ３３１：Ｎｏ）、これらピーク相関値Ｐｅａｋ（ｉ）群は、クロスコリレーションによるものだと判定する（Ｓ３２７）。
　このような２つの最大ノイズ値の設定方法の組合せを用いる場合で、母数が多い場合に、２番目に高い値のピーク相関値を用いるため、最大ノイズ値の設定に対して、母数の多さの影響を受けにくい。さらに、母数が多いため、最大ノイズ値を２番目に高い値のピーク相関値で設定しても、その値は十分な信頼性を有する。一方、母数が少ない場合には、平均値Ｅと標準偏差σとを用いる統計学的手法により、最大ノイズ値を設定するので、母数の少なさによる最大ノイズ値の信頼性低下を抑制できる。さらに、母数が少ないために演算処理負荷の増加に余り影響を及ぼさない。したがって、２つの最大ノイズ値の設定方法の組合せを用いる方法とすることで、最大ノイズ値の信頼性を維持し、且つ演算処理負荷をあまり増加させることなく、クロスコリレーションの判定を行うことができる。
　以上のような処理は、次に示す構成のＧＰＳ信号受信装置で実現できる。図８は、本発明の実施形態に係るＧＰＳ信号受信装置１の構成を示すブロック図である。
　ＧＰＳ信号受信装置１は、ＧＰＳ受信アンテナ１０、ＲＦ処理部２０、ベースバンド処理部３０、および測位演算部４０を備える。
　ＧＰＳ受信アンテナ１０は、各ＧＰＳ衛星から放送（送信）されるＧＰＳ信号を受信し、ＲＦ処理部２０へ出力する。ＲＦ処理部２０は、受信したＧＰＳ信号をダウンコンバートして、中間周波数信号（ＩＦ信号）生成し、ベースバンド処理部３０へ出力する。
　ベースバンド処理部３０は、本発明の「ピーク相関値取得部」および「判定部」を備える「信号サーチ装置」に対応する。また、ベースバンド処理部３０は、本発明の「捕捉追尾部」にも対応する。なお、ベースバンド処理部３０は、「判定部」に相当するハードウェアと、「ピーク相関値取得部」および「捕捉追尾部」に相当するハードウェアとを個別に実現してもよく、一体のハードウェアで実現してもよい。ベースバンド処理部３０は、ＩＦ信号にキャリア周波数信号を乗算することでベースバンド信号を生成し、当該ベースバンド信号によるＧＰＳ信号の捕捉処理、および追尾処理を行う。この際、捕捉処理に、上述の信号サーチ方法を用いる。これにより、クロスコリレーションの誤捕捉を抑制し、目的とするＧＰＳ信号を確実に捕捉できる。
　このように捕捉されたＧＰＳ信号に対する捕捉処理は、追尾処理に移行する。この追尾で得られるコード相関結果やキャリア相関結果、さらにはコード相関結果から得られる擬似距離は、測位演算部４０へ出力される。
　測位演算部４０は、コード相関結果に基づいて航法メッセージを復調するとともに、コード相関結果、キャリア位相結果、擬似距離から、ＧＰＳ信号受信装置１の測位を行う。
　このような構成を用い、上述の信号サーチ方法を用いることで、誤捕捉が抑制され且つ高速な捕捉を行えるので、ＧＰＳ信号の追尾精度が向上し、結果的に測位結果の精度を向上させることもできる。
　なお、上述の信号サーチ方法を実行するベースバンド処理部３０は、各処理を実行するハードウェア群で実現してもよく、上述の信号サーチ方法の各処理をプログラム化した状態で記憶媒体に記憶しておき、コンピュータで当該プログラムを読み出して実行する態様によって実現してもよい。
　また、このようなＧＰＳ信号受信装置１やＧＰＳ信号受信機能は、図９に示すような情報機器端末１００に利用される。図９は、本実施形態のＧＰＳ信号受信装置１を備えた情報機器端末１００の主要構成を示すブロック図である。
　図９に示すような情報機器端末１００は、例えば携帯電話機、カーナビゲーション装置、ＰＮＤ、カメラ、時計等であり、アンテナ１０、ＲＦ処理部２０、ベースバンド処理部３０、測位演算部４０、アプリケーション処理部１３０を備える。アンテナ１０、ＲＦ処理部２０、ベースバンド処理部３０、測位演算部４０は、上述の構成のものであり、これらにより上述のようにＧＰＳ信号受信装置１が構成されている。
　アプリケーション処理部１３０は、ＧＰＳ信号受信装置１から出力された測位結果に基づいて、自装置位置や自装置速度を表示したり、ナビゲーション等に利用するための処理を実行する。
　このような構成において、上述のように高精度な測位結果を得られることで、高精度な位置表示やナビゲーション等を実現することができる。
　なお、上述の説明では、ＧＰＳ信号のクロスコリレーションの場合を例に説明したが、他のＧＮＳＳ信号の捕捉にも同様に適用できる。さらには、所定周波数間隔で相関値にピークが現れるような無線通信信号の捕捉に対しても、同様に適用できる。
　なお、上述の説明では、ＧＰＳ信号（Ｌ１波）の捕捉を行う場合を示したが、他のＧＮＳＳ信号に対しても、同様に上述の信号サーチ方法および構成を適用することができる。この場合、目的とするＧＮＳＳ信号に用いられる拡散コードのビットレートに準じて周波数間隔を設定すればよい。
１：ＧＰＳ信号受信装置、１０：ＧＰＳ受信アンテナ、２０：ＲＦ処理部、３０：ベースバンド処理部、４０：測位演算部、１００：情報機器端末、１３０：アプリケーション処理部

Claims (13)

  1. 　目的の信号を捕捉するための信号サーチ方法であって、
    　受信信号と複数の周波数で生成された前記目的の信号のレプリカ信号との相関処理を行って相関値を取得する相関値取得工程と、
    　各周波数での最大相関値であるピーク相関値を取得するピーク相関値取得工程と、
    　前記各周波数のピーク相関値の大小関係に基づいて、前記目的の信号を捕捉したかどうかを判定する判定工程と、
    　を有する信号サーチ方法。
  2. 　請求項１に記載の信号サーチ方法であって、
    　前記判定工程は、
    　前記各周波数のピーク相関値の最大値となる最大ピーク値を検出する工程と、
    　前記最大ピーク値を除く前記各周波数のピーク相関値から最大ノイズ値を決定する工程と、
    　前記最大ピーク値と前記最大ノイズ値との比を算出する工程と、を有し、
    　前記比に基づいて前記目的の信号を捕捉したかどうかを判定する、信号サーチ方法。
  3. 　請求項２に記載の信号サーチ方法であって、
    　前記最大ノイズ値を決定する工程は、
    　前記最大ピーク値を除く前記各周波数のピーク相関値における最大のピーク相関値を用いる、信号サーチ方法。
  4. 　請求項２に記載の信号サーチ方法であって、
    　前記最大ノイズ値を決定する工程は、
    　前記複数の周波数の相関値の平均値および標準偏差を算出する工程と、
    　コード位相のサーチ範囲に基づいて設定された定数を決定する工程と、を有し、
    　前記標準偏差に前記定数を乗算し、この乗算結果に前記平均値を加算することで、前記最大ノイズ値を決定する、信号サーチ方法。
  5. 　請求項２に記載の信号サーチ方法であって、
    　前記最大ノイズ値を決定する工程は、
    　前記最大ピーク値を除く前記各周波数のピーク相関値における最大のピーク相関値を前記最大ノイズ値に設定する第１の決定工程と、
    　前記複数の周波数の相関値の平均値および標準偏差を算出する工程と、コード位相のサーチ範囲に基づいて設定された定数を決定する工程と、を有し、前記標準偏差に前記定数を乗算し、この乗算結果に前記平均値を加算することで、前記最大ノイズ値を決定する第２の決定工程と、を有し、
    　前記相関値を算出するコード位相範囲が広い場合に、前記第１の決定工程で最大ノイズ値を決定し、
    　前記相関値を算出するコード位相範囲が狭い場合に、前記第２の決定工程で最大ノイズ値を決定する、信号サーチ方法。
  6. 　請求項１乃至請求項５のいずれかに記載の信号サーチ方法であって、
    　前記判定工程によって目的の信号を捕捉できていないと判定された時、前記複数の周波数を、前記周波数間隔よりも短い一定の周波数幅で周波数シフトさせる、信号サーチ方法。
  7. 　請求項１乃至請求項６のいずれかに記載の信号サーチ方法であって、
    　前記目的の信号は、ＧＰＳ衛星毎に放送されるＧＰＳ信号である、信号サーチ方法。
  8. 　目的の信号を捕捉するための信号サーチ処理をコンピュータに実行させるための信号サーチプログラムであって、
    　所定の周波数間隔に設定された複数の周波数で、前記目的の信号のレプリカ信号と受信信号との相関処理を行って相関値を取得する相関値取得処理と、
    　各周波数での最大相関値であるピーク相関値を取得するピーク相関値取得処理と、
    　前記各周波数のピーク相関値の大小関係に基づいて、前記目的の信号を捕捉したかどうかを判定する判定処理と、
    　が含まれる信号サーチプログラム。
  9. 　請求項８に記載の信号サーチプログラムであって、
    　前記判定処理には、
    　前記各周波数のピーク相関値の最大値となる最大ピーク値を検出する処理と、
    　前記最大ピーク値を除く前記各周波数のピーク相関値から最大ノイズ値を決定する処理と、
    　前記最大ピーク値と前記最大ノイズ値との比を算出する処理と、が含まれ、
    　前記比に基づいて前記目的の信号を捕捉したかどうかが判定される、信号サーチプログラム。
  10. 　目的の信号を捕捉するための信号サーチ装置であって、
    　所定の周波数間隔に設定された複数の周波数で、前記目的の信号のレプリカ信号と受信信号との相関処理を行って相関値を取得し、各周波数での最大相関値であるピーク相関値を取得するピーク相関値取得部と、
    　前記各周波数のピーク相関値の大小関係に基づいて、前記目的の信号を捕捉したかどうかを判定する判定部と、
    　を備える信号サーチ装置。
  11. 　請求項１０に記載の信号サーチ装置であって、
    　前記判定部は、
    　前記各周波数のピーク相関値の最大値となる最大ピーク値と、該最大ピーク値を除く前記各周波数のピーク相関値から決定される最大ノイズ値を取得し、前記最大ピーク値と前記最大ノイズ値との比に基づいて前記目的の信号を捕捉したかどうかを判定する、
    　信号サーチ装置。
  12. 　請求項１０乃至請求項１１のいずれかに記載の信号サーチ装置と、
    　前記相関値から前記目的の信号を捕捉、追尾する捕捉追尾部と、
    　追尾結果に基づいて測位を行う測位演算部と、を備えるＧＮＳＳ信号受信装置。
  13. 　請求項１２に記載のＧＮＳＳ信号受信装置と、
    　前記測位演算部の測位演算結果を用いて所定のアプリケーションを実行するアプリケーション処理部と、を備える情報機器端末。
PCT/JP2013/053889 2012-03-22 2013-02-18 信号サーチ方法、信号サーチプログラム、信号サーチ装置、ｇｎｓｓ信号受信装置、および情報機器端末 WO2013140910A1 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2012064633 2012-03-22
JP2012-064633 2012-03-22

Publications (1)

Publication Number Publication Date
WO2013140910A1 true WO2013140910A1 (ja) 2013-09-26

Family

ID=49222375

Family Applications (1)

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

Country Status (1)

Country Link
WO (1) WO2013140910A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112924956A (zh) * 2021-02-05 2021-06-08 中国人民解放军国防科技大学 一种降低低频噪声级的背景均衡方法
CN114760319A (zh) * 2020-12-28 2022-07-15 武汉联影智融医疗科技有限公司 Hbc数据同步方法、计算机设备和存储介质

Citations (7)

* 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 偽信号相互相関検出方法、送信源選択制限方法及び衛星選択制限方法
JP2003110462A (ja) * 2001-09-28 2003-04-11 Japan Radio Co Ltd 相関値識別方法
WO2004005954A2 (en) * 2002-07-10 2004-01-15 Qualcomm, Incorporated Cross-correlation mitigation method and apparatus for use in a global positioning system receiver
JP2008107353A (ja) * 2006-10-26 2008-05-08 Nemerix Sa 相互相関拒否を有するｇｎｓｓ受信機
JP2008160642A (ja) * 2006-12-26 2008-07-10 Seiko Epson Corp 相関演算制御回路及び相関演算制御方法
JP2009525006A (ja) * 2006-01-26 2009-07-02 クゥアルコム・インコーポレイテッド 位置特定受信機のための相互相関抑制技術
WO2011065559A1 (ja) * 2009-11-30 2011-06-03 古野電気株式会社 不要信号判別装置、不要信号判別方法、不要信号判別プログラム、ｇｎｓｓ受信装置および移動端末

Patent Citations (7)

* 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 偽信号相互相関検出方法、送信源選択制限方法及び衛星選択制限方法
JP2003110462A (ja) * 2001-09-28 2003-04-11 Japan Radio Co Ltd 相関値識別方法
WO2004005954A2 (en) * 2002-07-10 2004-01-15 Qualcomm, Incorporated Cross-correlation mitigation method and apparatus for use in a global positioning system receiver
JP2009525006A (ja) * 2006-01-26 2009-07-02 クゥアルコム・インコーポレイテッド 位置特定受信機のための相互相関抑制技術
JP2008107353A (ja) * 2006-10-26 2008-05-08 Nemerix Sa 相互相関拒否を有するｇｎｓｓ受信機
JP2008160642A (ja) * 2006-12-26 2008-07-10 Seiko Epson Corp 相関演算制御回路及び相関演算制御方法
WO2011065559A1 (ja) * 2009-11-30 2011-06-03 古野電気株式会社 不要信号判別装置、不要信号判別方法、不要信号判別プログラム、ｇｎｓｓ受信装置および移動端末

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114760319A (zh) * 2020-12-28 2022-07-15 武汉联影智融医疗科技有限公司 Hbc数据同步方法、计算机设备和存储介质
CN114760319B (zh) * 2020-12-28 2023-05-30 武汉联影智融医疗科技有限公司 Hbc数据同步方法、计算机设备和存储介质
CN112924956A (zh) * 2021-02-05 2021-06-08 中国人民解放军国防科技大学 一种降低低频噪声级的背景均衡方法
CN112924956B (zh) * 2021-02-05 2022-05-13 中国人民解放军国防科技大学 一种降低低频噪声级的背景均衡方法

Similar Documents

Publication Publication Date Title
RU2472172C2 (ru) Обнаружение многолучевого распространения для принимаемого sps-сигнала
KR101070211B1 (ko) 충분한 의사범위 측정치들을 결정하는 시스템 및/또는 방법
KR101102283B1 (ko) 수신된 ｓｐｓ 신호에서의 모호성을 감소시키기 위한 시스템 및/또는 방법
KR101120444B1 (ko) Ｇｎｓｓ 신호의 포착 시스템 및/또는 방법
US8571088B2 (en) Suppression of multipath effects for received SPS signals
EP2510627B1 (en) Methods and apparatuses for identifying and mitigating interference in a wireless signal
EP2755049A1 (en) Method for performing consistency checks for multiple signals received from a transmitter
US20120319898A1 (en) Method of detecting multipath, gnss receiving apparatus , and mobile terminal
KR20150131020A (ko) 도달 시간 결정을 개선하기 위한 방법 및 그 시스템
EP2793050A1 (en) Gnss signal processing method, positioning method, gnss signal processing program, positioning program, gnss signal processing device, positioning device, and mobile terminal
JP4828308B2 (ja) 位相変調系列再生装置
JP5918351B2 (ja) 信号サーチ方法、信号サーチプログラム、信号サーチ装置、ｇｎｓｓ信号受信装置、および情報機器端末
JP6061316B2 (ja) Ｇｎｓｓ受信機の位置を決定する方法及び装置
JP5483750B2 (ja) 不要信号判別装置、不要信号判別方法、不要信号判別プログラム、ｇｎｓｓ受信装置および移動端末
US9891323B2 (en) GNSS signal processing method, positioning method, GNSS signal processing program, positioning program, GNSS signal processing device, positioning apparatus and mobile terminal
EP3112905A1 (en) Gnss receiver calculating a non-ambiguous discriminator to resolve subcarrier tracking ambiguities
CN107820212B (zh) 一种基于移动多媒体广播系统实现的定位方法及定位终端
JP2013053972A (ja) 信号捕捉方法、通信信号受信方法、ｇｎｓｓ信号受信方法、信号捕捉プログラム、通信信号受信プログラム、ｇｎｓｓ信号受信プログラム、信号捕捉装置、通信信号受信装置、ｇｎｓｓ信号受信装置、および移動端末
WO2013140910A1 (ja) 信号サーチ方法、信号サーチプログラム、信号サーチ装置、ｇｎｓｓ信号受信装置、および情報機器端末
JP2008209287A (ja) 衛星航法受信機
JP5521284B2 (ja) 測位方法及び測位装置
WO2013140909A1 (ja) 信号サーチ方法、信号サーチプログラム、信号サーチ装置、ｇｎｓｓ信号受信装置、および情報機器端末
CN116953748A (zh) 多路径环境中的直接及反射导航卫星信号参数的检测及估计
JP4116916B2 (ja) スペクトラム拡散信号処理装置
JP2011153891A (ja) 衛星信号受信機

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13763649

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 13763649

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: JP