JP5783693B2 - レーダ装置 - Google Patents

Description

この発明は、目標の高分解能画像を生成するレーダ装置に関するものである。

レーダにおける距離分解能は占有帯域幅に依存する。例えば３センチメートルの分解能を実現するためには５ギガヘルツの帯域幅を必要とするが、周波数資源の有効利用の観点から、これだけの帯域幅を占有することは難しい。この問題に対処する方法として、一対のレーダを用い、例えば第１のレーダでは５．０から５．５ギガヘルツの周波数帯域、第２のレーダ装置では９．５から１０．０ギガヘルツの周波数帯域を使用してそれぞれ観測を行い、その受信信号スペクトルを補間することによって、あたかも５から１０ギガヘルツの周波数帯域で観測したスペクトルを推定しようとする試みがなされている。この場合、２つのレーダから目標までの距離は一般に異なるので、それぞれで観測された受信信号スペクトルの位相変化は連続でなく、あるいはコヒーレントではない。そこで、従来のレーダ装置では、第１の周波数帯域で観測する第１のレーダと第２の周波数帯域で観測する第２のレーダの空間スペクトルを、レンジ方向について位相を一致させた上で、これらの空間スペクトルに基づいて帯域を拡張することによりレンジの高分解能化を図っていた（例えば特許文献１参照）。

特許文献１に掲載された従来のレーダ装置について説明する。
図１は従来のレーダ装置の観測のジオメトリを示す説明図である。図２は同レーダ装置の機能構成を示すブロック図である。図３は図２に示す相互コヒーレント化手段の詳細な構成を示すブロック図である。また、図４は同レーダ装置の信号スペクトルを示す模式図である。
図１に示すように、略同一位置に設置された２つのレーダが同一の目標１００を観測することを前提としている。すなわち、第１のレーダの送受信アンテナ１ａと第２のレーダの送受信アンテナ１ｂにより、目標１００に対するレーダ波の送受信を行う。ただし、２つのレーダの周波数帯域は異なる。このレーダ波のスペクトルを図４に示すが、横軸はレンジ方向の空間周波数あるいはレーダの瞬時送信周波数を表している。図４（ａ）において、１０１ａは第１のレーダの受信信号、１０１ｂは第２のレーダの受信信号を表す。図４（ｂ）において、１０２ａは相互コヒーレント化手段３の出力における第１のレーダの受信信号を表し、１０２ｂは相互コヒーレント化手段３の出力における第２のレーダの受信信号を表す。なお、以下では相互コヒーレント化手段３は第１のレーダの受信信号のみを補償することを想定して説明しているので、信号１０２ｂは信号１０１ｂと等しい。また、図４（ｃ）において、１０３は帯域内挿・拡張手段４の出力信号を表している。

このレーダ装置の場合、第１および第２のレーダによるレンジ分解能は、信号１０１ａのスペクトル帯域幅と信号１０１ｂのスペクトル帯域幅でそれぞれ規定され、これらの帯域幅が広いほど高いレンジ分解能を実現できる。そこで、２つのレーダで観測されたこれらの信号のスペクトルを内挿または外挿してスペクトルの帯域幅を拡大する。そのためには信号１０１ａと信号１０１ｂの位相の連続性（コヒーレンシー）が保たれている必要があるが、図４（ａ）に示されるように、この条件は一般に成立しない。その理由は２つのレーダのアンテナ１ａ，１ｂの位相中心から目標１００までの距離が、一般に完全には一致していないからである。この問題を解決するために、相互コヒーレント化手段３を用いて信号１０１ａと信号１０１ｂの位相のスペクトルが連続になるようにスペクトルを補償するようにしている。信号１０１ａのスペクトルを補償すると、図４（ｂ）に示される信号１０２ａのスペクトルのようになる。さらに、帯域内挿・拡張手段４により、スペクトルを内挿あるいは外挿してスペクトルの帯域幅を拡大する。拡張されたスペクトルは図４（ｃ）に示される信号１０３のスペクトルのようになる。次に、レンジ圧縮手段５により、拡張された信号１０３のスペクトルをフーリエ変換し、レンジを圧縮して高いレンジ分解能を実現する。

