WO2006030834A1

WO2006030834A1 - 信号到来方向推定装置、信号到来方向推定方法、および信号到来方向推定用プログラム

Info

Publication number: WO2006030834A1
Application number: PCT/JP2005/016956
Authority: WO
Inventors: Akira Tanaka
Original assignee: Hokkaido University NUC
Current assignee: Hokkaido University NUC
Priority date: 2004-09-14
Filing date: 2005-09-14
Publication date: 2006-03-23
Anticipated expiration: 2007-03-14
Also published as: US7436358B2; JP4660773B2; US20080001821A1; JPWO2006030834A1

Abstract

　入力信号間の相関行列や雑音相関行列が特異である場合でも、信号到来方向を推定することができる信号到来方向推定装置。本装置では、まず、観測手段（１）で、雑音のみが存在する状況と到来方向推定対象の音が入力されている状態とのそれぞれを観測し、短時間フーリエ変換手段（２）で、短時間フーリエ変換を行った後、相関行列算出手段（３）で、雑音相関行列および入力信号相関行列を算出する。得られた相関行列は、雑音相関行列保持手段（４）および入力信号相関行列保持手段（５）にそれぞれ保持される。次に、真固有値・真固有ベクトル・偽固有ベクトル算出手段（６）で、これらの相関行列を用いて、入力信号相関行列に対する雑音相関行列の真固有値と真固有ベクトルと偽固有ベクトルを算出し、補空間成分行列算出手段（７）で、得られた真固有値と真固有ベクトルと偽固有ベクトルとに基づいて、信号部分空間の補空間成分を求める行列を算出した後、到来方向探索手段（８）で、その信号部分空間の補空間成分を求める行列を用いて、音の到来方向を探索する。

Description

明細書

信号到来方向推定装置、信号到来方向推定方法、および信号到来方向推定用プログラム

技術分野

[0001] 本発明は、電波や音等の信号の到来方向を推定する信号到来方向推定装置、信号到来方向推定方法、および信号到来方向推定用プログラムに関する。

背景技術

[0002] 従来の信号到来方向推定方式として、 MUSIC (Multiple Signal Classification)法や ESPRIT (Estimation of signal Parameters via Rotational Invariance Techniques) 法等がよく知られている。

[0003] MUSIC法は、アレーアンテナ等に同時に到来する複数の電波等の信号の方向を測定するアルゴリズムであり、到来信号とは関係ない雑音の成分を用いて到来信号のパラメータを推定する方法である。 MUSIC法による電波到来方向測定では、例えば、受信信号力も得られる相関行列の固有値を求め、得られた固有値が到来波数によって信号固有値と雑音固有値とに分けられるため、雑音固有値に対応する雑音固有ベクトルを求め、得られた雑音固有ベクトル力角度スペクトルを求め、それらを角度（つまり、周波数)で平均して MUSICスペクトルを求めることにより、到来信号の到来方向を推定する。また、 ESPRIT法は、基本原理は MUSIC法と同様である力平行なアレー面を有する複数のアレー対を配置することにより、 MUSIC法において要請されるアレーの幾何学的な情報や計算量を緩和する手法である。

[0004] また、その改良手法として、例えば、特許文献 1、特許文献 2、および特許文献 3に記載の技術が提案されて、る。

[0005] 特許文献 1記載の技術は、固有値展開を用いる代わりに、入力信号相関行列の逆行列を利用することにより、少ない計算量での到来方向推定を実現している。また、特許文献 2記載の技術は、一つ目の受信器とその他の受信器の相関ベクトルのみからテブリッツ (To印 litz)行列を構成し、これを入力信号相関行列の代わりに用いることにより、計算量の削減を実現している。また、特許文献 3記載の技術は、アレーを物理的に回転させることにより、アレー間隔決定時に想定していな力つた周波数の信号源が存在する場合でも到来方向推定を実現している。これらの技術は、いずれも、入力信号相関行列が正則 (逆行列を持つこと)であることを前提として、計算量の削減や性能の向上を実現するものである。

特許文献 1 :特開 2002— 148324号公報

特許文献 2：特開平 11 133130号公報

特許文献 3：特開 2001— 108734号公報

発明の開示

発明が解決しょうとする課題

[0006] し力しながら、特許文献 1〜3に記載された信号到来方向推定方式にあっては、入力信号相関行列が正則であることを前提としているため、例えば、測定空間内にマイク数よりも少ない数の雑音源が存在する場合など、入力信号相関行列や雑音相関行列が特異である場合には、適用することができない、という問題がある。雑音に関する情報を利用せずに、単に入力信号相関行列の固有値展開を用いた場合は、求めた固有値が雑音に由来するもの力信号に由来するものかを判別できない可能性があり、また、入力信号相関行列と雑音相関行列を用いようとした場合は、原理的に雑音相関行列の逆行列が必要となるため、雑音相関行列が特異であった場合には、そもそも当該逆行列が存在しない、つまり、計算すらできないからである。

[0007] 本発明の目的は、入力信号相関行列や雑音相関行列が特異である場合でも、信号到来方向を推定することができる信号到来方向推定装置、信号到来方向推定方法、および信号到来方向推定用プログラムを提供することである。

課題を解決するための手段

[0008] 本発明は、到来方向推定対象の信号が存在しない雑音のみが存在する状況下における雑音相関行列と、到来方向推定対象の信号および雑音が存在する状況下における入力信号相関行列とを算出し、前記雑音相関行列および前記入力信号相関行列を用、て、前記入力信号相関行列に対する前記雑音相関行列の真固有値と真固有ベクトルと偽固有ベクトルとを算出し、前記真固有値と前記真固有ベクトルと前記偽固有ベクトルとに基づいて、前記信号の到来方向を推定する、構成を採る。発明の効果

[0009] 本発明によれば、入力信号相関行列に対する雑音相関行列の真固有値と真固有ベクトルと偽固有ベクトルとを求め、この真固有値と真固有ベクトルと偽固有ベクトルとに基づき信号の到来方向を推定するようにしたので、入力信号相関行列や雑音相関行列が特異な場合でも、信号到来方向を推定することができる。

図面の簡単な説明

[0010] [図 1]本発明の一実施の形態に係る信号到来方向推定装置の構成を示すブロック図発明を実施するための最良の形態

[0011] 以下、本発明の実施の形態について、図面を参照しながら説明する。なお、本実施の形態では、信号到来推定方式として、 MUSIC法を基礎として、その改良手法という形で説明する。

[0012] ここで、 MUSIC法とは、上記のように、アレーアンテナ等に同時に到来する複数の電波等の信号の方向を測定するアルゴリズムであり、到来信号とは関係ない雑音の成分を用いて、到来信号のパラメータを推定する方法である。 MUSIC法による電波到来方向測定では、上記のように、例えば、受信信号から得られる相関行列の固有値を求め、得られた固有値が到来波数によって信号固有値と雑音固有値とに分けられるため、雑音固有値に対応する雑音固有ベクトルを求め、得られた雑音固有べタトルカ角度スペクトルを求め、それらを角度（つまり、周波数)で平均して MUSICスベクトルを求めることにより、到来信号の到来方向を推定する。

[0013] 本実施の形態では、このような MUSIC法にぉ、て、雑音相関行列と入力信号相関行列とを求め、入力信号相関行列に対する雑音相関行列の真固有値と真固有べタトルと偽固有ベクトルとを算出して、信号部分空間の補空間成分を求める行列を算出し、得られた行列に基づ、て信号の到来方向を推定する。

[0014] まず、本発明の原理を説明する。

[0015] ここで、 i番目の信号源から j番目の受信器までの伝達関数を並べた行列を Aとしたときに、この行列 Aの列べ外ルが張る空間を「信号部分空間」と呼ぶ。信号部分空間は、信号が存在する部分空間である。これは、雑音がない場合は、入力信号相関行列の列ベクトルが張る空間と同じものであり、入力信号相関行列の固有ベクトルのうち、 0よりも大きい固有値に対応するものが張る部分空間と一致する。直感的には、受信器アレーの指向性を信号源の方向に向けたときの受信器アレーのウェイトべタトルが作る線形部分空間と一致する。一方、既存の手法で「雑音部分空間」と呼んでいるのは、信号部分空間の補空間である。雑音部分空間は、雑音のみが存在する部分空間である。例えば、対象とする空間が 2次元空間で信号部分空間が y軸の場合、 X軸が補空間のひとつの選び方になる。

[0016] 今、説明を簡単にするために、雑音が無相関であるとすると、この二つの部分空間、つまり、信号部分空間と雑音部分空間は直交する (まさに y軸と X軸のように)。

[0017] 到来方向の推定を行うに当たり、ある方向に指向性を持つウェイトベクトルを作ったとする。このとき、そのベクトルがたまたま信号源の方向を向いていたとすると、そのべタトルは信号部分空間上のベクトルになる。これは、上記の例ではベクトルが y軸上にのっていることを意味する。したがって、そのベクトルを補空間（つまり、 X軸）に正射影すると、そのベクトルはゼロベクトルに近いものになる。 MUSICスペクトルは、このベクトルの長さの 2乗の逆数であるため、結果として、ウェイトベクトルが信号の方向を向いているとき、 MUSICスペクトルは、非常に大きな値をとることになる。よって、 M USICスペクトルが大きくなる角度をもって信号の到来方向にするの力 MUSIC法の基本的な考え方である。

[0018] 一方、雑音が無相関でない場合は、一般に補空間は信号部分空間と直交しない。

例えば、 y=—xという直線が補空間になる。しかし、この場合でも、雑音相関行列の情報を用いて、正しい信号部分空間 (y軸)を特定し、 y軸に沿って y=—xという直線に斜めに射影すれば、直交している場合と同じ議論を行うことができる。これを数学的に実現しているのが、入力信号相関行列の雑音相関行列に対する一般ィ匕固有値展開である。これは、具体的には、 y軸と y=—xとが直交するように空間自体を特定の方向につぶしてから、直交している場合と同様の操作を行うことに対応する。このとき、既存の手法で非常に重要な点は、この補空間が雑音相関行列で特定できると仮定していることである。すなわち、雑音相関行列が補空間に関するすべての情報を持っていると仮定している。これは、雑音相関行列が正則である場合に相当する。逆に、雑音相関行列が正則でない場合は、入力信号相関行列の雑音相関行列に対する一般ィ匕固有値展開というものを計算することすらできない。なぜなら、この操作には、雑音相関行列の逆行列が必要となるからである。

[0019] 本発明で問題として、るのは、雑音相関行列によって補空間（上記の例では y= - X)が特定できない場合である。具体的には、アレーを構成する受信器数よりも少ない雑音源が存在する場合などに、このような現象が生じる。この場合は、上記のように、入力信号相関行列の雑音相関行列に対する一般化固有値展開を計算することすらできない。これに対しては、無理やり補空間が直交していると仮定して X軸に正射影してしまうという方法も考えられる。この方法では、勝手に X軸に射影したベクトル力本来の補空間に射影したベクトルの長さと同じになるとは限らず、そもそも、信号部分空間を特定する段階で、大きい固有値に対応する固有ベクトルが信号部分空間となる力どうかも保証されない、つまり、信号部分空間すら正しく求められない可能性がある

[0020] このような問題が生じる原因は、ウェイトベクトルが信号でも雑音でもない方向を向いている場合を全く考慮していないためである。本来、ウェイトベクトルがそのような方向を向いている場合には、その方向も補空間であると考えなくては、 MUSIC法の原理が機能しなくなる。ここで、信号でも雑音でもない方向の規定の仕方については、入力信号相関行列で張られる部分空間でもなくかつ雑音相関行列で張られる部分空間でもない部分空間ということになる。これを数学的に記述したのが、雑音相関行列の入力信号相関行列に対する「偽固有ベクトル」である。したがって、補空間として、入力信号相関行列が張る部分空間のうち、雑音に起因する部分 (既存の手法で「雑音部分空間」と呼んでいる空間）と偽固有ベクトルによって張られる部分空間とを用いることにより、正しい「信号部分空間の補空間」が得られ、ウェイトべ外ルを当該補空間へ射影することによって、正しい MUSICスペクトルを得ることができる。

[0021] なお、数学的な定義を与えると、例えば、 Aをエルミート行列、 Bを非負定値ェルミート行列とし、スカラー λおよびベクトル w力 Aw= λ Bw、 Bw≠0を満たすとき、 λ および wを、 Aの Bに対する真固有値および真固有ベクトルという。また、 Aw=Bw= 0となるときの wを、偽固有ベクトルという。

[0022] 要するに、 MUSIC法の本質は、「信号部分空間以外の空間」を適切に推定することであり、最も単純なケースの場合、つまり、雑音同士が無相関である場合には、信号部分空間と雑音部分空間が直交することが理論的にわ力つているため、信号相関行列の固有値展開によってそれぞれを特定することができる。そして、信号部分空間と雑音部分空間が直交しない場合には、雑音の相関を打ち消すような操作を加えて、直交する場合と同様の問題に帰着させている。これを数学的に実現するのが、信号相関行列の雑音相関行列に対する一般ィ匕固有値展開である。そして、既存の手法では、どちらのケースであれ、そもそも信号部分空間でない空間はすべて雑音部分空間であると仮定している (これが、雑音相関行列が正則という条件に対応する）ため、「信号部分空間以外の空間」を雑音部分空間と呼んでいる。ところが、一般的には、「信号部分空間以外の空間」には「雑音部分空間」と「信号でも雑音でもな！ヽ部分空間」が存在しうる (これが、雑音相関行列が特異な場合に対応する。正則な場合は非常に限定された条件である)。既存の手法では、この二つを正しく特定することができない。これに対して、本発明の枠組みでは、偽固有ベクトルによって「信号でも雑音でもない部分空間」を表現できるため、雑音相関行列が正則である場合を含むすべての問題を包含する一般的な枠組みを与えることができる。なお、雑音相関行列が正則である場合、本発明の枠組みでは、偽固有ベクトルは無くなってしまうため、結果として信号相関行列の雑音相関行列に対する一般ィ匕固有値展開と同じものになる。

[0023] 図 1は、本発明の実施の形態に係る信号到来方向推定装置の構成を示すブロック図である。

[0024] 図 1に示す信号到来方向推定装置は、観測手段 1、短時間フーリエ変換手段 2、相関行列算出手段 3、雑音相関行列保持手段 4、入力信号相関行列保持手段 5、真固有値 ·真固有ベクトル ·偽固有ベクトル算出手段 6、補空間成分行列算出手段 7、到来方向探索手段 8、および到来方向出力手段 9を有する。

[0025] 観測手段 1は、到来方向推定対象である音等の入力信号を複数のマイク等により観測する。

[0026] 短時間フーリエ変換手段 2は、観測手段 1で観測した音等の入力信号を短時間フ一リエ変換する。 [0027] 相関行列算出手段 3は、短時間フーリエ変換手段 2で得られた短時間フーリエ変換された系列力相関行列を算出する。

[0028] 雑音相関行列保持手段 4は、雑音のみの環境で相関行列算出手段 3により算出された相関行列 (雑音相関行列)を保持する。

[0029] 入力信号相関行列保持手段 5は、到来方向推定対象音が存在する状況下で相関行列算出手段 3により算出された相関行列 (入力信号相関行列)を保持する。

[0030] 真固有値 ·真固有ベクトル ·偽固有ベクトル算出手段 6は、雑音相関行列保持手段 4および入力信号相関行列保持手段 5にそれぞれ保持されてヽる相関行列 (雑音相関行列、入力信号相関行列)を用いて、入力信号相関行列に対する雑音相関行列の真固有値と、真固有ベクトルと、偽固有ベクトルとを算出する。なお、具体的な算出手順の一例については、後で詳述する。

[0031] 補空間成分行列算出手段 7は、入力信号相関行列に対する雑音相関行列の真固有値と、真固有ベクトルと、偽固有ベクトルとを用いて、信号部分空間の補空間成分を求めるために使用する行列を算出する。なお、具体的な算出手順の一例については、後で詳述する。

[0032] 到来方向探索手段 8は、補空間成分行列算出手段 7で算出された補空間成分行列を用いて音の到来方向を探索する。

[0033] 到来方向出力手段 9は、到来方向探索手段 8により探索された結果得られる音の到来方向を出力する。

[0034] なお、この例では、観測手段 1として、音等の観測手段を用いているが、電波等をァレーアンテナにより観測する観測手段等、その他の信号の観測手段であってももちろん良い。

[0035] 次に、上記構成を有する本実施の形態に係る音の到来方向推定装置について、その動作を説明する。

[0036] まず、前段階として、到来方向推定対象の音が存在しない状況、つまり、雑音のみが存在する状況で観測手段 1を用いて観測を行い、短時間フーリエ変換手段 2によつてその短時間フーリエ変換を行い、相関行列算出手段 3によって雑音相関行列を算出し、雑音相関行列保持手段 4に保持する。 [0037] 具体的には、周波数インデックスを ω、時間フレームインデックスを t、観測数を nとして、短時間フーリエ変換された各観測を n次元ベクトル表記したものを χ(ω ,t)とすると、 Eを時間フレームに関する期待値、「*」を共役転置を表す記号として、雑音相関行列 Q(co)は、次の（式 1)、つまり、

[数 1]

Qico)

