JP2009098127A - Navigator - Google Patents

Navigator Download PDF

Info

Publication number
JP2009098127A
JP2009098127A JP2008218126A JP2008218126A JP2009098127A JP 2009098127 A JP2009098127 A JP 2009098127A JP 2008218126 A JP2008218126 A JP 2008218126A JP 2008218126 A JP2008218126 A JP 2008218126A JP 2009098127 A JP2009098127 A JP 2009098127A
Authority
JP
Japan
Prior art keywords
acceleration
value
ins
offset
sensor
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2008218126A
Other languages
Japanese (ja)
Inventor
Hideki Sato
秀樹 佐藤
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.)
Yamaha Corp
Original Assignee
Yamaha 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 Yamaha Corp filed Critical Yamaha Corp
Priority to JP2008218126A priority Critical patent/JP2009098127A/en
Priority to US12/235,226 priority patent/US8108140B2/en
Priority to EP08016699.4A priority patent/EP2042832B1/en
Publication of JP2009098127A publication Critical patent/JP2009098127A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Navigation (AREA)
  • Traffic Control Systems (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a navigator which can detect positions with sufficient accuracy, even if it is the case where an offset error of an acceleration sensor exists. <P>SOLUTION: The navigator includes: an acceleration determination means which compares the acceleration obtained by the acceleration sensor with a predetermined threshold value; and a location calculation means which calculates the present location based on the acceleration and past location and moving velocity which are memorized by a memory means when it is determined that an absolute value of the acceleration is larger than the threshold value by the acceleration determination means and which calculates the present location based on an azimuth through the measurement by an azimuth sensor assuming as a constant velocity motion and based on the past location and moving velocity which are memorized in the memory means when it is determined that the absolute value of acceleration is smaller than the threshold value by the acceleration determination means. <P>COPYRIGHT: (C)2009,JPO&INPIT

Description

本発明は、現在位置の案内を行うナビゲーション装置に関する。   The present invention relates to a navigation apparatus for guiding a current position.

車両の位置を計測し走行案内を行うナビゲーション装置には、車内に配線等を行って設置する据え置き型のものと、PND(Personal Navigation Device)と呼ばれる持ち運び自在の携帯型のものとがある。近年、据え置き型のナビゲーション装置(例えば、特許文献1参照)では、GPS信号に基づき取得される位置情報の他、内蔵するジャイロセンサまたは地磁気センサより得られる方位や車両が備える車速計より得られる車速の情報も併用することによって、車両の現在位置を高精度に求める機能を搭載したものが実用化されている。このような機能は、特に、GPS信号が取得できないか、又はGPS信号の精度が良くない場合に有用となる。   There are two types of navigation devices that measure the position of a vehicle and provide travel guidance: a stationary type that is installed by wiring in the vehicle, and a portable type that is called a PND (Personal Navigation Device). In recent years, in a stationary navigation apparatus (see, for example, Patent Document 1), in addition to position information acquired based on GPS signals, a vehicle speed obtained from a built-in gyro sensor or a geomagnetic sensor and a vehicle speedometer provided in the vehicle. By using this information in combination, a device equipped with a function for obtaining the current position of the vehicle with high accuracy has been put into practical use. Such a function is particularly useful when the GPS signal cannot be acquired or the accuracy of the GPS signal is not good.

一方、携帯型のナビゲーション装置では、車速の情報を車両の計器盤等の内部に設置された車速計から得ることが簡単でないため、GPS信号が取得できない等の場合のための補助的手段として、従来、内蔵する加速度センサを用いて計測した加速度を積分することにより車速を得、これにより位置を求めるという方法が採用されている。
特開平9−42979号公報
On the other hand, in a portable navigation device, it is not easy to obtain vehicle speed information from a vehicle speedometer installed inside a vehicle instrument panel, etc., so as an auxiliary means for cases where GPS signals cannot be obtained, etc. Conventionally, a method has been adopted in which vehicle speed is obtained by integrating acceleration measured using a built-in acceleration sensor, thereby obtaining a position.
JP-A-9-42979

しかし、加速度センサは一般にオフセット誤差を有するため、位置の計測において十分に良好な精度が得られないという重大な問題がある。ここで、オフセット誤差とは、本来加速度運動による加速度がゼロであるにもかかわらずある有限の加速度値が加速度センサから出力されてしまう誤差のことである。こうしたオフセット誤差は、例えば、加速度センサが車両に対して水平に取り付けられていない場合や、車両が坂道など水平面から傾いた場所を走行する場合等に、重力加速度を車両に加えられた加速度と誤認してしまうことによって生じ得る。特にこの後者の場合は、たとえ加速度センサが車両に対して水平に取り付けてあったとしてもオフセット誤差の発生が避けられないため、大きな問題となる。   However, since an acceleration sensor generally has an offset error, there is a serious problem that a sufficiently good accuracy cannot be obtained in position measurement. Here, the offset error is an error in which a certain finite acceleration value is output from the acceleration sensor even though the acceleration due to the acceleration motion is originally zero. For example, when the acceleration sensor is not mounted horizontally with respect to the vehicle, or when the vehicle travels on a slope or other location inclined from a horizontal plane, the gravitational acceleration is mistaken as the acceleration applied to the vehicle. It can be caused by doing so. In particular, in the latter case, even if the acceleration sensor is mounted horizontally with respect to the vehicle, the occurrence of an offset error is inevitable, which is a serious problem.

本発明は上記の点に鑑みてなされたものであり、その目的は、加速度センサのオフセット誤差が存在する場合であっても精度良く位置を求めることが可能なナビゲーション装置を提供することにある。   The present invention has been made in view of the above points, and an object of the present invention is to provide a navigation device that can accurately obtain a position even when an offset error of an acceleration sensor exists.

本発明は上記の課題を解決するためになされたものであり、現在位置を求め位置案内を行うナビゲーション装置において、加速度センサから得られる加速度と所定の閾値とを比較する加速度判定手段と、前記加速度判定手段により加速度の絶対値が前記閾値より大きいと判定された場合は、該加速度と記憶手段に記憶されている過去の位置および移動速度とに基づいて現在の位置を算出し、前記加速度判定手段により加速度の絶対値が前記閾値より小さいと判定された場合は、等速運動であるとみなして方位センサによる計測に基づく方位と記憶手段に記憶されている過去の位置および移動速度とに基づいて現在の位置を算出する位置算出手段と、を備えることを特徴とする。   The present invention has been made to solve the above-described problem. In a navigation device that obtains a current position and performs position guidance, acceleration determination means that compares acceleration obtained from an acceleration sensor with a predetermined threshold value, and the acceleration When the determination means determines that the absolute value of the acceleration is greater than the threshold, the current position is calculated based on the acceleration and the past position and moving speed stored in the storage means, and the acceleration determination means If the absolute value of acceleration is determined to be smaller than the threshold value, it is assumed that the motion is constant velocity and based on the azimuth based on the measurement by the azimuth sensor and the past position and moving speed stored in the storage means. And a position calculation means for calculating the current position.

この発明によれば、測定された加速度の絶対値が所定の閾値より小さい場合は、その加速度は車両が傾いている(例えば坂道を走行している等)ことに起因するものであると想定して加速度をゼロとみなし、車両の動きを等速運動として現在の位置を求めるようにした。したがって、車両が傾くことにより加速度センサにオフセット誤差が生じていても、精度良く位置を求めることができる。   According to the present invention, when the absolute value of the measured acceleration is smaller than the predetermined threshold, it is assumed that the acceleration is caused by the vehicle leaning (for example, traveling on a slope). The acceleration is assumed to be zero, and the current position is obtained by assuming that the vehicle moves at a constant speed. Therefore, even if an offset error occurs in the acceleration sensor due to the vehicle tilting, the position can be obtained with high accuracy.

また、本発明は、上記のナビゲーション装置において、GPSにより位置を求めるとともに該求めた異なる時刻における複数の位置から移動速度を求める計測手段と、前記計測手段により求められた位置および移動速度と、前記位置算出手段により求められた位置および該位置に基づき計算される移動速度と、をカルマンフィルタに入力して現在の最も確からしい位置に関する最良値を推定する最良値推定手段と、前記最良値推定手段によって推定された最良値に基づき、前記位置算出手段により求められた位置の誤差が前記加速度センサのオフセット誤差に起因したものであるとして該オフセット誤差の分だけ加速度を補正するための加速度オフセット値を算出する加速度オフセット算出手段と、を更に備え、前記加速度判定手段は、前記加速度センサから得られる加速度を前記加速度オフセット算出手段により算出された加速度オフセット値を用いて補正し、該補正された加速度と前記閾値とを比較することを特徴とする。   Further, the present invention provides a navigation device that obtains a position by GPS and obtains a moving speed from a plurality of positions at the obtained different times, a position and a moving speed obtained by the measuring means, A best value estimating means for inputting a position obtained by the position calculating means and a moving speed calculated based on the position to a Kalman filter to estimate a best value relating to the current most probable position; and the best value estimating means Based on the estimated best value, an acceleration offset value for correcting the acceleration by the amount of the offset error is calculated on the assumption that the position error obtained by the position calculation means is caused by the offset error of the acceleration sensor. An acceleration offset calculating means for performing the acceleration determination, Acceleration obtained from the velocity sensor is corrected using the acceleration offset calculated by the acceleration offset calculation means, characterized by comparing the with the corrected acceleration threshold.

この発明によれば、加速度オフセット値を求めて加速度センサから得られる加速度を補正するようにしたので、オフセット誤差が補正され、更に精度良く位置を求めることができる。   According to the present invention, since the acceleration offset value is obtained and the acceleration obtained from the acceleration sensor is corrected, the offset error is corrected and the position can be obtained with higher accuracy.

また、本発明は、上記のナビゲーション装置において、前記加速度センサから得られる加速度の所定時間における変化が予め定められた値以下である場合に、前記加速度センサから得られた加速度の値に応じて加速度オフセット値を算出する加速度オフセット算出手段を更に備え、前記加速度判定手段は、前記加速度センサから得られる加速度を前記加速度オフセット算出手段により算出された加速度オフセット値を用いて補正し、該補正された加速度と前記閾値とを比較することを特徴とする。   According to the present invention, in the above navigation device, when the change in the acceleration obtained from the acceleration sensor at a predetermined time is equal to or less than a predetermined value, the acceleration is obtained according to the acceleration value obtained from the acceleration sensor. The apparatus further comprises acceleration offset calculating means for calculating an offset value, wherein the acceleration determining means corrects the acceleration obtained from the acceleration sensor using the acceleration offset value calculated by the acceleration offset calculating means, and the corrected acceleration And the threshold value are compared.

この発明によれば、加速度センサから得られた加速度の変化から加速度オフセット値を求めて加速度センサから得られる加速度を補正するようにしたので、GPS信号を取得できない、もしくはGPS信号の信頼性が十分でない場合においてもオフセット誤差が補正され、更に精度良く位置を求めることができる。   According to the present invention, since the acceleration offset value is obtained from the change in acceleration obtained from the acceleration sensor and the acceleration obtained from the acceleration sensor is corrected, the GPS signal cannot be acquired or the reliability of the GPS signal is sufficient. Even in the case where it is not, the offset error is corrected, and the position can be obtained with higher accuracy.

本発明によれば、加速度センサのオフセット誤差が存在する場合であっても精度良く位置を求めることが可能である。   According to the present invention, it is possible to obtain the position with high accuracy even when there is an offset error of the acceleration sensor.

以下、図面を参照しながら本発明の実施形態について詳しく説明する。
図1は、本発明の第1の実施形態によるナビゲーション装置1の構成図である。ナビゲーション装置1は、CPU10と、GPS部11と、加速度センサ12と、磁気センサ13と、メモリ14と、外部記憶装置15と、通信部16と、表示装置17と、音声出力装置18と、を含んで構成され、各部はバスライン19によって接続されている。ナビゲーション装置1のこれら各部は一つの筐体内に収容されており、ナビゲーション装置1が携帯可能であるとともに車両内の任意の位置に取付け可能になっている。以降、ナビゲーション装置1は車両内に取り付けられているものとする。
Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
FIG. 1 is a configuration diagram of a navigation apparatus 1 according to the first embodiment of the present invention. The navigation device 1 includes a CPU 10, a GPS unit 11, an acceleration sensor 12, a magnetic sensor 13, a memory 14, an external storage device 15, a communication unit 16, a display device 17, and an audio output device 18. Each part is connected by a bus line 19. These parts of the navigation device 1 are accommodated in a single casing, and the navigation device 1 can be carried and attached to any position in the vehicle. Henceforth, the navigation apparatus 1 shall be attached in the vehicle.

CPU(中央処理装置)10は、メモリ14に格納されたプログラムを読み込み、該読み込んだプログラムに従って、ナビゲーション装置1の各部を統括し制御するとともに、GPS部11や加速度センサ12、磁気センサ13より取得されるデータから現在の位置を算出する後述(図2〜図4)の演算処理や、ユーザに現在位置を案内するナビゲーション機能を実現する情報処理を実行する。   The CPU (central processing unit) 10 reads a program stored in the memory 14, controls and controls each part of the navigation device 1 according to the read program, and obtains it from the GPS unit 11, the acceleration sensor 12, and the magnetic sensor 13. The calculation processing described later (FIGS. 2 to 4) for calculating the current position from the data to be processed and the information processing for realizing the navigation function for guiding the current position to the user are executed.

GPS部11は、GPS(全地球測位システム)方式によって車両の位置(緯度経度)を測位するとともに、測位された最新の位置とその直前の位置とに基づいて車両の速度および進行方向を算出し、得られた位置、速度、進行方向を表す各情報をCPU10へ出力する。GPS方式では、GPS衛星からのGPS信号を受信し、該GPS信号を受信した受信時刻(図示しない計時手段により計時される)のデータを用いた演算処理によって、位置を測位する。   The GPS unit 11 measures the position (latitude and longitude) of the vehicle by a GPS (Global Positioning System) method, and calculates the speed and traveling direction of the vehicle based on the latest position and the position immediately before the position. Each information indicating the obtained position, speed, and traveling direction is output to the CPU 10. In the GPS system, a GPS signal from a GPS satellite is received, and the position is measured by arithmetic processing using data of a reception time (timed by a time measuring unit (not shown)) when the GPS signal is received.

加速度センサ12は、車両(即ちナビゲーション装置1)の加速度を測定し、その測定結果を表す加速度情報をCPU10へ出力する。測定される加速度は、車両が加速運動を行ったことによる加速度以外に、車両が坂道などで傾いたことによる重力加速度も含んでいる。加速度センサ12としては、互いに直角な異なる二方向の加速度成分を検出する2軸加速度センサを用いるものとし、本実施形態においては、加速度センサ12は、その一方の軸(これをy軸とし、他方をx軸とする)が次に説明する磁気センサ13の基準軸の向きと一致するようにして配置するものとする。CPU10へ出力される加速度情報は、各軸方向の(即ち、磁気センサ13の基準軸を基準としたxy座標系における)加速度成分である。各成分には、重力加速度が測定されてしまうことに起因するオフセット誤差が含まれている。   The acceleration sensor 12 measures the acceleration of the vehicle (that is, the navigation device 1) and outputs acceleration information representing the measurement result to the CPU 10. The measured acceleration includes gravitational acceleration due to the vehicle leaning on a slope or the like in addition to acceleration due to the vehicle performing acceleration motion. As the acceleration sensor 12, a two-axis acceleration sensor that detects acceleration components in two different directions perpendicular to each other is used. In the present embodiment, the acceleration sensor 12 has one axis (this is the y-axis and the other is the other). (X is the x-axis) is aligned with the orientation of the reference axis of the magnetic sensor 13 described below. The acceleration information output to the CPU 10 is an acceleration component in each axis direction (that is, in the xy coordinate system with reference to the reference axis of the magnetic sensor 13). Each component includes an offset error resulting from the gravitational acceleration being measured.

磁気センサ13は、車両(即ちナビゲーション装置1)における地磁気を測定し、磁気センサ13に定められた所定の基準軸が向けられている方位を算出してその方位を表す方位情報をCPU10へ出力する。本実施形態において、磁気センサ13は、この基準軸が車両正面方向に対してある角度だけずれた方向(これは、携帯可能なナビゲーション装置1を車内に置いたときにその置き方が任意であることに基づく)を向くようにして配置される。磁気センサ13としては、互いに直角な異なる二方向の磁気成分を検出する2軸磁気センサ、または互いに直角な異なる三方向の磁気成分を検出する3軸磁気センサ、のいずれかを用いるものとする。例えば2軸磁気センサの場合には、磁気の大きさを検出する2つの素子(磁気抵抗素子)が異なる二方向を向いて配置されており、それぞれの素子が検出する地磁気の大きさに基づいて上記基準軸の方位を算出する。   The magnetic sensor 13 measures the geomagnetism in the vehicle (that is, the navigation device 1), calculates the direction in which the predetermined reference axis defined by the magnetic sensor 13 is directed, and outputs the direction information indicating the direction to the CPU 10. . In the present embodiment, the magnetic sensor 13 has a direction in which the reference axis is deviated by a certain angle with respect to the vehicle front direction (this is arbitrary when the portable navigation device 1 is placed in the vehicle). Arranged based on the As the magnetic sensor 13, either a two-axis magnetic sensor that detects magnetic components in two different directions perpendicular to each other or a three-axis magnetic sensor that detects magnetic components in three different directions perpendicular to each other is used. For example, in the case of a biaxial magnetic sensor, two elements (magnetoresistive elements) for detecting the magnitude of magnetism are arranged in two different directions, and based on the magnitude of geomagnetism detected by each element. The direction of the reference axis is calculated.

メモリ14は、ROM(読出し専用メモリ)とRAM(随時書込み読出しメモリ)からなる。ROMには、CPU10で実行されるプログラムが記憶されている。RAMは、プログラム実行時に一時データを記憶する記憶領域として使用される。
外部記憶装置15は、例えばハードディスクドライブなどの大容量記憶装置であり、ナビゲーションに必要な地図情報等を記憶している。
通信部16は、例えば無線通信等によってインターネットと接続し、最新の地図情報等をインターネットを介して取得する。取得した地図情報等は、外部記憶装置15に記憶される。
The memory 14 includes a ROM (read only memory) and a RAM (anytime read / write memory). In the ROM, a program executed by the CPU 10 is stored. The RAM is used as a storage area for storing temporary data during program execution.
The external storage device 15 is a large-capacity storage device such as a hard disk drive, and stores map information and the like necessary for navigation.
The communication unit 16 is connected to the Internet, for example, by wireless communication, and acquires the latest map information and the like via the Internet. The acquired map information and the like are stored in the external storage device 15.

表示装置17は、CPU10によって実現するナビゲーション機能による走行案内を、例えば地図上に車両の現在位置と方位が分かるような形でグラフィカルに表示する。
音声出力装置18は、必要なタイミングで例えば交差点を曲がる等の音声案内用の音声を発生してスピーカから音声案内を行う。
The display device 17 graphically displays the travel guidance by the navigation function realized by the CPU 10 in a form that shows the current position and direction of the vehicle on a map, for example.
The voice output device 18 generates voice guidance voice such as turning an intersection at a necessary timing and performs voice guidance from a speaker.

次に、ナビゲーション装置1がGPS部11、加速度センサ12、および磁気センサ13より取得されるデータから現在の位置を算出する処理について説明する。図2は、ナビゲーション装置1の当該処理にかかる機能ブロックを示す図であり、当該処理は、INS系位置算出部101(加速度センサ12と磁気センサ13をINS(Inertial Navigation System)系と称する)と、最良値推定部102と、位置補正部103と、加速度オフセット算出部104と、加速度判定部105とによって実現される。なお、これら各部の機能は、前述したように、メモリ14に格納されたプログラムをCPU10が読み出して実行することにより実現される。   Next, processing in which the navigation apparatus 1 calculates the current position from data acquired from the GPS unit 11, the acceleration sensor 12, and the magnetic sensor 13 will be described. FIG. 2 is a diagram showing functional blocks related to the process of the navigation device 1, and this process is an INS system position calculation unit 101 (the acceleration sensor 12 and the magnetic sensor 13 are referred to as an INS (Internal Navigation System) system). The best value estimation unit 102, the position correction unit 103, the acceleration offset calculation unit 104, and the acceleration determination unit 105 are realized. Note that the functions of these units are realized by the CPU 10 reading and executing the program stored in the memory 14 as described above.

INS系位置算出部101へは、加速度センサ12からオフセット誤差を含んだ加速度a=(ax,ay)が入力され(axはx軸成分、ayはy軸成分)、磁気センサ13から上記基準軸の示す方位θが入力され、GPS部11から車両の進行方向φが入力される。更に、INS系測位部101へは、加速度オフセット算出部104から、加速度センサ12からの加速度aに含まれるオフセット誤差を補正するための加速度オフセット値b=(bx,by)が入力される(bxはx軸成分、byはy軸成分)。
INS系位置算出部101は、加速度判定部105からの指示に応じて、上記の入力された各データに基づき、車両の位置PINSおよび速度VINSを計算する。また、その計算結果を最良値推定部102と位置補正部103へ出力する。
An acceleration a = (ax, ay) including an offset error is input from the acceleration sensor 12 to the INS system position calculation unit 101 (ax is an x-axis component, ay is a y-axis component), and the reference axis Is input, and the traveling direction φ of the vehicle is input from the GPS unit 11. Further, the acceleration offset value b = (bx, by) for correcting the offset error included in the acceleration a from the acceleration sensor 12 is input from the acceleration offset calculation unit 104 to the INS system positioning unit 101 (bx Is the x-axis component, and by is the y-axis component).
In response to an instruction from the acceleration determination unit 105, the INS system position calculation unit 101 calculates a vehicle position P INS and a speed V INS based on each of the input data. The calculation result is output to the best value estimation unit 102 and the position correction unit 103.

具体的に、INS系位置算出部101は、加速度判定部105において加速度a+bの絶対値が所定の閾値以上であると判定された場合は、位置PINSおよび速度VINSを次式に従って計算する。
aN=−(ax+bx)・sinθ+(ay+by)・cosθ ……(1A)
aE= (ax+bx)・cosθ+(ay+by)・sinθ ……(1B)
VNINS=VN0INS+aN・Δt ……(1C)
VEINS=VE0INS+aE・Δt ……(1D)
PNINS=PN0INS+sin−1[(VNINS×Δt)/R] ……(1E)
PEINS=PE0INS+sin−1[(VEINS×Δt)/{R×cos(PN0INS)}] ……(1F)
但し、aN,aEはそれぞれ加速度a+bの緯度および経度方向成分、VNINS,VEINSはそれぞれ速度VINSの緯度および経度方向成分、PNINS,PEINSはそれぞれ位置PINSの緯度および経度成分である。また、VN0INS,VE0INS,PN0INS,PE0INSはそれぞれ前回計算されたVNINS,VEINS,PNINS,PEINSの値、Δtは前回の計算から今回の計算までの経過時間、Rは地球の半径である。
Specifically, if the acceleration determination unit 105 determines that the absolute value of the acceleration a + b is equal to or greater than a predetermined threshold, the INS system position calculation unit 101 calculates the position P INS and the velocity V INS according to the following equations.
aN = − (ax + bx) · sin θ + (ay + by) · cos θ (1A)
aE = (ax + bx) · cos θ + (ay + by) · sin θ (1B)
VN INS = VN0 INS + aN · Δt (1C)
VE INS = VE0 INS + aE · Δt (1D)
PN INS = PN0 INS + sin −1 [(VN INS × Δt) / R] (1E)
PE INS = PE0 INS + sin −1 [(VE INS × Δt) / {R × cos (PN0 INS )}] (1F)
Where aN and aE are the latitude and longitude components of acceleration a + b, VN INS and VE INS are the latitude and longitude components of velocity V INS , and PN INS and PE INS are the latitude and longitude components of position P INS , respectively. . Also, VN0 INS , VE0 INS , PN0 INS , PE0 INS are the previously calculated values of VN INS , VE INS , PN INS , PE INS , Δt is the elapsed time from the previous calculation to this calculation, and R is the earth Is the radius.

一方、INS系位置算出部101は、加速度判定部105において加速度a+bの絶対値が上記所定の閾値より小さいと判定された場合は、位置PINSおよび速度VINSを次式に従って計算する(PNINSとPEINSの式は上記と共通)。
V=√(VN0INS +VE0INS ) ……(2A)
VNINS=V×cosθ’ ……(2B)
VEINS=V×sinθ’ ……(2C)
PNINS=PN0INS+sin−1[(VNINS×Δt)/R ……(2D)
PEINS=PE0INS+sin−1[(VEINS×Δt)/{R×cos(PN0INS)}] ……(2E)
但し、θ’は磁気センサ13が車両正面方向からずれていることを考慮し方位θを後述(図4のフローチャートのステップS52を参照)のようにして補正した値である。
On the other hand, when the acceleration determination unit 105 determines that the absolute value of the acceleration a + b is smaller than the predetermined threshold, the INS system position calculation unit 101 calculates the position P INS and the velocity V INS according to the following equations (PN INS And PE INS are the same as above).
V = √ (VN0 INS 2 + VE0 INS 2 ) (2A)
VN INS = V × cos θ ′ (2B)
VE INS = V × sin θ ′ (2C)
PN INS = PN0 INS + sin −1 [(VN INS × Δt) / R (2D)
PE INS = PE0 INS + sin −1 [(VE INS × Δt) / {R × cos (PN0 INS )}] (2E)
However, θ ′ is a value obtained by correcting the azimuth θ as described later (see step S52 in the flowchart of FIG. 4) in consideration of the fact that the magnetic sensor 13 is deviated from the vehicle front direction.

このように、加速度の絶対値が所定の閾値以上の場合には、加速度aN,aEを用い車両の運動を加速度運動として速度VINSを求めてこれにより位置PINSを算出し、加速度の絶対値が所定の閾値より小さい場合には、前回の速度VN0INS,VE0INSを用い車両の運動を等速運動とみなして速度VINSを求めてこれにより位置PINSを算出する。したがって、車両が実際には等速運動を行っているが走行している場所が坂道などであるために加速度センサ12が出力する加速度aにオフセット誤差が含まれてしまう場合でも、INS系位置算出部101は後者の方式により、即ち等速運動であるとして位置を算出するので、オフセット誤差の影響を受けることなく精度良く位置を求めることができる。 Thus, when the absolute value of the acceleration is equal to or higher than the predetermined threshold value, the acceleration aN, thereby calculates the position P INS seeking speed V INS as acceleration motion the motion of the vehicle using aE, the absolute value of the acceleration Is smaller than the predetermined threshold value, the previous speeds VN0 INS and VE0 INS are used to regard the motion of the vehicle as a constant speed motion, and the speed V INS is obtained, thereby calculating the position P INS . Therefore, even when the vehicle is actually moving at a constant speed but the traveling place is on a slope or the like, an offset error is included in the acceleration a output from the acceleration sensor 12. Since the unit 101 calculates the position by the latter method, that is, assuming that it is a constant velocity motion, the position can be accurately obtained without being affected by the offset error.

なお、等速運動とみなして位置PINSを算出する処理は、GPS部11によりGPS信号を取得できない、もしくはGPS信号の信頼性が十分でない場合にのみ、行うようにしてもよい。GPS信号を取得できている間は後述する最良値推定部102および加速度オフセット算出部104が有効に機能して加速度オフセット値bを算出でき、これにより加速度センサ12のオフセット誤差が補正されるから、加速度aN,aEの値が小さくてもそれを信頼して上記の式(1A)〜(1F)により正しく位置PINSを求めることができる。これに対し、例えばトンネル内のようにGPS信号が取得できなくなると、最良値推定部102および加速度オフセット算出部104が有効に機能しないから加速度センサ12のオフセット誤差を補正できず、加速度aN,aEの値の信頼性がなくなる。このような場合に、加速度の絶対値が小さいとオフセット誤差の影響が大きく出てしまうので、等速運動とみなして位置PINSを算出する上記の式(2A)〜(2E)による処理に切り替えることの大きな意味がある。 It should be noted that the process of calculating the position P INS on the assumption that the motion is constant speed may be performed only when the GPS signal cannot be acquired by the GPS unit 11 or when the reliability of the GPS signal is not sufficient. While the GPS signal can be acquired, the best value estimation unit 102 and the acceleration offset calculation unit 104 described later function effectively to calculate the acceleration offset value b, thereby correcting the offset error of the acceleration sensor 12. Even if the values of the accelerations aN and aE are small, the position P INS can be correctly obtained by the above equations (1A) to (1F) with confidence. On the other hand, if the GPS signal cannot be acquired as in a tunnel, for example, the best value estimation unit 102 and the acceleration offset calculation unit 104 do not function effectively, so that the offset error of the acceleration sensor 12 cannot be corrected, and the accelerations aN, aE The reliability of the value is lost. In such a case, if the absolute value of the acceleration is small, the influence of the offset error is large. Therefore, the processing is switched to the processing according to the above formulas (2A) to (2E) for calculating the position P INS assuming that the motion is constant velocity. It has a big meaning.

最良値推定部102へは、INS系位置算出部101から上記計算された位置PINSおよび速度VINSが入力されるとともに、GPS部11から位置PGPSおよび速度VGPSが入力される。INS系位置算出部101によって計算される位置PINSおよび速度VINSと、GPS部11から得られる位置PGPSおよび速度VGPSとは、それぞれの測定原理が違うので互いの値には差があり、また各値は固有の誤差を持っている。最良値推定部102はこのような入力値から最も確からしい値を推定するために用いられる。 The calculated value P INS and velocity V INS are input from the INS system position calculation unit 101 to the best value estimation unit 102, and the position P GPS and velocity V GPS are input from the GPS unit 11. The position P INS and velocity V INS calculated by the INS system position calculation unit 101 and the position P GPS and velocity V GPS obtained from the GPS unit 11 are different from each other because their measurement principles are different. In addition, each value has an inherent error. The best value estimation unit 102 is used to estimate the most probable value from such input values.

最良値推定部102は、GPS部11からの入力がある度に、入力された位置の差分δP’=PGPS−PINSと入力された速度の差分δV’=VGPS−VINSを計算し、この差分δP’およびδV’に基づいてそれらの最も確からしい値(最良値)であるδPとδVを推定する。このδPとδVの推定は、カルマンフィルタを用いた演算によって行うものとする。演算の詳細については後述する。最良値推定部102は、推定により得られた最良値δPとδVを位置補正部103および加速度オフセット算出部104へ出力する。 The best value estimation unit 102 calculates the input position difference δP ′ = P GPS −P INS and the input speed difference δV ′ = V GPS −V INS every time there is an input from the GPS unit 11. Based on the differences δP ′ and δV ′, δP and δV which are the most probable values (best values) are estimated. It is assumed that δP and δV are estimated by calculation using a Kalman filter. Details of the calculation will be described later. Best value estimation section 102 outputs best values δP and δV obtained by the estimation to position correction section 103 and acceleration offset calculation section 104.

位置補正部103は、INS系位置算出部101から入力される車両の位置PINSと最良値推定部102から入力される位置の差分の最良値δPとに基づいて、次式
P=PINS+δP
に従って位置PINSを補正することにより現在の最も確からしい位置Pを計算し、得られた位置Pをナビゲーションソフト(CPU10)への出力とする。
Based on the vehicle position P INS input from the INS system position calculation unit 101 and the best value δP of the position difference input from the best value estimation unit 102, the position correction unit 103
P = P INS + δP
The current most probable position P is calculated by correcting the position P INS according to, and the obtained position P is used as an output to the navigation software (CPU 10).

加速度オフセット算出部104は、最良値推定部102から入力される位置の差分の最良値δPに基づいて、加速度センサ12から出力される加速度aに含まれるオフセット誤差を補正するための補正値である加速度オフセット値bを算出し、算出した加速度オフセット値bをINS系位置算出部101へ出力する。ここで、位置の差分の最良値δPから加速度オフセット値bを求めるための算出式は、次のような仮定に基づいて定めることとする。   The acceleration offset calculation unit 104 is a correction value for correcting an offset error included in the acceleration a output from the acceleration sensor 12 based on the best value δP of the position difference input from the best value estimation unit 102. The acceleration offset value b is calculated, and the calculated acceleration offset value b is output to the INS system position calculation unit 101. Here, a calculation formula for obtaining the acceleration offset value b from the best position difference δP is determined based on the following assumptions.

即ち、INS系位置算出部101において前回、加速度センサ12からの加速度aが完全に(つまりオフセット誤差がなくなるように)補正されていたとすれば、最良値推定部102によって得られる位置の差分の最良値δPはゼロになるはずである。一方、加速度aの補正が不完全である(つまりオフセット誤差が残っている)ならば、その不完全な分の加速度(これをb’とする)によってINS系位置算出部101で計算される位置PINSに誤差が生じ、この誤差に対応する値を持った位置の差分の最良値δPが最良値推定部102から得られることになるはずである。故に、最良値推定部102の出力δPは、加速度b’での等加速度運動によって生じる位置の誤差であり、次式
δP=(1/2)・b’・ΔT
が成り立つ。
That is, if the INS system position calculation unit 101 previously corrected the acceleration a from the acceleration sensor 12 completely (that is, so as to eliminate the offset error), the best position difference obtained by the best value estimation unit 102 is obtained. The value δP should be zero. On the other hand, if the correction of the acceleration a is incomplete (that is, an offset error remains), the position calculated by the INS system position calculation unit 101 by the incomplete amount of acceleration (referred to as b ′). An error occurs in P INS, and the best value δP of the position difference having a value corresponding to this error should be obtained from the best value estimation unit 102. Therefore, the output δP of the best value estimation unit 102 is a position error caused by the uniform acceleration motion at the acceleration b ′.
δP = (1/2) · b ′ · ΔT 2
Holds.

したがって、加速度オフセット算出部104が加速度オフセット値bを算出するための算出式を
b=b0+2・δPxy/ΔT
とする。但し、ここでb0はINS系位置算出部101での前回の加速度の補正に用いられた加速度オフセット値、ΔTは前回の加速度の補正から今回までの経過時間である。また、δPxyは、最良値推定部102から得られる緯度経度で表された最良値δPを車両に固定された上記のxy座標系に変換(方位θの回転)したものである。
Therefore, the calculation formula for the acceleration offset calculation unit 104 to calculate the acceleration offset value b is
b = b0 + 2 · δPxy / ΔT 2
And Here, b0 is an acceleration offset value used for the previous acceleration correction in the INS system position calculation unit 101, and ΔT is an elapsed time from the previous acceleration correction to the current time. Further, δPxy is obtained by converting the best value δP expressed by the latitude and longitude obtained from the best value estimation unit 102 into the xy coordinate system fixed to the vehicle (rotation of the azimuth θ).

こうして算出された加速度オフセット値bがINS系位置算出部101へ入力され、上述の式(1A),(1B)に示されるように加速度aと加算されることによって、加速度センサ12から得られる加速度aの補正が行われる。
なお、ナビゲーション装置1が加速度オフセット算出部104により加速度オフセット値bを算出し、この加速度オフセット値bを用いてINS系位置算出部101において加速度aを補正する処理は、必要に応じて行えばよく、省略することも可能である(省略する場合、次に述べる加速度判定部105が判定する加速度はa+bではなくaとなる)。但し、当該処理を行うと、加速度センサ12からのオフセット誤差を含んだ加速度が補正された上で更にその補正された加速度が上述した閾値より小さい場合に、等速運動として位置が算出されることになるので、より一層精度良く位置を求めることが可能である。
The acceleration offset value b calculated in this way is input to the INS system position calculation unit 101 and is added to the acceleration a as shown in the above formulas (1A) and (1B), whereby the acceleration obtained from the acceleration sensor 12 is obtained. Correction of a is performed.
The navigation device 1 may calculate the acceleration offset value b by the acceleration offset calculation unit 104 and correct the acceleration a by the INS system position calculation unit 101 using the acceleration offset value b as necessary. (The acceleration determined by the acceleration determination unit 105 described below is a instead of a + b). However, when this processing is performed, the position including the offset error from the acceleration sensor 12 is corrected, and the position is calculated as a constant velocity motion when the corrected acceleration is smaller than the threshold value described above. Therefore, it is possible to obtain the position with higher accuracy.

次に、ナビゲーション装置1の動作の具体例を説明する。図3は、ナビゲーション装置1の動作例を示すフローチャートである。   Next, a specific example of the operation of the navigation device 1 will be described. FIG. 3 is a flowchart showing an operation example of the navigation device 1.

まず、CPU10は、GPS部11から車両の位置PGPSと速度VGPSと進行方向φとを取得してメモリ14に格納し(ステップS1)、更に、加速度センサ12からオフセット誤差を含んだ加速度aを取得するとともに磁気センサ13から基準軸の示す方位θを取得して、それぞれをメモリ14に格納する(ステップS2)。次いでCPU10は、上記取得した加速度aとメモリ14に保持している加速度オフセット値bとの和a+bを計算し、その絶対値が所定の閾値以上であるか否かを判断する(ステップS3)。加速度オフセット値bは、ステップS3の初回実行時は初期値b=(0,0)であり、以下説明するステップS14を実行した後は当該ステップS14で算出された値であり、適宜、メモリ14にこれらが記憶されている。 First, the CPU 10 acquires the position P GPS , the speed V GPS, and the traveling direction φ of the vehicle from the GPS unit 11 and stores them in the memory 14 (step S1). Further, the acceleration a including an offset error from the acceleration sensor 12 is obtained. Are acquired from the magnetic sensor 13 and stored in the memory 14 (step S2). Next, the CPU 10 calculates the sum a + b of the acquired acceleration a and the acceleration offset value b held in the memory 14, and determines whether or not the absolute value is equal to or greater than a predetermined threshold (step S3). The acceleration offset value b is an initial value b = (0, 0) when the step S3 is executed for the first time, and is a value calculated at the step S14 after executing the step S14 described below. These are memorized.

和a+bの絶対値が所定の閾値以上である場合、INS系位置算出部101は、上述の式(1A)〜(1F)に従って車両の位置PINSおよび速度VINSを計算する(ステップS4)。和a+bの絶対値が所定の閾値より小さい場合、INS系位置算出部101は、上述の式(2A)〜(2E)に従って車両の位置PINSおよび速度VINSを計算する(ステップS5)。
なお、上記のステップS3の前に、GPS部11による計測の精度を判断するステップを設け、例えばGPS信号を取得できない場合やGPS信号を取得したがその精度が良好でない場合にのみステップS3を実行し、精度が十分に良好である場合にはステップS4により位置PINSおよび速度VINSを計算するようにしてもよい。これは、GPS部11からの計測結果の精度が良好であれば、以下のステップS14で加速度オフセット値bが算出され、この加速度オフセット値bにより補正された加速度を用いて上記の式(1A)〜(1F)から計算される位置PINSは、比較的信頼性が高いと期待できるからである。
When the absolute value of the sum a + b is equal to or greater than a predetermined threshold, the INS system position calculation unit 101 calculates the vehicle position P INS and the speed V INS according to the above-described equations (1A) to (1F) (step S4). When the absolute value of the sum a + b is smaller than the predetermined threshold, the INS system position calculation unit 101 calculates the vehicle position P INS and the speed V INS in accordance with the above-described equations (2A) to (2E) (step S5).
Before step S3, a step for determining the accuracy of measurement by the GPS unit 11 is provided. For example, step S3 is executed only when the GPS signal cannot be acquired or when the GPS signal is acquired but the accuracy is not good. If the accuracy is sufficiently good, the position P INS and the speed V INS may be calculated in step S4. If the accuracy of the measurement result from the GPS unit 11 is good, the acceleration offset value b is calculated in the following step S14, and the above equation (1A) is used using the acceleration corrected by the acceleration offset value b. This is because the position P INS calculated from (1F) can be expected to be relatively reliable.

次に、CPU10は、位置PGPSと速度VGPSが更新されたか否かを判断する(ステップS6)。後述するように上記ステップS1は時間ΔtGPS間隔で繰り返され、上記ステップS2は時間ΔtINS(<ΔtGPS)間隔で繰り返されるが、まだ時間ΔtGPSが経過していない時は位置PGPSと速度VGPSは更新されないので、ステップS6の判断はNOとなる。 Next, the CPU 10 determines whether or not the position P GPS and the speed V GPS have been updated (step S6). As will be described later, step S1 is repeated at time Δt GPS intervals, and step S2 is repeated at time Δt INS (<Δt GPS ) intervals. When time Δt GPS has not yet elapsed, position P GPS and speed are repeated. Since V GPS is not updated, the determination in step S6 is NO.

この場合、位置補正部103は、ステップS4またはステップS5で計算された位置PINSをメモリ14に保持している最新の最良値δPを用いて補正することにより現在の最も確からしい車両の位置Pを計算する(ステップS7)。なお、最良値δPは、以下説明するステップS13の前回の実行時に算出されてメモリ14に記憶されている。こうして、位置補正部103から、INS系位置算出部101で計算された位置PINSを補正した位置Pがナビゲーションソフトへ出力される(ステップS8)。 In this case, the position correcting unit 103 corrects the position P INS calculated in step S4 or step S5 using the latest best value δP held in the memory 14, thereby correcting the current most probable vehicle position P. Is calculated (step S7). The best value δP is calculated and stored in the memory 14 at the previous execution of step S13 described below. In this way, the position correction unit 103 outputs the position P obtained by correcting the position P INS calculated by the INS system position calculation unit 101 to the navigation software (step S8).

続いて、CPU10は、時間ΔtGPSが経過したか否かを判断し(ステップS9)、時間ΔtGPSが経過すると再びステップS1からの処理を繰り返す。時間ΔtGPSがまだ経過していなければ、CPU10は、次に時間ΔtINSが経過したか否かを判断し(ステップS10)、時間ΔtINSが経過すると再びステップS2からの処理を繰り返す。このようにして、GPS部11からのデータは時間ΔtGPS間隔で更新されるとともに、INS系のデータは時間ΔtINS間隔で更新されていく。なお、一例として、ΔtGPS=1秒、ΔtINS=0.2秒とする。 Subsequently, CPU 10 determines whether the elapsed time Delta] t GPS (step S9), and the process is repeated from step S1 once again the time Delta] t GPS has elapsed. If not elapsed time Delta] t GPS is still, CPU 10 determines whether or not next time Delta] t INS has elapsed (step S10), and repeats the process from step S2 again when the time Delta] t INS has elapsed. In this way, the data from the GPS unit 11 is updated at the time Δt GPS interval, and the INS data is updated at the time Δt INS interval. As an example, Δt GPS = 1 second and Δt INS = 0.2 seconds.

