JP2000310538A - Position detecting apparatus - Google Patents

Position detecting apparatus

Info

Publication number
JP2000310538A
JP2000310538A JP11986199A JP11986199A JP2000310538A JP 2000310538 A JP2000310538 A JP 2000310538A JP 11986199 A JP11986199 A JP 11986199A JP 11986199 A JP11986199 A JP 11986199A JP 2000310538 A JP2000310538 A JP 2000310538A
Authority
JP
Japan
Prior art keywords
position data
positioning
speed
estimated
current position
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP11986199A
Other languages
Japanese (ja)
Other versions
JP3601354B2 (en
Inventor
Takayuki Kamiya
隆行 神谷
Hiroyuki Kitagawa
弘之 北川
Yoshitaka Ozaki
義隆 尾崎
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.)
Denso Corp
Original Assignee
Denso 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 Denso Corp filed Critical Denso Corp
Priority to JP11986199A priority Critical patent/JP3601354B2/en
Publication of JP2000310538A publication Critical patent/JP2000310538A/en
Application granted granted Critical
Publication of JP3601354B2 publication Critical patent/JP3601354B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Abstract

PROBLEM TO BE SOLVED: To enhance the detecting accuracy during the stop of a position detecting apparatus which makes use of an artificial satellite. SOLUTION: At every computing timing, positioning position data Pg(i) which are found on the basis of received radio waves from a GPS satellite, a speed vector Vg(i) and a speed-error estimated value Vth are read out (S110). Estimated position data Pe(i) which is found at a previous computing timing and the positioning position data Pg(i) are weighted and averaged, and present- position data P(i) is found (S160). The speed vector Vg(i) is added to the present- position data P(i), and estimated position data Pe(i+1) at a next computing timing is found (S190). However, when the magnitude Vg(i) of the speed vector is smaller than the speed-error estimated value Vth (S180-Y), the estimated position data Pe(i) at this time is set as the estimated position data Pe(i+1) in the next time as it is.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、人工衛星から電波
を受信して、その受信地点の位置を検出する位置検出装
置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a position detecting device for receiving a radio wave from an artificial satellite and detecting the position of the receiving point.

【0002】[0002]

【従来の技術】従来より、この種の装置として、GPS
(Global Positioning System) 用に設けられた複数の
人工衛星(以下、GPS衛星という)を利用して、3或
いは4個のGPS衛星から同時に電波を受信し、電波に
重畳されたC/Aコードを用いて各GPS衛星との擬似
距離を求め、この擬似距離に基づき測位計算を行って現
在位置を求めるGPS受信装置が知られている。
2. Description of the Related Art Conventionally, as this kind of device, GPS
(Global Positioning System) Using a plurality of artificial satellites (hereinafter referred to as “GPS satellites”) provided for the purpose of receiving radio waves simultaneously from three or four GPS satellites and transmitting a C / A code superimposed on the radio waves 2. Description of the Related Art There is known a GPS receiving apparatus that calculates a pseudo distance with each GPS satellite using the calculated distance and calculates a positioning based on the pseudo distance to obtain a current position.

【0003】なお、擬似距離は、GPS衛星から決めら
れたタイミングで送信されるC/AコードがGPS受信
装置に到達するまでに要する時間を、GPSの基準時間
に従って動作する時計にて測定し、この測定時間に電波
の伝搬速度を乗じることで簡単に求めることができる。
[0003] The pseudorange is obtained by measuring the time required for a C / A code transmitted from a GPS satellite at a predetermined timing to reach a GPS receiver by a clock operating according to a GPS reference time. The measurement time can be easily obtained by multiplying the propagation time of the radio wave.

【0004】ところで、GPSでは、検出精度をコント
ロールするために、精度劣化(SA:Selective Availa
bility)が導入されており、意図的に、搬送波の周波数
等がある範囲内でゆらぐようにされているため、このゆ
らぎに応じて各種測定誤差が生じ、ひいては測位計算に
よって求めた測位位置がばらついてしまう。
Incidentally, in GPS, in order to control the detection accuracy, accuracy deterioration (SA: Selective Availa
abilities) are introduced, and the frequency of the carrier wave is intentionally fluctuated within a certain range. Therefore, various measurement errors occur according to the fluctuations, and the positioning position obtained by the positioning calculation varies. Would.

【0005】このような測位位置のばらつきを減少させ
る方法の一つとして、最新の測位位置を過去の測位位置
により平均化(例えば、前回の測位位置と今回の測位位
置により平均化)することが考えられるが、この方法を
適用すると、平均化により求められた測位位置が実際の
位置より時間的に遅れたものとなってしまうという問題
があった。
[0005] As one of the methods for reducing the variation of the positioning positions, the latest positioning position is averaged by past positioning positions (for example, by averaging the previous positioning position and the current positioning position). However, when this method is applied, there is a problem that the measured position obtained by averaging is delayed in time from the actual position.

【0006】これに対して、例えば、特開平8−686
51号公報には、従来通りC/Aコードに基づいて算出
される擬似距離から測位位置を求めるだけでなく、GP
S衛星との相対速度に応じて搬送波に生じるドップラシ
フトを利用して、GPS受信装置(又は該装置を搭載し
た移動体)の移動速度及び移動方向を示す速度ベクトル
を求め、この速度ベクトルを用いることにより、時間遅
れを生じることなく、GPSによる測位位置のバラツキ
を抑える装置が開示されている。
On the other hand, for example, Japanese Patent Application Laid-Open No. 8-686
No. 51 discloses not only a method of calculating a positioning position from a pseudo distance calculated based on a C / A code as in the past, but also a GP.
Using a Doppler shift generated in a carrier wave in accordance with a relative speed with respect to the S satellite, a speed vector indicating a moving speed and a moving direction of the GPS receiving device (or a moving body equipped with the device) is obtained, and this speed vector is used. Thus, a device that suppresses variations in the positioning position by GPS without causing a time delay is disclosed.

【0007】この装置では、既に算出されている現在位
置と速度ベクトルとに基づいて次回の測定タイミングで
の推定位置を求め、次回の測定タイミングではこの推定
位置と実測された測位位置との加重平均値を現在位置と
して算出するようにされている。
In this device, an estimated position at the next measurement timing is obtained based on the already calculated current position and the velocity vector, and a weighted average of the estimated position and the actually measured positioning position is obtained at the next measurement timing. The value is calculated as the current position.

【0008】つまり、過去(前回)の測位位置を、その
まま用いるのではなく、推定位置を算出するために用
い、この推定位置と実測された測位位置との間で平均化
を行っているため、従来装置のような時間遅れを生じる
ことがなく、また、GPS受信機内では搬送波の方が、
ビットレートが1MbpsのC/Aコードより分解能が
高いため、搬送波を用いて求めた速度ベクトルを組み合
わせて現在位置を求めた場合、C/Aコードを用いて求
めた擬似距離だけで現在位置を求めるよりも、バラツキ
を抑えることができ、測定精度が向上するのである。
In other words, the past (previous) measured position is not used as it is, but is used for calculating an estimated position, and averaging is performed between the estimated position and the actually measured position. There is no time lag unlike the conventional device, and the carrier wave in the GPS receiver is
Since the bit rate has a higher resolution than the C / A code of 1 Mbps, when the current position is obtained by combining the velocity vectors obtained by using the carrier, the current position is obtained only by the pseudo distance obtained by using the C / A code. Rather, the variation can be suppressed, and the measurement accuracy is improved.

【0009】[0009]

【発明が解決しようとする課題】しかし、上述のSAに
よる搬送波周波数のゆらぎは、ドップラシフトによる周
波数変位と区別できないため、速度ベクトルには、この
ゆらぎの大きさに応じた誤差成分が重畳されてしまう。
However, since the carrier frequency fluctuation due to SA cannot be distinguished from the frequency displacement due to Doppler shift, an error component corresponding to the magnitude of the fluctuation is superimposed on the velocity vector. I will.

【0010】そして、高速移動時には、その移動により
生じる速度ベクトル成分が十分に大きいため、誤差成分
の影響は小さいが、停止時や低速走行時には、移動によ
り生じる速度ベクトル成分が小さくなり、誤差成分が相
対的に大きくなるため、その影響が大きくなる。
[0010] When moving at high speed, the speed vector component generated by the movement is sufficiently large, so that the effect of the error component is small. The effect is greater because the distance is relatively large.

【0011】その結果、上述のGPS受信機を、車両用
ナビゲーション装置に適用している場合には、実際には
停車しているにも関わらず、速度ベクトルが検出されて
しまうため、地図上で現在位置が常時変化し、表示が見
難いという問題があった。また、例えば、交差点にて一
旦停車してから再度動き始める場合、誤差成分によって
停車中の現在位置がばらつくと、停車した位置とは異な
った位置から動き始めることになるため、図5に示すよ
うに、クランクした道路を走行しているような表示とな
る。その結果、GPS受信機にて求めた現在位置と、地
図上の道路とのマッチングを行っていると、場合によっ
ては、間違った道路上に現在位置が表示されてしまうと
いう問題があった。
As a result, when the above-mentioned GPS receiver is applied to a navigation device for a vehicle, a speed vector is detected in spite of the fact that the vehicle is actually stopped, so that the speed vector is detected on a map. There is a problem that the current position constantly changes and the display is difficult to see. Also, for example, when the vehicle stops once at an intersection and starts moving again, if the current position during the stop varies due to an error component, the vehicle will start moving from a position different from the stopped position, as shown in FIG. Is displayed as if the vehicle is running on a cranked road. As a result, if the current position obtained by the GPS receiver is matched with the road on the map, there is a problem that the current position is displayed on an incorrect road in some cases.

