JP2012154769A - Acceleration detection method, position calculation method and acceleration detection device - Google Patents

Acceleration detection method, position calculation method and acceleration detection device Download PDF

Info

Publication number
JP2012154769A
JP2012154769A JP2011013749A JP2011013749A JP2012154769A JP 2012154769 A JP2012154769 A JP 2012154769A JP 2011013749 A JP2011013749 A JP 2011013749A JP 2011013749 A JP2011013749 A JP 2011013749A JP 2012154769 A JP2012154769 A JP 2012154769A
Authority
JP
Japan
Prior art keywords
acceleration
vector
posture
detected
moving body
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.)
Withdrawn
Application number
JP2011013749A
Other languages
Japanese (ja)
Inventor
Kenji Onda
健至 恩田
Kumar Anand
アナン クマー
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Seiko Epson Corp
Original Assignee
Seiko Epson Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Seiko Epson Corp filed Critical Seiko Epson Corp
Priority to JP2011013749A priority Critical patent/JP2012154769A/en
Publication of JP2012154769A publication Critical patent/JP2012154769A/en
Withdrawn legal-status Critical Current

Links

Images

Abstract

PROBLEM TO BE SOLVED: To provide a new method for more appropriately detecting acceleration of a moving body.SOLUTION: An attitude adjustment section 20 adjusts a detection attitude detected by an attitude sensor 4 installed on a moving body using a detection acceleration vector detected by an acceleration sensor 3 installed on the moving body as well as a given reference vector relevant to a movement of the moving body. A coordinate conversion section 30 obtains a conversion coefficient related to coordinate conversion from a local coordinate system of the moving body to an absolute coordinate system using a detection attitude adjusted by the attitude adjustment section 20 and converts the detection acceleration vector detected by the acceleration sensor 3 in the local coordinate system into the acceleration vector in the absolute coordinate system.

Description

本発明は、加速度検出方法、位置算出方法及び加速度検出装置に関する。   The present invention relates to an acceleration detection method, a position calculation method, and an acceleration detection device.

いわゆるシームレス測位やモーションセンシング、姿勢制御など様々な分野において、慣性センサーの活用が注目されている。慣性センサーとしては、加速度センサーやジャイロセンサー、圧力センサー、地磁気センサーなどが広く知られている。慣性センサーの検出結果を利用して慣性航法演算を行って位置算出を行う技術も考案されている。   In various fields such as so-called seamless positioning, motion sensing, and attitude control, the use of inertial sensors is attracting attention. As an inertial sensor, an acceleration sensor, a gyro sensor, a pressure sensor, a geomagnetic sensor, and the like are widely known. There has also been devised a technique for calculating the position by performing inertial navigation calculation using the detection result of the inertial sensor.

慣性航法演算においては、慣性センサーの検出結果に含まれ得る種々の誤差成分に起因して位置算出の正確性が低下するという問題があり、位置算出の正確性を向上させるための様々な技術が考案されている。例えば、特許文献1には、慣性航法演算を行うことで測定したデータを、GPS(Global Positioning System)を利用して補正する技術が開示されている。   In inertial navigation calculation, there is a problem that the accuracy of position calculation is reduced due to various error components that can be included in the detection result of the inertial sensor, and there are various techniques for improving the accuracy of position calculation. It has been devised. For example, Patent Document 1 discloses a technique for correcting data measured by performing an inertial navigation calculation using a GPS (Global Positioning System).

特開2004−239643号公報JP 2004-239634 A

特許文献1に開示されている技術では、慣性航法演算を行うことで測定した移動体の位置や速度、加速度、姿勢角、方位角といった慣性航法演算に係る諸量それぞれについて、カルマンフィルターを利用して、各諸量に含まれ得る誤差を推定している。この場合、慣性航法演算に係る全ての諸量についてカルマンフィルターを利用して誤差推定演算を行う必要があるため、演算量は膨大になりがちである。   In the technique disclosed in Patent Document 1, a Kalman filter is used for each of various quantities related to the inertial navigation calculation such as the position, velocity, acceleration, posture angle, and azimuth of the moving body measured by performing the inertial navigation calculation. Thus, the error that can be included in each quantity is estimated. In this case, since it is necessary to perform error estimation calculation using the Kalman filter for all the quantities related to the inertial navigation calculation, the calculation amount tends to be enormous.

慣性航法演算では、例えば移動体に設置された加速度センサーを用いるのであれば、加速度センサーが検出した検出加速度を積分することで移動体の速度を算出し、さらに、その速度を積分することで位置を算出する。   In inertial navigation calculation, for example, if an acceleration sensor installed on a moving body is used, the speed of the moving body is calculated by integrating the detected acceleration detected by the acceleration sensor, and the position is calculated by integrating the speed. Is calculated.

この場合に算出される速度や位置に含まれ得る誤差は、そもそも加速度センサーの検出加速度に含まれ得る誤差に起因するものである。つまり、検出加速度に対して積分を繰り返すために、検出加速度に含まれていた誤差が、時間経過に伴って、速度や位置に誤差として重畳・累積されていくことによるものである。そのため、移動体の加速度を正しく検出することができれば、特許文献1のように慣性航法演算に係る全ての諸量の誤差を推定せずとも、慣性航法演算の精度を向上させることが可能であると考えられる。   In this case, the error that can be included in the calculated speed and position originates from an error that can be included in the acceleration detected by the acceleration sensor. That is, in order to repeat the integration with respect to the detected acceleration, the error included in the detected acceleration is superimposed and accumulated as an error in the speed and position as time elapses. Therefore, if the acceleration of the moving body can be detected correctly, it is possible to improve the accuracy of the inertial navigation calculation without estimating all the errors related to the inertial navigation calculation as in Patent Document 1. it is conceivable that.

本発明は上述した課題に鑑みて為されたものであり、その目的とするところは、移動体の加速度をより適切に検出するための新しい手法を提案することにある。   The present invention has been made in view of the above-described problems, and an object of the present invention is to propose a new method for more appropriately detecting the acceleration of a moving body.

以上の課題を解決するための第1の形態は、移動体に設置された姿勢センサーの検出姿勢を、前記移動体に設置された加速度センサーの検出加速度ベクトルと前記移動体の移動に係る所与の参照ベクトルとを用いて調整することと、前記調整された前記検出姿勢を用いて、前記加速度センサーに対応付けられたローカル座標系から絶対座標系への座標変換に係る変換係数を求めることと、前記変換係数を用いて、前記検出加速度ベクトルを前記絶対座標系における加速度ベクトルに座標変換することと、を含む加速度検出方法である。   In a first embodiment for solving the above-described problems, a detected attitude of an attitude sensor installed on a moving body is determined based on a detected acceleration vector of an acceleration sensor installed on the moving body and the movement of the moving body. Adjusting using the reference vector, and using the adjusted detection posture, obtaining a conversion coefficient related to coordinate conversion from a local coordinate system associated with the acceleration sensor to an absolute coordinate system; And converting the detected acceleration vector into an acceleration vector in the absolute coordinate system using the conversion coefficient.

また、他の形態として、移動体に設置された姿勢センサーの検出姿勢を、前記移動体に設置された加速度センサーの検出加速度ベクトルと前記移動体の移動に係る所与の参照ベクトルとを用いて調整する調整部と、前記調整された前記検出姿勢を用いて、前記加速度センサーに対応付けられたローカル座標系から絶対座標系への座標変換に係る変換係数を求める変換係数算出部と、前記変換係数を用いて、前記検出加速度ベクトルを前記絶対座標系における加速度ベクトルに座標変換する変換部と、を備えた加速度検出装置を構成してもよい。   As another form, the detected posture of the posture sensor installed on the moving body is detected using the detected acceleration vector of the acceleration sensor installed on the moving body and a given reference vector related to the movement of the moving body. An adjustment unit for adjusting, a conversion coefficient calculation unit for obtaining a conversion coefficient related to coordinate conversion from a local coordinate system to an absolute coordinate system associated with the acceleration sensor, using the adjusted detection posture; and the conversion You may comprise the acceleration detection apparatus provided with the conversion part which carries out coordinate conversion of the said detected acceleration vector into the acceleration vector in the said absolute coordinate system using a coefficient.

この第1の形態等によれば、移動体に設置された姿勢センサーの検出姿勢を、移動体に設置された加速度センサーの検出加速度ベクトルと移動体の移動に係る所与の参照ベクトルとを用いて調整する。そして、調整された検出姿勢を用いて、加速度センサーに対応付けられたローカル座標系から絶対座標系への座標変換に係る変換係数を求める。そして、求めた変換係数を用いて、加速度センサーの検出加速度ベクトルを絶対座標系における加速度ベクトルに座標変換する。   According to the first aspect and the like, the detected attitude of the attitude sensor installed on the moving body is determined using the detected acceleration vector of the acceleration sensor installed on the moving body and the given reference vector related to the movement of the moving body. Adjust. Then, using the adjusted detection posture, a conversion coefficient related to coordinate conversion from the local coordinate system associated with the acceleration sensor to the absolute coordinate system is obtained. Then, using the obtained conversion coefficient, the detected acceleration vector of the acceleration sensor is coordinate-converted to an acceleration vector in the absolute coordinate system.

加速度センサーが検出する加速度ベクトルは、加速度センサーに対応付けられたローカルな座標系における加速度ベクトルである。それに対して、例えば慣性航法演算においては、移動体の移動空間を定める絶対的な座標系を基準として演算を行うことが一般的である。姿勢センサーの検出姿勢には誤差が含まれ得るため、姿勢センサーの検出姿勢を、加速度センサーの検出加速度ベクトルと移動体の移動に係る所与の参照ベクトルとを用いて調整する。その上で、ローカル座標系から絶対座標系への変換係数を求めて座標変換を行うことで、絶対座標系における加速度ベクトルを適切に検出することができる。このようにして検出された絶対座標系における加速度ベクトルは、上記のように慣性航法演算等に好適である。   The acceleration vector detected by the acceleration sensor is an acceleration vector in a local coordinate system associated with the acceleration sensor. On the other hand, for example, in inertial navigation calculation, the calculation is generally performed based on an absolute coordinate system that defines the moving space of the moving body. Since an error may be included in the detected posture of the posture sensor, the detected posture of the posture sensor is adjusted using the detected acceleration vector of the acceleration sensor and a given reference vector related to the movement of the moving body. Then, by obtaining a conversion coefficient from the local coordinate system to the absolute coordinate system and performing coordinate conversion, the acceleration vector in the absolute coordinate system can be detected appropriately. The acceleration vector in the absolute coordinate system detected in this manner is suitable for inertial navigation calculation and the like as described above.

また、第2の形態として、第1の形態の加速度検出方法であって、前記変換係数を求めることは、前記検出加速度ベクトルに含まれ得る誤差を小さくすることが可能な前記変換係数を求めることである、加速度検出方法を構成してもよい。   Further, as a second form, the acceleration detection method according to the first form, wherein the conversion coefficient is obtained by obtaining the conversion coefficient capable of reducing an error that can be included in the detected acceleration vector. An acceleration detection method may be configured.

この第2の形態によれば、検出加速度ベクトルに含まれ得る誤差を小さくすることが可能な変換係数を求める。そのため、当該変換係数を用いて加速度センサーの検出加速度ベクトルを絶対座標系における加速度ベクトルに変換することで、加速度センサーの検出加速度ベクトルに含まれ得る誤差を低減させることができる。   According to the second embodiment, a conversion coefficient capable of reducing an error that can be included in the detected acceleration vector is obtained. Therefore, by converting the acceleration vector detected by the acceleration sensor into an acceleration vector in the absolute coordinate system using the conversion coefficient, errors that can be included in the acceleration vector detected by the acceleration sensor can be reduced.

また、第3の形態として、第1又は第2の形態の加速度検出方法であって、前記調整することと、前記変換係数を求めることと、前記座標変換することと、を前記移動体の移動中に行う、加速度検出方法を構成してもよい。   Further, as a third form, the acceleration detection method according to the first or second form, wherein the adjustment, the conversion coefficient is obtained, and the coordinate conversion is performed, and the moving body is moved. You may comprise the acceleration detection method performed inside.

この第3の形態によれば、姿勢センサーの検出姿勢の調整と、変換係数を求めることと、検出加速度ベクトルの座標変換と、を移動体の移動中に行うことで、移動中に移動体の姿勢が変化した場合であっても、移動体の加速度ベクトルを適切に検出することができる。   According to the third embodiment, the adjustment of the detected posture of the posture sensor, the conversion coefficient, and the coordinate conversion of the detected acceleration vector are performed while the moving body is moving, so that the moving body is moved during the movement. Even when the posture changes, the acceleration vector of the moving body can be detected appropriately.

また、第4の形態として、第1〜第3の何れかの形態の加速度検出方法であって、前記調整することは、前記検出姿勢のうちのヨー成分については調整を行わず、ロール成分及びピッチ成分の少なくとも一方を調整することを含む、加速度検出方法を構成してもよい。   Further, as a fourth mode, in the acceleration detection method according to any one of the first to third modes, the adjustment is performed without adjusting the yaw component of the detected posture, and the roll component and You may comprise the acceleration detection method including adjusting at least one of a pitch component.

この第4の形態によれば、検出姿勢のうちのヨー成分については調整を行わず、ロール成分及びピッチ成分の少なくとも一方を調整することで、検出姿勢の調整を簡易に行うことができる。   According to the fourth aspect, the yaw component in the detected posture is not adjusted, and the adjustment of the detected posture can be easily performed by adjusting at least one of the roll component and the pitch component.

また、第5の形態として、第1〜第4の何れかの形態の加速度検出方法であって、前記参照ベクトルは、前記移動体に設置された所定の速度計測システムによって計測された前記移動体の速度ベクトルを含み、前記調整することは、前記速度ベクトルを用いて、前記移動体の加速度ベクトルを算出することと、前記検出加速度ベクトルと前記算出した加速度ベクトルとを用いて、前記検出姿勢を調整することと、を含む、加速度検出方法を構成してもよい。   Further, as a fifth mode, the acceleration detecting method according to any one of the first to fourth modes, wherein the reference vector is measured by a predetermined speed measurement system installed in the mobile unit. The adjusting includes calculating the acceleration vector of the moving body using the velocity vector, and using the detected acceleration vector and the calculated acceleration vector to adjust the detected posture. And adjusting the acceleration.

この第5の形態によれば、参照ベクトルには、移動体に設置された所定の速度計測システムによって計測された移動体の速度ベクトルが含まれる。この場合において、速度計測システムによって計測された移動体の速度ベクトルを用いて、移動体の加速度ベクトルを算出する。そして、加速度センサーの検出加速度ベクトルと算出した加速度ベクトルとを用いて、姿勢センサーの検出姿勢を調整する。速度計測システムによって計測された移動体の速度ベクトルを利用することで、姿勢センサーの検出姿勢の調整を適切に行うことができる。   According to the fifth aspect, the reference vector includes the velocity vector of the moving object measured by the predetermined velocity measuring system installed on the moving object. In this case, the acceleration vector of the moving object is calculated using the velocity vector of the moving object measured by the speed measurement system. Then, the detected posture of the posture sensor is adjusted using the detected acceleration vector of the acceleration sensor and the calculated acceleration vector. By using the velocity vector of the moving body measured by the velocity measurement system, the detection posture of the posture sensor can be adjusted appropriately.

より具体的には、第6の形態のように、第5の形態における速度計測システムが、車速検出システム、又は、衛星測位システムである、加速度検出方法を構成することも可能である。   More specifically, as in the sixth embodiment, it is possible to configure an acceleration detection method in which the speed measurement system in the fifth embodiment is a vehicle speed detection system or a satellite positioning system.

この第6の形態によれば、車速検出システム、又は、衛星測位システムを利用することで、移動体の速度ベクトルを精度良く計測することが可能となり、移動体の加速度ベクトルの算出に供することができる。   According to the sixth embodiment, by using the vehicle speed detection system or the satellite positioning system, it becomes possible to measure the velocity vector of the moving body with high accuracy, and it is possible to use for the calculation of the acceleration vector of the moving body. it can.