次に、図３に示す相互コヒーレント化手段３の詳細動作を説明する。
ここでは、信号１０１ａと信号１０１ｂはＰ個の指数関数の和で表すものと仮定する。すなわち、次式のモデルで信号を表現する。

まず、ＭＵＳＩＣ手段７ａにより、信号１０１ａの極（pole）を求める。次に、all-pole推定手段８ａにおいて、極に基づいて（１）式の係数ａ_k 、ｐ_k を求め、モデルＭ₁ を決定する。信号１０１ｂについても同様の処理を行って係数ａ_k 、ｐ_k を求め、モデルＭ₂ を決定する。次に、非線形最小二乗手段９において、（２）式で表す指標Ｃを最小化するΔθとＡを求める。すなわち、２つのモデルＭ₁ とＭ₂ の位相の連続性を最大化するパラメータを探索する。ただし、Ａは２つのモデルの振幅比を調整する係数、Δθは２つのモデルの位相差を調整する係数、Ｎはスペクトルの標本点数（図４の横軸のサンプル数）とする。またｊは√（−１）である。

最後に、スペクトル補償手段１０において、非線形最小二乗手段９で求めたパラメータΔθと係数Ａを用いて信号１０１ａを補償し、図４（ｂ）に示されるように位相の連続性が保たれた信号１０２ａと１０２ｂの組が得られる。

次に、帯域内挿・拡張手段４の詳細動作を説明する。
なお、ここでも信号１０２ａと信号１０２ｂを合成した信号は指数関数の和で表されるものと仮定し、（３）式のモデルで信号を表現する。（３）式は（１）式と同じ形をしているが、信号１０２ａと信号１０２ｂを合成した信号に対するモデルであるので、こちらのモデルＭを global-pole と呼ぶことにする。

まず、ＭＵＳＩＣ手段７と同様の手段により、信号１０２ａと信号１０２ｂを合成した信号の極を求める。次に、all-pole推定手段８と同様の手段で、極に基づいて（３）式の係数ａ_k 、ｐ_k を求め、モデルＭの初期値を決定する。次に、非線形最小二乗手段９と同様の手段で、（４）式に示される指標Ｊを最小化するモデルＭを求めるためにａ_k 、ｐ_k をパラメータとして探索する。ただし、ｓ_n は信号１０２ａと信号１０２ｂを合成した信号のスペクトル、ｗ_n は適当な重み係数、＜ｎ＞は信号１０２ａと信号１０２ｂが値をもつ標本点である。

最後に、推定されたパラメータａ_k 、ｐ_k を（３）式に当てはめて帯域拡張された信号１０３のスペクトルを得る。以上のように動作することによって、従来のレーダ装置ではレンジについて高分解能化された目標の画像を得ていた。

米国特許第５９４５９４０号公報

従来のレーダ装置は、帯域の内挿・拡張を実施するために、送信帯域の異なる複数のレーダの信号を合成する前に、相互コヒーレント化手段によって観測信号の位置あわせおよび位相合わせを実施していた。また、各レーダで観測された信号の位置ずれと位相ずれの量を推定するために、信号合成前に各レーダの信号に対してＭＵＳＩＣ法のような超解像アルゴリズムを適用していた。この従来のレーダ装置において、精度良く位置あわせと位相あわせを実施し、帯域の内挿・拡張が成功するためには、各レーダで観測された信号の位置ずれと位相ずれを正確に推定して補償する必要がある。すなわち、信号合成前の各レーダの信号に対して適用した超解像アルゴリズムにより、精度の良い信号推定が達成されている必要があった。このことは、この従来のレーダ装置においては、各レーダの信号を合成する前に、各レーダで観測された信号をそれぞれ別々に利用しても十分に分解能が上げることが出来る環境でないと帯域の内挿・拡張が成功しないことを示唆している。したがって、この従来のレーダ装置による帯域の内挿・拡張の性能は、例えば、信号対雑音電力比（ＳＮＲ：Signal to Noise Ratio）が低い環境では著しく低下してしまう問題があった。

この発明は、上記問題点を解決するためになされたもので、各レーダの信号に超解像アルゴリズムを適用せずに位置ずれを推定して補正することが可能なレーダ装置を得ることを目的とする。