GPS部11からのデータが更新されると、ステップS6の判断はYESとなる。この場合、CPU10は、GPS部11による計測の精度が十分に良好であるか否かを判定する(ステップS11)。例えば、GPS部11から得られる車両の速度VGPSが所定値以上であれば精度が十分に良好であると判定するようにする。精度が十分に良好であると判定された場合、CPU10は、ステップS1で取得した進行方向φとステップS2で取得した磁気センサ13の基準軸の方位θの差分δ=θ−φを計算し、得られた差分δの値を配列Aに追加する(ステップS12)。この差分δは、磁気センサ13から得られる方位θが車両の進行方向φからどれだけずれているかを表している。 When the data from the GPS unit 11 is updated, the determination in step S6 is YES. In this case, the CPU 10 determines whether or not the accuracy of measurement by the GPS unit 11 is sufficiently good (step S11). For example, if the vehicle speed V GPS obtained from the GPS unit 11 is greater than or equal to a predetermined value, it is determined that the accuracy is sufficiently good. When it is determined that the accuracy is sufficiently good, the CPU 10 calculates a difference δ = θ−φ between the traveling direction φ acquired in step S1 and the orientation θ of the reference axis of the magnetic sensor 13 acquired in step S2, The obtained value of difference δ is added to array A (step S12). This difference δ represents how much the azimuth θ obtained from the magnetic sensor 13 deviates from the traveling direction φ of the vehicle.