【0012】更に、特に市街地の走行時には、建物等に
反射した電波を受信してしまうマルチパスが生じる。こ
のマルチパスが生じた場合、検出される擬似距離が長く
なるため、正しい測位が行われなくなる。但し、走行中
であれば、周囲の状況が短時間の間に変化するため、マ
ルチパス状態の継続時間も短く、その影響はほとんどな
いが、信号で一時停止したり、渋滞等により低速走行し
ている場合には、マルチパス状態が継続するため、その
影響は大きく、ナビゲーション装置の地図上に表示され
る現在位置が大きく誤ったものとなってしまうという問
題があった。
Further, especially when traveling in an urban area, a multipath occurs in which radio waves reflected on buildings and the like are received. When this multipath occurs, the pseudorange to be detected becomes longer, so that correct positioning cannot be performed. However, if the vehicle is traveling, the surrounding conditions change in a short period of time, so the duration of the multipath state is short, and there is almost no effect. In such a case, since the multipath state continues, the influence is large, and there is a problem that the current position displayed on the map of the navigation device is greatly incorrect.

【0013】そこで本発明は、上記問題点を解決するた
めに、人工衛星を利用した位置検出装置において、停止
中に検出される位置の精度を向上させることを目的とす
る。
Accordingly, an object of the present invention is to improve the accuracy of a position detected during stoppage in a position detecting device using an artificial satellite in order to solve the above problems.

【0014】[0014]

【課題を解決するための手段】上記目的を達成するため
になされた発明である請求項1に記載の位置検出装置で
は、受信手段が、複数の人工衛星から電波を受信し、そ
の受信電波に重畳されたタイミング情報に基づいて、擬
似距離検出手段が、人工衛星のそれぞれとの擬似距離を
検出し、更に測位演算手段が、検出された擬似距離に基
づき、予め設定された演算タイミング毎に測位位置を繰
り返し求める。
According to a first aspect of the present invention, there is provided a position detecting apparatus for receiving radio waves from a plurality of artificial satellites, and receiving the received radio waves from the plurality of artificial satellites. Based on the superimposed timing information, the pseudo-range detecting means detects a pseudo-range with each of the artificial satellites, and further, the positioning calculating means determines the positioning at each predetermined calculation timing based on the detected pseudo-distance. Find the position repeatedly.

【0015】これと共に、相対速度検出手段が、受信手
段が受信する電波の周波数変位に基づいて、当該装置と
人工衛星のそれぞれとの相対速度を検出し、速度ベクト
ル演算手段が、検出された相対速度に基づいて、演算タ
イミング毎に当該装置の移動速度及び移動方向を示す速
度ベクトルを繰り返し求める。即ち、電波の周波数は、
当該装置と人工衛星との相対速度に応じてドップラシフ
トするため、このドップラシフトによって生じた周波数
変位に基づいて、相対速度を求めることができるのであ
る。
At the same time, the relative velocity detecting means detects the relative velocity between the device and each of the artificial satellites based on the frequency displacement of the radio wave received by the receiving means, and the velocity vector calculating means detects the relative velocity. Based on the speed, a speed vector indicating the moving speed and moving direction of the device is repeatedly obtained for each calculation timing. That is, the frequency of the radio wave is
Since the Doppler shift is performed according to the relative speed between the device and the artificial satellite, the relative speed can be obtained based on the frequency displacement caused by the Doppler shift.

【0016】そして、位置推定手段が、速度ベクトル演
算手段にて求められた速度ベクトル、及び予め設定され
た当該装置の現在位置に基づいて、次回の演算タイミン
グでの当該装置の推定位置を求め、現在位置更新手段
が、測位演算手段にて求められた測位位置と、前回の演
算タイミング時に位置推定手段にて求められた推定位置
との加重平均値により、現在位置を更新する。
The position estimating means obtains an estimated position of the device at the next calculation timing based on the speed vector obtained by the speed vector calculating means and the preset current position of the device. The current position updating means updates the current position with a weighted average value of the positioning position obtained by the positioning calculating means and the estimated position obtained by the position estimating means at the previous calculation timing.

【0017】但し、停止判定手段にて停止状態と判定さ
れた場合、禁止手段が、位置推定手段による推定位置の
更新を禁止して、今回の演算タイミングにて使用された
推定位置を、そのまま次回の演算タイミングでの推定位
置として設定する。このように構成された本発明の位置
検出装置によれば、タイミング情報を用いて求めた測位
位置と、速度ベクトルを用いて求めた推定位置とに基づ
いて現在位置を求めているので、過去の測位位置を用い
てばらつきを抑える従来装置のような時間遅れを生じさ
せることなく現在位置の精度を向上させることができ
る。
However, if the stop determining means determines that the vehicle is in the stop state, the prohibiting means prohibits the update of the estimated position by the position estimating means, and replaces the estimated position used at the current calculation timing with the next time. Is set as the estimated position at the calculation timing. According to the position detection device of the present invention configured as described above, the current position is obtained based on the measured position obtained using the timing information and the estimated position obtained using the velocity vector. The accuracy of the current position can be improved without causing a time delay as in a conventional device that suppresses variations using the measured position.

【0018】これに加えて、本発明の位置検出装置によ
れば、当該装置(或いは当該装置を搭載する移動体)の
移動速度に基づいて停止状態を判定しており、停止状態
にある時には、速度ベクトルによる推定位置の更新を行
わないため、停止時や低速走行時に相対的に大きくなる
速度ベクトルの誤差成分によって、現在位置がばらつい
てしまうことを確実に抑えることができる。
In addition, according to the position detecting device of the present invention, the stop state is determined based on the moving speed of the device (or the moving body on which the device is mounted). Since the estimated position is not updated based on the speed vector, it is possible to reliably prevent the current position from fluctuating due to an error component of the speed vector that becomes relatively large when the vehicle is stopped or running at low speed.

【0019】なお、速度ベクトル演算手段は、前回の演
算タイミングで求めた速度ベクトルと今回の演算タイミ
ングで求めた速度ベクトルとを平均化処理したものを、
位置推定手段に供給する速度ベクトルとしてもよい。こ
の場合、速度ベクトルを平均化処理した場合、速度ベク
トルの加速度が考慮された値となるため、より正確に推
定位置を求めるとができ、ひいては、より高精度に現在
位置を求めることができる。
The speed vector calculating means averages the speed vector obtained at the previous calculation timing and the speed vector obtained at the current calculation timing,
The velocity vector may be supplied to the position estimating means. In this case, when the speed vector is averaged, the value takes into account the acceleration of the speed vector, so that the estimated position can be obtained more accurately, and the current position can be obtained with higher accuracy.

【0020】なお、停止判定手段にて停止状態を判定す
る場合、具体的には、例えば請求項2記載のように、速
度ベクトル演算手段にて求められる速度ベクトルの大き
さを、予め設定された判定しきい値と比較することによ
り判定することができる。この場合、速度ベクトルの大
きさが判定しきい値より小さければ、停止状態と判定す
ればよい。
When the stop state is judged by the stop judging means, specifically, the magnitude of the speed vector obtained by the speed vector calculating means is set in advance. The determination can be made by comparing with a determination threshold value. In this case, if the magnitude of the speed vector is smaller than the determination threshold, it may be determined that the vehicle is in the stopped state.

【0021】ところで、渋滞等により速度ベクトルの大
きさが判定しきい値以下となるような低速走行を行って
いる場合、実際には移動しているにも関わらず、推定位
置の更新が行われないため、推定位置が実際の位置から
大きく異なったものとなり、その結果、推定位置と測位
位置とに基づいて算出される現在位置の精度が低下して
しまう可能性がある。
When the vehicle is traveling at a low speed such that the magnitude of the speed vector is equal to or smaller than the determination threshold value due to traffic congestion or the like, the estimated position is updated even though the vehicle is actually moving. Therefore, the estimated position is greatly different from the actual position, and as a result, the accuracy of the current position calculated based on the estimated position and the measured position may be reduced.

【0022】そこで、請求項3記載の位置検出装置で
は、停止状態中に、測位位置と現在位置とが予め設定さ
れた上限距離以上離れている状態が、予め設定された上
限時間以上継続した場合、推定位置初期化手段が、測位
位置に基づいて推定位置を初期化する。
Therefore, in the position detecting device according to a third aspect of the present invention, when the measured position and the current position are separated from each other by a predetermined upper limit distance during the stop state for more than a predetermined upper limit time, The estimated position initialization means initializes the estimated position based on the measured position.

【0023】つまり、当該装置(当該装置を搭載した移
動体)が判定しきい値以下の低速で移動していれば、測
位位置と現在位置とが上限距離以上離れている状態が一
旦検出されると、判定しきい値以上での移動が開始され
ない限り、この状態が継続する。このため、上限時間以
上継続した場合には、当該装置は、実際には低速移動し
ているものとして、測位位置に基づいた推定位置の初期
化を行うことにより、低速移動を停止状態とみなしたこ
とにより生じた誤差を解消しているのである。
That is, if the device (moving body on which the device is mounted) is moving at a low speed equal to or less than the determination threshold, a state where the measured position is separated from the current position by the upper limit distance is once detected. This state continues unless movement at or above the determination threshold is started. For this reason, when continuing for more than the upper limit time, the apparatus considers that the low-speed movement is in the stopped state by initializing the estimated position based on the positioning position, assuming that the apparatus is actually moving at a low speed. The error caused by this is eliminated.

【0024】また、停止時に、マルチパスが生じた場
合、同様に推定位置と測位位置とが大きく異なったもの
となるが、ある程度時間が経過すると、衛星の移動によ
り受信状態が変化して、マルチパスが解消されることが
多い。このため、マルチパスが解消される可能性が高く
なるような十分な長さに上限時間を設定しておけば、こ
のマルチパスの影響が現在位置に蓄積されることがな
く、移動を再開した時には、マルチパスの影響が排除さ
れた精度のよい現在位置を得ることができる。
When a multipath occurs at the time of stoppage, the estimated position and the positioning position are also greatly different from each other. Paths are often eliminated. For this reason, if the upper limit time is set to a length long enough to increase the possibility that the multipath will be eliminated, the movement is restarted without the effect of the multipath being accumulated at the current position. At times, an accurate current position from which the influence of multipath has been eliminated can be obtained.