… （式 1 )

で表される。この雑音相関行列 Q(co)の行列群が、雑音相関行列保持手段 4に保持されること〖こなる。

[0038] 続いて、実運用状態、つまり、到来方向推定対象の音と雑音が入力されている状態で、上記と同様の操作により、相関行列算出手段 3が、入力信号相関行列を算出し、入力信号相関行列保持手段 5に保持させる。

[0039] この時、雑音と到来方向推定対象音とに由来する成分が無相関であれば、入力信号相関行列 Rx(co)は、次の（式 2)、つまり、

[数 2]

Εχ(ω) =∑[χ(ω, t)x{co, t)* J = (ω) + 0(ω) … （式 2)

と書くことができる。

[0040] ここで、 R(co)は、到来方向推定対象音のみ、つまり、入力信号のみに由来する成分の相関行列である。すなわち、雑音と到来方向推定対象音とに由来する成分が無相関であれば、入力信号相関行列 Rx(co)は、雑音相関行列 Q(co)と、到来方向推定対象音のみに由来する相関行列 R )との和になる。なお、入力信号相関行列 Rx( ω)と、雑音相関行列 Q(co)とは、実際に観測したデータ力も算出可能であるため、例えば、その差によって、到来方向推定対象音のみに由来する成分の相関行列 R(co) を求めることができる。また、 R(co)の階数は、 r(co)であるとしておく。

[0041] 次に、真固有値 ·真固有ベクトル ·偽固有べ外ル算出手段 6が、雑音相関行列保持手段 4および入力信号相関行列保持手段 5に保持されたこれらの雑音相関行列 Q (ω)、入力信号相関行列 Rx(co)を用いて、入力信号相関行列 Rx(co)に対する雑音相関行列 Q(co)の真固有ベクトル ·偽固有ベクトルを算出する。以下、この固有べタトル算出処理について説明する。なお、以下の動作は、全ての ωに対して同様に適用するので、以下、 ωは省略して記述することとする。

[0042] 具体的には、真固有値 ·真固有べ外ル'偽固有ベクトル算出手段 6は、雑音相関行列 Qと、入力信号のみに由来する相関行列 Rとを共傾変換で同時対角化する正則行列 Τを算出する。この結果、雑音相関行列 Qと、到来方向推定対象音のみ、つまり、入力信号のみに由来する相関行列 Rとは、次の（式 3)、つまり、

[数 3]

Q ^ ( }^li_{n r}r R = ( ⁱA ⁱ ··· (式 3 )

と書くことができる。

[0043] ここで、 I は、対角成分のうち最初の r個のみが 1であり、その他はすべて 0である n,r

行列を表し、 Λは対角行列となる。なお、 Tを求めるアルゴリズムは、例えば、ラオとミトラによる書籍「一般逆行列とその応用」（東京図書株式会社)の第 6章に記載されている。

[0044] この Tを用いることにより、雑音も含めた入力信号の相関行列 Rxは、雑音と到来方向推定対象音とに由来する成分が無相関であれば、上記 (式 2)で説明したように、入力信号相関行列 Rxは、雑音相関行列 Qと、到来方向推定対象音のみに由来する入力信号相関行列 Rとの和になるため、次の（式 4)、つまり、

画

Bx = {ry{A + I_{n r}) ^l … (式 4 )

であることがわかる。

[0045] なお、雑音相関行列 Qと入力信号相関行列 Rxを共傾変換で同時対角化した場合も、同じ Tが得られるため、入力信号のみに由来する相関行列 Rを求めなくても、（式 4)を得ることちでさる。

[0046] これらの結果から、雑音相関行列 Qと、雑音も含めた入力信号の相関行列 Rxは、次の（式 5)、つまり、

[数 5]

QT = RxT[\ + I_{n r}}l_{n r} … （式 5 ) を満たすことがわかる。

[0047] なお、「+」は行列のムーア ·ペンローズ型一般逆行列を表す記号とする。

[0048] そして、（Λ +1 )+の行列の対角成分のうち、 0の値がある個所を k番目とすると、 T

n,r 1

の行列の k列目の縦ベクトルが偽固有ベクトルである。この（Λ +1 )+の行列は、偽

1 n,r

固有ベクトルを特定するために使用するものである。なお、偽固有ベクトルは、複数ある可能性がある。

[0049] 一方、真固有値、つまり、（Λ +Ι )+1 の行列の対角成分に並んでいる値のうち、前

n,r n,r

記 k番目以外の対角成分で、その値が 1である個所を k番目であるとする。この (Λ

1 2

+ 1 )+1 の行列は、真固有ベクトルを特定するために使用するものである。なお、（Λ n,r n,r

+ 1 )+I の行列の対角成分が 1である箇所は、複数ある可能性がある。

n,r n,r

[0050] すなわち、 Tの行列の k列目の縦ベクトルが雑音部分空間の基底をなす真固有べ

2

タトルに対応し、それ以外の真固有ベクトルが信号部分空間の基底をなし、雑音部分空間の基底も、信号部分空間の基底も、複数個ある可能性がある

[0051] 例えば、 Λの行列の対角成分が、例えば、 [1, 0, 4, 0]であり、 I の行列の対角成

n,r

分が、例えば、 [1, 1, 0, 0]であったとする。

[0052] この場合、（Λ +Ι )+の行列の非対角成分は、理論的に 0になる一方、（Λ +Ι )+の

n,r n,r 対角成分は、 [0.5, 1, 0.25, 0]となり、（Λ +Ι )+の行列の対角成分のうち 0の値が

n,r

ある個所は 4番目であり、 k =4となる。すなわち、 Tの行列の 4列目の縦ベクトルが偽固有ベクトルに対応することになる。

[0053] これに対し、（Λ +Ι )+ X I の行列の対角成分は、 [0.5, 1, 0, 0]となるため、前記

n,r n,r

k番目以外の対角成分で、その値が 1である個所を k番目であるとすると、 k = 2と

1 2 2 なる。すなわち、 Tの行列の 2列目の縦ベクトルが、雑音部分空間に対応する真固有ベクトルに対応することになる。

[0054] そして、 Tの行列の残りの 1列目と、 3列目の縦ベクトルは、信号部分空間に対応する真固有ベクトルに対応することになる。

[0055] ここで、雑音部分空間に対応する真固有ベクトルを求めるのに、一見、偽固有べクトルを求めるために使用する前者の (Λ +Ι )+の行列を使っても良いように見える力

n,r

Λの中の 4が 1だった場合は、雑音に由来する 1なの力、信号に由来する 1なのかが判別できなくなる。そのため、雑音部分空間に対応する真固有ベクトルを求める場合に、 (Λ +Ι )+の行列を使用するのは不適当である。

n,r

[0056] このようにして、真固有値'真固有ベクトル '偽固有ベクトル算出手段 6は、雑音相関行列 Q( ω )および入力信号相関行列 Rx( ω )を用、て、入力信号相関行列 Rx( ω )に対する雑音相関行列 Q( co )の真固有値と、真固有ベクトルと、偽固有ベクトルとを求めると共に、その真固有値に基づき真固有ベクトルを切り分けて、雑音のみに依存する部分空間である雑音部分空間に対応する真固有ベクトルと、入力信号に依存する部分空間である信号部分空間に対応する真固有ベクトルとを算出する。

[0057] 次に、補空間成分行列算出手段 7は、信号部分空間に対応する真固有ベクトルを、次の（式 6)、つまり、

[数 6]

u_j， (i = l₇.., _s ) … (式 6 )

とおき、雑音部分空間に対応する真固有ベクトルを、次の (式 7)、つまり、 [数 7]

u,^N, = 1,..,»J … （式 7 )

とおき、偽固有ベクトルを、次の（式 8)、つまり、

[数 8]

1, … （式 8 )

とおき、次の（式 9)、つまり、

[数 9]

W =

■■■ (式 9 ) を構成し、保持する。

[0058] なお、n=n +n +nである。

S N I

[0059] 続いて、補空間成分行列算出手段 7は、上記 Wの随伴行列を算出し、それを、次の（式 10)、つまり、

[数 10]

* =卜 f,… , '…，， …， … （式 1 0 ) とおいたものから、次の（式 11)、つまり、

[数 11] = ..., ".. 」* … （式 1 1 )

のみを抽出した行列を、信号部分空間の補空間成分を求める行列として算出する。

[0060] そして、到来方向探索手段 8が、音の到来方向を探索する。具体的には、長さ 1であり、信号の到来方向 Θで決まるマイクロホンアレーに対するウェイトベクトルを a( Θ ) とし、次の（式 12)、つまり、

[数 12] … （式 1 2 )

力 Sピークを持つ Θを探索する。

[0061] 最後に、到来方向出力手段 9が、到来方向探索手段 8により探索された ΡΜ ( Θ )がピークを有するときの Θを、入力信号の到来方向として出力する。

[0062] なお、到来方向探索手段 8における PM ( Θ )はあくまで一例であり、 X、つまり、信号部分空間に対応する真固有ベクトルと、雑音部分空間に対応する真固有ベクトルと、偽固有ベクトルとを用いて信号の到来方向 Θを一意に決定することができる到来方向探索手段 8であれば何を用いても良、。

[0063] このように、本実施の形態によれば、入力信号相関行列に対する雑音相関行列の真固有値と真固有ベクトルと偽固有ベクトルとを求め、その真固有値と真固有べタトルと偽固有ベクトルとに基づき信号部分空間の補空間成分を求める行列を算出し、その行列に基づ!、て信号の到来方向を推定するため、入力信号相関行列や雑音相関行列が特異である場合でも、信号到来方向を推定することができる。

[0064] なお、本実施の形態では、図 1に示すように、信号到来方向推定装置の構成をプロック図により示してハードウェア的に説明したが、本発明は、これに限らず、この信号到来方向推定装置を、 CPUと、その CPUを上記のように実行させるためのプロダラムを記憶したハードディスクや、メモリ等の記憶装置を有する PC等の汎用のコンビュータがソフトウェア的に実行するようにしても勿論良い。この場合、上記の信号到来方向推定装置としての機能を果たすためのプログラムは、 CD等の記録媒体に記録されていて、その記録媒体力読み出してコンピュータ内の記憶装置にインストールしたり、またはインターネット等のネットワークを介してサーバ等力ダウンロードして記憶装置に記憶させる。このようにすれば、 PC等の汎用のコンピュータでも、上記の信号到来方向推定装置としての機能を果たすためのプログラムを実行することにより、上記の信号到来方向推定を実行することができる。

[0065] 本明細書は、 2004年 9月 14日出願の特願 2004— 267380に基づく。この内容はすべてここに含めておく。

産業上の利用可能性

[0066] 本発明に係る信号到来方向推定装置、信号到来方向推定方法、信号到来方向推定用プログラムは、雑音相関行列の入力信号相関行列に対する真固有値、真固有ベクトル、偽固有ベクトルを用いることにより、従来法では対応できなかった、相関行列が特異な場合についても、信号の到来方向を推定することができる、という効果を有し、移動通信や室内無線通信 (無線 LAN)などにおいて電波の到来方向を推定する、例えば、 MUSIC法を採用した信号到来方向推定装置、信号到来方向推定方法、信号到来方向推定用プログラムに有用である。

Claims

請求の範囲

[1] 信号の到来方向を推定する信号到来方向推定装置であって、

到来方向推定対象の信号が存在しない雑音のみが存在する状況下における雑音相関行列と、到来方向推定対象の信号および雑音が存在する状況下における入力信号相関行列とを算出する相関行列算出手段と、

前記雑音相関行列および前記入力信号相関行列を用、て、前記入力信号相関行列に対する前記雑音相関行列の真固有値と真固有ベクトルと偽固有ベクトルとを算出する真固有値'真固有ベクトル '偽固有ベクトル算出手段と、

前記真固有値と前記真固有ベクトルと前記偽固有ベクトルとに基づいて、前記信号の到来方向を推定する到来方向推定手段と、

を有する信号到来方向推定装置。

[2] 前記真固有値 ·真固有ベクトル ·偽固有ベクトル算出手段は、

前記真固有ベクトルから、前記真固有値に基づいて、信号部分空間に対応する真固有ベクトルと、雑音部分空間に対応する真固有ベクトルとを算出し、

当該信号到来方向推定装置は、

前記偽固有ベクトルと、前記信号部分空間に対応する真固有ベクトルと、前記雑音部分空間に対応する真固有ベクトルとを用いて、前記信号部分空間の補空間成分を求める行列を算出する補空間成分行列算出手段をさらに有し、

前記到来方向推定手段は、

前記信号部分空間の補空間成分を求める行列を用いて、前記真固有値と前記真固有べクトルと前記偽固有べクトルとに基づ!/、て、前記信号の到来方向を推定する、請求項 1記載の信号到来方向推定装置。

[3] 前記信号部分空間の補空間を求める行列を、 MUSIC (Multiple Signal Classificati on)法における、雑音部分空間の基底力もなる行列の代わりに用いる、

請求項 2記載の信号到来方向推定装置。

[4] 前記到来方向推定対象の信号は、電波または音である、請求項 1記載の信号到来方向推定装置。

[5] 信号の到来方向を推定する信号到来方向推定方法であって、到来方向推定対象の信号が存在しない雑音のみが存在する状況下における雑音相関行列と、到来方向推定対象の信号および雑音が存在する状況下における入力信号相関行列とを算出する工程と、

前記雑音相関行列および前記入力信号相関行列を用、て、前記入力信号相関行列に対する前記雑音相関行列の真固有値と真固有ベクトルと偽固有ベクトルとを算出する工程と、

前記真固有値と前記真固有ベクトルと前記偽固有ベクトルとに基づいて、前記信号の到来方向を推定する工程と、

を有する信号到来方向推定方法。

信号の到来方向を推定させるための信号到来方向推定用プログラムであって、コンピュータに、

到来方向推定対象の信号が存在しない雑音のみが存在する状況下における雑音相関行列と、到来方向推定対象の信号および雑音が存在する状況下における入力信号相関行列とを算出する工程と、

前記真固有ベクトルと前記偽固有ベクトルとに基づいて、前記信号の到来方向を推定する工程と、

を実行させるための信号到来方向推定用プログラム。