JP2005308428A - 位置推定装置 - Google Patents
位置推定装置 Download PDFInfo
- Publication number
- JP2005308428A JP2005308428A JP2004122505A JP2004122505A JP2005308428A JP 2005308428 A JP2005308428 A JP 2005308428A JP 2004122505 A JP2004122505 A JP 2004122505A JP 2004122505 A JP2004122505 A JP 2004122505A JP 2005308428 A JP2005308428 A JP 2005308428A
- Authority
- JP
- Japan
- Prior art keywords
- phase difference
- wake
- physical quantity
- calculated
- tracking
- 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
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
【課題】 マルチパスによる影響を受ける環境下で移動体が送信する電波を受信してこの移動体の位置を精度よく推定する。
【解決手段】 対象物からの電波を複数の受信器により受信して得た受信波間の位相差が整数値バイアスの異なる複数の受信波位相差から、追尾処理によって真の位相差を選択し対象物の位置を推定するために、複数の受信波位相差から第1及び第2の位相差を選択しこれらの基づく航跡を維持して平滑値を算出する追尾処理器22−1、22−2と、追尾処理器22−1、22−2が算出した平滑値から位相差とは異なる物理量であって対象物の運動状態を表す物理量を算出する物理量算出器23と、物理量算出器23により算出された物理量の連続性に基づいて第1及び第2位相差に基づく航跡とを入れ替える航跡入れ替え器24と、を備えた。
【選択図】 図1
【解決手段】 対象物からの電波を複数の受信器により受信して得た受信波間の位相差が整数値バイアスの異なる複数の受信波位相差から、追尾処理によって真の位相差を選択し対象物の位置を推定するために、複数の受信波位相差から第1及び第2の位相差を選択しこれらの基づく航跡を維持して平滑値を算出する追尾処理器22−1、22−2と、追尾処理器22−1、22−2が算出した平滑値から位相差とは異なる物理量であって対象物の運動状態を表す物理量を算出する物理量算出器23と、物理量算出器23により算出された物理量の連続性に基づいて第1及び第2位相差に基づく航跡とを入れ替える航跡入れ替え器24と、を備えた。
【選択図】 図1
Description
この発明は、電波を用いて対象物の位置を推定する装置に係るものであり、特に対象物が発する電波を複数の受信器で受信し、これら受信器間の受信波位相差を用いて対象物の位置を推定する技術に関する。
電波を用いて対象物の位置を計測する従来の方法としては、対象物から発信される電波を複数の受信器で受信し、これらの受信器の位置と受信器間の位相差、対象物の位置が満たす関係から対象物の位置を算出する、という方法がある(例えば、特許文献1)。
ここで、位相差はラジアン値で表されるので、2πのアンビギュイティが起こりうる。すなわち正しい位相差をΔφとした場合に、Δφ+2π×n(nは整数)となる位相差も同時に算出されうる。このように複数の位相差が算出される中で、対象物の位置に対応する位相差は一つだけである。複数の位相差から正しい位相差を見つけるために追尾技術を用いる方法が提案されている(例えば、非特許文献1)。
追尾技術は、過去に得られた値から予測値と平滑値とを外装法によって求め、現実の観測値との相関を評価して観測値の取捨選択を行う技術である。このような追尾技術によれば、各時刻(サンプル時刻と呼ぶ)において、複数の位相差を採りうる場合に、対象物に対応する最尤の位相差を選択することができる。
しかし、従来の追尾技術による位相差の選択方法は、受信器が対象物により送信される電波を直接波として受信する理想的な電波伝播環境を前提にしている。したがって、屋内や、屋外であっても壁面に囲まれた通路がある場合、位置測定の対象となる対象物が移動体である場合にこの対象物以外に併走している物体がある場合、そして地面や海面などによる乱反射によってマルチパスが発生する状況においては、上述のような選択方法を用いることは困難である。より具体的にいえば、マルチパス波の干渉によって位相差に誤差が生じ、正しい位相差と適切な相関が得られない状況が起こって、本来の位相差とは異なる位相差を選択してしまうおそれが強い。このような場合、対象物の位置に対応する位相差でない位相差に基づいて対象物の位置を計算することとなるため、位置の推定精度は著しく劣化することとなる。
この発明はこのような課題を解決するためになされたもので、マルチパス波等の干渉が生じる状況においても、その影響を軽減し、精度よく対象物の位置を推定することを目的とする。
対象物からの電波を複数の受信器により受信して得た受信波間の位相差が整数値バイアスの異なる複数の受信波位相差を含む場合に、これらの位相差に基づく追尾処理を行って上記対象物の位置推定に用いる平滑値を算出し、この平滑値から上記対象物の位置を推定する位置推定装置において、
上記複数の受信波位相差から第1位相差を選択しこの第1位相差に基づく航跡を維持して平滑値を算出する第1追尾処理器と、
上記第1位相差に隣接する第2位相差を選択してその第2位相差に基づく航跡を維持する第2追尾処理器と、
上記第1追尾処理器が算出した平滑値から位相差とは異なる物理量であって上記対象物の運動状態を表す物理量を算出する物理量算出器と、
上記物理量算出器により算出された物理量の連続性に基づいて上記第1位相差に基づく航跡と上記第2位相差に基づく航跡とを入れ替える航跡入れ替え器と、
を備えたものである。
上記複数の受信波位相差から第1位相差を選択しこの第1位相差に基づく航跡を維持して平滑値を算出する第1追尾処理器と、
上記第1位相差に隣接する第2位相差を選択してその第2位相差に基づく航跡を維持する第2追尾処理器と、
上記第1追尾処理器が算出した平滑値から位相差とは異なる物理量であって上記対象物の運動状態を表す物理量を算出する物理量算出器と、
上記物理量算出器により算出された物理量の連続性に基づいて上記第1位相差に基づく航跡と上記第2位相差に基づく航跡とを入れ替える航跡入れ替え器と、
を備えたものである。
ここで、追尾処理とは一般に、所定のサンプリング間隔で得られる観測値と前回サンプル時以前(前回を含む)の観測値及び平滑値から得られた今回サンプル用の予測値との相関を評価し、相関が認められた観測値に基づいて今回サンプル時の平滑値と次回サンプル時の予測値を算出する演算をサンプル毎に繰り返す処理をいう。そして、追尾処理においては、これまでのサンプル時に得られた観測値、予測値と平滑値は次回以降の予測値及び平滑値の算出に影響を及ぼすことが多いので、このような観測値、予測値と平滑値の時系列を、一般に航跡と呼んでいる。
また上述の位置推定装置において、位相差とは異なる物理量であって対象物の運動状態を表す物理量とは、例えば時間、位置、速度、角速度、加速度及び角加速度など、対象物の運動状態を表現する上で用いることのできる何れかの状態変数を指す。
このように、この発明に係る位置推定装置によれば、位相差にアンビギュイティが生じる状況において、位相差に基づく追尾処理のみでは対象物の位置に対応しない誤った位相差を選択してしまう状況であっても、位相差とは異なる対象物の物理量を算出しその物理量の整合性を評価することで、誤った位相差を選択することを抑制し、精度の高い位置推定処理を実現することができるのである。
さらにこの発明に係る位置推定装置においては、位相差にアンビギュイティがあることに基づいて上述のような構成としている。しかしながら、方位角を算出することのできるアレーアンテナを用いて構成した場合においては、位相差に替えて方位角のアンビギュイティに対しても同一の技術的思想を適用することができるのである。すなわち、かかる位置推定装置は、対象物からの電波を複数の受信器により受信して得た受信波から整数値バイアスの異なる複数の方位角が得られる場合に、これらの方位角に基づく追尾処理を行って上記対象物の位置推定に用いる主方位角を算出し、この主方位角から上記対象物の位置を推定する位置推定装置において、
上記複数の方位角から第1方位角を選択しこの第1方位角に基づく航跡を維持して上記主方位角を算出する第1追尾処理器と、
上記第1方位角に隣接する第2方位角を選択してこの第2方位角に基づく航跡を維持する第2追尾処理器と、
上記主方位角から方位角とは異なる物理量であって上記対象物の運動状態を表す物理量を算出する物理量算出器と、
上記物理量算出器により算出された物理量の連続性に基づいて上記第1方位角に基づく航跡と上記第2方位角に基づく航跡とを入れ替える航跡入れ替え器と、
を備えたもの、とすることもできる。
上記複数の方位角から第1方位角を選択しこの第1方位角に基づく航跡を維持して上記主方位角を算出する第1追尾処理器と、
上記第1方位角に隣接する第2方位角を選択してこの第2方位角に基づく航跡を維持する第2追尾処理器と、
上記主方位角から方位角とは異なる物理量であって上記対象物の運動状態を表す物理量を算出する物理量算出器と、
上記物理量算出器により算出された物理量の連続性に基づいて上記第1方位角に基づく航跡と上記第2方位角に基づく航跡とを入れ替える航跡入れ替え器と、
を備えたもの、とすることもできる。
このように、この発明に係る位置推定装置によれば、方位角にアンビギュイティが生じる状況において、方位角に基づく追尾処理のみでは対象物の位置に対応しない誤った方位角を選択してしまう状況であっても、方位角とは異なる対象物の物理量を算出しその物理量の整合性を評価することで、誤った方位角を選択することを抑制し、精度の高い位置推定処理を実現することができるのである。
以下、この発明の実施の形態を図を用いて説明する。
実施の形態1.
図1は、この発明の実施の形態1による位置推定装置の構成を示すブロック図であある。図において、対象物1には電波送信機が装着されている。位置推定装置2はこの発明の対象となる位置推定装置であって、M個の受信器21−1〜M、追尾処理器22−1〜3、物理量算出器23、航跡入れ替え器24を備えている。
実施の形態1.
図1は、この発明の実施の形態1による位置推定装置の構成を示すブロック図であある。図において、対象物1には電波送信機が装着されている。位置推定装置2はこの発明の対象となる位置推定装置であって、M個の受信器21−1〜M、追尾処理器22−1〜3、物理量算出器23、航跡入れ替え器24を備えている。
受信器21−iと受信器21−jはM個(Mは2以上の自然数)ある受信器21−1〜Mの中から選択された2個の受信器であって(i,j=1,2,…,M)、例えば対象物1に装着された電波送信機が送信する電波を受信する受信アンテナ、局部発信機、ダウンコンバータ、A/D変換器などを用いて構成され、最終的に受信した電波信号を所定のサンプリング間隔でデジタル信号として出力するようになっている。ここで、受信器21−iと21−jによって出力されるデジタル信号riとrjとの位相差Δφij、対象物1の送信する電波信号の波長をλ、対象物1の未知なる位置pをp=(x y z)、そして受信器21−iと21−jの位置を(xi yi zi)、(xj yj zj)とするならば、受信器21−i、21−jと対象物1の位置p、位相差Δφijは式(1)を満たす。
このようにすることで、対象物1が移動しても、一定のサンプリング間隔ごとに対象物1の位置が算出されるのである。なお、式(1)において、lijは位相差Δφijのアンビギュイティを表す整数バイアス値である。すなわち、位相差は角度値として表されるため、角度値としてはlijによらずΔφij+2πlijは同じ角を表すこととなる。ところが対象物1の位置は1つしかないので、式(1)によれば、対象物1の位置に対応するlijはたった一つしか存在しないはずである。このことから、誤ったlijを選択して式(1)に代入してしまうと、誤った位置を算出する結果を招くこととなる。
この問題を解決するための方法として、この発明の実施の形態1の位置推定装置は追尾処理により複数の位相差の中から正しい位相差を選択する。いま受信器21−iと21−jとの受信信号の位相差として得られる複数の位相差の集合を式(2)で示すUiで表すこととする。
追尾処理器22−1〜3は位相差をそれぞれ追尾処理する回路又は素子である。ここで追尾処理器22−1は追尾処理によって正しい位相差であると推定された位相差(第1位相差)の航跡を維持する追尾処理回路または素子であって、最終的に位相差を平滑した値である位相差平滑値を算出するようになっている。また追尾処理器22−2及び22−3は追尾処理器22−1が選択する第1位相差に隣接する第2位相差の航跡を維持する追尾処理回路または素子である。
図2は、追尾処理器22−1を、例えばカルマンフィルタによる追尾処理を行うものとして構成した場合の詳細な構成を示したブロック図である。図に示すように、追尾処理器22−1は位相差検出器25、相関判定器26、平滑器27、相関予測器28、位相差ゲイン算出器29、範囲算出器30を備えている。
このうち、位相差検出器25は、受信器21−i、21−jの出力するデジタル信号から位相差を検出する回路または素子である。また相関判定器26は、位相差検出器25が検出した位相差が後述する相関範囲算出器30によって設定される相関範囲に含まれるか否かを判定し、含まれると判断した場合は新たな観測値として追尾処理に採用し、含まれないと判断した場合はこの位相差を棄却して、代わりに予測器28が算出した予測値を出力する回路又は素子である。平滑器27は相関判定器26が出力するサンプル毎の位相差又は予測値に基づいてそのサンプルに対する位相差平滑値を算出する回路又は素子である。ここで算出される位相差平滑値は、追尾処理器22−nの出力値として外部に供給される。予測器28は、前回サンプル時の位相差平滑値に基づいて今回サンプルにおいて観測される位相差の予測値を算出する回路又は素子である。ただし前回サンプル時の位相差平滑値を直接の入力値とするのではなく、追尾処理器22−nの外部にある航跡入れ替え器24(後述)から取得した航跡の位相差平滑値を用いるようになっている。位相差ゲイン算出器29は、平滑器27が位相差平滑値の算出に用いる位相差ゲインを算出する回路又は素子である。相関範囲算出器30は、予測器28が算出した予測値から許容される今回サンプルの位相差の観測値の相関範囲を算出する回路又は素子である。
なお追尾処理器22−2及び22−3も追尾処理器22−1と同等の構成要素を備えているのでその構成については特に示さないこととする。
引き続き、図1に示した位置推定装置2の構成について説明する。物理量算出器23は、追尾処理器22−1〜3が算出した各位相差平滑値に基づいて、対象物1に関する物理量を算出する回路又は素子である。ただしここで算出される物理量とは位相差とは異なる物理量であって、対象物1の運動状態を示す状態量である。
図3は、物理量算出器23の例として対象物1の速度を算出する構成とした場合の詳細な構成を示すブロック図である。図において、測位演算器31は追尾処理器22−1〜3が算出した位相差平滑値に基づいて対象物1の位置を算出する回路又は素子である。なおこの発明の実施の形態1による位置推定装置2の出力結果である対象物の位置は、この部位の演算結果として求められる。速度算出器32は、測位演算器31が算出した対象物1の位置の時間変位から速度を算出する回路又は素子である。速度判定器33は、速度算出器32によって算出された対象物1の速度が所定の値域に含まれるかどうかを判断する部位である。
航跡入れ替え器24は、追尾処理器22−1〜3が使用する予測値を調整する回路又は素子である。すなわち、これらの予測値を調整することによって追尾処理器22−1〜3の航跡の入れ替えを行うようになっている。この部位は、前回サンプルまでに算出された予測値を1サンプルの間記憶しておき、さらに物理量算出器23の速度判定器33の判定結果に基づいて追尾処理器22−1〜3に従来の航跡をそのまま指定したり、入れ替えて指定する。
次に、位置推定装置2の動作について説明する。まず初めに位置推定装置2の動作原理について説明し、続いて具体的な動作の説明を行うこととする。対象物1から発信される電波は、位置推定装置2の受信器21−1〜21−Mによって受信され、検波処理を経てA/D変換されてデジタル信号として出力される。デジタル信号はサンプル毎に出力されるので、このデジタル信号から位相差を検出し、追尾処理に供するのである。ここで、ある受信器21−iと21−jに着目し、この受信器21−iおよび21−jの受信デジタル信号riとrjの位相差をΔΦijとし、この位相差の追尾処理について説明する。
図4は、追尾処理に基づいて複数の位相差から対象物1の位置に対応する位相差を選択する状況を説明する図である。図において、*(アスタリスク)は対象物1の位置に対応する本来の位相差を表す点である。また×で表される点は、本来の位相差と整数バイアス値分だけ異なる位相差である。さらに矩形で囲まれた位相差は、追尾処理に基づいて算出した位相差平滑値を表している。追尾処理によれば、一時的に対象物1の位置に対応する本来の位相差を選択しうるが、常に本来の位相差を選択することは保証されない。つまり、例えば位相差として、前回の位相差とは近い位相差を選択するため、本来の位相差ではなく、実は誤った位相差を選択してしまうことが生じうる。図4の例で言えば、区間101において本来の位相差と離れて隣接する整数値バイアスを有する位相差に航跡が遷ってしまう。このような結果は、位相差のみ相関処理で位相差を選択していたことに起因するものである。
そこで、位相差のみの相関処理ではなく、位相差から得られた対象物1の運動状態、例えば速度の相関を考慮して位相差を選択するようにする。図5は、図4の位相差平滑値から算出される対象物1の速度の時間変位を表すグラフである。図が示すように航跡が誤った位相差に乗り移る時間帯101において、速度が大きく変化している。仮に対象物1の速度が取りうる値の範囲を定めることができる場合、選択した位相差から算出した速度がこの範囲を逸脱するような場合には、不適切な位相差を選択したことが判断できるのである。
例えば、対象物1が競走馬の場合、速度の上限値はおよそ時速60km程度であり、さらに下限値はおよそ時速40km程度であることが知られている。したがって例えば対象物1の速度が時速40km〜60kmを逸脱するような位相差は取り得ないはずであり、そのような条件に基づく位相差の相関処理を行うことで、適切な位相差を選択する可能性を大きく向上し、結果として位置の推測精度を高めることにつながるのである。
以上の原理を詳細に説明する。追尾処理器22−nは位相差の追尾処理を行うが、ここでの追尾処理の運動モデルは式(4)で表されるものとする。
ここで、x kはサンプリング時刻tkにおける位相差の真値θkとその時間微分の真値θ'k(θ'kはθkの時間による一回微分)を表す状態ベクトルであり、式(5)で表される。
また、サンプリング時刻tk(k=1,2,…,n)は、位相差を検出するタイミングを表す離散時刻であり、以降、時刻tkと呼ぶ。
wkは時刻tkにおける駆動雑音であり、Γ1(k)は時刻tkにおける駆動雑音の変換行列である。例えば、位相差の運動モデルを等速直線運動と仮定したことによる打ち切り誤差項をΓ1(k−1)wk−1とみれば、wkは位相差の加速度ベクトルに相当する項であり、またΓ1(k−1)は式(7)で表される。なお、Tはサンプリング間隔である。平均を表す記号としてEを用いると、wkは平均0、分散Q kの正規分布の白色雑音とする。
一方、位相差の観測モデルは式(8)で表されるものとする。
ここで、ukは位相差の観測値、Hは観測行列で、式(9)で表される。
なお、式(8)において、ν kは時刻tkにおいて検出された位相差の観測雑音である。観測雑音は平均0、分散Rkの正規分布白色雑音に従うものと仮定する。また、時刻tkまでの間に位相差の追尾に用いた位相差の全体をUkとする。
以上の運動モデルに基づき、カルマンフィルタを用いた場合を例に追尾フィルタを構成する。1サンプリング前の位相差航跡の平滑ベクトルをx^k−1(+)(x^は数式において文字xの上に^を付した変数を表すものとする)、平滑誤差共分散行列をPk−1(+)とすると、時刻tkにおける位相差の状態ベクトルx kの予測ベクトルx^(−)及び予測誤差共分散行列Pk(−)はそれぞれ条件付平均ベクトル、条件付共分散行列で定義され、式(10)及び式(11)で表される。
次に、時刻tkに計測された複数の位相差Ukから対象物1の位置に対応する位相差を選択する方法について述べる。偽の位相差の含まれる複数の位相差から正しい位相差を判定するため、式(12)で表される位相差の予測値u k(−)を利用する。
過去の位相差情報Uk−1より得られる位相差の観測値の予測確率分布P[uk|Uk−1]は、式(13)に示す条件付確率密度関数で表される。
すなわち、位相差の観測値は、式(11)で与えられる位相差の予測値uk(−)を平均とし、式(14)で与えられるSkを分散とする1次元正規分布g(uk;uk(−),Sk)に従うものとする。
そして、式(15)の左辺値fiが最小となるような位相差un kを正しい位相差(第1位相差)と判断して選択するのである。さらにここでは、選択された位相差の前後の整数値バイアス値に対応する位相差un+1 kのいずれかを第2位相差とする。またun-1 kを第3位相差としてもよい。なお、式(15)において、dは目標との相関範囲を決めるパラメータであり、自由度1のχ自乗分布により算出する値である。
続いて、第1位相差の航跡と第2位相差及び第3位相差の航跡に対応する平滑ベクトルの算出方法について述べる。ゲイン行列Kk、位相差航跡の平滑ベクトルx^k(+)及び平滑誤差共分散行列Pk(+)は、カルマンフィルタの理論により式(16)〜式(18)で与えられる。
同様に、正しい位相差の前後の整数値バイアス値に対応する位相差un-1 k及びun+1 kをそれぞれ用いた位相差裏航跡の平滑ベクトルをx^a k(+)及びx^b k(+)とすると、これらは式(19)及び式(20)で与えられる。
次に、位相差航跡の位置または速度から、追尾が偽の位相差に乗り移った事象を判定し、位相差航跡を軌道修正する方法を、ここでは速度を利用した場合を例にとって説明する。第1位相差の航跡として選択された位相差un kを用いて算出した移動体の速度をVn k、第2位相差、第3位相差の航跡として選択された位相差un-1 k及びun+1 kを用いて算出した移動体の速度をそれぞれ、Vn-1 k及びVn+1 kとする。
乗り移りがあったと判定された場合、Vn kの代わりにVn-1 k及びVn+1 kについても式(21)及び式(22)の判定を行い、これらの条件を満たさない速度に対応する位相差で算出した航跡の平滑ベクトルを、位相差航跡として採用する。例えばVn-1 kが式(21)及び式(22)の条件を満たさない場合は、式(23)としてx^a k(+)を新たな位相差航跡とする。
そして、この時点でこれまでの第2位相差及び第3位相差の航跡を削除して、次サンプル時刻からは新たな位相差航跡を基準に、これまで同様、前後の整数値バイアスに対する位相差を用いて第2位相差及び第3位相差の航跡を生成するのである。
以上の処理方法に基づいて追尾処理器22−nと追尾処理器22−n+1、さらには物理量算出器23と航跡入れ替え器24の具体的な動作を示すと次のようになる。
追尾処理器22−1において、相関判定器26は、位相差検出器25から位相差を受け取ると、整数値バイアスを考慮した複数の位相差(偽の位相差を含む)に対して、相関範囲算出器30で算出された位相差の観測値の予測確率分布の平均値uk(−)と分散Sk、予め設定されたパラメータdを用いて、式(15)を満たす位相差のうち、fiが最小となる位相差un kを正しい位相差と判断して選択する。また追尾処理器22−2及び22−3は、追尾処理器22−1が選択した位相差un kの前後の整数値バイアス値に対応する位相差un-1 k及びun+1 kをそれぞれの位相差の航跡を生成するためのデータとして選択する。
追尾処理器22−1において、位相差ゲイン算出器29は、予測器28から取得した予測誤差共分散行列Pk(−)と予め設定された受信器21−1〜21−Mの観測雑音の分散Rkとから、式(16)に従い位相差のゲイン行列Kkを算出する。そして平滑器27は、位相差ゲイン算出器29から取得した位相差のゲイン行列Kk、予測ベクトルx^k(−)、予測誤差共分散行列Pk(−)と、相関判定器26から取得した位相差航跡に対応する位相差un kと、位相差の予測値uk(−)とを用いて、式(17)と式(18)に従い、位相差航跡の平滑ベクトルx^k(+)及び平滑誤差共分散行列Pk(+)を算出する。
同様に、追尾処理器22−2及び22−3においても、相関判定器26から取得した位相差の航跡に対応する位相差un-1 k及びun+1 k、位相差の予測値uk(−)と、位相差ゲイン算出器29から取得した位相差のゲイン行列Kk、予測ベクトルx^(−)とを用いて、式(19)と式(20)に従い、それぞれの位相差の航跡の平滑ベクトルx^a k(+)及びx^b k(+)を算出する。
以上のようにして位相差検出単位毎に選択された第1位相差の航跡と第2及び第3位相差の航跡に対応する位相差un k、un-1 k、un+1 kは、物理量算出器23に出力される。物理量算出器23の測位演算器31は、各位相差毎に式(1)の連立方程式を解き、測位解(位置)を算出する。
続いて速度算出器32は、測位演算器31で算出した各位相差に対応する位置と過去の位置情報から、各位相差に対応する速度Vn k、Vn-1 k、Vn+1 kを算出する。そして速度判定器33は、予め設定された移動体の速度の上限Vmax、下限Vmin及びその誤差σとから、速度算出器32から取得した各位相差に対応する速度Vn k、Vn-1 k、Vn+1 kが式(21)及び式(22)を満たすか否かを判定する。
航跡入れ替え器24は、速度判定器33において、位相差航跡の位相差un kに対応する速度Vn kが式(21)及び式(22)を満たした場合には、位相差航跡の位相差un kに対応する測位解を現サンプリング時刻における測位結果とし、満たさない場合は、式(21)及び式(22)を満たす位相差裏航跡の位相差が正しい位相差であると判定し、これに対応する測位解を現サンプリング時刻における測位解とする。航跡入れ替え器24は、速度判定器33が第2及び第3位相差の航跡の位相差が正しい位相差であると判定された場合にのみ、式(23)により位相差航跡の平滑ベクトルx^k(+)を位相差裏航跡の平滑ベクトルに置き換えるとともに、これまで算出してきた位相差裏航跡の平滑ベクトルを初期化する。
予測器28は、航跡入れ替え器24から取得した位相差航跡の平滑ベクトルx^k(+)、平滑誤差共分散行列Pk(+)と、予め設定された駆動雑音の誤差共分散行列Qkとから、式(10)及び式(11)に従い、次サンプリング時刻における位相差の予測ベクトルx^k+1(−)及び予測誤差共分散行列P^k+1(−)を算出する。相関範囲算出器30は、予測器28で算出された予測誤差共分散行列Pk+1(−)と、予め設定された受信器21−1〜Mの観測雑音の分散Rkとから、式(14)に従い位相差の観測値の予測確率分布の分散Skを算出する。
以上から明らかなように、この発明の実施の形態1の位置推定装置によれば、追尾処理器22−1だけでなく、追尾処理器22−2あるいは22−3を設け、第1位相差の航跡と、この位相差の前後の整数値バイアスに対応する第2の位相差とを航跡として維持して平滑することとした。そして、速度判定器26において、位相差航跡として選択された位相差により算出された測位解より算出された対象物1の速度が、移動体の取りうる速度範囲内であるかを判定して、範囲を超えている場合に、航跡入れ替え器24において、速度が正常な位相差を用いて算出した測位解を最終的な測位結果とすることとした。かかる構成を採用した結果、マルチパス波の干渉により位相差の誤差が発生し、位相差の追尾が偽の位相差に乗り移っても、追尾を真の位相差に修正するため、測位精度の劣化を防ぐことができる。
上述の処理では、最良の実施の形態として速度の上限と下限が定められる構成について説明した。しかし速度の値域については上限と下限のいずれか一方しか定めることができない場合であっても式(21)と式(22)のいずれか一方のみを選択して、速度を評価することで、不適切な位相差を検出することができるので、上限と下限の双方が定められる構成に限定されるものではない。
また上述の処理では、同じく最良の実施の形態として、第1位相差の前後の位相差の航跡として第2位相差及び第3位相差を2つとも維持する構成について説明した。しかし対象物1の位置測定環境における電波の伝播特性から、マルチパスによる電波の干渉が位相差に及ぼす影響が予測可能となって、第1位相差よりも大きな位相差あるいは第1位相差よりも小さな位相差に限定できる場合もある。このような場合においては、第2位相差と第3位相差のいずれか一方だけを追尾すれば足りるのである。したがってこのような場合には、追尾処理器22−1及び22−2のみで十分である。
実施の形態2.
実施の形態1による位置推定装置では、正しいと判定された位相差un-1 kで算出した速度に矛盾があった場合に、位相差un kの前後の整数値バイアス値に対応する位相差un-1 k及びun+1 kのうち、速度が正常となる位相差を正しい位相差として選択し直す構成とした。しかしながら、道路やトラックを周回する移動体のように、移動体の移動経路が限定されている場合は、対象物1がある基準点を通過する時間を予測して、この予測時間の整合性に基づいて位置推定結果の矛盾を判定するようにしてもよい。この発明の実施の形態2による位置推定装置は、かかる特徴を有するものである。
実施の形態1による位置推定装置では、正しいと判定された位相差un-1 kで算出した速度に矛盾があった場合に、位相差un kの前後の整数値バイアス値に対応する位相差un-1 k及びun+1 kのうち、速度が正常となる位相差を正しい位相差として選択し直す構成とした。しかしながら、道路やトラックを周回する移動体のように、移動体の移動経路が限定されている場合は、対象物1がある基準点を通過する時間を予測して、この予測時間の整合性に基づいて位置推定結果の矛盾を判定するようにしてもよい。この発明の実施の形態2による位置推定装置は、かかる特徴を有するものである。
この発明の実施の形態2による位置推定装置の構成は実施の形態1と同様に図1によって表される。ただし実施の形態2による位置推定装置は実施の形態1と対比した場合に、物理量算出器23の構成が異なっている。図6は、実施の形態2による位置推定装置2の物理量算出器23の詳細な構成を示すブロック図である。図において、通過判定器34は、測位演算器31の出力する対象物1の位置から対象物1が基準点を通過したかどうかを判定する回路又は装置である。通過時刻予測器35は、測位演算器31の出力する対象物1の位置の時間変位から、基準位置の通過時刻を算出する回路又は素子である。通過時刻判定器36は、通過時刻算出器34が算出した通過時刻を評価し、測位演算器31が位置の算出に用いた位相差平滑値の妥当性を判定する回路又は素子である。
ここで、対象物1の基準点通過が判定された時刻をt~とし、その直前の追尾処理において、位相差航跡として選択された第1位相差un kを用いて算出した対象物1の位置及び速度をそれぞれpn k及びVn k、第2位相差un-1 kに対応する位置及び速度をそれぞれpn-1 k及びVn-1 k、第3位相差un+1 kに対応する位置及び速度をそれぞれpn+1 k及びVn+1 k、基準点位置をPrefとすると、時刻t~における第1〜第3位相差の航跡に対応する移動体の予測位置は式(24)〜式(26)で表される。
このとき、予測位置と基準点位置の差は式(27)〜式(29)で表されるので、この差が最も小さいものが、対象物1の真の位相差航跡と判定できる。
以上のような原理に基づくこの発明の実施の形態2による位置推定装置2の具体的な動作を以下に説明する。実施の形態1と同様に、位置推定装置2においても追尾処理器22−1〜3において第1位相差〜第3位相差の追尾処理を行う。その結果、各位相差平滑値un k、un-1 k、un+1 kがそれぞれの追尾処理器により算出され、物理量算出器23に出力される。物理量演算器23の測位演算器31は、取得した位相差平滑値のそれぞれについて式(1)の連立方程式を解き、測位解(位置)を算出する。
通過判定器34は、移動体が起点を通過した時刻を記録し、その時刻を通過時刻予測手段に出力する。そして通過時刻予測器35は、通過判定器34から取得した移動体が基準点を通過した時刻t~と、測位演算器31から取得した対象物1の位置pn k及び速度Vn k、位置pn-1 k及び速度Vn-1 k、位置pn+1 k及び速度Vn+1 kとを用いて、式(24)〜(26)に従って、時刻t~における位相差航跡及び位相差裏航跡に対応する移動体の予測位置を算出する。通過時刻判定器36は、通過時刻予測器35から取得した時刻t~における第1〜第3位相差の航跡に対応する対象物1の予測位置と、予め設定された基準点位置prefを用いて、式(27)〜式(29)に従い予測位置と基準点位置の差を算出し、この中で最も小さい値を持つ予測位置に対応する位相差航跡または位相差裏航跡を、対象物1の真の位相差航跡と判定する。航跡入れ替え手段24は、通過時刻判定器36によって対象物1の真の位相差航跡と判定された予測位置が第2及び第3の位相差の航跡である場合には式(23)により航跡の平滑ベクトルx^k(+)を位相差裏航跡の平滑ベクトルに置き換えるとともに、これまで算出してきた第2及び第3の位相差の航跡の平滑ベクトルを初期化する。
予測器28は、航跡入れ替え器24から取得した位相差航跡の平滑ベクトルx^k(+)、平滑誤差共分散行列Pk(+)と、予め設定された駆動雑音の誤差共分散行列Qkとから、式(10)及び式(11)に従い、次サンプリング時刻における位相差の予測ベクトルx^k+1(−)及び予測誤差共分散行列Pk+1(−)を算出する。相関範囲算出器30は、予測器28で算出された予測誤差共分散行列Pk+1(−)と、予め設定された受信器21−1〜Mの観測雑音の分散Rkとから、式(14)に従い位相差の観測値の予測確率分布の分散Skを算出する。
以上から明らかなように、この発明の実施の形態2の位置推定装置によれば、通過時刻予測器35において算出された第1〜第3の位相差航跡の基準点通過時刻における予測位置と基準点位置の差から、対象物1の真の位相差を判定し、航跡を入れ替えることとしたので、マルチパス波の干渉により位相差の誤差が発生し、位相差の追尾が偽の位相差に乗り移っても、追尾を真の位相差に修正するため、測位精度の劣化を防ぐことができる。
なお、この発明の実施の形態2においても、電波の伝播条件を仮定できる場合においては、第1位相差の前後の位相差の航跡の一方のみを維持すれば足りる点で、実施の形態1と同様である。
また、実施の形態1と2では速度に上限値や下限値を設定したり、基準位置の通過時間に上限値や下限値を設定して、位置推定の対象となる対象物の位相差の航跡の乗り移りを検出しているが、速度や通過時間といった条件は一例に過ぎない。すなわち対象物の運動状態を示すある物理量が位相差に基づいて導出可能であり、さらにある対象物についてその物理量の時間変化が何らかの整合性や連続性を仮定できるのであれば、如何なる物理量であってもこの発明による技術的思想に包含されるのである。
実施の形態3.
実施の形態1及び2による位置推定装置では、位相差について追尾処理を行い、平滑値と予測値とを算出している。しかしながら、位相差に替えて対象物の直交座標による位置を追尾するような構成を採用してもよい。この実施の形態3による位置推定装置はかかる特徴を有するものである。
実施の形態1及び2による位置推定装置では、位相差について追尾処理を行い、平滑値と予測値とを算出している。しかしながら、位相差に替えて対象物の直交座標による位置を追尾するような構成を採用してもよい。この実施の形態3による位置推定装置はかかる特徴を有するものである。
この発明の実施の形態3による位置推定装置の構成は実施の形態1及び2と同様に図1によって表される。ただし実施の形態3による位置推定装置は追尾処理器22−1〜3及び物理量算出器23の詳細な構成が異なっている。図7は、実施の形態3による位置推定装置2の追尾処理器22−1の詳細な構成を示すブロック図である。図において、直交座標ゲイン算出器37は直交座標のゲインを算出する回路又は素子である。他の部位については実施の形態1及び2と同様であるが、後述するように、この追尾処理器22−1は対象物1との距離差を直交座標によって追尾するように構成されているので、処理に供されるデータが異なっている。
なお、追尾処理器21−2及び22−3の構成については、追尾処理器21−1と同様であるので説明を省略する。
図8は、この発明の実施の形態3における位置推定装置の物理量算出器23の詳細な構成を示すブロック図である。図において、速度判定器38は追尾処理器22−1〜22−3が算出した直交座標の平滑値の時間変位から速度を求めて判定する回路又は素子である。
ここで、この発明の実施の形態3における位置推定装置の追尾処理における運動モデルは式(30)で表される。
ここで、x kはサンプリング時刻tkにおける直交座標の位置の真値とその時間微分の真値を表す状態ベクトルであり、式(31)で表される。
またサンプリング時刻tk(k=1,2,…,n)は、位相差を検出するタイミングを表す離散時刻であり、以降、時刻tkと呼ぶ。Φk−1は時刻tk−1よりtkへの状態ベクトルの推移行列であり、式(32)で表される。
式(30)において、wkは時刻tkにおける駆動雑音であり、Γ1(k)は時刻tkにおける駆動雑音の変換行列である。例えば、移動体の運動モデルを等速直線運動と仮定したことによる打ち切り誤差項をΓ1(k−1)wk−1とみれば、wkは移動体の加速度ベクトルに相当するものである。また、Γ1(k−1)は式(33)で表される。
なお、Tはサンプリング間隔、Iは3行3列の単位行列である。平均を表す記号としてEを用いると、wkは平均0、共分散Qkの正規分布の白色雑音とする。
また、対象物1の観測モデルは式(34)で表される。
ここで、ukは距離差の観測値、h(x k)は観測行列で、式(35)で表される。
ただし、ν kは時刻tkにおいて検出された距離差の観測雑音である。観測雑音は平均0、分散Rkの正規分布白色雑音に従うものと仮定する。また、時刻tkまでの間に移動体の追尾に用いた距離差の全体をUkとする。
以上の運動モデルに基づき、拡張カルマンフィルタを用いて追尾フィルタを構成する。1サンプリング前の直交座標航跡の平滑ベクトルをx^k−1(+)、平滑誤差共分散行列をPk−1(+)とすると、時刻tkにおける移動体の状態ベクトルxkの予測ベクトルx^k(−)及び予測誤差共分散行列Pk(−)はそれぞれ条件付平均ベクトル、条件付共分散行列で定義され、式(36)及び式(37)で表される。
さらに、このようにして選択された距離差un kを算出した第1位相差の前後の整数値バイアス値をとる第2及び第3位相差から算出される距離差un-1 k及びun+1 kについても直交座標を航跡を生成するためのデータとして選択するのである。なお、以後の説明においては、距離差un kを第1距離差と呼び、距離差un-1 kを第2距離差、距離差un+1 kを第3距離差と呼ぶこととする。式(41)において、dは目標との相関範囲を決めるパラメータであり、自由度3のχ自乗分布により算出される。
また、予測ベクトル値x^k(−)は、式(42)と式(43)とを満たす。 続いて、第1〜第3位相差に対応する直交座標の航跡における平滑ベクトルの算出方法について述べる。第1距離差のゲイン行列Kk、直交座標航跡の平滑ベクトルx^k(+)及び平滑誤差共分散行列Pk(+)は拡張カルマンフィルタの理論により、式(44)〜式(46)で与えられる。
次に、直交座標航跡の位置または速度から、追尾が偽の位相差に乗り移った事象を判定し、直交座標航跡を軌道修正する方法を、ここでは速度を利用した場合を例にとって説明する。第1位相差の直交座標航跡の速度をVn k、第2及び第3位相差の直交座標航跡に対応する速度をそれぞれVn-1 k及びVn+1 kとする。
その後、第2及び第3位相差についての直交座標航跡を削除して、次サンプル時刻からは新たな直交座標航跡を基準に、これまで同様、前後の整数値バイアスに対する距離差を用いて直交座標航跡を生成する。
以上の処理方法に基づいて、この発明の実施の形態3における追尾処理器22−nと追尾処理器22−n+1、さらには物理量算出器23と航跡入れ替え器24の具体的な動作を示すと次のようになる。
追尾処理器22−1において、相関判定器26は、位相差検出器25から位相差を受け取ると、整数値バイアスを考慮した複数の距離差(偽の距離差を含む)に対して、相関範囲算出器30で算出された距離差の観測値の予測確率分布の平均値uk(−)と分散Sk、予め設定されたパラメータdとを用いて、式(41)を満たす距離差のうち、fiが最小となる距離差un kを正しい距離差と判断して選択する。また追尾処理器22−2及び22−3は、その位相差の前後の整数値バイアス値をとる位相差に対応する距離差un-1 k及びun+1 kを直交座標航跡を生成するためのデータとして選択する。
直交座標ゲイン算出器37は、予測手段31より入力した予測誤差共分散行列Pk(−)と予め設定された受信器21−1〜21−Mの観測雑音の分散Rkとから、式(44)に従い直交座標のゲイン行列Kkを算出する。平滑器27は、直交座標ゲイン算出器37から取得した直交座標のゲイン行列Kk、予測ベクトルx^k(−)、予測誤差共分散行列Pk(−)と、相関判定器26から取得した直交座標航跡に対応する距離差un kと、距離差の予測値uk(−)とを用いて、式(45)〜(46)に従い、直交座標航跡の平滑ベクトルx^k(+)及び平滑誤差共分散行列Pk(+)を算出する。
同様に、追尾処理器22−2及び22−3においても、平滑器27は、相関判定器26から取得した直交座標航跡に対応する距離差un-1 k及びun+1 k、距離差の予測値uk(−)と、直交座標ゲイン算出器37から取得した直交座標のゲイン行列Kk、予測ベクトルx^k(−)とを用いて、式(47)〜(48)に従い、それぞれの直交座標の航跡の平滑ベクトルx^a k(+)及びx^b k(+)を算出する。
以上のようにして各位相差検出単位毎に選択された直交座標航跡に対応する距離差un k、un-1 k、un+1 kは、物理量算出器23の速度判定器38に出力される。速度判定器38は、予め設定された移動体の速度の上限Vmax、下限Vmin及びその誤差σとから、追尾処理器22−1〜22−3の平滑器27から取得した各距離差に対応する速度Vn k、Vn-1 k、Vn+1 kが式(21)及び式(22)を満たすか否かを判定する。航跡入れ替え器24は、速度判定器38において、直交座標航跡の距離差un kに対応する速度Vn kが式(21)及び式(22)を満たす場合は、直交座標航跡を現サンプリング時刻における移動体の状態とし、満たさない場合は、式式(21)及び式(22)を満たす直交座標裏航跡を現サンプリング時刻における移動体の状態とする。航跡入れ替え器24は、速度判定器38において直交座標裏航跡の距離差が正しい距離差であると判定された場合にのみ、式(23)により直交座標航跡の平滑ベクトルx^k(+)を直交座標裏航跡の平滑ベクトルに置き換えるとともに、これまで算出してきた直交座標裏航跡の平滑ベクトルを初期化する。
予測器28は、航跡入れ替え器24から取得した直交座標航跡の平滑ベクトルx^k(+)、平滑誤差共分散行列Pk(+)と、予め設定された駆動雑音の誤差共分散行列Qkとから、式(36)及び式(37)に従い、次サンプリング時刻における直交座標の予測ベクトルx^k+1(−)及び予測誤差共分散行列Pk+1(−)を算出する。相関範囲算出器30は、予測器28が算出した予測誤差共分散行列Pk+1(−)と、予め設定された受信機21の観測雑音の分散Rkとから、式(40)を用いて距離差の観測値の予測確率分布の分散Skを算出する。
以上のように、この発明の実施の形態3の位置推定装置によれば、距離差に基づいて直交座標の航跡を維持し、直交座標の平滑値から直接的に対象物1の速度を算出することとした。この結果、実施の形態1の効果に加えて、実施の形態1及び2のように測位演算を別途行わなくてもよい、という効果が生じる。
また、これにより、追尾処理器22−1〜22−3における位相差検出器25の処理が同期していなくてもよいとの利点がある。
なお、この発明の実施の形態3による位置推定装置は、航跡間の乗り移りの発生を判断基準として、対象物の速度が所定の値域(式(21)及び式(22)で示される)に含まれるかどうかを調べた。しかしながら実施の形態2で説明したように運動状態を表す物理量としては、速度に限定されないのであり、例えば位置、通過時間、加速度、角速度等を用いてもよいことはいうまでもない。
また、この発明の実施の形態3から、位相差の乗り移りを判定するために必要となる追尾処理は、位相差自体の追尾処理に限られるものではなく、位相差に基づく直交座標の追尾処理でもよく、また他の物理量に関する追尾処理であってもよいことが理解されるはずである。さらに追尾処理のアルゴリズムとしてはカルマンフィルタに限定されるものではなく、例えば上述したように拡張カルマンフィルタを用いてもよく、またα−βフィルタその他の外挿法を採用してもこの発明の特徴が失われることはないのである。
実施の形態4.
実施の形態1〜3で示した位置推定装置は、マルチパス環境下における位相差のアンビギュイティを解決するために、対象物の運動状態を表す物理量であって、位相差から導出可能な物理量の時間的な整合性、連続性に着目し、この物理量が所定の条件を満たす範囲で位相差を選択することを特徴としていた。
実施の形態1〜3で示した位置推定装置は、マルチパス環境下における位相差のアンビギュイティを解決するために、対象物の運動状態を表す物理量であって、位相差から導出可能な物理量の時間的な整合性、連続性に着目し、この物理量が所定の条件を満たす範囲で位相差を選択することを特徴としていた。
しかしこの原理は、位相差のアンビギュイティを解決するためだけでなく、方位角のアンビギュイティを解決する上でも役立つ。この発明の実施の形態4では、かかる特徴を有する位置推定装置を示す。
図9は、この発明の実施の形態4における位置推定装置の構成を示すブロック図である。図において、図1と同一の符号を付した構成要素については、実施の形態1と同様であるので説明を省略する。ただしこの位置推定装置2においては、受信器21−1〜21−Mはアレーアンテナを構成するアレー素子であるものとする。方位角算出器40は、アレーアンテナ21−1〜21−Mから得られた受信信号から送信波を送出した対象物1の方位角を算出する回路又は素子である。また、図9における追尾処理器22−1〜22−3は方位角に関する追尾処理を行うように構成されている。
図10は、位置推定装置2の追尾処理器22−1の詳細な構成を示すブロック図である。図において、方位角ゲイン算出器41は方位角のゲインを算出する回路又は素子である。他の部位については実施の形態1及び2と同様であるが、後述するように、この追尾処理器22−1は対象物1の方位角を追尾するように構成されているので、処理に供されるデータが異なっている。
なお、追尾処理器21−2及び22−3の構成については、追尾処理器21−1と同様であるので説明を省略する。またこの発明の実施4では、物理量算出器23として実施の形態1と同様に図3に示す構成を有するものとする。
ここで、対象物1の方位角の運動モデルは式(4)によって示されるものとする。なお、x kはサンプリング時刻tkにおける方位角の真値θkとその時間微分の真値θ'kを表す状態ベクトルであり、式(5)で表される。ここで、サンプリング時刻tk(k=1,2,…,n)は、方位角を検出するタイミングを表す離散時刻であり、以降、時刻tkと呼ぶ。また、Φk−1は時刻tk−1よりtkへの状態ベクトルの推移行列であり、式(6)で表される。また、wkは時刻tkにおける駆動雑音であり、Γ1(k)は時刻tkにおける駆動雑音の変換行列である。例えば、位相差の運動モデルを等速直線運動と仮定したことによる打ち切り誤差項をΓ1(k−1)wk−1とみれば、wkは位相差の加速度ベクトル相当であり、Γ1(k−1)は式(7)で表される。なお、Tはサンプリング間隔である。平均を表す記号としてEを用いると、wkは平均0、分散Qkの正規分布の白色雑音とする。
方位角の観測モデルは式(8)で表される。ここで、ukは方位角の観測値、Hは観測行列で、式(9)で表される。ν kは時刻tkにおいて検出された方位角の観測雑音である。観測雑音は平均0、分散Rkの正規分布白色雑音に従うものと仮定する。また、時刻tkまでの間に方位角の追尾に用いた方位角の全体をUkとする。
以上の運動モデルに基づき、カルマンフィルタを用いた場合を例に追尾フィルタを構成する。1サンプリング前の方位角航跡の平滑ベクトルをx^k−1(+)、平滑誤差共分散行列をPk−1(+)とすると、時刻tkにおける方位角の状態ベクトルx kの予測ベクトルx^k(−)及び予測誤差共分散行列Pk(−)はそれぞれ条件付平均ベクトル、条件付共分散行列で定義され、式(10)及び式(11)で表される。
次に、時刻tkに計測された複数の方位角Ukの選択方法について述べる。偽の方位角の含まれる複数の方位角から正しい方位角を判定するため、式(12)で表される方位角の予測値u k(−)を利用する。
次に、時刻tkに計測された複数の方位角Ukの選択方法について述べる。偽の方位角の含まれる複数の方位角から正しい方位角を判定するため、式(12)で表される方位角の予測値u k(−)を利用する。
過去の方位角情報Uk-1より得られる方位角の観測値の予測確率分布P[uk|Uk-1]は、式(13)に示す条件付確率密度関数で表される。すなわち、方位角の観測値は、式(11)で与えられる方位角の予測値uk(−)を平均とし、式(14)で与えられるSkを分散とする1次元正規分布g(uk;uk(−),Sk)に従うものとする。そして、式(15)を満たす方位角のうち、fiが最小となる方位角un kを正しい方位角と判断し、選択する。また、選択された方位角の前後の整数値バイアス値に対応する方位角un-1 k及びun+1 kを方位角裏航跡を生成するためのデータとする。なお、dは目標との相関範囲を決めるパラメータであり、自由度1のχ自乗分布により算出する。
次に、方位角航跡と方位角裏航跡に対応する平滑ベクトルの算出方法について述べる。ゲイン行列Kk、方位角航跡の平滑ベクトルx^k(+)及び平滑誤差共分散行列Pk(+)はカルマンフィルタの理論により、式(16)〜式(18)で与えられる。同様に、正しい方位角の前後の整数値バイアス値に対応する方位角un-1 k及びun+1 kをそれぞれ用いた方位角裏航跡の平滑ベクトルをx^a k(+)及びx^b k(+)とすると、これらは式(19)及び式(20)で与えられる。
以上のように、方位角についても実施の形態1〜3で扱った位相差と全く同じ追尾処理を適用できることが分かる。そしてここで求めた方位角平滑値を複数組み合わせれば、対象物1の位置を推定できる。またその時間変位を算出すれば対象物1の速度を求めることもできる。
このことから各整数値バイアスをとる方位角について追尾処理を行い、その平滑値から対象物1の運動状態を示す物理量を算出して、式(21)及び式(22)、あるいは式(27)〜式(29)を最小にする方位角を算出することで、方位角のアンビギュイティを抑制し、対象物1の位置を精度よく推定することができるのである。
なお、アレーアンテナを用いて方位角を算出する方位角算出器40の処理については公知の方位角算出方法を用いて構成すればよい。このような方位角の算出方法としては、例えば、隣接するビーム間の受信信号の振幅や位相などの差(Δ値)と和(Σ値)の比(Δ/Σ値)を求め、このΔ/Σ値から反射波の入射方向を求める方法などがある。
以上から明らかなように、位相差に替えて方位角から対象物の位置を測定する機器においても、この発明の原理を適用することが可能である。
以上のように、この発明に係る位置推定装置は、例えば、マルチパスによる干渉波が発生する環境下において対象物から送信される電波を受信して、その対象物の位置を推定するシステムに適用することが可能である。
22−1〜22−3 追尾処理器、
23 物理量演算器、
24 航跡入れ替え器。
23 物理量演算器、
24 航跡入れ替え器。
Claims (8)
- 対象物からの電波を複数の受信器により受信して得た受信波間の位相差が整数値バイアスの異なる複数の受信波位相差を含む場合に、これらの位相差に基づく追尾処理を行って上記対象物の位置推定に用いる平滑値を算出し、この平滑値から上記対象物の位置を推定する位置推定装置において、
上記複数の受信波位相差から第1位相差を選択しこの第1位相差に基づく航跡を維持して平滑値を算出する第1追尾処理器と、
上記第1位相差に隣接する第2位相差を選択してその第2位相差に基づく航跡を維持する第2追尾処理器と、
上記第1追尾処理器が算出した平滑値から位相差とは異なる物理量であって上記対象物の運動状態を表す物理量を算出する物理量算出器と、
上記物理量算出器により算出された物理量の連続性に基づいて上記第1位相差に基づく航跡と上記第2位相差に基づく航跡とを入れ替える航跡入れ替え器と、
を備えたことを特徴とする位置推定装置。 - 上記物理量算出器は、上記平滑値から上記対象物の位置を推定するとともに、この位置の時間変位から上記対象物の速度を算出し、
上記航跡入れ替え器は、上記物理量算出器により算出された速度が所定の値域外にある場合に上記第1位相差に基づく航跡と上記第2位相差に基づく航跡とを入れ替えることを特徴とする請求項1に記載の位置推定装置。 - 上記物理量算出器は、上記平滑値から上記対象物の位置を推定するとともに、この位置の時間変位から上記対象物が基準地点を通過する予測時刻を算出し、
上記航跡入れ替え器は、上記物理量算出器により算出された予測時刻が所定の地域外にある場合に上記第1位相差に基づく航跡と上記第2位相差に基づく航跡とを入れ替えることを特徴とする請求項1に記載の位置推定装置。 - 上記第1及び第2追尾処理器は、上記第1及び第2位相差に基づく航跡としてそれらの位相差に基づいて算出される上記対象物の直交座標の航跡をそれぞれ維持することを特徴とする請求項1に記載の位置推定装置。
- 対象物からの電波を複数の受信器により受信して得た受信波から整数値バイアスの異なる複数の方位角が得られる場合に、これらの方位角に基づく追尾処理を行って上記対象物の位置推定に用いる平滑値を算出し、この平滑値から上記対象物の位置を推定する位置推定装置において、
上記複数の方位角から第1方位角を選択しこの第1方位角に基づく航跡を維持して上記平滑値を算出する第1追尾処理器と、
上記第1方位角に隣接する第2方位角を選択してその第2方位角に基づく航跡を維持する第2追尾処理器と、
上記平滑値から方位角とは異なる物理量であって上記対象物の運動状態を表す物理量を算出する物理量算出器と、
上記物理量算出器により算出された物理量の連続性に基づいて上記第1方位角に基づく航跡と上記第2方位角に基づく航跡とを入れ替える航跡入れ替え器と、
を備えたことを特徴とする位置推定装置。 - 上記物理量算出器は、上記平滑値から上記対象物の位置を推定するとともに、この位置の時間変位から上記対象物の速度を算出し、
上記航跡入れ替え器は、上記物理量算出器により算出された速度が所定の値域外にある場合に上記第1方位角に基づく航跡と上記第2方位角に基づく航跡とを入れ替えることを特徴とする請求項5に記載の位置推定装置。 - 上記物理量算出器は、上記平滑値から上記対象物の位置を推定するとともに、この位置の時間変位から上記対象物が基準地点を通過する予測時刻を算出し、
上記航跡入れ替え器は、上記物理量算出器により算出された予測時刻が所定の地域外にある場合に上記第1方位角に基づく航跡と上記第2方位角に基づく航跡とを入れ替えることを特徴とする請求項5に記載の位置推定装置。 - 上記第1及び第2追尾処理器は、上記第1及び第2方位角に基づく航跡としてそれらの方位角に基づいて算出される上記対象物の直交座標の航跡をそれぞれ維持することを特徴とする請求項5に記載の位置推定装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004122505A JP2005308428A (ja) | 2004-04-19 | 2004-04-19 | 位置推定装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004122505A JP2005308428A (ja) | 2004-04-19 | 2004-04-19 | 位置推定装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2005308428A true JP2005308428A (ja) | 2005-11-04 |
Family
ID=35437385
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004122505A Withdrawn JP2005308428A (ja) | 2004-04-19 | 2004-04-19 | 位置推定装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2005308428A (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009264968A (ja) * | 2008-04-25 | 2009-11-12 | Fujitsu Ten Ltd | レーダ装置、レーダ装置の信号処理方法及び、車両制御システム |
JP2012002809A (ja) * | 2010-06-17 | 2012-01-05 | Honeywell Internatl Inc | 干渉レーダーの曖昧さを除去するために証拠グリッドを使用するためのシステムと方法 |
JP2017067610A (ja) * | 2015-09-30 | 2017-04-06 | 三菱電機株式会社 | 信号処理装置、合成開口レーダ装置および信号処理方法 |
JP2021509717A (ja) * | 2018-03-28 | 2021-04-01 | 三菱電機株式会社 | 測位システム、方法及び媒体 |
JP2021524590A (ja) * | 2018-07-10 | 2021-09-13 | フリードリヒ−アレクサンダー−ウニベルジテート・エアランゲン−ニュルンベルク | 波動に基づいた信号を用いて少なくとも1つの物体の位置特定を行うための位置特定方法及び位置特定システム |
CN116753961A (zh) * | 2023-08-16 | 2023-09-15 | 中国船舶集团有限公司第七〇七研究所 | 基于状态观测的动力定位船舶高速循迹导航方法 |
-
2004
- 2004-04-19 JP JP2004122505A patent/JP2005308428A/ja not_active Withdrawn
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009264968A (ja) * | 2008-04-25 | 2009-11-12 | Fujitsu Ten Ltd | レーダ装置、レーダ装置の信号処理方法及び、車両制御システム |
JP2012002809A (ja) * | 2010-06-17 | 2012-01-05 | Honeywell Internatl Inc | 干渉レーダーの曖昧さを除去するために証拠グリッドを使用するためのシステムと方法 |
JP2017067610A (ja) * | 2015-09-30 | 2017-04-06 | 三菱電機株式会社 | 信号処理装置、合成開口レーダ装置および信号処理方法 |
JP2021509717A (ja) * | 2018-03-28 | 2021-04-01 | 三菱電機株式会社 | 測位システム、方法及び媒体 |
JP2021524590A (ja) * | 2018-07-10 | 2021-09-13 | フリードリヒ−アレクサンダー−ウニベルジテート・エアランゲン−ニュルンベルク | 波動に基づいた信号を用いて少なくとも1つの物体の位置特定を行うための位置特定方法及び位置特定システム |
JP7326422B2 (ja) | 2018-07-10 | 2023-08-15 | フリードリヒ-アレクサンダー-ウニベルジテート・エアランゲン-ニュルンベルク | 波動に基づいた信号を用いて少なくとも1つの物体の位置特定を行うための位置特定方法及び位置特定システム |
CN116753961A (zh) * | 2023-08-16 | 2023-09-15 | 中国船舶集团有限公司第七〇七研究所 | 基于状态观测的动力定位船舶高速循迹导航方法 |
CN116753961B (zh) * | 2023-08-16 | 2023-10-31 | 中国船舶集团有限公司第七〇七研究所 | 基于状态观测的动力定位船舶高速循迹导航方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10175348B2 (en) | Use of range-rate measurements in a fusion tracking system via projections | |
EP3688495B1 (en) | System for determining a physical metric such as position | |
DK3123197T3 (en) | METHODS AND APPARATUS FOR DETERMINING THE ARRANGEMENT (AOA) OF A RADAR WARNING RECEIVER | |
JP5972259B2 (ja) | レーダ装置 | |
CN110088642B (zh) | 用于检测目标对象的前进方向和速度的***和方法 | |
US10429500B2 (en) | Tracking apparatus, tracking method, and computer-readable storage medium | |
US10567918B2 (en) | Radio-location method for locating a target device contained within a region of space | |
US11486966B2 (en) | Method and apparatus for tracking target from radar signal using artificial intelligence | |
US9075142B2 (en) | Device and method for selecting signal, and radar apparatus | |
JP4116898B2 (ja) | 目標追尾装置 | |
JPH11118903A (ja) | 位置検出装置 | |
JP2021524590A (ja) | 波動に基づいた信号を用いて少なくとも1つの物体の位置特定を行うための位置特定方法及び位置特定システム | |
JP2011247596A (ja) | レーダ装置 | |
JP6095899B2 (ja) | 物標運動推定装置、物標運動推定方法、およびレーダ装置 | |
JP4766604B2 (ja) | 電波発射源可視化装置 | |
JP2005308428A (ja) | 位置推定装置 | |
JP5163765B2 (ja) | 測角装置、レーダ装置、測角方法及び測角プログラム | |
JP2008304329A (ja) | 測定装置 | |
JP2010237087A (ja) | レーダ装置及びそれを用いた電波到来方向の計測方法 | |
EP4119973A1 (en) | Positioning in a mobile communication system | |
JPH10246778A (ja) | 目標探知追尾装置 | |
JP5208408B2 (ja) | 相対位置推定システム | |
KR20140088683A (ko) | 차량용 레이더를 이용한 물체 감지 장치, 방법 및 컴퓨터 판독 가능한 기록 매체 | |
KR101052050B1 (ko) | 다수의 레이더를 가간섭적으로 결합하여 이동 표적의 감지율을 향상시키는 방법 및 그 방법이 적용된 레이더 시스템 | |
JP5599200B2 (ja) | 目標検出・追尾装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070129 |
|
A761 | Written withdrawal of application |
Effective date: 20100126 Free format text: JAPANESE INTERMEDIATE CODE: A761 |