JP4808984B2 - 到来波の方向推定装置 - Google Patents
到来波の方向推定装置 Download PDFInfo
- Publication number
- JP4808984B2 JP4808984B2 JP2005100314A JP2005100314A JP4808984B2 JP 4808984 B2 JP4808984 B2 JP 4808984B2 JP 2005100314 A JP2005100314 A JP 2005100314A JP 2005100314 A JP2005100314 A JP 2005100314A JP 4808984 B2 JP4808984 B2 JP 4808984B2
- Authority
- JP
- Japan
- Prior art keywords
- calculation
- matrix
- householder
- conversion calculation
- conversion
- 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.)
- Active
Links
Images
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
-
- 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/80—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 ultrasonic, sonic or infrasonic waves
- G01S3/8006—Multi-channel systems specially adapted for direction-finding, i.e. having a single aerial system capable of giving simultaneous indications of the directions of different signals
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Radio Transmission System (AREA)
Description
複数のアンテナ(15)で受信した、複数の到来波の受信データから、その相関行列(R)を演算する相関行列演算手段(6)、
前記相関行列演算手段で演算された前記相関行列に対して、ハウスホルダ変換演算(図3のステップS5からステップS8など)を実施して、前記相関行列を、固有ベクトル(Q)と固有値(λ)とに分解処理する固有値分解演算手段(7)、
前記固有値分解演算手段に設けられ、前記相関行列について、固有値を求めるための行列(QT)と固有ベクトルを求めるための行列(A)を生成する行列生成手段(7,図3のステップS1、ステップS2、ステップS10など)、
前記固有値分解演算手段に設けられ、前記固有値を求める行列において、前記ハウスホルダ変換演算を行うたび毎に、前記ハウスホルダ変換演算の対象となる行列に対して、各列の成分の列ノルムを演算する列ノルム演算手段(7,ステップS3など)、
前記固有値分解演算手段に設けられ、前記ハウスホルダ変換演算を行うたび毎に、前記列ノルム演算手段により演算されたハウスホルダ変換演算の対象となる行列の列ノルムの最大値が、所定の第1の閾値(Th)以下であるか否かを判定する、列ノルム最大値判定手段(7,図3のステップS4など)、
前記固有値分解演算手段に設けられ、前記列ノルム最大値判定手段により、前記列ノルム演算手段により演算されたハウスホルダ変換演算の対象となる行列の列ノルムの最大値が、前記所定の第1の閾値以下でないものと判定された場合に限り、前記固有値を求めるための行列に対してハウスホルダ変換演算を、繰り返し実行するハウスホルダ変換演算実行手段(7,図3のステップS3乃至ステップS8)、
前記固有値分解演算手段に設けられ、前記列ノルム最大値判定手段により、前記列ノルム演算手段により演算されたハウスホルダ変換演算の対象となる行列の列ノルムの最大値が、前記所定の第1の閾値以下であるものと判定された場合に、前記ハウスホルダ変換演算実行手段によるハウスホルダ変換演算の実行を中断すると共に、前記ハウスホルダ変換演算実行手段によるハウスホルダ変換演算の実行回数を所定のメモリに格納する変換演算実行制御手段(7,図3のステップS4など)、
前記固有値分解演算手段に設けられ、前記変換演算実行制御手段によりハウスホルダ変換演算の実行が中断された前記固有値を求めるための行列に対して、前記ハウスホルダ変換演算実行手段及び変換演算実行制御手段によるハウスホルダ変換演算を、所定回数に達するまで更に実行させて、その結果から、前記相関行列の前記固有ベクトルと固有値を求める、収束演算実行手段(7.図3のステップS11、S12,S14など)、及び、
前記求められた固有値を構成する値と所定の第2の閾値を比較し、当該第2の閾値よりも大きな値については到来波の信号の固有値とし、小さなものは雑音の固有値と判定して、到来波の数を推定演算すると共に、該推定演算された到来波の数と前記メモリ手段に格納されたハウスホルダ変換演算の実行回数を比較し、その少ない方の値を到来波の数と判定する、信号ランク判定手段(9,10)、
から構成される。
前記変換演算実行制御手段は、前記列ノルムの最大値が負の場合に、前記ハウスホルダ変換演算の更なる実行を中断させる、ことを特徴として構成される。
複数のアンテナ(15)で受信した、複数の到来波の受信データから、その相関行列(R)を演算する相関行列演算手段(6)、
前記相関行列演算手段で演算された前記相関行列に対して、ハウスホルダ変換演算(図3のステップS5からステップS8など)を実施して、前記相関行列を、固有ベクトル(Q)と固有値(λ)とに分解処理する固有値分解演算手段(7)、
前記固有値分解演算手段に設けられ、前記相関行列について、固有値を求めるための行列(QT)と固有ベクトルを求めるための行列(A)を生成する行列生成手段(7,図3のステップS1、ステップS2、ステップS10など)、
前記固有値分解演算手段に設けられ、前記固有値を求める行列において、前記ハウスホルダ変換演算を行うたび毎に、前記ハウスホルダ変換演算の対象となる行列に対して、各列の成分の列ノルムを演算する列ノルム演算手段(7,ステップS3など)、
所定の第1の閾値を、前記受信データにおいて雑音のみが捕捉された際の、前記固有値を求めるための行列の列ノルムよりも大きく、前記受信データの到来波の信号についての、前記固有値を求めるための行列の列ノルムよりも小さな値となるように設定する閾値設定手段、
前記固有値分解演算手段に設けられ、前記ハウスホルダ変換演算を行うたび毎に、前記列ノルム演算手段により演算されたハウスホルダ変換演算の対象となる行列の列ノルムの最大値が、前記所定の第1の閾値(Th)以下であるか否かを判定する、列ノルム最大値判定手段(7,図3のステップS4など)、
前記固有値分解演算手段に設けられ、前記列ノルム最大値判定手段により、前記列ノルム演算手段により演算されたハウスホルダ変換演算の対象となる行列の列ノルムの最大値が、前記所定の第1の閾値以下でないものと判定された場合に限り、前記固有値を求めるための行列に対してハウスホルダ変換演算を、繰り返し実行するハウスホルダ変換演算実行手段(7,図3のステップS3乃至ステップS8)、
前記固有値分解演算手段に設けられ、前記列ノルム最大値判定手段により、前記列ノルム演算手段により演算されたハウスホルダ変換演算の対象となる行列の列ノルムの最大値が、前記所定の第1の閾値以下であるものと判定された場合に、前記ハウスホルダ変換演算実行手段によるハウスホルダ変換演算の実行を中断すると共に、前記ハウスホルダ変換演算実行手段によるハウスホルダ変換演算の実行回数を所定のメモリに格納する変換演算実行制御手段(7,図3のステップS4など)、
前記固有値分解演算手段に設けられ、前記変換演算実行制御手段によりハウスホルダ変換演算の実行が中断された前記固有値を求めるための行列に対して、前記ハウスホルダ変換演算実行手段及び変換演算実行制御手段によるハウスホルダ変換演算を、所定回数に達するまで更に実行させて、その結果から、前記相関行列の前記固有ベクトルと固有値を求める、収束演算実行手段(7.図3のステップS11、S12,S14など)、及び、
前記メモリ手段に格納されたハウスホルダ変換演算の実行回数を到来波の数と判定する、信号ランク判定手段、
を有することを特徴として構成される。
前記ハウスホルダ変換演算実行手段によるハウスホルダ変換演算を、前記演算回数設定手段に設定されたハウスホルダ変換演算の実行回数の最大数以内で行うことを特徴として構成される。
前記変換演算回数制御手段は、N回目の前記ハウスホルダ変換演算実行手段及び変換演算実行制御手段によるハウスホルダ変換演算の実行回数を、それ以前の回におけるハウスホルダ変換演算の実行回数以下に制御することを特徴として構成される。
サンプリング時刻mΔT1(ΔT1はサンプリング間隔、mは自然数)において、各アンテナ素子15-1から15-K(Kはアンテナ素子数)で得られる複素データをx1(m),x2(m),…,xK(m)とすると、信号受信部では、式(1)で表されるサイズKのデータベクトルx(m)を得る。上付きのTは転置を表す。
相関行列演算部では式(2)を用いて、時刻nΔT2(ΔT2は相関行列を算出する時間間隔)における相関行列Rxx(n)を求める。Rxxは(K×K)の半正定値あるいは正定値エルミート行列である。あとでも触れるが、半正定値エルミート行列の場合、固有値はすべて0以上の実数となる(正定値の場合、固有値は0より大きい実数となる)。上付きのHは複素共役転置を表す。SSNはスナップショット数で、相関行列演算に用いるデータベクトルの数を表す。βは忘却係数で0〜1の範囲の値を取る。Rxx(n-1)は時刻(n-1)ΔT2における相関行列である。
固有値分解演算は、ハウスホルダ変換演算を利用したQR分解を繰りかえし利用することで行う。1回の固有値分解演算は複数回のQR分解演算を含み、1回のQR分解演算は複数回のハウスホルダ変演算換を含む。ここでは、1回の固有値分解演算中に行うQR分解演算の回数をNQR[回]とする。また、1回のQR分解演算においては、最大(K-1)回(Kは素子数)のハウスホルダ変換演算が実行される。このハウスホルダ変換演算の実行回数は適宜なメモリ手段に格納される。
列ノルムの最大値>=しきい値 ハウスホルダ変換を実行
列ノルムの最大値< しきい値 ハウスホルダ変換を中止
通常、しきい値Thは列ノルムの最大値が負であることをチェックすることのできる値、即ち、「0」または「0」に近い値が予め設定されている。先にも述べたが,Ryy(あるいはRxx)が半正定値のエルミート行列であるので、得られる固有値は必ず0以上の実数となる(正定値の場合0より大きな実数)。固有値と列ノルムは相互に強く関連している。列ノルムが負または0に近い値となった場合には、桁落ちなどの演算誤差が生じていることが考えられ、以後、更なるハウスホルダ変換演算を行っても求められる固有値の信頼性が低下してしまい、演算負荷ばかりが増加してしまうことから、ハウスホルダ変換を中止しステップ図14-S500へ入る。
例としてK=4の場合のj(=1,2,3)回目のハウスホルダ変換終了後に得られる行列A(j)を示す。ここで×印は不特定の値が入る。Hは(K×K)のハウスホルダ行列である。A(0)は初期値であり、最初は相関行列Ryyが入る。
こうして固有値分解演算部7において、相関行列Ryy(あるいはRxx)の固有値λが求められたところで、当該求められた固有値λを、図1の信号ランク推定演算部9に出力し、信号ランク推定演算部9は、固有値λから、到来波数(信号ランク)を推定する。推定法としては、AICやMDL、しきい値法などがある。しきい値法では、例えば予め設定さえた所定しきい値より大きな固有値については到来波の信号の固有値とし、小さなものは雑音の固有値と判定して、信号ランクLt、即ち到来波の数Ltを推定する。
信号ランク判定部10は、信号ランク推定演算部9で推定された信号ランク、即ち到電波の数Ltと、固有値分解演算部7で実行されたハウスホルダ変換実行回数Nh(i)(i=1,2,..,NQR)のうち、一番小さい値を信号ランク、即ち到来波の数Lとする演算を行い、求められた到来波(信号)の数Lを、スペクトル演算部11に対して出力する。
こうして、信号ランク推定部10により、ハウスホルダ変換演算回数Nhが考慮された形で、信号ランクLが判定され、その値Lが、図1に示す、スペクトル演算部11に出力されると、スペクトル演算部11は、公知の手法により、図2に示すようにMUSICスペクトラムSPEを演算して求め、当該求められたMUSICスペクトラムSPEを、到来波方位演算部12に出力する。
到来波方位演算部12は、当該MUSICスペクトラムSPEから、公知の手法で、到来波の到来する方位θ1、θ2を演算することが出来るが、MUSICスペクトラムSPEは、既に述べたように、信号の数Lが、相関行列Ryy(あるいはRxx)の固有値分解演算に際したハウスホルダ変換演算回数Nhに基づいて補正され、信頼性の低いデータが排除されているので、実態のない所謂「偽ピーク」を排除した形で、正確な到来波の方位を算出することができる。ただし「偽ピーク」を完全に排除できるわけではないので、より正確な推定には電力推定演算などを利用する必要がある(関連参考特許文献[5])。
文献[3]Gene H Golub et al,‘MATRIX COMPUTATION’(3rd edition),Jhon Hopkins Univ.Press,ISBN0-8018-5414-8
参考特許文献[4] 特願平10−152629
参考特許文献[5] 特願平10−290766
ここで、列ノルムccの最大値を示す列を抽出し、ind_kに代入して、その値をRAM17に格納する。図4(4)の場合、2列目の列ノルムccが最大値:2.9260を取るので、ind_k=2を代入する。
6……相関行列演算部
7……固有値分解演算部
9……信号ランク推定演算部
10……信号ランク判定部
15……アンテナ
Th……第1の閾値
Claims (7)
- 複数のアンテナで受信した複数の到来波の受信データから、前記到来波の到来方向を推定する、到来波の方向推定装置において、
前記複数のアンテナで受信した、複数の到来波の受信データから、その相関行列を演算する相関行列演算手段、
前記相関行列演算手段で演算された前記相関行列に対して、ハウスホルダ変換演算を実施して、前記相関行列を、固有ベクトルと固有値とに分解処理する固有値分解演算手段、
前記固有値分解演算手段に設けられ、前記相関行列について、固有値を求めるための行列と固有ベクトルを求めるための行列を生成する行列生成手段、
前記固有値分解演算手段に設けられ、前記固有値を求めるための行列において、前記ハウスホルダ変換演算を行うたび毎に、前記ハウスホルダ変換演算の対象となる行列に対して、各列の成分の列ノルムを演算する列ノルム演算手段、
前記固有値分解演算手段に設けられ、前記ハウスホルダ変換演算を行うたび毎に、前記列ノルム演算手段により演算されたハウスホルダ変換演算の対象となる行列の列ノルムの最大値が、所定の第1の閾値以下であるか否かを判定する、列ノルム最大値判定手段、
前記固有値分解演算手段に設けられ、前記列ノルム最大値判定手段により、前記列ノルム演算手段により演算されたハウスホルダ変換演算の対象となる行列の列ノルムの最大値が、前記所定の第1の閾値以下でないものと判定された場合に限り、前記固有値を求めるための行列に対してハウスホルダ変換演算を、繰り返し実行するハウスホルダ変換演算実行手段、
前記固有値分解演算手段に設けられ、前記列ノルム最大値判定手段により、前記列ノルム演算手段により演算されたハウスホルダ変換演算の対象となる行列の列ノルムの最大値が、前記所定の第1の閾値以下であるものと判定された場合に、前記ハウスホルダ変換演算実行手段によるハウスホルダ変換演算の実行を中断すると共に、前記ハウスホルダ変換演算実行手段によるハウスホルダ変換演算の実行回数を所定のメモリに格納する変換演算実行制御手段、
前記固有値分解演算手段に設けられ、前記変換演算実行制御手段によりハウスホルダ変換演算の実行が中断された前記固有値を求めるための行列に対して、前記ハウスホルダ変換演算実行手段及び変換演算実行制御手段によるハウスホルダ変換演算を、所定回数に達するまで更に実行させて、その結果から、前記相関行列の前記固有ベクトルと固有値を求める、収束演算実行手段、及び、
前記求められた固有値を構成する値と所定の第2の閾値を比較し、当該第2の閾値よりも大きな値については到来波の信号の固有値とし、小さなものは雑音の固有値と判定して、到来波の数を推定演算すると共に、該推定演算された到来波の数と前記メモリ手段に格納されたハウスホルダ変換演算の実行回数を比較し、その少ない方の値を到来波の数と判定する、信号ランク判定手段、
を有する、到来波の方向推定装置。 - 前記第1の閾値として、0又は0に近い値を格納するメモリ手段を有し、
前記変換演算実行制御手段は、前記列ノルムの最大値が負の場合に、前記ハウスホルダ変換演算の更なる実行を中断させる、ことを特徴とする、請求項1記載の到来波の方向推定装置。 - 複数のアンテナで受信した複数の到来波の受信データから、前記到来波の到来方向を推定する、到来波の方向推定装置において、
前記複数のアンテナで受信した、複数の到来波の受信データから、その相関行列を演算する相関行列演算手段、
前記相関行列演算手段で演算された前記相関行列に対して、ハウスホルダ変換演算を実施して、前記相関行列を、固有ベクトルと固有値とに分解処理する固有値分解演算手段、
前記固有値分解演算手段に設けられ、前記相関行列について、固有値を求めるための行列と固有ベクトルを求めるための行列を生成する行列生成手段、
前記固有値分解演算手段に設けられ、前記固有値を求めるための行列において、前記ハウスホルダ変換演算を行うたび毎に、前記ハウスホルダ変換演算の対象となる行列に対して、各列の成分の列ノルムを演算する列ノルム演算手段、
所定の第1の閾値を、前記受信データにおいて雑音のみが捕捉された際の、前記固有値を求めるための行列の列ノルムよりも大きく、前記受信データの到来波の信号についての、前記固有値を求めるための行列の列ノルムよりも小さな値となるように設定する閾値設定手段、
前記固有値分解演算手段に設けられ、前記ハウスホルダ変換演算を行うたび毎に、前記列ノルム演算手段により演算されたハウスホルダ変換演算の対象となる行列の列ノルムの最大値が、前記所定の第1の閾値以下であるか否かを判定する、列ノルム最大値判定手段、
前記固有値分解演算手段に設けられ、前記列ノルム最大値判定手段により、前記列ノルム演算手段により演算されたハウスホルダ変換演算の対象となる行列の列ノルムの最大値が、前記所定の第1の閾値以下でないものと判定された場合に限り、前記固有値を求めるための行列に対してハウスホルダ変換演算を、繰り返し実行するハウスホルダ変換演算実行手段、
前記固有値分解演算手段に設けられ、前記列ノルム最大値判定手段により、前記列ノルム演算手段により演算されたハウスホルダ変換演算の対象となる行列の列ノルムの最大値が、前記所定の第1の閾値以下であるものと判定された場合に、前記ハウスホルダ変換演算実行手段によるハウスホルダ変換演算の実行を中断すると共に、前記ハウスホルダ変換演算実行手段によるハウスホルダ変換演算の実行回数を所定のメモリに格納する変換演算実行制御手段、
前記固有値分解演算手段に設けられ、前記変換演算実行制御手段によりハウスホルダ変換演算の実行が中断された前記固有値を求めるための行列に対して、前記ハウスホルダ変換演算実行手段及び変換演算実行制御手段によるハウスホルダ変換演算を、所定回数に達するまで更に実行させて、その結果から、前記相関行列の前記固有ベクトルと固有値を求める、収束演算実行手段、及び、
前記メモリ手段に格納されたハウスホルダ変換演算の実行回数を到来波の数と判定する、信号ランク判定手段、
を有することを特徴とする、到来波の方向推定装置。 - 前記閾値設定手段は、前記所定の第1の閾値を、周囲の環境や観測距離に対応させて変動的に設定することを特徴とする、請求項3記載の到来波の方向推定装置。
- 前記ハウスホルダ変換演算実行手段によるハウスホルダ変換演算の実行回数の最大数を設定する演算回数設定手段を有し、
前記ハウスホルダ変換演算実行手段によるハウスホルダ変換演算を、前記演算回数設定手段に設定されたハウスホルダ変換演算の実行回数の最大数以内で行うことを特徴とする、請求項1又は3記載の到来波の方向推定装置。 - 前記演算回数設定手段は、前記ハウスホルダ変換演算の実行回数の最大数を、周囲の環境や観測距離に対応させて変動的に設定することを特徴とする、請求項5記載の到来波の方向推定装置。
- 前記収束演算実行手段による、前記ハウスホルダ変換演算実行手段及び変換演算実行制御手段による前記ハウスホルダ変換演算の、各回の実行回数を制御する変換演算回数制御手段を設け、
前記変換演算回数制御手段は、N回目の前記ハウスホルダ変換演算実行手段及び変換演算実行制御手段によるハウスホルダ変換演算の実行回数を、それ以前の回におけるハウスホルダ変換演算の実行回数以下に制御することを特徴とする、請求項1または3記載の到来波の方向推定装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005100314A JP4808984B2 (ja) | 2005-03-31 | 2005-03-31 | 到来波の方向推定装置 |
DE102006013448.6A DE102006013448B4 (de) | 2005-03-31 | 2006-03-20 | Vorrichtung zur Schätzung der Richtung einer ankommenden Welle |
US11/390,866 US7372404B2 (en) | 2005-03-31 | 2006-03-27 | Apparatus for estimating direction of arrival wave |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005100314A JP4808984B2 (ja) | 2005-03-31 | 2005-03-31 | 到来波の方向推定装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2006284180A JP2006284180A (ja) | 2006-10-19 |
JP4808984B2 true JP4808984B2 (ja) | 2011-11-02 |
Family
ID=37071873
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2005100314A Active JP4808984B2 (ja) | 2005-03-31 | 2005-03-31 | 到来波の方向推定装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US7372404B2 (ja) |
JP (1) | JP4808984B2 (ja) |
DE (1) | DE102006013448B4 (ja) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
US7839835B2 (en) * | 2006-08-22 | 2010-11-23 | Nec Laboratories America, Inc. | Quantized precoding over a set of parallel channels |
JP4531738B2 (ja) * | 2006-11-21 | 2010-08-25 | 日本電信電話株式会社 | 行列における数値分解方法 |
KR100932789B1 (ko) * | 2007-12-15 | 2009-12-21 | 한국전자통신연구원 | 다중입력 다중출력 시스템에서 qr 분해 장치 및 그 방법 |
JP5600866B2 (ja) * | 2008-03-04 | 2014-10-08 | 富士通株式会社 | 探知測距装置および探知測距方法 |
JP5519132B2 (ja) * | 2008-07-28 | 2014-06-11 | 株式会社デンソー | レーダ装置 |
JP5102165B2 (ja) * | 2008-09-22 | 2012-12-19 | 株式会社デンソー | レーダ装置 |
JP5511479B2 (ja) * | 2010-04-16 | 2014-06-04 | 三菱電機株式会社 | 測定装置 |
US8224627B2 (en) | 2010-04-28 | 2012-07-17 | Mitov Iliya P | Technique for determination of the signal subspace dimension |
US20110286306A1 (en) * | 2010-05-21 | 2011-11-24 | Leo Eisner | Determining origin and mechanism of microseismic events in the earth's subsurface by deviatoric moment inversion |
US9294310B2 (en) * | 2010-10-08 | 2016-03-22 | Blackberry Limited | Method and apparatus for LTE channel state information estimation |
KR20120071851A (ko) * | 2010-12-23 | 2012-07-03 | 한국전자통신연구원 | 통신 시스템에서 신호 도래 방향 추정 장치 및 방법 |
JP5628732B2 (ja) * | 2011-04-04 | 2014-11-19 | 富士通テン株式会社 | レーダ装置用の演算装置、レーダ装置、レーダ装置用の演算方法およびプログラム |
JP6195490B2 (ja) * | 2013-08-23 | 2017-09-13 | 日本放送協会 | Mimo受信装置 |
NL2014562A (en) * | 2014-04-04 | 2015-10-13 | Asml Netherlands Bv | Control system, positioning system, lithographic apparatus, control method, device manufacturing method and control program. |
CN104297594A (zh) * | 2014-10-13 | 2015-01-21 | 电子科技大学 | 一种智能家居的家电设备运行状态监测方法 |
RU2604871C2 (ru) * | 2015-04-15 | 2016-12-20 | федеральное государственное автономное образовательное учреждение высшего профессионального образования "Южный федеральный университет" (Южный федеральный университет") | Способ определения местоположения объекта навигации |
US11067659B2 (en) | 2016-12-15 | 2021-07-20 | Bae Systems Information And Electronic Systems Integration Inc. | System and method for rank estimation of electromagnetic emitters |
US10473746B1 (en) * | 2016-12-15 | 2019-11-12 | Bae Systems Information And Electronic Systems Integration Inc. | System and method for rank estimation of electromagnetic emitters |
US11255955B2 (en) * | 2018-12-28 | 2022-02-22 | Panasonic Intellectual Property Management Co., Ltd. | Estimation method, estimation device, and recording medium |
CN111381228A (zh) * | 2018-12-28 | 2020-07-07 | 松下知识产权经营株式会社 | 推测方法以及推测装置 |
US10986509B1 (en) * | 2020-05-14 | 2021-04-20 | At&T Intellectual Property I, L.P. | Placement of antennas for fifth generation (5G) or other next generation networks |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09212489A (ja) * | 1996-01-31 | 1997-08-15 | Fujitsu Ltd | 対称行列の固有値問題を解く並列処理装置および方法 |
JP3081522B2 (ja) * | 1996-02-14 | 2000-08-28 | 株式会社エイ・ティ・アール光電波通信研究所 | 受信信号処理装置 |
JPH11344517A (ja) | 1998-06-02 | 1999-12-14 | Nec Corp | 電波環境分析装置 |
JP2000121716A (ja) * | 1998-10-13 | 2000-04-28 | Anritsu Corp | 電波到来方向推定装置 |
JP2005176020A (ja) * | 2003-12-12 | 2005-06-30 | Rikogaku Shinkokai | 復号方法および復号装置 |
-
2005
- 2005-03-31 JP JP2005100314A patent/JP4808984B2/ja active Active
-
2006
- 2006-03-20 DE DE102006013448.6A patent/DE102006013448B4/de active Active
- 2006-03-27 US US11/390,866 patent/US7372404B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2006284180A (ja) | 2006-10-19 |
DE102006013448B4 (de) | 2022-06-09 |
DE102006013448A1 (de) | 2007-01-18 |
US7372404B2 (en) | 2008-05-13 |
US20060224655A1 (en) | 2006-10-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4808984B2 (ja) | 到来波の方向推定装置 | |
Bazzi et al. | Detection of the number of superimposed signals using modified MDL criterion: A random matrix approach | |
Athley | Threshold region performance of maximum likelihood direction of arrival estimators | |
JP4138825B2 (ja) | ウェイト算出方法、ウェイト算出装置、アダプティブアレーアンテナ、及びレーダ装置 | |
US6567034B1 (en) | Digital beamforming radar system and method with super-resolution multiple jammer location | |
US7576682B1 (en) | Method and system for radar target detection and angle estimation in the presence of jamming | |
EP1425607A1 (en) | Radar system and method including superresolution raid counting | |
CN107870315B (zh) | 一种利用迭代相位补偿技术估计任意阵列波达方向方法 | |
WO2006067869A1 (ja) | 到来方向推定装置及びプログラム | |
JP2008096137A (ja) | レーダ装置及び測角装置 | |
CN109471063B (zh) | 基于延迟快拍的均匀线列阵高分辨波达方向估计方法 | |
CN109901103B (zh) | 基于非正交波形的mimo雷达doa估算方法及设备 | |
CN116881385B (zh) | 轨迹平滑方法、装置、电子设备及可读存储介质 | |
CN117092585B (zh) | 单比特量化DoA估计方法、***和智能终端 | |
Yan et al. | Computationally efficient direction of arrival estimation with unknown number of signals | |
WO2010066306A1 (en) | Apparatus and method for constructing a sensor array used for direction of arrival (doa) estimation | |
CN106844886B (zh) | 基于主分量分析的目标波达方向获取方法 | |
JP5152949B2 (ja) | ウェイト算出方法、ウェイト算出装置、アダプティブアレーアンテナ、及びレーダ装置 | |
CN115421098A (zh) | 嵌套面阵下降维求根music的二维doa估计方法 | |
JP2008032437A (ja) | ウェイト算出方法、ウェイト算出装置、アダプティブアレーアンテナ、及びレーダ装置 | |
CN108051773A (zh) | 基于盖式圆盘准则估计信源数目的epuma方法 | |
CN109683128B (zh) | 冲击噪声环境下的单快拍测向方法 | |
Roemer et al. | Higher order SVD based subspace estimation to improve multi-dimensional parameter estimation algorithms | |
CN113281698A (zh) | 一种嵌套阵中基于级联的非高斯信源测向方法 | |
CN111880167A (zh) | 一种基于先随机后优化的波达方向估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070614 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100402 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20101026 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101214 |
|
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: 20110726 |
|
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: 20110818 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140826 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4808984 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140826 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |