JP2020125979A - 目標検知装置、目標推定方法およびプログラム - Google Patents

Abstract

【課題】レーダにおいて、複数の目標の検出が高い精度で行えるようにする。【解決手段】送信信号の反射波を受信して得た受信信号から、距離方向の情報についての処理を行って、距離方向（ファストタイム方向）の次元と、パルスヒット方向と周波数方向とアンテナ素子方向の少なくともいずれか一つの方向（スロータイム方向）の別の次元とを有する多次元信号を取得する。そして、多次元信号について、距離方向の次元以外の少なくとも１つ以上の次元の信号を処理し、その処理された信号から、目標速度と目標距離と角度の少なくともいずれか１つを検出する。さらに、検出した目標から、距離ビン毎の信号複素振幅を求め、検出した目標の情報と、距離ビン毎の信号複素振幅をもとに、減算信号を生成し、生成した減算信号を、多次元信号から減算して、別の目標を検出するようにした。【選択図】図１

Description

本発明は、送信信号の反射波を受信して目標を検知するレーダ技術を適用した目標検知装置、目標推定方法およびプログラムに関する。

近年、ミリ波帯の利用による車載レーダ装置やインフラとして設置されたレーダ装置の高分解能化や高信頼性化が求められている。すなわち、既存の車載レーダ装置であっても、車両などの比較的大きな物体の検知は可能であるが、歩行者やさらに小さな対象物の検知を可能とする必要があり、ミリ波レーダのさらなる高分解能化や高信頼性化が求められている。

特許文献１には、レーダ装置において、信号強度の強いターゲットの検出結果から、それを再現する受信信号を推定し、元の受信信号から推定した受信信号を減算して、減算で残った信号に対して検出を繰り返すことで、微弱な信号のターゲットについての検出を可能とする技術についての記載がある。
特許文献１に記載された減算処理を繰り返すことで、複数の対象物を検知することが可能になる。

特開２０１５−４９０７４号公報

ところで、レーダ装置において、受信信号から速度や距離を検出する際には、アナログ系の受信回路が備えるアナログフィルタの特性やケーブルの特性に受信信号が大きく影響を受けてしまう。
一方、元の受信信号から推定した受信信号を減算する際の減算信号は、速度や距離の検出部で得られたピークとなる周波数ビンにおける「速度、距離、複素振幅」から信号である。

したがって、特許文献１に記載された技術では、減算信号を生成して、その減算信号を受信信号から減算したとしても、検出したターゲットの信号を十分に減算したものとはならず、減算残りが生じてしまい、複数のターゲットの検出精度が低下してしまう。特に微弱な信号のターゲットの検出精度が悪くなって、誤検出につながってしまう。

例えば、図１９（Ａ）に示すように、車載レーダの前方に、車とバイクと人が存在した場合、レーダで得られる反射波としては、図１９（Ｂ）に示すように、最も大きなターゲットである車の成分ｍ１と、それより小さなターゲットであるバイクの成分ｍ２と、最も小さなターゲットである人の成分ｍ３とが重なった状態である。
ここで、バイクや人を精度よく検出するためには、車の成分ｍ１やバイクの成分ｍ２を正確に検出する必要があるが、従来の技術では高精度化に限りがあった。

本発明の目的は、複数の目標の検出が高い精度で行うことができる目標検知装置、目標推定方法およびプログラムを提供することにある。

本発明の目標検知装置は、送信信号の反射波を受信して得た受信信号から、距離方向の情報についての処理を行って、距離方向（実施の形態例で述べるファストタイム処理の処理方向）の次元と、パルスヒット方向と周波数方向とアンテナ素子方向の少なくともいずれか一つの方向の別の次元（実施の形態例で述べるスロータイム処理の処理方向）とを有する多次元信号を取得する取得部と、多次元信号について、距離方向の次元以外の少なくとも１つ以上の次元の信号を処理する多次元信号処理部と、多次元信号処理部で処理された信号を使って、目標速度と目標距離と角度の少なくともいずれか１つを検出する目標検出部と、目標検出部が検出した目標速度と目標距離と角度の少なくともいずれか１つから、距離ビン毎の信号複素振幅を求める距離方向複素信号振幅生成部と、目標速度と目標距離と角度の少なくともいずれか１つと、距離ビン毎の信号複素振幅をもとに、減算信号を生成する減算信号生成部と、減算信号生成部が生成した減算信号を、取得部で得られた多次元信号から減算する減算処理部と、を備え、減算処理部で減算された多次元信号を、多次元信号処理部に供給するようにしたものである。

また、本発明の目標推定方法は、送信信号の反射波を受信する受信処理と、受信処理で得た受信信号から、距離方向の情報についての処理を行って、距離方向の次元（実施の形態例で述べるファストタイム処理の処理方向）と、パルスヒット方向と周波数方向とアンテナ素子方向の少なくともいずれか一つの方向（実施の形態例で述べるスロータイム処理の処理方向）の別の次元とを有する多次元信号を取得する信号を取得する取得処理と、多次元信号について、距離方向の次元以外の少なくとも１つ以上の次元の信号を処理する多次元信号処理と、多次元信号処理により処理された信号を使って、目標速度と目標距離と角度の少なくともいずれか１つを検出する目標検出処理と、目標検出処理により検出した目標速度と目標距離と角度の少なくともいずれか１つから、距離ビン毎の信号複素振幅を求める距離方向複素信号振幅生成処理と、目標速度と目標距離と角度の少なくともいずれか１つと、距離ビン毎の信号複素振幅をもとに、減算信号を生成する減算信号生成処理と、減算信号生成処理により生成した減算信号を、取得部で得られた多次元信号から減算して、多次元信号処理を行うための信号とする減算処理と、を含む。

また、本発明のプログラムは、上記目標推定方法の各処理をステップ化して、コンピュータに実行させるようにしたものである。

本発明によると、信号電力の大きい目標のサイドローブに埋もれるような微弱な電力の目標を検知することができ、複数の目標の正確な検出が可能になる。

本発明の第１の実施の形態例による構成例を示すブロック図である。 本発明の第２の実施の形態例による構成例を示すブロック図である。 本発明の第２の実施の形態例による処理手順の例を示すフローチャートである。 本発明の第３の実施の形態例による構成例を示すブロック図である。 本発明の第４の実施の形態例による構成例を示すブロック図である。 本発明の第４の実施の形態例によるアンテナ構成の例を示す図である。 本発明の第５の実施の形態例による構成例を示すブロック図である。 本発明の第５の実施の形態例による処理手順の例を示すフローチャートである。 本発明の第６の実施の形態例による構成例を示すブロック図である。 本発明の第７の実施の形態例による構成例を示すブロック図である。 本発明の第８の実施の形態例による構成例を示すブロック図である。 本発明の各実施の形態例での処理状態の例を示す図である。 本発明の各実施の形態例でのスロータイム処理の例を示す図である。 本発明の各実施の形態例での多次元信号の例を示す図である。 本発明の各実施の形態例での減算処理の例を示す波形図である。 本発明の各実施の形態例での減算後の多次元信号の例を示す図である。 本発明の各実施の形態例での複数の目標（第１目標、第２目標、第３目標）を検出する例を示す図である。 本発明の各実施の形態例でのスロータイム信号の例を示す図である。 車載レーダが複数の目標を検出する例を示す図である。

以下、本発明の実施の形態例を、添付図面を参照して説明する。以下の説明では、まず図１〜図１１を参照して、各実施の形態例（第１〜第８の実施の形態例）のレーダ装置１０ａ〜１０ｈの基本的な構成を順に説明し、その後、各実施の形態例のレーダ装置１０ａ〜１０ｈの各構成要素が行う処理の詳細を、図１２〜図１８を参照して説明する。

［１．第１の実施の形態例］
図１は、本発明の第１の実施の形態例のレーダ装置１０ａの構成例を示す。
第１の実施の形態例のレーダ装置１０ａは、受信部１１、アナログ／デジタル変換器１２、取得部１３、減算器１４、スロータイム処理部１５、目標検出部１６、距離方向複素信号振幅生成部１７、減算信号生成部１８、減算処理部１９、および送信部２１を備える。

送信部２１は、符号変調したパルス信号を、パルス繰り返し間隔（ＰＲＩ：Pulse Repetition Interval）で繰り返し連続して送信する。パルス信号を繰り返し送信する場合、ここでは、後述する図１８（Ｂ）に示すように、パルス信号を送信する毎に周波数を変化させた、時分割多周波ステップ方式とする。時分割多周波ステップ方式の場合、周期的に周波数を繰り返し変化させる。送信周波数帯域としては、例えば７９ＧＨｚ帯などの高い周波数帯を使用する。

受信部１１は、送信部２１で送信した信号の反射波を受信する受信処理を行う。反射波には、目標で反射した成分が含まれる。ここでの目標には、複数の目標が存在する。

受信部１１で得られた受信信号は、符号変調された時分割多周波ステップ受信信号を受信するため、［距離ビン方向と、距離ビン毎のパルスヒット方向の受信信号］×周波数ステップ数分の信号になる。受信部１１で得た受信信号は、アナログフィルタの影響受けた信号（ファストタイム信号）となっており、複数の目標の信号が重畳されている。
なお、以下の説明で複数の目標のそれぞれを、第１目標、第２目標、第３目標、・・・と個別に述べたとき、第１目標の受信信号が最も受信強度が高く、第２目標、第３目標、・・・と順に受信強度が低下するものとする。

受信部１１で得られた受信信号は、アナログ／デジタル変換器１２でサンプリングされる。アナログ／デジタル変換器１２でサンプリングされた受信信号は、取得部１３に供給される。

取得部１３は、受信信号から、距離方向の情報についての処理を行って、距離方向の次元と、それとは別の次元（ここではパルスビット方向の次元）とを有する多次元信号を取得する。なお、取得部１３は、ファストタイム信号を取得処理するため、ファストタイム信号処理部として機能する。
取得部１３は、パルス圧縮処理部１３１と、パルスドップラ処理部１３２とを備える。
パルス圧縮処理部１３１は、受信信号をパルス圧縮処理する。すなわち、パルス圧縮処理部１３１は、符号変調された多周波ステップ受信信号をパルス波形化する。パルス波形化した信号は、アナログフィルタの影響を受けた信号である。
パルス圧縮処理部１３１でパルス波形化した多周波ステップ受信信号は、パルスドップラ処理部１３２に供給される。

パルスドップラ処理部１３２は、パルスドップラフィルタ処理により、距離ビン毎のパルスヒット方向の受信信号をドップラ信号に変換する処理が行われる。ここでのドップラ信号は、目標速度に対応した信号であり、距離方向の次元と速度方向の次元（パルスビット方向の次元）とを有する２次元信号（多次元信号）が生成される。

取得部１３で得られた２次元信号は、減算器１４に供給される。減算器１４では、減算処理部１９から減算信号が供給されるとき、取得部１３から供給される２次元信号から、減算信号を減算処理する。減算信号が供給されないタイミングでは、減算器１４での減算は行われない。減算器１４で減算処理が行われた２次元信号（または減算処理が行われていない２次元信号）は、スロータイム処理部１５に供給される。
なお、減算器１４は、減算後の信号を記憶する記憶部を備え、後述する繰り返しでの減算時に、記憶部が記憶した信号から減算する処理が行われる。

スロータイム処理部１５は、供給される多次元信号（２次元信号）の内で、距離方向の次元以外の少なくとも１つの次元の信号を処理する多次元信号処理部である。スロータイム処理部１５の具体例については、第２〜第４の実施の形態例で説明する。

スロータイム処理部１５で、距離方向の次元以外の次元であるパルスビット方向の次元について処理された信号は、目標検出部１６に供給される。目標検出部１６では、スロータイム処理部１５が出力する２次元信号から、信号振幅が最大となる目標速度と目標距離を検出する。検出した目標速度と目標距離の情報は、出力部１６１に得られる。

そして、目標検出部１６が検出した目標距離の情報は、距離ビン毎の信号複素振幅を求める距離方向複素信号振幅生成部１７に供給される。
距離方向複素振幅生成部１７では、目標速度における距離ビン毎の信号複素振幅を求める距離方向複素信号振幅生成処理を行う。例えば、目標検出部１６が第１目標の目標速度と目標距離を検出したとき、その第１目標の目標速度におけるアナログフィルタの影響受けた複素信号波形を求める。

距離方向複素振幅生成部１７で求めた、目標速度における距離ビン毎の信号複素振幅は、減算信号生成部１８に供給される。
減算信号生成部１８は、目標検出部１６で得られた目標速度および目標距離と、距離方向信複素振幅生成部１７で得られた距離ビン毎の信号複素振幅をもとにパルスヒット方向に亘る減算信号（例えば第１目標の信号）を生成する減算信号生成処理を行う。減算信号生成部１８が生成した減算信号は、減算処理部１９に供給される。

