JPWO2006067874A1 - 測位装置および測位装置の位置特定方法 - Google Patents

Abstract

未知電波を用いる場合にも、各基地局での直接波受信時刻の計測に基づく電波源の位置の時間差測位を行うことができる測位装置を得る。電波源１から発信された電波を受信し、参照信号と受信電波との相関から電波受信時刻を計測する複数の受信時刻算出基地局２と、複数の受信時刻算出基地局２からのそれぞれの電波受信時刻に基づいて電波源１の位置を特定する時間差測位器３とを備えた測位装置において、電波源１から発信された未知電波を推定して基準時刻における参照信号を生成する特定基地局２（ｋ0）をさらに備え、複数の受信時刻算出基地局２は、特定基地局２（ｋ0）からの参照信号と受信した未知電波との相関から電波受信時刻を計測する。また、電波源としてＲＦＩＤタグ、特定基地局および受信時刻算出基地局としてＲＦＩＤリーダ／ライタを使用できる。

Description

本発明は、ＲＦＩＤ（Radio Frequency IDentification）タグ等の電波源から発射された電波をＲＦＩＤリーダ等の複数の基地局で受信し、各基地局での電波受信時刻差から電波源の位置を特定する時間差測位方式に関するものである。

複数の基地局により受信された直接波の受信時刻から、その直接波の送信元である電波源の位置を求めることのできる時間差測位方法がある（例えば、特許文献１参照）。基準となる基地局における直接波受信時刻と、各基地局自身における直接波受信時刻との差から、それぞれの基地局から電波の送信元である電波源までの距離を求め、複数の基地局からの距離が一致する場所を特定することにより電波源の位置を推定するものである。

特開２００４−５３２３１４号公報（第１頁、図３）

しかしながら、従来技術には次のような課題がある。従来の方式では、各基地局で受信される電波が既知であることが前提であり、未知電波を用いる場合には時間計測を行うことができないという問題がある。

本発明は上述のような課題を解決するためになされたもので、未知電波を用いる場合にも、各基地局での直接波受信時刻の計測に基づく電波源の位置の時間差測位を行うことができる測位装置および測位装置の位置特定方法を得ることを目的とする。

本発明に係る測位装置は、電波源から発信された電波を受信し、参照信号と受信電波との相関から電波受信時刻を計測する複数の受信時刻算出基地局と、複数の受信時刻算出基地局からのそれぞれの電波受信時刻に基づいて電波源の位置を特定する時間差測位器とを備えた測位装置において、電波源から発信された未知電波を推定して基準時刻における参照信号を生成する特定基地局をさらに備え、複数の受信時刻算出基地局は、特定基地局からの参照信号と受信した未知電波との相関から電波受信時刻を計測するものである。

本発明に係る測位装置の位置特定方法は、特定基地局で、電波源から発信された未知電波を推定して基準時刻における参照信号を生成するステップと、複数の受信時刻算出基地局で、特定基地局からの参照信号と受信した未知電波との相関から電波受信時刻を計測するステップと、時間差測位器で、複数の受信時刻算出基地局からのそれぞれの電波受信時刻に基づいて電波源の位置を特定するステップとを備えたものである。

さらに、本発明に係る測位装置の位置特定方法は、複数のＲＦＩＤリーダ／ライタあるいは特定ＲＦＩＤリーダ／ライタの何れかからＲＦＩＤタグに電波を送信するステップと、ＲＦＩＤタグで、電波の受信電力を検知することにより未知電波を発信するステップと、特定ＲＦＩＤリーダ／ライタで、ＲＦＩＤタグから発信された未知電波を推定して基準時刻における参照信号を生成するステップと、複数のＲＦＩＤリーダ／ライタで、特定ＲＦＩＤリーダ／ライタからの参照信号と受信した未知電波との相関から電波受信時刻を計測するステップと、時間差測位器で、複数の複数のＲＦＩＤリーダ／ライタからのそれぞれの電波受信時刻に基づいてＲＦＩＤタグの位置を特定するステップとを備えたものである。

本発明によれば、特定基地局において未知電波の推定を行い、それぞれの基地局は特定基地局による推定結果に基づいて直接波受信時刻を計測することにより、未知電波を用いる場合にも、各基地局での直接波受信時刻の計測に基づく電波源の位置の時間差測位を可能とした測位装置および測位装置の位置特定方法を得ることができる。

本発明の実施の形態１における測位装置の構成図である。 本発明の実施の形態１の特定基地局における未知電波の推定を行う構成を示した図である。 本発明の実施の形態１における参照信号生成器の構成図である。 本発明の実施の形態１における未知電波の推定結果を用いて直接波受信時刻を計測する基地局の構成図である。 本発明の実施の形態２における未知電波の推定結果を用いて直接波の受信時刻を計測する基地局の構成図である。 本発明の実施の形態２におけるフィルタバンクの構成図である。 本発明の実施の形態３において時間差測位器が電波源の位置を特定するために行う処理手順を示したフローチャートである。 本発明の実施の形態５の特定基地局における未知電波の推定を行う構成を示した図である。 本発明の実施の形態５における複素振幅比対応型参照信号生成器の構成図である。 本発明の実施の形態６において車両位置の測位に位置特定方法を適用した応用例の説明図である。 本発明の実施の形態６の応用例における時間差測位の説明図である。 本発明の実施の形態６の応用例における車両に搭載されたＥＴＣ情報送信アンテナを含むＥＴＣの構成図である。 本発明の実施の形態７において車両位置の測位に位置特定方法を適用した応用例の説明図である。 本発明の実施の形態７の応用例における車両に搭載されたスマートプレート送信アンテナを含むスマートプレートの構成図である。 時間差測位により電波源の位置を特定する測位装置の基本構成を示した図である。 直接波受信時刻の計測を行うための基本構成を示す図である。 ｋ₀番目の基地局から電波源までの距離をｒ_k0と仮定したときに特定される各基地局から電波源までの距離の関係を示す図である。 ある値のｒ_k0において、それぞれの円周がある１点で交わった状態を示す図である。

以下、本発明の測位装置の好適な実施の形態につき図面を用いて説明する。
本発明の測位装置は、未知電波を用いて電波源の位置を測位するものであり、特定基地局により未知電波の推定を行って基準時刻の参照信号を生成し、この参照信号に基づいてそれぞれの基地局が直接波受信時刻を求めることを特徴としている。

実施の形態１．
まず始めに、時間差測位の基本原理について説明する。図１５は、時間差測位により電波源の位置を特定する測位装置の基本構成を示した図であり、電波源１、Ｋ個の基地局１５（１）〜１５（Ｋ）（ただし、Ｋは３以上の整数）、および時間差測位器３で構成される。図１５は、３個の基地局を示しており、Ｋ＝３の場合に相当する。

電波源１は、各基地局１５（１）〜１５（Ｋ）に対して電波を発信する。これに対して、各基地局１５（１）〜１５（Ｋ）は、電波源１から発信されそれぞれの基地局に直接到達する直接波と、マルチパスにより直接波と異なる経路を経由して到達する間接波とを含む電波を受信し、その受信した電波から直接波受信時刻を計測する。

時間差測位器３は、それぞれの基地局１５（１）〜１５（Ｋ）により計測された直接波受信時刻を受け取り、それぞれの基地局の位置および直接波受信時刻を基に電波源１の位置を時間差測位処理により特定する。

また、図１６は、直接波受信時刻の計測を行うための基地局の基本構成を示す図であり、中原他“ＦＦＴ−ＭＵＳＩＣ法とＦＦＴ演算型相関法の多重伝播遅延時間分解能の比較検討”、信学総大、Ｂ−２７、１９９５に示されたものである。それぞれの基地局１５（１）〜１５（Ｋ）は、図１６に示す基本構成を有しており、自らが受信した電波に基づいて直接波受信時刻を求めることができる。

この基地局は、パッシブアンテナ４、受信機５、Ａ／Ｄ変換器６、ＦＦＴ処理器７ａ、７ｂ、参照信号生成器１６、相関処理器１２、ＭＵＳＩＣ処理器１３で構成される。ここで、ＦＦＴ処理器７ａ、７ｂは周波数成分算出手段に、相関処理器１２は相関信号生成手段に、ＭＵＳＩＣ処理器１３は電波受信時刻推定手段にそれぞれ相当している。

電波源１から発信された電波は、パッシブアンテナ４により受信される。そして、パッシブアンテナ４で受信された電波は、受信機５により帯域制限、位相検波等が行われた後に、Ａ／Ｄ変換器５によりアナログ信号からディジタル信号に変換される。

このようにして得られたディジタル信号は、さらに、ＦＦＴ処理器７ａによりＦＦＴ（Fast Fourier transform）処理が行われ、時間領域の信号から周波数領域の信号に変換される。一方、参照信号生成器１６で生成された参照信号も、同様に、ＦＦＴ処理器７ｂによりＦＦＴ（Fast Fourier transform）処理が行われ、時間領域の信号から周波数領域の信号に変換される。

ＦＦＴ処理器７ａ、７ｂから出力されるそれぞれの周波数領域の信号は、相関処理器１２により相関演算が行われる。さらに、相関演算後の信号は、ＭＵＳＩＣ処理器１３によりＭＵＳＩＣ（Multiple Signal classification）処理が施され、高分解能に計測した直接波受信時刻が求められる。このようにして計測された直接波受信時刻は、それぞれの基地局１５（１）〜１５（Ｋ）から時間差測位器３に伝達され、時間差測位器３は、それぞれの直接波受信時刻に基づいて電波源の位置を特定することとなる。

次に、測位装置の一連の動作、および測位のための具体的な演算式について説明する。それぞれの基地局１５（１）〜１５（Ｋ）は、電波源１から発信された電波が直接到達する直接波と、マルチパスにより直接波と異なる経路を経由して到達する間接波とを含む電波を、パッシブアンテナ４で受信する。受信時刻は、全ての基地局共通のあらかじめ設定された基準時刻を零として計られる。

ここで、あるｋ番目（ただし、ｋは１≦ｋ≦Ｋの整数）の基地局を例にして、その処理内容を説明する。ｋ番目の基地局が直接波を受信した時刻をτ_1,kとし、間接波を受信した時刻をτ_2,kとする。パッシブアンテナ４で受信された電波は、受信機５に伝達され、帯域制限、位相検波が行われる。さらに、受信機５からの出力信号は、Ａ／Ｄ変換器６に伝達される。

Ａ／Ｄ変換器６は、受信電波をサンプルしてディジタル信号に変換する。ここで、Ａ／Ｄ変換器６において、サンプル周期をＴ_samp、サンプル数をＮ（ただし、Ｎは１以上の整数）としてディジタル化された受信信号をｚ_1,k・・・ｚ_N,kとする。次に、ＦＦＴ処理器７ａは、ディジタル化された受信信号ｚ_1,k・・・ｚ_N,kにＦＦＴ処理を施して受信信号の周波数成分ｙ_1,k・・・ｙ_N,kを算出する。

ここで、受信信号の周波数成分ｙ_1,k・・・ｙ_N,kに含まれる電波源１から発信された直接波成分と間接波成分とを加えた周波数成分ｙ'_1,k・・・ｙ'_N,kは、次式（１）、（２）により表される。この式（１）、（２）において、Γ_i,iは参照信号（詳細は後述）の周波数成分を表している。また、ｃ_1,kは直接波の振幅値、ｃ_2,kは間接波の振幅値をそれぞれ表している。

参照信号生成器１６は、各基地局で受信される電波が既知であることを前提として、電波源１から発信された電波が時刻零で受信されたことを想定して参照信号ｚ０₁,・・・,ｚ０_Nを生成する。参照信号ｚ０₁,・・・,ｚ０_Nは、ＦＦＴ処理器７ｂによりＦＦＴ処理が施され、参照信号の周波数成分Γ_1,1,・・・,Γ_N,Nが生成される。相関処理器１２は、次式（３）を用いて、参照信号の周波数成分Γ_1,1,・・・,Γ_N,Nと受信信号の周波数成分ｙ_1,k ・・・ ｙ_N,kとの相関演算を行い、周波数領域における相関信号ｘ_1,k・・・ｘ_N,kを生成する。

周波数領域における相関信号ｘ_i,kに含まれる直接波成分と間接波成分とを加えた相関信号ｘ'_i,kは、具体的には次式（４）により表される。

実際に相関処理器１２から出力される相関信号ｘ_i,kは、直接波成分と間接波成分とを加えた相関信号ｘ'_i,kを用いて、次式（５）のように表される。式（５）において、ｎ_i,kは、ＭＵＳＩＣ処理入力時において式（４）の相関信号ｘ'_i,kに重畳している受信機雑音成分を表している。

ＭＵＳＩＣ処理器１３は、相関信号ｘ_i,k（ただし、ｉは１≦ｉ≦Ｎとなる整数）を入力信号として取り込む。ベクトルｘ_iをｘ_i＝［ｘ_i,k,ｘ_i+1,k・・・ｘ_i+M-1,k］^Tとして定義するとき、ベクトルｘ_iは、次式（６）のように表される。ＭＵＳＩＣ処理器１３は、ａ（τ）をステアリングベクトルとしてＭＵＳＩＣ処理を行い、直接波の観測時刻の計測値に相当する直接波受信時刻の推定値τ_1,k（チルダ）を算出し、時間差測位器３に伝達する。なお、ここで用いたＭは、ステアリングベクトルの次数であり、サンプル数Ｎよりも小さい整数である。

ここで、ＭＵＳＩＣ処理を行った際に、直接波と間接波の２種類の受信時刻が計測されることとなるが、直接波は、間接波に比較して速く到達するので、小さい受信時刻の方を直接波の到達時刻として計測値τ（チルダ）を推定する。

時間差測位器３は、各基地局から伝達されたτ（チルダ）に基づき、電波源１の位置を時間差測位する。時間差測位器３は、特定基地局（ｋ₀番目の基地局）を設定し、そのｋ₀番目の基地局から電波源１までの距離をｒ_k0として仮定する。このとき、時間差測位器３は、ｋ₀番目の基地局における直接波受信時刻と、他の基地局における直接波受信時刻との差から、次式（１０）により、一意的に、ｋ₀番目の基地局以外の各基地局から電波源１までの距離を特定できる。ここで、ｃは、光速を表している。

図１７は、ｋ₀番目の基地局から電波源１までの距離をｒ_k0と仮定したときに特定される各基地局から電波源１までの距離の関係を示す図である。ｒ_k0の仮定のもとに、あるｋ番目の基地局の位置を中心として半径ｒ_kの円をそれぞれ描いた状態を示している。もし、距離ｒ_kの位置に電波源１が存在するとき、半径ｒ_kの円周上のいずれかの地点に電波源１が存在することとなる。

そこで、最初に仮定したｒ_k0を変えていき、そのときのそれぞれの半径ｒ_kを求めると、ｒ_k0がある値となったときに、それぞれの円周がある１点で交わる状態が発生することとなる。図１８は、ある値のｒ_k0において、それぞれの円周がある１点で交わった状態を示す図である。従って、時間差測位器３は、各基地局における直接波受信時刻に基づいて、図１８のような状態となるｒ_k0を求めることにより、電波源１の位置を測位することができることとなる。

実際には、直接波受信時刻の計測誤差等により、各基地局の位置を中心とするそれぞれの円周は、ある１点で完全には交わらない。そこで、時間差測位器３は、最小二乗法等を用いて距離の二乗誤差が最小となるような位置を算出することにより、電波源１の位置を推定することとなる。

上述した方式では、各基地局で受信される電波が既知であることを前提とし、既知電波に基づいて参照信号を生成し、直接波受信時刻を計測している。これに対して、本発明は、未知電波を受信した際にも時間計測に基づく測位を行うことができる点を特徴としており、その詳細について、次に説明する。

図１は、本発明の実施の形態１における測位装置の構成図であり、電波源１、Ｋ個の未知電波対応型基地局２（１）〜２（Ｋ）、および時間差測位器３で構成される。電波源１および時間差測位器３は、図１５の基本構成で示したものと同一であり、未知電波対応型基地局２（１）〜２（Ｋ）を中心に説明する。

未知電波対応型基地局２（ｋ₀）（ただし、ｋ₀は１≦ｋ₀≦Ｋの整数）は、未知電波を推定して、その推定結果として基準時刻の参照信号を生成して他の未知電波対応型基地局２（１）、２（２）、２（Ｋ）に対して送信する特定基地局である。

一方、未知電波対応型基地局２（１）、２（２）、２（Ｋ）は、未知電波対応型基地局２（ｋ₀）からの参照信号を用いて、電波源１から発信された直接波の受信時刻を計測する受信時刻算出基地局である。

このように、本発明の測位装置は、ある特定基地局（未知電波対応型基地局２（ｋ₀）に相当）により未知電波の推定を行って基準時刻の参照信号を生成するとともに、他の基地局（未知電波対応型基地局２（１）、２（２）、２（Ｋ）に相当）において、この参照信号に基づいて直接波受信時刻を計測し、電波源の測位を行う。

図２は、本発明の実施の形態１の特定基地局における未知電波の推定を行う構成を示した図であり、特定基地局である未知電波対応型基地局２（ｋ₀）の構成を示すものである。図２におけるパッシブアンテナ４、受信機５、Ａ／Ｄ変換器６、ＦＦＴ処理器７は、図１６で示したものと同一である。

参照信号生成手段である参照信号生成器８は、ＦＦＴ処理された未知電波の周波数成分の信号を推定し、直接波受信時刻を算出するための参照信号を生成する。図３は、本発明の実施の形態１における参照信号生成器８の構成図である。この参照信号生成器８は、間接波受信時刻推定器９、振幅比推定器１０、初期位相推定器１１で構成される。

間接波受信時刻推定器９は、ＦＦＴ処理された未知電波の周波数成分の信号から間接波の受信時刻を推定する。振幅比推定器１０は、時刻ゼロで受信される直接波の振幅と、推定された間接波の受信時刻における間接波の振幅との振幅比を推定する。さらに、初期位相推定器１１は、電波源から発信された電波の初期位相を推定する。

図４は、本発明の実施の形態１における未知電波の推定結果を用いて直接波受信時刻を計測する基地局の構成図であり、受信時刻算出基地局である未知電波対応型基地局２（１）、２（２）、２（Ｋ）の構成を示すものである。図４におけるパッシブアンテナ４、受信機５、Ａ／Ｄ変換器６、ＦＦＴ処理器７、相関処理器１２、ＭＵＳＩＣ処理器１３は、図１６で示したものと同一である。

図４の構成を有する受信時刻算出基地局は、図１６で示した基地局の構成と比較すると、参照信号生成器１６およびＦＦＴ処理器７ｂが不要となっており、その代わりに参照信号を特定基地局から取り込む構成となっている。

次に、一連の動作、および測位のための具体的な演算式について説明する。それぞれの未知電波対応型基地局２（１）、２（２）、２（Ｋ）、２（ｋ₀）は、電波源１から発信された電波が直接到達する直接波と、マルチパスにより直接波と異なる経路を経由して到達する間接波とを含む電波をパッシブアンテナ４で受信する。

まず始めに、特定基地局である未知電波対応型基地局２（ｋ₀）は、受信した電波に基づいて、直接波の受信時刻をゼロとした基準時刻における参照信号を生成する。その後、受信時刻算出基地局である未知電波対応型基地局２（１）、２（２）、２（Ｋ）は、未知電波対応型基地局２（ｋ₀）で生成された参照信号と、自らが生成したＦＦＴ処理後の受信信号の周波数成分との相関演算を行うことにより、直接波受信時刻の計測を行う。

未知電波対応型基地局２（ｋ₀）は、先の図２に示したような未知電波の推定を行うための構成を有している。電波源１から発信された電波は、マルチパスを発生する。そして、パッシブアンテナ４は、直接波と間接波とを含む電波を受信することとなる。その後、受信機５、Ａ／Ｄ変換器６、ＦＦＴ処理器７での処理は、図１６で説明したものと同一であり、ＦＦＴ処理器７から受信信号の周波数成分ｙ_1,k0 ・・・ｙ_N,k0が出力される。

