JP2006148457A

JP2006148457A - 測距システム，送信端末，受信端末，測距方法，およびコンピュータプログラム

Info

Publication number: JP2006148457A
Application number: JP2004334632A
Authority: JP
Inventors: Hiroaki Takano; 裕昭高野
Original assignee: Sony Corp
Current assignee: Sony Corp
Priority date: 2004-11-18
Filing date: 2004-11-18
Publication date: 2006-06-08
Anticipated expiration: 2024-11-18
Also published as: US20060104198A1; US7860014B2; CN1777158A; JP4561329B2

Abstract

【課題】マルチキャリア伝送方式において測距の誤差を軽減することが可能な，測距システム送信端末，受信端末，測距方法，およびコンピュータプログラムを提供する。
【解決手段】上記測距システム送信端末，受信端末，測距方法，およびコンピュータプログラムのうち，マルチキャリア伝送方式における端末間の距離を測定する測距システムは，パケットの送信時及び／又は受信時のデータが高速離散フーリエ変換され，さらに所定処理後の周波数領域のデータに対して，スペクトル推定アルゴリズムを適用することで，インパルスレスポンスを取得する周波数領域データ処理部（２１２）とを備えている。
【選択図】図２

Description

本発明はマルチキャリア伝送方式における測距システム，送信端末，受信端末，測距方法，およびコンピュータプログラムに関する。

従来から，複数のコンピュータを情報通信手段で接続し，ファイルやデータ等の情報を無線通信で送受信したり，共有したりすることが可能な無線ＬＡＮ（無線ローカルエリアネットワーク）を構成することが一般的に知られている。

上記無線ＬＡＮには，例えば，ＩＥＥＥ８０２．１１ａ規格等が存在している。上記ＩＥＥＥ８０２．１１ａの規格では，最大で５４Ｍｂｐｓの通信速度を達成する５ＧｈｚのＯＦＤＭ変調方式をサポートし，５２本のサブキャリアを利用して無線通信が行われる。

上記ＩＥＥＥ８０２．１１ａ等の無線ＬＡＮを利用して装置間の距離を測定する測距システムが開発されている。この無線ＬＡＮで送信器と受信器の装置間の距離を測る場合，まず相互でパケットの受信位置を特定する必要がある。

上記ＩＥＥＥ８０２．１１ａ等の無線ＬＡＮでパケットの受信位置を特定するには，同期位置でパケットの受信位置を特定する方法が考えられる。例えば，受信器に備わるＡＤ変換機による動作後，８０ＭＨｚのクロックで受信器に備わるディジタル回路が動作しているとする。

上記パケットの先頭には，既知の信号系列が埋め込まれており，受信機のＡＤ変換機による処理後，この既知の信号を用いて，受信信号との相関が計算される。パケットの既知信号の位置では，この相関値が大きくなるため，相関値のピークを探索（サーチ）することにより，どの位置でパケットを受信できたかを判別することができ，装置間の距離を測ることができる。

また一方で，無線ＬＡＮにおける測距とは全く異なる分野であるが，例えば，いわゆるアダプティブアレーで伝播解析する際に利用される技術としてスペクトル推定アルゴリズムのうちのＭＵＳＩＣアルゴリズムが存在する（非特許文献１，参照）。

菊間信良著，「アレーアンテナによる適応信号処理」，科学技術出版，１９９８年９月２０日，ｐ．１９１−２０３

しかしながら，例えば，マルチキャリア伝送方式である５ＧＨｚ系のＯＦＤＭ変調方式では，システムクロックが８０ＭＨｚ程度であって，そのシステムクロック単位でしか，送信器と受信器は同期の位置を判別することができないため，パケットの受信時刻に基づいて距離を測定する場合には，誤差が大きかった。

例えば，８０ＭＨｚのクロックでは１クロックあたり，１／８０ＭＨｚ＝１２ｎｓ（ナノ秒）となる。つまり，８０ＭＨｚのクロックにおける１クロックの誤差は距離にすると，光の速度は３．０×１０^８ｍ／ｓであることから，３．０×１０^８ｍ／ｓ×１２ｎｓ＝３．６ｍとなり，メートルオーダーの精度の測距しか行えず，誤差が大きいことが分かる。

したがって，家庭内や小規模なオフィス内等で上記ＩＥＥＥ８０２．１１ａなどの無線ＬＡＮを利用し，装置間の測距を実行したとしても，誤差が大きすぎるため有効な活用が望めなかった。

ところで，近年においては，ＭＵＳＩＣアルゴリズムが無線ＬＡＮに活用することが期待されている。従来におけるＭＵＳＩＣアルゴリズムは，アダプティブアレーの伝播解析に主として利用されるものであったが，かかるＭＵＳＩＣアルゴリズムを，各種無線ＬＡＮ分野において活用されることが期待されている。

本発明は，上記問題点に鑑みてなされたものであり，本発明の目的は，マルチキャリア伝送方式において測距の誤差を軽減することが可能な，新規かつ改良された測距システム，送信端末，受信端末，測距方法，およびコンピュータプログラムを提供することである。

上記課題を解決するため，本発明の第１の観点によれば，マルチキャリア伝送方式における端末間の距離を測定する測距システムが提供される。上記測距システムは，パケットの送信時及び／又は受信時のデータが高速離散フーリエ変換され，さらに所定処理後の周波数領域のデータに対して，スペクトル推定アルゴリズムを適用することで，インパルスレスポンスを取得する周波数領域データ処理部とを備えることを特徴としている。

本発明によれば，端末間の距離を測定する測距システムでは，端末間の距離を測定するために，パケットの送信時及び／又は受信時のデータを高速離散フーリエ変換し，さらに所定処理を施した後の周波数領域のデータに対して，スペクトル推定アルゴリズムを適用し，インパルスレスポンスを求めている。

上記インパルスレスポンスは，パケットの受信時刻及び／又は送信時刻を高精度に示す時刻情報であるように構成してもよい。かかる構成により，正確なパケットの受信時刻及び／又は送信時刻を取得することが可能となり，端末間の距離をより高精度に測距することが可能となる。

上記測距システムには，周波数領域のデータに存在するキャリアホールを補間する補間処理部と；上記補間後の周波数領域のデータに対して周波数平均処理する周波数平均処理部とを備え，上記所定処理には，補間処理部による補間処理と，周波数平均処理部による周波数平均処理とが含まれるように構成してもよい。かかる構成により，パケット送受信時のデータを高速離散フーリエ変換した周波数領域のデータにキャリアホールが存在したとしてもスペクトル推定アルゴリズムを適用し，適正なインパルスレスポンスを求めることが可能である。

上記補間処理部は，キャリアホールに相当する部分のデータを隣のキャリアからのデータで補間するように構成してもよい。

上記課題を解決するために，本発明の別の観点によれば，マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末が提供される。上記送信端末は，第１のパケットを通信相手の受信端末に送信するパケット送信部と；第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成部とを備え，さらに，上記第１のパケット送信時刻情報と，受信端末により第１のパケットの受信とともに送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，受信端末と送信端末との間のパケット往復時間を概算する往復時間概算部と；受信端末からの第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理部と；第１の周波数領域のデータに所定処理を施した後，該所定処理後の第１の周波数領域のデータに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求める周波数領域データ処理部と；第３のパケットを受信端末から所定のタイミングで受信し，パケット往復時間から第３のパケットに含む内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに第１のインパルスレスポンスと第３のパケットに含む第２のインパルスレスポンスとに基づいて，電波の伝搬往復時間を補正後伝搬往復時間に補正する伝搬時間補正部と；補正後伝搬往復時間と電波の速度とから送信端末と受信端末との間の距離を求める距離算出部とを備えることを特徴としている。

本発明によれば，送信端末は，第１のパケット送信時刻と第２のパケット受信時刻とからパケットが受信端末に行って戻ってくるターンアラウンドタイム等に相当するパケット往復時間を概算し，さらに内部処理遅延時間と，第１のインパルスレスポンスと，第２のインパルスレスポンスから，電波が受信端末に行って送信端末に再び戻ってくる高精度な伝搬往復時間を求めている。かかる構成によれば，高精度な伝搬往復時間を求めることが可能であるので，たとえマルチキャリア伝送方式であっても測距による誤差を極小化し，端末間の距離を正確に測定することが可能となる。

上記第２のインパルスレスポンスは，受信端末によって第１のパケット受信時のデータから第２の周波数領域のデータに高速離散フーリエ変換され，所定処理が施された後，その第２の周波数領域のデータにスペクトル推定アルゴリズムが適用されることで求められ，内部処理遅延時間は，該受信端末の内部処理の遅延時間であるように構成してもよい。

上記第１のインパルスレスポンスは，上記第２のパケットの受信時刻を高精度に示す時刻情報であって，上記第２のインパルスレスポンスは，第１のパケットの受信時刻を高精度に示す時刻情報であるように構成してもよい。

上記送信端末は，周波数領域データに存在するキャリアホールを補間する補間処理部と；上記補間後の周波数領域データに対して周波数平均処理する周波数平均処理部とを備え，上記所定処理には，補間処理部による補間処理と，周波数平均処理部による周波数平均処理とが含まれるように構成してもよい。

上記課題を解決するために，本発明の別の観点によれば，マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末の測距方法が提供される。上記測距方法は，第１のパケットを通信相手の受信端末に送信するパケット送信ステップと；第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成ステップとを含み，さらに，第１のパケット送信時刻情報と，受信端末により第１のパケットの受信とともに送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，受信端末と送信端末との間のパケット往復時間を概算する往復時間概算ステップと；受信端末からの第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換ステップと；第１の周波数領域のデータに所定処理を施した後，該所定処理後の第１の周波数領域のデータに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求めるステップと；第３のパケットを受信端末から所定のタイミングで受信し，パケット往復時間から第３のパケットに含む内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに第１のインパルスレスポンスと第３のパケットに含む第２のインパルスレスポンスとに基づいて，電波の伝搬往復時間を補正後伝搬往復時間に補正する伝搬時間補正ステップと；補正後伝搬往復時間と電波の速度とから送信端末と受信端末との間の距離を求める距離算出ステップとを含むことを特徴としている。

上記課題を解決するために，本発明の別の観点によれば，コンピュータをして，マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末として機能させるコンピュータプログラムが提供される。上記コンピュータプログラムは，第１のパケットを通信相手の受信端末に送信するパケット送信手段と；第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成手段とを含み，さらに，第１のパケット送信時刻情報と，受信端末により第１のパケットの受信とともに送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，受信端末と送信端末との間のパケット往復時間を概算する往復時間概算手段と；受信端末からの第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理手段と；第１の周波数領域のデータに所定処理を施した後，該所定処理後の第１の周波数領域のデータに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求める周波数領域データ処理手段と；第３のパケットを受信端末から所定のタイミングで受信し，パケット往復時間から第３のパケットに含む内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに第１のインパルスレスポンスと第３のパケットに含む第２のインパルスレスポンスとに基づいて，電波の伝搬往復時間を補正後伝搬往復時間に補正する伝搬時間補正手段と；補正後伝搬往復時間と電波の速度とから送信端末と受信端末との間の距離を求める距離算出手段とを含むことを特徴としている。

上記課題を解決するために，本発明の別の観点によれば，マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる受信端末が提供される。上記受信端末は，通信相手の送信端末からの第１のパケットを受信するとともに，第２のパケットを生成し，上記送信端末に第２のパケットを送信する第１のパケット送信部と；送信端末からの第１のパケット受信時のデータを周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理部と；周波数領域のデータに所定処理を施した後，上記所定処理後の周波数領域のデータに対してスペクトル推定アルゴリズムを適用し，インパルスレスポンスを求める周波数領域データ処理部と；受信端末自体の内部処理による処理遅延時間として内部処理遅延時間を求める処理遅延算出部と；上記インパルスレスポンスと内部処理遅延時間とが少なくとも含まれた第３のパケットを生成し，所定のタイミングで，送信端末に送信する第２のパケット送信部とを備えることを特徴としている。なお，本発明に係る所定のタイミングは，上記第１のパケットを受信した後であれば，いかなるタイミングでもよい。また，本発明にかかる第１のパケットを受信するとともに，第２のパケットを生成するタイミングは，例えば，第１のパケットを受信するのと略同時または第１のパケットの受信直後等を例示できる。

インパルスレスポンスは，前記第１のパケットの受信時刻を高精度に示す時刻情報であるように構成してもよい。

上記受信端末は，周波数領域データに存在するキャリアホールを補間する補間処理部と；上記補間後の周波数領域データに対して周波数平均処理する周波数平均処理部とを備え，上記所定処理には，補間処理部による補間処理と，前記周波数平均処理部による周波数平均処理とが含まれるように構成してもよい。

補間処理部は，キャリアホールに相当する部分のデータを隣のキャリアからのデータで補間するように構成してもよい。

上記課題を解決するために，本発明の別の観点によれば，マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末が提供される。上記送信端末は，第１のパケットを通信相手の受信端末に送信するパケット送信部と；第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成部とを備え，さらに，第１のパケット送信時刻情報と，受信端末により第１のパケットの受信とともに送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，受信端末と該送信端末との間のパケット往復時間を概算する往復時間概算部と；受信端末からの第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理部と；第３のパケットを受信端末から所定のタイミングで受信し，該第３のパケットに含まれる第２の周波数領域データに所定処理を施した後，該所定処理後の第２の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求めるとともに，第１の周波数領域データに所定処理を施した後，該所定処理後の第１の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第２のインパルスレスポンスを求める周波数領域データ処理部と；パケット往復時間から第３のパケットに含む内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに第１のインパルスレスポンスと第２のインパルスレスポンスとに基づいて，電波の伝搬往復時間を補正後伝搬往復時間に補正する伝搬時間補正部と；補正後伝搬往復時間と電波の速度とから送信端末と受信端末との間の距離を求める距離算出部とを備えることを特徴としている。

本発明によれば，第１のパケット送信時刻と第２のパケット受信時刻とからパケットが受信端末に行って戻ってくるターンアラウンドタイム等に相当するパケット往復時間を概算し，さらに内部処理遅延時間と，第１のインパルスレスポンスと，第２のインパルスレスポンスから，電波が受信端末に行って送信端末に再び戻ってくる高精度な伝搬往復時間を求めている。かかる構成によれば，高精度な伝搬往復時間を求めることが可能であるので，たとえマルチキャリア伝送方式であっても測距による誤差を極小化し，端末間の距離を正確に測定することが可能となる。また，送信端末が上記第１のインパルスレスポンスと，第２のインパルスレスポンスとを一元的に求めているため測距処理を集約化することができる。さらに受信端末の内部構成を送信端末よりも簡略化できる。

上記第２の周波数領域のデータは，受信端末によって第１のパケット受信時のデータから第２の周波数領域のデータに高速離散フーリエ変換され，内部処理遅延時間は，該受信端末の内部処理の遅延時間であるように構成してもよい。

上記第１のインパルスレスポンスは，第１のパケットの受信時刻を高精度に示す時刻情報であって，上記第２のインパルスレスポンスは，第２のパケットの受信時刻を高精度に示す時刻情報であるように構成してもよい。

上記課題を解決するために，本発明の別の観点によれば，マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末の測距方法が提供される。上記測距方法は，第１のパケットを通信相手の受信端末に送信するパケット送信ステップと；第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成ステップとを含み，さらに，上記第１のパケット送信時刻情報と，受信端末により第１のパケットの受信とともに送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，受信端末と該送信端末との間のパケット往復時間を概算する往復時間概算ステップと；受信端末からの第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理ステップと；第３のパケットを受信端末から所定のタイミングで受信し，該第３のパケットに含まれる第２の周波数領域データに所定処理を施した後，該所定処理後の第２の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求めるとともに，第１の周波数領域データに所定処理を施した後，該所定処理後の第１の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第２のインパルスレスポンスを求めるステップと；上記パケット往復時間から第３のパケットに含む内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに第１のインパルスレスポンスと第２のインパルスレスポンスとに基づいて，電波の伝搬往復時間を補正後伝搬往復時間に補正する伝搬時間補正ステップと；補正後伝搬往復時間と電波の速度とから送信端末と受信端末との間の距離を求める距離算出ステップとを含むことを特徴としている。

また，上記課題を解決するために，本発明の別の観点によれば，コンピュータをして，マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末として機能させるコンピュータプログラムが提供される。第１のパケットを通信相手の受信端末に送信するパケット送信手段と；上記第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成手段とを含み，さらに，第１のパケット送信時刻情報と，受信端末により第１のパケットの受信とともに送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，受信端末と該送信端末との間のパケット往復時間を概算する往復時間概算手段と；受信端末からの第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理手段と；第３のパケットを受信端末から所定のタイミングで受信し，該第３のパケットに含まれる第２の周波数領域データに所定処理を施した後，該所定処理後の第２の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求めるとともに，第１の周波数領域データに所定処理を施した後，該所定処理後の第１の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第２のインパルスレスポンスを求める周波数領域データ処理手段と；パケット往復時間から第３のパケットに含む内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに第１のインパルスレスポンスと第２のインパルスレスポンスとに基づいて，電波の伝搬往復時間を補正後伝搬往復時間に補正する伝搬時間補正手段と；補正後伝搬往復時間と電波の速度とから送信端末と受信端末との間の距離を求める距離算出手段とを含むことを特徴としている。

本発明によれば，コンピュータプログラムは，送信端末に対して，第１のパケット送信時刻と第２のパケット受信時刻とからパケットが受信端末に行って戻ってくるターンアラウンドタイム等に相当するパケット往復時間を概算させ，さらに内部処理遅延時間と，第１のインパルスレスポンスと，第２のインパルスレスポンスから，電波が受信端末に行って送信端末に再び戻ってくる高精度な伝搬往復時間を求めさせている。かかる構成によれば，高精度な伝搬往復時間を算出させることが可能であるので，送信端末にマルチキャリア伝送方式であっても測距による誤差を極小化させて，端末間の距離を正確に測定させることが可能となる。また，例えば，インパルスレスポンスを求める処理など，送信端末に測距に関する処理を一元的に実行させることが可能であるため，受信端末に実行させる処理を簡略化し，処理負荷を軽減させることができる。

上記課題を解決するために，本発明の別の観点によれば，マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる受信端末が提供される。上記受信端末は，通信相手の送信端末からの第１のパケットを受信するとともに，第２のパケットを生成し，送信端末に該第２のパケットを送信する第１のパケット送信部と；送信端末からの第１のパケット受信時のデータを周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理部と；受信端末自体の内部処理による処理遅延時間として内部処理遅延時間を求める処理遅延算出部と；周波数領域のデータと内部処理遅延時間とが少なくとも含まれた第３のパケットを生成し，所定のタイミングで，送信端末に送信する第２のパケット送信部とを備えることを特徴としている。

上記課題を解決するために，本発明の別の観点によれば，マルチキャリア伝送方式における端末間の距離を測定する測距システムが提供される。上記測距システムにおける送信端末は，第１のパケットを受信端末に送信する第１のパケット送信部と；第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成部とを備えており，また受信端末は，第１のパケットを受信するとともに，第２のパケットを生成し，送信端末に該第２のパケットを送信する第２のパケット送信部と；送信端末からの第１のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する第１の高速離散フーリエ変換処理部と；第１の周波数領域データに所定処理を施した後，該所定処理後の第１の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求める第１の周波数領域データ処理部と；受信端末自体の内部処理による処理遅延時間として内部処理遅延時間を求める処理遅延算出部と；第１のインパルスレスポンスと内部処理遅延時間とが少なくとも含まれた第３のパケットを生成し，所定のタイミングで，送信端末に送信する第３のパケット送信部とを備えており，さらに送信端末は，第１のパケット送信時刻情報と，時刻情報生成部により生成される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報とに基づき，受信端末と該送信端末間のパケット往復時間を概算する往復時間概算部と；受信端末からの第２のパケット受信時のデータを第２の周波数領域のデータに高速離散フーリエ変換する第２の高速離散フーリエ変換処理部と；第２の周波数領域データに所定処理を施した後，該所定処理後の第２の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第２のインパルスレスポンスを求める第２の周波数領域データ処理部と；パケット往復時間から内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに第１のインパルスレスポンスと第２のインパルスレスポンスとに基づいて電波の伝搬往復時間を補正することで，補正後伝搬往復時間を求める伝搬時間補正部と，補正後伝搬往復時間と電波の速度とから送信端末と受信端末との間の距離を求める距離算出部とを備えることを特徴としている。なお，上記測距システムでは，送信端末および受信端末がインパルスレスポンスを各々求めている。

上記課題を解決するために，本発明の別の観点によれば，マルチキャリア伝送方式における端末間の距離を測定する測距システムが提供される。上記測距システムにおける送信端末は，第１のパケットを受信端末に送信する第１のパケット送信部と；第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成部とを備えており，受信端末は，第１のパケットを受信するとともに，第２のパケットを生成し，送信端末に該第２のパケットを送信する第２のパケット送信部と；送信端末からの第１のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する第１の高速離散フーリエ変換処理部と；受信端末自体の内部処理による処理遅延時間として内部処理遅延時間を求める処理遅延算出部と；第１の周波数領域のデータと内部処理遅延時間とが少なくとも含まれた第３のパケットを生成し，所定のタイミングで，送信端末に送信する第３のパケット送信部とを備えており，さらに送信端末は，第１のパケット送信時刻情報と，時刻情報生成部により生成される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報とに基づき，受信端末と該送信端末間のパケット往復時間を概算する往復時間概算部と；受信端末からの第２のパケット受信時のデータを第２の周波数領域のデータに高速離散フーリエ変換する第２の高速離散フーリエ変換処理部と；第３のパケットに含まれる第１の周波数領域データに所定処理を施した後，該所定処理後の第１の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求めるとともに，第２の周波数領域データに所定処理を施した後，該所定処理後の第２の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第２のインパルスレスポンスを求める第１の周波数領域データ処理部と；パケット往復時間から内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに第１のインパルスレスポンスと第２のインパルスレスポンスとに基づいて電波の伝搬往復時間を補正することで，補正後伝搬往復時間を求める伝搬時間補正部と，補正後伝搬往復時間と電波の速度とから送信端末と受信端末との間の距離を求める距離算出部とを備えることを特徴としている。なお，上記測距システムでは，送信端末のみがインパルスレスポンスを全て一元的に求めている。

上記課題を解決するために，本発明の別の観点によれば，マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末が提供される。上記送信端末は，第１のパケットを受信端末に送信するパケット送信部と；第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成部とを備え，さらに送信端末は，第１のパケット送信時刻情報と，受信端末により第１のパケットの受信とともに該送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，受信端末と該送信端末間のパケット往復時間を概算する往復時間概算部と；受信端末からの第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理部と；第２の周波数領域データに所定処理を施した後，該所定処理後の第２の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求める周波数領域データ処理部と；第１パケット受信時のデータから第２の周波数領域のデータに高速離散フーリエ変換され，所定処理が施された後，その第２の周波数領域のデータにスペクトル推定アルゴリズムが適用されることで求められる第２のインパルスレスポンスと，該受信端末の内部処理の遅延時間である内部処理遅延時間とが含まれる第３のパケットを受信端末から所定のタイミングで受信し，パケット往復時間から内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに第１のインパルスレスポンスと第２のインパルスレスポンスとに基づいて電波の伝搬往復時間を補正することで，補正後伝搬往復時間を求める伝搬時間補正部と，補正後伝搬往復時間と電波の速度とから送信端末と受信端末との間の距離を求める距離算出部とを備えることを特徴としている。

以上説明したように，本発明によれば，マルチキャリア伝送方式における端末間の測距の精度を飛躍的に向上させ，端末間の測距による距離の誤差を軽減させることが可能となる。

以下，本発明の好適な実施の形態について，添付図面を参照しながら詳細に説明する。なお，以下の説明及び添付図面において，略同一の機能及び構成を有する構成要素については，同一符号を付することにより，重複説明を省略する。

（第１の実施の形態について）
まず，図１を参照しながら，第１の実施の形態にかかる測距システム１００について説明する。図１は，第１の実施の形態にかかる測距システムの概略的な構成の一例を示す説明図である。

図１に示すように，測距システム１００には，距離測定をする端末１０１ａ（以下，端末Ａと記載する場合もある）と，距離測定をされる端末１０１ｂ（以下，端末Ｂと記載する場合もある）とが少なくとも存在する。なお，端末１０１ａが距離測定をされてもよいし，端末１０１ｂが距離測定をしてもよく，さらに測距システム１００に３台以上の端末１０１が存在する場合でもよい。

また，第１の実施の形態にかかる端末１０１は，例えば，パーソナルコンピュータ（ＰＣ）等を例示することができるが，無線ＬＡＮでデータを送受信可能な装置であれば，かかる例に限定されず，如何なる装置でも実施可能である。

図１に示す測距システム１００における測距とは，端末Ａから端末Ｂまでの距離を測定することである。これは，端末Ａから送信したパケットが端末Ｂを経由して再びパケットが戻ってくるまでのターンアランドタイムの測定値（パケット往復時間）から端末１０１特有の遅延時間等を引き算することにより，ある空間での電波の遅延時間がまず求められる。その遅延時間から電波の速さ等を考慮して距離を測定するというものである。

次に，図２を参照しながら，第１の実施の形態にかかる端末１０１ａについて説明する。図２は，第１の実施の形態にかかる距離測定をする側の端末の概略的な構成の一例を示すブロック図である。

図２に示すように，第１の実施の形態にかかる端末１０１ａは，距離を測定する機能を有する距離測定ブロック２０１ａと，無線ＬＡＮでデータを送受信する機能を有する無線通信ブロック２０３とを備える。

上記距離測定ブロック２０１ａには，さらに，キャリアホールの補間処理を実行する補間処理部２１０と，移動平均（周波数平均）処理を実行する移動平均処理部（周波数平均処理部）２１１と，ＭＵＳＩＣアルゴリズムを適用してパケットの高精度な受信／送信位置を特定するＭＵＳＩＣアルゴリズム処理部（周波数領域データ処理部）２１２と，例えばパケット１〜パケット３等の種々のデータを記憶可能なメモリ２１３（２１３ａ，２１３ｂ，…２１３ｅ）と，端末１００間の距離を算出する距離計算処理部（往復時間概算部，伝搬時間補正部，距離算出部）２１４と，時刻設定部（時刻情報生成部）２１８とが備わっている。