減算処理部１９は、減算器１４で、取得部１３から供給される多次元信号（２次元信号）から、減算信号生成部１８で生成された減算信号を減算する減算処理を行う。
この減算処理を行うことで、取得部１３から出力されたアナログフィルタの影響受けた多次元信号（２次元信号）から、アナログフィルタの影響が考慮された減算信号（例えば最初の段階では第１目標の信号）を減算する処理が行われる。したがって、減算器１４で減算した後の信号には、第２目標以降の目標の成分が含まれることになる。

減算器１４で減算された信号は、第１目標を検出する場合と同様に、スロータイム処理部１５、目標検出部１６で処理され、第２目標についての目標距離、目標速度が検出される。さらに、第３目標が存在する場合には、距離方向複素信号振幅生成部１７と減算信号生成部１８と減算処理部１９での処理が繰り返され、減算器１４で、第２目標以降の目標の成分が含まれる信号（減算器１４が記憶した信号）から、第２目標の信号が減算され、第３目標以降の目標の成分が含まれる信号が減算器１４で得られる。

このようにして、減算器１４、スロータイム処理部１５、目標検出部１６、距離方向信複素振幅生成部１７、減算信号生成部１８、および減算処理部１９での処理が、必要な目標が検出されるまで繰り返される。

以上説明したように、本実施の形態例のレーダ装置１０ａによると、受信信号に含まれる複数の目標（第１目標、第２目標、第３目標、・・・）を、個別に検出することができるようになる。すなわち、振幅の大きな目標に埋もれた振幅の小さな目標の検知が高精度に行うことが可能となり、誤検知を抑制できるようになる。

［２．第２の実施の形態例］
図２は、本発明の第２の実施の形態例のレーダ装置１０ｂの構成例を示す。
第２の実施の形態例のレーダ装置１０ｂは、第２の実施の形態例のレーダ装置１０ａのスロータイム処理部１５として、合成帯域処理部１５ａを適用したものである。

合成帯域処理部１５ａでは、取得部１３から減算器１４を経由して供給される多周波ステップ受信信号を合成して、広帯域化する処理が行われる。合成帯域処理部１５ａで広帯域化した信号は、目標検出部１６に供給される。
図２に示すレーダ装置１０ｂのその他の箇所は、図１に示すレーダ装置１０ａと同じ構成である。

図３は、合成帯域処理部１５ａを備えたレーダ装置１０ｂにおける目標検出処理の流れを示すフローチャートである。この図３のフローチャートでは、レーダ装置１０ｂにおける、減算器１４以降の処理を示す。

まず、減算処理部１９が目標数を更新する（ステップＳ１１）。
そして、減算処理部１９で得られた減算信号により、減算器１４で減算処理を行う（ステップＳ１２）。ここでの減算処理は、ステップＳ１１で設定した目標数が、１以上であるときに行われ、目標数１のときには、ここでの減算処理が省略される。

その後、合成帯域処理部１５ａにより、多周波ステップ受信信号を合成して広帯域化する合成帯域処理が行われる（ステップＳ１３）。さらに、合成帯域処理が行われた速度方向と距離方向の２次元信号から、目標検出部１６で、信号振幅が最大となる目標速度と目標距離の検出処理が行われる（ステップＳ１４）。

そして、距離方向複素振幅生成部１７が、目標速度における距離ビン毎の信号複素振幅を算出する（ステップＳ１５）。この距離方向複素振幅生成部１７で算出された目標速度における距離ビン毎の信号複素振幅と、目標検出部１６で得られた目標速度および目標距離とをもとに、減算信号生成部１８が、パルスヒット方向に亘る減算信号を生成し、生成した減算信号を減算器１４で減算する（ステップＳ１６）。
その後、減算器１４は、減算信号を減算して得られた信号を記憶する（ステップＳ１７）。

そして、減算処理部１９が終了判定の条件を満たすか否かを判断する（ステップＳ１８）。ここでの終了判定の条件を満たす場合としては、例えば減算後の信号波形が、雑音レベルであると判断したとき、終了と判定する。この雑音レベルの判定の一例については、［数８］式および［数９］式を用いて後述する。
このステップＳ１８での判断で、終了判定の条件を満たしていない場合（ステップＳ１８のＮＯ）、ステップＳ１１の目標数の更新に戻る。

また、ステップＳ１８での判断で、終了判定の条件を満たしたと判断したとき（ステップＳ１８のＹＥＳ）、取得部１３で得られた信号についての処理を終了する。

本発明の第２の実施の形態例のレーダ装置１０ｂによると、スロータイム処理部１５として、合成帯域処理部１５ａを適用することで、広帯域化により目標距離の検出精度が向上する。すなわち、複数ステップの受信信号を合成することにより、目標検出部１６に得られる信号のＳＮ比が向上し、目標距離の検出精度が向上する。

ここで、第２の実施の形態例のレーダ装置１０ｂで行われる処理を、数式を用いて説明する。
なお、本明細書で説明する各式に用いられる変数は、以下の通りである。
ｎ：周波数ステップ番号、ｍ：シーケンス繰り返し番号、ｓ：距離ビン番号、ｌ：ドップラビン番号、ｋ：レンジサンプル番号、ｈ：アレー素子番号、ｑ：角度番号、tgtn：目標番号、tgtＮ：目標数、ｃ：光速、ｆ_０：搬送波周波数、ｆ_ｓ：サンプリング周波数、Ｎ：周波数ステップ数、Δｆ：周波数ステップ幅、Ｔ_ｐｒｉ：パルス送信間隔、τ：目標までの電波の往復時間、ｆ_ｄ：目標のドップラ周波数、Ｒ：目標の距離、θ：目標の角度、ΔＶ：目標速度探索刻み幅、ΔＲ：目標距離探索刻み幅、Δθ：目標距離探索刻み幅

まず、取得部１３の出力信号を、次の［数１］式に示す。

この取得部１３の出力は、合成帯域処理部１５ａにて、次に示すように、ｎ方向フーリエ変換がなされる。

合成帯域処理部１５ａの出力は、目標検出部１６にて、しきい値処理がなされ、目標数個の速度と距離が紐づけされて検出される。距離方向複素信号振幅演算部１７により、目標検出部１６にて検出された目標速度と目標距離を用いて距離方向の複素信号振幅を下式より生成する。

距離方向複素信号振幅演算部７の出力より、減算信号生成部８により減算信号が生成される。

減算処理部１９では、減算信号生成部１８により生成された減算信号を記憶する。減算器１４では、取得部１３の出力から、記憶した減算信号が減算される。

