(57)【要約】【課題】ＭＭＴは、受信信号が直接波１波とマルチパス１波とからなる場合において、マルチパスの影響を低減できる。しかしながら、実際の環境において、マルチパス波の数は１波とは限らない。受信信号に複数のマルチパス波が含まれている場合には、マルチパスの影響を除去しきれない。反対にマルチパス波が含まれない場合には、直接波１波を直接波１波とマルチパス１波として推定してしまうため、誤差を生じてしまう。また、時間領域信号を対象に最尤推定を行うことは計算上の困難を伴う。【解決手段】本発明に係る測位装置は、周波数領域での最尤推定により信号モデルのパラメータを推定し、情報量基準により信号モデルを推定することによって、マルチパス誤差を抑制する。【選択図】図２
2008216111.tif 000002
【特許請求の範囲】
【請求項１】
送信機から送信された信号を受信する信号受信手段と、
前記信号受信手段で受信した前記信号を推定する信号推定手段と、
前記信号推定手段によって推定した前記信号から前記信号の伝播遅延時間を算出する伝播遅延時間算出手段と、
前記伝播遅延時間から前記送信機と前記信号受信手段との測距をする測距手段とを備え、
前記信号推定手段は、信号モデルのパラメータを周波数領域で最尤推定したことを特徴とする測距装置。
【請求項２】
信号推定手段は、受信信号を複数のマルチパス波を含む信号として推定する場合には１波少ないマルチパス波を含む信号として推定した結果を利用したことを特徴とする請求項１に記載の測距装置。
【請求項３】
信号推定手段は、受信信号を複数のマルチパス波を含む信号として推定する場合には受信信号に含まれる各信号の到来時刻初期値を離散時点上から算出したことを特徴とする請求項２に記載の測距装置。
【請求項４】
信号推定手段は、情報量基準を利用して受信信号に含まれる信号の数を推定したことを特徴とする請求項１に記載の測距装置。
【請求項５】
信号推定手段は、最初に到来した信号を直接波として伝播遅延時間を算出したことを特徴とする請求項４に記載の測距装置。
【請求項６】
信号推定手段は、信号強度が所定の閾値を越えた最初の到来信号を直接波として伝播遅延時間を算出したことを特徴とする請求項４に記載の測距装置。
【請求項７】
信号推定手段は、直接波の信号強度がマルチパス波よりも十分に大きいと想定される場合には、受信信号に含まれる信号の数を所定の数とし、信号強度が一番大きな前記信号を直接波として伝播遅延時間を算出したことを特徴とする請求項１に記載の測距装置。
【請求項８】
信号受信手段は、少なくとも３個の送信機から送信された信号を受信し、
請求項１から７のいずれか１項に記載の測距装置を用いて信号受信手段の位置を測位したことを特徴とする測位装置。
【請求項９】
送信機から送信された信号を受信機で受信する信号受信工程と、
前記信号受信工程で受信した前記信号を推定する信号推定工程と、
前記信号推定工程によって推定した前記信号から前記信号の伝播遅延時間を算出する伝播遅延時間算出工程と、
前記伝播遅延時間から前記送信機と前記受信機との測距をする測距工程とを備え、
前記信号推定工程は、信号モデルのパラメータを周波数領域で最尤推定したことを特徴とする測距方法。
【請求項１０】
受信機は少なくとも３個の送信機から送信された信号を受信し、
請求項９に記載の測距方法を用いて前記受信機の位置を測位する測位工程を備えたことを特徴とする測位方法。
【発明の詳細な説明】
【技術分野】
【０００１】
無線局（送信機）が送信した信号の伝搬遅延時間に基づいて受信機との距離を測定する測距装置およびその測距方法に関するものである。また、当該測距装置および測距方法を用いて受信機の位置を測位する測位装置およびその測位方法に関するものである。
【背景技術】
【０００２】
ＧＰＳ（Ｇｌｏｂａｌ  Ｐｏｓｉｔｉｏｎｉｎｇ  Ｓｙｓｔｅｍ）衛星を用いた測位が、幅広い分野で利用されている。都心部などにおいても需要が高まっているが、マルチパスの影響が測位における大きな誤差要因の一つとなっている。これまで、マルチパス誤差を低減するために、ＭＭＴ（Ｍｕｌｔｉｐａｔｈ  Ｍｉｔｉｇａｔｉｏｎ  Ｔｅｃｈｎｏｌｏｇｙ）（特許文献１）やＮａｒｒｏｗ  Ｃｏｒｒｅｌａｔｏｒ（非特許文献１の１２０ページ）が開発されてきた。
ＭＭＴは、直接波１波とマルチパス１波の信号モデルにおいて、最尤推定法により信号モデルのパラメータ推定を行うことにより、従来困難であった遅延距離の短いマルチパスの影響を低減できることで知られている。
【０００３】
【特許文献１】米国特許番号  ６３７０２０７
【非特許文献１】Global Positioning Systems, Inertial Navigation, and Integration Mohinder S. Grewal, Lawrence R. Weill, Angus P. Andrews John Wiley & Sons, Inc., 2001
【発明の開示】
【発明が解決しようとする課題】
【０００４】
ＭＭＴは、受信信号が直接波１波とマルチパス１波とからなる場合において、マルチパスの影響を低減できる。しかしながら、実際の環境において、マルチパス波の数は１波とは限らない。受信信号に複数のマルチパス波が含まれている場合には、マルチパスの影響を除去しきれない。反対にマルチパス波が含まれない場合には、直接波１波を直接波１波とマルチパス１波として推定してしまうため、誤差を生じてしまう。また、時間領域信号を対象に最尤推定を行うことは計算上の困難を伴う。そこで、この発明は少ない計算量で測距、測位をすることを目的としたものである。
【課題を解決するための手段】
【０００５】
この発明の測距装置は、送信機から送信された信号を受信する信号受信手段と、信号受信手段で受信した信号を推定する信号推定手段と、信号推定手段によって推定した信号から信号の伝播遅延時間を算出する伝播遅延時間算出手段と、伝播遅延時間から送信機と信号受信手段との測距をする測距手段とを備え、信号推定手段は、信号モデルのパラメータを周波数領域で最尤推定したことを特徴とするものである。
【０００６】
また、この発明の測位装置は、信号受信手段は、少なくとも３個の送信機から送信された信号を受信し、前記測距装置を用いて信号受信手段の位置を測位したことを特徴とするものである。
【０００７】
さらに、この発明の測距方法は、送信機から送信された信号を受信機で受信する信号受信工程と、信号受信工程で受信した信号を推定する信号推定工程と、信号推定工程によって推定した信号から信号の伝播遅延時間を算出する伝播遅延時間算出工程と、伝播遅延時間から送信機と受信機との測距をする測距工程とを備え、信号推定工程は、信号モデルのパラメータを周波数領域で最尤推定したことを特徴とするものである。
【０００８】
また、この発明の測位方法は、受信機は少なくとも３個の送信機から送信された信号を受信し、前記測距方法を用いて受信機の位置を測位する測位工程を備えたことを特徴とするものである。
【発明の効果】
【０００９】
この発明の測距装置は、送信機から送信された信号を受信する信号受信手段と、信号受信手段で受信した信号を推定する信号推定手段と、信号推定手段によって推定した信号から信号の伝播遅延時間を算出する伝播遅延時間算出手段と、伝播遅延時間から送信機と信号受信手段との測距をする測距手段とを備え、信号推定手段は、信号モデルのパラメータを周波数領域で最尤推定したので、少ない計算量で測距することができる。
【００１０】
また、この発明の測位装置は、信号受信手段は、少なくとも３個の送信機から送信された信号を受信し、前記測距装置を用いて信号受信手段の位置を測位したので、少ない計算量で測位することができる。
【００１１】
さらに、この発明の測距方法は、送信機から送信された信号を受信機で受信する信号受信工程と、信号受信工程で受信した信号を推定する信号推定工程と、信号推定工程によって推定した信号から信号の伝播遅延時間を算出する伝播遅延時間算出工程と、伝播遅延時間から送信機と受信機との測距をする測距工程とを備え、信号推定工程は、信号モデルのパラメータを周波数領域で最尤推定したので、少ない計算量で測距することができる。
【００１２】
また、この発明の測位方法は、受信機は少なくとも３個の送信機から送信された信号を受信し、前記測距方法を用いて受信機の位置を測位する測位工程を備えた、少ない計算量で測位することができる。
【発明を実施するための最良の形態】
【００１３】
実施の形態１．
図１は本発明の実施の形態１における測位および測距を行うフローチャートである。複数の送信機（例えば、人工衛星）から送信された信号を信号受信手段によって受信し（ＳＴ１０）、受信した信号を信号推定手段により推定する。さらに、推定した信号から伝播遅延時間算出手段により信号の伝播遅延時間を算出する（ＳＴ１４）。算出した伝播遅延時間から位置算出手段により受信機の位置を算出する（ＳＴ１５）。
なお、信号推定手段は、初期値算出手段と信号モデルパラメータ推定手段と信号モデル推定手段とからなる。初期値算出手段が信号モデルのパラメータ初期値を算出し（ＳＴ１１）、信号モデルパラメータ推定手段が信号モデルのパラメータを周波数領域で推定する（ＳＴ１２）。信号モデル推定手段は情報量基準を利用して、信号モデルすなわち受信信号に含まれる信号の数を推定する（ＳＴ１２）。信号推定手段は、信号モデル推定手段によって推定された信号モデルにおける、信号モデルパラメータ推定手段によって推定された信号モデルパラメータを伝播遅延時間算出手段に出力する（ＳＴ１３）。また、信号モデル推定として妥当な推定結果が得られるまでＳＴ１１、ＳＴ１２を繰り返す（ＳＴ１３）。
【００１４】
ここでは、複数の送信機から送信される信号の伝播遅延時間を用いて受信機位置を算出する方法についてまとめて述べたが、受信機と１個の送信機との距離を同様の方法で個々に算出することができる。また、測位をする場合には、送信機となる人工衛星は、少なくとも３個必要となる。４個の送信機から送信された信号を受信した場合には、受信機が備える内部時計の時間のオフセットを調整することができ、正確な測位が可能となる。これに対して、３個の送信機から送信された信号を受信した場合には、例えば、受信機側で地表面データを別途保有することで、正確な測位が可能となる。
【００１５】
図２は本発明の実施の形態１における測位装置及び測距装置のブロック図である。実施の形態１では、ＧＰＳ衛星を用いた測位において実施の形態を説明する。本発明の測位装置及び測距装置は、複数のＧＰＳ衛星から送信されるＧＰＳ信号を受信機のアンテナ１で受信する。受信した信号は、ＲＦ(Ｒａｄｉｏ  Ｆｒｅｑｕｅｎｃｙ)モジュール２で中間周波数信号に周波数変換され、Ａ／Ｄ変換機３において所定周期でデジタル信号としてサンプリングされる。サンプリングされた信号は、信号処理部４においてベースバンド信号に変換され、航法データが復調される。ベースバンド信号および航法データはＲＡＭ７に保存される。初期値算出手段と信号モデルパラメータ推定手段と信号モデル推定手段を備えた信号推定手段、伝播遅延時間算出手段、および位置算出手段は、プログラムとしてＲＯＭ６に保存され、ＣＰＵ５で実行される。
【００１６】
サンプリング間隔Ｔでサンプリングされたベースバンド信号モデルを次式で表す。
2008216111.tif 000003
【００１７】
ここで、ｍ（ｔ）は、時間ｔの関数であり、信号帯域幅に合わせて帯域制限されたＣ−Ａコード、Ｐはマルチパスによる到来信号数を表す。また、各到来信号の振幅をαp、初期位相をθp、コード遅延量をτpで表す。eiθｐは各到来信号の搬送波の位相シフトに対応する複素係数を表し、jはサンプリング時点のインデックス(j番目)を表す。さらに、
α＝（α1，・・・，αP）T
θ＝（θ1，・・・，θP）T
τ＝（τ1，・・・，τP）T
とする。
【００１８】
また、ｉは虚数単位である。計算を容易にするため、実際には等価な次式を用いる。
2008216111.tif 000004
すなわち、αpeiθｐ＝ａp＋ｉｂpとなるよう変数変換を行う。以降、
ａ＝（ａ1，・・・，ａP）T
ｂ＝（ｂ1，・・・，ｂP）T
とする。
【００１９】
（２）式を離散フーリエ変換した式は次式となる。
2008216111.tif 000005
ここで、Ｍ（ω）は、ｍ（ｊＴ）を離散フーリエ変換したものである。
【００２０】
受信したベースバンド信号をｒ（ｊ）、その離散フーリエ変換をＲ（ω）とする。
受信したベースバンド信号は、ノイズを含むので、ｒ（ｊ）＝ｑ（ｊ）＋ｎ（ｊ）と仮定する。ここで、ｎ（ｊ）は複素数のホワイトノイズである。
【００２１】
時間領域信号を対象とした最尤推定においては、次式を最小とするａ、ｂ、τを求める。
2008216111.tif 000006
ここで、最尤推定とは、尤度を手持ちの観測データのもとで、あるパラメータ値が得られる確率とみなして (つまり尤度が未知パラメータの関数とみなして)、尤度を最大化するようなパラメータ値を探索する推定方法をいう。ここでは、ｒ（ｊ）−ｑ（ｊ）を複素数のホワイトノイズとして、その生起確率を尤度として、受信信号に対する信号モデルパラメータａ、ｂ、τを推定する。
【００２２】
しかしながら、（４）式において、コード遅延量τpをｋＴ（ｋは整数）ではない値として求めることは、ｍ（ｊＴ−τp）の算出に計算量を要する。また、τpをｋＴとして（４）式の最小化を行う場合、計算誤差が生じるとともに、各到来信号のτpの組み合わせを探索しようとすると、到来信号数Ｐが大きくなると計算量の爆発を招く。τpをｋＴに丸める処理を入れて非線形最小化手法を用いようとしても、計算の不安定性を招く可能性がある。
【００２３】
そこで、本発明における信号モデルパラメータ推定手段では、上記問題を解決するため、受信信号に対する信号モデルパラメータの推定を周波数領域での最尤推定により行う。すなわち、次式を最小化する。
2008216111.tif 000007
ここで、標本数をＮとするとＮΛ＝Λ´となる。
（５）式を展開すると、次式となる。
2008216111.tif 000008
ここで、Ｒｅ［・］は・の実部を、Ｍ*（ω）はＭ（ω）の共役複素数を表す。
【００２４】
（６）式を最小化するため、（７）、（８）、（９）式を充たすａ、ｂ、τを求める。
2008216111.tif 000009
なお、（７）式〜（９）式において、ｋは１〜Ｐまでの値をとる。
【００２５】
（７）、（８）、（９）式から、それぞれ（１０）、（１１）、（１２）式が導かれる。
2008216111.tif 000010
ここで、Ｉｍ［・］は・の虚部を表す。
【００２６】
（１０）式および（１１）式は、ａk、ｂkに関して線形式となっているため、τの値が決まれば、連立一次方程式を解くことによりａ、ｂを算出することができる。そのため、本発明では、図３のフローチャートのように、信号モデルパラメータを最尤推定する。
まず、τk（ｋ＝１，・・・，Ｐ）の初期値を設定する（ＳＴ２０）。次に、連立一次方程式を解くことによりａk、ｂk（ｋ＝１，…，Ｐ）を算出する（ＳＴ２１）。τkの更新を行う（ＳＴ２２）。更新されたτkが収束したかの収束判定を行う（ＳＴ２３）。τkが収束するまでＳＴ２１とＳＴ２２を繰り返す。τkが収束した時点でのａk、ｂkを算出する（ＳＴ２４）。
【００２７】
τの更新には、ａ、ｂを定数として、Ｎｅｗｔｏｎ法と同様の方法を用いる。
具体的には、まず（１２）式を、ｆk(τ)と置く。ｆk(τ)をτk、τlで偏微分するとそれぞれ次式となる（ｋ、ｌは１〜Ｐまでの値）。
2008216111.tif 000011
【００２８】
ｆ(τ)＝（ｆ1(τ)，・・・，ｆP(τ)）Tとすると、そのＪａｃｏｂｉ行列は次式となる。
2008216111.tif 000012
【００２９】
τの更新値を、τ(new)とすると、τ(new)は次式で算出される（ＳＴ２２）。
2008216111.tif 000013
すべてのτkの変化量が所定の閾値以下となった場合に、τが収束したと判定することができる（ＳＴ２３）。
【００３０】
次に、初期値算出手段におけるτの初期値の算出方法を説明する。Ｐ＝１の場合は、通常のＧＰＳ受信機において利用されている相関器によって算出された、コード遅延量を利用することができる。Ｐ＞１の場合は、Ｐ−１の際に算出したコード遅延量τ1，・・・，τP-1を利用する。まず、τ1，・・・，τP-1をサンプリング時刻ｊＴ上に丸める。複数のτpが同一のサンプリング時刻に丸められた場合は、次のサンプリング時刻に割り当てるようにする。よって、すべてのτpがτp＝τp+1とならないようにする。次に、τ1，・・・，τP-1を定数とし、τ1，・・・，τP-1と異なるサンプリング時刻ｊＴ上から、（６）式の値が最小となるτPを探索する。こうすることにより、ｒ（ｊ）とｍ（ｊＴ）相関関数およびｍ（ｊＴ）の自己相関関数を利用して高速に（６）式を評価することができるため、コード遅延量の初期値を高速に算出できる。また、局所解に陥ることなく、信号モデルパラメータの最尤推定を高速に実行することができる。初期値の算出方法においては、受信したベースバンド信号ｒ（ｊ）とＣ−Ａコードｍ（ｊＴ）とを高分解能化して利用しても良い。
【００３１】
初期値算出手段と信号モデルパラメータ推定手段を以上のように構成することにより、受信信号に複数の信号が含まれている場合に、到来信号数をn波として推定する場合には、（ｎ−１）波の場合での推定結果を利用することができる。また、到来信号数をn波として推定する場合に、（ｎ−１）波の場合での推定結果を利用する際、各信号の到来時刻初期値を離散時点上から算出することができる。
【００３２】
信号モデル推定手段を説明する。本発明の信号モデル推定手段は、情報量基準を利用して、信号モデルすなわち信号の数を推定する。なお、情報量基準とは、モデルの将来の値の分布を予測するための基準であり、真の分布の標本分布に関するエントロピーを最大にする（最大の情報量を得る）ために、モデルのパラメータ白由度を決定する手法である。
【００３３】
ここでは、ＢＩＣ（ベイズ情報量規準）を用いて説明する。ＢＩＣにおいては、次式を最小とするモデルを良いモデルとする。
2008216111.tif 000014
ここで、Θは最大尤度、Ｎは標本数、ｓは独立変数の数を表す。
【００３４】
（１）式の信号の数がＰのモデルにおいては、（１７）式は次式となる。
2008216111.tif 000015
ここで、σは、最小化した（５）式の値をＮの二乗で除算し、平方根をとることにより算出される残差の標準偏差である。（１８）式の右辺第一項は、σを複素数のホワイトノイズに関する標準偏差として用いて、複素数のホワイトノイズの生起確率を算出し、その値を最大尤度Θとすることにより算出される。（１８）式の右辺第二項は、（１）式の信号モデルにおいて、信号１波に対して、独立変数がαk，θk，τkの３つ含まれることから導出される。
信号モデル推定手段は、初期値算出手段と信号モデルパラメータ推定手段をＰ＝１から順に繰り返し、ＢＩＣ（Ｐ）＜ＢＩＣ（Ｐ＋１）となるＰを受信したベースバンド信号ｒ（ｊ）に含まれる信号の数（信号モデル）とする。
【００３５】
情報量基準を利用して受信信号に含まれる信号の数を推定する信号モデル推定手段を備えることで、適切な信号数の信号モデルで、信号モデルパラメータを推定することによって、信号モデルパラメータを正確に推定することができる。
【００３６】
伝播遅延時間算出手段は、信号推定手段によって推定された信号モデルにおける、信号モデルパラメータ推定手段によって推定された信号モデルパラメータとＲＡＭ７に保存された航法データを使って直接波の伝播遅延時間を算出する。ここでは、推定した信号モデルにおいて、推定された信号モデルパラメータから、最初に到来した信号か、あるいは信号強度が所定の閾値を越えた最初の到来信号を直接波として伝播遅延時間を算出する。
位置算出手段は、複数のＧＰＳ衛星からの信号の伝播遅延時間と、ＲＡＭ７に保存された航法データを用いて、通常のＧＰＳ受信機と同様に受信機位置を算出する。
【００３７】
このように、送信機から送信された信号を受信する信号受信手段と、信号受信手段で受信した信号を推定する信号推定手段と、信号推定手段によって推定した信号から信号の伝播遅延時間を算出する伝播遅延時間算出手段と、伝播遅延時間算出手段によって算出した伝播遅延時間から受信機の位置を算出する位置算出手段とを備え、信号推定手段は信号モデルのパラメータ初期値を算出する初期値算出手段と、初期値算出手段によって算出された初期値を利用して、信号モデルのパラメータを周波数領域で推定する信号モデルパラメータ推定手段で、受信信号に対する信号モデルパラメータの最尤推定を周波数領域で行うことによって、複数のマルチパス波を含む受信信号から少ない計算量で信号モデルパラメータを推定することができる。
【００３８】
また、信号モデルパラメータ推定手段において、受信信号に対する信号モデルパラメータの最尤推定を１波少ないマルチパス波を含む信号として推定した結果を利用して周波数領域で行う。さらに、初期値算出手段においては、１波少ないマルチパス波を含む信号として推定した結果を利用して、受信信号に含まれる各信号の到来時刻初期値を離散時点上から算出することによって、受信信号に対する信号モデルパラメータを少ない計算量で安定して推定することができる。
【００３９】
さらに、信号モデル推定手段において、情報量基準を利用して信号モデル、すなわち受信信号に含まれる信号の数を推定する。また、推定した信号モデルにおいて、推定された信号モデルパラメータから最初に到来した信号を直接波として、直接波の伝播遅延時間を算出することができる。よって、適切な信号数の信号モデルで推定した信号モデルパラメータを利用することにより、伝播遅延時間を正確に算出することができる。
【００４０】
また、推定した信号モデルにおいて、推定された信号モデルパラメータから信号強度が所定の閾値を越えた最初の到来信号を直接波として伝播遅延時間を算出する。これによって、直接波よりも早く到来した信号として誤って推定された信号が含まれている場合においても、直接波の伝播遅延時間を正確に算出することができる。
【００４１】
本発明の実施の形態１における直接波の伝播遅延時間推定のシミュレーション結果を図４に示す。信号には、信号強度−１２９ｄＢｍの直接波１波と信号強度−１３５ｄＢｍのマルチパス波が１波含まれている。信号帯域幅は４.０９２ＭＨｚ、直接波とマルチパス波の相対位相差は０°としている。図の横軸は直接波に対するマルチパス波の相対遅延を表し、縦軸は推定誤差を２乗平均平方根誤差であるＲＭＳＥ（Ｒｏｏｔ  Ｍｅａｎ  Ｓｑｕａｒｅ  Ｅｒｒｏｒ）で表している。本発明によるミュレーション結果は、四角のマーカーで示されている。図には、不偏推定量の分散の下限であるＣｒａｍｅｒ−Ｒａｏの下限から算出した、推定誤差の下限も図示している。２波の信号として推定した場合における推定誤差の下限を実線で、本来２波の信号を１波として推定した場合における推定誤差の下限を点破線で示している。本発明の結果は、ほぼ推定誤差の下限を達成していることが分かる。
【００４２】
次に、本発明の実施の形態１における、実測データに対する第１到来波と第２到来波のコード遅延推定量のグラフを図５に示す。また、第１到来波と第２到来波の位相差推定量のグラフを図６に示す。図５および図６の横軸は時間、図５の縦軸はコード遅延量、図６の縦軸は位相差である。図５では、第１到来波のコード遅延推定量における１次成分は消去し、第１到来波の平均コード遅延量を０としている。図５から、およそ５０ｍ遅延したマルチパス波を識別できていることが分かるが、これは計測条件から妥当な値となっている。また、図６から、衛星の移動に伴い、遅延量が次第に変化していることも確認できる。
【００４３】
実施の形態２．
実施の形態２は、実施の形態１に比べ、信号モデル推定手段と伝播遅延時間算出手段とが異なる。
直接波の信号強度がマルチパス波よりも十分に大きいと想定される場合には、信号モデル推定手段が所定の信号数の信号モデルを選択し、伝播遅延時間算出手段が、その信号モデルにおける、信号モデルパラメータ推定手段が推定した信号の中で、信号強度が一番大きな信号を直接波として、伝播遅延時間を算出する。これによって、受信機位置を更に正確に算出することができる。実施の形態２は、信号帯域幅が狭く、推定できる信号の数が少ない場合においても好適である。特に、直接波の信号強度がマルチパス波よりも十分大きい場合に、直接波に対する相対遅延時間が小さいマルチパス波が含まれている場合に有効である。
【００４４】
本発明の実施の形態２における直接波の伝播遅延時間推定のシミュレーション結果を図７に示す。信号には、信号強度−１２９ｄＢｍの直接波１波と信号強度−１３５ｄＢｍのマルチパス波が１波含まれている。信号帯域幅２.０４６ＭＨｚ、直接波とマルチパス波の相対位相差は０°としている。図の横軸は直接波に対するマルチパス波の相対遅延を表し、縦軸は推定誤差をＲＭＳＥで表している。実施の形態１によるミュレーション結果を四角のマーカーで、実施の形態２によるミュレーション結果を点のマーカーで示している。図には、不偏推定量の分散の下限であるＣｒａｍｅｒ−Ｒａｏの下限から算出した、推定誤差の下限も図示している。２波の信号として推定した場合における推定誤差の下限を実線で、本来２波の信号を１波として推定した場合における推定誤差の下限を点破線で示している。実施の形態２の結果は、実施の形態１の結果よりも推定誤差が小さく優れていることが分かる。なお、不偏推定量の分散の下限であるＣｒａｍｅｒ−Ｒａｏの下限から算出した推定誤差の下限よりも、実施の形態２の結果が優れているのは、実施の形態２の推定量が不偏推定量ではない（推定量に偏りがある）からである。
【図面の簡単な説明】
【００４５】
【図１】本発明の実施の形態１における測位および測距を行うフローチャートである。
【図２】本発明の実施の形態１における測位装置および測距装置のブロック図である。
【図３】本発明の実施の形態１における最尤推定のフローチャートである。
【図４】本発明の実施の形態１における伝播遅延時間推定のシミュレーション結果を表す図である。
【図５】本発明の実施の形態１における実測データに対する第１到来波と第２到来波のコード遅延推定量の図である。
【図６】本発明の実施の形態１における実測データに対する第１到来波と第２到来波の位相差推定量の図である。
【図７】本発明の実施の形態２における伝播遅延時間推定のシミュレーション結果を表す図である。
【符号の説明】
【００４６】
１  アンテナ、２  ＲＦモジュール、３  Ａ／Ｄ変換機、４  信号処理部、５  ＣＰＵ、６  ＲＯＭ、７  ＲＡＭ。
【図１】
2008216111.tif 000016
【図２】
2008216111.tif 000017
【図３】
2008216111.tif 000018
【図４】
2008216111.tif 000019
【図５】
2008216111.tif 000020
【図６】
2008216111.tif 000021
【図７】
2008216111.tif 000022