補間処理部２１０は，ＦＦＴ（高速離散フーリエ変換処理部）２２４によるＦＦＴ処理後の周波数領域のデータには存在するキャリアホールが存在するが，そのキャリアホールに相当する部分のデータを例えば，該当部分の隣のキャリアからのデータで補間する。

上記ＭＵＳＩＣアルゴリズム処理部２１２は，パケットの受信／送信位置を正確に特定するために，ＭＵＳＩＣアルゴリズムを適用し，インパルスレスポンスを求める。なお，本実施の形態にかかるＭＵＳＩＣアルゴリズムの詳細については，一般的なＭＵＳＩＣアルゴリズムとほぼ同様であり，詳細には，例えば，菊間信良著の「アレーアンテナによる適応信号処理」，科学技術出版，１９９８年９月２０日，ｐ．１９１−２０３等に記載の通りである。

なお，第１の実施の形態にかかるメモリ２１３は，ＲＡＭ等の記憶手段に限定されず，例えば，フラッシュメモリまたはＨＤＤ（ハードディスクドライブ）等の記憶手段でも良い。

また，上記無線通信ブロック２０３は，ＲＦ処理部２２０と，信号をＡＤ変換するＡＤ変換処理部２２１と，パケット同期処理部２２２と，ＧＴ除去処理部２２３と，ＦＦＴ（ＦａｓｔＦｏｕｒｉｅｒＴｒａｎｓｆｏｒｍｓ：高速離散フーリエ変換）処理を実行するＦＦＴ（高速離散フーリエ変換処理部）２２４と，デコーダー２２５と，エンコーダー２２６と，ＩＦＦＴ（逆高速離散フーリエ変換）処理を実行するＩＦＦＴ２２７と，ＧＴ付与処理部２２８と，ＤＡ変換処理部２２９とを備えている。

上記ＲＦ処理部２２０は，少なくとも５ＧＨｚのＲＦ（Ｒａｄｉｏ−Ｆｒｅｑｕｅｎｃｙ）を，アンテナを介して送出／受信する。

パケット同期処理部２２２は，例えばパケットの先頭に存在する既知信号を用いて相関値のピークサーチをすることによって受信位置を特定するための同期位置を検出する。上記同期位置をもとに時刻設定部２１８によって，パケットの受信時刻／送信時刻が求められる。

メモリ２１３ａは，距離測定用のパケット１（第１のパケット）を送信した時刻を示す送信時刻情報（パケット１送信時刻情報（第１のパケット送信時刻情報））を記憶し，メモリ２１３ｂは，距離測定用のパケット２を受信した時刻を示す受信時刻情報（パケット２受信時刻情報）を記憶し，メモリ２１３ｃは，距離測定用のパケット２（第２のパケット）の受信時刻よりも一層正確性が高い補正後の時刻を示す補正値（パケット２受信時刻補正値（第２のインパルスレスポンス））を記憶し，メモリ２１３ｄは，通信相手の端末１０１ｂから受信したパケット３に含む端末１０１ｂの内部処理遅延値（内部処理で生じる遅延時間）を記憶し，メモリ２１３ｅは，通信相手の端末１０１ｂからのパケット３（第３のパケット）の受信時刻よりも高精度な時刻を示す補正値（パケット３受信時刻補正値（第３のインパルスレスポンス））と，通信相手の端末１０１ｂがパケット１を受信した受信時刻よりも正確性が高い補正後の時刻を示す補正値（パケット１受信時刻補正値（第１のインパルスレスポンス））を記憶する。なお，上記パケット１〜パケット３については後述する
なお，距離測定用のパケット１の受信時刻，パケット２の受信時刻，またはパケット３の受信時刻は，端末１０１によって，既知であるパケットの先頭ブロック信号とシステムクロックの同期位置に基づいて求められるため，大まかな受信時刻となる。より受信時刻の正確性を高めるため端末１０１は，例えば，上記ＭＵＳＩＣアルゴリズムを適用することで，上記パケット１またはパケット２の受信時刻よりも高精度な受信時刻（パケット１受信時刻補正値／パケット２受信時刻補正値）を取得することが可能となる。

上記時刻設定部２１８は，タイマー機能を有し，時々刻々と計時しながら時刻を示す時刻情報を生成している。さらに，パケット１〜パケット３のうちいずれかのパケットを送信又は受信したタイミングの時刻を示す時刻情報をパケット受信時刻情報又はパケット送信時刻情報としてメモリ２１３（２１３ａ，２１３ｂ）各々に格納する。

次に，図３を参照しながら，第１の実施の形態にかかる端末１０１ｂについて説明する。図３は，第１の実施の形態にかかる距離測定される側の端末の概略的な構成の一例を示すブロック図である。

図３に示すように，第１の実施の形態にかかる端末１０１ｂは，距離を測定するために所定処理を実行する機能を有する距離測定ブロック２０１ｂと，無線ＬＡＮでデータを送受信する機能を有する無線通信ブロック２０３とを備える。

上記距離測定ブロック２０１ｂには，補間処理部２１０と，移動平均処理部２１１と，ＭＵＳＩＣアルゴリズム処理部２１２と，メモリ２１３（２１３ｆ）と，パケットを受信するパケット受信部２１５と，パケットを生成又は送信するパケット生成／送信部２１６（２１６ａ，２１６ｂ）と，遅延演算部２１７とが備わっている。

メモリ２１３ｆには，端末１０１ｂの内部の処理遅延値（内部処理遅延値）が格納され，保持される。内部処理遅延値は，例えば，パケットを受信し，新たなパケットを生成し，そのパケットを送信する一連の端末１０１ｂの内部処理で生じる遅延時間である。

遅延演算部２１７は，端末１０１ｂの内部に備わる例えばアナログ回路やディジタル回路で実行される処理で発生する遅延時間を合計し，上記内部処理遅延値を取得する。求められた上記内部処理遅延値は，遅延演算部２１７によって，メモリ２１３ｆに格納される。

パケット受信部２１５は，パケットを受信し，そのうちの距離測定用のパケット１を検出する。なお，パケット１については後程説明する。

パケット生成／送信部２１６ａは，ＭＵＳＩＣアルゴリズム処理部２１２で高分解能のインパルスレスポンスが求められた結果から距離測定用のパケット１の受信時刻補正値（パケット１受信時刻補正値）を求める。なお，パケット１の受信時刻補正値は，端末１０１ｂがパケット１を受信した時刻の精度，正確さを高めるための補正値である。

さらに，パケット生成／送信部２１６ａは，メモリ２１３ｆに格納された内部処理遅延値を含んだパケット３を生成し，無線通信ブロック２０３に送信する。

パケット生成／送信部２１６ｂは，上記パケット受信部２１５でパケット１を受信したタイミングと同時に，すぐにパケット２を生成し，無線通信ブロック２０３に送信する。なお，上記パケット２が端末１０１ａに到達することで，大まかなターンアラウンドタイムが求まり，端末１０１間の大まかな距離が求められる。

また，第１の実施の形態にかかる端末１０１ａには複数のメモリ２１３（２１３ａ〜２１３ｅ）が備わる場合を例に挙げて説明したが，かかる例に限定されず，例えば，端末１０１ａには単一のメモリ２１３のみが備わる場合でも実施可能である。

また，第１の実施の形態にかかる端末１０１ｂに備わる無線通信ブロック２０３は，上記説明の端末１０１ａに備わる無線通信ブロック２０３と実質的に同一であるため詳細な説明は省略する。

なお，第１の実施の形態にかかる端末１０１ａには，距離測定ブロック２０１ｂが備わらない場合を例に挙げて説明したが，かかる例に限定されず，例えば，端末１０１ａにさらに距離測定ブロック２０１ｂが備わる場合でも実施可能である。

なお，第１の実施の形態にかかる端末１０１ｂには，距離測定ブロック２０１ａが備わらない場合を例に挙げて説明したが，かかる例に限定されず，例えば，端末１０１ｂにさらに距離測定ブロック２０１ａが備わる場合でも実施可能である。

また，第１の実施の形態にかかる端末１０１ｂには，パケット生成／送信部２１６ａとパケット生成／送信部２１６ｂとが別体として備わる場合を例に挙げて説明したが，かかる例に限定されず，例えば，端末１０１ｂに一体化された単一のパケット生成／送信部２１６が備わる場合等でも実施することができる。

上記説明したように，図２及び図３に示す端末１０１の構成によれば，距離測定する側（端末１０１ａ）と，距離測定される側（端末１０１ｂ）の双方にＭＵＳＩＣアルゴリズム処理部２１２を搭載することで，高精度なパケット位置の特定を行うことができる。

上記高精度なパケット位置を特定するために，端末１０１ｂに備わるＦＦＴ２２４によるＦＦＴ処理後の周波数領域のデータに対して，高分解能のスペクトル推定アルゴリズムの一つであるＭＵＳＩＣアルゴリズムを処理するＭＵＳＩＣアルゴリズム２１２が高分解能のインパルスレスポンスを得ることによって特定される。

周波数領域のデータは，通常，互いに相関が高いため，周波数平均という前処理を行う必要がある。移動平均処理部２１１が上記前処理に相当する。これは，ＭＵＳＩＣアルゴリズムを適用する場合の前処理において一般的なものである。

無線ＬＡＮの場合，周波数領域のデータには，さらにＤＣ成分を使用しない目的のために，通常，キャリアホールが存在する。当該キャリアホールのキャリアでは，エネルギーが０である。つまり，ＭＵＳＩＣアルゴリズム処理部２１２による処理が実行される前の前処理の周波数平均（移動平均）を処理実行するには，予めこのキャリアホールの位置の周波数データを隣のキャリアデータから補間する必要がある。

したがって，図２及び図３に示すように，端末１０１には，補間処理部２１０と移動平均処理部２１１が備わり，ＭＵＳＩＣアルゴリズム処理部２１２の処理前に，補間処理部２１０と移動平均処理部２１１による処理が実行される。

以上から，従来に係るＭＵＳＩＣアルゴリズムは主としてアダプティブアレーの伝播解析処理等に適用されていたが，そのまま単にＭＵＳＩＣアルゴリズムを無線ＬＡＮにおける測距システムに適用したとしても高精度の測距を実現することができなかった。

次に，図４を参照しながら，第１の実施の形態にかかる測距システム１００の一連の処理について説明する。図４は，第１の実施の形態にかかる測距システムの一連の処理の概略の一例を示すシーケンス図である。なお，距離測定を行う端末を端末１０１ａ（送信機），距離測定をされる端末を端末１０１ｂ（受信機）とするが，かかる例に限定されない。

図４に示すように，まず，端末１０１ａは，端末１０１ａと端末１０１ｂとの間の距離を測定するため距離測定用のパケット１を送信すると同時に，そのパケット１の送信時刻を示すパケット１送信時刻情報をメモリ２１３ａに格納する（Ｓ４０１）。

次に，端末１０１ｂは，パケットを受信し，そのうち距離測定用のパケット１を認識すると同時に，パケット２の生成処理を開始し，可能な限りすぐに，距離測定用のパケット２を送信する（Ｓ４０３）。パケット２の送信処理の際，ＤＡ変換前等のディジタル処理による処理遅延や，ＤＡ変換後等のアナログ処理による処理遅延を合計した遅延時間を示す値（内部処理遅延値ｔ_１）をメモリ２１３ｆに記憶する（Ｓ４０３）。

次に，端末１０１ｂは，端末１０１ａからのパケット１を受信し，そのパケット１は無線通信ブロック２０３におけるＦＦＴ２２４によってＦＦＴ処理される。次にＦＦＴ処理後の周波数データのキャリアホールを補間処理し，その補間処理後のデータについて周波数平均処理を行った後に，ＭＵＳＩＣアルゴリズムを適用して高分解能のインパルスレスポンスを求める。さらに，上記インパルスレスポンスに基づいてパケット１受信時刻補正値が求められる（Ｓ４０５）。上記パケット１受信時刻補正値は，メモリ２１３に記憶する。なお，パケット１受信時刻補正値は，メモリ２１３ｆに格納されてもよく，または端末１０１ｂに備わるその他の記憶手段（図示せず。）に格納される場合でもよい。

次に，端末１０１ｂに備わるパケット生成／送信部２１６ａは，上記内部処理遅延値とパケット１受信時刻補正値とを含んだパケット３を生成し，そのパケット３を無線通信ブロック２０３に送信する（Ｓ４０７）。上記無線通信ブロック２０３は，所定処理実行後，ＲＦ処理部２２０から無線で端末１０１ａに送信される。

次に，端末１０１ａがパケット２を受信すると，そのパケット２を受信した時刻を示すパケット２受信時刻情報がメモリ２１３ｂに記憶される（Ｓ４０９）。

次に，端末１０１ａは，パケット２を受信し，そのパケット２はＦＦＴ２２４によってＦＦＴ処理される。次にＦＦＴ処理後の周波数データのキャリアホールを補間処理し，その補間処理後のデータについて周波数平均処理を行った後に，ＭＵＳＩＣアルゴリズムを適用し，高分解能のインパルスレスポンスが求められる。さらに上記インパルスレスポンスに基づいてパケット２受信時刻補正値が求められ，メモリ２１３ｃに記憶される（Ｓ４１１）。

次に，端末１０１ａは，端末１０１ｂにより送信したパケット３を受信する。そして端末１０１ａはパケット３から，端末１０１ｂで求められたパケット１受信時刻補正値と，内部処理遅延値ｔ_１を取得し，パケット２受信時刻補正値保持と内部処理遅延値ｔ_１とはメモリ２１３ｅと，メモリ２１３ｄに各々格納される（Ｓ４１３）。

次に，端末１０１ａに備わる距離計算処理部（往復時間概算部）２１４は，メモリ２１３ａに記憶されたパケット１送信時刻情報と，メモリ２１３ｂに記憶されたパケット２受信時刻情報とを取得し，その双方の情報から，ターンアランウドタイムの概算をする（Ｓ４１５）。

なお，第１の実施の形態にかかるターンアラウンドタイムの概算は，パケット２を受信した受信時刻からパケット１を送信した送信時刻を引き算することで求められる。つまり，パケット２受信時刻情報−パケット送信時刻情報＝ターンアラウンドタイムの概算値（パケット往復時間）となる。

次に，端末１０１ａに備わる距離計算処理部２１４は，上記ターンアラウンドタイムの概算値から，さらに端末１０１ｂ（通信相手）の内部処理遅延値ｔ_１を引き算することにより，往復の電波の伝搬時間（伝搬往復時間）の概算値を求める（Ｓ４１７）。なお，伝搬往復時間の概算値は，“伝搬往復時間＝ターンアラウンドタイムの概算値−通信相手の内部処理遅延値”の式から求まる。

次に，端末１０１ａに備わる距離計算処理部（伝搬時間補正部）２１４は，まずメモリ２１３ｃに記憶されたパケット２受信時刻補正値と，メモリ２１３ｅに記憶されたパケット１受信時刻補正値を取得し，上記双方の補正値に基づいて，Ｓ４１７で求められた伝搬往復時間の概算値を補正し，高精度な伝搬往復時間（補正後伝搬往復時間）を得る。さらに，距離計算処理部（距離算出部）２１４は，上記高精度な伝搬往復時間と，電波の速度（３．０×１０^８ｍ／ｓ）とから端末１０１間の距離を計算する（Ｓ４１９）。以上で，第１の実施の形態にかかる測距システム１００の一連の動作の説明が終了する。

次に，第１の実施の形態にかかる測距システム１００の実施例を示す。なお，以下の実施例では，５ＧＨｚのＯＦＤＭ（直交周波数分割多重）変調方式のうち，４０ＭＨｚ帯域幅を使用したＯＦＤＭ変調方式の場合を例に挙げて説明するが，マルチキャリア伝送方式であれば，かかる例に限定されない。

さらに，測距システム１００の実施例において，ＦＦＴ２２４のポイント数は１２８ポイントのものを使用し，サブキャリアは１１４本とする。またＤＣ成分の近辺には３本のキャリアホールが存在する。上記の条件下において，以下に測距システム１００の実施例を示す。

まず，上記ＤＣ成分の近辺に３本のキャリアホールが存在するため，上記説明した通り，３本のキャリアホールは隣の又は近辺のキャリア値を用いることで補間する。したがって，合計で１１７本の周波数データに対してＭＵＳＩＣアルゴリズムを適用する。なお，ＭＵＳＩＣアルゴリズムを適用する前に，上記説明の通り１１７本の周波数データの平均処理を事前に行っておく。

なお，ＭＵＳＩＣアルゴリズムを適用することで求められるインパルスレスポンスの対象を図５に示すように直接波が５ｎｓにくるようなインパルスレスポンスとする。図５は，第１の実施の形態にかかる直接波が５ｎｓになる場合のインパルスレスポンスの概略を示す説明図である。

端末１０１に備わるパケット同期処理部２２２によって判別されたパケットの同期位置を使用し，ＦＥＴ処理後のデータの切り出し位置を図６に示すように実行した場合，ＭＵＳＩＣアルゴリズムを適用した結果のＭＵＳＩＣスペクトラムは図７に示すようになる。なお，図６は，データの切り出し位置の概略を示す説明図であり，図７は，ＭＵＳＩＣスペクトラムの概略的な構成の一例を示す説明図である。

なお，図６に示すガードインターバルを設けることで，例えば，乱反射などによって受信地点に時間的なずれを持った信号（ゴースト）が到来しても，影響（マルチパス障害）が出ないようにすることができる。

図７に示すように，ＭＵＳＩＣスペクトラムの最大のピークは略５ｎｓであることが分かる。上記図６に示すデータの切り出し位置では，ＭＵＳＩＣスペクトラムは５ｎｓでピークとなっているのが分かるが，かかる切り出し位置で切り出すのは必ずしも容易であるとは限らず，切り出し位置がずれる場合は多々ある。

そこで，データの切り出しに多少余裕を持たせるために，例えば，ＦＥＴ処理後のデータの切り出し位置を１サンプル分（２５ｎｓ）前にずらした場合，データの切り出し位置は図８に示すようになる。なお，図８は，データを２５ｎｓだけ前にずらして切り出した場合の切り出し位置の概略を示す説明図である。

その図８に示すデータの切り出し位置で，データを切り出し，ＭＵＳＩＣアルゴリズムを適用した場合のＭＵＳＩＣスペクトラムは図９に示すようになる。なお，図９は，ＭＵＳＩＣスペクトラムの概略的な構成の一例を示す説明図である。

図９に示すように，ＭＵＳＩＣスペクトラムのピークは，略３０ｎｓであり，これは図７に示すＭＵＳＩＣスペクトラムのピーク（５ｎｓ）に切り出し位置をずらした分（２５ｎｓ）を加算した値となる。つまり，切り出し位置のずらす程度と同じだけ，ＭＵＳＩＣスペクトラムのピークがずれることとなる。なお，上記ずれる分が予め分かれば，ＭＵＳＩＣスペクトラムを補正し，例えば５ｎｓにすることができ，高精度の測距を実現できる。以上から，ＯＦＤＭ変調方式においてＦＦＴ処理後のデータにＭＵＳＩＣアルゴリズムを適用することによって，パケットの送信／受信位置をシステムクロックの幅よりも細かい精度で正確に特定することができる。つまり，端末１０１間の距離をより正確に測定することができる。

一方，ＵＷＢ（ウルトラワイドバンド通信方式）の無線方式では，高速なパルスを使用するため，例えば，２ＧＨｚでのパケットの同期をとることができる。２ＧＨｚのクロック幅は，１／２ＧＨｚ＝０．５ｎｓであるため，３．０×１０^８ｍ／ｓ×０．５ｎｓ＝０．１５ｍの精度で測距を行うことができる。本実施例では，このＵＷＢまで高速なクロックではないにも関わらず，高精度の測距を実現している。

最後に，図１０に測距システム１００で距離を測定した結果の精度について示す。なお，ＳＮ＝２０ｄＢとする。図１０を参照すれば，測距した結果の「距離の誤差」は，試行回数がいずれの場合も０．０９ｃｍ〜０．２７ｃｍの範囲内（平均は，０．１４４ｃｍ）であり，極めて距離の誤差を最小化することができる。つまり，メートルオーダーではなくセンチメートルオーダーの測距を実行することができる。なお設定時刻は，想定するインパルスレスポンスの対象となる値（５ｎｓ）である。

なお，従来にかかるＭＵＳＩＣアルゴリズムについてもアダプティブアレーの伝播解析処理等には適用されていた。しかしながら，無線ＬＡＮを利用した測距システム１００に組み合わせて適用することは，容易に想到することはできなかった。それは，ＦＦＴ処理後のデータをそのままＭＵＳＩＣアルゴリズムに適用したとしても，上記補間処理部２１０による補間処理等を実行しない限り，高精度の測距は実現できないためである。したがって，無線ＬＡＮを利用した測距システム１００に上記ＭＵＳＩＣアルゴリズムを適用することは当業者が容易に想到できたものではない。

（第２の実施の形態について）
次に，図１１，図１２を参照しながら，第２の実施の形態にかかる端末１０１について説明する。図１１は，第２の実施の形態にかかる距離を測定する側の端末（端末１０１ｃ）の概略的な構成を示すブロック図であり，図１２は，第２の実施の形態にかかる距離を測定される側の端末（端末１０１ｄ）の概略的な構成を示すブロック図である。

なお，第２の実施の形態にかかる測距システムは，図１を参照しながら説明した第１の実施の形態にかかる測距システム１００と実質的に同一であるため詳細な説明は省略する。また，第２の実施の形態にかかる端末１０１についての説明は，上記説明した第１の実施の形態にかかる端末１０１と相違する点について詳細に説明し，その他の点については実質的に同一であるため詳細な説明は省略する。

第２の実施の形態にかかる測距システムでは，距離測定する側の端末（端末１０１ｃ）だけがＭＵＳＩＣアルゴリズムを適用することで端末１０１間の距離を測定する。なお，第１の実施の形態にかかる測距システムでは，距離測定する側及び距離測定される側の端末１０１がＭＵＳＩＣアルゴリズムを適用し，端末１０１間の距離を測定する。

第２の実施の形態にかかる距離測定をする側の端末１０１ｃは，図２に示す第１の実施の形態にかかる距離測定をする側の端末１０１ａと比較し，距離測定をする側の端末１０１ｄから受信するパケット３にはパケット１受信時刻補正値が含まれておらず，上記端末１０１ｄからパケット１を受信時の周波数領域のデータを受信し，端末１０１ｃがパケット１受信時刻補正値を求める点で相違する。

より具体的には，図１１に示すように，第２の実施の形態にかかる端末１０１ｃには，端末１０１ｃが求めたパケット１受信時刻補正値を記憶するメモリ２１３ｇと，メモリ２１３ｅの代わりに端末１０１ｄでパケット１を受信した時の周波数領域のデータを記憶するメモリ２１３ｈが備わっている点で，第１の実施の形態にかかる端末１０１ａと相違する。

なお，第２の実施の形態にかかる端末１０１ｃに備わるメモリ２１３ｇ及び／又はメモリ２１３ｈは，他のメモリ２１３（２１３ａ，２１３ｂ，２１３ｃ，２１３ｄ）と別体である場合を例に挙げて説明したが，かかる例に限定されず，例えば，メモリ２１３ｇ及び／又はメモリ２１３ｈは，他のメモリ２１３（２１３ａ，２１３ｂ，２１３ｃ，２１３ｄ）のうち少なくとも一つのメモリ２１３と一体化される場合等でも実施可能である。

次に，第２の実施の形態にかかる距離測定される側の端末１０１ｄは，第１の実施の形態にかかる端末１０１ｂと比較して，ＭＵＳＩＣアルゴリズム処理部２１２と，移動平均処理部２１１と，補間処理部２１０とを備えず，パケット生成／送信部２１６ａの代わりにパケット生成／送信部２１６ｃを備える点で相違する。

上記パケット生成／送信部２１６ｃは，端末１０１ｄがパケット１を受信した時の周波数領域のデータと，端末１０１ｄの内部遅延処理値とが少なくとも含まれるパケット３を生成し，そのパケット３を無線通信ブロック２０３に送信する。その他の点については，第１の実施の形態にかかる端末１０１ｂと実質的に同じである。

つまり，図１２に示す第２の実施の形態にかかる端末１０１ｄは，パケット３で端末１０１ｃに送信する情報としては，パケット１受信時刻補正値ではなく，パケット１を受信した時の周波数領域のデータをそのまま加工しないで送信する。

また，図１１に示す第２の実施の形態にかかる端末１０１ｃは，端末１０１ｄから上記パケット３を受信すると，パケット３に含まれるパケット１を受信した時の周波数領域のデータと，内部処理遅延値を取得し，上記周波数領域のデータを用いてＭＵＳＩＣアルゴリズムを適用することで，パケット１受信時刻補正値を求める。

したがって，第１の実施の形態にかかる端末１０１ｂでは，パケット１受信時刻補正値をＭＵＳＩＣアルゴリズムを適用することで求めていたが，第２の実施の形態にかかる端末１０１ｄでは，パケット１受信時刻補正値を求めずに，パケット１の受信時の周波数データを送信すればよいため，端末１０１ｄに対する処理負荷が軽減され，端末１０１ｄに備わる回路等を省略化できる。

また，第２の実施の形態にかかる測距システムの一連の動作は，第１の実施の形態にかかる測距システムの一連の動作と比較して，距離を測定する側の端末１０１ｃがパケット１受信時刻補正値とパケット２受信時刻補正値のパケット受信時刻補正値全てを求めている点で相違するが，その他の点については第１の実施の形態にかかる測距システムとほぼ同様であるため詳細な説明を省略する。以上で，第２の実施の形態にかかる測距システムの説明を終了する。

次に，第１の実施の形態及び第２の実施の形態で説明した一連の測距システムに関する処理は，専用のハードウェアにより行うこともできるし，ソフトウェアにより行うこともできる。一連の処理をソフトウェアによって行う場合には，そのソフトウェアを構成するプログラムが，汎用のコンピュータやマイクロコンピュータ等に相当する端末１０１に備わるメモリ２１３等の記憶手段にインストールされる。

なお，測距システムの動作を端末１０１に実行させるためのプログラムは，コンピュータに内蔵されている記録手段としてのＨＤＤ（ハードディスクドライブ）やＲＯＭ，ＲＡＭ等に予め記録しておくことができる。

あるいはまた，プログラムは，ＨＤＤに限らず，フレキシブルディスク，ＣＤ−ＲＯＭ（ＣｏｍｐａｃｔＤｉｓｃＲｅａｄＯｎｌｙＭｅｍｏｒｙ），ＭＯ（ＭａｇｎｅｔｏＯｐｔｉｃａｌ）ディスク，ＤＶＤ（ＤｉｇｉｔａｌＶｅｒｓａｔｉｌｅＤｉｓｃ），磁気ディスク，半導体メモリなどのリムーバブル記録媒体に，一時的あるいは永続的に格納（記録）しておくことができる。このようなリムーバブル記録媒体は，いわゆるパッケージソフトウエアとして提供することができる。

なお，プログラムは，上述したようなリムーバブル記録媒体からコンピュータにインストールする他，ダウンロードサイトから，ディジタル衛星放送用の人工衛星を介して，端末１０１に無線で転送したり，ＬＡＮ（ＬｏｃａｌＡｒｅａＮｅｔｗｏｒｋ），インターネットといったネットワークを介して，端末１０１に有線で転送してもよい。端末１０１では，上記のように転送されてくるプログラムを，受信し，内蔵するＨＤＤ等にインストールされる。

ここで，本明細書において，端末１０１に各種の処理を行わせるためのプログラムを記述する処理ステップは，必ずしも図４に示すシーケンス図として記載された順序に沿って時系列に処理する必要はなく，並列的あるいは個別に実行される処理等も含むものである。

最後に，従来にかかるＭＵＳＩＣアルゴリズムについてもアダプティブアレーの伝播解析処理等には適用されていた。しかしながら，無線ＬＡＮを利用した本明細書に記載の測距システムに組み合わせて適用することは，容易に想到することはできなかった。それは，ＦＦＴ処理後のデータをそのままＭＵＳＩＣアルゴリズムに適用したとしても，上記補間処理部２１０による補間処理等を実行しない限り，高精度の測距は実現できないためである。したがって，無線ＬＡＮを利用した測距システム１００に上記ＭＵＳＩＣアルゴリズムを適用することは当業者が容易に想到できたものではない。

以上，添付図面を参照しながら本発明の好適な実施形態について説明したが，本発明はかかる例に限定されない。当業者であれば，特許請求の範囲に記載された技術的思想の範疇内において各種の変更例または修正例を想定し得ることは明らかであり，それらについても当然に本発明の技術的範囲に属するものと了解される。

上記実施形態においては，測距システム１００ではＭＵＳＩＣアルゴリズムを適用して端末間の距離を測定する場合を例にあげて説明したが，スペクトル推定アルゴリズムであれば，本発明はかかる例に限定されず如何なるアルゴリズムでも実施することができる。

また上記実施形態においては，メモリ２１３は，複数のメモリから構成される場合を例に挙げて説明したが，かかる例に限定されず，例えば，メモリ単体から構成される場合や，ハードディスクドライブ単体から構成される場合等でもよい。

また上記実施形態においては，端末１０１ｂには，パケット生成／送信部２１６ａとパケット生成／送信部２１６ｂとが別体として備わる場合を例に挙げて説明したが，かかる例に限定されず，例えば，端末１０１ｂに一体化された単一のパケット生成／送信部２１６が備わる場合等でも実施することができる。

上記実施形態においては，端末１０１に備わる各部はハードウェアからなる場合を例にあげて説明したが，本発明はかかる例に限定されない。例えば，上記各部は，１又は２以上のモジュールまたはコンポーネントから構成されるプログラムの場合であってもよい。

本発明は，マルチキャリア伝送方式で測距可能な測距システム，送信端末，受信端末，測距方法，およびコンピュータプログラムに適用可能である。

第１の実施の形態にかかる測距システムの概略的な構成の一例を示す説明図である。第１の実施の形態にかかる距離測定をする側の端末の概略的な構成の一例を示すブロック図である。第１の実施の形態にかかる距離測定される側の端末の概略的な構成の一例を示すブロック図である。第１の実施の形態にかかる測距システムの一連の処理の概略の一例を示すシーケンス図である。第１の実施の形態にかかる直接波が５ｎｓになる場合のインパルスレスポンスの概略を示す説明図である。データの切り出し位置の概略を示す説明図である。ＭＵＳＩＣスペクトラムの概略的な構成の一例を示す説明図である。データを２５ｎｓだけ前にずらして切り出した場合の切り出し位置の概略を示す説明図である。ＭＵＳＩＣスペクトラムの概略的な構成の一例を示す説明図である。測距システムで距離を測定した結果の精度の概略を示す説明図である。第２の実施の形態にかかる距離を測定する側の端末（端末１０１ｃ）の概略的な構成を示すブロック図である。第２の実施の形態にかかる距離を測定される側の端末（端末１０１ｄ）の概略的な構成を示すブロック図である。

符号の説明

１０１端末
２１０補間処理部
２１１移動平均処理部
２１２ＭＵＳＩＣアルゴリズム処理部
２１３メモリ
２１４距離計算処理部
２１８時刻設定部
２２４ＦＦＴ

Claims