減算処理部１９では、以下の式による終了判定を実施する。終了判定では、設定値εより小さいとき終了する。

終了判定で、設定値εより小さいという条件を満たさないときは、目標数の更新に戻り、目標数に１を加えたものが新たに設定される。

［３．第３の実施の形態例］
図４は、本発明の第３の実施の形態例の構成例を示す。
第３の実施の形態例のレーダ装置１０ｃは、第１の実施の形態例のレーダ装置１０ａのスロータイム処理部１５として、パルスドップラ処理部１５ｂを適用したものである。また、取得部１３は、パルス圧縮処理部１３１と合成帯域処理部１３２ａとを備える。

パルスドップラ処理部１５ｂでは、取得部１３から減算器１４を経由して供給される多周波ステップ受信信号について、パルスドップラフィルタ処理により、距離ビン毎のパルスヒット方向の受信信号をドップラ信号に変換する処理が行われる。

パルスドップラ処理部１５ｂでパルスドップラフィルタ処理が行われた信号は、目標検出部１６に供給される。
図４に示すレーダ装置１０ｃのその他の箇所は、図１に示すレーダ装置１０ａと同じ構成である。

本発明の第３の実施の形態例のレーダ装置１０ｃによると、スロータイム処理部１５として、パルスドップラ処理部１５ｂを適用することで、複数の目標の検出を精度よく行うことができる。

ここで、第３の実施の形態例のレーダ装置１０ｃで行われる処理を、数式を用いて説明する。
まず、取得部１３の出力信号を、次の［数８］式に示す。

取得部１３の出力は、ドップラパルス処理部１５ｂにて、ｍ方向のフーリエ変換がなされ、以下の［数９］式、［数１０］式で示すことができる。

そして、ドップラパルス処理部１５ｂの出力は、目標検出部１６にて、しきい値処理がなされ、目標数個の速度と距離が紐づけされて検出される。距離方向複素信号振幅演算部１７により、目標検出部１６にて検出された目標速度と目標距離を用いて距離方向の複素信号振幅を以下の［数１１］式より生成する。

距離方向複素信号生成部１７の出力より、減算信号生成部１８により次に示す減算信号が生成される。

減算器１４では、元の取得部１３の出力から減算信号生成部１８の出力が減算される。

なお、減算処理部１９での減算処理の終了判定では、次式の設定値εより小さいとき、終了と判定する。

終了判定で、設定値εより小さいという条件を満たさないときは、目標数の更新に戻り、目標数に１を加えたものが新たに設定される。

［４．第４の実施の形態例］
図５は、本発明の第４の実施の形態例のレーダ装置１０ｄの構成例を示す。
第４の実施の形態例のレーダ装置１０ｄは、第１の実施の形態例のレーダ装置１０ａのスロータイム処理部１５として、ビームフォーミング処理部１５ｃを適用したものである。
また、レーダ装置１０ｄは、図６に示すように、受信部１１に複数の受信アンテナ１１１，１１２，１１３，１１４，・・・が接続されている。受信アンテナ１１１，１１２，１１３，１１４，・・・は、例えば受信角度（方位角）をずらして配置する。

受信部１１に複数の受信アンテナ１１１，１１２，１１３，１１４，・・・が接続されていることで、取得部１３で得られる多次元信号に含まれる１つの次元には、目標の角度の次元の信号が含まれるようになる。

そして、スロータイム処理部１５であるビームフォーミング処理部１５ｃが、ビームフォーミング処理で角度ごとの受信状態の情報を取得する。
ビームフォーミング処理部１５ｃで得られた角度ごとの受信状態の情報は、目標検出部１６に供給される。目標検出部１６では、角度ごとの受信状態の情報から、目標の角度（方位角）を検出する。目標検出部１６が検出した目標の角度の情報は、出力部１６１から出力される。
図５に示すレーダ装置１０ｂのその他の箇所は、図１に示すレーダ装置１０ａと同じ構成である。

本発明の第４の実施の形態例のレーダ装置１０ｄによると、スロータイム処理部１５として、ビームフォーミング処理部１５ｃを適用することで、複数の目標の角度の検出を精度よく行うことができる。

［５．第５の実施の形態例］
図７は、本発明の第５の実施の形態例のレーダ装置１０ｅの構成例を示す。
第５の実施の形態例のレーダ装置１０ｅは、減算処理部１９が、目標数推定部１９１を備えるようにしたものである。
目標数推定部１９１は、レーダ装置１０ｅでの過去や現在の目標の検出状況や、他の装置（例えばカメラなど）での検出状況から、レーダ装置１０ｅの周囲の目標の数を推定する。目標数推定部１９１が推定する目標の数には、上限を決めてもよい。
図７に示すレーダ装置１０ｅのその他の箇所は、図１に示すレーダ装置１０ａと同じ構成である。

図８は、減算処理部１９が、目標数推定部１９１を備えた場合の処理の流れを示すフローチャートである。

まず、減算処理部１９が目標数を更新する（ステップＳ２１）。こ
そして、減算処理部１９で得られた減算信号により、減算器１４で減算処理を行う（ステップＳ２２）。ここでの減算処理は、ステップＳ２１で設定した目標数が、１以上であるときに行われ、目標数１のときには、ここでの減算処理が省略される。

その後、スロータイム処理部１５でのスロータイム処理が行われる（ステップＳ２３）。さらに、スロータイム処理が行われた速度方向と距離方向の２次元信号から、目標検出部１６で、信号振幅が最大となる目標速度と目標距離の検出処理が行われる（ステップＳ２４）。

そして、距離方向複素振幅生成部１７が、目標速度における距離ビン毎の信号複素振幅を算出する（ステップＳ２５）。この距離方向複素振幅生成部１７で算出された目標速度における距離ビン毎の信号複素振幅と、目標検出部１６で得られた目標速度および目標距離とをもとに、減算信号生成部１８が、パルスヒット方向に亘る減算信号を生成し、生成した減算信号を減算器１４で減算する（ステップＳ２６）。
その後、減算器１４は、減算信号を減算して得られた信号を記憶する（ステップＳ２７）。