【0025】但し、衛星からの電波を正常に受信してい
る場合でも、個々に求められる測位位置は誤差が比較的
大きいため、請求項4記載のように、推定位置初期化手
段が推定位置の初期化を行う場合には、停止状態中に求
められた複数の測位位置を平均化して得られる位置を用
いることが望ましい。
However, even if the radio waves from the satellites are normally received, the positioning positions individually obtained have relatively large errors, so that the estimated position initializing means sets the estimated position as the estimated position. When performing initialization, it is desirable to use a position obtained by averaging a plurality of positioning positions obtained during the stop state.

【0026】なお、平均化とは、一定期間内の測位位置
を単純平均してもよいし、時系列で重み付け平均しても
よいし、その時々の測位精度に応じて重み付け平均して
もよい。ところで、測位位置の精度は、測位に使用する
人工衛星の組合せや、人工衛星自体が持つ精度のばらつ
きによって様々に変化し、これに応じて速度ベクトルの
誤差がばらつく範囲も様々に変化する。また、当該装置
が停止状態にあるか否かを判定する判定しきい値は、速
度ベクトルの誤差のばらつきより小さく設定しても、有
効な判定を行うことはできない。
The averaging may be a simple averaging of the positioning positions within a certain period, a weighted averaging in a time series, or a weighted averaging according to the positioning accuracy at each time. . By the way, the accuracy of the positioning position changes variously depending on the combination of artificial satellites used for positioning and the variation in accuracy of the artificial satellite itself, and accordingly, the range in which the error of the velocity vector varies also changes. Further, even if the determination threshold value for determining whether or not the device is in the stopped state is set to be smaller than the variation in the error of the speed vector, an effective determination cannot be made.

【0027】そこで、請求項5記載のように、判定しき
い値設定手段を設け、受信手段が受信する電波に重畳さ
れた精度情報に基づいて、停止判定手段が用いる判定し
きい値を可変設定することが望ましい。これにより、測
位に使用する人工衛星の組合せや配置により、刻々と変
化する誤差に応じて、常に最適な判定しきい値を設定で
きるため、システムの性能を最大限に引き出すことがで
きる。
Therefore, a determination threshold value setting means is provided, and the determination threshold value used by the stop determination means is variably set based on the accuracy information superimposed on the radio wave received by the receiving means. It is desirable to do. Thus, the optimum determination threshold can always be set according to the error that changes every moment depending on the combination and arrangement of the satellites used for positioning, so that the system performance can be maximized.

【0028】但し、判定しきい値を可変設定した場合、
この値が過大になると、移動を開始しているにも関わら
ず、ある程度大きな速度にならなければ、速度ベクトル
を用いた推定位置の更新が行われないため、推定位置の
更新が開始されると、現在位置が急激に変化し、これに
基づいてナビゲーション装置などの地図上に現在位置の
表示を行った場合には、現在位置の表示が不自然な動き
をすることになる。
However, when the judgment threshold is variably set,
If this value is excessive, the update of the estimated position using the speed vector is not performed unless the speed is increased to a certain degree even though the movement is started, so that the update of the estimated position is started. If the current position suddenly changes and the current position is displayed on a map such as a navigation device based on the sudden change, the display of the current position moves unnaturally.

【0029】そこで、請求項6記載のように、判定しき
い値設定手段は、制度情報を所定の演算式に適用して判
定しきい値を算出する場合、算出値が予め設定された上
限値より大きければ、判定しきい値を上限値に制限する
ことが望ましい。
According to a sixth aspect of the present invention, when the determination threshold value is calculated by applying the accuracy information to a predetermined arithmetic expression, the determination threshold value setting means sets the calculated value to a predetermined upper limit value. If it is larger, it is desirable to limit the determination threshold to the upper limit.

【0030】[0030]

【発明の実施の形態】以下に本発明の実施例を図面と共
に説明する。図1は、本実施例のGPS受信機2を用い
て構成された車両用ナビゲーション装置の全体構成を表
すブロック図である。
Embodiments of the present invention will be described below with reference to the drawings. FIG. 1 is a block diagram illustrating an overall configuration of a vehicular navigation device configured using the GPS receiver 2 of the present embodiment.

【0031】図1に示すように、車両用ナビゲーション
装置は、GPS衛星からの電波を受信して車両の現在位
置を示す現在位置データPを出力するGPS受信機2
と、道路地図の地図データを記憶する地図データ記憶媒
体4と、この地図データ記憶媒体4からの地図データお
よびGPS受信機2からの現在位置データPに基づい
て、車両の現在位置を中心にした車両の走行領域の道路
地図を表示部6に表示させるナビ演算部8とを備えてい
る。
As shown in FIG. 1, the vehicular navigation device receives a radio wave from a GPS satellite and outputs a current position data P indicating the current position of the vehicle.
And a map data storage medium 4 for storing map data of a road map, and a current position of the vehicle centered on the basis of the map data from the map data storage medium 4 and the current position data P from the GPS receiver 2. The navigation unit 8 includes a navigation unit 8 that displays a road map of the traveling area of the vehicle on the display unit 6.

【0032】ここで、本実施例のGPS受信機2は、G
PS衛星からの電波を受波するアンテナ10と、アンテ
ナ10からの受信信号に基づき、衛星信号の捕捉/追跡
の他、搬送波の抽出,C/Aコードの復調,C/Aコー
ドにてスペクトル拡散されたデータの復号,搬送波の周
波数変位量の抽出,C/Aコードの伝達時間の測定等を
行う受信回路部12と、CPU,ROM,RAMを中心
に構成されたマイクロコンピュータからなる演算部14
とを備えている。
Here, the GPS receiver 2 of the present embodiment
An antenna 10 for receiving a radio wave from a PS satellite, and based on a signal received from the antenna 10, capture / track a satellite signal, extract a carrier, demodulate a C / A code, and spread a spectrum using a C / A code. A receiving circuit unit 12 for decoding decoded data, extracting a frequency displacement amount of a carrier wave, measuring a transmission time of a C / A code, and an arithmetic unit 14 including a microcomputer mainly including a CPU, a ROM, and a RAM.
And

【0033】また、図1には、演算部14の機能ブロッ
クが示されており、図示されているように、演算部14
は、受信回路部12にて復号されたデータから精度情報
を抽出し、この精度情報に基づいて判定しきい値として
の速度誤差見積値Vthを算出する速度誤差推定部20
と、受信回路部12にて各GPS衛星毎に測定されたC
/Aコードの伝達時間に基づき、当該ナビゲーション装
置を搭載した車両と各GPS衛星との擬似距離を求め、
この擬似距離に基づいての車両の測位位置を表す測位位
置データPgを算出する測位演算部22と、受信回路部
12にて抽出された搬送波の周波数変位量に基づき、測
位演算に用いる各GPS衛星と車両との相対速度を求
め、この相対速度に基いて車両の移動速度及び移動方向
を表す速度ベクトルVgを求める速度ベクトル演算部2
4と、測位演算部22及び速度ベクトル演算部24にて
求められた測位位置データPg及び速度ベクトルVgに
基づき、現在位置データPを算出するフィルタリング部
26とからなる。
FIG. 1 shows functional blocks of the operation unit 14, and as shown in FIG.
A speed error estimating unit 20 that extracts accuracy information from data decoded by the receiving circuit unit 12 and calculates an estimated speed error value Vth as a determination threshold based on the accuracy information.
And C measured by the receiving circuit unit 12 for each GPS satellite.
The pseudo distance between the vehicle equipped with the navigation device and each GPS satellite is determined based on the transmission time of the / A code,
A positioning operation unit 22 that calculates positioning position data Pg representing a positioning position of the vehicle based on the pseudo distance, and each GPS satellite used for positioning operation based on the frequency displacement amount of the carrier extracted by the receiving circuit unit 12 Vector calculation unit 2 for obtaining a relative speed between the vehicle and the vehicle and obtaining a speed vector Vg representing the moving speed and the moving direction of the vehicle based on the relative speed.
4 and a filtering unit 26 that calculates the current position data P based on the positioning position data Pg and the speed vector Vg obtained by the positioning calculation unit 22 and the speed vector calculation unit 24.

【0034】このうち、速度誤差推定部20では、復号
データから抽出する精度情報として詳細軌道情報(エフ
ェメリス)に含まれるURA(User Range Accuracy)
を用いると共に、測位演算に用いるGPS衛星の配置に
よって決まる精度劣化係数(DOP:Dilution Of Prec
ision) を用いて、次の(1)式に従って速度誤差見積
値Vth[m/s]を算出する。
The speed error estimator 20 includes a user range accuracy (URA) included in detailed trajectory information (ephemeris) as accuracy information extracted from the decoded data.
And the accuracy degradation coefficient (DOP: Dilution Of Prec) determined by the arrangement of GPS satellites used for positioning calculation.
ision), the speed error estimated value Vth [m / s] is calculated according to the following equation (1).

【0035】 Vth=k×PMura ×PMhdop (1) なお、PMura は、URAに基づいて設定されるパラメ
ータであり、[表1]に示すように、URAに対応した
値域の平均値、例えば、URA=7であればPMura =
36が用いられる。また、PMhdopは、DOPから求め
られる水平方向の精度劣化係数(HDOP)であり、更
にkは、GPS受信機2の内部雑音によって発生するド
ップラーシフトの測定誤差を速度換算したものである。
Vth = k × PMura × PMhdop (1) Here, PMura is a parameter set based on URA, and as shown in [Table 1], an average value of a range corresponding to URA, for example, URA = 7 if PMura =
36 are used. PMhdop is a horizontal accuracy degradation coefficient (HDOP) obtained from the DOP, and k is a value obtained by converting the measurement error of the Doppler shift generated by the internal noise of the GPS receiver 2 into speed.

【0036】[0036]

【表1】 [Table 1]

【0037】但し、速度誤差推定部20は、(1)式に
て算出される速度誤差見積値Vthが、上限値(本実施例
では1[m/s])より大きい場合には、速度誤差見積
値Vthを、この上限値に制限(即ち、0<Vth≦1とな
る)し、また、測位演算に用いるGPS衛星の組合せや
配置が変化する毎に、速度誤差見積値Vthを繰り返し算
出するように構成されている。
However, if the estimated speed error Vth calculated by the equation (1) is larger than the upper limit value (1 [m / s] in this embodiment), the speed error estimating unit 20 determines that the speed error The estimated value Vth is limited to this upper limit (that is, 0 <Vth ≦ 1), and the speed error estimated value Vth is repeatedly calculated every time the combination or arrangement of the GPS satellites used for the positioning calculation changes. It is configured as follows.

【0038】また、測位演算部22及び速度ベクトル演
算部24は、受信回路部12にて同じタイミングで受信
された信号に基づき、測位位置データPg及び速度ベク
トルVgを、同じタイミングで繰り返し算出するように
構成されている。次に、測位演算部22では、少なくと
も3つ(好ましくは4つ以上)のGPS衛星との擬似距
離を同時に検出し、この擬似距離を公知の測位演算に適
用することにより、測位位置データPgを算出する。
The positioning calculation unit 22 and the speed vector calculation unit 24 repeatedly calculate the positioning position data Pg and the speed vector Vg at the same timing based on the signals received by the reception circuit unit 12 at the same timing. Is configured. Next, the positioning calculation unit 22 simultaneously detects at least three (preferably four or more) pseudoranges with the GPS satellites, and applies the pseudorange to a known positioning calculation to convert the positioning position data Pg. calculate.

【0039】また、速度ベクトル演算部24では、
(2)式を用いて求められる北方向速度Vn,及び東方
向速度Veに基づき、これらを合成することで速度ベク
トルVgを算出する。
In the speed vector calculation unit 24,
Based on the north direction velocity Vn and the east direction velocity Ve obtained by using the equation (2), the velocity vector Vg is calculated by synthesizing these.

【0040】[0040]

【数1】 (Equation 1)

【0041】但し、Vuは上方向速度、VLFは受信機の
基準発振器のオフセットの速度換算値、ni,ei,u
iはGPS衛星iに向かう単位ベクトルの北,東,上方
向余弦、Vdui は衛星iからの電波(搬送波)のドップ
ラーシフトの予測値(計算値)と受信回路部12にて検
出された搬送波の周波数変位の実測値との差を示してい
る。なお、この速度ベクトル演算部24における速度ベ
クトル演算は、測位演算部22における測位演算と同様
に公知のものである。
Where Vu is the upward speed, VLF is the speed-converted value of the offset of the reference oscillator of the receiver, ni, ei, u
i is the north, east, and upward cosine of the unit vector heading for the GPS satellite i, and Vdui is the predicted value (calculated value) of the Doppler shift of the radio wave (carrier wave) from the satellite i and the carrier wave detected by the receiving circuit unit 12. The difference from the measured value of the frequency displacement is shown. Note that the speed vector calculation in the speed vector calculation unit 24 is a known one, similarly to the positioning calculation in the positioning calculation unit 22.

【0042】ここで、速度誤差推定部20からの速度誤
差見積値Vth,測位演算部22からの測位位置データP
g,速度ベクトル演算部24からの速度ベクトルVgに
基づいて、演算部14が行うフィルタリング処理を、図
2に示すフローチャートに沿って説明する。
Here, the speed error estimated value Vth from the speed error estimating unit 20 and the positioning position data P from the positioning calculating unit 22
g, the filtering process performed by the calculation unit 14 based on the speed vector Vg from the speed vector calculation unit 24 will be described with reference to the flowchart shown in FIG.

【0043】本処理は車両の運転開始後、周期的(本実
施例では1[s]周期)に繰り返し起動される。本処理
が起動されると、まずS110では、速度誤差推定部2
0,測位演算部22,速度ベクトル演算部24から、速
度誤差見積値Vth,測位位置データPg,速度ベクトル
Vgを入力し、続くS120では、本処理の今回の起動
が、車両の運転開始後、最初の起動であるか否かを判断
する。
This process is started periodically (1 [s] cycle in the present embodiment) repeatedly after the start of the operation of the vehicle. When this processing is started, first, in S110, the speed error estimating unit 2
0, the estimated speed error value Vth, the measured position data Pg, and the speed vector Vg are input from the positioning calculation unit 22 and the speed vector calculation unit 24. In the subsequent S120, the current start of this process is performed after the vehicle starts running. It is determined whether or not this is the first activation.

【0044】そして、最初の起動であると判定された場
合には、S130に移行して、S110にて読み込んだ
測位位置データPgにより現在位置データP,及び推定
位置データPeを初期化し、続くS140では、演算タ
イミングを表すパラメータiを零に初期化(i←0)し
て、S170に移行する。
If it is determined that this is the first activation, the process proceeds to S130 to initialize the current position data P and the estimated position data Pe with the positioning position data Pg read in S110, and then to S140. Then, the parameter i representing the operation timing is initialized to zero (i ← 0), and the routine goes to S170.

【0045】以後、本処理の最初の起動を0回目とし
て、本処理のi回目の起動時(以下、i回目の演算タイ
ミングという)にS110にて読み込まれるデータをP
g(i),Vg(i)、及び同じi回目の演算タイミン
グにて算出される現在位置データをP(i)、更にP
(i)の算出に用いる推定位置データをPe(i)にて
表すものとする。
Thereafter, assuming that the first activation of this processing is the 0th activation, the data read in S110 at the i-th activation of this processing (hereinafter referred to as the i-th computation timing) is P
g (i), Vg (i), and the current position data calculated at the same i-th calculation timing are P (i) and P (i).
Assume that the estimated position data used for the calculation of (i) is represented by Pe (i).

【0046】一方、先のS120にて、本処理の今回の
起動が最初の起動ではないと判定された場合にはS15
0に移行し、ビル等による電波遮断によって捕捉中のG
PS衛星数が2以下となることにより、測位演算部22
及び速度ベクトル演算部24による演算処理が中断され
ているか否かを判断する。
On the other hand, if it is determined in S120 that the current activation of this process is not the first activation, the process proceeds to S15.
Go to 0, and G is being captured by blocking radio waves by buildings, etc.
When the number of PS satellites becomes 2 or less, the positioning operation unit 22
Then, it is determined whether or not the calculation processing by the speed vector calculation unit 24 is interrupted.

【0047】そして、測位演算が中断されることなく行
われていれば、S160に移行して、次の(3)式に従
って現在位置データP(i)の算出を行う。即ち、先の
S120にて読み込んだ測位位置データPg(i)と、
前回の演算タイミングで後述のS190,S240,S
250,S260のいずれかにて算出された推定位置デ
ータPe(i)との加重平均値を求めた後、S170に
移行する。
If the positioning calculation has been performed without interruption, the flow shifts to S160 to calculate the current position data P (i) according to the following equation (3). That is, the positioning position data Pg (i) read in S120, and
At the previous calculation timing, S190, S240, S
After calculating the weighted average value with the estimated position data Pe (i) calculated in any of 250 and S260, the process proceeds to S170.

【0048】[0048]

【数2】 (Equation 2)

【0049】但し、m>0,n>0であり、例えば、m
=9,n=1に設定する。なお、mの値をnより大きく
するほど、算出される現在位置データP(i)のばらつ
きは小さくなる。この加重平均値を求める処理は、測位
位置データPgに対する速度ベクトルVgによるフィル
タリングにて、現在位置データPの軌跡が速度ベクトル
Vgにて規定される所定の範囲に入るように現在位置デ
ータPを更新していくものであると考えることができ
る。
Here, m> 0 and n> 0. For example, m
= 9, n = 1. Note that, as the value of m is larger than n, the variation in the calculated current position data P (i) becomes smaller. In the process of obtaining the weighted average value, the current position data P is updated by filtering the positioning position data Pg with the speed vector Vg so that the locus of the current position data P falls within a predetermined range defined by the speed vector Vg. Can be thought of as something to do.

【0050】そして、S170では、S130又はS1
60にて設定,算出された現在位置データP(i)を、
ナビ演算部8に出力する。この時、ナビ演算部8では、
GPS受信機2からの現在位置データP(i)に基づ
き、車両の現在位置および現在位置付近の道路地図を表
示部6に表示させる。
Then, in S170, S130 or S1
The current position data P (i) set and calculated at 60 is
Output to the navigation operation unit 8. At this time, the navigation calculation unit 8
Based on the current position data P (i) from the GPS receiver 2, the display unit 6 displays the current position of the vehicle and a road map near the current position.

【0051】続くS180では、S110にて読み込ん
だ速度誤差見積値Vthを判定しきい値として、同じくS
110にて読み込んだ速度ベクトルの大きさ|Vg
(i)|、即ち車両の移動速度が、速度誤差見積値Vth
より小さいか否かを判断する。そして、速度ベクトルの
大きさ|Vg(i)|が速度誤差見積値Vth以上である
と判定された場合には、車両は移動中であるものとして
S190に移行し、移動中の車両の速度・方向は短時間
では急激に変化することはないことから、(4)式に示
すように、現在位置データP(i)と速度ベクトルVg
(i)とを加算することにより、次回の演算タイミング
での推定位置データPe(i+1)を求める。
At S180, the estimated speed error value Vth read at S110 is used as a judgment threshold value.
The magnitude of the velocity vector read at 110 | Vg
(I) |, that is, the moving speed of the vehicle is the speed error estimated value Vth
It is determined whether it is smaller than. If it is determined that the magnitude | Vg (i) | of the speed vector is equal to or greater than the estimated speed error value Vth, the process proceeds to S190 assuming that the vehicle is moving, and the speed / speed of the moving vehicle is determined. Since the direction does not change rapidly in a short time, the current position data P (i) and the speed vector Vg are calculated as shown in the equation (4).
(I) is added to obtain estimated position data Pe (i + 1) at the next calculation timing.

【0052】 Pe(i+1)=P(i)+Vg(i) (4) 一方、先のS180にて、速度ベクトルの大きさ|Vg
(i)|が速度誤差見積値Vthより小さいと判定された
場合には、車両は停止中であるものとしてS200に移
行し、今度は、測位位置と現在位置との距離|Pg
(i)−P(i)|が、予め設定された上限距離Pthよ
り大きいか否かを判断し、肯定判定された場合、即ち、
測位位置と現在位置との距離が許容範囲(上限距離Pt
h)を越えている場合にはS210に移行する。
Pe (i + 1) = P (i) + Vg (i) (4) On the other hand, at S180, the magnitude of the velocity vector | Vg
(I) If it is determined that | is smaller than the estimated speed error value Vth, the process proceeds to S200 assuming that the vehicle is stopped, and then the distance | Pg between the positioning position and the current position is determined.
(I) -P (i) | is determined to be greater than or equal to a preset upper limit distance Pth, and if a positive determination is made, that is,
The distance between the positioning position and the current position is within the allowable range (upper limit distance Pt).
If h) is exceeded, the process moves to S210.

【0053】S210では、予め設定された上限時間T
thを計時するタイマが起動済みであるか否かを判断し、
未だ起動していなければS220に移行して、タイマを
起動後、S250に移行する。但し、タイマは、S15
0にて肯定判定された場合、及びS180にて否定判定
された場合には、動作が停止され計時値もリセットされ
る。
At S210, a preset upper limit time T
judge whether the timer for measuring th has been started,
If the timer has not been started, the process proceeds to S220, and after the timer is started, the process proceeds to S250. However, the timer is S15
If a positive determination is made at 0 and if a negative determination is made at S180, the operation is stopped and the clock value is reset.

【0054】一方、S210にてタイマが起動済みであ
ると判定された場合には、S230に移行して、今度は
タイマがタイムアウトしているか否かを判断し、タイム
アウトしていればS240に移行して、(5)式に示す
ように、S110にて読み込んだ測位位置データPg
(i)を、次回の推定位置データPe(i+1)として
設定する。
On the other hand, if it is determined in S210 that the timer has been started, the process proceeds to S230, and it is determined whether the timer has timed out. If so, the process proceeds to S240. Then, as shown in equation (5), the positioning position data Pg read in S110
(I) is set as the next estimated position data Pe (i + 1).

【0055】 Pe(i+1)=Pg(i) (5) また、先のS200にて、測位位置と現在位置との距離
|Pg(i)−P(i)|が上限距離Pth以下(即ち許
容範囲内)である場合、又は、上限距離Pthより大きく
てもタイムアウトしていない場合には、S250に移行
して、(6)式に示すように、今回の演算タイミングで
の推定位置データPe(i)を、次回の演算タイミング
での推定位置データPe(i+1)として設定する。
Pe (i + 1) = Pg (i) (5) In S200, the distance | Pg (i) −P (i) | between the measured position and the current position is equal to or smaller than the upper limit distance Pth (that is, allowable). If it is within the range), or if the timeout has not occurred even if it is greater than the upper limit distance Pth, the flow shifts to S250, and as shown in Expression (6), the estimated position data Pe ( i) is set as the estimated position data Pe (i + 1) at the next calculation timing.

【0056】 Pe(i+1)=Pe(i) (6) また、先のS150にて、測位中断中であると判定され
た場合には、S260に移行して、(7)式に示すよう
に、測位が中断される直前に最後に読み込まれた速度ベ
クトルVg(i0)を、現在位置データP(i)に加算
することにより、推定位置データPe(i+1)を算出
する。
Pe (i + 1) = Pe (i) (6) If it is determined in S150 that the positioning is suspended, the process proceeds to S260, and as shown in Expression (7). Then, the estimated position data Pe (i + 1) is calculated by adding the velocity vector Vg (i0) read last immediately before the positioning is interrupted to the current position data P (i).

【0057】 Pe(i+1)=P(i)+Vg(iO) (7) ここで、i0 は、最後に測位した時刻を示すパラメータ
を示す。つまり、測位中断中に車両の速度・方位が急激
に変化しなければ、このような速度ベクトルV(i0)に
基づいて、次回以降の演算タイミングでの推定位置デー
タPe(i+1)を、ある程度正確に求めることがで
き、測位が再開された時には、その推定位置データPe
によって現在位置をある程度正確に補正することができ
るのである。
Pe (i + 1) = P (i) + Vg (iO) (7) Here, i0 indicates a parameter indicating the time of last positioning. In other words, if the speed / azimuth of the vehicle does not change abruptly during the suspension of the positioning, the estimated position data Pe (i + 1) at the next and later calculation timings can be corrected to some extent based on the speed vector V (i0). And when the positioning is resumed, the estimated position data Pe
Thus, the current position can be corrected to some extent accurately.

【0058】そして、S190,S240,S250,
S260のいずれかにて次回の推定位置データPe(i
+1)が算出されると、S270に移行して、パラメー
タiをインクリメント(i←i+1)して本処理を終了
する。即ち、本処理においては、各演算タイミング毎
に、読み込んだ(S110)測位位置データPg(i)
と、前回の演算タイミングにて求めた推定位置データP
e(i)とにより、現在位置データP(i)を求めてい
る(S160)。
Then, S190, S240, S250,
In any of S260, the next estimated position data Pe (i
When +1) is calculated, the flow shifts to S270, where the parameter i is incremented (i ← i + 1), and this processing ends. That is, in this processing, the read (S110) positioning position data Pg (i) at each calculation timing.
And the estimated position data P obtained at the previous calculation timing
The current position data P (i) is obtained from e (i) (S160).

【0059】また、本処理においては、次回の演算タイ
ミングで用いる推定位置データPe(i+1)の算出
を、場合に応じて次のように行っている。即ち、車両が
移動中(S180−NO)の場合は、車両の速度・方向
が急激に変化することはないものとして、現在位置デー
タP(i)と速度ベクトルVg(i)とを加算すること
により推定位置データPe(i+1)を求め(S19
0)、一方、車両が停止中(S180−YES)の場合
は、基本的には、推定位置データの更新を行わず、今回
の推定位置データPe(i)をそのまま次回の推定位置
データPe(i+1)として設定(S250)してい
る。
In this processing, the calculation of the estimated position data Pe (i + 1) used at the next calculation timing is performed as follows, as the case may be. That is, when the vehicle is moving (S180-NO), the current position data P (i) and the speed vector Vg (i) are added assuming that the speed and direction of the vehicle do not change rapidly. To obtain the estimated position data Pe (i + 1) (S19).
0) On the other hand, when the vehicle is stopped (S180-YES), basically, the estimated position data Pe (i) is not updated and the next estimated position data Pe (i) is not updated. (i + 1)) (S250).

【0060】但し、本処理では、車両が停止しているか
否かの判断を、速度ベクトルの大きさ|Vg(i)|が
速度誤差見積値Vthより小さいか否かにより行っている
ため、車両が停止中であると判断しても、実際には、速
度誤差見積値Vthより小さい速度で低速移動している場
合が考えられる。
However, in this processing, whether the vehicle is stopped or not is determined based on whether the magnitude | Vg (i) | of the speed vector is smaller than the estimated speed error value Vth. Is determined to be stopped, the vehicle may actually be moving at a lower speed at a speed smaller than the estimated speed error value Vth.

【0061】この場合、推定位置データPeの更新が行
われないことから、推定位置データPe(i)に基づい
て算出される現在位置データP(i)と、測位演算部2
2にて算出される測位位置データPg(i)との間の誤
差は、時間の経過に従って増大する。
In this case, since the estimated position data Pe is not updated, the current position data P (i) calculated based on the estimated position data Pe (i) and the positioning operation unit 2
The error from the positioning position data Pg (i) calculated in 2 increases with the passage of time.

【0062】そのため、この誤差、即ち、現在位置と測
位位置との距離|Pg(i)−P(i)|が、許容範囲
(上限距離Pth)を越えた場合(S200−YES)に
は、推定位置データPe(i+1)を測位位置データP
g(i)により再初期化(S240)して、両データ間
の誤差が小さくなるようにしている。
Therefore, when this error, that is, the distance | Pg (i) -P (i) | between the current position and the measured position exceeds the allowable range (upper limit distance Pth) (S200-YES), The estimated position data Pe (i + 1) is converted to the positioning position data P
Reinitialization is performed using g (i) (S240) so that the error between the two data is reduced.

【0063】なお、現在位置と測位位置との距離|Pg
(i)−P(i)|が許容範囲を越えるような状況は、
車両の低速移動時に限らず、マルチパスが生じた場合に
も発生する。但し、マルチパスは、GPS衛星の移動に
より、一定時間が経過すると解消される可能性が高い。
このため、上記状況が上限時間Tth以上継続した場合
(S210−YES,S230−YES)にのみ、マル
チパスの影響ではないものとして、上記再初期化(S2
40)を行うことにより、推定位置データPe,ひいて
は現在位置データPから、マルチパスの影響が排除され
るようにしている。
The distance | Pg between the current position and the measured position
A situation where (i) -P (i) |
It occurs not only when the vehicle moves at low speed but also when a multipath occurs. However, there is a high possibility that the multipath will be canceled after a certain period of time due to the movement of the GPS satellite.
Therefore, only when the above situation continues for the upper limit time Tth or more (S210-YES, S230-YES), the re-initialization (S2
By performing step 40), the influence of multipath is eliminated from the estimated position data Pe and, consequently, the current position data P.

【0064】ここで、図3は、フィルタリング処理によ
り求められる現在位置データP(i)および推定位置デ
ータPe(i)が変化する様子を表す説明図である。な
お、図3(a)〜(c)のいずれも、演算タイミングi
=2までは、車両が速度誤差見積値Vth以上の速度で走
行している場合を示している。
FIG. 3 is an explanatory diagram showing how the current position data P (i) and the estimated position data Pe (i) obtained by the filtering process change. 3A to 3C, the calculation timing i
= 2 indicates a case where the vehicle is running at a speed equal to or higher than the estimated speed error value Vth.