マルチキャリア伝送方式における端末間の距離を測定する測距システムであって：
パケットの送信時及び／又は受信時のデータが高速離散フーリエ変換され，さらに所定処理後の周波数領域のデータに対して，スペクトル推定アルゴリズムを適用することで，インパルスレスポンスを取得する周波数領域データ処理部と；
を備えることを特徴とする，測距システム。
前記インパルスレスポンスは，前記パケットの受信時刻及び／又は送信時刻を高精度に示す時刻情報であることを特徴とする，請求項１に記載の測距システム。
前記測距システムには，前記周波数領域のデータに存在するキャリアホールを補間する補間処理部と；
前記補間後の周波数領域のデータに対して周波数平均処理する周波数平均処理部とを備え，
前記所定処理には，前記補間処理部による補間処理と，前記周波数平均処理部による周波数平均処理とが含まれることを特徴とする，請求項１に記載の測距システム。
前記補間処理部は，前記キャリアホールに相当する部分のデータを隣のキャリアからのデータで補間することを特徴とする，請求項３に記載の測距システム。
マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末であって：
第１のパケットを通信相手の受信端末に送信するパケット送信部と；
前記第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成部とを備え，
さらに，前記第１のパケット送信時刻情報と，前記受信端末により前記第１のパケットの受信とともに前記送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，前記受信端末と送信端末との間のパケット往復時間を概算する往復時間概算部と；
前記受信端末からの第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理部と；
前記第１の周波数領域のデータに所定処理を施した後，該所定処理後の第１の周波数領域のデータに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求める周波数領域データ処理部と；
第３のパケットを前記受信端末から所定のタイミングで受信し，前記パケット往復時間から前記第３のパケットに含む内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに前記第１のインパルスレスポンスと前記第３のパケットに含む第２のインパルスレスポンスとに基づいて，前記電波の伝搬往復時間を補正後伝搬往復時間に補正する伝搬時間補正部と；
前記補正後伝搬往復時間と電波の速度とから前記送信端末と受信端末との間の距離を求める距離算出部とを備えることを特徴とする，送信端末。
前記第２のインパルスレスポンスは，前記受信端末によって前記第１のパケット受信時のデータから第２の周波数領域のデータに高速離散フーリエ変換され，所定処理が施された後，その第２の周波数領域のデータにスペクトル推定アルゴリズムが適用されることで求められ，前記内部処理遅延時間は，該受信端末の内部処理の遅延時間であることを特徴とする，請求項５に記載の送信端末。
前記第１のインパルスレスポンスは，前記第２のパケットの受信時刻を高精度に示す時刻情報であって，前記第２のインパルスレスポンスは，前記第１のパケットの受信時刻を高精度に示す時刻情報であることを特徴とする，請求項５に記載の送信端末。
前記送信端末は，前記周波数領域データに存在するキャリアホールを補間する補間処理部と；
前記補間後の周波数領域データに対して周波数平均処理する周波数平均処理部とを備え，
前記所定処理には，前記補間処理部による補間処理と，前記周波数平均処理部による周波数平均処理とが含まれることを特徴とする，請求項５に記載の送信端末。
前記補間処理部は，前記キャリアホールに相当する部分のデータを隣のキャリアからのデータで補間することを特徴とする，請求項８に記載の送信端末。
マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末の測距方法であって：
第１のパケットを通信相手の受信端末に送信するパケット送信ステップと；
前記第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成ステップとを含み，
さらに，前記第１のパケット送信時刻情報と，前記受信端末により前記第１のパケットの受信とともに前記送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，前記受信端末と送信端末との間のパケット往復時間を概算する往復時間概算ステップと；
前記受信端末からの第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換ステップと；
前記第１の周波数領域のデータに所定処理を施した後，該所定処理後の第１の周波数領域のデータに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求めるステップと；
第３のパケットを前記受信端末から所定のタイミングで受信し，前記パケット往復時間から前記第３のパケットに含む内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに前記第１のインパルスレスポンスと前記第３のパケットに含む第２のインパルスレスポンスとに基づいて，前記電波の伝搬往復時間を補正後伝搬往復時間に補正する伝搬時間補正ステップと；
前記補正後伝搬往復時間と電波の速度とから前記送信端末と受信端末との間の距離を求める距離算出ステップとを含むことを特徴とする，測距方法。
コンピュータをして，マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末として機能させるコンピュータプログラムであって：
第１のパケットを通信相手の受信端末に送信するパケット送信手段と；
前記第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成手段とを含み，
さらに，前記第１のパケット送信時刻情報と，前記受信端末により前記第１のパケットの受信とともに前記送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，前記受信端末と送信端末との間のパケット往復時間を概算する往復時間概算手段と；
前記受信端末からの第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理手段と；
前記第１の周波数領域のデータに所定処理を施した後，該所定処理後の第１の周波数領域のデータに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求める周波数領域データ処理手段と；
第３のパケットを前記受信端末から所定のタイミングで受信し，前記パケット往復時間から前記第３のパケットに含む内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに前記第１のインパルスレスポンスと前記第３のパケットに含む第２のインパルスレスポンスとに基づいて，前記電波の伝搬往復時間を補正後伝搬往復時間に補正する伝搬時間補正手段と；
前記補正後伝搬往復時間と電波の速度とから前記送信端末と受信端末との間の距離を求める距離算出手段とを含むことを特徴とする，コンピュータプログラム。
マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる受信端末であって：
通信相手の送信端末からの第１のパケットを受信するとともに，第２のパケットを生成し，前記送信端末に該第２のパケットを送信する第１のパケット送信部と；
前記送信端末からの第１のパケット受信時のデータを周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理部と；
前記周波数領域のデータに所定処理を施した後，該所定処理後の周波数領域のデータに対してスペクトル推定アルゴリズムを適用し，インパルスレスポンスを求める周波数領域データ処理部と；
前記受信端末自体の内部処理による処理遅延時間として内部処理遅延時間を求める処理遅延算出部と；
前記インパルスレスポンスと前記内部処理遅延時間とが少なくとも含まれた第３のパケットを生成し，所定のタイミングで，前記送信端末に送信する第２のパケット送信部とを備えることを特徴とする，受信端末。
前記インパルスレスポンスは，前記第１のパケットの受信時刻を高精度に示す時刻情報であることを特徴とする，請求項１２に記載の受信端末。
前記受信端末は，前記周波数領域データに存在するキャリアホールを補間する補間処理部と；
前記補間後の周波数領域データに対して周波数平均処理する周波数平均処理部とを備え，
前記所定処理には，前記補間処理部による補間処理と，前記周波数平均処理部による周波数平均処理とが含まれることを特徴とする，請求項１２に記載の受信端末。
前記補間処理部は，前記キャリアホールに相当する部分のデータを隣のキャリアからのデータで補間することを特徴とする，請求項１４に記載の受信端末。
マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末であって：
第１のパケットを通信相手の受信端末に送信するパケット送信部と；
前記第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成部とを備え，
さらに，前記第１のパケット送信時刻情報と，前記受信端末により前記第１のパケットの受信とともに前記送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，前記受信端末と該送信端末との間のパケット往復時間を概算する往復時間概算部と；
前記受信端末からの前記第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理部と；
第３のパケットを前記受信端末から所定のタイミングで受信し，該第３のパケットに含まれる第２の周波数領域データに所定処理を施した後，該所定処理後の第２の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求めるとともに，前記第１の周波数領域データに所定処理を施した後，該所定処理後の第１の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第２のインパルスレスポンスを求める周波数領域データ処理部と；
前記パケット往復時間から前記第３のパケットに含む内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに前記第１のインパルスレスポンスと前記第２のインパルスレスポンスとに基づいて，前記電波の伝搬往復時間を補正後伝搬往復時間に補正する伝搬時間補正部と；
前記補正後伝搬往復時間と電波の速度とから前記送信端末と受信端末との間の距離を求める距離算出部とを備えることを特徴とする，送信端末。
前記第２の周波数領域のデータは，前記受信端末によって前記第１のパケット受信時のデータから第２の周波数領域のデータに高速離散フーリエ変換され，前記内部処理遅延時間は，該受信端末の内部処理の遅延時間であることを特徴とする，請求項１６に記載の送信端末。
前記第１のインパルスレスポンスは，前記第１のパケットの受信時刻を高精度に示す時刻情報であって，前記第２のインパルスレスポンスは，前記第２のパケットの受信時刻を高精度に示す時刻情報であることを特徴とする，請求項１６に記載の送信端末。
前記送信端末は，前記周波数領域データに存在するキャリアホールを補間する補間処理部と；
前記補間後の周波数領域データに対して周波数平均処理する周波数平均処理部とを備え，
前記所定処理には，前記補間処理部による補間処理と，前記周波数平均処理部による周波数平均処理とが含まれることを特徴とする，請求項１６に記載の送信端末。
前記補間処理部は，前記キャリアホールに相当する部分のデータを隣のキャリアからのデータで補間することを特徴とする，請求項１９に記載の送信端末。
マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末の測距方法であって：
第１のパケットを通信相手の受信端末に送信するパケット送信ステップと；
前記第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成ステップとを含み，
さらに，前記第１のパケット送信時刻情報と，前記受信端末により前記第１のパケットの受信とともに前記送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，前記受信端末と該送信端末との間のパケット往復時間を概算する往復時間概算ステップと；
前記受信端末からの前記第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理ステップと；
第３のパケットを前記受信端末から所定のタイミングで受信し，該第３のパケットに含まれる第２の周波数領域データに所定処理を施した後，該所定処理後の第２の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求めるとともに，前記第１の周波数領域データに所定処理を施した後，該所定処理後の第１の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第２のインパルスレスポンスを求めるステップと；
前記パケット往復時間から前記第３のパケットに含む内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに前記第１のインパルスレスポンスと前記第２のインパルスレスポンスとに基づいて，前記電波の伝搬往復時間を補正後伝搬往復時間に補正する伝搬時間補正ステップと；
前記補正後伝搬往復時間と電波の速度とから前記送信端末と受信端末との間の距離を求める距離算出ステップとを含むことを特徴とする，測距方法。
コンピュータをして，マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる送信端末として機能させるコンピュータプログラムであって：
第１のパケットを通信相手の受信端末に送信するパケット送信手段と；
前記第１のパケットを送信した時刻を示す第１のパケット送信時刻情報を少なくとも生成する時刻情報生成手段とを含み，
さらに，前記第１のパケット送信時刻情報と，前記受信端末により前記第１のパケットの受信とともに前記送信端末に送信される第２のパケットを受信した時刻を示す第２のパケット受信時刻情報と，に基づき，前記受信端末と該送信端末との間のパケット往復時間を概算する往復時間概算手段と；
前記受信端末からの前記第２のパケット受信時のデータを第１の周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理手段と；
第３のパケットを前記受信端末から所定のタイミングで受信し，該第３のパケットに含まれる第２の周波数領域データに所定処理を施した後，該所定処理後の第２の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第１のインパルスレスポンスを求めるとともに，前記第１の周波数領域データに所定処理を施した後，該所定処理後の第１の周波数領域データに対してスペクトル推定アルゴリズムを適用し，第２のインパルスレスポンスを求める周波数領域データ処理手段と；
前記パケット往復時間から前記第３のパケットに含む内部処理遅延時間を差し引くことで電波の往復の伝搬時間を概算し，さらに前記第１のインパルスレスポンスと前記第２のインパルスレスポンスとに基づいて，前記電波の伝搬往復時間を補正後伝搬往復時間に補正する伝搬時間補正手段と；
前記補正後伝搬往復時間と電波の速度とから前記送信端末と受信端末との間の距離を求める距離算出手段とを含むことを特徴とする，コンピュータプログラム。
マルチキャリア伝送方式における端末間の距離を測定する測距システムに備わる受信端末であって：
通信相手の送信端末からの第１のパケットを受信するとともに，第２のパケットを生成し，前記送信端末に該第２のパケットを送信する第１のパケット送信部と；
前記送信端末からの第１のパケット受信時のデータを周波数領域のデータに高速離散フーリエ変換する高速離散フーリエ変換処理部と；
前記受信端末自体の内部処理による処理遅延時間として内部処理遅延時間を求める処理遅延算出部と；
前記周波数領域のデータと前記内部処理遅延時間とが少なくとも含まれた第３のパケットを生成し，所定のタイミングで，前記送信端末に送信する第２のパケット送信部とを備えることを特徴とする，受信端末。