この発明に係るレーダ装置は、
それぞれ近傍に配設され、異なる周波数帯域で目標を観測する複数のレーダと、
前記各レーダの信号に対して、空間平均の手法によって信号空間の共分散行列を推定する処理を行う共分散行列推定手段と、
共分散行列推定手段で推定された各レーダの信号に対する共分散行列を合成する処理を行う共分散行列合成手段と、
前記共分散行列合成手段によって算出された合成共分散行列を用いて、超解像アルゴリズムにより高分解能レンジプロフィールを推定する超解像レンジ圧縮手段とを備え、
前記共分散行列合成手段は、
前記各レーダで観測された信号に対して、Ｄ(τ)（なお、Ｄ(τ)は、未知の位置ずれ量τを補償するための項であり、exp(j2πτm)（ただし、ｍ＝0,1,…,Ｍ−１)を対角成分とする、ＭｘＭの対角行列である。）と、Ｄ(τ) ^H とを用いて、複数の仮定した位置ずれ量のそれぞれを補償するための位相補償を施す共分散行列補償手段と、
前記共分散行列補償手段において位相補償を施した各レーダの信号に対応する共分散行列の平均を算出する共分散行列平均手段と、
前記共分散行列平均手段において算出した平均の共分散行列を固有値分解し、その固有値スペクトルの広がりを測る指標であるエントロピーを算出するエントロピー算出手段と、
前記複数の仮定した位置ずれ量の中から、前記エントロピー算出手段において算出されたエントロピーが最小となる位置ずれ量を、位置ずれ量の推定値として出力する位置ずれ判定手段と、
前記位置ずれ判定手段において出力された前記位置ずれ量の推定値に基づいて位相補償を施した各レーダの信号に対応する共分散行列の平均を、前記合成共分散行列として、前記超解像レンジ圧縮手段に出力する合成共分散行列選択手段とを備える。

この発明に係るレーダ装置によれば、各レーダの受信信号に対して、空間平均の手法によって信号空間の共分散行列を共分散行列推定手段で推定し、共分散行列合成手段により各レーダの信号に対する共分散行列を合成する処理を行う過程で各レーダの信号の位置ずれ量を補償しながら共分散行列の合成をし、この共分散行列合成手段によって算出された合成共分散行列を用いて、超解像アルゴリズムにより高分解能レンジプロフィールを超解像レンジ圧縮手段で推定することにより、各レーダの信号に超解像アルゴリズムを適用せず、各レーダの信号に対する共分散行列の合成共分散行列に超解像アルゴリズムを適用することで、位置ずれを推定して補正し高分解能レンジプロフィールを推定することが可能である効果を奏する。

従来のレーダ装置の観測のジオメトリを示す説明図である。 従来のレーダ装置の機能構成を示すブロック図である。 図２に示す相互コヒーレント化手段の詳細な構成を示すブロック図である。 従来のレーダ装置の信号スペクトルを示す模式図である。 この発明の実施の形態１によるレーダ装置の機能構成を示すブロック図である。 この発明の実施の形態１によるレーダ装置における共分散行列合成手段の詳細な機能構成を示すブロック図である。 実施の形態１によるレーダ装置の観測信号の空間スペクトルの概念図である。 この発明の実施の形態２によるレーダ装置の機能構成を示すブロック図である。 この発明の実施の形態３によるレーダ装置の機能構成を示すブロック図である。

実施の形態１．
図５はこの発明の実施の形態１によるレーダ装置の機能構成を示すブロック図である。
図５において、送受信機２ａ，２ｂは、送受信アンテナ１ａ，１ｂを用いて異なる周波数帯域で目標を観測する第１、第２のレーダを構成する。共分散行列推定手段１１ａと１１ｂは各レーダの信号に対して、空間平均の手法によって信号空間の共分散行列を推定する処理を行う手段である。共分散行列合成手段１２は、共分散行列推定手段１１ａと１１ｂで推定された各レーダの信号に対する共分散行列を合成する処理を行う手段である。超解像レンジ圧縮手段１３は、共分散行列合成手段１２で合成された共分散行列をＭＵＳＩＣ法などの超解像アルゴリズムにより高分解能レンジプロフィールを推定する手段である。表示手段６は、超解像レンジ圧縮手段１３で処理して得られた目標の画像を表示する手段である。

図６は図５に示した前記共分散行列合成手段１２の詳細な機能構成を示すブロック図である。前記共分散行列合成手段１２は、共分散行列補償手段１４と共分散行列平均手段１５とエントロピー算出手段１６と位置ずれ判定手段１７と合成共分散行列選択手段１８を備えており、その詳細については後述する。

次に、動作について説明する。本実施の形態１によるレーダ装置の観測のジオメトリは、図１に示す従来のレーダ装置の観測のジオメトリと同様である。すなわち、送信帯域の異なる複数のレーダによって、目標物を同時観測しているとする。図１には２つのレーダを図示しているが、２つ以上のいくつでもかまわない。ただし、目標とレーダの間の距離は、異なるレーダ間の距離に比べて十分に大きいものとする。このようなジオメトリで観測を行った場合、各レーダから目標までの距離は若干異なる。各レーダの位置は事前に計測しておくことは可能であるが、レーダの分解能以上のオーダで精密に位置あわせをすることは困難である。したがって、複数のレーダの信号を合成するにあたって、観測信号の位置ずれを補償する必要がある。なお、本明細書においては、以下では、特に断らない限り、レーダの数が２つの場合について記述するが、３つの場合に拡張することは容易である。

図７は本実施の形態１によるレーダ装置の観測信号の空間スペクトルの概念図である。本実施の形態１によるレーダ装置の方式は、この空間スペクトル信号に対して、Ｃａｐｏｎ法やＭＵＳＩＣ（MULTIPLE SIGNAL CLASSIFICATION）法などの超解像アルゴリズムを適用して、２つのレーダで使用している送信信号帯域幅の合計で決まる分解能よりも高い分解能を得るものである。

Ｃａｐｏｎ法やＭＵＳＩＣ法などの超解像アルゴリズムによってレーダの距離分解能を向上する場合、まず、観測された空間スペクトル信号に対して、空間スムージングの手段を適用して、共分散行列を推定する。本実施の形態１によるレーダ装置においてもこれは同様であるが、各レーダで観測された信号には位置ずれが発生していることを考慮するために、まず、前記共分散行列推定手段１１ａと１１ｂにおいて、第１のレーダと第２のレーダによって観測された空間スペクトル信号より、それぞれ共分散行列を推定する。

式(５)において、R₁ と R₂は、それぞれ、第１のレーダと第２のレーダによって観測された空間スペクトル信号から推定された共分散行列である。ｘ_kと ｙ₁ は、第１のレーダと第２のレーダによって観測された空間スペクトル信号の中から連続するＭ個のサンプルを取り出して構成したＭ次元ベクトルである。ここで、Ｍ個のサンプルの取り出し方は図７に示すように、スライディングウィンドウをかける要領で取り出せばよい。ＫとＬは、このようにして取り出したＭ次元ベクトルの数である。第１のレーダと第２のレーダの送信信号帯域は必ずしも一致している必要は無いので、ＫとＬは必ずしも同じ値ではない。

次に、前記共分散行列平均手段１５において、R₁ と R₂の共分散行列を合成する。前記の通り、各レーダで観測された信号には位置ずれが発生していることを考慮して、次式によって R₁ と R₂ の共分散行列を合成する。

ここで、Ｒ(τ) が合成された共分散行列である。以下では、これを合成共分散行列と呼ぶ。また、Ｄ(τ) は、未知の位置ずれ量τを補償するための項であり、exp(j2πτm) （ただし、ｍ＝0,1,…,Ｍ−１)を対角成分とする、ＭｘＭの対角行列である。
図６のブロック図において、前記共分散行列補償手段１４は、式（６）のＤ(τ)Ｒ₂Ｄ(τ)^Hの演算を実施する。前記共分散行列平均手段１５において、式（６）の残りの演算を実施している。

なお、共分散行列を推定するために使用するサンプル数を等価的に増加させて、推定精度を向上する方式として、Forward-Backward空間スムージング法が知られているが、ここでもこの手法を適用することは可能である。その場合、式（５）および式（６）は、それぞれ、式（７）および式（８）で置き換える。