受信信号の周波数成分ｙ_1,k0 ・・・ｙ_N,k0は、参照信号生成器８に伝達される。ここで、参照信号生成器８は、図３に示す処理構成となっている。参照信号生成器８に入力した受信信号の周波数成分ｙ_1,k0 ・・・ｙ_N,k0は、間接波受信時刻推定器９に伝達される。間接波受信時刻推定器９は、最初に、次式（１１）（１２）により受信信号周波数成分の電力値Ｐｙ（ｉ）を算出する。

さらに、次式（１３）により信号ベクトルＰｙ_i,k0を設定する。

このとき、信号ベクトルＰｙ_i,k0は、次式（１４）〜（１６）により表すことができる。

さらに、間接波受信時刻推定器９は、次式（１７）により信号ベクトルＰｙ_i,k0の相関行列Ｒｐ_k0を算出し、ＭＵＳＩＣ処理により、直接波の受信時刻をゼロとしたときの間接波の受信時刻の推定値τ_2,k0（チルダ）を高分解能で算出する。

間接波受信時刻推定器９から出力される受信時刻の推定値τ_2,k0（チルダ）は、振幅比推定器１０に伝達される。振幅比推定器１０は、受信時刻の推定値τ_2,k0（チルダ）に基づいて直接波と間接波との振幅比を推定する。最初に、次式（１８）により振幅比Ｃ_k0（ｉ）の推定値Ｃ_k0（ｉ）（チルダ）を算出する。

次に、Ｃ_k0（ｉ）（チルダ）を平均して、次式（１９）により直接波と間接波との振幅比の推定値Ｃ_k0（チルダ）を算出する。

振幅比推定器１０は、間接波受信時刻の推定値τ_2,k0（チルダ）、および直接波と間接波との振幅比の推定値Ｃ_k0（チルダ）を初期位相推定器１１に伝達する。次に、初期位相推定器１１は、次式（２０）で表される位相検出用信号Ｈ_iをまず始めに算出する。

次に、初期位相推定器１１は、上式（２０）で算出した位相検出用信号Ｈ_iの実部と虚部から参照信号の周波数成分Γ_i,iの初期位相の推定値φ_i（チルダ）を算出する。さらに、初期位相推定器１１は、次式（２１）により、参照信号の周波数成分の推定値Γ_i,i（チルダ）を算出する。

初期位相推定器１１は、このようにして算出された参照信号の周波数成分の推定値Γ_i,i（チルダ）を各未知電波対応型基地局２（１）、２（２）、２（Ｋ）に伝達する。各未知電波対応型基地局２（１）、２（２）、２（Ｋ）は、式（３）における参照信号の周波数成分Γ_i,iに換えて、未知電波対応型基地局２（ｋ₀）から送られてきた参照信号の周波数成分の推定値Γ_i,i（チルダ）を用いて周波数領域での相関信号ｘ_i,kを算出する。

さらに、ＭＵＳＩＣ処理器１３は、各基地局で受信される電波が既知であることを前提としてすでに図１６を用いて説明したものと同一の手順により、直接波の観測時刻の計測値として直接波受信時刻の推定値τ_1,k（チルダ）を算出する。

なお、電波源１から発信する未知電波として、周波数成分の絶対値が周波数に依らずに一定となるものを採用することにより、より高精度な測位が可能となる。すなわち、特定基地局は、周波数領域の信号の振幅（絶対値）が一定のものを用いれば、未知パラメータが位相だけとなるため、初期位相の検出を確実に行うことができる。

また、直接波と間接波の受信時刻の差に応じて高分解能の処理が必要な場合には、上述したようなＭＵＳＩＣ処理が用いられる。この場合にも、周波数成分の絶対値が周波数に依らずに一定となる電波を未知電波として採用することにより、より高精度な測位が可能となる。

実施の形態１によれば、特定基地局は、受信した未知電波に基づいて直接波の受信時刻をゼロとした基準時刻における参照信号を生成することができる。さらに、受信時刻算出基地局は、特定基地局により生成された参照信号に基づいて、自らが受信した未知電波から直接波受信時刻を計測することができる。これにより、未知電波を用いた場合にも直接波受信時刻の計測が可能となり、電波源の位置の時間差測位を行うことができる。

実施の形態２．
本実施の形態２では、未知電波の推定結果を用いて直接波の受信時刻を計測する別の方法を説明する。図５は、本発明の実施の形態２における未知電波の推定結果を用いて直接波の受信時刻を計測する基地局の構成図であり、受信時刻算出基地局である未知電波対応型基地局２（１）、２（２）、２（Ｋ）の構成を示すものである。

図５におけるパッシブアンテナ４、受信機５、Ａ／Ｄ変換器６、相関処理器１２、ＭＵＳＩＣ処理器１３は、実施の形態１における図４と同じである。本実施の形態２では、ＦＦＴ処理器７の代わりにフィルタバンク１４を用いており、このフィルタバンク１４の機能を中心に説明する。

図６は、本発明の実施の形態２におけるフィルタバンク１４の構成図である。フィルタバンク１４は、Ｎ個（ただし、Ｎは２以上の整数であり、サンプル数Ｎと同じ数）のＦＩＲ（Finite Impulse Response）型帯域通過フィルタ１４（１）〜１４（Ｎ）で構成される。帯域通過フィルタ１４（１）〜１４（Ｎ）は、入力電波の周波数成分の存在範囲Ｗｐに対して、それぞれ異なる通過帯域を有している。

次に、一連の動作、および測位のための具体的な演算式について説明する。それぞれの未知電波対応型基地局２（１）、２（２）、２（Ｋ）、２（ｋ₀）は、電波源１から発信された電波が直接到達する直接波と、マルチパスにより直接波と異なる経路を経由して到達する間接波とを含む電波をパッシブアンテナ４で受信する。

まず始めに、特定基地局である未知電波対応型基地局２（ｋ₀）は、受信した電波に基づいて基準時刻における参照信号を生成する。その後、受信時刻算出基地局である未知電波対応型基地局２（１）、２（２）、２（Ｋ）は、未知電波対応型基地局２（ｋ₀）で生成された参照信号と、自らが生成したフィルタバンクによる処理後の受信信号の周波数成分との相関演算を行うことにより、直接波受信時刻の計測を行う。

未知電波対応型基地局２（ｋ₀）は、先の図２に示したような未知電波の推定を行うための構成を有している。電波源１から発信された電波は、マルチパスを発生する。そして、パッシブアンテナ４は、直接波と間接波とを含む電波を受信することとなる。その後、受信機５、Ａ／Ｄ変換器６、での処理は、図１６で説明したものと同一であり、Ａ／Ｄ変換器６からディジタル化された受信信号ｚ_1,k・・・ｚ_N,kが出力される。

ディジタル化された受信信号ｚ_1,k・・・ｚ_N,kは、フィルタバンク１４に伝達される。フィルタバンク１４は、Ｎ個の帯域通過フィルタ１４（１）〜１４（Ｎ）から構成されており、ｉ番目（ただし、ｉは１≦ｉ≦Ｎの整数）の帯域通過フィルタ１４（ｉ）は、次式（２２）により受信信号のｉ番目の周波数成分の信号ｙ_i,kを出力する。

