WO2018138830A1 - 信号捕捉装置 - Google Patents

信号捕捉装置 Download PDF

Info

Publication number
WO2018138830A1
WO2018138830A1 PCT/JP2017/002730 JP2017002730W WO2018138830A1 WO 2018138830 A1 WO2018138830 A1 WO 2018138830A1 JP 2017002730 W JP2017002730 W JP 2017002730W WO 2018138830 A1 WO2018138830 A1 WO 2018138830A1
Authority
WO
WIPO (PCT)
Prior art keywords
signal
code
code replica
signal acquisition
unit
Prior art date
Application number
PCT/JP2017/002730
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 三菱電機株式会社
Priority to JP2017532196A priority Critical patent/JP6253855B1/ja
Priority to US16/471,516 priority patent/US11169275B2/en
Priority to PCT/JP2017/002730 priority patent/WO2018138830A1/ja
Publication of WO2018138830A1 publication Critical patent/WO2018138830A1/ja

Links

Images

Classifications

    • 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
    • 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/246Acquisition or tracking or demodulation of signals transmitted by the system involving long acquisition integration times, extended snapshots of signals or methods specifically directed towards weak signal acquisition
    • 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/35Constructional details or hardware or software details of the signal processing chain
    • G01S19/37Hardware or software details of the signal processing chain
    • 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/7073Synchronisation aspects
    • H04B1/7075Synchronisation aspects with code phase acquisition
    • 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/7073Synchronisation aspects
    • H04B1/7085Synchronisation aspects using a code tracking loop, e.g. a delay-locked loop
    • 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

