JP4722132B2 - 到来波数推定方法、到来波数推定装置及び無線装置 - Google Patents
到来波数推定方法、到来波数推定装置及び無線装置 Download PDFInfo
- Publication number
- JP4722132B2 JP4722132B2 JP2007524485A JP2007524485A JP4722132B2 JP 4722132 B2 JP4722132 B2 JP 4722132B2 JP 2007524485 A JP2007524485 A JP 2007524485A JP 2007524485 A JP2007524485 A JP 2007524485A JP 4722132 B2 JP4722132 B2 JP 4722132B2
- Authority
- JP
- Japan
- Prior art keywords
- matrix
- correlation
- estimation
- incoming
- antenna
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/74—Multi-channel systems specially adapted for direction-finding, i.e. having a single antenna system capable of giving simultaneous indications of the directions of different signals
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
- H04B7/0842—Weighted combining
- H04B7/086—Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Radar Systems Or Details Thereof (AREA)
Description
信号の到来方向推定問題に関して、推定精度と計算演算量などの立場から信号部分空間と雑音部分空間の直交性を利用した部分空間手法(Subspace-based method)がよく知られている。その代表例として、MUSIC (Multiple signal classification) である(非特許文献1参照)。また、完全な相関性をもつ多重波の到来方向推定問題への対応策として、空間スムージングを用いた部分空間手法(Subspace-based method with spatial smoothing) は空間スムージングMUSIC (Spatial smoothing based MUSIC)としてよく知られている(非特許文献2および非特許文献3参照)。これら従来の部分空間手法は、アレー共分散行列の固有値分解(Eigenvalue decomposition: EVD) 或いはアレーデータ行列の特異値分解(Singular value decomposition: SVD)により信号部分空間或いは雑音部分空間を得るため、アレーに入射する信号の個数に関する情報を必要とする。このため、固有値分解や特異値分解に限らずアレーの受信データから到来信号の個数を推定することは高分解能をもつ到来方向推定手法の絶対必要な条件である。
length)という情報理論基準(information theoretic criteria)に基づく推定手法(非特許文献4参照) がよく知られている。また、完全な相関性をもつ多重波の個数推定問題への対応策として、空間スムージング(spatial smoothing(SS))を用いた情報理論基準に基づく推定手法もよく知られている。その代表例として、SS-AIC やSS-MDLがある (非特許文献4及び非特許文献5を参照)。
無相関信号の個数を推定するAIC やMDL 手法では、アレーアンテナの素子の受信データからアレー共分散行列を求め、この共分散行列の固有値分解を行い、小さい固有値の数が信号の個数に関連するという特性を利用して信号数を推定する。また、相関性をもつ信号(完全な相関性を持つ多重波を含む)の到来方向に関しては、入射信号間の相関を抑圧するために線形等間隔アレーをサブアレー化し、各サブアレーの共分散行列の平均操作を行い、空間的に平均された共分散行列の小さい固有値の数が信号の個数に関連するという特性を利用して相関信号の個数を推定する。
ここで,M 個アレー素子数をもつ線形等間隔アレーにp 個の多重波狭帯域信号[sk(n)] が角度[θk ] からアレーアンテナに入射しているとする。各素子のアレー受信信号は次式 のように表せる。
Rs=E[s(n)sH(n)]は入射する多重波の共分散行列であり、IMはM×Mの単位行列である。さらに、観測データyi(n),yk(n)の相関rikをrik=E{yi(n) y*k(n)}で定義する。ただし、rik=r*kiという関係が存在する。(・)*は複素共役を表す。また、(2)式 内のアレーの共分散行列Rは次式で明確的に表現できる。
のサブアレーの受信ベクトル
は(4)式で表現できる。
である。また、am(θk) とAmはサブアレーの応答ベクトルと応答行列である。従って、
のサブアレーの共分散行列は(5)式 で与えられる。
を空間的に平均すると,(6)式 のような共分散行列が得られる。
の固有ベクトルと固有値であり、Eは[ei]を列とする行列、
は[λi]を要素とする対角行列である。また、λ1≧λ2≧・・・≧λp>λp+1=・・・=λm=σ2である。従って,最小固有値の数a(=m−p )から入射信号の個数p(=m−a)を推定できる.このことを利用して,サンプル時刻n=1, 2,・・・ , N までの受信ベクトル
を用いて、各サブアレーのサンプル共分散行列を次式で求める。すると、N=∞である(5)式は次式で表現できる。
の固有値分解を(9)式 のように計算できる。
を用いて、信号数を推定するAIC とMDL基準は次式与えられる。
を得るために空間的に平均されたアレー共分散行列
の固有値分解を行う必要性がある。
しかし、従来のAIC やMDL 手法では、アレー素子の数が多い時また実時間処理で変化する到来方向を推定する時、必須の固有値分解処理や特異値分解処理が複雑となって計算量が大量となり、計算時間がかなりかかる。このため、従来の固有値分解に基づく入射信号数推定手法の実際への応用は、計算の負担となる固有値分解に制限されてしまい、高速かつ高精度でアレーに入射する信号の個数と到来方向を推定できない問題が生じる。
また、従来のAIC やMDL 手法ではアレーアンテナの受信データ長が短い場合や信号対雑音比(Signal-to-noise ratio: SNR)が低い場合に高精度でアレーに入射する信号の個数と到来方向を推定できない問題がある。
また、信号到来方向を正確に推定できなくなると、精確な基地局の受信/送信ビーム形成できなくなり、基地局の受信及び送信システムの性能が劣化する。
従って、本発明の目的は、時空間的に無相関な白色雑音環境または空間的に相関する雑音環境において、線形等間隔アレー(Uniform linear array, ULA) に入射する狭帯域信号に対して、固有値分解を利用しない少ない演算量で到来信号の個数を正確に推定できる手法を提供することである。なお、対象とする狭帯域信号は、例えば、無相関或いは相関或いは完全相関(即ち多重波)、または一部相関または一部完全相関の狭帯域信号である。
本発明の別の目的は、アレーアンテナの受信データ長が短い場合でも、また信号対雑音比SNRが低い場合でも、固有値分解を利用しない少ない演算量で到来信号の個数を正確に推定できる手法を提供することである。
本発明の別の目的は、空間的に角度が近い複数の信号が入射しても高い検出確率で到来信号の個数を正確に推定できる手法を提供することである。
R.O. Schmidt, "Multiple emitter location and signal parameter estimation," IEEE Trans. Antenna Propagation, vol. 34, no. 3, pp. 276-280 (1986) J. Shan, M. Wax and T. Kailath, "On spatial smoothing for direction-of-arrival estimation of coherent signals," IEEE Trans. Acoust., Speech, Signal Processing, vol. 33, no. 4, pp. 806-811 (1985) S.U. Pillai and B.H. Kwon, "Forward/backward spatial smoothing techniques for coherent signals identification," IEEE Trans. Acoust., Speech, Signal, vol. 37, no. 1, pp. 8-15(1989) Wax and T. Kailath、 "Detection of signals by information theoretic criteria, " IEEE Trans. Acoust. Speech, Signal Processing, vol. 33、 no. 2, pp. 387-392 (1985) T.-J. Shan, M. Wax and T. Kailath, "On spatial smoothing for direction-of-arrival estimation of coherent signals," IEEE Trans. Acoust., Speech, Signal Processing, vol. 33, no. 4, pp. 806-811 (1985)
本発明の到来波数推定方法は、前記アレーアンテナの各アンテナ素子の受信信号と他アンテナ素子の受信信号間の相関を算出するステップ、各アンテナ素子の受信信号と他アンテナ素子の受信信号間の相関を行列に配列してなるアレー共分散行列を構成する所定の行あるいは列より対角要素を除き、該対角要素が除かれた(M−1)個の相関より所定数
個の相関を順次1要素づつずらしながら取り出して行列に配列してなる相関行列を作成するステップ、該相関行列を用いて到来波数を推定するための推定行列を作成するステップ、該推定行列をQR分解して得られた上三角行列因子の各行要素に基づいて到来電波の個数を推定するステップ、を備えている。
上記の相関行列を作成するステップでは、q(0≦q≦M−1)を雑音の空間的相関の長さとするとき、前記アレー共分散行列を構成する所定の行あるいは列より対角要素および該対角要素に連鎖的に隣接するq個の要素を除き、(q+1)個の要素が除かれた(M−q−1)個の相関より、所定数の相関を取り出し、順次1要素づつずらしながら前記所定数の相関を取り出して行列に配列してなる相関行列を作成する。
本発明の到来波数推定装置は、アレーアンテナの各アンテナ素子の受信信号と他アンテナ素子の受信信号間の相関を算出する相関演算部、各アンテナ素子の受信信号と他アンテナ素子の受信信号間の相関を行列に配列してなるアレー共分散行列を構成する所定の行あるいは列より対角要素を除き、該対角要素が除かれた(M−1)個の相関より所定数の相関を順次1要素づつずらしながら取り出して行列に配列してなる相関行列を作成する相関行列作成部、 該相関行列を用いて到来波数を推定するための推定行列を作成する推定行列作成部、該推定行列をQR分解するQR分解部、QR分解により得られた上三角行列因子の各行要素に基づいて到来電波の個数を決定する到来波数決定部を備えている。
上記の相関行列作成部は、q(0≦q≦M−1)を雑音の空間的相関の長さとするとき、前記アレー共分散行列を構成する所定の行あるいは列より対角要素および該対角要素に連鎖的に隣接するq個の要素を除き、(q+1)個の要素が除かれた(M−q−1)個の相関より、所定数の相関を順次1要素づつずらしながら取り出して行列に配列して前記相関行列を作成する。
また、前記相関行列作成部は、前記アレー共分散行列の先頭行及び最終行並びに先頭列及び最終列のうち1以上について前記相関行列を作成し、前記推定行列作成部は、該作成された相関行列を並べることにより前記推定行列を作成する。
前記到来波数決定部は、前記上三角行列因子の第i行目の各行列要素の絶対値の総和と第(i+1)行目の各行列要素の絶対値の総和の比を計算する計算部、該比が最小のiを到来電波の個数と決定する決定部を有している。
本発明の無線装置は、到来電波の個数を推定する到来波数推定部、アンテナ素子毎の受信データと到来電波の個数を用いて信号の到来方向を推定する到来方向推定部、到来方向の推定値を用いて信号源方向にピークを有するように受信ビームフォーミング処理を行う受信ビームフォーマを備え、前記到来波数推定部は、
前記アレーアンテナの各アンテナ素子の受信信号と他アンテナ素子の受信信号間の相関を算出する相関演算部、各アンテナ素子の受信信号と他アンテナ素子の受信信号間の相関を行列に配列してなるアレー共分散行列を構成する所定の行あるいは列より対角要素を除き、該対角要素が除かれた(M−1)個の相関より所定数の相関を順次1要素づつずらしながら取り出して行列に配列してなる相関行列を作成する相関行列作成部、該相関行列を用いて到来波数を推定するための推定行列を作成する推定行列作成部、該推定行列をQR分解するQ
R分解部、QR分解により得られた上三角行列因子の各行要素に基づいて到来電波の個数を決定する到来波数決定部を備えている。
本発明はアレーアンテナを用いて電波到来方向を推定する基地局における到来電波の個数推定方法及び推定装置に関するものであり、図面に従って第1実施例の到来電波の個数推定方法及び推定装置について説明する。なお、以降の図において、概略同じ物あるいは同じ機能を有するものについては同じ符号を付する。
図1は距離dの間隔で直線的にM個のアンテナ素子を配列したアレーアンテナの構成図である。図2は送信源10と基地局受信アンテナ(アレーアンテナ)30との配置関係図である。アレーアンテナ30は図1に示すように等間隔直線アレーアンテナ構成を有し、多重波到来方向推定システムを構成する。図2において、送信源10からアレーアンテナ30にまっすぐ入射するものは直接波11であり、建物BL1,BL2などによって反射されてからアレーアンテナ30に入射するものは反射波12である。図2では一つ例として、二つ反射波を示すが、以下では送信源10からの直接波と反射波の総個数はp(実際には未知数)とする。また、M>2pである。さらに、直接波と反射波の関係は、次式
図4は到来波数推定部50の構成図である。この到来波数推定部50は、アレーデータ間の相関を計算する相関計算手段51、Hankel相関行列形成部52、到来電波の個数を推定するための推定行列の形成手段53、該推定行列をQR分解するQR 分解手段54,QR 分解により得られた上三角行列(upper triangular matrix) 因子より評価基準を計算する評価基準計算手段55及び到来波の個数を決定する個数決定手段56で構成されている
一般には、アレーアンテナ30により受信した信号から電波到来方向を推定する際に、アンテナ素子受信信号ベクトルy(n)(=[y1(n),y2(n),…,yM(n)]T)の各受信信号間の相関r11〜rMMを演算して行列に配列したアレー共分散行列Rが用いられる。このアレー共分散行列Rは、受信信号ベクトルy(n)の複素共役をyH(n)とすると、無相関白色雑音環境において次式で与えられる。
yi(n)=xi(n)+wi(n)
E[ wi(n) wj*(n)]=σ2 (i=j)
E[ wi(n) wj*(n) ]=0 (i≠j)
である。すなわち、無相関白色雑音環境ではアレー共分散行列Rの対角要素r11,r22,...rMMに雑音が含まれている。
さて、アレー共分散行列Rは共役対称になっているから、到来波の個数や到来方向の推定に際して任意の行、列、例えば図5に示すようにその1列目と最終列、あるいは図6に示すように1行目と最終行を計算するだけで十分である。ただし、対角要素には前述のように雑音が含まれているから図5、図6に示すように各列、各行から対角要素r11,rMMを除外する。
を用いて、(14)式により信号y(n)とy*M(n),およびy(n)とy*1(n)の相関ベクトル
のHankel相関行列
を形成する。すなわち、相関計算手段51は、アレーアンテナにおけるM番目のアンテナ素子の受信信号と第1,2,…M−1番目のアンテナ素子の受信信号間の (M−1)個の相関
(図5の最終列参照)を算出する。なお、Mはアレーアンテナの素子数、Lはサブアレーの数、
はサブアレーのサイズ(即ち,サブアレーの素子数)で、M/2より大きくない最大整数であるとする。すなわち、
である。
個の相関を取り出し、順次上から下方向に1個づつずらしながら
組の相関を取り出し、取り出した各組の相関を1行目から順番に行列に配列してなる
の相関行列
を作成する。
同様に、相関計算手段51は、前記アレーアンテナにおける1番目のアンテナ素子の受信信号と第2,3,…M番目のアンテナ素子の受信信号間の (M−1)個の相関
(図5の第1列参照)を算出する。ついで、Hankel相関行列形成部52は、図5の右下に示すように、第1列目の(M−1)個の相関よりサブアレーサイズ
個の相関を取り出し、順次上から下方向に1個づつずらしながら
組の相関を取り出し、取り出した各組の相関行列を1行目から順番に行列に配列してなる
の相関行列
を作成する。
(図6の第1行目参照)を算出する。ついで、Hankel相関行列形成部52は、図6の右上に示すように、第1行目の(M−1)個の相関よりサブアレーサイズ
個の相関を取り出し、順次右から左に1個づつずらしながら
組の相関を取り出し、取り出した各組の相関行列を1行目から順番に行列に配列してなる
の相関行列
を作成する。
同様に、相関計算手段51は、前記アレーアンテナにおけるM番目のアンテナ素子の受信信号と第1,2,…M−1番目のアンテナ素子の受信信号間の (M−1)個の相関
(図6の最終行参照)を算出する。ついで、Hankel相関行列形成部52は、図6の右下に示すように最終行の(M−1)個の相関よりサブアレーサイズ
個の相関を取り出し、順次右から左に1個づつずらしながら
組の相関を取り出し、取り出した各組の相関行列を1行目から順番に行列に配列してなる
の相関行列
を作成する。
の行列である。(14a),(14b)式より、データの長さN が十分大きければ(A)に示した4つの各Hankel相関行が雑音の影響を受けないことは明らかである。また、
の相関行列
を求めて次式
で表現される。ここで直交行列とは任意の2つの行の対応する要素を掛け合わせた総和が0となる行列であり、上三角行列とは対角線より下の行列要素がすべて0となる行列である。したがって、(17)式において、
はそれぞれ
のユニタリー行列(直交行列)、上三角行列となり、上三角行列の各要素次のような値をとる。
のランク(階数)は到来波の個数
となる。したがって、雑音のない場合または受信データの長さNが十分大きい理想的な場合には上三角行列の下から
個の行の要素が0となる。
次に、評価基準計算手段55 は、QR 上三角行列因子の第i行の要素を用いて補助量ζ(i)
を次式
個の補助量ζ(i)は0となるから、ζ(i)=0となる行を求め、その数より到来波の個数
を求めることができる。しかし、雑音が存在する場合には厳密にζ(i)=0とならないため、以下の処理を行う。すなわち、到来波数決定手段56は、評価基準ξ(i) を次式
であり,この評価基準ξ(i) を最大にするiを到来波数の推定値
とする。すなわち、
以下、計算機シミュレーションの具体例を通してさらに説明する。シミュレーションにおいて、アンテナ素子数M は10であり、また,同じパワーをもつ2 つ多重波(p=2 は、到来方向θ1=50 とθ2=120でアレーアンテナに入射するものとする。
まず,信号対雑音比(SNR) に対する本発明に基づく多重波到来波の個数推定性能を考察する。ここで、受信データの長さNを128とし、SNR を‐10 から25dB までの範囲に変化させる。各SNR に対して1000 回の計算を行う。図7 の曲線Aは到来波の個数を正確に検出した本発明の確率(検出確率)を示す特性図であり、横軸はSNR(dB)、縦軸は検出確率である。なお、比較するために従来のSS-AIC (曲線B)とSS-MDL (曲線C)を用いた結果をプロットして示している。
図7から明らかなように本発明に基づく到来波の個数推定手法は、相関行列
を構成するアレーデータの相関
を適切に選択することにより,観測雑音の影響を除去できるので、本発明の手法の検出性能は固有値分解を用いたSS-AIC とSS-MDL 法よりはるかに良くなっている。また、本発明の手法の演算量は固有値分解処理が不要のためかなり低減となっている.
次に、受信データの長さNに対する本発明に基づく多重波到来方向の推定性能を示す。ここで、SNR を5dB とし、受信データの長さN を10 〜1000 の範囲に変化させる。図8 の曲線Aは本発明において1000回のシミュレーション計算で得られた個数推定の検出確率を示す。図8 から明らかなように、受信データが少ない場合でも,本発明の方法は従来のSS-AIC(曲線B) とSS-MDL (曲線C)より高い検出確率を示している。
図9は多重波の個数推定の性能を示すシミュレーション結果で、到来波の個数を正確に検出した確率(検出確率)を示している。図において、横軸はサブアレーサイズ
、縦軸は到来波の個数を正確に検出した確率である。ただし、アンテナ素子数M=10、受信データ長N=128 、SNR =5 dBとし、サブアレーサイズを2 から7までの範囲で変化させ、各サブアレーサイズ に対して1000 回の計算を行ったときの到来波の個数を正確に検出した確率を示している。なお、実際の到来波の個数は2である。図9より、
と仮定して個数推定すれば正確に到来波の個数検出することが可能となる。
上述の第1実施例では、次式
により推定行列を形成して、到来波の個数推定方法を説明したが、(16)式に関連して説明したようにHankel相関行列
の任意の組み合わせで推定行列
を形成し、無相関白色雑音環境における完全相関信号すなわち多重波、または部分相関信号、または無相関信号の到来個数を推定することができる。
上記(A)における4つの相関行列のうち、任意の1つを選んで、推定行列を形成する。すなわち、次式
(b) 任意の二組の相関行列を用いる実施例
上記(A)における4つの組のうち、任意の2つの相関行列を選んで、推定行列を形成する。例えば、次式
上記(A)における4つの組のうち、任意の3つの相関行列を選んで、推定行列を形成する。例えば、次式
第1、第2実施例は、無相関白色雑音環境における到来電波の個数推定の実施例であり、第iアンテナ受信素子と第jアンテナ受信素子の受信信号に含まれる雑音wi(n)、 wj(n)の相関が次式
E[ wi(n) wj*(n)]=σ2 (i=j)
E[ wi(n) wj*(n)]=0 (i≠j)
で表現される場合である。すなわち、雑音の空間的相関の長さが1の場合である。第3実施例では雑音の空間的相関の長さがq(>1)の場合の実施例である。雑音の空間的相関の長さがqの場合、第iアンテナ受信素子と第jアンテナ受信素子の受信信号に含まれる雑音wi(n)、 wj(n)の相関は次式
E[ wi(n) wj*(n)]≠0 (|i−j|≦q)
E[ wi(n) wj*(n) ]=0 (|i−j|>q)
で表現される
すなわち、雑音の空間的相関の長さがqの場合、アレー共分散行列Rの対角要素r11,r22,…,rMMおよび該対角要素からqの距離にある相関要素には雑音が含まれている。
アレー共分散行列Rは共役対称になっているから、到来方向の推定には該アレー共分散行列Rの任意の行、列、例えば、図10に示すようにその1列目と最終列、あるいは図11に示すように1行目と最終行を計算するだけで十分である。ただし、対角要素および該対角要素からqの距離にある相関要素には雑音が含まれているから図10に示すように各列から要素
を除外する。また、図11に示すように、各行から要素
を除外する。
(図10の最終列参照)を算出する。ただし、
である。ついで、Hankel相関行列形成部52は図10の右上に示すように、最終列の(M−q−1)個の相関よりサブアレーサイズ
個の相関を取り出し、順次上から下方向に1個づつずらしながら
組取り出し、取り出した各組の相関を1行目から順番に行列に配列してなる
の相関行列
を作成する。
同様に、相関計算手段51は、前記アレーアンテナにおける1番目のアンテナ素子の受信信号と第q+2,q+3,…M番目のアンテナ素子の受信信号間の (M−q−1)個の相関
(図10の第1列参照)を算出する。ついで、Hankel相関行列形成部52は、図10の右下に示すように、第1列目の(M−q−1)個の相関よりサブアレーサイズ
個の相関を取り出し、以後、順次上から下方向に1個づつずらしながら
組の相関を取り出し、取り出した各組の相関を1行目から順番に行列に配列してなる
の相関行列
を作成する。
しかる後、行についても同様の操作を行なって、図11に示すように相関行列
を作成する。
4つのHankel相関行列が求まれば、推定行列作成部53は次式
また、第2実施例と同様に、4つのHankel相関行列
の任意の組み合わせで推定行列
を形成し、空間相関雑音環境境における多重波の到来個数を推定することができる。
以上の実施例では、(17)式に示したQR分解にしたがって到来波の個数を推定したが、(17)式QR分解に代わって列回転QR分解(QR decomposition with column pivoting)することもできる。列回転QR分解では、
の置換行列(permutation matrix)Πを用いて次式
(26)式の列回転QR 分解において置換行列Πを単位行列にすれば、たとえば、
とすれば(17)式のQR 分解となる。この列回転QR 分解では、上三角行列の対角要素の絶対値を減少するように置換行列Π を決定できる(QRP と略称する)。また、行列
の列指数最大差異二分法(column index maximum-difference bisection)で事前決定した置換行列Π を使用することができる(QRPP と略称する)。例えば、
の場合、
とする。ただし、eiは
の単位ベクトルであり、その第i番目の要素=1、他の要素=0である。
図12は列回転QR 分解によるシミュレーション結果であり、測定条件は図7のシミュレーションと同一であり、Aは第1実施例の検出確率特性、BはQRPによる検出確率特性、CはQRPPによる検出確率特性である。
・基地局受信装置
図13は本発明を適用した基地局受信装置の構成図である。アレーアンテナ30は信号を受信し、べースバンド及びディジタル処理部40に入力する。ディジタル処理部40はアンテナ素子毎に信号処理して複素ディジタル受信データを出力する。到来波数推定部50は第1実施例〜第4実施例で説明した方法により、到来電波の個数pを推定する。到来方向推定部60は各アンテナ素子毎の複素ディジタル受信データと到来電波の個数p用いて信号の到来方向を推定する。なお、到来方向推定方法は例えばWO2005/001504 A1およびPCT/JP04/011598(出願日2004.8月12日)を参照されたい。
ビーム形成器(受信ビームフォーマ)70は到来方向推定部60から取得した信号の到来方向の推定値を用いて信号源方向にピークを有するようにビームを形成する(受信ビームフォーミング)。すなわち、ビーム形成器70は干渉や雑音などを抑圧しながら希望信号を抽出してチャネル受信部80に送る。チャネル受信部80は周知の方法で受信処理を行って受信データを復調、出力する。
なお、ビーム形成器70として種々の構成が可能であるが、例えば、O.L. Frost、 "An algorithm for linearly constrained adaptive array processing、" Proc. IEEE、 vol.
60、no. 8、 pp. 926-935 (1975)及びJ. Xin、 H. Tsuji、 Y. Hase、 and A. Sano、 "Array beamforming based on cyclic signal detection、" Proc. IEEE 48th Vehicular Technology Conference、 pp. 890-894、 Ottawa、Canada (1998)などに記載されたビーム形成手法を活用して、希望の信号到来方向にビームを向けて受信することが可能である。
図14は基地局送信装置の構成図である。なお、図14には基地局受信装置も図示している。
送信ビームフォーマ90は、送信部100から送信データが入力されると、到来方向推定部60により推定された方向にピークが向くように送信ビームパターンを形成し、複素ディジタル送信信号をべースバンド及びデジタル信号処理部40′に入力する。信号処理部40′は複素ディジタル送信データを無線信号に変換してアレーアンテナ30′の各アンテナ素子に入
力する。この結果、受信局に向けてビームが発射され、誤り率を低下できる。図14のアレーアンテナ30,30′を共通化することができる。
なお、本発明は多重波或いは部分相関或いは無相関信号の到来方向をオンラインで推定/追尾する際に適用できるものである。
以上説明したように、本発明によれば、固有値分解を利用せず少ない演算量で空間的に角度が近い到来波の個数推定を行うことができる。また、計算機シミュレーションの具体例で明らかにしたように、本発明によれば、短い受信データ長であっても、また、低い信号対雑音比(SNR)であっても、基地局におけるアレーアンテナに入射する到来波の個数を高い正確検出確率で推定することができる。
また、本発明によれば、固有値分解を用いた従来のSS-AIC やSS-MDL 手法に比較して少ない演算量で優れた推定性能が得られる。すなわち、従来技術に比べて、多重波の到来方向の個数と推定精度の向上を図ることができる。
更に、本発明によれば、信号源方向にビームが向くようにビームフォーミングができ、受信時における品質を向上でき、しかも、送信時において希望方向に指向性をもつビームを形成できるため可能な限り他の干渉とならないようにできる。
Claims (10)
- 複数個(=M)のアンテナ素子を同一素子間隔で直線状に配列したアレーアンテナで到来電波を受信し、該到来電波の個数を推定する到来波数推定方法において、
前記アレーアンテナの各アンテナ素子の受信信号と他アンテナ素子の受信信号間の相関を算出し、
各アンテナ素子の受信信号と他アンテナ素子の受信信号間の相関を行列に配列してなるアレー共分散行列を構成する所定の行あるいは列より対角要素を除き、該対角要素が除かれた(M−1)個の相関より所定数
個の相関を順次1要素づつずらしながら取り出して行列に配列してなる相関行列を作成し、
該相関行列を用いて到来波数を推定するための推定行列を作成し、
該推定行列をQR分解して得られた上三角行列因子の各行要素に基づいて到来電波の個数を推定する、
ことを特徴とする到来波数推定方法。 - q(0≦q≦M−1)を雑音の空間的相関の長さとするとき、前記アレー共分散行列を構成する所定の行あるいは列より対角要素および該対角要素に連鎖的に隣接するq個の要素を除き、(q+1)個の要素が除かれた(M−q−1)個の相関より、所定数の相関を順次1要素づつずらしながら前記所定数
個の相関を取り出して行列に配列してなる前記相関行列を作成する、
ことを特徴とする請求項1記載の到来波数推定方法。 - 複数個(=M)のアンテナ素子を同一素子間隔で直線状に配列したアレーアンテナで到来電波を受信し、該到来電波の個数を推定する到来波数推定装置において、
前記アレーアンテナの各アンテナ素子の受信信号と他アンテナ素子の受信信号間の相関を算出する相関演算部、
各アンテナ素子の受信信号と他アンテナ素子の受信信号間の相関を行列に配列してなるアレー共分散行列を構成する所定の行あるいは列より対角要素を除き、該対角要素が除かれた(M−1)個の相関より所定数の相関を順次1要素づつずらしながら取り出して行列に配列してなる相関行列を作成する相関行列作成部、
該相関行列を用いて到来波数を推定するための推定行列を作成する推定行列作成部、
該推定行列をQR分解するQR分解部、
QR分解により得られた上三角行列因子の各行要素に基づいて到来電波の個数を決定する到来波数決定部、
を備えたことを特徴とする到来波数推定装置。 - 前記相関行列作成部は、q(0≦q≦M−1)を雑音の空間的相関の長さとするとき、前記アレー共分散行列を構成する所定の行あるいは列より対角要素および該対角要素に連鎖的に隣接するq個の要素を除き、(q+1)個の要素が除かれた(M−q−1)個の相関より、所定数の相関を順次1要素づつずらしながら前記所定数
個の相関を取り出して行列に配列してなる前記相関行列を作成する、
ことを特徴とする請求項3記載の到来波数推定装置。 - 前記相関行列作成部は、前記アレー共分散行列の先頭行及び最終行並びに先頭列及び最終列のうち1以上について前記相関行列を作成し、
前記推定行列作成部は、該作成された相関行列を並べることにより前記推定行列を作成する、
ことを特徴とする請求項3または4記載の到来波数推定装置。 - 前記相関行列作成部は、前記所定数
をM/2より小さい最大の整数とする、
ことを特徴とする請求項3または4記載の到来波数推定装置。 - 前記到来波数決定部は、前記上三角行列因子の階数を到来電波の個数とする、ことを特徴とする請求項3または4記載の到来波数推定方法。
- 前記到来波数決定部は、前記上三角行列因子の第i行目の各行列要素の絶対値の総和と第(i+1)行目の各行列要素の絶対値の総和の比を計算する計算部、
該比が最小のiを到来電波の個数と決定する決定部、
を有することを特徴とする請求項3または4記載の到来波数推定装置。 - 前記 QR分解部は、QR分解として列回転QR分解を行う、
ことを特徴とする請求項3記載の到来波数推定装置。 - 複数個(=M)のアンテナ素子を同一素子間隔で直線状に配列したアレーアンテナで到来電波を受信し、該電波の到来方向を推定する無線装置において、
到来電波の個数を推定する到来波数推定部、
アンテナ素子毎の受信データと到来電波の個数を用いて信号の到来方向を推定する到来方向推定部、
到来方向の推定値を用いて信号源方向にピークを有するように受信ビームフォーミング処理を行う受信ビームフォーマを備え、
前記到来波数推定部は、
前記アレーアンテナの各アンテナ素子の受信信号と他アンテナ素子の受信信号間の相関を算出する相関演算部、
各アンテナ素子の受信信号と他アンテナ素子の受信信号間の相関を行列に配列してなるアレー共分散行列を構成する所定の行あるいは列より対角要素を除き、該対角要素が除かれた(M−1)個の相関より所定数の相関を順次1要素づつずらしながら取り出して行列に配列してなる相関行列を作成する相関行列作成部、
該相関行列を用いて到来波数を推定するための推定行列を作成する推定行列作成部、
該推定行列をQR分解するQR分解部、
QR分解により得られた上三角行列因子の各行要素に基づいて到来電波の個数を決定する到来波数決定部、
を備えたことを特徴とする無線装置。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2005/012758 WO2007007390A1 (ja) | 2005-07-11 | 2005-07-11 | 到来波数推定方法、到来波数推定装置及び無線装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2007007390A1 JPWO2007007390A1 (ja) | 2009-01-29 |
JP4722132B2 true JP4722132B2 (ja) | 2011-07-13 |
Family
ID=37636798
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007524485A Expired - Fee Related JP4722132B2 (ja) | 2005-07-11 | 2005-07-11 | 到来波数推定方法、到来波数推定装置及び無線装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US7450067B2 (ja) |
EP (1) | EP1903348B1 (ja) |
JP (1) | JP4722132B2 (ja) |
WO (1) | WO2007007390A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9588221B2 (en) | 2012-05-16 | 2017-03-07 | Denso Corporation | Radar apparatus and method of processing incoming waves for the radar apparatus |
Families Citing this family (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3895228B2 (ja) * | 2002-05-07 | 2007-03-22 | 松下電器産業株式会社 | 無線通信装置および到来方向推定方法 |
EP1903348B1 (en) * | 2005-07-11 | 2015-12-30 | Fujitsu Ltd. | Number-of-arriving-waves estimating method, number-of-arriving-waves estimating device, and radio device |
JP4708224B2 (ja) * | 2006-03-06 | 2011-06-22 | 日本電信電話株式会社 | 無線信号分離方法および受信装置並びにそのプログラムと記録媒体 |
US7948436B2 (en) * | 2007-03-09 | 2011-05-24 | Telefonaktiebolaget L M Ericsson (Publ) | Array antenna arrangement |
US8446312B2 (en) * | 2007-12-25 | 2013-05-21 | Honda Elesys Co., Ltd. | Electronic scanning type radar device, estimation method of direction of reception wave, and program estimating direction of reception wave |
JP5311832B2 (ja) * | 2008-01-15 | 2013-10-09 | 三菱電機株式会社 | 信号処理装置 |
WO2009105418A1 (en) | 2008-02-20 | 2009-08-27 | Hobbit Wave | Beamforming devices and methods |
JP5600866B2 (ja) * | 2008-03-04 | 2014-10-08 | 富士通株式会社 | 探知測距装置および探知測距方法 |
JP5320792B2 (ja) | 2008-03-28 | 2013-10-23 | 富士通株式会社 | 到来方向推定装置、到来方向推定方法および到来方向推定プログラム |
KR101514816B1 (ko) | 2008-06-02 | 2015-04-23 | 삼성전자주식회사 | 알려진 신호를 이용하여 간섭을 인지하는 인지 무선 통신시스템 |
CN101436890B (zh) * | 2008-12-10 | 2012-11-21 | ***通信集团广东有限公司 | 一种天线阵列误差参数估计方法及其设备 |
JP5675285B2 (ja) * | 2010-11-10 | 2015-02-25 | 富士通テン株式会社 | レーダ装置 |
US8948718B2 (en) | 2012-03-07 | 2015-02-03 | Hobbit Wave, Inc. | Devices and methods using the Hermetic Transform |
US9154353B2 (en) | 2012-03-07 | 2015-10-06 | Hobbit Wave, Inc. | Devices and methods using the hermetic transform for transmitting and receiving signals using OFDM |
FR3011641B1 (fr) * | 2013-10-04 | 2015-12-11 | Thales Sa | Procede de goniometrie avec un systeme ayant moins de voies de reception que de capteurs |
WO2015061709A1 (en) | 2013-10-25 | 2015-04-30 | HOBBIT, Wave | Devices and methods employing hermetic transforms for encoding and decoding digital information in spread-spectrum communication systems |
WO2015105592A2 (en) | 2013-11-22 | 2015-07-16 | Hobbit Wave | Radar using hermetic transforms |
US11304661B2 (en) | 2014-10-23 | 2022-04-19 | VertoCOMM, Inc. | Enhanced imaging devices, and image construction methods and processes employing hermetic transforms |
US9871684B2 (en) | 2014-11-17 | 2018-01-16 | VertoCOMM, Inc. | Devices and methods for hermetic transform filters |
GB2539722B (en) | 2015-06-25 | 2021-10-13 | Airspan Ip Holdco Llc | Bearing calculation |
GB2539727B (en) | 2015-06-25 | 2021-05-12 | Airspan Ip Holdco Llc | A configurable antenna and method of operating such a configurable antenna |
GB2539734A (en) | 2015-06-25 | 2016-12-28 | Airspan Networks Inc | An antenna apparatus and method of performing spatial nulling within the antenna apparatus |
GB2539731B (en) | 2015-06-25 | 2021-08-04 | Airspan Ip Holdco Llc | Quality of service in wireless backhauls |
GB2539733A (en) | 2015-06-25 | 2016-12-28 | Airspan Networks Inc | An antenna apparatus and method of configuring a transmission beam for the antenna apparatus |
WO2016207603A1 (en) | 2015-06-25 | 2016-12-29 | Airspan Networks Inc. | Managing external interference in a wireless network |
GB2539732A (en) | 2015-06-25 | 2016-12-28 | Airspan Networks Inc | A configurable antenna and method of operating such a configurable antenna |
GB2539735A (en) * | 2015-06-25 | 2016-12-28 | Airspan Networks Inc | Sub-sampling antenna elements |
GB2539730B (en) | 2015-06-25 | 2021-04-07 | Airspan Ip Holdco Llc | Node role assignment in networks |
GB2539736A (en) | 2015-06-25 | 2016-12-28 | Airspan Networks Inc | Wireless network configuration using path loss determination between nodes |
US10324166B2 (en) * | 2015-09-28 | 2019-06-18 | Rockwell Collins, Inc. | Affordable combined pulsed/FMCW radar AESA |
US10305717B2 (en) | 2016-02-26 | 2019-05-28 | VertoCOMM, Inc. | Devices and methods using the hermetic transform for transmitting and receiving signals using multi-channel signaling |
JP6826873B2 (ja) * | 2016-12-07 | 2021-02-10 | 三菱重工業株式会社 | 測距・測角装置、測距・測角方法及び受信アレイセンサの設計方法 |
KR101956449B1 (ko) * | 2017-04-25 | 2019-03-08 | 서강대학교산학협력단 | 행렬의 거듭제곱기법을 이용한 도래방향각 추정방법, 상기 방법을 이용한 도래방향각 추정장치, 및 상기 방법을 이용한 도래방향각 추정시스템 |
CN112136057B (zh) * | 2018-05-07 | 2023-05-23 | 三菱电机株式会社 | 到达波数推定装置及到达波数到达方向推定装置 |
CN109901111A (zh) * | 2019-02-25 | 2019-06-18 | 西安电子科技大学 | 基于偏最小二乘回归的近场声源定位方法 |
CN110736958B (zh) * | 2019-09-18 | 2021-09-28 | 中国人民解放军国防科技大学 | 一种基于多策略矩阵重构的自适应信源数估计方法 |
CN110967664B (zh) * | 2019-11-28 | 2024-01-23 | 宁波大学 | 基于cold阵列增强四元数esprit的doa估计方法 |
CN112014792B (zh) * | 2020-08-28 | 2023-03-31 | 陕西理工大学 | 一种改进的二维互质阵列doa估计方法 |
DE102021107220B4 (de) * | 2021-03-23 | 2022-12-22 | Airbus Defence and Space GmbH | System und Verfahren zum Detektieren und Lokalisieren einer Signalquelle |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11281725A (ja) * | 1998-03-26 | 1999-10-15 | Nec Corp | 多重伝搬波パラメータ計測方法及び装置並びにプログラムを記録した機械読み取り可能な記録媒体 |
US6144711A (en) * | 1996-08-29 | 2000-11-07 | Cisco Systems, Inc. | Spatio-temporal processing for communication |
JP2002243826A (ja) * | 2000-12-12 | 2002-08-28 | Matsushita Electric Ind Co Ltd | 電波到来方向推定装置、電波到来方向推定方法及び指向性可変送受信装置 |
WO2005001504A1 (ja) * | 2003-06-25 | 2005-01-06 | Fujitsu Limited | 電波到来方向推定方法及び装置 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4750147A (en) * | 1985-11-06 | 1988-06-07 | Stanford University | Method for estimating signal source locations and signal parameters using an array of signal sensor pairs |
DE19858951C2 (de) * | 1998-12-10 | 2001-06-13 | Hertz Inst Heinrich | Verfahren und Anordnung zur Erzeugung vorgegebener Richtcharakteristiken |
EP1387180A1 (en) * | 2000-12-12 | 2004-02-04 | Matsushita Electric Industrial Co., Ltd. | Radio-wave arrival-direction estimating apparatus and directional variable transceiver |
WO2004104620A1 (ja) * | 2003-05-22 | 2004-12-02 | Fujitsu Limited | 固有値分解を利用しない信号到来方向推定手法および受信ビーム形成装置 |
JP4545150B2 (ja) * | 2004-08-12 | 2010-09-15 | 富士通株式会社 | 電波到来方向の適応推定追尾方法および装置 |
JP4805591B2 (ja) * | 2005-03-17 | 2011-11-02 | 富士通株式会社 | 電波到来方向の追尾方法及び電波到来方向追尾装置 |
JP4808984B2 (ja) * | 2005-03-31 | 2011-11-02 | 株式会社デンソーアイティーラボラトリ | 到来波の方向推定装置 |
EP1903348B1 (en) * | 2005-07-11 | 2015-12-30 | Fujitsu Ltd. | Number-of-arriving-waves estimating method, number-of-arriving-waves estimating device, and radio device |
-
2005
- 2005-07-11 EP EP05758310.6A patent/EP1903348B1/en not_active Expired - Fee Related
- 2005-07-11 WO PCT/JP2005/012758 patent/WO2007007390A1/ja not_active Application Discontinuation
- 2005-07-11 JP JP2007524485A patent/JP4722132B2/ja not_active Expired - Fee Related
-
2008
- 2008-01-08 US US11/970,645 patent/US7450067B2/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6144711A (en) * | 1996-08-29 | 2000-11-07 | Cisco Systems, Inc. | Spatio-temporal processing for communication |
JPH11281725A (ja) * | 1998-03-26 | 1999-10-15 | Nec Corp | 多重伝搬波パラメータ計測方法及び装置並びにプログラムを記録した機械読み取り可能な記録媒体 |
JP2002243826A (ja) * | 2000-12-12 | 2002-08-28 | Matsushita Electric Ind Co Ltd | 電波到来方向推定装置、電波到来方向推定方法及び指向性可変送受信装置 |
WO2005001504A1 (ja) * | 2003-06-25 | 2005-01-06 | Fujitsu Limited | 電波到来方向推定方法及び装置 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9588221B2 (en) | 2012-05-16 | 2017-03-07 | Denso Corporation | Radar apparatus and method of processing incoming waves for the radar apparatus |
Also Published As
Publication number | Publication date |
---|---|
WO2007007390A1 (ja) | 2007-01-18 |
EP1903348A4 (en) | 2014-06-25 |
JPWO2007007390A1 (ja) | 2009-01-29 |
EP1903348A1 (en) | 2008-03-26 |
EP1903348B1 (en) | 2015-12-30 |
US20080129596A1 (en) | 2008-06-05 |
US7450067B2 (en) | 2008-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4722132B2 (ja) | 到来波数推定方法、到来波数推定装置及び無線装置 | |
JP4805591B2 (ja) | 電波到来方向の追尾方法及び電波到来方向追尾装置 | |
JP4339801B2 (ja) | 固有値分解を利用しない信号到来方向推定手法および受信ビーム形成装置 | |
US7084812B2 (en) | Method and device for tracking the directions-of-arrival of radio waves | |
US7391370B2 (en) | Method and apparatus for adaptive direction estimation of radio waves | |
Varade et al. | Robust algorithms for DOA estimation and adaptive beamforming for smart antenna application | |
US8248304B2 (en) | Method for measuring incoming angles of coherent sources using space smoothing on any sensor network | |
CN107255796B (zh) | 一种非均匀噪声下窄带近场信号源定位方法 | |
JP2007159130A (ja) | 分散型アンテナ移動通信システムにおける上り受信方法及び装置 | |
JP2011014980A (ja) | 電波伝搬パラメータ推定装置、電波伝搬パラメータ推定方法 | |
Ramezanpour et al. | Two-stage beamforming for rejecting interferences using deep neural networks | |
RU2659608C1 (ru) | Способ синтеза многолучевой самофокусирующейся адаптивной антенной решетки с использованием параметрической модели корреляционной матрицы принимаемого сигнала | |
JP5408966B2 (ja) | 受信信号の到着時間(toa)を推定するための方法 | |
JP4660562B2 (ja) | 移動局方向推定方法及び装置 | |
Tayem et al. | Capon root-MUSIC-like direction of arrival estimation based on real data | |
Suleesathira | Direction of arrival identification using MUSIC method and NLMS beamforming | |
CN111337879B (zh) | 一种基于多簇的加权music声源直接定位方法 | |
Lin et al. | Efficient doa estimation method employing unitary improved polynomial rooting | |
Le et al. | Fast 1D Direction Finding Method In Multipath Environment | |
Aminu et al. | Performance evaluation of combined methods for the estimation of fading coefficients of uncorrelated signal sources in multipath propagation | |
Tayem et al. | A fast algorithm for direction of arrival estimation in multipath environments | |
Jovanović et al. | Performance analysis of Capon and Capon-like algorithm for smart antenna system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20101026 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101125 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20110405 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110405 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140415 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |