JP6017983B2 - Vehicle position measuring method and vehicle position measuring system - Google Patents

Vehicle position measuring method and vehicle position measuring system Download PDF

Info

Publication number
JP6017983B2
JP6017983B2 JP2013020572A JP2013020572A JP6017983B2 JP 6017983 B2 JP6017983 B2 JP 6017983B2 JP 2013020572 A JP2013020572 A JP 2013020572A JP 2013020572 A JP2013020572 A JP 2013020572A JP 6017983 B2 JP6017983 B2 JP 6017983B2
Authority
JP
Japan
Prior art keywords
satellite
epoch
calculated
doppler shift
calculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2013020572A
Other languages
Japanese (ja)
Other versions
JP2014153084A (en
Inventor
春生 山本
春生 山本
知二 高須
知二 高須
久保 信明
信明 久保
明生 安田
明生 安田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Railway Technical Research Institute
Original Assignee
Railway Technical Research Institute
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 Railway Technical Research Institute filed Critical Railway Technical Research Institute
Priority to JP2013020572A priority Critical patent/JP6017983B2/en
Publication of JP2014153084A publication Critical patent/JP2014153084A/en
Application granted granted Critical
Publication of JP6017983B2 publication Critical patent/JP6017983B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Description

本発明は、車両の位置を計測するための車両位置計測方法等に関する。   The present invention relates to a vehicle position measurement method for measuring the position of a vehicle.

鉄道の安全な運行には列車の位置の把握が不可欠である。
位置検出方式としては、例えば軌道回路を用いる方式、速度発電機を利用する方式(例えば、特許文献1)、GPS(Global Positioning System)による測位情報を利用する方式(例えば、特許文献2、特許文献3)等が知られるところである。
Understanding the location of trains is essential for the safe operation of railways.
As the position detection method, for example, a method using a track circuit, a method using a speed generator (for example, Patent Document 1), a method using positioning information by GPS (Global Positioning System) (for example, Patent Document 2, Patent Document) 3) etc. are known.

特開2012−20734号公報JP 2012-20734 A 特開2006−240478号公報Japanese Patent Laid-Open No. 2006-240478 特開2002−335607号公報JP 2002-335607 A

軌道回路等の地上設備を必要としない、或いは、最低限に抑える車両位置の計測は、例えば閑散線区などにおいて要望が高い。軌道回路に代わってGPSによる測位情報を利用する場合には、計測誤差をいかに小さくするかが重要とされる。GPSの計測誤差の一つとして「マルチパス(多重反射)」が知られる。マルチパスは、地表面や樹木、建物などから反射された信号が、GPS衛星からの直線経路外から届き、直線経路で到達した信号に干渉することで生じる誤差である。なお、GPS衛星以外の衛星測位システムを利用する場合にもマルチパスの問題がある。また、自動車等の列車以外の車両にもこの問題がある。   The measurement of the vehicle position that does not require or minimizes ground facilities such as track circuits is highly demanded, for example, in a quiet line area. When using GPS positioning information instead of the track circuit, it is important to reduce the measurement error. As one of GPS measurement errors, “multipath” is known. Multipath is an error that occurs when a signal reflected from the ground surface, a tree, a building, or the like arrives from outside a straight path from a GPS satellite and interferes with a signal that has reached the straight path. Note that there is a problem of multipath even when a satellite positioning system other than GPS satellites is used. Moreover, vehicles other than trains, such as a motor vehicle, have this problem.

本発明がなされた目的は、衛星測位システムを利用して車両位置を計測する場合に、マルチパスの影響を低減させることにある。   An object of the present invention is to reduce the influence of multipath when measuring a vehicle position using a satellite positioning system.

以上の課題を解決するための第1の発明は、衛星から受信した衛星信号に基づいて車両の位置を計測する車両位置計測方法であって、各候補衛星の衛星信号を捕捉する捕捉ステップ(例えば、図7のステップS2)と、
前記捕捉した衛星信号に係る擬似距離を算出する擬似距離算出ステップ(例えば、図7のステップS8)と、
前記捕捉した衛星信号に係るドップラシフト周波数を算出するドップラシフト周波数算出ステップ(例えば、図7のステップS10)と、
前記候補衛星それぞれについて、前記算出した擬似距離に基づく擬似距離変化量と、前記算出したドップラシフト周波数に基づく当該候補衛星と当該車両との相対移動量との乖離量を算出する乖離量算出ステップ(例えば、図7のステップS24)と、
前記算出した乖離量を用いて、前記候補衛星の中から測位に使用する衛星を選択する衛星選択ステップ(例えば、図7のステップS40)と、
前記選択した衛星の衛星信号に基づいて測位する測位ステップ(例えば、図7のステップS46)と、を含む車両位置計測方法である。
A first invention for solving the above problems is a vehicle position measuring method for measuring the position of a vehicle based on a satellite signal received from a satellite, and a capturing step (for example, capturing a satellite signal of each candidate satellite) , Step S2) of FIG.
A pseudo distance calculating step (for example, step S8 in FIG. 7) for calculating a pseudo distance related to the captured satellite signal;
A Doppler shift frequency calculating step (for example, Step S10 in FIG. 7) for calculating a Doppler shift frequency related to the captured satellite signal;
For each of the candidate satellites, a deviation amount calculating step for calculating a deviation amount between the pseudo distance change based on the calculated pseudo distance and a relative movement amount between the candidate satellite and the vehicle based on the calculated Doppler shift frequency ( For example, step S24) in FIG.
A satellite selection step (for example, step S40 in FIG. 7) for selecting a satellite to be used for positioning from the candidate satellites using the calculated deviation amount;
And a positioning step (for example, step S46 in FIG. 7) for positioning based on the satellite signal of the selected satellite.

また別形態として、衛星から受信した衛星信号に基づいて車両の位置を測位する車両位置計測システムであって、各候補衛星の衛星信号を捕捉する捕捉部(例えば、図3のCPU24、ベースバンド処理部30、捕捉部301、図7のステップS2)と、
前記捕捉した衛星信号に係る擬似距離を算出する擬似距離算出部(例えば、図3のCPU24、ベースバンド処理部30、擬似距離算出部303、図7のステップS8)と、
前記捕捉した衛星信号に係るドップラシフト周波数を算出するドップラシフト周波数算出部(例えば、図3のCPU24、ベースバンド処理部30、ドップラシフト周波数算出部305、図7のステップS10)と、
前記候補衛星それぞれについて、前記算出した擬似距離に基づく擬似距離変化量と、前記算出したドップラシフト周波数に基づく当該候補衛星と当該車両との相対移動量との乖離量を算出する乖離量算出部(例えば、図3のCPU24、乖離量算出部32、図7のステップS24)と、
前記算出した乖離量を用いて、前記候補衛星の中から測位に使用する衛星を選択する衛星選択部(例えば、図3のCPU24、衛星選択部36、図7のステップS40)と、
前記選択した衛星の衛星信号に基づいて測位する測位部(例えば、図3のCPU24、測位演算部38、図7のステップS46))と、を備えた車両位置計測システムを構成できる。
In another embodiment, the vehicle position measurement system measures the position of the vehicle based on the satellite signal received from the satellite, and includes a capturing unit that captures the satellite signal of each candidate satellite (for example, the CPU 24 in FIG. 3, baseband processing). Unit 30, capture unit 301, step S2) of FIG.
A pseudo distance calculation unit (for example, the CPU 24 in FIG. 3, the baseband processing unit 30, the pseudo distance calculation unit 303, and step S8 in FIG. 7) that calculates a pseudo distance related to the captured satellite signal;
A Doppler shift frequency calculation unit for calculating a Doppler shift frequency related to the captured satellite signal (for example, the CPU 24 in FIG. 3, the baseband processing unit 30, the Doppler shift frequency calculation unit 305, Step S10 in FIG. 7);
For each of the candidate satellites, a deviation amount calculation unit that calculates a deviation amount between the pseudo distance change based on the calculated pseudo distance and a relative movement amount between the candidate satellite and the vehicle based on the calculated Doppler shift frequency ( For example, the CPU 24 in FIG. 3, the divergence amount calculation unit 32, step S24 in FIG.
A satellite selection unit (for example, CPU 24 in FIG. 3, satellite selection unit 36, step S40 in FIG. 7) for selecting a satellite to be used for positioning from the candidate satellites using the calculated deviation amount;
A vehicle position measurement system including a positioning unit (for example, the CPU 24, the positioning calculation unit 38, and step S46 in FIG. 7) that performs positioning based on the satellite signal of the selected satellite can be configured.

第1の発明等によれば、衛星信号を用いて車両の位置を計測することができる。そして、位置計測に当たっては、候補衛星それぞれについて、擬似距離に基づく擬似距離変化量と、ドップラシフト周波数に基づく当該候補衛星と車両との相対移動量との乖離量を算出できる。擬似距離変化量は原理的にマルチパスの影響を受け易いが、ドップラシフト周波数はマルチパスの影響を受け難い。よって、これら両者の差である乖離量は、マルチパスの影響度合を示す指数と言える。従って、算出した乖離量を用いて候補衛星の中から測位に使用する衛星を選択することで、マルチパスの影響を低減することができる。   According to the first invention and the like, the position of the vehicle can be measured using the satellite signal. In the position measurement, for each candidate satellite, a deviation amount between the pseudo distance change based on the pseudo distance and the relative movement amount between the candidate satellite and the vehicle based on the Doppler shift frequency can be calculated. The pseudo-range change amount is easily affected by multipath in principle, but the Doppler shift frequency is hardly affected by multipath. Therefore, the amount of divergence, which is the difference between the two, can be said to be an index indicating the degree of multipath influence. Therefore, the influence of multipath can be reduced by selecting a satellite to be used for positioning from the candidate satellites using the calculated deviation amount.

より具体的には、第2の発明として、前記乖離量算出ステップは、第1のエポックで算出した擬似距離と第2のエポックで算出した擬似距離の差から前記擬似距離変化量を算出するステップ(例えば図7のステップS12)と、
前記第1のエポックで算出したドップラシフト周波数と前記第2のエポックで算出したドップラシフト周波数の平均周波数、及び、前記第1のエポックと前記第2のエポック間の時間間隔を用いて前記相対移動量を算出するステップ(例えば、図7のステップS20)と、を含む第1の発明の車両位置計測方法を構成することができる。
More specifically, as the second invention, the divergence amount calculating step calculates the pseudo distance change amount from a difference between the pseudo distance calculated in the first epoch and the pseudo distance calculated in the second epoch. (For example, step S12 in FIG. 7);
The relative movement using the Doppler shift frequency calculated in the first epoch, the average frequency of the Doppler shift frequency calculated in the second epoch, and the time interval between the first epoch and the second epoch. The vehicle position measuring method according to the first aspect of the present invention can include the step of calculating the quantity (for example, step S20 in FIG. 7).

そして、第1の発明及び第2の発明を、鉄道に適用するならば、第3の発明として、前記車両を、軌道を走行する列車とし、前記軌道に沿った区域別に定めた前記乖離量を判定するための判定基準値の中から、前記測位の結果に対応する判定基準値を選択する判定基準値選択ステップ(例えば、図7のステップS26)を更に含み、前記衛星選択ステップを、前記選択した判定基準値と、前記算出した乖離量とを比較し、当該比較結果に基づいて、測位に使用する衛星を選択するステップとした第1又は第2の発明の車両位置計測方法を構成すると良い。   Then, if the first invention and the second invention are applied to a railway, as the third invention, the vehicle is a train traveling on a track, and the deviation amount determined for each area along the track The method further includes a determination reference value selection step (for example, step S26 in FIG. 7) for selecting a determination reference value corresponding to the positioning result from the determination reference values for determination, and the satellite selection step includes the selection. The vehicle position measurement method according to the first or second aspect of the invention may be configured such that the determined reference value is compared with the calculated divergence amount and the satellite used for positioning is selected based on the comparison result. .

第3の発明によれば、第1又は第2の発明と同様の効果が得られるとともに、軌道沿いの区域別の判定基準値を用いることで、場所によってマルチパスの影響が変化するような軌道であっても、マルチパスの影響を適切に低減できる。例えば、比較的マルチパスの影響が高くなる区域では、他の区域よりも判定基準値を大きめに設定するならば、判定基準を厳しくし過ぎた結果、測位に使用できる衛星数が限られてしまい測位できなくなることを避けつつ、マルチパスの影響が大きい衛星の使用を避けることができる。   According to the third invention, the same effect as the first or second invention is obtained, and the trajectory in which the influence of the multipath changes depending on the location by using the criterion value for each area along the trajectory. Even so, the effects of multipath can be reduced appropriately. For example, in an area where the influence of multipath is relatively high, if the criterion value is set larger than that in other regions, the criterion becomes too strict, resulting in a limited number of satellites that can be used for positioning. While avoiding the inability to perform positioning, it is possible to avoid the use of satellites that are highly affected by multipath.

車両位置計測システムの構成例を示す図。The figure which shows the structural example of a vehicle position measurement system. 車両位置計測システムがマルチパスの影響を低減する原理図。The principle diagram in which a vehicle position measurement system reduces the influence of multipath. 車両位置計測システムの機能構成例を示す機能ブロック図。The functional block diagram which shows the function structural example of a vehicle position measurement system. 判定基準値設定データのデータ構成の一例を示す図。The figure which shows an example of the data structure of determination reference value setting data. 擬似距離履歴データのデータ構成の一例を示す図。The figure which shows an example of a data structure of pseudo distance history data. ドップラシフト周波数履歴データのデータ構成の一例を示す図。The figure which shows an example of a data structure of Doppler shift frequency historical data. 一回の測位を行うための処理の流れを説明するためのフローチャート。The flowchart for demonstrating the flow of the process for performing positioning once. 乖離量履歴データのデータ構成の一例を示す図。The figure which shows an example of a data structure of deviation | shift amount log | history data. 変形例における、一回の測位を行うための処理の流れを説明するためのフローチャート。The flowchart for demonstrating the flow of the process for performing positioning once in the modification.

〔第1実施形態〕
本発明を適用した車両位置計測システムについて説明する。尚、本実施形態ではGPSを利用するものとして説明するが、他の衛星測位システムに適宜置き換えることができる。また、車両位置の計測タイミングを以降では「エポック」と呼ぶこととする。
[First Embodiment]
A vehicle position measurement system to which the present invention is applied will be described. In the present embodiment, the GPS is described as being used, but can be appropriately replaced with another satellite positioning system. In addition, the measurement timing of the vehicle position is hereinafter referred to as “epoch”.

[システム構成の説明]
図1は、本実施形態における車両位置計測システムの構成例を示す図である。尚、車両位置計測システム2へ電力を供給する回路等についての説明は省略している。
[Description of system configuration]
FIG. 1 is a diagram illustrating a configuration example of a vehicle position measurement system in the present embodiment. In addition, the description about the circuit etc. which supply electric power to the vehicle position measurement system 2 is abbreviate | omitted.

本実施形態の車両位置計測システム2は、軌道3を走行する車両である列車4の位置を計測するシステムであって列車4に搭載される。
車両位置計測システム2は、GPS衛星6から送信される衛星信号7を受信するGPSアンテナ21と、RF(Radio Frequency)受信部22と、CPU(Central Processing Unit)24と、ICメモリ26と、RTC(Real Time Clock)27と、通信I/F(Interface)28とを有する。RF受信部22、CPU24、ICメモリ26、RTC27および通信I/F28は、バス29によりデータ送受可能に接続されている。
The vehicle position measurement system 2 of the present embodiment is a system that measures the position of a train 4 that is a vehicle traveling on a track 3 and is mounted on the train 4.
The vehicle position measurement system 2 includes a GPS antenna 21 that receives a satellite signal 7 transmitted from a GPS satellite 6, an RF (Radio Frequency) receiver 22, a CPU (Central Processing Unit) 24, an IC memory 26, an RTC (Real Time Clock) 27 and a communication I / F (Interface) 28. The RF receiver 22, the CPU 24, the IC memory 26, the RTC 27, and the communication I / F 28 are connected via a bus 29 so that data can be transmitted and received.

GPSアンテナ21は、GPS衛星6から発信される電波を電気信号に変換してRF受信部22へ出力する。
RF受信部22は、GPSアンテナ21から入力した信号を中間周波数に変換し、デジタル信号に変換してデジタルIF信号として出力する。RF受信部22は、例えば、各種フィルタ回路や低雑音増幅器(LNA:Low Noise Amplifier)の他、中間周波数のIF信号(IF:Intermediate Frequency)へ変換するダウンコンバータ、ダウンコンバートされたIF信号をデジタル信号に変換するA/D変換器、等を有する。尚、RF受信部22は、公知のGPS用のフロントエンドアナログ回路や、フロントエンドICを流用可能である。
The GPS antenna 21 converts a radio wave transmitted from the GPS satellite 6 into an electrical signal and outputs the electrical signal to the RF receiving unit 22.
The RF receiver 22 converts the signal input from the GPS antenna 21 into an intermediate frequency, converts it into a digital signal, and outputs it as a digital IF signal. For example, the RF receiving unit 22 digitally converts various down-converted IF signals into down-converters that convert to IF signals (IF: Intermediate Frequency) in addition to various filter circuits and low-noise amplifiers (LNA). An A / D converter that converts the signal into a signal is included. The RF receiver 22 can use a well-known GPS front-end analog circuit or front-end IC.

CPU24は、(1)RF受信部22から入力されるデジタルIF信号をもとに衛星信号の捕捉や追尾、擬似距離の算出、アルマナック(Almanac)やエフェメリス(Ephemeris)等の航法メッセージのデコードなどを含む所謂ベースバンド処理、(2)測位演算処理、(3)測位結果である位置座標を通信I/F28へ出力する外部出力制御処理、を実行することができる。   The CPU 24 (1) captures and tracks satellite signals based on the digital IF signal input from the RF receiver 22, calculates pseudoranges, decodes navigation messages such as Almanac and Ephemeris, etc. In other words, so-called baseband processing, (2) positioning calculation processing, and (3) external output control processing for outputting position coordinates as positioning results to the communication I / F 28 can be executed.

ICメモリ26は、CPU24による演算処理に必要なプログラムやデータを一時的或いは不揮発に記憶する情報記憶媒体である。RAMやROMにより実現される。また、代替の記憶部として、ハードディスクを備えることとしてもよい。   The IC memory 26 is an information storage medium that temporarily or nonvolatilely stores programs and data necessary for arithmetic processing by the CPU 24. It is realized by RAM or ROM. Moreover, it is good also as providing a hard disk as an alternative memory | storage part.

RTC27は、車両位置計測システム2における基準時刻を出力する。適宜、図示されないバッテリーに接続された構成とすると良い。   The RTC 27 outputs a reference time in the vehicle position measurement system 2. The configuration may be appropriately connected to a battery (not shown).

通信I/F28は、当該計測システム外の装置とのデータ通信を実現する。いわゆるインタフェースIC、通信端子、無線装置等を適宜流用可能である。例えば、車両位置情報(本実施形態では、GPS座標値)等を、当該計測システム外へ出力することができる。   The communication I / F 28 realizes data communication with a device outside the measurement system. A so-called interface IC, communication terminal, wireless device, or the like can be appropriately used. For example, vehicle position information (in this embodiment, a GPS coordinate value) or the like can be output outside the measurement system.

尚、CPU24とICメモリ26は、1チップのシステムLSIとして実現されるとしても良い。また、GPSアンテナ21やRF受信部22を更に含めた一つのモジュールとして実現される構成でも良い。   The CPU 24 and the IC memory 26 may be realized as a one-chip system LSI. Moreover, the structure implement | achieved as one module further including the GPS antenna 21 and RF receiving part 22 may be sufficient.

[原理の説明]
図2は、車両位置計測システム2がマルチパスの影響を低減する原理を説明するための図である。車両位置計測システム2は、エポックt1においてGPSアンテナ21でGPS衛星6(6a)から衛星信号7(7a)を受信する。このとき、受信した衛星信号7(7a)に基づいて、GPS衛星星6(6a)から列車4(4a)までの擬似距離L1と、両者の相対移動に係るドップラシフト周波数fd1とを求めることができる。
[Description of Principle]
FIG. 2 is a diagram for explaining the principle that the vehicle position measurement system 2 reduces the influence of multipath. The vehicle position measurement system 2 receives the satellite signal 7 (7a) from the GPS satellite 6 (6a) with the GPS antenna 21 at the epoch t1. At this time, based on the received satellite signal 7 (7a), the pseudo distance L1 from the GPS satellite star 6 (6a) to the train 4 (4a) and the Doppler shift frequency fd1 related to the relative movement of both are obtained. it can.

同様に、時間間隔△t後のエポックt2において、車両位置計測システム2は、時間間隔△tの間に移動したGPS衛星6(6b)から衛星信号7(7b)を受信する。そして、この時受信した衛星信号に基づいて、GPS衛星6(6b)から列車4(4b)までの擬似距離L2と、両者の相対移動速度に基づくドップラシフト周波数fd2とを求めることができる。   Similarly, at the epoch t2 after the time interval Δt, the vehicle position measurement system 2 receives the satellite signal 7 (7b) from the GPS satellite 6 (6b) moved during the time interval Δt. Based on the satellite signal received at this time, the pseudo distance L2 from the GPS satellite 6 (6b) to the train 4 (4b) and the Doppler shift frequency fd2 based on the relative moving speed of both can be obtained.

ここで、時間間隔△tの間に移動した列車4の移動距離に関連して二つのパラメータを考える。1)擬似距離変化量△Lと、2)GPS衛星6と列車4との相対移動量△Dである。   Here, two parameters are considered in relation to the moving distance of the train 4 moved during the time interval Δt. 1) pseudo-range change ΔL, and 2) relative movement ΔD between the GPS satellite 6 and the train 4.