次いで、最良値推定部102は、ステップS1,S4,S5で得られた位置および速度から計算される差分δP’=PGPS−PINSとδV’=VGPS−VINSに基づいて、それらの最良値δPとδVを推定する(ステップS13)。この推定結果はメモリ14に記憶される。加速度オフセット算出部104は、上記得られた最良値δPから加速度オフセット値bを算出する(ステップS14)。算出された加速度オフセット値bは、メモリ14に記憶されステップS3において加速度a+bの閾値判定のために使用される。ステップS14の後、位置補正部103はステップS13で推定された最良値δPを用いて現在の最も確からしい車両の位置Pを計算する(上記ステップS7)。 Next, the best value estimator 102 calculates the difference δP ′ = P GPS −P INS and δV ′ = V GPS −V INS calculated from the positions and velocities obtained in steps S1, S4 and S5. The best values δP and δV are estimated (step S13). This estimation result is stored in the memory 14. The acceleration offset calculation unit 104 calculates the acceleration offset value b from the obtained best value δP (step S14). The calculated acceleration offset value b is stored in the memory 14 and used for threshold value determination of the acceleration a + b in step S3. After step S14, the position correcting unit 103 calculates the current most probable vehicle position P using the best value δP estimated in step S13 (step S7).

一方、ステップS11においてGPS部11による計測の精度が十分に良好でないと判定された場合には、ステップS13における最良値δPとδVの推定およびステップS14における加速度オフセット値bの算出を行うことなく、ステップS7へ処理を移す。これは、GPS部11からの精度の悪いデータに基づいて精度の悪い最良値δPが推定されてしまう結果、位置補正部103が出力する位置Pや加速度オフセット算出部104が出力する加速度オフセット値bの精度が悪くなってしまうことがないようにするためである。   On the other hand, if it is determined in step S11 that the accuracy of measurement by the GPS unit 11 is not sufficiently good, the best values δP and δV in step S13 and the acceleration offset value b in step S14 are not calculated. The process moves to step S7. This is because the inaccurate best value δP is estimated based on inaccurate data from the GPS unit 11, and as a result, the position P output from the position correcting unit 103 and the acceleration offset value b output from the acceleration offset calculating unit 104. This is to prevent the accuracy of the process from deteriorating.