Definitions

  • the present invention relates to a signal acquisition device that performs signal detection and initial synchronization by obtaining a correlation between a received signal for one period of a ranging code and a code replica.
  • GNSS Global Navigation Satellite System
  • FPGA Field-Programmable Gate Array
  • the detection and detection of signals are performed by correlating the received signal for one period of the ranging code and the code replica
  • the initial synchronization of the code phase and the Doppler frequency is performed.
  • Cyclic convolution is widely used for correlation processing in signal acquisition (see Non-Patent Document 1, for example). Correlation processing using cyclic convolution is established because the distance measurement code included in the GNSS signal is a periodic signal.
  • a feature of this method is that the correlation processing can be realized with low computation.
  • Signal acquisition by cyclic convolution assumes the complete periodicity of the ranging code.
  • the polarity may change in the middle.
  • GNSS signals in which codes called overlay codes and secondary codes are superimposed on distance measurement codes have increased in recent years (for example, E5a signal of Galileo system).
  • the sign of the superposition code determines the polarity of the distance measurement code, so that the polarity change of the distance measurement code frequently occurs.
  • the conventional method based on cyclic convolution is greatly affected by the capture performance deterioration due to the polarity change. Therefore, as a countermeasure for this problem, there has been a method of obtaining a received signal having a distance ranging code period or more and obtaining a correlation with the code replica by linear convolution (see, for example, Non-Patent Document 2).
  • the method using the linear convolution as described above does not cause a decrease in capture probability or a Doppler frequency error due to a change in the polarity of the ranging code.
  • a received signal having a distance ranging code period or more is required, and the processing scale increases.
  • the present invention has been made to solve such a problem, and an object of the present invention is to provide a signal capturing device capable of preventing performance deterioration due to polarity change while suppressing an increase in processing scale.
  • the signal acquisition device includes a RF front end that performs frequency conversion and A / D conversion on a reception signal received from a signal transmitted from a satellite, and a polarity change in a ranging code with respect to an output from the RF front end
  • a signal acquisition unit that performs signal detection and initial synchronization by cyclic convolution operation processing using a first code replica corresponding to the case where there is no polarity and a second code replica corresponding to a case where there is a polarity change, and a signal acquisition unit
  • a signal tracking unit that performs synchronous tracking using the output signal acquisition result as an initial value.
  • the signal acquisition device performs cyclic convolution operation processing using a first code replica corresponding to a case where there is no polarity change in the ranging code and a second code replica corresponding to a case where there is a polarity change. It is what I do. As a result, it is possible to suppress the length of the received signal, suppress an increase in processing scale, and prevent performance deterioration due to polarity change.
  • 6B is an explanatory diagram showing a case where there is a change in polarity. It is explanatory drawing which shows the correlation value after the low-pass filter process in the signal capture device of Embodiment 2 of this invention. It is a block diagram of the calculator in the signal acquisition apparatus of Embodiment 3 of this invention.
  • FIG. 1 is a configuration diagram of a signal capturing device according to the present embodiment.
  • the illustrated signal capturing apparatus includes a receiving antenna 1, an RF front end 2, a signal capturing unit 3, and a signal tracking unit 4.
  • the receiving antenna 1 is an antenna that receives a signal transmitted from a satellite.
  • the RF front end 2 is a processing unit that performs frequency conversion and A / D conversion on a reception signal received by the reception antenna 1.
  • the signal acquisition unit 3 For the signal from the RF front end 2, the signal acquisition unit 3 includes a first code replica corresponding to the case where there is no polarity change in the distance measurement code and a second code replica corresponding to the case where the polarity change occurs. It is a processing unit that performs signal detection and initial synchronization by the used cyclic convolution operation processing.
  • the signal tracking unit 4 is a processing unit that performs synchronization tracking using the signal acquisition result output from the signal acquisition unit 3 as an initial value.
  • the signal capturing apparatus when a GNSS signal is received by the receiving antenna 1, frequency conversion and A / D conversion are performed by the RF front end 2 and transferred to the signal capturing unit 3.
  • the signal acquisition unit 3 detects the GNSS signal and performs initial synchronization, and outputs the result to the signal tracking unit 4.
  • the signal tracking unit 4 performs synchronous tracking of the GNSS signal using the signal acquisition result as an initial value.
  • the signal acquisition unit 3 includes an input interface 5 that inputs a received signal after A / D conversion, an output interface 6 that outputs a signal acquisition result to the signal tracking unit 4, a calculator 7 that performs a signal acquisition calculation, and a code replica And the like.
  • the computing unit 7 for example, a CPU (Central Processing Unit) or an FPGA (Field Programmable Gate Array) can be used.
  • FIG. 2 is a functional configuration diagram of the computing unit 7.
  • the calculator 7 includes a received signal input unit 9, an oscillator 10, a code replica generator 11, a multiplier 12, a convolution calculator 13, and a signal detector 14.
  • the reception signal input unit 9 acquires reception signals for one period of the ranging code.
  • the oscillator 10 generates a reference signal for frequency conversion.
  • the code replica generator 11 generates a distance measurement code replica and stores it in the memory 8.
  • the multiplier 12 mixes the received signal acquired by the received signal input unit 9 with the reference signal generated by the oscillator 10 and converts the frequency of the received signal.
  • the convolution calculator 13 calculates a correlation value between the frequency-converted signal and the code replica generated by the code replica generator 11 and held in the memory 8 by using cyclic convolution.
  • the convolution operation unit 13 repeatedly performs these processes while changing the frequency of the reference signal generated by the oscillator 10.
  • the signal detection unit 14 determines whether or not the GNSS signal is received from the maximum intensity of the obtained correlation value, and if it is determined that the GNSS signal is received, the Doppler frequency and code phase of the GNSS signal from the peak position of the correlation value. Get.
  • FIG. 3 shows a flowchart of processing in the computing unit 7.
  • the feature of this embodiment is that, in addition to the first code replica that has been conventionally used in the code replica generator 11, a correlation value can be obtained without loss when there is no change in the polarity of the ranging code in the received signal.
  • a second code replica that can obtain a correlation value without loss is generated, and cyclic convolution between the two types of code replicas and the received signal is performed. Yes (steps ST3 and ST4).
  • the first code replica r 1 [n] is represented by the following equation using the ranging code c (t) ( ⁇ 1 or +1 binary signal) transmitted from the GNSS satellite.
  • a second code replica represented by the following formula is used.
  • M is an arbitrary integer.
  • the second code replica r 2 [n] is obtained by applying a unique phase rotation to the ranging code for each sampling point.
  • a reception signal for the code period is acquired by the reception signal input unit 9 (step ST1), and the multiplication unit 12 multiplies the reception signal by the reference signal from the oscillator 10 (step ST2) and gives the result to the convolution operation unit 13.
  • A is the signal amplitude
  • N ⁇ is the sample position (code phase) at the beginning of the ranging code period
  • p [n] is the polarity (+1 or ⁇ 1) of the ranging code
  • f IF is the IF center frequency.
  • F d is the Doppler frequency
  • f is the oscillation frequency of the oscillator 10
  • is the initial phase.
  • the convolution operation unit 13 performs a cyclic convolution operation on such an input signal using the first code replica and the second code replica stored in the memory 8 (step ST3 and step ST4).
  • polar p [n] can only be varied at the boundary position N tau of ranging code (. See bit change in an arrow showing an example in FIG. 4). Therefore, the following equation holds.
  • p 1 and P 2 are scalar values that take values of +1 or ⁇ 1.
  • equation (8) By substituting equation (8) into equations (6) and (7), the following equation is obtained.
  • the correlation value when there is no change in polarity in the received signal used for calculating the correlation value, the correlation value can be obtained without any loss by cyclic convolution using the first code replica. (
  • AN))
  • the correlation value is obtained without any loss by cyclic convolution using the second code replica. It can be seen that the value is obtained (
  • AN)).
  • the convolution operation unit 13 holds the maximum value of the cyclic convolution output using the first code replica and the second code replica (step ST5), so that the presence or absence of the polarity change can be determined while using the cyclic convolution. Regardless of this, a correlation value without loss can be obtained.
  • the convolution operation unit 13 changes the reference signal generated by the oscillator 10 for such a circular convolution operation (step ST7), and searches for all Doppler frequencies (step ST6).
  • the signal detection unit 14 compares the maximum value obtained by the convolution operation unit 13 with a predetermined detection threshold, and the maximum value satisfies the detection threshold. If so, it is determined that the GNSS signal has been received (step ST8).
  • the search range of the Doppler frequency is a value determined from the assumed satellite speed and receiver speed, and for example, about ⁇ several kHz is used.
  • the maximum value is held in the above-described step ST5
  • a value that is equal to or larger than a predetermined set value including the maximum value may be held, and one of the values may be selected.
  • the RF front end that performs frequency conversion and A / D conversion on the received signal that has received the signal transmitted from the satellite, and the output from the RF front end are used.
  • signal detection and initial synchronization are performed by cyclic convolution operation processing using the first code replica corresponding to the case where there is no polarity change in the distance measurement code and the second code replica corresponding to the case where there is a polarity change. Since it has a signal acquisition unit and a signal tracking unit that performs synchronous tracking using the signal acquisition result output from the signal acquisition unit as an initial value, the received signal length is suppressed to prevent an increase in processing scale, and due to polarity change Performance degradation can be prevented.
  • the first code replica and the second code replica are generated by giving a unique phase rotation to the ranging code for each sampling point. Therefore, the first code replica and the second code replica can be easily obtained.
  • the signal acquisition unit has a set value equal to or greater than a set value among the correlation values obtained by the cyclic convolution operation process using the first code replica and the second code replica. Since any value is selected and signal detection is performed, a correlation value without loss can be obtained regardless of the presence or absence of polarity change.
  • Embodiment 2 the maximum value of each output of cyclic convolution using the first code replica and the second code replica is selected and the signal is detected, but after adding both outputs (non-coherent integration) It is also possible to detect a signal, which will be described below as a second embodiment.
  • the signal acquisition unit 3 of the second embodiment is configured to perform signal detection by adding the correlation values obtained by the cyclic convolution operation process using the first code replica and the second code replica. . Since the receiving antenna 1, the RF front end 2 and the signal tracking unit 4 are the same as those in the first embodiment, description thereof is omitted here.
  • FIG. 5 is a flowchart showing the operation of the arithmetic unit 7 in the second embodiment. Steps ST1 to ST4 are the same as in the first embodiment, and a description thereof is omitted here. Next, in step ST5 of the first embodiment, only the maximum value is held for the correlation values obtained in step ST3 and step ST4. However, in the second embodiment, the convolution operation unit 13 holds all correlation values. (Step ST9).
  • the cyclic convolution using the first code replica when there is a polarity change in the received signal and the cyclic convolution using the second code replica when there is no polarity change in the received signal are both true Doppler frequencies. Correlation power is scattered around and the correlation value is not maximized at the true Doppler frequency.
  • FIG. 6 shows the image of the correlation output when the first code replica is used.
  • FIG. 6A shows a case where there is no change in polarity
  • FIG. 6B shows a state where there is a change in polarity.
  • the first code replica is used, as shown in FIG. 6B, a phenomenon occurs in which the correlation output is scattered around the true Doppler frequency.
  • the second code replica is used, the states of FIGS. 6A and 6B are reversed.
  • p 1 and p 2 are scalar values that take a value of +1 or ⁇ 1 in equation (8).
  • the above equation (12) is as follows.
  • step ST10 A low pass filter is applied in the Doppler frequency direction (step ST10).
  • the correlation value peak becomes a true Doppler frequency.
  • FIG. 7 is an explanatory diagram showing this processing. Since the low-pass filter is applied in the Doppler frequency direction, the correlation value peak is in the vicinity of the true Doppler frequency (0 Hz) as shown by the curve 100 in the figure. This process also has an effect of reducing the side lobe of the Doppler frequency.
  • the convolution calculation unit 13 performs non-coherent integration of both correlation values (step ST11), and the signal detection unit 14 performs signal acquisition by comparing the maximum value with a detection threshold (step ST8).
  • the correlation value peak is accumulated at the true signal position, while the dispersion other than the true signal position is reduced. Therefore, the detection probability is improved as compared with the signal detection using the single correlation value.
  • a moving average filter etc. can also be used as a low-pass filter of said step ST10.
  • the signal acquisition unit uses the correlation values obtained by the cyclic convolution operation processing using the first code replica and the second code replica. Since addition is performed to perform signal detection, the detection probability in signal detection can be improved.
  • the signal acquisition unit adds the correlation values after applying a low-pass filter in the Doppler frequency direction, so that the correlation value peak can be obtained more accurately. Furthermore, the detection probability in signal detection can be improved.
  • Embodiment 3 A signal acquisition processing configuration in which cyclic convolution operation in GNSS signal acquisition is performed using DFT (Discrete Fourier Transform) or FFT (Fast Fourier Transform) is called PCS (Parallel Code-Phase Search). Therefore, it is widely used in GNSS signal acquisition devices (for example, see Non-Patent Document 2).
  • the signal capturing apparatus according to the first embodiment and the second embodiment can be applied to a PCS configuration, which will be described below as a third embodiment.
  • FIG. 8 is a configuration diagram illustrating a main part of the signal capturing device according to the third embodiment.
  • This signal acquisition device is a configuration example by PCS using FFT, and shows the configuration of the third embodiment corresponding to the arithmetic unit 7 shown in FIG.
  • the oscillator 10 of the third embodiment is configured to output reference signals f 1 to f K
  • the multiplier 12 includes multipliers 121-1 to 121-K, and these multipliers 121-1 to 121-K.
  • the convolution operation unit 13 is configured to perform a convolution operation on each of the plurality of signals output from the multiplication unit 12.
  • the FFT units 131-1 to 131-K corresponding to the plurality of output signals from the multiplication unit 12 and the multipliers 132a-1 to 132a-K for multiplying the first code replica stored in the memory 8, respectively.
  • Multipliers 132b-1 to 132b-K that respectively multiply the second code replicas stored in the memory 8
  • IFFT units 133a-1 to 133a- corresponding to the outputs of the multipliers 132a-1 to 132a-K.
  • K and IFFT units 133b-1 to 133b-K corresponding to the outputs of the multipliers 132b-1 to 132b-K.
  • the first code replica and the second code replica stored in the memory 8 are code replicas after the code replica generated by the code replica generator 11 is subjected to FFT processing by the FFT unit 15.
  • the IFFT units 133a-1 to 133a-K and 133b-1 to 133b-K are processing units that perform inverse FFT (Inverse FFT) in the FFT units 131-1 to 131-K.
  • the arithmetic unit 7 obtains the correlation value of all code phases using the FFT for each oscillator frequency from the oscillator 10 and outputs the correlation value to the signal detection unit 14.
  • a cyclic convolution operation can be realized with a low calculation amount and high speed.
  • the signal acquisition unit performs the cyclic convolution operation processing using the fast Fourier transform or the discrete Fourier transform. It can be realized in quantity and at high speed.
  • the signal acquisition device relates to a configuration for performing signal acquisition using a cyclic convolution operation, and is suitable for use in a GNSS signal acquisition device.

Abstract

信号捕捉部(3)は、RFフロントエンド(2)からの出力に対して、測距コードに極性変化がない場合に対応した第1のコードレプリカと極性変化がある場合に対応した第2のコードレプリカとを用いた巡回畳み込み演算処理により信号検出及び初期同期を行う。信号追尾部(4)は、信号捕捉部(3)が出力した信号捕捉結果を初期値として同期追尾を行う。

Description

信号捕捉装置
 本発明は、測距コード1周期分の受信信号とコードレプリカとの相関を取ることで信号検出と初期同期を行う信号捕捉装置に関する。
 ソフトウェアやFPGA(Field-Programmable Gate Array)で構成されたGNSS(Global Navigation Satellite System)信号捕捉装置では、測距コード1周期分の受信信号とコードレプリカとの相関を取ることで、信号の検出と、コード位相及びドップラ周波数の初期同期を行う。信号捕捉における相関処理には、巡回畳み込みが広く利用されている(例えば、非特許文献1参照)。巡回畳み込みを用いた相関処理は、GNSS信号に含まれる測距コードが周期性のある信号であるため成り立つ。本手法の特徴として、相関処理が低演算で実現可能なことが挙げられる。
 巡回畳み込みによる信号捕捉は、測距コードの完全な周期性を前提としている。しかし、GNSS信号の測距コードは、周期的ではあるものの途中で極性が変化することがある。測距コードの極性変化が発生している受信信号を用いて巡回畳み込みを行った場合、相関値低下や相関ピーク位置のずれが発生し、検出確率低下やドップラ周波数誤りを引き起こす。
 従来から用いられてきた民生用コードであるC/Aコードでは、航法メッセージのデータ符号が測距コードの極性を決定することから、航法メッセージビット境界(1航法メッセージビット長=測距コード20周期分)でのみ極性変化が起こり得る。そのため、極性変化の発生頻度が低く、巡回畳み込みを用いた信号捕捉に対する影響を無視することも可能であった。
 一方、測距コードに対し、オーバーレイコードやセカンダリコードと呼ばれるコードが重畳されたGNSS信号が近年増加している(例えば、GalileoシステムのE5a信号)。このようなGNSS信号においては、航法メッセージに加えて、重畳コードの符号が測距コードの極性を決定するため、測距コードの極性変化が頻繁に発生する。そのため、巡回畳み込みによる従来手法では、極性変化による捕捉性能劣化の影響を大きく受けてしまう。そこで、この問題の対策として、測距コード周期以上の受信信号を取得し、コードレプリカとの相関を線形畳み込みで求める手法があった(例えば、非特許文献2参照)。