擬似距離変化量△Lは、最新のエポックt2(第2のエポック)における擬似距離L2と、最も時間的に近い過去のエポックt1(第1のエポック)における擬似距離L1との差から求められる(図2中の数式A)。   The pseudo distance change amount ΔL is obtained from the difference between the pseudo distance L2 in the latest epoch t2 (second epoch) and the pseudo distance L1 in the past epoch t1 (first epoch) closest in time (first epoch) ( Formula A) in FIG.

相対移動距離△Dは、a)最新のエポックt2(第2のエポック)におけるドップラシフト周波数fd2と最も時間的に近い過去のエポックt1(第1のエポック)におけるドップラシフト周波数fd1との平均値と、b)衛星信号の搬送波波長λと、c)最新のエポックt2から最も時間的に近い過去のエポックt1までのエポックの時間間隔△tとの積を求め、その正負符合を逆転することで得られる(図2中の数式B)。   The relative movement distance ΔD is a) the average value of the Doppler shift frequency fd2 at the latest epoch t2 (second epoch) and the Doppler shift frequency fd1 at the past epoch t1 (first epoch) closest in time. B) Obtain the product of the carrier wave wavelength λ of the satellite signal and c) the time interval Δt of the epoch from the latest epoch t2 to the nearest epoch t1 in time, and reverse the sign (Formula B in FIG. 2).

そして、ドップラシフト周波数fdi(i=1,2,…)が原理的に擬似距離Li(i=1,2,…)よりもマルチパスの影響を受け難い点に着目すると、もしマルチパスの影響が皆無で且つその他の誤差要因を考慮しないと仮定すると、擬似距離変化量△Lと相対移動量△Dとが一致するはずである。反対に、マルチパスの影響を受けていれば擬似距離変化量△Lと相対移動量△Dとの乖離が大きくなる。   If the Doppler shift frequency fdi (i = 1, 2,...) Is in principle less susceptible to multipath than the pseudorange Li (i = 1, 2,...). Assuming that there is no error and other error factors are not taken into account, the pseudo distance change amount ΔL and the relative movement amount ΔD should match. On the other hand, if it is affected by the multipath, the difference between the pseudo distance change amount ΔL and the relative movement amount ΔD becomes large.

位置座標を算出するためには少なくとも4機のGPS衛星が必要であり、数が多ければそれだけ高精度且つ高速に位置座標を算出できる。天空には4機を超えるGPS衛星が常に配置されているため、測位に用いるGPS衛星を選択することが可能である。
そこで、本実施形態では擬似距離変化量△Lと相対移動量△Dの差すなわち両者の乖離量Qe(図2中の数式C)を、候補となるGPS衛星(候補衛星)別に求め、それぞれの乖離量Qeを用いて測位に使用する衛星を選択することでマルチパスの影響を低減する。具体的には、本実施形態では、マルチパス誤差を含む許容できる誤差の上限として「判定基準値k」を設定し、乖離量Qeが当該判定基準値kより大きければマルチパスの影響を強く受けていると見なし、測位に使用する衛星の候補から除外することとする。
In order to calculate the position coordinates, at least four GPS satellites are necessary. If the number is large, the position coordinates can be calculated with high accuracy and high speed. Since four or more GPS satellites are always arranged in the sky, it is possible to select a GPS satellite to be used for positioning.
Therefore, in the present embodiment, the difference between the pseudo distance change amount ΔL and the relative movement amount ΔD, that is, the deviation amount Qe (formula C in FIG. 2) is obtained for each candidate GPS satellite (candidate satellite). The influence of multipath is reduced by selecting a satellite to be used for positioning using the deviation amount Qe. Specifically, in this embodiment, “determination reference value k” is set as the upper limit of an allowable error including multipath error, and if the deviation amount Qe is larger than the determination reference value k, it is strongly influenced by multipath. Therefore, it is excluded from the satellite candidates used for positioning.

