[go: up one dir, main page]

JPH11142497A - Wave bearing detector - Google Patents

Wave bearing detector

Info

Publication number
JPH11142497A
JPH11142497A JP9305547A JP30554797A JPH11142497A JP H11142497 A JPH11142497 A JP H11142497A JP 9305547 A JP9305547 A JP 9305547A JP 30554797 A JP30554797 A JP 30554797A JP H11142497 A JPH11142497 A JP H11142497A
Authority
JP
Japan
Prior art keywords
incident
signal processing
processing device
array antenna
wave number
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.)
Granted
Application number
JP9305547A
Other languages
Japanese (ja)
Other versions
JP3647621B2 (en
Inventor
Kazufumi Hirata
和史 平田
Tetsuo Kirimoto
哲郎 桐本
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP30554797A priority Critical patent/JP3647621B2/en
Publication of JPH11142497A publication Critical patent/JPH11142497A/en
Application granted granted Critical
Publication of JP3647621B2 publication Critical patent/JP3647621B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a wave bearing detector capable of increasing accuracy by reducing the amount of calculation when the number of incident waves is small. SOLUTION: This wave bearing detector comprises a alley antenna 1; an incident wave estimating means 10 which estimates the number of incident waves; an algorithm selector means 20 which selects an algorithm to obtain a wave input direction according to the number of incident waves; a signal processing device 30 which obtains the input direction from the correlation values of received signals when the number of incident waves is K1 or less; a signal processing device 40 which calculates the correlated matrixes and inter-correlated vectors of the received signals when the number of incident waves is K1 or larger and K2 or less, calculates a load value from the calculated values, and solves an equation with the load value as a factor so as to obtain an input direction; and a signal processing device 50 which calculates the load value with an adaptive filter from the received signals when the number of incident waves is K2 or larger, calculates an antenna pattern from the load value, and detects the null point of the antenna pattern so as to obtain the input direction. Thus, when the number of incident waves is less, accuracy can be increased and the amount of calculation can be decreased.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】この発明は、アレーアンテナ
の受信信号から、複数の入射波の到来方向を測定する電
波方位探知装置に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a radio azimuth detecting device for measuring arrival directions of a plurality of incident waves from a signal received by an array antenna.

【0002】[0002]

【従来の技術】従来の電波方位探知装置について図5を
参照しながら説明する。図5は、例えばR.O.Schmid
t、”Multiple Emitter Location and Signal Paramete
r Estimation”、IEEE Trans.Antennas Propag.、Vol.A
P-34、No.3、PP.276-280(March 1986)に示された従来
のミュージック(MUSIC:multiple signal classi
fication)アルゴリズムを用いた電波方位探知装置を示
す図である。
2. Description of the Related Art A conventional radio direction finder will be described with reference to FIG. FIG. 5, for example, ROSchmid
t, “Multiple Emitter Location and Signal Paramete
r Estimation ”, IEEE Trans.Antennas Propag., Vol.A
Conventional music (MUSIC: multiple signal classi) shown in P-34, No. 3, and PP.276-280 (March 1986)
FIG. 3 is a diagram illustrating a radio wave direction detecting device using an algorithm.

【0003】図5において、1はアレーアンテナであり
その受信信号をMUSIC処理器2に入力する。
[0005] In FIG. 5, reference numeral 1 denotes an array antenna which inputs a received signal to a MUSIC processor 2.

【0004】また、同図において、3〜7はMUSIC
処理器2内の処理フローを示している。3はアレーアン
テナ1の出力を入力して、各アレー素子間の相関を求め
て「相関行列R」を計算するステップ、4は相関行列R
の「固有値λ」、及び「固有ベクトルe」を計算するス
テップ、5はステップ4で求められた固有値のうちノイ
ズに対応した固有値に対応する固有ベクトルを算出する
ステップ、6は測定する各方位に対する各アレーアンテ
ナ1の受信位相差や、利得の情報を持つ「ステアリング
ベクトル」とステップ5で得られる固有ベクトルの内積
の絶対値の2乗の逆数を「方位評価関数P(θ)」とし
て算出するステップ、7はステップ6で得られる方位評
価関数P(θ)に現れるピークから入射波の到来方向θ
を推定するステップである。
Further, in FIG.
2 shows a processing flow in the processor 2. 3 is a step of inputting the output of the array antenna 1 and calculating a correlation matrix R by obtaining a correlation between the array elements, and 4 is a correlation matrix R
Calculating the eigenvalue λ and eigenvector e of step 5; calculating the eigenvector corresponding to the eigenvalue corresponding to the noise among the eigenvalues obtained in step 4; A step of calculating the reciprocal of the square of the absolute value of the inner value of the inner product of the “steering vector” having the reception phase difference and gain of the antenna 1 and the eigenvector obtained in step 5 as the “azimuth evaluation function P (θ)”, 7 Is the arrival direction θ of the incident wave from the peak appearing in the azimuth evaluation function P (θ) obtained in step 6.
Is the step of estimating.

【0005】つぎに、前述した従来のMUSICアルゴ
リズムを用いた電波方位探知装置の動作について説明す
る。
Next, the operation of the radio wave direction detecting device using the conventional MUSIC algorithm will be described.

【0006】MUSICアルゴリズムは、M個の素子ア
ンテナを有するアレーアンテナ1に、Mより少ないK波
の信号が入射する場合に、アレーアンテナ1の受信信号
の相関行列Rのノイズに対応した固有値に対応する固有
ベクトルと電波の入射方向に対応するステアリングベク
トルが直交する性質を利用した方探アルゴリズムであ
る。
[0006] The MUSIC algorithm corresponds to an eigenvalue corresponding to noise of the correlation matrix R of a received signal of the array antenna 1 when a signal of K waves smaller than M is incident on the array antenna 1 having M element antennas. This is a direction finding algorithm that utilizes the property that the eigenvector and the steering vector corresponding to the incident direction of the radio wave are orthogonal.

【0007】アレーアンテナ1を構成する素子アンテナ
は、無指向性のアンテナとし、m番目のアレー素子アン
テナの受信信号を「xm」、ノイズを「nm」、k番目の
到来電波「sk」、到来方向を「θk」とする。このと
き、アレーアンテナ1の「受信信号列ベクトル」は、次
の式(1)のように表される。
[0007] element antennas constituting an array antenna 1, an omnidirectional antenna, "x m" received signal of the m-th array element antenna, the noise "n m", k-th arrival wave "s k ”, And the arrival direction is“ θ k ”. At this time, the “received signal column vector” of the array antenna 1 is represented by the following equation (1).

【0008】[0008]

【数1】 (Equation 1)

【0009】上記の式(1)のa(θk)は、次の式
(2)で示される到来方向θkに対応したステアリング
ベクトルである。
A (θ k ) in the above equation (1) is a steering vector corresponding to the arrival direction θ k shown in the following equation (2).

【0010】[0010]

【数2】 (Equation 2)

【0011】上記の式(3)において、gm(θk)は素
子アンテナmの到来方向θkの「複素ゲイン」、exp
(j(2π/λ)pk・dm)は受信位相差を表す複素数
である。また、pkはθk方向の「単位ベクトル」、dm
は素子アンテナmの「位置ベクトル」である。
In the above equation (3), g mk ) is the “complex gain” of the arrival direction θ k of the element antenna m, exp
(J (2π / λ) p k · d m) is a complex number representing a reception phase difference. Also, pk is a “unit vector” in the θ k direction, d m
Is the “position vector” of the element antenna m.

【0012】上記の式(1)において、左辺の受信信号
ベクトルを「X」、ステアリングベクトルで構成される
行列を「A」、到来信号ベクトルを「S」、ノイズベク
トルを「N」とすると、次の式(4)のように書き換え
られる。
In the above equation (1), if the received signal vector on the left side is “X”, the matrix composed of the steering vector is “A”, the incoming signal vector is “S”, and the noise vector is “N”, It can be rewritten as the following equation (4).

【0013】[0013]

【数3】 (Equation 3)

【0014】相関行列Rは、上記の式(4)から、XX
Tとして求められる。以上の処理が図5のステップ3に
おいて行われる処理である。アレーアンテナ数M>到来
信号数Kの時、相関行列Rの最小固有値は受信機雑音に
対応し、その固有ベクトルEminと電波到来方向θkに対
応したステアリングベクトルa(θk)とは直交する性
質がある。そこで、MUSICアルゴリズムでは次の式
(5)で示されるようなP(θ)を方位評価関数として
用いる。
From the above equation (4), the correlation matrix R is given by XX
Required as T. The above processing is the processing performed in step 3 of FIG. When the number of array antennas M> the number of arriving signals K, the smallest eigenvalue of the correlation matrix R corresponds to the receiver noise, and the eigenvector E min and the steering vector a (θ k ) corresponding to the radio wave arrival direction θ k are orthogonal. There is nature. Therefore, in the MUSIC algorithm, P (θ) as shown in the following equation (5) is used as the azimuth evaluation function.

【0015】[0015]

【数4】 (Equation 4)

【0016】なお、相関行列Rの固有値及び固有ベクト
ルは図5のステップ4により求められ、図5のステップ
5において、求められたM個の固有値のうち最小な固有
値λminに対応する固有ベクトルEminが選定され、式
(5)に示す処理は、図5のステップ6において行われ
る。式(5)に示した方位評価関数より、固有ベクトル
minと電波到来方向θに対応したステアリングベクト
ルa(θ)が直交する時に、方位評価関数P(θ)は極
大となり、この方向θが電波の到来方向として求められ
る。図5のステップ7では、この方位評価関数に現れる
ピークから入射波の到来方向を推定する。
The eigenvalues and eigenvectors of the correlation matrix R are obtained in step 4 of FIG. 5. In step 5 of FIG. 5, the eigenvector E min corresponding to the smallest eigenvalue λ min among the obtained M eigenvalues is calculated. The process selected and shown in equation (5) is performed in step 6 of FIG. According to the azimuth evaluation function shown in equation (5), when the eigenvector E min and the steering vector a (θ) corresponding to the radio wave arrival direction θ are orthogonal, the azimuth evaluation function P (θ) becomes maximum, and this direction θ becomes equal to the radio wave. Is required as the direction of arrival. In step 7 in FIG. 5, the arrival direction of the incident wave is estimated from the peak appearing in the azimuth evaluation function.

【0017】MUSICアルゴリズムでは、全方向θに
対し、方位評価関数P(θ)を求めることが理想である
が、現実には離散的なθ列に対してしか、方位評価関数
P(θ)を計算することができない。このθ列のポイン
ト数をiとする。このとき、(測定する全方位角)/i
の分解能で方位評価関数を計算することができ、ポイン
ト数iが大きい程、より高分解能な方探が可能となる
が、それだけ演算量が増大する欠点がある。
In the MUSIC algorithm, it is ideal to find the azimuth evaluation function P (θ) for all directions θ, but in reality, the azimuth evaluation function P (θ) is obtained only for discrete θ columns. Cannot calculate. The number of points in the θ column is defined as i. At this time, (all azimuths to be measured) / i
The azimuth evaluation function can be calculated with the following resolution. The larger the number i of points, the higher the direction of the search. However, there is a disadvantage that the amount of calculation increases accordingly.

【0018】[0018]

【発明が解決しようとする課題】上述したような従来の
MUSICアルゴリズムを用いた電波方位探知装置で
は、入射波数が少ない場合でも演算量が多い処理を行う
必要があり、精度を上げるためにはさらに測定時間がか
かるという問題点があった。
In the conventional radio azimuth detecting device using the MUSIC algorithm as described above, it is necessary to perform a large amount of processing even when the number of incident waves is small. There was a problem that measurement time was required.

【0019】この発明は、前述した問題点を解決するた
めになされたもので、入射波数が少ない場合には演算量
を少なくでき、精度を高めることができる電波方位探知
装置を得ることを目的とする。
SUMMARY OF THE INVENTION The present invention has been made to solve the above-described problems, and has as its object to provide a radio azimuth detecting device capable of reducing the amount of calculation and increasing the accuracy when the number of incident waves is small. I do.

【0020】[0020]

【課題を解決するための手段】この発明に係る電波方位
探知装置は、電波を受信するアレーアンテナと、前記ア
レーアンテナの受信信号から入射波数を推定する入射波
数推定手段と、前記入射波数推定手段により推定された
入射波数に応じて電波の入射方向を求める信号処理装置
を切り替える信号処理装置切替手段と、前記入射波数が
正の整数である第1の定数以下である場合は、前記アレ
ーアンテナの受信信号の相関値に基づいて電波の入射方
向を求める第1の信号処理装置と、前記入射波数が前記
第1の定数より大きく、正の整数であって前記第1の定
数より大きい値の第2の定数より小さい場合は、前記ア
レーアンテナの受信信号の相関行列及び相互相関ベクト
ルを計算し、前記相関行列及び相互相関ベクトルから荷
重値を計算し、前記荷重値を係数とする所定の方程式を
解いて電波の入射方向を求める第2の信号処理装置と、
前記入射波数が前記第2の定数より大きい場合には、前
記アレーアンテナの受信信号から適応フィルタにより荷
重値を計算し、前記荷重値からアンテナパターンを計算
し、前記アンテナパターンのヌル点を検出して電波の入
射方向を求める第3の信号処理装置とを備えたものであ
る。
A radio azimuth detecting device according to the present invention comprises an array antenna for receiving radio waves, an incident wave number estimating means for estimating an incident wave number from a signal received by the array antenna, and an incident wave number estimating means. A signal processing device switching means for switching a signal processing device that determines an incident direction of a radio wave in accordance with the number of incident waves estimated by the above, and when the number of incident waves is equal to or less than a first constant that is a positive integer, the array antenna A first signal processing device for determining an incident direction of a radio wave based on a correlation value of a received signal; and a first signal processing device in which the number of incident waves is larger than the first constant and is a positive integer and larger than the first constant. If it is smaller than the constant of 2, the correlation matrix and the cross-correlation vector of the received signal of the array antenna are calculated, and the weight value is calculated from the correlation matrix and the cross-correlation vector. A second signal processing device for determining the direction of incidence of radio waves by solving predetermined equations for the coefficient a load value,
When the incident wave number is larger than the second constant, a weight value is calculated by an adaptive filter from a reception signal of the array antenna, an antenna pattern is calculated from the weight value, and a null point of the antenna pattern is detected. And a third signal processing device for determining the incident direction of the radio wave.

【0021】また、この発明に係る電波方位探知装置
は、前記第1の信号処理装置が、前記アレーアンテナの
2つの素子アンテナの受信信号の相関値を計算してその
結果から2つの素子アンテナの受信位相差を計算して入
射方向を求め、前記第2の信号処理装置が、前記アレー
アンテナの一部の受信信号の相関処理を行って部分的な
相関行列と相互相関ベクトルを計算し、前記部分的な相
関行列と相互相関ベクトルから信号の平均電力を最小化
する荷重値を計算し、前記荷重値を係数とする前記入射
波数と同じ次数の方程式を解いて入射方向を求め、前記
第3の信号処理装置が、信号の平均電力を最小化する荷
重値を適応フィルタにより漸化的に計算し、前記アレー
アンテナの各素子に荷重値を乗じた場合のアンテナパタ
ーンを計算し、前記アンテナパターンのヌル点を検出し
て電波の入射方向を求めるものである。
Further, in the radio azimuth detecting device according to the present invention, the first signal processing device calculates a correlation value of a reception signal of the two element antennas of the array antenna and calculates a correlation value of the two element antennas from the result. Calculating the reception phase difference to determine the incident direction, the second signal processing device calculates a partial correlation matrix and a cross-correlation vector by performing a correlation process on a part of the received signal of the array antenna, A weight value that minimizes the average power of the signal is calculated from the partial correlation matrix and the cross-correlation vector, and the incident direction is obtained by solving an equation of the same order as the incident wave number using the weight value as a coefficient. The signal processing device, recursively calculates a weight value that minimizes the average power of the signal by an adaptive filter, calculates an antenna pattern when multiplying each element of the array antenna by a weight value, And requests the incident direction of the radio wave by detecting a null point of the antenna pattern.

【0022】また、この発明に係る電波方位探知装置
は、前記入射波数推定手段が、入射波が周波数の異なる
狭帯域な信号である場合には受信信号のスペクトルから
入射波数を推定するものである。
In the radio azimuth detecting apparatus according to the present invention, the incident wave number estimating means estimates the number of incident waves from the spectrum of a received signal when the incident waves are narrow band signals having different frequencies. .

【0023】また、この発明に係る電波方位探知装置
は、前記入射波数推定手段が、前記アレーアンテナの受
信信号の相関処理を行い相関行列を計算し、前記相関行
列の固有値を計算し、前記固有値の大きさから入射波数
を推定するものである。
In the radio azimuth detecting apparatus according to the present invention, the incident wave number estimating means calculates a correlation matrix by performing a correlation process on a signal received by the array antenna, calculates an eigenvalue of the correlation matrix, and calculates an eigenvalue of the correlation matrix. Is used to estimate the number of incident waves from the magnitude of.

【0024】さらに、この発明に係る電波方位探知装置
は、電波を受信するアレーアンテナと、前記アレーアン
テナの受信信号の相関処理を行い相関行列を計算し、前
記相関行列の固有値を計算し、前記固有値の大きさから
入射波数を推定する入射波数推定手段と、前記入射波数
推定手段により推定された入射波数に応じて電波の入射
方向を求める信号処理装置を切り替える信号処理装置切
替手段と、前記入射波数が正の整数である第1の定数以
下である場合は、前記アレーアンテナの2つの素子アン
テナの受信信号の相関値を計算してその結果から2つの
素子アンテナの受信位相差を計算して入射方向を求める
第1の信号処理装置と、前記入射波数が前記第1の定数
より大きく、正の整数であって前記第1の定数より大き
い値の第2の定数より小さい場合は、前記相関行列から
信号の平均電力を最小化する荷重値を計算し、前記荷重
値を係数とする前記入射波数と同じ次数の方程式を解い
て入射方向を求める第2の信号処理装置と、前記入射波
数が前記第2の定数より大きい場合には、信号の平均電
力を最小化する荷重値を適応フィルタにより漸化的に計
算し、前記アレーアンテナの各素子に荷重値を乗じた場
合のアンテナパターンを計算し、前記アンテナパターン
のヌル点を検出して電波の入射方向を求める第3の信号
処理装置とを備えたものである。
Further, the radio azimuth detecting device according to the present invention calculates a correlation matrix by performing correlation processing of an array antenna receiving radio waves and a signal received by the array antenna, calculates an eigenvalue of the correlation matrix, Incident wave number estimating means for estimating the number of incident waves from the magnitude of the eigenvalue, signal processing device switching means for switching a signal processing device for obtaining the incident direction of the radio wave according to the incident wave number estimated by the incident wave number estimating means, When the wave number is equal to or smaller than a first constant that is a positive integer, the correlation value of the reception signals of the two element antennas of the array antenna is calculated, and the reception phase difference of the two element antennas is calculated from the result. A first signal processing device for determining an incident direction, and a second constant, wherein the number of incident waves is larger than the first constant, is a positive integer, and is larger than the first constant. If it is smaller, the second signal processing for calculating a load value that minimizes the average power of the signal from the correlation matrix and solving the equation of the same order as the incident wave number using the load value as a coefficient to obtain the incident direction. The apparatus, if the incident wave number is larger than the second constant, recursively calculates a weight value for minimizing the average power of the signal by an adaptive filter, and multiplies each element of the array antenna by the weight value. And a third signal processing device for calculating an antenna pattern in the case of the above, detecting a null point of the antenna pattern and obtaining an incident direction of a radio wave.

【0025】さらに、この発明に係る電波方位探知装置
は、電波を受信するアレーアンテナと、前記アレーアン
テナの受信信号の相関処理を行い相関行列を計算し、前
記相関行列の固有値を計算し、前記固有値の大きさから
入射波数を推定する入射波数推定手段と、前記入射波数
推定手段により推定された入射波数に応じて電波の入射
方向を求める信号処理装置を切り替える信号処理装置切
替手段と、前記入射波数が1である場合は、前記固有値
のうち最も大きな固有値に対応する固有ベクトルを算出
し、前記固有ベクトルをステアリングベクトルとして入
射方向を求める第1の信号処理装置と、前記入射波数が
1より大きく、正の整数であって1より大きい値の定数
より小さい場合は、前記相関行列から信号の平均電力を
最小化する荷重値を計算し、前記荷重値を係数とする前
記入射波数と同じ次数の方程式を解いて入射方向を求め
る第2の信号処理装置と、前記入射波数が前記定数より
大きい場合には、信号の平均電力を最小化する荷重値を
適応フィルタにより漸化的に計算し、前記アレーアンテ
ナの各素子に荷重値を乗じた場合のアンテナパターンを
計算し、前記アンテナパターンのヌル点を検出して電波
の入射方向を求める第3の信号処理装置とを備えたもの
である。
Further, the radio azimuth detecting device according to the present invention calculates an correlation antenna by performing a correlation process on an array antenna receiving radio waves and a signal received by the array antenna, calculates an eigenvalue of the correlation matrix, Incident wave number estimating means for estimating the number of incident waves from the magnitude of the eigenvalue, signal processing device switching means for switching a signal processing device for obtaining the incident direction of the radio wave according to the incident wave number estimated by the incident wave number estimating means, When the wave number is 1, a first signal processing device that calculates an eigenvector corresponding to the largest eigenvalue among the eigenvalues and determines an incident direction using the eigenvector as a steering vector, If the integer is less than a constant greater than 1, a weight value that minimizes the average power of the signal from the correlation matrix Calculate, solve the equation of the same order as the incident wave number with the load value as a coefficient, solve the equation of the incident direction, and, if the incident wave number is larger than the constant, the average power of the signal The weight value to be minimized is recursively calculated by the adaptive filter, the antenna pattern when each element of the array antenna is multiplied by the weight value is calculated, the null point of the antenna pattern is detected, and the incident direction of the radio wave is detected. And a third signal processing device for obtaining

【0026】さらに、この発明に係る電波方位探知装置
は、電波を受信するアレーアンテナと、前記アレーアン
テナの受信信号の相関処理を行い相関行列を計算し、前
記相関行列の固有値を計算し、前記固有値の大きさから
入射波数を推定する入射波数推定手段と、前記入射波数
推定手段により推定された入射波数に応じて電波の入射
方向を求める信号処理装置を切り替える信号処理装置切
替手段と、前記入射波数が1である場合は、前記固有値
のうち最も大きな固有値に対応する固有ベクトルを算出
し、前記固有ベクトルをステアリングベクトルとして入
射方向を求める第1の信号処理装置と、前記入射波数が
1より大きい場合には、ノイズに対応した固有値に対応
する固有ベクトルを求め、前記アレーアンテナの各方位
に対する位相特性や振幅特性を要素として持つステアリ
ングベクトルと前記固有ベクトルから方位評価関数を算
出し、前記方位評価関数のピークを検出して入射方向を
求める第2の信号処理装置とを備えたものである。
Further, the radio azimuth detecting device according to the present invention calculates a correlation matrix by performing correlation processing of an array antenna for receiving radio waves and a signal received by the array antenna, and calculates an eigenvalue of the correlation matrix. Incident wave number estimating means for estimating the number of incident waves from the magnitude of the eigenvalue, signal processing device switching means for switching a signal processing device for obtaining the incident direction of the radio wave according to the incident wave number estimated by the incident wave number estimating means, When the wave number is 1, a first signal processing device that calculates an eigenvector corresponding to the largest eigenvalue among the eigenvalues and obtains an incident direction using the eigenvector as a steering vector; Obtains an eigenvector corresponding to an eigenvalue corresponding to noise, and calculates a phase characteristic for each direction of the array antenna. Calculating the orientation evaluation function from the steering vector eigenvector having amplitude characteristics as an element, in which a second signal processing device for determining the direction of incidence by detecting the peak of the orientation evaluation function.

【0027】[0027]

【発明の実施の形態】実施の形態1.この発明の実施の
形態1に係る電波方位探知装置について図1を参照しな
がら説明する。図1は、この発明の実施の形態1に係る
電波方位探知装置の構成を示す図である。なお、各図
中、同一符号は同一又は相当部分を示す。
DESCRIPTION OF THE PREFERRED EMBODIMENTS Embodiment 1 The radio wave direction detecting device according to Embodiment 1 of the present invention will be described with reference to FIG. FIG. 1 is a diagram showing a configuration of a radio azimuth detecting device according to Embodiment 1 of the present invention. In the drawings, the same reference numerals indicate the same or corresponding parts.

【0028】図1において、1はアレーアンテナ、10
はアレーアンテナ1の受信信号から入射波数を推定する
入射波数推定手段、20は入射波数推定手段10で求め
られた入射波数に応じて電波方向探知アルゴリズムを切
り替えるアルゴリズム切替手段(信号処理装置切替手
段)、30は入射波数Kがある正の整数K1以下である
場合(K≦K1)にアレーアンテナ1の受信信号の相関
値より干渉計の原理で入射方位を求める信号処理装置
(アルゴリズム1)、40は入射波数KがK1より大き
く正の整数K2より小さい場合(K1<K<K2)にアレ
ーアンテナ1の受信信号の相関行列及び相互相関ベクト
ルを算出し相関行列の逆行列と相互相関ベクトルの積か
ら荷重値を計算してこれを係数として持つ方程式の解か
ら入射方位を算出する信号処理装置(アルゴリズム
2)、50は入射波数KがK2以上である場合(K2
K)に適応フィルタにより荷重値を計算し、この荷重値
から算出されるアンテナパターンのヌルから入射方位を
算出する信号処理装置(アルゴリズム3)である。
In FIG. 1, 1 is an array antenna, 10
Is an incident wave number estimating means for estimating the number of incident waves from the received signal of the array antenna 1, and 20 is an algorithm switching means (signal processing device switching means) for switching a radio wave direction finding algorithm according to the number of incident waves obtained by the incident wave number estimating means 10. If 30 is a positive integer K 1 below there is a incident wave number K (K ≦ K 1) to the signal processor for determining the incident orientation on the principle of the interferometer than the correlation value of the received signal of the array antenna 1 (algorithm 1) , 40 inverse matrix of the incident wave number K may largely positive integer K 2 smaller than K 1 (K 1 <K <K 2) calculating a correlation matrix and cross correlation vector of the received signal of the array antenna 1 in the correlation matrix A signal processing device (algorithm 2) for calculating a load value from the product of the cross-correlation vector and calculating the incident direction from the solution of an equation having the coefficient as a coefficient. If it is 2 or more (K 2
K) is a signal processing device (algorithm 3) that calculates a load value by an adaptive filter and calculates an incident azimuth from null of an antenna pattern calculated from the load value.

【0029】つぎに、この実施の形態1の動作について
説明する。まず、アレーアンテナ1の受信信号を入射波
数推定手段10に入力し、受信信号に含まれる入射波の
数を推定する。このような手段として、入射波が周波数
の異なる狭帯域な信号である場合には受信信号のスペク
トルから入射波数を推定することが出来る。
Next, the operation of the first embodiment will be described. First, the received signal of the array antenna 1 is input to the incident wave number estimating means 10 to estimate the number of incident waves included in the received signal. As such means, when the incident wave is a narrow band signal having a different frequency, the number of incident waves can be estimated from the spectrum of the received signal.

【0030】次に、アルゴリズム切替手段20によっ
て、入射波数推定手段10で求められた入射波数に応じ
て電波方位探知に用いるアルゴリズムを切り替える。ア
ルゴリズムを切り替える境界となる入射波数K1、K
2は、各アルゴリズムの計算量や精度により決定する。
ここでは、例として、K1=1、K2=5とする。また、
各素子アンテナは無指向性でありgm(θk)=1である
とする。
Next, the algorithm switching means 20 switches the algorithm used for radio wave direction detection according to the number of incident waves obtained by the number of incident waves estimation means 10. The incident wave numbers K 1 and K that are the boundaries for switching the algorithm
2 is determined by the calculation amount and accuracy of each algorithm.
Here, as an example, it is assumed that K 1 = 1 and K 2 = 5. Also,
It is assumed that each element antenna is omnidirectional and g mk ) = 1.

【0031】まず、入射波数が「1」である場合は、ア
ルゴリズム1を選択する。アルゴリズム1では、処理ス
テップ30により、2つの素子アンテナの受信信号の相
関値を計算してその結果から2つの素子アンテナの受信
位相差を計算して入射方向θを求める。
First, when the incident wave number is "1", the algorithm 1 is selected. In the algorithm 1, in the processing step 30, the correlation value of the reception signals of the two element antennas is calculated, and the reception phase difference of the two element antennas is calculated from the calculation result to obtain the incident direction θ.

【0032】例えば、入射波s0が入射するときアレー
アンテナ1の0番目および1番目のアンテナの受信信号
の相関値は式(1)〜(3)より次の式(6)のように
なる。
For example, when the incident wave s 0 is incident, the correlation values of the received signals of the 0th and 1st antennas of the array antenna 1 are as shown in the following expression (6) from the expressions (1) to (3). .

【0033】[0033]

【数5】 (Equation 5)

【0034】この相関値(複素数)の偏角から、入射方
向ベクトルp0を求めることができる。このように入射
波が1波である場合には、従来のMUSICアルゴリズ
ムによる多くの処理ステップを行う必要なく簡単に入射
方向を知ることが出来る。また、MUSICアルゴリズ
ムでは、評価方向ごとにノイズに対応する固有ベクトル
とステアリングベクトルの内積から評価関数を求めるた
め、入射方向と評価方向が一致せず精度が悪くなること
が考えられるが、この方法では式(6)から直接、入射
方向ベクトルp0を求めるので精度の面でも有効であ
る。
The incident direction vector p 0 can be obtained from the argument of the correlation value (complex number). As described above, when the number of incident waves is one, it is possible to easily know the incident direction without performing many processing steps by the conventional MUSIC algorithm. Further, in the MUSIC algorithm, since the evaluation function is obtained from the inner product of the eigenvector corresponding to noise and the steering vector for each evaluation direction, the incident direction and the evaluation direction do not match, which may degrade the accuracy. Since the incident direction vector p 0 is obtained directly from (6), it is also effective in terms of accuracy.

【0035】次に、入射波数が2以上4以下である場合
について説明する。この場合には、アルゴリズム2によ
る信号処理装置40を選択する。MUSICアルゴリズ
ムの場合と同様にアレーアンテナ1の一部(例えば、3
素子のアレーアンテナであれば2素子、5素子のアレー
アンテナであれば4素子)の素子アンテナの相関行列、
および相互相関ベクトルを計算する。アレーアンテナ1
の素子数をMとし、各素子アンテナの受信信号をx
0(t)、x1(t)、・・・、xM-1(t)とする。こ
のとき、x1(t)〜xM-1(t)の相関行列R、及びx
0(t)とx1(t)〜xM-1(t)の相互相関ベクトル
Pをステップ41において次の式(7)及び式(8)に
示すように求める。
Next, the case where the number of incident waves is 2 or more and 4 or less will be described. In this case, the signal processing device 40 based on the algorithm 2 is selected. As in the case of the MUSIC algorithm, a part of the array antenna 1 (for example, 3
A correlation matrix of two element antennas for an element antenna, four elements for a five element array antenna,
And calculate the cross-correlation vector. Array antenna 1
Is the number of elements of M, and the received signal of each element antenna is x
0 (t), x 1 ( t), ···, and x M-1 (t). At this time, the correlation matrix R of x 1 (t) to x M-1 (t) and x
0 (t) and obtaining a cross-correlation vector P of x 1 (t) ~x M- 1 (t) as shown in the following equation (7) and (8) in step 41.

【0036】[0036]

【数6】 (Equation 6)

【0037】次に、式(9)により計算される信号e
(t)の平均電力を最小化する荷重w1〜wM-1をステッ
プ42において式(10)のように求める。
Next, the signal e calculated by the equation (9)
In step 42 the load w 1 ~w M-1 that minimizes the average power of (t) determined by the equation (10).

【0038】[0038]

【数7】 (Equation 7)

【0039】次に、ステップ43において、式(10)
で求められた荷重を係数として持つ次の式(11)に示
すようなuについてのK次方程式を解く。
Next, in step 43, equation (10)
The K-order equation for u as shown in the following equation (11) having the load obtained as a coefficient as a coefficient is solved.

【0040】[0040]

【数8】 (Equation 8)

【0041】上記の式(11)の解(複素数)の偏角か
ら、入射角が求められる。入射波が1波の場合と同様に
入射角を解析的に解くので、精度と演算量の面で有効で
ある。しかし、式(10)における逆行列の計算、及び
式(11)のK次方程式を解く計算は入射波の数が多
く、次数が増加すると極度に難しくなる。このため、あ
る次数以上では、パラメータを振って解を見つけなくて
はならない。
The incident angle is obtained from the argument of the solution (complex number) of the above equation (11). Since the incident angle is analytically solved as in the case of one incident wave, it is effective in terms of accuracy and calculation amount. However, the calculation of the inverse matrix in equation (10) and the calculation to solve the K-th order equation in equation (11) involve a large number of incident waves, and becomes extremely difficult as the order increases. Therefore, above a certain order, it is necessary to find a solution by varying parameters.

【0042】上記のような理由から、入射波数が5以上
となる場合にはアルゴリズム3による信号処理装置50
を選択する。まず、ステップ51において、式(9)の
平均電力を最小化する荷重w1〜wM-1を適応フィルタに
より漸化的に求める。
For the above reason, when the number of incident waves is 5 or more, the signal processing device 50 based on the algorithm 3
Select First, in step 51, recurrently seeking load w 1 ~w M-1 that minimizes the average power of the formula (9) by the adaptive filter.

【0043】次に、ステップ52において、式(12)
で計算されるようなアレーアンテナ1の各素子に荷重w
1〜wM-1を乗じた場合のアンテナパターンを計算する。
Next, in step 52, equation (12)
The load w is applied to each element of the array antenna 1 as calculated by
Calculating the antenna pattern when multiplied by 1 ~w M-1.

【0044】[0044]

【数9】 (Equation 9)

【0045】適応フィルタは、入射波の方向のゲインを
最小化するように荷重を制御するので、上記の式(1
2)で計算されるパターンには入射方向にヌルが形成さ
れる。この結果、入射波の入射方向θを知ることができ
る。
The adaptive filter controls the load so as to minimize the gain in the direction of the incident wave.
Nulls are formed in the incident direction in the pattern calculated in 2). As a result, the incident direction θ of the incident wave can be known.

【0046】以上説明したように、この実施の形態1は
入射波の数に応じて方向探知アルゴリズム(信号処理装
置)を切り替えるので、特に入射波が少ない場合に、高
精度で演算量の少ない方向探知を行うことが出来る。
As described above, in the first embodiment, the direction detection algorithm (signal processing device) is switched according to the number of incident waves. Can detect.

【0047】実施の形態2.この発明の実施の形態2に
係る電波方位探知装置について図2を参照しながら説明
する。図2は、この発明の実施の形態2に係る電波方位
探知装置の構成を示す図である。
Embodiment 2 Embodiment 2 A radio azimuth detecting device according to Embodiment 2 of the present invention will be described with reference to FIG. FIG. 2 is a diagram showing a configuration of a radio wave direction detecting device according to Embodiment 2 of the present invention.

【0048】この実施の形態2では、上記の実施の形態
1の処理のうち、入射波数推定手段10を以下のような
ステップ処理により構成する。
In the second embodiment, of the processing of the first embodiment, the incident wave number estimating means 10 is constituted by the following step processing.

【0049】まず、ステップ11において、アレーアン
テナ1の受信信号を入力して次の式(13)に示すよう
な相関行列Raを算出する。
[0049] First, in step 11, and inputs the received signal of the array antenna 1 calculates a correlation matrix R a as shown in the following equation (13).

【0050】[0050]

【数10】 (Equation 10)

【0051】次に、ステップ12〜13において、相関
行列の固有値を算出し、固有値の大きさから入射波数を
推定する。固有値の大きさは、受信機雑音に対応するも
のは小さく一般にばらつきが小さい。入射信号に対応す
る固有値はノイズに対応するものよりも大きいので、入
射波数を判別することが出来る。
Next, in steps 12 and 13, the eigenvalue of the correlation matrix is calculated, and the number of incident waves is estimated from the magnitude of the eigenvalue. The magnitude of the eigenvalue corresponding to the receiver noise is small and generally has a small variation. Since the eigenvalue corresponding to the incident signal is larger than that corresponding to noise, the number of incident waves can be determined.