次に、上述したステップS5の具体的な処理を説明する。図4は、加速度a+bの絶対値が所定の閾値より小さい場合にINS系位置算出部101が位置および速度を算出する処理の例を示すフローチャートである。
INS系位置算出部101は、ステップS12で差分δの値が格納された配列A内の各差分δ(異なる時刻に対応する複数の差分δ)の平均δaveを計算し(ステップS51)、この平均値δaveを用いて、ステップS2で取得された磁気センサ13からの方位θを次式
θ’=θ−δave
に従って補正する(ステップS52)。これにより、方位θは車両の進行方向を表すθ’に補正され、磁気センサ13の出力θから車両の進行方向が求まることになる。この補正後の方位θ’を用いて、INS系位置算出部101は上述の式(2B)および(2C)に従って速度VINSを計算し(ステップS53)、更に式(2D)および(2E)に従って位置PINSを計算する(ステップS54)。
Next, the specific process of step S5 mentioned above is demonstrated. FIG. 4 is a flowchart illustrating an example of processing in which the INS system position calculation unit 101 calculates the position and speed when the absolute value of the acceleration a + b is smaller than a predetermined threshold.
The INS system position calculation unit 101 calculates an average δave of each difference δ (a plurality of differences δ corresponding to different times) in the array A in which the value of the difference δ is stored in step S12 (step S51). Using the value δave, the azimuth θ from the magnetic sensor 13 obtained in step S2 is
θ ′ = θ−δave
(Step S52). As a result, the bearing θ is corrected to θ ′ representing the traveling direction of the vehicle, and the traveling direction of the vehicle is obtained from the output θ of the magnetic sensor 13. Using this corrected azimuth θ ′, INS system position calculation unit 101 calculates speed V INS according to the above-described equations (2B) and (2C) (step S53), and further according to equations (2D) and (2E). The position P INS is calculated (step S54).