尚、第1のエポックと第2のエポックとは、1エポック差とすると最も精度良くマルチパスの影響を低減できる。   The first epoch and the second epoch can reduce the influence of the multipath with the highest accuracy when the difference is one epoch.

[機能構成の説明]
図3は、本実施形態における車両位置計測システム2の機能構成例を示すブロック図である。車両位置計測システム2のCPU24は、ベースバンド処理部30と、乖離量算出部32と、判定基準値選択部34と、衛星選択部36と、測位演算部38と、を有する。
[Description of functional configuration]
FIG. 3 is a block diagram illustrating a functional configuration example of the vehicle position measurement system 2 according to the present embodiment. The CPU 24 of the vehicle position measurement system 2 includes a baseband processing unit 30, a deviation amount calculation unit 32, a determination reference value selection unit 34, a satellite selection unit 36, and a positioning calculation unit 38.

ベースバンド処理部30は、衛星信号の捕捉や追跡、航法メッセージのデコードを行う。また、ベースバンド処理部30は、捕捉部301と、擬似距離算出部303と、ドップラシフト周波数算出部305とを有する。尚、デコードされた航法メッセージは、ICメモリ26に記憶される(航法メッセージ44)。   The baseband processing unit 30 captures and tracks satellite signals and decodes navigation messages. The baseband processing unit 30 includes a capturing unit 301, a pseudo distance calculating unit 303, and a Doppler shift frequency calculating unit 305. The decoded navigation message is stored in the IC memory 26 (navigation message 44).

捕捉部301は、測位に使用する候補となるGPS衛星(候補衛星)の衛星信号を捕捉する。例えば、取得済の航法メッセージ44に基づいて各GPS衛星の位置を予測計算し、可視範囲内のGPS衛星のC/Aコード(レプリカコード)を生成し、デジタルIF信号との相関値がピークとなるようにC/Aコード位相制御する。   The capturing unit 301 captures satellite signals of GPS satellites (candidate satellites) that are candidates for use in positioning. For example, the position of each GPS satellite is predicted and calculated based on the acquired navigation message 44, a C / A code (replica code) of the GPS satellite within the visible range is generated, and the correlation value with the digital IF signal is a peak. The C / A code phase control is performed as follows.

擬似距離算出部303は、捕捉できたGPS衛星別に擬似距離Li(i=1,2,…)を算出する。擬似距離の算出方法は従来公知の方法と同様である。例えば、受信したGPS衛星信号に示されている発信時刻と、RTC27の計時時刻に基づく受信時刻との差から擬似距離を求めてもよい。算出された擬似距離Li(i=1,2,…)は、ICメモリ26の擬似距離履歴データ48に衛星IDと、エポックID(例えば、エポックの示す時刻、又は起動から連番付与される識別情報)と対応づけて格納される。   The pseudo distance calculation unit 303 calculates a pseudo distance Li (i = 1, 2,...) For each GPS satellite that has been captured. The calculation method of the pseudo distance is the same as a conventionally known method. For example, the pseudo distance may be obtained from the difference between the transmission time indicated in the received GPS satellite signal and the reception time based on the time measured by the RTC 27. The calculated pseudo distance Li (i = 1, 2,...) Is identified in the pseudo distance history data 48 of the IC memory 26 by the satellite ID and the epoch ID (for example, the time indicated by the epoch or the serial number from the start. Information) and stored.

ドップラシフト周波数算出部305は、捕捉できたGPS衛星別にドップラシフト周波数fdi(i=1,2,…)を算出する。ドップラシフト周波数の算出方法は従来公知の方法と同様である。例えば、捕捉した衛星信号の受信周波数から搬送波周波数(1.57542[GHz])を減算して求めても良い。算出されたドップラシフト周波数fdi(i=1,2,…)は、ICメモリ26のドップラシフト周波数履歴データ50に衛星IDと、エポックIDと対応づけて格納される。   The Doppler shift frequency calculation unit 305 calculates a Doppler shift frequency fdi (i = 1, 2,...) For each GPS satellite that has been captured. The calculation method of the Doppler shift frequency is the same as a conventionally known method. For example, the carrier frequency (1.57542 [GHz]) may be subtracted from the reception frequency of the captured satellite signal. The calculated Doppler shift frequency fdi (i = 1, 2,...) Is stored in the Doppler shift frequency history data 50 of the IC memory 26 in association with the satellite ID and the epoch ID.

尚、本実施形態ではベースバンド処理部30をCPU24の演算処理によって実現する構成としたが、公知のベースバンド処理ICなどのハードウェアを流用し、ベースバンド処理部30をCPU24とは独立したプロセッサで実現することもできる。   In this embodiment, the baseband processing unit 30 is configured to be realized by arithmetic processing of the CPU 24. However, hardware such as a known baseband processing IC is used, and the baseband processing unit 30 is a processor independent of the CPU 24. Can also be realized.

乖離量算出部32は、捕捉されたGPS衛星別に、(1)算出された擬似距離Li(i=1,2,…)に基づく擬似距離変化量△Lと、(2)算出されたドップラシフト周波数fdi(i=1,2,…)に基づく当該GPS衛星と車両との相対移動量△Dと、(3)擬似距離変化量△Lと相対移動量△Dとの差である乖離量Qeと、を算出する。   The divergence amount calculation unit 32 includes (1) a pseudo distance change amount ΔL based on the calculated pseudo distance Li (i = 1, 2,...) And (2) a calculated Doppler shift for each captured GPS satellite. A deviation amount Qe which is a difference between the relative movement amount ΔD between the GPS satellite and the vehicle based on the frequency fdi (i = 1, 2,...) And (3) the pseudo-range change amount ΔL and the relative movement amount ΔD. And are calculated.

判定基準値選択部34は、乖離量Qeを判定するための判定基準値の中から、判定基準値を選択する。本実施形態では、車両の例として軌道3を走行する列車4を挙げている。そこで、軌道3に沿って予め「区域」を定めておき、その区域別の判定基準値kを定義する判定基準値設定データ42をICメモリ26に記憶させておく。そして、判定基準値選択部34は、最後に算出した列車4の位置座標に該当する区域の判定基準値kを判定基準値設定データ42から選択する。   The determination reference value selection unit 34 selects a determination reference value from among the determination reference values for determining the deviation amount Qe. In the present embodiment, the train 4 traveling on the track 3 is cited as an example of the vehicle. Therefore, a “zone” is determined in advance along the trajectory 3, and judgment reference value setting data 42 that defines a judgment reference value k for each zone is stored in the IC memory 26. Then, the determination reference value selection unit 34 selects, from the determination reference value setting data 42, the determination reference value k of the area corresponding to the position coordinate of the train 4 calculated last.

衛星選択部36は、乖離量Qeを用いて、捕捉できるGPS衛星(候補衛星)の中から測位に使用する衛星を選択する。具体的には、候補衛星毎に得られた乖離量Qeと、判定基準値選択部34で選択した判定基準値kとを比較し、乖離量Qeが判定基準値k以上の候補衛星を選択しない。
より具体的には、本実施形態では、乖離量Qeが判定基準値k以上となる候補衛星の衛星IDを除外衛星リスト60に登録する。一方、乖離量Qeが判定基準値k未満の候補衛星の衛星IDが除外衛星リスト60に登録されていた場合には登録を抹消する。衛星選択部36は、この除外衛星リスト60に登録されている候補衛星を除いて測位に使用する衛星を選択する。
The satellite selection unit 36 selects a satellite to be used for positioning from GPS satellites (candidate satellites) that can be captured using the deviation amount Qe. Specifically, the deviation amount Qe obtained for each candidate satellite is compared with the determination reference value k selected by the determination reference value selection unit 34, and a candidate satellite having the deviation amount Qe equal to or larger than the determination reference value k is not selected. .
More specifically, in this embodiment, the satellite IDs of candidate satellites whose divergence amount Qe is equal to or greater than the determination reference value k are registered in the excluded satellite list 60. On the other hand, when the satellite ID of the candidate satellite having the deviation amount Qe less than the determination reference value k is registered in the excluded satellite list 60, the registration is deleted. The satellite selection unit 36 selects satellites to be used for positioning excluding candidate satellites registered in the excluded satellite list 60.