また、第7の形態として、第1〜第6の何れかの形態の加速度検出方法を用いて前記絶対座標系における加速度ベクトルを検出することと、所与の移動基準位置と前記検出された加速度ベクトルとを用いて、前記絶対座標系における前記移動体の位置を算出することと、を含む位置算出方法を構成してもよい。   Further, as a seventh aspect, the acceleration detection method of any one of the first to sixth aspects is used to detect an acceleration vector in the absolute coordinate system, and a given movement reference position and the detected acceleration. You may comprise the position calculation method including calculating the position of the said mobile body in the said absolute coordinate system using a vector.

この第7の形態によれば、所与の移動基準位置と上記の形態の加速度検出方法を用いて検出した加速度ベクトルとを用いることで、移動体の位置を高い正確性で算出することができる。   According to the seventh aspect, the position of the moving object can be calculated with high accuracy by using the given movement reference position and the acceleration vector detected using the acceleration detection method of the above aspect. .

加速度検出システムのシステム構成の説明図。Explanatory drawing of the system configuration | structure of an acceleration detection system. 第1の加速度検出システムのシステム構成の説明図。Explanatory drawing of the system configuration | structure of a 1st acceleration detection system. 第2の加速度検出システムのシステム構成の説明図。Explanatory drawing of the system configuration | structure of a 2nd acceleration detection system. 第3の加速度検出システムのシステム構成の説明図。Explanatory drawing of the system configuration | structure of a 3rd acceleration detection system. 絶対座標速度ベクトル補正の原理の説明図。Explanatory drawing of the principle of absolute coordinate velocity vector correction | amendment. ナビゲーションシステムのシステム構成の説明図。Explanatory drawing of the system configuration | structure of a navigation system. カーナビゲーション装置の機能構成の説明図。Explanatory drawing of a function structure of a car navigation apparatus. 第1の位置算出システムのシステム構成の説明図。Explanatory drawing of the system configuration | structure of a 1st position calculation system. ナビゲーション処理の流れを示すフローチャート。The flowchart which shows the flow of a navigation process. 第1の位置算出システムによる位置算出結果の説明図。Explanatory drawing of the position calculation result by a 1st position calculation system. 第1の位置算出システムによる姿勢算出結果の説明図。Explanatory drawing of the attitude | position calculation result by a 1st position calculation system. 第2の位置算出システムのシステム構成の説明図。Explanatory drawing of the system configuration | structure of a 2nd position calculation system. 第2の位置算出システムによる位置算出結果の説明図。Explanatory drawing of the position calculation result by a 2nd position calculation system. 第2の位置算出システムによる姿勢算出結果の説明図。Explanatory drawing of the attitude | position calculation result by a 2nd position calculation system. 第3の位置算出システムのシステム構成の説明図。Explanatory drawing of the system configuration | structure of a 3rd position calculation system.

以下、図面を参照して、本発明の好適な実施形態の一例について説明する。但し、本発明を適用可能な実施形態が以下説明する実施形態に限定されるわけでないことは勿論である。   Hereinafter, an example of a preferred embodiment of the present invention will be described with reference to the drawings. However, it goes without saying that embodiments to which the present invention can be applied are not limited to the embodiments described below.

1.原理
1−1.システム構成
図1は、本実施形態における加速度検出方法を実現するための主要な機能ブロックで構成される加速度検出システム1のシステム構成図である。加速度検出システム1は、各種の移動体に適用されるシステムであり、移動体の加速度(詳細には加速度ベクトル)を検出して出力するシステムである。
1. Principle 1-1. System Configuration FIG. 1 is a system configuration diagram of an acceleration detection system 1 configured by main functional blocks for realizing the acceleration detection method according to the present embodiment. The acceleration detection system 1 is a system that is applied to various types of moving objects, and is a system that detects and outputs accelerations (specifically, acceleration vectors) of the moving objects.

加速度検出システム1は、慣性センサーの一種である加速度センサー3及び姿勢センサー4を備えて構成される。図1では、慣性センサーを二重線で図示し、慣性センサーの検出結果を利用して各種処理を行う処理ブロックを一重線で図示することで、両者を区別している。以降の図面についても同様である。   The acceleration detection system 1 includes an acceleration sensor 3 and an attitude sensor 4 that are a kind of inertial sensors. In FIG. 1, the inertial sensor is illustrated by a double line, and a processing block for performing various processes using the detection result of the inertial sensor is illustrated by a single line to distinguish the two. The same applies to the subsequent drawings.

加速度センサー3は、移動体に設置されて移動体の加速度を検出するセンサーであり、例えばMEMS(Micro Electro Mechanical Systems)の技術を利用したMEMSセンサーが用いられる。本実施形態において、加速度センサー3は、移動体を基準として3軸のローカル座標系における3軸のローカル座標加速度ベクトルを検出して出力するように構成されている。つまり、加速度センサー3に対応付けられたローカル座標系が移動体を基準とする座標系と一致するような姿勢で、加速度センサー3が移動体に設置されているものとして説明する。   The acceleration sensor 3 is a sensor that is installed on a moving body and detects the acceleration of the moving body. For example, a MEMS sensor using a technology of MEMS (Micro Electro Mechanical Systems) is used. In the present embodiment, the acceleration sensor 3 is configured to detect and output a triaxial local coordinate acceleration vector in a triaxial local coordinate system with the moving body as a reference. That is, it is assumed that the acceleration sensor 3 is installed on the moving body in such a posture that the local coordinate system associated with the acceleration sensor 3 matches the coordinate system based on the moving body.

姿勢センサー4は、移動体の姿勢を検出するセンサーであり、例えばジャイロセンサーのみで構成してもよいし、磁気センサーや加速度センサーをさらに組み合わせて構成してもよい。姿勢センサー4は移動体に固定され、直交3軸でなるローカル座標系における移動体の姿勢を検出する。すなわち、姿勢センサー4は移動体の姿勢角を計測するための計測器であり、姿勢センサー4が計測した値は、絶対座標系で表した移動体の向きとなる。   The attitude sensor 4 is a sensor that detects the attitude of the moving body. For example, the attitude sensor 4 may be configured by only a gyro sensor, or may be configured by further combining a magnetic sensor and an acceleration sensor. The posture sensor 4 is fixed to the moving body and detects the posture of the moving body in a local coordinate system having three orthogonal axes. That is, the posture sensor 4 is a measuring instrument for measuring the posture angle of the moving body, and the value measured by the posture sensor 4 is the direction of the moving body expressed in an absolute coordinate system.

移動体の加速度及び速度は、方向及び大きさを持っている。そのため、本明細書では、スカラーとベクトルとを適宜区別して説明する。原則的に、単に加速度や速度と言ったときは加速度や速度の大きさ(スカラー量)を表すものとし、加速度ベクトルや速度ベクトルと言ったときは方向及び大きさを考慮した加速度及び速度を表すものとする。また、ローカル座標系において定義されるベクトルと、絶対座標系において定義されるベクトルとを区別するため、各変量を表す文言の先頭に座標系の種類を付して説明する。   The acceleration and speed of the moving body have a direction and a magnitude. For this reason, in this specification, a scalar and a vector are appropriately distinguished and described. In principle, simply referring to acceleration and speed means the magnitude of the acceleration and speed (scalar amount), and referring to acceleration vector and speed vector means the acceleration and speed considering the direction and magnitude. Shall. Further, in order to distinguish between a vector defined in the local coordinate system and a vector defined in the absolute coordinate system, a description will be given with the type of the coordinate system added to the head of a word representing each variable.

座標変換の際には、加速度センサー3の検出加速度ベクトルのうち、加速度ベクトルの向きを座標変換し、加速度ベクトルの大きさは維持する。言い換えると、加速度センサー3の相対的な検出方向を、絶対座標系における方向(絶対方向)に座標変換する。   In the coordinate conversion, the direction of the acceleration vector among the detected acceleration vectors of the acceleration sensor 3 is coordinate-converted, and the magnitude of the acceleration vector is maintained. In other words, the relative detection direction of the acceleration sensor 3 is coordinate-converted into a direction (absolute direction) in the absolute coordinate system.

また、本明細書では、姿勢センサー4から出力される姿勢角を「検出姿勢角」と称し、この検出姿勢角によって定まる移動体の姿勢を「検出姿勢」と称する。また、検出姿勢角に対して姿勢調整部20が調整を行うことで取得した姿勢角を「調整姿勢角」と称し、この調整姿勢角によって定まる移動体の姿勢を「調整姿勢」と称する。   Further, in this specification, the posture angle output from the posture sensor 4 is referred to as “detected posture angle”, and the posture of the moving body determined by the detected posture angle is referred to as “detected posture”. In addition, the posture angle acquired by adjusting the detected posture angle by the posture adjustment unit 20 is referred to as an “adjustment posture angle”, and the posture of the moving body determined by the adjustment posture angle is referred to as an “adjustment posture”.

加速度検出システム1は、主要な処理ブロックとして、姿勢調整部20と、座標変換部30とを有する。   The acceleration detection system 1 includes an attitude adjustment unit 20 and a coordinate conversion unit 30 as main processing blocks.

姿勢調整部20は、姿勢センサー4の検出姿勢を、加速度センサー3の検出加速度ベクトルであるローカル座標加速度ベクトルと、移動体の移動に係る所与の参照ベクトルとを用いて調整する調整部である。   The posture adjusting unit 20 is an adjusting unit that adjusts the detected posture of the posture sensor 4 using a local coordinate acceleration vector that is a detected acceleration vector of the acceleration sensor 3 and a given reference vector related to the movement of the moving body. .

座標変換部30は、姿勢調整部20により調整された姿勢を用いて、加速度センサー3の検出加速度ベクトルであるローカル座標加速度ベクトルを絶対座標加速度ベクトルに座標変換する変換部である。具体的には、姿勢調整部20により調整された検出姿勢である調整姿勢を用いて、移動体のローカル座標系から絶対座標系への座標変換に係る変換係数(詳細には座標変換行列)を求める。そして、求めた変換係数を用いて、ローカル座標加速度ベクトルを絶対座標加速度ベクトルに座標変換する。詳細後述するように、座標変換行列は、座標変換に係る変換係数を成分とする行列である。   The coordinate conversion unit 30 is a conversion unit that converts a local coordinate acceleration vector, which is a detected acceleration vector of the acceleration sensor 3, into an absolute coordinate acceleration vector using the posture adjusted by the posture adjustment unit 20. Specifically, using the adjustment posture that is the detection posture adjusted by the posture adjustment unit 20, a conversion coefficient (specifically, a coordinate conversion matrix) for coordinate conversion from the local coordinate system to the absolute coordinate system of the moving object is obtained. Ask. Then, using the obtained conversion coefficient, the local coordinate acceleration vector is coordinate-converted to an absolute coordinate acceleration vector. As will be described in detail later, the coordinate transformation matrix is a matrix having a transformation coefficient related to coordinate transformation as a component.

座標変換部30は、座標変換に係る変換係数(座標変換行列)を求める変換係数算出部である。座標変換により絶対座標系における加速度ベクトルが算出されるため、座標変換部30は、絶対座標加速度ベクトル算出部であるとも言える。   The coordinate conversion unit 30 is a conversion coefficient calculation unit that calculates a conversion coefficient (coordinate conversion matrix) related to coordinate conversion. Since the acceleration vector in the absolute coordinate system is calculated by the coordinate conversion, it can be said that the coordinate conversion unit 30 is an absolute coordinate acceleration vector calculation unit.

1−2.加速度検出の原理
次に、本実施形態における加速度検出の原理について説明する。最初に、姿勢センサー4の検出姿勢を、加速度センサー3の検出加速度ベクトルであるローカル座標加速度ベクトルと移動体の移動に係る所与の参照ベクトルとを用いて調整することが、本実施形態の特徴の1つである。次いで、調整された検出姿勢を用いて、加速度センサー3に対応付けられたローカル座標系から絶対座標系への座標変換に係る変換係数(座標変換行列)を求める。そして、求めた変換係数を用いて、加速度センサー3の検出加速度ベクトルを絶対座標系における加速度ベクトルに座標変換することで、移動体の加速度を検出する。
1-2. Next, the principle of acceleration detection in this embodiment will be described. First, the detection posture of the posture sensor 4 is adjusted using a local coordinate acceleration vector that is a detection acceleration vector of the acceleration sensor 3 and a given reference vector related to the movement of the moving body. It is one of. Next, a conversion coefficient (coordinate conversion matrix) related to coordinate conversion from the local coordinate system associated with the acceleration sensor 3 to the absolute coordinate system is obtained using the adjusted detection posture. Then, using the obtained conversion coefficient, the acceleration of the moving body is detected by converting the detected acceleration vector of the acceleration sensor 3 into the acceleration vector in the absolute coordinate system.

移動体の移動に係る参照ベクトルとしては、例えば、移動体の加速度ベクトルや速度ベクトルを適用することができる。参照ベクトルは、種々の手法によって取得することが可能である。例えば、移動体に設置された所定の速度計測システムから取得することができる。速度計測システムは、車速検出システム、又は、衛星測位システムである。   As the reference vector related to the movement of the moving body, for example, an acceleration vector or a velocity vector of the moving body can be applied. The reference vector can be obtained by various methods. For example, it can be acquired from a predetermined speed measurement system installed on the moving body. The speed measurement system is a vehicle speed detection system or a satellite positioning system.

車速検出システムでは、例えば車速パルスによって移動体の速度の大きさを検出することができる。磁気センサー等の方位センサーを組み合わせて車速検出システムを構成した場合には、移動体の移動方向も併せて検出することができる。この場合、車速検出システムは、移動体の速度の大きさ及び移動方向を検出する。このため、車速検出システムから、参照ベクトルとして移動体の速度ベクトルを取得することができる。   In the vehicle speed detection system, for example, the magnitude of the speed of the moving body can be detected by a vehicle speed pulse. When a vehicle speed detection system is configured by combining a direction sensor such as a magnetic sensor, the moving direction of the moving body can also be detected. In this case, the vehicle speed detection system detects the magnitude and moving direction of the moving body. For this reason, the speed vector of the moving body can be acquired as the reference vector from the vehicle speed detection system.

衛星測位システムは、代表的にはGPS(Global Positioning System)である。GPSでは、複数のGPS衛星の位置や各GPS衛星からGPS受信機までの擬似距離等の情報に基づいて、GPS受信機の位置座標及び時計誤差を求める位置計算を行う。なお、GPSを利用した位置計算は従来公知であるため、ここでは詳細な説明を省略する。   The satellite positioning system is typically a GPS (Global Positioning System). In GPS, based on information such as the positions of a plurality of GPS satellites and pseudo distances from the respective GPS satellites to the GPS receiver, position calculation for obtaining the GPS receiver position coordinates and clock error is performed. In addition, since the position calculation using GPS is conventionally well-known, detailed description is abbreviate | omitted here.

位置計算の他に、GPSでは、GPS衛星及びGPS受信機の相対的な位置関係の変化により生ずるドップラーを利用して、GPS受信機の速度ベクトルを算出することが可能である。すなわち、GPS衛星から受信したGPS衛星信号のドップラー周波数を測定することで、GPS受信機とGPS衛星との間の相対速度ベクトルを算出する。GPS衛星の速度ベクトル(移動速度及び移動方向)は、衛星軌道データ(エフェメリスやアルマナック)から求めることができる。そのため、相対速度ベクトルとGPS衛星の速度ベクトルとを用いれば、移動体の速度ベクトル(移動速度及び移動方向)を算出することができる。このため、GPSから、参照ベクトルとして移動体の速度ベクトルを取得することができる。   In addition to the position calculation, the GPS can calculate the velocity vector of the GPS receiver by using Doppler generated by the change in the relative positional relationship between the GPS satellite and the GPS receiver. That is, the relative velocity vector between the GPS receiver and the GPS satellite is calculated by measuring the Doppler frequency of the GPS satellite signal received from the GPS satellite. The velocity vector (movement speed and direction) of the GPS satellite can be obtained from satellite orbit data (ephemeris and almanac). Therefore, if the relative velocity vector and the velocity vector of the GPS satellite are used, the velocity vector (movement velocity and movement direction) of the moving body can be calculated. For this reason, the velocity vector of the moving object can be acquired from the GPS as the reference vector.

移動体の速度ベクトルを時間微分すれば、移動体の加速度ベクトルが求まる。そのため、上記の速度計測システムから参照ベクトルとして取得した速度ベクトルを時間微分することで、移動体の加速度ベクトルを算出することができる。   If the velocity vector of the moving body is differentiated with respect to time, the acceleration vector of the moving body can be obtained. Therefore, the acceleration vector of the moving body can be calculated by differentiating the velocity vector acquired as a reference vector from the velocity measurement system with respect to time.

そこで、本実施形態の姿勢調整方法では、上記の参照ベクトルを用いて、姿勢調整に用いる移動体の加速度ベクトル(以下、「姿勢調整用加速度ベクトル」と称す。)を算出する。そして、加速度センサー3の検出加速度ベクトルであるローカル座標加速度ベクトルと算出した姿勢調整用加速度ベクトルとを用いて、検出姿勢角に含まれる誤差(以下、「姿勢角誤差」と称す。)を推定する。この場合、姿勢調整用加速度ベクトルとして、絶対座標系における加速度ベクトルを算出すると好適である。最終的に、推定した姿勢角誤差を用いて検出姿勢角を調整する。   Therefore, in the posture adjustment method of the present embodiment, the acceleration vector of the moving body used for posture adjustment (hereinafter referred to as “posture adjustment acceleration vector”) is calculated using the reference vector. Then, an error included in the detected posture angle (hereinafter referred to as “posture angle error”) is estimated using the local coordinate acceleration vector that is the detected acceleration vector of the acceleration sensor 3 and the calculated posture adjustment acceleration vector. . In this case, it is preferable to calculate an acceleration vector in the absolute coordinate system as the posture adjustment acceleration vector. Finally, the detected posture angle is adjusted using the estimated posture angle error.

説明を分かり易くするために、座標系を具体的に例示して説明する。ここでは、ローカル座標系として機体座標系(Body Frame)(以下、「Bフレーム」と称する。)を適用し、絶対座標系として北東下座標系として知られるNED(North East Down)座標系(以下、「Nフレーム」と称する。)を適用した場合を一例として説明する。   In order to make the description easy to understand, a coordinate system will be specifically exemplified and described. Here, an aircraft coordinate system (Body Frame) (hereinafter referred to as “B frame”) is applied as a local coordinate system, and an NED (North East Down) coordinate system (hereinafter referred to as a northeast lower coordinate system) is known as an absolute coordinate system. , Referred to as “N frame”) will be described as an example.

Bフレームは、例えば、移動体の前方を正とする前後方向をX軸(ロール軸)、右方を正とする左右方向をY軸(ピッチ軸)、鉛直下方を正とする上下方向をZ軸(ヨー軸)とする三次元直交座標系である。それに対して、Nフレームは、例えば、北方を正とする南北方向をN軸、東方を正とする東西方向をE軸、鉛直下方を正とする高度方向をD軸とする三次元直交座標系である。   The B frame is, for example, an X-axis (roll axis) with the front / rear direction positive with respect to the front of the moving body, a Y-axis (pitch axis) with the right side positive, and a vertical direction with the vertical downward direction positive with Z It is a three-dimensional orthogonal coordinate system with an axis (yaw axis). On the other hand, the N frame is, for example, a three-dimensional orthogonal coordinate system in which the north-south direction with the north as positive is the N axis, the east-west direction with the east as positive is the E axis, and the altitude direction with the vertical downward is the D axis. It is.

Bフレーム加速度ベクトルからNフレーム加速度ベクトルへの変換は、次式(1)で定式化される。

Figure 2012154769
The conversion from the B frame acceleration vector to the N frame acceleration vector is formulated by the following equation (1).
Figure 2012154769

式(1)における“fb”及び“fn”は、それぞれBフレーム加速度ベクトル及びNフレーム加速度ベクトルであり、次式(2)及び(3)で与えられる。

Figure 2012154769
Figure 2012154769
“F b ” and “f n ” in equation (1) are a B frame acceleration vector and an N frame acceleration vector, respectively, and are given by the following equations (2) and (3).
Figure 2012154769
Figure 2012154769

また、式(1)における“Cb n”はBフレームからNフレームへの座標変換行列であり、次式(4)で与えられる。

Figure 2012154769
Further, “C b n ” in the equation (1) is a coordinate conversion matrix from the B frame to the N frame, and is given by the following equation (4).
Figure 2012154769

式(4)において、(φ,θ,Ψ)は、それぞれBフレームにおけるX軸(ロール軸)回りの回転角(ロール角)、Y軸(ピッチ軸)回りの回転角(ピッチ角)、Z軸(ヨー軸)回りの回転角(ヨー角)である。また、座標変換行列“Cb n”の各成分は、BフレームからNフレームへの座標変換に係る変換係数である。 In Expression (4), (φ, θ, Ψ) respectively represents a rotation angle (roll angle) around the X axis (roll axis), a rotation angle (pitch angle) around the Y axis (pitch axis), and Z The rotation angle (yaw angle) about the axis (yaw axis). Each component of the coordinate transformation matrix “C b n ” is a transformation coefficient related to coordinate transformation from the B frame to the N frame.

式(1)〜(4)より、次式(5)が導かれる。

Figure 2012154769
From the formulas (1) to (4), the following formula (5) is derived.
Figure 2012154769

ここで、式(5)をロール角“φ”及びピッチ角“θ”の関数とみなし、次式(6)のように定式化する。

Figure 2012154769
Here, the expression (5) is regarded as a function of the roll angle “φ” and the pitch angle “θ”, and is formulated as the following expression (6).
Figure 2012154769

式(6)において、ハット「^」付きの姿勢角は、それぞれロール角“φ”及びピッチ角“θ”の実測値を示し、ハット「^」無しの姿勢角は、それぞれロール角“φ”及びピッチ角“θ”の真値を示す。また、デルタ「Δ」付きの姿勢角は、それぞれロール角“φ”及びピッチ角“θ”の実測値に含まれ得る誤差(姿勢角誤差)を示す。“Δφ”はロール角誤差であり、“Δθ”はピッチ角誤差である。   In Expression (6), the posture angles with a hat “^” indicate the actual measurement values of the roll angle “φ” and the pitch angle “θ”, respectively, and the posture angles without the hat “^” are the roll angles “φ”. And the true value of the pitch angle “θ”. The posture angle with a delta “Δ” indicates an error (posture angle error) that can be included in the actual measurement values of the roll angle “φ” and the pitch angle “θ”. “Δφ” is a roll angle error, and “Δθ” is a pitch angle error.

式(6)を一次の項までテイラー展開すると、次式(7)が導かれる。

Figure 2012154769
When the equation (6) is Taylor-expanded to the first order term, the following equation (7) is derived.
Figure 2012154769

式(6)及び式(7)から、次式(8)が導かれる。

Figure 2012154769
From the equations (6) and (7), the following equation (8) is derived.
Figure 2012154769

式(8)の左辺について、第1項の“fn=(fx n,fy n,fz n)”は、Nフレーム加速度ベクトルの真値を示し、第2項のハット^付きの“Fn(φ,θ)=(Fx n(φ,θ),Fy n(φ,θ),Fz n(φ,θ))”は、Nフレーム加速度ベクトルの実測値を示す。 The left side of Equation (8), "f n = (f x n, f y n, f z n)" in the first term represents the true value of N frame acceleration vector, the second term hat ^ with the “F n (φ, θ) = (F x n (φ, θ), F y n (φ, θ), F z n (φ, θ))” indicates an actually measured value of the N frame acceleration vector.

Nフレーム加速度ベクトルの実測値であるハット^付きの“Fn(φ,θ)”は、加速度センサー3の検出加速度ベクトルであるBフレーム加速度ベクトル“fb=(fx b,fy b,fz b)”と、姿勢センサー4の検出姿勢角(φ,θ,Ψ)とを、それぞれ式(5)に代入することで算出することができる。つまり、Nフレーム加速度ベクトルの実測値“Fn(φ,θ)”は既知である。 N frame the actual measurement value is hat ^ with the acceleration vector "F n (φ, θ) " is, B-frame acceleration vector "f b is a detected acceleration vector of the acceleration sensor 3 = (f x b, f y b, f z b ) ”and the detected attitude angles (φ, θ, Ψ) of the attitude sensor 4 can be calculated by substituting them into the equation (5), respectively. That is, the actual measurement value “F n (φ, θ)” of the N frame acceleration vector is known.

それに対して、Nフレーム加速度ベクトルの真値“fn=(fx n,fy n,fz n)”は未知であり、その値を知ることはできない。しかし、本実施形態では、移動体の移動に係る参照ベクトルが外部から与えられる。そこで、参照ベクトルを用いて、絶対座標加速度ベクトルであるNフレーム加速度ベクトルを算出し、これをNフレーム加速度ベクトルの真値“fn=(fx n,fy n,fz n)”とみなして計算する。参照ベクトルが十分信頼できるものであれば、参照ベクトルから算出されるNフレーム加速度ベクトルを真値とみなして計算しても、十分な精度が保証される。 In contrast, the true value of N frame acceleration vector "f n = (f x n , f y n, f z n)" is unknown, it is impossible to know the value. However, in this embodiment, a reference vector related to the movement of the moving object is given from the outside. Therefore, by using a reference vector, calculates the N frame acceleration vector is an absolute coordinate acceleration vector, which the true value of N frame acceleration vector "f n = (f x n , f y n, f z n)" and Consider and calculate. If the reference vector is sufficiently reliable, sufficient accuracy is guaranteed even if the N frame acceleration vector calculated from the reference vector is regarded as a true value.

式(8)の右辺について、偏微分行列(ヤコビ行列)の各成分の値は、それぞれ式(5)を、ロール角及びピッチ角の実測値(ハット^付きの(φ,θ))で偏微分することで求められる。姿勢角の実測値は、姿勢センサー4により検出された姿勢角(検出姿勢角)である。   Regarding the right side of Equation (8), the values of each component of the partial differentiation matrix (Jacobi matrix) are obtained by substituting Equation (5) with the measured values of roll angle and pitch angle ((φ, θ) with a hat ^). It is obtained by differentiating. The measured value of the posture angle is the posture angle (detected posture angle) detected by the posture sensor 4.

式(8)において、例えば、左辺を“Δρ”と置き、右辺の偏微分行列を“H”と置き、右辺の姿勢角誤差を成分とする行列(姿勢角誤差行列)を“Δx”と置くと、式(8)は式(9)のように書き換えることができる。

Figure 2012154769
In Expression (8), for example, the left side is set as “Δρ”, the partial differential matrix on the right side is set as “H”, and the matrix (posture angle error matrix) including the posture angle error on the right side as a component is set as “Δx”. Equation (8) can be rewritten as Equation (9).
Figure 2012154769

式(9)における未知数は、姿勢角誤差行列“Δx”のみである。そのため、式(9)は、次式(10)のように解くことができる。

Figure 2012154769
The only unknown in equation (9) is the attitude angle error matrix “Δx”. Therefore, the equation (9) can be solved as the following equation (10).
Figure 2012154769

コンピューター(計算機)を用いた計算では、姿勢角誤差行列“Δx”を近似的に求めることになる。例えば、求根アルゴリズムの一種であり、反復法として知られるニュートン法を用いることで、姿勢角誤差行列“Δx”の近似解を求めることができる。   In the calculation using a computer (computer), the posture angle error matrix “Δx” is approximately obtained. For example, an approximate solution of the attitude angle error matrix “Δx” can be obtained by using a Newton method which is a kind of root finding algorithm and is known as an iterative method.

姿勢角誤差行列“Δx”が求まったならば、姿勢角誤差(ロール角誤差“Δφ”及びピッチ角誤差“Δθ”)を姿勢センサー4の検出姿勢角(ハット^付きの(φ,θ,Ψ))に加算することで、検出姿勢角を調整する。   When the posture angle error matrix “Δx” is obtained, the posture angle errors (roll angle error “Δφ” and pitch angle error “Δθ”) are detected by the posture angle detected by the posture sensor 4 ((φ, θ, Ψ with a hat ^). )) Is added to adjust the detected posture angle.

なお、本実施形態では、ロール角“φ”及びピッチ角“θ”に含まれる誤差を近似的に求め、これらの誤差を用いて、検出姿勢角のうちのロール角“φ”及びピッチ角“θ”を調整する。つまり、移動体の検出姿勢のうちのヨー成分については調整を行わず、ロール成分及びピッチ成分のみを調整することになる。   In this embodiment, errors included in the roll angle “φ” and the pitch angle “θ” are approximately obtained, and the roll angle “φ” and the pitch angle “of the detected attitude angles are used by using these errors. Adjust θ ”. That is, the yaw component in the detected posture of the moving body is not adjusted, and only the roll component and the pitch component are adjusted.

上記のようにして求めた調整姿勢角を用いてローカル座標系から絶対座標系への座標変換行列(変換係数)を求め、これを用いて加速度センサー3の検出加速度ベクトルを座標変換すると、検出加速度ベクトルに含まれ得る誤差が座標変換により低減される。つまり、元々検出加速度ベクトルに含まれていた誤差よりも誤差が小さくなる。その理由は、以下の通りである。   When the coordinate transformation matrix (transformation coefficient) from the local coordinate system to the absolute coordinate system is obtained using the adjustment attitude angle obtained as described above, and the detected acceleration vector of the acceleration sensor 3 is transformed using this, the detected acceleration is obtained. Errors that may be included in the vector are reduced by coordinate transformation. That is, the error is smaller than the error originally included in the detected acceleration vector. The reason is as follows.

加速度センサー3により検出されたローカル座標加速度ベクトルには、移動体の移動による誤差が含まれ得る。この誤差を含み得るローカル座標加速度ベクトルを座標変換した場合は、変換後の絶対座標加速度ベクトルにも当然に誤差が重畳されることとなる。   The local coordinate acceleration vector detected by the acceleration sensor 3 may include an error due to movement of the moving body. When the local coordinate acceleration vector that may contain this error is subjected to coordinate conversion, the error is naturally superimposed on the converted absolute coordinate acceleration vector.

しかし、上記の原理で説明したように、本実施形態では、加速度センサー3により検出されたローカル座標加速度ベクトルそのものを用いて移動体の検出姿勢を調整し、その調整結果である調整姿勢角を用いて座標変換行列を求めている。これは、加速度センサー3の検出加速度ベクトルに誤差が含まれていることを前提とし、言ってみれば検出加速度ベクトルと辻褄が合うような変換係数を求めることを意味する。つまり、検出加速度ベクトルの誤差を検出姿勢の誤差に含めて考えることで、検出加速度ベクトルに含まれる誤差を低減させようとする考え方である。   However, as described in the above principle, in this embodiment, the detected posture of the moving body is adjusted using the local coordinate acceleration vector itself detected by the acceleration sensor 3, and the adjustment posture angle that is the adjustment result is used. To obtain a coordinate transformation matrix. This presupposes that an error is included in the detected acceleration vector of the acceleration sensor 3, and in other words, it means obtaining a conversion coefficient that matches the detected acceleration vector. That is, it is an idea to reduce the error included in the detected acceleration vector by including the error in the detected acceleration vector in the error in the detected posture.

上記の姿勢調整部20による検出姿勢(検出姿勢角)の調整と、座標変換部30による変換係数(座標変換行列)の算出と、座標変換部30による座標変換とは、移動体の移動中に行うと好適である。勿論、移動体の停止時に行うこととしてもよいが、本実施形態の手法は、移動体の移動中に特に効果を発揮する。   The adjustment of the detected posture (detected posture angle) by the posture adjusting unit 20, the calculation of the conversion coefficient (coordinate conversion matrix) by the coordinate converting unit 30, and the coordinate conversion by the coordinate converting unit 30 are performed while the moving body is moving. It is preferred to do so. Of course, it may be performed when the moving body is stopped, but the method of the present embodiment is particularly effective during the movement of the moving body.