次に、減算処理部１９は、目標数推定部１９１が推定した目標数に対応した試行回数だけ減算を繰り返したか否かを判断する（ステップＳ２８）。このステップＳ２８での判断で、目標検出の試行回数が、推定した目標数に対応した数でない場合（ステップＳ２８のＮＯ）、ステップＳ２２の減算処理に戻る。
また、ステップＳ２８での判断で、目標検出の試行回数が、推定した目標数に対応した数であるとき（ステップＳ２８のＹＥＳ）、目標検出についての終了条件を満たすか否かを判断する（ステップＳ２９）。
ここで、目標数推定部１９１の推定方法として、目標検出部１９１で検出したピークに対応する波形を減算処理した後のサイドローブレベルと減算前のサイドローブレベルを比較して、有意な差があれば検知したピークを目標とみなし、目標数をカウントアップする。有意な差が見られなければピークを目標と見なさず目標数のカウントアップを行わない。目標検出についての終了条件を満たしたときのカウントを目標数として出力する。

このステップＳ２９で終了条件を満たさない場合（ステップＳ２９のＮＯ）、ステップＳ２１の目標数の更新に戻る。また、ステップＳ２９での判断で、終了条件を満たす場合（ステップＳ２９のＹＥＳ）、取得部１３で得られた信号についての処理を終了する。

本発明の第５の実施の形態例のレーダ装置１０ｅによると、目標数推定部１９１を設けて、目標数推定部１９１が推定した目標数に応じた回数だけ、目標を検出する処理を繰り返すことで、適切な数の目標数を検出できるようになる。

［６．第６の実施の形態例］
図９は、本発明の第６の実施の形態例のレーダ装置１０ｆの構成例を示す。
第６の実施の形態例のレーダ装置１０ｆは、取得部１３ａが、パルス圧縮処理部１３１と、ビームフォーミング処理部１３３と、パルスドップラ処理部１３２とを備えるようにしたものである。
また、レーダ装置１０ｆは、図６に示すように、受信部１１に複数の受信アンテナ１１１，１１２，１１３，１１４，・・・が接続されている。受信アンテナ１１１，１１２，１１３，１１４，・・・は、例えば受信角度（方位角）をずらして配置する。

受信部１１に複数の受信アンテナ１１１，１１２，１１３，１１４，・・・が接続されていることで、取得部１３のビームフォーミング処理部１３３で得られる多次元信号に含まれる１つの次元には、目標の角度の次元の信号が含まれるようになる。ビームフォーミング処理部１３３は、パルス圧縮処理部１３１の出力から、目標の角度の次元の信号を得る。
そして、ビームフォーミング処理部１３３で目標角度の次元の信号が含まれた信号が、パルスドップラ処理部１３２に供給される。
図９に示すレーダ装置１０ｆのその他の箇所は、図１に示すレーダ装置１０ａと同じ構成である。

［７．第７の実施の形態例］
図１０は、本発明の第７の実施の形態例のレーダ装置１０ｇの構成例を示す。
第７の実施の形態例のレーダ装置１０ｇは、スロータイム処理部１５ｄが、パルスドップラ処理部１５１と、帯域合成処理部１５２と、ビームフォーミング処理部１５３とを備えるようにしたものである。ここで、レーダ装置１０ｇのスロータイム処理部１５ｄは、パルスドップラ処理部１５１、帯域合成処理部１５２、ビームフォーミング処理部１５３の順で処理を行うようにした。
また、レーダ装置１０ｇは、図６に示すように、受信部１１に複数の受信アンテナ１１１，１１２，１１３，１１４，・・・が接続されている。受信アンテナ１１１，１１２，１１３，１１４，・・・は、例えば受信角度（方位角）をずらして配置する。なお、取得部１３は、パルス圧縮処理部１３１を備える。

受信部１１に複数の受信アンテナ１１１，１１２，１１３，１１４，・・・が接続されていることで、ビームフォーミング処理部１５３は、目標角度の次元の信号についてのスロータイム処理を行うようになる。
図１０に示すレーダ装置１０ｇのその他の箇所は、図１に示すレーダ装置１０ａと同じ構成である。

ここで、第７の実施の形態例のレーダ装置１０ｇで行われる処理を、数式を用いて説明する。
まず、取得部１３の出力信号を、次の［数１５］式に示す。

この取得部１３の出力は、パルスドップラ処理部１５１にて、ｍ方向のフーリエ変換がなされる。

パルスドップラ処理部１５１の出力は、合成帯域処理部１５２にて、ｎ方向フーリエ変換がなされる。

合成帯域処理部１５２の出力は、目標検出部１６にて、しきい値処理がなされ、目標数個の速度と距離が紐づけされて検出される。そして、距離方向複素信号振幅演算部１７により、目標検出部１６にて検出された目標速度と目標距離を用いて距離方向の複素信号振幅を下式より生成する。

さらに、距離方向複素信号振幅演算部１７の出力より、減算信号生成部１８により減算信号が生成される。

減算処理部１９では、減算信号生成部１８により生成された減算信号を記憶する。
減算器１４では、元の取得部１３の出力から減算処理部１９に記憶された減算信号が減算される。

ここで、減算処理部１９では、目標数の更新により目標数が１として設定される。その後、目標検出の試行回数が設定回数に達すると、以下の式による終了判定を実施する。終了判定では、設定値εより小さいとき終了する。

終了判定で、設定値εより小さいという条件を満たさないときは、目標数の更新に戻り、目標数に１を加えたものが新たに設定される。

［８．第８の実施の形態例］
図１１は、本発明の第８の実施の形態例のレーダ装置１０ｈの構成例を示す。
第８の実施の形態例のレーダ装置１０ｈは、スロータイム処理部１５ｅが、ビームフォーミング処理部１５４と、パルスドップラ処理部１５５と、帯域合成処理部１５６とを備えるようにしたものである。ここで、レーダ装置１０ｈのスロータイム処理部１５ｅは、ビームフォーミング処理部１５４、パルスドップラ処理部１５５、帯域合成処理部１５６の順で処理を行うようにした。なお、取得部１３は、パルス圧縮処理部１３１を備える。
また、レーダ装置１０ｈは、図６に示すように、受信部１１に複数の受信アンテナ１１１，１１２，１１３，１１４，・・・が接続されている。受信アンテナ１１１，１１２，１１３，１１４，・・・は、例えば受信角度（方位角）をずらして配置する。

受信部１１に複数の受信アンテナ１１１，１１２，１１３，１１４，・・・が接続されていることで、ビームフォーミング処理部１５４は、目標角度の次元の信号についてのスロータイム処理を行うようになる。
図１１に示すレーダ装置１０ｈのその他の箇所は、図１に示すレーダ装置１０ａと同じ構成である。