【0065】この時、測位位置データPg(i)と、よ
りバラツキの小さい速度ベクトルVg(i−1)に基づ
く推定位置データPe(i)との加重平均により算出さ
れる現在位置データP(i)は、測位位置データPg
(i)そのものよりもばらつきが小さくなる。
At this time, the current position data P (i) calculated by the weighted average of the positioning position data Pg (i) and the estimated position data Pe (i) based on the speed vector Vg (i-1) having less variation. ) Indicates positioning position data Pg
(I) Variation is smaller than that of itself.

【0066】そして、図3(a)において、演算タイミ
ングi=3以降は、車両が実際に停止している場合を示
しており、この場合、推定位置データPe(i)が固定
され、速度ベクトルVg(i)に含まれる誤差の影響を
全く受けないため、現在位置データP(i)のばらつき
が、より一層低減されることになる。
FIG. 3A shows a case where the vehicle is actually stopped after the calculation timing i = 3. In this case, the estimated position data Pe (i) is fixed and the speed vector Since there is no influence of the error included in Vg (i), the variation in the current position data P (i) is further reduced.

【0067】また、図3(b)において、演算タイミン
グi=3以降は、車両が速度誤差見積値Vthより小さい
速度で低速走行している場合を示している。ここでは、
演算タイミングi=4の時に、現在位置と測位位置との
距離|Pg(4)−P(4)|が上限距離Pthより大き
くなり、その後、上限時間Tthを経過した演算タイミン
グi=X−1の時に、推定位置データPeの再初期化が
行われている。
FIG. 3B shows a case in which the vehicle is traveling at a low speed at a speed smaller than the estimated speed error value Vth after the calculation timing i = 3. here,
When the calculation timing i = 4, the distance | Pg (4) −P (4) | between the current position and the measured position becomes larger than the upper limit distance Pth, and thereafter, the calculation timing i = X−1 after the upper limit time Tth has elapsed. At the time of, the estimated position data Pe is re-initialized.

【0068】このように、低速走行を停止状態と判定し
てしまうと、測位位置データPg(i)と、推定位置デ
ータPe(i)との誤差が大きくなるが、再初期化を行
うことで、現在位置に蓄積されるこの誤差が解消される
ことがわかる。更に、図3(c)において、演算タイミ
ングi=3以降は、車両は実際に停止しているがマルチ
パスが発生し、その後、演算タイミングi=6にて、マ
ルチパスが解消した場合を示している。つまり、速度ベ
クトルVgが速度誤差見積値Vthより小さければ、推定
位置データPeの更新が行われないため、マルチパスが
発生して測位位置データPgが大きく変化したとして
も、これに追従してしまうことがなく、マルチパスが解
消された時に算出される現在位置データP(6)は、マ
ルチパスが発生する以前に求められた現在位置データP
(2)の近傍に復帰する。つまり、その後、車両が走行
を再開した時に、マルチパスの影響が蓄積されていない
現在位置データP(i)が算出される。
As described above, if the low-speed running is determined to be in the stop state, the error between the positioning position data Pg (i) and the estimated position data Pe (i) increases. It can be seen that this error accumulated at the current position is eliminated. Further, FIG. 3C shows a case where the vehicle is actually stopped after the calculation timing i = 3, but a multipath occurs, and thereafter, the multipath is resolved at the calculation timing i = 6. ing. In other words, if the speed vector Vg is smaller than the estimated speed error value Vth, the estimated position data Pe is not updated. Therefore, even if a multipath occurs and the positioning position data Pg greatly changes, the positioning position data Pg follows the change. And the current position data P (6) calculated when the multipath is canceled is the current position data P (6) obtained before the occurrence of the multipath.
It returns to the vicinity of (2). That is, when the vehicle resumes running, the current position data P (i) in which the influence of the multipath is not accumulated is calculated.

【0069】そして、このようにして算出された現在位
置データPは、図2のフローチャートで説明したよう
に、S170でナビ演算部8に出力される。そして、車
両の現在位置が、その付近の道路地図と共に表示部6に
表示される。図4に、実際に走行した場合の現在位置の
走行軌跡を示す。図5に示した従来装置、即ち速度誤差
見積値Vthによる停止判定を行わない場合は、車両の停
止時に算出される現在位置データPのばらつきが大きい
が、本実施例のGPS受信機2、即ち停止判定を行い、
停止状態にあると判定すると推定位置データPeの更新
を禁止する場合は、停止時における現在位置データPの
ばらつきが大幅に減少し、安定した軌跡となっているこ
とがわかる。
The current position data P calculated in this manner is output to the navigation operation section 8 in S170 as described in the flowchart of FIG. Then, the current position of the vehicle is displayed on the display unit 6 together with a nearby road map. FIG. 4 shows a traveling locus of the current position when the vehicle actually travels. In the case of not performing the stop determination based on the conventional device shown in FIG. 5, that is, when the stop determination based on the estimated speed error value Vth is performed, the dispersion of the current position data P calculated when the vehicle stops is large. Make a stop decision,
When the update of the estimated position data Pe is prohibited when it is determined that the vehicle is in the stop state, it is understood that the variation in the current position data P at the time of the stop is greatly reduced, and the trajectory is stable.

【0070】以上説明したように、本実施例のGPS受
信機2によれば、現在位置データP(i)と速度ベクト
ルVg(i)とに基づいて次の演算タイミングでの推定
位置データPe(i+1)を算出し、次の演算タイミン
グでは、この推定位置データPe(i+1)と測位位置
データPg(i+1)とに基づいて、過去の測位位置デ
ータPg(i)をそのまま用いることなく、現在位置デ
ータP(i+1)を算出しているので、現在位置データ
Pを、時間遅れを生じさせることなく精度よく算出でき
る。
As described above, according to the GPS receiver 2 of this embodiment, the estimated position data Pe () at the next calculation timing is calculated based on the current position data P (i) and the velocity vector Vg (i). i + 1) is calculated, and at the next calculation timing, based on the estimated position data Pe (i + 1) and the positioning position data Pg (i + 1), the current positioning position data Pg (i) is not used as it is, and the current position is determined. Since the data P (i + 1) is calculated, the current position data P can be calculated accurately without causing a time delay.

【0071】また、本実施例のGPS受信機2によれ
ば、速度ベクトルの大きさ|Vg(i)|が、測位に用
いるGPS衛星の状態によって決まる速度誤差見積値V
thより小さく、速度ベクトルVg中の誤差成分による影
響が大きい場合には、車両が停止中であると判定して、
推定位置データPe(i)を固定し、速度ベクトルVg
(i)による値の更新を行わないようにしているので、
速度ベクトルVg(i)の誤差成分によって、停止中の
現在位置データP(i)がばらついてしまうことを確実
に抑えることができる。
According to the GPS receiver 2 of the present embodiment, the magnitude | Vg (i) | of the velocity vector is equal to the estimated velocity error value V determined by the state of the GPS satellite used for positioning.
If it is smaller than th and the influence of the error component in the speed vector Vg is large, it is determined that the vehicle is stopped,
With the estimated position data Pe (i) fixed, the velocity vector Vg
Since the value is not updated according to (i),
Variations in the stopped current position data P (i) due to the error component of the speed vector Vg (i) can be reliably suppressed.

【0072】更に、本実施例のGPS受信機2では、車
両が停止中である間に、測位位置と現在位置との距離|
Pg(i)−P(i)|が、上限距離Pthより離れた状
態が、上限時間Tth以上継続した場合に、次回の演算タ
イミングでの推定位置データP(i+1)を、測位位置
データPg(i)で再初期化するようにされている。
Further, in the GPS receiver 2 of the present embodiment, while the vehicle is stopped, the distance between the positioning position and the current position |
If Pg (i) -P (i) | is longer than the upper limit distance Pth for more than the upper limit time Tth, the estimated position data P (i + 1) at the next calculation timing is converted to the positioning position data Pg ( Re-initialization is performed in i).

【0073】このため、停止中であると判定されている
にも関わらず、実際には車両が速度誤差見積値Vthより
小さな速度で低速移動している場合に、これを検出し
て、低速移動により生じた誤差を解消できるだけでな
く、マルチパスの影響により、測位位置と現在位置との
距離が上限距離以上離れた状態が発生したとしても、上
限時間Tth以上継続しなければ、測位位置データPg
(i)による現在位置データP(i)の再初期化を行わ
ないため、マルチパスによる測位位置データPg(i)
の異常が現在位置データP(i)に蓄積されてしまうこ
とを防止できる。
For this reason, when the vehicle is actually moving at a low speed at a speed smaller than the estimated speed error value Vth, it is detected that the vehicle is moving at a low speed, although it is determined that the vehicle is stopped. In addition to eliminating the error caused by the above, even if a state in which the distance between the positioning position and the current position is longer than the upper limit distance due to the influence of multipath occurs, if the distance does not continue for the upper limit time Tth, the positioning position data Pg
Since the re-initialization of the current position data P (i) by (i) is not performed, the positioning position data Pg (i) by multipath is not used.
Can be prevented from being accumulated in the current position data P (i).

【0074】また、本実施例のGPS受信機2では、停
止中であるか否かを判定するための速度誤差見積値Vth
を可変設定し、しかも、その上限値を制限するようにさ
れている。従って、測位に使用する人工衛星の組合せや
配置により、刻々と変化する誤差に応じて、常に最適な
速度誤差見積値Vthを設定できるため、システムの性能
を最大限に引き出すことができると共に、速度誤差見積
値Vthが必要以上に大きな値に設定されることにより、
地図上に表示される現在位置の表示が、車両が動き始め
ているにも関わらずなかなか動き出さないといったよう
な、不自然な動きとなることを確実に防止できる。
In the GPS receiver 2 of the present embodiment, the estimated speed error value Vth for determining whether or not the vehicle is stopped.
Is variably set, and the upper limit is limited. Therefore, the optimum speed error estimated value Vth can always be set according to the error that changes every moment depending on the combination and arrangement of the satellites used for positioning, so that the system performance can be maximized and the speed can be maximized. By setting the error estimation value Vth to a value larger than necessary,
It is possible to reliably prevent the display of the current position displayed on the map from becoming an unnatural movement, such as a situation in which the vehicle does not readily start moving even though it is starting to move.

