JP2004226374A - 位置計測システム、位置計測方法、プログラム及び航走体運用システム - Google Patents
位置計測システム、位置計測方法、プログラム及び航走体運用システム Download PDFInfo
- Publication number
- JP2004226374A JP2004226374A JP2003017741A JP2003017741A JP2004226374A JP 2004226374 A JP2004226374 A JP 2004226374A JP 2003017741 A JP2003017741 A JP 2003017741A JP 2003017741 A JP2003017741 A JP 2003017741A JP 2004226374 A JP2004226374 A JP 2004226374A
- Authority
- JP
- Japan
- Prior art keywords
- information
- time
- mother ship
- position information
- speed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 82
- 238000000691 measurement method Methods 0.000 title description 13
- 238000004364 calculation method Methods 0.000 claims description 44
- 238000000034 method Methods 0.000 claims description 29
- 230000005540 biological transmission Effects 0.000 claims description 28
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 27
- 230000033001 locomotion Effects 0.000 claims description 22
- 239000006185 dispersion Substances 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 11
- 230000005856 abnormality Effects 0.000 description 5
- 230000014509 gene expression Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000013535 sea water Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 230000032258 transport Effects 0.000 description 1
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Navigation (AREA)
- Traffic Control Systems (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
【課題】母船から離れて航行する航走体の位置を、誤差やばらつきを少なくし、正確に把握する。
【解決手段】母船1の位置としての第1位置情報を出力する測位装置31と、その第1位置情報と、母船1から送出された航走体2の母船1に対する相対位置としての相対位置情報と、航走体2の対地速度としての速度情報と、航走体2の向かう方位としての方位情報とに基づいて、航走体2の位置としての第2位置情報を推定する位置計測装置4とを具備する位置計測システムを用いる。これにより、母船から離れて航行する航走体の位置を正確に把握することができる。
【選択図】 図1
【解決手段】母船1の位置としての第1位置情報を出力する測位装置31と、その第1位置情報と、母船1から送出された航走体2の母船1に対する相対位置としての相対位置情報と、航走体2の対地速度としての速度情報と、航走体2の向かう方位としての方位情報とに基づいて、航走体2の位置としての第2位置情報を推定する位置計測装置4とを具備する位置計測システムを用いる。これにより、母船から離れて航行する航走体の位置を正確に把握することができる。
【選択図】 図1
Description
【0001】
【発明の属する技術分野】
本発明は、位置計測システム及び位置計測方法等に係り、特に、母船と母船に制御されながら航行する航走体とを有する系における位置計測システム、位置計測方法、当該方法をコンピューターに実行させるプログラム及び航走体運用システムに関する。
【0002】
【従来の技術】
水中で行なう各種探査においては、基地となる母船により制御された無人の航走体が、実際の探査作業を行なうことが多い。その場合、母船は、航走体の正確な絶対位置(以下、単に「位置」ともいう)を把握している必要がある。
【0003】
そのような位置把握を行なう方法として、以下のような方法が知られている。すなわち、まずGPS(Global Positioning System:衛星航法システム)のような測位方法を用いて母船の位置を把握する。同時に、航走体のピンガー(海中の物標や地形の状況などの探査のための音波発信器)のような水中音響装置が発する音響(ピンガー音など)の受信状態により、航走体の母船に対する相対位置を把握する。そして、母船の位置と母船に対する相対位置とに基づいて、航走体の(絶対)位置を把握する方法である。
【0004】
上記方法では、航走体の位置計測の誤差が大きいという問題があった。例えば、水中音響装置を用いた場合には、海中の状態により突然に遠くの場所に位置計算結果が飛んでしまうことが起こる。それに対応するために、位置計算結果が飛ばないように、現在位置から一定距離内への位置更新のみを認めるロジックを入れる。そうすることで、位置計測結果が遠くの場所に飛ぶことをある程度回避できる。
【0005】
しかし、初期位置が正しく計測されていないと、ずっと正しくない位置を出力しつづけることになる。また、ピンガーのような水中音響による母船と航走体間の位置計測が途絶する(航走体からの信号を母船が受信できない)と、航走体の位置が求められない。
【0006】
このため、航走体の位置を、誤差やばらつきを少なく、正確に把握することが可能な技術が求められている。また、母船の測位装置に異常が生じた場合や、水中音響による相対位置の計測が困難な状況でも航走体の位置を正確に把握することが可能な技術が求められている。また、水中のようなGPS信号が届かない場所にいる航走体の位置を正確に把握する技術が求められている。
【0007】
【発明が解決しようとする課題】
従って、本発明の目的は、母船から離れて航行する航走体の位置を正確に把握可能な位置計測システム、位置計測方法、当該方法をコンピューターに実行させるプログラム及び航走体運用システムを提供することである。
【0008】
また、本発明の他の目的は、母船に搭載された測位装置に問題が生じても航走体の位置を正確に把握可能な位置計測システム、位置計測方法、当該方法をコンピューターに実行させるプログラム及び航走体運用システムを提供することである。
【0009】
本発明の更に他の目的は、水中音響による航走体の相対位置の計測が困難な場合でも航走体の位置を正確に把握可能な位置計測システム、位置計測方法、当該方法をコンピューターに実行させるプログラム及び航走体運用システムを提供することである。
【0010】
本発明の別の目的は、GPS信号の届かない水中においても、航走体の位置を正確に求めることが可能な位置計測システム、位置計測方法、当該方法をコンピューターに実行させるプログラム及び航走体運用システムを提供することである。
【0011】
【課題を解決するための手段】
以下に、[発明の実施の形態]で使用される番号・符号を用いて、課題を解決するための手段を説明する。これらの番号・符号は、[特許請求の範囲]の記載と[発明の実施の形態]との対応関係を明らかにするために付加されたものである。ただし、それらの番号・符号を、[特許請求の範囲]に記載されている発明の技術的範囲の解釈に用いてはならない。
【0012】
従って、上記課題を解決するために、本発明の位置計測システムは、母船(1)の位置としての第1位置情報を出力する測位装置(31)と、その第1位置情報と、母船(1)から送出された航走体(2)の母船(1)に対する相対位置としての相対位置情報と、航走体(2)の対地速度としての速度情報と、航走体(2)の向かう方位としての方位情報とに基づいて、航走体(2)の位置としての第2位置情報を推定する位置計測装置(4)とを具備する。
【0013】
また、本発明の位置計測システムは、位置計測装置(4)が、航走体(2)から送信されるその速度情報とその方位情報とを受信する送受信部(8)を具備する。
【0014】
また、本発明の位置計測システムは、位置計測装置(4)が、航走体(2)から発せられる探査信号を受信し、その探査信号の到着時刻としての第2時刻情報と、その探査信号の母船(1)に対する方向とを出力する受波器(6)と、航走体(2)から発せられるその探査信号の発信時刻としての第1時刻情報と、その第2時刻情報と、その方向とに基づいて、その相対位置情報を算出する制御部(5)とを更に具備し、送受信部(8)は、更に、その第1時刻情報を受信する。
【0015】
また、本発明の位置計測システムは、制御部(5)が、更に、その第1位置情報と、その相対位置情報と、その速度情報と、その方位情報とに基づいて、その第2位置情報を推定する。
【0016】
更に、本発明の位置計測システムは、制御部(5)が、カルマンフィルタを用いてその第2位置情報を推定する。
【0017】
更に、本発明の位置計測システムは、制御部(5)が、予め設定された演算時間ΔT毎のその第2位置情報を下記式(1)乃至式(3)に基づいて計算し、下記式(4)によりその第2位置情報を推定する。ここで、
【数3】
ただし、添え字k:k番目の演算時間、Kk:カルマンフィルタゲイン、XS:第1位置情報と相対位置情報とに基づいて求める航走体位置計測値、X^:第2位置情報、Xk^:演算時間ΔT毎の第2位置情報、φ:方位情報、A,B,C,D:航走体の運動モデルの状態空間表現、u:速度情報、Q:システム雑音wの分散、R:測定誤差nの分散、N:システム雑音wと測定誤差nの共分散である。
【0018】
更に、本発明の位置計測システムは、制御部(5)が、そのΔT毎のその第2位置情報を、前記式(1)及び式(2)を用いず、前記式(3)と予め設定されたカルマンフィルタゲインKとを用いて計算し、計算されたそのΔT毎のその第2位置情報に基づいて、前記式(4)でその第2位置情報を推定する。
【0019】
更に、本発明の位置計測システムは、母船(1)から送出された航走体(2)の対地速度としての速度情報と、航走体(2)の向かう方位としての方位情報とを航走体(2)から受信する送受信部(8)と、航走体(2)から発せられる探査信号を受信し、その探査信号の到着時刻としての第2時刻情報と、その探査信号の母船(1)に対する方向とを出力する受波器(6)と、航走体(2)から発せられる探査信号の発信時刻としての第1時刻情報と、その第2時刻情報と、その方向とに基づいて、その相対位置情報を算出する制御部(5)とを有する位置計測装置(4)を具備し、送受信部(8)は、更に、その第1時刻情報を受信し、制御部(5)は、更に、その相対位置情報と、その速度情報と、その方位情報とに基づいて、航走体(2)の位置としての第2位置情報を推定する。
【0020】
更に、本発明の位置計測システムは、制御部(5)が、カルマンフィルタを用いてその第2位置情報を推定する。
【0021】
更に、本発明の位置計測システムは、制御部(5)が、予め設定された演算時間ΔT毎のその第2位置情報を下記式(5)乃至式(7)に基づいて計算し、下記式(8)によりその第2位置情報を推定する。ここで、
【数4】
ただし、添え字k:k番目の演算時間、Kk:カルマンフィルタゲイン、XS:相対位置情報に基づいて求める航走体位置計測値、X^:第2位置情報、Xk^:演算時間ΔT毎の第2位置情報、φ:方位情報、A,B,C,D:航走体の運動モデルの状態空間表現、u:速度情報、Q:システム雑音wの分散、R:測定誤差nの分散、N:システム雑音wと測定誤差nの共分散である。
【0022】
更に、本発明の位置計測システムは、その速度情報が、航走体(2)の対地速度検出器(11)及び対水速度検出器(13)の少なくとも一方の計測結果に基づいている。
【0023】
更に、本発明の位置計測システムは、母船(1)は船舶であり、航走体(2)は水中を探査する潜水艇である。
【0024】
上記課題を解決するために、本発明の航走体運用システムは、上記各項のいずれか一項に記載の位置計測システムと、その位置計測システムを有する母船(1)と、母船(1)に制御され運動する航走体(2)とを具備する。
【0025】
上記課題を解決するために、本発明の位置計測方法は、母船(1)の位置としての第1位置情報を取得するステップと、母船(1)から送出された航走体(2)の対地速度としての速度情報と、航走体(2)の向かう方位としての方位情報と、航走体(2)から発せられるその探査信号の発信時刻としての第1時刻情報とを航走体(2)から取得するステップと、その探査信号に基づいて、その探査信号の到着時刻としての第2時刻情報と、その探査信号の母船(1)に対する方向とを取得するステップと、その第1時刻情報と、その第2時刻情報と、その方向とに基づいて、航走体(2)の母船(1)に対する位置である相対位置情報を算出するステップと、その第1位置情報と、その相対位置情報と、その速度情報と、その方位情報とに基づいて、カルマンフィルタを用いて航走体(2)の位置としての第2位置情報を推定するステップとを具備する。
【0026】
上記課題を解決するために、本発明に関わるプログラムは、母船(1)の位置としての第1位置情報を取得するステップと、母船(1)から送出された航走体(2)の対地速度としての速度情報と、航走体(2)の向かう方位としての方位情報と、航走体(2)から発せられるその探査信号の発信時刻としての第1時刻情報とを航走体(2)から取得するステップと、その探査信号の到着時刻としての第2時刻情報と、その探査信号の母船(1)に対する方向とを取得するステップと、その第1時刻情報と、その第2時刻情報と、その方向とに基づいて、航走体(2)の母船(1)に対する位置である相対位置情報を算出するステップと、その第1位置情報と、その相対位置情報と、その速度情報と、その方位情報とに基づいて、カルマンフィルタを用いて航走体(2)の位置としての第2位置情報を推定するステップとをコンピューターに実行させる。
【0027】
【発明の実施の形態】
以下、本発明である位置計測システム、位置計測方法、当該方法をコンピューターに実行させるプログラム、及び航走体運用システムの実施の形態に関して、添付図面を参照して説明する。
本実施の形態において、母船に制御され、水中探査に使用される航走体を例に示して説明する。ただし、他の母船に制御され、母船を離れて動作を行なう航走体(水中に限らない)についても適用可能である。
なお、各実施の形態において同一又は相当部分には同一の符号を付して説明する。
【0028】
本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態の構成について説明する。
図1は、本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態に関わる全体構成を示す図である。
位置計測装置4及び測位装置としてのGPS受信装置31を有する母船1と、航走体制御システム10を有する航走体2と、ケーブル3とを有するとともに、GPS用衛星群30を使用する。
【0029】
母船1は、航走体2を探査海域まで輸送し、探査海域において航走体2を送出し、航走体2の運行を制御する。位置計測システムである位置計測装置4及びGPS受信装置31を有する。
測位装置としてのGPS受信装置31は、母船1に属し、GPS用衛星群30を利用して、母船1の位置としての第1位置情報を出力する。第1位置情報は、(緯度、経度)で例示される。
位置計測装置4は、母船1に属し、母船1に対する航走体2の相対位置としての相対位置情報(例示:母船1からの距離、方向)を算出すると共に、航走体2の対地速度としての速度情報(例示:航行速度)と、航走体2の向かう方位としての方位情報(例示:経線を基準とする緯度方向及び水深方向の角度)とを取得する。そして、前述の第1位置情報と、相対位置情報と、速度情報と、方位情報とに基づいて、航走体2の位置としての第2位置情報を算出する。第2位置情報は、(緯度、経度、水深)で例示される。
【0030】
航走体2は、母船1で探査海域まで輸送され、探査海域において母船1から送出され、母船1に制御されながら水中を航行する。航走体2は、無人水中探査船で例示される。航走体制御システム10を有する。
航走体制御システム10は、航走体2に属する。母船1からの航走体2の相対位置を母船1において算出するために必要な探査信号(例示:ピンガー音)を母船1に向けて発信する。同時に、探査信号の発信時刻としての第1時刻情報を母船1へ出力する。また、航走体2の対地速度としての速度情報(例示:航行速度)と、航走体2の向かう方位としての方位情報(例示:緯線を基準とする経度方向及び水深方向の角度)を計測する。更に、航走体2の対水速度として対水速度情報(例示:海水に対する相対速度)を計測する。計測された速度情報、方位情報及び対水速度情報は、それぞれ必要に応じて、母船1へ出力される。
【0031】
ケーブル3は、母船1と航走体2とを電気的に接続する。航走体2(の航走体制御システム10)から出力される速度情報、方位情報及び対水速度情報を母船1(の位置計測装置4)へ伝達する。また、母船1から出力される航走体2を制御するための制御信号を航走体2へ伝達する。ケーブル3を用いず、無線を用いて行なっても良い。
【0032】
GPS用衛星群30は、一般に用いられるGPSの運用に使用される複数の衛星である。その衛星群からの情報により、GPS受信装置31は、母船1の第1位置情報を出力(算出)する。
【0033】
次に、位置計測装置4について、更に説明する。
図2は、本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態において適用される位置計測装置4の構成を示す図である。位置計測装置4は、制御部5、受波器6、記憶部7、送受信部8及び入出力部9を具備する。
【0034】
制御部5は、位置計測装置4に属する。位置計測装置4全般の制御及び航走体2全般の制御を行なう。また、探査信号の母船1への到着時刻としての第2時刻情報と、第1時刻情報とに基づいて、母船1からの航走体2の相対位置情報を算出する。そして、第1位置情報と、相対位置情報と、速度情報と、方位情報(必要に応じて対水速度情報)とに基づいて、第2位置情報を算出する。算出方法は、後述する。
送受信部8は、位置計測装置4に属する。航走体2から送信される速度情報と方位情報とを、ケーブル3を介して受信する。そして、航走体2から発せられる探査信号の発信時刻としての第1時刻情報を、ケーブル3を介して受信する。また、航走体2を制御するための信号を制御部5から受け取り、ケーブル3を介して航走体2へ出力する。ただし、ケーブル3を用いない場合には、海中無線装置が送受信部8となる。
受波器6は、位置計測装置4に属し、航走体2からの探査信号を受信する。そして、探査信号に基づいて、探査信号の発信方向と、探査信号の母船1への到着時刻としての第2時刻情報とを制御部5へ出力する。
記憶部7は、位置計測装置4に属し、第2位置情報を算出するために利用されるカルマンフィルタゲイン(後述)のテーブルが記憶されている。必要に応じて、既述の各種情報を保持する。
入出力部9は、位置計測装置4に属し、制御部5での制御に必要な情報の入出力を行なう。ケーブル3を伝送する情報もここを経由する。海中無線装置を用いる場合には、入出力部9はアンテナを含む。
【0035】
次に、航走体制御システム10について、更に説明する。
図3は、本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態において適用される航走体制御システム10の構成を示す図である。航走体制御システム10は、対地速度検出器11、記憶部12、対水速度検出器13、方位センサ14、送受信部15、船体制御部16及びピンガー17を具備する。
【0036】
対地速度検出器11は、航走体制御システム10に属し、航走体2の(対地)速度を計測し、速度情報として出力する。対地速度検出器11は、ドップラーソナーに例示される。
記憶部12は、計測した各種情報(速度情報、対水速度情報、方位情報)や航行スケジュール(例示:時刻と航走体2の位置を記したテーブル)、探査信号発信時刻などを記憶する。
対水速度検出器13は、航走体制御システム10に属し、航走体2の(対水)速度を計測し、対水速度情報として出力する。対水速度検出器13は、電磁ログ方式の検出器に例示される。
方位センサ14は、航走体制御システム10に属し、航走体2の進行する方位を計測し、方位情報として出力する。方位センサ14は、方位磁石に例示される。
送受信部15は、航走体制御システム10に属する。速度情報と方位情報とを、ケーブル3を介して母船1へ送信する。また、探査信号の発信時刻としての第1時刻情報を、ケーブル3を介して母船1へ送信する。更に、必要に応じて対水速度情報を、ケーブル3を介して母船1へ送信する。また、母船1から送信されてくる航走体2を制御するための信号(電力、制御指令など)を、航走体2へ伝達する。
船体制御部16は、航走体制御システム10に属する。母船1から送信されてくる航走体2を制御するための信号に基づいて、航走体2を制御する。また、航走体2から出力される信号・情報の入出力を管理する。
ピンガー17は、航走体制御システム10に属する。探査信号を母船1に向けて発信する。探査信号は、ピンガー音に例示される。また、探査信号発信と同時に、発信時刻を示す第1時刻情報を出力する。
【0037】
次に、本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態の動作について、添付図面を参照して説明する。
図5は、本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態の動作に関わる計算モデルに用いる座標系を示す図である。本例では、2次元平面での計算について示し、深さ(水深)方向に付いては省略している。しかし、通常の3次元への数学的拡張方法を用いれば、3次元空間における運動に応用することが可能である。
絶対座標系xyを設定し、母船1の座標を(X,Y)、航走体2の計測した座標を(xS,yS)、航走体2の(対地)速度を(ug,vg)、航走体2の方位をφとする。
【0038】
本発明においては、母船1にて第1位置情報(計測値)と相対位置情報(計測値)とから得られた航走体2位置(xS,yS)(計測値)に基づいて、カルマンフィルタ(Kalman Filter)を用いて、航走体2位置(x^,y^)を正確に推定する。この推定により、第1位置情報と相対位置情報とから得られる航走体2の位置に比較して、より精度が高く、ばらつきが少ない位置(第2位置情報)を得ることが可能となる。
ただし、文章中の右肩添字“^”は、推定値を示す(明細書中で同じ)。数式表現では真上添字で表現される場合も有る。
【0039】
次に、位置計算方法のアルゴリズムについて説明する。
位置計算のステップkにおけるアルゴリズムに用いる各式(9)〜(12)を以下に示す。
【数5】
ただし、
添え字k:k番目の演算時間に関する
Kk:カルマンフィルタゲイン
XS:航走体2位置(計測値):第1位置情報と相対位置情報とに基づいて求める航走体位置計測値
(3次元の場合には、(xS,yS,zS)である。)
X^:航走体2位置(計算値):第2位置情報(目的とする最終的な推定値)
(3次元の場合には、(x^,y^,z^)である。)
Xk^:演算時間(ΔT)ごとの航走体2位置(計算値)(推定値)
(3次元の場合には、(xk^,yk^,zk^)である。)
φ:航走体2方位(計測値):方位情報
A,B,C,D:航走体2の運動モデルの状態空間表現(後述の式(17)及び式(18)参照)
u:航走体2の運動モデルへの入力(絶対座標系での対地速度、計測値):速度情報
Q:システム雑音wの分散
R:測定誤差nの分散
N:そのシステム雑音wとその測定誤差nの共分散
である。
【0040】
本例では、2次元における航走体の運動について推定するので、以下の式(13)〜(16)を用いる。すなわち、
【数6】
ただし、
添え字k:k番目の演算時間に関する
Kk:カルマンフィルタゲイン
(xS,yS):航走体2位置(計測値):第1位置情報と相対位置情報とに基づいて求める
(x^,y^):航走体2位置(計算値):第2位置情報(目的とする最終的な推定値)
(xk^,yk^):演算時間(ΔT)ごとの航走体2位置(計算値)(推定値)
φ:航走体2方位(計測値):方位情報
A,B,C,D:航走体2の運動モデルの状態空間表現(後述の式(17)及び式(18)参照)
u:航走体2の運動モデルへの入力(絶対座標系での対地速度、計測値):速度情報(後述の式(17)参照)
Q:システム雑音wの分散(後述の式(19)参照)
R:測定誤差nの分散(式(19)参照)
N:システム雑音wと測定誤差nの共分散(式(18)参照)
【0041】
ここで、離散時間系での航走体2の運動推定モデルを次式(17),(18)で表す。
【数7】
ここで、Dは、0行列であり、式(19)は下記で示される。
【数8】
ただし、
ΔT:演算時間刻み
(ug,vg)=u:航走体2の(対地)速度(絶対座標系での対地速度、計測値):速度情報(前記式(17)]参照)
(uw,vw):対地速度計測誤差
(nx,ny):測位装置(GPS)及び音響測位による位置の計測誤差
q:対地速度計測誤差の分散
r:測位装置(GPS)及び音響測位による位置の計測誤差の分散
κ:カルマンフィルタゲイン(演算時間ΔT毎に変化)
【0042】
演算手順(I)について説明する。
(1)k=0である初期値P0が与えられる(k=1であるP1でも良い)。
同時に、計測値{(xS,yS):航走体2位置、φ:航走体2方位:方位情報、u:航走体2の運動モデルへの入力(絶対座標系での対地速度):速度情報}が与えられる。
(2)k=k+1として、前記式(13)及び(14)(漸化式)からPkが計算される(式(14)では、Pk+1と表示)。
(3)Pkに基づいて、式(13)からKkが計算される。
(4)Kkに基づいて、式(15)及び式(16)から(xk^,yk^)が計算(推定)される。
(5)|Pk−Pk−1|<δ(δ:予め設定された基準値){ただし、kΔT<TS(TS:予め設定された許容時間)}、の状況になった場合、計算を終了し、(6)に進む。そうでない場合には、(2)へ戻り、計算を繰り返す。
(6)航走体2の位置(x^,y^)は、(xk^,yk^)となり、そのときの時刻は、(直前の推定値算出時刻)+kΔT、で示される。
【0043】
ΔT毎の航走体2の位置(xk^,yk^)が演算(推定)され、前記式(15)の右辺第3項の因子(中括弧{ }内)に基づいて、ΔT毎に、測定誤差(nx,ny)を漸化的に除去する数学的フィルタリングが実行される。この数学的フィルタリングは、高精度な数学的推定方法である。
【0044】
ただし、上記演算方法(I)においては、計算が発散する{(5)の条件をkΔT<TSで満足できない}危険性を有している。従って、演算の発散の危険性を回避するためには、各ステップkでのカルマンフィルタのフィルタゲインKkの更新は行なわず、定常解を用いることとする。定常解は次式(20)で与えられる。
【数9】
ただし、Pは次のRiccati方程式の正定対称な解である。
【数10】
ただし、各符号は式(13)〜(19)での説明を参照のこと。
【0045】
演算手順(II)について説明する。
(1)計測値{(xS,yS):航走体2位置、φ:航走体2方位:方位情報、u:航走体2の運動モデルへの入力(絶対座標系での対地速度):速度情報}が与えられる。
(2)式(20)及び(21)からPが計算される。
(3)Pに基づいて、式(20)からKが計算される。
(4)Kに基づいて、式(13)及び(14)から(xk^,yk^)が計算(推定)される。
(5)航走体2の位置(x^,y^)は、(xk^,yk^)となり、そのときの時刻は、(直前の推定値算出時刻)+ΔT、で示される。
【0046】
上記演算方法(II)においては、上記の式(20)に基づいて、(2)、(3)でP及びKについて演算を行なっている。しかし、P及びKについて事前に演算しておき、実際の制御では、(1)、(4)、(5)のみを実行するようにすることも可能である。
【0047】
その場合、カルマンフィルタの定常解を予め演算して図4に示すようなテーブルを用意しておく。そして、位置計算のステップにおいては、テーブルから適切なフィルタゲインを選択し、それを用いて演算を行なう。そして、フィルタゲインの更新を行なわないこととする。それにより、カルマンフィルタの発散により航走体2の位置演算に失敗し、航走体2を見失う危険を回避することが出来る。
【0048】
図4は、フィルタゲインを有するのフィルタゲインテーブルを示す表である。
フィルタゲインテーブル20は、パターン21、q22、r23及びKk24とからなる。
パターン21は、気象や海象の状況に基づいて分類されたパターンである。あるいは、対地速度検出器11に異常が生じて対水速度検出器13に切り替えた場合、及び、GPS受信装置31による位置計測に異常が生じて第1位置情報を用いないで制御する場合に対応したパターンもある。航走体2が動作を行なう海域の気象や海象、各機器の状況に対応したパターンを選択する。
q22は、対地速度検出器11による対地速度計測の誤差の分散(上記q)である。パターン21に対応した値である。
r23は、GPS受信装置31による第1位置情報計測及びピンガー17−受波器6による相対位置情報計測の誤差の分散(上記r)である。パターン21に対応した値である。
Kk24は、カルマンフィルタにおけるフィルタゲインである。式(20)及び(21)においては、“K”と記されている。式(20)及び(21)に基づいて計算される。
【0049】
次に、本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態の動作の詳細について説明する。
i)母船1は、搭載していた航走体2を、探査海域において海中へ送出する。母船1と航走体2とは、ケーブル3を介した電気通信が可能である(海中無線通信でも良い)。
ii)母船1の位置計測装置4の制御部5は、入出力部9からの航走体2への移動指令(例示:緯度、経度、水深など)又は記憶部7の移動スケジュールによる移動指令に基づいて、その移動指令を送受信部8から船走体2の航走体制御システム10へケーブル3を介して送信する。
iii)航走体2は、送受信部15を介して移動指令を受信する。そして、その移動指令に基づいて移動を行なう。
iv)航走体2では、航走体制御システム10の船体制御部16からの指令により、対地速度検出器11が、対地速度を計測している。また、方位センサが、進行方向の方位を計測している。更に、対水速度検出器13が、海水に対する速度を計測している。それら計測された値は、記憶部12に記憶されている。
v)航走体2では、適当な時間毎に、航走体制御システム10の船体制御部16からの指令により、ピンガー17が、船外へピンガー音を発する。それと共に、ピンガー音を発した時刻(第1時刻)である第1時刻情報、第1時刻あるいは第1時刻に最も近い時刻での速度情報(対地速度)及び方位情報(進行方向の方位)を、送受信部15から母船1の位置計測装置4へケーブル3を介して送信する。
vi)母船1のGPS受信装置31は、GPS衛星群30からのGPS信号に基づいて、母船1の位置を算出する。そして、その位置を第1位置情報として位置計測装置4へ出力する。第1位置情報は、その時刻と共に記憶部12に保持される。
vi)母船1の位置計測装置4は、受波器6により、ピンガー音を受信する。そしてそのピンガー音に基づいて、受信した時刻(到着時刻)である第2時刻情報と、ピンガー音の発信方向を記憶部7に記憶する。
viii)母船1の位置計測装置4は、送受信部8にて、第1時刻情報、速度情報及び方位情報を受信する。
ix)母船1の位置計測装置4は、制御部5において、第1時刻情報と、第2時刻情報と、ピンガー音(の発信方向)とに基づいて、航走体2の母船に対する相対位置である相対位置情報を求める。すなわち、ピンガー音の音速×(第2時刻情報−第1時刻情報)=相対距離である。また、発信方向は、水平方向の角度(例示:緯線からの角度)及び垂直方向の角度(例示:水平面からの角度)である。
x)母船1の位置計測装置4は、制御部5において、第1位置情報と相対位置情報とに基づいて、航走体2の位置(xS,yS)(計測値)を求める。そして、その航走体位置(計測地)と、速度情報と、方位情報とに基づいて、航走体2の位置としての第2位置情報を推定する。推定方法は、既述の演算手順(I)あるいは(II)を用いる。
なお、演算手順(II)を用いる場合には、位置計測装置4の記憶部7に保持されているフィルタゲインテーブル20を用いる。
xi)推定された航走体2の第2位置情報を用いて、母船1の位置計測装置4の制御部5は、適宜移動指令を発し、その運動を制御する。
xii)適宜(ii){又は(iv)}〜(xi)を繰り返す。
【0050】
以上の動作により、カルマンフィルタを用いた新規の位置計測システム、位置計測方法及び当該方法をコンピューターに実行させるプログラムを使用して、航走体2の位置を正確に把握することが可能となる。
特に、短い時間毎に連続して位置推定を繰り返すことにより、常時正確な位置を把握しつづけることが出来る。
【0051】
また、GPS受信装置31による位置計測に異常が生じ、第1位置情報に問題が生じても、第1位置情報を用いないで、航走体2の位置を推定することが出来る。例えば、演算手順(II)における(xS,yS)について、それまでの第1位置情報及び相対位置情報とに基づいて計算してしていたものを、相対位置情報である航走体2の相対位置(xS’,yS’)とする。そして、式(15)において、(xS,yS)を(xS’,yS’)に置き換え、予めフィルタゲインテーブル(図4)に記憶された対応するカルマンフィルタゲインKを用いて、第2位置情報(x^,y^)を推定する。
すなわち、対地速度検出器11のみの演算に切り替えることにより、位置計算を継続することができる。すなわち、航走体2を見失うことが無い。
【0052】
また、対地速度検出器11に異常が生じ、速度情報に問題が生じても、対水速度検出器13に切り替えることにより対応することが可能である。例えば、母船1において海流の速度を計測しておき、海流速度と対水速度とから航走体の対地速度を求めることが出来る。そして、その対地速度を速度情報とすれば良い。
【0053】
あるいは、演算手順(II)における速度情報uについて、それまでの対地速度検出器11で測定された対地速度を用いていたものを、対水速度検出器13で測定された対水速度を用いることとする(速度情報u’)。そして、式(15)〜(17)において、uをu’に置き換え、予めフィルタゲインテーブル(図4)に記憶された対応するカルマンフィルタゲインKを用いて、第2位置情報(x^,y^)を推定する。
すなわち、対水速度検出器13による速度情報を用いることにより、位置計算を計測することが出来、位置計算を継続することができる。すなわち、航走体2を見失うことが無い。
【0054】
上記の位置計測システム及び位置計測方法を実施した実験例を以下に示す。ここでは、航走体2の実際の位置におけるyの値を固定し、xの位置を変化させて、位置を比較している。
図6は、航走体2の実際の位置(x0,y0):航走体位置(実際)、GPS受信装置31による第1位置情報及びピンガー17−受波器6による相対位置情報とから得られる航走体の位置(xS,yS):航走体位置(計測)、カルマンフィルタを用いて計算した航走体2の位置(x^,y^):航走体位置(計算)のそれぞれの航走体2の位置を比較した図である。ただし、縦軸及び横軸は、図5で示した軸である。
図から明らかなように、航走体位置(計算)は、航走体位置(実際)とほぼ同じ位置を示している。そして、航走体位置(計測)に比較して、誤差及びばらつきがほとんど発生していないことがわかる。
【0055】
図7は、母船1の実際の位置(X0,Y0):母船位置(実際)、母船1のGPS受信装置31による第1位置情報:母船位置(計測)のそれぞれの母船1の位置を比較した図である。ただし、縦軸及び横軸は、図6と同様である。また、母船1の実際の位置におけるYの値を固定し、Xの位置を変化させて、各位置を比較している。
図7から明らかなように、母船位置(計測)は、母船位置(実際)を常に正確に表しているわけではなく、ある程度の誤差とばらつきを有することがわかる。
【0056】
そして、図6及び図7から明らかなように、図7のように母船位置が必ずしも正確に把握されていない状態にもかかわらず、図6のような正確な航走体2の位置を把握することが可能となる。
【0057】
【発明の効果】
本発明により、母船から離れて航行する航走体の位置を、誤差やばらつきが少なく、正確に把握することが可能となる。
【図面の簡単な説明】
【図1】本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態に関わる全体構成を示す説明図である。
【図2】前記実施の形態において適用される位置計測装置の構成を示すブロック図である。
【図3】前記実施の形態において適用される航走体制御システムの構成を示すブロック図である。
【図4】前記実施の形態において使用されるカルマンフィルタゲインを示す表である。
【図5】前記実施の形態の動作に関わる計算モデルに用いる座標系を示す説明図である。
【図6】航走体の実際の位置、航走体の計測による位置及び航走体の推測による位置を比較した説明図である。
【図7】母船の実際の位置、母船の計測による位置を比較した説明図である。
【符号の説明】
1 母船
2 航走体
3 ケーブル
4 位置計測装置
5 制御部
6 受波器
7 記憶部
8 送受信部
9 入出力部
10 航走体制御システム
11 対地速度検出器
12 記憶部
13 対水速度検出器
14 方位センサ
15 送受信部
16 船体制御部
17 ピンガー
20 フィルタゲインテーブル
21 パターン
22 q
23 r
24 Kk
【発明の属する技術分野】
本発明は、位置計測システム及び位置計測方法等に係り、特に、母船と母船に制御されながら航行する航走体とを有する系における位置計測システム、位置計測方法、当該方法をコンピューターに実行させるプログラム及び航走体運用システムに関する。
【0002】
【従来の技術】
水中で行なう各種探査においては、基地となる母船により制御された無人の航走体が、実際の探査作業を行なうことが多い。その場合、母船は、航走体の正確な絶対位置(以下、単に「位置」ともいう)を把握している必要がある。
【0003】
そのような位置把握を行なう方法として、以下のような方法が知られている。すなわち、まずGPS(Global Positioning System:衛星航法システム)のような測位方法を用いて母船の位置を把握する。同時に、航走体のピンガー(海中の物標や地形の状況などの探査のための音波発信器)のような水中音響装置が発する音響(ピンガー音など)の受信状態により、航走体の母船に対する相対位置を把握する。そして、母船の位置と母船に対する相対位置とに基づいて、航走体の(絶対)位置を把握する方法である。
【0004】
上記方法では、航走体の位置計測の誤差が大きいという問題があった。例えば、水中音響装置を用いた場合には、海中の状態により突然に遠くの場所に位置計算結果が飛んでしまうことが起こる。それに対応するために、位置計算結果が飛ばないように、現在位置から一定距離内への位置更新のみを認めるロジックを入れる。そうすることで、位置計測結果が遠くの場所に飛ぶことをある程度回避できる。
【0005】
しかし、初期位置が正しく計測されていないと、ずっと正しくない位置を出力しつづけることになる。また、ピンガーのような水中音響による母船と航走体間の位置計測が途絶する(航走体からの信号を母船が受信できない)と、航走体の位置が求められない。
【0006】
このため、航走体の位置を、誤差やばらつきを少なく、正確に把握することが可能な技術が求められている。また、母船の測位装置に異常が生じた場合や、水中音響による相対位置の計測が困難な状況でも航走体の位置を正確に把握することが可能な技術が求められている。また、水中のようなGPS信号が届かない場所にいる航走体の位置を正確に把握する技術が求められている。
【0007】
【発明が解決しようとする課題】
従って、本発明の目的は、母船から離れて航行する航走体の位置を正確に把握可能な位置計測システム、位置計測方法、当該方法をコンピューターに実行させるプログラム及び航走体運用システムを提供することである。
【0008】
また、本発明の他の目的は、母船に搭載された測位装置に問題が生じても航走体の位置を正確に把握可能な位置計測システム、位置計測方法、当該方法をコンピューターに実行させるプログラム及び航走体運用システムを提供することである。
【0009】
本発明の更に他の目的は、水中音響による航走体の相対位置の計測が困難な場合でも航走体の位置を正確に把握可能な位置計測システム、位置計測方法、当該方法をコンピューターに実行させるプログラム及び航走体運用システムを提供することである。
【0010】
本発明の別の目的は、GPS信号の届かない水中においても、航走体の位置を正確に求めることが可能な位置計測システム、位置計測方法、当該方法をコンピューターに実行させるプログラム及び航走体運用システムを提供することである。
【0011】
【課題を解決するための手段】
以下に、[発明の実施の形態]で使用される番号・符号を用いて、課題を解決するための手段を説明する。これらの番号・符号は、[特許請求の範囲]の記載と[発明の実施の形態]との対応関係を明らかにするために付加されたものである。ただし、それらの番号・符号を、[特許請求の範囲]に記載されている発明の技術的範囲の解釈に用いてはならない。
【0012】
従って、上記課題を解決するために、本発明の位置計測システムは、母船(1)の位置としての第1位置情報を出力する測位装置(31)と、その第1位置情報と、母船(1)から送出された航走体(2)の母船(1)に対する相対位置としての相対位置情報と、航走体(2)の対地速度としての速度情報と、航走体(2)の向かう方位としての方位情報とに基づいて、航走体(2)の位置としての第2位置情報を推定する位置計測装置(4)とを具備する。
【0013】
また、本発明の位置計測システムは、位置計測装置(4)が、航走体(2)から送信されるその速度情報とその方位情報とを受信する送受信部(8)を具備する。
【0014】
また、本発明の位置計測システムは、位置計測装置(4)が、航走体(2)から発せられる探査信号を受信し、その探査信号の到着時刻としての第2時刻情報と、その探査信号の母船(1)に対する方向とを出力する受波器(6)と、航走体(2)から発せられるその探査信号の発信時刻としての第1時刻情報と、その第2時刻情報と、その方向とに基づいて、その相対位置情報を算出する制御部(5)とを更に具備し、送受信部(8)は、更に、その第1時刻情報を受信する。
【0015】
また、本発明の位置計測システムは、制御部(5)が、更に、その第1位置情報と、その相対位置情報と、その速度情報と、その方位情報とに基づいて、その第2位置情報を推定する。
【0016】
更に、本発明の位置計測システムは、制御部(5)が、カルマンフィルタを用いてその第2位置情報を推定する。
【0017】
更に、本発明の位置計測システムは、制御部(5)が、予め設定された演算時間ΔT毎のその第2位置情報を下記式(1)乃至式(3)に基づいて計算し、下記式(4)によりその第2位置情報を推定する。ここで、
【数3】
ただし、添え字k:k番目の演算時間、Kk:カルマンフィルタゲイン、XS:第1位置情報と相対位置情報とに基づいて求める航走体位置計測値、X^:第2位置情報、Xk^:演算時間ΔT毎の第2位置情報、φ:方位情報、A,B,C,D:航走体の運動モデルの状態空間表現、u:速度情報、Q:システム雑音wの分散、R:測定誤差nの分散、N:システム雑音wと測定誤差nの共分散である。
【0018】
更に、本発明の位置計測システムは、制御部(5)が、そのΔT毎のその第2位置情報を、前記式(1)及び式(2)を用いず、前記式(3)と予め設定されたカルマンフィルタゲインKとを用いて計算し、計算されたそのΔT毎のその第2位置情報に基づいて、前記式(4)でその第2位置情報を推定する。
【0019】
更に、本発明の位置計測システムは、母船(1)から送出された航走体(2)の対地速度としての速度情報と、航走体(2)の向かう方位としての方位情報とを航走体(2)から受信する送受信部(8)と、航走体(2)から発せられる探査信号を受信し、その探査信号の到着時刻としての第2時刻情報と、その探査信号の母船(1)に対する方向とを出力する受波器(6)と、航走体(2)から発せられる探査信号の発信時刻としての第1時刻情報と、その第2時刻情報と、その方向とに基づいて、その相対位置情報を算出する制御部(5)とを有する位置計測装置(4)を具備し、送受信部(8)は、更に、その第1時刻情報を受信し、制御部(5)は、更に、その相対位置情報と、その速度情報と、その方位情報とに基づいて、航走体(2)の位置としての第2位置情報を推定する。
【0020】
更に、本発明の位置計測システムは、制御部(5)が、カルマンフィルタを用いてその第2位置情報を推定する。
【0021】
更に、本発明の位置計測システムは、制御部(5)が、予め設定された演算時間ΔT毎のその第2位置情報を下記式(5)乃至式(7)に基づいて計算し、下記式(8)によりその第2位置情報を推定する。ここで、
【数4】
ただし、添え字k:k番目の演算時間、Kk:カルマンフィルタゲイン、XS:相対位置情報に基づいて求める航走体位置計測値、X^:第2位置情報、Xk^:演算時間ΔT毎の第2位置情報、φ:方位情報、A,B,C,D:航走体の運動モデルの状態空間表現、u:速度情報、Q:システム雑音wの分散、R:測定誤差nの分散、N:システム雑音wと測定誤差nの共分散である。
【0022】
更に、本発明の位置計測システムは、その速度情報が、航走体(2)の対地速度検出器(11)及び対水速度検出器(13)の少なくとも一方の計測結果に基づいている。
【0023】
更に、本発明の位置計測システムは、母船(1)は船舶であり、航走体(2)は水中を探査する潜水艇である。
【0024】
上記課題を解決するために、本発明の航走体運用システムは、上記各項のいずれか一項に記載の位置計測システムと、その位置計測システムを有する母船(1)と、母船(1)に制御され運動する航走体(2)とを具備する。
【0025】
上記課題を解決するために、本発明の位置計測方法は、母船(1)の位置としての第1位置情報を取得するステップと、母船(1)から送出された航走体(2)の対地速度としての速度情報と、航走体(2)の向かう方位としての方位情報と、航走体(2)から発せられるその探査信号の発信時刻としての第1時刻情報とを航走体(2)から取得するステップと、その探査信号に基づいて、その探査信号の到着時刻としての第2時刻情報と、その探査信号の母船(1)に対する方向とを取得するステップと、その第1時刻情報と、その第2時刻情報と、その方向とに基づいて、航走体(2)の母船(1)に対する位置である相対位置情報を算出するステップと、その第1位置情報と、その相対位置情報と、その速度情報と、その方位情報とに基づいて、カルマンフィルタを用いて航走体(2)の位置としての第2位置情報を推定するステップとを具備する。
【0026】
上記課題を解決するために、本発明に関わるプログラムは、母船(1)の位置としての第1位置情報を取得するステップと、母船(1)から送出された航走体(2)の対地速度としての速度情報と、航走体(2)の向かう方位としての方位情報と、航走体(2)から発せられるその探査信号の発信時刻としての第1時刻情報とを航走体(2)から取得するステップと、その探査信号の到着時刻としての第2時刻情報と、その探査信号の母船(1)に対する方向とを取得するステップと、その第1時刻情報と、その第2時刻情報と、その方向とに基づいて、航走体(2)の母船(1)に対する位置である相対位置情報を算出するステップと、その第1位置情報と、その相対位置情報と、その速度情報と、その方位情報とに基づいて、カルマンフィルタを用いて航走体(2)の位置としての第2位置情報を推定するステップとをコンピューターに実行させる。
【0027】
【発明の実施の形態】
以下、本発明である位置計測システム、位置計測方法、当該方法をコンピューターに実行させるプログラム、及び航走体運用システムの実施の形態に関して、添付図面を参照して説明する。
本実施の形態において、母船に制御され、水中探査に使用される航走体を例に示して説明する。ただし、他の母船に制御され、母船を離れて動作を行なう航走体(水中に限らない)についても適用可能である。
なお、各実施の形態において同一又は相当部分には同一の符号を付して説明する。
【0028】
本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態の構成について説明する。
図1は、本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態に関わる全体構成を示す図である。
位置計測装置4及び測位装置としてのGPS受信装置31を有する母船1と、航走体制御システム10を有する航走体2と、ケーブル3とを有するとともに、GPS用衛星群30を使用する。
【0029】
母船1は、航走体2を探査海域まで輸送し、探査海域において航走体2を送出し、航走体2の運行を制御する。位置計測システムである位置計測装置4及びGPS受信装置31を有する。
測位装置としてのGPS受信装置31は、母船1に属し、GPS用衛星群30を利用して、母船1の位置としての第1位置情報を出力する。第1位置情報は、(緯度、経度)で例示される。
位置計測装置4は、母船1に属し、母船1に対する航走体2の相対位置としての相対位置情報(例示:母船1からの距離、方向)を算出すると共に、航走体2の対地速度としての速度情報(例示:航行速度)と、航走体2の向かう方位としての方位情報(例示:経線を基準とする緯度方向及び水深方向の角度)とを取得する。そして、前述の第1位置情報と、相対位置情報と、速度情報と、方位情報とに基づいて、航走体2の位置としての第2位置情報を算出する。第2位置情報は、(緯度、経度、水深)で例示される。
【0030】
航走体2は、母船1で探査海域まで輸送され、探査海域において母船1から送出され、母船1に制御されながら水中を航行する。航走体2は、無人水中探査船で例示される。航走体制御システム10を有する。
航走体制御システム10は、航走体2に属する。母船1からの航走体2の相対位置を母船1において算出するために必要な探査信号(例示:ピンガー音)を母船1に向けて発信する。同時に、探査信号の発信時刻としての第1時刻情報を母船1へ出力する。また、航走体2の対地速度としての速度情報(例示:航行速度)と、航走体2の向かう方位としての方位情報(例示:緯線を基準とする経度方向及び水深方向の角度)を計測する。更に、航走体2の対水速度として対水速度情報(例示:海水に対する相対速度)を計測する。計測された速度情報、方位情報及び対水速度情報は、それぞれ必要に応じて、母船1へ出力される。
【0031】
ケーブル3は、母船1と航走体2とを電気的に接続する。航走体2(の航走体制御システム10)から出力される速度情報、方位情報及び対水速度情報を母船1(の位置計測装置4)へ伝達する。また、母船1から出力される航走体2を制御するための制御信号を航走体2へ伝達する。ケーブル3を用いず、無線を用いて行なっても良い。
【0032】
GPS用衛星群30は、一般に用いられるGPSの運用に使用される複数の衛星である。その衛星群からの情報により、GPS受信装置31は、母船1の第1位置情報を出力(算出)する。
【0033】
次に、位置計測装置4について、更に説明する。
図2は、本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態において適用される位置計測装置4の構成を示す図である。位置計測装置4は、制御部5、受波器6、記憶部7、送受信部8及び入出力部9を具備する。
【0034】
制御部5は、位置計測装置4に属する。位置計測装置4全般の制御及び航走体2全般の制御を行なう。また、探査信号の母船1への到着時刻としての第2時刻情報と、第1時刻情報とに基づいて、母船1からの航走体2の相対位置情報を算出する。そして、第1位置情報と、相対位置情報と、速度情報と、方位情報(必要に応じて対水速度情報)とに基づいて、第2位置情報を算出する。算出方法は、後述する。
送受信部8は、位置計測装置4に属する。航走体2から送信される速度情報と方位情報とを、ケーブル3を介して受信する。そして、航走体2から発せられる探査信号の発信時刻としての第1時刻情報を、ケーブル3を介して受信する。また、航走体2を制御するための信号を制御部5から受け取り、ケーブル3を介して航走体2へ出力する。ただし、ケーブル3を用いない場合には、海中無線装置が送受信部8となる。
受波器6は、位置計測装置4に属し、航走体2からの探査信号を受信する。そして、探査信号に基づいて、探査信号の発信方向と、探査信号の母船1への到着時刻としての第2時刻情報とを制御部5へ出力する。
記憶部7は、位置計測装置4に属し、第2位置情報を算出するために利用されるカルマンフィルタゲイン(後述)のテーブルが記憶されている。必要に応じて、既述の各種情報を保持する。
入出力部9は、位置計測装置4に属し、制御部5での制御に必要な情報の入出力を行なう。ケーブル3を伝送する情報もここを経由する。海中無線装置を用いる場合には、入出力部9はアンテナを含む。
【0035】
次に、航走体制御システム10について、更に説明する。
図3は、本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態において適用される航走体制御システム10の構成を示す図である。航走体制御システム10は、対地速度検出器11、記憶部12、対水速度検出器13、方位センサ14、送受信部15、船体制御部16及びピンガー17を具備する。
【0036】
対地速度検出器11は、航走体制御システム10に属し、航走体2の(対地)速度を計測し、速度情報として出力する。対地速度検出器11は、ドップラーソナーに例示される。
記憶部12は、計測した各種情報(速度情報、対水速度情報、方位情報)や航行スケジュール(例示:時刻と航走体2の位置を記したテーブル)、探査信号発信時刻などを記憶する。
対水速度検出器13は、航走体制御システム10に属し、航走体2の(対水)速度を計測し、対水速度情報として出力する。対水速度検出器13は、電磁ログ方式の検出器に例示される。
方位センサ14は、航走体制御システム10に属し、航走体2の進行する方位を計測し、方位情報として出力する。方位センサ14は、方位磁石に例示される。
送受信部15は、航走体制御システム10に属する。速度情報と方位情報とを、ケーブル3を介して母船1へ送信する。また、探査信号の発信時刻としての第1時刻情報を、ケーブル3を介して母船1へ送信する。更に、必要に応じて対水速度情報を、ケーブル3を介して母船1へ送信する。また、母船1から送信されてくる航走体2を制御するための信号(電力、制御指令など)を、航走体2へ伝達する。
船体制御部16は、航走体制御システム10に属する。母船1から送信されてくる航走体2を制御するための信号に基づいて、航走体2を制御する。また、航走体2から出力される信号・情報の入出力を管理する。
ピンガー17は、航走体制御システム10に属する。探査信号を母船1に向けて発信する。探査信号は、ピンガー音に例示される。また、探査信号発信と同時に、発信時刻を示す第1時刻情報を出力する。
【0037】
次に、本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態の動作について、添付図面を参照して説明する。
図5は、本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態の動作に関わる計算モデルに用いる座標系を示す図である。本例では、2次元平面での計算について示し、深さ(水深)方向に付いては省略している。しかし、通常の3次元への数学的拡張方法を用いれば、3次元空間における運動に応用することが可能である。
絶対座標系xyを設定し、母船1の座標を(X,Y)、航走体2の計測した座標を(xS,yS)、航走体2の(対地)速度を(ug,vg)、航走体2の方位をφとする。
【0038】
本発明においては、母船1にて第1位置情報(計測値)と相対位置情報(計測値)とから得られた航走体2位置(xS,yS)(計測値)に基づいて、カルマンフィルタ(Kalman Filter)を用いて、航走体2位置(x^,y^)を正確に推定する。この推定により、第1位置情報と相対位置情報とから得られる航走体2の位置に比較して、より精度が高く、ばらつきが少ない位置(第2位置情報)を得ることが可能となる。
ただし、文章中の右肩添字“^”は、推定値を示す(明細書中で同じ)。数式表現では真上添字で表現される場合も有る。
【0039】
次に、位置計算方法のアルゴリズムについて説明する。
位置計算のステップkにおけるアルゴリズムに用いる各式(9)〜(12)を以下に示す。
【数5】
ただし、
添え字k:k番目の演算時間に関する
Kk:カルマンフィルタゲイン
XS:航走体2位置(計測値):第1位置情報と相対位置情報とに基づいて求める航走体位置計測値
(3次元の場合には、(xS,yS,zS)である。)
X^:航走体2位置(計算値):第2位置情報(目的とする最終的な推定値)
(3次元の場合には、(x^,y^,z^)である。)
Xk^:演算時間(ΔT)ごとの航走体2位置(計算値)(推定値)
(3次元の場合には、(xk^,yk^,zk^)である。)
φ:航走体2方位(計測値):方位情報
A,B,C,D:航走体2の運動モデルの状態空間表現(後述の式(17)及び式(18)参照)
u:航走体2の運動モデルへの入力(絶対座標系での対地速度、計測値):速度情報
Q:システム雑音wの分散
R:測定誤差nの分散
N:そのシステム雑音wとその測定誤差nの共分散
である。
【0040】
本例では、2次元における航走体の運動について推定するので、以下の式(13)〜(16)を用いる。すなわち、
【数6】
ただし、
添え字k:k番目の演算時間に関する
Kk:カルマンフィルタゲイン
(xS,yS):航走体2位置(計測値):第1位置情報と相対位置情報とに基づいて求める
(x^,y^):航走体2位置(計算値):第2位置情報(目的とする最終的な推定値)
(xk^,yk^):演算時間(ΔT)ごとの航走体2位置(計算値)(推定値)
φ:航走体2方位(計測値):方位情報
A,B,C,D:航走体2の運動モデルの状態空間表現(後述の式(17)及び式(18)参照)
u:航走体2の運動モデルへの入力(絶対座標系での対地速度、計測値):速度情報(後述の式(17)参照)
Q:システム雑音wの分散(後述の式(19)参照)
R:測定誤差nの分散(式(19)参照)
N:システム雑音wと測定誤差nの共分散(式(18)参照)
【0041】
ここで、離散時間系での航走体2の運動推定モデルを次式(17),(18)で表す。
【数7】
ここで、Dは、0行列であり、式(19)は下記で示される。
【数8】
ただし、
ΔT:演算時間刻み
(ug,vg)=u:航走体2の(対地)速度(絶対座標系での対地速度、計測値):速度情報(前記式(17)]参照)
(uw,vw):対地速度計測誤差
(nx,ny):測位装置(GPS)及び音響測位による位置の計測誤差
q:対地速度計測誤差の分散
r:測位装置(GPS)及び音響測位による位置の計測誤差の分散
κ:カルマンフィルタゲイン(演算時間ΔT毎に変化)
【0042】
演算手順(I)について説明する。
(1)k=0である初期値P0が与えられる(k=1であるP1でも良い)。
同時に、計測値{(xS,yS):航走体2位置、φ:航走体2方位:方位情報、u:航走体2の運動モデルへの入力(絶対座標系での対地速度):速度情報}が与えられる。
(2)k=k+1として、前記式(13)及び(14)(漸化式)からPkが計算される(式(14)では、Pk+1と表示)。
(3)Pkに基づいて、式(13)からKkが計算される。
(4)Kkに基づいて、式(15)及び式(16)から(xk^,yk^)が計算(推定)される。
(5)|Pk−Pk−1|<δ(δ:予め設定された基準値){ただし、kΔT<TS(TS:予め設定された許容時間)}、の状況になった場合、計算を終了し、(6)に進む。そうでない場合には、(2)へ戻り、計算を繰り返す。
(6)航走体2の位置(x^,y^)は、(xk^,yk^)となり、そのときの時刻は、(直前の推定値算出時刻)+kΔT、で示される。
【0043】
ΔT毎の航走体2の位置(xk^,yk^)が演算(推定)され、前記式(15)の右辺第3項の因子(中括弧{ }内)に基づいて、ΔT毎に、測定誤差(nx,ny)を漸化的に除去する数学的フィルタリングが実行される。この数学的フィルタリングは、高精度な数学的推定方法である。
【0044】
ただし、上記演算方法(I)においては、計算が発散する{(5)の条件をkΔT<TSで満足できない}危険性を有している。従って、演算の発散の危険性を回避するためには、各ステップkでのカルマンフィルタのフィルタゲインKkの更新は行なわず、定常解を用いることとする。定常解は次式(20)で与えられる。
【数9】
ただし、Pは次のRiccati方程式の正定対称な解である。
【数10】
ただし、各符号は式(13)〜(19)での説明を参照のこと。
【0045】
演算手順(II)について説明する。
(1)計測値{(xS,yS):航走体2位置、φ:航走体2方位:方位情報、u:航走体2の運動モデルへの入力(絶対座標系での対地速度):速度情報}が与えられる。
(2)式(20)及び(21)からPが計算される。
(3)Pに基づいて、式(20)からKが計算される。
(4)Kに基づいて、式(13)及び(14)から(xk^,yk^)が計算(推定)される。
(5)航走体2の位置(x^,y^)は、(xk^,yk^)となり、そのときの時刻は、(直前の推定値算出時刻)+ΔT、で示される。
【0046】
上記演算方法(II)においては、上記の式(20)に基づいて、(2)、(3)でP及びKについて演算を行なっている。しかし、P及びKについて事前に演算しておき、実際の制御では、(1)、(4)、(5)のみを実行するようにすることも可能である。
【0047】
その場合、カルマンフィルタの定常解を予め演算して図4に示すようなテーブルを用意しておく。そして、位置計算のステップにおいては、テーブルから適切なフィルタゲインを選択し、それを用いて演算を行なう。そして、フィルタゲインの更新を行なわないこととする。それにより、カルマンフィルタの発散により航走体2の位置演算に失敗し、航走体2を見失う危険を回避することが出来る。
【0048】
図4は、フィルタゲインを有するのフィルタゲインテーブルを示す表である。
フィルタゲインテーブル20は、パターン21、q22、r23及びKk24とからなる。
パターン21は、気象や海象の状況に基づいて分類されたパターンである。あるいは、対地速度検出器11に異常が生じて対水速度検出器13に切り替えた場合、及び、GPS受信装置31による位置計測に異常が生じて第1位置情報を用いないで制御する場合に対応したパターンもある。航走体2が動作を行なう海域の気象や海象、各機器の状況に対応したパターンを選択する。
q22は、対地速度検出器11による対地速度計測の誤差の分散(上記q)である。パターン21に対応した値である。
r23は、GPS受信装置31による第1位置情報計測及びピンガー17−受波器6による相対位置情報計測の誤差の分散(上記r)である。パターン21に対応した値である。
Kk24は、カルマンフィルタにおけるフィルタゲインである。式(20)及び(21)においては、“K”と記されている。式(20)及び(21)に基づいて計算される。
【0049】
次に、本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態の動作の詳細について説明する。
i)母船1は、搭載していた航走体2を、探査海域において海中へ送出する。母船1と航走体2とは、ケーブル3を介した電気通信が可能である(海中無線通信でも良い)。
ii)母船1の位置計測装置4の制御部5は、入出力部9からの航走体2への移動指令(例示:緯度、経度、水深など)又は記憶部7の移動スケジュールによる移動指令に基づいて、その移動指令を送受信部8から船走体2の航走体制御システム10へケーブル3を介して送信する。
iii)航走体2は、送受信部15を介して移動指令を受信する。そして、その移動指令に基づいて移動を行なう。
iv)航走体2では、航走体制御システム10の船体制御部16からの指令により、対地速度検出器11が、対地速度を計測している。また、方位センサが、進行方向の方位を計測している。更に、対水速度検出器13が、海水に対する速度を計測している。それら計測された値は、記憶部12に記憶されている。
v)航走体2では、適当な時間毎に、航走体制御システム10の船体制御部16からの指令により、ピンガー17が、船外へピンガー音を発する。それと共に、ピンガー音を発した時刻(第1時刻)である第1時刻情報、第1時刻あるいは第1時刻に最も近い時刻での速度情報(対地速度)及び方位情報(進行方向の方位)を、送受信部15から母船1の位置計測装置4へケーブル3を介して送信する。
vi)母船1のGPS受信装置31は、GPS衛星群30からのGPS信号に基づいて、母船1の位置を算出する。そして、その位置を第1位置情報として位置計測装置4へ出力する。第1位置情報は、その時刻と共に記憶部12に保持される。
vi)母船1の位置計測装置4は、受波器6により、ピンガー音を受信する。そしてそのピンガー音に基づいて、受信した時刻(到着時刻)である第2時刻情報と、ピンガー音の発信方向を記憶部7に記憶する。
viii)母船1の位置計測装置4は、送受信部8にて、第1時刻情報、速度情報及び方位情報を受信する。
ix)母船1の位置計測装置4は、制御部5において、第1時刻情報と、第2時刻情報と、ピンガー音(の発信方向)とに基づいて、航走体2の母船に対する相対位置である相対位置情報を求める。すなわち、ピンガー音の音速×(第2時刻情報−第1時刻情報)=相対距離である。また、発信方向は、水平方向の角度(例示:緯線からの角度)及び垂直方向の角度(例示:水平面からの角度)である。
x)母船1の位置計測装置4は、制御部5において、第1位置情報と相対位置情報とに基づいて、航走体2の位置(xS,yS)(計測値)を求める。そして、その航走体位置(計測地)と、速度情報と、方位情報とに基づいて、航走体2の位置としての第2位置情報を推定する。推定方法は、既述の演算手順(I)あるいは(II)を用いる。
なお、演算手順(II)を用いる場合には、位置計測装置4の記憶部7に保持されているフィルタゲインテーブル20を用いる。
xi)推定された航走体2の第2位置情報を用いて、母船1の位置計測装置4の制御部5は、適宜移動指令を発し、その運動を制御する。
xii)適宜(ii){又は(iv)}〜(xi)を繰り返す。
【0050】
以上の動作により、カルマンフィルタを用いた新規の位置計測システム、位置計測方法及び当該方法をコンピューターに実行させるプログラムを使用して、航走体2の位置を正確に把握することが可能となる。
特に、短い時間毎に連続して位置推定を繰り返すことにより、常時正確な位置を把握しつづけることが出来る。
【0051】
また、GPS受信装置31による位置計測に異常が生じ、第1位置情報に問題が生じても、第1位置情報を用いないで、航走体2の位置を推定することが出来る。例えば、演算手順(II)における(xS,yS)について、それまでの第1位置情報及び相対位置情報とに基づいて計算してしていたものを、相対位置情報である航走体2の相対位置(xS’,yS’)とする。そして、式(15)において、(xS,yS)を(xS’,yS’)に置き換え、予めフィルタゲインテーブル(図4)に記憶された対応するカルマンフィルタゲインKを用いて、第2位置情報(x^,y^)を推定する。
すなわち、対地速度検出器11のみの演算に切り替えることにより、位置計算を継続することができる。すなわち、航走体2を見失うことが無い。
【0052】
また、対地速度検出器11に異常が生じ、速度情報に問題が生じても、対水速度検出器13に切り替えることにより対応することが可能である。例えば、母船1において海流の速度を計測しておき、海流速度と対水速度とから航走体の対地速度を求めることが出来る。そして、その対地速度を速度情報とすれば良い。
【0053】
あるいは、演算手順(II)における速度情報uについて、それまでの対地速度検出器11で測定された対地速度を用いていたものを、対水速度検出器13で測定された対水速度を用いることとする(速度情報u’)。そして、式(15)〜(17)において、uをu’に置き換え、予めフィルタゲインテーブル(図4)に記憶された対応するカルマンフィルタゲインKを用いて、第2位置情報(x^,y^)を推定する。
すなわち、対水速度検出器13による速度情報を用いることにより、位置計算を計測することが出来、位置計算を継続することができる。すなわち、航走体2を見失うことが無い。
【0054】
上記の位置計測システム及び位置計測方法を実施した実験例を以下に示す。ここでは、航走体2の実際の位置におけるyの値を固定し、xの位置を変化させて、位置を比較している。
図6は、航走体2の実際の位置(x0,y0):航走体位置(実際)、GPS受信装置31による第1位置情報及びピンガー17−受波器6による相対位置情報とから得られる航走体の位置(xS,yS):航走体位置(計測)、カルマンフィルタを用いて計算した航走体2の位置(x^,y^):航走体位置(計算)のそれぞれの航走体2の位置を比較した図である。ただし、縦軸及び横軸は、図5で示した軸である。
図から明らかなように、航走体位置(計算)は、航走体位置(実際)とほぼ同じ位置を示している。そして、航走体位置(計測)に比較して、誤差及びばらつきがほとんど発生していないことがわかる。
【0055】
図7は、母船1の実際の位置(X0,Y0):母船位置(実際)、母船1のGPS受信装置31による第1位置情報:母船位置(計測)のそれぞれの母船1の位置を比較した図である。ただし、縦軸及び横軸は、図6と同様である。また、母船1の実際の位置におけるYの値を固定し、Xの位置を変化させて、各位置を比較している。
図7から明らかなように、母船位置(計測)は、母船位置(実際)を常に正確に表しているわけではなく、ある程度の誤差とばらつきを有することがわかる。
【0056】
そして、図6及び図7から明らかなように、図7のように母船位置が必ずしも正確に把握されていない状態にもかかわらず、図6のような正確な航走体2の位置を把握することが可能となる。
【0057】
【発明の効果】
本発明により、母船から離れて航行する航走体の位置を、誤差やばらつきが少なく、正確に把握することが可能となる。
【図面の簡単な説明】
【図1】本発明である位置計測システム、位置計測方法、プログラム及び航走体運用システムの実施の形態に関わる全体構成を示す説明図である。
【図2】前記実施の形態において適用される位置計測装置の構成を示すブロック図である。
【図3】前記実施の形態において適用される航走体制御システムの構成を示すブロック図である。
【図4】前記実施の形態において使用されるカルマンフィルタゲインを示す表である。
【図5】前記実施の形態の動作に関わる計算モデルに用いる座標系を示す説明図である。
【図6】航走体の実際の位置、航走体の計測による位置及び航走体の推測による位置を比較した説明図である。
【図7】母船の実際の位置、母船の計測による位置を比較した説明図である。
【符号の説明】
1 母船
2 航走体
3 ケーブル
4 位置計測装置
5 制御部
6 受波器
7 記憶部
8 送受信部
9 入出力部
10 航走体制御システム
11 対地速度検出器
12 記憶部
13 対水速度検出器
14 方位センサ
15 送受信部
16 船体制御部
17 ピンガー
20 フィルタゲインテーブル
21 パターン
22 q
23 r
24 Kk
Claims (14)
- 母船の位置としての第1位置情報を出力する測位装置と、
前記第1位置情報と、前記母船から送出された航走体の前記母船に対する相対位置としての相対位置情報と、前記航走体の対地速度としての速度情報と、前記航走体の向かう方位としての方位情報とに基づいて、前記航走体の位置としての第2位置情報を推定する位置計測装置とを具備することを特徴とする位置計測システム。 - 前記位置計測装置は、前記航走体から送信される前記速度情報と前記方位情報とを受信する送受信部を具備する請求項1記載の位置計測システム。
- 前記位置計測装置は、前記航走体から発せられる探査信号を受信し、前記探査信号の到着時刻としての第2時刻情報と、前記探査信号の前記母船に対する方向とを出力する受波器と、
前記航走体から発せられる前記探査信号の発信時刻としての第1時刻情報と前記第2時刻情報と前記方向とに基づいて前記相対位置情報を算出する制御部とを更に具備し、
前記送受信部は、更に、前記第1時刻情報を受信する請求項2記載の位置計測システム。 - 前記制御部は、更に、前記第1位置情報と、前記相対位置情報と、前記速度情報と、前記方位情報とに基づいて、前記第2位置情報を推定する、請求項3記載の位置計測システム。
- 前記制御部は、カルマンフィルタを用いて前記第2位置情報を推定する請求項4記載の位置計測システム。
- 前記制御部は、予め設定された演算時間ΔT毎の前記第2位置情報を下記式(1)乃至式(3)に基づいて計算するか、
又は下記式(1)及び式(2)を用いず、下記式(3)と予め設定されたカルマンフィルタゲインKとを用いて計算し、
計算された前記ΔT毎の前記第2位置情報に基づいて、下記式(4)により前記第2位置情報を推定する、
添え字k:k番目の演算時間
Kk:カルマンフィルタゲイン
XS:前記第1位置情報と前記相対位置情報とに基づいて求める航走体位置計測値
X^:前記第2位置情報
Xk^:演算時間ΔT毎の前記第2位置情報
φ:前記方位情報
A,B,C,D:前記航走体の運動モデルの状態空間表現
u:前記速度情報
Q:システム雑音wの分散
R:測定誤差nの分散
N:前記システム雑音wと前記測定誤差nの共分散
である)
請求項5記載の位置計測システム。 - 母船から送出された航走体の対地速度としての速度情報と、前記航走体の向かう方位としての方位情報とを前記航走体から受信する送受信部と、
前記航走体から発せられる探査信号を受信し、前記探査信号の到着時刻としての第2時刻情報と、前記探査信号の前記母船に対する方向とを出力する受波器と、
前記航走体から発せられる前記探査信号の発信時刻としての第1時刻情報と、前記第2時刻情報と、前記方向とに基づいて、前記相対位置情報を算出する制御部とを有する位置計測装置を具備し、
前記送受信部は、更に、前記第1時刻情報を受信し、
前記制御部は、更に、前記相対位置情報と、前記速度情報と、前記方位情報とに基づいて、前記航走体の位置としての第2位置情報を推定することを特徴とする位置計測システム。 - 前記制御部は、カルマンフィルタを用いて前記第2位置情報を推定する請求項7記載の位置計測システム。
- 前記制御部は、予め設定された演算時間ΔT毎の前記第2位置情報を下記式(5)乃至式(7)に基づいて計算し、計算された前記ΔT毎の前記第2位置情報に基づいて、下記式(8)により前記第2位置情報を推定する、
添え字k:k番目の演算時間
Kk:カルマンフィルタゲイン
XS:前記相対位置情報に基づいて求める航走体位置計測値
X^:前記第2位置情報
Xk^:演算時間ΔT毎の前記第2位置情報
φ:前記方位情報
A,B,C,D:前記航走体の運動モデルの状態空間表現
u:前記速度情報
Q:システム雑音wの分散
R:測定誤差nの分散
N:前記システム雑音wと前記測定誤差nの共分散
である)
請求項8記載の位置計測システム。 - 前記速度情報は、前記航走体の対地速度検出器及び対水速度検出器の少なくとも一方の計測結果に基づいている請求項7,8又は9記載の位置計測システム。
- 前記母船は船舶であり、前記航走体は水中を探査する潜水艇である請求項1乃至10のいずれか一項に記載の位置計測システム。
- 請求項1乃至11のいずれか一項に記載の位置計測システムと、前記位置計測システムを有する前記母船と、前記母船に制御され運動する前記航走体とを具備することを特徴とする航走体運用システム。
- 母船の位置としての第1位置情報を取得するステップと、
前記母船から送出された航走体の対地速度としての速度情報と、前記航走体の向かう方位としての方位情報と、前記航走体から発せられる前記探査信号の発信時刻としての第1時刻情報とを前記航走体から取得するステップと、
前記探査信号に基づいて、前記探査信号の到着時刻としての第2時刻情報と、前記探査信号の前記母船に対する方向とを取得するステップと、
前記第1時刻情報と、前記第2時刻情報と、前記方向とに基づいて、前記航走体の前記母船に対する位置である相対位置情報を算出するステップと、
前記第1位置情報と、前記相対位置情報と、前記速度情報と、前記方位情報とに基づいて、カルマンフィルタを用いて前記航走体の位置としての第2位置情報を推定するステップとを具備することを特徴とする位置計測方法。 - 母船の位置としての第1位置情報を取得するステップと、
前記母船から送出された航走体の対地速度としての速度情報と、前記航走体の向かう方位としての方位情報と、前記航走体から発せられる前記探査信号の発信時刻としての第1時刻情報とを前記航走体から取得するステップと、
前記探査信号の到着時刻としての第2時刻情報と、前記探査信号の前記母船に対する方向とを取得するステップと、
前記第1時刻情報と、前記第2時刻情報と、前記方向とに基づいて、前記航走体の前記母船に対する位置である相対位置情報を算出するステップと、
前記第1位置情報と、前記相対位置情報と、前記速度情報と、前記方位情報とに基づいて、カルマンフィルタを用いて前記航走体の位置としての第2位置情報を推定するステップとをコンピューターに実行させるためのプログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003017741A JP2004226374A (ja) | 2003-01-27 | 2003-01-27 | 位置計測システム、位置計測方法、プログラム及び航走体運用システム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003017741A JP2004226374A (ja) | 2003-01-27 | 2003-01-27 | 位置計測システム、位置計測方法、プログラム及び航走体運用システム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2004226374A true JP2004226374A (ja) | 2004-08-12 |
Family
ID=32904821
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003017741A Pending JP2004226374A (ja) | 2003-01-27 | 2003-01-27 | 位置計測システム、位置計測方法、プログラム及び航走体運用システム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2004226374A (ja) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008308077A (ja) * | 2007-06-15 | 2008-12-25 | Nec Corp | 制御装置、および速度制御方法 |
KR101026187B1 (ko) | 2008-06-23 | 2011-03-31 | 삼성중공업 주식회사 | 예인선의 부유체 예인 제어장치 |
JP2011116188A (ja) * | 2009-12-01 | 2011-06-16 | Nec Corp | 航走体位置報告装置、航走体位置報告方法、及びプログラム |
WO2013012218A1 (ko) * | 2011-07-15 | 2013-01-24 | 삼성중공업 주식회사 | 수중 이동체 위치측정장치 및 그 방법 |
WO2015129337A1 (ja) * | 2014-02-25 | 2015-09-03 | 古野電気株式会社 | 表層流推定装置、表層流推定システム、海洋モデル推定装置、及び危険度判定装置 |
JP2018072244A (ja) * | 2016-11-01 | 2018-05-10 | 三菱重工業株式会社 | 移動体観測システム、移動体観測方法及びプログラム |
WO2022054170A1 (en) * | 2020-09-09 | 2022-03-17 | Volvo Truck Corporation | Method and device for estimating coupling angle, tractor, computer program, and computer-readable medium storing computer program |
-
2003
- 2003-01-27 JP JP2003017741A patent/JP2004226374A/ja active Pending
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008308077A (ja) * | 2007-06-15 | 2008-12-25 | Nec Corp | 制御装置、および速度制御方法 |
KR101026187B1 (ko) | 2008-06-23 | 2011-03-31 | 삼성중공업 주식회사 | 예인선의 부유체 예인 제어장치 |
JP2011116188A (ja) * | 2009-12-01 | 2011-06-16 | Nec Corp | 航走体位置報告装置、航走体位置報告方法、及びプログラム |
US9341717B2 (en) | 2009-12-01 | 2016-05-17 | Nec Corporation | Apparatus for position notification of vehicle, method and computer-readable medium |
WO2013012218A1 (ko) * | 2011-07-15 | 2013-01-24 | 삼성중공업 주식회사 | 수중 이동체 위치측정장치 및 그 방법 |
KR101304579B1 (ko) * | 2011-07-15 | 2013-09-05 | 삼성중공업 주식회사 | 수중 이동체 위치측정장치 및 그 방법 |
US9335172B2 (en) | 2011-07-15 | 2016-05-10 | Samsung Heavy Ind. Co., Ltd | Apparatus for measuring location of underwater vehicle and method thereof |
AU2012284778B2 (en) * | 2011-07-15 | 2016-05-12 | Samsung Heavy Ind. Co., Ltd | Device for measuring location of underwater vehicle and method thereof |
WO2015129337A1 (ja) * | 2014-02-25 | 2015-09-03 | 古野電気株式会社 | 表層流推定装置、表層流推定システム、海洋モデル推定装置、及び危険度判定装置 |
JPWO2015129337A1 (ja) * | 2014-02-25 | 2017-03-30 | 古野電気株式会社 | 表層流推定装置、表層流推定システム、海洋モデル推定装置、及び危険度判定装置 |
JP2018072244A (ja) * | 2016-11-01 | 2018-05-10 | 三菱重工業株式会社 | 移動体観測システム、移動体観測方法及びプログラム |
WO2022054170A1 (en) * | 2020-09-09 | 2022-03-17 | Volvo Truck Corporation | Method and device for estimating coupling angle, tractor, computer program, and computer-readable medium storing computer program |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10983206B2 (en) | Low cost high precision GNSS systems and methods | |
Alcocer et al. | Underwater acoustic positioning systems based on buoys with GPS | |
US10261176B2 (en) | Rotating attitude heading reference systems and methods | |
Kepper et al. | A navigation solution using a MEMS IMU, model-based dead-reckoning, and one-way-travel-time acoustic range measurements for autonomous underwater vehicles | |
Kinsey et al. | In situ alignment calibration of attitude and Doppler sensors for precision underwater vehicle navigation: Theory and experiment | |
US20190204430A1 (en) | Submerged Vehicle Localization System and Method | |
Whitcomb et al. | Combined Doppler/LBL based navigation of underwater vehicles | |
Webster et al. | Preliminary deep water results in single-beacon one-way-travel-time acoustic navigation for underwater vehicles | |
CN103697910B (zh) | 自主水下航行器多普勒计程仪安装误差的校正方法 | |
Kinsey et al. | Adaptive identification on the group of rigid-body rotations and its application to underwater vehicle navigation | |
KR100906362B1 (ko) | 2개의 기준점에 대한 거리정보와 저정밀도 관성센서를 이용한 무인잠수정 선단의 의사 lbl 수중항법시스템 | |
CN113311388B (zh) | 一种水下机器人的超短基线定位系统 | |
Xu et al. | A novel self-adapting filter based navigation algorithm for autonomous underwater vehicles | |
Bo et al. | Cooperative localisation of AUVs based on huber-based robust algorithm and adaptive noise estimation | |
JP2008128968A (ja) | 水中測位システムおよび水中測位方法 | |
JP4250764B2 (ja) | 誘導制御システム、母船、誘導制御方法及びプログラム | |
JP2018084445A (ja) | 水中音響測位システム | |
Stateczny et al. | Precise bathymetry as a step towards producing bathymetric electronic navigational charts for comparative (terrain reference) navigation | |
JP2017165332A (ja) | 航走管理装置及び航走管理方法 | |
Jalving et al. | Terrain referenced navigation of AUVs and submarines using multibeam echo sounders | |
JP6716127B2 (ja) | 水中航走体の航走支援装置、航走体及び水中航走体の航走支援方法 | |
JP6946865B2 (ja) | 航走制御装置及び航走制御方法 | |
JP2004226374A (ja) | 位置計測システム、位置計測方法、プログラム及び航走体運用システム | |
CN111708008B (zh) | 一种基于imu和tof的水下机器人单信标导航方法 | |
Troni et al. | Field performance evaluation of new methods for in-situ calibration of attitude and Doppler sensors for underwater vehicle navigation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050216 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050713 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050912 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20051013 |