(57)【要約】【課題】抽象狭帯域測距を使用して第１の無線信号送受信機と第２の無線信号送受信機との間の距離を決定する方法を提供する。【解決手段】前記方法は、信号位相及び信号強度の測定（２）に基づいて、一方向周波数領域チャネル応答に比例する値の予備推定値を計算するステップ（４）と、隣接周波数対に対して、第１の無線信号送受信機によって使用されるクロックと第２の無線信号送受信機によって使用されるクロックとの間での時間同期オフセットに比例する値の推定値を計算するステップ（６）と、時間同期オフセットに比例する値の予備推定値及び隣接推定値に基づいて、一方向周波数領域チャネル応答に比例する値の最終推定値を決定するステップ（６）とを含み、最終的な距離決定（１０）には最終的な推定値が使用される。【選択図】図３
2019113529.tif 000002
【特許請求の範囲】
【請求項１】
第１の無線信号送受信機と第２の無線信号送受信機との間の距離を決定する方法であって、前記方法は、
第１の測定結果セットと第２の測定結果セットとを受信するステップ（２）を含み、
前記第１の測定結果セットは、前記第２の無線信号送受信機から送信された信号に基づいて、前記第１の無線信号送受信機によって取得され、
前記第２の測定結果セットは、前記第１の無線信号送受信機から送信された信号に基づいて、前記第２の無線信号送受信機によって取得され、
前記第１の測定結果セットは、複数の周波数の各々について、位相測定値と信号強度測定の測定対を含み、
前記第２の測定結果セットは位相測定値を含み、又はオプションで位相測定値及び信号強度測定値の測定対を含み、
前記方法は、
前記複数の周波数の各周波数に対して、前記第１の測定結果セットと前記位相測定値からの測定対に基づいて、又はオプションで前記第２の測定結果セットからの測定対に基づいて、一方向周波数領域チャネル応答に比例する値の予備推定値を計算するステップ（４）と、
隣接周波数対に対して、前記第１の無線信号送受信機によって使用されるクロックと第２の無線信号送受信機によって使用されるクロックとの間での、時間同期オフセットに比例する値の推定値を計算するステップ（６）とを含み、
前記時間同期オフセットに比例する値の推定値は、隣接周波数対における各周波数について、一方向周波数領域チャネル応答に比例する値の計算された予備推定値の位相及び位相測定値に基づいており、又はオプションで第１の測定結果セットからの測定対に基づいており、及びオプションで、隣接周波数対における周波数間の周波数差に基づいており、
前記方法は、
前記複数の周波数の各周波数に対して、前記時間同期オフセットに比例する値の予備推定値及び隣接推定値に基づいて、前記一方向周波数領域チャネル応答に比例する値の最終推定値を決定するステップ（８）と、
前記最終推定値に基づいて前記第１の無線信号送受信機と前記第２の無線信号送受信機との間の距離を決定するステップ（１０）とを含む方法。
【請求項２】
前記一方向周波数領域チャネル応答に比例する値の予備推定値の計算するステップ（４）は、
前記第１の測定結果セットと前記位相測定値からの測定対に基づいて、又はオプションで前記第２の測定結果セットからの測定対に基づいて、双方向周波数領域チャネル応答に比例する値の推定値を計算するステップと、
前記双方向周波数領域チャネル応答に比例する値の推定値に基づいて、前記一方向周波数領域チャネル応答に比例する前記予備推定値を計算するステップとを含む請求項１に記載の方法。
【請求項３】
前記測定対の各測定対は、複素数として表し、
前記複素数の絶対値は信号強度測定値に対応する振幅を表し、前記複素数の偏角は位相測定値を表し、
前記予備推定値及び前記最終推定値はそれぞれ複素数で表し、前記複素数の絶対値が振幅応答を表し、前記複素数の前記偏角が位相応答を表す請求項１又は２に記載の方法。
【請求項４】
前記双方向周波数領域チャネル応答に比例する値の推定値を計算するステップは、前記第１の測定結果セットからの前記測定対を表す複素数に、前記第２の測定結果セットからの前記測定対を表す複素数を乗算することを含み、又は含むものとして表す請求項３に記載の方法。
【請求項５】
前記双方向周波数領域チャネル応答に比例する値の推定値に基づいて、前記一方向周波数領域チャネル応答に比例する値の予備推定値を計算するステップ（４）は、前記双方向周波数領域チャネル応答に比例した推定値の複素数平方根をとるステップを含み、又は含むものとして表す請求項３又は４に記載の方法。
【請求項６】
前記平方根を取っているときに、−π／２とπ／２との間の位相を有する解が選択される請求項５に記載の方法。
【請求項７】
隣接周波数対に対する時間同期オフセットに比例する値の推定値を計算するステップ（６）は、
前記隣接周波数対における各周波数について、前記第１の測定結果セットからの前記測定対を表す複素数と、予備推定値を表す複素数の複素共役との積の偏角を取るステップと、
前記偏角の差を計算するステップと、
前記偏角の差を、前記隣接周波数対における周波数間の周波数差の２π倍で除算するステップとを含む請求項３〜６のうちのいずれか１つに記載の方法。
【請求項８】
前記一方向周波数領域チャネル応答に比例する値の最終推定値を決定するステップ（８）は、
前記時間同期オフセットに比例する値の各推定値について、前記時間同期オフセットに比例する値の隣接推定値に関して、所定の値によって決定される度合いまで、前記時間同期オフセットに比例する値の推定値が異なるかどうかを判定ステップと、
そのような差が判定されたときに、前記予備推定値の位相を反転させる前記最終推定値を形成するステップとを含む請求項１〜７のうちのいずれか１つに記載の方法。
【請求項９】
前記所定の値は、時間同期オフセットに関連して表されるときに、前記周波数差の逆数の半分である請求項８に記載の方法。
【請求項１０】
前記第１の無線信号送受信機と前記第２の無線信号送受信機との間の距離を決定するステップは、ＩＦＦＴ及び／又は超分解能アルゴリズムに基づいて決定される請求項１〜８のうちのいずれか１つに記載の方法。
【請求項１１】
処理ユニット上で実行されるとコンピュータプログラム製品は、処理ユニットに対して、請求項１〜１０のうちのいずれか１つに記載の方法を実行させるコンピュータ可読媒体。
【請求項１２】
第２の無線信号送受信機との距離を決定するように構成された無線信号送受信機であって、
前記無線信号送受信機は、第２の無線信号送受信機から送信された信号に基づいて、複数の周波数のそれぞれについて、位相測定値と信号強度測定値との測定対を含む第１の測定結果セットを取得する（２）ように構成された測定ユニットと、
第１の無線信号送受信機から送信された信号に基づいて、前記第２の無線信号送受信機によって取得された第２の測定結果セットを受信するように構成された受信機とを備え、
前記第２の測定結果セットは、複数の周波数のそれぞれについて、位相測定値を含み、又はオプションで位相測定値及び信号強度測定値の測定対を含み、
前記無線信号送受信機は、
処理ユニットを備え、
前記処理ユニットは、
前記複数の周波数の各周波数に対して、前記位相測定値に基づいて、又はオプションで第１の測定結果セットからの測定対と、第２の測定結果セットからの測定対とに基づいて、一方向周波数領域チャネル応答に比例する値の予備推定値を計算し（４）、
隣接周波数対に対して、第１の無線信号送受信機によって使用されるクロックと第２の無線信号送受信機によって使用されるクロックとの間での時間同期オフセットに比例する値の推定値を計算し（６）、
前記時間同期オフセットに比例する値の推定値は、隣接周波数対における各周波数について、一方向周波数領域チャネル応答に比例する値の計算された予備推定値の位相及び位相測定値に基づくものであり、又はオプションで第１の測定結果セットからの測定対に基づくものであり、及びオプションで隣接周波数対における周波数間の周波数差に基づくものであり、
前記処理ユニットは、
前記複数の周波数の各周波数に対して、前記時間同期オフセットに比例する値の予備推定値及び隣接推定値に基づいて、前記一方向周波数領域チャネル応答に比例する値の最終推定値を決定し（８）、
前記最終推定値に基づいて前記第１の無線信号送受信機と前記第２の無線信号送受信機との間の距離を決定する（１０）無線信号送受信機。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の概念は、第１の無線信号送受信機と第２の無線信号送受信機との間の距離を決定する方法に関する。本発明はまた、処理装置上で実行されるときコンピュータプログラム製品は処理装置に方法を実行させるコンピュータ可読命令を格納するコンピュータ可読媒体を備えるコンピュータプログラム製品に関し、さらに、第２の無線信号送受信機との距離を決定するように構成された無線信号送受信機に関する。
【背景技術】
【０００２】
２つの無線送受信機（トランシーバ）間の距離を決定するための狭帯域測距方法は、従来技術において知られている。多くの場合、狭帯域システムは受信信号強度を使用する。しかし、同様に、あるいは代替的に、いくつかの異なるキャリア（マルチキャリア位相差）の位相測定値を備える方法が例えば非特許文献１に記載されている。
【０００３】
最終的な距離（範囲）決定のために、このような測定の結果は、例えば、非特許文献２，３及び４に記載されているように、逆高速フーリエ変換（ＩＦＦＴ）アルゴリズム又はいわゆる「超分解能アルゴリズム」に供給されてもよい。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Kluge: Ranging with IEEE 802.15.4 Narrow-Band PHY, Response to Call for Final Proposals, IEEE P802.15, November 2009.
【非特許文献２】Schmidt: IEEE Transactions on Antennas and Propagation, Vol AP-34, No. 3, pp. 276-280, March 1986.
【非特許文献３】Sakar: IEEE Antennas and Propagation Magazine, Vol. 37, No. 1, pp. 48-55, February 1995.
【非特許文献４】Roy: IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. 37, No. 7, July 1989.
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明の概念の目的は、２つの無線送受信機間の距離を決定するための改良された狭帯域測距方法を提供することである。
【０００６】
本発明のこれらの目的及び他の目的は、独立請求項に定義される本発明によって少なくとも部分的に満たされる。好ましい実施形態は、従属請求項に記載されている。
【課題を解決するための手段】
【０００７】
本発明の概念によれば、第１の無線信号送受信機と第２の無線信号送受信機との間の距離を決定する方法であって、前記方法は、
第１の測定結果セット及び第２の測定結果セットを受信するステップを含み、
前記第１の測定結果セットは、前記第２の無線信号送受信機から送信された信号に基づいて、前記第１の無線信号送受信機によって取得され、
前記第２の測定結果セットは、前記第１の無線信号送受信機から送信された信号に基づいて、前記第２の無線信号送受信機によって取得され、
前記第１の測定結果セットは、複数の周波数の各々について、位相測定値と信号強度測定値の第１の測定対とを含み、
前記第２の測定値セットは位相測定値を含み、又はオプションとして、位相測定値及び信号強度測定値の第２測定対を含み、
前記方法は、
前記複数の周波数の各々に対して、第１の測定結果セットと位相測定値からの測定対に基づいて、又はオプションとして前記第２の測定結果セットから前記測定対に基づいて、一方向周波数領域チャネル応答に比例する値の予備推定値を計算するステップと、
隣接周波数対に対して、前記第１の無線信号送受信機によって使用されるクロックと前記第２の無線信号送受信機によって使用されるクロックとの間での、時間同期オフセットに比例する値の推定値を計算するステップとを含み、
前記時間同期オフセットに比例する値の推定値は、隣接周波数対における各周波数について、一方向周波数領域チャネル応答に比例する値の計算された予備推定値の位相と、前記位相測定値とに基づくものであり、又はオプションとして、第１の測定結果セットからの測定対に基づくものであり、及びオプションとして、隣接周波数対における周波数間の周波数差に基づくものであり、
前記方法は、
前記複数の周波数の各周波数に対して、前記時間同期オフセットに比例する値に対する、予備推定値及び隣接推定値に基づいて、一方向周波数領域チャネル応答に比例する値の最終推定値を決定するステップと、
前記最終推定値に基づいて前記第１無線信号送受信機と前記第２無線信号送受信機との間の距離を決定するステップとを含む。
【０００８】
信号強度測定という用語は、受信信号の電力又は振幅のいずれかに比例する任意の測定をカバーすると理解されるべきである。
【０００９】
測定結果の受信は、別の装置から送信された結果を受信する装置、又は装置にすでにローカルに保存されている測定結果を使用する装置のいずれかと理解する必要がある。
【００１０】
この方法は、距離決定が一方向周波数領域チャネル応答に基づくことを可能にする。そうでなければ、一方向周波数領域チャネル応答は、２つの送受信機でのクロック間のタイミング同期によってもたらされる位相の曖昧さのために、明確に再構成することが困難である。この曖昧さは、上記のように時間同期オフセットに比例する値の推定値を計算することによって解決される。
【００１１】
２つのクロック間のタイミング同期オフセットは、測定間でゆっくりと変化すると予想される。従って、タイミング同期オフセットに隣接推定値を使用すると、上記のように一方向周波数領域チャネル応答に比例した値の予備推定値を補正するために使用できる余分な情報が得られる。時間同期オフセットに比例した値の推定値の計算すること、とりわけ予備推定値に基づいて計算することで、予備推定値の位相曖昧性に関する誤った選択に敏感になる時間同期オフセットに比例する値の推定値を生成し、最終的な推定のために予備的推定の位相を反転させるべきであることを通知することができる。
【００１２】
一方向周波数領域チャネル応答に基づいて距離判定を実行することにより、双方向周波数領域チャネル応答は一方向伝播モードの多数の交差積を含むので、マルチパス伝搬からの干渉に対して敏感でなくなる結果を生成し、伝播モードの数の２乗で成長する。従って、感度の低い距離測定方法が提供され、すなわち、マルチパス効果に対してよりロバストである。
【００１３】
第１の測定結果セットと第２の測定結果セットの両方を使用して距離を計算するために、他の送受信機からの測定結果を受信する方法は、送受信機のうちの１つで実行されてもよい。しかし、この方法は、任意の装置で実行されてもよく、例えば、場合によっては第１及び第２の送受信機よりも処理能力が高い外部装置で実行されてもよい。次いで、外部装置は、それぞれの送受信機から第１及び第２の測定結果セットを受信し、それらの間の距離を決定することができる。外部装置はさらに、送受信機がそれらの間の距離を知ることができるように、決定された距離を送受信機に通信で通知することができる。
【００１４】
一実施形態によれば、前記一方向周波数領域チャネル応答に比例する値の予備推定値の前記計算は、
第１の測定結果セットからの測定対と第２の測定結果セットからの測定対とに基づいて、双方向周波数領域チャネル応答に比例する値の推定値を計算するステップと、
前記双方向周波数領域チャネル応答に比例する値の推定値に基づいて、前記一方向周波数領域チャネル応答に比例する予備推定値を計算するステップとを含む。
【００１５】
前記測定対に基づいて双方向周波数領域チャネル応答に比例する値をまず計算することは、ある２つの送受信機間の時間同期オフセットが自然に相殺されるという利点を有する。このことは、所定周波数において、一方の方向におけるこのオフセットに起因する測定された位相シフトが、他方の方向と反対の符号で作用するからである。
【００１６】
一実施形態によれば、前記測定対の各測定対は複素数として表すことができ、前記複素数の絶対値は信号強度測定値に対応する振幅を表し、前記複素数の偏角は位相測定値を表し、そして前記予備推定値及び前記最終推定値はそれぞれ複素数で表すことができ、前記複素数の絶対値が振幅応答を表し、前記複素数の偏角が位相応答を表す。
【００１７】
当該分野での慣習として確立されているように、複素数表現は、周期的に変化する信号に関連する便利な表記法及び便利な計算を提供し、ここでもまた、確立された慣例に従って、実際の物理的実数値信号は、対応する複素数の実数部によって表される。しかしながら、当業者によって容易に理解されるように、実際の計算を実行するときは任意の他の適切な表現が使用されてもよい。特に、計算が１つ以上の複素数を含む演算として「表現されてもよい」と述べられている場合、使用された実際の表現に関係なく、数学的に等価な計算をカバーすることが理解されるであろう。
【００１８】
一実施形態によれば、前記双方向周波数領域チャネル応答に比例する値の推定値の前記計算は、前記第１の測定結果セットからの前記測定対を表す複素数に、前記第２の測定結果セットからの前記測定対を表す複素数を乗算することを含み、又は含むものとして表してもよい。
【００１９】
一実施形態によれば、双方向周波数領域チャネル応答に比例する値の推定値に基づいており、一方向周波数領域チャネル応答に比例する値の予備推定値の前記計算は、双方向周波数領域チャネル応答に比例する推定値の複素数平方根を取ることを含み、又は含むものとして表してもよい。複素数平方根をとることは、固有のπ（１８０度）の位相曖昧性を持ち、２つの可能な解のうちの１つを選択する必要がある。例えば、一実施形態によれば、前記平方根をとるとき、−π／２とπ／２との間の位相を有する解が選択されてもよく、すなわち、実数部が正の解が選択されてもよい。
【００２０】
一実施形態によれば、隣接周波数対の時間同期オフセットに比例する値の前記推定値の前記計算は、前記隣接周波数対における各周波数について、前記第１の測定結果セットからの前記測定対を表す複素数と前記予備推定を表す前記複素数の前記複素共役との積の偏角を取るステップと、前記偏角の差を計算するステップと、前記偏角の前記差を、前記隣接周波数対における前記周波数間の前記周波数差の２π倍で割るステップとを含む。
【００２１】
一実施形態によれば、前記一方向周波数領域チャネル応答に比例する値の最終推定値の決定は、前記時間同期オフセットに比例する値の各推定値について、
前記時間同期オフセットに比例する値の前記推定値が、前記時間同期オフセットに比例する値の隣接推定値に対して、所定の値によって決定される程度に異なるかどうかを判定するステップと、
そのような差が決定されたときに、前記予備的推定値の位相を反転させるステップと、
及びオプションとして、前記最終推定値を形成するステップとを含む。
【００２２】
時間同期オフセットに比例する値の隣接推定値に関して、所定の値によって決定される程度まで、時間同期オフセットに比例する値の推定値が異なるかどうかを判定することは、厳密には値の減算を含むものとして理解されるべきではなく、減算、除算、又は比較を含む任意の適切な方法を含むことが理解されるべきである。次いで、所定の値が、使用されるスキームと互換性があるように選択され、全てが本発明の範囲内である。
【００２３】
一実施形態によれば、前記所定の値は、時間同期オフセットに関連して表される場合、前記周波数差の逆数の半分である。これは、仮推定のために間違った位相解が選択された場合、推定された時間同期オフセットが予想されるジャンプである。逆に、そのようなジャンプが観察されない場合、予備的な推定の修正は必要ない。従って、タイミングオフセットの隣接推定値の間の差に関連して表され、ゼロに近い差、又は周波数差の逆数の半分のオーダーの絶対値との差が予想され、測定の不確実性を考慮する。例えば、代わりに、タイミングオフセットの隣接推定値を除算することに関して表現される場合、１に近い値又は１／２又は２に近い値のいずれかが予想される。
【００２４】
一実施形態によれば、第１の無線信号送受信機と第２の無線信号送受信機との間の距離の決定は、ＩＦＦＴ及び／又は超分解能アルゴリズムに基づくアルゴリズムを使用する。
【００２５】
本発明の第２の態様によれば、コンピュータ可読命令を格納するコンピュータ可読媒体を含むコンピュータプログラム製品であって、処理ユニット上で実行されると、コンピュータプログラム製品は、処理ユニットに上記の方法を実行させるようにする。
【００２６】
この第２の態様の効果及び特徴は、第１の態様に関連して上述したものと概ね類似している。第１の態様に関して述べた実施形態は、第２の態様と大きく互換性がある。
【００２７】
本発明の第３の態様によれば、第２の無線信号送受信機との距離を決定するように構成された無線信号送受信機であって、
前記無線信号送受信機は、
第２の無線信号送受信機から送信された信号に基づいて第１の測定結果セットを取得するように構成された測定ユニットを備え、
第１の測定結果セットは、複数の周波数の各々について、位相測定値と信号強度測定値との測定対とを含み、
前記無線信号送受信機は、
前記第１の無線信号送受信機から送信された信号に基づいて、第２の無線信号送受信機によって取得された第２の測定結果セットを受信するように構成された受信機を備え、
前記第２の測定結果セットは、複数の周波数の各々について、位相測定値と信号強度測定値との測定対とを含み、
前記無線信号送受信機は処理ユニットを備え、
前記処理ユニットは、
前記複数の周波数の各周波数に対して、第１の測定結果セットからの測定対と第２の測定結果セットからの測定対とに基づいて、一方向周波数領域チャネル応答に比例する値の予備推定値を計算し、
隣接周波数対に対して、時間同期オフセットに比例する値の推定値を計算する第１の無線信号送受信機によって使用されるクロックと、第２の無線信号送受信機によって使用されるクロックとの間での、時間同期オフセットの推定値を計算し、
前記時間同期オフセットの推定値は、隣接周波数対における各周波数について、一方向性周波数領域のチャネル応答に比例する値の計算された予備推定値と、第１の測定結果セットからの測定対に基づいており、
及びオプションとして、隣接周波数対における周波数間の周波数差に基づき、
前記処理ユニットは、
前記複数の周波数の各周波数に対して、前記時間同期オフセットに比例する値の予備推定値及び隣接推定値に基づいて、一方向周波数領域チャネル応答に比例した値の最終推定値を決定し、
前記最終推定値に基づいて前記第１の無線信号送受信機と前記第２の無線信号送受信機との間の距離を決定する。
【００２８】
この第３の態様の効果及び特徴は、第１及び第２の態様に関連して上述したものと概ね類似している。第１及び第２の態様に関連して述べた実施形態は、第３の態様と大きく互換性がある。
【００２９】
従って、無線信号送受信機は、ロバストな方法で第２の無線信号送受信機との距離を決定することができる。測定を実行する一対の無線信号送受信機において、無線信号送受信機の１つは、送受信機間の距離を決定するために計算を実行することができる。次いで、この送受信機は、両方の送受信機がそれらの間の距離を知るように、決定された距離を他の送受信機に通信することができる。
【００３０】
本発明の概念の上記の目的、追加の目的、特徴及び利点は、添付の図面を参照して、以下の例示的及び非限定的な詳細な説明によってよりよく理解されるであろう。図面において、特に断らない限り、類似の要素には同様の参照番号が使用される。
【図面の簡単な説明】
【００３１】
【図１】局部発振器（ＬＯ）周波数及び送受信機Ａ及びＢの状態を時間で示した図である。
【図２】適用されるメソッドの例である。
【図３】この方法のステップを要約したフローチャートである。
【発明を実施するための形態】
【００３２】
この方法は、第１の無線信号送受信機、装置Ａと第２の無線信号送受信機、装置Ｂとの間の距離の決定のための、マルチキャリア位相差（ＭＣＰＤ）測距の原理を使用して、第１の測定結果セットと第２の測定結果セットを入力として使用する。第１の測定結果セットは、装置Ｂのような第２の無線信号送受信装置から送信された信号に基づいて、装置Ａのような第１の無線信号送受信機によって取得され、第２の測定結果セットは、第１の無線信号送受信機から送信された信号に基づいて、装置Ｂのような第２の無線信号送受信機によって取得される。各測定結果セットは、複数の周波数のそれぞれについて、位相測定値と信号強度測定値との測定対を生成する。
【００３３】
測定結果の取得は、２つの装置Ａ及びＢが測距パラメータにおいて一致し、それらの周波数を整列調整させることから開始し（例えば、キャリア周波数オフセット（ＣＦＯ）推定及び較正を使用して）、粗い時間同期を実現し、すなわち、ＡとＢの両方がクロック（クロック）例えば、開始フレームデリミタ（Ｓｔａｒｔ  Ｆｒａｍｅ  Ｄｅｌｉｍｉｔｅｒ；ＳＦＤ）の送受信時に、両方の装置Ａ及びＢはローカルステートマシンを制御するために使用される。ステートマシンは、どの送受信機が何をしているかを制御する。
【００３４】
図１に示すように、測定は、以下のステップで実行されてもよい。
０：装置Ａ及び装置Ｂは、局部発振器ＬＯを所定の周波数に設定し、ループカウンタｋ＝０を設定する。
１：装置ＡはＬＯを送信し、装置Ｂは位相測定（φＢ［ｋ］）を実行する。さらに、装置Ｂは、受信信号強度指示測定ＲＳＳＩＢ［ｋ］（図示せず）を実行する。
２：装置Ａ及びＢは送信方向を変更し、ＬＯを安定させるためのガード時間を可能にする。
３：装置ＢはＬＯを送信し、装置Ａは位相測定を実行する（φＡ［ｋ］）。さらに、装置Ａは、受信信号強度指示測定ＲＳＳＩＡ［ｋ］（図示せず）を実行する。
４：装置Ａ及び装置ＢはＬＯ周波数を所定の間隔Δｆだけ増加させ、ステップ１に戻る。このループは所定回数（Ｋｆ）繰り返される。
【００３５】
装置Ａと装置Ｂの両方は、それぞれの局部発振器（ＬＯ）信号を生成するために連続的に動作する位相ロックループ（ＰＬＬ）を備える。送信から受信又はその逆に切り替えるとき、ＰＬＬはオンのままであり、連続した位相信号を可能にする。
【００３６】
図３は、この方法のステップを要約している。図３のブロック２において、一旦測定が実行されると、距離を決定する方法はあまり時間的に重要ではない。従って、より多くの処理能力を有する第３の装置／エンティティ上で計算することができ、これは、例えば、クラウドの中で、エンティティが両方の送受信機からの測定データにアクセスできると仮定する。従って、この方法は、装置Ａ及び／又はＢ上で実行されてもよいが、第３の装置Ｃ上で収集されてもよく、ＡとＢの間の距離を計算することができ、装置Ｃがクラウド内に存在する可能性がある。装置が方法を実行しない場合、方法は、その方法を実行する装置に測定結果を送信するか、送信させる可能性がある。従って、例えば、装置Ｃが装置上で実行される場合、装置Ｂは、装置Ｃ（φＢ［０：Ｋｆ−１］）にすべての位相測定値を有するフレームを送信し、装置Ａは、すべての位相測定値を装置Ｃに送信する（φＡ［０：Ｋｆ−１］）。同様に、ＲＳＳＩ測定値ＲＳＳＩＡ［０：Ｋｆ−１］及び、オプションでＲＳＳＩＢ［０：Ｋｆ−１］は、この方法を実行する装置、例えば装置Ｃに送信される。従って、本方法を実施する装置は、φＡ［０：Ｋｆ−１］及びＲＳＳＩＡ［０：Ｋｆ−１］を含む第１の測定結果セットと、φＢ［０：Ｋｆ−１］及びオプションでＲＳＳＩＢ［０：Ｋｆ−１］を含む第２の測定結果セットとを含む。
【００３７】
とって代わって、装置Ａは方法を実行し、上記のように、装置Ｂのような第２の無線信号送受信装置から送信された信号に基づいて、第１の測定結果セットを取得するように構成された測定ユニットを備える。装置Ａは、装置Ａのような第１無線信号送受信機から送信された信号に基づいて、装置Ｂなどの第２の無線信号送受信機によって取得された第２の測定結果セットを受信するように構成された受信機をさらに備える。さらに、装置Ａは、以下に説明するように、本方法のステップを実行するための処理ユニットを備える。
【００３８】
各周波数及び測定のセットについて、一方向周波数領域応答に比例した複素数が形成され、ここで、その絶対値は信号強度測定値に対応する振幅を表し、前記複素数の偏角は位相測定値を表す。
【００３９】
ＨＡ［ｋ］＝ＡＡ［ｋ］ｅｘｐ（ｊφＡ［ｋ］）
【００４０】
ＨＢ［ｋ］＝ＡＢ［ｋ］ｅｘｐ（ｊφＢ［ｋ］）
【００４１】
ここで、ＡＡ［ｋ］及びＡＢ［ｋ］は、信号振幅に比例する値であり、例えば、対応するＲＳＳＩ値の平方根をとることによって計算される。
【００４２】
熱雑音又は位相雑音がない場合、ｋ番目の周波数におけるこれらの測定された大きさ及び位相は、以下のように実際のチャネル応答Ｈ［ｋ］に関連する。
【００４３】
ＨＡ［ｋ］∝Ｈ［ｋ］ｅｘｐ（ｊ２πｆ［ｋ］Δｔ［ｋ］）
【００４４】
ここで、Δｔ［ｋ］は、ｋ番目の周波数の測定中にＡとＢとの間の時間オフセットを示す記号∝は比例関係を示し、すなわち、ａ［ｋ］∝ｂ［ｋ］は、ｋのすべての値に対してａ［ｋ］＝ｃ＊ｂ［ｋ］であり、ここで、ｃは未知の複素値であり、すべてのｋに対して同じである。タイミングオフセットΔｔ［ｋ］は、クロックドリフト及びジッタのために測定ごとに変化することがあり、オフセットは測定から測定までゆっくりとしか変化しない。同様に、Ｂにおいて、我々は以下を測定する。
【００４５】
ＨＢ［ｋ］∝Ｈ［ｋ］ｅｘｐ（−ｊ２πｆ［ｋ］Δｔ［ｋ］）
【００４６】
図３のブロック４において、双方向周波数領域のチャネル応答に比例する値の推定値Ｘ［ｋ］を形成することができる２つの値を乗算することによって、タイミングオフセットに関連する要因を相殺する。このことは、ＡとＢにおける測定値が互いに直後に取られるので、デルタｔ［ｋ］は変化しないと仮定することができるからである。
【００４７】
Ｘ［ｋ］＝ＨＡ［ｋ］ＨＢ［ｋ］∝（Ｈ［ｋ］）２
【００４８】
従って、双方向周波数領域チャネル応答に比例する値の推定値の計算は、第１の測定結果セットからの測定対と第２の測定結果セットからの測定対とに基づいて、さらに、前記計算は、前記第１の測定結果セットからの前記測定対を表す複素数に、前記第２の測定結果セットからの前記測定対を表す複素数を乗算することを含み、又は含むものとして表してもよい。
【００４９】
とって代わって、Ｘ［ｋ］は振幅に関して、次式のように、Ａのみの測定値に基づいて計算することができる。
【００５０】
【数１】
2019113529.tif 000003
【００５１】
ここで、
【数２】
2019113529.tif 000004
は絶対二乗演算子を表す。
【数３】
2019113529.tif 000005
はＲＳＳＩＡ［ｋ］と等しいことに注意してください。
【００５２】
従って、ここでは、双方向周波数領域チャネル応答に比例する値の推定値を計算することは、第１の測定結果セットからの測定対と、第２の測定結果セットからの位相測定値とに基づいている。
【００５３】
以下では、一方向周波数領域チャネル応答Ｈ［ｋ］は、捕捉された位相測定値φＡ［ｋ］及びφＢ［ｋ］からａ）Ｘ［ｋ］及びｂ）サイド情報を使用して再構成される。
【００５４】
一方向周波数領域のチャネル応答Ｈ［ｋ］の予備推定値Ｈｓｑｒｔ［ｋ］は、双方向周波数領域チャネル応答Ｘ［ｋ］に比例した推定値の平方根をとることによって、次式のように計算される。
【００５５】
【数４】
2019113529.tif 000006
【００５６】
前記式は、上記の比例関係に基づく真の一方向周波数領域チャネル応答に関連し、ここで、ｃ［ｋ］は＋１又は１のいずれかであり、複雑な平方根をとる固有の位相の曖昧さに起因する。ｃ［ｋ］の値を推定するために、時間オフセットが急激に変化することはないという事実を利用する。予備的な推定のために、−π／２とπ／２との間の位相を有する解を選択することができ、すなわち、正の実数部を選択することができる。曖昧さを解決するために、サイド情報を使用することが提案される。すなわち、ＡとＢとの間の時間オフセットは、以下に説明するように、手順中にゆっくりと変化するだけである。
【００５７】
従って、各周波数について、第１の測定結果セットと位相測定値とからの測定対に基づいて、又はオプションで第２の測定結果セットからの測定対に基づいて、一方向周波数領域チャネル応答に比例する値の予備推定値が計算される。
【００５８】
図３のブロック６において、Δｆ−１に等しい曖昧性を有する時間測定オフセットΔｔの値を推定することができる。指数ｋを持つ周波数に対するΔｔの推定値は、次のように計算される。
【００５９】
【数５】
2019113529.tif 000007
【００６０】
ここで、∠は複素数の偏角を取ることを示す。推定値
【数６】
2019113529.tif 000008
は、すべてのｋについて近似的に同じでなければならない。
【００６１】
ＨＡの偏角／位相のみが使用されるので、この値の大きさは、ＨＡの偏角を使用して計算から除外することができる。例えば、周波数間隔が一定である場合には、除算を省略することもできる。従って、隣接周波数対に対して、第１の無線信号送受信機によって使用されるクロックと第２の無線信号送受信機によって使用されるクロックとの間での、時間同期オフセットに比例する値の推定値は、隣接周波数対における各周波数について、一方向周波数領域チャネル応答に比例する値の計算された予備推定値の位相及び位相測定値に基づいて、又はオプションでは、第１の測定結果セットからの測定対に基づいて、及びオプションで隣接周波数対における周波数間の周波数差に基づいている。さらに、前記推定値は、第１の測定結果セットから前記測定対を表す複素数と、前記予備推定値を表す複素数の複素共役との積の偏角を取り、前記偏角の差を計算することで計算され、及びオプションで、前記偏角の差を、隣接周波数対における周波数間の周波数差を２π倍した値で除算することで計算される。
【００６２】
図３のブロック８において、推定値
【数７】
2019113529.tif 000009
は、隣接周波数間で１／２Δｆ−１のオーダーでジャンプするときに、このことは、予備推定値Ｈｓｑｒｔ［ｋ］における不正確な符号変化を示す。従って、ｋ≧ｍであるすべての値Ｈｓｑｒｔ［ｋ］は−１を乗算しなければならない。２つのジャンプがｍとｎでｍ＜ｎで観察されるときに、ｍ≦ｋ＜ｎの間の値には−１を乗算し、剰余には（−１）２＝１を乗算する。３つ以上のジャンプの場合、手順は単純に拡張される。乗算の回数を制限するために、マスクを作成して、どの値を符号反転させ、どの値を反転させないかを追跡することができる。マスクは
【数８】
2019113529.tif 000010
で示される。一方向周波数領域チャネル応答の最終推定値は次式になる。
【００６３】
【数９】
2019113529.tif 000011
【００６４】
従って、各周波数について、一方向周波数領域チャネル応答に比例する値の最終推定値は、時間同期オフセットに比例する値の予備推定値及び隣接推定値に基づいて決定される。さらに、時間同期オフセットに比例する値の各推定値について、それが異なるかどうかが判定され、時間同期オフセットに比例する値の隣接推定値に関して、所定の値によって決定される程度までそのような差異が判定された場合には、予備推定値と、オプションで後続の予備推定値とは反転されて、一方向周波数領域チャネル応答の最終推定値を形成する。
【００６５】
最後に、図３のブロック１０において、第１の無線信号送受信機と第２の無線信号送受信機との間の距離は、最終推定値に基づいて決定されてもよい。例えば、それ自体が知られている逆高速フーリエ変換（ＩＦＦＴ）が使用されることができるだけでなく、より進んだ信号処理技術も、それ自体は周知のように、一般に超分解能アルゴリズムと呼ばれる。再構成された一方向周波数領域チャネル応答Ｈ［ｋ］は、ほとんどの測距アルゴリズムがマルチパスからより多くの干渉を軽減することを可能にし、問題の順序／コンポーネントの数が減少するため、問題が発生する。マルチパスが存在する場合、見通し内（ＬＯＳ）成分の遅延の推定に干渉する成分の数が減少し、測距及び位置特定がより正確になる。
【００６６】
コンピュータ可読媒体を含むコンピュータプログラム製品は、コンピュータ可読命令処理ユニット上で実行されると、コンピュータプログラム製品は、処理ユニットに上記の方法を実行させるようにする。
【００６７】
この方法は、処理ユニット内で実行されてもよく、上記のように装置Ａ、Ｂ又はＣに配置することができる。
【００６８】
処理ユニットは、ハードウェアで、又はソフトウェアとハードウェアの任意の組み合わせとして実装することができる。処理ユニットの機能の少なくとも一部は、例えば、汎用コンピュータ上で実行されるソフトウェアとして実現されてもよい。従って、システムは、１つ又は複数の処理ユニットを備えることができ、中央演算処理装置（ＣＰＵ）所望の機能を実現するために１つ又は複数のコンピュータプログラムの命令を実行することができる。
【００６９】
あるいは、処理ユニットは、例えば組み込みシステムにおいて、又は特別に設計された処理ユニットとして、特定用途向け集積回路（ＡＳＩＣ）又はフィールドプログラマブルゲートアレイ（ＦＰＧＡ）のような、特定用途向け集積回路を含む。
【００７０】
以下の議論は、双方向周波数領域チャネル応答（ＴＷＦＤＣＲ）Ｘ［ｋ］が、一方向周波数領域チャネル応答であるマルチパス効果に対してより敏感である理由を示すのに役立つ。Ｘ［ｋ］の問題は、それが無線コンポーネントだけではなく、多くのクロス積（外積）もある。この議論では、Ｘ［ｋ］は以下のように定義される。
【００７１】
Ｘ［ｋ］＝（Ｈ［ｋ］）２
【００７２】
ここで、Ｈ［ｋ］は次式で表される。
【００７３】
【数１０】
2019113529.tif 000012
【００７４】
すなわち、ＴＷＦＤＣＲ  Ｘ［ｋ］は一方向周波数チャネル応答Ｈ［ｋ］の２乗である。
【００７５】
言い換えると、無線チャネルは、Ｌ成分（ＬＯＳ及び反射／マルチパス）の合計であり、各成分は、ある複素振幅ａ１及び遅延τｌを有する。一般性を失うことなく、ｌ＝０はＬＯＳ成分を表す。最終的には、ＬＯＳ構成要素の遅延を推定することが目的であり、他の構成要素は推定に干渉する。これをＸ［ｋ］の方程式に代入すると、次式を得る。
【００７６】
【数１１】
2019113529.tif 000013
【００７７】
言い換えれば、Ｘ［ｋ］はＫ＜Ｌ２での複素指数の和である。ある程度単純化することはできるが、組み合わせの数は依然として約１／２Ｌ２である。換言すれば、推定に干渉する成分の数は、パス／成分の数と共に２の累乗で増加する。
【００７８】
方法などのアルゴリズムの動作を説明するために、それは測定されたチャネルに適用された。チャネルは、周波数指数ｋに対応する１ＭＨｚの周波数ステップサイズを有する８０個の異なる周波数で測定された。結果を図２に示す。真の一方向周波数領域チャネル応答（１ＷＦＤＣＲ）Ｈ［ｋ］は実線で示され、予備推定値Ｈｓｑｒｔ［ｋ］が示されている。周波数インデックス（５−１９，３７−４９，５３，６１−６２，８０）に対して、Ｈ［ｋ］とＨｓｑｒｔ［ｋ］は異なる位相（右上の分図）を持つことに注意してください。左下の分図には、推定された時間オフセットΔｔ，ｅｓｔが示されている。推定された時間オフセットは、位相が正しくない区間が開始又は終了する場合を除き、一定であることに注意する。また、時間オフセットジャンプは、予想どおり０．５マイクロ秒に等しい大きさを有することにも留意されたい。
【００７９】
上記の本発明の概念は、限られた数の例を参照して主に説明されている。しかし、当業者には容易に理解されるように、上述したもの以外の他の例も、添付の請求項によって定義される本発明の概念の範囲内で等しく可能である。
【図１】
2019113529.tif 000014
【図２】
2019113529.tif 000015
【図３】
2019113529.tif 000016
【外国語明細書】
2019113529.tif 000017
2019113529.tif 000018
2019113529.tif 000019
2019113529.tif 000020
2019113529.tif 000021
2019113529.tif 000022
2019113529.tif 000023
2019113529.tif 000024
2019113529.tif 000025
2019113529.tif 000026
2019113529.tif 000027
2019113529.tif 000028
2019113529.tif 000029
2019113529.tif 000030
2019113529.tif 000031
2019113529.tif 000032
2019113529.tif 000033
2019113529.tif 000034
2019113529.tif 000035
2019113529.tif 000036
2019113529.tif 000037
2019113529.tif 000038
2019113529.tif 000039
2019113529.tif 000040