共分散行列の最終的な推定値を得るためには、上記未知の位置ずれ量τを推定する必要がある。本実施の形態１によるレーダ装置においては、式(６)において，τが真の位置ずれに一致しない場合、点目標数増加と同等の効果があるため、Ｒ(τ）の固有値λ_m(τ）（ただし、ｍ＝1,2,…,Ｍ）のスペクトルが広がる現象があることに着目する。また、固有値スペクトルの広がりは，次式で定義されるエントロピーＨ(τ）で測ることができる。

式（１０）で定義されるエントロピーＨ(τ）を指標として、エントロピーＨ(τ）を最小にするτを探索することにより、固有値スペクトルの広がりを最小にするτ、すなわち、真の位置ずれ量に良く一致するτを見つけることが出来る。

以上を踏まえ、前記エントロピー算出手段１６は、位置ずれ量τを変化させながら、式（７）によってエントロピーＨ(τ）を算出する。次いで、前記位置ずれ判定手段１７は、前記エントロピー算出手段１６によって算出されたエントロピーＨ(τ）の値が最小となる位置ずれ量τを、位置ずれ量の推定値 τ’として出力する。前記合成共分散行列選択手段１８は、前記位置ずれ判定手段１７によって出力された位置ずれ量の推定値τ’を式（６）に代入して得られる合成共分散行列を出力する。

なお、前記共分散行列推定手段１１ａと１１ｂにおいて、Forward-Ｂackward空間スムージング法を用いた場合は、前記位置ずれ判定手段１７によって出力された位置ずれ量の推定値τ’を式（８）に代入して得られる合成共分散行列を出力する。

次いで、前記超解像レンジ圧縮手段１３では、前記合成共分散行列選択手段１８によって出力された合成共分散行列R（τ'）を用いて、Ｃａｐｏｎ法やＭＵＳＩＣ（MULTIPLE SIGNAL CLASSIFICATION）法などの超解像アルゴリズムを適用して、２つのレーダで使用している送信信号帯域幅の合計で決まる分解能よりも高い分解能を有するプロフィールを得る。なお、超解像アルゴリズムの処理手順は公知であり、この実施の形態１では公知の処理手順を用いる。

以上のように、この実施の形態１によれば、未知の位置ずれ量を様々に仮定しながら合成共分散行列を生成し、この合成共分散行列の固有値スペクトルの広がりをエントロピーを指標として評価し、エントロピーの値が最も小さくなる位置ずれ量を、未知の位置ずれ量の推定値として用いるように構成したので、各レーダの信号に超解像アルゴリズムを適用せずに位置ずれを推定して補正することが可能である効果を奏する。

実施の形態２．
図８はこの発明の実施の形態２によるレーダ装置の機能構成を示すブロック図である。
発明の実施の形態２によるレーダ装置は、相互相関関数算出手段１９と位置ずれ粗推定手段２０を備えている。

次に、動作について説明する。本実施の形態２によるレーダ装置においては、まず、前記相互相関関数算出手段１９において、第１のレーダと第２のレーダで観測したレンジプロフィールの相互相関関数を算出する。次いで、位置ずれ粗推定手段２０において、相互相関関数のピークを検索することにより、第１のレーダと第２のレーダの信号の位置ずれを大まかに推定することが可能である。この推定結果を用いて、前記共分散行列平均手段１５において、未知の位置ずれ量として設定する値の範囲を決定する。この範囲は、例えば、相互相関関数のピークを含み、各レーダの距離分解能程度の範囲に設定すれば良い。また、各レーダの距離分解能が異なる場合は、もっとも低い距離分解能に相当する範囲に設定すれば良い。

以上のように、この実施の形態２によれば、位置ずれ量の候補として、第１のレーダと第２のレーダで観測したレンジプロフィールの相互相関関数のピーク周辺に絞り込むことができるので、演算量を低減できる効果を奏する。

実施の形態３．
図９はこの発明の実施の形態３によるレーダ装置の機能構成を示すブロック図である。図９は、図５に示した実施の形態１によるレーダ装置の前記共分散行列合成手段１２の構成図の中で、前記エントロピー算出手段１６を統計モデル適合度指標算出手段２１に置き換えたものである。

前記統計モデル適合度指標算出手段２１においては、Ｒ(τ)の固有値λ_m(τ)（ただし、ｍ＝1,2,…,Ｍ）のスペクトルに対して、信号に含まれる点目標の数をパラメータとする統計モデルを当てはめたときの適合度を測る指標を算出する。このような指標としては、例えば、赤池情報基準（ＡＩＣ：Akaike Information Criteria）やＭＤＬ（Minimum Description Length）が良く知られている。

以下では、ＡＩＣを用いる場合を例に説明する。前記統計モデル適合度指標算出手段２１においては、位置ずれ量τを変化させながら、各ずれ量τにおける固有値スペクトルに対して、仮定する点目標の数を様々に変化させた場合のＡＩＣを算出する。次いで、各ずれ量τにおいて、ＡＩＣを最小値とする点目標の数を決定し、このときのＡＩＣを各ずれ量τにおけるＡＩＣの評価値として出力する。

次いで、位置ずれ判定手段１７では、これらのＡＩＣを最小値とする点目標の数が最も少ないずれ量τを選択する。このとき、点目標の数が最小となるずれ量τが一意に決まらない場合がある。この場合は、点目標の数が最小となるずれ量τのうち、ＡＩＣの評価値が最小となるずれ量τを位置ずれの推定値τ’として出力する。

以上のように、この実施の形態３によれば、統計モデルを当てはめたときの適合度を測る指標を用いてずれ量を推定するように構成したので、各レーダの信号に超解像アルゴリズムを適用せずに位置ずれを推定して補正することが可能である効果を奏する。また、前記超解像レンジ圧縮手段１３において、ＭＵＳＩＣ法のように点目標数の推定を必要とするアルゴリズムを採用する場合、位置ずれの推定値τ’におけるＡＩＣの最小値を与える点目標の数を、点目標数の推定値としてそのまま使用できる効果を奏する。

この発明によるレーダ装置は、目標の高分解能画像が生成できるので、目標画像、特に目標の高分解能画像を必要とする場合に用いるレーダ装置に利用される可能性がある。

１ａ，１ｂ；送受信アンテナ、２ａ，２ｂ；送受信機、１１ａ，１１ｂ；共分散行列推定手段、１２；共分散行列合成手段、１３；超解像レンジ圧縮手段、１４；共分散行列補償手段、１５；共分散行列平均手段、１６；エントロピー算出手段、１７；位置ずれ判定手段、１８；合成共分散行列選択手段、１９；相互相関関数算出手段、２０；位置ずれ粗推定手段、２１；統計モデル適合度指標算出手段。

Claims (6)

1. それぞれ近傍に配設され、異なる周波数帯域で目標を観測する複数のレーダと、
   前記各レーダの信号に対して、空間平均の手法によって信号空間の共分散行列を推定する処理を行う共分散行列推定手段と、
   共分散行列推定手段で推定された各レーダの信号に対する共分散行列を合成する処理を行う共分散行列合成手段と、
   前記共分散行列合成手段によって算出された合成共分散行列を用いて、超解像アルゴリズムにより高分解能レンジプロフィールを推定する超解像レンジ圧縮手段と
   を備え、
   前記共分散行列合成手段は、
   前記各レーダで観測された信号に対して、Ｄ(τ)（なお、Ｄ(τ)は、未知の位置ずれ量τを補償するための項であり、exp(j2πτm)（ただし、ｍ＝0,1,…,Ｍ−１)を対角成分とする、ＭｘＭの対角行列である。）と、Ｄ(τ)^Hとを用いて、複数の仮定した位置ずれ量のそれぞれを補償するための位相補償を施す共分散行列補償手段と、
   前記共分散行列補償手段において位相補償を施した各レーダの信号に対応する共分散行列の平均を算出する共分散行列平均手段と、
   前記共分散行列平均手段において算出した平均の共分散行列を固有値分解し、その固有値スペクトルの広がりを測る指標であるエントロピーを算出するエントロピー算出手段と、
   前記複数の仮定した位置ずれ量の中から、前記エントロピー算出手段において算出されたエントロピーが最小となる位置ずれ量を、位置ずれ量の推定値として出力する位置ずれ判定手段と、
   前記位置ずれ判定手段において出力された前記位置ずれ量の推定値に基づいて位相補償を施した各レーダの信号に対応する共分散行列の平均を、前記合成共分散行列として、前記超解像レンジ圧縮手段に出力する合成共分散行列選択手段と
   を備えることを特徴とするレーダ装置。
2. それぞれ近傍に配設され、異なる周波数帯域で目標を観測する複数のレーダと、
   前記各レーダの信号に対して、空間平均の手法によって信号空間の共分散行列を推定する処理を行う共分散行列推定手段と、
   共分散行列推定手段で推定された各レーダの信号に対する共分散行列を合成する処理を行う共分散行列合成手段と、
   前記共分散行列合成手段によって算出された合成共分散行列を用いて、超解像アルゴリズムにより高分解能レンジプロフィールを推定する超解像レンジ圧縮手段と
   を備え、
   前記共分散行列合成手段は、
   前記各レーダで観測された信号に対して、Ｄ(τ)（なお、Ｄ(τ)は、未知の位置ずれ量τを補償するための項であり、exp(j2πτm)（ただし、ｍ＝0,1,…,Ｍ−１)を対角成分とする、ＭｘＭの対角行列である。）と、Ｄ(τ)^Hとを用いて、複数の仮定した位置ずれ量のそれぞれを補償するための位相補償を施す共分散行列補償手段と、
   前記共分散行列補償手段において位相補償を施した各レーダの信号に対応する共分散行列の平均を算出する共分散行列平均手段と、
   前記共分散行列平均手段において算出した平均の共分散行列を固有値分解し、その固有値スペクトルに対して、信号に含まれる点目標数をパラメータとする統計モデルを当てはめたときの適合度を測る指標である統計モデル適合度指標を、前記複数の仮定した位置ずれ量のそれぞれにおいて仮定する点目標数ごとに算出し、前記複数の仮定した位置ずれ量のそれぞれについて前記統計モデル適合度指標が最小となる点目標数を決定し、該決定された点目標数に対応する統計モデル適合度指標の最小値を前記仮定した位置ずれ量の統計モデル適合度指標として出力する統計モデル適合度指標算出手段と、
   前記統計モデル適合度指標算出手段において前記複数の仮定した位置ずれ量のそれぞれについて決定された点目標数の中で最も少ない点目標数に対応する位置ずれ量を、該最も少ない点目標数が一意に決まらない場合には該点目標数が最も少ない位置ずれ量のうち前記統計モデル適合度指標算出手段から出力された統計モデル適合度指標が最小となる位置ずれ量を、位置ずれ量の推定値として出力する位置ずれ判定手段と、
   前記位置ずれ判定手段において出力された前記位置ずれ量の推定値に基づいて位相補償を施した各レーダの信号に対応する共分散行列の平均を、前記合成共分散行列として、前記超解像レンジ圧縮手段に出力する合成共分散行列選択手段とを備え、
   前記統計モデル適合度指標算出手段は、統計モデルの適合度を測る指標として、赤池情報基準（ＡＩＣ：Akaike Information Criteria）またはＭＤＬ（Minimum Description Length）を採用する
   ことを特徴とするレーダ装置。
3. 前記各レーダで観測したレンジプロフィールの相互相関関数を算出する相互相関関数算出手段と、
   相互相関関数のピークを検索することにより、前記各レーダの信号の間の位置ずれを大まかに推定し、前記共分散行列合成手段において様々に仮定する位置ずれ量の設定範囲を絞る位置ずれ粗推定手段と
   を備えることを特徴とする請求項１または２に記載のレーダ装置。
4. 前記位置ずれ粗推定手段は、前記共分散行列合成手段において様々に仮定する位置ずれ量の設定範囲として、相互相関関数のピークに対して、前記各レーダの距離分解能程度あるいは各レーダの距離分解能のうち最も低い距離分解能の範囲とすることを特徴とする請求項３に記載のレーダ装置。
5. 前記超解像レンジ圧縮手段において、超解像アルゴリズムとしてＣａｐｏｎ法を用いることを特徴とする請求項１乃至４の何れか１項に記載のレーダ装置。
6. 前記超解像レンジ圧縮手段において、超解像アルゴリズムとしてＭＵＳＩＣ法を用いることを特徴とする請求項１乃至４の何れか１項に記載のレーダ装置。

Images