(57)【要約】【課題】アレイアンテナでの受信信号数を正確に検出する。【解決手段】方位検出回路４における、アダプティブアレーによる到来方位検出の位ベクトルに基づく式（４）による計算を評価式計算回路５で行い、その計算された方位方向を予備知識としたＤＣＭＰのアルゴリズムにより最適ウェイトを算出し、乗算回路７で干渉波抑圧信号を出力するとき、式（４）においては、到来波数を正確に設定してはじめて、乗算回路７における干渉波抑圧が適正に行なわれる。  そこで、式（４）における信号数を変化させつつ算出される特徴量算出回路８からの干渉波抑圧信号の特徴量を評価回路９において相対的に評価し、真の受信信号数（到来波数）を適正かつ正確に判定し検出する。【選択図】図１
2008008781.tif 000002
【特許請求の範囲】
【請求項１】
アレイ状に配列され、高周波信号を受信する複数のアンテナと、
これら複数のアンテナにそれぞれ接続され、受信された前記高周波信号をＡ／Ｄ変換する複数のＡ／Ｄ変換器と、
これら複数のＡ／Ｄ変換器に共通接続され、アダプティブアレーにより前記受信された高周波信号の到来方位を検出する方位検出手段と、
この方位検出手段により到来方位が検出された受信高周波信号の方向ベクトルに基づく、コレクトステアリングベクトルを複数回計算する計算手段と、
この計算手段により計算された前記コレクトステアリングベクトルの方位方向を予備知識とした方向拘束付出力電力最小化法（ＤＣＭＰ）のアルゴリズムにより、前記複数のアンテナのアンテナパターンに対する最適ウェイトを算出するウェイト算出手段と、
このウェイト算出手段で算出した最適ウェイトと前記複数のＡ／Ｄ変換器から出力された受信高周波信信号とを乗算して干渉波抑圧信号を出力する乗算手段と、
この乗算手段から出力された前記干渉波抑圧信号の特徴量を算出する特徴量算出手段と、
この特徴量算出手段で算出された前記特徴量を評価する評価手段と
を具備することを特徴とする到来波数検出装置。
【請求項２】
前記方位検出手段は、アダプティブヌルステアリングにより、前記受信された高周波信号の到来方位を検出することを特徴とする請求項１に記載の到来波数検出装置。
【請求項３】
前記計算手段は、前記コレクトステアリングベクトルは下記式に基づいて算出することを特徴とする請求項１または請求項２に記載の到来波数検出装置。
ａ（θｄ）＝（Ｉ-ＶＮＶＮＨ）ｂ（θｄ）
但し、
Ｉ      ：  単位行列
ＶN     ：  雑音空間部分の固有ベクトル
Ｈ      ：  ＶN の複素共役転置
ａ(θd)：  補正前のステアリングベクトル
【請求項４】
前記計算手段は、前記方位検出手段で検出される受信高周波信号の受信数を変化させて前記コレクトステアリングベクトルを計算することを特徴とする請求項１ないし請求項３のうちのいずれか１項に記載の到来波数検出装置。
【請求項５】
前記特徴量算出手段は、前記干渉波抑圧信号の電力、または周波数スペクトル、あるいは信号対雑音比（Ｓ／Ｎ）を特徴量として算出することを特徴とする請求項１ないし請求項４のうちのいずれか１項に記載の到来波数検出装置。
【請求項６】
前記特徴量算出手段は、高速フーリエ変換（ＦＦＴ）回路で構成され、前記干渉波抑圧信号の波形全体の特徴量を算出することを特徴とする請求項１ないし請求項５のうちのいずれか１項に記載の到来波数検出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、アレイアンテナを介した高周波信号の受信において、受信方位の検出に採用して好適な高周波信号の到来波数検出装置の改良に関する。
【背景技術】
【０００２】
アレイアンテナを備えた３次元レーダ等では、移相制御によるビーム走査を行い、送信パルスを一定周期で送信し、移動目標の運動に対応したドップラ周波数の目標反射信号を受信して目標位置の検出が行われる。
【０００３】
その３次元レーダでは、目標反射信号から、目標の位置する方位を検出するために、図２に示す方位検出装置が組み込まれ、受信信号の到来方位の角度測定、すなわち測角が行われる。
【０００４】
図２に示したように、従来の方位検出装置は、アレイ状に配列されたｎ（複数）個のアンテナ１（１１〜１ｎ）で目標反射信号を受信し、その受信高周波信号は、各アンテナ１（１１〜１ｎ）に対応して接続された受信回路のＡ／Ｄ変換器２（２１〜２ｎ）にそれぞれ供給される。
【０００５】
Ａ／Ｄ変換され、ディジタル化された受信高周波信号は、相関行列演算器３１に供給される。
【０００６】
相関行列演算器３１は、各アンテナ１（１１〜１ｎ）における受信高周波信号に関する相関行列の演算から固有ベクトルと固有値とを算出し、算出された固有ベクトルは測角演算器３２に供給され、また同じく算出された固有値は到来波数検出器３３に供給される（例えば、非特許文献３参照。）。
【０００７】
測角演算器３２では、下記（１）式で示されるＭＵＳＩＣ（ MUltiple SIgnal Classification ）法のアルゴリズムによる測角演算により、受信された受信高周波信号毎にその到来方位を検出することができる。（例えば、同じく非特許文献３参照。）。
【数１】
2008008781.tif 000003
【０００８】
但し、
ＰＭＵ        ＭＵＳＩＣスペクトルの評価関数
α（θ）    方向ベクトル
ＥＮ        固有ベクトルと到来信号数から決定される雑音部分空間
Ｅ＊Ｎ      ＥＮ  複素共役転置
α＊(θ)    α（θ）の複素共役転置
すなわち、測角演算器３２は、相関行列演算器３１から供給された固有ベクトルと各アンテナ１（１１〜１ｎ）の配置位置により定まる方向ベクトルとで定義されるＭＵＳＩＣスペクトルの評価関数ＰＭＵ の値が、各アンテナ１（１１〜１ｎ）の配置位置により定まる方向ベクトルα（θ）と実際に受信された高周波信号の到来方向とが一致したときに、評価関数のスペクトラム（信号強度）が到来（角度）方向に鋭い極値（ピーク）を示すことを利用して到来方向の検出、すなわち測角が行われる。
【０００９】
測角演算器３２における、ＭＵＳＩＣ法のアルゴリズムによる受信高周波信号（到来波）の方位検出では、上記（１）式から分かるように、固有ベクトルと受信信号数（到来波数）から決定される雑音部分空間ＥＮ を信号部分空間から分離する必要がある。
【００１０】
そこで、雑音部分空間ＥＮ を信号部分空間から分離するために、到来波数検出器３３において検出された受信高周波信号数（すなわち、到来波数）が測角演算器３２に供給される。
【００１１】
到来波数検出器３３では、相関行列演算器３１で算出された固有値と、下記（２）及び（３）式で示される因子分析という手法が用いられたＭＤＬ（Minimum Description Length）規範やＡＩＣ（Akaike Information Criterion）規範の評価関数を用いて、ほぼ同時に受信される到来波数の推定または検出が行われる。
【数２】
2008008781.tif 000004
【数３】
2008008781.tif 000005
【００１２】
但し、
ＭＤＬ（ｋ）  ＭＤＬ規範の評価関数
ＡＩＣ（ｋ）  ＡＩＣ規範の評価関数
Ｎ      相関行列計算に使用したサンプル数
ｌｉ    相関行列の固有値を大きさの大きい方から順に並べたときのｉ番目の固有値
ｐ      相関行列の固有値の個数
【非特許文献１】福江, 畝田, 外園, 「レーダ信号処理におけるフーリエ係数を用いたＭＵＳＩＣ法による分離・測角性能の改善」  信学論   Vol.J85-B, no.12, pp.2380-2388 Dec. 2002.
【非特許文献２】M. Wax and T. Kailath,「 Detection of signals by information theoretic criteria」 IEEE Trans. Acoustics, Speech and Signal Processing, Vol. ASSP-33, no.2, pp.387-392, Apr. 1985.
【非特許文献３】R. O. Schmidt, 「Multiple Emitter Location and Signal Parameter Estimation」 IEEE Trans. Antennas and Propagation, Vol. AP-34, no3, pp276-280, Mar. 1986.
【発明の開示】
【発明が解決しようとする課題】
【００１３】
ＭＵＳＩＣ法のアルゴリズムを採用して受信高周波信号の到来方位を検出する従来の測角演算回路では、到来波数検出器３３からの到来波数の検出を得て、信号部分空間と雑音部分空間とを分離して測角操作が行われる。
【００１４】
そこで到来波数検出器３３では、上記のように、因子分析という手法が用いられたＭＤＬ規範等の評価関数は、アレイ状のアンテナ間の相関に基づくので、ホワイトノイズ（白色雑音）のように、複数のアンテナ間において無相関な雑音環境のもとではともかく、空電雑音や干渉波等のように複数のアンテナ間で相関を有したノイズが受信される環境のもとでは、目的とする到来波数を正確に検出できないという欠点があった。
【００１５】
そこで本発明は、干渉波等の雑音がアレイアンテナ間で相関を有して受信される環境下においても、到来数をより正確に検出可能な到来波数検出装置を提供することを目的とする。
【課題を解決するための手段】
【００１６】
本発明の到来波数検出装置は、アレイ状に配列され、高周波信号を受信する複数のアンテナと、これら複数のアンテナにそれぞれ接続され、受信された前記高周波信号をＡ／Ｄ変換する複数のＡ／Ｄ変換器と、これら複数のＡ／Ｄ変換器に共通接続され、アダプティブアレーにより前記受信された高周波信号の到来方位を検出する方位検出手段と、この方位検出手段により到来方位が検出された受信高周波信号の方向ベクトルに基づき、求まるコレクトステアリングベクトルを信号数に応じて計算する計算手段と、この計算手段により計算された前記コレクトステアリングベクトルの方位方向を予備知識とした方向拘束付出力電力最小化法（ＤＣＭＰ）のアルゴリズムにより、前記複数のアンテナのアンテナパターンに対する最適ウェイトを算出するウェイト算出手段と、このウェイト算出手段で算出した最適ウェイトＷoptと前記複数のＡ／Ｄ変換器から出力された受信高周波信信号とを乗算して干渉波抑圧信号を出力する乗算手段と、この乗算手段から出力された前記干渉波抑圧信号の特徴量を算出する特徴量算出手段と、この特徴量算出手段で算出された前記特徴量を評価する評価手段とを具備することを特徴とする。
【発明の効果】
【００１７】
アレイ状に配置されたアンテナにより高周波信号を受信し、アダプティブヌルステアリングにより干渉波抑圧信号を出力するとき、実形成される受信アンテナパターンと理想とされる設計上の受信アンテナパターンとの間に誤差を生ずるのは避けられず、その両者間の誤差に起因して、干渉波抑圧に誤差が発生する。
【００１８】
そこで、現実の受信アンテナパターンと設計上の受信アンテナパターンとの間の誤差を補正し、正しい測角を行なうべく、到来方向を決めるコレクトステアリングベクトルａ(θd)を表した下記（４）式の評価式を用いる手法が知られている。（例えば、［非特許文献４］「IEEE TRANSACTIONS ON SIGNAL PROCESSING 」VOL.42.NO.6.JUNE 1994 pp1543-1547 参照。）。
【００１９】
ａ（θｄ）＝（Ｉ-ＶＮＶＮＨ）ｂ（θｄ）      ・・・（４）
但し、
Ｉ      ：  単位行列
ＶＮ    ：  雑音空間部分の固有ベクトル
Ｈ      ：  ＶＮ の複素共役転置
ｂ(θd)：  補正前のステアリングベクトル
上記（４）式において、仮の到来波数を設定し、アンテナ１の数ｎを（例えば、ｎ＝１０）とし、到来波数Ｓを３とすると、雑音空間の固有ベクトルＶＮ における行列の次元（N）は (ｎ−Ｓ)＝７となるから、雑音空間の固有ベクトルＶＮ における行列の次元（N）は、（１０−３＝）７次元の行列式となる。
【００２０】
また、方向拘束付出力電力最小化法（ＤＣＭＰ：Directionally Constrained Minimization of Power、以下ＤＣＭＰと称する）によるアレイ信号処理アルゴリズムにより、最適ウェイト（ベクトル）Ｗoptを得るアダプティブビームフォーミングでは、受信電波の到来方向を予備知識として持つ必要があることが知られている。（例えば、［非特許文献５］菊間 信良著 「アレイアンテナによる適応信号処理」科学技術出版  1999 、参照）。
【００２１】
そこで本発明は、ＤＣＭＰのアルゴリズムにより最適ウェイトＷoptを算出して、干渉波抑圧信号を出力するとき、上記（４）式の共分散行列を指針として、現実の受信アンテナパターンと理想上の受信アンテナパターンとの間の誤差が正しく補正されたときにはじめて干渉波の除去ないしは抑圧が適正に行なわれ、かつコレクトステアリングａ(θd)は、設定した到来波数、すなわち行列の次元を規定する受信信号数に応じて求められることに着目してなされたもので、その行列の次元を規定した設定受信信号数が正しいか否かを干渉波抑圧信号の特徴点から判定するものである。
【００２２】
すなわち、上記（４）式で表されたコレクトステアリングベクトルａ(θd)では、到来方位のアンテナパターン、すなわち方向ベクトルＣに当て嵌め補正され、到来波数が正しいか否かによって特徴量に差が生ずる（例えば、上記
【非特許文献４】ｐ９０．参照。）。
【００２３】
そこで、本発明の到来波数検出装置は、アダプティブアレーにより到来方位を検出する方位検出手段と、検出された受信高周波信号の方位ベクトルに基づくコレクトステアリングベクトルを計算する計算手段と、計算されたコレクトステアリングベクトルの方位方向を予備知識としたＤＣＭＰのアルゴリズムにより最適ウェイトを算出するウェイト算出手段と、算出した最適ウェイトと受信高周波信信号とを乗算して干渉波抑圧信号を出力する乗算手段と、出力された干渉波抑圧信号の特徴量を算出する特徴量算出手段とを具備し、干渉波抑圧信号の評価手段における特徴量の評価から、受信信号数、すなわち到来波数を正確に判定する。
【発明を実施するための最良の形態】
【００２４】
以下、本発明による到来波数検出装置の一実施例を、図１を参照して詳細に説明する。なお、到来波数検出装置は方位測定装置への組み込みに適したものであるから、図２に示した従来の方位測定装置と同一構成には同一符号を付して、詳細な説明は省略する。
【００２５】
すなわち、この実施例の到来波数検出装置は、目標反射信号等の高周波信号を受信するアレイ状に配列されたｎ（複数）個のアンテナ１（１１〜１ｎ）と、各アンテナ１（１１〜１ｎ）に対応して接続されたアナログ／ディジタル（Ａ／Ｄ）変換回路２（２１〜２ｎ）と、各Ａ／Ｄ変換回路２（２１〜２ｎ）においてディジタル化された受信高周波信号を導入し、アダプティブヌルステアリングにより、受信高周波信号の到来方位を検出する方位検出回路４と、この方位検出回路４により検出される受信信号の到来方位に基づき、アンテナパターンにおける到来方向のコレクトステアリングベクトルを上記（４）式により示される評価式を指針として計算する計算回路５と、この計算回路５により計算された固有ベクトルを有するアンテナパターンに関し、ＤＣＭＰのアルゴリズムにより最適ウェイトＷoptを算出するウェイト算出回路６と、このウェイトベクトル算出回路６で算出した最適ウェイトＷoptと複数のＡ／Ｄ変換器２（２１〜２ｎ）からのディジタル化された受信高周波信号とを乗算して干渉波抑圧信号を出力する乗算回路７と、この乗算回路７から出力された干渉波抑圧信号の波形全体の特徴量、すなわち例えば前記干渉波抑圧信号の電力、または周波数スペクトル、あるいは信号対雑音比（Ｓ／Ｎ）を算出すべく高速フーリエ変換（ＦＦＴ）回路で構成された特徴量算出回路８と、この特徴量算出回路８で算出された特徴量を評価する評価回路９とを備え、仮信号数設定回路１０で適宜設定された受信信号数に基づいて出力される特徴量算出回路８からの出力を評価回路９が評価するように構成されている。
【００２６】
上記構成により、方位検出回路４において形成された方位アンテナパターンにより、例えばひとつの受信高周波信号（到来波）を受信したときに、乗算回路７は最適ウェイトＷoptと受信高周波信号との乗算により、コレクトステアリングａ(θ)で表される所望方向以外の信号（干渉波）が抑圧された信号波形を出力し、ＦＦＴからなる特徴量算出回路８は、出力信号波形全体の特徴量を算出する。
【００２７】
そこで、評価式計算回路５における（４）式の計算を順次、到来波数Ｓを変えて行い、特徴量算出回路８は、その到来波数の変化に基づき乗算回路７から出力されたコレクトステアリングａ(θd)で表される干渉波抑圧信号波形全体の特徴量を算出して、評価回路９に供給する。
【００２８】
すなわち、評価式計算回路５での上記（４）式で表されたコレクトステアリングベクトルａ(θd)では、到来方位のアンテナパターン、すなわち方向ベクトルＣに当て嵌め補正されるから、到来波（受信信号）数設定回路１０で設定された到来波数が正しいか否かによって、特徴量算出回路８から出力の特徴量に差が生ずる。
【００２９】
従って、本実施例において、到来波数設定回路１０で設定された到来波数Ｓにおける特徴量算出回路８の出力特徴量を評価回路９で評価し、特徴量が到来波数Ｓが正しく設定された場合の特徴量と、正しい到来波数Ｓとは異なって設定された場合の特徴量とを相対評価することで正しい到来波数を検知できる。
【００３０】
上記説明のように、本実施例に係る到来波数検出装置によれば、到来波数が正しく設定されてはじめて、上記（４）式における補正が正しく行なわれ、それに伴い、評価回路９における特徴量の評価では、例えば特徴量を表す出力電力が最も大であるときの設定到来波数、あるいは周波数スペクトラムの信号強度のピークが最も大となったときの設定到来波数、あるいはＳ／Ｎが最も良好なときの設定到来波数を検出することで、到来波数を正確に判定し検出できる。
【００３１】
従って、本実施例に係る到来波数検出装置を方位検出に採用したとき、方位検出装置における到来波の方位を分解能良く高精度に検出できる。
【図面の簡単な説明】
【００３２】
【図１】本発明による到来波数検出装置の一実施例を示した構成図である。
【図２】従来の到来波数検出器を組み込んだ方位測定装置の構成図である。
【符号の説明】
【００３３】
１(１１〜１ｎ)  アンテナ
２  Ａ／Ｄ変換器
４  方位検出回路（ＤＢＦ）
５  評価式計算回路（計算手段）
６  ウェイト算出回路（ウェイト算出手段）
７  乗算回路（乗算手段）
８  特徴量算出回路（特徴量算出手段）
９  評価回路（評価手段）
１０  到来波数（受信信号数）設定回路（設定手段）
【図１】
2008008781.tif 000006
【図２】
2008008781.tif 000007