JP2011137802A - Received signal integrating method and receiving device - Google Patents

Received signal integrating method and receiving device Download PDF

Info

Publication number
JP2011137802A
JP2011137802A JP2010220172A JP2010220172A JP2011137802A JP 2011137802 A JP2011137802 A JP 2011137802A JP 2010220172 A JP2010220172 A JP 2010220172A JP 2010220172 A JP2010220172 A JP 2010220172A JP 2011137802 A JP2011137802 A JP 2011137802A
Authority
JP
Japan
Prior art keywords
period
received signal
signal
code
received
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2010220172A
Other languages
Japanese (ja)
Other versions
JP2011137802A5 (en
Inventor
Yoshitaka Yamagata
整功 山形
Norio Teruuchi
則生 照内
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Seiko Epson Corp
Original Assignee
Seiko Epson Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Seiko Epson Corp filed Critical Seiko Epson Corp
Priority to JP2010220172A priority Critical patent/JP2011137802A/en
Publication of JP2011137802A publication Critical patent/JP2011137802A/en
Publication of JP2011137802A5 publication Critical patent/JP2011137802A5/en
Pending legal-status Critical Current

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To provide an appropriate received signal integrating method in consideration of periodic deviation. <P>SOLUTION: A received signal of GPS satellite signal is sampled for each one clock as sampling time interval to acquire a plurality pieces of sampling data. Then, using sampling data sets with different code periods when the received signal is time-divided by an assumed period of CA code as diffusion sign of the GPS satellite signal, a periodic deviation coefficient representing a periodic deviation between a true period of CA code and the assumed period is calculated for each code period. Then the received signals are integrated using the periodic deviation coefficient. <P>COPYRIGHT: (C)2011,JPO&INPIT

Description

本発明は、受信信号積算方法及び受信装置に関する。   The present invention relates to a received signal integration method and a receiving apparatus.

測位用信号を利用した測位システムとしては、GPS(Global Positioning System)が広く知られており、携帯型電話機やカーナビゲーション装置等に内蔵されたGPS受信装置に利用されている。GPS受信装置は、複数のGPS衛星の位置や各GPS衛星から受信装置までの擬似距離等の情報に基づいて受信装置の位置を示す3次元の座標値と時計誤差とを求める位置算出処理を行う。   A GPS (Global Positioning System) is widely known as a positioning system using positioning signals, and is used in a GPS receiver built in a mobile phone or a car navigation device. The GPS receiver performs position calculation processing for obtaining a three-dimensional coordinate value indicating the position of the receiver and a clock error based on information such as the positions of a plurality of GPS satellites and pseudo distances from each GPS satellite to the receiver. .

GPS衛星信号は、スペクトラム拡散変調方式として古くから知られるCDMA(CodeDivision Multiple Access)方式で拡散変調された通信信号の一種である。受信信号からGPS衛星信号を捕捉する際には、受信信号と、GPS衛星信号の拡散符号であるCA(Coarse and Acquisition)コードのレプリカ信号との相関処理を、周波数及びコード位相を変化させながら行って(いわゆる周波数方向及び位相方向の相関演算。周波数サーチや位相サーチとも呼ばれる。)決定する手法が一般的である(例えば特許文献1)。   A GPS satellite signal is a type of communication signal that is spread-modulated by a CDMA (Code Division Multiple Access) method that has long been known as a spread spectrum modulation method. When capturing a GPS satellite signal from a received signal, correlation processing between the received signal and a CA (Coarse and Acquisition) code replica signal, which is a spread code of the GPS satellite signal, is performed while changing the frequency and code phase. (So-called correlation calculation in frequency direction and phase direction, also called frequency search or phase search) is generally determined (for example, Patent Document 1).

特開2007−256111号公報JP 2007-256111 A

GPS衛星信号の受信信号が弱電界となる環境(例えばインドア環境、以下「弱電界環境」と称す。)では受信信号が微弱であるため、相関処理を行うことで得られる相関値に差が現れず、相関値のピークが埋もれてしまう場合がある。そこで、弱電界環境等においては、受信信号を所定期間に亘って積算し、積算した信号に対する相関処理を行うことで、相関値のピークの判別を容易にする手法が用いられる。   Since the received signal is weak in an environment where the received signal of the GPS satellite signal is a weak electric field (for example, indoor environment, hereinafter referred to as “weak electric field environment”), a difference appears in the correlation value obtained by performing correlation processing. In some cases, the correlation value peak may be buried. Therefore, in a weak electric field environment or the like, a technique is used that makes it easy to determine the peak of the correlation value by integrating received signals over a predetermined period and performing correlation processing on the integrated signals.

ところで、GPS衛星信号の搬送波周波数は1.57542[GHz]である。また、GPS衛星信号の拡散符号であるCAコードは、コード長1023チップを1PNフレームとする繰返し周期1msの擬似ランダム雑音符号であり、チップレートは1.023[MHz]である。従って、理論上は、1チップ当たりの搬送波の周期数は1540周期であり、CAコードの1コード周期当たりでは1540×1023=1,575,420周期である。   By the way, the carrier wave frequency of the GPS satellite signal is 1.57542 [GHz]. The CA code, which is a spreading code for GPS satellite signals, is a pseudo-random noise code with a repetition period of 1 ms using a code length of 1023 chips as a 1PN frame, and the chip rate is 1.023 [MHz]. Therefore, theoretically, the number of carrier wave cycles per chip is 1540 cycles, and 1540 × 1023 = 1,575,420 cycles per CA code cycle.

しかし、GPS衛星信号を実際に受信した際の受信周波数には、いわゆるドップラー周波数や、ローカルクロックの誤差(時計誤差)に起因する周波数誤差が含まれる。これらの周波数誤差の存在により、本来繰り返し周期1msであるはずのCAコード1周期分を受信機側で推測したコード周期(以下、この周期を「想定周期」と称する。)は、真の周期からずれることとなる。搬送波が1,575,420周期する毎にCAコードは1コード周期となり、この周期が「真の周期」となる。しかし、受信機は搬送波の周期をカウントしてコード周期を決定しているわけではない。コード周期は、搬送波の周期をカウントせずに、想定して決定している。より詳細には、コード周期を直接想定するのではなく、周波数を決定することで、等価的にコード周期を想定している。   However, the reception frequency when the GPS satellite signal is actually received includes a so-called Doppler frequency and a frequency error caused by a local clock error (clock error). Due to the presence of these frequency errors, a code cycle (hereinafter, this cycle is referred to as an “assumed cycle”) in which one CA code cycle, which should originally be a repetition cycle of 1 ms, is estimated from the true cycle. It will shift. Each time the carrier wave is 1,575,420 cycles, the CA code becomes one code cycle, and this cycle becomes the “true cycle”. However, the receiver does not determine the code period by counting the period of the carrier wave. The code period is determined on the assumption without counting the period of the carrier wave. More specifically, instead of directly assuming the code period, the code period is equivalently assumed by determining the frequency.

このため、想定周期で受信信号を時分割した場合、搬送波が1,575,420周期ちょうどで時分割されるとは限らず、いくらかの位相のずれが生じる。すなわち、コード周期を想定して時分割しているために、あるコード周期の開始時点における搬送波の位相と、次のコード周期の開始時点における搬送波の位相との間にはずれが生じる(以下、このずれのことを「周期ずれ」と称する。)。周期ずれは、CAコードの真の周期と想定周期とのずれと等価である。上記のように、コード周期の想定とは、受信周波数を求めることとも言える。受信周波数に周波数誤差が含まれることにより、周期ずれが生じるのである。   For this reason, when the received signal is time-divided at the assumed period, the carrier wave is not always time-divided at the 1,575,420 period, and some phase shift occurs. In other words, since time division is performed assuming a code period, there is a difference between the phase of the carrier at the start of a certain code period and the phase of the carrier at the start of the next code period (hereinafter referred to as this The deviation is referred to as “periodic deviation”.) The period deviation is equivalent to a deviation between the true period of the CA code and the assumed period. As described above, it can be said that the assumption of the code period is to obtain the reception frequency. The frequency deviation is included in the reception frequency, resulting in a period shift.

弱電界環境において、相関値のピークの判別を容易にするために受信信号を積算しようとしても、周期ずれが生じた状態で受信信号を積算した場合には、受信信号の振幅が却って小さくなる場合がある。周期を正確に把握できずに積算していくと、積算する信号の位相がずれていき、振幅の正負が入れ替わった信号を積算してしまう可能性があるためである。また、周期ずれが生じた状態で受信信号を積算し、相関処理を行った場合、判別された相関値のピークが正しい結果でない可能性もある。   In a weak electric field environment, even if you try to integrate received signals to make it easier to determine the correlation value peak, if the received signals are integrated with a period shift, the amplitude of the received signal will be smaller There is. This is because if the period is accumulated without accurately grasping the period, the phase of the signal to be accumulated is shifted, and a signal in which the amplitude is changed may be accumulated. In addition, when the received signals are integrated and correlation processing is performed in a state where a period shift has occurred, the peak of the determined correlation value may not be a correct result.

本発明は、上述した課題に鑑みて為されたものであり、その目的とするところは、周期ずれを考慮した適切な受信信号積算手法を提案することにある。   The present invention has been made in view of the above-described problems, and an object of the present invention is to propose an appropriate received signal integration method in consideration of a period shift.

以上の課題を解決するための第1の形態は、衛星信号を受信する際に前記衛星信号の拡散符号の符号周期時間を推測した想定周期で前記衛星信号の受信信号を時分割した場合に、前記受信信号のうちの第1受信信号部分と前記第1受信信号部分とは想定周期が異なる第2受信信号部分とを用いて、前記拡散符号の真の周期と前記想定周期との周期ずれを表す周期ずれ係数を算出することと、前記受信信号を前記周期ずれ係数を用いて積算することと、を含む受信信号積算方法である。   The first mode for solving the above problems is when the received signal of the satellite signal is time-divided with an assumed period in which the code period time of the spreading code of the satellite signal is estimated when the satellite signal is received. The first received signal portion and the first received signal portion of the received signal have a second received signal portion having a different assumed cycle, and the period deviation between the true cycle of the spread code and the assumed cycle is changed. A received signal integration method comprising: calculating a periodic shift coefficient to be expressed; and integrating the received signal using the periodic shift coefficient.

この第1の形態によれば、衛星信号の拡散符号の想定周期で衛星信号の受信信号を時分割した場合に、受信信号のうちの第1受信信号部分と、想定周期が異なる第2受信信号部分とを用いて、拡散符号の真の周期と想定周期との周期ずれを表す周期ずれ係数を算出する。そして、受信信号を周期ずれ係数を用いて積算する。   According to the first aspect, when the satellite signal reception signal is time-divided in the assumed cycle of the satellite signal spreading code, the second received signal is different from the first received signal portion in the received signal. Is used to calculate a period deviation coefficient representing a period deviation between the true period of the spread code and the assumed period. Then, the received signals are integrated using a period shift coefficient.

周期ずれは、衛星信号の拡散符号の真の周期と想定周期との差である。詳細は後述するが、本願発明者は、拡散符号の想定周期で衛星信号の受信信号を時分割することで得られる受信信号部分のうち、想定周期が異なる受信信号部分を用いることで、周期ずれを表す係数(指数)が算出可能であることを発見した。この周期ずれ係数を用いて受信信号を積算することで、周期ずれを考慮した適切な積算受信信号を得ることができる。   The period shift is a difference between the true period of the spread code of the satellite signal and the assumed period. Although details will be described later, the inventor of the present application uses a received signal portion having a different assumed cycle among received signal portions obtained by time-division of a received signal of a satellite signal with an assumed cycle of a spread code, thereby causing a period shift. It was discovered that a coefficient (index) representing can be calculated. By integrating the reception signals using this period deviation coefficient, it is possible to obtain an appropriate integrated reception signal in consideration of the period deviation.

また、第2の形態として、第1の形態の受信信号積算方法であって、前記受信信号を積算することは、前記第1受信信号部分及び前記第2受信信号部分と前記周期ずれ係数とを積和することで前記受信信号を積算した信号を生成することである、受信信号積算方法を構成してもよい。   Further, as a second mode, in the received signal integration method according to the first mode, the integration of the received signals includes the first received signal portion, the second received signal portion, and the period deviation coefficient. You may comprise the received signal integration method which is producing | generating the signal which integrated | accumulated the said received signal by multiply-adding.

この第2の形態によれば、第1受信信号部分及び第2受信信号部分と周期ずれ係数とを積和することで、受信信号を積算した信号を生成する。受信信号部分に周期ずれ係数を乗算することで、周期ずれに起因する誤差成分を除去することができる。そのため、受信信号を積算した信号は、周期ずれに起因する信号劣化の生じていない高品質な信号となる。   According to the second aspect, the sum of the received signals is generated by multiplying and summing the first received signal portion and the second received signal portion and the period shift coefficient. By multiplying the received signal portion by a period deviation coefficient, an error component due to the period deviation can be removed. For this reason, the signal obtained by integrating the received signals is a high-quality signal in which no signal deterioration due to the period shift occurs.

また、第3の形態として、第1又は第2の形態の受信信号積算方法であって、前記周期ずれ係数を算出することは、前記受信信号のうち、前記想定周期がn個分(nは自然数)異なる前記第2受信信号部分それぞれについての前記周期ずれ係数を算出することを含み、前記受信信号を積算することは、前記想定周期がn個分異なる前記第2受信信号部分それぞれを積算する際に、対応する前記周期ずれ係数を用いて積算することを含む受信信号積算方法を構成してもよい。   Further, as a third mode, in the received signal integration method according to the first or second mode, the calculation of the period deviation coefficient may include n assumption cycles of the received signal (n is (Natural number) calculating the period deviation coefficient for each of the different second received signal parts, and integrating the received signals integrates the second received signal parts having n different assumed periods. In this case, a reception signal integration method including integration using the corresponding period deviation coefficient may be configured.

この第3の形態によれば、拡散符号の想定周期がn個分異なる第2受信信号部分それぞれについての周期ずれ係数を算出する。そして、想定周期がn個分異なる第2受信信号部分それぞれを積算する際に、対応する周期ずれ係数を用いて積算する。複数の受信信号部分を周期ずれ係数を用いて積算することで、より品質の高い積算信号を得ることができる。   According to the third embodiment, the period deviation coefficient is calculated for each of the second received signal portions having n different spread code assumption periods. Then, when integrating each of the second received signal portions having different assumed periods by n, the corresponding period deviation coefficients are used for integration. By integrating a plurality of received signal portions using a period shift coefficient, an integrated signal with higher quality can be obtained.

また、第4の形態として、第1〜第3の何れかの形態の受信信号積算方法であって、前記周期ずれ係数を算出することは、前記第1及び前記第2受信信号部分において、前記想定周期内で同一タイミングとなる信号部分を用いて、前記周期ずれ係数を算出することを含む、受信信号積算方法を構成してもよい。   Further, as a fourth mode, in the received signal integration method according to any one of the first to third modes, the calculation of the period shift coefficient may be performed in the first and second received signal portions. You may comprise the received signal integration method including calculating the said period shift coefficient using the signal part which becomes the same timing within an assumption period.

この第4の形態によれば、第1及び第2受信信号部分において、想定周期内で同一タイミングとなる信号部分を用いて、周期ずれ係数を算出する。想定周期内でタイミングが一致する信号部分を用いることで、周期ずれ係数を適切に算出することができる。   According to this 4th form, a period shift coefficient is calculated using the signal part which becomes the same timing within an assumption period in a 1st and 2nd received signal part. By using a signal portion having the same timing within the assumed period, the period deviation coefficient can be calculated appropriately.

また、第5の形態として、第1〜第4の何れかの形態の受信信号積算方法であって、前記周期ずれ係数を算出することは、前記第1受信信号部分と、前記第2受信信号部分の複素共役とを乗算することで前記周期ずれ係数を算出することを含む、受信信号積算方法を構成してもよい。   Further, as a fifth mode, in the received signal integration method according to any one of the first to fourth modes, calculating the period deviation coefficient includes the first received signal portion and the second received signal. You may comprise the received signal integration method including calculating the said period shift coefficient by multiplying with the complex conjugate of a part.

この第5の形態によれば、第1受信信号部分と、第2受信信号部分の複素共役とを乗算するといった簡易な演算により、周期ずれ係数を算出することができる。   According to the fifth embodiment, the period deviation coefficient can be calculated by a simple calculation such as multiplying the first received signal portion by the complex conjugate of the second received signal portion.

また、第6の形態として、第4の形態の受信信号積算方法であって、前記周期ずれ係数を算出することは、前記第1受信信号部分のうちの前記想定周期内で異なるタイミングとなる信号部分と、第n番目(nは自然数)の前記第2受信信号部分のうちの対応する同一タイミングとなる信号部分の複素共役とを乗算することと、第n番目の前記第2受信信号部分に関して得られた前記各タイミングの前記乗算結果の平均を、前記第1受信信号部分と第n番目の前記第2受信信号部分との周期ずれ係数とすることと、を含む、受信信号積算方法を構成してもよい。   Further, as a sixth aspect, in the received signal integration method according to the fourth aspect, the calculation of the period deviation coefficient is a signal having a different timing within the assumed period of the first received signal portion. Multiplying the portion by the complex conjugate of the corresponding signal portion of the nth (n is a natural number) second received signal portion and the nth second received signal portion. The average of the obtained multiplication results of the respective timings is used as a period shift coefficient between the first received signal portion and the nth second received signal portion, and a received signal integrating method is configured. May be.

この第6の形態によれば、第1受信信号部分のうちの想定周期内で異なるタイミングとなる信号部分と、第n番目の第2受信信号部分のうちの対応する同一タイミングとなる信号部分の複素共役とを乗算する。そして、第n番目の第2受信信号部分に関して得られた各タイミングの乗算結果の平均を、第1受信信号部分と第n番目の第2受信信号部分との周期ずれ係数とする。このような処理を行うことで、より正確な周期ずれ係数を算出することができる。   According to the sixth aspect, the signal portion having the different timing within the assumed period in the first received signal portion and the corresponding signal portion in the nth second received signal portion having the same timing. Multiply by complex conjugate. Then, the average of the multiplication results at the respective timings obtained for the nth second received signal portion is taken as the period shift coefficient between the first received signal portion and the nth second received signal portion. By performing such processing, a more accurate periodic deviation coefficient can be calculated.

また、第7の形態として、第1〜第6の何れかの形態の受信信号積算方法であって、前記受信信号を積算することは、前記衛星信号の搬送波が除去されていない状態の受信信号を前記周期ずれ係数を用いて積算することである、受信信号積算方法を構成してもよい。   Further, as a seventh aspect, the reception signal integration method according to any one of the first to sixth aspects, wherein the integration of the reception signals is a reception signal in a state where a carrier wave of the satellite signal is not removed. May be configured to integrate the received signal using the period deviation coefficient.

この第7の形態によれば、衛星信号の搬送波が除去されていない状態の受信信号を周期ずれ係数を用いて積算する。受信信号から衛星信号の搬送波を除去する必要がないため、衛星信号の受信回路の簡素化が実現される。   According to the seventh aspect, the received signals in a state where the carrier wave of the satellite signal is not removed are integrated using the period shift coefficient. Since it is not necessary to remove the carrier wave of the satellite signal from the received signal, simplification of the satellite signal receiving circuit is realized.

また、第8の形態として、衛星信号を受信する際に前記衛星信号の拡散符号の符号周期時間を推測した想定周期で前記衛星信号の受信信号を時分割した場合に、前記受信信号のうちの第1受信信号部分と前記第1受信信号部分とは前記想定周期が異なる第2受信信号部分とを用いて、前記拡散符号の真の周期と前記想定周期との周期ずれを表す周期ずれ係数を算出する周期ずれ係数算出部と、前記受信信号を前記周期ずれ係数を用いて積算する受信信号積算部と、前記受信信号積算部により積算された信号に対する相関処理を行う相関処理部と、前記相関処理の結果に基づいて前記衛星信号を捕捉する捕捉部と、を備えた受信装置を構成してもよい。   Further, as an eighth aspect, when the satellite signal reception signal is time-divided with an assumed period in which the code period time of the spreading code of the satellite signal is estimated when the satellite signal is received, The first received signal portion and the first received signal portion use a second received signal portion having a different assumed period, and a period deviation coefficient representing a period deviation between the true period of the spread code and the assumed period A periodic deviation coefficient calculating unit to calculate, a received signal integrating unit that integrates the received signal using the periodic deviation coefficient, a correlation processing unit that performs correlation processing on the signal integrated by the received signal integrating unit, and the correlation You may comprise the receiver provided with the acquisition part which acquires the said satellite signal based on the result of a process.

この第8の形態によれば、周期ずれ係数算出部により、拡散符号の真の周期と想定周期との周期ずれを表す周期ずれ係数が算出される。そして、受信信号積算部により周期ずれ係数を用いて受信信号が積算され、積算された信号に対する相関処理が、相関処理部により行われる。そして、捕捉部により、相関処理の結果に基づいて衛星信号が捕捉される。かかる構成により、第1の形態と同様の効果が発揮されるとともに、周期ずれを考慮した積算受信信号に対する相関処理を行うことで、衛星信号を捕捉するための相関値を適切に求めることができる。   According to the eighth embodiment, the period deviation coefficient calculating unit calculates a period deviation coefficient representing the period deviation between the true period of the spread code and the assumed period. Then, the received signal is integrated by the received signal integrating unit using the period shift coefficient, and the correlation processing for the integrated signal is performed by the correlation processing unit. Then, the capturing unit captures the satellite signal based on the correlation processing result. With this configuration, the same effect as that of the first embodiment is exhibited, and the correlation value for capturing the satellite signal can be appropriately obtained by performing the correlation process on the accumulated reception signal in consideration of the period shift. .

受信信号積算の原理の説明図。Explanatory drawing of the principle of reception signal integration. 受信信号積算処理の流れを示すフローチャート。The flowchart which shows the flow of a received signal integration process. 周期ずれ係数算出の原理の説明図。Explanatory drawing of the principle of period shift coefficient calculation. 受信信号積算の原理の説明図。Explanatory drawing of the principle of reception signal integration. 携帯型電話機の機能構成を示すブロック図。The block diagram which shows the function structure of a portable telephone. ベースバンド処理回路部の回路構成の一例を示す図。The figure which shows an example of the circuit structure of a baseband process circuit part. 周期ずれ係数算出処理の流れを示すフローチャート。The flowchart which shows the flow of a period shift coefficient calculation process. 受信信号積算処理の流れを示すフローチャート。The flowchart which shows the flow of a received signal integration process. ベースバンド処理の流れを示すフローチャート。The flowchart which shows the flow of a baseband process. コード位相検出処理の流れを示すフローチャート。The flowchart which shows the flow of a code phase detection process.

以下、図面を参照して、本発明の好適な実施形態について説明する。以下では、GPS(Global Positioning System)衛星から発信されているGPS衛星信号を受信・捕捉するGPS受信装置に本発明を適用した場合について説明する。尚、本発明を適用可能な実施形態が以下説明する実施形態に限定されるわけでないことは勿論である。   Hereinafter, preferred embodiments of the present invention will be described with reference to the drawings. Below, the case where this invention is applied to the GPS receiver which receives and captures the GPS satellite signal transmitted from the GPS (Global Positioning System) satellite is demonstrated. Of course, embodiments to which the present invention is applicable are not limited to the embodiments described below.

1.原理
先ず、本実施形態における受信信号積算の原理について説明する。
GPS衛星は、測位用衛星の一種であり、6つの地球周回軌道面それぞれに4機以上ずつ配置され、原則、地球上のどこからでも常時4機以上の衛星が幾何学的配置のもとで観測できるように運用されている。
1. Principle First, the principle of received signal integration in this embodiment will be described.
GPS satellites are a kind of positioning satellites, and four or more satellites are arranged on each of the six orbiting surfaces of the earth. In principle, four or more satellites are always observed from anywhere on the earth in a geometrical arrangement. It is operated as possible.

GPS衛星は、アルマナックやエフェメリス等の航法メッセージを、測位用信号の一種であるGPS衛星信号に含めて発信している。GPS衛星信号は、拡散符号の一種であるCA(Coarse and Acquisition)コードによって、スペクトラム拡散方式として知られるCDMA(Code Division Multiple Access)方式によって変調された1.57542[GHz]の通信信号である。CAコードは、コード長1023チップを1PNフレームとする繰返し周期1msの擬似ランダム雑音符号であり、GPS衛星毎に異なるものである。   GPS satellites transmit navigation messages such as almanac and ephemeris in GPS satellite signals that are a type of positioning signals. The GPS satellite signal is a 1.57542 [GHz] communication signal modulated by a CDMA (Code Division Multiple Access) system known as a spread spectrum system by a CA (Coarse and Acquisition) code which is a kind of spreading code. The CA code is a pseudo-random noise code having a repetition period of 1 ms with a code length of 1023 chips as 1 PN frame, and is different for each GPS satellite.

GPS衛星がGPS衛星信号を発信する際の周波数(搬送波周波数)は、1.57542[GHz]と予め規定されているが、GPS衛星やGPS受信装置の移動により生ずるドップラーの影響等により、GPS受信装置がGPS衛星信号を受信する際の周波数は、必ずしも搬送波周波数とは一致しない。そのため、従来のGPS受信装置は、受信信号の中からGPS衛星信号を捕捉するための周波数方向の相関演算である周波数サーチ及び位相方向の相関演算である位相サーチを行ってGPS衛星信号を捕捉する。   The frequency (carrier frequency) at which a GPS satellite transmits a GPS satellite signal is specified in advance as 1.57542 [GHz]. However, GPS reception is possible due to the influence of Doppler caused by the movement of the GPS satellite and the GPS receiver. The frequency at which the device receives GPS satellite signals does not necessarily match the carrier frequency. Therefore, a conventional GPS receiver captures a GPS satellite signal by performing a frequency search that is a correlation calculation in the frequency direction and a phase search that is a correlation calculation in the phase direction for capturing the GPS satellite signal from the received signal. .

周波数サーチ及び位相サーチでは、受信信号とGPS衛星信号の拡散符号のレプリカであるレプリカコードの発生信号との相関処理を行って相関値を取得し、相関値が最大となった周波数及び位相を検出する。GPS衛星信号の受信信号が強電界の信号となる環境(例えば屋外環境、以下「強電界環境」と称す。)では、受信信号に対して相関処理を行った場合に相関値に明確な差が現れやすいため、相関値のピーク検出は比較的容易である。   In frequency search and phase search, the correlation value is acquired by performing correlation processing between the received signal and the generated signal of the replica code that is a replica of the spread code of the GPS satellite signal, and the frequency and phase at which the correlation value is maximized are detected. To do. In an environment where the received signal of the GPS satellite signal is a strong electric field signal (for example, an outdoor environment, hereinafter referred to as “strong electric field environment”), there is a clear difference in the correlation value when correlation processing is performed on the received signal. Since it appears easily, the peak detection of the correlation value is relatively easy.

しかし、GPS衛星信号の受信信号が弱電界の信号となる環境(例えば屋内環境、以下「弱電界環境」と称す。)では、受信信号に対して相関処理を行った場合に相関値に明確な差が現れず、相関値のピークの判別が容易ではない場合が多い。そのため、弱電界環境等では、受信信号を所定期間に亘って積算し、積算した受信信号(以下、「積算受信信号」と称す。)に対して相関処理を行って相関値を取得する手法が用いられる。   However, in an environment where the received signal of the GPS satellite signal is a weak electric field signal (for example, indoor environment, hereinafter referred to as “weak electric field environment”), the correlation value is clear when the correlation processing is performed on the received signal. In many cases, the difference does not appear and it is not easy to determine the peak of the correlation value. Therefore, in a weak electric field environment or the like, there is a technique for accumulating received signals over a predetermined period and performing correlation processing on the accumulated received signals (hereinafter referred to as “integrated received signals”) to obtain a correlation value. Used.

図1は、従来における受信信号積算の概念の説明図である。GPS衛星信号の拡散符号であるCAコードは周期性を有している。具体的には、コード長1023チップを1PNフレームとし、繰返し周期を1msとして、GPS衛星から繰り返し送信されている。そのため、GPS衛星信号の受信信号をCAコードの周期時間間隔で足し上げていけば、振幅(パワー)が大きな積算信号が得られるはずである。   FIG. 1 is an explanatory diagram of the concept of conventional received signal integration. The CA code, which is a spreading code for GPS satellite signals, has periodicity. Specifically, it is repeatedly transmitted from a GPS satellite with a code length of 1023 chips as a 1PN frame and a repetition period of 1 ms. Therefore, if the received signal of the GPS satellite signal is added at the CA code cycle time interval, an integrated signal with a large amplitude (power) should be obtained.

具体的に説明すると、GPS衛星信号の受信信号を所定のサンプリング時間間隔でサンプリングすることで、受信信号のサンプリングデータを取得する。本実施形態では、受信信号のサンプリング単位を「クロック」と称し、1クロックの経過時間を「T」で表す。サンプリングの時間間隔(すなわち時間T)は、CAコードの1チップ分の時間間隔としてもよいし、1チップをより細かく細分化した時間間隔としてもよい。   More specifically, sampling data of the received signal is acquired by sampling the received signal of the GPS satellite signal at a predetermined sampling time interval. In this embodiment, the sampling unit of the received signal is referred to as “clock”, and the elapsed time of one clock is represented by “T”. The sampling time interval (that is, time T) may be a time interval for one chip of the CA code, or may be a time interval obtained by finely subdividing one chip.

ここで、時刻「t」における受信信号「r(t)」は、次式(1)のように表現することができる。

Figure 2011137802
式(1)において、「I(t)」と「Q(t)」はそれぞれ受信信号「r(t)」のIQ成分を示している。すなわち「I(t)」は受信信号「r(t)」の同相成分(実部)を示し、「Q(t)」は「r(t)」の直交成分(虚部)を示す。「CA(t)」はGPS衛星信号のCAコードを示しており、「+1」と「−1」の何れかの値である。また、「exp(iωt)」は、GPS衛星信号の搬送波を表す項である。 Here, the received signal “r (t)” at time “t” can be expressed as the following equation (1).
Figure 2011137802
In Expression (1), “I (t)” and “Q (t)” indicate IQ components of the received signal “r (t)”, respectively. That is, “I (t)” represents the in-phase component (real part) of the received signal “r (t)”, and “Q (t)” represents the quadrature component (imaginary part) of “r (t)”. “CA (t)” indicates a CA code of the GPS satellite signal, and is a value of “+1” or “−1”. “Exp (iωt)” is a term representing a carrier wave of a GPS satellite signal.

式(1)において「ω」は受信信号の周波数であり、次式(2)で表される。

Figure 2011137802
但し、「ωc」はGPS衛星信号の搬送波周波数であり、「ωd」は周波数誤差(例えば、ドップラー周波数やローカルクロックの誤差(時計誤差))である。 In equation (1), “ω” is the frequency of the received signal, and is represented by the following equation (2).
Figure 2011137802
However, “ω c ” is the carrier frequency of the GPS satellite signal, and “ω d ” is a frequency error (for example, an error of the Doppler frequency or a local clock (clock error)).

受信信号「r(t)」を1クロック毎にサンプリングしていき、1つのCAコードの周期につき、例えばM+1(m=0,1,2,・・・,M)個のサンプリングデータを取得する。そして、例えばN+1(n=0,1,2,・・・,N)周期分サンプリングを行うと、全部で(M+1)×(N+1)個のサンプリングデータが得られる。   The received signal “r (t)” is sampled every clock and, for example, M + 1 (m = 0, 1, 2,..., M) sampling data is acquired for each CA code period. . For example, if sampling is performed for N + 1 (n = 0, 1, 2,..., N) periods, a total of (M + 1) × (N + 1) sampling data is obtained.

注意しなければならないのは、ここで言うCAコードの周期は、GPS受信装置がCAコードの周期時間を推測した周期(想定周期)であって、CAコードの真の周期とは異なることである。前述したように、GPS受信装置は、GPS衛星信号の周波数を決定することで、等価的にCAコードの周期を決定している。しかし、ドップラーの影響によりGPS受信装置がGPS衛星信号を受信する際の周波数は搬送波周波数とは完全に一致しない。そのため、見かけ上のCAコードの周期はCAコードの真の周期と乖離する場合がある。   It should be noted that the CA code period mentioned here is a period (assumed period) in which the GPS receiver estimates the period time of the CA code, and is different from the true period of the CA code. . As described above, the GPS receiver determines the frequency of the GPS satellite signal to determine the CA code period equivalently. However, the frequency at which the GPS receiver receives a GPS satellite signal does not completely match the carrier frequency due to the influence of Doppler. For this reason, the apparent CA code period may deviate from the true period of the CA code.

また、GPS受信装置内のローカルクロックの誤差(時計誤差)により、GPS受信装置が内部で計測した1msの周期が正確ではなく、想定周期と真の周期とが乖離する場合もある。そのため、通常は、真の周期と想定周期との間には「周期ずれ」が生ずる。この周期ずれは、前述したように、あるコード周期の開始時点における搬送波の位相と、別のコード周期の開始時点における搬送波の位相とのずれと等価である。   In addition, due to an error (clock error) of the local clock in the GPS receiver, the 1 ms period internally measured by the GPS receiver may not be accurate, and the assumed period may deviate from the true period. Therefore, normally, a “period shift” occurs between the true period and the assumed period. As described above, this period deviation is equivalent to a deviation between the phase of the carrier at the start of a certain code period and the phase of the carrier at the start of another code period.

本実施形態では、CAコードの想定周期「TCA」でGPS衛星信号の受信信号「r(t)」を時分割した場合における各想定周期の番号「n」を「コード周期番号」と称し、各コード周期における受信信号のサンプリングの番号「m」を「サンプリング番号」と称する。但し、「n」と「m」は共に自然数である。 In the present embodiment, the number “n” of each assumed period when the received signal “r (t)” of the GPS satellite signal is time-divided by the assumed period “T CA ” of the CA code is referred to as “code period number”. The sampling number “m” of the received signal in each code period is referred to as “sampling number”. However, “n” and “m” are both natural numbers.

そして、コード周期番号「n」におけるサンプリング番号「m」のサンプリングデータを「rn,m(t)」と表記する。すなわち、コード周期番号及びサンプリング番号をこの順で下付きの添え字で表し、対応する時刻を括弧書きで表す。また、対応するCAコードを「CAn,m(t)」と表記する。 The sampling data of the sampling number “m” in the code cycle number “n” is expressed as “r n, m (t)”. That is, the code cycle number and the sampling number are represented by subscripts in this order, and the corresponding time is represented by parentheses. The corresponding CA code is expressed as “CA n, m (t)”.

図1では、第0コード周期(n=0)については{r0,0(t),r0,1(t+T),r0,2(t+2T),・・・,r0,M(t+MT)}の「M+1個」のサンプリングデータが得られる。同様に、第1コード周期(n=1)については{r1,0(t+TCA),r1,1(t+T+TCA),r1,2(t+2T+TCA),・・・,r1,M(t+MT+TCA)}が,・・・,第Nコード周期(n=N)については{rN,0(t+NTCA),rN,1(t+T+NTCA),rN,2(t+2T+NTCA),・・・,rN,M(t+MT+NTCA)}のサンプリングデータが得られる。すなわち、各コード周期「n」について「M+1個」のサンプリングデータが得られる。 In FIG. 1, for the 0th code period (n = 0), {r 0,0 (t), r 0,1 (t + T), r 0,2 (t + 2T),..., R 0, M (t + MT )} “M + 1” sampling data is obtained. Similarly, for the first code period (n = 1), {r 1,0 (t + T CA ), r 1,1 (t + T + T CA ), r 1,2 (t + 2T + T CA ),..., R 1, M (T + MT + T CA )}, for the Nth code period (n = N), {r N, 0 (t + NT CA ), r N, 1 (t + T + NT CA ), r N, 2 (t + 2T + NT CA ), .., R N, M (t + MT + NT CA )} sampling data is obtained. That is, “M + 1” sampling data is obtained for each code period “n”.

次いで、各サンプリング番号「m」それぞれについて、各コード周期(n=0〜N)のサンプリングデータを合算して、積算サンプリングデータ「Rm」を取得する。具体的には、次式(3)に従って積算サンプリングデータ「Rm(t)」を算出する。

Figure 2011137802
但し、CAコードの周期性から「CAn,m(t)=CAn,m(t+nTCA)」となる性質を利用している。 Next, for each sampling number “m”, the sampling data of each code period (n = 0 to N) is added to obtain integrated sampling data “R m ”. Specifically, the integrated sampling data “R m (t)” is calculated according to the following equation (3).
Figure 2011137802
However, the property of “CA n, m (t) = CA n, m (t + nT CA )” is used from the periodicity of the CA code.

例えば、サンプリング番号「m=1」に着目して説明する。図1に示すように、第0コード周期のサンプリングデータ「r0,1(t+T)」と、第1コード周期のサンプリングデータ「r1,1(t+T+TCA)」と、第2コード周期のサンプリングデータ「r2,1(t+T+2TCA)」と,・・・,第Nコード周期のサンプリングデータ「rN,1(t+T+NTCA)」とを合算して、積算サンプリングデータ「R1(t+T)」を算出する。 For example, the description will be given focusing on the sampling number “m = 1”. As shown in FIG. 1, sampling data “r 0,1 (t + T)” in the 0th code period, sampling data “r 1,1 (t + T + T CA )” in the first code period, and sampling in the second code period The data “r 2,1 (t + T + 2T CA )”,..., The sampling data “r N, 1 (t + T + NT CA )” of the Nth code period are added together, and the integrated sampling data “R 1 (t + T)” Is calculated.

ここで、式(3)を見ると、積算サンプリングデータ「Rm(t)」には、「Σexp(iωnTCA)」で表される項が含まれている。この項の絶対値は、次式(4)で表されるように、最大で「N+1」、最小で「0」となる項である。

Figure 2011137802
Here, looking at equation (3), the integrated sampling data “R m (t)” includes a term represented by “Σexp (iωnT CA )”. The absolute value of this term is a term that is “N + 1” at the maximum and “0” at the minimum, as represented by the following equation (4).
Figure 2011137802

従って、「Σexp(iωnTCA)」で表される項の値が「1」よりも小さければ、積算サンプリングデータ「Rm(t)」は、元のサンプリングデータ「rn,m(t)」と比べて値が小さくなる。 Therefore, if the value of the term represented by “Σexp (iωnT CA )” is smaller than “1”, the integrated sampling data “R m (t)” is converted into the original sampling data “r n, m (t)”. The value is smaller than

この意味を考察する。式(4)の「exp(iωnTCA)」は、第0コード周期の開始時点における搬送波の位相と、第nコード周期の開始時点における搬送波の位相とのずれ、すなわち第nコード周期における周期ずれを表していると考えられる。従って、各コード周期における周期ずれの大きさによっては、受信信号を積算した場合に積算した信号の振幅(パワー)が小さくなる場合があることを意味している。 Consider this meaning. “Exp (iωnT CA )” in Expression (4) is a shift between the phase of the carrier at the start of the 0th code period and the phase of the carrier at the start of the nth code period, that is, a period shift in the nth code period. It is thought that it represents. Therefore, depending on the magnitude of the period deviation in each code period, it means that the amplitude (power) of the accumulated signal may be reduced when the received signal is accumulated.

上記の計算式からわかるように、このような問題が生ずる原因は、搬送波「exp(iωt)」の存在によるものである。すなわち、受信信号「r(t)」から搬送波「exp(iωt)」を除去せずに受信信号を積算していくと、周期ずれの存在により、信号が強め合うどころか逆に弱め合う場合があるということである。   As can be seen from the above calculation formula, the cause of such a problem is due to the presence of the carrier wave “exp (iωt)”. That is, if the received signal is integrated without removing the carrier wave “exp (iωt)” from the received signal “r (t)”, the signal may be weakened rather than strengthened due to the presence of a period shift. That's what it means.

この問題を解決するため、本願発明者は、各コード周期における周期ずれを表す指数として「周期ずれ係数」と呼ぶ係数を算出し、この周期ずれ係数を用いて受信信号を積算する全く新しい受信信号の積算手法を考案した。   In order to solve this problem, the inventor of the present application calculates a coefficient called “period deviation coefficient” as an index representing the period deviation in each code period, and uses this period deviation coefficient to integrate the received signals. We have devised an integration method.

図2は、本実施形態における受信信号積算処理を示すフローチャートである。
最初に、受信信号のサンプリング処理及び蓄積処理を行う(ステップA1)。具体的には、受信信号を所定のサンプリング時間間隔(1クロック毎のタイミング)でサンプリングし、そのサンプリングデータ「rn,m(t)」を記憶部に蓄積する。
FIG. 2 is a flowchart showing the received signal integration process in this embodiment.
First, a sampling process and an accumulation process of the received signal are performed (step A1). Specifically, the received signal is sampled at a predetermined sampling time interval (timing for each clock), and the sampling data “r n, m (t)” is accumulated in the storage unit.

次いで、周期ずれ係数算出処理を行う(ステップA3)。周期ずれ係数算出処理では、各コード周期(n=0〜N)それぞれについて、異なるコード周期の受信信号部分を用いて、周期ずれを表す周期ずれ係数を算出する。具体的には、例えば第0コード周期(n=0)を基準コード周期とした場合に、他のコード周期(n=1〜N)の周期ずれ係数「Sn」を、次式(5)に従って算出する。

Figure 2011137802
但し、「CA0,m(t)×CAn,m(t+nTCA)=1」となる性質を利用している。また、上付きの添え字の「*」は複素共役を示している。 Next, a period deviation coefficient calculation process is performed (step A3). In the period deviation coefficient calculation process, for each code period (n = 0 to N), a period deviation coefficient representing the period deviation is calculated using the received signal portions having different code periods. Specifically, for example, when the 0th code period (n = 0) is set as the reference code period, the period deviation coefficient “S n ” of other code periods (n = 1 to N) is expressed by the following equation (5). Calculate according to
Figure 2011137802
However, the property of “CA 0, m (t) × CA n, m (t + nT CA ) = 1” is used. The superscript “*” indicates a complex conjugate.

図3は、周期ずれ係数算出の原理の説明図である。基準コード周期(第0コード周期)のサンプリングデータ(第1受信信号部分)と、他のコード周期(第nコード周期)のサンプリングデータ(第2受信信号部分)とを用いて周期ずれ係数を算出する。具体的には、図3に示すように、第0コード周期で異なるタイミングとなる信号部分、すなわち各サンプリング番号「m」(m=0〜M)のサンプリングデータ「r0,m」と、第nコード周期内で対応する同一タイミングとなる信号部分、すなわち対応するサンプリング番号「m」(m=0〜M)のサンプリングデータの複素共役「{rn,m」とを、サンプリング番号「m」を揃えてそれぞれ乗算する。 FIG. 3 is an explanatory diagram of the principle of calculating the period deviation coefficient. A period deviation coefficient is calculated using sampling data (first received signal portion) of a reference code cycle (0th code cycle) and sampling data (second received signal portion) of another code cycle (nth code cycle). To do. Specifically, as shown in FIG. 3, signal portions having different timings in the 0th code period, that is, sampling data “r 0, m ” of each sampling number “m” (m = 0 to M), A signal portion corresponding to the same timing within the n code period, that is, a complex conjugate “{r n, m } * ” of the sampling data corresponding to the sampling number “m” (m = 0 to M) is set to the sampling number “ m ”are aligned and multiplied.

そして、M+1個のサンプリングデータそれぞれについて得られた乗算結果「r0,m・{rn,m」の平均を、基準コード周期(第0コード周期)と第nコード周期との周期ずれ係数「Sn」とする。すなわち、受信信号のうち、各コード周期における受信信号部分としてのサンプリングデータの組のうち、コード周期内で同一タイミングとなるサンプリングデータを用いて、各コード周期それぞれについて周期ずれ係数を算出する。 Then, the average of the multiplication results “r 0, m · {r n, m } * ” obtained for each of the M + 1 pieces of sampling data is calculated as the period deviation between the reference code period (0th code period) and the nth code period. The coefficient is “S n ”. That is, among the sets of sampling data as the received signal portion in each code period in the received signal, the period deviation coefficient is calculated for each code period using sampling data having the same timing within the code period.

図2の受信信号積算処理に戻って、周期ずれ係数算出処理を行った後、受信信号積算処理を行う(ステップA5)。受信信号積算処理では、受信信号のサンプリングデータ「rn,m」と、周期ずれ係数算出処理で算出された周期ずれ係数「Sn」とを用いて、次式(6)に従って積算サンプリングデータ「Rm」を算出する。

Figure 2011137802
但し、「CAn,m(t+nTCA)=CAn,m(t)」となる性質を利用している。 Returning to the reception signal integration process of FIG. 2, after performing the period deviation coefficient calculation process, the reception signal integration process is performed (step A5). In the received signal integration process, the integrated sampling data “r n, m ” and the period deviation coefficient “S n ” calculated in the period deviation coefficient calculation process are used according to the following equation (6). to calculate the R m ".
Figure 2011137802
However, the property of “CA n, m (t + nT CA ) = CA n, m (t)” is used.

図4は、受信信号積算の原理の説明図である。各コード周期それぞれにおけるサンプリングデータ「rn,m」と、対応するコード周期の周期ずれ係数「Sn」とを積和することで積算サンプリングデータ「Rm」を算出する。具体的には、図4に示すように、各サンプリング番号「m」(m=0〜M)それぞれについて、第nコード周期のサンプリングデータ「rn,m」と、第nコード周期の周期ずれ係数「Sn」とを乗算して、乗算値「rn,m・Sn」を算出する。そして、N+1個のコード周期それぞれについて得られた乗算値「rn,m・Sn」を合算することで、各サンプリング番号「m」それぞれについて積算サンプリングデータ「Rm」を算出する。最終的に、各サンプリング番号「m」それぞれについて得られた1周期分の積算サンプリングデータ「Rm」の組が、積算された受信信号(以下、「積算受信信号」と称す。)となる。 FIG. 4 is an explanatory diagram of the principle of reception signal integration. The integrated sampling data “R m ” is calculated by multiplying the sampling data “r n, m ” in each code period by the period deviation coefficient “S n ” of the corresponding code period. Specifically, as shown in FIG. 4, for each sampling number “m” (m = 0 to M), sampling data “rn , m ” in the nth code period and a period shift in the nth code period The product “r n, m · S n ” is calculated by multiplying the coefficient “S n ”. Then, the integrated sampling data “R m ” is calculated for each sampling number “m” by adding the multiplication values “r n, m · S n ” obtained for each of the N + 1 code periods. Finally, a set of integrated sampling data “R m ” for one period obtained for each sampling number “m” becomes an integrated received signal (hereinafter referred to as “integrated received signal”).

式(5)からわかるように、周期ずれ係数「Sn」は、第0コード周期のサンプリングデータ「r0,m(t)」と、第nコード周期のサンプリングデータの複素共役(複素数の虚部の符号を反転させたもの)「{rn,m(t+nT)}」とを乗算することで算出され、「Sn=exp(−iωnTCA)」で表される。これを見ると、周期ずれ係数「Sn」は、式(3)及び(4)で問題となった「exp(iωnTCA)」の部分の複素共役となっていることがわかる。従って、サンプリングデータ「rn,m」に周期ずれ係数「Sn」を乗算することで「exp(iωnTCA)」の部分を消去することができる。 As can be seen from the equation (5), the period deviation coefficient “S n ” is the complex conjugate of the sampling data “r 0, m (t)” of the 0th code period and the sampling data of the nth code period (imaginary complex number). It is calculated by multiplying by “{r n, m (t + nT)} * ” and represented by “S n = exp (−iωnT CA )”. From this, it can be seen that the period deviation coefficient “S n ” is a complex conjugate of the part of “exp (iωnT CA )” that is a problem in the equations (3) and (4). Therefore, the portion of “exp (iωnT CA )” can be eliminated by multiplying the sampling data “r n, m ” by the period shift coefficient “S n ”.

このようにして、各コード周期「n」について「exp(iωnTCA)」の部分を消去し、その上で、各コード周期nについて得られた乗算値「rn,m・Sn」を足し上げることで、「Σexp(iωnTCA)」の項を含まない積算受信信号を得ることができる。 In this way, the part of “exp (iωnT CA )” is deleted for each code period “n”, and then the multiplication value “r n, m · S n ” obtained for each code period n is added. By increasing the value, an integrated reception signal that does not include the term “Σexp (iωnT CA )” can be obtained.

式(3)と式(6)とを比較すると、積算サンプリングデータ「Rm」の式中の「Σexp(iωnTCA)」の項が、コード周期の総数である「N+1」に置き換わっていることがわかる。「N+1」は定数であり、値が変化することはない。従って、式(6)に従って算出される積算サンプリングデータ「Rm」は、利得が増した強い信号となる。 Comparing equation (3) with equation (6), the term “Σexp (iωnT CA )” in the equation of accumulated sampling data “R m ” is replaced with “N + 1”, which is the total number of code periods. I understand. “N + 1” is a constant and the value does not change. Therefore, the integrated sampling data “R m ” calculated according to the equation (6) is a strong signal with an increased gain.

以上により、周期ずれ係数「S」を用いて受信信号「r(t)」を積算することで、GPS衛星信号の受信信号から搬送波「exp(iωt)」を除去することなく、利得が増した高品質な積算受信信号を得ることができる。   As described above, the gain is increased without removing the carrier wave “exp (iωt)” from the received signal of the GPS satellite signal by integrating the received signal “r (t)” using the period shift coefficient “S”. A high-quality integrated reception signal can be obtained.

尚、式(2)の「ωd」をドップラー周波数及びGPS受信装置内のローカルクロックの誤差(時計誤差)に起因する周波数誤差としたが、その他の周波数誤差を更に含めて考えることができる。「ωd」がいかなる値であっても、上記の受信信号積算手法により相関処理に適した信号を得ることができる。 In addition, although “ω d ” in Expression (2) is a frequency error caused by an error (clock error) of the Doppler frequency and the local clock in the GPS receiver, other frequency errors can be considered. Regardless of the value of “ω d ”, a signal suitable for correlation processing can be obtained by the received signal integration method.

2.実施例
次に、上述した原理を適用したGPS受信装置の実施例について説明する。ここでは、GPS受信装置を搭載した電子機器の一種である携帯型電話機1を具体例として説明する。
2. Embodiment Next, an embodiment of a GPS receiver to which the above-described principle is applied will be described. Here, a mobile phone 1 which is a kind of electronic device equipped with a GPS receiver will be described as a specific example.

2−1.構成
図5は、携帯型電話機1の機能構成を示すブロック図である。携帯型電話機1は、GPSアンテナ9と、GPS受信部10と、ホストCPU(Central Processing Unit)30と、操作部40と、表示部50と、携帯電話用アンテナ60と、携帯電話用無線通信回路部70と、記憶部80とを備えて構成される。
2-1. Configuration FIG. 5 is a block diagram showing a functional configuration of the mobile phone 1. The mobile phone 1 includes a GPS antenna 9, a GPS receiving unit 10, a host CPU (Central Processing Unit) 30, an operation unit 40, a display unit 50, a mobile phone antenna 60, and a mobile phone radio communication circuit. Unit 70 and storage unit 80.

GPSアンテナ9は、GPS衛星から発信されているGPS衛星信号を含むRF(Radio Frequency)信号を受信するアンテナであり、受信信号をGPS受信部10に出力する。   The GPS antenna 9 is an antenna that receives an RF (Radio Frequency) signal including a GPS satellite signal transmitted from a GPS satellite, and outputs a received signal to the GPS receiver 10.

GPS受信部10は、GPSアンテナ9から出力された信号に基づいて携帯型電話機1の位置を計測する位置算出回路であり、いわゆるGPS受信装置に相当する機能ブロックである。GPS受信部10は、RF(Radio Frequency)受信回路部11と、ベースバンド処理回路部20とを備えて構成される。尚、RF受信回路部11と、ベースバンド処理回路部20とは、それぞれ別のLSI(Large Scale Integration)として製造することも、1チップとして製造することも可能である。   The GPS receiver 10 is a position calculation circuit that measures the position of the mobile phone 1 based on a signal output from the GPS antenna 9, and is a functional block corresponding to a so-called GPS receiver. The GPS receiving unit 10 includes an RF (Radio Frequency) receiving circuit unit 11 and a baseband processing circuit unit 20. The RF receiving circuit unit 11 and the baseband processing circuit unit 20 can be manufactured as separate LSIs (Large Scale Integration) or can be manufactured as one chip.

RF受信回路部11は、RF信号の処理回路ブロックであり、所定の発振信号を分周或いは逓倍することで、RF信号乗算用の発振信号を生成する。そして、生成した発振信号を、GPSアンテナ9から出力されたRF信号に乗算することで、RF信号を中間周波数の信号(以下、「IF(Intermediate Frequency)信号」と称す。)にダウンコンバートし、IF信号を増幅等した後、A/D変換器でデジタル信号に変換して、ベースバンド処理回路部20に出力する。   The RF receiving circuit unit 11 is an RF signal processing circuit block, and generates an oscillation signal for RF signal multiplication by dividing or multiplying a predetermined oscillation signal. Then, by multiplying the generated oscillation signal by the RF signal output from the GPS antenna 9, the RF signal is down-converted to an intermediate frequency signal (hereinafter referred to as an "IF (Intermediate Frequency) signal"), After the IF signal is amplified, it is converted into a digital signal by an A / D converter and output to the baseband processing circuit unit 20.

ベースバンド処理回路部20は、RF受信回路部11から出力されたIF信号に対して相関演算処理等を行ってGPS衛星信号を捕捉・抽出し、データを復号して航法メッセージや時刻情報等を取り出す回路部である。   The baseband processing circuit unit 20 performs correlation calculation processing or the like on the IF signal output from the RF receiving circuit unit 11 to acquire and extract GPS satellite signals, decodes the data, and displays navigation messages, time information, and the like. It is a circuit part to take out.

図6は、ベースバンド処理回路部20の回路構成の一例を示す図である。ベースバンド処理回路部20は、衛星信号捕捉部21と、CPU25と、記憶部27とを備えて構成される。   FIG. 6 is a diagram illustrating an example of a circuit configuration of the baseband processing circuit unit 20. The baseband processing circuit unit 20 includes a satellite signal capturing unit 21, a CPU 25, and a storage unit 27.

衛星信号捕捉部21は、RF受信回路部11から出力されたIF信号である受信信号からGPS衛星信号を捕捉する回路部であり、受信信号積算処理回路部211と、レプリカ信号発生部213と、相関処理部215とを備えて構成される。   The satellite signal capturing unit 21 is a circuit unit that captures a GPS satellite signal from a reception signal that is an IF signal output from the RF reception circuit unit 11, and includes a reception signal integration processing circuit unit 211, a replica signal generation unit 213, And a correlation processing unit 215.

受信信号積算処理回路部211は、RF受信回路部11から出力されたIF信号である受信信号「r(t)」を積算する処理を行う回路部であり、積算受信信号「R(t)」を相関処理部215に出力する。本実施形態では、受信信号積算処理回路部211は、デジタルシグナルプロセッサー(DSP(Digital Signal Processor))等のプロセッサーとメモリとを備え、図2で説明したフローチャートに従って、受信信号積算処理をデジタル信号処理として実行することとして説明する。   The reception signal integration processing circuit unit 211 is a circuit unit that performs processing to integrate the reception signal “r (t)” that is an IF signal output from the RF reception circuit unit 11, and the integration reception signal “R (t)”. Is output to the correlation processing unit 215. In this embodiment, the received signal integration processing circuit unit 211 includes a processor such as a digital signal processor (DSP) and a memory, and performs the received signal integration processing according to the flowchart described in FIG. Will be described as being executed.

受信信号積算処理回路部211は、各種のデータを格納するメモリとしての記憶部212を備えている。記憶部212には、例えば、受信信号をサンプリングすることで得られる受信信号のサンプリングデータ2121と、各コード周期について算出された周期ずれ係数のデータである周期ずれ係数データ2122と、積算された受信信号のデータである積算受信信号データ2123とが格納される。   The reception signal integration processing circuit unit 211 includes a storage unit 212 as a memory for storing various data. In the storage unit 212, for example, received signal sampling data 2121 obtained by sampling the received signal, period deviation coefficient data 2122 which is data of the period deviation coefficient calculated for each code period, and the integrated reception The accumulated reception signal data 2123 which is signal data is stored.

受信信号積算処理回路部211は、受信信号のサンプリングデータ2121を用いて周期ずれ係数を算出する周期ずれ係数算出部として機能するとともに、受信信号のサンプリングデータ2121を周期ずれ係数を用いて積算する受信信号積算部として機能する。   The reception signal integration processing circuit unit 211 functions as a period deviation coefficient calculation unit that calculates a period deviation coefficient using the sampling data 2121 of the reception signal, and receives the sampling data 2121 of the reception signal using the period deviation coefficient. Functions as a signal integration unit.

従来のGPS受信装置では、受信信号「r(t)」から搬送波「exp(iωt)」を除去してから受信信号「r(t)」を積算する必要があったため、例えば受信信号積算処理回路部211に搬送波を除去するための検波部(キャリア再生部)を設ける必要があった。しかし、本実施形態では、原理で説明したように周期ずれ係数を用いて受信信号を積算するため、受信信号積算処理回路部211に検波部を設ける必要がない。   In the conventional GPS receiver, it is necessary to integrate the received signal “r (t)” after removing the carrier wave “exp (iωt)” from the received signal “r (t)”. It is necessary to provide a detection unit (carrier reproduction unit) for removing the carrier wave in the unit 211. However, in the present embodiment, as described in the principle, the reception signals are integrated using the period shift coefficient, so that it is not necessary to provide a detection unit in the reception signal integration processing circuit unit 211.

レプリカ信号発生部213は、GPS衛星信号のCAコードの拡散符号レプリカの発生信号であるレプリカ信号を生成する回路部である。レプリカ信号発生部213は、CPU25から出力されるCAコード指示信号(捕捉対象衛星の指示信号)に従ったレプリカ信号「CAR(t)」を生成して、相関処理部215に出力する。 The replica signal generation unit 213 is a circuit unit that generates a replica signal that is a generation signal of a spread code replica of a CA code of a GPS satellite signal. The replica signal generation unit 213 generates a replica signal “CA R (t)” according to the CA code instruction signal (capture target satellite instruction signal) output from the CPU 25 and outputs the replica signal “CA R (t)” to the correlation processing unit 215.

相関処理部215は、受信信号積算処理回路部211から入力した積算受信信号「R(t)」と、レプリカ信号発生部213から入力したレプリカ信号「CAR(t)」との相関処理を行う回路部である。相関処理部215は、CPU25から入力した位相指示信号に従って、レプリカ信号の位相「Δt」を変化させながら「R(t)」と「CAR(t+Δt)」との相関を計算し、その相関値「P(Δt)」をCPU25に出力する。 Correlation processing unit 215 performs the cumulative received signal input from reception signal integration processing circuit unit 211 "R (t)", the correlation between the replica signal input from the replica signal generation unit 213 "CA R (t)" It is a circuit part. Correlation processing unit 215, according to the phase instruction signal input from the CPU 25, to calculate the correlation while changing the replica signal phase "Delta] t" and "R (t)" and "CA R (t + Δt)", the correlation value “P (Δt)” is output to the CPU 25.

CPU25は、記憶部27に記憶されているシステムプログラム等の各種プログラムに従ってベースバンド処理回路部20の各部を統括的に制御するプロセッサーである。CPU25は、各捕捉対象衛星について、相関処理部215から出力される相関値「P(Δt)」に基づいてコード位相を検出する処理を行う。そして、検出したコード位相を用いて当該捕捉対象衛星と携帯型電話機1間の擬似距離を算出し、算出した擬似距離を利用した位置算出計算を行って、携帯型電話機1の位置を算出する。   The CPU 25 is a processor that comprehensively controls each unit of the baseband processing circuit unit 20 according to various programs such as a system program stored in the storage unit 27. The CPU 25 performs a process of detecting the code phase for each capture target satellite based on the correlation value “P (Δt)” output from the correlation processing unit 215. Then, a pseudo distance between the capture target satellite and the mobile phone 1 is calculated using the detected code phase, and a position calculation calculation using the calculated pseudo distance is performed to calculate the position of the mobile phone 1.

CPU25は、捕捉対象衛星のCAコード(捕捉対象衛星のPRN番号)を指示するためのCAコード指示信号をレプリカ信号発生部213に出力し、捕捉対象衛星のレプリカ信号「CAR(t)」をレプリカ信号発生部213に生成させる。また、レプリカ信号「CAR(t)」の位相「Δt」を指示するための位相指示信号を相関処理部215に出力し、相関処理部215に、レプリカ信号「CAR(t)」の位相「Δt」を変化させながら相関処理を実行させる。 CPU25 outputs a CA code instruction signal for instructing the CA code acquisition target satellite (PRN number of the acquisition target satellite) to the replica signal generation unit 213, a replica signal of the acquisition target satellite "CA R (t)" The replica signal generator 213 generates the signal. Also it outputs a phase instruction signal for instructing the phase "Δt" of the replica signal "CA R (t)" to the correlation processing unit 215, a correlation processing unit 215, the phase of the replica signal "CA R (t)" Correlation processing is executed while changing “Δt”.

記憶部27は、ROM(Read Only Memory)やフラッシュROM、RAM(Random Access Memory)等の記憶装置によって構成され、CPU25がベースバンド処理回路部20を制御するためのシステムプログラムや、位置算出機能を実現するための各種プログラムやデータ等を記憶している。また、CPU25により実行されるシステムプログラム、各種処理プログラム、各種処理の処理中データ、処理結果などを一時的に記憶するワークエリアを形成している。   The storage unit 27 is configured by a storage device such as a ROM (Read Only Memory), a flash ROM, or a RAM (Random Access Memory), and has a system program for the CPU 25 to control the baseband processing circuit unit 20 and a position calculation function. Various programs, data, etc. for realizing are stored. In addition, a work area for temporarily storing a system program executed by the CPU 25, various processing programs, data being processed in various processing, processing results, and the like is formed.

ホストCPU30は、記憶部80に記憶されているシステムプログラム等の各種プログラムに従って携帯型電話機1の各部を統括的に制御するプロセッサーである。ホストCPU30は、ベースバンド処理回路部20から入力した位置情報を表示部50に表示させる処理を行ったり、当該位置情報を利用した各種のアプリケーション処理を行う。   The host CPU 30 is a processor that comprehensively controls each unit of the mobile phone 1 according to various programs such as a system program stored in the storage unit 80. The host CPU 30 performs processing for displaying the position information input from the baseband processing circuit unit 20 on the display unit 50, and performs various application processes using the position information.

操作部40は、例えばタッチパネルやボタンスイッチ等により構成される入力装置であり、押下されたキーやボタンの信号をホストCPU30に出力する。この操作部40の操作により、通話要求やメール送受信要求、位置算出要求等の各種指示入力がなされる。   The operation unit 40 is an input device configured by, for example, a touch panel, a button switch, or the like, and outputs a pressed key or button signal to the host CPU 30. By operating the operation unit 40, various instructions such as a call request, a mail transmission / reception request, and a position calculation request are input.

表示部50は、LCD(Liquid Crystal Display)等により構成され、ホストCPU30から入力される表示信号に基づいた各種表示を行う表示装置である。表示部50には、位置表示画面や時刻情報等が表示される。   The display unit 50 is configured by an LCD (Liquid Crystal Display) or the like, and is a display device that performs various displays based on display signals input from the host CPU 30. The display unit 50 displays a position display screen, time information, and the like.

携帯電話用アンテナ60は、携帯型電話機1の通信サービス事業者が設置した無線基地局との間で携帯電話用無線信号の送受信を行うアンテナである。   The cellular phone antenna 60 is an antenna that transmits and receives cellular phone radio signals to and from a radio base station installed by a communication service provider of the cellular phone 1.

携帯電話用無線通信回路部70は、RF変換回路、ベースバンド処理回路等によって構成される携帯電話の通信回路部であり、携帯電話用無線信号の変調・復調等を行うことで、通話やメールの送受信等を実現する。   The cellular phone wireless communication circuit unit 70 is a cellular phone communication circuit unit configured by an RF conversion circuit, a baseband processing circuit, and the like, and performs modulation and demodulation of the cellular phone radio signal, thereby enabling communication and mailing. Realize transmission / reception and so on.

記憶部80は、ホストCPU30が携帯型電話機1を制御するためのシステムプログラムや、位置算出機能を実現するための各種プログラムやデータ等を記憶する記憶装置である。   The storage unit 80 is a storage device that stores a system program for the host CPU 30 to control the mobile phone 1 and various programs and data for realizing a position calculation function.

2−2.処理の流れ
(1)受信信号積算処理回路部211の処理
図7は、受信信号積算処理回路部211が、図2の受信信号積算処理のステップA3において実行する周期ずれ係数算出処理の流れを示すフローチャートである。特に説明しないが、受信信号積算処理の実行中は、GPSアンテナ9によるRF信号の受信や、RF受信回路部11によるRF信号のIF信号へのダウンコンバージョンが行われ、IF信号に変換された受信信号「r(t)」がベースバンド処理回路部20に随時出力される状態にあるものとする。
2-2. Processing Flow (1) Processing of Received Signal Integration Processing Circuit Unit 211 FIG. 7 shows a flow of periodic shift coefficient calculation processing executed by the received signal integration processing circuit unit 211 in step A3 of the received signal integration processing of FIG. It is a flowchart. Although not specifically described, during reception signal integration processing, reception of an RF signal received by the GPS antenna 9 or down-conversion of the RF signal to an IF signal by the RF reception circuit unit 11 is performed and converted into an IF signal. It is assumed that the signal “r (t)” is being output to the baseband processing circuit unit 20 as needed.

先ず、受信信号積算処理回路部211は、RF受信回路部11から出力された受信信号「r(t)」を所定のサンプリング時間間隔でサンプリングして受信信号のサンプリングデータ「rn,m」を取得し、記憶部212にサンプリングデータ2121として記憶させる(ステップB1)。 First, the reception signal integration processing circuit unit 211 samples the reception signal “r (t)” output from the RF reception circuit unit 11 at a predetermined sampling time interval to obtain the sampling data “r n, m ” of the reception signal. Acquired and stored in the storage unit 212 as sampling data 2121 (step B1).

次いで、受信信号積算処理回路部211は、基準コード周期である第0コード周期(n=0)の周期ずれ係数「S0」を「1」に設定して、記憶部212の周期ずれ係数データ2122に記憶させる(ステップB3)。 Next, the reception signal integration processing circuit unit 211 sets the period deviation coefficient “S 0 ” of the 0th code period (n = 0), which is the reference code period, to “1”, and the period deviation coefficient data stored in the storage unit 212. 2122 (step B3).

その後、受信信号積算処理回路部211は、基準コード周期を除く他のコード周期番号(n=1〜N)それぞれについてループAの処理を実行する(ステップB5〜B19)。また、ループAの処理では、各サンプリング番号「m」(m=0〜M)それぞれについてループBの処理を実行する(ステップB7〜B13)。   Thereafter, the reception signal integration processing circuit unit 211 executes the process of loop A for each of the other code cycle numbers (n = 1 to N) excluding the reference code cycle (steps B5 to B19). In the loop A process, the loop B process is executed for each sampling number “m” (m = 0 to M) (steps B7 to B13).

ループBの処理では、受信信号積算処理回路部211は、サンプリングデータ「rn,m」の複素共役「{rn,m」を演算する(ステップB9)。そして、基準コード周期のサンプリングデータ「r0,m」と、演算した複素共役「{rn,m」とを乗算する(ステップB11)。そして、次のサンプリング番号へと処理を移行する。 In the process of loop B, the received signal integration processing circuit unit 211 calculates the complex conjugate “{r n, m } * ” of the sampling data “r n, m ” (step B9). Then, the sampling data “r 0, m ” of the reference code period is multiplied by the calculated complex conjugate “{r n, m } * ” (step B11). Then, the processing shifts to the next sampling number.

全てのサンプリング番号「m」についてステップB9及びB11の処理を行った後、受信信号積算処理回路部211は、ループBの処理を終了する(ステップB13)。そして、各サンプリング番号「m」それぞれについて得られた乗算値「r0,m・{rn,m」を合算する(ステップB15)。 After performing the processing of steps B9 and B11 for all the sampling numbers “m”, the received signal integration processing circuit unit 211 ends the processing of loop B (step B13). Then, the multiplication values “r 0, m · {r n, m } * ” obtained for each sampling number “m” are added together (step B15).

そして、受信信号積算処理回路部211は、その合算値をサンプリング番号の総数M+1で除算することで第nコード周期の周期ずれ係数「Sn」を算出し、記憶部212の周期ずれ係数データ2122に記憶させる(ステップB17)。そして、次のコード周期番号へと処理を移行する。 The received signal integration processing circuit unit 211 then calculates the period deviation coefficient “S n ” of the nth code period by dividing the sum by the total number of sampling numbers M + 1, and the period deviation coefficient data 2122 in the storage unit 212. (Step B17). Then, the processing shifts to the next code cycle number.

全てのコード周期番号「n」についてステップB7〜B17の処理を行って周期ずれ係数「Sn」を算出したら、受信信号積算処理回路部211は、ループAの処理を終了して(ステップB19)、周期ずれ係数算出処理を終了する。 When the processing of steps B7 to B17 is performed for all the code cycle numbers “n” to calculate the cycle shift coefficient “S n ”, the received signal integration processing circuit unit 211 ends the processing of loop A (step B19). Then, the period deviation coefficient calculation process is terminated.

図8は、受信信号積算処理回路部211が、図2の受信信号積算処理のステップA5において実行する受信信号積算処理の流れを示すフローチャートである。   FIG. 8 is a flowchart showing the flow of the reception signal integration process executed by the reception signal integration processing circuit unit 211 in step A5 of the reception signal integration process of FIG.

先ず、受信信号積算処理回路部211は、各サンプリング番号(m=0〜M)それぞれについてループCの処理を実行する(ステップC1〜C11)。また、ループCの処理では、各コード周期番号(n=0〜N)それぞれについてループDの処理を実行する(ステップC3〜C7)。   First, the received signal integration processing circuit unit 211 executes processing of loop C for each sampling number (m = 0 to M) (steps C1 to C11). In the process of loop C, the process of loop D is executed for each code cycle number (n = 0 to N) (steps C3 to C7).

ループDの処理では、受信信号積算処理回路部211は、サンプリングデータ「rn,m」と周期ずれ係数「Sn」とを乗算する(ステップC5)。そして、次のコード周期番号へと処理を移行する。 In the process of loop D, the received signal integration processing circuit unit 211 multiplies the sampling data “r n, m ” by the period deviation coefficient “S n ” (step C5). Then, the processing shifts to the next code cycle number.

全てのコード周期番号「n」についてステップC5の処理を行った後、受信信号積算処理回路部211は、ループDの処理を終了する(ステップC7)。その後、各コード周期番号「n」それぞれについて得られた乗算結果「rn,m・Sn」を合算することで、積算受信信号「R」のうちのサンプリング番号「m」の部分を算出し、記憶部212の積算受信信号データ2123に記憶させる(ステップC9)。そして、次のサンプリング番号へと処理を移行する。 After performing the process of step C5 for all code cycle numbers “n”, the received signal integration processing circuit unit 211 ends the process of loop D (step C7). After that, by multiplying the multiplication results “r n, m · S n ” obtained for each code cycle number “n”, the portion of the sampling number “m” in the integrated received signal “R” is calculated. Then, the accumulated reception signal data 2123 of the storage unit 212 is stored (step C9). Then, the processing shifts to the next sampling number.

全てのサンプリング番号「m」についてステップC3〜C9の処理を行うことで積算受信信号「R(t)」を算出した後、受信信号積算処理回路部211は、ループCの処理を終了して(ステップC11)、受信信号積算処理を終了する。   After calculating the integrated reception signal “R (t)” by performing the processing of steps C3 to C9 for all sampling numbers “m”, the reception signal integration processing circuit unit 211 ends the processing of the loop C ( Step C11), the received signal integration process is terminated.

(2)CPU25の処理
図9は、ベースバンド処理回路部20のCPU25が実行するベースバンド処理の流れを示すフローチャートである。
(2) Processing of CPU 25 FIG. 9 is a flowchart showing a flow of baseband processing executed by the CPU 25 of the baseband processing circuit unit 20.

先ず、CPU25は、捕捉対象衛星判定処理を行う(ステップD1)。具体的には、不図示の時計部で計時されている現在時刻において、所与の基準位置の天空に位置するGPS衛星を、アルマナックやエフェメリス等の衛星軌道データを用いて判定して、捕捉対象衛星とする。基準位置は、例えば、電源投入後の初回の位置算出の場合は、いわゆるサーバーアシストによって携帯型電話機1の基地局から取得した位置とし、2回目以降の位置算出の場合は、前回の位置算出で求めた最新のGPS算出位置とする等の方法で設定できる。   First, the CPU 25 performs capture target satellite determination processing (step D1). Specifically, at the current time measured by a clock unit (not shown), a GPS satellite located in the sky at a given reference position is determined using satellite orbit data such as an almanac or an ephemeris, and is captured. Satellite. The reference position is, for example, the position acquired from the base station of the mobile phone 1 by so-called server assist in the case of the first position calculation after power-on, and in the case of the second or subsequent position calculation, It can be set by a method such as obtaining the latest calculated GPS position.

次いで、CPU25は、ステップD1で判定した各捕捉対象衛星について、ループEの処理を実行する(ステップD3〜D11)。ループEの処理では、CPU25は、当該GPS衛星信号に含まれる航法メッセージに基づいて、当該捕捉対象衛星の衛星位置、衛星移動速度及び衛星移動方向等の衛星情報を算出する(ステップD5)。そして、CPU25は、コード位相検出処理を行う(ステップD7)。   Next, the CPU 25 executes a process of loop E for each capture target satellite determined in step D1 (steps D3 to D11). In the process of loop E, the CPU 25 calculates satellite information such as the satellite position, satellite moving speed, and satellite moving direction of the capture target satellite based on the navigation message included in the GPS satellite signal (step D5). Then, the CPU 25 performs code phase detection processing (step D7).

図10は、コード位相検出処理の流れを示すフローチャートである。
先ず、CPU25は、当該捕捉対象衛星のCAコードの指示信号をレプリカ信号発生部213に出力する(ステップE1)。そして、CPU25は、位相のサーチ範囲及びサーチ間隔を設定して、位相サーチに使用するサーチ位相を決定する(ステップE3)。
FIG. 10 is a flowchart showing the flow of the code phase detection process.
First, the CPU 25 outputs an instruction signal of the CA code of the capture target satellite to the replica signal generator 213 (step E1). Then, the CPU 25 sets a search range and a search interval for the phase, and determines a search phase to be used for the phase search (step E3).

次いで、CPU25は、ステップE3で設定した各サーチ位相について、ループFの処理を実行する(ステップE5〜E11)。ループFの処理では、CPU25は、当該サーチ位相「Δt」の指示信号を相関処理部215に出力する(ステップE7)。   Next, the CPU 25 executes loop F processing for each search phase set in step E3 (steps E5 to E11). In the process of loop F, the CPU 25 outputs an instruction signal of the search phase “Δt” to the correlation processing unit 215 (step E7).

ステップE7が実行されると、前述したように相関処理部215は、受信信号積算処理回路部211が図1〜図4で説明した原理に基づいて算出した積算受信信号「R(t)」と、レプリカ信号発生部213から入力したレプリカ信号「CAR(t)」との相関処理を行う。相関処理部215は、CPU25から入力した位相指示信号に従って、レプリカ信号の位相「Δt」を変化させながら「R(t)」と「CAR(t+Δt)」との相関を計算し、その相関値「P(Δt)」をCPU25に出力する。 When step E7 is executed, as described above, the correlation processing unit 215 determines that the reception signal integration processing circuit unit 211 calculates the integrated reception signal “R (t)” calculated based on the principle described with reference to FIGS. performs correlation processing with the replica signal input from the replica signal generation unit 213 "CA R (t)". Correlation processing unit 215, according to the phase instruction signal input from the CPU 25, to calculate the correlation while changing the replica signal phase "Delta] t" and "R (t)" and "CA R (t + Δt)", the correlation value “P (Δt)” is output to the CPU 25.

CPU25は、相関処理部215から相関値「P(Δt)」を入力すると、当該相関値「P(Δt)」を記憶部27に記憶させる(ステップE9)。そして、CPU25は、次のサーチ位相へと処理を移行する。   When the correlation value “P (Δt)” is input from the correlation processing unit 215, the CPU 25 stores the correlation value “P (Δt)” in the storage unit 27 (step E9). Then, the CPU 25 shifts the processing to the next search phase.

全てのサーチ位相についてステップE7及びE9の処理を行った後、CPU25は、ループFの処理を終了する(ステップE11)。そして、CPU25は、記憶部27に記憶されている相関値「P(Δt)」が最大となったサーチ位相「Δt」をコード位相に決定する(ステップE13)。そして、CPU25は、コード位相検出処理を終了する。   After performing the processing of steps E7 and E9 for all search phases, the CPU 25 ends the processing of loop F (step E11). Then, the CPU 25 determines the search phase “Δt” at which the correlation value “P (Δt)” stored in the storage unit 27 is maximized as the code phase (step E13). Then, the CPU 25 ends the code phase detection process.

図9のベースバンド処理に戻って、コード位相検出処理を終了した後、CPU25は、ステップD5で算出した衛星情報と、ステップD7で検出したコード位相とを用いて、当該捕捉対象衛星と携帯型電話機1間の擬似距離を算出する(ステップD9)。擬似距離の整数部分は、例えば最新のGPS算出位置と衛星位置とを用いて算出することができ、擬似距離の端数部分は、コード位相を用いて算出することができる。擬似距離を算出した後、CPU25は、次の捕捉対象衛星へと処理を移行する。   Returning to the baseband process of FIG. 9, after the code phase detection process is completed, the CPU 25 uses the satellite information calculated in step D5 and the code phase detected in step D7 to determine the acquisition target satellite and the portable type. The pseudo distance between the telephones 1 is calculated (step D9). The integer part of the pseudorange can be calculated using, for example, the latest GPS calculation position and the satellite position, and the fractional part of the pseudorange can be calculated using the code phase. After calculating the pseudorange, the CPU 25 shifts the processing to the next capture target satellite.

全ての捕捉対象衛星についてステップD5〜D9の処理を行った後、CPU25は、ループEの処理を終了する(ステップD11)。その後、CPU25は、ステップD9で各捕捉対象衛星について算出された擬似距離を利用したGPS位置算出処理を行って携帯型電話機1の位置を算出する(ステップD13)。尚、擬似距離を利用した位置算出計算の詳細については従来公知であるため、詳細な説明を省略する。   After performing the processing of steps D5 to D9 for all the capture target satellites, the CPU 25 ends the processing of loop E (step D11). Thereafter, the CPU 25 calculates the position of the mobile phone 1 by performing GPS position calculation processing using the pseudorange calculated for each capture target satellite in step D9 (step D13). The details of the position calculation calculation using the pseudo distance are well known in the art and will not be described in detail.

次いで、CPU25は、GPS位置算出処理で算出した位置をホストCPU30に出力する(ステップD15)。そして、CPU25は、位置算出を終了するか否かを判定し(ステップD17)、まだ終了しないと判定した場合は(ステップD17;No)、ステップD1に戻る。また、位置算出を終了すると判定した場合は(ステップD17;Yes)、ベースバンド処理を終了する。   Next, the CPU 25 outputs the position calculated by the GPS position calculation process to the host CPU 30 (step D15). Then, the CPU 25 determines whether or not to end the position calculation (step D17). When it is determined that the position calculation has not ended yet (step D17; No), the CPU 25 returns to step D1. If it is determined that the position calculation is to be ended (step D17; Yes), the baseband processing is ended.

2−3.作用効果
ベースバンド処理回路部20の衛星信号捕捉部21において、RF受信回路部11から出力された受信信号が受信信号積算処理回路部211により積算される。すなわち、受信信号がサンプリング時間間隔である1クロック毎にサンプリングされて複数のサンプリングデータが取得される。そして、GPS衛星信号の拡散符号であるCAコードの想定周期で受信信号を時分割した場合のコード周期を異にするサンプリングデータ組を用いて、CAコードの真の周期と想定周期との周期ずれを表す周期ずれ係数が算出される。そして、受信信号が周期ずれ係数を用いて積算される。
2-3. Action In the satellite signal acquisition unit 21 of the baseband processing circuit unit 20, the reception signal output from the RF reception circuit unit 11 is integrated by the reception signal integration processing circuit unit 211. That is, the received signal is sampled every clock that is a sampling time interval, and a plurality of sampling data is acquired. Then, using a sampling data set having different code periods when the received signal is time-divided with an assumed period of the CA code that is a spreading code of the GPS satellite signal, a period deviation between the true period of the CA code and the assumed period Is calculated. Then, the received signal is integrated using a period shift coefficient.

周期ずれは、CAコードの真の周期と想定周期との差である。CAコードの想定周期でGPS衛星信号の受信信号を時分割することで得られるサンプリングデータの組のうち、想定周期を異にするサンプリングデータの組を用いることで、周期ずれを表す係数を算出することができる。そして、この周期ずれ係数を用いて受信信号を積算することで、周期ずれを考慮した適切な積算受信信号を得ることができる。   The period shift is a difference between the true period of the CA code and the assumed period. A coefficient representing a period shift is calculated by using a sampling data set having a different assumed period from among a set of sampling data obtained by time-division of the received signal of the GPS satellite signal with an assumed period of the CA code. be able to. Then, by integrating the reception signals using this period deviation coefficient, it is possible to obtain an appropriate integrated reception signal considering the period deviation.

すなわち、受信信号をサンプリング時間間隔である1クロック毎にサンプリングして、1つのコード周期につきM+1個のサンプリングデータを取得する。そして、例えば第0コード周期を基準コード周期として、基準コード周期のサンプリングデータ(第1受信信号部分)と、他のN個のコード周期のサンプリングデータ(第2受信信号部分)とを用いて、各コード周期それぞれについて周期ずれ係数を算出する。そして、基準コード周期及び他コード周期のN+1個のコード周期のサンプリングデータを積算する際に、対応するコード周期の周期ずれ係数を用いて積算する。   That is, the received signal is sampled every clock that is a sampling time interval, and M + 1 pieces of sampling data are acquired for one code period. Then, for example, using the 0th code period as a reference code period, using sampling data (first received signal part) of the reference code period and sampling data (second received signal part) of other N code periods, A period deviation coefficient is calculated for each code period. Then, when sampling data of N + 1 code periods of the reference code period and the other code periods are integrated, they are integrated using the period deviation coefficient of the corresponding code period.

サンプリングデータを積算する際に、予めサンプリングデータに周期ずれ係数を乗算することで、周期ずれによる誤差成分を除去することができる。そして、周期ずれによる誤差成分が除去された信号を積算することで、振幅(パワー)が定数倍された相関処理に適した信号を得ることができる。これは言い換えるならば、搬送波が除去されていない状態のGPS衛星信号の受信信号を周期ずれ係数を用いて積算することで、周期ずれに起因する誤差成分を含まない高品質な積算受信信号を得ることができることを意味する。搬送波の除去が不用であるため、ベースバンド処理回路部20にキャリア再生部を設ける必要がなく、GPS受信回路の簡素化が実現される。   When integrating sampling data, an error component due to the period deviation can be removed by multiplying the sampling data by a period deviation coefficient in advance. Then, by integrating the signals from which the error component due to the period deviation is removed, a signal suitable for correlation processing in which the amplitude (power) is multiplied by a constant can be obtained. In other words, a high-quality integrated received signal that does not include an error component due to the period deviation is obtained by integrating the GPS satellite signal reception signals from which carrier waves have not been removed using the period deviation coefficient. Means that you can. Since it is unnecessary to remove the carrier wave, it is not necessary to provide a carrier reproducing unit in the baseband processing circuit unit 20, and the GPS receiving circuit can be simplified.

3.変形例
3−1.適用システム
上述した実施形態では、GPS衛星信号の捕捉を例に挙げて説明したが、GPS衛星信号以外の信号を受信する受信装置についても、本発明を同様に適用可能である。すなわち、拡散符号で拡散変調された衛星信号の受信信号を積算し、積算した信号に対する相関処理を行って衛星信号を捕捉する受信装置であれば本発明を適用可能である。
3. Modification 3-1. Application System In the above-described embodiment, the acquisition of a GPS satellite signal has been described as an example. However, the present invention can be similarly applied to a receiving apparatus that receives a signal other than a GPS satellite signal. In other words, the present invention can be applied to any receiving apparatus that integrates received signals of satellite signals that have been spread-modulated with spreading codes, and performs correlation processing on the integrated signals to capture satellite signals.

3−2.電子機器
また、上述した実施形態では、電子機器の一種である携帯型電話機に本発明を適用した場合を例に挙げて説明したが、本発明を適用可能な電子機器はこれに限られるわけではない。例えば、カーナビゲーション装置や携帯型ナビゲーション装置、パソコン、PDA(Personal Digital Assistant)、腕時計といった他の電子機器についても同様に適用することが可能である。
3-2. In the above-described embodiment, the case where the present invention is applied to a mobile phone that is a kind of electronic device has been described as an example. However, electronic devices to which the present invention can be applied are not limited thereto. Absent. For example, the present invention can be similarly applied to other electronic devices such as a car navigation device, a portable navigation device, a personal computer, a PDA (Personal Digital Assistant), and a wristwatch.

3−3.衛星位置算出システム
上述した実施形態では、衛星位置算出システムとしてGPSを例に挙げて説明したが、WAAS(Wide Area Augmentation System)、QZSS(Quasi Zenith Satellite System)、GLONASS(GLObal NAvigation Satellite System)、GALILEO等の他の衛星位置算出システムであってもよい。
3-3. In the above-described embodiment, the GPS is described as an example of the satellite position calculation system, but WAAS (Wide Area Augmentation System), QZSS (Quasi Zenith Satellite System), GLONASS (GLObal NAvigation Satellite System), GALILEO Other satellite position calculation systems may be used.

3−4.周期ずれ係数算出処理
上述した実施形態では、サンプリングデータ「r0,m」とサンプリングデータの複素共役「{rn,m」とをサンプリング番号「m=0〜M」を揃えてそれぞれ乗算することでM+1個の乗算結果「r0,m・{rn,m」を取得し、それらを平均して第nコード周期の周期ずれ係数「Sn」とするものとして説明したが、次のようにして周期ずれ係数を算出してもよい。すなわち、各コード周期の全てのサンプリングデータの組を用いるのではなく、各コード周期について任意に選択した1又は2以上の受信信号部分としてのサンプリングデータを用いて周期ずれ係数「Sn」を算出する。
3-4. In the embodiment described above, the sampling data “r 0, m ” and the complex conjugate “{r n, m } * ” of the sampling data are multiplied by the sampling numbers “m = 0 to M”, respectively. In the above description, M + 1 multiplication results “r 0, m · {r n, m } * ” are acquired and averaged to obtain a period deviation coefficient “S n ” of the nth code period. The period deviation coefficient may be calculated as follows. That is, instead of using a set of all sampling data for each code period, the period deviation coefficient “S n ” is calculated using sampling data as one or more received signal portions arbitrarily selected for each code period. To do.

任意に選択したL個(1≦L≦M+1)のサンプリング番号のサンプリングデータを用いる場合は、選択したサンプリング番号それぞれについての乗算結果をL個で平均することで、周期ずれ係数「Sn」を同様に算出することができる。例えば、任意に選択した1個のサンプリング番号「m」のサンプリングデータを用いる場合は、次式(7)に従って周期ずれ係数「Sn」を算出すればよい。

Figure 2011137802
In the case of using sampling data with arbitrarily selected L (1 ≦ L ≦ M + 1) sampling numbers, the multiplication result for each of the selected sampling numbers is averaged with L to obtain the period deviation coefficient “S n ”. It can be calculated similarly. For example, when sampling data of one sampling number “m” arbitrarily selected is used, the period deviation coefficient “S n ” may be calculated according to the following equation (7).
Figure 2011137802

3−5.受信信号積算処理
上述した実施例では、受信信号積算処理回路部211が受信信号の積算をデジタル信号処理としてソフトウェア的に行うものとして説明した。しかし、ソフトウェア的に行うのではなく、論理回路等の回路素子を用いたデジタル回路で構成することも可能である。
3-5. Reception signal integration processing In the above-described embodiments, the reception signal integration processing circuit unit 211 has been described as performing software integration of reception signals as digital signal processing. However, it is also possible to configure with a digital circuit using circuit elements such as a logic circuit instead of using software.

1 携帯型電話機、 10 GPS受信部、 11 RF受信回路部、
20 ベースバンド処理回路部、 21 衛星信号捕捉部、 25 CPU、
27 記憶部、 30 ホストCPU、 40 操作部、 50 表示部、
60 携帯電話用アンテナ、 70 携帯電話用無線通信回路部、 80 記憶部、
211 受信信号積算処理回路部、 213 レプリカ信号発生部、
215 相関処理部
1 mobile phone, 10 GPS receiver, 11 RF receiver,
20 baseband processing circuit section, 21 satellite signal capturing section, 25 CPU,
27 storage unit, 30 host CPU, 40 operation unit, 50 display unit,
60 cellular phone antenna, 70 wireless communication circuit unit for cellular phone, 80 storage unit,
211 reception signal integration processing circuit unit, 213 replica signal generation unit,
215 Correlation processing unit

Claims (8)

衛星信号を受信する際に前記衛星信号の拡散符号の符号周期時間を推測した想定周期で前記衛星信号の受信信号を時分割した場合に、前記受信信号のうちの第1受信信号部分と前記第1受信信号部分とは前記想定周期が異なる第2受信信号部分とを用いて、前記拡散符号の真の周期と前記想定周期との周期ずれを表す周期ずれ係数を算出することと、
前記受信信号を前記周期ずれ係数を用いて積算することと、
を含む受信信号積算方法。
When the satellite signal reception signal is time-divided with an assumed period in which the code period time of the spread code of the satellite signal is estimated when the satellite signal is received, the first reception signal portion of the reception signal and the first Calculating a period deviation coefficient representing a period deviation between the true period of the spreading code and the assumed period using a second received signal part having a different assumed period from the one received signal part;
Integrating the received signal using the periodic shift coefficient;
A received signal integrating method including:
前記受信信号を積算することは、前記第1受信信号部分及び前記第2受信信号部分と前記周期ずれ係数とを積和することで前記受信信号を積算した信号を生成することである、
請求項1に記載の受信信号積算方法。
The integration of the reception signals is to generate a signal obtained by integrating the reception signals by multiplying and summing the first reception signal portion and the second reception signal portion and the period deviation coefficient.
The received signal integration method according to claim 1.
前記周期ずれ係数を算出することは、前記受信信号のうち、前記想定周期がn個分(nは自然数)異なる前記第2受信信号部分それぞれについての前記周期ずれ係数を算出することを含み、
前記受信信号を積算することは、前記想定周期がn個分異なる前記第2受信信号部分それぞれを積算する際に、対応する前記周期ずれ係数を用いて積算することを含む、
請求項1又は2に記載の受信信号積算方法。
Calculating the period deviation coefficient includes calculating the period deviation coefficient for each of the second received signal portions of the received signal, the assumption period being different by n (n is a natural number);
Accumulating the received signals includes integrating using the corresponding period deviation coefficients when integrating each of the second received signal parts having different n assumed periods.
The received signal integration method according to claim 1 or 2.
前記周期ずれ係数を算出することは、前記第1及び前記第2受信信号部分において、前記想定周期内で同一タイミングとなる信号部分を用いて、前記周期ずれ係数を算出することを含む、
請求項1〜3の何れか一項に記載の受信信号積算方法。
Calculating the period deviation coefficient includes calculating the period deviation coefficient using a signal part having the same timing within the assumed period in the first and second received signal parts.
The received signal integration method according to any one of claims 1 to 3.
前記周期ずれ係数を算出することは、前記第1受信信号部分と、前記第2受信信号部分の複素共役とを乗算することで前記周期ずれ係数を算出することを含む、
請求項1〜4の何れか一項に記載の受信信号積算方法。
Calculating the period deviation coefficient includes calculating the period deviation coefficient by multiplying the first received signal portion by a complex conjugate of the second received signal part.
The received signal integration method according to any one of claims 1 to 4.
前記周期ずれ係数を算出することは、
前記第1受信信号部分のうちの前記想定周期内で異なるタイミングとなる信号部分と、第n番目(nは自然数)の前記第2受信信号部分のうちの対応する同一タイミングとなる信号部分の複素共役とを乗算することと、
第n番目の前記第2受信信号部分に関して得られた前記各タイミングの前記乗算結果の平均を、前記第1受信信号部分と第n番目の前記第2受信信号部分との周期ずれ係数とすることと、
を含む、
請求項4に記載の受信信号積算方法。
Calculating the period deviation coefficient is
A complex of a signal portion having a different timing within the assumed period of the first received signal portion and a corresponding signal portion of the nth (n is a natural number) second received signal portion corresponding to the same timing. Multiplying with the conjugate,
An average of the multiplication results of the respective timings obtained for the nth second received signal portion is used as a period shift coefficient between the first received signal portion and the nth second received signal portion. When,
including,
The received signal integration method according to claim 4.
前記受信信号を積算することは、前記衛星信号の搬送波が除去されていない状態の受信信号を前記周期ずれ係数を用いて積算することである、
請求項1〜6の何れか一項に記載の受信信号積算方法。
Accumulating the received signal is to integrate the received signal in a state in which the carrier wave of the satellite signal is not removed, using the period deviation coefficient.
The received signal integration method according to any one of claims 1 to 6.
衛星信号を受信する際に前記衛星信号の拡散符号の符号周期時間を推測した想定周期で前記衛星信号の受信信号を時分割した場合に、前記受信信号のうちの第1受信信号部分と前記第1受信信号部分とは前記想定周期が異なる第2受信信号部分とを用いて、前記拡散符号の真の周期と前記想定周期との周期ずれを表す周期ずれ係数を算出する周期ずれ係数算出部と、
前記受信信号を前記周期ずれ係数を用いて積算する受信信号積算部と、
前記受信信号積算部により積算された信号に対する相関処理を行う相関処理部と、
前記相関処理の結果に基づいて前記衛星信号を捕捉する捕捉部と、
を備えた受信装置。
When the satellite signal reception signal is time-divided with an assumed period in which the code period time of the spread code of the satellite signal is estimated when the satellite signal is received, the first reception signal portion of the reception signal and the first A period deviation coefficient calculating unit that calculates a period deviation coefficient representing a period deviation between a true period of the spread code and the assumed period using a second received signal part having a different assumed period from the one received signal part; ,
A received signal integrating unit that integrates the received signal using the periodic deviation coefficient;
A correlation processing unit that performs correlation processing on the signal integrated by the received signal integration unit;
A capturing unit that captures the satellite signal based on a result of the correlation processing;
A receiving device.
JP2010220172A 2010-09-30 2010-09-30 Received signal integrating method and receiving device Pending JP2011137802A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010220172A JP2011137802A (en) 2010-09-30 2010-09-30 Received signal integrating method and receiving device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010220172A JP2011137802A (en) 2010-09-30 2010-09-30 Received signal integrating method and receiving device

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2009294831A Division JP2011133424A (en) 2009-12-25 2009-12-25 Received signal integration method and receiving device

Publications (2)

Publication Number Publication Date
JP2011137802A true JP2011137802A (en) 2011-07-14
JP2011137802A5 JP2011137802A5 (en) 2012-05-17

Family

ID=44349344

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010220172A Pending JP2011137802A (en) 2010-09-30 2010-09-30 Received signal integrating method and receiving device

Country Status (1)

Country Link
JP (1) JP2011137802A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011133424A (en) * 2009-12-25 2011-07-07 Seiko Epson Corp Received signal integration method and receiving device
CN110109150A (en) * 2019-04-22 2019-08-09 长沙翼盾电子科技有限公司 A kind of high-precision array signal simulator and method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003014841A (en) * 2001-06-27 2003-01-15 Mitsubishi Electric Corp Radar apparatus and coherent integration method
JP2011133424A (en) * 2009-12-25 2011-07-07 Seiko Epson Corp Received signal integration method and receiving device

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003014841A (en) * 2001-06-27 2003-01-15 Mitsubishi Electric Corp Radar apparatus and coherent integration method
JP2011133424A (en) * 2009-12-25 2011-07-07 Seiko Epson Corp Received signal integration method and receiving device

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011133424A (en) * 2009-12-25 2011-07-07 Seiko Epson Corp Received signal integration method and receiving device
CN110109150A (en) * 2019-04-22 2019-08-09 长沙翼盾电子科技有限公司 A kind of high-precision array signal simulator and method

Similar Documents

Publication Publication Date Title
US9065549B2 (en) Correlation calculating method, satellite signal capturing method, and correlation calculating circuit
JP5524226B2 (en) Method for processing combined navigation signals
JP5974487B2 (en) Satellite signal capturing method and satellite signal capturing apparatus
KR20080112123A (en) Coherent integration enhancement method, positioning method, storage medium, coherent integration enhancement circuit, positioning circuit, and electronic instrument
JP5321500B2 (en) Signal acquisition method
JP2009097897A (en) Positioning method, program, positioning device, and electronic device
JP5321499B2 (en) Signal acquisition method
US9229113B2 (en) Receiving device and correlation integrating method
JP4911220B2 (en) Satellite signal capturing method and satellite signal receiving apparatus
US20110235687A1 (en) Signal acquisition method, signal acquisition apparatus and electronic device
JP2011137802A (en) Received signal integrating method and receiving device
JP2011133424A (en) Received signal integration method and receiving device
JP5614029B2 (en) Acquisition frequency determination method
JP2011203234A (en) Signal acquisition method, signal acquisition apparatus and electronic device
JP5488380B2 (en) Satellite signal capturing method and satellite signal receiving apparatus
US8761232B2 (en) Bit change determination method and bit change determination device
JP2009300204A (en) Search frequency correction method, positioning method, program, and receiving device
JP5494394B2 (en) Receiver
JP2011174912A (en) Method and device for signal acquisition
JP6326767B2 (en) Satellite signal receiving method for positioning, satellite signal receiving device for positioning, and electronic device
US8022866B2 (en) Method and system for doppler estimation
JP2011174913A (en) Method and device for signal acquisition
JP2011078052A (en) Capture frequency-determining method and receiving apparatus
JP2008051574A (en) Program, positioning circuit, and electronic device

Legal Events

Date Code Title Description
A521 Written amendment

Effective date: 20120328

Free format text: JAPANESE INTERMEDIATE CODE: A523

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20121218

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130319

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130520

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140311

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140430

A02 Decision of refusal

Effective date: 20141021

Free format text: JAPANESE INTERMEDIATE CODE: A02