【0075】以上、本発明の一実施例について説明した
が、本発明は上記実施例に限定されるものではなく、本
発明の要旨を逸脱しない範囲において、様々な態様にて
実施することができる。例えば、上記実施例では、測位
位置データPg(i)と推定位置データPe(i)とか
ら現在位置データP(i)を算出する際に、加重平均の
重み付け係数m,nを固定して用いているが、例えば、
測位に使用するGPS衛星の組合せが変化する等して、
測位精度が変化する毎に、この測位精度に応じて重み付
け係数m,nを変化させるようにしてもよい。具体的に
は、測位精度が高い時には、測位位置データPg(i)
の重み付け係数nを大きくし、逆に測位精度が低い時に
は、推定位置データPe(i)の重み付け係数mを大き
くすればよい。この場合、現在位置データP(i)の精
度及び安定性を向上させることができる。
As described above, one embodiment of the present invention has been described. However, the present invention is not limited to the above-described embodiment, and can be implemented in various modes without departing from the gist of the present invention. . For example, in the above embodiment, when calculating the current position data P (i) from the positioning position data Pg (i) and the estimated position data Pe (i), the weighting coefficients m and n of the weighted average are fixed and used. But, for example,
As the combination of GPS satellites used for positioning changes,
Each time the positioning accuracy changes, the weighting coefficients m and n may be changed according to the positioning accuracy. Specifically, when the positioning accuracy is high, the positioning position data Pg (i)
When the positioning accuracy is low, the weighting coefficient m of the estimated position data Pe (i) may be increased. In this case, the accuracy and stability of the current position data P (i) can be improved.

【0076】また、実際の車両では加速度や回転が生じ
ているため、車両の運動を一定速度と仮定して推定位置
データPe(i)を算出している上記実施例では、推定
位置データPe(i)に誤差が生じる。従って、一定速
度から外れる車両の運動分を推定位置データPe(i)
の誤差として考慮し、現在位置データP(i)を算出す
る際の重み付け係数m,nを決定するようにすれば、よ
り精度良く現在位置データP(i)を算出することがで
きる。
Further, since acceleration and rotation occur in the actual vehicle, the estimated position data Pe (i) is calculated in the above embodiment in which the motion of the vehicle is assumed to be constant and the estimated position data Pe (i) is calculated. An error occurs in i). Therefore, the movement amount of the vehicle deviating from the constant speed is estimated position data Pe (i).
If the weighting factors m and n for calculating the current position data P (i) are determined in consideration of the error of the current position data P (i), the current position data P (i) can be calculated more accurately.

【0077】また、上記実施例では、S190におい
て、車両の運動を一定速度運動としてモデル化し、
(4)式を用いて推定位置データPe(i+1)を算出
しているが加速度運動まで考慮し、(8)式を用いて推
定位置データPe(i+1)を算出するようにしてもよ
い。
In the above embodiment, in S190, the motion of the vehicle is modeled as a constant speed motion.
Although the estimated position data Pe (i + 1) is calculated using the expression (4), the estimated position data Pe (i + 1) may be calculated using the expression (8) in consideration of the acceleration motion.

【0078】 Pe(i+1)=P(i)+{V(i−1)+V(i)}/2 (8) 更に、上記実施例では、S240において、推定位置デ
ータPe(i+1)の再初期化を行う際に、測位位置デ
ータPg(i)をそのまま用いているが、タイマが起動
してからタイムアウトするまでの間に、演算タイミング
毎にS110にて読み込まれた複数の測位位置データP
gに基づいて得られる平均化された測位位置データを用
いるようにしてもよい。この場合の平均化処理として
は、単純平均してもよいし、時系列に従って重み付け平
均してもよいし、その時々の測位精度に応じて重み付け
平均してもよい。
Pe (i + 1) = P (i) + {V (i−1) + V (i)} / 2 (8) Further, in the above embodiment, in S240, the re-initialization of the estimated position data Pe (i + 1) is performed. When performing the conversion, the positioning position data Pg (i) is used as it is, but a plurality of positioning position data Pg read in S110 at each calculation timing from the start of the timer to the time-out.
The averaged positioning position data obtained based on g may be used. In this case, the averaging process may be simple averaging, weighted averaging according to a time series, or weighted averaging according to the positioning accuracy at each time.

【0079】また、上記実施例では、停止状態を、速度
ベクトルの大きさが予め設定された値より小さいか否か
により判定しているが、車速パルス等、車速を測定する
ために車載されている他の装置からの測定信号等に基づ
いて、直接停止状態を判定してもよい。
In the above-described embodiment, the stop state is determined based on whether or not the magnitude of the speed vector is smaller than a predetermined value. Alternatively, the stop state may be directly determined based on a measurement signal or the like from another device.

【0080】上記実施例において、GPS受信機2が位
置検出装置、アンテナ10及び受信回路部12が受信手
段、速度誤差推定部20が判定しきい値設定手段、測位
演算部22が疑似距離検出手段及び測位演算手段、速度
ベクトル演算部24が相対速度検出手段及び速度ベクト
ル演算手段、S190が位置推定手段、S160が現在
位置更新手段、S180が停止判定手段、S250が禁
止手段、S210〜S240が推定位置初期化手段に相
当する。
In the above embodiment, the GPS receiver 2 is a position detecting device, the antenna 10 and the receiving circuit unit 12 are a receiving unit, the speed error estimating unit 20 is a determination threshold setting unit, and the positioning calculating unit 22 is a pseudo distance detecting unit. And speed calculation means 24, relative speed detection means and speed vector calculation means, S190 is position estimation means, S160 is current position update means, S180 is stop determination means, S250 is inhibition means, and S210 to S240 are estimation. It corresponds to position initialization means.

【図面の簡単な説明】[Brief description of the drawings]

【図1】 実施例のGPS受信機を用いて構成されたナ
ビゲーション装置の全体構成を表すブロック図である。
FIG. 1 is a block diagram illustrating an overall configuration of a navigation device configured using a GPS receiver according to an embodiment.

【図2】 フィルタリング部にて実行されるフィルタリ
ング処理の内容を表すフローチャートである。
FIG. 2 is a flowchart illustrating the content of a filtering process performed by a filtering unit.

【図3】 フィルタリング処理により算出される現在位
置データ及び推定位置データが変化する様子を表す説明
図である。
FIG. 3 is an explanatory diagram showing how the current position data and the estimated position data calculated by the filtering process change.

【図4】 地図上に表示された現在位置データの軌跡を
示す説明図である。
FIG. 4 is an explanatory diagram showing a locus of current position data displayed on a map.

【図5】 従来装置による現在位置データの軌跡を示す
説明図である。
FIG. 5 is an explanatory diagram showing a locus of current position data by a conventional device.

【符号の説明】[Explanation of symbols]

2…GPS受信機 4…地図データ記憶媒体
6…表示部 8…ナビ演算部 10…アンテナ 1
2…受信回路部 14…演算部 20…速度誤差推定部
22…測位演算部 24…速度ベクトル演算部 26…フィルタリング
2 ... GPS receiver 4 ... Map data storage medium
6 display unit 8 navigation operation unit 10 antenna 1
2 ... Receiving circuit unit 14 ... Calculation unit 20 ... Speed error estimation unit
22 positioning measurement unit 24 speed vector calculation unit 26 filtering unit

─────────────────────────────────────────────────────
────────────────────────────────────────────────── ───

【手続補正書】[Procedure amendment]

【提出日】平成11年4月28日(1999.4.2
8)
[Submission date] April 28, 1999 (1999.4.2
8)

【手続補正1】[Procedure amendment 1]

【補正対象書類名】明細書[Document name to be amended] Statement

【補正対象項目名】0048[Correction target item name] 0048

【補正方法】変更[Correction method] Change

【補正内容】[Correction contents]

【0048】[0048]

【数2】 (Equation 2)

───────────────────────────────────────────────────── フロントページの続き (72)発明者 尾崎 義隆 愛知県刈谷市昭和町1丁目1番地 株式会 社デンソー内 Fターム(参考) 2C032 HB22 HC08 HD03 2F029 AA02 AB07 AC02 AC06 AC12 AC14 AD03 5H180 AA01 BB15 FF05 FF22 FF27 FF32 9A001 GG01 JJ78  ────────────────────────────────────────────────── ─── Continuing on the front page (72) Inventor Yoshitaka Ozaki 1-1-1, Showa-cho, Kariya-shi, Aichi F-term in DENSO Corporation (reference) 2C032 HB22 HC08 HD03 2F029 AA02 AB07 AC02 AC06 AC12 AC14 AD03 5H180 AA01 BB15 FF05 FF22 FF27 FF32 9A001 GG01 JJ78

Claims (6)