［９．受信部が受信する信号の例］
図１２は、各実施の形態例の受信部１１での受信信号の例を示す。
図１２（Ａ）は、送信信号ＴＸの一例を示す。
図１２（Ａ）に示す例では、パルス繰り返し周期ＰＲＩごとにパルス信号を送信する状態を示す。図１２では、１周期ＰＲＩごとの周波数の相違は無視する。

このような送信信号ＴＸを送信部２１が送信したとき、受信部１１では、図１２（Ｂ）に示すように、送信信号ＴＸの反射波である受信信号ＲＸが得られる。この受信信号ＲＸの受信タイミングは、目標との距離に応じて送信タイミングから遅れたものになる。
図１２（Ｂ）に示す受信信号ＲＸは、目標が１つである場合の例を示し、目標が２つの場合には、図１２（Ｃ）に示す状態になる。すなわち、図１２（Ｃ）に示すように、第１目標に対応した強い受信信号ＲＸ１と、第２目標に対応した弱い受信信号ＲＸ２とが重畳された状態になる。第１目標と第２目標とに距離が相違する場合、それぞれの信号ＲＸ１，ＲＸ２の受信タイミングが異なる。

この図１２（Ｃ）に示す受信信号ＲＸ１，ＲＸ２が受信される状態は理想的な場合を示し、実際には受信部１１でのアナログフィルタ処理により、図１２（Ｄ）に示すようなアナログ的に強度が変化する受信信号ＲＸａの波形が得られる。
本発明の各実施の形態例のレーダ装置１０ａ〜１０ｈでは、図１２（Ｄ）に示すような受信信号ＲＸａから、第１目標の成分と第２目標の成分を取りだして、それぞれの目標を高精度に検出する処理が行われる。

［１０．ファストタイム処理とスロータイム処理の例］
図１３は、本発明の各実施の形態例で行われるファストタイム処理とスロータイム処理の概要を示す。
各実施の形態例では、取得部１３が、パルス圧縮処理などのファストタイム処理を行い、スロータイム処理部１５が、合成帯域処理などのスロータイム処理を行うようにした。
まず、受信部１１で図１３（Ａ）に示す繰り返し間隔ＰＲＩごとの受信信号ＲＸａが得られるとする。この受信信号は、図１２（Ｄ）で説明した受信信号ＲＸａと同じ波形である。この受信信号ＲＸａには、図１２で説明した第１目標の成分と第２目標の成分が含まれている。

ここで、取得部１３では、パルス圧縮処理を行う際に、図１３（Ｂ）に示す参照信号を取得し、受信信号をサンプリングした信号と、参照信号との相互相関を取る。ここでの参照信号は、図１２（Ａ）に示す送信信号に対応する。
参照信号との相互相関を取ることで、図１３（Ｃ）に示すパルス圧縮（相互相関）された信号Ｐ１が、パルス繰り返し周期ＰＲＩ毎（ファストタイム方向）に得られる。このパルス圧縮された信号Ｐ１は、パルスの終了時刻にピークを持つ波形である。
なお、パルス圧縮された信号Ｐ１には、図１３（Ｃ）に破線で示す第２目標の信号Ｐ２が含まれるが、信号Ｐ１に埋もれている。
ここまでのパルス圧縮処理は、本実施の形態例で述べるファストタイム処理である。

一方、スロータイム処理部１５でのスロータイム処理では、複数の繰り返し間隔ＰＲＩでの信号についての処理が行われる。すなわち、図１３（Ｄ）に示すように、それぞれの間隔ＰＲＩの信号の内で、同一レンジビンに対して、パルスビット方向にフーリエ変換する処理が行われる。図１３（Ｄ）に示す状態は、パルス圧縮された信号Ｐ１のピーク位置の各ＰＲＩのレンジビン（第１目標が得られるレンジビン）ｆａに対して、パルスビット方向にフーリエ変換する状態を示す。また、図１３（Ｅ）に示す状態は、パルス圧縮された信号Ｐ１の第２目標が得られる箇所のレンジビンｆｂに対して、パルスビット方向にフーリエ変換する状態を示す。

図１４は、ファストタイム処理後のレンジビンの信号（横軸）と、スロータイム処理後のドップラビンの信号（縦軸）とを、２次元に展開したものである。
減算器１４での減算が行われていない最初の状態では、図１４（Ａ）に示すように、レンジビンとドップラビンのいずれの方向でも、第１目標の位置が最も高く、その第１目標の位置から徐々にレベルが低下する状態を示す。
図１４（Ｂ）は、スロータイム処理後のドップラビンの方向に、フーリエ変換した状態を矢印で示す。例えば、最も高いレベルのレンジビンｆａが、第１目標が含まれる信号になる。

そして、距離方向複素信号振幅生成部１７での距離方向複素信号振幅の生成で、レンジビン方向に減算波形を生成することで、図１５（Ａ）に示すように、第１目標についての減算波形α１が生成される。この減算波形α１を元の波形から減算することで、図１５（Ｂ）に示すように、第１目標の成分が除去された第２目標（およびそれより信号強度が弱い目標）の成分の波形α２が得られる。
この減算後の第２目標の成分の波形α２から、距離、速度、角度などの検出が可能になる。

図１６は、図１５（Ｂ）に示す減算後の信号から、第２目標を検出した例を示す。図１６は、図１４と同様に、ファストタイム処理後のレンジビンの信号（横軸）と、スロータイム処理後のドップラビンの信号（縦軸）とを、２次元に展開したものである。
この図１６に示すように、減算後の信号からは、第２目標の信号がピークとして明確に現れ、第２目標についての情報（速度、距離、方位など）を高精度に検出できるようになる。

図１７は、速度と距離に展開した信号（左側）と、信号の電力（右側）とを示す。図１７（Ａ）は、受信信号を示す。図１７（Ａ）に示すように、受信信号には、第１目標ｍ１の情報と、第２目標の情報ｍ２と、第３目標の情報ｍ３とが含まれている。
図１７（Ｂ）は、この受信信号から第１目標を得る信号を示す。第１目標を得る際には減算処理を行っていないので、図１７（Ａ）の信号と図１７（Ｂ）の信号は同じであり、第１目標ｍ１がピークになる。このピーク値より第１目標を検出することができる。そして、第１目標の目標速度における距離ビン毎の信号複素振幅を求める処理が行われる。

図１７（Ｃ）は、１回目の減算処理を行った信号を示す。すなわち、第１目標を検出した際の目標速度における距離ビン毎の信号複素振幅から、第１目標についての減算信号を生成し、減算処理を行うことで、振幅が小さな第２目標と第３目標が重畳した信号が、図１７（Ｃ）に示すように残る。この信号では第２目標ｍ２がピークになり、第２目標ｍ２の距離や速度を検出することができる。そして、第２目標の目標速度における距離ビン毎の信号複素振幅を求める処理が行われる。

図１７（Ｄ）は、２回目の減算処理を行った信号を示す。すなわち、第２目標を検出した際の目標速度における距離ビン毎の信号複素振幅から、第２目標についての減算信号を生成し、減算処理を行うことで、振幅が小さな第３目標の信号が、図１７（Ｄ）に示すように残る。この信号では第３目標ｍ３がピークになり、第３目標ｍ３の距離や速度を検出することができる。
さらに第４目標以降が存在する場合には、減算処理を繰り返す。

［１１．スロータイム信号の例］
上述した各実施の形態例では、符号変調された時分割多周波ステップ信号を受信して、ファストタイム処理とスロータイム処理を行うようにした。
時分割多周波ステップ信号について、スロータイム処理を行うのは一例であり、その他の受信信号を、スロータイム信号として処理することもできる。

例えば、図１８（Ａ）は、符号変調したパルス信号を、同じ周波数を使って、パルス繰り返し間隔（ＰＲＩ）で繰り返し送信する場合に、一定数の繰り返し間隔で、スロータイム信号を生成する例を示す。

図１８（Ｂ）は、符号変調したパルス信号を、順に周波数を変化させて、パルス繰り返し間隔（ＰＲＩ）で繰り返し送信する場合に、その周波数変化を繰り返す単位で、スロータイム信号を生成する例を示す。つまり、図１８（Ｂ）は、時分割多周波ステップ信号の例である。

この時分割多周波ステップ信号の場合、図１８（Ｃ）に示すように、例えば周波数がＮ段階（Ｎは２以上の任意の整数）の変化であり、周波数ｆ_０からｆ_Ｎ−１まで変化するものを、ｍ＝０のスロータイム信号とし、以下、同じ周期で、ｍ＝１、ｍ＝２、・・・とスロータイム信号を生成する。
なお、図１８（Ｂ）や図１８（Ｃ）に示す例では、周波数を順に変化させる例を示したが、図示の例とは異なる順序で周波数を変化させてもよい。

図１８（Ｄ）は、複数の受信アンテナによるアンテナアレイを備えて、その複数の受信アンテナの受信信号をビームフォーミング処理する場合におけるスロータイム信号の例を示す。この例は、図５および図６に示す第４の実施の形態例に相当する。
それぞれのアンテナの受信信号がパルス信号であるとき、全てのアンテナの受信信号で、スロータイム信号を生成する。この場合のスロータイム信号では、距離や速度の情報の他に、目標の方向（方位角）の情報を取得することができる。

［１２．具体的な信号処理の例］
次に、第１の実施の形態例の構成で、信号処理を行う場合の具体的な信号処理の詳細を、説明する。
まず、受信部１１で受信される信号は、送信部２１から送信した信号が目標で反射し、目標までの往復時間に相当する時間遅延τの後、受信波として受信アンテナに入射し、受信部１１で受信され、アナログ／デジタル変換器１２でサンプリングされる。
サンプリングされた受信波は、［数２４］式で示される。

ここで、ｔ_ｎ，ｍは周波数ステップ番号ｎのｍ回目の時刻、ｓは距離ゲート番号、codeは、ＣＰＣ符号列番号（code=0,1）、ｆ_ｄ，ｎは周波数ステップ番号ｎの周波数におけるドップラ周波数である。
複数目標が存在する環境下におけるサンプリング信号は、［数２４］式の線形和として表される。
アナログ／デジタル変換器１２でサンプリングされた信号は、パルス圧縮処理部１３１でパルス圧縮処理され、その出力は、次の［数２５］式で示される。

ここで、sinc_ｃｏｄｅは各ＣＰＣ符号の自己相関関数を示す。パルス圧縮出力は、パルスドップラ処理部１３２のパルスドップラフィルタ（ＰＤＦ）処理により、周波数ステップごとにｍ方向のフーリエ変換が施され、その出力は、次の［数２６］式で示される。

フーリエ変換が施された信号は、加算処理により、各codeのパルスドップラフィルタ出力が加算処理される。加算処理された出力は、次の［数２７］式で示される。

そして、パルスドップラ処理部１３２の出力は、スロータイム処理部１５に供給され、スロータイム処理が行われる。スロータイム処理として、第２の実施の形態例のように合成帯域処理部１５ａでの合成帯域処理を行った場合、ｎ方向のフーリエ変換が施される。
ｎ方向のフーリエ変換が施された信号は、次の［数２８］式で示される。

合成帯域処理部１５ａの出力は，目標検出部１６にて、しきい値処理がなされ、最大電力に対応する速度と距離が紐づけされて検出される。
そして、減算信号生成部１８において、減算する目標数がＫのときのパルス圧縮後の距離ビン毎の複素振幅は、次の［数２９］式で示される。

ここで、tgtn（＝０，１，…Ｋ−１）は、減算する目標番号を示す。距離ビン毎の距離方向複素信号振幅生成部１７の出力を用いて、距離ビン毎の減算信号波形が、次の［数３０］式で生成される。

この減算信号の減算処理により、元のパルス圧縮出力から減算信号を減算する。減算後の波形は、次の［数３１］式で示される。

ここで、減算処理部１９は、例えば収束判定により、Ｖ＾およびＲ＾（＾は本来は数式に示すように符号の上に付与される）が収束するまで、減算処理を繰り返す。収束したら、目標／サイドローブ判定に移行する。例えば、減算する目標数がＫのときと、Ｋ−１のときの合成帯域後のサイドローブレベルを比較して、サイドローブレベルに低減が見られないときは、最後に検知した目標はサイドローブの誤検知であるとして、検知対象から除外し、次に大きなピーク電力をＫ番目の目標として検知して処理を続ける。サイドローブレベルに低減が確認されたら、目標であるとして、終了判定に移行する。

例えば、［数３１］式の減算後の波形が、［数３２］式に示すように雑音レベルになったとき、全体の処理を終了する。

減算後の波形が雑音レベルにならないとき、減算目標数更新により減算する目標に新たに推定した目標を加え、減算する目標数をＫ＋１として、減算処理を繰り返す。
このような終了判定処理は、例えば図３のフローチャートのステップＳ１８や、図８のフローチャートのステップＳ２９での終了判定に相当する。
なお、このような雑音レベルの判断に基づいて終了判定を行うのは一例であり、その他の処理で終了判定してもよい。