本実施形態の手法とは別に、移動体が停止している場合に加速度センサー3により検出される加速度ベクトルを用いて、移動体の姿勢を推定する手法が考案されている。具体的には、次式(11)及び(12)に従って、ロール角“φ”及びピッチ角“θ”を推定する。

Figure 2012154769
Figure 2012154769
Apart from the method of the present embodiment, a method has been devised that estimates the posture of the moving body using the acceleration vector detected by the acceleration sensor 3 when the moving body is stopped. Specifically, the roll angle “φ” and the pitch angle “θ” are estimated according to the following equations (11) and (12).
Figure 2012154769
Figure 2012154769

移動体が停止しているのであれば、移動体の姿勢は通常は変化しないため、式(11)及び(12)を用いることで、移動体の姿勢を正しく推定することができる。よって、移動体の停止時においては、式(11)及び(12)に従って推定したロール角“φ”及びピッチ角“θ”を用いて座標変換行列を求め、加速度センサー3の検出加速度ベクトルを座標変換することで、絶対座標加速度ベクトルを適切に求めることができる。   If the moving body is stopped, since the posture of the moving body does not normally change, the posture of the moving body can be correctly estimated by using the equations (11) and (12). Therefore, when the moving body is stopped, a coordinate transformation matrix is obtained using the roll angle “φ” and the pitch angle “θ” estimated according to the equations (11) and (12), and the detected acceleration vector of the acceleration sensor 3 is coordinated. By converting, the absolute coordinate acceleration vector can be obtained appropriately.

しかし、移動体の移動中は事態が異なる。移動体の移動により、移動体の姿勢は随時変化するようになる。すると、停止時に推定した姿勢角は正しくなくなり、停止時に推定した姿勢角から求めた座標変換行列は、加速度センサー3の検出加速度ベクトルとの関係で辻褄が合わなくなる。かかる座標変換行列を用いて検出加速度ベクトルを座標変換すると、検出加速度ベクトルに含まれる誤差を低減させることはできず、絶対座標加速度ベクトルには大きな誤差が残存してしまう。   However, the situation is different while the moving body is moving. Due to the movement of the moving body, the posture of the moving body changes at any time. Then, the posture angle estimated at the time of stopping becomes incorrect, and the coordinate transformation matrix obtained from the posture angle estimated at the time of stopping becomes inconsistent with the detected acceleration vector of the acceleration sensor 3. When the detected acceleration vector is coordinate-converted using such a coordinate conversion matrix, the error included in the detected acceleration vector cannot be reduced, and a large error remains in the absolute coordinate acceleration vector.

その点、本実施形態の手法では、加速度センサー3の検出加速度ベクトルと辻褄が合うような座標変換行列を求めるため、移動中に移動体の姿勢が随時変化しても問題はなく、検出加速度ベクトルに含まれる誤差を低減させて、誤差の影響が抑制された絶対座標加速度ベクトルを得ることができる。   In that respect, in the method of the present embodiment, since a coordinate transformation matrix that matches the detected acceleration vector of the acceleration sensor 3 is obtained, there is no problem even if the posture of the moving body changes at any time during movement. Thus, an absolute coordinate acceleration vector in which the influence of the error is suppressed can be obtained.

図1に示した加速度検出システム1は、本実施形態における加速度検出を実現するための最低限の構成要素を示したものである。実際には、加速度検出システムの種々のバリエーションが考えられる。そこで、加速度検出システムのバリエーションについて以下説明する。なお、加速度検出システム1と同一の機能ブロックについては同一の符号を付して説明を省略し、加速度検出システム1とは異なる部分を中心に説明する。   The acceleration detection system 1 shown in FIG. 1 shows the minimum components for realizing acceleration detection in the present embodiment. Actually, various variations of the acceleration detection system are conceivable. Therefore, variations of the acceleration detection system will be described below. In addition, about the functional block same as the acceleration detection system 1, the same code | symbol is attached | subjected and description is abbreviate | omitted and it demonstrates centering on a different part from the acceleration detection system 1. FIG.

図2は、第1の加速度検出システム1Aの説明図である。第1の加速度検出システム1Aは、加速度検出システム1の各機能ブロックの他に、姿勢調整用加速度ベクトル算出部40を備えて構成される。姿勢調整用加速度ベクトル算出部40は、移動体の移動に係る参照ベクトルを用いて、姿勢角調整用の加速度ベクトル(姿勢角調整用加速度ベクトル)を算出する。   FIG. 2 is an explanatory diagram of the first acceleration detection system 1A. The first acceleration detection system 1 </ b> A includes a posture adjustment acceleration vector calculation unit 40 in addition to the functional blocks of the acceleration detection system 1. The posture adjustment acceleration vector calculation unit 40 calculates a posture angle adjustment acceleration vector (posture angle adjustment acceleration vector) using a reference vector related to the movement of the moving object.

具体的には、姿勢調整用加速度ベクトル算出部40は、参照ベクトルとして外部入力した移動体の参照速度ベクトルを用いて絶対座標加速度ベクトルを算出する。例えば、参照速度ベクトルとして絶対座標速度ベクトルが入力される場合は、当該絶対座標速度ベクトルを時間微分することで、絶対座標加速度ベクトルを算出する。そして、算出した絶対座標加速度ベクトルを、姿勢調整用加速度ベクトルとして姿勢調整部20に出力する。   Specifically, the posture adjustment acceleration vector calculation unit 40 calculates the absolute coordinate acceleration vector using the reference velocity vector of the moving body externally input as the reference vector. For example, when an absolute coordinate velocity vector is input as the reference velocity vector, the absolute coordinate acceleration vector is calculated by differentiating the absolute coordinate velocity vector with respect to time. Then, the calculated absolute coordinate acceleration vector is output to the posture adjustment unit 20 as a posture adjustment acceleration vector.

図3は、第2の加速度検出システム1Bの説明図である。第2の加速度検出システム1Bは、処理ブロックとして、姿勢調整部20と、座標変換部30と、姿勢調整用加速度ベクトル算出部40と、速度ベクトル算出部50と、速度ベクトル補正部60とを有する。   FIG. 3 is an explanatory diagram of the second acceleration detection system 1B. The second acceleration detection system 1B includes an attitude adjustment unit 20, a coordinate conversion unit 30, an attitude adjustment acceleration vector calculation unit 40, a velocity vector calculation unit 50, and a velocity vector correction unit 60 as processing blocks. .

速度ベクトル算出部50は、座標変換部30から入力した絶対座標加速度ベクトルを積分及び継ぎ足し処理することで、絶対座標系における移動体の速度ベクトル(絶対座標速度ベクトル)を算出する。そして、算出した絶対座標速度ベクトルを速度ベクトル補正部60に出力する。ここで、「積分」とは、所定時間分の値を累積的に加算することを意味する。また、「継ぎ足す」とは、積分により得られた結果を最後に更新した結果に足し合わせることを意味する。   The velocity vector calculation unit 50 calculates the velocity vector (absolute coordinate velocity vector) of the moving body in the absolute coordinate system by integrating and adding the absolute coordinate acceleration vector input from the coordinate conversion unit 30. Then, the calculated absolute coordinate velocity vector is output to the velocity vector correction unit 60. Here, “integration” means that values for a predetermined time are cumulatively added. “Adding” means adding the result obtained by integration to the result of the last update.

具体的には、速度ベクトル算出部50は、所定時間分の絶対座標加速度ベクトルを積分することで、その所定時間の間における移動体の速度ベクトルの変化分を算出する。そして、算出した速度ベクトルの変化分を最後に更新した絶対座標速度ベクトルに継ぎ足すことで、絶対座標速度ベクトルを新たに算出・更新する。   Specifically, the velocity vector calculation unit 50 integrates the absolute coordinate acceleration vector for a predetermined time to calculate the change in the velocity vector of the moving body during the predetermined time. Then, the absolute coordinate velocity vector is newly calculated and updated by adding the calculated change in the velocity vector to the last updated absolute coordinate velocity vector.

速度ベクトル補正部60は、速度ベクトル算出部50から入力した絶対座標速度ベクトルを補正する。具体的には、移動体の移動に係る参照ベクトルとして外部入力した参照速度ベクトルを利用して、絶対座標速度ベクトルを補正する。参照速度ベクトルとしては、例えば絶対座標系における速度ベクトルを取得することができる。そのため、例えば、2種類の絶対座標速度ベクトルの平均処理を行うことで、絶対座標速度ベクトルを補正する。   The velocity vector correction unit 60 corrects the absolute coordinate velocity vector input from the velocity vector calculation unit 50. Specifically, the absolute coordinate velocity vector is corrected using a reference velocity vector input externally as a reference vector related to the movement of the moving object. As the reference velocity vector, for example, a velocity vector in an absolute coordinate system can be acquired. For this reason, for example, the absolute coordinate velocity vector is corrected by averaging two types of absolute coordinate velocity vectors.

この場合における平均処理としては、単純な算術平均を適用することもできるし、加重平均を適用することもできる。参照速度ベクトルを信頼して速度ベクトルを補正するのであれば、参照速度ベクトルの重みを速度ベクトル算出部50により算出された速度ベクトルの重みよりも高くして加重平均すればよい。逆に、速度ベクトル算出部50により算出された速度ベクトルを信頼して速度ベクトルを補正するのであれば、当該速度ベクトルの重みを参照速度ベクトルの重みよりも高くして加重平均すればよい。   As an averaging process in this case, a simple arithmetic average can be applied, or a weighted average can be applied. If the reference speed vector is trusted and the speed vector is corrected, the weight of the reference speed vector may be set higher than the weight of the speed vector calculated by the speed vector calculation unit 50 and weighted. Conversely, if the speed vector calculated by the speed vector calculation unit 50 is trusted and the speed vector is corrected, the weight of the speed vector may be set higher than the weight of the reference speed vector and weighted averaged.

速度ベクトル補正部60は、補正結果として得られる絶対座標速度ベクトル(以下、「絶対座標補正速度ベクトル」と称す。)を外部出力し、且つ、姿勢調整用加速度ベクトル算出部40にフィードバックする。姿勢調整用加速度ベクトル算出部40は、速度ベクトル補正部60から入力した絶対座標補正速度ベクトルを用いて姿勢調整用加速度ベクトルを算出し、姿勢調整部20に出力する。   The velocity vector correction unit 60 externally outputs an absolute coordinate velocity vector (hereinafter referred to as “absolute coordinate correction velocity vector”) obtained as a correction result, and feeds it back to the posture adjustment acceleration vector calculation unit 40. The posture adjustment acceleration vector calculation unit 40 calculates a posture adjustment acceleration vector using the absolute coordinate correction velocity vector input from the velocity vector correction unit 60, and outputs it to the posture adjustment unit 20.

図4は、第3の加速度検出システム1Cの説明図である。第3の加速度検出システム1Cは、第2の加速度検出システム1Bの変形例である。具体的には、第2の加速度検出システム1Bとは異なり、第3の加速度検出システム1Cでは、姿勢センサー4の検出姿勢が、姿勢調整部20と速度ベクトル補正部60とにそれぞれ入力されるように構成されている。   FIG. 4 is an explanatory diagram of the third acceleration detection system 1C. The third acceleration detection system 1C is a modification of the second acceleration detection system 1B. Specifically, unlike the second acceleration detection system 1B, in the third acceleration detection system 1C, the detected attitude of the attitude sensor 4 is input to the attitude adjustment unit 20 and the velocity vector correction unit 60, respectively. It is configured.

速度ベクトル補正部60は、移動体の移動に係る参照情報として外部入力した参照速度と、姿勢センサー4の検出姿勢とを用いて、速度ベクトル算出部50から入力した絶対座標速度ベクトルを補正する。例えば、車速検出システムから参照情報を取得可能な場合において、車速検出システムが移動体の速度の大きさのみを検出可能に構成されている場合は、参照情報として速度の大きさを取得し、姿勢センサー4の検出姿勢と併用して、絶対座標速度ベクトルを補正する。   The velocity vector correction unit 60 corrects the absolute coordinate velocity vector input from the velocity vector calculation unit 50 using the reference velocity externally input as the reference information related to the movement of the moving body and the detected posture of the posture sensor 4. For example, when reference information can be acquired from the vehicle speed detection system, if the vehicle speed detection system is configured to detect only the speed of the moving body, the speed information is acquired as the reference information, and the attitude The absolute coordinate velocity vector is corrected in combination with the detection posture of the sensor 4.

図5は、絶対座標速度ベクトルの補正の原理を説明するための図である。ここでは、絶対座標系としてNフレームを考え、東西方向(E軸)及び南北方向(N軸)でなる2次元平面において速度ベクトルを補正する場合について図示・説明する。図5において、横軸は東西方向の軸(E軸)を示しており、縦軸は南北方向の軸(N軸)を示している。   FIG. 5 is a diagram for explaining the principle of correcting the absolute coordinate velocity vector. Here, N frames are considered as an absolute coordinate system, and a case where a velocity vector is corrected in a two-dimensional plane composed of an east-west direction (E axis) and a north-south direction (N axis) will be illustrated and described. In FIG. 5, the horizontal axis indicates the east-west axis (E axis), and the vertical axis indicates the north-south axis (N axis).

Nフレーム加速度ベクトルの東西方向成分及び南北方向成分を積分して最後に更新したNフレーム速度ベクトルに継ぎ足すと、Nフレーム速度ベクトルの東西方向成分及び南北方向成分が得られる。図5において、E軸上に記載した矢印がNフレーム速度ベクトルの東西方向成分を示し、N軸上に記載した矢印がNフレーム速度ベクトルの南北方向成分を示す。   When the east-west direction component and the north-south direction component of the N frame acceleration vector are integrated and added to the last updated N frame velocity vector, the east-west direction component and the north-south direction component of the N frame velocity vector are obtained. In FIG. 5, the arrow described on the E axis indicates the east-west direction component of the N frame rate vector, and the arrow described on the N axis indicates the north-south direction component of the N frame rate vector.

この東西方向成分及び南北方向成分を合成することで、大きさが“|V1|”で、向きが“θ1”のNフレーム速度ベクトルが得られる。すなわち、動径を“|V1|”、偏角を“θ1”とする円座標でNフレーム速度ベクトルが定義される。但し、図5では、Nフレーム速度ベクトルの向きを、速度ベクトルとE軸との成す角度で表現している。上記のNフレーム速度ベクトルが、速度ベクトル算出部50により算出されるNフレーム速度ベクトルである。   By combining the east-west direction component and the north-south direction component, an N frame velocity vector having a magnitude of “| V1 |” and a direction of “θ1” is obtained. That is, the N frame velocity vector is defined by a circular coordinate having a moving radius of “| V1 |” and a declination of “θ1”. However, in FIG. 5, the direction of the N frame velocity vector is represented by an angle formed by the velocity vector and the E axis. The N frame rate vector is the N frame rate vector calculated by the rate vector calculating unit 50.

このNフレーム速度ベクトルの大きさ及び向きを、例えば次のように補正する。先ず、Nフレーム速度ベクトルの大きさ“|V1|”(第1の速度)と、外部入力した参照速度(第2の速度)とを平均処理することで合成する。そして、その合成結果として得られる速度を、Nフレーム補正速度ベクトルの大きさ“|V2|”とする。   The magnitude and direction of the N frame rate vector are corrected as follows, for example. First, the N frame rate vector magnitude “| V1 |” (first rate) and the externally input reference rate (second rate) are averaged and combined. The speed obtained as a result of the synthesis is assumed to be the magnitude “| V2 |” of the N frame correction speed vector.