測位演算部38は、衛星選択部36により測位に使用する衛星として選択されたGPS衛星の衛星信号から位置座標を算出し、ICメモリ26に位置座標62を記憶する。   The positioning calculation unit 38 calculates position coordinates from the satellite signals of the GPS satellites selected as the satellites used for positioning by the satellite selection unit 36 and stores the position coordinates 62 in the IC memory 26.

ICメモリ26は、揮発性及び不揮発性のICメモリで実現できる。本実施形態のICメモリ26は、測位プログラム40と、判定基準値設定データ42と、航法メッセージ44と、擬似距離履歴データ48と、ドップラシフト周波数履歴データ50と、除外衛星リスト60と、位置座標62とを記憶する。   The IC memory 26 can be realized by a volatile and non-volatile IC memory. The IC memory 26 of the present embodiment includes a positioning program 40, determination reference value setting data 42, navigation message 44, pseudorange history data 48, Doppler shift frequency history data 50, excluded satellite list 60, position coordinates. 62 is stored.

測位プログラム40は、CPU24が読み出して実行することにより、ベースバンド処理部30、乖離量算出部32、判定基準値選択部34、衛星選択部36、測位演算部38としての機能を実現させるためのプログラムである。尚、ベースバンド処理部30をハードウェア回路として実現する構成の場合には、測位プログラム40を、乖離量算出部32、判定基準値選択部34、衛星選択部36、測位演算部38としての機能を実現させるためのプログラムとしたり、専用回路として実現することができる。   The positioning program 40 is read and executed by the CPU 24, thereby realizing functions as a baseband processing unit 30, a deviation amount calculation unit 32, a determination reference value selection unit 34, a satellite selection unit 36, and a positioning calculation unit 38. It is a program. In the case where the baseband processing unit 30 is configured as a hardware circuit, the positioning program 40 functions as a deviation amount calculation unit 32, a determination reference value selection unit 34, a satellite selection unit 36, and a positioning calculation unit 38. Can be realized as a program for realizing the above, or as a dedicated circuit.

判定基準値設定データ42は、乖離量Qeを判定する判定基準値を定義する。例えば、図4に示すように、区域421別に、衛星方位条件423と判定基準値425とを対応付けて格納する。区域421は、例えば、ある駅(α駅)からの距離や位置座標の範囲を定義する。衛星方位条件423は、車両の進行方向を基準とする方位範囲(例えば、−30°〜+30°)として定義してもよいし、東西南北を基準として定義してもよい。   The determination reference value setting data 42 defines a determination reference value for determining the deviation amount Qe. For example, as shown in FIG. 4, the satellite azimuth condition 423 and the determination reference value 425 are stored in association with each area 421. The area 421 defines, for example, a distance from a certain station (α station) and a range of position coordinates. The satellite azimuth condition 423 may be defined as an azimuth range (for example, −30 ° to + 30 °) based on the traveling direction of the vehicle, or may be defined based on east, west, south, and north.

航法メッセージ44は、受信したGPS衛星信号をデコードすることで得られるデータであり、各GPS衛星の衛星軌道情報を含む。なお、測位を開始する前の当初の航法メッセージ44、或いは、それに含まれる衛星軌道情報のみを、通信I/F28を介して外部から取得することとしてもよい。   The navigation message 44 is data obtained by decoding the received GPS satellite signal, and includes satellite orbit information of each GPS satellite. In addition, it is good also as acquiring the initial navigation message 44 before starting a positioning, or only the satellite orbit information included in it from the exterior via communication I / F28.

擬似距離履歴データ48は、GPS衛星別の擬似距離Lの履歴を格納する。例えば図5に示すように、衛星ID481と、エポックID483に対応付けられた擬似距離値485とを格納する。エポックID483には、例えば計測時刻を利用できる。   The pseudo distance history data 48 stores the history of the pseudo distance L for each GPS satellite. For example, as shown in FIG. 5, the satellite ID 481 and the pseudo distance value 485 associated with the epoch ID 483 are stored. For the epoch ID 483, for example, a measurement time can be used.

ドップラシフト周波数履歴データ50は、GPS衛星別のドップラシフト周波数fdの履歴を格納する。例えば図6に示すように、衛星ID501と、エポックID503に対応づけられたドップラシフト周波数値505とを格納する。   The Doppler shift frequency history data 50 stores the history of the Doppler shift frequency fd for each GPS satellite. For example, as shown in FIG. 6, the satellite ID 501 and the Doppler shift frequency value 505 associated with the epoch ID 503 are stored.

除外衛星リスト60は、測位に使用するのに不適と判断された衛星IDを登録する。
本実施形態では、GPS衛星別に乖離量Qeと判定基準値kとが比較され、判定基準値k以上であれば当該GPS衛星の衛星IDが登録される。判定基準値k未満であれば登録抹消される。なお、判定基準値kを超える場合に衛星IDを登録し、判定基準値k以下の場合に衛星IDの登録を抹消することとしてもよい。要は、判定基準値kが閾値となるということである。
The excluded satellite list 60 registers satellite IDs determined to be unsuitable for use in positioning.
In the present embodiment, the divergence amount Qe and the determination reference value k are compared for each GPS satellite, and if it is equal to or greater than the determination reference value k, the satellite ID of the GPS satellite is registered. If it is less than the criterion value k, the registration is deleted. Note that the satellite ID may be registered when the determination reference value k is exceeded, and the satellite ID registration may be deleted when the determination reference value k or less. The point is that the criterion value k is a threshold value.

位置座標62は、車両(本実施形態では列車4)の位置座標である。
その他、ICメモリ26にはベースバンド処理に必要な各種情報を格納することができるものとする。
The position coordinates 62 are position coordinates of the vehicle (the train 4 in this embodiment).
In addition, it is assumed that the IC memory 26 can store various information necessary for baseband processing.

[動作の説明]
次に、本実施形態における車両位置計測システム2の動作説明として、CPU24における測位関連の処理の流れに付いて説明する。尚、航法メッセージ44のデコード・更新についてはGPSを用いた公知の測位技術と同様に実現できるのでここでの説明は省略する。また、電源投入直後の、航法メッセージ44の再取得と、可視可能なGPS衛星6の位置の予測計算と、当該衛星からの衛星信号を探すサーチを行うが、それらも公知の測位技術と同様に実現できるので説明は省略する。
[Description of operation]
Next, as a description of the operation of the vehicle position measurement system 2 in the present embodiment, the flow of positioning related processing in the CPU 24 will be described. The decoding / updating of the navigation message 44 can be realized in the same manner as a known positioning technique using GPS, and therefore the description thereof is omitted here. In addition, immediately after the power is turned on, the navigation message 44 is reacquired, the position of the visible GPS satellite 6 is predicted, and a search for a satellite signal from the satellite is performed. Since it is realizable, explanation is omitted.

図7は、CPU24における一回の測位を行うための処理の流れを説明するためのフローチャートであって、車両位置計測システム2が起動している間、エポックの時間差と同じ周期で繰返し実行される。   FIG. 7 is a flowchart for explaining the flow of processing for performing positioning once in the CPU 24, and is repeatedly executed at the same cycle as the epoch time difference while the vehicle position measurement system 2 is activated. .

同処理において、CPU24は先ず、GPS衛星6の捕捉処理を実行する(ステップS2)。例えば、RF受信部22から入力されたデジタルIF信号(受信信号)とレプリカコードとの相関演算を、GPS衛星別に定められたレプリカコードを変更しながら行い、相関が取れた場合に、該当するGPS衛星の信号を捕捉したと判定する。   In the process, the CPU 24 first executes a GPS satellite 6 capturing process (step S2). For example, when the correlation between the digital IF signal (received signal) input from the RF receiver 22 and the replica code is performed while changing the replica code determined for each GPS satellite, and the correlation is obtained, the corresponding GPS It is determined that the satellite signal has been captured.

次に、CPU24は、捕捉した各GPS衛星6それぞれについてループAを実行する(ステップS4〜S36)。
ループAでは、先ず処理対象としているGPS衛星6の今回のエポックにおける擬似距離Lを算出し、擬似距離履歴データ48に格納する(ステップS8;図5参照)。次いで、処理対象のGPS衛星6の今回のエポックにおけるドップラシフト周波数fdを算出し、ドップラシフト周波数履歴データ50に格納する(ステップS10;図6参照)。
Next, the CPU 24 executes a loop A for each captured GPS satellite 6 (steps S4 to S36).
In the loop A, first, the pseudo distance L at the current epoch of the GPS satellite 6 to be processed is calculated and stored in the pseudo distance history data 48 (step S8; see FIG. 5). Next, the Doppler shift frequency fd in the current epoch of the GPS satellite 6 to be processed is calculated and stored in the Doppler shift frequency history data 50 (step S10; see FIG. 6).