次に、最良値推定部102が最良値δPとδVを推定する演算処理の詳細について説明する。
最良値推定部102の演算処理にはカルマンフィルタを適用し、本実施形態では、カルマンフィルタへの入力となる測定ベクトルz、その出力である状態ベクトルx、状態ベクトルxの時間変化を表す状態方程式を次のように定義する。
Next, details of the calculation process in which the best value estimation unit 102 estimates the best values δP and δV will be described.
A Kalman filter is applied to the arithmetic processing of the best value estimation unit 102. In the present embodiment, a measurement vector z that is an input to the Kalman filter, a state vector x that is an output thereof, and a state equation that represents a time change of the state vector x are Define as follows.

Figure 2009098127
Figure 2009098127

ここで、PNGPS,PEGPSはそれぞれGPS部11から得られる位置PGPSの緯度および経度成分、VNGPS,VEGPSはそれぞれGPS部11から得られる速度VGPSの緯度および経度方向成分である。また、δPN,δPEはそれぞれ最良値δPの緯度および経度成分、δVN,δVEはそれぞれ最良値δVの緯度および経度方向成分である。なお添字INSの変数はINS系のものである(既述)。 Here, PN GPS and PE GPS are the latitude and longitude components of the position P GPS obtained from the GPS unit 11, and VN GPS and VE GPS are the latitude and longitude components of the velocity V GPS obtained from the GPS unit 11, respectively. Further, δPN and δPE are the latitude and longitude components of the best value δP, respectively, and δVN and δVE are the latitude and longitude components of the best value δV, respectively. Note that the variable of the subscript INS is of the INS system (described above).