［１３．変形例］
なお、上述した各実施の形態例で示す構成は、それぞれの構成を行う処理回路などのハードウェアを設けて、専用の回路で構成してもよいが、処理回路の一部または全てを、コンピュータで構成してもよい。コンピュータで構成する場合には、それぞれの処理を行うプログラムを用意し、そのプログラムのコンピュータへの実装で、各処理回路での処理が実現される。
ここでのコンピュータには、ＣＰＵ（Central Processing Unit：中央処理装置）とその周辺機器をバスで接続した一般的な計算機の他に、ＦＰＧＡ（Field Programmable Gate Array）やＡＳＩＣ（Application Specific Integrated Circuit）等の半導体チップで構成されたコンピュータによって実現してもよい。

１０ａ，１０ｂ，１０ｃ，１０ｄ，１０ｅ，１０ｆ，１０ｇ，１０ｈ…レーダ装置、１１…受信部、１２…アナログ／デジタル変換器、１３…取得部、１４…減算器、１５…スロータイム処理部、１５ａ…合成帯域処理部、１５ｂ…パルスドップラ処理部、１５ｃ…ビームフォーミング処理部、１５ｄ，１５ｅ…スロータイム処理部、１６…目標検出部、１７…距離方向複素信号振幅生成部、１８…減算信号生成部、１９…減算処理部、２１…送信部、１１１，１１２，１１３，１１４…アンテナ、１３１…パルス圧縮処理部、１３２…パルスドップラ処理部、１３２ａ…帯域合成処理部、１３３…ビームフォーミング処理部、１５１…パルスドップラ処理部、１５２…帯域合成処理部、１５３…ビームフォーミング処理部、１５４…ビームフォーミング処理部、１５５…パルスドップラ処理部、１５６…帯域合成処理部、１６１…目標検出信号出力部、１９１…目標数推定部

Claims (8)

1. 送信信号の反射波を受信して得た受信信号から、距離方向の情報についての処理を行って、距離方向の次元と、パルスヒット方向と周波数方向とアンテナ素子方向の少なくともいずれか一つの方向の別の次元とを有する多次元信号を取得する取得部と、
   前記多次元信号について、前記距離方向の次元以外の少なくとも１つ以上の次元の信号を処理する多次元信号処理部と、
   前記多次元信号処理部で処理された信号を使って、目標速度と目標距離と角度の少なくともいずれか１つを検出する目標検出部と、
   前記目標検出部が検出した前記目標速度と前記目標距離と前記角度の少なくともいずれか１つから、距離ビン毎の信号複素振幅を求める距離方向複素信号振幅生成部と、
   前記目標速度と前記目標距離と前記角度の少なくともいずれか１つと、前記距離ビン毎の信号複素振幅をもとに、減算信号を生成する減算信号生成部と、
   前記減算信号生成部が生成した減算信号を、前記取得部で得られた多次元信号から減算する減算処理部と、を備え、
   前記減算処理部で減算された多次元信号を、前記多次元信号処理部に供給するようにした
   目標検知装置。
2. 前記送信信号はパルス信号であり、
   前記取得部は、パルス圧縮処理部である
   請求項１に記載の目標検知装置。
3. 前記多次元信号処理部は、パルスドップラ処理を行うパルスドップラ処理部であり、
   前記目標検出部は、前記パルスドップラ処理部の出力に基づいて、目標速度と目標距離と角度の少なくともいずれか１つを検出する
   請求項２に記載の目標検知装置。
4. 前記多次元信号処理部は、合成帯域処理を行う合成帯域処理部であり、
   前記目標検出部は、前記合成帯域処理部の出力に基づいて、目標速度と目標距離と角度の少なくともいずれか１つ検出する
   請求項２に記載の目標検知装置。
5. 前記多次元信号処理部は、ビームフォーミング処理を行うビームフォーミング処理部であり、
   前記目標検出部は、前記ビームフォーミング処理部の出力に基づいて、目標速度と目標距離と角度の少なくともいずれか１つ検出する
   請求項２に記載の目標検知装置。
6. 前記減算処理部は、目標数推定部を含む
   請求項２に記載の目標検知装置。
7. 送信信号の反射波を受信する受信処理と、
   前記受信処理で得た受信信号から、距離方向の情報についての処理を行って、距離方向の次元と、パルスヒット方向と周波数方向とアンテナ素子方向の少なくともいずれか一つの方向の別の次元とを有する多次元信号を取得する信号を取得する取得処理と、
   前記多次元信号について、前記距離方向の次元以外の少なくとも１つ以上の次元の信号を処理する多次元信号処理と、
   前記多次元信号処理により処理された信号を使って、目標速度と目標距離と角度の少なくともいずれか１つを検出する目標検出処理と、
   前記目標検出処理により検出した前記目標速度と前記目標距離と前記角度の少なくともいずれか１つから、距離ビン毎の信号複素振幅を求める距離方向複素信号振幅生成処理と、
   前記目標速度と前記目標距離と前記角度の少なくともいずれか１つと、前記距離ビン毎の信号複素振幅をもとに、減算信号を生成する減算信号生成処理と、
   前記減算信号生成処理により生成した減算信号を、取得処理で得られた多次元信号から減算して、前記多次元信号処理を行うための信号とする減算処理と、を含む
   目標推定方法。
8. 送信信号の反射波を受信して得た受信信号から、距離方向の情報についての処理を行って、距離方向の次元と、パルスヒット方向と周波数方向とアンテナ素子方向の少なくともいずれか一つの方向の別の次元とを有する多次元信号を取得する取得ステップと、
   前記多次元信号について、前記距離方向の次元以外の少なくとも１つ以上の次元の信号を処理する多次元信号処理ステップと、
   前記多次元信号処理ステップで処理された信号を使って、目標速度と目標距離と角度の少なくともいずれか１つを検出する目標検出ステップと、
   前記目標検出ステップにより検出した前記目標速度と前記目標距離と前記角度の少なくともいずれか１つから、距離ビン毎の信号複素振幅を求める距離方向複素信号振幅生成ステップと、
   前記目標速度と前記目標距離と前記角度の少なくともいずれか１つと、前記距離ビン毎の信号複素振幅をもとに、減算信号を生成する減算信号生成ステップと、
   前記減算信号生成ステップにより生成した減算信号を、前記取得ステップで得られた多次元信号から減算して、前記多次元信号処理ステップでの処理を行うための信号とする減算ステップと、
   をコンピュータに実行させる
   プログラム。