(57)【要約】      （修正有）【課題】高精度かつ小さい計算負荷で到来波の方向と強度を推定する。【解決手段】アレーアンテナの第Ｋ（Ｋ＝１〜Ｍ）素子の信号をそれ以外の素子の信号によって予測することにより複数のウエイトベクトルを算出し、それぞれのウエイトベクトルから複数の角度スペクトラムを計算する。複数のウエイトベクトルから算出される角度スペクトラムを平均処理によって合成することにより、真の到来方向でないものはそのスペクトラム強度が下がり、真の到来方向にあたるスペクトラム強度のみが高いままに残り、真の到来方向を推定することが可能となる。【選択図】          図４
2004198218.tif 000002
【特許請求の範囲】
【請求項１】
到来波から得られる入力信号から相関行列を算出する相関行列算出手段と、
前記相関行列の逆行列を計算する逆行列計算手段と、
前記相関行列の逆行列から複数のウエイトベクトルを計算するウエイトベクトル算出手段と、
前記複数のウエイトベクトルを用いて、推定される到来波の方向と信号強度の分布を算出する到来波分布算出手段と、
前記推定される到来波の信号強度に基づいて到来波の方向を判定する到来波判定手段と
を有することを特徴とする到来波推定装置。
【請求項２】
到来波を受信する複数のアンテナから構成されるアレーアンテナと、
前記アレーアンテナで受信した信号を前記入力信号へ変換する受信手段と
をさらに有することを特徴とする請求項１に記載の到来波推定装置。
【請求項３】
前記到来波分布算出手段は、
前記複数のウエイトベクトルの各々から角度スペクトラムを算出する角度スペクトラム算出手段と、
前記算出した複数の角度スペクトラムを合成して信号強度を推定する信号強度推定手段と
を有することを特徴とする請求項１又は請求項２に記載の到来波推定装置。
【請求項４】
前記到来波分布算出手段は、
前記複数のウエイトベクトルのいずれか１つを用いて角度スペクトラムを算出する角度スペクトラム算出手段と、
前記算出した角度スペクトラムの極大値を検出して到来波候補とする到来波候補検出手段と、
前記複数のウエイトベクトルを用いて前記到来波候補の信号強度を推定する信号強度推定手段と
を有することを特徴とする請求項１又は請求項２に記載の到来波推定装置。
【請求項５】
空間平均演算手段をさらに有する
ことを特徴とする請求項１乃至請求項４のいずれか１項に記載の到来波推定装置。
【請求項６】
到来波から得られる入力信号から相関行列を算出する相関行列算出工程と、
前記相関行列の逆行列を計算する逆行列計算工程と、
前記相関行列の逆行列から複数のウエイトベクトルを計算するウエイトベクトル算出工程と、
前記複数のウエイトベクトルを用いて、推定される到来波の方向と信号強度の分布を算出する到来波分布算出工程と、
前記推定される到来波の信号強度に基づいて到来波の方向を判定する到来波判定工程と
を有することを特徴とする到来波推定方法。
【請求項７】
前記到来波分布算出工程は、
前記複数のウエイトベクトルの各々から角度スペクトラムを算出する角度スペクトラム算出工程と、
前記算出した複数の角度スペクトラムを合成して信号強度を推定する信号強度推定工程と
を含むことを特徴とする請求項６に記載の到来波推定方法。
【請求項８】
前記到来波分布算出工程は、
前記複数のウエイトベクトルのいずれか１つを用いて角度スペクトラムを算出する角度スペクトラム算出工程と、
前記算出した角度スペクトラムの極大値を検出して到来波候補とする到来波候補検出工程と、
前記複数のウエイトベクトルを用いて前記到来波候補の信号強度を推定する信号強度推定工程と
を含むことを特徴とする請求項６に記載の到来波推定方法。
【請求項９】
波源から得られる入力信号から相関行列を算出する相関行列算出手段と、
前記相関行列の逆行列を計算する逆行列計算手段と、
前記相関行列の逆行列から複数のウエイトベクトルを計算するウエイトベクトル算出手段と、
前記複数のウエイトベクトルを用いて、推定される波源の空間分布または時間分布を算出する信号強度分布算出手段と、
前記推定される波源の信号強度分布に基づいて波源の空間的位置または時間的位置を推定する波源位置判定手段と
を有することを特徴とする波源位置推定装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、電波の到来波を推定する到来波推定装置に関し、複数のアンテナで構成されるアレーアンテナで受信された信号に基づいて、その到来波の方向と強度を推定する到来波推定装置及び到来波推定方法に関するものである。また、電波の波源の空間的位置または時間的位置を推定する波源位置推定装置に関する。
【０００２】
【従来の技術】
従来、複数のアンテナで構成されるアレーアンテナを用いて到来波分布を把握するために多くの到来方向推定手法が提案されている。古典的な手法としてはビームフォーマ法(beamformer)、Ｃａｐｏｎ法、線形予測法などがある。また、高分解能到来方向推定手法としては、ＭＵＳＩＣ(Multiple Signal Classification)やＥＳＰＲＩＴ(Estimation of Signal Parameters via Rotational Invariance Techniques)などがある。これらは非特許文献１に詳しく記載されている。
【０００３】
中でもＭＵＳＩＣやＥＳＰＲＩＴは高分解能かつ高精度に到来方向を推定できることから広く利用されている。しかし、ＭＵＳＩＣやＥＳＰＲＩＴでは、その信号処理過程において、固有値・固有ベクトルを計算する必要があり、計算負荷が非常に大きいため処理を高速化することが難しい。これに対して、線形予測法は固有値計算が必要ないために、小さい計算負荷で到来波の方向を検出することが可能である。しかし、ビームフォーマ法やＣａｐｏｎ法とは異なり、算出される角度スペクトラムにおいてスペクトラム強度は到来波の受信強度に比例しない。そのため、検出した方向が真の到来方向であるかどうかを判定するためには、到来波の受信信号強度を別手段によって推定しなければならない。
【０００４】
特許文献１に記載の電波到来方向推定装置によれば、線形予測法により到来方向候補を推定し、それら候補の信号強度を推定する部分には電力推定精度の高いＣａｐｏｎ法を用いて、到来方向候補のアンテナ出力電力を算出している。
【０００５】
【特許文献１】
特開２００２−１４８３２４号公報
【非特許文献１】
菊間信良著、「アレーアンテナによる適応信号処理」科学技術出版、１９９８年１１月、ｐ．１７３−２４５
【０００６】
【発明が解決しようとする課題】
しかしながら、Ｃａｐｏｎ法による電力推定演算も必ずしも計算負荷が小さいわけではない。Ｃａｐｏｎ法による角度θにおける電力Ｐcp（θ）を求める演算は式（１）で表される。
【数１】
2004198218.tif 000003
ここで、ａはモードベクトル、Ｒxxは到来波信号の相関行列である。ａ、Ｒxxは複素ベクトル及び複素行列であるが、複素数演算は計算負荷が高い。そのため高速な処理を実現させるためには複素数演算を削減する必要がある。
【０００７】
また、線形予測法とＣａｐｏｎ法では精度や分解能に差がある。そのため真の到来方向であったとしても、両者の角度スペクトラムにおいてはスペクトラム強度のピーク位置がずれることがある。特に近接した方向から複数の到来波が受信される場合には、線形予測法では分解可能でもＣａｐｏｎ法では分解不可能であるような状況がある。このような場合には、特許文献１に記載の線形予測法とＣａｐｏｎ法を組み合わせた手法によっても、到来方向候補の信号強度を正しく推定することができなくなる。
【０００８】
本発明は上記のような課題を解決するものであり、複数のアンテナで構成されるアレーアンテナで受信した信号から、高精度かつ小さい計算負荷で到来波の方向（つまり波源の空間的位置）を推定することが可能な到来波推定装置及び到来波推定方法、並びに波源位置推定装置を提供することを目的とする。
【０００９】
また、移動通信や無線ＬＡＮにおいてはフェージングによって通信品質が劣化することがある。多重波の伝播遅延時間はこのフェージングを特徴付ける重要なパラメータであり、この伝播遅延時間を推定することが重要である（非特許文献、Ｐ２６９−２７４）。本発明は、波源から得られる信号を基に、高精度かつ小さい計算負荷で多重波の伝播遅延時間（つまり波源の時間的位置）を推定することが可能な波源位置推定装置を提供することを目的とする。
【００１０】
【課題を解決するための手段】
本願の請求項１に記載の到来波推定装置は、到来波から得られる入力信号から相関行列を算出する相関行列算出手段と、前記相関行列の逆行列を計算する逆行列計算手段と、前記相関行列の逆行列から複数のウエイトベクトルを計算するウエイトベクトル算出手段と、前記複数のウエイトベクトルを用いて、推定される到来波の方向と信号強度の分布を算出する到来波分布算出手段と、前記推定波源の信号強度に基づいて到来波の方向を判定する到来波判定手段とを有することを特徴とする。
【００１１】
さらに本願の請求項２に記載の到来波推定装置は、到来波を受信する複数のアンテナから構成されるアレーアンテナと、前記アレーアンテナで受信した信号を前記入力信号へ変換する受信手段とをさらに有することを特徴とする。
【００１２】
また本願の請求項３に記載の到来波推定装置は、到来波分布算出手段が、前記複数のウエイトベクトルの各々から角度スペクトラムを算出する角度スペクトラム算出手段と、前記算出した角度スペクトラムを合成して信号強度を推定する信号強度推定手段とを有することを特徴とする。
【００１３】
また本願の請求項４に記載の到来波推定装置は、到来波分布算出手段が、前記複数のウエイトベクトルのいずれか１つを用いて角度スペクトラムを算出する角度スペクトラム算出手段と、前記角度スペクトラムの極大値を検出して到来波候補とする到来波候補検出手段と、前記複数のウエイトベクトルを用いて、前記到来波候補の信号強度を推定する信号強度推定手段とを有することを特徴とする。
【００１４】
本願の請求項５に記載の到来波推定装置は、空間平均演算手段をさらに有することを特徴とする。
【００１５】
また、本願の請求項６に記載の到来波推定方法は、到来波から得られる入力信号から相関行列を算出する相関行列算出工程と、前記相関行列の逆行列を計算する逆行列計算工程と、前記相関行列の逆行列から複数のウエイトベクトルを計算するウエイトベクトル算出工程と、前記複数のウエイトベクトルを用いて推定される到来波の方向と信号強度の分布を算出する到来波分布算出工程と、前記推定波源の信号強度に基づいて到来波の方向を判定する到来波判定工程とを有することを特徴とする。
【００１６】
また、本願の請求項７に記載の到来波推定方法は、前記到来波分布算出工程に、前記複数のウエイトベクトルの各々から角度スペクトラムを算出する角度スペクトラム算出工程と、前記算出した角度スペクトラムを合成して信号強度を推定する信号強度推定工程とを含むことを特徴とする。
【００１７】
また、本願の請求項８に記載の到来波推定方法は、前記到来波分布算出工程に、前記複数のウエイトベクトルのいずれか１つを用いて角度スペクトラムを算出する角度スペクトラム算出工程と、前記算出した角度スペクトラムの極大値を検出して到来波候補とする到来波候補検出工程と、前記複数のウエイトベクトルを用いて前記到来波候補の信号強度を推定する信号強度推定工程とを含むことを特徴とする。
【００１８】
また、本願の請求項９に記載の波源位置推定装置は、波源から得られる入力信号から相関行列を算出する相関行列算出手段と、前記相関行列の逆行列を計算する逆行列計算手段と、前記相関行列の逆行列から複数のウエイトベクトルを計算するウエイトベクトル算出手段と、前記複数のウエイトベクトルを用いて、推定される波源の空間分布または時間分布を算出する信号強度分布算出手段と、前記推定される波源の信号強度分布に基づいて波源の空間的位置または時間的位置を推定する波源位置判定手段とを有することを特徴とする。
【００１９】
また、本願の発明は、多重波の伝播遅延時間を推定する遅延時間推定装置であって、多重波の入力信号から相関行列を算出する相関行列算出手段と、前記算出した相関行列の逆行列を計算する逆行列算出手段と、前記計算した逆行列から複数のウエイトベクトルを計算するウエイトベクトル計算手段と、前記複数のウエイトベクトルを用いて、多重波の伝播遅延時間に対する信号強度の分布を算出する遅延時間分布算出手段と、前記信号強度に基づいて多重波の伝播遅延時間を算出する遅延時間算出手段とを有することを特徴とする。
【００２０】
【発明の作用及び効果】
本発明の到来波推定装置は、到来波から得られる入力信号から相関行列を算出し、算出した相関行列の逆行列を計算する。計算された相関行列の逆行列から複数のウエイトベクトルを計算し、その複数のウエイトベクトルを用いて、推定される到来波の方向と信号強度の分布を算出し、その結果、推定波源の信号強度に基づいて到来波の数とその方向とを判定するものである。
【００２１】
一般に固有値計算を必要としない線形予測法はＭＵＳＩＣやＥＳＰＩＲＩＴに比べて計算負荷が小さいことはよく知られている。従来の線形予測法では相関行列の第１列から導かれるウエイトベクトルにより角度スペクトラムを算出し、スペクトラム強度のピークを検出することで到来方向を推定する。しかし、真でない方向、つまり実際には到来しない方向にもピークが検出されることがある。一方、ウエイトベクトルは相関行列の第１列に限らず、どの列成分からでも導出可能であり、相関行列の次数分だけ算出することができる。当然それぞれのウエイトベクトルに対応した角度スペクトラムが算出可能であるが、真でない方向に突出するピークは使用するウエイトベクトルによって出現したり消失したりする。したがって、複数のウエイトベクトルを使用することにより、真でない到来方向のピークを判断することが可能である。本発明はこのことを利用することで、高精度かつ小さい計算負荷で到来波の方向と強度を推定することができる。
【００２２】
以上のように本発明は線形予測法を基にしており、ＭＵＳＩＣやＥＳＰＲＩＴで行われる固有値計算を必要としないため高速演算が可能である。また、従来の線形予測法では直接推定することができない受信信号強度を少ない演算量で推定できるため、真の到来方向を高速に検出することができる。
【００２３】
ここで、従来の線形予測法と本発明との違いを説明する。
線形予測法では式（２）によってウエイトベクトルを計算し、それを用いて角度スペクトラム（式（３））が計算される。このウエイトベクトルはアレイアンテナの第１素子の信号を第２素子から第Ｍ素子の信号を用いて予測し、その予測誤差を最小化することによって導かれたものである。図１に線形予測法による角度スペクトラムの例を示す。図１において真の到来方向は０°と１０°であるが、それ以外の方向にもピークが存在することが見て取れる。
【数２】
2004198218.tif 000004
【００２４】
本発明では線形予測法の考え方を拡張して、第Ｋ（Ｋ＝１〜Ｍ）素子の信号をそれ以外の素子の信号によって予測することにより複数のウエイトベクトルを算出する。したがって、それぞれのウエイトベクトルから角度スペクトラムを計算することができる。図２に角度スペクトラムを示す。これらの角度スペクトラムには次のような性質がある。
（ａ）いずれの角度スペクトラムにおいても真の到来方向はピークとして現れる。
（ｂ）真の到来方向でないピークはその強度や位置が角度スペクトラムによって様々に変化する。
したがって、複数のウエイトベクトルから算出される角度スペクトラムを例えば調和平均によって合成すると（図３）、真の到来方向でないものはそのピークが下がり、真の到来方向にあたるピークのみが高いままに残る。また、図３のような合成後の角度スペクトラムには相関行列の持つ情報がほぼすべて含まれているため、従来の線形予測法による角度スペクトラムに比べると、ピークの高さが到来波の電力に対応したものとなっている。
【００２５】
また、本発明では算出した複数のウエイトベクトルの全てを使わずとも、例えば１つのウエイトベクトルを用いて角度スペクトラムを求めるようにしてもよい。この場合、求めた角度スペクトラムの極大値を到来波の候補として、候補となった到来波方向の信号強度について複数のウエイトベクトルを用いて求めるようにする。このため、複数の角度スペクトラムを算出することなく、信号強度を到来波候補についてのみ計算するだけでよく、さらに計算負荷を低減することができる。
【００２６】
また、算出した相関行列に対して空間平均演算（例えば、Forward-Backward空間平均演算等）を行うことによって、算出すべきウエイトベクトルの個数をほぼ半分にすることができるため、さらに計算負荷を低減することができる。
【００２７】
本発明は上記のように、複数のウエイトベクトルを用いることで、高精度かつ小さい計算負荷で真の到来方向を判断することができるという効果が得られるのである。
【００２８】
また本発明の波源位置推定装置は、波源から得られる多重波の信号、例えば周波数領域のアレーデータから相関行列を算出し、算出した相関行列の逆行列を計算する。計算された相関行列の逆行列から複数のウエイトベクトルを計算し、その複数のウエイトベクトルを用いて、推定される波源から得られる多重波の信号の遅延時間と信号強度の分布（すなわち時間分布）を算出し、その結果、推定波源の信号強度に基づいて多重波の伝播遅延時間（時間的位置）を推定するものである。
【００２９】
本発明の到来波推定装置は、複数のアンテナから構成されるアレーアンテナで受信した信号を入力信号とし、角度スペクトラムを算出しているのに対して、本発明の波源位置推定装置は、一つのアンテナで受信した多重波を入力信号とし、遅延時間スペクトラムを算出している。本発明の波源位置推定装置は、本発明の到来波推定装置と同様にして、複数のウエイトベクトルを用いることで、高精度かつ小さい計算負荷で多重波の伝播遅延時間（時間的位置）を推定することができるという効果が得られるのである。
【００３０】
【発明の実施の形態】
本発明の到来波推定装置の実施の形態を説明する。
（第１の実施の形態）
図４は本発明の第１の実施の形態の構成要素を示すものである。本発明の到来波推定装置は、到来波を受信するＭ個のアンテナから構成されるアレーアンテナ１０、アンテナからの信号が入力される受信機２０、受信信号から相関行列を算出する相関行列算出器３０、相関行列の逆行列を算出する逆行列算出器４０、ウエイトベクトルを算出するウエイトベクトル算出器５０、ウエイトベクトルを用いて角度スペクトラムを算出する角度スペクトラム算出器６０、複数の角度スペクトラムを合成して到来波の信号強度分布を推定する信号強度推定器７０、推定された信号強度によって到来波の方向を判定する到来波判定器８０からなる。
【００３１】
次に到来波検出の流れを説明する。まず、各アレーアンテナで受信された信号ｘ1〜ｘMを元に、相関行列算出器３０において式（４）により相関行列Ｒxxが算出される。ただし、Ｌはスナップショット数である。さらに逆行列算出器４０により相関行列の逆行列が算出される。
【数３】
2004198218.tif 000005
【００３２】
次にウエイトベクトル算出器５０によって、相関行列の逆行列を基にウエイトベクトルが算出される。本実施の形態ではＮ個のウエイトベクトル算出器が用意され、式（５）に基づいてｉ番目のウエイトベクトル算出器がウエイトベクトルＷiを算出する。ここでＮは相関行列Ｒxxの次数である。またＴiはｉ番目の成分だけ１で他の成分は０であるベクトルを表す。なお、本実施の形態では複数のウエイトベクトル算出器を用いているが、必ずしも複数である必要はなく１つのウエイトベクトル算出器であってもよい。複数のウエイトベクトルが算出できればよい。
【数４】
2004198218.tif 000006
【００３３】
角度スペクトラム算出器６０ではそれぞれウエイトベクトルＷiを用いて式（６）により角度スペクトラムＰi（θ）が算出される。信号強度推定器７０では式（７）に基づき、それぞれの角度スペクトラムＰi（θ）が合成される。式（７）は複数の角度スペクトラムを調和平均するものである。これにより、Ｐi（θ）のいずれか一つで突発的に高い値をとるものが存在しても平均値Ｐ（θ）にはあまり影響を及ぼさず、すべてのＰi（θ）が高いときに平均値Ｐ（θ）も高くなる。したがって、平均値Ｐ（θ）が高くなるのはすべてのＰi（θ）においてピークをとる真の到来方向のみである。
【数５】
2004198218.tif 000007
【００３４】
到来波判定器８０では、信号強度推定器７０で算出された平均スペクトラムから信号強度の強いピークを検出し、到来波の方向と信号強度を算出する。
【００３５】
（第２の実施の形態）
図５に本発明の第２の実施の形態の到来波推定装置の構成を示す。この構成において角度スペクトラム算出器６１は一つであり、新たに到来波候補検出器６２が導入される。
【００３６】
複数のウエイトベクトルを算出するところまでは第１の実施の形態と同じであるので省略する。角度スペクトラム算出器６１において、算出された複数のウエイトベクトルの一つ（ここではＷ1）を用いて角度スペクトラムＰLP（θ）が算出される。到来波候補検出器６２は算出された角度スペクトラムの極大値を検出して到来波候補θjとする。信号強度推定器７０では検出された到来波候補θjについてのみ、式（７）を使用して信号強度が計算される。到来波判定器８０では推定された信号強度を基に到来波候補の中から真の到来波が判定される。
【００３７】
本実施の形態では、Ｎ個すべての角度スペクトラムを算出する必要がなく、信号強度も到来波候補についてのみ計算するだけでよいので、第１の実施の形態に比べ演算量を削減することができる。
【００３８】
尚、本実施の形態では、角度スペクトラムを算出するために、ウエイトベクトルＷ1を用いているが、他のウエイトベクトルであっても良い。また、一つのウエイトベクトルだけでなく、到来波の候補を検出できればそれ以上の複数のウエイトベクトルを用いるようにしても構わない。
【００３９】
（第３の実施の形態）
第３の実施の形態の到来波推定装置の構成を図６に示す。上記実施の形態との違いは、相関行列演算器３０と逆行列演算器４０との間に空間平均演算器９０が加えられることである。
【００４０】
この空間平均演算器９０では相関行列にForward-Backward空間平均法が施される。Forward-Backward空間平均法は、入力ベクトル（Forwardベクトル：式（８））から計算される相関行列（式（１０））と、入力ベクトルの複素共役を取り成分を逆順に並べたベクトル（Backwardベクトル：式（９））から計算される相関行列（式（１１））を式（１２）に示すように平均する手法である。
【数６】
2004198218.tif 000008
【００４１】
この操作により、信号強度推定器７０における式（６）の演算結果には式（１３）に示すような関係が成り立つ。そのため、Ｎ個すべてのウエイトベクトルを算出する必要はなく、Ｎが偶数のときにＮ／２番目まで、Ｎが奇数のとき（Ｎ＋１）／２番目までのウエイトベクトルを算出すればよいので、ウエイトベクトル算出器の個数をほぼ半減することができる。さらに、信号強度推定器７０における演算も式（１４）のように変更できるため、複素数演算もほぼ半減することができるようになる。
【数７】
2004198218.tif 000009
【００４２】
（第４の実施の形態）
図７に本発明の第４の実施の形態に係る波源位置推定装置の構成を示す。本実施の形態の波源位置推定装置は、波源からの信号を受信するアンテナ１１、アンテナからの信号が入力される受信機２１、受信信号から相関行列を算出する相関行列算出器３１、相関行列の逆行列を算出する逆行列算出器４１、ウエイトベクトルを算出するウエイトベクトル算出器５１、ウエイトベクトルを用いて遅延時間スペクトラムを算出する遅延時間スペクトラム算出器１６０、複数の遅延時間スペクトラムを合成して信号強度分布を算出する信号強度推定器１７０、推定された信号強度によって多重波の伝播遅延時間を推定する遅延時間推定器１８０からなる。
【００４３】
ある波源からは、式（１５）に示すようなΔｆの一定間隔で掃引されたｍ個の周波数の信号が順次送信される。波源から送信された信号は様々な構造物による反射や回折によって、多重波となってアンテナ１１に到来する。アンテナ１１は、波源から送信されたｍ個の周波数の信号を受信する。アンテナ１１で受信されたｍ個の信号ｘ（ｆi）（ｉ＝１，…，ｍ）が、受信機２１から、相関行列算出器３１に渡され、式（１６）により相関行列Ｒxxが算出される。
【数８】
2004198218.tif 000010
【００４４】
以降の手順は第１の実施の形態と同様である。第１の実施の形態では、角度スペクトラム算出器６０において角度スペクトラム（横軸−角度、縦軸−スペクトラム強度）を算出したのに対し、本実施の形態では遅延時間スペクトラム算出器１６０において遅延時間スペクトラム（横軸−遅延時間、縦軸−スペクトラム強度）を算出する。そして、算出した各遅延時間スペクトラムを信号強度推定器１７０において合成（調和平均）する。遅延時間推定器１８０では、算出された平均スペクトラムから信号強度の強いピークを算出し、遅延時間と信号強度を算出する。このように、複数のウエイトベクトルを用いて、遅延スペクトラムを算出することで、高精度かつ小さい計算負荷で多重波の伝送遅延時間を推定することが出来る。
【００４５】
なお、本実施の形態の波源位置推定装置は、多重波の入力信号から相関行列を算出する相関行列算出手段と、前記算出した相関行列の逆行列を計算する逆行列算出手段と、前記計算した逆行列から複数のウエイトベクトルを計算するウエイトベクトル計算手段と、前記複数のウエイトベクトルを用いて、多重波の伝播遅延時間に対する信号強度の分布を算出する遅延時間分布算出手段と、前記信号強度に基づいて多重波の伝播遅延時間を算出する遅延時間算出手段とを有することを特徴とする多重波の伝播遅延時間を推定する遅延時間推定装置とも言い換えることができる。
【００４６】
また、本発明の到来波推定装置及び到来波推定方法、並びに波源位置推定装置は、上述の図示例にのみ限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変更を加え得ることは勿論である。
【図面の簡単な説明】
【図１】従来技術の線形予測法による角度スペクトラムの例を示す図。
【図２】本発明の第１の実施の形態の複数ウエイトベクトルによる角度スペクトラムを示す図。
【図３】本発明の第１の実施の形態の複数ウエイトベクトルによる角度スペクトラムを合成した角度スペクトラムを示す図。
【図４】本発明の第１の実施の形態の構成を表すブロック図。
【図５】本発明の第２の実施の形態の構成を表すブロック図。
【図６】本発明の第３の実施の形態の構成を表すブロック図。
【図７】本発明の第４の実施の形態の構成を表すブロック図。
【図１】
2004198218.tif 000011
【図２】
2004198218.tif 000012
【図３】
2004198218.tif 000013
【図４】
2004198218.tif 000014
【図５】
2004198218.tif 000015
【図６】
2004198218.tif 000016
【図７】
2004198218.tif 000017