このような定義の下、カルマンフィルタの出力である状態ベクトルxは、次式のように表される。但し、k(=1,2,…)は最良値の推定がk回目であることを示す添字であり、時間に対応している。   Under such a definition, the state vector x, which is the output of the Kalman filter, is expressed as follows. However, k (= 1, 2,...) Is a subscript indicating that the best value is estimated for the kth time, and corresponds to time.

Figure 2009098127
Figure 2009098127

ここで、Kはカルマンゲインであり、次式で計算される。 Here, K k is a Kalman gain and is calculated by the following equation.

Figure 2009098127
Figure 2009098127

但し、Δtは前回(k−1回目)推定した最良値δPに基づき算出される加速度オフセット値bによって加速度の補正を行ってから今回(k回目)までの経過時間であり、各行列R,P,Qの初期値は次のとおりとする。   However, Δt is an elapsed time from the correction of the acceleration by the acceleration offset value b calculated based on the best value δP estimated at the previous time (k−1th time) to the current time (k time), and each matrix R, P The initial values of Q are as follows.

Figure 2009098127
Figure 2009098127

最良値推定部102は、これらの式に従って測定ベクトルzから状態ベクトルxを計算する。この計算により求められる状態ベクトルxが、最良値δP(の緯度および経度成分)とδV(の緯度および経度方向成分)である。   The best value estimation unit 102 calculates the state vector x from the measurement vector z according to these equations. The state vector x obtained by this calculation is the best value δP (its latitude and longitude components) and δV (its latitude and longitude components).

以上、図面を参照してこの発明の一実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、この発明の要旨を逸脱しない範囲内において様々な設計変更等をすることが可能である。
例えば、本発明は携帯型のナビゲーション装置に用いて好適であるが、据え置き型のナビゲーション装置にも同様に適用可能である。また、上記実施形態において、車両とは自動車、自動二輪車、自転車等を含むものであるとする。
As described above, the embodiment of the present invention has been described in detail with reference to the drawings. However, the specific configuration is not limited to the above, and various design changes and the like can be made without departing from the scope of the present invention. It is possible to
For example, the present invention is suitable for use in a portable navigation device, but is also applicable to a stationary navigation device. In the above embodiment, the vehicle includes an automobile, a motorcycle, a bicycle, and the like.

また、加速度オフセット算出部104が加速度オフセット値bを算出するための算出式は、次式
b=b0+2・δPxy/ΔT・r
のようにしてもよい。ここで、rは加速度オフセット値を前回の値b0から更新する際の変更の度合を調節するための係数であり、例えばr=0.03のように適宜設定する。
また、加速度オフセット算出部104が加速度オフセット値bを算出するための算出式は、最良値推定部102から得られる速度の最良値δVを用いて、次式
b=b0+δVxy/ΔT
のようにしてもよい。但し、δVxyは緯度経度で表された最良値δVを車両に固定された上記のxy座標系に変換(方位θの回転)したものである。
The calculation formula for the acceleration offset calculation unit 104 to calculate the acceleration offset value b is as follows:
b = b0 + 2 · δPxy / ΔT 2 · r
It may be as follows. Here, r is a coefficient for adjusting the degree of change when the acceleration offset value is updated from the previous value b0, and is set as appropriate, for example, r = 0.03.
The calculation formula for the acceleration offset calculation unit 104 to calculate the acceleration offset value b is expressed by the following equation using the best speed value δV obtained from the best value estimation unit 102:
b = b0 + δVxy / ΔT
It may be as follows. However, δVxy is obtained by converting the best value δV represented by latitude and longitude into the xy coordinate system fixed to the vehicle (rotation of the azimuth θ).

また、最良値推定部102は、位置および速度の差分δP’=PGPS−PINS,δV’=VGPS−VINSを計算することなく、これらの位置PGPSおよびPINSと速度VGPSおよびVINSから直接、位置の最良値Pと速度の最良値Vを推定するようにしてもよい。この場合、推定を行うためのカルマンフィルタの演算において、入力となる測定ベクトルzおよび出力である状態ベクトルxを適宜修正するとともに、計算された状態ベクトルxにより与えられる最良値Pをそのままナビゲーションソフト(CPU10)への出力とすればよい。
また、磁気センサ13に代えてジャイロセンサなどの他の方位センサを用いてもよい。
The best value estimation unit 102 calculates the position P GPS and P INS and the speed V GPS without calculating the position and speed differences δP ′ = P GPS −P INS and δV ′ = V GPS −V INS. The best position value P and the best speed value V may be estimated directly from V INS . In this case, in the calculation of the Kalman filter for estimation, the measurement vector z to be input and the state vector x to be output are appropriately corrected, and the best value P given by the calculated state vector x is directly used as navigation software (CPU 10 ).
Further, instead of the magnetic sensor 13, another direction sensor such as a gyro sensor may be used.

また、GPS部11がGPS信号を取得できない、もしくはGPS信号の信頼性が十分でない場合は、図5および図6に示す方法によって加速度オフセット値を算出することも可能である。
すなわち、車両の加速や減速は通常ある短い期間(例えば、数秒間)で完結するため、それによる加速度は短い期間で急激に変化するのに対し、道路の傾斜は急激に変わらないため、それによる加速度は、長い期間ほぼ一定値を示すことに着目すると、加速度センサから出力される加速度が、車両の加速や減速によるものか、道路の傾斜等によるものかを判断することができるので、その判断結果に基づいて加速度オフセット値を算出することができる。
If the GPS unit 11 cannot acquire a GPS signal or if the reliability of the GPS signal is not sufficient, the acceleration offset value can be calculated by the method shown in FIGS.
That is, the acceleration and deceleration of the vehicle is usually completed in a short period (for example, a few seconds), and the acceleration caused by the acceleration changes abruptly in a short period, while the slope of the road does not change abruptly. Focusing on the fact that acceleration shows a constant value over a long period of time, it is possible to determine whether the acceleration output from the acceleration sensor is due to acceleration or deceleration of the vehicle, or due to road inclination, etc. An acceleration offset value can be calculated based on the result.