次に、擬似距離変化量△Lを算出する(ステップS12)。
図2の数式Aを参照して述べるならば、ステップS8で新たに算出した擬似距離Lを第2のエポックにおける「擬似距離L2」とし、擬似距離履歴データ48に登録されている「時間的に最も近い過去のエポック」における擬似距離値485を第1のエポックにおける「擬似距離L1」として、擬似距離変化量△Lを算出する。
Next, the pseudo distance change amount ΔL is calculated (step S12).
2, the pseudo distance L newly calculated in step S8 is referred to as “pseudo distance L2” in the second epoch, and the “temporal distance” registered in the pseudo distance history data 48 is changed. The pseudo distance change ΔL is calculated by setting the pseudo distance value 485 in the “nearest past epoch” as the “pseudo distance L1” in the first epoch.

もし、一つ前のエポックにおいて擬似距離が算出できていれば、実質的に直前のエポックにおける擬似距離に基づいて擬似距離変化量△Lが算出される。仮に、一つ前のエポックにおいて跨線橋をくぐり抜けた時など、一時的に当該GPS衛星を捕捉できていなかったとしても、更にそれより過去のエポックにおいて擬似距離を算出できていれば、擬似距離変化量△Lを算出できる。しかし、コールドスタート時や長いトンネル通過などGPS衛星6を長時間に渡って捕捉できない状況から抜けた直後などでは、「時間的に最も近い過去のエポック」に該当する擬似距離485が無いケースがある。その場合はステップS12で擬似距離変化量△Lは算出できないことになる。そこで、擬似距離変化量△Lが算出できなければ(ステップS14のNO)、当該処理対象のGPS衛星6の衛星IDを除外衛星リスト60に登録し(ステップS16)、ループAを終了する(ステップS36)。   If the pseudo distance can be calculated in the immediately preceding epoch, the pseudo distance change ΔL is calculated based on the pseudo distance in the immediately preceding epoch. Even if the GPS satellite could not be captured temporarily, such as when passing through the overpass in the previous epoch, if the pseudorange can be calculated in the epoch in the past, the pseudorange change amount ΔL can be calculated. However, there is a case where there is no pseudo distance 485 corresponding to “the past epoch closest in time” immediately after exiting from a situation where the GPS satellite 6 cannot be captured for a long time, such as during a cold start or through a long tunnel. . In this case, the pseudo distance change amount ΔL cannot be calculated in step S12. Therefore, if the pseudo distance change amount ΔL cannot be calculated (NO in step S14), the satellite ID of the GPS satellite 6 to be processed is registered in the excluded satellite list 60 (step S16), and the loop A is terminated (step S16). S36).

もし、ステップS12で擬似距離変化量△Lを算出できたならば(ステップS14のYES)、CPU24は次に相対移動量△Dを算出する(ステップS20)。
図2の数式Bを参照して述べるならば、ステップS10で新たに算出したドップラシフト周波数fdを第2のエポックにおける「ドップラシフト周波数fd2」とし、ドップラシフト周波数履歴データ50に登録されている「時間的に最も近い過去のエポック」におけるドップラシフト周波数値505を第1のエポックにおける「ドップラシフト周波数fd1」として、エポック時間差△tを今エポックと時間的に最も近い過去のエポックとの時間差として、相対移動距離△Dを算出する。
If the pseudo distance change amount ΔL can be calculated in step S12 (YES in step S14), the CPU 24 then calculates the relative movement amount ΔD (step S20).
2, the Doppler shift frequency fd newly calculated in Step S10 is set to “Doppler shift frequency fd2” in the second epoch, and is registered in the Doppler shift frequency history data 50. The Doppler shift frequency value 505 in the past epoch closest in time is set as the “Doppler shift frequency fd1” in the first epoch, and the epoch time difference Δt is set as the time difference between the current epoch and the past epoch closest in time. The relative movement distance ΔD is calculated.

先ほどの擬似距離変化量△Lの算出と同様の理由で相対移動距離△Dが算出できなければ(ステップS22のNO)、CPU24は当該処理対象のGPS衛星6の衛星IDを除外衛星リスト60に登録し(ステップS16)、ループAを終了する(ステップS36)。   If the relative movement distance ΔD cannot be calculated for the same reason as the calculation of the pseudo distance change ΔL (NO in step S22), the CPU 24 adds the satellite ID of the GPS satellite 6 to be processed to the excluded satellite list 60. Registration is performed (step S16), and loop A is terminated (step S36).

もし、相対移動距離△Dが算出できたならば(ステップS22のYES)、CPU24は、次に処理対象のGPS衛星6についての乖離量Qeを算出する(ステップS24;図2の数式C参照)。   If the relative movement distance ΔD can be calculated (YES in step S22), the CPU 24 next calculates a deviation amount Qe for the GPS satellite 6 to be processed (step S24; see formula C in FIG. 2). .

次いで、判定基準値設定データ42を参照して、最後に算出した位置座標62から該当する区域421を割り出すとともに、航法メッセージ44に含まれる衛星軌道情報に基づいて処理対象のGPS衛星6の方位に該当する衛星方位条件423を選択する。そして、選択した判定基準値425を、今回の判定に使用する判定基準値kとする(ステップS26;図4参照)。
そして、この判定基準値kと乖離量Qeとを比較し、乖離量Qeが判定基準値k以上であれば(ステップS28のYES)、CPU24は当該処理対象のGPS衛星6の衛星IDを除外衛星リスト60に登録し(ステップS16)、ループAを終了する(ステップS36)。もし、乖離量Qeが判定基準値k未満であれば(ステップS28のNO)、CPU24は当該処理対象のGPS衛星6の衛星IDを除外衛星リスト60から登録抹消し(ステップS30)、ループAを終了する(ステップS36)。
Next, referring to the determination reference value setting data 42, the corresponding area 421 is determined from the position coordinate 62 calculated last, and the direction of the GPS satellite 6 to be processed is determined based on the satellite orbit information included in the navigation message 44. A corresponding satellite orientation condition 423 is selected. Then, the selected determination reference value 425 is set as a determination reference value k used for the current determination (step S26; see FIG. 4).
Then, the determination reference value k is compared with the deviation amount Qe. If the deviation amount Qe is equal to or larger than the determination reference value k (YES in step S28), the CPU 24 excludes the satellite ID of the GPS satellite 6 to be processed as an excluded satellite. It registers in the list 60 (step S16) and ends the loop A (step S36). If the deviation amount Qe is less than the determination reference value k (NO in step S28), the CPU 24 deregisters the satellite ID of the GPS satellite 6 to be processed from the excluded satellite list 60 (step S30), and loops A. The process ends (step S36).

ステップS2にて捕捉された全てのGPS衛星6についてループAを実行したならば、CPU24は、捕捉されたGPS衛星6のうち除外衛星リスト60に登録されている衛星を除外した中から、所定数の測位に使用する衛星を選択する(ステップS40)。もし、ここで測位に必要な所定数の衛星を選択できなければ(ステップS44のNO)、測位演算は行わずにステップS2に戻る。もし、所定数の衛星を選択できれば(ステップS44のYES)、それら選択した衛星を使用して測位演算を実行して今回のエポックにおける位置座標を算出して、ICメモリ26に位置座標62に記憶する(ステップS46)。   If the loop A is executed for all the GPS satellites 6 acquired in step S2, the CPU 24 excludes the satellites registered in the excluded satellite list 60 from the acquired GPS satellites 6, and then the predetermined number. A satellite to be used for positioning is selected (step S40). If a predetermined number of satellites necessary for positioning cannot be selected here (NO in step S44), the process returns to step S2 without performing positioning calculation. If a predetermined number of satellites can be selected (YES in step S44), the position calculation in the current epoch is calculated by using the selected satellites and stored in the IC memory 26 in the position coordinates 62. (Step S46).

以上、本実施形態によれば、原理的にマルチパスの影響を受け易い擬似距離変化量と、マルチパスの影響を受けにくいドップラシフト周波数に基づく相対移動量との乖離量から、マルチパスの影響を受けているかどうかを検定することが可能となる。   As described above, according to the present embodiment, the influence of the multipath is obtained from the divergence between the pseudo distance change amount that is easily affected by the multipath and the relative movement amount based on the Doppler shift frequency that is not easily affected by the multipath. It is possible to test whether or not you have received

また、本実施形態では、擬似距離とドップラシフト周波数の履歴データを保持するので、擬似距離変化量△Lを、最新のエポックにおける擬似距離と、当該エポックと最も時間的に近い過去のエポックにおける擬似距離から算出することができる。また、当該過去のエポックにおけるドップラシフト周波数を取得して相対移動量△Dを算出できる。
よって、常時は時間的に連続するエポック間での擬似距離及びドップラシフト周波数に基づいて精度良くマルチパスの影響を検定できる。仮に跨線橋をくぐり抜けた時など、一時的にGPS衛星からの電波が正しく受信できなかったエポックが単数回又は複数回連続して有ったとしても、それらよりも過去のエポックに遡って擬似距離やドップラシフト周波数を参照し、当該過去のエポックまでの時間差△tに基づいて相対移動距離△Dを算出できるので、マルチパスの影響の検定を継続的に実現することができる。
Further, in the present embodiment, since the history data of the pseudo distance and the Doppler shift frequency is held, the pseudo distance change ΔL is set to the pseudo distance in the latest epoch and the pseudo epoch in the past epoch closest to the epoch. It can be calculated from the distance. Further, the relative movement amount ΔD can be calculated by acquiring the Doppler shift frequency in the past epoch.
Therefore, it is possible to accurately test the influence of multipath based on the pseudorange and Doppler shift frequency between epochs that are continuous in time. Even if there are epochs that were temporarily unable to receive radio waves from GPS satellites singly or multiple times, such as when passing through an overpass, the pseudo- Since the relative movement distance ΔD can be calculated based on the time difference Δt up to the past epoch with reference to the Doppler shift frequency, it is possible to continuously realize the multipath effect test.

〔変形例〕
以上、本発明を適用した実施形態について説明したが、本発明の形態はこれに限定されるものではなく、適宜構成要素の追加・省略・変更を施すことができる。
[Modification]
As described above, the embodiment to which the present invention is applied has been described. However, the embodiment of the present invention is not limited to this, and additions, omissions, and changes of components can be appropriately made.

[その1]
例えば、上記実施形態では、乖離量Qeと判定基準値kとを比較することで、測位演算に使用しない衛星をふるい落としているが、判定基準値kの設定を省略することもできる。例えば、乖離量算出部32が、ICメモリ26にGPS衛星別に算出した乖離量Qeを記憶する構成とし、記憶された乖離量Qeの小さいほうから順に測位演算に使用する衛星を選択する構成も可能である。
[Part 1]
For example, in the above-described embodiment, the satellites that are not used for the positioning calculation are screened out by comparing the deviation amount Qe with the determination reference value k. However, the setting of the determination reference value k may be omitted. For example, the divergence amount calculation unit 32 may be configured to store the divergence amount Qe calculated for each GPS satellite in the IC memory 26, and may be configured to select satellites used for positioning calculation in order from the smaller stored divergence amount Qe. It is.

具体的には、乖離量算算出部32は、図8に示すように、衛星別の乖離量履歴データ64をICメモリ26に記憶する。乖離量履歴データ64は、衛星ID641別に、エポックID643と、当該エポックにおいて算出された乖離量値645とを対応付けて格納する。   Specifically, the divergence amount calculation calculation unit 32 stores the divergence amount history data 64 for each satellite in the IC memory 26 as shown in FIG. The divergence amount history data 64 stores the epoch ID 643 and the divergence amount value 645 calculated in the epoch in association with each satellite ID 641.

そして、図9に示すように、ステップS16、ステップS26〜S30に代えて、(1)ステップS24で算出した乖離量QeをループAでの処理対象GPS衛星に対応する乖離量履歴データ64に記憶するステップS32と、(2)ステップS14で否定判定の場合、及びステップS22で否定判定の場合に、ループAでの処理対象GPS衛星に対応する乖離量履歴データ64に所定値(例えばNULL)を格納するステップS34と、を実行する。   As shown in FIG. 9, instead of steps S16 and S26 to S30, (1) the deviation amount Qe calculated in step S24 is stored in the deviation amount history data 64 corresponding to the GPS satellite to be processed in loop A. (2) When a negative determination is made at step S14 and when a negative determination is made at step S22, a predetermined value (for example, NULL) is set in the deviation amount history data 64 corresponding to the processing target GPS satellite in loop A. Step S34 of storing is executed.

次いで、ループAの処理の後、最新(今回)のエポックにおける乖離量Qeが、小さい方から順に所定数(4個以上であり、例えば5個でもよいし、4個でもよい。)のGPS衛星を選択する(ステップS42)。そして、選択したGPS衛星の衛星信号を用いて位置座標を算出する(ステップS46)。
この場合、乖離量Qeは、GPS衛星を選択するか否かの指標値となる。
Next, after the processing of loop A, the GPS satellites of a predetermined number (four or more, for example, five or four) may be used in order from the smallest divergence amount Qe in the latest (current) epoch. Is selected (step S42). Then, the position coordinates are calculated using the satellite signal of the selected GPS satellite (step S46).
In this case, the deviation amount Qe is an index value indicating whether or not to select a GPS satellite.

[その2]
また、上述した実施形態では、衛星測位システムとしてGPSを例に挙げて説明したが、GPS以外の衛星測位システムに本発明を適用することも可能である。
また、車両を列車としたが、自動車に適用することも可能である。
[Part 2]
In the above-described embodiment, the GPS has been described as an example of the satellite positioning system. However, the present invention can be applied to a satellite positioning system other than the GPS.
Moreover, although the vehicle is a train, it can also be applied to an automobile.

2…車両位置計測システム
3…軌道
4…列車
6…GPS衛星
7…衛星信号
21…GPSアンテナ
22…RF受信部
24…CPU
26…ICメモリ
28…通信I/F
30…ベースバンド処理部
301…捕捉部
303…擬似距離算出部
305…ドップラシフト周波数算出部
32…乖離量算出部
34…判定基準値選択部
36…衛星選択部
38…測位演算部
40…測位プログラム
42…判定基準値設定データ
44…航法メッセージ
48…擬似距離履歴データ
50…ドップラシフト周波数履歴データ
60…除外衛星リスト
62…位置座標
64…乖離量履歴データ
2 ... Vehicle position measurement system 3 ... Track 4 ... Train 6 ... GPS satellite 7 ... Satellite signal 21 ... GPS antenna 22 ... RF receiver 24 ... CPU
26 ... IC memory 28 ... Communication I / F
DESCRIPTION OF SYMBOLS 30 ... Baseband process part 301 ... Acquisition part 303 ... Pseudo distance calculation part 305 ... Doppler shift frequency calculation part 32 ... Deviation amount calculation part 34 ... Judgment reference value selection part 36 ... Satellite selection part 38 ... Positioning calculation part 40 ... Positioning program 42 ... Determination reference value setting data 44 ... Navigation message 48 ... Pseudo-range history data 50 ... Doppler shift frequency history data 60 ... Excluded satellite list 62 ... Position coordinates 64 ... Deviation amount history data

Claims (2)

