(57)【特許請求の範囲】
【請求項１】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された推定値に対して、上記フィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算する第一のフィルタ係数計算手段と、
上記縦続接続されたフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された上記クラッタの中心周波数の推定値に対して、オフセット周波数を持たせて自フィルタの零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算手段と、
該第二〜第Ｎのフィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎのフィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタ手段と、
上記受信信号及び上記第一〜第Ｎのフィルタ手段の出力信号から、上記第一〜第Ｎのフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項２】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
該フィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段から送られるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタ手段と、
上記受信信号、又は上記フィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記フィルタ手段へ出力する切替手段と、
該切替手段へ送られる上記フィルタ手段の出力信号を記憶しておく出力信号記憶手段と、
上記受信信号及び上記フィルタ手段の出力信号から、上記フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項３】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
該フィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタ手段と、
上記受信信号、又は上記第１のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第１のフィルタ手段へ出力する切替手段と、
上記フィルタ係数計算手段から出力されるフィルタ係数を記憶しておくフィルタ係数記憶手段と、
上記受信信号及び上記第１のフィルタ手段の出力信号から、上記第１のフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタ手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項４】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された推定値に対して、上記フィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算する第一のフィルタ係数計算手段と、
上記従続接続されたフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された上記クラッタの中心周波数の推定値に対して、オフセット周波数を持たせて自フィルタの零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算手段と、
上記縦続接続されたフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された上記クラッタの中心周波数の推定値に対して、オフセット周波数を持たせて自フィルタの零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算手段と、
該第二〜第Ｎのフィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記受信信号中の複数のクラッタに対してフィルタのノッチが形成され、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記オフセット周波数に記憶されているオフセット周波数を設定するオフセット周波数選択手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎのフィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタ手段と、
上記受信信号及び上記第一〜第Ｎのフィルタ手段の出力信号から、上記第一〜第Ｎのフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項５】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
上記受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記フィルタ係数計算手段へ送られるオフセット周波数を選択するオフセット周波数設定手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタ手段と、
上記受信信号、又は上記フィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記フィルタ手段へ出力する切替手段と、
該切替手段に送られる上記フィルタ手段の出力信号を記憶しておく出力信号記憶手段と、
上記受信信号及び上記フィルタ手段の出力信号から、上記フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項６】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
上記受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記フィルタ係数計算手段へ送られるオフセット周波数を選択するオフセット周波数設定手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタ手段と、
上記受信信号、又は上記第１のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第１のフィルタ手段へ出力する切替手段と、
上記フィルタ係数計算手段から出力されるフィルタ係数を記憶しておくフィルタ係数記憶手段と、
上記受信信号及び上記第１のフィルタ手段の出力信号から、上記第１のフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタ手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項７】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する多次フィルタ手段と、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算手段へ送られた以外の推定値に対して、上記フィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算する第一のフィルタ係数計算手段と、
上記縦続接続されたフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定されたクラッタ中心周波数の推定値うち、上記多次フィルタ係数計算手段へ送られた以外の推定値に対して、オフセット周波数を持たせて自フィルタの零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算手段と、
該第二〜第Ｎのフィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎのフィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタ手段と、
上記受信信号及び上記第一〜第Ｎのフィルタ手段の出力信号から、上記第一〜第Ｎのフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
上記フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項８】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一のフィルタ手段と、
上記クラッタ中心周波数推定手段により推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算手段へ送られた以外の推定値に対して、オフセット周波数を持たせて自フィルタの零点を割り当てるようにフィルタ係数を計算するフィルタ係数計算手段と、
該フィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段から出力されるフィルタ係数を用いて切替手段から出力される信号中のクラッタを抑圧する第二のフィルタ手段と、
上記第一のフィルタ手段の出力信号、又は上記第二のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第二のフィルタ手段へ出力する切替手段と、
該切替手段へ送られる上記第二のフィルタ手段の出力信号を記憶しておく出力信号記憶手段と、
上記第一及び第二のフィルタ手段の出力信号から、フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、  該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項９】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一のフィルタ手段と、
上記クラッタ中心周波数推定手段により推定された推定値のうち、上記多次フィルタ係数計算手段へ送られた推定値以外の推定値に対して、オフセット周波数を持たせて自フィルタの零点を割り当てるようにフィルタ係数を計算するフィルタ係数計算手段と、
該フィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段から出力されるフィルタ係数を用いて切替手段から出力される信号中のクラッタを抑圧する第二のフィルタ手段と、
上記第一のフィルタ手段の出力信号、又は上記第二のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第二のフィルタ手段へ出力する切替手段と、
上記フィルタ係数計算手段から出力されるフィルタ係数を記憶しておくフィルタ係数記憶手段と、
上記第一及び第二のフィルタ手段の出力信号から、フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、  該フィルタ次数選定手段から出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行うフィルタ手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項１０】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
上記クラッタ中心周波数推定ステップにより推定された推定値に対して、複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算する第一のフィルタ係数計算ステップと、
上記クラッタ中心周波数推定ステップにより推定された上記クラッタの中心周波数の推定値に対して、オフセット周波数を持たせて上記フィルタ群の零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算ステップと、
該第二〜第Ｎのフィルタ係数計算ステップにおいて使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎのフィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタステップと、
上記受信信号及び上記第一〜第Ｎのフィルタステップの出力信号から、上記第一〜第Ｎのフィルタステップにおいて入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項１１】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
該フィルタ係数計算ステップで使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップにより出力されるフィルタ係数を用いて受信信号中のクラッタを抑圧するフィルタステップと、
上記受信信号、又は上記フィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記フィルタステップへ出力する切替ステップと、
上記切替ステップにおいて使用される上記フィルタステップの出力信号を記憶しておく出力信号記憶ステップと、
上記受信信号及び上記フィルタステップの出力信号から、上記フィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項１２】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
該フィルタ係数計算ステップにおいて使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタステップと、
上記受信信号、又は上記第１のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記フィルタステップへ出力する切替ステップと、
上記フィルタ係数計算ステップにより出力されるフィルタ係数を記憶しておくフィルタ係数記憶ステップと、
上記受信信号及び上記第１のフィルタステップの出力信号から、上記第１のフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項１３】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
上記クラッタ中心周波数推定ステップにより推定された推定値に対して、複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算する第一のフィルタ係数計算ステップと、
上記クラッタ中心周波数推定ステップにより推定された上記クラッタの中心周波数の推定値に対して、オフセット周波数を持たせて上記フィルタ群の零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算ステップと、
該第二〜第Ｎのフィルタ係数計算ステップにおいて使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
上記受信信号中の複数のクラッタに対してフィルタのノッチが形成され、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記オフセット周波数に記憶されているオフセット周波数を設定するオフセット周波数選択ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎのフィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタステップと、
上記受信信号及び上記第一〜第Ｎのフィルタステップの出力信号から、上記第一〜第Ｎのフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項１４】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
上記受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記フィルタ係数計算ステップにおいて使用されるオフセット周波数を選択するオフセット周波数設定ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタステップと、
上記受信信号、又は上記フィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記フィルタステップへ出力する切替ステップと、
該切替ステップにおいて使用される上記フィルタステップの出力信号を記憶しておく出力信号記憶ステップと、
上記受信信号及び上記フィルタステップによる出力信号から、上記フィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項１５】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
上記受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記フィルタ係数計算ステップにおいて使用されるオフセット周波数を選択するオフセット周波数設定ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタステップと、
上記受信信号、又は上記第１のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記第１のフィルタステップへ出力する切替ステップと、
上記フィルタ係数計算ステップにより出力されるフィルタ係数を記憶しておくフィルタ係数記憶ステップと、
上記受信信号及び上記第１のフィルタステップの出力信号から、上記第１のフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項１６】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する多次フィルタステップと、
上記クラッタ中心周波数推定ステップにより推定された推定値のうち、上記多次フィルタ係数計算ステップにおいて使用された以外の推定値に対して、複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算する第一のフィルタ係数計算ステップと、
上記クラッタ中心周波数推定ステップにより推定された上記クラッタの中心周波数の推定値のうち、上記多次フィルタ係数計算ステップにおいて使用された以外の推定値に対して、オフセット周波数を持たせて上記フィルタ群の零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算ステップと、
該第二〜第Ｎのフィルタ係数計算ステップにおいて使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎのフィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタステップと、
上記受信信号及び上記第一〜第Ｎのフィルタステップの出力信号から、上記第一〜第Ｎのフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
上記フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項１７】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタステップと、
上記クラッタ中心周波数推定ステップにより推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算ステップにおいて使用された以外の推定値に対して、オフセット周波数を持たせて縦続接続されたフィルタ群の零点を割り当てるようにフィルタ係数を計算するフィルタ係数計算ステップと、
該フィルタ係数計算ステップにおいて使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップにより出力されるフィルタ係数を用いて切替ステップにより出力される信号中のクラッタを抑圧する第二のフィルタステップと、
上記第一のフィルタステップの出力信号、又は上記第二のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記第二のフィルタステップへ出力する切替ステップと、
該切替ステップにおいて使用される上記第二のフィルタステップの出力信号を記憶しておく出力信号記憶ステップと、
上記第一及び第二のフィルタステップによる出力信号から、フィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項１８】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタステップと、
上記クラッタ中心周波数推定ステップにより推定された推定値のうち、上記多次フィルタ係数計算ステップにおいて使用された推定値以外の推定値に対して、オフセット周波数を持たせて縦続接続されたフィルタ群の零点を割り当てるようにフィルタ係数を計算するフィルタ係数計算ステップと、
該フィルタ係数計算ステップにおいて使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップにより出力されるフィルタ係数を用いて切替ステップにより出力される信号中のクラッタを抑圧する第二のフィルタステップと、
上記第一のフィルタステップの出力信号、又は上記第二のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記第二のフィルタステップへ出力する切替ステップと、
上記フィルタ係数計算ステップにより出力されるフィルタ係数を記憶しておくフィルタ係数記憶ステップと、
上記第一及び第二のフィルタステップの出力信号から、フィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行うフィルタステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項１９】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された推定値に対して、上記フィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
上記縦続接続されたフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された上記クラッタの中心周波数の推定値が、上記フィルタ係数計算手段に送られたクラッタ中心周波数推定値と同じであれば、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する第一〜第Ｎ-1の多重零点処理手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段及び上記第一〜第Ｎ-1の多重零点処理手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタ手段と、
上記受信信号及び上記第一〜第Ｎのフィルタ手段の出力信号から、上記第一〜第Ｎのフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項２０】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、フィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
２回以上同じ周波数に対して零点を形成する場合に、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段及び上記多重零点処理手段から送られるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタ手段と、
上記受信信号、又は上記フィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記フィルタ手段へ出力する切替手段と、
該切替手段へ送られる上記フィルタ手段の出力信号を記憶しておく出力信号記憶手段と、
上記受信信号及び上記フィルタ手段の出力信号から、上記フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項２１】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、フィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
２回以上同じ周波数に対して零点を形成する場合に，多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段及び上記多重零点処理手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタ手段と、
上記受信信号、又は上記第１のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第１のフィルタ手段へ出力する切替手段と、
上記フィルタ係数計算手段及び上記多重零点処理手段から出力されるフィルタ係数を記憶しておくフィルタ係数記憶手段と、
上記受信信号及び上記第１のフィルタ手段の出力信号から、上記第１のフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタ手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項２２】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する多次フィルタ手段と、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算手段へ送られた以外の推定値に対して、上記フィルタ群全体の周波数特性の多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する第一〜第Ｎ-1の多重零点処理手段と、
受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記第一〜第Ｎ-1の多重零点処理手段へ送られる零点周波数を選択する多重零点割当選択手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎ−１の多重零点処理手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎ−１のフィルタ手段と、
上記受信信号及び上記多次フィルタ手段及び上記第一〜第Ｎ−１のフィルタ手段の出力信号から、上記多次フィルタ手段及び上記第一〜第Ｎ−１のフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項２３】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一のフィルタ手段と、
上記クラッタ中心周波数推定手段により推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算手段へ送られた以外の推定値に対して、２回以上同じ周波数に対して零点を形成する場合に，多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多重零点処理手段から出力されるフィルタ係数を用いて切替手段から出力される信号中のクラッタを抑圧する第二のフィルタ手段と、
上記第一のフィルタ手段の出力信号、又は上記第二のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第二のフィルタ手段へ出力する切替手段と、
該切替手段へ送られる上記第二のフィルタ手段の出力信号を記憶しておく出力信号記憶手段と、
上記第一及び第二のフィルタ手段の出力信号から、該第二のフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項２４】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一のフィルタ手段と、
上記クラッタ中心周波数推定手段により推定された推定値のうち、上記多次フィルタ係数計算手段へ送られた推定値以外の推定値に対して、２回以上同じ周波数に対して零点を形成する場合に、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多重零点処理手段から出力されるフィルタ係数を用いて切替手段から出力される信号中のクラッタを抑圧する第二のフィルタ手段と、
上記第一のフィルタ手段の出力信号、又は上記第二のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第二のフィルタ手段へ出力する切替手段と、
上記多重零点処理手段から出力されるフィルタ係数を記憶しておくフィルタ係数記憶手段と、
上記第一及び第二のフィルタ手段の出力信号から、該第二のフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第三のフィルタ手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項２５】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
上記クラッタ中心周波数推定ステップにより推定された推定値に対して、複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
上記クラッタ中心周波数推定ステップにより推定された上記クラッタの中心周波数の推定値が、上記フィルタ係数計算ステップにおいて使用されたクラッタ中心周波数推定値と同じであれば、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する第一〜第Ｎ-1の多重零点処理ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップ及び上記第一〜第Ｎ-1の多重零点処理ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタステップと、
上記受信信号及び上記第一〜第Ｎのフィルタステップの出力信号から、上記第一〜第Ｎのフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項２６】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、フィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
２回以上同じ周波数に対して零点を形成する場合に、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップ及び上記多重零点処理ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタステップと、
上記受信信号、又は上記フィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記フィルタステップへ出力する切替ステップと、
該切替ステップにおいて使用される上記フィルタステップの出力信号を記憶しておく出力信号記憶ステップと、
上記受信信号及び上記フィルタステップの出力信号から、上記フィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項２７】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、フィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
２回以上同じ周波数に対して零点を形成する場合に，多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理ステップと
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップ及び上記多重零点処理ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタステップと、
上記受信信号、又は上記第１のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記第１のフィルタステップへ出力する切替ステップと、
上記フィルタ係数計算ステップ及び上記多重零点処理ステップにより出力されるフィルタ係数を記憶しておくフィルタ係数記憶ステップと、
上記受信信号及び上記第１のフィルタステップの出力信号から、上記第１のフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項２８】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する多次フィルタステップと、
上記クラッタ中心周波数推定ステップにより推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算ステップへ送られた以外の推定値に対して、複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群全体の周波数特性の多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する第一〜第Ｎ-1の多重零点処理ステップと、
受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記第一〜第Ｎ-1の多重零点処理ステップへ送られる零点周波数を選択する多重零点割当選択ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎ−１の多重零点処理ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎ−１のフィルタステップと、
上記受信信号及び上記多次フィルタステップ及び上記第一〜第Ｎ−１のフィルタステップの出力信号から、上記多次フィルタステップ及び上記第一〜第Ｎ−１のフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項２９】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一のフィルタステップと、
上記クラッタ中心周波数推定ステップにより推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算ステップにおいて使用された以外の推定値に対して、２回以上同じ周波数に対して零点を形成する場合に，多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多重零点処理ステップにより出力されるフィルタ係数を用いて切替ステップにより出力される信号中のクラッタを抑圧する第二のフィルタステップと、
上記第一のフィルタステップの出力信号、又は上記第二のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記第二のフィルタステップへ出力する切替ステップと、
上記切替ステップで利用される上記第二のフィルタステップの出力信号を記憶しておく出力信号記憶ステップと、
上記第一及び第二のフィルタステップの出力信号から、該第二のフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項３０】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一のフィルタステップと、
上記クラッタ中心周波数推定ステップにより推定された推定値のうち、上記多次フィルタ係数計算ステップにおいて使用された推定値以外の推定値に対して、２回以上同じ周波数に対して零点を形成する場合に、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多重零点処理ステップにより出力されるフィルタ係数を用いて切替ステップにより出力される信号中のクラッタを抑圧する第二のフィルタステップと、
上記第一のフィルタステップの出力信号、又は上記第二のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記第二のフィルタステップへ出力する切替ステップと、
上記多重零点処理ステップにより出力されるフィルタ係数を記憶しておくフィルタ係数記憶ステップと、
上記第一及び第二のフィルタステップの出力信号から、該第二のフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第三のフィルタステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項３１】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された推定値に対して、多重零点を形成できるようにを形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
上記縦続接続されたフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された上記クラッタの中心周波数の推定値に対して、フィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する第一〜第Ｎ-1の多重零点処理手段と、
受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記第一〜第Ｎ-1の多重零点処理手段へ送られる零点周波数を選択する多重零点割当選択手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段及び上記第一〜第Ｎ-1の多重零点処理手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタ手段と、
上記受信信号及び上記第一〜第Ｎのフィルタ手段の出力信号から、上記第一〜第Ｎのフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項３２】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理手段と、
受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記多重零点処理手段へ送られる零点周波数を選択する多重零点割当選択手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段及び上記多重零点処理手段から送られるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタ手段と、
上記受信信号、又は上記フィルタ手段の出力信号を、上記多重零点処理手段及び上記フィルタ手段へ出力する切替手段と、
該切替手段へ送られる上記フィルタ手段の出力信号を記憶しておく出力信号記憶手段と、
上記受信信号及び上記フィルタ手段の出力信号から、上記フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項３３】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理手段と、
受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記多重零点処理手段へ送られる零点周波数を選択する多重零点割当選択手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段及び上記多重零点処理手段から送られるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタ手段と、
上記受信信号、又は上記第１のフィルタ手段の出力信号を、上記多重零点処理手段及び上記第１のフィルタ手段へ出力する切替手段と、
上記多重零点処理手段から出力されるフィルタ係数を記憶しておくフィルタ係数記憶手段と、
上記受信信号及び上記フィルタ手段の出力信号から、上記フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタ手段と、を備えたことを特徴とするクラッタ抑圧装置。
【請求項３４】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定ステップにより推定された推定値に対して、多重零点を形成できるようにを形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
上記縦続接続されたフィルタ群の一部を構成し、上記クラッタ中心周波数推定ステップにより推定された上記クラッタの中心周波数の推定値に対して、フィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する第一〜第Ｎ-1の多重零点処理ステップと、
受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記第一〜第Ｎ-1の多重零点処理ステップへ送られる零点周波数を選択する多重零点割当選択ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップ及び上記第一〜第Ｎ-1の多重零点処理ステップから出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタステップと、
上記受信信号及び上記第一〜第Ｎのフィルタステップの出力信号から、上記第一〜第Ｎのフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップから出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項３５】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理ステップと、
受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記多重零点処理ステップへ送られる零点周波数を選択する多重零点割当選択ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップ及び上記多重零点処理ステップから送られるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタステップと、
上記受信信号、又は上記フィルタステップの出力信号を、上記多重零点処理ステップ及び上記フィルタステップへ出力する切替ステップと、
該切替ステップへ送られる上記フィルタステップの出力信号を記憶しておく出力信号記憶ステップと、
上記受信信号及び上記フィルタステップの出力信号から、上記フィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップから出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたことを特徴とするクラッタ抑圧方法。
【請求項３６】
スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理ステップと、
受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記多重零点処理ステップへ送られる零点周波数を選択する多重零点割当選択ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップ及び上記多重零点処理ステップから送られるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタステップと、
上記受信信号、又は上記第１のフィルタステップの出力信号を、上記多重零点処理ステップ及び上記第１のフィルタステップへ出力する切替ステップと、
上記多重零点処理ステップから出力されるフィルタ係数を記憶しておくフィルタ係数記憶ステップと、
上記受信信号及び上記フィルタ手段の出力信号から、上記フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップから出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタステップと、を備えたことを特徴とするクラッタ抑圧方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この発明は、例えばレーダにより受信された信号のうち目標信号以外のクラッタ等の不要反射エコーをフィルタを利用して除去するクラッタ抑圧装置およびクラッタ抑圧方法に関するものである。
【０００２】
【従来の技術】
レーダで移動目標（航空機など）を見ようとする場合、アンテナの形状や目標の位置によっては、どうしても同一ビームの中で、同一距離のところに固定目標が含まれることが多い。このような場合、固定目標からの反射波は、移動目標からの反射波に比べ、はるかに大きいことがある。従って、反射波の強度だけから移動目標と固定目標を区別することは困難である。そこで、パルスレーダの機能を損なうことなく、移動目標の反射波のみをドプラ効果を利用して取り出すことが考えられる。このような方式を移動目標識別（Ｍｏｖｉｎｇ  Ｔａｒｇｅｔ  Ｉｎｄｉｃａｔｉｏｎ；ＭＴＩ）処理という。すなわち、レーダ送信機からの送信パルスが移動目標から反射されると、その反射波の周波数はドプラ効果によって変化するが、固定目標からの反射波は周波数変位を受けない。従って、ドプラ効果による周波数変位あるいはこれに対応する位相変位だけを検出するかあるいはそれに伴う振幅変化を検出すれば移動目標だけを検出することができる。これを実現するため、受信信号を中間周波数信号に変換した後、この中間周波数信号を基準信号発生器からの基準信号で位相検波してビデオ信号を生成する。さらに、このビデオ信号を１周期遅らせた信号と次の受信信号から生成したビデオ信号との差をとる。これによりクラッタなどの固定目標をキャンセルすることができる。
【０００３】
このようにＭＴＩ処理ではビデオ信号を１周期前のビデオ信号との差をとることにより固定目標と移動目標が混在している状態から移動目標のみを抽出することができる。しかしながら、差をとることにより、得られる信号の振幅は正弦成分を含んでいるため、ドプラ周波数がパルス繰り返し周波数の整数倍となる場合には、速度レスポンスが零となるいわゆるブラインド速度を生じ、その速度の目標の検出が不可能となるという問題がある。そこで、ＭＴＩ処理において、ブラインド速度領域を狭くするために周期的な不等間隔パルスを送信するスタガトリガ方式が用いられる。これは、２種以上のパルス繰り返し周期を連続的に切り替えることにより、常にはブラインド速度の条件が成立しないようにするものである。
【０００４】
従来、スタガトリガ方式を用いたパルスレーダにより受信されたクラッタを除去するためのクラッタ抑圧装置として、電子情報通信学会  信学技報ＳＡＮＥ９９−１１５  ｐｐ．５７−６２（２０００年２月）の関口らによる「スタガトリガ方式において２つの周波数に零点を持つ単一フィルタによる複峰性クラッタの抑圧」に示されたものが知られている。
なお、一般に、レーダの受信信号は、パルス繰り返し間隔の受信信号を当該間隔よりかなり短い所定の時間で標本化したときの１標本化時間に対応する電波伝播距離を示すレンジビンと、送信パルスに対して受信されるパルスを示すパルスヒットの２次元データとして扱われるが、ここでは特に断りがない限り、任意の一レンジビン、或いは同一の処理を行う複数レンジビン単位の処理について記述することにする。
【０００５】
図２８は、従来のクラッタ抑圧装置を示すブロック図である。図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、６２はフィルタ係数を算出するフィルタ係数計算手段、６３はフィルタ係数計算手段６２からのフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段である。
【０００６】
次に、図２８に従って、動作について説明する。例えば捜索レーダのように離散フーリエ変換等による周波数解析が行えないほどパルスヒット数が少ないレーダにおいて、受信するクラッタ数を推定する方法としては、以下で説明するクラッタ中心周波数推定手段６で計算されるＡＲ（Auto Regressive）モデルの係数aijを利用する方法等が考えられるが、ここではクラッタ数は既知であるものとする。
まず、クラッタ中心周波数推定手段６において、受信信号からクラッタの中心周波数を推定する。クラッタ中心周波数は、パルスヒット数を十分確保することができれば、離散フーリエ変換等を利用して周波数解析を行い、ピーク検出処理をすることで推定できる。しかしながら、スタガトリガ方式を用いる必要があるようなパルス繰り返し周波数が高くないレーダでは、離散フーリエ変換によって十分な周波数解析を実行できるほどのパルスヒット数を確保することができない。そこで、従来のクラッタ抑圧装置では、少ないデータ数でもスペクトル推定を行う際に有効な最大エントロピー法を使用して、クラッタ中心周波数の推定を行っている。この最大エントロピー法については、S.Haykin,"Nonlinear Methods of Spectral Analysis ",Springer-Verlag(1983)に説明されている。
【０００７】
クラッタスペクトルピークが１つの単峰性クラッタの場合、１次のＡＲモデルで近似でき、このＡＲモデルの係数をａ１１とすれば、１次の極を計算で求めることにより、式（１）でクラッタ中心周波数（パルス繰り返し周波数で規格化された値）f0を推定することができる。
【０００８】
【数１】
0003932957.tif 000002
【０００９】
クラッタスペクトルピークが２つの双峰性クラッタの場合は、２次のＡＲモデルで近似でき、このＡＲモデルの係数をａｉｉ（ｉ＝１，２）、０次前向き予測誤差をe0(n)、０次後向き予測誤差をb0(n)、i次前向き予測誤差をei(n)、i次後向き予測誤差をbi(n)とすれば、２次の極を計算で求めることにより、式（３）でクラッタ中心周波数f01、f02を推定することができる。
【００１０】
【数２】
0003932957.tif 000003
【００１１】
次に、以上の式によって推定されたクラッタ中心周波数は、フィルタ係数計算手段６２に送られる。フィルタ係数計算手段６２では、スタガトリガ方式の場合においてもフィルタの深いノッチが形成できるように、フィルタの係数が時変になる。すなわち、パルス間隔ごとにフィルタの係数が変わることになる。なお、フィルタはＦＩＲ（Ｆｉｎｉｔｅ  Ｉｍｐｕｌｓｅ  Ｒｅｓｐｏｎｓｅ）形とするので、インパルス応答がフィルタ係数に対応する。
クラッタが単峰性の場合は双峰性の場合に含まれるので、ここではクラッタがスペクトルピークを２つ有する双峰性であると仮定する。スタガトリガ方式の等価振幅２乗特性を次式のように定義する。
【００１２】
【数３】
0003932957.tif 000004
【００１３】
ここで、Nはフィルタのインパルス応答長、hlnはフィルタのインパルス応答、Ｌはスタガ数、τi（i=1,2,…,L）はパルス間隔である。
次式のように時間δTを定義し、遅延演算子ｚ-1をδTだけの遅延を与える素子とする。このとき、δTはスタガ比１を示すため、スタガ比Ｒの遅延に対する遅延演算子はｚ-Ｒと表される。従って、伝達関数Cl(z)（ｌ＝０，１，…，Ｌ−１）は式（１２）のように表される。
【００１４】
【数４】
0003932957.tif 000005
【００１５】
推定した２つのクラッタ中心周波数f01、f02に、それぞれk1重零点、k2重零点を持つフィルタの係数を生成するものとする。lを固定して考えて、Cl(z)、及びそのｚに関する１階から（k1-1）階導関数に対して、z1=exp[j2πf01δT]のときに０となるようにする。
【００１６】
【数５】
0003932957.tif 000006
【００１７】
かつ、Cl(z)、及びそのｚに関する１階から（k2-1）階導関数に対して、z2=exp[j2πf02δT]のときに０となるようにする。
【００１８】
【数６】
0003932957.tif 000007
【００１９】
以上より、k1+ k2=(N-1)個の方程式ができる。係数に対して方程式が１つ少ないが、hl0を任意に決めれば良い。このようにしてできた(N-1)元1次連立方程式を解けば、特定のlに対するフィルタ係数が得られる。これをl =0,1,…,L-1に対して繰り返して実行すれば、全てのフィルタ係数を求めることができる。
【００２０】
フィルタ係数計算手段６２で生成されたフィルタ係数hlnは、フィルタ手段６３へ送られる。フィルタ手段６３では、次式に示す畳込み演算が行われて、出力信号y(tk)では、受信信号x(tk)中のクラッタが抑圧される。
【００２１】
【数７】
0003932957.tif 000008
【００２２】
以上の説明では双峰性クラッタの場合を示したが、単峰性クラッタの場合は、Cl(z)の導関数に関する条件を、例えば式（１３）についてのみ考慮すればよい。
【００２３】
【発明が解決しようとする課題】
従来のスタガトリガ方式によるレーダにおけるクラッタ抑圧装置は、推定したクラッタ中心周波数に対してフィルタのノッチを形成するようにフィルタ係数を調整するが、クラッタを抑圧するための時変ノッチフィルタの次数を予め決めておく必要があるため、受信するクラッタの状況によってはフィルタ次数の過不足が発生して、クラッタの消え残りが生じたり、フィルタのノッチの阻止域がクラッタに対して過剰になって目標信号の減衰を誘発する、という問題点があった。
【００２４】
この発明は上記のような問題点を解決するためになされたもので、クラッタ抑圧フィルタを複数のフィルタの縦続接続構成にして、各フィルタの入出力信号の電力比を用いてクラッタ抑圧に必要なフィルタ次数を推定することで、自動的に最適なフィルタ次数で抑圧処理を実施することが可能なクラッタ抑圧装置を得ることを目的とする。
【００２５】
【課題を解決するための手段】
第１の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧装置において、  上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された推定値に対して、上記フィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算する第一のフィルタ係数計算手段と、
上記縦続接続されたフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された上記クラッタの中心周波数の推定値に対して、オフセット周波数を持たせて自フィルタの零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算手段と、
該第二〜第Ｎのフィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎのフィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタ手段と、
上記受信信号及び上記第一〜第Ｎのフィルタ手段の出力信号から、上記第一〜第Ｎのフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたものである。
【００２６】
また、第２の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
該フィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段から送られるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタ手段と、
上記受信信号、又は上記フィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記フィルタ手段へ出力する切替手段と、
該切替手段へ送られる上記フィルタ手段の出力信号を記憶しておく出力信号記憶手段と、
上記受信信号及び上記フィルタ手段の出力信号から、上記フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたものである。
【００２７】
また、第３の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
該フィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタ手段と、
上記受信信号、又は上記第１のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第１のフィルタ手段へ出力する切替手段と、
上記フィルタ係数計算手段から出力されるフィルタ係数を記憶しておくフィルタ係数記憶手段と、
上記受信信号及び上記第１のフィルタ手段の出力信号から、上記第１のフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタ手段と、を備えたものである。
【００２８】
また、第４の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された推定値に対して、上記フィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算する第一のフィルタ係数計算手段と、
上記縦続接続されたフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された上記クラッタの中心周波数の推定値に対して、オフセット周波数を持たせて自フィルタの零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算手段と、
該第二〜第Ｎのフィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記受信信号中の複数のクラッタに対してフィルタのノッチが形成され、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記オフセット周波数に記憶されているオフセット周波数を設定するオフセット周波数選択手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎのフィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタ手段と、
上記受信信号及び上記第一〜第Ｎのフィルタ手段の出力信号から、上記第一〜第Ｎのフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたものである。
【００２９】
また、第５の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
上記受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記フィルタ係数計算手段へ送られるオフセット周波数を選択するオフセット周波数設定手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタ手段と、  上記受信信号、又は上記フィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記フィルタ手段へ出力する切替手段と、
該切替手段に送られる上記フィルタ手段の出力信号を記憶しておく出力信号記憶手段と、
上記受信信号及び上記フィルタ手段の出力信号から、上記フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたものである。
【００３０】
また、第６の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
上記受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記フィルタ係数計算手段へ送られるオフセット周波数を選択するオフセット周波数設定手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタ手段と、
上記受信信号、又は上記第１のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第１のフィルタ手段へ出力する切替手段と、
上記フィルタ係数計算手段から出力されるフィルタ係数を記憶しておくフィルタ係数記憶手段と、
上記受信信号及び上記第１のフィルタ手段の出力信号から、上記第１のフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタ手段と、を備えたものである。
【００３１】
また、第７の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する多次フィルタ手段と、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算手段へ送られた以外の推定値に対して、上記フィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算する第一のフィルタ係数計算手段と、
上記縦続接続されたフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定されたクラッタ中心周波数の推定値うち、上記多次フィルタ係数計算手段へ送られた以外の推定値に対して、オフセット周波数を持たせて自フィルタの零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算手段と、
該第二〜第Ｎのフィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎのフィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタ手段と、
上記受信信号及び上記第一〜第Ｎのフィルタ手段の出力信号から、上記第一〜第Ｎのフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
上記フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたものである。
【００３２】
また、第８の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一のフィルタ手段と、
上記クラッタ中心周波数推定手段により推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算手段へ送られた以外の推定値に対して、オフセット周波数を持たせて自フィルタの零点を割り当てるようにフィルタ係数を計算するフィルタ係数計算手段と、
該フィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段から出力されるフィルタ係数を用いて切替手段から出力される信号中のクラッタを抑圧する第二のフィルタ手段と、
上記第一のフィルタ手段の出力信号、又は上記第二のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第二のフィルタ手段へ出力する切替手段と、
該切替手段へ送られる上記第二のフィルタ手段の出力信号を記憶しておく出力信号記憶手段と、
上記第一及び第二のフィルタ手段の出力信号から、フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、  該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたものである。
【００３３】
また、第９の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一のフィルタ手段と、
上記クラッタ中心周波数推定手段により推定された推定値のうち、上記多次フィルタ係数計算手段へ送られた推定値以外の推定値に対して、オフセット周波数を持たせて自フィルタの零点を割り当てるようにフィルタ係数を計算するフィルタ係数計算手段と、
該フィルタ係数計算手段へ送られるオフセット周波数を記憶するオフセット周波数記憶手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段から出力されるフィルタ係数を用いて切替手段から出力される信号中のクラッタを抑圧する第二のフィルタ手段と、
上記第一のフィルタ手段の出力信号、又は上記第二のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第二のフィルタ手段へ出力する切替手段と、
上記フィルタ係数計算手段から出力されるフィルタ係数を記憶しておくフィルタ係数記憶手段と、
上記第一及び第二のフィルタ手段の出力信号から、フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、  該フィルタ次数選定手段から出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行うフィルタ手段と、を備えたものである。
【００３４】
また、第１０の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
上記クラッタ中心周波数推定ステップにより推定された推定値に対して、複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算する第一のフィルタ係数計算ステップと、
上記クラッタ中心周波数推定ステップにより推定された上記クラッタの中心周波数の推定値に対して、オフセット周波数を持たせて上記フィルタ群の零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算ステップと、
該第二〜第Ｎのフィルタ係数計算ステップにおいて使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎのフィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタステップと、
上記受信信号及び上記第一〜第Ｎのフィルタステップの出力信号から、上記第一〜第Ｎのフィルタステップにおいて入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたものである。
【００３５】
また、第１１の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
該フィルタ係数計算ステップで使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップにより出力されるフィルタ係数を用いて受信信号中のクラッタを抑圧するフィルタステップと、
上記受信信号、又は上記フィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記フィルタステップへ出力する切替ステップと、
上記切替ステップにおいて使用される上記フィルタステップの出力信号を記憶しておく出力信号記憶ステップと、
上記受信信号及び上記フィルタステップの出力信号から、上記フィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたものである。
【００３６】
また、第１２の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
該フィルタ係数計算ステップにおいて使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタステップと、
上記受信信号、又は上記第１のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記フィルタステップへ出力する切替ステップと、  上記フィルタ係数計算ステップにより出力されるフィルタ係数を記憶しておくフィルタ係数記憶ステップと、
上記受信信号及び上記第１のフィルタステップの出力信号から、上記第１のフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタステップと、を備えたものである。
【００３７】
また、第１３の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
上記クラッタ中心周波数推定ステップにより推定された推定値に対して、複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算する第一のフィルタ係数計算ステップと、
上記クラッタ中心周波数推定ステップにより推定された上記クラッタの中心周波数の推定値に対して、オフセット周波数を持たせて上記フィルタ群の零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算ステップと、
該第二〜第Ｎのフィルタ係数計算ステップにおいて使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
上記受信信号中の複数のクラッタに対してフィルタのノッチが形成され、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記オフセット周波数に記憶されているオフセット周波数を設定するオフセット周波数選択ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎのフィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタステップと、
上記受信信号及び上記第一〜第Ｎのフィルタステップの出力信号から、上記第一〜第Ｎのフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたものである。
【００３８】
また、第１４の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
上記受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記フィルタ係数計算ステップにおいて使用されるオフセット周波数を選択するオフセット周波数設定ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタステップと、
上記受信信号、又は上記フィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記フィルタステップへ出力する切替ステップと、
該切替ステップにおいて使用される上記フィルタステップの出力信号を記憶しておく出力信号記憶ステップと、
上記受信信号及び上記フィルタステップによる出力信号から、上記フィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたものである。
【００３９】
また、第１５の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、オフセット周波数を加えてフィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
上記受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記フィルタ係数計算ステップにおいて使用されるオフセット周波数を選択するオフセット周波数設定ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタステップと、
上記受信信号、又は上記第１のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記第１のフィルタステップへ出力する切替ステップと、
上記フィルタ係数計算ステップにより出力されるフィルタ係数を記憶しておくフィルタ係数記憶ステップと、
上記受信信号及び上記第１のフィルタステップの出力信号から、上記第１のフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタステップと、を備えたものである。
【００４０】
また、第１６の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する多次フィルタステップと、
上記クラッタ中心周波数推定ステップにより推定された推定値のうち、上記多次フィルタ係数計算ステップにおいて使用された以外の推定値に対して、複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算する第一のフィルタ係数計算ステップと、
上記クラッタ中心周波数推定ステップにより推定された上記クラッタの中心周波数の推定値のうち、上記多次フィルタ係数計算ステップにおいて使用された以外の推定値に対して、オフセット周波数を持たせて上記フィルタ群の零点を割り当てるようにフィルタ係数を計算する第二〜第Ｎのフィルタ係数計算ステップと、
該第二〜第Ｎのフィルタ係数計算ステップにおいて使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎのフィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタステップと、
上記受信信号及び上記第一〜第Ｎのフィルタステップの出力信号から、上記第一〜第Ｎのフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
上記フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたものである。
【００４１】
また、第１７の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタステップと、
上記クラッタ中心周波数推定ステップにより推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算ステップにおいて使用された以外の推定値に対して、オフセット周波数を持たせて縦続接続されたフィルタ群の零点を割り当てるようにフィルタ係数を計算するフィルタ係数計算ステップと、
該フィルタ係数計算ステップにおいて使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップにより出力されるフィルタ係数を用いて切替ステップにより出力される信号中のクラッタを抑圧する第二のフィルタステップと、
上記第一のフィルタステップの出力信号、又は上記第二のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記第二のフィルタステップへ出力する切替ステップと、
該切替ステップにおいて使用される上記第二のフィルタステップの出力信号を記憶しておく出力信号記憶ステップと、
上記第一及び第二のフィルタステップによる出力信号から、フィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたものである。
【００４２】
また、第１８の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタステップと、
上記クラッタ中心周波数推定ステップにより推定された推定値のうち、上記多次フィルタ係数計算ステップにおいて使用された推定値以外の推定値に対して、オフセット周波数を持たせて縦続接続されたフィルタ群の零点を割り当てるようにフィルタ係数を計算するフィルタ係数計算ステップと、
該フィルタ係数計算ステップにおいて使用されるオフセット周波数を記憶するオフセット周波数記憶ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップにより出力されるフィルタ係数を用いて切替ステップにより出力される信号中のクラッタを抑圧する第二のフィルタステップと、
上記第一のフィルタステップの出力信号、又は上記第二のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記第二のフィルタステップへ出力する切替ステップと、
上記フィルタ係数計算ステップにより出力されるフィルタ係数を記憶しておくフィルタ係数記憶ステップと、
上記第一及び第二のフィルタステップの出力信号から、フィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行うフィルタステップと、を備えたものである。
【００４３】
また、第１９の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された推定値に対して、上記フィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
上記縦続接続されたフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された上記クラッタの中心周波数の推定値が、上記フィルタ係数計算手段に送られたクラッタ中心周波数推定値と同じであれば、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する第一〜第Ｎ-1の多重零点処理手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段及び上記第一〜第Ｎ-1の多重零点処理手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタ手段と、
上記受信信号及び上記第一〜第Ｎのフィルタ手段の出力信号から、上記第一〜第Ｎのフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたものである。
【００４４】
また、第２０の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、フィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
２回以上同じ周波数に対して零点を形成する場合に、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段及び上記多重零点処理手段から送られるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタ手段と、
上記受信信号、又は上記フィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記フィルタ手段へ出力する切替手段と、
該切替手段へ送られる上記フィルタ手段の出力信号を記憶しておく出力信号記憶手段と、
上記受信信号及び上記フィルタ手段の出力信号から、上記フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたものである。
【００４５】
また、第２１の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、フィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
２回以上同じ周波数に対して零点を形成する場合に，多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理手段と
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段及び上記多重零点処理手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタ手段と、
上記受信信号、又は上記第１のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第１のフィルタ手段へ出力する切替手段と、
上記フィルタ係数計算手段及び上記多重零点処理手段から出力されるフィルタ係数を記憶しておくフィルタ係数記憶手段と、
上記受信信号及び上記第１のフィルタ手段の出力信号から、上記第１のフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタ手段と、を備えたものである。
【００４６】
また、第２２の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する多次フィルタ手段と、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算手段へ送られた以外の推定値に対して、上記フィルタ群全体の周波数特性の多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する第一〜第Ｎ-1の多重零点処理手段と、
受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記第一〜第Ｎ-1の多重零点処理手段へ送られる零点周波数を選択する多重零点割当選択手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎ−１の多重零点処理手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎ−１のフィルタ手段と、
上記受信信号及び上記多次フィルタ手段及び上記第一〜第Ｎ−１のフィルタ手段の出力信号から、上記多次フィルタ手段及び上記第一〜第Ｎ−１のフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたものである。
【００４７】
また、第２３の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一のフィルタ手段と、
上記クラッタ中心周波数推定手段により推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算手段へ送られた以外の推定値に対して、２回以上同じ周波数に対して零点を形成する場合に，多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多重零点処理手段から出力されるフィルタ係数を用いて切替手段から出力される信号中のクラッタを抑圧する第二のフィルタ手段と、
上記第一のフィルタ手段の出力信号、又は上記第二のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第二のフィルタ手段へ出力する切替手段と、
該切替手段へ送られる上記第二のフィルタ手段の出力信号を記憶しておく出力信号記憶手段と、
上記第一及び第二のフィルタ手段の出力信号から、該第二のフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたものである。
【００４８】
また、第２４の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一のフィルタ手段と、
上記クラッタ中心周波数推定手段により推定された推定値のうち、上記多次フィルタ係数計算手段へ送られた推定値以外の推定値に対して、２回以上同じ周波数に対して零点を形成する場合に、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多重零点処理手段から出力されるフィルタ係数を用いて切替手段から出力される信号中のクラッタを抑圧する第二のフィルタ手段と、
上記第一のフィルタ手段の出力信号、又は上記第二のフィルタ手段の出力信号を、上記フィルタ係数計算手段及び上記第二のフィルタ手段へ出力する切替手段と、
上記多重零点処理手段から出力されるフィルタ係数を記憶しておくフィルタ係数記憶手段と、
上記第一及び第二のフィルタ手段の出力信号から、該第二のフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、  該フィルタ次数選定手段から出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第三のフィルタ手段と、を備えたものである。
【００４９】
また、第２５の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
上記クラッタ中心周波数推定ステップにより推定された推定値に対して、複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群全体の周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
上記クラッタ中心周波数推定ステップにより推定された上記クラッタの中心周波数の推定値が、上記フィルタ係数計算ステップにおいて使用されたクラッタ中心周波数推定値と同じであれば、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する第一〜第Ｎ-1の多重零点処理ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップ及び上記第一〜第Ｎ-1の多重零点処理ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタステップと、
上記受信信号及び上記第一〜第Ｎのフィルタステップの出力信号から、上記第一〜第Ｎのフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたものである。
【００５０】
また、第２６の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、フィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
２回以上同じ周波数に対して零点を形成する場合に、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップ及び上記多重零点処理ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタステップと、
上記受信信号、又は上記フィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記フィルタステップへ出力する切替ステップと、
該切替ステップにおいて使用される上記フィルタステップの出力信号を記憶しておく出力信号記憶ステップと、
上記受信信号及び上記フィルタステップの出力信号から、上記フィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたものである。
【００５１】
また、第２７の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである単峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された推定値に対して、フィルタ周波数特性の零点を形成するようにフィルタ係数を計算するフィルタ係数計算ステップと、
２回以上同じ周波数に対して零点を形成する場合に，多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理ステップと
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算ステップ及び上記多重零点処理ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第１のフィルタステップと、
上記受信信号、又は上記第１のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記第１のフィルタステップへ出力する切替ステップと、
上記フィルタ係数計算ステップ及び上記多重零点処理ステップにより出力されるフィルタ係数を記憶しておくフィルタ係数記憶ステップと、
上記受信信号及び上記第１のフィルタステップの出力信号から、上記第１のフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第２のフィルタステップと、を備えたものである。
【００５２】
また、第２８の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する多次フィルタステップと、
上記クラッタ中心周波数推定ステップにより推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算ステップへ送られた以外の推定値に対して、複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群全体の周波数特性の多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する第一〜第Ｎ-1の多重零点処理ステップと、
受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記第一〜第Ｎ-1の多重零点処理ステップへ送られる零点周波数を選択する多重零点割当選択ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記第一〜第Ｎ−１の多重零点処理ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎ−１のフィルタステップと、
上記受信信号及び上記多次フィルタステップ及び上記第一〜第Ｎ−１のフィルタステップの出力信号から、上記多次フィルタステップ及び上記第一〜第Ｎ−１のフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたものである。
【００５３】
また、第２９の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一のフィルタステップと、
上記クラッタ中心周波数推定ステップにより推定されたクラッタ中心周波数の推定値のうち、上記多次フィルタ係数計算ステップにおいて使用された以外の推定値に対して、２回以上同じ周波数に対して零点を形成する場合に，多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多重零点処理ステップにより出力されるフィルタ係数を用いて切替ステップにより出力される信号中のクラッタを抑圧する第二のフィルタステップと、
上記第一のフィルタステップの出力信号、又は上記第二のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記第二のフィルタステップへ出力する切替ステップと、
上記切替ステップで利用される上記第二のフィルタステップの出力信号を記憶しておく出力信号記憶ステップと、
上記第一及び第二のフィルタステップの出力信号から、該第一及び第二のフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて出力信号を切り替える出力信号選択ステップと、を備えたものである。
【００５４】
また、第３０の発明に係るクラッタ抑圧方法は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧方法において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定ステップと、
該クラッタ中心周波数推定ステップにより推定された一つの推定値に対して、二次以上の次数の時変係数フィルタである多次フィルタの周波数特性の零点を形成するように上記多次フィルタの係数を計算する多次フィルタ係数計算ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多次フィルタ係数計算ステップにより出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一のフィルタステップと、
上記クラッタ中心周波数推定ステップにより推定された推定値のうち、上記多次フィルタ係数計算ステップにおいて使用された推定値以外の推定値に対して、２回以上同じ周波数に対して零点を形成する場合に、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理ステップと、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記多重零点処理ステップにより出力されるフィルタ係数を用いて切替ステップにより出力される信号中のクラッタを抑圧する第二のフィルタステップと、
上記第一のフィルタステップの出力信号、又は上記第二のフィルタステップの出力信号を、上記フィルタ係数計算ステップ及び上記第二のフィルタステップへ出力する切替ステップと、
上記多重零点処理ステップにより出力されるフィルタ係数を記憶しておくフィルタ係数記憶ステップと、
上記第一及び第二のフィルタステップの出力信号から、該第二のフィルタステップの入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定ステップと、
該フィルタ次数選定ステップにより出力される選定結果に基づいて、選定されたフィルタ係数を使用してクラッタ抑圧処理を行う第三のフィルタステップと、を備えたものである。
【００５５】
また、第３１の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
複数Ｎ個の時変係数フィルタを縦続接続して成るフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された推定値に対して、多重零点を形成できるようにを形成するようにフィルタ係数を計算するフィルタ係数計算手段と、
上記縦続接続されたフィルタ群の一部を構成し、上記クラッタ中心周波数推定手段により推定された上記クラッタの中心周波数の推定値に対して、フィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する第一〜第Ｎ-1の多重零点処理手段と、
受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記第一〜第Ｎ-1の多重零点処理手段へ送られる零点周波数を選択する多重零点割当選択手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段及び上記第一〜第Ｎ-1の多重零点処理手段から出力されるフィルタ係数を用いて上記受信信号中のクラッタを抑圧する第一〜第Ｎのフィルタ手段と、
上記受信信号及び上記第一〜第Ｎのフィルタ手段の出力信号から、上記第一〜第Ｎのフィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたものである。
【００５６】
また、第３２の発明に係るクラッタ抑圧装置は、スタガトリガ方式のパルスレーダにおける受信信号中に含まれる不要反射エコーである双峰性クラッタを抑圧するクラッタ抑圧装置において、
上記クラッタの中心周波数を推定するクラッタ中心周波数推定手段と、
該クラッタ中心周波数推定手段により推定された推定値に対して、多重零点を形成できるようにフィルタ伝達関数の高階導関数を用いてフィルタ係数を計算する多重零点処理手段と、
受信信号中のクラッタ数に応じて、上記クラッタ中心周波数に対してフィルタの零点を交互に設定するように、上記多重零点処理手段へ送られる零点周波数を選択する多重零点割当選択手段と、
上記スタガトリガ方式によって異なるパルス間隔ごとに上記フィルタ係数計算手段及び上記多重零点処理手段から送られるフィルタ係数を用いて上記受信信号中のクラッタを抑圧するフィルタ手段と、
上記受信信号、又は上記フィルタ手段の出力信号を、上記多重零点処理手段及び上記フィルタ手段へ出力する切替手段と、
該切替手段へ送られる上記フィルタ手段の出力信号を記憶しておく出力信号記憶手段と、
上記受信信号及び上記フィルタ手段の出力信号から、上記フィルタ手段の入出力信号電力比を計算して電力比変化に基づいてフィルタ全体の次数を決定するフィルタ次数選定手段と、
該フィルタ次数選定手段から出力される選定結果に基づいて出力信号を切り替える出力信号選択手段と、を備えたものである。
【００５７】
【発明の実施の形態】
実施の形態１．
図１はこの発明に係るクラッタ抑圧装置の実施の形態１を示す構成図である。図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、１ａはこのクラッタ中心周波数推定手段６から出力される推定値に対して、フィルタのノッチを形成するようにフィルタ係数を調整する第一のフィルタ係数計算手段、２ａはこのフィルタ係数計算手段１ａで算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第一のフィルタ手段、５はクラッタ中心周波数推定手段６の推定値に対して、クラッタ抑圧性能を確保するために若干の周波数オフセットを与えるためのオフセット周波数記憶手段、１ｂ〜１Ｎはクラッタ中心周波数推定手段６から出力される推定値に、オフセット周波数記憶手段５から出力されるオフセット周波数を加えた周波数値に対して、フィルタのノッチを形成するようにフィルタ係数を調整する第二〜第Ｎのフィルタ係数計算手段、２ｂ〜２Ｎはこのフィルタ係数計算手段１ｂ〜１Ｎで算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第二〜第Ｎのフィルタ手段、３は受信信号及び第一〜第Ｎのフィルタ手段２ａ〜２Ｎの出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、４はフィルタ次数選定手段３から出力される選定結果に基づいて、クラッタ抑圧フィルタとしての出力信号を選択する出力信号選択手段である。
図１４はこの発明に係るクラッタ抑圧装置の実施の形態１の動作を説明するためのブロック図である。
図２はこの発明に係るクラッタ抑圧装置の実施の形態１の動作を説明するためのフローチャートである。
図３はこの発明に係るクラッタ抑圧装置の実施の形態１の動作を説明するための説明図である。
【００５８】
次に図１〜３及び図１４を参照して動作について説明する。離散フーリエ変換等による周波数解析が行えないほどパルスヒット数が少ないレーダにおいて、受信するクラッタ数を推定する方法としては、以下で説明するクラッタ中心周波数推定手段１で計算されるＡＲ係数aijを利用する方法等が考えられるが、本発明の実施の形態では、クラッタ数は既知であるものとする。本実施の形態１では、クラッタは単峰性であるとする。
まず、図１において受信信号がクラッタ中心周波数推定手段６に送られて、従来のクラッタ抑圧装置と同様な方法でクラッタ中心周波数が推定される。
【００５９】
次に、クラッタ中心周波数推定手段６から出力されるクラッタ中心周波数の推定値に従ってクラッタ抑圧フィルタの係数が計算される。従来のクラッタ抑圧装置と異なる点は、最適なフィルタ次数を自動的に設定するために、従来単一のフィルタで実施していたクラッタ抑圧処理を、複数の１次フィルタの縦続接続構成にし、且つ、同じ周波数に複数のフィルタでノッチを形成する場合に２段目以降のフィルタに入力されるべきクラッタ中心周波数の推定値にオフセット周波数を加えるようにしたことである。本実施の形態では縦続接続したフィルタの次数は１次としたが、２次以上に設定してもかまわない。
【００６０】
スタガトリガ方式では不等間隔パルスになるため、ノッチ周波数が０以外のフィルタを実現するには、フィルタ係数がパルス間隔に応じて変わる時変係数フィルタが必要である。このような時変係数フィルタを縦続接続した場合のインパルス応答は、２つのフィルタのインパルス応答の単純な畳み込みにはならない。縦続接続した後段のフィルタ処理では、後段のフィルタの入力信号における現時点でのサンプルとそれ以前のサンプルとでは、それらのサンプルを算出するために用いた前段のフィルタ係数が異なるからである。以下に、このような時変係数フィルタを縦続接続した場合のインパルス応答を示す。ここでは説明の都合上、図１４に示すように縦続接続するフィルタ数は２、すなわち図１の時変フィルタ＃１と時変フィルタ＃２について取り扱い、また、次数をそれぞれK1、K2とした一般的な場合について説明する。
パルスヒット番号に対応する時刻をtkとする。第一のフィルタ手段２ａの入力信号をu(tk)、第二のフィルタ手段２ｂの入力信号をx(tk)、第二のフィルタ手段２ｂの出力信号をy(tk)とする。今、スタガ数Lを３、第一のフィルタ係数計算手段１ａで計算されるフィルタ係数をhk0、hk1、hk2、hk3、第二の係数計算手段１ｂで計算されるフィルタ係数をgk0、gk1、gk2とする。ある時刻の出力信号y(tk)は次式で表すことができる。
【００６１】
【数８】
0003932957.tif 000009
【００６２】
これより、＊を[  ]で囲んだ２つのベクトルの畳み込み演算とすると、２つの時変フィルタを縦続接続した場合のインパルス応答は、次式のようになる。なお、同式中、＋は畳み込みの結果得られた数列をベクトルとみなしたときのベクトル加算である。
【００６３】
【数９】
0003932957.tif 000010
【００６４】
第一のフィルタ係数計算手段１ａの次数をK1、第二のフィルタ係数計算手段１ｂの次数をK2として、上式を一般形に拡張する。得られる畳み込み演算の結果をelk（l=0,1,…,L-1、k=0,1,…,K1+ K2）とすると、次式のように表すことができる。
【００６５】
【数１０】
0003932957.tif 000011
【００６６】
上式において、パルス間隔はスタガ数Ｌで周期的に変化するため、時刻t0,tL,t2L,…には、g00,g01,…,g0,K2が対応し、時刻t1,tL+1,t2L+1,…には、g10,g11,…,g1,K2が対応し、時刻t2,tL+2,t2L+2,…には、g20,g21,…,g2,K2が対応する（以降繰り返し）。
このときの等価振幅２乗特性E(f)は、次式のようになる。
【００６７】
【数１１】
0003932957.tif 000012
【００６８】
ノッチ周波数が０のフィルタのノッチ周波数を任意の周波数に移動させる場合、第一のフィルタ係数計算手段１ａで使用される、ノッチ周波数０のフィルタ係数をhlk'、移動先のノッチ周波数をf01とする。また、仮に第二のフィルタ係数計算手段１ｂで使用される、ノッチ周波数０のフィルタ係数をglk'、移動先のノッチ周波数をf02とする。このとき、フィルタ単体で考えると、hlk、glkは次式のように再計算される。
【００６９】
【数１２】
0003932957.tif 000013
【００７０】
しかし、時変フィルタの縦続接続構成の場合、第二のフィルタ係数計算手段１ｂで再計算される係数glkは、上式のような単純なノッチフィルタのノッチ周波数の移動にはならない。第二のフィルタ係数計算手段１ｂで計算される係数が、第一のフィルタ係数計算手段１ａで計算される係数に依存するからである。
以下に、第二のフィルタ係数計算手段１ｂにおいて実施される、フィルタ係数計算方法について示す。式（１０）に示したように時間δTを定義し、遅延演算子z-1をδTだけの遅延を与える素子とする。このとき、時変フィルタ縦続接続時の伝達関数は、次式のように表すことができる。
【００７１】
【数１３】
0003932957.tif 000014
【００７２】
式（２４）の指数部分を式（２６）のように置き換えると、式（２５）のようになる。
【００７３】
【数１４】
0003932957.tif 000015
【００７４】
式（２５）に式（１９）を代入して整理すると、式（２７）のようになる。
【００７５】
【数１５】
0003932957.tif 000016
【００７６】
上記伝達関数El(z)において、zに関する導関数を求め、式（２８）のときに０となる、すなわち、式（２９）が成り立つようにすると、f=f02にK2重零点を持たせることができる。
【００７７】
【数１６】
0003932957.tif 000017
【００７８】
上式におけるEl(z)の導関数は、式（３０）で求めることができる。
【００７９】
【数１７】
0003932957.tif 000018
【００８０】
このようにして、フィルタ係数の数より１つ少ないK2個の方程式ができるが、gl0を任意に決めればよい。また、縦続接続された２つのフィルタにより、フィルタ入出力前後の受信機雑音電力が変化しないようにするためには、式（３１）の定数clをgl0，gl1，…，gl,K2に乗じてやればよい。
【００８１】
【数１８】
0003932957.tif 000019
【００８２】
以上の関係を行列表現するために、ベクトルを式（３２）のように定義すると、式（３３）に示す連立方程式が得られる。
【００８３】
【数１９】
0003932957.tif 000020
【００８４】
上式において、K2×K2行列のAl  、ベクトルblは次式で表される。
【００８５】
【数２０】
0003932957.tif 000021
【００８６】
なお、式（３５）は任意の値ｇｌ０を１とおいて計算したものである。式（３３）を解くことで、特定のlに関する最適な係数を求めることができる。従って、式（３３）をl=0,1,…,L-1についてL回解くことにより、f=f02にK2重零点を持つ２段目のフィルタ係数を求めることができる。
このようにして、第一のフィルタ係数計算手段１ａで計算された係数を第一のフィルタ手段２ａに適用し、第二のフィルタ係数計算手段１ｂで計算された係数を第二のフィルタ手段２ｂに適用することで、フィルタ全体の特性としてf=f01とf=f02にノッチを形成することができる。
【００８７】
以上より、２つの時変フィルタを縦続接続しても、２個所の周波数にノッチを形成できることが示された。従って、縦続接続した時変フィルタ＃１と時変フィルタ＃２を一つのフィルタとみなして、上記フィルタ係数計算方法により、このフィルタと時変フィルタ＃３を縦続接続した時のフィルタ係数を第三のフィルタ係数計算手段１ｃで計算するようにすれば、縦続接続したフィルタの数に等しい数の周波数に対してノッチを形成するようなクラッタ抑圧フィルタを構成することが可能である。
本実施の形態１では、単峰性のクラッタを想定しているため、第一のフィルタ係数計算手段１ａに送られるクラッタ中心周波数推定値と、第二〜第のフィルタ係数計算手段１ｂ〜１Ｎに送られるクラッタ中心周波数推定値は同一の周波数になる。つまり、クラッタ中心周波数の推定値をf0とすると、f01=f02=…=f0として各フィルタ係数が計算される。
しかしながら、以上に記述した方法では、複数のフィルタで零点を割り当てる周波数を重複して設定する場合、後段のフィルタで割り当てた零点の多重度のみが支配的になり、縦続接続したフィルタ次数に相当するフィルタ特性が得られないことが経験的に知られている。具体的には、フィルタのノッチの帯域幅が狭くなるため、従来のクラッタ抑圧装置に比べてクラッタ抑圧性能が劣化してしまう。
【００８８】
そこで、本実施の形態１では、図１に示すようにオフセット周波数記憶手段５を設け、第一のフィルタ係数計算手段１ａで零点が割り当てられた、単峰性クラッタの中心周波数推定値f0に対して、第二〜第Ｎのフィルタ係数計算手段１ｂ〜１Ｎで零点を割り当てる周波数がf0に一致しないようにしている。第二〜第Ｎのフィルタ係数計算手段１ｂ〜１Ｎでは、クラッタ中心周波数推定手段６から出力される推定値f0と、オフセット周波数計算手段５から出力されるオフセット周波数値δfとを次式に示すように加算している。
fm+1=f0+δfm（m=1,2,…,N-1）  （３６）
ただし、fmは縦続接続されたm番目の時変フィルタで零点を設定する周波数、δfmはm番目のオフセット周波数である。Ｎ個のフィルタを縦続接続した場合は、式（３６）に示したようにオフセット周波数δfがＮ−１種類必要になる。上式により計算された周波数fmに対して零点を割り当てるようにし、かつfm=f0にはならないため、従来のクラッタ抑圧装置におけるフィルタと同等な等価振幅２乗特性を得ることができる。
【００８９】
図１のように、１次の時変係数フィルタを多数縦続接続すると、各フィルタの入力信号中のクラッタ電力が大きく、フィルタのノッチがクラッタの中心周波数近辺に形成されていれば、各フィルタ処理後の出力信号電力は入力信号電力よりも小さくなる。クラッタがある程度抑圧されてしまえば、新たにフィルタ処理をしても出力信号電力が大きく変化することはないので、この電力変化を検出してフィルタの次数を選定する。各フィルタの入出力信号電力比はフィルタの番号２ａ〜２Ｎに向かって大きくなるに連れて次第に小さくなり、その変化も次第に小さくなる。従って、フィルタの番号を横軸にし、入出力電力比を縦軸にすると、電力比の頂点は下に凸で遠方はほぼ平らになるなめらかな曲線上に乗る。従って、判定方法は、任意に設定できるしきい値と各フィルタの入出力信号電力比を比較し、電力比がしきい値を下回った状態はフィルタ次数が過剰と判断できるので、電力比がしきい値を初めて上回った時のフィルタの出力信号をクラッタ抑圧フィルタの出力信号にすることで、必要十分なフィルタ次数でクラッタ抑圧処理を実施することができる。
以上のフィルタ次数選定処理を、図２のフローチャートに従って説明する。まず、ステップ２０でフィルタ次数のカウンタを初期化する。ステップ２１では、フィルタ次数がフィルタ次数の上限値NFIL以下であればステップ２２へ処理が移行し、そうでなければステップ２５へ処理が移行する。ステップ２２、２３では、それぞれ各フィルタの入出力信号の平均電力、電力比を計算する。このとき、レンジビン毎の電力比のばらつきをおさえるため、複数のレンジビンで平均化処理を行う。以上の処理で各フィルタの入出力電力比が計算される。ステップ２５では仮のフィルタ次数Ｍに上限値NFILを代入する。ステップ２６ではフィルタ次数が１次と判定された時点で処理が終了するように判定する。ステップ２７ではステップ２３で計算した各フィルタの入出力信号電力比としきい値を比較して、電力比がしきい値を上回った時点で処理がステップ２９に移行する。そうでなければ、ステップ２８で仮のフィルタ次数Ｍを更新して、ステップ２６へ移る。即ち、電力比がしきい値を下回っている間は仮りのフィルタ次数Ｍを1つずつ減らし、電力比が初めてしきい値を上回ったときにステップ２９へ移る。ステップ２９では仮のフィルタ次数Ｍに基づいて、最適なフィルタ次数を選定して処理が終了する。
【００９０】
図３に計算例を示す。クラッタ抑圧フィルタへの入力信号（受信信号）は、疑似的に発生させた単峰性クラッタ、目標信号と受信機雑音を加算した複素信号とした。スタガトリガ数を3、パルス繰り返し周期を400,300,500μs（平均400μs）、目標信号のドップラ周波数を500Hz、目標信号が存在するレンジビンを210レンジビン、目標信号電力を35dB、単峰性クラッタの中心周波数を-600Hz、ガウス型のスペクトルを持つクラッタの帯域幅を平均50Hz、クラッタ電力を140-40log10(100+r)dB（rはレンジビン番号）、受信機雑音電力を20dB（平均0の正規分布）、フィルタの最大数を４として、計算機シミュレーションを行った。(a)は受信信号及び各フィルタの出力信号の変化を示すグラフで、横軸がレンジビン、縦軸が平均電力を表す。フィルタ＃３以降のフィルタでは、出力信号電力がほとんど雑音電力付近であることがわかる。(b)は各フィルタの入出力信号電力比の変化を示すグラフで、横軸がレンジビン、縦軸がフィルタ入出力信号電力比を表す。ここでは２５レンジビン単位で電力比を平均しているので、階段状に電力比が変化している。フィルタ＃４では電力比が0dB付近であり、入力信号と出力信号がほとんど変化していないことがわかる。フィルタ次数選定のしきい値を1.0dBに設定して、この(b)の電力比と比較すると、フィルタ次数はレンジビン１〜５０までが３次、５０以降は２次と選定された。その時のクラッタ抑圧フィルタの出力信号が(c)であり、横軸がレンジビン、縦軸が平均電力を表す。受信信号中のクラッタが抑圧されて目標信号が検出可能になることがわかる。目標信号が存在するレンジビン付近以外のレンジビンでは、信号電力はほぼ受信機雑音電力程度であり、クラッタの消え残りはないことがわかる。
【００９１】
以上のように、１つのクラッタ中心周波数推定値に対して、複数の１次フィルタを縦続接続した構成にして、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。
【００９２】
実施の形態２．
図４はこの発明に係るクラッタ抑圧装置の実施の形態２を示す構成図である。図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、８はこのクラッタ中心周波数推定手段６から出力されるクラッタ中心周波数の推定値に対して、フィルタのノッチを形成するようにフィルタ係数を調整するフィルタ係数計算手段、９はこのフィルタ係数計算手段８で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、５はクラッタ中心周波数推定手段６から出力される推定値に対して、クラッタ抑圧性能を確保するために若干の周波数オフセットを与えるためのオフセット周波数を記憶するオフセット周波数記憶手段、７はフィルタ手段９の出力信号を記憶しておく出力信号記憶手段、１０は受信信号及びフィルタ手段９の出力信号を切り替えるサーキュレータ、３はフィルタ手段９の出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、４はフィルタ次数選定手段３から出力される選定結果に基づいて、出力信号記憶手段７からクラッタ抑圧フィルタとしての出力信号を選択する出力信号選択手段である。
【００９３】
次に図４を参照して動作について説明する。クラッタは単峰性であるとする。本実施の形態２は、基本的な動作は実施の形態１と同じであるので、動作が異なる点のみ説明する。
本実施の形態２では、サーキュレータ１０を設けてフィルタ係数計算手段８、及びフィルタ手段９への入力信号を切り替えることで、ハードウェア規模を小さくすることができる。サーキュレータ１０は、まず、受信信号をフィルタ係数計算手段８、及びフィルタ手段９へ出力する。受信信号のフィルタ処理が完了した時点で、フィルタ手段９の出力信号がサーキュレータ１０に送られるので、サーキュレータ１０は、フィルタ処理完了と同期させながらフィルタ手段９から出力される信号をフィルタ係数計算手段８、及びフィルタ手段９へ出力する。以降、フィルタ次数の上限値まで同様の処理を繰り返す。このとき、出力信号記憶手段７では、フィルタ手段９の出力信号を記憶しておく。フィルタ次数選定手段３で選定されたフィルタ次数は、出力信号選択手段４へ送られるので、出力信号選択手段４はこの結果に基づいて出力信号記憶手段７から該当する出力信号を選択して、クラッタ抑圧フィルタの出力信号として出力する。
また、フィルタ次数選定手段３の処理をフィルタ手段９の出力信号が得られる毎に実施すれば、フィルタ次数の上限値を設けることなくフィルタ次数の選定処理を行うことができる。
【００９４】
以上のように、１つのクラッタ中心周波数推定値に対して、１次フィルタでの処理を巡回的に実施する構成にして、フィルタ処理毎の入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができる。
【００９５】
実施の形態３．
図５はこの発明に係るクラッタ抑圧装置の実施の形態３を示す構成図である。図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、８はこのクラッタ中心周波数推定手段６から出力されるクラッタ中心周波数の推定値に対して、フィルタのノッチを形成するようにフィルタ係数を調整するフィルタ係数計算手段、９はこのフィルタ係数計算手段８で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、５はクラッタ中心周波数推定手段６から出力される推定値に対して、クラッタ抑圧性能を確保するために若干の周波数オフセットを与えるためのオフセット周波数記憶手段、１１はフィルタ係数計算手段８で算出される係数を記憶しておくフィルタ係数記憶手段、１０は受信信号及びフィルタ手段９の出力信号を切り替えるサーキュレータ、３はフィルタ手段９の出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、
１２はフィルタ次数選定手段３から出力される選定結果に基づいて、フィルタ係数記憶手段１１が選択したフィルタ係数を用いてクラッタ抑圧処理を行う第２のフィルタ手段である。
【００９６】
次に図５を参照して動作について説明する。クラッタは単峰性であるとする。本実施の形態３は、基本的な動作は実施の形態１及び実施の形態２と同じであるので、動作が異なる点のみ説明する。
実施の形態２では、フィルタ処理を巡回的に実施する際、出力信号を記憶していたが、本実施の形態３では、フィルタ係数記憶手段１１において、使用したフィルタの係数を記憶するようにしている。フィルタ次数選定手段３で選定されたフィルタ次数に対応するフィルタ係数を第２のフィルタ手段１２へ出力し、フィルタ処理を実施することで出力信号を得る。
【００９７】
以上のように、１つのクラッタ中心周波数推定値に対して、１次フィルタでの処理を巡回的に実施する構成にして、フィルタ処理毎の入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができる。
【００９８】
実施の形態４．
図６はこの発明に係るクラッタ抑圧装置の実施の形態４を示す構成図である。図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、１ａはこのクラッタ中心周波数推定手段６から出力されるクラッタ中心周波数の推定値に対して、フィルタのノッチを形成するようにフィルタ係数を調整する第一のフィルタ係数計算手段、２ａはこのフィルタ係数計算手段１ａで算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第一のフィルタ手段、５はクラッタ中心周波数推定手段６の推定値に対して、クラッタ抑圧性能を確保するために若干の周波数オフセットを与えるためのオフセット周波数を記憶するオフセット周波数記憶手段、１３はオフセット周波数を設定するフィルタ番号に応じて、オフセット周波数値をオフセット周波数記憶手段５から選択して出力するオフセット周波数選択手段、１４は、オフセット周波数記憶手段５およびオフセット周波数選択手段１３を含むオフセット周波数設定手段、１ｂ〜１Ｎはクラッタ中心周波数推定手段６から出力されるクラッタ中心周波数の推定値に、オフセット周波数記憶手段５から出力されるオフセット周波数を加えた周波数値に対して、フィルタのノッチを形成するようにフィルタ係数を調整する第二〜第Ｎのフィルタ係数計算手段、２ｂ〜２Ｎはこのフィルタ係数計算手段１ｂ〜１Ｎで算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第二〜第Ｎのフィルタ手段、１９は受信信号及び第一〜第Ｎのフィルタ手段２ａ〜２Ｎの出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、４はフィルタ次数選定手段３から出力される選定結果に基づいて、クラッタ抑圧フィルタとしての出力信号を選択する出力信号選択手段である。
図７はこの発明の実施の形態４の動作を説明するためのフローチャートである。
図８はこの発明の実施の形態４の動作を説明するための図である。
【００９９】
次に図６〜８を参照して動作について説明する。クラッタは双峰性であるとする。本実施の形態４において、クラッタ中心周波数推定手段６、第一〜第Ｎのフィルタ係数計算手段１ａ〜１Ｎ、第一〜第Ｎのフィルタ手段２ａ〜２Ｎ、出力信号選択手段４の動作は実施の形態１と同様であるため、ここでは動作の異なるオフセット周波数設定手段１４、及びフィルタ次数選定手段１９の動作を中心に説明する。
【０１００】
スタガトリガ方式のレーダでは、各フィルタは時変係数フィルタになるので、実施の形態１と同様に前段のフィルタ係数を考慮しながら設計する必要がある。図６において、フィルタ＃1〜フィルタ＃Ｎは全て１次フィルタであるとする。ここでは受信クラッタとして、帯域幅の差が大きくない２つの異なるスペクトルピークを有する双峰性クラッタを想定する。本実施の形態４で実施するフィルタ次数選定方法は、各フィルタの出力信号電力を基準として用いるため、連続するフィルタで１種類のクラッタを集中的に抑圧しても電力変化を検出できない可能性がある。つまり、２つのクラッタが等電力で受信されたとすると、複数のフィルタの零点を１つのクラッタの中心周波数に集中して設定してこのクラッタを完全に抑圧しても、他方のクラッタは抑圧されずそのまま残るため、受信信号中のクラッタの総電力は高々半分にしかならないためである。そこで、第一〜第Ｎのフィルタ係数計算手段１ａ〜１Ｎでは、２つのクラッタの中心周波数付近に対して零点を交互に割り当てるように動作する。例えば、奇数番号のフィルタはクラッタ１の中心周波数fc1付近に対して零点を設定し、偶数番号のフィルタはクラッタ２の中心周波数fc2付近に対して零点を設定する。このようにすると、２つのクラッタの相対電力比によって多少の変化はあるものの、処理する信号中にクラッタ成分が存在していれば、フィルタ１つおきに入出力信号電力の変化が大きくなる。この変化を検出すれば良い。これに伴い、オフセット周波数設定手段１４は偶数番号フィルタ用のオフセット周波数と、奇数番号フィルタ用のオフセット周波数とを分けて管理する必要がある。オフセット周波数は１つのクラッタ中心周波数に対して重複しなければよいので、Ｎ個の時変フィルタが縦続接続されている場合、オフセット周波数はｍ＝（Ｎ−２）／２種類用意すればよい。具体的には、オフセット周波数をフィルタ＃３，＃５，…，フィルタ＃Ｎ−1に対して、それぞれδ1，δ2，…，δmを設定し，同様に、フィルタ＃４，＃６，…，フィルタ＃Ｎに対しても、それぞれδ1，δ2，…，δmを設定する。
【０１０１】
次にフィルタ次数選定手段１９の動作を図７のフローチャートに従って説明する。
まず、ステップ３０でフィルタ次数のカウンタを初期化する。ステップ３１では、フィルタ次数がフィルタ次数の上限値NFIL以下であればステップ３２へ処理が移行し、そうでなければステップ３５へ処理が移行する。ステップ３２〜３４では、それぞれ各フィルタの入出力信号の平均電力、電力比を計算する。このとき、レンジビン毎の電力比のばらつきをおさえるため、複数のレンジビンで平均化処理を行う。以上の処理で各フィルタの入出力電力比が計算される。ステップ３５では仮のフィルタ次数Ｍに初期値１を代入する。ステップ３６ではフィルタ次数がフィルタ次数上限値NFILになると、処理がステップ４０に移行するように判定する。そうでない場合、ステップ３７〜３９ではステップ３３で計算した各フィルタの入出力信号電力比と一つ目のしきい値THLD1を比較して、電力比がしきい値を上回る最大の次数をordに代入する。ステップ４０以降は、受信したクラッタの相対電力比が大きい場合に発生する可能性があるクラッタの消え残りを回避するための処理である。まず、ステップ４０で仮の追加次数Lに上限値２を代入し、現時点での選定次数ordをFordに代入する。ステップ４１〜４４では、追加次数を２、１と変化させた時のフィルタ入出力電力比と２つ目のしきい値THLD2を比較して、電力比がTHLD2を上回ることがないかどうかをチェックする。電力比がTHLD2を上回る場合には、ステップ４３で選定するフィルタ次数を追加してFordに代入する。ステップ４５では、Fordの結果から選定するフィルタ次数を決定して、出力信号選択手段４へ出力する。
出力信号選択手段４では、出力されたフィルタ次数を基に出力信号を選択して、クラッタ抑圧フィルタの出力信号として出力する。
【０１０２】
図８に計算例を示す。クラッタ抑圧フィルタへの入力信号（受信信号）は，疑似的に発生させた双峰性クラッタ，目標信号と受信機雑音を加算した複素信号とした。スタガトリガ数を3、パルス繰り返し周期を400,300,500μs（平均400μs）、目標信号のドップラ周波数を500Hz、目標信号が存在するレンジビンを210レンジビン、目標信号電力を15dB、双峰性クラッタの中心周波数をそれぞれ0Hz、-600Hz、ガウス型のスペクトルを持つクラッタの帯域幅を双方とも平均30Hz、クラッタ電力を双方とも140-40log10(100+r)dB（rはレンジビン番号）、受信機雑音電力を0dB（平均0の正規分布）、フィルタの最大数を７として、計算機シミュレーションを行った。(a)は受信信号及び各フィルタの出力信号の変化を示すグラフで、横軸がレンジビン、縦軸が平均電力を表す。フィルタ＃５以降のフィルタでは、出力信号電力がほとんど雑音電力付近であることがわかる。(b)は各フィルタの入出力信号電力比の変化を示すグラフで、横軸がレンジビン、縦軸がフィルタ入出力信号電力比を表す。ここでは２５レンジビン単位で電力比を平均しているので、階段状に電力比が変化している。フィルタ＃５、６では電力比が0dB付近であり、入力信号と出力信号がほとんど変化していないことがわかる。フィルタ次数選定のしきい値を1.0dBに設定して、この(b)の電力比変化と比較すると、フィルタ次数はレンジビン１〜１５０までが６次、１５０以降は４次と選定された。その時のクラッタ抑圧フィルタの出力信号が(c)であり、横軸がレンジビン、縦軸が平均電力を表す。クラッタ抑圧処理により、受信信号中のクラッタが抑圧されて目標信号が検出可能になることがわかる。目標信号が存在するレンジビン付近以外のレンジビンでは、信号電力はほぼ受信機雑音電力程度であり、クラッタの消え残りはないことがわかる。
また、以上の説明では、本実施の形態４のクラッタ抑圧装置は双峰性クラッタに対して動作するように説明したが、Ｋ個のクラッタを受信した場合は、縦続接続したフィルタで設定する零点周波数をＫ個のクラッタ中心周波数推定値ごとに順番に割り振ること、及びオフセット周波数の組合せをＫ種類用意することで、３つ以上のクラッタにも対処することが可能である。
【０１０３】
以上のように、複数の１次フィルタを縦続接続した構成にして、縦続接続されたフィルタの零点を交互に２つのクラッタ中心周波数推定値に対して設定するようにし、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。
【０１０４】
実施の形態５．
図９はこの発明に係るクラッタ抑圧装置の実施の形態５を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、８はこのクラッタ中心周波数推定手段６から出力される２つのクラッタ中心周波数の推定値に対して、交互にフィルタの零点を設定するようにフィルタ係数を調整するフィルタ係数計算手段、９はこのフィルタ係数計算手段８で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、１４はクラッタ中心周波数推定手段６から出力されるクラッタ中心周波数の推定値に対して、クラッタ抑圧性能を確保するようにフィルタ毎に若干の周波数オフセットを選択して与えるためのオフセット周波数設定手段、７はフィルタ手段９の出力信号を記憶しておく出力信号記憶手段、１０は受信信号及びフィルタ手段９の出力信号を切り替えるサーキュレータ、１９はフィルタ手段９の出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、４はフィルタ次数選定手段３から出力される選定結果に基づいて、出力信号記憶手段７からクラッタ抑圧フィルタとしての出力信号を選択する出力信号選択手段である。
【０１０５】
次に図９を参照して動作について説明する。クラッタは双峰性であるとする。本実施の形態５は、基本的な動作は実施の形態２及び４と同じであるため、動作が異なる点を中心に説明する。
本実施の形態５では、サーキュレータ１０を設けてフィルタ係数計算手段８、及びフィルタ手段９への入力信号を切り替えることで、ハードウェア規模を小さくすることができる。サーキュレータ１０は、まず、受信信号をフィルタ係数計算手段８、及びフィルタ手段９へ出力する。受信信号のフィルタ処理が完了した時点で、フィルタ手段９の出力信号がサーキュレータ１０に出力されるので、サーキュレータ１０は、フィルタ処理完了と同期させながらフィルタ手段９から出力される信号をフィルタ係数計算手段８、及びフィルタ手段９へ出力する。以降、フィルタ次数の上限値まで同様の処理を繰り返す。ただし、実施の形態４に示したように、双峰性クラッタを受信した場合を想定しているので、２つのクラッタ中心周波数に対してフィルタの零点を交互に設定する必要がある。出力信号記憶手段７では、フィルタ手段９の出力信号を記憶しておく。オフセット周波数設定手段１４は、実施の形態４と同様な動作を行い、適用するフィルタごとにオフセット周波数を適宜選択してフィルタ係数計算手段８へ出力する。フィルタ次数選定手段１９では、実施の形態４に示した方法と同様な方法でフィルタ次数が選定される。選定されたフィルタ次数は、出力信号選択手段４へ出力されるので、出力信号選択手段４はこの結果に基づいて出力信号記憶手段７から該当する出力信号を選択して、クラッタ抑圧フィルタの出力信号として出力する。
また、フィルタ次数選定手段１９の処理をフィルタ手段９の出力信号が得られる毎に実施すれば、フィルタ次数の上限値を設けることなくフィルタ次数の選定処理を行うことができる。
【０１０６】
以上のように、双峰性クラッタの２つの中心周波数推定値のうち、一方のクラッタ中心周波数推定値に対しては次数を固定したフィルタで抑圧処理を行い、その出力信号に対してサーキュレータと１個の１次フィルタにより複数の１次フィルタを縦続接続した構成と同等の構成を構築して、縦続接続されたフィルタの零点を交互に２つのクラッタ中心周波数推定値に対して設定するようにし、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができる。
【０１０７】
実施の形態６．
図１０はこの発明に係るクラッタ抑圧装置の実施の形態６を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、８はこのクラッタ中心周波数推定手段６から出力される２つのクラッタ中心周波数の推定値に対して、交互にフィルタの零点を設定するようにフィルタ係数を調整するフィルタ係数計算手段、９はこのフィルタ係数計算手段８で算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第１のフィルタ手段、１４はクラッタ中心周波数推定手段６の推定値に対して、クラッタ抑圧性能を確保するようにフィルタ毎に若干の周波数オフセットを選択して与えるためのオフセット周波数設定手段、１１はフィルタ係数計算手段８で算出される係数を記憶しておくフィルタ係数記憶手段、１０は受信信号及び第１のフィルタ手段９の出力信号を切り替えるサーキュレータ、１９は第１のフィルタ手段９の出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、１２はフィルタ次数選定手段３から出力される選定結果に基づいて、フィルタ係数記憶手段１１から該当するフィルタ係数を選択し、その係数を用いてクラッタ抑圧処理を行う第２のフィルタ手段である。
【０１０８】
次に図１０を参照して動作について説明する。クラッタは双峰性であるとする。本実施の形態６は、基本的な動作は実施の形態５と同じであるため、動作が異なる点を中心に説明する。
実施の形態５では、フィルタ処理を巡回的に実施する際、出力信号を記憶していたが、本実施の形態６では、フィルタ係数記憶手段１１において、使用したフィルタの係数を記憶するようにしている。フィルタ次数選定手段１９で選定されたフィルタ次数に対応するフィルタ係数を第２のフィルタ手段１２へ出力し、フィルタ処理を実施することで出力信号を得る。
【０１０９】
以上のように、双峰性クラッタの２つの中心周波数推定値のうち、一方のクラッタ中心周波数推定値に対しては次数を固定したフィルタで抑圧処理を行い、その出力信号に対してサーキュレータと１個の１次フィルタにより複数の１次フィルタを縦続接続した構成と同等の構成を構築して、縦続接続されたフィルタの零点を交互に２つのクラッタ中心周波数推定値に対して設定するようにし、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができる。
【０１１０】
実施の形態７．
図１１はこの発明に係るクラッタ抑圧装置の実施の形態７を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、５０はクラッタ中心周波数推定手段から出力されるクラッタ中心周波数推定値の中の一つの周波数に対して、多重零点を設定するようにフィルタ係数を計算する多次フィルタ係数計算手段、５１は多次フィルタ係数計算手段５０で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、１ａはこのクラッタ中心周波数推定手段６から出力される推定値のうち、多次フィルタ係数計算手段５０へ出力された推定値以外の推定値に対して、フィルタのノッチを形成するようにフィルタ係数を調整する第一のフィルタ係数計算手段、２ａはこのフィルタ係数計算手段１ａで算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第一のフィルタ手段、５はクラッタ中心周波数推定手段６の推定値に対して、クラッタ抑圧性能を確保するために若干の周波数オフセットを与えるためのオフセット周波数記憶手段、１ｂ〜１Ｎはクラッタ中心周波数推定手段６から出力される推定値に、オフセット周波数記憶手段５から出力されるオフセット周波数を加えた周波数値に対して、フィルタのノッチを形成するようにフィルタ係数を調整する第二〜第Ｎのフィルタ係数計算手段、２ｂ〜２Ｎはこのフィルタ係数計算手段１ｂ〜１Ｎで算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第二〜第Ｎのフィルタ手段、３は受信信号及び第一〜第Ｎのフィルタ手段２ａ〜２Ｎの出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、４はフィルタ次数選定手段３から出力される選定結果に基づいて、クラッタ抑圧フィルタとしての出力信号を選択する出力信号選択手段である。
【０１１１】
次に図１１を参照して動作について説明する。クラッタは双峰性であるとする。
まず、図１１において他の実施例と同様に受信信号がクラッタ中心周波数推定手段６へ出力されて、従来のクラッタ抑圧装置と同様な方法でクラッタ中心周波数が推定される。
【０１１２】
クラッタ中心周波数推定手段６から出力される２つのクラッタ中心周波数の推定値のうち、一方が多次フィルタ係数計算手段５０へ出力される。双峰性クラッタの場合、海面からの不要反射エコーであるシークラッタか、地面からの不要反射エコーであるグランドクラッタと雲や雨からの不要反射エコーであるウェザクラッタが受信されることが多い。この中で、シークラッタとグランドクラッタに関しては、送信ビームの放射角度や気象情報等からその電力や帯域幅を経験的に概ね推定可能である場合が多く、これらを抑圧するフィルタの次数を予め見積もっておくことが可能である。本実施の形態７では、双峰性クラッタを受信した場合、このタイプのクラッタに対する抑圧処理はフィルタの次数が固定のフィルタで行い、残ったクラッタに対しては次数選定処理を行う。一般に、グランドクラッタやシークラッタはほとんどドップラ周波数が発生しないため、クラッタ中心周波数推定手段６で推定される周波数のうち、零に近い推定値を多次フィルタ係数計算手段５０へ出力する。
多次フィルタ係数計算手段５０では、クラッタ中心周波数推定手段６から出力された推定値に多重ノッチを形成するように時変フィルタ係数を計算する。算出された係数はフィルタ手段５１へ出力され、クラッタ抑圧処理が実施される。この処理により双峰性クラッタのうち一方が抑圧されるので、フィルタ手段５１の出力信号中のクラッタは、単峰性クラッタとして扱うことができる。
【０１１３】
クラッタ中心周波数推定手段６から第一〜第Ｎのフィルタ係数計算手段１ａ〜１Ｎへ送られる推定値が、多次フィルタ係数計算手段へ出力された推定値以外の推定値である点が異なるだけで、以降の動作は、実施の形態１と同様である。
【０１１４】
以上のように、双峰性クラッタの２つの中心周波数推定値のうち、一方のクラッタ中心周波数推定値に対しては次数を固定したフィルタで抑圧処理を行い、その出力信号に対して複数の１次フィルタを縦続接続した構成にして、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、双峰性クラッタを受信した場合でも、単峰性クラッタ受信時と同等な処理でクラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。
【０１１５】
実施の形態８．
図１２はこの発明に係るクラッタ抑圧装置の実施の形態８を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、５０はクラッタ中心周波数推定手段から出力されるクラッタ中心周波数の推定値の中の一つの周波数に対して、多重零点を設定するようにフィルタ係数を計算する多次フィルタ係数計算手段、５１は多次フィルタ係数計算手段５０で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、８はクラッタ中心周波数推定手段６から出力される他のクラッタ中心周波数の推定値に対して、フィルタのノッチを形成するようにフィルタ係数を調整するフィルタ係数計算手段、９はこのフィルタ係数計算手段８で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、５はクラッタ中心周波数推定手段６から出力されるクラッタ中心周波数の推定値に対して、クラッタ抑圧性能を確保するために若干の周波数オフセットを与えるためのオフセット周波数を記憶するオフセット周波数記憶手段、７はフィルタ手段９の出力信号を記憶しておく出力信号記憶手段、１０はフィルタ手段５１の出力信号及びフィルタ手段９の出力信号を切り替えるサーキュレータ、３はフィルタ手段９の出力信号から各フィルタ手段の入出力信号電力比を計算し、その入出力信号電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、４はフィルタ次数選定手段３から出力される選定結果に基づいて、出力信号記憶手段７からクラッタ抑圧フィルタとしての出力信号を選択する出力信号選択手段である。
【０１１６】
次に図１２を参照して動作について説明する。クラッタは双峰性であるとする。本実施の形態８は、基本的な動作は実施の形態２及び７と同等であるため、動作が異なる点を中心に説明する。
クラッタ中心周波数推定手段６から出力される２つの推定値のうち、一方が多次フィルタ係数計算手段５０に入力される。多次フィルタ係数計算手段５０では、入力された推定値に多重ノッチを形成するように時変フィルタ係数を計算する。算出されたフィルタ係数はフィルタ手段５１に入力され、クラッタ抑圧処理が実施される。この処理により双峰性クラッタのうち一方が抑圧されるので、フィルタ手段５１の出力信号中のクラッタは、単峰性クラッタとして扱うことができる。
【０１１７】
クラッタ中心周波数推定手段６からフィルタ係数計算手段８へ送出されるクラッタ中心周波数の推定値が、多次フィルタ係数計算手段へ送出されたクラッタ中心周波数の推定値以外の推定値である点、及びサーキュレータ１０の入力信号の一つが受信信号ではなく、フィルタ手段５１からの出力信号である点が異なるだけで、以降の動作は、実施の形態２と同様である。
【０１１８】
以上のように、双峰性クラッタの２つの中心周波数推定値のうち、一方のクラッタ中心周波数推定値に対しては次数を固定したフィルタで抑圧処理を確実に行い、この一方のクラッタ中心周波数が抑圧された出力信号に対してサーキュレータと１個の１次フィルタにより複数の１次フィルタを縦続接続した構成と同等の構成を構築して、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、双峰性クラッタを受信した場合でも、単峰性クラッタ受信時と同等な処理でクラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができる。
【０１１９】
実施の形態９．
図１３はこの発明に係るクラッタ抑圧装置の実施の形態９を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、５０はクラッタ中心周波数推定手段６から出力されるクラッタ中心周波数推定値の中の一つの周波数に対して、多重零点を設定するようにフィルタ係数を計算する多次フィルタ係数計算手段、５１は多次フィルタ係数計算手段５０で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、８はこのクラッタ中心周波数推定手段６から出力される推定値に対して、フィルタのノッチを形成するようにフィルタ係数を調整するフィルタ係数計算手段、９はこのフィルタ係数計算手段８で算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第一のフィルタ手段、５はクラッタ中心周波数推定手段６の推定値に対して、クラッタ抑圧性能を確保するために若干の周波数オフセットを与えるためのオフセット周波数記憶手段、１１はフィルタ係数計算手段８で算出される係数を記憶しておくフィルタ係数記憶手段、１０はフィルタ手段５１の出力信号及び第一のフィルタ手段９の出力信号を切り替えるサーキュレータ、３は第一のフィルタ手段９の出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、１２はフィルタ次数選定手段３から出力される選定結果に基づいて、フィルタ係数記憶手段１１から該当するフィルタ係数を選択し、その係数を用いてクラッタ抑圧処理を行う第２のフィルタ手段である。
【０１２０】
次に図１３を参照して動作について説明する。クラッタは双峰性であるとする。本実施の形態９は、基本的な動作は実施の形態３及び７と同等であるため、動作が異なる点を中心に説明する。
クラッタ中心周波数推定手段６から出力される２つのクラッタ中心周波数の推定値のうち、一方が多次フィルタ係数計算手段５０に入力される。多次フィルタ係数計算手段５０では、入力されたクラッタ中心周波数の推定値に多重ノッチを形成するように時変フィルタ係数を計算する。算出された係数はフィルタ手段５１に入力され、クラッタ抑圧処理が実施される。この処理により双峰性クラッタのうち一方が抑圧されるので、フィルタ手段５１の出力信号中のクラッタは、単峰性クラッタとして扱うことができる。
クラッタ中心周波数推定手段６からフィルタ係数計算手段８へ送られるクラッタ中心周波数の推定値が、多次フィルタ係数計算手段へ出力されたクラッタ中心周波数の推定値以外の推定値である点、及びサーキュレータ１０の入力信号の一つが受信信号ではなく、フィルタ手段５１の出力信号である点が異なるだけで、以降の動作は、実施の形態３と同様である。
【０１２１】
以上のように、双峰性クラッタの２つの中心周波数の推定値のうち、一方のクラッタ中心周波数の推定値に対しては次数を固定したフィルタで抑圧処理を確実に行い、この一方のクラッタ中心周波数が抑圧された出力信号に対してサーキュレータと１個の１次フィルタにより複数の１次フィルタを縦続接続した構成と同等の構成を構築して、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、双峰性クラッタを受信した場合でも、単峰性クラッタ受信時と同等な処理でクラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができる。
実施の形態１０．
図１５はこの発明に係るクラッタ抑圧装置の実施の形態１０を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、１ａはこのクラッタ中心周波数推定手段６から出力される推定値に対して、フィルタのノッチを形成するようにフィルタ係数を調整する第一のフィルタ係数計算手段、２ａはこのフィルタ係数計算手段１ａで算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第一のフィルタ手段、７１〜７０＋Ｎ−１はクラッタ中心周波数推定手段６の推定値に対して、クラッタ抑圧性能を確保するために多重零点が形成されるようにフィルタ係数を計算する第１〜第Ｎ−１の多重零点処理手段、２ｂ〜２Ｎは上記第１〜第Ｎ−１の多重零点処理手段で算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第二〜第Ｎのフィルタ手段、３は受信信号及び第一〜第Ｎのフィルタ手段２ａ〜２Ｎの出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、４はフィルタ次数選定手段３から出力される選定結果に基づいて、クラッタ抑圧フィルタとしての出力信号を選択する出力信号選択手段である。
図１６はこの発明に係るクラッタ抑圧装置の実施の形態１０の動作を説明するためのフローチャートである。
図１７はこの発明に係るクラッタ抑圧装置の実施の形態１０の動作を説明するための説明図である。
【０１２２】
次に図１５〜１７を参照して動作について説明する。離散フーリエ変換等による周波数解析が行えないほどパルスヒット数が少ないレーダにおいて、受信するクラッタ数を推定する方法としては、以下で説明するクラッタ中心周波数推定手段６で計算されるＡＲ係数aijを利用する方法等が考えられるが、本実施の形態１０では、クラッタ数は既知であるものとする。また、本実施の形態１０では、クラッタは単峰性であるとする。
まず、図１５において受信信号がクラッタ中心周波数推定手段６に送られて、従来のクラッタ抑圧装置と同様な方法でクラッタ中心周波数が推定される。
【０１２３】
次に、クラッタ中心周波数推定手段６から出力されるクラッタ中心周波数の推定値に従ってクラッタ抑圧フィルタの係数が計算される。従来のクラッタ抑圧装置と異なる点は、最適なフィルタ次数を自動的に設定するために、従来単一のフィルタで実施していたクラッタ抑圧処理を、複数の１次フィルタの縦続接続構成にし、且つ、同じ周波数に複数のフィルタでノッチを形成する場合に２段目以降のフィルタ計算時に零点多重度が確保できるようにしたことである。本実施の形態１０では縦続接続したフィルタの次数は１次としたが、２次以上に設定してもかまわない。
【０１２４】
スタガトリガ方式では不等間隔パルスになるため、ノッチ周波数が０以外のフィルタを実現するには、フィルタ係数がパルス間隔に応じて変わる時変係数フィルタが必要である。このような時変係数フィルタを縦続接続した場合のインパルス応答は、２つのフィルタのインパルス応答の単純な畳み込みにはならない。縦続接続した後段のフィルタ処理では、後段のフィルタの入力信号における現時点でのサンプルとそれ以前のサンプルとでは、それらのサンプルを算出するために用いた前段のフィルタ係数が異なるからである。
本実施の形態１０では実施の形態１と同様に縦続接続したフィルタの数に等しい数の周波数に対してノッチを形成するようなクラッタ抑圧フィルタを構成することが可能である。
本実施の形態１０では、単峰性のクラッタを想定しているため、フィルタ係数計算手段１ａに送られるクラッタ中心周波数推定値と、第１〜第Ｎ−１の多重零点処理手段７１〜７０＋Ｎ−１に送られるクラッタ中心周波数推定値は同一の周波数になる。つまり、クラッタ中心周波数の推定値をf0とすると、f01=f02=…=f0として各フィルタ係数が計算される。
しかしながら、以上に記述した方法では、複数のフィルタで零点を割り当てる周波数を重複して設定する場合、後段のフィルタで割り当てた零点の多重度のみが支配的になり、縦続接続したフィルタ次数に相当するフィルタ特性が得られないことが経験的に知られている。具体的には、フィルタのノッチの帯域幅が狭くなるため、従来のクラッタ抑圧装置に比べてクラッタ抑圧性能が劣化してしまう。
【０１２５】
そこで、本実施の形態１０では、図１５に示すように多重零点処理手段７１〜７０＋Ｎ−１を設け、フィルタ係数計算手段１ａで零点が割り当てられた、単峰性クラッタの中心周波数推定値f0に対して、後段の第一の多重零点処理手段７１以降において零点多重度が確保できるように、フィルタ係数を計算するようにしている。
以下に零点多重度を確保する手順を具体的に示す。
時変フィルタでは、同一の零点周波数を持つフィルタを縦続接続すると零点多重度が保存されない。フィルタの零点多重度は、設計時に式（２８）〜（２９）に示したフィルタ伝達関数の導関数をどこまで使用するかが影響する。例えば、５次の時変フィルタを２次と３次の時変フィルタの縦続接続で実現しようとすると、５次のフィルタでは伝達関数に関する０階〜４階の導関数を使用してフィルタ係数を計算するが、縦続接続構成では、０階〜２階の導関数しか使用されない。伝達関数に関する高次の微係数が０になるようにフィルタ係数を計算することは、その周波数に対して周波数特性がフラットに近くなる、即ちノッチフィルタのノッチ幅が広くなる方向へ作用するため、このままではフィルタ次数の合計は同じであっても、単一フィルタと複数フィルタの縦続構成とでは周波数特性は一致しない。
このような現象に対応するために、図１５の多重零点処理手段７１〜７０＋Ｎ−１では、クラッタ中心周波数推定手段６から送られる推定値が、前段のフィルタで零点を設定した周波数と一致しているかどうかを判定し、一致していなければ零点多重度が保存されるので、時変フィルタを縦続接続した場合の通常の設計方法によりフィルタ係数が計算される。一致している場合は零点多重度が保存さないので、零点多重度を確保する必要がある。そこで、零点多重度を確保するため、以下の計算方法に基づいてフィルタ係数計算が行われる。
（１）フィルタ#1、フィルタ#2の次数をそれぞれK1、K2とすると、フィルタ#2の係数計算においては、フィルタ伝達関数のK1から（K1+K2-1）階導関数を使用する。
（２）基本的な設計方法は、実施の形態１（但し、オフセット処理は除く）と同様である。
従って、式（３４）の行列A、式（３５）のベクトルbを、それぞれ式（３７）、式（３８）に示すように置き換え、式（３９）に示す一次連立方程式解くことでフィルタ#2の係数glを決定する。
【０１２６】
【数２１】
0003932957.tif 000022
【０１２７】
以上の計算方法により、周波数f0に対して等価的に（K1+K2）重の零点が形成されるため、従来のクラッタ抑圧装置におけるフィルタと同等な振幅２乗特性を得ることができる。
【０１２８】
図１５に示すように、１次の時変係数フィルタを多数縦続接続すると、各フィルタの入力信号中のクラッタ電力が大きく、フィルタのノッチがクラッタの中心周波数近辺に形成されていれば、各フィルタ処理後の出力信号電力は入力信号電力よりも小さくなる。クラッタがある程度抑圧されてしまえば、新たにフィルタ処理をしても出力信号電力が大きく変化することはないので、この電力変化を検出してフィルタの次数を選定する。各フィルタの入出力信号電力比はフィルタの番号が２ａ〜２Ｎに向かって大きくなるに連れて次第に小さくなり、その変化も次第に小さくなる。従って、フィルタの番号を横軸にし、入出力電力比を縦軸にすると、電力比の頂点は下に凸で遠方はほぼ平らになるなめらかな曲線上に乗る。従って、判定方法は、任意に設定できるしきい値と各フィルタの入出力信号電力比を比較し、電力比がしきい値を下回った状態はフィルタ次数が過剰と判断できるので、電力比がしきい値を初めて上回った時のフィルタの出力信号をクラッタ抑圧フィルタの出力信号にすることで、必要十分なフィルタ次数でクラッタ抑圧処理を実施することができる。
以上のフィルタ次数選定処理を、図１６のフローチャートに従って説明する。まず、ステップ２０でフィルタ次数のカウンタを初期化する。ステップ２１では、フィルタ次数がフィルタ次数の上限値NFIL以下であればステップ２２へ処理が移行し、そうでなければステップ２５へ処理が移行する。ステップ２２、２３では、それぞれ各フィルタの入出力信号の平均電力、電力比を計算する。このとき、レンジビン毎の電力比のばらつきをおさえるため、複数のレンジビンで平均化処理を行う。以上の処理で各フィルタの入出力電力比が計算される。ステップ２５では仮のフィルタ次数Ｍに上限値NFILを代入する。ステップ２６ではフィルタ次数が１次と判定された時点で処理が終了するように判定する。ステップ２７ではステップ２３で計算した各フィルタの入出力信号電力比としきい値を比較して、電力比がしきい値を上回った時点で処理がステップ２９に移行する。そうでなければ、ステップ２８で仮のフィルタ次数Ｍを更新して、ステップ２６へ移る。即ち、電力比がしきい値を下回っている間は仮りのフィルタ次数Ｍを1つずつ減らし、電力比が初めてしきい値を上回ったときにステップ２９へ移る。ステップ２９では仮のフィルタ次数Ｍに基づいて、最適なフィルタ次数を選定して処理が終了する。
【０１２９】
図１７に計算例を示す。クラッタ抑圧フィルタへの入力信号（受信信号）は，疑似的に発生させた単峰性クラッタ，目標信号と受信機雑音を加算した複素信号とした。スタガトリガ数を3、パルス繰り返し周期を400,300,500μs（平均400μs）、目標信号のドップラ周波数を500Hz、目標信号が存在するレンジビンを210レンジビン、目標信号電力を35dB、単峰性クラッタの中心周波数を-600Hz、ガウス型のスペクトルを持つクラッタの帯域幅を平均50Hz、クラッタ電力を140-40log10(100+r)dB（rはレンジビン番号）、受信機雑音電力を20dB（平均0の正規分布）、フィルタの最大数を４として、計算機シミュレーションを行った。(a)は受信信号及び各フィルタの出力信号の変化を示すグラフで、横軸がレンジビン、縦軸が平均電力を表す。フィルタ＃３以降のフィルタでは、出力信号電力がほとんど雑音電力付近であることがわかる。(b)は各フィルタの入出力信号電力比の変化を示すグラフで、横軸がレンジビン、縦軸がフィルタ入出力信号電力比を表す。ここでは２５レンジビン単位で電力比を平均しているので、階段状に電力比が変化している。フィルタ＃４では電力比が0dB付近であり、入力信号と出力信号がほとんど変化していないことがわかる。フィルタ次数選定のしきい値を1.0dBに設定して、この(b)の電力比と比較すると、フィルタ次数はレンジビン１〜５０までが３次、５０以降は２次と選定された。その時のクラッタ抑圧フィルタの出力信号が(c)であり、横軸がレンジビン、縦軸が平均電力を表す。受信信号中のクラッタが抑圧されて目標信号が検出可能になることがわかる。目標信号が存在するレンジビン付近以外のレンジビンでは、信号電力はほぼ受信機雑音電力程度であり、クラッタの消え残りはないことがわかる。
【０１３０】
以上のように、１つのクラッタ中心周波数推定値に対して、複数の１次フィルタを縦続接続した構成にして、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。
【０１３１】
実施の形態１１．
図１８はこの発明に係るクラッタ抑圧装置の実施の形態１１を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、５はこのクラッタ中心周波数推定手段６から出力されるクラッタ中心周波数の推定値に対して、フィルタのノッチを形成するようにフィルタ係数を調整するフィルタ係数計算手段、９はこのフィルタ係数計算手段８で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、１８はクラッタ中心周波数推定手段６から出力される推定値に対して、クラッタ抑圧性能を確保するために多重零点を形成するようにフィルタ係数を計算する多重零点処理手段、７はフィルタ手段９の出力信号を記憶しておく出力信号記憶手段、１０は受信信号及びフィルタ手段９の出力信号を切り替えるサーキュレータ、３はフィルタ手段９の出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、４はフィルタ次数選定手段３から出力される選定結果に基づいて、出力信号記憶手段７からクラッタ抑圧フィルタとしての出力信号を選択する出力信号選択手段である。
【０１３２】
次に図１８を参照して動作について説明する。クラッタは単峰性であるとする。本実施の形態１１は、基本的な動作は実施の形態１０と同じであるので、動作が異なる点のみ説明する。
本実施の形態１１では、サーキュレータ１０を設けてフィルタ係数計算手段５、及びフィルタ手段９への入力信号を切り替えることで、ハードウェア規模を小さくすることができる。サーキュレータ１０は、まず、受信信号をフィルタ係数計算手段５、及びフィルタ手段９へ出力する。受信信号のフィルタ処理が完了した時点で、フィルタ手段９の出力信号がサーキュレータ１０に送られるので、サーキュレータ１０は、フィルタ処理完了と同期させながらフィルタ手段９から出力される信号をフィルタ係数計算手段５、及びフィルタ手段９へ出力する。このとき、フィルタ係数計算手段５において使用されたクラッタ中心周波数に対してフィルタ処理を実施する場合は、多重零点処理手段１８においてフィルタ係数が計算され、フィルタ手段９へ送られる。以降、フィルタ次数の上限値まで同様の処理を繰り返す。このとき、出力信号記憶手段７では、フィルタ手段９の出力信号を記憶しておく。フィルタ次数選定手段３で選定されたフィルタ次数は、出力信号選択手段４へ送られるので、出力信号選択手段４はこの結果に基づいて出力信号記憶手段７から該当する出力信号を選択して、クラッタ抑圧フィルタの出力信号として出力する。
また、フィルタ次数選定手段３の処理をフィルタ手段９の出力信号が得られる毎に実施すれば、フィルタ次数の上限値を設けることなくフィルタ次数の選定処理を行うことができる。
【０１３３】
以上のように、１つのクラッタ中心周波数推定値に対して、１次フィルタでの処理を巡回的に実施する構成にして、フィルタ処理毎の入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができる。
【０１３４】
実施の形態１２．
図１９はこの発明に係るクラッタ抑圧装置の実施の形態１２を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、５はこのクラッタ中心周波数推定手段６から出力されるクラッタ中心周波数の推定値に対して、フィルタのノッチを形成するようにフィルタ係数を調整するフィルタ係数計算手段、９はこのフィルタ係数計算手段５で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、１８はクラッタ中心周波数推定手段６から出力される推定値に対して、クラッタ抑圧性能を確保するために多重零点を形成するようにフィルタ係数を計算する多重零点処理手段、１１はフィルタ係数計算手段５及び多重零点処理手段１８で算出される係数を記憶しておくフィルタ係数記憶手段、１０は受信信号及びフィルタ手段９の出力信号を切り替えるサーキュレータ、３はフィルタ手段９の出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、１２はフィルタ次数選定手段３から出力される選定結果に基づいて、フィルタ係数記憶手段１１が選択したフィルタ係数を用いてクラッタ抑圧処理を行う第２のフィルタ手段である。
【０１３５】
次に図１９を参照して動作について説明する。クラッタは単峰性であるとする。本実施の形態１２は、基本的な動作は実施の形態１０及び実施の形態１１と同じであるので、動作が異なる点のみ説明する。
実施の形態１１では、フィルタ処理を巡回的に実施する際、出力信号を記憶していたが、本実施の形態１２では、フィルタ係数記憶手段１１において、使用したフィルタの係数を記憶するようにしている。フィルタ次数選定手段３で選定されたフィルタ次数に対応するフィルタ係数を第２のフィルタ手段１２へ出力し、フィルタ処理を実施することで出力信号を得る。
【０１３６】
以上のように、１つのクラッタ中心周波数推定値に対して、１次フィルタでの処理を巡回的に実施する構成にして、フィルタ処理毎の入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができる。
【０１３７】
実施の形態１３．
図２０はこの発明に係るクラッタ抑圧装置の実施の形態１３を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、１ａはこのクラッタ中心周波数推定手段６から出力されるクラッタ中心周波数の推定値に対して、フィルタのノッチを形成するようにフィルタ係数を調整するフィルタ係数計算手段、２ａはこのフィルタ係数計算手段１ａで算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第一のフィルタ手段、２３は受信クラッタ数に応じてフィルタに割り当てる零点周波数を選択する多重零点割当選択手段、７１〜７０＋Ｎ−１はクラッタ中心周波数推定手段６から出力される推定値に対して、クラッタ抑圧性能を確保するために多重零点が形成されるようにフィルタ係数を計算する第一〜第Ｎ−１の多重零点処理手段、２ｂ〜２Ｎはこの第一〜第Ｎ−１の多重零点処理手段７１〜７０＋Ｎ−１で算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第二〜第Ｎのフィルタ手段、１９は受信信号及び第一〜第Ｎのフィルタ手段２ａ〜２Ｎの出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、４はフィルタ次数選定手段３から出力される選定結果に基づいて、クラッタ抑圧フィルタとしての出力信号を選択する出力信号選択手段である。
図２１はこの発明の実施の形態１３の動作を説明するためのフローチャートである。
図２２はこの発明の実施の形態１３の動作を説明するための図である。
【０１３８】
次に図２０〜２２を参照して動作について説明する。クラッタは双峰性であるとする。本実施の形態１３において、クラッタ中心周波数推定手段６、第一〜第Ｎ−１の多重零点処理手段７１〜７０＋Ｎ−１、第一〜第Ｎのフィルタ手段２ａ〜２Ｎ、出力信号選択手段４の動作は実施の形態１０と同様であるため、ここでは動作の異なる多重零点割当選択手段２３、及びフィルタ次数選定手段１９の動作を中心に説明する。
【０１３９】
スタガトリガ方式のレーダでは、各フィルタは時変係数フィルタになるので、実施の形態１０と同様に前段のフィルタ係数を考慮しながら設計する必要がある。図２０において、フィルタ＃1〜フィルタ＃Ｎは全て１次フィルタであるとする。ここでは受信クラッタとして、帯域幅の差が大きくない２つの異なるスペクトルピークを有する双峰性クラッタを想定する。本実施の形態１３で実施するフィルタ次数選定方法は、各フィルタの出力信号電力を基準として用いるため、連続するフィルタで１種類のクラッタを集中的に抑圧しても電力変化を検出できない可能性がある。つまり、２つのクラッタが等電力で受信されたとすると、複数のフィルタの零点を１つのクラッタの中心周波数に集中して設定してこのクラッタを完全に抑圧しても、他方のクラッタは抑圧されずそのまま残るため、受信信号中のクラッタの総電力は高々半分にしかならないためである。そこで、フィルタ係数計算手段１ａ及び第一〜第Ｎ−１の多重零点処理手段７１〜７０＋Ｎ−１では、２つのクラッタの中心周波数付近に対して零点を交互に割り当てるように動作する。例えば、奇数番号のフィルタはクラッタ１の中心周波数fc1付近に対して零点を設定し、偶数番号のフィルタはクラッタ２の中心周波数fc2付近に対して零点を設定する。このようにすると、２つのクラッタの相対電力比によって多少の変化はあるものの、処理する信号中にクラッタ成分が存在していれば、フィルタ１つおきに入出力信号電力の変化が大きくなる。この変化を検出すれば良い。
【０１４０】
次にフィルタ次数選定手段１９の動作を図２１のフローチャートに従って説明する。
まず、ステップ３０でフィルタ次数のカウンタを初期化する。ステップ３１では、フィルタ次数がフィルタ次数の上限値NFIL以下であればステップ３２へ処理が移行し、そうでなければステップ３５へ処理が移行する。ステップ３２〜３４では、それぞれ各フィルタの入出力信号の平均電力、電力比を計算する。このとき、レンジビン毎の電力比のばらつきをおさえるため、複数のレンジビンで平均化処理を行う。以上の処理で各フィルタの入出力電力比が計算される。ステップ３５では仮のフィルタ次数Ｍに初期値１を代入する。ステップ３６ではフィルタ次数がフィルタ次数上限値NFILになると、処理がステップ４０に移行するように判定する。そうでない場合、ステップ３７〜３９ではステップ３３で計算した各フィルタの入出力信号電力比と一つ目のしきい値THLD1を比較して、電力比がしきい値を上回る最大の次数をordに代入する。ステップ４０以降は、受信したクラッタの相対電力比が大きい場合に発生する可能性があるクラッタの消え残りを回避するための処理である。まず、ステップ４０で仮の追加次数Lに上限値２を代入し、現時点での選定次数ordをFordに代入する。ステップ４１〜４４では、追加次数を２、１と変化させた時のフィルタ入出力電力比と２つ目のしきい値THLD2を比較して、電力比がTHLD2を上回ることがないかどうかをチェックする。電力比がTHLD2を上回る場合には、ステップ４３で選定するフィルタ次数を追加してFordに代入する。ステップ４５では、Fordの結果から選定するフィルタ次数を決定して、出力信号選択手段４へ出力する。
出力信号選択手段４では、出力されたフィルタ次数を基に出力信号を選択して、クラッタ抑圧フィルタの出力信号として出力する。
【０１４１】
図２２に計算例を示す。クラッタ抑圧フィルタへの入力信号（受信信号）は、疑似的に発生させた双峰性クラッタ、目標信号と受信機雑音を加算した複素信号とした。スタガトリガ数を3、パルス繰り返し周期を400,300,500μs（平均400μs）、目標信号のドップラ周波数を500Hz、目標信号が存在するレンジビンを210レンジビン、目標信号電力を15dB、双峰性クラッタの中心周波数をそれぞれ0Hz、-600Hz、ガウス型のスペクトルを持つクラッタの帯域幅を双方とも平均30Hz、クラッタ電力を双方とも140-40log10(100+r)dB（rはレンジビン番号）、受信機雑音電力を0dB（平均0の正規分布）、フィルタの最大数を７として、計算機シミュレーションを行った。(a)は受信信号及び各フィルタの出力信号の変化を示すグラフで、横軸がレンジビン、縦軸が平均電力を表す。フィルタ＃５以降のフィルタでは、出力信号電力がほとんど雑音電力付近であることがわかる。(b)は各フィルタの入出力信号電力比の変化を示すグラフで、横軸がレンジビン、縦軸がフィルタ入出力信号電力比を表す。ここでは２５レンジビン単位で電力比を平均しているので、階段状に電力比が変化している。フィルタ＃５、６では電力比が0dB付近であり、入力信号と出力信号がほとんど変化していないことがわかる。フィルタ次数選定のしきい値を1.0dBに設定して、この(b)の電力比変化と比較すると、フィルタ次数はレンジビン１〜１５０までが６次、１５０以降は４次と選定された。その時のクラッタ抑圧フィルタの出力信号が(c)であり、横軸がレンジビン、縦軸が平均電力を表す。クラッタ抑圧処理により、受信信号中のクラッタが抑圧されて目標信号が検出可能になることがわかる。目標信号が存在するレンジビン付近以外のレンジビンでは、信号電力はほぼ受信機雑音電力程度であり、クラッタの消え残りはないことがわかる。
また、以上の説明では、本実施の形態１３のクラッタ抑圧装置は双峰性クラッタに対して動作するように説明したが、Ｋ個のクラッタを受信した場合は、縦続接続したフィルタで設定する零点周波数をＫ個のクラッタ中心周波数推定値ごとに順番に割り振ることで、３つ以上のクラッタにも対処することが可能である。
【０１４２】
以上のように、複数の１次フィルタを縦続接続した構成にして、縦続接続されたフィルタの零点を交互に２つのクラッタ中心周波数推定値に対して設定するようにし、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。
【０１４３】
実施の形態１４．
図２３はこの発明に係るクラッタ抑圧装置の実施の形態１４を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、１８はこのクラッタ中心周波数推定手段６から出力される２つのクラッタ中心周波数の推定値に対して、多重零点が形成されるようにフィルタ係数を調整する多重零点処理手段、９はこの多重零点処理手段１８で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、２３は受信クラッタ数に応じてフィルタに割り当てる零点周波数を選択する多重零点割当選択手段、７はフィルタ手段９の出力信号を記憶しておく出力信号記憶手段、１０は受信信号及びフィルタ手段９の出力信号を切り替えるサーキュレータ、１９はフィルタ手段９の出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、４はフィルタ次数選定手段３から出力される選定結果に基づいて、出力信号記憶手段７からクラッタ抑圧フィルタとしての出力信号を選択する出力信号選択手段である。
【０１４４】
次に図２３を参照して動作について説明する。クラッタは双峰性であるとする。本実施の形態１４は、基本的な動作は実施の形態１１及び実施の形態１３と同じであるため、動作が異なる点を中心に説明する。
本実施の形態１４では、サーキュレータ１０を設けて多重零点処理手段１８、及びフィルタ手段９への入力信号を切り替えることで、ハードウェア規模を小さくすることができる。サーキュレータ１０は、まず、受信信号を多重零点処理手段１８、及びフィルタ手段９へ出力する。受信信号のフィルタ処理が完了した時点で、フィルタ手段９の出力信号がサーキュレータ１０に出力されるので、サーキュレータ１０は、フィルタ処理完了と同期させながらフィルタ手段９から出力される信号を多重零点処理手段１８、及びフィルタ手段９へ出力する。以降、フィルタ次数の上限値まで同様の処理を繰り返す。ただし、実施の形態１３に示したように、双峰性クラッタを受信した場合を想定しているので、２つのクラッタ中心周波数に対してフィルタの零点を交互に設定する必要がある。出力信号記憶手段７では、フィルタ手段９の出力信号を記憶しておく。多重零点割当選択手段２３は、実施の形態１３と同様な動作を行い、受信クラッタ数に応じて多重零点処理を行う零点周波数を選択する。フィルタ次数選定手段１９では、実施の形態１３に示した方法と同様な方法でフィルタ次数が選定される。選定されたフィルタ次数は、出力信号選択手段４へ出力されるので、出力信号選択手段４はこの結果に基づいて出力信号記憶手段７から該当する出力信号を選択して、クラッタ抑圧フィルタの出力信号として出力する。
また、フィルタ次数選定手段１９の処理をフィルタ手段９の出力信号が得られる毎に実施すれば、フィルタ次数の上限値を設けることなくフィルタ次数の選定処理を行うことができる。
【０１４５】
以上のように、双峰性クラッタの２つの中心周波数推定値のうち、一方のクラッタ中心周波数推定値に対しては次数を固定したフィルタで抑圧処理を行い、その出力信号に対してサーキュレータと１個の１次フィルタにより複数の１次フィルタを縦続接続した構成と同等の構成を構築して、縦続接続されたフィルタの零点を交互に２つのクラッタ中心周波数推定値に対して設定するようにし、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができる。
【０１４６】
実施の形態１５．
図２４はこの発明に係るクラッタ抑圧装置の実施の形態１５を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、１８は多重零点が形成されるようにフィルタ係数を調整する多重零点処理手段、９はこの多重零点処理手段１８で算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第一のフィルタ手段、２３は受信クラッタ数に応じてフィルタに割り当てる零点周波数を選択する多重零点割当選択手段、１１は多重零点処理手段１８で算出される係数を記憶しておくフィルタ係数記憶手段、１０は受信信号及び第一のフィルタ手段９の出力信号を切り替えるサーキュレータ、１９は第一のフィルタ手段９の出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、１２はフィルタ次数選定手段３から出力される選定結果に基づいて、フィルタ係数記憶手段１１から該当するフィルタ係数を選択し、その係数を用いてクラッタ抑圧処理を行う第２のフィルタ手段である。
【０１４７】
次に図２４を参照して動作について説明する。クラッタは双峰性であるとする。本実施の形態１５は、基本的な動作は実施の形態１４と同じであるため、動作が異なる点を中心に説明する。
実施の形態１４では、フィルタ処理を巡回的に実施する際、出力信号を記憶していたが、本実施の形態１５では、フィルタ係数記憶手段１１において、使用したフィルタの係数を記憶するようにしている。フィルタ次数選定手段１９で選定されたフィルタ次数に対応するフィルタ係数を第２のフィルタ手段１２へ出力し、フィルタ処理を実施することで出力信号を得る。
【０１４８】
以上のように、双峰性クラッタの２つの中心周波数推定値のうち、一方のクラッタ中心周波数推定値に対しては次数を固定したフィルタで抑圧処理を行い、その出力信号に対してサーキュレータと１個の１次フィルタにより複数の１次フィルタを縦続接続した構成と同等の構成を構築して、縦続接続されたフィルタの零点を交互に２つのクラッタ中心周波数推定値に対して設定するようにし、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができる。
【０１４９】
実施の形態１６．
図２５はこの発明に係るクラッタ抑圧装置の実施の形態１６を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、５０はクラッタ中心周波数推定手段から出力されるクラッタ中心周波数推定値の中の一つの周波数に対して、多重零点を設定するようにフィルタ係数を計算する多次フィルタ係数計算手段、５１は多次フィルタ係数計算手段５０で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、２３は受信クラッタ数に応じてフィルタに割り当てる零点周波数を選択する多重零点割当選択手段、７１〜７０＋Ｎ−１はクラッタ中心周波数推定手段６の推定値に対して、クラッタ抑圧性能を確保するために多重零点が形成されるようにフィルタ係数を計算する第一〜第Ｎ−１の多重零点処理手段、２ａ〜２Ｎはこの第一〜第Ｎ−１の多重零点処理手段７１〜７０＋Ｎ−１で算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第一〜第Ｎ−１のフィルタ手段、３は受信信号及び第一〜第Ｎ−１のフィルタ手段２ａ〜２Ｎの出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、４はフィルタ次数選定手段３から出力される選定結果に基づいて、クラッタ抑圧フィルタとしての出力信号を選択する出力信号選択手段である。
【０１５０】
次に図２５を参照して動作について説明する。クラッタは双峰性であるとする。
まず、図２５において他の実施例と同様に受信信号がクラッタ中心周波数推定手段６へ出力されて、従来のクラッタ抑圧装置と同様な方法でクラッタ中心周波数が推定される。
【０１５１】
クラッタ中心周波数推定手段６から出力される２つのクラッタ中心周波数の推定値のうち、一方が多次フィルタ係数計算手段５０へ出力される。双峰性クラッタの場合、海面からの不要反射エコーであるシークラッタか、地面からの不要反射エコーであるグランドクラッタと雲や雨からの不要反射エコーであるウェザクラッタが受信されることが多い。この中で、シークラッタとグランドクラッタに関しては、送信ビームの放射角度や気象情報等からその電力や帯域幅を経験的に概ね推定可能である場合が多く、これらを抑圧するフィルタの次数を予め見積もっておくことが可能である。本実施の形態１６では、双峰性クラッタを受信した場合、このタイプのクラッタに対する抑圧処理はフィルタの次数が固定のフィルタで行い、残ったクラッタに対しては次数選定処理を行う。一般に、グランドクラッタやシークラッタはほとんどドップラ周波数が発生しないため、クラッタ中心周波数推定手段６で推定される周波数のうち、零に近い推定値を多次フィルタ係数計算手段５０へ出力する。
多次フィルタ係数計算手段５０では、クラッタ中心周波数推定手段６から出力された推定値に多重ノッチを形成するように時変フィルタ係数を計算する。算出された係数はフィルタ手段５１へ出力され、クラッタ抑圧処理が実施される。この処理により双峰性クラッタのうち一方が抑圧されるので、フィルタ手段５１の出力信号中のクラッタは、単峰性クラッタとして扱うことができる。
【０１５２】
クラッタ中心周波数推定手段６から第一〜第Ｎ−１の多重零点処理手段７１〜７０＋Ｎ−１へ送られる推定値が、多次フィルタ係数計算手段へ出力された推定値以外の推定値である点が異なるだけで、以降の動作は、実施の形態１０と同様である。
【０１５３】
以上のように、双峰性クラッタの２つの中心周波数推定値のうち、一方のクラッタ中心周波数推定値に対しては次数を固定したフィルタで抑圧処理を行い、その出力信号に対して複数の１次フィルタを縦続接続した構成にして、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、双峰性クラッタを受信した場合でも、単峰性クラッタ受信時と同等な処理でクラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。
【０１５４】
実施の形態１７．
図２６はこの発明に係るクラッタ抑圧装置の実施の形態１７を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、５０はクラッタ中心周波数推定手段から出力されるクラッタ中心周波数の推定値の中の一つの周波数に対して、多重零点を設定するようにフィルタ係数を計算する多次フィルタ係数計算手段、５１は多次フィルタ係数計算手段５０で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、１８はこのクラッタ中心周波数推定手段６から出力される２つのクラッタ中心周波数の推定値に対して、多重零点が形成されるようにフィルタ係数を調整する多重零点処理手段、９はこの多重零点処理手段１８で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、７はフィルタ手段９の出力信号を記憶しておく出力信号記憶手段、１０はフィルタ手段５１の出力信号及びフィルタ手段９の出力信号を切り替えるサーキュレータ、３はフィルタ手段９の出力信号から各フィルタ手段の入出力信号電力比を計算し、その入出力信号電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、４はフィルタ次数選定手段３から出力される選定結果に基づいて、出力信号記憶手段７からクラッタ抑圧フィルタとしての出力信号を選択する出力信号選択手段である。
【０１５５】
次に図２６を参照して動作について説明する。クラッタは双峰性であるとする。本実施の形態１７は、基本的な動作は実施の形態１１及び実施の形態１６と同等であるため、動作が異なる点を中心に説明する。
クラッタ中心周波数推定手段６から出力される２つの推定値のうち、一方が多次フィルタ係数計算手段５０に入力される。多次フィルタ係数計算手段５０では、入力された推定値に多重ノッチを形成するように時変フィルタ係数を計算する。算出されたフィルタ係数はフィルタ手段５１に入力され、クラッタ抑圧処理が実施される。この処理により双峰性クラッタのうち一方が抑圧されるので、フィルタ手段５１の出力信号中のクラッタは、単峰性クラッタとして扱うことができる。
【０１５６】
クラッタ中心周波数推定手段６から多重零点処理手段１８へ送出されるクラッタ中心周波数の推定値が、多次フィルタ係数計算手段へ送出されたクラッタ中心周波数の推定値以外の推定値である点、及びサーキュレータ１０の入力信号の一つが受信信号ではなく、フィルタ手段５１からの出力信号である点が異なるだけで、以降の動作は、実施の形態１１と同様である。
【０１５７】
以上のように、双峰性クラッタの２つの中心周波数推定値のうち、一方のクラッタ中心周波数推定値に対しては次数を固定したフィルタで抑圧処理を確実に行い、この一方のクラッタ中心周波数が抑圧された出力信号に対してサーキュレータと１個の１次フィルタにより複数の１次フィルタを縦続接続した構成と同等の構成を構築して、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、双峰性クラッタを受信した場合でも、単峰性クラッタ受信時と同等な処理でクラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができる。
【０１５８】
実施の形態１８．
図２７はこの発明に係るクラッタ抑圧装置の実施の形態１８を示す構成図である。
図において、６は受信信号から受信信号中のクラッタ中心周波数を推定するクラッタ中心周波数推定手段、５０はクラッタ中心周波数推定手段６から出力されるクラッタ中心周波数推定値の中の一つの周波数に対して、多重零点を設定するようにフィルタ係数を計算する多次フィルタ係数計算手段、５１は多次フィルタ係数計算手段５０で算出されたフィルタ係数を用いてクラッタ抑圧処理を行うフィルタ手段、１８はこのクラッタ中心周波数推定手段６から出力される２つのクラッタ中心周波数の推定値に対して、多重零点が形成されるようにフィルタ係数を調整する多重零点処理手段、９はこの多重零点処理手段１８で算出されたフィルタ係数を用いてクラッタ抑圧処理を行う第一のフィルタ手段、１１はフィルタ係数計算手段１８で算出される係数を記憶しておくフィルタ係数記憶手段、１０は受信信号及び第一のフィルタ手段９の出力信号を切り替えるサーキュレータ、３は第一のフィルタ手段９の出力信号から各フィルタ手段の入出力信号電力比を計算し、その電力比を用いて最適なフィルタ次数を選定するフィルタ次数選定手段、１２はフィルタ次数選定手段３から出力される選定結果に基づいて、フィルタ係数記憶手段１１から該当するフィルタ係数を選択し、その係数を用いてクラッタ抑圧処理を行う第２のフィルタ手段である。
【０１５９】
次に図２７を参照して動作について説明する。クラッタは双峰性であるとする。本実施の形態１８は、基本的な動作は実施の形態１２及び実施の形態１６と同等であるため、動作が異なる点を中心に説明する。
クラッタ中心周波数推定手段６から出力される２つのクラッタ中心周波数の推定値のうち、一方が多次フィルタ係数計算手段５０に入力される。多次フィルタ係数計算手段５０では、入力されたクラッタ中心周波数の推定値に多重ノッチを形成するように時変フィルタ係数を計算する。算出された係数はフィルタ手段５１に入力され、クラッタ抑圧処理が実施される。この処理により双峰性クラッタのうち一方が抑圧されるので、フィルタ手段５１の出力信号中のクラッタは、単峰性クラッタとして扱うことができる。
クラッタ中心周波数推定手段６から多重零点処理手段１８へ送られるクラッタ中心周波数の推定値が、多次フィルタ係数計算手段へ出力されたクラッタ中心周波数の推定値以外の推定値である点、及びサーキュレータ１０の入力信号の一つが受信信号ではなく、フィルタ手段５１の出力信号である点が異なるだけで、以降の動作は、実施の形態１２と同様である。
【０１６０】
以上のように、双峰性クラッタの２つの中心周波数の推定値のうち、一方のクラッタ中心周波数の推定値に対しては次数を固定したフィルタで抑圧処理を確実に行い、この一方のクラッタ中心周波数が抑圧された出力信号に対してサーキュレータと１個の１次フィルタにより複数の１次フィルタを縦続接続した構成と同等の構成を構築して、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、双峰性クラッタを受信した場合でも、単峰性クラッタ受信時と同等な処理でクラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができる。
【０１６１】
なお、実施の形態１６〜１８に多重零点割り当て手段を適用してもよい。これにより、実施の形態１３と同様の効果を奏する。
【０１６２】
【発明の効果】
この発明によれば、１つのクラッタ中心周波数推定値に対して、複数の１次フィルタを縦続接続した構成にして、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができるという効果を奏する。
【０１６３】
また、この発明によれば、１つのクラッタ中心周波数推定値に対して、１次フィルタでの処理を巡回的に実施する構成にして、フィルタ処理毎の入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができるという効果を奏する。
【０１６４】
また、この発明によれば、１つのクラッタ中心周波数推定値に対して、１次フィルタでの処理を巡回的に実施する構成にして、フィルタ処理毎の入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができるという効果を奏する。
【０１６５】
また、この発明によれば、複数の１次フィルタを縦続接続した構成にして、縦続接続されたフィルタの零点を交互に２つのクラッタ中心周波数推定値に対して設定するようにし、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。
という効果を奏する。
【０１６６】
また、この発明によれば、双峰性クラッタの２つの中心周波数推定値のうち、一方のクラッタ中心周波数推定値に対しては次数を固定したフィルタで抑圧処理を行い、その出力信号に対してサーキュレータと１個の１次フィルタにより複数の１次フィルタを縦続接続した構成と同等の構成を構築して、縦続接続されたフィルタの零点を交互に２つのクラッタ中心周波数推定値に対して設定するようにし、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができるという効果を奏する。
【０１６７】
また、この発明によれば、双峰性クラッタの２つの中心周波数推定値のうち、一方のクラッタ中心周波数推定値に対しては次数を固定したフィルタで抑圧処理を行い、その出力信号に対してサーキュレータと１個の１次フィルタにより複数の１次フィルタを縦続接続した構成と同等の構成を構築して、縦続接続されたフィルタの零点を交互に２つのクラッタ中心周波数推定値に対して設定するようにし、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、クラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができるという効果を奏する。
【０１６８】
また、この発明によれば、双峰性クラッタの２つの中心周波数推定値のうち、一方のクラッタ中心周波数推定値に対しては次数を固定したフィルタで抑圧処理を行い、その出力信号に対して複数の１次フィルタを縦続接続した構成にして、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、双峰性クラッタを受信した場合でも、単峰性クラッタ受信時と同等な処理でクラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができるという効果を奏する。
【０１６９】
また、この発明によれば、双峰性クラッタの２つの中心周波数推定値のうち、一方のクラッタ中心周波数推定値に対しては次数を固定したフィルタで抑圧処理を確実に行い、この一方のクラッタ中心周波数が抑圧された出力信号に対してサーキュレータと１個の１次フィルタにより複数の１次フィルタを縦続接続した構成と同等の構成を構築して、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、双峰性クラッタを受信した場合でも、単峰性クラッタ受信時と同等な処理でクラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができるという効果を奏する。
【０１７０】
また、この発明によれば、双峰性クラッタの２つの中心周波数の推定値のうち、一方のクラッタ中心周波数の推定値に対しては次数を固定したフィルタで抑圧処理を確実に行い、この一方のクラッタ中心周波数が抑圧された出力信号に対してサーキュレータと１個の１次フィルタにより複数の１次フィルタを縦続接続した構成と同等の構成を構築して、各フィルタの入出力信号電力比をパラメータにしてクラッタを抑圧するのに最低限必要なフィルタ次数を選定するようにしているので、双峰性クラッタを受信した場合でも、単峰性クラッタ受信時と同等な処理でクラッタの受信状況に応じて最適なフィルタ次数を自動的に設定することが可能になり、フィルタ次数が不足することによるクラッタの消え残りや過剰なフィルタ処理による目標信号の減衰等を防ぐことができる。また、フィルタ処理を巡回的に実施することにより、ハードウェア規模を小さくすることができるという効果を奏する。
【０１７１】
また、この発明によれば、受信信号中のクラッタ数に応じて、多重零点処理を行うので、クラッタ抑圧性能の劣化を防ぐことができるという項かを奏する。
【図面の簡単な説明】
【図１】  この発明に係るクラッタ抑圧装置の実施の形態１を示す構成図である。
【図２】  この発明の実施の形態１の動作を説明するためのフローチャートである。
【図３】  この発明の実施の形態１の動作を説明するための説明図である。
【図４】  この発明に係るクラッタ抑圧装置の実施の形態２を示す構成図である。
【図５】  この発明に係るクラッタ抑圧装置の実施の形態３を示す構成図である。
【図６】  この発明に係るクラッタ抑圧装置の実施の形態４を示す構成図である。
【図７】  この発明の実施の形態４の動作を説明するためのフローチャートである。
【図８】  この発明の実施の形態４の動作を説明するための図である。
【図９】  この発明に係るクラッタ抑圧装置の実施の形態５を示す構成図である。
【図１０】  この発明に係るクラッタ抑圧装置の実施の形態６を示す構成図である。
【図１１】  この発明に係るクラッタ抑圧装置の実施の形態７を示す構成図である。
【図１２】  この発明に係るクラッタ抑圧装置の実施の形態８を示す構成図である。
【図１３】  この発明に係るクラッタ抑圧装置の実施の形態９を示す構成図である。
【図１４】  この発明に係るクラッタ抑圧装置の実施の形態１の動作を説明するためのブロック図である。
【図１５】  この発明に係るクラッタ抑圧装置の実施の形態１０を示す構成図である。
【図１６】  この発明に係るクラッタ抑圧装置の実施の形態１０の動作を説明するためのフローチャートである。
【図１７】  この発明に係るクラッタ抑圧装置の実施の形態１０の動作を説明するための説明図である。
【図１８】  この発明に係るクラッタ抑圧装置の実施の形態１１を示す構成図である。
【図１９】  この発明に係るクラッタ抑圧装置の実施の形態１２を示す構成図である。
【図２０】  この発明に係るクラッタ抑圧装置の実施の形態１３を示す構成図である。
【図２１】  この発明の実施の形態１３の動作を説明するためのフローチャートである。
【図２２】  この発明の実施の形態１３の動作を説明するための図である。
【図２３】  この発明に係るクラッタ抑圧装置の実施の形態１４を示す構成図である。
【図２４】  この発明に係るクラッタ抑圧装置の実施の形態１５を示す構成図である。
【図２５】  この発明に係るクラッタ抑圧装置の実施の形態１６を示す構成図である。
【図２６】  この発明に係るクラッタ抑圧装置の実施の形態１７を示す構成図である。
【図２７】  この発明に係るクラッタ抑圧装置の実施の形態１８を示す構成図である。
【図２８】  従来のクラッタ抑圧装置を示すブロック図である。
【符号の説明】
１ａ〜１Ｎ  フィルタ係数計算手段、１ｂ  フィルタ係数計算手段、２ａ〜２ｘフィルタ手段、７１〜７０＋Ｎ−１  第１〜第Ｎ−１の多重零点処理手段、２ｂ〜２Ｎ  第二〜第Ｎのフィルタ手段、３  フィルタ次数選定手段、４  出力信号選択手段、５  オフセット周波数記憶手段、６  クラッタ中心周波数推定手段、７  出力信号記憶手段、８  フィルタ係数計算手段、１８  多重零点処理手段、９  フィルタ手段、１０  サーキュレータ、１１  フィルタ係数記憶手段、１２  フィルタ手段、１３  オフセット周波数選択手段、２３  多重零点割当選択手段、１４  オフセット周波数設定手段、１９  フィルタ次数選定手段、５０  多次フィルタ係数計算手段、５１  フィルタ手段、６２  フィルタ係数計算手段、６３  フィルタ手段。
【図１】
0003932957.tif 000023
【図２】
0003932957.tif 000024
【図３】
0003932957.tif 000025
【図４】
0003932957.tif 000026
【図５】
0003932957.tif 000027
【図６】
0003932957.tif 000028
【図７】
0003932957.tif 000029
【図８】
0003932957.tif 000030
【図９】
0003932957.tif 000031
【図１０】
0003932957.tif 000032
【図１１】
0003932957.tif 000033
【図１２】
0003932957.tif 000034
【図１３】
0003932957.tif 000035
【図１４】
0003932957.tif 000036
【図１５】
0003932957.tif 000037
【図１６】
0003932957.tif 000038
【図１７】
0003932957.tif 000039
【図１８】
0003932957.tif 000040
【図１９】
0003932957.tif 000041
【図２０】
0003932957.tif 000042
【図２１】
0003932957.tif 000043
【図２２】
0003932957.tif 000044
【図２３】
0003932957.tif 000045
【図２４】
0003932957.tif 000046
【図２５】
0003932957.tif 000047
【図２６】
0003932957.tif 000048
【図２７】
0003932957.tif 000049
【図２８】
0003932957.tif 000050