【0052】このような方法で入射波数を判別した場合
には、上記の実施の形態1におけるステップ41の処理
は不要となる。なぜならば、相関行列Rは相関行列Ra
の2〜(M−1)行、2〜(M−1)列の要素で構成さ
れる行列であり、相互相関ベクトルPは相関行列Ra
2〜(M−1)行、1列の要素で構成される行列である
からである。このような構成にすることにより、上記の
実施の形態1と同様な効果を得ることが出来る。
When the number of incident waves is determined by such a method, the processing of step 41 in the first embodiment becomes unnecessary. Because the correlation matrix R is the correlation matrix R a
Of 2 (M-1) row, 2 (M-1) is a matrix composed of the column, the cross-correlation vector P is 2 to the correlation matrix R a (M-1) row, the first column This is because the matrix is composed of elements. With such a configuration, the same effect as in the first embodiment can be obtained.

【0053】実施の形態3.この発明の実施の形態3に
係る電波方位探知装置について図3を参照しながら説明
する。図3は、この発明の実施の形態3に係る電波方位
探知装置の構成を示す図である。
Embodiment 3 Embodiment 3 A radio azimuth detecting device according to Embodiment 3 of the present invention will be described with reference to FIG. FIG. 3 is a diagram showing a configuration of a radio wave direction detecting apparatus according to Embodiment 3 of the present invention.

【0054】この実施の形態3では、入射波数推定手段
10の推定結果が1であった場合に、上記の実施の形態
1及び2のアルゴリズム1(ステップ30)の処理に代
わってステツプ60の処理(信号処理装置)により入射
波の方向探知を行う。
In the third embodiment, when the estimation result of the incident wave number estimating means 10 is 1, the processing of the step 60 replaces the processing of the algorithm 1 (step 30) of the first and second embodiments. (Signal processing device) The direction of the incident wave is detected.

【0055】このステップ60(アルゴリズム1)は、
入射波が1波である場合にステップ12で求められた固
有値のなかで最も大きな固有値に対応する固有ベクトル
を求める。この固有ベクトルは入射波のステアリングベ
クトルa(θ0)であるので、これより入射方向θを求
めることができる。このような構成にすることにより、
上記の実施の形態1及び2と同様な効果を得ることが出
来る。
This step 60 (algorithm 1)
When the number of incident waves is one, an eigenvector corresponding to the largest eigenvalue among the eigenvalues obtained in step 12 is obtained. Since this eigenvector is the steering vector a (θ 0 ) of the incident wave, the incident direction θ can be obtained from this. With such a configuration,
The same effects as in the first and second embodiments can be obtained.

【0056】実施の形態4.この発明の実施の形態4に
係る電波方位探知装置について図4を参照しながら説明
する。図4は、この発明の実施の形態4に係る電波方位
探知装置の構成を示す図である。
Embodiment 4 Fourth Embodiment A radio azimuth detecting device according to a fourth embodiment of the present invention will be described with reference to FIG. FIG. 4 is a diagram showing a configuration of a radio wave direction detecting device according to Embodiment 4 of the present invention.

【0057】この実施の形態4では、アルゴリズム切替
手段20において、入射波数KがK1より大きい場合
(K1<K)にアルゴリズム2として、従来のMUSI
Cアルゴリズムと同様の処理を信号処理装置70におい
て行う。すなわち、ステップ5〜7において、ノイズに
対応する固有ベクトルを算出し、この固有ベクトルと評
価する方位のステアリングベクトルとから式(5)に示
すように方位評価関数P(θ)を算出する。この方位評
価関数P(θ)に現れるピークから入射波の入射方向θ
を求める。このような構成にすることにより、上記実施
の形態1ないし3と同様な効果を得ることが出来る。
[0057] In the fourth embodiment, in the algorithm switching means 20, as an algorithm 2 when the incident wave number K is greater than K 1 (K 1 <K), the conventional MUSI
Processing similar to the C algorithm is performed in the signal processing device 70. That is, in steps 5 to 7, the eigenvector corresponding to the noise is calculated, and the azimuth evaluation function P (θ) is calculated from the eigenvector and the steering vector of the azimuth to be evaluated as shown in Expression (5). From the peak appearing in this azimuth evaluation function P (θ), the incident direction θ of the incident wave
Ask for. With such a configuration, the same effect as in the first to third embodiments can be obtained.

【0058】[0058]

【発明の効果】この発明に係る電波方位探知装置は、以
上説明したとおり、電波を受信するアレーアンテナと、
前記アレーアンテナの受信信号から入射波数を推定する
入射波数推定手段と、前記入射波数推定手段により推定
された入射波数に応じて電波の入射方向を求める信号処
理装置を切り替える信号処理装置切替手段と、前記入射
波数が正の整数である第1の定数以下である場合は、前
記アレーアンテナの受信信号の相関値に基づいて電波の
入射方向を求める第1の信号処理装置と、前記入射波数
が前記第1の定数より大きく、正の整数であって前記第
1の定数より大きい値の第2の定数より小さい場合は、
前記アレーアンテナの受信信号の相関行列及び相互相関
ベクトルを計算し、前記相関行列及び相互相関ベクトル
から荷重値を計算し、前記荷重値を係数とする所定の方
程式を解いて電波の入射方向を求める第2の信号処理装
置と、前記入射波数が前記第2の定数より大きい場合に
は、前記アレーアンテナの受信信号から適応フィルタに
より荷重値を計算し、前記荷重値からアンテナパターン
を計算し、前記アンテナパターンのヌル点を検出して電
波の入射方向を求める第3の信号処理装置とを備えたの
で、特に入射波数が少ない場合には精度を高くできると
ともに演算量を少なくすることができるという効果を奏
する。
As described above, the radio azimuth detecting device according to the present invention comprises: an array antenna for receiving radio waves;
Incident wave number estimating means for estimating the number of incident waves from the received signal of the array antenna, and a signal processing device switching means for switching a signal processing device that determines the incident direction of a radio wave according to the number of incident waves estimated by the incident wave number estimating means, When the number of incident waves is equal to or less than a first constant that is a positive integer, a first signal processing device that determines an incident direction of a radio wave based on a correlation value of a reception signal of the array antenna; If the value is greater than the first constant, is a positive integer, and is smaller than the second constant having a value larger than the first constant,
A correlation matrix and a cross-correlation vector of the received signal of the array antenna are calculated, a load value is calculated from the correlation matrix and the cross-correlation vector, and a predetermined equation using the load value as a coefficient is solved to obtain an incident direction of a radio wave. A second signal processing device, when the number of incident waves is larger than the second constant, calculate a weight value by an adaptive filter from a reception signal of the array antenna, calculate an antenna pattern from the weight value, Since a third signal processing device for detecting the null point of the antenna pattern to obtain the incident direction of the radio wave is provided, the accuracy can be improved and the amount of calculation can be reduced particularly when the number of incident waves is small. To play.

【0059】また、この発明に係る電波方位探知装置
は、以上説明したとおり、前記第1の信号処理装置が、
前記アレーアンテナの2つの素子アンテナの受信信号の
相関値を計算してその結果から2つの素子アンテナの受
信位相差を計算して入射方向を求め、前記第2の信号処
理装置が、前記アレーアンテナの一部の受信信号の相関
処理を行って部分的な相関行列と相互相関ベクトルを計
算し、前記部分的な相関行列と相互相関ベクトルから信
号の平均電力を最小化する荷重値を計算し、前記荷重値
を係数とする前記入射波数と同じ次数の方程式を解いて
入射方向を求め、前記第3の信号処理装置が、信号の平
均電力を最小化する荷重値を適応フィルタにより漸化的
に計算し、前記アレーアンテナの各素子に荷重値を乗じ
た場合のアンテナパターンを計算し、前記アンテナパタ
ーンのヌル点を検出して電波の入射方向を求めるので、
特に入射波数が少ない場合には精度を高くできるととも
に演算量を少なくすることができるという効果を奏す
る。
Further, as described above, the radio azimuth detecting device according to the present invention includes the first signal processing device,
The correlation value of the reception signals of the two element antennas of the array antenna is calculated, and the reception phase difference of the two element antennas is calculated from the calculation result to determine the incident direction. A partial correlation matrix and a cross-correlation vector are calculated by performing a correlation process on a part of the received signal, and a weight value that minimizes the average power of the signal from the partial correlation matrix and the cross-correlation vector is calculated, The incident direction is determined by solving the equation of the same order as the incident wave number having the load value as a coefficient, and the third signal processing device gradually reduces the load value for minimizing the average power of the signal by an adaptive filter. Calculate, calculate the antenna pattern when multiplying each element of the array antenna by a load value, and detect the null point of the antenna pattern to determine the incident direction of the radio wave,
In particular, when the number of incident waves is small, there is an effect that the accuracy can be increased and the amount of calculation can be reduced.

【0060】また、この発明に係る電波方位探知装置
は、以上説明したとおり、前記入射波数推定手段が、入
射波が周波数の異なる狭帯域な信号である場合には受信
信号のスペクトルから入射波数を推定するので、特に入
射波数が少ない場合には精度を高くできるとともに演算
量を少なくすることができるという効果を奏する。
As described above, in the radio azimuth detecting device according to the present invention, when the incident wave number estimating means determines that the number of incident waves is a narrow band signal having a different frequency, the number of incident waves is determined from the spectrum of the received signal. Since the estimation is performed, there is an effect that the accuracy can be improved and the amount of calculation can be reduced particularly when the number of incident waves is small.

【0061】また、この発明に係る電波方位探知装置
は、以上説明したとおり、前記入射波数推定手段が、前
記アレーアンテナの受信信号の相関処理を行い相関行列
を計算し、前記相関行列の固有値を計算し、前記固有値
の大きさから入射波数を推定するので、特に入射波数が
少ない場合には精度を高くできるとともに演算量を少な
くすることができるという効果を奏する。
In the radio azimuth detecting device according to the present invention, as described above, the incident wave number estimating means calculates a correlation matrix by performing a correlation process on the reception signal of the array antenna, and calculates an eigenvalue of the correlation matrix. Since the calculation is performed and the number of incident waves is estimated from the magnitude of the eigenvalue, there is an effect that the accuracy can be increased and the amount of calculation can be reduced particularly when the number of incident waves is small.

【0062】さらに、この発明に係る電波方位探知装置
は、以上説明したとおり、電波を受信するアレーアンテ
ナと、前記アレーアンテナの受信信号の相関処理を行い
相関行列を計算し、前記相関行列の固有値を計算し、前
記固有値の大きさから入射波数を推定する入射波数推定
手段と、前記入射波数推定手段により推定された入射波
数に応じて電波の入射方向を求める信号処理装置を切り
替える信号処理装置切替手段と、前記入射波数が正の整
数である第1の定数以下である場合は、前記アレーアン
テナの2つの素子アンテナの受信信号の相関値を計算し
てその結果から2つの素子アンテナの受信位相差を計算
して入射方向を求める第1の信号処理装置と、前記入射
波数が前記第1の定数より大きく、正の整数であって前
記第1の定数より大きい値の第2の定数より小さい場合
は、前記相関行列から信号の平均電力を最小化する荷重
値を計算し、前記荷重値を係数とする前記入射波数と同
じ次数の方程式を解いて入射方向を求める第2の信号処
理装置と、前記入射波数が前記第2の定数より大きい場
合には、信号の平均電力を最小化する荷重値を適応フィ
ルタにより漸化的に計算し、前記アレーアンテナの各素
子に荷重値を乗じた場合のアンテナパターンを計算し、
前記アンテナパターンのヌル点を検出して電波の入射方
向を求める第3の信号処理装置とを備えたので、特に入
射波数が少ない場合には精度を高くできるとともに演算
量を少なくすることができるという効果を奏する。
Further, as described above, the radio azimuth detecting device according to the present invention calculates the correlation matrix by performing the correlation processing of the array antenna for receiving the radio wave and the signal received by the array antenna, and calculates the eigenvalue of the correlation matrix. Signal processing device switching means for calculating an incident wave number estimating means for estimating an incident wave number from the magnitude of the eigenvalue, and a signal processing device for obtaining an incident direction of a radio wave according to the incident wave number estimated by the incident wave number estimating means. Means for calculating the correlation value of the reception signals of the two element antennas of the array antenna when the number of incident waves is equal to or less than a first constant that is a positive integer, and calculating the reception position of the two element antennas from the result. A first signal processing apparatus for calculating an incident direction by calculating a phase difference, wherein the number of incident waves is larger than the first constant and is a positive integer and is larger than the first constant. If it is smaller than the second constant of the threshold value, a load value that minimizes the average power of the signal is calculated from the correlation matrix, and an equation of the same order as the incident wave number using the load value as a coefficient is solved to obtain the incident direction. And a second signal processing device for calculating the weighting value for minimizing the average power of the signal when the incident wave number is larger than the second constant. Calculate the antenna pattern when multiplying each element by the load value,
Since a third signal processing device for detecting the null point of the antenna pattern to determine the incident direction of the radio wave is provided, especially when the number of incident waves is small, the accuracy can be improved and the amount of calculation can be reduced. It works.

【0063】さらに、この発明に係る電波方位探知装置
は、以上説明したとおり、電波を受信するアレーアンテ
ナと、前記アレーアンテナの受信信号の相関処理を行い
相関行列を計算し、前記相関行列の固有値を計算し、前
記固有値の大きさから入射波数を推定する入射波数推定
手段と、前記入射波数推定手段により推定された入射波
数に応じて電波の入射方向を求める信号処理装置を切り
替える信号処理装置切替手段と、前記入射波数が1であ
る場合は、前記固有値のうち最も大きな固有値に対応す
る固有ベクトルを算出し、前記固有ベクトルをステアリ
ングベクトルとして入射方向を求める第1の信号処理装
置と、前記入射波数が1より大きく、正の整数であって
1より大きい値の定数より小さい場合は、前記相関行列
から信号の平均電力を最小化する荷重値を計算し、前記
荷重値を係数とする前記入射波数と同じ次数の方程式を
解いて入射方向を求める第2の信号処理装置と、前記入
射波数が前記定数より大きい場合には、信号の平均電力
を最小化する荷重値を適応フィルタにより漸化的に計算
し、前記アレーアンテナの各素子に荷重値を乗じた場合
のアンテナパターンを計算し、前記アンテナパターンの
ヌル点を検出して電波の入射方向を求める第3の信号処
理装置とを備えたので、特に入射波数が少ない場合には
精度を高くできるとともに演算量を少なくすることがで
きるという効果を奏する。
Further, as described above, the radio azimuth detecting device according to the present invention calculates the correlation matrix by performing correlation processing of an array antenna for receiving radio waves and a signal received by the array antenna, and calculates an eigenvalue of the correlation matrix. Signal processing device switching means for calculating an incident wave number estimating means for estimating an incident wave number from the magnitude of the eigenvalue, and a signal processing device for obtaining an incident direction of a radio wave according to the incident wave number estimated by the incident wave number estimating means. A means for calculating an eigenvector corresponding to the largest eigenvalue among the eigenvalues when the number of incident waves is 1, and a first signal processing device for determining an incident direction using the eigenvector as a steering vector; If it is greater than 1 and less than a constant that is a positive integer and greater than 1, the average power of the signal from the correlation matrix Calculate the load value to minimize the, the second signal processing device to determine the incident direction by solving the equation of the same order as the incident wave number with the load value as a coefficient, when the incident wave number is greater than the constant Calculates the weight value that minimizes the average power of the signal recursively by an adaptive filter, calculates the antenna pattern when each element of the array antenna is multiplied by the weight value, and calculates the null point of the antenna pattern. Since the third signal processing device for detecting the incident direction of the radio wave by detecting the signal is provided, it is possible to increase the accuracy and reduce the amount of calculation, particularly when the number of incident waves is small.

【0064】さらに、この発明に係る電波方位探知装置
は、以上説明したとおり、電波を受信するアレーアンテ
ナと、前記アレーアンテナの受信信号の相関処理を行い
相関行列を計算し、前記相関行列の固有値を計算し、前
記固有値の大きさから入射波数を推定する入射波数推定
手段と、前記入射波数推定手段により推定された入射波
数に応じて電波の入射方向を求める信号処理装置を切り
替える信号処理装置切替手段と、前記入射波数が1であ
る場合は、前記固有値のうち最も大きな固有値に対応す
る固有ベクトルを算出し、前記固有ベクトルをステアリ
ングベクトルとして入射方向を求める第1の信号処理装
置と、前記入射波数が1より大きい場合には、ノイズに
対応した固有値に対応する固有ベクトルを求め、前記ア
レーアンテナの各方位に対する位相特性や振幅特性を要
素として持つステアリングベクトルと前記固有ベクトル
から方位評価関数を算出し、前記方位評価関数のピーク
を検出して入射方向を求める第2の信号処理装置とを備
えたので、特に入射波数が少ない場合には精度を高くで
きるとともに演算量を少なくすることができるという効
果を奏する。
Further, as described above, the radio azimuth detecting device according to the present invention calculates the correlation matrix by performing the correlation processing of the array antenna receiving the radio wave and the signal received by the array antenna, and calculates the eigenvalue of the correlation matrix. Signal processing device switching means for calculating an incident wave number estimating means for estimating an incident wave number from the magnitude of the eigenvalue, and a signal processing device for obtaining an incident direction of a radio wave according to the incident wave number estimated by the incident wave number estimating means. A means for calculating an eigenvector corresponding to the largest eigenvalue among the eigenvalues when the number of incident waves is 1, and a first signal processing device for determining an incident direction using the eigenvector as a steering vector; If it is greater than 1, an eigenvector corresponding to the eigenvalue corresponding to the noise is obtained, and each of the array antennas A second signal processing device for calculating an azimuth evaluation function from a steering vector having phase characteristics and amplitude characteristics with respect to the position as elements and the eigenvector and detecting a peak of the azimuth evaluation function to obtain an incident direction. In particular, when the number of incident waves is small, there is an effect that the accuracy can be increased and the amount of calculation can be reduced.

【図面の簡単な説明】[Brief description of the drawings]

【図1】 この発明の実施の形態1に係る電波方位探知
装置の構成を示す図である。
FIG. 1 is a diagram showing a configuration of a radio azimuth detecting device according to Embodiment 1 of the present invention.

【図2】 この発明の実施の形態2に係る電波方位探知
装置の構成を示す図である。
FIG. 2 is a diagram showing a configuration of a radio azimuth detecting device according to Embodiment 2 of the present invention.

【図3】 この発明の実施の形態3に係る電波方位探知
装置の構成を示す図である。
FIG. 3 is a diagram showing a configuration of a radio azimuth detecting device according to Embodiment 3 of the present invention.

【図4】 この発明の実施の形態4に係る電波方位探知
装置の構成を示す図である。
FIG. 4 is a diagram showing a configuration of a radio wave direction detecting apparatus according to Embodiment 4 of the present invention.

【図5】 従来の電波方位探知装置の構成を示す図であ
る。
FIG. 5 is a diagram showing a configuration of a conventional radio wave direction detecting device.