次に、Nフレーム速度ベクトルの向き“θ1”を、姿勢センサー4の検出姿勢に基づいて、移動体の移動方向“θ2”に合わせるように補正する。Nフレーム速度ベクトルの向き“θ1”は、Nフレーム補正加速度ベクトルを積分及び継ぎ足し処理して間接的に求められる向きであるのに対し、移動方向“θ2”は、姿勢センサー4の検出姿勢から直接的に定まる移動体の向きである。そのため、上記の補正により、Nフレーム補正速度ベクトルの向きは、移動体の移動方向に沿ったより正確な向きとなる。   Next, the direction “θ1” of the N frame velocity vector is corrected based on the detected posture of the posture sensor 4 so as to match the moving direction “θ2” of the moving body. The direction “θ1” of the N frame velocity vector is a direction obtained indirectly by integrating and adding the N frame corrected acceleration vector, whereas the moving direction “θ2” is directly determined from the detected posture of the posture sensor 4. This is the direction of the moving object. Therefore, by the above correction, the direction of the N frame correction speed vector becomes a more accurate direction along the moving direction of the moving body.

2.実施例
次に、上記の加速度検出システムを適用した位置算出システム及びナビゲーションシステムの実施例について説明する。但し、本発明を適用可能な実施例が以下説明する実施例に限定されるわけではないことは勿論である。
2. Example Next, an example of a position calculation system and a navigation system to which the above-described acceleration detection system is applied will be described. However, it goes without saying that the embodiments to which the present invention can be applied are not limited to the embodiments described below.

2−1.システム構成
図6は、本実施例におけるナビゲーションシステムのシステム構成の説明図である。ナビゲーションシステムは、移動体の一種である四輪自動車(以下、単に「自動車」と称す。)に、加速度検出装置及び位置算出装置を具備した電子機器の一種であるカーナビゲーション装置1000が搭載されたシステムである。
2-1. System Configuration FIG. 6 is an explanatory diagram of the system configuration of the navigation system in the present embodiment. In the navigation system, a car navigation device 1000, which is a type of electronic device including an acceleration detection device and a position calculation device, is mounted on a four-wheeled vehicle (hereinafter simply referred to as "automobile") that is a type of mobile body. System.

カーナビゲーション装置1000は、自動車に搭載される電子機器であり、自動車の加速度ベクトルを算出する加速度検出装置及び自動車の位置を算出する位置算出装置を備えて構成される。また、カーナビゲーション装置1000は、加速度センサー3及びジャイロセンサー5を有するセンサーユニットとしてIMU(Inertial Measurement Unit)600を具備している。   The car navigation apparatus 1000 is an electronic device mounted on a car, and includes an acceleration detection device that calculates an acceleration vector of the car and a position calculation device that calculates the position of the car. In addition, the car navigation apparatus 1000 includes an IMU (Inertial Measurement Unit) 600 as a sensor unit having the acceleration sensor 3 and the gyro sensor 5.

IMU600は、慣性計測ユニットとして知られるセンサーユニットであり、自動車に対する相対的な座標系の加速度ベクトル及び角速度を検出して出力するように設計されている。本実施例では、この相対的な座標系がBフレームと一致するような姿勢でカーナビゲーション装置1000が自動車に設置・固定されているものとして説明する。   The IMU 600 is a sensor unit known as an inertial measurement unit, and is designed to detect and output an acceleration vector and an angular velocity of a coordinate system relative to the automobile. In the present embodiment, the car navigation apparatus 1000 will be described as being installed and fixed in an automobile in such a posture that the relative coordinate system coincides with the B frame.

また、カーナビゲーション装置1000は、衛星測位システムの一種であるGPSを利用したセンサーユニットとしてGPSユニット200を具備している。GPSユニット200は、GPS衛星から発信されているGPS衛星信号を利用して、例えばNフレームにおける自動車の位置及び速度ベクトルを計測して出力可能に構成されている。   In addition, the car navigation apparatus 1000 includes a GPS unit 200 as a sensor unit using GPS, which is a kind of satellite positioning system. The GPS unit 200 is configured to be able to measure and output, for example, the position and speed vector of an automobile in an N frame using a GPS satellite signal transmitted from a GPS satellite.

カーナビゲーション装置1000は、IMU600の検出結果を利用して慣性航法演算処理を行うことで、自動車の位置を算出する。そして、当該算出位置に対するマップマッチング処理を行うなどしてナビゲーション画面を生成し、生成したナビゲーション画面を表示部500であるディスプレイに表示させる。   The car navigation apparatus 1000 calculates the position of the automobile by performing inertial navigation calculation processing using the detection result of the IMU 600. Then, a navigation screen is generated by performing a map matching process on the calculated position, and the generated navigation screen is displayed on the display as the display unit 500.

2−2.機能構成
図7は、カーナビゲーション装置1000の機能構成の一例を示すブロック図である。カーナビゲーション装置1000は、主要な機能構成として、GPSアンテナ100と、GPSユニット200と、処理部300と、操作部400と、表示部500と、IMU600と、記憶部700とを備えて構成される。
2-2. Functional Configuration FIG. 7 is a block diagram illustrating an example of a functional configuration of the car navigation apparatus 1000. The car navigation apparatus 1000 includes a GPS antenna 100, a GPS unit 200, a processing unit 300, an operation unit 400, a display unit 500, an IMU 600, and a storage unit 700 as main functional configurations. .

GPSアンテナ100は、測位用衛星の一種であるGPS衛星から発信されているGPS衛星信号を含むRF(Radio Frequency)信号を受信するアンテナであり、受信信号をGPSユニット200に出力する。GPS衛星信号は、拡散符号の一種であるCA(Coarse and Acquisition)コードによって、スペクトラム拡散方式として知られるCDMA(Code Division Multiple Access)方式によって変調された1.57542[GHz]の通信信号である。CAコードは、コード長1023チップを1PNフレームとする繰返し周期1msの擬似ランダム雑音符号であり、GPS衛星毎に異なるものである。   The GPS antenna 100 is an antenna that receives an RF (Radio Frequency) signal including a GPS satellite signal transmitted from a GPS satellite that is a type of positioning satellite, and outputs a received signal to the GPS unit 200. The GPS satellite signal is a 1.57542 [GHz] communication signal modulated by a CDMA (Code Division Multiple Access) system known as a spread spectrum system by a CA (Coarse and Acquisition) code which is a kind of spreading code. The CA code is a pseudo-random noise code having a repetition period of 1 ms with a code length of 1023 chips as 1 PN frame, and is different for each GPS satellite.

GPSユニット200は、GPSアンテナ100から出力された信号に基づいて自動車の位置、移動速度及び移動方向を計測する回路或いは装置であり、いわゆるGPS受信装置に相当する機能ブロックである。図示を省略するが、GPSユニット200は、RF受信回路部と、ベースバンド処理回路部とを備えて構成される。なお、RF受信回路部と、ベースバンド処理回路部とは、それぞれ別のLSI(Large Scale Integration)として製造することも、1チップとして製造することも可能である。   The GPS unit 200 is a circuit or device that measures the position, moving speed, and moving direction of an automobile based on a signal output from the GPS antenna 100, and is a functional block corresponding to a so-called GPS receiving device. Although not shown, the GPS unit 200 includes an RF receiving circuit unit and a baseband processing circuit unit. Note that the RF receiving circuit unit and the baseband processing circuit unit can be manufactured as separate LSIs (Large Scale Integration) or can be manufactured as one chip.

ベースバンド処理回路部は、RF受信回路部から出力された受信信号に対して相関処理等を行ってGPS衛星信号を捕捉・追尾する。そして、GPS衛星信号から取り出した衛星軌道データや時刻データ等に基づいて、所定の位置算出計算を行って自動車の位置を算出する。また、ベースバンド処理回路部は、GPS衛星から受信したGPS衛星信号のドップラー周波数に基づいて、所定の速度算出計算を行って自動車の速度ベクトルを算出する。本実施例では、GPSユニット200は、自動車の位置(位置座標)及び速度ベクトルをNフレームにおいて算出して出力可能に構成されているものとして説明する。   The baseband processing circuit unit performs correlation processing or the like on the received signal output from the RF receiving circuit unit to capture and track the GPS satellite signal. Then, based on satellite orbit data, time data, and the like extracted from the GPS satellite signal, a predetermined position calculation is performed to calculate the position of the automobile. Further, the baseband processing circuit unit performs a predetermined speed calculation calculation based on the Doppler frequency of the GPS satellite signal received from the GPS satellite to calculate the speed vector of the automobile. In this embodiment, the GPS unit 200 will be described as being configured to be able to calculate and output the position (position coordinates) and speed vector of an automobile in N frames.

処理部300は、記憶部700に記憶されているシステムプログラム等の各種プログラムに従ってカーナビゲーション装置1000の各部を統括的に制御する制御装置であり、CPU(Central Processing Unit)等のプロセッサーを有して構成される。   The processing unit 300 is a control device that comprehensively controls each unit of the car navigation apparatus 1000 according to various programs such as a system program stored in the storage unit 700, and includes a processor such as a CPU (Central Processing Unit). Composed.

処理部300は、IMU600の検出結果を利用した慣性航法演算処理を行って自動車の位置(位置座標)を算出する。そして、算出した位置をもとに、自動車の現在位置を指し示した地図を表示部500に表示させる処理を行う。   The processing unit 300 performs inertial navigation calculation processing using the detection result of the IMU 600 to calculate the position (position coordinates) of the automobile. And based on the calculated position, the process which displays the map which pointed to the present position of a motor vehicle on the display part 500 is performed.

図8は、原理で説明した第1の加速度検出システム1Aを適用した位置算出システムである第1の位置算出システム2Aのシステム構成の一例を示す図であり、処理部300を機能ブロックとして表した概念図である。   FIG. 8 is a diagram illustrating an example of a system configuration of a first position calculation system 2A that is a position calculation system to which the first acceleration detection system 1A described in the principle is applied, and the processing unit 300 is represented as a functional block. It is a conceptual diagram.

処理部300は、姿勢算出部310と、姿勢調整部320と、座標変換部330と、姿勢調整用加速度ベクトル算出部340と、速度ベクトル算出部350と、位置算出部380とを機能部として有する。姿勢調整部320〜速度ベクトル算出部350は、原理で説明した加速度検出システムの姿勢調整部20〜速度ベクトル算出部50にそれぞれ対応する。   The processing unit 300 includes a posture calculation unit 310, a posture adjustment unit 320, a coordinate conversion unit 330, a posture adjustment acceleration vector calculation unit 340, a velocity vector calculation unit 350, and a position calculation unit 380 as functional units. . The posture adjustment unit 320 to the velocity vector calculation unit 350 correspond to the posture adjustment unit 20 to the velocity vector calculation unit 50 of the acceleration detection system described in the principle.

姿勢算出部310は、ジャイロセンサー5により検出されたBフレーム角速度を積分及び継ぎ足していくことで、自動車の姿勢角を算出する。すなわち、ジャイロセンサー5及び姿勢算出部310によって、原理で説明した加速度検出システムにおける姿勢センサー4が構成される。   The attitude calculation unit 310 calculates the attitude angle of the automobile by integrating and adding the B frame angular velocity detected by the gyro sensor 5. In other words, the gyro sensor 5 and the attitude calculation unit 310 constitute the attitude sensor 4 in the acceleration detection system described in the principle.

具体的には、姿勢算出部310は、ジャイロセンサー5により検出された所定時間分のBフレーム角速度を積分することで、その所定時間の間の移動体の姿勢角の変化方向及び変化分を算出する。所定時間はBフレーム角速度の検出時間間隔よりも長く設定され、通常は一定時間とされる。   Specifically, the posture calculation unit 310 calculates the change direction and the change amount of the posture angle of the moving body during the predetermined time by integrating the B frame angular velocity for a predetermined time detected by the gyro sensor 5. To do. The predetermined time is set longer than the detection time interval of the B frame angular velocity, and is usually set to a fixed time.

例えば、所定時間を100ミリ秒とし、Bフレーム角速度の検出時間間隔(より正確にはジャイロセンサー5から出力された信号のサンプリング時間間隔)が1〜10ミリ秒の間で適宜調整されるとした場合には、検出時間間隔に関わらず、100ミリ秒毎に移動体の姿勢角の変化方向及び変化分を算出する。そして、100ミリ秒毎に算出される姿勢角の変化方向及び変化分を最後に更新した姿勢角に継ぎ足すことで、最新の姿勢角を算出する。   For example, the predetermined time is set to 100 milliseconds, and the detection time interval of the B frame angular velocity (more precisely, the sampling time interval of the signal output from the gyro sensor 5) is appropriately adjusted between 1 to 10 milliseconds. In this case, the change direction and change amount of the posture angle of the moving body are calculated every 100 milliseconds regardless of the detection time interval. Then, the latest posture angle is calculated by adding the change direction and change amount of the posture angle calculated every 100 milliseconds to the last updated posture angle.

以下の説明では、姿勢算出部310が算出した姿勢角のことを「算出姿勢角」と称し、この算出姿勢角によって定まる自動車の姿勢のことを「算出姿勢」と称する。姿勢算出部310は、算出姿勢(算出姿勢角)を姿勢調整部320に出力する。   In the following description, the attitude angle calculated by the attitude calculation unit 310 is referred to as “calculated attitude angle”, and the attitude of the vehicle determined by the calculated attitude angle is referred to as “calculated attitude”. The posture calculation unit 310 outputs the calculated posture (calculated posture angle) to the posture adjustment unit 320.

姿勢調整用加速度ベクトル算出部340は、GPSユニット200からNフレームGPS速度ベクトルを参照速度ベクトルとして取得して、姿勢調整用加速度ベクトルを算出する。そして、算出した姿勢調整用加速度ベクトルを姿勢調整部320に出力する。   The posture adjustment acceleration vector calculation unit 340 acquires an N frame GPS velocity vector from the GPS unit 200 as a reference velocity vector, and calculates a posture adjustment acceleration vector. Then, the calculated attitude adjustment acceleration vector is output to attitude adjustment section 320.

位置算出部380は、速度ベクトル算出部350により算出された自動車のNフレーム速度ベクトルを用いて、Nフレームにおける自動車の位置を算出する。具体的には、Nフレーム速度ベクトルを積分及び継ぎ足し処理することで、自動車の位置変化分を算出する。そして、自動車の所与の移動基準位置(初期位置)を基準として、算出した位置変化分を累積的に加算していくことで、自動車の位置を算出・更新する。移動基準位置(初期位置)は、例えばナビゲーションの開始位置としてユーザーにより入力された位置とすることができる。   The position calculation unit 380 calculates the position of the vehicle in the N frame using the vehicle N frame speed vector calculated by the speed vector calculation unit 350. Specifically, the position change of the vehicle is calculated by integrating and adding the N frame speed vector. Then, the position of the automobile is calculated / updated by cumulatively adding the calculated position changes with the given movement reference position (initial position) of the automobile as a reference. The movement reference position (initial position) can be, for example, a position input by the user as a navigation start position.

図7の説明に戻り、操作部400は、例えばタッチパネルやボタンスイッチ等により構成される入力装置であり、押下されたキーやボタンの信号を処理部300に出力する。この操作部400の操作により、目的地の入力等の各種指示入力がなされる。   Returning to the description of FIG. 7, the operation unit 400 is an input device configured by, for example, a touch panel or a button switch, and outputs a pressed key or button signal to the processing unit 300. By operating the operation unit 400, various instructions such as destination input are input.

表示部500は、LCD(Liquid Crystal Display)等により構成され、処理部300から入力される表示信号に基づいた各種表示を行う表示装置である。表示部500には、ナビゲーション画面等が表示される。   The display unit 500 is a display device that includes an LCD (Liquid Crystal Display) or the like and performs various displays based on a display signal input from the processing unit 300. A navigation screen or the like is displayed on the display unit 500.

IMU600は、加速度センサー3及びジャイロセンサー5を備え、Bフレームの各軸における加速度ベクトル及び各軸の軸回りの角速度を検出可能に構成されている。なお、加速度センサー3とジャイロセンサー5とは、それぞれが独立したセンサーであってもよいし、一体型のセンサーであってもよい。   The IMU 600 includes an acceleration sensor 3 and a gyro sensor 5, and is configured to be able to detect an acceleration vector in each axis of the B frame and an angular velocity around each axis. The acceleration sensor 3 and the gyro sensor 5 may be independent sensors or may be integrated sensors.

記憶部700は、ROM(Read Only Memory)やフラッシュROM、RAM(Random Access Memory)等の記憶装置によって構成され、カーナビゲーション装置1000のシステムプログラムや、ナビゲーション機能等の各種機能を実現するための各種プログラム、データ等を記憶している。また、各種処理の処理中データ、処理結果などを一時的に記憶するワークエリアを有する。   The storage unit 700 is configured by a storage device such as a ROM (Read Only Memory), a flash ROM, or a RAM (Random Access Memory), and various types for realizing various functions such as a system program of the car navigation device 1000 and a navigation function. Stores programs, data, etc. In addition, it has a work area for temporarily storing data being processed and results of various processes.

2−3.データ構成
図7に示すように、記憶部700には、プログラムとして、処理部300により読み出され、ナビゲーション処理(図9参照)として実行されるナビゲーションプログラム701が記憶されている。ナビゲーション処理については、フローチャートを用いて詳細に後述する。
2-3. Data Configuration As shown in FIG. 7, the storage unit 700 stores a navigation program 701 that is read by the processing unit 300 and executed as a navigation process (see FIG. 9) as a program. The navigation process will be described later in detail using a flowchart.

また、記憶部700には、データとして、センサーデータ703と、参照データ705と、算出姿勢角データ707と、姿勢角誤差データ709と、調整姿勢角データ711と、移動ベクトルデータ713と、算出位置データ715とが記憶される。   The storage unit 700 also includes sensor data 703, reference data 705, calculated attitude angle data 707, attitude angle error data 709, adjusted attitude angle data 711, movement vector data 713, and calculated positions. Data 715 is stored.

センサーデータ703は、IMU600により検出されたBフレーム加速度ベクトル及びBフレーム角速度が記憶されたデータである。   The sensor data 703 is data in which the B frame acceleration vector and the B frame angular velocity detected by the IMU 600 are stored.

参照データ705は、自動車の移動に係る参照値が記憶されたデータである。本実施例では、GPSユニット200により計測された自動車のNフレーム速度ベクトルが参照データ705として記憶される。   The reference data 705 is data in which a reference value related to the movement of the automobile is stored. In this embodiment, the N frame speed vector of the automobile measured by the GPS unit 200 is stored as the reference data 705.

算出姿勢角データ707は、姿勢算出部310により算出された自動車の姿勢角が記憶されたデータである。ロール成分、ピッチ成分及びヨー成分の姿勢角がそれぞれ算出されて記憶される。   The calculated attitude angle data 707 is data in which the attitude angle of the automobile calculated by the attitude calculation unit 310 is stored. The attitude angles of the roll component, pitch component, and yaw component are calculated and stored.

姿勢角誤差データ709は、姿勢調整部320が姿勢角誤差の推定処理を行って推定した姿勢角誤差が記憶されたデータである。例えば、ロール成分及びピッチ成分の誤差が算出されて記憶される。   The posture angle error data 709 is data in which the posture angle error estimated by the posture adjustment unit 320 performing the posture angle error estimation process is stored. For example, the roll component and pitch component errors are calculated and stored.

調整姿勢角データ711は、姿勢調整部320が算出姿勢角の調整を行うことで取得した自動車の調整姿勢角が記憶されたデータである。   The adjusted attitude angle data 711 is data in which the adjusted attitude angle of the automobile acquired by the attitude adjustment unit 320 adjusting the calculated attitude angle is stored.

移動ベクトルデータ713は、自動車の移動状態を示すベクトルのデータであり、座標変換部330により算出されたNフレーム加速度ベクトルや、速度ベクトル算出部350により算出されたNフレーム速度ベクトルといったデータが記憶される。   The movement vector data 713 is vector data indicating the moving state of the automobile, and stores data such as an N frame acceleration vector calculated by the coordinate conversion unit 330 and an N frame velocity vector calculated by the velocity vector calculation unit 350. The

算出位置データ715は、位置算出部380により算出された自動車の位置が記憶されたデータである。   The calculated position data 715 is data in which the position of the automobile calculated by the position calculation unit 380 is stored.

2−4.処理の流れ
図9は、記憶部700に記憶されているナビゲーションプログラム701が処理部300により読み出されて実行されることで、カーナビゲーション装置1000において実行されるナビゲーション処理の流れを示すフローチャートである。以下のナビゲーション処理では、IMU600により検出された加速度ベクトル及び角速度が、記憶部700のセンサーデータ703に随時記憶されるものとする。
2-4. Processing Flow FIG. 9 is a flowchart showing the flow of navigation processing executed in the car navigation device 1000 when the navigation program 701 stored in the storage unit 700 is read and executed by the processing unit 300. . In the following navigation processing, it is assumed that the acceleration vector and angular velocity detected by the IMU 600 are stored in the sensor data 703 of the storage unit 700 as needed.

最初に、姿勢算出部310は、姿勢角算出処理を行う(ステップA1)。具体的には、ジャイロセンサー5の検出結果を積分及び継ぎ足し処理して自動車の姿勢角を算出し、算出姿勢角データ707として記憶部700に記憶させる。   First, the posture calculation unit 310 performs a posture angle calculation process (step A1). Specifically, the attitude angle of the vehicle is calculated by integrating and adding the detection results of the gyro sensor 5 and stored in the storage unit 700 as calculated attitude angle data 707.

次いで、処理部300は、GPSユニット200から参照データ705を取得したか否かを判定し(ステップA3)、参照データ705を取得しなかったと判定した場合は(ステップA3;No)、ステップA9へと処理を移行する。   Next, the processing unit 300 determines whether or not the reference data 705 has been acquired from the GPS unit 200 (step A3). If it is determined that the reference data 705 has not been acquired (step A3; No), the processing unit 300 proceeds to step A9. And migrate the process.

また、参照データ705を取得したと判定した場合は(ステップA3;Yes)、姿勢調整用加速度ベクトル算出部340が、姿勢調整用加速度ベクトル算出処理を行う(ステップA5)。すなわち、GPSユニット200から取得したNフレームGPS速度ベクトルを時間微分することで、Nフレームにおける姿勢調整用加速度ベクトルを算出する。   If it is determined that the reference data 705 has been acquired (step A3; Yes), the posture adjustment acceleration vector calculation unit 340 performs posture adjustment acceleration vector calculation processing (step A5). That is, the N-frame GPS velocity vector acquired from the GPS unit 200 is time-differentiated to calculate the attitude adjustment acceleration vector in the N frame.

次いで、姿勢調整部320は、姿勢角誤差推定処理を行う(ステップA7)。具体的には、加速度センサー3により検出されたBフレーム加速度ベクトルと、ステップA5において算出した姿勢調整用加速度ベクトルと、ステップA1において算出した姿勢角とを用いて、例えばニュートン法を用いた反復計算を行って姿勢角誤差を推定する。そして、推定した姿勢角誤差を姿勢角誤差データ709として記憶部700に記憶させる。   Next, the posture adjustment unit 320 performs posture angle error estimation processing (step A7). Specifically, using the B frame acceleration vector detected by the acceleration sensor 3, the posture adjustment acceleration vector calculated in step A5, and the posture angle calculated in step A1, for example, iterative calculation using the Newton method. To estimate the attitude angle error. Then, the estimated posture angle error is stored in the storage unit 700 as posture angle error data 709.

その後、姿勢調整部320は、姿勢角誤差データ709に記憶された最新の姿勢角誤差を用いて、算出姿勢角データ707に記憶された最新の算出姿勢角を調整し、調整姿勢角データ711として記憶部700に記憶させる(ステップA9)。   Thereafter, the posture adjustment unit 320 adjusts the latest calculated posture angle stored in the calculated posture angle data 707 using the latest posture angle error stored in the posture angle error data 709, and serves as adjusted posture angle data 711. It memorize | stores in the memory | storage part 700 (step A9).

次いで、座標変換部330は、Nフレーム加速度ベクトル算出処理を行う(ステップA11)。具体的には、調整姿勢角データ711に記憶された最新の調整姿勢角を用いてBフレームからNフレームへの座標変換行列を求める。そして、当該座標変換行列を用いて加速度センサー3により検出されたBフレーム加速度ベクトルをNフレーム加速度ベクトルに変換し、移動ベクトルデータ713として記憶部700に記憶させる。   Next, the coordinate conversion unit 330 performs an N frame acceleration vector calculation process (step A11). Specifically, a coordinate conversion matrix from the B frame to the N frame is obtained using the latest adjustment attitude angle stored in the adjustment attitude angle data 711. Then, the B frame acceleration vector detected by the acceleration sensor 3 is converted into an N frame acceleration vector using the coordinate conversion matrix, and is stored in the storage unit 700 as movement vector data 713.

その後、速度ベクトル算出部350は、Nフレーム速度ベクトル算出処理を行う(ステップA13)。すなわち、所定時間分のNフレーム加速度ベクトルを積分し、最新のNフレーム速度ベクトルに継ぎ足すことで、Nフレーム速度ベクトルを新たに算出する。そして、移動ベクトルデータ713として記憶部700に記憶させる。   Thereafter, the speed vector calculation unit 350 performs an N frame speed vector calculation process (step A13). That is, the N frame acceleration vector for a predetermined time is integrated and added to the latest N frame velocity vector to newly calculate the N frame velocity vector. Then, it is stored in the storage unit 700 as movement vector data 713.

次いで、位置算出部380は、位置算出処理を行う(ステップA15)。具体的には、移動ベクトルデータ713に記憶されている所定時間分のNフレーム速度ベクトルを積分及び継ぎ足し処理することで、移動体の位置変化分を算出する。そして、算出した位置変化分を、算出位置データ715に記憶されている最新の算出位置に加算することで、自動車の位置(以下、「慣性航法演算位置」と称す。)を新たに算出・更新する。   Next, the position calculation unit 380 performs position calculation processing (step A15). Specifically, the position change of the moving body is calculated by integrating and adding N frame velocity vectors for a predetermined time stored in the movement vector data 713. The calculated position change is added to the latest calculated position stored in the calculated position data 715 to newly calculate / update the position of the automobile (hereinafter referred to as “inertial navigation calculation position”). To do.

なお、ステップA15の位置算出処理において、GPSユニット200からGPS算出位置を取得し、取得したGPS算出位置を加味して最終的な自動車の位置(最終位置)を決定するように構成することも可能である。例えば、慣性航法演算位置とGPS算出位置との平均処理を行うこととしてもよい。また、定期的なタイミングにおいて、GPS算出位置で慣性航法演算位置をリセットすることとしてもよい。この場合は、GPS算出位置を移動基準位置として、新たに慣性航法演算を開始することとなる。   In addition, in the position calculation process of step A15, it is also possible to acquire the GPS calculated position from the GPS unit 200, and to determine the final position of the automobile (final position) in consideration of the acquired GPS calculated position. It is. For example, an average process between the inertial navigation calculation position and the GPS calculation position may be performed. Further, the inertial navigation calculation position may be reset at the GPS calculation position at regular timing. In this case, the inertial navigation calculation is newly started with the GPS calculated position as the movement reference position.

その後、位置算出部380は、ステップA15で算出した最終位置に対してマップマッチング処理等を行い、表示部500のナビゲーション画面を更新する(ステップA17)。なお、マップマッチングに係る処理については従来公知であるため、ここでは詳細な説明を省略する。   Thereafter, the position calculation unit 380 performs a map matching process on the final position calculated in step A15, and updates the navigation screen of the display unit 500 (step A17). Since the process related to map matching is conventionally known, detailed description thereof is omitted here.

次いで、処理部300は、ナビゲーション処理を終了するか否かを判定する(ステップA19)。例えば、操作部400を介してユーザーによりナビゲーションの終了指示操作がなされた場合に、ナビゲーション処理を終了すると判定する。まだ処理を終了しないと判定した場合は(ステップA19;No)、処理部300は、ステップA1に戻る。また、終了すると判定した場合は(ステップA19;Yes)、処理部300は、ナビゲーション処理を終了する。   Next, the processing unit 300 determines whether or not to end the navigation process (step A19). For example, when a navigation end instruction operation is performed by the user via the operation unit 400, it is determined that the navigation process is to be ended. When it determines with not complete | finishing a process yet (step A19; No), the process part 300 returns to step A1. If it is determined that the process is to be ended (step A19; Yes), the processing unit 300 ends the navigation process.

2−5.実験結果
図10は、上述した第1の位置算出システム2Aにおいて実際に位置算出を行った実験結果の一例を示す図である。自動車を予め定められた周回経路に沿って1周させ、その場合に算出される位置を東西南北の2次元平面にプロットした。
2-5. Experimental Results FIG. 10 is a diagram illustrating an example of experimental results obtained by actually performing position calculation in the above-described first position calculation system 2A. The automobile was caused to make one turn along a predetermined circulation route, and the positions calculated in that case were plotted on a two-dimensional plane of east, west, south, and north.

図10において、横軸は東西方向、縦軸は南北方向(単位はメートル)をそれぞれ示している。東西方向0m、南北方向0mの位置がスタート地点であり、出発地点から北方向に向かい、大きく反時計回り方向に1周周回するような経路を辿った。ゴール地点はスタート地点と同じである。また、自動車の真の軌跡を点線で、姿勢角の調整を行わずに自動車の位置を求めた軌跡を細実線で、姿勢角の調整を行って自動車の位置を求めた軌跡を太実線でそれぞれ示している。   In FIG. 10, the horizontal axis indicates the east-west direction, and the vertical axis indicates the north-south direction (unit: meters). The starting point was 0m in the east-west direction and 0m in the north-south direction, and the route followed the route from the departure point to the north direction and one round in the counterclockwise direction. The goal point is the same as the starting point. In addition, the true trajectory of the car is indicated by a dotted line, the trajectory obtained by determining the car position without adjusting the attitude angle is indicated by a thin solid line, and the trajectory obtained by adjusting the attitude angle and obtaining the car position is indicated by a thick solid line. Show.

この図を見ると、姿勢角の調整を行わずに自動車の位置を求めた場合は、最初のうちはほぼ真の軌跡に沿った位置が得られているものの、時間経過に伴って徐々に真の軌跡から外れていき、最終的には大きな位置ずれが発生していることがわかる。これは、姿勢算出部310により算出された自動車の姿勢角に誤差が含まれているために、当該姿勢角を用いてBフレーム加速度ベクトルを座標変換すると、Nフレーム加速度ベクトルに誤差が重畳されてしまうためである。   As can be seen from this figure, when the position of the car is obtained without adjusting the attitude angle, a position along the true trajectory is obtained at first, but it gradually increases with time. It can be seen that a large positional deviation has finally occurred. This is because an error is included in the attitude angle of the vehicle calculated by the attitude calculation unit 310, and when the B frame acceleration vector is coordinate-converted using the attitude angle, the error is superimposed on the N frame acceleration vector. It is because it ends.

それに対して、姿勢角の調整を行って自動車の位置を求めた場合は、ほぼ真の軌跡に沿った結果が得られていることがわかる。時間が経過しても真の軌跡から大きく外れることはなく、時間経過に伴う誤差の累積が防止されていることがわかる。   On the other hand, when the position of the automobile is obtained by adjusting the attitude angle, it can be seen that a result along a substantially true locus is obtained. Even if time passes, it does not deviate greatly from the true trajectory, and it can be seen that accumulation of errors with the passage of time is prevented.

図11は、上記の位置算出における姿勢角の時間変化を示すグラフである。図11(A)はロール角“φ”の時間変化を示し、図11(B)はピッチ角“θ”の時間変化を示し、図11(C)はヨー角“Ψ”の時間変化を示す。   FIG. 11 is a graph showing the time change of the posture angle in the position calculation. 11A shows the time change of the roll angle “φ”, FIG. 11B shows the time change of the pitch angle “θ”, and FIG. 11C shows the time change of the yaw angle “Ψ”. .

各図において、自動車の真の姿勢角を点線で示す。ロール角“φ”及びピッチ角“θ”の真値はゼロである。また、調整を行っていない姿勢角(算出姿勢角)を細実線で示し、調整を行った姿勢角(調整姿勢角)を太実線で示す。なお、本実施形態ではヨー成分については調整を行っていないため、ロール成分及びピッチ成分にのみ調整結果が示されている。   In each figure, the true attitude angle of the automobile is indicated by a dotted line. The true values of the roll angle “φ” and the pitch angle “θ” are zero. In addition, a posture angle that has not been adjusted (calculated posture angle) is indicated by a thin solid line, and a posture angle that has been adjusted (adjusted posture angle) is indicated by a thick solid line. In the present embodiment, since the yaw component is not adjusted, the adjustment results are shown only for the roll component and the pitch component.

この結果を見ると、調整を行っていない姿勢角については、ロール角“φ”及びピッチ角“θ”の何れも、時間経過に伴い、真値から値がずれていく傾向があることがわかる。時間経過に伴って、姿勢角の真値からのずれ量(姿勢角誤差)の振幅が大きくなっていくことが特徴的である。図10に示したように、この算出姿勢角の誤差は、検出加速度ベクトルの誤差と相まって、位置算出の正確性を大きく低下させる要因となっている。   As can be seen from the results, for the posture angle that is not adjusted, both the roll angle “φ” and the pitch angle “θ” tend to deviate from the true value with time. . It is characteristic that the amplitude of the deviation amount (posture angle error) from the true value of the posture angle increases with time. As shown in FIG. 10, the error in the calculated attitude angle is coupled with the error in the detected acceleration vector and is a factor that greatly reduces the accuracy of position calculation.

それに対して、調整姿勢角については、ロール角“φ”及びピッチ角“θ”の何れも、時間が経過しても、ほぼ真値に沿った値が得られていることがわかる。図10の結果から明らかなように、本実施形態の手法で調整を施した姿勢角を用いて検出加速度を座標変換し、慣性航法演算を行って位置を算出することで、位置算出の正確性が向上する。   On the other hand, with regard to the adjustment posture angle, it can be seen that both the roll angle “φ” and the pitch angle “θ” have values that are substantially in line with the true value even if time passes. As is clear from the results of FIG. 10, the accuracy of position calculation is obtained by performing coordinate conversion of detected acceleration using the attitude angle adjusted by the method of the present embodiment, and calculating the position by performing inertial navigation calculation. Will improve.

2−6.作用効果
加速度検出システム1Aにおいて、姿勢算出部310は、自動車に設置されたジャイロセンサー5が検出したBフレーム角速度を用いて、自動車の姿勢(姿勢角)を算出する。姿勢調整部320は、姿勢算出部310により算出された自動車の姿勢を、自動車に設置された加速度センサー3が検出したBフレーム加速度ベクトルとGPSユニット200が計測した参照速度ベクトルとを用いて調整する。すなわち、姿勢調整用加速度ベクトル算出部340により、参照速度ベクトルを用いて自動車のNフレーム加速度ベクトルが算出される。姿勢調整部320は、Bフレーム加速度ベクトルとNフレーム加速度ベクトルとを入力値として、姿勢角誤差を推定する所定の推定演算を行って姿勢角誤差を推定する。そして、姿勢調整部320は、推定した姿勢角誤差を用いて姿勢角を調整する。
2-6. Effects In the acceleration detection system 1A, the posture calculation unit 310 calculates the posture (posture angle) of the vehicle using the B frame angular velocity detected by the gyro sensor 5 installed in the vehicle. The attitude adjustment unit 320 adjusts the attitude of the vehicle calculated by the attitude calculation unit 310 using the B frame acceleration vector detected by the acceleration sensor 3 installed in the vehicle and the reference velocity vector measured by the GPS unit 200. . That is, the attitude adjustment acceleration vector calculation unit 340 calculates the N frame acceleration vector of the automobile using the reference speed vector. The posture adjustment unit 320 estimates a posture angle error by performing a predetermined estimation calculation for estimating a posture angle error using the B frame acceleration vector and the N frame acceleration vector as input values. Then, posture adjustment unit 320 adjusts the posture angle using the estimated posture angle error.

加速度センサー3が検出する加速度ベクトルは、加速度センサー3に対応付けられたローカルな座標系(例えばBフレーム)における加速度ベクトルである。それに対し、例えば慣性航法演算においては、移動体の移動空間を定める絶対的な座標系(例えばNフレーム)を基準として演算を行うことが一般的である。ジャイロセンサー5の検出結果には誤差が含まれ得るため、これを用いて算出される自動車の姿勢角にも誤差が含まれ得る。そこで、ジャイロセンサー5の検出結果を利用して算出した自動車の姿勢角を、加速度センサー3の検出加速度ベクトルとGPSユニット200の計測結果から求めた加速度ベクトルとを用いて調整する。その上で、ローカル座標系(Bフレーム)から絶対座標系(Nフレーム)への座標変換行列を求めることで、絶対座標系における加速度ベクトルを適切に求めることができる。   The acceleration vector detected by the acceleration sensor 3 is an acceleration vector in a local coordinate system (for example, B frame) associated with the acceleration sensor 3. On the other hand, for example, in inertial navigation calculation, the calculation is generally performed based on an absolute coordinate system (for example, N frame) that defines the moving space of the moving body. Since the detection result of the gyro sensor 5 may include an error, an error may also be included in the attitude angle of the automobile calculated using this. Therefore, the attitude angle of the vehicle calculated using the detection result of the gyro sensor 5 is adjusted using the detected acceleration vector of the acceleration sensor 3 and the acceleration vector obtained from the measurement result of the GPS unit 200. Then, by obtaining a coordinate transformation matrix from the local coordinate system (B frame) to the absolute coordinate system (N frame), the acceleration vector in the absolute coordinate system can be obtained appropriately.

座標変換部330は、姿勢調整部320により調整された姿勢角を用いて、BフレームからNフレームへの座標変換行列を求める。この座標変換行列は、加速度センサー3により検出されたBフレーム加速度ベクトルに含まれ得る誤差を座標変換で低減可能な座標変換行列である。つまり、調整姿勢角に基づき定義される座標変換行列を用いてBフレーム加速度ベクトルを座標変換することで、Bフレーム加速度ベクトルに含まれ得る誤差を低減させ、誤差の影響が抑制されたNフレーム加速度ベクトルを得ることができる。   The coordinate conversion unit 330 uses the posture angle adjusted by the posture adjustment unit 320 to obtain a coordinate conversion matrix from the B frame to the N frame. This coordinate transformation matrix is a coordinate transformation matrix that can reduce an error that can be included in the B frame acceleration vector detected by the acceleration sensor 3 by coordinate transformation. That is, by converting the B frame acceleration vector using the coordinate conversion matrix defined based on the adjustment attitude angle, the error that can be included in the B frame acceleration vector is reduced, and the influence of the error is suppressed. You can get a vector.

速度ベクトル算出部350は、座標変換部330により算出されたNフレーム加速度ベクトルを用いてNフレーム速度ベクトルを算出する。誤差が低減されたNフレーム加速度ベクトルを用いてNフレーム速度ベクトルを求めるため、加速度センサー3が検出したBフレーム加速度ベクトルに含まれていた誤差がNフレーム速度ベクトルに重畳することが防止され、時間経過に伴う誤差の累積を防止することができる。   The velocity vector calculation unit 350 calculates an N frame velocity vector using the N frame acceleration vector calculated by the coordinate conversion unit 330. Since the N frame velocity vector is obtained using the N frame acceleration vector in which the error is reduced, it is possible to prevent the error included in the B frame acceleration vector detected by the acceleration sensor 3 from being superimposed on the N frame velocity vector. Accumulation of errors with progress can be prevented.

本実施形態の加速度検出方法では、加速度センサー3の検出加速度ベクトルに含まれ得る誤差を座標変換で低減可能な変換係数(変換行列)を求めて検出加速度ベクトルを座標変換するため、検出加速度ベクトルそれ自体の補正を必要とせずに、誤差の影響が抑制された絶対座標系の加速度ベクトルを検出することができる。そして、所与の移動基準位置と、検出した絶対座標系における加速度ベクトルとを用いることで、移動体の位置をより正確に求めることができる。   In the acceleration detection method of the present embodiment, the detected acceleration vector is converted into coordinates by converting a detected acceleration vector by obtaining a conversion coefficient (conversion matrix) that can reduce errors that can be included in the detected acceleration vector of the acceleration sensor 3 by coordinate conversion. The acceleration vector of the absolute coordinate system in which the influence of the error is suppressed can be detected without requiring correction of itself. And the position of a moving body can be calculated | required more correctly by using the given movement reference position and the acceleration vector in the detected absolute coordinate system.

3.変形例
本発明を適用可能な実施例は、上記の実施例に限定されることなく、本発明の趣旨を逸脱しない範囲で適宜変更可能であることは勿論である。以下、変形例について説明するが、上記の実施例と同一の構成要素については同一の符号を付して説明を省略し、上記の実施例とは異なる部分を中心に説明する。
3. Modifications Embodiments to which the present invention can be applied are not limited to the above-described embodiments, and can be changed as appropriate without departing from the spirit of the present invention. Hereinafter, although a modification is demonstrated, the same code | symbol is attached | subjected about the component same as said Example, description is abbreviate | omitted, and it demonstrates centering on a different part from said Example.

3−1.位置算出システム
原理で説明した各種のバリエーションの加速度検出システムを適用して位置算出システムを構築することも勿論可能である。
3-1. Position calculation system Of course, it is possible to construct a position calculation system by applying the acceleration detection system of various variations described in the principle.

図12は、第2の加速度検出システム1Bを適用した第2の位置算出システム2Bのシステム構成の一例を示す図である。第2の位置算出システム2Bでは、処理部300は、姿勢算出部310と、姿勢調整部320と、座標変換部330と、姿勢調整用加速度ベクトル算出部340と、速度ベクトル算出部350と、速度ベクトル補正部360と、位置算出部380とを機能ブロックとして有する。   FIG. 12 is a diagram illustrating an example of a system configuration of a second position calculation system 2B to which the second acceleration detection system 1B is applied. In the second position calculation system 2B, the processing unit 300 includes an attitude calculation unit 310, an attitude adjustment unit 320, a coordinate conversion unit 330, an attitude adjustment acceleration vector calculation unit 340, a speed vector calculation unit 350, a speed A vector correction unit 360 and a position calculation unit 380 are provided as functional blocks.

速度ベクトル補正部360は、GPSユニット200から参照速度ベクトルとしてNフレームGPS速度ベクトルを取得し、速度ベクトル算出部350により算出されたNフレーム速度ベクトルを補正する。そして、その補正結果として得られるNフレーム補正速度ベクトルを位置算出部380に出力する。また、Nフレーム補正速度ベクトルを姿勢調整用加速度ベクトル算出部340にフィードバックする。   The speed vector correction unit 360 acquires an N frame GPS speed vector as a reference speed vector from the GPS unit 200 and corrects the N frame speed vector calculated by the speed vector calculation unit 350. Then, the N frame correction speed vector obtained as the correction result is output to position calculation section 380. Further, the N frame correction speed vector is fed back to the attitude adjustment acceleration vector calculation unit 340.

図13は、第2の位置算出システム2Bを利用して実際に位置算出を行った実験結果の一例を示す図である。図の見方は図10と同じである。この図を見ると、第2の位置算出システム2Bを利用して位置算出を行った場合も、自動車の真の軌跡に沿った軌跡が得られており、誤差の累積が抑えられていることがわかる。   FIG. 13 is a diagram illustrating an example of an experimental result of actual position calculation using the second position calculation system 2B. The way of viewing the figure is the same as in FIG. As can be seen from the figure, even when the position calculation is performed using the second position calculation system 2B, a trajectory along the true trajectory of the automobile is obtained, and accumulation of errors is suppressed. Recognize.

図14は、図13の位置算出結果に対応する姿勢角の時間変化を示すグラフである。図の見方は図12と同じである。第2の位置算出システム2Bにより姿勢角の調整を行った場合も、時間経過に伴う姿勢角誤差の増大が抑えられ、真の姿勢角に近い姿勢角が得られていることがわかる。   FIG. 14 is a graph showing the temporal change of the posture angle corresponding to the position calculation result of FIG. The way of viewing the figure is the same as in FIG. It can be seen that even when the posture angle is adjusted by the second position calculation system 2B, an increase in the posture angle error with the passage of time is suppressed, and a posture angle close to the true posture angle is obtained.

図15は、第3の加速度検出システム1Cを適用した第3の位置算出システム2Cのシステム構成の一例を示す図である。第3の位置算出システム2Cでは、姿勢算出部310により算出された自動車の姿勢(算出姿勢)が速度ベクトル補正部360に入力され、且つ、車速検出システム800により検出された自動車の速度が速度ベクトル補正部360に入力されるように構成されている。   FIG. 15 is a diagram illustrating an example of a system configuration of a third position calculation system 2C to which the third acceleration detection system 1C is applied. In the third position calculation system 2C, the vehicle attitude (calculated attitude) calculated by the attitude calculation unit 310 is input to the speed vector correction unit 360, and the vehicle speed detected by the vehicle speed detection system 800 is the speed vector. It is configured to be input to the correction unit 360.

車速検出システム800は、自動車の速度を検出するシステムであり、例えば車速パルスを備えて構成される。車速検出システム800は、車軸の回転数に応じた電圧又はパルスを発生させることで車速を検出する。   The vehicle speed detection system 800 is a system that detects the speed of an automobile, and is configured to include vehicle speed pulses, for example. The vehicle speed detection system 800 detects the vehicle speed by generating a voltage or pulse corresponding to the rotational speed of the axle.

速度ベクトル補正部360は、姿勢算出部310から入力した算出姿勢と、車速検出システム800から入力した車速とを用いて、速度ベクトル算出部350により算出されたNフレーム速度ベクトルを補正する。そして、その補正結果であるNフレーム補正速度ベクトルを位置算出部380に出力する。   The speed vector correction unit 360 corrects the N frame speed vector calculated by the speed vector calculation unit 350 using the calculated posture input from the posture calculation unit 310 and the vehicle speed input from the vehicle speed detection system 800. Then, an N frame correction speed vector, which is the correction result, is output to position calculation section 380.

3−2.姿勢センサー
上記の実施例では、ジャイロセンサーを姿勢センサーとして例示したが、姿勢センサーはこれに限られない。例えば、地磁気センサーとして知られる方位センサーを利用して移動体の姿勢を検出することとしてもよい。また、ジャイロセンサーに方位センサーや加速度センサーを組み合わせて姿勢センサーを構成してもよい。
3-2. Attitude Sensor In the above embodiment, the gyro sensor is exemplified as the attitude sensor, but the attitude sensor is not limited to this. For example, the orientation of the moving body may be detected using an orientation sensor known as a geomagnetic sensor. Further, a posture sensor may be configured by combining a gyro sensor with an orientation sensor or an acceleration sensor.

3−3.座標系
上記の実施例では、ローカル座標系の一例として機体座標系(Bフレーム)を挙げ、絶対座標系の一例として北東下座標系(Nフレーム)を挙げて説明した。しかし、これらの座標系はあくまでも一例であり、他の座標系を適用することとしてもよいことは勿論である。例えば、絶対座標系として、東北上座標系として知られるENU(East North Up)座標系や、地球中心地球固定座標系として知られるECEF(Earth Centered Earth Fixed)座標系を適用してもよい。
3-3. Coordinate System In the above embodiment, the machine coordinate system (B frame) is given as an example of the local coordinate system, and the northeast lower coordinate system (N frame) is given as an example of the absolute coordinate system. However, these coordinate systems are merely examples, and it is needless to say that other coordinate systems may be applied. For example, an ENU (East North Up) coordinate system known as the Tohoku upper coordinate system or an ECEF (Earth Centered Earth Fixed) coordinate system known as the earth-centered earth fixed coordinate system may be applied as the absolute coordinate system.