以下、具体的に車両の加速や減速と道路の傾斜を区別する方法を説明する。図5は、車両の走行によって得られる加速度センサの出力例を示したグラフである。図5において、加速度波形a1は車両の加速や減速による加速度センサ12の出力の変化を表しており、この加速度波形a1では、期間t1において加速度がaからbまで大きく変化している。また、加速度波形a2は道路の傾斜による加速度センサ12の出力の変化を表しており、この加速度波形a2では、期間t2において加速度がcからdまでほとんど変化していない。このように、ある期間内の加速度の変化量は、車両の加速や減速の場合と、道路の傾斜の場合とで異なっている。そこで、加速度の最大値と最小値の差(変化量)が予め定められた所定値より大きければ車両の加速や減速と判断し、予め定められた所定値よりも小さければ道路の傾斜と判断することで、両者を区別することができる。   Hereinafter, a method for distinguishing between vehicle acceleration and deceleration and road inclination will be described in detail. FIG. 5 is a graph showing an output example of the acceleration sensor obtained by running the vehicle. In FIG. 5, an acceleration waveform a1 represents a change in the output of the acceleration sensor 12 due to acceleration or deceleration of the vehicle. In the acceleration waveform a1, the acceleration changes greatly from a to b in the period t1. Further, the acceleration waveform a2 represents a change in the output of the acceleration sensor 12 due to the road inclination. In the acceleration waveform a2, the acceleration hardly changes from c to d in the period t2. Thus, the amount of change in acceleration within a certain period differs between the case of acceleration or deceleration of the vehicle and the case of slope of the road. Therefore, if the difference (change amount) between the maximum value and the minimum value of acceleration is larger than a predetermined value, it is determined that the vehicle is accelerating or decelerating. If it is smaller than the predetermined value, it is determined that the road is inclined. Thus, the two can be distinguished.

次に、加速度センサから出力される加速度に従って、上記方法を利用して加速度オフセット値を算出する具体例を説明する。図6は、加速度センサの出力から加速度オフセット値を算出する方法の例を示したグラフである。まず、加速度センサ12から図6(A)に示す加速度波形A1のような加速度が出力された場合、低域通過フィルタを用いることよって車両の振動成分を示す短周期のノイズを抑圧し、図6(B)に示すフィルタ出力波形B1を得る。   Next, a specific example of calculating an acceleration offset value using the above method according to the acceleration output from the acceleration sensor will be described. FIG. 6 is a graph showing an example of a method for calculating the acceleration offset value from the output of the acceleration sensor. First, when acceleration such as the acceleration waveform A1 shown in FIG. 6A is output from the acceleration sensor 12, a short-cycle noise indicating a vehicle vibration component is suppressed by using a low-pass filter, and FIG. A filter output waveform B1 shown in (B) is obtained.

次に、加速度オフセット算出部104は、図6(B)のフィルタ出力波形B1から、予め定められた所定の期間内における平均値と、該期間内の最大値および最小値とを算出する。加速度オフセット算出部104は、算出した最大値と最小値との差が予め定められた所定の値以上である場合は、車両の加速や減速と判断し、算出した最大値と最小値との差が予め定められた所定の値以下である場合は、車両が加速や減速を行っていない、すなわち、道路の傾斜等による加速度であると判断する。図6(B)のフィルタ出力波形B1においては、期間aと期間cの期間が加速または減速をしていると判断され、期間bが道路の傾斜と判断される。   Next, the acceleration offset calculation unit 104 calculates an average value in a predetermined period, and a maximum value and a minimum value in the period, from the filter output waveform B1 in FIG. The acceleration offset calculation unit 104 determines that the vehicle is accelerated or decelerated when the difference between the calculated maximum value and the minimum value is equal to or greater than a predetermined value, and the difference between the calculated maximum value and the minimum value is determined. Is less than or equal to a predetermined value, it is determined that the vehicle is not accelerating or decelerating, that is, acceleration due to road inclination or the like. In the filter output waveform B1 of FIG. 6B, it is determined that the period a and the period c are accelerating or decelerating, and the period b is determined as a road slope.

加速度オフセット算出部104は、加速または減速を行っていないと判断した期間(期間b)は、算出した平均値を加速度センサ12のオフセット値とする。また、加速度オフセット算出部104は、加速または減速をしていると判断した期間(期間aと期間c)のオフセット値は、直前の値を維持して変化させない。つまり、時刻“0”ではまだオフセット値が求まっていないため、期間aではオフセット値=“0”を維持し、期間cでは期間bの最後のオフセット値を維持する。加速度オフセット算出部104は、このようにして算出した加速度オフセット値をINS系位置算出部101へ出力する。この加速度オフセット算出部104が算出した加速度オフセット値をグラフに表したものが、図6(C)のオフセット波形C1である。このように、GPS部11がGPS信号を取得できない、もしくはGPS信号の信頼性が十分でない場合には、GPS信号を用いることなく、加速度オフセット値を算出することが可能となる。   The acceleration offset calculation unit 104 sets the calculated average value as the offset value of the acceleration sensor 12 during the period (period b) in which it is determined that acceleration or deceleration is not performed. Further, the acceleration offset calculation unit 104 maintains the previous value and does not change the offset value of the period (period a and period c) determined as accelerating or decelerating. That is, since the offset value has not yet been obtained at time “0”, offset value = “0” is maintained in period a, and the last offset value in period b is maintained in period c. The acceleration offset calculation unit 104 outputs the acceleration offset value calculated in this way to the INS system position calculation unit 101. A graph representing the acceleration offset value calculated by the acceleration offset calculation unit 104 is an offset waveform C1 in FIG. As described above, when the GPS unit 11 cannot acquire the GPS signal or the reliability of the GPS signal is not sufficient, the acceleration offset value can be calculated without using the GPS signal.

また、車両の走行においては、カーブを走行中の遠心力によっても加速度センサ12から出力される加速度が変化する。このカーブによる加速度を考慮して車両の加速や減速による正しい加速度を求めるため、上述したステップS52によって求められた車両の進行方向θ’に基づき、加速度センサ12から出力される加速度のθ’方向の成分と、θ’に垂直な方向の成分を求め、前者を車両の加速や減速と判断し、後者をカーブにおける遠心力と判断する。そして、進行方向θ’の成分のみを車両の加速や減速による加速度の値として使用することも可能である。このことにより、より精度良く車両の加速度を求めることが可能となる。   Further, when the vehicle is traveling, the acceleration output from the acceleration sensor 12 is also changed by the centrifugal force while traveling along the curve. In order to obtain the correct acceleration due to acceleration or deceleration of the vehicle in consideration of the acceleration due to this curve, the acceleration in the θ ′ direction of the acceleration output from the acceleration sensor 12 is based on the traveling direction θ ′ of the vehicle obtained in step S52 described above. A component and a component in a direction perpendicular to θ ′ are obtained, and the former is determined as acceleration or deceleration of the vehicle, and the latter is determined as a centrifugal force on the curve. It is also possible to use only the component in the traveling direction θ ′ as the acceleration value due to acceleration or deceleration of the vehicle. As a result, the acceleration of the vehicle can be obtained with higher accuracy.

本発明の一実施形態によるナビゲーション装置の構成図である。It is a block diagram of the navigation apparatus by one Embodiment of this invention. ナビゲーション装置の機能ブロックを示す図である。It is a figure which shows the functional block of a navigation apparatus. ナビゲーション装置の動作例を示すフローチャートである。It is a flowchart which shows the operation example of a navigation apparatus. 加速度の絶対値が所定の閾値より小さい場合に位置および速度を算出する処理の例を示すフローチャートである。It is a flowchart which shows the example of a process which calculates a position and speed when the absolute value of an acceleration is smaller than a predetermined threshold value. 車両の走行によって得られる加速度センサの出力例を示したグラフである。It is the graph which showed the output example of the acceleration sensor obtained by driving | running | working of a vehicle. 加速度センサの出力から加速度オフセット値を算出する方法の例を示したグラフである。It is the graph which showed the example of the method of calculating an acceleration offset value from the output of an acceleration sensor.

符号の説明Explanation of symbols

10…CPU 11…GPS部 12…加速度センサ 13…磁気センサ 14…メモリ 15…外部記憶装置 16…通信部 17…表示装置 18…音声出力装置 19…バスライン 101…INS系位置算出部 102…最良値推定部 103…位置補正部 104…加速度オフセット算出部 105…加速度判定部   DESCRIPTION OF SYMBOLS 10 ... CPU 11 ... GPS part 12 ... Acceleration sensor 13 ... Magnetic sensor 14 ... Memory 15 ... External storage device 16 ... Communication part 17 ... Display device 18 ... Audio output device 19 ... Bus line 101 ... INS system position calculation part 102 ... Best Value estimating unit 103 ... Position correcting unit 104 ... Acceleration offset calculating unit 105 ... Acceleration determining unit

Claims (3)