【符号の説明】[Explanation of symbols]

1 アレーアンテナ、10 入射波数推定手段、20
アルゴリズム切替手段、30 アルゴリズム1による信
号処理装置、40 アルゴリズム2による信号処理装
置、50 アルゴリズム3による信号処理装置、60
アルゴリズム1による信号処理装置、70 アルゴリズ
ム2による信号処理装置。
1 array antenna, 10 incident wave number estimation means, 20
Algorithm switching means, 30 a signal processing device using algorithm 1, 40 a signal processing device using algorithm 2, 50 a signal processing device using algorithm 3, 60
Signal processing device according to algorithm 1, 70 Signal processing device according to algorithm 2.

Claims (7)

【特許請求の範囲】[Claims] 【請求項1】 電波を受信するアレーアンテナと、 前記アレーアンテナの受信信号から入射波数を推定する
入射波数推定手段と、 前記入射波数推定手段により推定された入射波数に応じ
て電波の入射方向を求める信号処理装置を切り替える信
号処理装置切替手段と、 前記入射波数が正の整数である第1の定数以下である場
合は、前記アレーアンテナの受信信号の相関値に基づい
て電波の入射方向を求める第1の信号処理装置と、 前記入射波数が前記第1の定数より大きく、正の整数で
あって前記第1の定数より大きい値の第2の定数より小
さい場合は、前記アレーアンテナの受信信号の相関行列
及び相互相関ベクトルを計算し、前記相関行列及び相互
相関ベクトルから荷重値を計算し、前記荷重値を係数と
する所定の方程式を解いて電波の入射方向を求める第2
の信号処理装置と、 前記入射波数が前記第2の定数より大きい場合には、前
記アレーアンテナの受信信号から適応フィルタにより荷
重値を計算し、前記荷重値からアンテナパターンを計算
し、前記アンテナパターンのヌル点を検出して電波の入
射方向を求める第3の信号処理装置とを備えたことを特
徴とする電波方位探知装置。
An array antenna for receiving a radio wave, an incident wave number estimating means for estimating an incident wave number from a signal received by the array antenna, and an incident direction of the radio wave according to the incident wave number estimated by the incident wave number estimating means. Signal processing device switching means for switching a signal processing device to be obtained; and, when the number of incident waves is equal to or less than a first constant that is a positive integer, calculating an incident direction of a radio wave based on a correlation value of a reception signal of the array antenna. A first signal processing device; a signal received by the array antenna when the incident wave number is larger than the first constant and smaller than a second constant that is a positive integer and larger than the first constant; Calculating a correlation matrix and a cross-correlation vector, calculating a weight value from the correlation matrix and the cross-correlation vector, solving a predetermined equation using the weight value as a coefficient, Second determining the elevation direction
When the number of incident waves is larger than the second constant, a weight value is calculated from a reception signal of the array antenna by an adaptive filter, an antenna pattern is calculated from the weight value, and the antenna pattern And a third signal processing device for detecting a null point of the radio wave to determine an incident direction of the radio wave.
【請求項2】 前記第1の信号処理装置は、前記アレー
アンテナの2つの素子アンテナの受信信号の相関値を計
算してその結果から2つの素子アンテナの受信位相差を
計算して入射方向を求め、 前記第2の信号処理装置は、前記アレーアンテナの一部
の受信信号の相関処理を行って部分的な相関行列と相互
相関ベクトルを計算し、前記部分的な相関行列と相互相
関ベクトルから信号の平均電力を最小化する荷重値を計
算し、前記荷重値を係数とする前記入射波数と同じ次数
の方程式を解いて入射方向を求め、 前記第3の信号処理装置は、信号の平均電力を最小化す
る荷重値を適応フィルタにより漸化的に計算し、前記ア
レーアンテナの各素子に荷重値を乗じた場合のアンテナ
パターンを計算し、前記アンテナパターンのヌル点を検
出して電波の入射方向を求めることを特徴とする請求項
1記載の電波方位探知装置。
2. The first signal processing device calculates a correlation value of reception signals of two element antennas of the array antenna, calculates a reception phase difference of the two element antennas from the calculation result, and changes an incident direction. The second signal processing device calculates a partial correlation matrix and a cross-correlation vector by performing a correlation process on a part of the received signal of the array antenna, and calculates the partial correlation matrix and the cross-correlation vector from the partial correlation matrix and the cross-correlation vector. A load value that minimizes the average power of the signal is calculated, and an incident direction is obtained by solving an equation of the same order as the incident wave number using the load value as a coefficient. The third signal processing device calculates an average power of the signal. Is calculated recursively by an adaptive filter, the antenna pattern when each element of the array antenna is multiplied by the load value is calculated, and a null point of the antenna pattern is detected to detect radio waves. Telecommunications orientation detection device of claim 1, wherein the determining the direction of incidence.
【請求項3】 前記入射波数推定手段は、入射波が周波
数の異なる狭帯域な信号である場合には受信信号のスペ
クトルから入射波数を推定することを特徴とする請求項
1又は2記載の電波方位探知装置。
3. The radio wave according to claim 1, wherein said incident wave number estimating means estimates the incident wave number from a spectrum of a received signal when the incident wave is a narrow band signal having a different frequency. Direction detection device.
【請求項4】 前記入射波数推定手段は、前記アレーア
ンテナの受信信号の相関処理を行い相関行列を計算し、
前記相関行列の固有値を計算し、前記固有値の大きさか
ら入射波数を推定することを特徴とする請求項1又は2
記載の電波方位探知装置。
4. The incident wave number estimating means calculates a correlation matrix by performing a correlation process on a reception signal of the array antenna,
The eigenvalue of the correlation matrix is calculated, and the number of incident waves is estimated from the magnitude of the eigenvalue.
Radio azimuth detection device as described.
【請求項5】 電波を受信するアレーアンテナと、 前記アレーアンテナの受信信号の相関処理を行い相関行
列を計算し、前記相関行列の固有値を計算し、前記固有
値の大きさから入射波数を推定する入射波数推定手段
と、 前記入射波数推定手段により推定された入射波数に応じ
て電波の入射方向を求める信号処理装置を切り替える信
号処理装置切替手段と、 前記入射波数が正の整数である第1の定数以下である場
合は、前記アレーアンテナの2つの素子アンテナの受信
信号の相関値を計算してその結果から2つの素子アンテ
ナの受信位相差を計算して入射方向を求める第1の信号
処理装置と、 前記入射波数が前記第1の定数より大きく、正の整数で
あって前記第1の定数より大きい値の第2の定数より小
さい場合は、前記相関行列から信号の平均電力を最小化
する荷重値を計算し、前記荷重値を係数とする前記入射
波数と同じ次数の方程式を解いて入射方向を求める第2
の信号処理装置と、 前記入射波数が前記第2の定数より大きい場合には、信
号の平均電力を最小化する荷重値を適応フィルタにより
漸化的に計算し、前記アレーアンテナの各素子に荷重値
を乗じた場合のアンテナパターンを計算し、前記アンテ
ナパターンのヌル点を検出して電波の入射方向を求める
第3の信号処理装置とを備えたことを特徴とする電波方
位探知装置。
5. An array antenna for receiving a radio wave, and a correlation matrix is calculated by performing a correlation process on a signal received by the array antenna, an eigenvalue of the correlation matrix is calculated, and the number of incident waves is estimated from the magnitude of the eigenvalue. An incident wave number estimating unit; a signal processing device switching unit that switches a signal processing device that obtains an incident direction of a radio wave according to the incident wave number estimated by the incident wave number estimating unit; and a first in which the incident wave number is a positive integer. If not greater than a constant, a first signal processing apparatus for calculating the correlation value of the received signals of the two element antennas of the array antenna and calculating the reception phase difference between the two element antennas from the result to determine the incident direction And if the incident wave number is greater than the first constant, is a positive integer and is smaller than a second constant having a value greater than the first constant, a signal from the correlation matrix Calculates a load value that minimizes the average power, the second to determine the incident direction by solving an equation of the same order as the incident wave number and coefficient of the load value
When the number of incident waves is larger than the second constant, a weight value for minimizing the average power of the signal is recursively calculated by an adaptive filter, and a weight is applied to each element of the array antenna. A third signal processing device that calculates an antenna pattern when multiplied by a value, detects a null point of the antenna pattern, and obtains an incident direction of the radio wave.
【請求項6】 電波を受信するアレーアンテナと、 前記アレーアンテナの受信信号の相関処理を行い相関行
列を計算し、前記相関行列の固有値を計算し、前記固有
値の大きさから入射波数を推定する入射波数推定手段
と、 前記入射波数推定手段により推定された入射波数に応じ
て電波の入射方向を求める信号処理装置を切り替える信
号処理装置切替手段と、 前記入射波数が1である場合は、前記固有値のうち最も
大きな固有値に対応する固有ベクトルを算出し、前記固
有ベクトルをステアリングベクトルとして入射方向を求
める第1の信号処理装置と、 前記入射波数が1より大きく、正の整数であって1より
大きい値の定数より小さい場合は、前記相関行列から信
号の平均電力を最小化する荷重値を計算し、前記荷重値
を係数とする前記入射波数と同じ次数の方程式を解いて
入射方向を求める第2の信号処理装置と、 前記入射波数が前記定数より大きい場合には、信号の平
均電力を最小化する荷重値を適応フィルタにより漸化的
に計算し、前記アレーアンテナの各素子に荷重値を乗じ
た場合のアンテナパターンを計算し、前記アンテナパタ
ーンのヌル点を検出して電波の入射方向を求める第3の
信号処理装置とを備えたことを特徴とする電波方位探知
装置。
6. An array antenna for receiving a radio wave, a correlation matrix is calculated by performing a correlation process on a signal received by the array antenna, an eigenvalue of the correlation matrix is calculated, and the number of incident waves is estimated from the magnitude of the eigenvalue. An incident wave number estimating unit; a signal processing device switching unit that switches a signal processing device that determines an incident direction of a radio wave according to the incident wave number estimated by the incident wave number estimating unit; and the eigenvalue when the incident wave number is 1. A first signal processing device that calculates an eigenvector corresponding to the largest eigenvalue and determines the incident direction using the eigenvector as a steering vector; and the incident wave number is greater than 1 and is a positive integer greater than 1. If the value is smaller than the constant, a load value for minimizing the average power of the signal is calculated from the correlation matrix, and the input value having the load value as a coefficient is calculated. A second signal processing device for solving the equation of the same order as the wave number to obtain the incident direction; and when the incident wave number is larger than the constant, a load value for minimizing the average power of the signal is gradually reduced by an adaptive filter. And a third signal processing device for calculating an antenna pattern when each element of the array antenna is multiplied by a load value, and detecting a null point of the antenna pattern to obtain an incident direction of a radio wave. A radio azimuth detecting device characterized by the above-mentioned.
【請求項7】 電波を受信するアレーアンテナと、 前記アレーアンテナの受信信号の相関処理を行い相関行
列を計算し、前記相関行列の固有値を計算し、前記固有
値の大きさから入射波数を推定する入射波数推定手段
と、 前記入射波数推定手段により推定された入射波数に応じ
て電波の入射方向を求める信号処理装置を切り替える信
号処理装置切替手段と、 前記入射波数が1である場合は、前記固有値のうち最も
大きな固有値に対応する固有ベクトルを算出し、前記固
有ベクトルをステアリングベクトルとして入射方向を求
める第1の信号処理装置と、 前記入射波数が1より大きい場合には、ノイズに対応し
た固有値に対応する固有ベクトルを求め、前記アレーア
ンテナの各方位に対する位相特性や振幅特性を要素とし
て持つステアリングベクトルと前記固有ベクトルから方
位評価関数を算出し、前記方位評価関数のピークを検出
して入射方向を求める第2の信号処理装置とを備えたこ
とを特徴とする電波方位探知装置。
7. An array antenna for receiving a radio wave, and a correlation matrix is calculated by performing a correlation process on a signal received by the array antenna, an eigenvalue of the correlation matrix is calculated, and the number of incident waves is estimated from the magnitude of the eigenvalue. An incident wave number estimating unit; a signal processing device switching unit that switches a signal processing device that determines an incident direction of a radio wave according to the incident wave number estimated by the incident wave number estimating unit; and the eigenvalue when the incident wave number is 1. A first signal processing device that calculates an eigenvector corresponding to the largest eigenvalue and determines an incident direction using the eigenvector as a steering vector, and, when the number of incident waves is greater than 1, corresponds to an eigenvalue corresponding to noise. An eigenvector is obtained, and the steering having the phase characteristic and the amplitude characteristic for each direction of the array antenna as an element. Calculating the orientation evaluation function from vector and the eigenvector, wave direction detection apparatus characterized by comprising a second signal processing apparatus detect and determine the direction of incidence of the peak of the orientation evaluation function.
JP30554797A 1997-11-07 1997-11-07 Radio direction detector Expired - Fee Related JP3647621B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP30554797A JP3647621B2 (en) 1997-11-07 1997-11-07 Radio direction detector

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP30554797A JP3647621B2 (en) 1997-11-07 1997-11-07 Radio direction detector

Publications (2)

Publication Number Publication Date
JPH11142497A true JPH11142497A (en) 1999-05-28
JP3647621B2 JP3647621B2 (en) 2005-05-18

Family

ID=17946482

Family Applications (1)

Application Number Title Priority Date Filing Date
JP30554797A Expired - Fee Related JP3647621B2 (en) 1997-11-07 1997-11-07 Radio direction detector

Country Status (1)

Country Link
JP (1) JP3647621B2 (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001029989A1 (en) * 1999-10-21 2001-04-26 Matsushita Electric Industrial Co., Ltd. Array-antenna wireless communication device, and method of generating weight coefficient
US6573859B2 (en) 2000-02-07 2003-06-03 Toyota Jidosha Kabushiki Kaisha Radar apparatus
WO2004104620A1 (en) * 2003-05-22 2004-12-02 Fujitsu Limited Technique for estimating signal arrival direction not by utilizing eigenvalue decomposition and reception beam shaper
JP2005121581A (en) * 2003-10-20 2005-05-12 Mitsubishi Electric Corp Radar device
JP2009036539A (en) * 2007-07-31 2009-02-19 Mitsubishi Electric Corp Radar signal processing apparatus and method
JP2009204501A (en) * 2008-02-28 2009-09-10 Mitsubishi Electric Corp Wavenumber estimating device
JP2010237087A (en) * 2009-03-31 2010-10-21 Hitachi Automotive Systems Ltd Radar apparatus and method for measuring radio wave arrival direction using the same
JP2012088236A (en) * 2010-10-21 2012-05-10 Honda Elesys Co Ltd Electron scanning type radar device, received wave direction estimation method, and received wave direction estimation program
JP2014044160A (en) * 2012-08-28 2014-03-13 Mitsubishi Electric Corp Positioning device

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001029989A1 (en) * 1999-10-21 2001-04-26 Matsushita Electric Industrial Co., Ltd. Array-antenna wireless communication device, and method of generating weight coefficient
US7110733B1 (en) 1999-10-21 2006-09-19 Matsushita Electric Industrial Co., Ltd. Array antenna radio communication apparatus and weight coefficient generating method
US6573859B2 (en) 2000-02-07 2003-06-03 Toyota Jidosha Kabushiki Kaisha Radar apparatus
WO2004104620A1 (en) * 2003-05-22 2004-12-02 Fujitsu Limited Technique for estimating signal arrival direction not by utilizing eigenvalue decomposition and reception beam shaper
US7068221B2 (en) 2003-05-22 2006-06-27 Fujitsu Limited Technique for direction-of-arrival estimation without eigendecomposition and its application to beamforming at base station
JP2005121581A (en) * 2003-10-20 2005-05-12 Mitsubishi Electric Corp Radar device
JP2009036539A (en) * 2007-07-31 2009-02-19 Mitsubishi Electric Corp Radar signal processing apparatus and method
JP2009204501A (en) * 2008-02-28 2009-09-10 Mitsubishi Electric Corp Wavenumber estimating device
JP2010237087A (en) * 2009-03-31 2010-10-21 Hitachi Automotive Systems Ltd Radar apparatus and method for measuring radio wave arrival direction using the same
JP2012088236A (en) * 2010-10-21 2012-05-10 Honda Elesys Co Ltd Electron scanning type radar device, received wave direction estimation method, and received wave direction estimation program
JP2014044160A (en) * 2012-08-28 2014-03-13 Mitsubishi Electric Corp Positioning device

Also Published As

Publication number Publication date
JP3647621B2 (en) 2005-05-18

Similar Documents

Publication Publication Date Title
JP2001305202A (en) MUSIC spectrum calculation method, apparatus and medium
US20060208947A1 (en) Apparatus and method for estimating direction of arrival of radio wave
JP2988463B2 (en) Direction finding device and measurement result processing device therefor
JPH08262121A (en) Position measuring instrument
US20190280740A1 (en) Wireless receiver, wireless receiving method, and wireless system
JP2000121716A (en) Radio wave propagation estimating equipment
KR20170027375A (en) Method of spatial interpolation for linear phased array antenna and appratus thereof
JPH11142497A (en) Wave bearing detector
JP3629195B2 (en) Arrival direction estimation apparatus and arrival direction estimation method
JPH08201498A (en) Angle measuring apparatus
JP4187985B2 (en) Angle measuring device, angle measuring method and program
JPH11166965A (en) Angle measurement signal processing device and angle measurement signal processing method in it
JP4977849B2 (en) Radio wave arrival direction detector
JP6573745B2 (en) Adaptive array antenna device
JP2948459B2 (en) Angle measuring device
Godara Limitations and capabilities of directions-of-arrival estimation techniques using an array of antennas: a mobile communications perspective
Featherstone et al. A novel method to improve the performance of Capon's minimum variance estimator
JP2008008781A (en) Device of detecting the number of incoming waves
WO2022185379A1 (en) Angle width estimation device, angle width estimation method, and target angle measurement device
JP4119719B2 (en) Mobile station direction estimation method and apparatus
Chargé et al. A direction finding method under sensor gain and phase uncertainties
JP2002107440A (en) Method and device for detecting direction
JP3817875B2 (en) Direction finding device
CN114690114A (en) Signal direction finding method, device, system and storage medium
JP7585870B2 (en) Adaptive phasing processing device, adaptive phasing system, adaptive phasing processing method, and program for adaptive phasing processing method

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20040907

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: 20050208

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050209

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080218

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090218

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100218

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100218

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110218

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120218

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130218

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130218

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140218

Year of fee payment: 9

LAPS Cancellation because of no payment of annual fees