その後、相関処理器１２およびＭＵＳＩＣ処理器１３は、実施の形態１と同様に動作して直接波受信時刻の推定値τ_1,k（チルダ）を計測する。

実施の形態２によれば、基地局のいずれか１つを特定して未知電波の推定を行うとともに、他の基地局において、この推定結果に基づいて直接波の受信時刻を計測することにより、各基地局で未知電波を受信した場合にも直接波受信時刻が計測でき、電波源の位置の時間差測位を行うことができる。さらに、受信時刻算出基地局である他の基地局においては、ＦＦＴ処理器の代わりに複数の帯域通過フィルタを有するフィルタバンクを用いることによっても、直接波受信時刻を計測することが可能となる。

実施の形態３．
実施の形態１および２では、１台の未知電波対応型基地局２（ｋ₀）を未知電波の推定を行う基地局として特定することにより直接波受信時刻の計測を行い、電波源１の位置を推定する場合について説明した。本実施の形態３では、直接波受信時刻の計測および電波源の位置の算出を、異なる基地局を未知電波の推定を行う基地局として特定した場合についてそれぞれ個別に行うことにより複数の推定結果を求め、それらの推定結果から最終的に電波源１の位置を特定する場合について説明する。実施の形態３における測位装置の構成は、図１と同様である。

本実施の形態３におけるそれぞれの未知電波対応型基地局２（１）、２（２）、２（ｋ₀）、２（Ｋ）は、それぞれが未知電波を推定して参照信号を生成し、その参照信号を他の未知電波対応型基地局に対して送信する機能を有している。さらに、それぞれの未知電波対応型基地局２（１）、２（２）、２（ｋ₀）、２（Ｋ）は、それぞれの参照信号に基づいて、直接波受信時刻を求め、時間差測位器３に送信する機能も兼ね備えている。

図７は、本発明の実施の形態３において時間差測位器３が電波源１の位置を特定するために行う処理手順を示したフローチャートである。ステップＳ７０１において、時間差測位器３は、未知電波を推定する特定基地局の番号をｋ₀＝１とする。ステップＳ７０２において、時間差測位器３は、特定基地局をｋ₀＝１としたときにそれぞれの未知電波対応型基地局で算出された直接波受信時刻に基づいて、実施の形態１と同様の手順で電波源１の位置ｐ１（チルダ）を算出する。

そして、ステップＳ７０３において、時間差測位器３は、特定基地局の番号ｋ₀が基地局の台数Ｋよりも小さいかを判断する。時間差測位器３は、ｋ₀＜Ｋであると判断した場合には、ステップＳ７０４において、特定基地局の番号ｋ₀に１を加算してｋ₀を更新する。そして、時間差測位器３は、ステップＳ７０２、Ｓ７０３の処理を繰り返すことにより、Ｋ台全ての基地局のそれぞれを特定基地局とした場合における電波源１の位置ｐｋ₀（ただし、ｋ₀は１≦ｋ₀≦Ｋの整数）を算出する。

ステップＳ７０３でｋ₀＜Ｋでないと判断した場合には、Ｋ台のそれぞれを特定基地局としたときの電波源１の位置ｐｋ₀がそろったこととなり、時間差測位器３は、ステップＳ７０５において、ｐｋ₀（チルダ）の平均値ｐ（チルダ）を算出する。このようにして、複数の基地局を未知電波推定のための特定基地局として設定し、電波源１の測位結果を複数求め、複数の測位結果を平均化することにより電波源１の位置を１つに特定することができる。

なお、図７の説明においては、Ｋ台全てを特定基地局とした場合のＫ通りの電波源１の位置を求めたが、これに限定されない。必ずしもＫ台全てではなく、２台以上の複数の基地局を特定基地局とする場合にも、複数の測定結果から電波源１の位置を１つに特定することができる。

また、単純に全ての測位結果を平均化せずに、例えば、最大値および最小値を削除して平均化することも可能であり、種々の統計的な手法により、複数の測位結果から１つの測位結果を求めることができる。

また、実施の形態３においては、複数の特定基地局が、参照信号を生成する本来の機能と、電波受信時刻を計測する機能（すなわち、本来は受信時刻算出基地局が有している機能）とを兼ね備えている場合について説明したが、これに限定されない。参照信号を生成する本来の機能のみを有する基地局を複数の特定基地局として用いることも可能である。

実施の形態３によれば、複数の基地局を未知電波推定のための特定基地局として設定し、複数の測位結果に基づいて電波源の位置を特定することが可能となる。これにより、１台の基地局のみを未知電波推定のための特定基地局として設定して求めた測位結果と比較して、測位誤差を抑えた測位結果を得ることができ、測定精度の向上を図ることができる。

実施の形態４．
本実施の形態４では、電波源としてＲＦＩＤタグを用い、基地局としてＲＦＩＤリーダ／ライタを用いる場合について説明する。基本的な構成及び動作は、実施の形態１〜３で記載したものと同一である。

ＲＦＩＤタグは、耐環境性に優れており、水、油、薬品等の汚れ、あるいは外乱光による影響を受けない。さらに、ＲＦＩＤリーダ／ライタ側からの非接触電力伝送を行うことにより、ＲＦＩＤタグは、電池を持たずに半永久的に利用可能とすることができる。

実施の形態１〜３で説明した電波源１は、自らが電池（電源）を有しており、自ら能動的に電波を発信するものであり、ＲＦＩＤタグも同じように使用可能である。さらに、電波源１としてＲＦＩＤタグを用いる場合には、上述の非接触電力伝送技術を活用することにより、ＲＦＩＤタグ自身は、電池を持っておらず、自ら能動的に電波を発信はしない代わりに、ＲＦＩＤリーダ／ライタから送信された電波（あるいは電磁波）の受信電力を検知することにより、ＲＦＩＤリーダ／ライタに対して電波（あるいは電磁波）を送信することができる。

従って、ＲＦＩＤタグは、任意のＲＦＩＤリーダ／ライタから電波を受信することにより、測位装置内の複数のＲＦＩＤリーダ／ライタに対して電波を送信することができる。一方、複数のＲＦＩＤリーダ／ライタの中の特定のＲＦＩＤリーダ／ライタは、実施の形態１〜３で説明した特定基地局として機能し、ＲＦＩＤタグからの電波を未知電波として受信し、この未知電波を推定して基準時刻における参照信号を生成することができる。

さらに、特定のＲＦＩＤリーダ／ライタ以外の複数のＲＦＩＤリーダ／ライタは、実施の形態１〜３で説明した受信時刻算出基地局として機能し、特定のＲＦＩＤリーダ／ライタから受信した基準信号と、ＲＦＩＤタグから受信した未知電波との相関から電波受信時刻を計測する。そして、時間差測位器３は、特定のＲＦＩＤリーダ／ライタ以外の複数のＲＦＩＤリーダ／ライタで計測されたそれぞれの電波受信時刻荷基づいて、ＲＦＩＤタグの位置を特定する。

実施の形態４によれば、電波源および基地局として、ＲＦＩＤタグおよびＲＦＩＤリーダ／ライタを用いた場合にも、ＲＦＩＤタグからの未知電波に基づいてＲＦＩＤタグの位置を特定できる測位装置を得ることができる。さらに、ＲＦＩＤリーダ／ライタ側からの非接触電力伝送技術を用いれば、ＲＦＩＤタグは、電池レス化され、耐環境性に優れるとともに、半永久的に利用可能となる。

実施の形態５．
本実施の形態５は、直接波と間接波の振幅比が複素数の場合にも、受信した未知電波に基づいて直接波の受信時刻をゼロとした基準時刻における参照信号を生成することができる特定基地局について説明する。

図８は、本発明の実施の形態５の特定基地局における未知電波の推定を行う構成を示した図であり、特定基地局である未知電波対応型基地局２（ｋ₀）の構成を示すものである。図８におけるパッシブアンテナ４、受信機５、Ａ／Ｄ変換器６、ＦＦＴ処理器７は、実施の形態１の構成を示した図２のものと同一である。

複素振幅比対応型参照信号生成器１８は、図２における参照信号生成器８に相当し、ＦＦＴ処理された未知電波の周波数成分の信号を推定し、直接波受信時刻を算出するための参照信号を生成する。ここで、複素振幅比対応型参照信号生成器１８は、直接波と間接波の振幅比が複素数のときにも参照信号が生成でき、この点が図２における参照信号生成器８と異なる。

図９は、本発明の実施の形態５における複素振幅比対応型参照信号生成器１８の構成図である。この複素振幅比対応型参照信号生成器１８は、間接波受信時刻推定器９、複素振幅比推定器１９、複素振幅比対応型初期位相推定器２０で構成される。

間接波受信時刻計測器９は、実施の形態１で説明したものと同一であり、ＦＦＴ処理された未知電波の周波数成分の信号から間接波の受信時刻を推定する。複素振幅比推定器１９は、時刻ゼロで受信される直接波の振幅と、推定された間接波の受信時刻における間接波の振幅との複素振幅比を推定する。さらに、複素振幅比対応型初期位相推定器２０は、間接波受信時刻計測器９の出力信号とＦＦＴ出力器７から出力される受信信号の周波数成分ｙ_i,k0から参照信号の初期位相を推定する。

次に、動作、および測位のための具体的な演算式について説明する。未知電波対応型基地局２（ｋ₀）は、先の図８に示したような未知電波の推定を行うための構成を有している。電波源１から発信された電波は、マルチパスを発生する。そして、パッシブアンテナ４は、直接波と間接波とを含む電波を受信することとなる。

その後、受信機５、Ａ／Ｄ変換器６、ＦＦＴ処理器７での処理は、図１６で説明したものと同一である。ＦＦＴ出力器７から出力される受信信号の周波数成分ｙ_i,k0は、電波源１から発信された直接波成分と間接波成分とを加えた周波数成分ｙ'_i,k0を含んでいる。この周波数成分ｙ'_i,k0は、次式（２３）で表される。

直接波の振幅値であるＣ’_1、k0と間接波の振幅値であるＣ’_2、k0との振幅比は、次式（２４）のように、Ｃ’_k0で表すことができる。ここで、Ｃ’_k0は、複素数となっている。

未知電波対応型基地局２（ｋ₀）による受信信号の周波数成分ｙ_i,k0は、図９に示した複素振幅比対応型参照信号生成器１８内の間接波受信時刻推定器９に伝達される。間接波受信時刻推定器９は、実施の形態１で説明したように、間接波の受信時刻の推定値τ_2,k0（チルダ）と、式（１３）で示した信号ベクトルＰｙ_i、k0を算出する。

間接波の受信時刻の推定値τ_2,k0（チルダ）と信号ベクトルＰｙ_i、k0は、複素振幅比推定器１９に入力される。複素振幅比推定器１９は、まず、次式（２５）〜（２７）により信号ベクトルＳ_i、k0（チルダ）を求める。

ここで、信号ベクトルＳ_i、k0（チルダ）は、次式（２８）に示すような成分で構成される。

次に、複素振幅比推定器１９は、次式（２９）により、複素振幅比Ｃ’_k0（ｉ）の偏角値の推定値Ａｒｇ［Ｃ’_k0（チルダ）（ｉ）］を求める。ここで、Ａｒｇ［α］は、複素数αの偏角を表している。

また、複素振幅比推定器１９は、次の二次方程式（３０）の解から、複素振幅比Ｃ’_k0（ｉ）の絶対値の推定値｜Ｃ’_k0（チルダ）（ｉ）｜を求める。

具体的には、この解は、次式（３１）により表される。

さらに、複素振幅比推定器１９は、式（２９）の偏角推定値Ａｒｇ［Ｃ’_k0（チルダ）（ｉ）］と、式（３１）の絶対値推定値｜Ｃ’_k0（チルダ）（ｉ）｜を用いて、次式（３２）により、複素振幅比推定値Ｃ’_k0（チルダ）を求める。

最終的に、複素振幅比推定器１９は、間接波の受信時刻の推定値τ_2、k0（チルダ）並びに複素振幅比推定値Ｃ’_k0（チルダ）を複素振幅比対応型初期位相推定器２０に伝達する。

これに対して、複素振幅比対応型初期位相推定器２０は、複素振幅比推定器１９から間接波の受信時刻の推定値τ_2、k0（チルダ）および複素振幅比推定値Ｃ’_k0（チルダ）を受信するとともに、ＦＦＴ処理器７から受信信号の周波数成分ｙ_i,k0を受信する。そして、複素振幅比対応型初期位相推定器２０は、まず始めに、次式（３３）で表される位相検出用信号Ｈ_iを算出する。

次に、複素振幅比対応型初期位相推定器２０は、上式（３３）で算出した位相検出用信号Ｈ_iの偏角から、参照信号周波数成分Γ_i、jの初期位相推定値φ_i（チルダ）を次式（３４）により算出する。

さらに、複素振幅比対応型初期位相推定器２０は、次式（３５）により、参照信号の周波数成分の推定値Γ_i、j（チルダ）を算出する。

このようにして、複素振幅比対応型参照信号生成器１８は、参照信号生成器８と同様に、受信信号の周波数成分ｙ_i,k0から参照信号の周波数成分の推定値Γ_i、j（チルダ）を算出することができる。

その後、各基地局は、実施の形態１と同様に動作し、複素振幅比対応型初期位相推定器２０から出力される参照信号の周波数成分の推定値Γ_i、i（チルダ）に基づいて、直接波受信時刻の推定値τ_1、k（チルダ）を計測する。さらに、時間差測位器３は、各基地局からの直接波受信時刻の推定値τ_1、k（チルダ）に基づいて電波源１の位置を求める。

実施の形態５によれば、特定基地局は、直接波と間接波の振幅比が複素数のときにも、受信した未知電波に基づいて直接波の受信時刻をゼロとした基準時刻における参照信号を生成することができる。

実施の形態６．
本実施の形態６では、本発明による測位装置および位置特定方法を車両位置の測位に適用する応用例について説明する。図１０は、本発明の実施の形態６において車両位置の測位に位置特定方法を適用した応用例の説明図である。フリーフローＥＴＣ（Electronic Toll Collection）システムを用いて、車両位置の測位を実現するものである。

図１０において、車両２１は、ＥＴＣからの情報を電波として送信するＥＴＣ情報送信アンテナ２２を備えている。一方、路側には、ＥＴＣ情報送信アンテナ２２からの送信電波を受信するＫ個の路側アンテナ２３（１）〜２３（Ｋ）が備えられている。

図１１は、本発明の実施の形態６の応用例における時間差測位の説明図である。時間差測位器３は、各路側アンテナ２３（１）〜２３（Ｋ）からの直接波受信時刻の推定値τ_1、k（チルダ）から、車両２１を時間差測位する。未知電波対応型基地局２と時間差測位器３の機能は、実施の形態１で説明したものと同じである。

図１２は、本発明の実施の形態６の応用例における車両２１に搭載されたＥＴＣ情報送信アンテナ２２を含むＥＴＣの構成図である。ＥＴＣ情報格納器２４は、ＥＴＣ処理を行うために必要な情報を格納している。送信機２５は、ＥＴＣ情報格納器２４に格納された情報を、ＥＴＣ情報送信アンテナ２２を介して電波として送信する。

次に、動作について説明する。ＥＴＣ情報送信アンテナ２２を介して送信された情報は、各路側アンテナ２３（１）〜２３（Ｋ）で受信される。路側アンテナ２３（１）〜２３（Ｋ）で受信された電波は、未知電波対応型基地局２（１）〜２（Ｋ）に伝達される。

実施の形態５と同様に、ｋ₀番目の未知電波対応型基地局２（ｋ₀）は、参照信号の周波数成分の推定値Γ_i、i（チルダ）を求めて、各未知電波対応型基地局２（１）〜２（Ｋ）に伝達する。各未知電波対応型基地局２（１）〜２（Ｋ）は、直接波受信時刻の推定値τ_1、k（チルダ）を計測し、時間差測位器３に伝達する。時間差測位器３は、直接波受信時刻の推定値τ_1、k（チルダ）を用いて車両２１の位置を時間差測位する。

ここで、ＥＴＣを搭載していない車両が路側アンテナ２３（１）〜２３（Ｋ）の下を通過すると、通行料金が徴収できない等の問題が発生する。そのため、ＥＴＣを搭載していない車両を判別することが必要となる。

本実施の形態６の応用例によれば、ＥＴＣを搭載している車両位置を特定することができる。したがって、ＥＴＣを搭載している車両とＥＴＣを搭載していない車両が同時に路側アンテナ２３（１）〜２３（Ｋ）の下を通過する際にも、特定される車両位置の情報を用いることにより、車両のＥＴＣ搭載の有無を判別することができる。

実施の形態６によれば、本発明の測位装置および位置特定方法をＥＴＣシステムに適用することにより、車両位置の測位を容易に実現することができる。

実施の形態７．
本実施の形態７では、本発明による測位装置および位置特定方法を車両位置の測位に適用する別の応用例について説明する。図１３は、本発明の実施の形態７において車両位置の測位に位置特定方法を適用した応用例の説明図である。スマートプレートを用いて、車両位置の測位を実現するものである。

ここで、スマートプレートとは、車両に関連する自動車登録情報等の種々の情報を記録したＩＣチップが、車両のナンバープレート上に搭載されており、種々の情報を送信する機能を有するもののことである。

図１３において、車両２１は、スマートプレートのＩＣチップからの情報を電波として送信するスマートプレート送信アンテナ２４を備えている。一方、路側には、スマートプレート送信アンテナ２４からの送信電波を受信するＫ個の路側アンテナ２３（１）〜２３（Ｋ）が備えられている。

図１４は、本発明の実施の形態７の応用例における車両２１に搭載されたスマートプレート送信アンテナ２４を含むスマートプレートの構成図である。車両情報格納器２６は、種々の情報を記録したＩＣチップに相当する。送信機２５は、車両情報格納器２６に格納された情報を、スマートプレート送信アンテナ２４を介して電波として送信する。

次に、動作について説明する。スマートプレート送信アンテナ２４を介して送信された情報は、各路側アンテナ２３（１）〜２３（Ｋ）で受信される。路側アンテナ２３（１）〜２３（Ｋ）で受信された電波は、未知電波対応型基地局２（１）〜２（Ｋ）に伝達される。

実施の形態６と同様に、ｋ₀番目の未知電波対応型基地局２（ｋ₀）は、参照信号の周波数成分の推定値Γ_i、i（チルダ）を求めて、各未知電波対応型基地局２（１）〜２（Ｋ）に伝達する。各未知電波対応型基地局２（１）〜２（Ｋ）は、直接波受信時刻の推定値τ_1、k（チルダ）を計測し、時間差測位器３に伝達する。時間差測位器３は、直接波受信時刻の推定値τ_1、k（チルダ）を用いて車両２１の位置を時間差測位する。

実施の形態７によれば、本発明の測位装置および位置特定方法をスマートプレートシステムに適用することにより、車両位置の測位を容易に実現することができる。スマートプレートを搭載した車両の位置が特定できることにより、立ち入り禁止区域に入った車両を発見できる等の利点がある。

Claims (16)

1. 電波源から発信された電波を受信し、参照信号と受信電波との相関から電波受信時刻を計測する複数の受信時刻算出基地局と、
   前記複数の受信時刻算出基地局からのそれぞれの前記電波受信時刻に基づいて前記電波源の位置を特定する時間差測位器と
   を備えた測位装置において、
   前記電波源から発信された未知電波を推定して基準時刻における参照信号を生成する特定基地局をさらに備え、
   前記複数の受信時刻算出基地局は、前記特定基地局からの前記参照信号と受信した未知電波との相関から電波受信時刻を計測する測位装置。
2. 請求項１に記載の測位装置において、
   前記電波源は、ＲＦＩＤタグであり、
   前記複数の受信時刻算出基地局は、複数のＲＦＩＤリーダ／ライタであり、
   前記特定基地局は、特定ＲＦＩＤリーダ／ライタである
   測位装置。
3. 請求項２に記載の測位装置において、
   前記ＲＦＩＤタグは、前記複数のＲＦＩＤリーダ／ライタあるいは前記特定ＲＦＩＤリーダ／ライタの何れかからの電波の受信電力を検知することにより電波を発信する測位装置。
4. 請求項２に記載の測位装置において、
   前記ＲＦＩＤタグは、未知電波として周波数成分の絶対値が周波数に依らずに一定となる電波を発信し、
   前記特定ＲＦＩＤリーダ／ライタは、前記未知電波の位相を求めることにより基準時刻における参照信号を生成する
   測位装置。
5. 請求項２に記載の測位装置において、
   前記特定ＲＦＩＤリーダ／ライタは、
   前記ＲＦＩＤタグから発信された電波を受信する受信機と、
   前記受信機からの受信電波を複数の周波数成分の周波数信号に変換する周波数成分算出手段と、
   前記周波数成分算出手段からの前記周波数信号に基づいてそれぞれの周波数成分に対応した参照信号を生成する参照信号生成手段と
   を備えた測位装置。
6. 請求項５に記載の測位装置において、
   前記参照信号生成手段は、
   前記周波数成分算出手段からの前記周波数信号に基づいて、直接波の受信時刻をゼロとしたときのマルチパスを経由して受信された間接波の受信時刻を推定する間接波受信時刻推定手段と、
   時刻ゼロにおける直接波の振幅と、推定された前記間接波の受信時刻における間接波の振幅との振幅比を推定する振幅比推定手段と、
   前記間接波の受信時刻および前記振幅比に基づいて前記ＲＦＩＤタグから発信された電波の周波数成分の初期位相を推定してそれぞれの周波数成分に対応した参照信号を生成する初期位相推定手段と
   を備えた測位装置。
7. 請求項５に記載の測位装置において、
   前記参照信号生成手段は、
   前記周波数成分算出手段からの前記周波数信号に基づいて、直接波の受信時刻をゼロとしたときのマルチパスを経由して受信された間接波の受信時刻を推定する間接波受信時刻推定手段と、
   時刻ゼロにおける直接波の振幅と、推定された前記間接波の受信時刻における間接波の振幅との複素振幅比を推定する複素振幅比推定手段と、
   前記間接波の受信時刻、前記複素振幅比、および前記周波数信号に基づいて前記ＲＦＩＤタグから発信された電波の周波数成分の初期位相を推定してそれぞれの周波数成分に対応した参照信号を生成する複素振幅比対応型初期位相推定手段と
   を備えた測位装置。
8. 請求項６または７に記載の測位装置において、
   前記間接波受信時刻推定手段は、ＭＵＳＩＣ処理を用いて前記間接波の受信時刻を推定する測位装置。
9. 請求項２に記載の測位装置において、
   前記ＲＦＩＤリーダ／ライタは、
   前記ＲＦＩＤタグから発信された電波を受信する受信機と、
   前記受信機からの受信電波を複数の周波数成分の周波数信号に変換する周波数成分算出手段と、
   前記特定ＲＦＩＤリーダ／ライタからの参照信号と前記周波数成分算出手段からの周波数信号とに基づいて周波数成分ごとの相関信号を生成する相関信号生成手段と、
   前記相関信号生成手段からの前記相関信号に基づいて未知電波の電波受信時刻を計測する電波受信時刻推定手段と
   を備えた測位装置。
10. 請求項９に記載の測位装置において、
    前記電波受信時刻推定手段は、ＭＵＳＩＣ処理を用いて前記電波受信時刻を推定する測位装置。
11. 請求項５ないし１０のいずれか１項に記載の測位装置において、
    前記周波数成分算出手段は、ＦＦＴ処理を行うことにより前記受信電波から前記複数の周波数成分の周波数信号を算出するＦＦＴ処理手段である測位装置。
12. 請求項５ないし１０のいずれか１項に記載の測位装置において、
    前記周波数成分算出手段は、帯域通過フィルタを用いることにより前記受信電波から前記複数の周波数成分の周波数信号を算出するフィルタバンクである測位装置。
13. 請求項１ないし１０のいずれか１項に記載の測位装置において、
    前記測位装置は、前記特定ＲＦＩＤリーダ／ライタとして複数の特定ＲＦＩＤリーダ／ライタを備え、
    前記複数のＲＦＩＤリーダ／ライタは、前記複数の特定ＲＦＩＤリーダ／ライタからのそれぞれの参照信号と受信した未知電波との相関から電波受信時刻をそれぞれの参照信号別に計測し、
    前記時間差測位器は、前記複数のＲＦＩＤリーダ／ライタからのそれぞれの前記電波受信時刻に基づいて参照信号別に前記ＲＦＩＤタグの位置を特定し、さらに参照信号別に特定された前記ＲＦＩＤタグの位置から前記ＲＦＩＤタグの位置を特定する
    測位装置。
14. 請求項１３に記載の測位装置において、
    前記複数の特定ＲＦＩＤリーダ／ライタのそれぞれは、他の特定ＲＦＩＤリーダ／ライタからのそれぞれの参照信号と受信した未知電波との相関から電波受信時刻をそれぞれの参照信号別に計測し、
    前記時間差測位器は、前記複数のＲＦＩＤリーダ／ライタおよび前記複数の特定ＲＦＩＤリーダ／ライタからのそれぞれの前記電波受信時刻に基づいて参照信号別に前記ＲＦＩＤタグの位置を特定し、さらに参照信号別に特定された前記ＲＦＩＤタグの位置から前記ＲＦＩＤタグの位置を特定する
    測位装置。
15. 特定基地局で、電波源から発信された未知電波を推定して基準時刻における参照信号を生成するステップと、
    複数の受信時刻算出基地局で、前記特定基地局からの前記参照信号と受信した未知電波との相関から電波受信時刻を計測するステップと、
    時間差測位器で、前記複数の受信時刻算出基地局からのそれぞれの前記電波受信時刻に基づいて前記電波源の位置を特定するステップと
    を備えた測位装置の位置特定方法。
16. 複数のＲＦＩＤリーダ／ライタあるいは特定ＲＦＩＤリーダ／ライタの何れかからＲＦＩＤタグに電波を送信するステップと、
    前記ＲＦＩＤタグで、前記電波の受信電力を検知することにより未知電波を発信するステップと、
    前記特定ＲＦＩＤリーダ／ライタで、前記ＲＦＩＤタグから発信された前記未知電波を推定して基準時刻における参照信号を生成するステップと、
    前記複数のＲＦＩＤリーダ／ライタで、前記特定ＲＦＩＤリーダ／ライタからの前記参照信号と受信した未知電波との相関から電波受信時刻を計測するステップと、
    時間差測位器で、前記複数の複数のＲＦＩＤリーダ／ライタからのそれぞれの前記電波受信時刻に基づいて前記ＲＦＩＤタグの位置を特定するステップと
    を備えた測位装置の位置特定方法。
    

Images