現在位置を求め位置案内を行うナビゲーション装置において、
加速度センサから得られる加速度と所定の閾値とを比較する加速度判定手段と、
前記加速度判定手段により加速度の絶対値が前記閾値より大きいと判定された場合は、該加速度と記憶手段に記憶されている過去の位置および移動速度とに基づいて現在の位置を算出し、前記加速度判定手段により加速度の絶対値が前記閾値より小さいと判定された場合は、等速運動であるとみなして方位センサによる計測に基づく方位と記憶手段に記憶されている過去の位置および移動速度とに基づいて現在の位置を算出する位置算出手段と、
を備えることを特徴とするナビゲーション装置。
In a navigation device that obtains the current position and provides position guidance,
Acceleration determining means for comparing the acceleration obtained from the acceleration sensor with a predetermined threshold;
If the acceleration determining means determines that the absolute value of the acceleration is greater than the threshold, the current position is calculated based on the acceleration and the past position and moving speed stored in the storage means, and the acceleration When the determination means determines that the absolute value of the acceleration is smaller than the threshold value, it is regarded as a constant velocity motion, and the azimuth based on the measurement by the azimuth sensor and the past position and movement speed stored in the storage means A position calculating means for calculating a current position based on;
A navigation device comprising:
GPSにより位置を求めるとともに該求めた異なる時刻における複数の位置から移動速度を求める計測手段と、
前記計測手段により求められた位置および移動速度と、前記位置算出手段により求められた位置および該位置に基づき計算される移動速度と、をカルマンフィルタに入力して現在の最も確からしい位置に関する最良値を推定する最良値推定手段と、
前記最良値推定手段によって推定された最良値に基づき、前記位置算出手段により求められた位置の誤差が前記加速度センサのオフセット誤差に起因したものであるとして該オフセット誤差の分だけ加速度を補正するための加速度オフセット値を算出する加速度オフセット算出手段と、
を更に備え、
前記加速度判定手段は、前記加速度センサから得られる加速度を前記加速度オフセット算出手段により算出された加速度オフセット値を用いて補正し、該補正された加速度と前記閾値とを比較する
ことを特徴とする請求項1に記載のナビゲーション装置。
Measuring means for obtaining a position by GPS and obtaining a moving speed from a plurality of positions at different times obtained;
The position and moving speed obtained by the measuring means and the position obtained by the position calculating means and the moving speed calculated based on the position are input to the Kalman filter to obtain the best value regarding the current most probable position. A best value estimating means for estimating;
Based on the best value estimated by the best value estimating means, the position error obtained by the position calculating means is assumed to be caused by the offset error of the acceleration sensor, and the acceleration is corrected by the offset error. Acceleration offset calculating means for calculating an acceleration offset value of
Further comprising
The acceleration determination means corrects the acceleration obtained from the acceleration sensor using the acceleration offset value calculated by the acceleration offset calculation means, and compares the corrected acceleration with the threshold value. Item 4. The navigation device according to Item 1.
前記加速度センサから得られる加速度の所定時間における変化が予め定められた値以下である場合に、前記加速度センサから得られた加速度の値に応じて加速度オフセット値を算出する加速度オフセット算出手段を更に備え、
前記加速度判定手段は、前記加速度センサから得られる加速度を前記加速度オフセット算出手段により算出された加速度オフセット値を用いて補正し、該補正された加速度と前記閾値とを比較する
ことを特徴とする請求項1に記載のナビゲーション装置。
Acceleration offset calculating means for calculating an acceleration offset value according to the acceleration value obtained from the acceleration sensor when a change in the acceleration obtained from the acceleration sensor at a predetermined time is equal to or less than a predetermined value. ,
The acceleration determination means corrects the acceleration obtained from the acceleration sensor using the acceleration offset value calculated by the acceleration offset calculation means, and compares the corrected acceleration with the threshold value. Item 4. The navigation device according to Item 1.
JP2008218126A 2007-09-25 2008-08-27 Navigator Pending JP2009098127A (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2008218126A JP2009098127A (en) 2007-09-25 2008-08-27 Navigator
US12/235,226 US8108140B2 (en) 2007-09-25 2008-09-22 Navigation device
EP08016699.4A EP2042832B1 (en) 2007-09-25 2008-09-23 Navigation device with GPS unit, acceleration and bearing sensor

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2007247783 2007-09-25
JP2008218126A JP2009098127A (en) 2007-09-25 2008-08-27 Navigator

Publications (1)

Publication Number Publication Date
JP2009098127A true JP2009098127A (en) 2009-05-07

Family

ID=40701259

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008218126A Pending JP2009098127A (en) 2007-09-25 2008-08-27 Navigator

Country Status (1)

Country Link
JP (1) JP2009098127A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013132968A (en) * 2011-12-26 2013-07-08 Denso Corp Vehicle attitude estimating device and theft informing device
JP2014528061A (en) * 2011-09-01 2014-10-23 クアルコム,インコーポレイテッド MEMS accelerometer takeoff and landing detection of passenger aircraft

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08285621A (en) * 1995-04-14 1996-11-01 Omron Corp Navigation device
JPH09189564A (en) * 1996-01-11 1997-07-22 Matsushita Electric Ind Co Ltd Traveling body position speed calculating device
JPH10153620A (en) * 1996-11-25 1998-06-09 Murata Mfg Co Ltd Signal-processing system of acceleration sensor
JP2002333332A (en) * 2001-05-08 2002-11-22 Pioneer Electronic Corp Hybrid processing method and device thereof, navigation system, and computer program
JP2005164590A (en) * 2004-11-25 2005-06-23 Pioneer Electronic Corp Navigation device
JP2005195395A (en) * 2004-01-05 2005-07-21 Mitsubishi Electric Corp Moving object acceleration/distance estimating circuit, pseudo-distance estimating circuit for positioning navigation, moving object positioning device, and moving object positioning method
WO2007015137A1 (en) * 2005-08-01 2007-02-08 Toyota Jidosha Kabushiki Kaisha Zero-point correction apparatus and method for an angular speed sensor
JP2007232443A (en) * 2006-02-28 2007-09-13 Yokogawa Electric Corp Inertia navigation system and its error correction method

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08285621A (en) * 1995-04-14 1996-11-01 Omron Corp Navigation device
JPH09189564A (en) * 1996-01-11 1997-07-22 Matsushita Electric Ind Co Ltd Traveling body position speed calculating device
JPH10153620A (en) * 1996-11-25 1998-06-09 Murata Mfg Co Ltd Signal-processing system of acceleration sensor
JP2002333332A (en) * 2001-05-08 2002-11-22 Pioneer Electronic Corp Hybrid processing method and device thereof, navigation system, and computer program
JP2005195395A (en) * 2004-01-05 2005-07-21 Mitsubishi Electric Corp Moving object acceleration/distance estimating circuit, pseudo-distance estimating circuit for positioning navigation, moving object positioning device, and moving object positioning method
JP2005164590A (en) * 2004-11-25 2005-06-23 Pioneer Electronic Corp Navigation device
WO2007015137A1 (en) * 2005-08-01 2007-02-08 Toyota Jidosha Kabushiki Kaisha Zero-point correction apparatus and method for an angular speed sensor
JP2008542782A (en) * 2005-08-01 2008-11-27 トヨタ自動車株式会社 Apparatus and method for correcting zero point of angular velocity sensor
JP2007232443A (en) * 2006-02-28 2007-09-13 Yokogawa Electric Corp Inertia navigation system and its error correction method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014528061A (en) * 2011-09-01 2014-10-23 クアルコム,インコーポレイテッド MEMS accelerometer takeoff and landing detection of passenger aircraft
JP2013132968A (en) * 2011-12-26 2013-07-08 Denso Corp Vehicle attitude estimating device and theft informing device

Similar Documents

Publication Publication Date Title
JP4466705B2 (en) Navigation device
US8185308B2 (en) Angular velocity correcting device, angular velocity correcting method, and navigation device
US9448250B2 (en) Detecting mount angle of mobile device in vehicle using motion sensors
EP2646777B1 (en) Inertial sensor aided heading and positioning for gnss vehicle navigation
US8108140B2 (en) Navigation device
JP4941199B2 (en) Navigation device
WO2002018873A2 (en) Calibration of multi-axis accelerometer in vehicle navigation system using gps data
JP2007033395A (en) Position correction method and navigation apparatus
US20080243384A1 (en) Azimuth determination apparatus, azimuth determination method and azimuth determination program
JP2004138553A (en) Device, method, program, and record medium for detecting position of moving object
KR100870091B1 (en) Method and apparatus for decide turn condition using sensor
CN109141470A (en) Electronic equipment, error calibration method and recording medium
JPH10253352A (en) Movement judging apparatus and navigator
JP2015004593A (en) Navigation device
JP2009098127A (en) Navigator
JP4376738B2 (en) Apparatus and method for detecting zero point error of angular velocity sensor
JP2008008653A (en) Mounting direction determination method, navigation device and program
US9791277B2 (en) Apparatus and method for measuring velocity of moving object in a navigation system
JP4824522B2 (en) In-vehicle device
JP4955131B1 (en) Wind correction device, wind correction method, wind correction program, and recording medium
JP4967724B2 (en) Location measuring apparatus and method
JP5916540B2 (en) Speed calculation device, speed calculation method, speed calculation program, and recording medium
JP3551661B2 (en) Navigation device
JP2006266842A (en) Method for calibrating acceleration, and navigation apparatus
JP4271541B2 (en) Gradient angle detection apparatus and method

Legal Events

Date Code Title Description
A621 Written request for application examination

Effective date: 20110620

Free format text: JAPANESE INTERMEDIATE CODE: A621

A131 Notification of reasons for refusal

Effective date: 20121002

Free format text: JAPANESE INTERMEDIATE CODE: A131

A977 Report on retrieval

Effective date: 20121003

Free format text: JAPANESE INTERMEDIATE CODE: A971007

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20130226