【特許請求の範囲】[Claims] 【請求項1】 複数の人工衛星から電波を受信する受信
手段と、 該受信手段が受信する電波に重畳されたタイミング情報
に基づき、前記人工衛星のそれぞれとの擬似距離を検出
する擬似距離検出手段と、 該擬似距離検出手段にて検出された擬似距離に基づき、
予め設定された演算タイミング毎に測位位置を繰り返し
求める測位演算手段と、 前記受信手段が受信する電波の周波数変位に基づき、当
該装置と前記人工衛星のそれぞれとの相対速度を検出す
る相対速度検出手段と、 該相対速度検出手段にて検出された相対速度に基づき、
前記演算タイミング毎に当該装置の移動速度及び移動方
向を示す速度ベクトルを繰り返し求める速度ベクトル演
算手段と、 該速度ベクトル演算手段にて求められた速度ベクトル、
及び予め設定された当該装置の現在位置に基づいて、次
回の演算タイミングでの当該装置の推定位置を求める位
置推定手段と、 前記測位演算手段にて求められた測位位置と、前回の演
算タイミング時に前記位置推定手段にて求められた推定
位置との加重平均値により、前記現在位置を更新する現
在位置更新手段と、 を備えた位置検出装置において、 当該位置検出装置の停止状態を判定する停止判定手段
と、 前記停止判定手段により停止状態と判定された場合に、
前記位置推定手段による前記推定位置の更新を禁止する
禁止手段と、 を備えることを特徴とする位置検出装置。
1. A receiving means for receiving radio waves from a plurality of artificial satellites, and a pseudo distance detecting means for detecting a pseudo distance to each of the artificial satellites based on timing information superimposed on the radio waves received by the receiving means. And, based on the pseudo distance detected by the pseudo distance detecting means,
Positioning calculation means for repeatedly obtaining a positioning position for each predetermined calculation timing; and relative speed detection means for detecting a relative speed between the device and the artificial satellite based on a frequency displacement of a radio wave received by the reception means. Based on the relative speed detected by the relative speed detection means,
Speed vector calculating means for repeatedly obtaining a speed vector indicating a moving speed and a moving direction of the device at each of the calculation timings; a speed vector obtained by the speed vector calculating means;
And a position estimating means for obtaining an estimated position of the device at the next calculation timing based on the preset current position of the device, and a positioning position obtained by the positioning calculation means, A current position updating means for updating the current position based on a weighted average value of the estimated position obtained by the position estimating means; and a stop determination for judging a stop state of the position detecting apparatus. Means, when it is determined that the stop state by the stop determination means,
Prohibiting means for prohibiting updating of the estimated position by the position estimating means.
【請求項2】 請求項1記載の位置検出装置において、 前記停止判定手段は、前記速度ベクトル演算手段にて求
められる速度ベクトルの大きさが予め設定された判定し
きい値より小さい場合に、停止状態と判定することを特
徴とする位置検出装置。
2. The position detecting device according to claim 1, wherein said stop judging means stops when a magnitude of a speed vector obtained by said speed vector calculating means is smaller than a predetermined judgment threshold value. A position detecting device for determining a state.
【請求項3】 請求項1又は請求項2記載の位置検出装
置において、 前記停止状態中に、前記測位位置と前記現在位置とが予
め設定された上限距離以上離れている状態が、予め設定
された上限時間以上継続した場合、前記測位位置に基づ
いて前記推定位置を初期化する推定位置初期化手段を設
けたことを特徴とする位置検出装置。
3. The position detection device according to claim 1, wherein a state in which the positioning position and the current position are separated by a predetermined upper limit distance or more during the stop state is set in advance. A position detecting device provided with an estimated position initializing means for initializing the estimated position based on the measured position, when the estimated time has continued for a predetermined time or more.
【請求項4】 請求項3記載の位置検出装置において、 前記推定位置初期化手段は、前記停止状態中に求められ
た複数の測位位置を平均化して得られる位置により、前
記推定位置を初期化することを特徴とする位置検出装
置。
4. The position detecting device according to claim 3, wherein the estimated position initializing means initializes the estimated position by a position obtained by averaging a plurality of positioning positions obtained during the stop state. A position detecting device.
【請求項5】 請求項1ないし請求項4いずれか記載の
位置検出装置において、 前記受信手段が受信する電波に重畳された精度情報に基
づいて、前記停止判定手段が用いる判定しきい値を可変
設定する判定しきい値設定手段を設けたことを特徴とす
る位置検出装置。
5. The position detecting device according to claim 1, wherein a judgment threshold value used by said stop judging means is variable based on accuracy information superimposed on a radio wave received by said receiving means. A position detecting device provided with a judgment threshold value setting means for setting.
【請求項6】 請求項5記載の位置検出装置において、 前記判定しきい値設定手段は、前記制度情報を所定の演
算式に適用して判定しきい値を算出し、該算出値が予め
設定された上限値より大きい場合には、前記判定しきい
値を前記上限値に制限することを特徴とする位置検出装
置。
6. The position detecting device according to claim 5, wherein the determination threshold value setting means calculates a determination threshold value by applying the accuracy information to a predetermined arithmetic expression, and the calculated value is set in advance. The position detecting device, wherein, when the value is larger than the upper limit, the determination threshold is limited to the upper limit.
JP11986199A 1999-04-27 1999-04-27 Position detection device Expired - Lifetime JP3601354B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP11986199A JP3601354B2 (en) 1999-04-27 1999-04-27 Position detection device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP11986199A JP3601354B2 (en) 1999-04-27 1999-04-27 Position detection device

Publications (2)

Publication Number Publication Date
JP2000310538A true JP2000310538A (en) 2000-11-07
JP3601354B2 JP3601354B2 (en) 2004-12-15

Family

ID=14772095

Family Applications (1)

Application Number Title Priority Date Filing Date
JP11986199A Expired - Lifetime JP3601354B2 (en) 1999-04-27 1999-04-27 Position detection device

Country Status (1)

Country Link
JP (1) JP3601354B2 (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005274506A (en) * 2004-03-26 2005-10-06 Hitachi Ltd Navigation device and present position estimation method
JP2006125876A (en) * 2004-10-26 2006-05-18 Alpine Electronics Inc Position calculation system and calculation method for moving body
JP2008134086A (en) * 2006-11-27 2008-06-12 Denso Corp On-vehicle navigation device
JP2009156734A (en) * 2007-12-27 2009-07-16 Seiko Epson Corp Positioning method, program, positioning device, and electronic device
JP2009300361A (en) * 2008-06-17 2009-12-24 Seiko Epson Corp Positioning technique, program, and positioning apparatus
US7688259B2 (en) 2006-05-29 2010-03-30 Seiko Epson Corporation Positioning device, method of controlling positioning device, and recording medium
JP2011122903A (en) * 2009-12-10 2011-06-23 Kyosan Electric Mfg Co Ltd Train control system
JP2014035218A (en) * 2012-08-07 2014-02-24 Seiko Epson Corp Stop continuation determination method and stop continuation determination device
JP2015225615A (en) * 2014-05-29 2015-12-14 トヨタ自動車株式会社 Drive assist system
US9799836B2 (en) 2014-02-28 2017-10-24 Seiko Epson Corporation Light emitting element, light emitting device, authentication device, and electronic device
JP7438718B2 (en) 2019-11-12 2024-02-27 ヤンマーパワーテクノロジー株式会社 Automated driving system for work vehicles

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005274506A (en) * 2004-03-26 2005-10-06 Hitachi Ltd Navigation device and present position estimation method
JP2006125876A (en) * 2004-10-26 2006-05-18 Alpine Electronics Inc Position calculation system and calculation method for moving body
JP4498095B2 (en) * 2004-10-26 2010-07-07 アルパイン株式会社 Moving body position calculating apparatus and calculating method
US8028014B2 (en) 2006-05-29 2011-09-27 Seiko Epson Corporation Positioning device, method of controlling positioning device, and recording medium
US7688259B2 (en) 2006-05-29 2010-03-30 Seiko Epson Corporation Positioning device, method of controlling positioning device, and recording medium
JP2008134086A (en) * 2006-11-27 2008-06-12 Denso Corp On-vehicle navigation device
JP2009156734A (en) * 2007-12-27 2009-07-16 Seiko Epson Corp Positioning method, program, positioning device, and electronic device
JP2009300361A (en) * 2008-06-17 2009-12-24 Seiko Epson Corp Positioning technique, program, and positioning apparatus
JP2011122903A (en) * 2009-12-10 2011-06-23 Kyosan Electric Mfg Co Ltd Train control system
JP2014035218A (en) * 2012-08-07 2014-02-24 Seiko Epson Corp Stop continuation determination method and stop continuation determination device
US9799836B2 (en) 2014-02-28 2017-10-24 Seiko Epson Corporation Light emitting element, light emitting device, authentication device, and electronic device
JP2015225615A (en) * 2014-05-29 2015-12-14 トヨタ自動車株式会社 Drive assist system
US10106154B2 (en) 2014-05-29 2018-10-23 Toyota Jidosha Kabushiki Kaisha Driving support apparatus
US10696297B2 (en) 2014-05-29 2020-06-30 Toyota Jidosha Kabushiki Kaisha Driving support apparatus
US11364901B2 (en) 2014-05-29 2022-06-21 Toyota Jidosha Kabushiki Kaisha Driving support apparatus
US11738744B2 (en) 2014-05-29 2023-08-29 Toyota Jidosha Kabushiki Kaisha Driving support apparatus
JP7438718B2 (en) 2019-11-12 2024-02-27 ヤンマーパワーテクノロジー株式会社 Automated driving system for work vehicles

Also Published As

Publication number Publication date
JP3601354B2 (en) 2004-12-15

Similar Documents

Publication Publication Date Title
US11441907B2 (en) Positioning device and positioning method
JP5270184B2 (en) Satellite navigation / dead reckoning integrated positioning system
US7136015B2 (en) Method and apparatus for satellite positioning
JP5521531B2 (en) Position calculation method and position calculation system
JP3448976B2 (en) Current position detection device for vehicles
US11079494B2 (en) Positioning device
JP2012093240A (en) Positioning device for mobile object and car navigation apparatus
JP2000310538A (en) Position detecting apparatus
US9423486B2 (en) Position calculating method and position calculating device
US20140104101A1 (en) Position calculating method and position calculating device
JP2009133702A (en) Positioning device
JP4905054B2 (en) Mobile satellite radio receiver
JP2008157705A (en) Navigation system and gps position accuracy determination method
JP4931113B2 (en) Own vehicle position determination device
JP2007225459A (en) On-board unit
JP2008232771A (en) Positioning device
JP2008051572A (en) Navigation apparatus, method therefor, and program therefor
JP2004150852A (en) Satellite signal receiver
JP5163511B2 (en) GNSS receiver and positioning method
JP3210502B2 (en) Positioning method by satellite navigation
JP3827598B2 (en) Moving body position measurement system
JP3731686B2 (en) Position calculation device
JP2011128138A (en) Method and system for calculating position
JP2004045126A (en) Satellite signal receiver
JP3648935B2 (en) GPS receiver

Legal Events

Date Code Title Description
A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20040622

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040720

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20040729

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20040913

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101001

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111001

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121001

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121001

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131001

Year of fee payment: 9

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term