3−4.検出姿勢の調整
上記の実施形態では、姿勢センサーの検出姿勢のうちのロール成分及びピッチ成分の両方について調整を行うものとして説明したが、ロール成分又はピッチ成分の何れか一方のみについて調整を行うこととしてもよい。例えば、ロール成分のみについて調整を行う場合は、式(8)をロール角“φ”の式に書き換えてロール角誤差“Δφ”の推定演算を行えばよい。また、ピッチ成分のみについて調整を行う場合は、式(8)をピッチ角“θ”の式に書き換えてピッチ角誤差“Δθ”の推定演算を行えばよい。
3-4. Adjustment of detected posture In the above-described embodiment, it has been described that both the roll component and the pitch component of the detected posture of the posture sensor are adjusted. However, only one of the roll component and the pitch component is adjusted. It is good. For example, when adjusting only the roll component, Equation (8) may be rewritten to the equation of the roll angle “φ” and the roll angle error “Δφ” may be estimated. Further, when adjusting only the pitch component, the formula (8) may be rewritten with the formula of the pitch angle “θ” and the pitch angle error “Δθ” may be estimated.

3−5.姿勢角誤差の推定演算
また、上記の実施形態で説明した姿勢角誤差の推定演算はあくまでも一例であり、適宜設計変更可能である。例えば、上記の実施形態では、式(8)〜式(10)に基づいて、求根アルゴリズムの一種であるニュートン法を用いて姿勢角誤差を推定するものとして説明した。求根アルゴリズムとしては、ニュートン法の他にも、割線法や挟み撃ち法といったアルゴリズムを適用することができることは勿論である。
3-5. Attitude Angle Error Estimation Calculation The attitude angle error estimation calculation described in the above embodiment is merely an example, and the design can be changed as appropriate. For example, in the above embodiment, the posture angle error is estimated based on the equations (8) to (10) using the Newton method which is a kind of root finding algorithm. Of course, as the root finding algorithm, in addition to the Newton method, algorithms such as the secant method and the pinch-shot method can be applied.

また、上記の実施形態とは異なり、カルマンフィルターを利用して姿勢角誤差を推定するように構成することも可能である。この場合は、例えば、姿勢角誤差を状態推定値(状態ベクトル)として設定し、加速度センサーの検出加速度ベクトルと移動体の移動に係る参照ベクトルから求めた姿勢調整用加速度ベクトルとの差分を観測値(観測ベクトル)として設定すればよい。そして、予測演算(時刻更新)と補正演算(観測更新)とを繰り返し行って、最適な状態推定値(姿勢角誤差)を推定する。   In addition, unlike the above-described embodiment, the posture angle error can be estimated using a Kalman filter. In this case, for example, the posture angle error is set as a state estimated value (state vector), and the difference between the detected acceleration vector of the acceleration sensor and the posture adjustment acceleration vector obtained from the reference vector related to the movement of the moving object is an observed value. What is necessary is just to set as (observation vector). Then, a prediction calculation (time update) and a correction calculation (observation update) are repeatedly performed to estimate an optimum state estimated value (attitude angle error).

3−6.電子機器
上記の実施例では、四輪自動車に搭載するナビゲーション装置に本発明を適用した場合を例に挙げて説明したが、本発明を適用可能な電子機器はこれに限られるわけではない。例えば、二輪自動車に搭載するナビゲーション装置に適用してもよいし、携帯型ナビゲーション装置に適用することとしてもよい。
3-6. Electronic Device In the above embodiment, the case where the present invention is applied to a navigation device mounted on a four-wheeled vehicle has been described as an example. However, an electronic device to which the present invention can be applied is not limited thereto. For example, the present invention may be applied to a navigation device mounted on a two-wheeled vehicle, or may be applied to a portable navigation device.

また、ナビゲーション以外の用途の電子機器についても本発明を同様に適用可能であることは勿論である。例えば、携帯型電話機やパソコン、PDA(Personal Digital Assistant)といった他の電子機器についても本発明を同様に適用して、当該電子機器の加速度検出及び位置算出を実現することが可能である。   Of course, the present invention can be similarly applied to electronic devices for purposes other than navigation. For example, the present invention can be similarly applied to other electronic devices such as a mobile phone, a personal computer, and a PDA (Personal Digital Assistant), and acceleration detection and position calculation of the electronic device can be realized.

3−7.衛星測位システム
上述した実施形態では、衛星測位システムとしてGPSを例に挙げて説明したが、WAAS(Wide Area Augmentation System)、QZSS(Quasi Zenith Satellite System)、GLONASS(GLObal NAvigation Satellite System)、GALILEO等の他の衛星測位システムであってもよいことは勿論である。
3-7. Satellite positioning system In the above-described embodiment, the GPS has been described as an example of the satellite positioning system. Of course, other satellite positioning systems may be used.

1 加速度検出システム、 2 位置算出システム、 3 加速度センサー、 4 姿勢センサー、 5 ジャイロセンサー、 20 姿勢調整部、 30 座標変換部、 40 姿勢調整用加速度ベクトル算出部、 50 速度ベクトル算出部、 60 速度ベクトル補正部、 100 GPSアンテナ、 200 GPSユニット、 300 処理部、 400 操作部、 500 表示部、 600 IMU、 700 記憶部 DESCRIPTION OF SYMBOLS 1 Acceleration detection system, 2 Position calculation system, 3 Acceleration sensor, 4 Attitude sensor, 5 Gyro sensor, 20 Attitude adjustment part, 30 Coordinate conversion part, 40 Attitude adjustment acceleration vector calculation part, 50 Speed vector calculation part, 60 Speed vector Correction unit, 100 GPS antenna, 200 GPS unit, 300 processing unit, 400 operation unit, 500 display unit, 600 IMU, 700 storage unit

Claims (8)

移動体に設置された姿勢センサーの検出姿勢を、前記移動体に設置された加速度センサーの検出加速度ベクトルと前記移動体の移動に係る所与の参照ベクトルとを用いて調整することと、
前記調整された前記検出姿勢を用いて、前記加速度センサーに対応付けられたローカル座標系から絶対座標系への座標変換に係る変換係数を求めることと、
前記変換係数を用いて、前記検出加速度ベクトルを前記絶対座標系における加速度ベクトルに座標変換することと、
を含む加速度検出方法。
Adjusting the detected attitude of the attitude sensor installed in the moving body using a detected acceleration vector of the acceleration sensor installed in the moving body and a given reference vector related to the movement of the moving body;
Using the adjusted detection attitude to obtain a conversion coefficient for coordinate conversion from a local coordinate system associated with the acceleration sensor to an absolute coordinate system;
Using the conversion coefficient, coordinate-transforming the detected acceleration vector into an acceleration vector in the absolute coordinate system;
An acceleration detection method including:
前記変換係数を求めることは、前記検出加速度ベクトルに含まれ得る誤差を小さくすることが可能な前記変換係数を求めることである、
請求項1に記載の加速度検出方法。
Obtaining the conversion coefficient is to obtain the conversion coefficient capable of reducing an error that may be included in the detected acceleration vector.
The acceleration detection method according to claim 1.
前記調整することと、前記変換係数を求めることと、前記座標変換することと、を前記移動体の移動中に行う、
請求項1又は2に記載の加速度検出方法。
Performing the adjustment, obtaining the conversion coefficient, and converting the coordinates during the movement of the moving body;
The acceleration detection method according to claim 1 or 2.
前記調整することは、前記検出姿勢のうちのヨー成分については調整を行わず、ロール成分及びピッチ成分の少なくとも一方を調整することを含む、
請求項1〜3の何れか一項に記載の加速度検出方法。
The adjusting includes adjusting at least one of a roll component and a pitch component without adjusting the yaw component of the detected posture,
The acceleration detection method as described in any one of Claims 1-3.
前記参照ベクトルは、前記移動体に設置された所定の速度計測システムによって計測された前記移動体の速度ベクトルを含み、
前記調整することは、
前記速度ベクトルを用いて、前記移動体の加速度ベクトルを算出することと、
前記検出加速度ベクトルと前記算出した加速度ベクトルとを用いて、前記検出姿勢を調整することと、
を含む、
請求項1〜4の何れか一項に記載の加速度検出方法。
The reference vector includes a velocity vector of the moving body measured by a predetermined speed measuring system installed on the moving body,
The adjustment is
Calculating an acceleration vector of the moving object using the velocity vector;
Adjusting the detected posture using the detected acceleration vector and the calculated acceleration vector;
including,
The acceleration detection method as described in any one of Claims 1-4.
前記速度計測システムは、車速検出システム、又は、衛星測位システムである、
請求項5に記載の加速度検出方法。
The speed measurement system is a vehicle speed detection system or a satellite positioning system.
The acceleration detection method according to claim 5.
請求項1〜6の何れか一項に記載の加速度検出方法を用いて前記絶対座標系における加速度ベクトルを検出することと、
所与の移動基準位置と前記検出された加速度ベクトルとを用いて、前記絶対座標系における前記移動体の位置を算出することと、
を含む位置算出方法。
Detecting an acceleration vector in the absolute coordinate system using the acceleration detection method according to claim 1;
Calculating a position of the moving body in the absolute coordinate system using a given movement reference position and the detected acceleration vector;
Position calculation method including
移動体に設置された姿勢センサーの検出姿勢を、前記移動体に設置された加速度センサーの検出加速度ベクトルと前記移動体の移動に係る所与の参照ベクトルとを用いて調整する調整部と、
前記調整された前記検出姿勢を用いて、前記加速度センサーに対応付けられたローカル座標系から絶対座標系への座標変換に係る変換係数を求める変換係数算出部と、
前記変換係数を用いて、前記検出加速度ベクトルを前記絶対座標系における加速度ベクトルに座標変換する変換部と、
を備えた加速度検出装置。
An adjustment unit that adjusts the detected posture of the posture sensor installed in the moving body using a detected acceleration vector of the acceleration sensor installed in the moving body and a given reference vector related to the movement of the moving body;
A conversion coefficient calculation unit for obtaining a conversion coefficient related to coordinate conversion from the local coordinate system to the absolute coordinate system associated with the acceleration sensor, using the adjusted detection posture;
A conversion unit that converts the detected acceleration vector into an acceleration vector in the absolute coordinate system using the conversion coefficient;
Acceleration detection apparatus comprising:
JP2011013749A 2011-01-26 2011-01-26 Acceleration detection method, position calculation method and acceleration detection device Withdrawn JP2012154769A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2011013749A JP2012154769A (en) 2011-01-26 2011-01-26 Acceleration detection method, position calculation method and acceleration detection device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011013749A JP2012154769A (en) 2011-01-26 2011-01-26 Acceleration detection method, position calculation method and acceleration detection device

Publications (1)

Publication Number Publication Date
JP2012154769A true JP2012154769A (en) 2012-08-16

Family

ID=46836643

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011013749A Withdrawn JP2012154769A (en) 2011-01-26 2011-01-26 Acceleration detection method, position calculation method and acceleration detection device

Country Status (1)

Country Link
JP (1) JP2012154769A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104316079A (en) * 2014-09-30 2015-01-28 北京航天控制仪器研究所 Drop point precision estimation method for inertia measurement system based on rocket sled test
CN104848857A (en) * 2015-04-30 2015-08-19 北京航天控制仪器研究所 Method for automatically distributing accuracy indexes of ballistic missile inertia measurement system
CN105334349A (en) * 2014-08-12 2016-02-17 上海新跃仪表厂 Closed-loop test system of liquid floated accelerometer and test method of test system
CN105894126A (en) * 2016-04-25 2016-08-24 王景国 Impact point prediction model generation and positioning methods for rocket debris
JP2016211906A (en) * 2015-05-01 2016-12-15 セイコーエプソン株式会社 Gradient measurement method and apparatus, electronic apparatus and program
JP2016211904A (en) * 2015-05-01 2016-12-15 セイコーエプソン株式会社 Error measurement method and apparatus, bias correction method and apparatus, electronic apparatus and program
CN111323043A (en) * 2020-03-26 2020-06-23 深圳市创客火科技有限公司 Sensor data processing method and system
CN113038821A (en) * 2018-10-30 2021-06-25 东京计器株式会社 Error correction device
CN113793509A (en) * 2018-03-12 2021-12-14 深圳鼎然信息科技有限公司 Road condition analysis method, device, equipment and medium based on motion sensor

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105334349A (en) * 2014-08-12 2016-02-17 上海新跃仪表厂 Closed-loop test system of liquid floated accelerometer and test method of test system
CN104316079A (en) * 2014-09-30 2015-01-28 北京航天控制仪器研究所 Drop point precision estimation method for inertia measurement system based on rocket sled test
CN104848857A (en) * 2015-04-30 2015-08-19 北京航天控制仪器研究所 Method for automatically distributing accuracy indexes of ballistic missile inertia measurement system
CN104848857B (en) * 2015-04-30 2017-11-28 北京航天控制仪器研究所 Ballistic missile inertial measurement system precision index auto-allocation method
JP2016211906A (en) * 2015-05-01 2016-12-15 セイコーエプソン株式会社 Gradient measurement method and apparatus, electronic apparatus and program
JP2016211904A (en) * 2015-05-01 2016-12-15 セイコーエプソン株式会社 Error measurement method and apparatus, bias correction method and apparatus, electronic apparatus and program
CN105894126A (en) * 2016-04-25 2016-08-24 王景国 Impact point prediction model generation and positioning methods for rocket debris
CN113793509A (en) * 2018-03-12 2021-12-14 深圳鼎然信息科技有限公司 Road condition analysis method, device, equipment and medium based on motion sensor
CN113793509B (en) * 2018-03-12 2024-04-16 深圳鼎然信息科技有限公司 Road condition analysis method, device, equipment and medium based on motion sensor
CN113038821A (en) * 2018-10-30 2021-06-25 东京计器株式会社 Error correction device
CN113038821B (en) * 2018-10-30 2022-12-06 东京计器株式会社 Error correction device
CN111323043A (en) * 2020-03-26 2020-06-23 深圳市创客火科技有限公司 Sensor data processing method and system

Similar Documents

Publication Publication Date Title
JP2012154769A (en) Acceleration detection method, position calculation method and acceleration detection device
US8781776B2 (en) Position calculating method and position calculating device
US9494428B2 (en) Attitude determination method, position calculation method, and attitude determination device
US9791575B2 (en) GNSS and inertial navigation system utilizing relative yaw as an observable for an ins filter
US10641625B2 (en) Method and apparatus for calibrating a magnetic sensor
JP5606656B2 (en) Positioning device
US8131462B2 (en) Vehicle guidance and sensor bias determination
JP5521531B2 (en) Position calculation method and position calculation system
JP6083279B2 (en) Movement status information calculation method and movement status information calculation device
US9759567B2 (en) Position calculation method and position calculation device
JP5742450B2 (en) Position calculation method and position calculation apparatus
US9121929B2 (en) Position output device using satellite navigation system
JP2012008096A (en) Bias estimation method, attitude estimation method, bias estimation device and attitude estimation device
JP2009236532A (en) Method for geolocation, program, and apparatus for geolocation
WO2016203744A1 (en) Positioning device
JP2013170903A (en) Measurement apparatus and measurement method
JP5892845B2 (en) Calibration quality determination apparatus and method
JP2012002678A (en) Acceleration calculation method, position calculation method, and acceleration calculation device
JP2014219340A (en) Offset correction method and offset correction device
JP4343581B2 (en) Moving body posture detection device
JP2015135349A (en) Bias estimation method, attitude estimation method, bias estimation device and attitude estimation device
JP2015102330A (en) Movement information calculation device, movement information calculation method, movement information calculation program, and mobile body
JP2011128138A (en) Method and system for calculating position
Jiang et al. Design of a Pedestrian Navigation System based on the Fusion INS with GPS
JP2013195325A (en) Positioning method and device for moving object

Legal Events

Date Code Title Description
A300 Withdrawal of application because of no request for examination

Free format text: JAPANESE INTERMEDIATE CODE: A300

Effective date: 20140401