J. B. Y. Tsui著、「Fundamentals of Global Positioning System Receivers: A Software Approach」John Wiley & Sons, Inc.、2000年発行、p.140~144 C. Yang, "FFT Acquisition of Periodic, Aperiodic, Puncture, and Overlaid Code Sequences in GPS," Proceedings of the ION GPS 2001, p.137-147, 2011.
 上記のような線形畳み込みを用いた手法では、測距コードの極性変化による捕捉確率低下やドップラ周波数誤りは発生しない。しかしながら、測距コード周期以上の受信信号が必要となり、処理規模も増大するという課題があった。
 この発明は、かかる問題を解決するためになされたもので、処理規模の増大を抑えながら、極性変化による性能劣化を防ぐことのできる信号捕捉装置を提供することを目的とする。
 この発明に係る信号捕捉装置は、衛星から送信された信号を受信した受信信号を周波数変換及びA/D変換するRFフロントエンドと、RFフロントエンドからの出力に対して、測距コードに極性変化がない場合に対応した第1のコードレプリカと極性変化がある場合に対応した第2のコードレプリカとを用いた巡回畳み込み演算処理により信号検出及び初期同期を行う信号捕捉部と、信号捕捉部が出力した信号捕捉結果を初期値として同期追尾を行う信号追尾部とを備えたものである。
 この発明に係る信号捕捉装置は、測距コードに極性変化がない場合に対応した第1のコードレプリカと、極性変化がある場合に対応した第2のコードレプリカとを用いて巡回畳み込み演算処理を行うようにしたものである。これにより、受信信号長を抑えて処理規模の増大を抑制し、かつ、極性変化による性能劣化を防ぐことができる。
この発明の実施の形態1による信号捕捉装置の構成図である。 この発明の実施の形態1の信号捕捉装置における演算器の構成図である。 この発明の実施の形態1の信号捕捉装置における演算器の処理を示すフローチャートである。 測距コードの極性変化を示す説明図である。 この発明の実施の形態2の信号捕捉装置における演算器の処理を示すフローチャートである。 図6Aは、この発明の実施の形態2の信号捕捉装置における第1のコードレプリカを用いたときの極性変化無しの場合、図6Bは極性変化ありの場合を示す説明図である。 この発明の実施の形態2の信号捕捉装置におけるローパスフィルタ処理後の相関値を示す説明図である。 この発明の実施の形態3の信号捕捉装置における演算器の構成図である。
 以下、この発明をより詳細に説明するために、この発明を実施するための形態について、添付の図面に従って説明する。
実施の形態1.
 図1は、本実施の形態による信号捕捉装置の構成図である。
 図示の信号捕捉装置は、受信アンテナ1とRFフロントエンド2と信号捕捉部3と信号追尾部4とを備える。受信アンテナ1は衛星から送信された信号を受信するアンテナである。RFフロントエンド2は、受信アンテナ1で受信された受信信号を周波数変換すると共にA/D変換する処理部である。信号捕捉部3は、RFフロントエンド2からの信号に対して、測距コードに極性変化がない場合に対応した第1のコードレプリカと極性変化がある場合に対応した第2のコードレプリカとを用いた巡回畳み込み演算処理により信号検出及び初期同期を行う処理部である。信号追尾部4は、信号捕捉部3が出力した信号捕捉結果を初期値として同期追尾を行う処理部である。
 すなわち、実施の形態1の信号捕捉装置では、受信アンテナ1でGNSS信号を受信した場合、RFフロントエンド2で周波数変換及びA/D変換され、信号捕捉部3に転送される。信号捕捉部3ではGNSS信号の検出及び初期同期を行い、結果を信号追尾部4に出力する。信号追尾部4は、信号捕捉結果を初期値として、GNSS信号の同期追尾を実施する。
 信号捕捉部3は、A/D変換後の受信信号を入力する入力インタフェース5と、信号捕捉結果を信号追尾部4に出力する出力インタフェース6と、信号捕捉演算を行う演算器7と、コードレプリカなどを保持するメモリ8で構成されている。演算器7には、例えばCPU(Central Processing Unit)やFPGA(Field Programmable Gate Array)を用いることができる。
 図2は、演算器7の機能構成図である。演算器7は、受信信号入力部9、オシレータ10、コードレプリカ生成器11、乗算部12、畳み込み演算部13、信号検出部14により構成される。
 受信信号入力部9では、測距コード1周期分の受信信号を取得する。オシレータ10では、周波数変換用の基準信号を生成する。コードレプリカ生成器11では、測距コードのレプリカを生成し、メモリ8に保存する。乗算部12では、受信信号入力部9で取得した受信信号と、オシレータ10で生成した基準信号とをミキシングし、受信信号を周波数変換する。畳み込み演算部13では、周波数変換後の信号と、コードレプリカ生成器11で生成しメモリ8に保持しておいたコードレプリカとの相関値を、巡回畳み込みを用いて算出する。畳み込み演算部13では、これらの処理を、オシレータ10で生成する基準信号の周波数を変化させながら繰り返し実施する。信号検出部14では、得られた相関値の最大強度からGNSS信号の受信有無を判定し、GNSS信号の受信を判定した場合であれば、相関値のピーク位置からGNSS信号のドップラ周波数及びコード位相を得る。
 演算器7における処理のフローチャートを図3に示す。本実施の形態の特徴は、コードレプリカ生成器11において、従来から用いられてきた、受信信号中に測距コードの極性変化がない場合に損失なく相関値が得られる第1のコードレプリカに加えて、受信信号中に測距コードの極性変化がある場合に損失なく相関値が得られる第2のコードレプリカを生成し、2種類のコードレプリカと受信信号との巡回畳み込みをそれぞれ実施することにある(ステップST3、ST4)。
 従来手法で用いられる第1のコードレプリカを、r[n](n=0~N-1、nはサンプル番号)とする。ここで、Nは測距コード1周期分のサンプル数である。すなわち、Tを測距コード周期、Tをサンプリング周期とした場合、N=T/Tと表される。このとき、第1のコードレプリカr[n]は、GNSS衛星から送信される測距コードc(t)(-1または+1の2値信号)を用いて、下式のとおり表される。
Figure JPOXMLDOC01-appb-I000001
 本実施の形態では、上式の第1のコードレプリカに加えて、下式で表される第2のコードレプリカを用いる。ここで、Mは任意の整数である。すなわち第2のコードレプリカr[n]は、測距コードに対して、サンプリング点毎にそれぞれ固有の位相回転を与えたものである。
Figure JPOXMLDOC01-appb-I000002
 先ず、受信信号入力部9でコード周期分の受信信号を取得し(ステップST1)、乗算部12で受信信号にオシレータ10からの基準信号を乗算して(ステップST2)、畳み込み演算部13に与える。
 畳み込み演算部13の入力信号となる測距コード1周期(=T)分の受信信号中のGNSS信号成分s(n,f)(n=0~N-1)は下式のとおり表すことができる。ここで、Aは信号振幅を、Nτは測距コード周期先頭のサンプル位置(コード位相)を、p[n]は測距コードの極性(+1または-1)を、fIFはIF中心周波数を、fはドップラ周波数を、fはオシレータ10の発振周波数を、Φは初期位相を示す。
Figure JPOXMLDOC01-appb-I000003
 畳み込み演算部13は、このような入力信号に対して、メモリ8に格納されている第1のコードレプリカと第2のコードレプリカとを用いて巡回畳み込み演算を行う(ステップST3及びステップST4)。
Figure JPOXMLDOC01-appb-I000004

Figure JPOXMLDOC01-appb-I000005

Figure JPOXMLDOC01-appb-I000006
Figure JPOXMLDOC01-appb-I000007

Figure JPOXMLDOC01-appb-I000008

Figure JPOXMLDOC01-appb-I000009
 ここで、極性p[n]は、測距コードの境界位置Nτでのみ変化し得る(図4に例を示す。図中矢印のビット変化参照)。そのため、下式が成り立つ。ただし、p及びPは+1または-1の値を取るスカラ値である。
Figure JPOXMLDOC01-appb-I000010
 式(6)及び式(7)に式(8)を代入することで、下式が求まる。下式は、極性変化がない場合(p=p)、及び極性変化がある場合(p=-p)それぞれの相関値を示している。
Figure JPOXMLDOC01-appb-I000011