衛星から受信した衛星信号に基づいて車両の位置を計測する車両位置計測方法であって、
各候補衛星の衛星信号を捕捉する捕捉ステップと、
前記捕捉した衛星信号に係る擬似距離を算出する擬似距離算出ステップと、
前記捕捉した衛星信号に係るドップラシフト周波数を算出するドップラシフト周波数算出ステップと、
前記候補衛星それぞれについて、1)第1のエポックで算出された擬似距離と第2のエポックで算出された擬似距離の差から擬似距離変化量を算出することと、2)前記第1のエポックで算出されたドップラシフト周波数と前記第2のエポックで算出されたドップラシフト周波数の平均周波数、及び、前記第1のエポックと前記第2のエポック間の時間間隔を用いて当該候補衛星と前記車両との相対移動量を算出することと、3)前記疑似距離変化量と前記相対移動量との差である乖離量を算出することと、を実行する乖離量算出ステップと、
前記候補衛星の中から前記乖離量の小さい順に所定数の前記候補衛星を測位に使用する衛星として選択する衛星選択ステップと、
前記選択した衛星の衛星信号に基づいて測位する測位ステップと、
を含む車両位置計測方法。
A vehicle position measurement method for measuring the position of a vehicle based on a satellite signal received from a satellite,
A capture step for capturing satellite signals of each candidate satellite;
A pseudorange calculating step for calculating a pseudorange related to the captured satellite signal;
A Doppler shift frequency calculating step for calculating a Doppler shift frequency related to the captured satellite signal;
For each of the candidate satellites, 1) calculating the pseudorange variation from the difference between the pseudorange calculated in the first epoch and the pseudorange calculated in the second epoch, and 2) in the first epoch Using the calculated Doppler shift frequency, the average frequency of the Doppler shift frequency calculated in the second epoch, and the time interval between the first epoch and the second epoch, the candidate satellite and the vehicle A divergence amount calculating step of performing a calculation of a relative movement amount, and 3) calculating a divergence amount that is a difference between the pseudo distance change amount and the relative movement amount ;
A satellite selection step of selecting a predetermined number of the candidate satellites as satellites to be used for positioning from the candidate satellites in ascending order of the deviation amount ;
A positioning step of positioning based on the satellite signal of the selected satellite;
A vehicle position measuring method including:
衛星から受信した衛星信号に基づいて車両の位置を測位する車両位置計測システムであって、
各候補衛星の衛星信号を捕捉する捕捉部と、
前記捕捉した衛星信号に係る擬似距離を算出する擬似距離算出部と、
前記捕捉した衛星信号に係るドップラシフト周波数を算出するドップラシフト周波数算出部と、
前記候補衛星それぞれについて、1)第1のエポックで算出された擬似距離と第2のエポックで算出された擬似距離の差から擬似距離変化量を算出することと、2)前記第1のエポックで算出されたドップラシフト周波数と前記第2のエポックで算出されたドップラシフト周波数の平均周波数、及び、前記第1のエポックと前記第2のエポック間の時間間隔を用いて当該候補衛星と前記車両との相対移動量を算出することと、3)前記疑似距離変化量と前記相対移動量との差である乖離量を算出することと、を実行する乖離量算出部と、
前記候補衛星の中から前記乖離量の小さい順に所定数の前記候補衛星を測位に使用する衛星として選択する衛星選択部と、
前記選択した衛星の衛星信号に基づいて測位する測位部と、
を備えた車両位置計測システム。
A vehicle position measurement system that measures the position of a vehicle based on a satellite signal received from a satellite,
A capture unit for capturing satellite signals of each candidate satellite;
A pseudorange calculator for calculating a pseudorange related to the captured satellite signal;
A Doppler shift frequency calculation unit for calculating a Doppler shift frequency related to the captured satellite signal;
For each of the candidate satellites, 1) calculating the pseudorange variation from the difference between the pseudorange calculated in the first epoch and the pseudorange calculated in the second epoch, and 2) in the first epoch Using the calculated Doppler shift frequency, the average frequency of the Doppler shift frequency calculated in the second epoch, and the time interval between the first epoch and the second epoch, the candidate satellite and the vehicle A divergence amount calculation unit that executes the following: 3) calculating a relative movement amount; and 3) calculating a divergence amount that is a difference between the pseudo distance change amount and the relative movement amount ;
A satellite selecting unit that selects a satellite to be used for positioning a predetermined number of said candidate satellites in the ascending order of the deviation amount from among the candidate satellites,
A positioning unit for positioning based on the satellite signal of the selected satellite;
A vehicle position measurement system.
JP2013020572A 2013-02-05 2013-02-05 Vehicle position measuring method and vehicle position measuring system Expired - Fee Related JP6017983B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013020572A JP6017983B2 (en) 2013-02-05 2013-02-05 Vehicle position measuring method and vehicle position measuring system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013020572A JP6017983B2 (en) 2013-02-05 2013-02-05 Vehicle position measuring method and vehicle position measuring system

Publications (2)

Publication Number Publication Date
JP2014153084A JP2014153084A (en) 2014-08-25
JP6017983B2 true JP6017983B2 (en) 2016-11-02

Family

ID=51575104

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013020572A Expired - Fee Related JP6017983B2 (en) 2013-02-05 2013-02-05 Vehicle position measuring method and vehicle position measuring system

Country Status (1)

Country Link
JP (1) JP6017983B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103148813B (en) * 2013-01-31 2015-09-09 湖南致力工程科技有限公司 For the treatment of the method for GPS deformation measurement data
JP6017984B2 (en) * 2013-02-05 2016-11-02 公益財団法人鉄道総合技術研究所 Vehicle position measuring method and vehicle position measuring system
JP2016188792A (en) * 2015-03-30 2016-11-04 独立行政法人交通安全環境研究所 Position measuring method and position measuring system
JP6946660B2 (en) * 2017-02-23 2021-10-06 株式会社デンソー Positioning device

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3811566B2 (en) * 1998-02-06 2006-08-23 関西電力株式会社 Moving body positioning method, moving body positioning apparatus, and unmanned moving body guidance apparatus
JP2001124840A (en) * 1999-10-28 2001-05-11 Matsushita Electric Ind Co Ltd Gps receiver and positioning method
JP2006133142A (en) * 2004-11-08 2006-05-25 Furuno Electric Co Ltd Positioning receiver
JP2006242911A (en) * 2005-03-07 2006-09-14 Denso Corp Position detector
JP2008145183A (en) * 2006-12-07 2008-06-26 Mitsubishi Electric Corp Measurement device, multipath determination method, and measurement method
JP5094344B2 (en) * 2007-11-20 2012-12-12 古野電気株式会社 Abnormal satellite detection device and positioning device
US7821454B2 (en) * 2007-11-20 2010-10-26 Sirf Technology, Inc. Systems and methods for detecting GPS measurement errors
JP5244490B2 (en) * 2008-07-28 2013-07-24 独立行政法人交通安全環境研究所 GPS positioning apparatus and method therefor
JP2010249620A (en) * 2009-04-15 2010-11-04 Japan Radio Co Ltd Positioning device
JP5413118B2 (en) * 2009-10-09 2014-02-12 トヨタ自動車株式会社 Positioning system
WO2011105445A1 (en) * 2010-02-26 2011-09-01 古野電気株式会社 Pseudo range estimation method, pseudo range estimation program, gnss receiver apparatus, and mobile terminal
US9121936B2 (en) * 2010-02-26 2015-09-01 Furuno Electric Company Limited Positioning method, GNSS receiving apparatus, and mobile terminal
WO2011105447A1 (en) * 2010-02-26 2011-09-01 古野電気株式会社 Multipath detection method, multipath detection program, gnss receiver apparatus, and mobile terminal
JP6017984B2 (en) * 2013-02-05 2016-11-02 公益財団法人鉄道総合技術研究所 Vehicle position measuring method and vehicle position measuring system

Also Published As

Publication number Publication date
JP2014153084A (en) 2014-08-25

Similar Documents

Publication Publication Date Title
CN102713674B (en) Centimeter positioning using low cost single frequency GNSS receivers
JP4655139B2 (en) Mobile positioning device
US10739471B2 (en) GNSS receiver with a capability to resolve ambiguities using an uncombined formulation
US11099276B2 (en) Anti-spoofing system for GNSS receivers
JP6017984B2 (en) Vehicle position measuring method and vehicle position measuring system
KR20100021382A (en) Method and system for calibrating group delay errors in a combined gps and glonass receiver
US9253751B2 (en) Method of calculating movement speed and device for calculating movement speed
JP6017983B2 (en) Vehicle position measuring method and vehicle position measuring system
US20130016007A1 (en) Gnss receiver and positioning method
JP2009133702A (en) Positioning device
JP5636552B2 (en) How to optimize the acquisition of spread spectrum signals from satellites by mobile receivers
KR20130111315A (en) A navigation bit boundary determination apparatus and a method thereof
JP2008232771A (en) Positioning device
JP5163474B2 (en) GNSS receiver and positioning method
JP2010164340A (en) Gnss-receiving device and positioning method
JP2009025045A (en) Positioning device for moving body
JP2006250681A (en) Setting method of search frequency in gps receiver, gps receiver and on-vehicle unit
JP2019168257A (en) Moving body information estimation device and program
JP2010145178A (en) Moving body position specification device
JP2010112759A (en) Mobile body positioning apparatus
JP2008089309A (en) Positional information acquisition device, positional information acquisition method, and program for terminal device
JP4738944B2 (en) GPS receiver
JP2001311768A (en) Gps receiver with multi-path judging function
JP2010139355A (en) Positioning method using gps
US10816675B2 (en) Coordinate output method and coordinate output device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150508

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160303

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160315

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160516

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160920

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160929

R150 Certificate of patent or registration of utility model

Ref document number: 6017983

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees