JP3817875B2 - Direction finding device - Google Patents

Direction finding device Download PDF

Info

Publication number
JP3817875B2
JP3817875B2 JP33693997A JP33693997A JP3817875B2 JP 3817875 B2 JP3817875 B2 JP 3817875B2 JP 33693997 A JP33693997 A JP 33693997A JP 33693997 A JP33693997 A JP 33693997A JP 3817875 B2 JP3817875 B2 JP 3817875B2
Authority
JP
Japan
Prior art keywords
antenna
calculator
arrival direction
evaluation function
antenna elements
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
Application number
JP33693997A
Other languages
Japanese (ja)
Other versions
JPH11166966A (en
Inventor
暁 村山
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 JP33693997A priority Critical patent/JP3817875B2/en
Publication of JPH11166966A publication Critical patent/JPH11166966A/en
Application granted granted Critical
Publication of JP3817875B2 publication Critical patent/JP3817875B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Description

【0001】
【発明の属する技術分野】
この発明は、複数の電波源から到来する少なくとも2以上の複数電波が同一周波数帯にある場合に、それぞれの到来方向を同時に測定し、電波の到来方向を推定する方向探知装置に関する。
【0002】
【従来の技術】
図11は従来の方向探知装置を示すもので、この図はSchmidt著、"Multiple Emitter Location and Signal Parameter Estimation",IEEE Trance,AP−34,3,pp.276−280(1986).に記述されている内容に基づいた方向探知装置の構成図である。
【0003】
また、図12は前記方向探知装置の運用図である。
【0004】
図11において、1−1、1−2、・・・、1−Mは電波源からの信号を受信するアンテナ素子、2は1−1、1−2、・・・、1−Mのアンテナ素子を取り付けたアンテナタワー、3は従来の信号処理を行う信号処理器A、4−1、4−2、・・・、4−Mは受信機、5−1、5−2、・・・、5−Mは受信機4の出力信号を2分配し、互いの位相差がπ/2ずれるように位相検波するディジタル位相検波器、6は受信信号共分散行列算出器、7は雑音部分空間行列算出器、8はステアリングベクトル算出器、9は到来方向評価関数算出器、10は到来方向検出器である。
【0005】
図12において、11は電波源の一例として無線送信機を持つ自動車を示している。
【0006】
次に動作について説明する。ここで用いる記号を以下に定義する。Kは到来波数、kkは到来波番号(kk≦K)、Mはアンテナ素子数、mはアンテナ素子番号(1≦m≦M)であり、K<Mとする。また、前記M個のアンテナ素子同士は、任意に配列されている。
【0007】
図11において、K個の電波源からアンテナ素子1−mに到来したK個の到来波は、受信機4−mに入力され、受信機4−mからの出力信号が前記ディジタル位相検波器5−mでA/D変換後、位相検波され、受信信号共分散行列算出器6において、ディジタル位相検波器5−1、5−2、・・・、5−Mの出力であるM個のディジタル信号を入力し受信信号共分散行列を求める。ここで、ディジタル位相検波器5−1、5−2、・・・、5−Mの出力であるM個のディジタル信号を"数1"に、受信信号共分散行列を"数2"に示す。
【0008】
【数1】

Figure 0003817875
【0009】
【数2】
Figure 0003817875
【0010】
次に、雑音部分空間行列算出器7において、受信信号共分散行列算出器6の出力である受信信号共分散行列の固有値・固有ベクトルを雑音固有値・雑音固有ベクトルとして求め、雑音固有ベクトルより雑音部分空間行列を算出する。ここで、受信信号共分散行列の固有値・固有ベクトルから求められる雑音部分空間行列を“数3”に示す。
【0011】
【数3】
Figure 0003817875
【0012】
またステアリングベクトル算出器8において、前記M個のアンテナ素子の配列から求められるアレー応答に対応したステアリングベクトルを算出し、これを“数4”に示す。
【0013】
【数4】
Figure 0003817875
【0014】
到来方向評価関数算出器9において、雑音部分空間行列算出器7の出力である雑音部分空間行列とステアリングベクトル算出器8の出力であるステアリングベクトルから"数5"に示す到来方向評価関数P(θ)を算出し、到来方向検出器10において、到来方向評価関数算出器9の出力である到来方向評価関数P(θ)を角度θについてサーチし、図13に示すようにP(θ)の最大ピークからK個を選ぶことによりそのピークに対応した角度を検出し、到来波の到来方向(θ1、θ2、・・・、θK)とする。これによって、ピークに対応した角度を検出し擬像を抑圧する。但し、図13では例として到来波数を2波としている。
【0015】
【数5】
Figure 0003817875
【0016】
【発明が解決しようとする課題】
ところが、上記のような従来の方向探知装置では、アンテナ素子数と同数の受信機、ディジタル位相検波器が必要で、製品の重量が重くなり、また受信機が高価なものであるため、製品価格が高価になりやすいという問題点があった。
【0017】
この発明は、上記のような問題点を解決するものであり、アンテナ素子の総数に関わらず選択された3個以上のアンテナ素子と、前記選択された3個以上のアンテナ素子に接続された3個以上の受信機、ディジタル位相検波器で信号処理を行うことにより粗到来方向評価関数を求め、全アンテナ素子からJ個のアンテナ素子を選択する全組み合わせそれぞれから求められる粗到来方向評価関数を積分処理することにより精到来方向評価関数を求め、到来方向を検出することを目的としている。
【0018】
【課題を解決するための手段】
この発明の第1の発明による方向探知装置は、電波源から到来する少なくとも2以上の電波を受信する少なくとも4個以上の複数のアンテナ素子と、前記複数のアンテナ素子のうち任意の3個以上のアンテナ素子を選択するアンテナ切替器と、前記アンテナ切替器で選択されるアンテナ素子の選択の全組み合わせを制御する切替制御器と、前記アンテナ切替器で選択されるそれぞれのアンテナ素子端に接続される3個以上の受信機と、前記それぞれの受信機によって供給される信号から受信信号共分散行列を算出する受信信号共分散行列算出器と、前記受信信号共分散行列算出器より供給される受信信号共分散行列から固有値/固有ベクトルを算出し、前記固有ベクトルから構成される雑音部分空間行列を算出する雑音部分空間行列算出器と、前記アンテナ切替器で選択されたアンテナ素子間の配列に基づいてステアリングベクトルを算出するステアリングベクトル算出器と、前記雑音部分空間行列算出器から供給される雑音部分空間行列と前記ステアリングベクトル算出器から供給されるステアリングベクトルに基づいて、粗到来方向評価関数を算出する粗到来方向評価関数算出器と、前記アンテナ切替器で選択されるアンテナ素子の全組み合わせに対応して前記粗到来方向評価関数算出器から供給されるそれぞれの全粗到来方向評価関数を積分し、精到来方向評価関数を算出する精到来方向評価関数算出器と、前記精到来方向評価関数算出器から供給される精到来方向評価関数から電波の到来方向を検出する到来方向検出器とを備えたものである。
【0019】
また、第2の発明による方向探知装置は、電波源から到来する少なくとも2以上の電波を受信する少なくとも4個以上の複数のアンテナ素子と、前記複数のアンテナ素子のうち任意の3個以上のアンテナ素子を選択するアンテナ切替器と、前記アンテナ切替器で選択されるアンテナ素子の選択の組み合わせを制御する切替制御器と、前記アンテナ切替器で選択されるそれぞれのアンテナ素子端に接続される3個以上の受信機と、前記それぞれの受信機によって供給される信号から受信信号共分散行列を算出する受信信号共分散行列算出器と、前記受信信号共分散行列算出器より供給される受信信号共分散行列から固有値/固有ベクトルを算出し、前記固有ベクトルから構成される雑音部分空間行列を算出する雑音部分空間行列算出器と、前記アンテナ切替器で選択されたアンテナ素子間の配列に基づいてステアリングベクトルを算出するステアリングベクトル算出器と、前記雑音部分空間行列算出器から供給される雑音部分空間行列と前記ステアリングベクトル算出器から供給されるステアリングベクトルに基づいて、粗到来方向評価関数を算出する粗到来方向評価関数算出器と、前記アンテナ切替器で選択されるアンテナ素子の全組み合わせに対応して前記粗到来方向評価関数算出器から供給されるそれぞれの粗到来方向評価関数から電波の到来方向を検出する粗到来方向検出器と、前記粗到来方向検出器から供給される全組み合わせ分の粗到来方向を平均し、精到来方向を算出する到来方向平均処理器とを備えた。
【0020】
また、第3の発明による方向探知装置は、第1の発明、もしくは第2の発明によるアンテナ切替器に、前記複数のアンテナ素子のうちの1つのアンテナ素子を直接前記受信機の一つに接続したまま、他の複数のアンテナ素子のうちの任意の2個以上のアンテナ素子を選択する手段を備えたものである。
【0021】
また、第4の発明による方向探知装置は、第1の発明、もしくは第2の発明の前記切替制御器において、それぞれのアンテナ素子を1回以上選択するよう制御する手段を備えたものである。
【0022】
また、第5の発明による方向探知装置は、第3の発明の前記切替制御器において、前記他の複数のアンテナ素子を1回以上選択するよう制御する手段を備えたものである。
【0023】
【発明の実施の形態】
実施の形態1.
図1はこの発明の実施の形態1を示す構成図である。従来例を示す図11では、アンテナ素子数と同数の受信機、ディジタル位相検波器を用いて信号処理を行い電波の到来方向を検出していたが、この実施の形態では価格低減のため、全アンテナ素子のうち全アンテナ素子数よりも少なく到来波の数より多い、例えば3個以上のアンテナ素子を選択するアンテナ切替器及びアンテナ切替器を制御する切替制御器を設け、3個以上の受信機、ディジタル位相検波器を用いて粗到来方向評価関数を算出する。この粗到来方向評価関数を用いた場合、真の到来方向以外に擬像の発生する場合があるが、前記アンテナ切替器によって選択されるそれぞれ異なる3個以上のアンテナ素子の全組み合わせから算出されるそれぞれの粗到来方向評価関数を積分処理し精到来方向評価関数を求めることによって擬像を抑圧し、従来例と同等の精度で電波源から到来する複数電波の到来方向を検出できる。このため、従来と比べて受信機の数を少なくでき、軽量化および低価格な方向探知装置を得ることが出来る。
【0024】
図1において、1−1、1−2、・・・、1−M(M≧4)は電波源からの信号を受信するアンテナ素子、2はアンテナ素子1−1、1−2、・・・、1−M(以下アンテナ素子1)を取り付けたアンテナタワー、12はアンテナタワー1から供給される信号の信号処理を行う信号処理装置B、13はアンテナタワー2における各アンテナ素子1との接続を選択するアンテナ切替器A、14はアンテナ切替器A13における接続を順次切替えて全組み合わせを制御する切替制御器A、4−1、4−2、・・・、4−Jはアンテナ切替器A13の切替に応じてアンテナ1との接続が切替わる受信機、5−1、5−2、・・・、5−Jは受信機4の出力信号を2分配し、互いの位相差がπ/2ずれるように位相検波するディジタル位相検波器、6は受信信号共分散行列算出器、7は雑音部分空間行列算出器、8はステアリングベクトル算出器、15は到来方向評価関数算出器、16は精到来方向評価関数算出器、10は到来方向検出器で、信号処理装置B12は4〜8、10、および13〜16によって信号処理を行う。
【0025】
次に図1の動作を説明する。ここで用いる記号を以下に定義する。Kは到来波個数、kkは到来波番号(kk≦K)、Mはアンテナ素子1の個数、mはアンテナ素子番号(1≦m≦M)、Jは選択されたアンテナ素子個数、jjは選択されたJ個のアンテナ素子のうちjj番目(1≦jj≦J)のアンテナ素子番号であり、K<J<Mかつ3≦Jである。前記M個のそれぞれのアンテナ素子同士は任意に配列されている。また、電波源は例えば数10km先の遠方を走行中の自動車、船などの移動体から発生されるものを考慮しており、測角処理を行っている間の処理時間内において電波源の移動はないと近似されるものとする。
【0026】
図1において、K個の電波源からアンテナ素子1−mに受信された受信信号のうち、アンテナ切替器A13で選択されたJ個のアンテナ素子1のうちアンテナ素子1−jjからの出力が、受信機4−jjに入力され、受信機4−jjからの出力信号が前記ディジタル位相検波器5−jjでA/D変換後位相検波され、受信信号共分散行列算出器6において、J個のアンテナ素子1のそれぞれに対応して検波されるディジタル位相検波器5−1、5−2、・・・、5−Jの出力であるJ個のディジタル信号が入力され、受信信号共分散行列が求められる。ここで、切替制御器A14で制御され、アンテナ切替器A13で切り換えられる選択されたアンテナ素子J個の全組み合わせ数をYとすると、“数6”のように示される。前記ディジタル位相検波器5−1、5−2、・・・、5−Jの出力であるY通りあるアンテナ組み合わせのy通り目のJ個のディジタル信号を“数7”に、受信信号共分散行列を“数8”に示す。
【0027】
【数6】
Figure 0003817875
【0028】
【数7】
Figure 0003817875
【0029】
【数8】
Figure 0003817875
【0030】
次に、雑音部分空間行列算出器7において、受信信号共分散行列算出器6の出力である受信信号共分散行列の固有値・固有ベクトルを雑音固有値、雑音固有ベクトルとして求め、雑音固有ベクトルより雑音部分空間行列を算出し、これを“数9”に示す。
【0031】
【数9】
Figure 0003817875
【0032】
またステアリングベクトル算出器8において、前記選択されたJ個のアンテナ素子の配列から求められるアレー応答に対応したJ列のステアリングベクトルを算出し、これを“数10”に示す。
【0033】
【数10】
Figure 0003817875
【0034】
粗到来方向評価関数算出器15において、雑音部分空間行列算出器7の出力である雑音部分空間行列とステアリングベクトル算出器8の出力であるステアリングベクトルから粗到来方向評価関数を算出するが、ここで、y番目の粗到来方向評価関数Py (θ)(y=1...Y)は“数11”に示される。但し、前記アンテナ切替器A13におけるアンテナの切り替え方(選択の仕方)は任意である。
【0035】
【数11】
Figure 0003817875
【0036】
次に、精到来方向評価関数算出器16において、切替制御器A14で制御されたアンテナ切替器A13で選択されたJ個のアンテナ素子の全組み合わせ分の粗到来方向評価関数Py (θ)を積分処理し精到来方向評価関数AP(θ)を算出する。ここで、精到来方向評価関数AP(θ)は“数12”に示される。1例として、図2に粗到来方向評価関数Py (θ)と、それを積分処理して求まる精到来方向評価関数AP(θ)を示す。これを求めることにより、受信機の数を減らしたことによって生じる擬像を抑圧することができるとともに、S/Nを改善出来る。
【0037】
【数12】
Figure 0003817875
【0038】
最後に、前記10の到来方向検出器において、精到来方向評価関数算出器16の出力である精到来方向評価関数AP(θ)からK個の到来波の到来方向(ψ1 、ψ2 、・・・、ψK )を検出する。この検出は従来の図11と同様、AP(θ)の最大ピークからK個を選ぶことによりそのピークに対応した角度を検出し、到来波の到来方向(ψ1 、ψ2 、・・・、ψK )を求める。
【0039】
実施の形態2.
図3はこの発明の実施の形態2を示す構成図である。実施の形態1では、アンテナ切替器13においてM個のアンテナ素子のうちJ個のアンテナ素子を選択したが、この実施の形態では、図3のように1個のアンテナ素子を1個の受信機に接続し、残りの(M−1)個のアンテナ素子のうちの任意の(J−1)個のアンテナ素子を選択するアンテナ切替器B17と、アンテナ切替器B17の全組み合わせを制御する切替制御器B18を用いた信号処理装置C19を用いることにより、実施の形態1と同様な効果を得ることができ、かつ、電波源の信号を聴音する場合に、復調時のアンテナ切り替えノイズを取り除く装置を付加せずに、常にアンテナ素子に固定されている受信機から電波源の信号を聴音することが可能となる。
【0040】
実施の形態3.
図4はこの発明の実施の形態3を示す構成図である。実施の形態1では、信号処理装置B12においてY個の粗到来方向評価関数Py (θ)(y=1...Y)から、積分処理によって精到来方向評価関数AP(θ)を求め、精到来方向(ψ1 、ψ2 、・・・、ψK )を検出した。この実施の形態では、図4のように20の信号処理装置Dを設け、精到来方向評価関数算出器16の代わりに、粗到来方向評価関数算出器15の出力である粗到来方向評価関数Py (θ)(y=1...Y)が入力され、粗到来方向を検出する21の粗到来方向検出器と、到来方向検出器10の代わりに、粗到来方向検出器21から出力されるアンテナ切替器13における全切替の組み合わせ分Y個の粗到来方向(θ(y) 1、θ(y) 2、・・・、θ(y) K)(y=1...Y)について、値の近いもの同士が同じ集団になるように分け、数の多い集団からK個の集団を選び、それぞれ平均処理することによって精到来方向(ψ1 、ψ2 、・・・、ψK )を算出する22の到来方向平均処理器を用いて実施の形態1と同様な効果を得ることができる。
【0041】
実施の形態4.
図5はこの発明の実施の形態4を示す構成図である。実施の形態3では、アンテナ切替器13においてM個のアンテナ素子のうちJ個のアンテナ素子を選択したが、この実施の形態では、図5のように1個のアンテナ素子を1個の受信機に接続し、残りの(M−1)個のアンテナ素子のうちの任意の(J−1)個の素子アンテナ素子を選択するアンテナ切替器B17と、アンテナ切替器B17の全組み合わせを制御する切替制御器B18を用いた信号処理装置D23を用いることにより、実施の形態3と同様な効果を得ることができ、かつ、電波源の信号を聴音する場合に、復調時のアンテナ切り替えノイズを取り除く装置を付加せずに、常にアンテナ素子に固定されている受信機から電波源の信号を聴音することが可能となる。
【0042】
実施の形態5.
図6はこの発明の実施の形態5を示す構成図である。実施の形態1では、切替制御器A14においてM個のアンテナ素子からJ個のアンテナ素子を選択する全組み合わせY個の制御を行ったが、アンテナの組み合わせによっては全体のアンテナ開口径が小さくなり測角精度が悪くなる。この実施の形態では、図6のようにM個のアンテナ素子からJ個のアンテナ素子を選択する組み合わせ数が全組み合わせ数でなく、それぞれのアンテナ素子を1回以上とり、かつ全体のアンテナ開口径が大きくとれる組み合わせ分選択し制御する24の切替制御器Cを用いて測角信号処理を行う25の信号処理装置Fを用いることにより、切替制御器A14の処理がより簡単になり、処理時間が早くなるとともに実施の形態1と同様な効果を得ることができる。
【0043】
実施の形態6.
図7はこの発明の実施の形態6を示す構成図である。実施の形態2では、切替制御器B18において(M−1)個のアンテナ素子から(J−1)個のアンテナ素子を選択する全組み合わせの制御を行ったが、アンテナの組み合わせによっては全体のアンテナ開口径が小さくなり測角精度が悪くなる。この実施の形態6は、図7のように(M−1)個のアンテナ素子から(J−1)個のアンテナ素子を選択する組み合わせ数が全組み合わせ数でなくても、それぞれのアンテナ素子を1回以上とり、かつ全体のアンテナ開口径が大きくとれる組み合わせ分選択し制御する切替制御器D26を用いて測角信号処理を行う信号処理装置G27を用いることにより、実施の形態2と同様な効果を得ることができる。1例としてアンテナ素子数7個、受信機3個の場合のアンテナ選択の組み合わせを図8に示す。
【0044】
実施の形態7.
図9はこの発明の実施の形態7を示す構成図である。実施の形態3では、切替制御器A14においてM個のアンテナ素子からJ個のアンテナ素子を選択する全組み合わせY個の制御を行ったが、アンテナの組み合わせによっては全体のアンテナ開口径が小さくなり測角精度が悪くなる。この実施の形態では、図9のようにM個のアンテナ素子からJ個のアンテナ素子を選択する組み合わせ数が全組み合わせ数でなくても、それぞれのアンテナ素子を1回以上とり、かつ全体のアンテナ開口径が大きくとれる組み合わせ分選択し制御する切替制御器C24を用いて測角信号処理を行う信号処理装置H28を用いることにより、実施の形態3と同様な効果を得ることができる。
【0045】
実施の形態8.
図10はこの発明の実施の形態8を示す構成図である。実施の形態4では、切替制御器B18において(M−1)個のアンテナ素子から(J−1)個のアンテナ素子を選択する全組み合わせの制御を行ったが、アンテナの組み合わせによっては全体のアンテナ開口径が小さくなり測角精度が悪くなる。この実施の形態では、図10のように(M−1)個のアンテナ素子から(J−1)個のアンテナ素子を選択する組み合わせ数が全組み合わせ数でなくても、それぞれのアンテナ素子を1回以上とり、かつ全体のアンテナ開口径が大きくとれる組み合わせ分選択し制御する切替制御器D26用いて測角信号処理を行う信号処理装置I29を用いても実施の形態2と同様な効果を得ることができる。
【0046】
【発明の効果】
この発明に係る第1から第5の発明によれば、アンテナ素子数と同数の受信機を用いる方向探知装置と同程度の検知精度を保持したまま、受信信号を処理する受信器の数の少ない、軽量で低価格な方向探知装置を得ることができる。
【0047】
また、第3、第5の発明によれば、M個のうちの1つのアンテナ素子を常に受信機の一つに接続し、他の(M−1)個のアンテナ素子のうちの(J−1)個のアンテナ素子を選択するアンテナ切替器を備えることにより、第1の発明と同様な効果を得ることができるとともに、アンテナ素子の切替時に常に一つのアンテナ素子が一つの受信機に接続されているため、電波源の信号を聴音する場合、復調時におけるアンテナ切り替え雑音を取り除く装置を新たに付加せずに、直接その受信機から復調し聴音することができる。
【0048】
さらにまた、第4、第5の発明によれば、アンテナ切替器におけるアンテナ素子の組み合わせ数を全て用いず、それぞれのアンテナ素子を1回以上とり、かつ全体のアンテナ開口径が大きくとれるような組み合わせ分で測角信号処理を行うことにより、第1から第3の発明と比べて演算処理時間を低減できる。
【図面の簡単な説明】
【図1】 この発明による方向探知装置の実施の形態1の一実施例を示す図である。
【図2】 積分処理によって偽像を抑圧していることを示す図である。
【図3】 この発明による方向探知装置の実施の形態2の一実施例を示す図である。
【図4】 この発明による方向探知装置の実施の形態3の一実施例を示す図である。
【図5】 この発明による方向探知装置の実施の形態4の一実施例を示す図である。
【図6】 この発明による方向探知装置の実施の形態5の一実施例を示す図である。
【図7】 この発明による方向探知装置の実施の形態6の一実施例を示す図である。
【図8】 アンテナの組み合わせの一例を示す図である。
【図9】 この発明による方向探知装置の実施の形態7の一実施例を示す図である。
【図10】 この発明による方向探知装置の実施の形態8の一実施例を示す図である。
【図11】 従来の方向探知装置の一実施例を示す図である。
【図12】 運用図である。
【図13】 到来方向検出器における検出方法を示した図である。
【符号の説明】
1 アンテナ素子、2 アンテナタワー、3 信号処理器A、4 受信機、5ディジタル位相検波器、6 受信信号共分散行列算出器、7 雑音部分空間行列算出器、8 ステアリングベクトル算出器、9 到来方向評価慣習算出器、10 到来方向検出器、11 自動車、12 信号処理装置B、13 アンテナ切替器A、14 切替制御器A、15 粗到来方向評価関数算出器、16 精到来方向評価関数算出器、17 アンテナ切替器B、18 切替制御器B、19 信号処理装置C、20 信号処理装置D、21 粗到来方向検出器、22 到来方向平均処理器、23 信号処理装置E、24 切替制御器C、25 信号処理装置F、26 切替制御器D、27 信号処理装置G、28 信号処理装置H、29 信号処理装置I。[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a direction detection apparatus that simultaneously measures the arrival directions of each radio wave and estimates the radio wave arrival directions when at least two or more radio waves coming from a plurality of radio wave sources are in the same frequency band.
[0002]
[Prior art]
FIG. 11 shows a conventional direction finding device, which is illustrated by Schmidt, “Multiple Emitter Location and Signal Parameter Estimation”, IEEE Transform, AP-34, 3, pp. 276-280 (1986). It is a block diagram of the direction detection apparatus based on the content described in (1).
[0003]
FIG. 12 is an operational diagram of the direction finding device.
[0004]
In FIG. 11 , 1-1, 1-2,..., 1-M are antenna elements that receive a signal from a radio wave source, 2 is an antenna of 1-1, 1-2,. Antenna tower to which elements are attached, 3 is a signal processor A that performs conventional signal processing, 4-1, 4-2,..., 4-M is a receiver, 5-1, 5-2,. , 5-M distributes the output signal of the receiver 4 in two and detects the phase so that the mutual phase difference is shifted by π / 2, 6 is a received signal covariance matrix calculator, and 7 is a noise subspace. A matrix calculator, 8 is a steering vector calculator, 9 is an arrival direction evaluation function calculator, and 10 is an arrival direction detector.
[0005]
In FIG. 12 , reference numeral 11 denotes an automobile having a wireless transmitter as an example of a radio wave source.
[0006]
Next, the operation will be described. The symbols used here are defined below. K is the number of incoming waves, kk is the incoming wave number (kk ≦ K), M is the number of antenna elements, m is the antenna element number (1 ≦ m ≦ M), and K <M. The M antenna elements are arbitrarily arranged.
[0007]
In FIG. 11 , K incoming waves arriving at the antenna element 1-m from K radio wave sources are input to the receiver 4-m, and an output signal from the receiver 4-m is the digital phase detector 5. -M, after A / D conversion, phase detection is performed, and the received signal covariance matrix calculator 6 outputs M digital outputs which are outputs of digital phase detectors 5-1, 5-2, ..., 5-M. Input a signal and obtain a received signal covariance matrix. Here, M digital signals which are outputs of the digital phase detectors 5-1, 5-2,... .
[0008]
[Expression 1]
Figure 0003817875
[0009]
[Expression 2]
Figure 0003817875
[0010]
Next, in the noise subspace matrix calculator 7, the eigenvalue / eigenvector of the received signal covariance matrix that is the output of the received signal covariance matrix calculator 6 is obtained as a noise eigenvalue / noise eigenvector, and the noise subspace matrix is obtained from the noise eigenvector. calculate. Here, the noise subspace matrix obtained from the eigenvalues and eigenvectors of the received signal covariance matrix is shown in “Expression 3”.
[0011]
[Equation 3]
Figure 0003817875
[0012]
Further, the steering vector calculator 8 calculates a steering vector corresponding to the array response obtained from the arrangement of the M antenna elements.
[0013]
[Expression 4]
Figure 0003817875
[0014]
In the direction-of-arrival evaluation function calculator 9, the direction-of-arrival evaluation function P (θ shown in “Expression 5” is calculated from the noise subspace matrix output from the noise subspace matrix calculator 7 and the steering vector output from the steering vector calculator 8. ) And the arrival direction evaluation function P (θ), which is the output of the arrival direction evaluation function calculation unit 9, is searched for the angle θ, and the maximum of P (θ) is obtained as shown in FIG. By selecting K pieces from the peaks, the angle corresponding to the peaks is detected and set as the arrival directions (θ1, θ2,..., ΘK) of the incoming waves. This detects the angle corresponding to the peak and suppresses the pseudo image. However, in FIG. 13 , the number of incoming waves is two as an example.
[0015]
[Equation 5]
Figure 0003817875
[0016]
[Problems to be solved by the invention]
However, the conventional direction finding device as described above requires the same number of receivers and digital phase detectors as the number of antenna elements, which increases the weight of the product and makes the receiver expensive. However, there was a problem that it was easy to become expensive.
[0017]
The present invention solves the above-described problems, and the three or more selected antenna elements regardless of the total number of antenna elements, and the three connected to the selected three or more antenna elements. The coarse arrival direction evaluation function is obtained by performing signal processing with more than one receiver and digital phase detector, and the coarse arrival direction evaluation function obtained from all combinations of selecting J antenna elements from all antenna elements is integrated. The purpose is to obtain a precision arrival direction evaluation function by processing and detect the arrival direction.
[0018]
[Means for Solving the Problems]
A direction detecting apparatus according to a first aspect of the present invention includes at least four or more antenna elements that receive at least two or more radio waves coming from a radio wave source, and any three or more of the plurality of antenna elements. An antenna switch for selecting an antenna element, a switching controller for controlling all combinations of antenna element selections selected by the antenna switch, and an antenna element end selected by the antenna switch. Three or more receivers, a received signal covariance matrix calculator for calculating a received signal covariance matrix from signals supplied by the respective receivers, and a received signal supplied from the received signal covariance matrix calculator A noise subspace matrix calculator for calculating eigenvalues / eigenvectors from a covariance matrix and calculating a noise subspace matrix composed of the eigenvectors; A steering vector calculator that calculates a steering vector based on an array between antenna elements selected by the antenna switch, a noise subspace matrix supplied from the noise subspace matrix calculator, and a steering vector calculator A rough arrival direction evaluation function calculator for calculating a rough arrival direction evaluation function based on the steering vector, and the rough arrival direction evaluation function calculator corresponding to all combinations of antenna elements selected by the antenna switcher Integrating each coarse arrival direction evaluation function supplied from each of them, and calculating a fine arrival direction evaluation function, and a fine arrival direction evaluation function supplied from the fine arrival direction evaluation function calculator And a direction-of-arrival detector for detecting the direction of arrival of radio waves.
[0019]
In addition, a direction finding device according to a second aspect of the present invention includes at least four or more antenna elements that receive at least two or more radio waves coming from a radio wave source, and any three or more antennas among the plurality of antenna elements. An antenna switch for selecting an element, a switch controller for controlling a combination of antenna element selections selected by the antenna switch, and three pieces connected to respective antenna element ends selected by the antenna switch The above receiver, a received signal covariance matrix calculator for calculating a received signal covariance matrix from signals supplied by the respective receivers, and a received signal covariance supplied from the received signal covariance matrix calculator A noise subspace matrix calculator for calculating an eigenvalue / eigenvector from a matrix and calculating a noise subspace matrix composed of the eigenvector; A steering vector calculator for calculating a steering vector based on an arrangement between antenna elements selected by the antenna switcher, a noise subspace matrix supplied from the noise subspace matrix calculator, and a steering vector calculator From the rough arrival direction evaluation function calculator for calculating the rough arrival direction evaluation function based on the steering vector, and from the rough arrival direction evaluation function calculator corresponding to all combinations of the antenna elements selected by the antenna switch The coarse arrival direction detector for detecting the arrival direction of the radio wave from each of the supplied coarse arrival direction evaluation functions and the coarse arrival directions for all combinations supplied from the coarse arrival direction detector are averaged to obtain the fine arrival direction. And a direction-of-arrival average processor for calculation.
[0020]
A direction finding device according to a third aspect of the present invention is the antenna switcher according to the first or second aspect, wherein one antenna element of the plurality of antenna elements is directly connected to one of the receivers. As it is, there is provided means for selecting any two or more of the other antenna elements.
[0021]
According to a fourth aspect of the present invention, there is provided a direction finder according to the first or second aspect of the present invention, further comprising means for controlling each antenna element to be selected one or more times.
[0022]
According to a fifth aspect of the present invention, there is provided a direction finder according to the third aspect of the present invention, further comprising means for controlling the plurality of other antenna elements to be selected one or more times in the switching controller.
[0023]
DETAILED DESCRIPTION OF THE INVENTION
Embodiment 1 FIG.
FIG. 1 is a block diagram showing Embodiment 1 of the present invention. In FIG. 11 showing the conventional example, signal processing is performed by using the same number of receivers and digital phase detectors as the number of antenna elements to detect the direction of arrival of radio waves. Provided with an antenna switch for selecting three or more antenna elements and a switching controller for controlling the antenna switch, for example, three or more receivers. A coarse arrival direction evaluation function is calculated using a digital phase detector. When this rough arrival direction evaluation function is used, a false image may be generated in addition to the true arrival direction, but it is calculated from all combinations of three or more different antenna elements selected by the antenna switcher. By integrating each coarse arrival direction evaluation function and obtaining a fine arrival direction evaluation function, the pseudo image is suppressed, and the arrival directions of a plurality of radio waves coming from the radio wave source can be detected with the same accuracy as the conventional example. For this reason, the number of receivers can be reduced as compared with the prior art, and a light weight and low cost direction finding device can be obtained.
[0024]
1, 1-1, 1-2,..., 1-M (M ≧ 4) are antenna elements that receive a signal from a radio wave source, 2 are antenna elements 1-1, 1-2,. The antenna tower to which 1-M (hereinafter referred to as antenna element 1) is attached, 12 is a signal processing device B that performs signal processing of signals supplied from the antenna tower 1, and 13 is a connection to each antenna element 1 in the antenna tower 2. The antenna switchers A and 14 for selecting the switching controllers A4, 4-1, 4-2,..., 4-J for sequentially switching the connections in the antenna switcher A13 and controlling all combinations are the antenna switcher A13. , 5-J distributes the output signal of the receiver 4 in two, and the phase difference between the receivers 5-1, 5-2,. Digital phase detector for phase detection so as to shift by two, 6 Received signal covariance matrix calculator, 7 is a noise subspace matrix calculator, 8 is a steering vector calculator, 15 is an arrival direction evaluation function calculator, 16 is a fine arrival direction evaluation function calculator, and 10 is an arrival direction detector. The signal processing device B12 performs signal processing using 4 to 8, 10 and 13 to 16.
[0025]
Next, the operation of FIG. 1 will be described. The symbols used here are defined below. K is the number of incoming waves, kk is the number of incoming waves (kk ≦ K), M is the number of antenna elements 1, m is the number of antenna elements (1 ≦ m ≦ M), J is the number of selected antenna elements, and jj is selected. Among the J antenna elements, the jjth (1 ≦ jj ≦ J) antenna element number, and K <J <M and 3 ≦ J. The M antenna elements are arbitrarily arranged. In addition, for example, a radio wave source that is generated from a moving body such as an automobile or a ship that is traveling several tens of kilometers away is considered. The radio wave source moves within the processing time during the angle measurement process. Approximate that there is no.
[0026]
In FIG. 1, among the received signals received by the antenna element 1-m from the K radio wave sources, the output from the antenna element 1-jj among the J antenna elements 1 selected by the antenna switch A13 is as follows. The digital signal is input to the receiver 4-jj, and the output signal from the receiver 4-jj is A / D-converted and phase-detected by the digital phase detector 5-jj. J digital signals that are output from the digital phase detectors 5-1, 5-2,..., 5-J detected corresponding to each of the antenna elements 1 are input, and a received signal covariance matrix is obtained. Desired. Here, when the total number of combinations of J selected antenna elements controlled by the switching controller A14 and switched by the antenna switching unit A13 is Y, it is expressed as "Equation 6". The digital signal detectors 5-1, 5-2,..., 5-J output Y-th J digital signals of the antenna combination into “Equation 7”, and the received signal covariance The matrix is shown in “Equation 8”.
[0027]
[Formula 6]
Figure 0003817875
[0028]
[Expression 7]
Figure 0003817875
[0029]
[Equation 8]
Figure 0003817875
[0030]
Next, in the noise subspace matrix calculator 7, eigenvalues / eigenvectors of the received signal covariance matrix, which is the output of the received signal covariance matrix calculator 6, are obtained as noise eigenvalues and noise eigenvectors, and the noise subspace matrix is obtained from the noise eigenvectors. This is calculated and shown in “Equation 9”.
[0031]
[Equation 9]
Figure 0003817875
[0032]
Further, the steering vector calculator 8 calculates a steering vector of J columns corresponding to the array response obtained from the array of the selected J antenna elements, and this is shown in “Equation 10”.
[0033]
[Expression 10]
Figure 0003817875
[0034]
The coarse arrival direction evaluation function calculator 15 calculates a coarse arrival direction evaluation function from the noise subspace matrix that is the output of the noise subspace matrix calculator 7 and the steering vector that is the output of the steering vector calculator 8. The y-th coarse arrival direction evaluation function P y (θ) (y = 1... Y) is expressed by “Equation 11”. However, the antenna switching method (selection method) in the antenna switch A13 is arbitrary.
[0035]
[Expression 11]
Figure 0003817875
[0036]
Next, in the fine arrival direction evaluation function calculator 16, the coarse arrival direction evaluation function P y (θ) for all combinations of the J antenna elements selected by the antenna switch A13 controlled by the switching controller A14 is obtained. Integration processing is performed to calculate the precision arrival direction evaluation function AP (θ). Here, the precision arrival direction evaluation function AP (θ) is expressed by “Equation 12”. As an example, FIG. 2 shows a rough arrival direction evaluation function P y (θ) and a fine arrival direction evaluation function AP (θ) obtained by integrating the function. By obtaining this, it is possible to suppress a false image caused by reducing the number of receivers and improve S / N.
[0037]
[Expression 12]
Figure 0003817875
[0038]
Finally, in the ten arrival direction detectors, the arrival directions (ψ 1 , ψ 2 ,...) Of K arrival waves are derived from the fine arrival direction evaluation function AP (θ) that is the output of the fine arrival direction evaluation function calculator 16. .., ψ K ) is detected. In this detection, similarly to the conventional FIG. 11, the angle corresponding to the peak is detected by selecting K peaks from the maximum peak of AP (θ), and the arrival direction of the incoming wave (ψ 1 , ψ 2 ,... ψ K ).
[0039]
Embodiment 2. FIG.
FIG. 3 is a block diagram showing Embodiment 2 of the present invention. In the first embodiment, the antenna switcher 13 selects J antenna elements among the M antenna elements. However, in this embodiment, one antenna element is replaced with one receiver as shown in FIG. Switching control for controlling all combinations of the antenna switch B17 and the antenna switch B17 that select any (J-1) antenna elements among the remaining (M-1) antenna elements. By using the signal processing device C19 using the device B18, it is possible to obtain an effect similar to that of the first embodiment, and to remove the antenna switching noise at the time of demodulation when listening to the signal of the radio wave source. Without adding, it is possible to always listen to the signal of the radio wave source from the receiver fixed to the antenna element.
[0040]
Embodiment 3 FIG.
4 is a block diagram showing Embodiment 3 of the present invention. In the first embodiment, the signal processing device B12 obtains the fine arrival direction evaluation function AP (θ) from the Y coarse arrival direction evaluation functions P y (θ) (y = 1... Y) by integration processing. The direction of fine arrival (ψ 1 , ψ 2 ,..., Ψ K ) was detected. In this embodiment, 20 signal processing devices D are provided as shown in FIG. 4, and instead of the fine arrival direction evaluation function calculator 16, the coarse arrival direction evaluation function P that is the output of the coarse arrival direction evaluation function calculator 15. y (θ) (y = 1... Y) is input, and is output from the coarse arrival direction detector 21 instead of the 21 coarse arrival direction detectors that detect the coarse arrival direction and the arrival direction detector 10. Assuming that there are Y coarse arrival directions (θ (y) 1 , θ (y) 2 ,..., Θ (y) K ) (y = 1... Y) for all combinations in the antenna switch 13. , Divide the ones with similar values into the same group, select K groups from a large number of groups, and perform average processing for each group (ψ 1 , ψ 2 ,..., Ψ K ). The same effect as in the first embodiment can be obtained by using the 22 arrival direction averaging processors that calculate
[0041]
Embodiment 4 FIG.
FIG. 5 is a block diagram showing Embodiment 4 of the present invention. In the third embodiment, J antenna elements are selected from among the M antenna elements in the antenna switch 13, but in this embodiment, one antenna element is replaced with one receiver as shown in FIG. And an antenna switch B17 that selects an arbitrary (J-1) element antenna element among the remaining (M-1) antenna elements, and a switch that controls all combinations of the antenna switch B17. By using the signal processing device D23 using the controller B18, it is possible to obtain the same effect as in the third embodiment, and to remove the antenna switching noise at the time of demodulation when listening to the signal of the radio wave source It is possible to always listen to the signal of the radio wave source from the receiver fixed to the antenna element without adding.
[0042]
Embodiment 5 FIG.
6 is a block diagram showing Embodiment 5 of the present invention. In the first embodiment, the switching controller A14 controls all Y combinations of selecting J antenna elements from M antenna elements. However, depending on the combination of antennas, the overall antenna aperture diameter may be reduced. Angular accuracy deteriorates. In this embodiment, as shown in FIG. 6, the number of combinations for selecting J antenna elements from M antenna elements is not the total number of combinations, but each antenna element is taken once or more, and the total antenna aperture diameter By using the 25 signal processing devices F that perform angle measurement signal processing using the 24 switching controllers C that select and control the combinations for which a large value can be obtained, the processing of the switching controller A14 becomes easier and the processing time The effect similar to that of the first embodiment can be obtained while being faster.
[0043]
Embodiment 6 FIG.
FIG. 7 is a block diagram showing Embodiment 6 of the present invention. In the second embodiment, the switching controller B18 controls all combinations of selecting (J-1) antenna elements from (M-1) antenna elements. However, depending on the combination of antennas, the entire antenna may be controlled. The aperture diameter becomes smaller and the angle measurement accuracy becomes worse. In the sixth embodiment, even if the number of combinations for selecting (J-1) antenna elements from (M-1) antenna elements is not the total number of combinations as shown in FIG. The effect similar to that of the second embodiment can be obtained by using the signal processing device G27 that performs the angle measurement signal processing using the switching controller D26 that takes one or more times and selects and controls the combination in which the entire antenna aperture diameter is large. Can be obtained. As an example, FIG. 8 shows a combination of antenna selections in the case of seven antenna elements and three receivers.
[0044]
Embodiment 7 FIG.
FIG. 9 is a block diagram showing Embodiment 7 of the present invention. In Embodiment 3, the switching controller A14 controls all Y combinations of selecting J antenna elements from M antenna elements. However, depending on the combination of antennas, the overall antenna aperture diameter may be reduced. Angular accuracy deteriorates. In this embodiment, even if the number of combinations for selecting J antenna elements from M antenna elements is not the total number of combinations as shown in FIG. The effect similar to that of the third embodiment can be obtained by using the signal processing device H28 that performs angle measurement signal processing using the switching controller C24 that selects and controls the combination having a large opening diameter.
[0045]
Embodiment 8 FIG.
10 is a block diagram showing Embodiment 8 of the present invention. In the fourth embodiment, all combinations of selecting (J-1) antenna elements from (M-1) antenna elements in the switching controller B18 are controlled. However, depending on the combination of antennas, the entire antenna may be controlled. The aperture diameter becomes smaller and the angle measurement accuracy becomes worse. In this embodiment, even if the number of combinations for selecting (J-1) antenna elements from (M-1) antenna elements is not the total number of combinations as shown in FIG. Even if the signal processing device I29 that performs angle measurement signal processing using the switching controller D26 that selects and controls the combination that can take a large number of times and that the entire antenna aperture diameter is large can be obtained. Can do.
[0046]
【The invention's effect】
According to the first to fifth aspects of the present invention, the number of receivers that process a received signal is small while maintaining detection accuracy comparable to that of a direction finding device that uses the same number of receivers as the number of antenna elements. A lightweight and low-cost direction finding device can be obtained.
[0047]
According to the third and fifth inventions, one of the M antenna elements is always connected to one of the receivers, and the other (M−1) antenna elements (J− 1) By providing an antenna switch for selecting one antenna element, the same effect as in the first invention can be obtained, and one antenna element is always connected to one receiver when switching the antenna element. Therefore, when listening to a signal from a radio wave source, it is possible to directly demodulate and listen to the sound from the receiver without newly adding a device for removing antenna switching noise during demodulation.
[0048]
Furthermore, according to the fourth and fifth inventions, a combination that does not use the total number of antenna elements in the antenna switching unit, takes each antenna element once or more, and can take a large antenna aperture diameter. By performing the angle measurement signal processing in minutes, the calculation processing time can be reduced as compared with the first to third inventions.
[Brief description of the drawings]
FIG. 1 is a diagram showing an example of a first embodiment of a direction finding device according to the present invention.
FIG. 2 is a diagram showing that a false image is suppressed by integration processing;
FIG. 3 is a diagram showing an example of a second embodiment of the direction finding device according to the present invention;
FIG. 4 is a diagram showing an example of a third embodiment of a direction finding apparatus according to the present invention.
FIG. 5 is a diagram showing an example of a fourth embodiment of the direction finding apparatus according to the present invention.
FIG. 6 is a diagram showing an example of a fifth embodiment of the direction finding apparatus according to the present invention.
FIG. 7 is a diagram showing an example of a sixth embodiment of the direction finding apparatus according to the present invention.
FIG. 8 is a diagram illustrating an example of a combination of antennas.
FIG. 9 is a diagram showing an example of a seventh embodiment of the direction finding apparatus according to the present invention.
FIG. 10 is a diagram showing an example of an embodiment 8 of the direction detecting apparatus according to the present invention.
FIG. 11 is a diagram illustrating an example of a conventional direction finding device.
FIG. 12 is an operation diagram.
FIG. 13 is a diagram showing a detection method in an arrival direction detector.
[Explanation of symbols]
1 antenna element, 2 antenna tower, 3 signal processor A, 4 receiver, 5 digital phase detector, 6 received signal covariance matrix calculator, 7 noise subspace matrix calculator, 8 steering vector calculator, 9 direction of arrival Evaluation customary calculator, 10 arrival direction detector, 11 automobile, 12 signal processor B, 13 antenna switcher A, 14 switching controller A, 15 coarse arrival direction evaluation function calculator, 16 fine arrival direction evaluation function calculator, 17 antenna switch B, 18 switch controller B, 19 signal processor C, 20 signal processor D, 21 coarse arrival direction detector, 22 arrival direction average processor, 23 signal processor E, 24 switch controller C, 25 Signal processor F, 26 Switching controller D, 27 Signal processor G, 28 Signal processor H, 29 Signal processor I

Claims (5)

電波源から到来する少なくとも2つ以上の電波を受信する少なくとも4個以上の複数のアンテナ素子と、前記複数のアンテナ素子のうち任意の3個以上のアンテナ素子を選択するアンテナ切替器と、前記アンテナ切替器で選択されるアンテナ素子の選択の組み合わせを制御する切替制御器と、前記アンテナ切替器で選択されるそれぞれのアンテナ素子端に接続される3個以上の受信機と、前記それぞれの受信機によって供給される信号から受信信号共分散行列を算出する受信信号共分散行列算出器と、前記受信信号共分散行列算出器より供給される受信信号共分散行列から固有値/固有ベクトルを算出し、前記固有ベクトルから構成される雑音部分空間行列を算出する雑音部分空間行列算出器と、前記アンテナ切替器で選択されたアンテナ素子間の配列に基づいてステアリングベクトルを算出するステアリングベクトル算出器と、前記雑音部分空間行列算出器から供給される雑音部分空間行列と前記ステアリングベクトル算出器から供給されるステアリングベクトルに基づいて、粗到来方向評価関数を算出する粗到来方向評価関数算出器と、前記アンテナ切替器で選択されるアンテナ素子の全組み合わせに対応して前記粗到来方向評価関数算出器から供給されるそれぞれの粗到来方向評価関数を積分し、精到来方向評価関数を算出する精到来方向評価関数算出器と、前記精到来方向評価関数算出器から供給される精到来方向評価関数から電波の到来方向を検出する到来方向検出器とを備えたことを特徴とする方向探知装置。  A plurality of antenna elements that receive at least two or more radio waves coming from a radio wave source, an antenna switch that selects any three or more antenna elements among the plurality of antenna elements, and the antenna A switching controller for controlling a combination of antenna element selections selected by the switch; three or more receivers connected to respective antenna element ends selected by the antenna switch; and the respective receivers. A received signal covariance matrix calculator for calculating a received signal covariance matrix from the signal supplied by the receiver, and calculating an eigenvalue / eigenvector from the received signal covariance matrix supplied from the received signal covariance matrix calculator. A noise subspace matrix calculator for calculating a noise subspace matrix comprising the antenna element selected by the antenna switch A steering vector calculator that calculates a steering vector based on an array between the noise subspace matrix, a noise subspace matrix supplied from the noise subspace matrix calculator, and a steering vector supplied from the steering vector calculator Coarse arrival direction evaluation function calculator for calculating a direction evaluation function and each coarse arrival direction evaluation function supplied from the coarse arrival direction evaluation function calculator corresponding to all combinations of antenna elements selected by the antenna switcher The arrival direction detection that detects the arrival direction of the radio wave from the fine arrival direction evaluation function calculator that integrates the function and calculates the fine arrival direction evaluation function and the fine arrival direction evaluation function supplied from the fine arrival direction evaluation function calculator A direction finding device comprising a device. 電波源から到来する少なくとも2つ以上の電波を受信する少なくとも4個以上の複数のアンテナ素子と、前記複数のアンテナ素子のうち任意の3個以上のアンテナ素子を選択するアンテナ切替器と、前記アンテナ切替器で選択されるアンテナ素子の選択の組み合わせを制御する切替制御器と、前記アンテナ切替器で選択されるそれぞれのアンテナ素子端に接続される3個以上の受信機と、前記それぞれの受信機によって供給される信号から受信信号共分散行列を算出する受信信号共分散行列算出器と、前記受信信号共分散行列算出器より供給される受信信号共分散行列から固有値/固有ベクトルを算出し、前記固有ベクトルから構成される雑音部分空間行列を算出する雑音部分空間行列算出器と、前記アンテナ切替器で選択されたアンテナ素子間の配列に基づいてステアリングベクトルを算出するステアリングベクトル算出器と、前記雑音部分空間行列算出器から供給される雑音部分空間行列と前記ステアリングベクトル算出器から供給されるステアリングベクトルに基づいて粗到来方向評価関数を算出する粗到来方向評価関数算出器と、前記アンテナ切替器で選択されるアンテナ素子の全組み合わせに対応して前記粗到来方向評価関数算出器から供給されるそれぞれの粗到来方向評価関数から電波の到来方向を検出する粗到来方向検出器と、前記粗到来方向検出器から供給される全組み合わせ分の粗到来方向を平均し、精到来方向を算出する到来方向平均処理器とを備えたことを特徴とする方向探知装置。  A plurality of antenna elements that receive at least two or more radio waves coming from a radio wave source, an antenna switch that selects any three or more antenna elements among the plurality of antenna elements, and the antenna A switching controller for controlling a combination of antenna element selections selected by the switch; three or more receivers connected to respective antenna element ends selected by the antenna switch; and the respective receivers. A received signal covariance matrix calculator for calculating a received signal covariance matrix from the signal supplied by the receiver, and calculating an eigenvalue / eigenvector from the received signal covariance matrix supplied from the received signal covariance matrix calculator. A noise subspace matrix calculator for calculating a noise subspace matrix comprising the antenna element selected by the antenna switch Steering vector calculator for calculating a steering vector based on the arrangement between, a noise subspace matrix supplied from the noise subspace matrix calculator, and a coarse arrival direction based on the steering vector supplied from the steering vector calculator Coarse arrival direction evaluation function calculator for calculating an evaluation function and each coarse arrival direction evaluation function supplied from the coarse arrival direction evaluation function calculator corresponding to all combinations of antenna elements selected by the antenna switch. A coarse arrival direction detector that detects the arrival direction of radio waves from the signal, and an arrival direction average processor that calculates the fine arrival direction by averaging the coarse arrival directions for all combinations supplied from the coarse arrival direction detector. A direction finding device characterized by that. 前記アンテナ切換器は、前記複数のアンテナ素子のうちの1つのアンテナ素子を直接前記受信機の一つに接続したまま、他の複数のアンテナ素子のうちの任意の2個以上のアンテナ素子を選択する手段を備えたことを特徴とする請求項1もしくは2記載の方向探知装置。  The antenna switching unit selects any two or more antenna elements of the plurality of other antenna elements while directly connecting one antenna element of the plurality of antenna elements to one of the receivers. The direction detecting device according to claim 1 or 2, further comprising means for performing the operation. 前記切替制御器は、それぞれのアンテナ素子を1回以上選択するよう制御する手段を備えたことを特徴とする請求項1もしくは2記載の方向探知装置。  3. The direction detecting apparatus according to claim 1, wherein the switching controller includes means for controlling each antenna element to be selected one or more times. 前記切替制御器は、前記他の複数のアンテナ素子を1回以上選択するよう制御するよう制御する手段を備えたことを特徴とする請求項3記載の方向探知装置。  4. The direction detecting apparatus according to claim 3, wherein the switching controller includes means for controlling to select the other plurality of antenna elements one or more times.
JP33693997A 1997-12-08 1997-12-08 Direction finding device Expired - Fee Related JP3817875B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP33693997A JP3817875B2 (en) 1997-12-08 1997-12-08 Direction finding device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP33693997A JP3817875B2 (en) 1997-12-08 1997-12-08 Direction finding device

Publications (2)

Publication Number Publication Date
JPH11166966A JPH11166966A (en) 1999-06-22
JP3817875B2 true JP3817875B2 (en) 2006-09-06

Family

ID=18304034

Family Applications (1)

Application Number Title Priority Date Filing Date
JP33693997A Expired - Fee Related JP3817875B2 (en) 1997-12-08 1997-12-08 Direction finding device

Country Status (1)

Country Link
JP (1) JP3817875B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002214318A (en) * 2001-01-19 2002-07-31 Futaba Corp Incoming wave estimation system
JP5047661B2 (en) * 2007-03-27 2012-10-10 三菱電機株式会社 Angle measuring device and angle measuring method
JP5955140B2 (en) * 2012-07-10 2016-07-20 三菱電機株式会社 Angle measuring device and angle measuring method
CN112034416B (en) * 2020-08-07 2023-07-07 浙江大学 Method for automatically detecting direction-finding unmanned aerial vehicle

Also Published As

Publication number Publication date
JPH11166966A (en) 1999-06-22

Similar Documents

Publication Publication Date Title
CN110018452B (en) Method and apparatus for estimating direction of arrival using generation of virtual received signals
WO2006067857A1 (en) Arrival direction estimating device and program
KR102422396B1 (en) Method of spatial interpolation for linear phased array antenna and appratus thereof
JP6556399B2 (en) Radar equipment
JP2001281316A (en) Radio wave arrival direction estimating device and directional variable transmitter-receiver
KR102183439B1 (en) Method and apparatus for estimating direction of arrival using combined beamspace music and tma
JP3817875B2 (en) Direction finding device
US11614535B2 (en) Radar apparatus for vehicle and method for controlling the same
JP3629195B2 (en) Arrival direction estimation apparatus and arrival direction estimation method
KR20140006275A (en) Hybrid direction finding antenna and its using receiver and method
JP2005003579A (en) Angle measurement system and positioning system
JP5047002B2 (en) Wave number estimation device
JPH08201498A (en) Angle measuring apparatus
WO2020261834A1 (en) Receiver, and radar apparatus, vehicle, and communication system provided with receiver
JP2948459B2 (en) Angle measuring device
JP3647621B2 (en) Radio direction detector
JP2002048853A (en) Radio wave arrival direction estimator and directivity variable transmitter-receiver
JP5025170B2 (en) Arrival wave number detector
JP5722026B2 (en) Direction of arrival estimation method
JP5377345B2 (en) Radio wave receiver and direction of arrival measurement method
JP4119719B2 (en) Mobile station direction estimation method and apparatus
JP4660562B2 (en) Mobile station direction estimation method and apparatus
JP7519817B2 (en) Direction of arrival estimation device, direction of arrival estimation method, and direction of arrival estimation program
JP2619144B2 (en) Direction finder
JP4991279B2 (en) Autonomous mobile device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040225

RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7421

Effective date: 20040720

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20060220

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060307

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060425

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060605

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

Free format text: PAYMENT UNTIL: 20100623

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20100623

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20110623

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20120623

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20130623

Year of fee payment: 7

LAPS Cancellation because of no payment of annual fees