Figure JPOXMLDOC01-appb-I000012
 上記の式(9)及び式(10)により、相関値算出に用いる受信信号中に極性変化がない場合は、第1のコードレプリカを用いた巡回畳み込みにより一切の損失なく相関値が求まることがわかる(|(R(Nτ,f)|=AN)。一方、受信信号中に極性変化が発生している場合は、第2のコードレプリカを用いた巡回畳み込みにより一切の損失なく相関値が求まることがわかる(|(R(Nτ,f)|=AN)。
 従って、畳み込み演算部13は、第1のコードレプリカ及び第2のコードレプリカを用いた巡回畳み込み出力の最大値を保持する(ステップST5)ことで、巡回畳み込みを利用しながら、極性変化の有無に依らず損失のない相関値を得ることができる。
 畳み込み演算部13は、このような巡回畳み込み演算を、オシレータ10で生成する基準信号を変更して(ステップST7)、全ドップラ周波数を探索する(ステップST6)。ステップST6において、全ドップラ周波数が探索済みであった場合、信号検出部14は、畳み込み演算部13で得られた最大値と、予め定めた検出閾値とを比較し、最大値が検出閾値を満たしていた場合が、これをGNSS信号の受信と判定する(ステップST8)。ここで、ドップラ周波数の探索範囲は、想定する衛星速度と受信機速度から決定される値であり、例えば±数kHz程度が用いられる。
 なお、第1のコードレプリカとして従来用いられてきた式(1)を挙げたが、第1のコードレプリカとして下式を用いても良い。ただし、Mは任意の整数である。
Figure JPOXMLDOC01-appb-I000013
 また、上記ステップST5において最大値を保持するとしたが、最大値を含む予め決められた設定値以上の値を保持し、その中からいずれかの値を選択するようにしてもよい。
 以上説明したように、実施の形態1の信号捕捉装置によれば、衛星から送信された信号を受信した受信信号を周波数変換及びA/D変換するRFフロントエンドと、RFフロントエンドからの出力に対して、測距コードに極性変化がない場合に対応した第1のコードレプリカと極性変化がある場合に対応した第2のコードレプリカとを用いた巡回畳み込み演算処理により信号検出及び初期同期を行う信号捕捉部と、信号捕捉部が出力した信号捕捉結果を初期値として同期追尾を行う信号追尾部とを備えたので、受信信号長を抑えて処理規模の増大を抑制し、かつ、極性変化による性能劣化を防ぐことができる。
 また、実施の形態1の信号捕捉装置によれば、第1のコードレプリカ及び第2のコードレプリカは、測距コードに対し、サンプリング点毎にそれぞれ固有の位相回転を与えて生成するようにしたので、容易に第1のコードレプリカ及び第2のコードレプリカを求めることができる。
 また、実施の形態1の信号捕捉装置によれば、信号捕捉部は、第1のコードレプリカと第2のコードレプリカとを用いた巡回畳み込み演算処理により得られた相関値のうち設定値以上のいずれかの値を選択し、信号検出を行うようにしたので、極性変化の有無に依らず損失のない相関値を得ることができる。
実施の形態2.
 実施の形態1では、第1のコードレプリカ及び第2のコードレプリカを用いた巡回畳み込みの各出力の最大値を選択し、信号検出していたが、両出力を加算(ノンコヒーレント積分)した後に信号検出することも可能であり、これを実施の形態2として次に説明する。
 実施の形態2における信号捕捉装置の図面上の構成については図1及び図2の構成と同様であるため、これらの図面の構成を用いて説明する。実施の形態2の信号捕捉部3は、第1のコードレプリカと第2のコードレプリカとを用いた巡回畳み込み演算処理により得られた相関値同士を加算し、信号検出を行うよう構成されている。受信アンテナ1、RFフロントエンド2及び信号追尾部4については実施の形態1と同様であるため、ここでの説明は省略する。
 図5は、実施の形態2における演算器7の動作を示すフローチャートである。ステップST1~ステップST4は、実施の形態1と同様であるため、ここでの説明は省略する。次に、実施の形態1のステップST5では、ステップST3及びステップST4で求めた相関値に対し最大値のみを保持していたが、実施の形態2では畳み込み演算部13が全相関値を保持しておく(ステップST9)。
 受信信号中に極性変化があるときの第1のコードレプリカを用いた巡回畳み込み及び受信信号中に極性変化がないときの第2のコードレプリカを用いた巡回畳み込みでは、いずれも真のドップラ周波数の周辺に相関電力が散らばり、真のドップラ周波数において相関値が最大とならない特徴を持つ。
 図6に第1のコードレプリカを用いたときの相関出力のイメージを示す。図6Aが極性変化なしの場合であり、図6Bが極性変化ありの状態を示す。第1のコードレプリカを用いた場合、図6Bに示すように,真のドップラ周波数周辺に相関出力が散らばる現象が発生する。第2のコードレプリカを用いた場合は、図6Aと図6Bの状態が逆になる。
Figure JPOXMLDOC01-appb-I000014

Figure JPOXMLDOC01-appb-I000015
 上式において、p及びpは,式(8)の+1または-1の値を取るスカラ値である。ここで、極性変化がある場合(p=-p)、上式(12)は以下の通りとなる。
Figure JPOXMLDOC01-appb-I000016

Figure JPOXMLDOC01-appb-I000017
 そこで、実施の形態2の畳み込み演算部13では、ステップST6で全ドップラ周波数を探索済みであった場合、次に、第1のコードレプリカによる相関値及び第2のコードレプリカによる相関値に対し、ドップラ周波数方向にローパスフィルタをかける(ステップST10)。この処理により、相関値ピークが真のドップラ周波数となる。図7はこの処理を示す説明図である。ドップラ周波数方向にローパスフィルタをかけたため、図中の曲線100に示すように、相関値ピークが真のドップラ周波数(0Hz)付近となる。また、この処理は、ドップラ周波数のサイドローブを低減させる効果もある。
 その後、畳み込み演算部13は、両相関値をノンコヒーレント積分し(ステップST11)、信号検出部14は、その最大値を検出閾値と比較することで信号捕捉を行う(ステップST8)。ノンコヒーレント積分により、相関値ピークを真の信号位置に積み上げる一方、真の信号位置以外の分散が低減されるため、単独相関値での信号検出と比較して検出確率を向上させる効果がある。
 なお、上記ステップST10のローパスフィルタとして、移動平均フィルタなどを用いることもできる。
 以上説明したように、実施の形態2の信号捕捉装置によれば、信号捕捉部は、第1のコードレプリカと第2のコードレプリカとを用いた巡回畳み込み演算処理により得られた相関値同士を加算し、信号検出を行うようにしたので、信号検出における検出確率を向上させることができる。
 また、実施の形態2の信号捕捉装置によれば、信号捕捉部は、ドップラ周波数方向にローパスフィルタをかけた後に相関値同士の加算を行うようにしたので、相関値ピークをより正確に求めることができ、さらに、信号検出における検出確率を向上させることができる。
実施の形態3.
 GNSS信号捕捉における巡回畳み込み演算を、DFT(Discrete Fourier Transform:離散フーリエ変換)やFFT(Fast Fourier Transform:高速フーリエ変換)を用いて実現した信号捕捉処理構成はPCS(Parallel Code-phase Search)と呼ばれ、GNSS信号捕捉装置に広く用いられている(例えば、非特許文献2参照)。上記実施の形態1及び実施の形態2の信号捕捉装置は、PCS構成に適用することも可能であり、これを実施の形態3として以下に説明する。
 図8は、実施の形態3の信号捕捉装置の要部を示す構成図である。この信号捕捉装置は、FFTを用いたPCSによる構成例であり、図2で示した演算器7に対応した実施の形態3の構成を示している。実施の形態3のオシレータ10は、基準信号f~fを出力するよう構成され、乗算部12は、乗算器121-1~121-Kを備え、これら乗算器121-1~121-Kによって、オシレータ10からの基準信号f~fと受信信号入力部9からの受信信号とをそれぞれ乗算するよう構成されている。畳み込み演算部13は、乗算部12から出力された複数の信号に対してそれぞれ畳み込み演算を行うよう構成されている。すなわち、乗算部12からの複数の出力信号に対応したFFT部131-1~131-Kと、メモリ8に格納されている第1のコードレプリカをそれぞれ乗算する乗算器132a-1~132a-K、メモリ8に格納されている第2のコードレプリカをそれぞれ乗算する乗算器132b-1~132b-Kと、乗算器132a-1~132a-Kの出力に対応したIFFT部133a-1~133a-Kと、乗算器132b-1~132b-Kの出力に対応したIFFT部133b-1~133b-Kとを備えている。なお、メモリ8に格納される第1のコードレプリカと第2のコードレプリカは、コードレプリカ生成器11で生成されたコードレプリカをFFT部15でFFT処理した後のコードレプリカである。また、IFFT部133a-1~133a-K,133b-1~133b-Kは、FFT部131-1~131-KにおけるFFTの逆変換(Inverse FFT)を行う処理部である。
 このように、実施の形態3の演算器7では、オシレータ10からのオシレータ周波数毎にFFTを用いて全コード位相の相関値を求め、信号検出部14に出力する。本構成により、巡回畳み込み演算が低演算量かつ高速に実現できる。
 なお、図6ではFFTを用いた構成としているが、もちろんDFTを用いた構成としてもかまわない。
 以上説明したように、実施の形態3の信号捕捉装置によれば、信号捕捉部は、高速フーリエ変換または離散フーリエ変換を用いて巡回畳み込み演算処理を行うようにしたので、巡回畳み込み演算を低演算量でかつ高速に実現することができる。
 なお、本願発明はその発明の範囲内において、各実施の形態の自由な組み合わせ、あるいは各実施の形態の任意の構成要素の変形、もしくは各実施の形態において任意の構成要素の省略が可能である。
 以上のように、この発明に係る信号捕捉装置は、巡回畳み込み演算を用いて信号捕捉を行う構成に関するものであり、GNSS信号捕捉装置に用いるのに適している。
 1 受信アンテナ、2 RFフロントエンド、3 信号捕捉部、4 信号追尾部、5 入力インタフェース、6 出力インタフェース、7 演算器、8 メモリ、9 受信信号入力部、10 オシレータ、11 コードレプリカ生成器、12 乗算部、13 畳み込み演算部、14 信号検出部、15,131-1~131-K FFT部、121-1~121-K,132a-1~132a-K,132b-1~132b-K 乗算器、133a-1~133a-K,133b-1~133b-K IFFT部。

Claims (6)

  1.  衛星から送信された信号を受信した受信信号を周波数変換及びA/D変換するRFフロントエンドと、
     前記RFフロントエンドからの出力に対して、測距コードに極性変化がない場合に対応した第1のコードレプリカと極性変化がある場合に対応した第2のコードレプリカとを用いた巡回畳み込み演算処理により信号検出及び初期同期を行う信号捕捉部と、
     前記信号捕捉部が出力した信号捕捉結果を初期値として同期追尾を行う信号追尾部とを備えたことを特徴とする信号捕捉装置。
  2.  前記第1のコードレプリカ及び第2のコードレプリカは、測距コードに対し、サンプリング点毎にそれぞれ固有の位相回転を与えて生成することを特徴とする請求項1記載の信号捕捉装置。
  3.  前記信号捕捉部は、前記第1のコードレプリカと前記第2のコードレプリカとを用いた巡回畳み込み演算処理により得られた相関値のうち設定値以上のいずれかの値を選択し、前記信号検出を行うことを特徴とする請求項1記載の信号捕捉装置。
  4.  前記信号捕捉部は、前記第1のコードレプリカと前記第2のコードレプリカとを用いた巡回畳み込み演算処理により得られた相関値同士を加算し、信号検出を行うことを特徴とする請求項1記載の信号捕捉装置。
  5.  前記信号捕捉部は、ドップラ周波数方向にローパスフィルタをかけた後に前記相関値同士の加算を行うことを特徴とする請求項4記載の信号捕捉装置。
  6.  前記信号捕捉部は、高速フーリエ変換または離散フーリエ変換を用いて巡回畳み込み演算処理を行うことを特徴とする請求項1記載の信号捕捉装置。
PCT/JP2017/002730 2017-01-26 2017-01-26 信号捕捉装置 WO2018138830A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2017532196A JP6253855B1 (ja) 2017-01-26 2017-01-26 信号捕捉装置
US16/471,516 US11169275B2 (en) 2017-01-26 2017-01-26 Signal acquisition device
PCT/JP2017/002730 WO2018138830A1 (ja) 2017-01-26 2017-01-26 信号捕捉装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2017/002730 WO2018138830A1 (ja) 2017-01-26 2017-01-26 信号捕捉装置

Publications (1)

Publication Number Publication Date
WO2018138830A1 true WO2018138830A1 (ja) 2018-08-02

Family

ID=60860167

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2017/002730 WO2018138830A1 (ja) 2017-01-26 2017-01-26 信号捕捉装置

Country Status (3)

Country Link
US (1) US11169275B2 (ja)
JP (1) JP6253855B1 (ja)
WO (1) WO2018138830A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013116408A1 (en) 2012-01-30 2013-08-08 California Institute Of Technology Filtration membranes and related compositions, methods and systems

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010060436A (ja) * 2008-09-04 2010-03-18 Furuno Electric Co Ltd Gnss受信装置および測位装置
US20100141520A1 (en) * 2007-08-21 2010-06-10 Giorgio Ghinamo Method for the acquisition of signals of a global navigation satellite system
JP2011202958A (ja) * 2010-03-24 2011-10-13 Seiko Epson Corp 信号捕捉方法、信号捕捉装置及び電子機器
JP2014216871A (ja) * 2013-04-26 2014-11-17 三菱電機株式会社 信号処理装置、信号処理方法及び信号処理プログラム

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4755920B2 (ja) * 2006-02-23 2011-08-24 古野電気株式会社 キャリア位相追尾装置および擬似雑音コード信号追尾装置
US8525727B2 (en) * 2009-12-29 2013-09-03 Texas Instruments Incorporated Position and velocity uncertainty metrics in GNSS receivers
US8855171B2 (en) * 2011-08-10 2014-10-07 Qualcomm Incorporated Satellite signal acquisition
WO2014171999A2 (en) * 2013-02-04 2014-10-23 Vanderbilt University Method and system for high-accuracy differential tracking of global positioning system (gps) receivers
US11372112B2 (en) * 2018-08-22 2022-06-28 Stmicroelectronics S.R.L. GNSS receiver apparatus with GNSS pseudo random noise delayed sequence generator

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100141520A1 (en) * 2007-08-21 2010-06-10 Giorgio Ghinamo Method for the acquisition of signals of a global navigation satellite system
JP2010060436A (ja) * 2008-09-04 2010-03-18 Furuno Electric Co Ltd Gnss受信装置および測位装置
JP2011202958A (ja) * 2010-03-24 2011-10-13 Seiko Epson Corp 信号捕捉方法、信号捕捉装置及び電子機器
JP2014216871A (ja) * 2013-04-26 2014-11-17 三菱電機株式会社 信号処理装置、信号処理方法及び信号処理プログラム

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JEON, SANGHOON ET AL.: "Bit Transition Cancellation Signal Acquisition Method for Modernized GPS and Galileo Signal", PROCEEDINGS OF THE 24TH INTERNATIONAL TECHNICAL MEETING OF THE SATELLITE DIVISION OF THE INSTITUTE OF NAVIGATION (ION GNSS 2011), 23 September 2011 (2011-09-23), pages 1028 - 1039, XP055554628 *
MASAYUKI HATA ET AL.: "A Noise Reduction Method in Calculating Correlation of FFT-GPS", IEICE TECHNICAL REPORT, vol. 106, no. 22, 14 April 2006 (2006-04-14), pages 55 - 60 *

Also Published As

Publication number Publication date
JPWO2018138830A1 (ja) 2019-02-07
US11169275B2 (en) 2021-11-09
US20200018862A1 (en) 2020-01-16
JP6253855B1 (ja) 2017-12-27

Similar Documents

Publication Publication Date Title
US20030227963A1 (en) Spread spectrum bit boundary correlation search acquisition system
US10042056B2 (en) Signal processing method for ultra-fast acquisition and tracking of severely attenuated spread spectrum signals with doppler frequency and apparatus thereof
JP4869022B2 (ja) 衛星信号追尾装置及びそれを備えた衛星信号受信機
US20120183023A1 (en) Implementations for wireless signal processing
JP2005031073A (ja) Gps相関ピーク信号の探索方法及びこれのためのシステム。
JP2009258107A (ja) Gnss信号を高速に取得するためのシステムおよび方法
JP2004340855A (ja) 衛星測位システム及び衛星測位方法
EP2829895B1 (en) Signal search method, signal search program, signal search device, global navigation satellite system (gnss) signal receiver, and information terminal
Leclere et al. Modified parallel code-phase search for acquisition in presence of sign transition
US9124356B2 (en) Reconfigurable correlator for a navigation system
JP6253855B1 (ja) 信号捕捉装置
KR100981458B1 (ko) 신호 수신 방법 및 수신기
WO2011001601A1 (ja) キャリア周波数同期検出回路及び相関演算器
KR101356239B1 (ko) Gnss 및 대역확산 신호의 신속한 신호 획득을 위한 압축 감지 방법 및 그 장치
US9612338B2 (en) Method to improve satellite signal detection
JP2008510426A (ja) 共通の復調テンプレートを使用する信号収集のための装置、方法およびコンピュータプログラム製品
CN113671547B (zh) 一种改进的高动态捕获方法、装置、设备及存储介质
US20120274512A1 (en) Signal Processing Method, Device and System
JP2006217601A (ja) 地理的局所化システムの位置決め信号を捕捉するための方法、地理的局所化システム用受信機及び上記の方法を実行するためのプログラム命令を含んでいるコンピュータデータ媒体
CN112764063A (zh) 一种实现捕获处理的方法及接收机
WO2013140910A1 (ja) 信号サーチ方法、信号サーチプログラム、信号サーチ装置、gnss信号受信装置、および情報機器端末
JP2014070960A (ja) 衛星信号受信装置
Iswariya et al. FFT based acquisition techniques of GPS L2C signals
CN101206256B (zh) 一种gps接收机信号非相干积分的数据处理方法
CN116347329B (zh) 定位信号捕获方法及装置、计算设备及计算机存储介质

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2017532196

Country of ref document: JP

Kind code of ref document: A

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

Ref document number: 17893958

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: 17893958

Country of ref document: EP

Kind code of ref document: A1