JPH10227851A - クラッタ抑圧装置 - Google Patents
クラッタ抑圧装置Info
- Publication number
- JPH10227851A JPH10227851A JP9033926A JP3392697A JPH10227851A JP H10227851 A JPH10227851 A JP H10227851A JP 9033926 A JP9033926 A JP 9033926A JP 3392697 A JP3392697 A JP 3392697A JP H10227851 A JPH10227851 A JP H10227851A
- Authority
- JP
- Japan
- Prior art keywords
- clutter
- filter coefficient
- notch filter
- level
- bandwidth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000001629 suppression Effects 0.000 title abstract description 32
- 238000012545 processing Methods 0.000 claims abstract description 40
- 230000005540 biological transmission Effects 0.000 claims abstract description 20
- 238000000034 method Methods 0.000 claims abstract description 14
- 238000001228 spectrum Methods 0.000 claims description 36
- 238000001914 filtration Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 27
- 238000010586 diagram Methods 0.000 description 17
- 239000013256 coordination polymer Substances 0.000 description 9
- 230000000694 effects Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 2
- 230000008034 disappearance Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001094 effect on targets Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
(57)【要約】
【課題】 ドップラ周波数がゼロでないクラッタでもそ
の影響を抑圧し、また、目標のドップラ周波数がクラッ
タのドップラ中心周波数に近い場合でも、クラッタ成分
のみを抑圧することを目的とする。 【解決手段】 レーダ送信周波数よりクラッタ中心周波
数を推定するクラッタ中心周波数算出手段と、この推定
されたクラッタ中心周波数と、予め設定したクラッタレ
ベルとクラッタ帯域幅に対応したノッチフィルタ係数と
からフィルタの係数を算出するノッチフィルタ係数発生
手段とを備え、受信した入力信号に上記フィルタ係数を
用いたノッチ・フィルタ処理によりクラッタを除去する
ようにした。
の影響を抑圧し、また、目標のドップラ周波数がクラッ
タのドップラ中心周波数に近い場合でも、クラッタ成分
のみを抑圧することを目的とする。 【解決手段】 レーダ送信周波数よりクラッタ中心周波
数を推定するクラッタ中心周波数算出手段と、この推定
されたクラッタ中心周波数と、予め設定したクラッタレ
ベルとクラッタ帯域幅に対応したノッチフィルタ係数と
からフィルタの係数を算出するノッチフィルタ係数発生
手段とを備え、受信した入力信号に上記フィルタ係数を
用いたノッチ・フィルタ処理によりクラッタを除去する
ようにした。
Description
【０００１】
【発明の属する技術分野】この発明は、短波帯のレーダ
の信号処理におけるシークラッタのクラッタ抑圧処理に
関するものである。
の信号処理におけるシークラッタのクラッタ抑圧処理に
関するものである。
【０００２】
【従来の技術】レーダにおいて地面等からの不要反射信
号であるクラッタは、目標の検出に対し大きな障害とな
る。クラッタの抑圧に対しては、クラッタが静止物から
の反射であり、クラッタのドップラ周波数がゼロである
ことが多いことから、反射信号のドップラ周波数で航空
機等の移動目標と弁別、抑圧する処理が一般的である。
図８は、例えばレーダ技術（吉田考監修、社団法人電子
通信学会編）に記載の単一消去器を用いた従来技術のク
ラッタ抑圧装置を用いた移動目標指示装置（Ｍｏｖｉｎ
ｇ Ｔａｒｇｅｔ Ｉｎｄｉｃａｔｏｒ；以下ＭＴＩと
いう）の機能系統図である。図において、１は位相検波
された入力信号である位相検波ビデオ信号、２は１サン
プル毎に位相検波ビデオ信号１を遅延させる遅延回路
（通常、遅延時間はレーダのパルス繰り返し周期に等し
い）、３は現在の入力信号と１サンプル前の位相検波ビ
デオ信号１との差をとる減算回路、４は位相検波ビデオ
信号１からクラッタを除去した出力信号である。
号であるクラッタは、目標の検出に対し大きな障害とな
る。クラッタの抑圧に対しては、クラッタが静止物から
の反射であり、クラッタのドップラ周波数がゼロである
ことが多いことから、反射信号のドップラ周波数で航空
機等の移動目標と弁別、抑圧する処理が一般的である。
図８は、例えばレーダ技術（吉田考監修、社団法人電子
通信学会編）に記載の単一消去器を用いた従来技術のク
ラッタ抑圧装置を用いた移動目標指示装置（Ｍｏｖｉｎ
ｇ Ｔａｒｇｅｔ Ｉｎｄｉｃａｔｏｒ；以下ＭＴＩと
いう）の機能系統図である。図において、１は位相検波
された入力信号である位相検波ビデオ信号、２は１サン
プル毎に位相検波ビデオ信号１を遅延させる遅延回路
（通常、遅延時間はレーダのパルス繰り返し周期に等し
い）、３は現在の入力信号と１サンプル前の位相検波ビ
デオ信号１との差をとる減算回路、４は位相検波ビデオ
信号１からクラッタを除去した出力信号である。
【０００３】従来の上記のように構成されたクラッタ抑
圧装置の動作を説明する。位相検波ビデオ信号１は、減
算回路３において、１サンプル前の位相検波ビデオ信号
１が減算される。電波が移動目標から反射してくると
き、その反射波の周波数ｆはドップラ効果によって、次
式（１）だけ偏移する。
圧装置の動作を説明する。位相検波ビデオ信号１は、減
算回路３において、１サンプル前の位相検波ビデオ信号
１が減算される。電波が移動目標から反射してくると
き、その反射波の周波数ｆはドップラ効果によって、次
式（１）だけ偏移する。
【０００４】
【数１】
【０００５】ここにｃは光速度、ｖは目標の電波進行方
向の分速度、ｆc は送信周波数を表す。一方、固定目標
からの反射波（クラッタ）は、周波数偏移を受けない。
従ってドップラ効果による周波数偏移、あるいはこれに
対応する位相偏移を検出すれば、移動目標だけを検出す
ることができる。パルスレーダにおいては、反射波を位
相検波し、１周期だけの遅延を与え、次の周期の検波出
力と逆極性にして加え合わせると固定目標からのクラッ
タは抑圧されることになる。簡単のため、入力の位相検
波ビデオ信号１を、ドップラ各周波数をωd ＝２πｆd
として Ｅ1 ＝Ａｓｉｎωd ｔ （２） と表すと、１周期前の遅延回路２の出力は、遅延時間を
Ｔとして Ｅ2 ＝Ａｓｉｎωd （ｔ−Ｔ） （３） となる。これをベクトル表示にすれば次式（４）とな
る。従って、減算回路２の出力信号４は式（５）とな
る。
向の分速度、ｆc は送信周波数を表す。一方、固定目標
からの反射波（クラッタ）は、周波数偏移を受けない。
従ってドップラ効果による周波数偏移、あるいはこれに
対応する位相偏移を検出すれば、移動目標だけを検出す
ることができる。パルスレーダにおいては、反射波を位
相検波し、１周期だけの遅延を与え、次の周期の検波出
力と逆極性にして加え合わせると固定目標からのクラッ
タは抑圧されることになる。簡単のため、入力の位相検
波ビデオ信号１を、ドップラ各周波数をωd ＝２πｆd
として Ｅ1 ＝Ａｓｉｎωd ｔ （２） と表すと、１周期前の遅延回路２の出力は、遅延時間を
Ｔとして Ｅ2 ＝Ａｓｉｎωd （ｔ−Ｔ） （３） となる。これをベクトル表示にすれば次式（４）とな
る。従って、減算回路２の出力信号４は式（５）とな
る。
【０００６】
【数２】
【０００７】これは遅延系をｅｘｐ（−ｊωd ｔ）なる
伝送特性を持つ素子と考えることができる。いま、Ｚを
式（６）とすると、遅延回路２と減算回路３で構成され
る消去器の伝送特性は式（７）となり、その速度レスポ
ンスは式（８）で表され、この単一消去ＭＴＩの特性を
図９に示す。
伝送特性を持つ素子と考えることができる。いま、Ｚを
式（６）とすると、遅延回路２と減算回路３で構成され
る消去器の伝送特性は式（７）となり、その速度レスポ
ンスは式（８）で表され、この単一消去ＭＴＩの特性を
図９に示す。
【０００８】
【数３】
【０００９】
【発明が解決しようとする課題】従来のクラッタ抑圧装
置は、以上のように構成されているので、短波帯のレー
ダではクラッタのドップラ効果周波数がゼロでない、例
えばＢｒａｇｇ・ｌｉｎｅとよばれるシークラッタを抑
圧できないという課題があった。また、クラッタの電力
スペクトルの形状に対してフィルタのカットオフ特性が
緩慢なため、目標のドップラ周波数がクラッタのドップ
ラ中心周波数の近傍に存在する場合は、クラッタととも
に目標信号も抑圧されるという課題があった。
置は、以上のように構成されているので、短波帯のレー
ダではクラッタのドップラ効果周波数がゼロでない、例
えばＢｒａｇｇ・ｌｉｎｅとよばれるシークラッタを抑
圧できないという課題があった。また、クラッタの電力
スペクトルの形状に対してフィルタのカットオフ特性が
緩慢なため、目標のドップラ周波数がクラッタのドップ
ラ中心周波数の近傍に存在する場合は、クラッタととも
に目標信号も抑圧されるという課題があった。
【００１０】この発明は上記の課題を解消するためにな
されたもので、ドップラ周波数がゼロでないクラッタで
もその影響を抑圧し、また目標信号のドップラ周波数と
クラッタのドップラ中心周波数が近い場合でも、クラッ
タ成分だけ抑圧するクラッタ抑圧装置を得ることを目的
とする。
されたもので、ドップラ周波数がゼロでないクラッタで
もその影響を抑圧し、また目標信号のドップラ周波数と
クラッタのドップラ中心周波数が近い場合でも、クラッ
タ成分だけ抑圧するクラッタ抑圧装置を得ることを目的
とする。
【００１１】
【課題を解決するための手段】この発明に係るクラッタ
抑圧装置は、レーダ送信周波数よりクラッタ中心周波数
を推定するクラッタ中心周波数算出手段と、この推定さ
れたクラッタ中心周波数と、予め設定したクラッタレベ
ルとクラッタ帯域幅に対応したノッチフィルタ係数とか
らフィルタの係数を算出するノッチフィルタ係数発生手
段とを備え、受信した入力信号に上記フィルタ係数を用
いたノッチフィルタ処理によりクラッタを除去するよう
にした。
抑圧装置は、レーダ送信周波数よりクラッタ中心周波数
を推定するクラッタ中心周波数算出手段と、この推定さ
れたクラッタ中心周波数と、予め設定したクラッタレベ
ルとクラッタ帯域幅に対応したノッチフィルタ係数とか
らフィルタの係数を算出するノッチフィルタ係数発生手
段とを備え、受信した入力信号に上記フィルタ係数を用
いたノッチフィルタ処理によりクラッタを除去するよう
にした。
【００１２】また更に、入力信号の電力スペクトルと推
定されたクラッタ中心周波数からクラッタのレベルを算
出するクラッタレベル推定手段と、予め各種のクラッタ
・レベルに相当する減衰量を設定したフィルタ係数から
推定されたクラッタ・レベルに相当する減衰量のフィル
タ係数を選択するノッチフィルタ係数選択手段を付加
し、推定されたクラッタ・レベル対応のノッチフィルタ
係数を発生させてノッチフィルタ処理によりクラッタを
除去するようにした。
定されたクラッタ中心周波数からクラッタのレベルを算
出するクラッタレベル推定手段と、予め各種のクラッタ
・レベルに相当する減衰量を設定したフィルタ係数から
推定されたクラッタ・レベルに相当する減衰量のフィル
タ係数を選択するノッチフィルタ係数選択手段を付加
し、推定されたクラッタ・レベル対応のノッチフィルタ
係数を発生させてノッチフィルタ処理によりクラッタを
除去するようにした。
【００１３】また更に、レーダ送信周波数よりクラッタ
中心周波数を推定するクラッタ中心周波数算出手段に換
えて、入力信号の電力スペクトルとレーダ送信周波数か
らクラッタ中心周波数を推定するクラッタ電力中心周波
数算出手段とし、入力信号の電力スペクトルと推定され
たクラッタ電力中心周波数からクラッタのレベルを算出
するクラッタレベル推定手段と、予め各種のクラッタ・
レベルに相当する減衰量を設定したフィルタ係数から上
記推定されたクラッタ・レベルに相当する減衰量のフィ
ルタ係数を選択するノッチフィルタ係数選択手段を付加
し、この推定されたクラッタ・レベル対応のノッチフィ
ルタ係数を発生させてノッチフィルタ処理によりクラッ
タを除去するようにした。
中心周波数を推定するクラッタ中心周波数算出手段に換
えて、入力信号の電力スペクトルとレーダ送信周波数か
らクラッタ中心周波数を推定するクラッタ電力中心周波
数算出手段とし、入力信号の電力スペクトルと推定され
たクラッタ電力中心周波数からクラッタのレベルを算出
するクラッタレベル推定手段と、予め各種のクラッタ・
レベルに相当する減衰量を設定したフィルタ係数から上
記推定されたクラッタ・レベルに相当する減衰量のフィ
ルタ係数を選択するノッチフィルタ係数選択手段を付加
し、この推定されたクラッタ・レベル対応のノッチフィ
ルタ係数を発生させてノッチフィルタ処理によりクラッ
タを除去するようにした。
【００１４】また更に、各種のクラッタ・レベルに相当
する減衰量を設定したフィルタ係数から対応フィルタ係
数を選択するノッチフィルタ係数選択手段に換えて、各
種のクラッタ帯域幅に相当する阻止帯域幅を設定したフ
ィルタ係数から対応帯域幅のフィルタ係数を選択するノ
ッチフィルタ係数選択手段とし、クラッタレベル推定手
段に換えて、入力信号の電力スペクトルと上記推定され
たクラッタ電力中心周波数からクラッタの帯域幅を推定
するクラッタ帯域幅推定手段とし、この推定されたクラ
ッタの帯域幅対応のノッチフィルタ係数を発生させてノ
ッチフィルタ処理によりクラッタを除去するようにし
た。
する減衰量を設定したフィルタ係数から対応フィルタ係
数を選択するノッチフィルタ係数選択手段に換えて、各
種のクラッタ帯域幅に相当する阻止帯域幅を設定したフ
ィルタ係数から対応帯域幅のフィルタ係数を選択するノ
ッチフィルタ係数選択手段とし、クラッタレベル推定手
段に換えて、入力信号の電力スペクトルと上記推定され
たクラッタ電力中心周波数からクラッタの帯域幅を推定
するクラッタ帯域幅推定手段とし、この推定されたクラ
ッタの帯域幅対応のノッチフィルタ係数を発生させてノ
ッチフィルタ処理によりクラッタを除去するようにし
た。
【００１５】また更に、各種のクラッタ・レベルに相当
する減衰量を設定したフィルタ係数から対応フィルタ係
数を選択するノッチフィルタ係数選択手段に換えて、各
種のクラッタ・レベルに相当する減衰量とクラッタ帯域
幅に相当する阻止帯域幅を設定したフィルタ係数から対
応レベルと帯域幅のフィルタ係数を選択するノッチフィ
ルタ係数選択手段とし、クラッタレベル推定手段に換え
て、入力信号の電力スペクトルと上記推定されたクラッ
タ電力中心周波数からクラッタ・レベルとクラッタ帯域
幅を推定するクラッタ・レベル帯域幅推定手段とし、こ
の推定されたクラッタ・レベルとクラッタ帯域幅対応の
ノッチフィルタ係数を発生させてノッチフィルタ処理に
よりクラッタを除去するようにした。
する減衰量を設定したフィルタ係数から対応フィルタ係
数を選択するノッチフィルタ係数選択手段に換えて、各
種のクラッタ・レベルに相当する減衰量とクラッタ帯域
幅に相当する阻止帯域幅を設定したフィルタ係数から対
応レベルと帯域幅のフィルタ係数を選択するノッチフィ
ルタ係数選択手段とし、クラッタレベル推定手段に換え
て、入力信号の電力スペクトルと上記推定されたクラッ
タ電力中心周波数からクラッタ・レベルとクラッタ帯域
幅を推定するクラッタ・レベル帯域幅推定手段とし、こ
の推定されたクラッタ・レベルとクラッタ帯域幅対応の
ノッチフィルタ係数を発生させてノッチフィルタ処理に
よりクラッタを除去するようにした。
【００１６】また更に、入力信号の電力スペクトルから
クラッタ対雑音比（Ｃ／Ｎ比）を算出するＣ／Ｎ比算出
手段と、この算出されたＣ／Ｎ比によりフィルタ処理を
行うかまたは入力信号を出力するかを切り換えるスイッ
チを付加し、算出されたＣ／Ｎ比によりノッチフィルタ
処理によるクラッタ除去の採否をするようにした。
クラッタ対雑音比（Ｃ／Ｎ比）を算出するＣ／Ｎ比算出
手段と、この算出されたＣ／Ｎ比によりフィルタ処理を
行うかまたは入力信号を出力するかを切り換えるスイッ
チを付加し、算出されたＣ／Ｎ比によりノッチフィルタ
処理によるクラッタ除去の採否をするようにした。
【００１７】
実施の形態１．図１はこの発明の実施の形態１における
クラッタ抑圧装置のブロック図である。図において、５
はレーダ送信周波数、６はレーダ送信周波数情報からシ
ークラッタの中心周波数算出式によりクラッタの中心周
波数を算出するクラッタ中心周波数算出部（ａ）、７は
クラッタ中心周波数と予め設定したノッチフィルタ係数
とから、ノッチフィルタのノッチがクラッタの中心とな
るフィルタの係数を算出するノッチフィルタ係数発生部
（ａ）、８は位相検波ビデオ信号１に対し、算出された
フィルタ係数から構成されるフィルタによりクラッタを
除去するフィルタ処理部である。１、４の入力信号は従
来装置と同等である。
クラッタ抑圧装置のブロック図である。図において、５
はレーダ送信周波数、６はレーダ送信周波数情報からシ
ークラッタの中心周波数算出式によりクラッタの中心周
波数を算出するクラッタ中心周波数算出部（ａ）、７は
クラッタ中心周波数と予め設定したノッチフィルタ係数
とから、ノッチフィルタのノッチがクラッタの中心とな
るフィルタの係数を算出するノッチフィルタ係数発生部
（ａ）、８は位相検波ビデオ信号１に対し、算出された
フィルタ係数から構成されるフィルタによりクラッタを
除去するフィルタ処理部である。１、４の入力信号は従
来装置と同等である。
【００１８】上記の構成のクラッタ抑圧装置の動作を説
明する。まず、クラッタ中心周波数算出部（ａ）６にお
けるクラッタの中心周波数の算出方法を説明する。図２
は短波帯のレーダ装置において海面からの反射信号に含
まれるシークラッタのドップラスペクトルを説明する図
である。図においてドップラスペクトルの正負に顕著な
ピーク２３が見られる。これはＢｒａｇｇ・ｌｉｎｅと
呼ばれるシークラッタである。２つのピーク２３のドッ
プラ周波数ｆｄは、次式（９）であることが知られてい
る。 ±√｛（ｇ・ｆ）／（π・ｃ）｝ （９） ｆ：レーダの送信周波数 ｇ：重力加速度 π：円周率 ｃ：光速 従ってレーダの送信周波数５によりＢｒａｇｇ・ｌｉｎ
ｅのピーク２３のドップラ周波数は一意に決まる。クラ
ッタ中心周波数算出部（ａ）６はドップラ周波数ｆd を
クラッタの中心周波数ｆdcとして次式（１０）のように
出力する。ｆdc＝ｆd （１０）
明する。まず、クラッタ中心周波数算出部（ａ）６にお
けるクラッタの中心周波数の算出方法を説明する。図２
は短波帯のレーダ装置において海面からの反射信号に含
まれるシークラッタのドップラスペクトルを説明する図
である。図においてドップラスペクトルの正負に顕著な
ピーク２３が見られる。これはＢｒａｇｇ・ｌｉｎｅと
呼ばれるシークラッタである。２つのピーク２３のドッ
プラ周波数ｆｄは、次式（９）であることが知られてい
る。 ±√｛（ｇ・ｆ）／（π・ｃ）｝ （９） ｆ：レーダの送信周波数 ｇ：重力加速度 π：円周率 ｃ：光速 従ってレーダの送信周波数５によりＢｒａｇｇ・ｌｉｎ
ｅのピーク２３のドップラ周波数は一意に決まる。クラ
ッタ中心周波数算出部（ａ）６はドップラ周波数ｆd を
クラッタの中心周波数ｆdcとして次式（１０）のように
出力する。ｆdc＝ｆd （１０）
【００１９】ノッチフィルタ係数発生部（ａ）７は、入
力されるクラッタ中心周波数と、予め設定したノッチフ
ィルタ係数とからフィルタ処理部８に出力するノッチフ
ィルタ係数を算出する。まず、クラッタを抑圧するノッ
チフィルタの係数を予め設計しておく。例えば文献「Ｄ
ＯＮＡＬＤ Ｅ．ＢＡＲＲＩＣＫ，ＪＡＭＥＳ Ｍ．Ｈ
ＥＡＤＲＩＣＫ，ＲＯＢＥＲＴ Ｗ．ＢＯＧＬＥ，ＤＯ
ＵＧＬＡＳＳ Ｄ．ＣＲＯＭＢＩＥ（１９７７），Ｓｅ
ａ Ｂａｃｋｓｃａｔｔｅｒ ａｔ ＨＦ：Ｉｎｔｅｒ
ｐｒｅｔａｔｉｏｎ ａｎｄ Ｕｔｉｌｉｚａｔｉｏｎ
ｏｆ ｔｈｅＥｃｈｏ」からクラッタの形状はガウス
分布で近似でき、指定される減衰量及び阻止帯域幅より
ノッチフィルタの係数を設計できる。この設計したノッ
チフィルタ係数はノッチフィルタ係数発生部（ａ）７に
予め記憶しておく。
力されるクラッタ中心周波数と、予め設定したノッチフ
ィルタ係数とからフィルタ処理部８に出力するノッチフ
ィルタ係数を算出する。まず、クラッタを抑圧するノッ
チフィルタの係数を予め設計しておく。例えば文献「Ｄ
ＯＮＡＬＤ Ｅ．ＢＡＲＲＩＣＫ，ＪＡＭＥＳ Ｍ．Ｈ
ＥＡＤＲＩＣＫ，ＲＯＢＥＲＴ Ｗ．ＢＯＧＬＥ，ＤＯ
ＵＧＬＡＳＳ Ｄ．ＣＲＯＭＢＩＥ（１９７７），Ｓｅ
ａ Ｂａｃｋｓｃａｔｔｅｒ ａｔ ＨＦ：Ｉｎｔｅｒ
ｐｒｅｔａｔｉｏｎ ａｎｄ Ｕｔｉｌｉｚａｔｉｏｎ
ｏｆ ｔｈｅＥｃｈｏ」からクラッタの形状はガウス
分布で近似でき、指定される減衰量及び阻止帯域幅より
ノッチフィルタの係数を設計できる。この設計したノッ
チフィルタ係数はノッチフィルタ係数発生部（ａ）７に
予め記憶しておく。
【００２０】このノッチフィルタ係数はノッチの中心が
周波数０Ｈｚにあるため、そのまま係数をノッチフィル
タの荷重値として用いるとクラッタ中心周波数からはず
れてクラッタを抑圧できない。そこで予め設定したノッ
チフィルタ係数とクラッタ中心周波数から次式により、
ノッチフィルタのノッチがクラッタ中心周波数ｆdcとな
るようにノッチフィルタ係数を調整する。
周波数０Ｈｚにあるため、そのまま係数をノッチフィル
タの荷重値として用いるとクラッタ中心周波数からはず
れてクラッタを抑圧できない。そこで予め設定したノッ
チフィルタ係数とクラッタ中心周波数から次式により、
ノッチフィルタのノッチがクラッタ中心周波数ｆdcとな
るようにノッチフィルタ係数を調整する。
【００２１】
【数４】
【００２２】ここで、ｈ0 （ｉ）は予め設定したノッチ
フィルタ係数、Ｌはノッチフィルタの段数を表す。上記
のノッチフィルタ係数ｈ（ｉ）はフィルタ処理部８へ出
力される。フィルタ処理部８ではノッチフィルタ係数ｈ
（ｉ）を用いて次式のディジタルフィルタ処理を行い、
次のように位相検波ビデオ信号からクラッタ信号を除去
する。 Ｚ（ｉ）＝Ｈ（ｉ）・Ｓ（ｉ） （１２） ここで、Ｈ（ｉ）はノッチフィルタ係数ｈ（ｉ）を高速
フーリエ変換した結果、Ｓ（ｉ）は位相検波ビデオ信号
ｓ（ｉ）を高速フーリエ変換した結果を表す。フィルタ
処理部８はこのＺ（ｉ）を逆高速フーリエ変換し、クラ
ッタ信号だけを除去した信号として出力する。このよう
に従来技術では抑圧できなかったドップラ周波数がゼロ
でないクラッタを抑圧できる。また、単一消去器ＭＴＩ
ではクラッタの電力スペクトルの形状よりカットオフ特
性が緩慢なため、目標のドップラ周波数がクラッタの中
心周波数の近傍に存在する場合は目標信号も抑圧する可
能性があるが、ノッチフィルタではクラッタの電力スペ
クトルの形状にあったフィルタを形成できるため、目標
信号の抑圧を極力少なくしてクラッタを抑圧できる。
フィルタ係数、Ｌはノッチフィルタの段数を表す。上記
のノッチフィルタ係数ｈ（ｉ）はフィルタ処理部８へ出
力される。フィルタ処理部８ではノッチフィルタ係数ｈ
（ｉ）を用いて次式のディジタルフィルタ処理を行い、
次のように位相検波ビデオ信号からクラッタ信号を除去
する。 Ｚ（ｉ）＝Ｈ（ｉ）・Ｓ（ｉ） （１２） ここで、Ｈ（ｉ）はノッチフィルタ係数ｈ（ｉ）を高速
フーリエ変換した結果、Ｓ（ｉ）は位相検波ビデオ信号
ｓ（ｉ）を高速フーリエ変換した結果を表す。フィルタ
処理部８はこのＺ（ｉ）を逆高速フーリエ変換し、クラ
ッタ信号だけを除去した信号として出力する。このよう
に従来技術では抑圧できなかったドップラ周波数がゼロ
でないクラッタを抑圧できる。また、単一消去器ＭＴＩ
ではクラッタの電力スペクトルの形状よりカットオフ特
性が緩慢なため、目標のドップラ周波数がクラッタの中
心周波数の近傍に存在する場合は目標信号も抑圧する可
能性があるが、ノッチフィルタではクラッタの電力スペ
クトルの形状にあったフィルタを形成できるため、目標
信号の抑圧を極力少なくしてクラッタを抑圧できる。
【００２３】実施の形態２．次に、クラッタのレベルの
変動によるクラッタの消え残りや過剰抑圧を押さえたク
ラッタ抑圧装置を説明する。図３は、この発明の本実施
の形態におけるクラッタ抑圧装置のブロック図である。
図において、図１と同等要素は同一符号を付してその説
明を省略する。新規な要素として、９は位相検波ビデオ
信号１の電力スペクトルを算出する電力スペクトル算出
部、１０は電力スペクトルからクラッタのレベルを算出
するクラッタレベル推定部、１１はクラッタのレベルに
相当する減衰量となるようなフィルタ係数を出力するノ
ッチフィルタ係数選択部（ａ）である。１２はクラッタ
中心周波数とフィルタ係数とからノッチフィルタのノッ
チがクラッタの中心となるフィルタ係数を算出するノッ
チフィルタ係数発生部（ｂ）である。
変動によるクラッタの消え残りや過剰抑圧を押さえたク
ラッタ抑圧装置を説明する。図３は、この発明の本実施
の形態におけるクラッタ抑圧装置のブロック図である。
図において、図１と同等要素は同一符号を付してその説
明を省略する。新規な要素として、９は位相検波ビデオ
信号１の電力スペクトルを算出する電力スペクトル算出
部、１０は電力スペクトルからクラッタのレベルを算出
するクラッタレベル推定部、１１はクラッタのレベルに
相当する減衰量となるようなフィルタ係数を出力するノ
ッチフィルタ係数選択部（ａ）である。１２はクラッタ
中心周波数とフィルタ係数とからノッチフィルタのノッ
チがクラッタの中心となるフィルタ係数を算出するノッ
チフィルタ係数発生部（ｂ）である。
【００２４】上記構成の装置の動作を説明する。図３に
おいて、電力スペクトル算出部９は位相検波ビデオ信号
１から高速フーリエ変換等の手法により入力信号の電力
スペクトルを算出する。クラッタレベル推定部１０はク
ラッタ中心周波数算出部（ａ）６より算出されたクラッ
タの中心周波数における電力スペクトルの値をクラッタ
のレベルとする。クラッタレベルは海面の状況により変
動するため実施の形態１で予め推定した減衰量ではクラ
ッタの過剰抑圧または、クラッタの消え残りを生じるお
それがある。フィルタ係数選択部（ａ）１１はクラッタ
レベルと予め設定した複数の減衰量とを比較し、クラッ
タレベルに相当する減衰量を選択してその減衰量に対応
したフィルタ係数をノッチフィルタ係数発生部（ｂ）１
２へ出力する。減衰量と減衰量に対応したフィルタ係数
はノッチフィルタ係数選択部（ａ）１１に予め記憶して
おく。
おいて、電力スペクトル算出部９は位相検波ビデオ信号
１から高速フーリエ変換等の手法により入力信号の電力
スペクトルを算出する。クラッタレベル推定部１０はク
ラッタ中心周波数算出部（ａ）６より算出されたクラッ
タの中心周波数における電力スペクトルの値をクラッタ
のレベルとする。クラッタレベルは海面の状況により変
動するため実施の形態１で予め推定した減衰量ではクラ
ッタの過剰抑圧または、クラッタの消え残りを生じるお
それがある。フィルタ係数選択部（ａ）１１はクラッタ
レベルと予め設定した複数の減衰量とを比較し、クラッ
タレベルに相当する減衰量を選択してその減衰量に対応
したフィルタ係数をノッチフィルタ係数発生部（ｂ）１
２へ出力する。減衰量と減衰量に対応したフィルタ係数
はノッチフィルタ係数選択部（ａ）１１に予め記憶して
おく。
【００２５】クラッタレベルに相当する減衰量を選択す
るのは、具体的には以下のように行う。次式（１３）よ
りクラッタレベルＣＰc と予め設定した減衰量ＣＰ
（ｊ）の差を求め、クラッタレベルを下回らない最小の
減衰量を選択する。 ＣＰdef ＝ＣＰc −ＣＰ（ｊ） （１３） Ｊ＝１，２，・・・，ＪＪ ここで、ＣＰc はクラッタレベル、ＣＰ（ｊ）は予め設
定した減衰量、ＪＪは予め設定した減衰量の数を表す。
ＣＰdef が最小となるＣＰ（ｊ）に対応するノッチフィ
ルタ係数をフィルタ処理部８へ出力する。
るのは、具体的には以下のように行う。次式（１３）よ
りクラッタレベルＣＰc と予め設定した減衰量ＣＰ
（ｊ）の差を求め、クラッタレベルを下回らない最小の
減衰量を選択する。 ＣＰdef ＝ＣＰc −ＣＰ（ｊ） （１３） Ｊ＝１，２，・・・，ＪＪ ここで、ＣＰc はクラッタレベル、ＣＰ（ｊ）は予め設
定した減衰量、ＪＪは予め設定した減衰量の数を表す。
ＣＰdef が最小となるＣＰ（ｊ）に対応するノッチフィ
ルタ係数をフィルタ処理部８へ出力する。
【００２６】ノッチフィルタ係数発生部（ｂ）１２はノ
ッチフィルタ係数選択部（ａ）１１から出力されたノッ
チフィルタ係数とクラッタ中心周波数とからノッチフィ
ルタのノッチがクラッタ中心周波数となるように調整す
る。ノッチフィルタ係数の算出方法は実施の形態１と同
様であるため動作の詳細記述を省略する。図３に示す装
置では、クラッタレベルが変動してもクラッタレベルに
あわせてクラッタ抑圧を変動させるので、消え残りや過
剰抑圧を押さえることができる。
ッチフィルタ係数選択部（ａ）１１から出力されたノッ
チフィルタ係数とクラッタ中心周波数とからノッチフィ
ルタのノッチがクラッタ中心周波数となるように調整す
る。ノッチフィルタ係数の算出方法は実施の形態１と同
様であるため動作の詳細記述を省略する。図３に示す装
置では、クラッタレベルが変動してもクラッタレベルに
あわせてクラッタ抑圧を変動させるので、消え残りや過
剰抑圧を押さえることができる。
【００２７】実施の形態３．上記実施の形態では、レー
ダの送信周波数から求まるドップラ周波数をクラッタの
中心周波数としている。しかし、実際のＢｒａｇｇ・ｌ
ｉｎｅのドップラ周波数は海面の状態によりレーダ送信
周波数から求まるクラッタ中心周波数からずれる。本実
施の形態ではこのずれを補正する場合について説明す
る。図４はこの発明の本実施の形態におけるクラッタ抑
圧装置のブロック図である。図において、先の実施の形
態での要素と同じ要素は同一符号を付してその説明を省
略する。新規な要素として、１３はレーダ送信周波数に
基いて得られたシークラッタの中心周波数と、電力スペ
クトルとからクラッタの中心周波数を算出するクラッタ
中心周波数算出部（ｂ）である。
ダの送信周波数から求まるドップラ周波数をクラッタの
中心周波数としている。しかし、実際のＢｒａｇｇ・ｌ
ｉｎｅのドップラ周波数は海面の状態によりレーダ送信
周波数から求まるクラッタ中心周波数からずれる。本実
施の形態ではこのずれを補正する場合について説明す
る。図４はこの発明の本実施の形態におけるクラッタ抑
圧装置のブロック図である。図において、先の実施の形
態での要素と同じ要素は同一符号を付してその説明を省
略する。新規な要素として、１３はレーダ送信周波数に
基いて得られたシークラッタの中心周波数と、電力スペ
クトルとからクラッタの中心周波数を算出するクラッタ
中心周波数算出部（ｂ）である。
【００２８】上記構成の装置の動作を説明する。クラッ
タ中心周波数算出部（ｂ）１３は位相検波ビデオの電力
スペクトルからクラッタ中心周波数を算出する。まずレ
ーダの送信周波数から求まるＢｒａｇｇ・ｌｉｎｅのド
ップラ周波数ｆd を中心として予め設定したドップラ範
囲内±ＩＦＧＡＴＥで電力の最大値を検出し、その周波
数をクラッタピーク周波数ｆpcとする。最大値は次式
（１４）より算出する。 Ｐc ＝Ｍａｘ（｜ｓ（ｆ）｜2 ， ｆ＝ｆd −ＩＦＧＡＴＥ〜ｆd ＋ＩＦＧＡＴＥ） （１４） ここで、Ｍａｘ（）は（）内の数式の最大値をとる関
数、ｓ（ｆ）は位相検波ビデオ信号の電力スペクトル、
ｆd は計算式より算出したクラッタのドップラ周波数、
ＩＦＧＡＴＥは予め設定したドップラ範囲を表す。電力
がＰc となる周波数をｆpcとする。
タ中心周波数算出部（ｂ）１３は位相検波ビデオの電力
スペクトルからクラッタ中心周波数を算出する。まずレ
ーダの送信周波数から求まるＢｒａｇｇ・ｌｉｎｅのド
ップラ周波数ｆd を中心として予め設定したドップラ範
囲内±ＩＦＧＡＴＥで電力の最大値を検出し、その周波
数をクラッタピーク周波数ｆpcとする。最大値は次式
（１４）より算出する。 Ｐc ＝Ｍａｘ（｜ｓ（ｆ）｜2 ， ｆ＝ｆd −ＩＦＧＡＴＥ〜ｆd ＋ＩＦＧＡＴＥ） （１４） ここで、Ｍａｘ（）は（）内の数式の最大値をとる関
数、ｓ（ｆ）は位相検波ビデオ信号の電力スペクトル、
ｆd は計算式より算出したクラッタのドップラ周波数、
ＩＦＧＡＴＥは予め設定したドップラ範囲を表す。電力
がＰc となる周波数をｆpcとする。
【００２９】上記で求めたｆpcを中心に、別に設定した
モーメント処理ゲート内でモーメント法を適用し、クラ
ッタ中心周波数ｆdcを算出する。まず１次モーメントＭ
1 は、次式（１５）で表される。ここで、ＭＯＭＷはモ
ーメント処理ゲートを表す。これを用いて式（１６）が
得られ、クラッタ中心周波数ｆdcが算出される。
モーメント処理ゲート内でモーメント法を適用し、クラ
ッタ中心周波数ｆdcを算出する。まず１次モーメントＭ
1 は、次式（１５）で表される。ここで、ＭＯＭＷはモ
ーメント処理ゲートを表す。これを用いて式（１６）が
得られ、クラッタ中心周波数ｆdcが算出される。
【００３０】
【数５】
【００３１】図４ではクラッタの中心周波数が変動する
場合でも、ノッチフィルタのノッチを常にクラッタの中
心周波数に形成することができる。
場合でも、ノッチフィルタのノッチを常にクラッタの中
心周波数に形成することができる。
【００３２】次に、クラッタの帯域幅の変動によるクラ
ッタの消え残りや過剰抑圧を押さえたクラッタ抑圧装置
を説明する。図５は本実施の形態における他のクラッタ
抑圧装置のブロック図である。図において１４は電力ス
ペクトルからクラッタの帯域幅を算出するクラッタ推定
部、１５はクラッタの帯域幅に相当する阻止帯域幅とな
るようなフィルタを係数出力するノッチフィルタ係数選
択部（ｂ）である。他の要素は上記実施の形態の対応番
号の要素と同じものである。
ッタの消え残りや過剰抑圧を押さえたクラッタ抑圧装置
を説明する。図５は本実施の形態における他のクラッタ
抑圧装置のブロック図である。図において１４は電力ス
ペクトルからクラッタの帯域幅を算出するクラッタ推定
部、１５はクラッタの帯域幅に相当する阻止帯域幅とな
るようなフィルタを係数出力するノッチフィルタ係数選
択部（ｂ）である。他の要素は上記実施の形態の対応番
号の要素と同じものである。
【００３３】図５において、クラッタ帯域幅推定部１４
ではクラッタ中心周波数と位相検波ビデオ信号の電力ス
ペクトルからクラッタの帯域幅を推定し、ノッチフィル
タ係数選択部（ｂ）１５に出力する。クラッタ帯域幅の
推定はクラッタの形状がガウス分布で近似できることを
利用する。クラッタのピーク周波数ｆpcを中心に、別に
設定したモーメント処理ゲート内でモーメント法を適用
し、クラッタの帯域幅ＢＷc を算出する。まず１次モー
メントＭ1 は、次式（１７）で表される。ここで、ｓ
（ｆ）は位相検波ビデオの電力スペクトル、ＭＯＭＷは
モーメント処理ゲートを表す。クラッタの帯域幅は２次
モーメントを用いて求める。２次モーメントを用いて式
（１８）が得られる。σはクラッタ分布の標準偏差を示
すから、クラッタの帯域幅ＢＷc とすると、これは式
（１９）となり、クラッタの帯域幅が推定される。
ではクラッタ中心周波数と位相検波ビデオ信号の電力ス
ペクトルからクラッタの帯域幅を推定し、ノッチフィル
タ係数選択部（ｂ）１５に出力する。クラッタ帯域幅の
推定はクラッタの形状がガウス分布で近似できることを
利用する。クラッタのピーク周波数ｆpcを中心に、別に
設定したモーメント処理ゲート内でモーメント法を適用
し、クラッタの帯域幅ＢＷc を算出する。まず１次モー
メントＭ1 は、次式（１７）で表される。ここで、ｓ
（ｆ）は位相検波ビデオの電力スペクトル、ＭＯＭＷは
モーメント処理ゲートを表す。クラッタの帯域幅は２次
モーメントを用いて求める。２次モーメントを用いて式
（１８）が得られる。σはクラッタ分布の標準偏差を示
すから、クラッタの帯域幅ＢＷc とすると、これは式
（１９）となり、クラッタの帯域幅が推定される。
【００３４】
【数６】
【００３５】ノッチフィルタ係数選択部（ｂ）１５は電
力スペクトルより推定したクラッタの帯域幅ＢＷc と予
め設定した複数の阻止帯域幅とを比較し、クラッタの帯
域幅に相当する阻止帯域幅を選択してその阻止帯域幅に
対応したフィルタ係数をノッチフィルタ係数発生部
（ｂ）１２へ出力する。阻止帯域幅と阻止帯域幅に対応
したフィルタ係数はノッチフィルタ係数選択部（ｂ）１
５に予め記憶しておく。
力スペクトルより推定したクラッタの帯域幅ＢＷc と予
め設定した複数の阻止帯域幅とを比較し、クラッタの帯
域幅に相当する阻止帯域幅を選択してその阻止帯域幅に
対応したフィルタ係数をノッチフィルタ係数発生部
（ｂ）１２へ出力する。阻止帯域幅と阻止帯域幅に対応
したフィルタ係数はノッチフィルタ係数選択部（ｂ）１
５に予め記憶しておく。
【００３６】クラッタの帯域幅に相当する阻止帯域幅を
選択するのは、具体的には以下のように行う。次式（２
０）よりクラッタ帯域幅ＢＷc と予め設定した阻止帯域
幅ＢＷ（ｊ）の差を求め、クラッタの帯域幅を下回らな
い最小の阻止帯域幅を選択する。 ＢＷdef ＝ＢＷc −ＢＷ（ｊ） （２０） Ｊ＝１，２，・・・，ＪＪ ここで、ＢＷc はクラッタの帯域幅、ＢＷ（ｊ）は予め
設定した阻止帯域幅、ＪＪは予め設定した阻止帯域幅の
数を表す。ＢＷdef が最小となるＢＷ（ｊ）に対応する
ノッチフィルタ係数をフィルタ処理部８へ出力する。図
５に示す装置では、クラッタの帯域幅は変動しても、ク
ラッタの帯域幅にあわせてクラッタ抑圧を変動させるの
で、クラッタの消え残りや過剰抑圧を押さえることがで
きる。
選択するのは、具体的には以下のように行う。次式（２
０）よりクラッタ帯域幅ＢＷc と予め設定した阻止帯域
幅ＢＷ（ｊ）の差を求め、クラッタの帯域幅を下回らな
い最小の阻止帯域幅を選択する。 ＢＷdef ＝ＢＷc −ＢＷ（ｊ） （２０） Ｊ＝１，２，・・・，ＪＪ ここで、ＢＷc はクラッタの帯域幅、ＢＷ（ｊ）は予め
設定した阻止帯域幅、ＪＪは予め設定した阻止帯域幅の
数を表す。ＢＷdef が最小となるＢＷ（ｊ）に対応する
ノッチフィルタ係数をフィルタ処理部８へ出力する。図
５に示す装置では、クラッタの帯域幅は変動しても、ク
ラッタの帯域幅にあわせてクラッタ抑圧を変動させるの
で、クラッタの消え残りや過剰抑圧を押さえることがで
きる。
【００３７】次に、クラッタのレベルと帯域幅の変動に
よるクラッタの消え残りや過剰抑圧を押さえたクラッタ
抑圧装置を説明する。図６は本実施の形態における他の
クラッタ抑圧装置のブロック図である。図において、１
６は電力スペクトルからクラッタのレベルと帯域幅を算
出するクラッタレベル帯域幅推定部、１７はクラッタの
レベルに相当する減衰量とクラッタの帯域幅に相当する
阻止帯域幅に対応したフィルタ係数を出力するノッチフ
ィルタ係数選択部（ｃ）である。他の要素は上記実施の
形態の対応番号の要素と同じものである。
よるクラッタの消え残りや過剰抑圧を押さえたクラッタ
抑圧装置を説明する。図６は本実施の形態における他の
クラッタ抑圧装置のブロック図である。図において、１
６は電力スペクトルからクラッタのレベルと帯域幅を算
出するクラッタレベル帯域幅推定部、１７はクラッタの
レベルに相当する減衰量とクラッタの帯域幅に相当する
阻止帯域幅に対応したフィルタ係数を出力するノッチフ
ィルタ係数選択部（ｃ）である。他の要素は上記実施の
形態の対応番号の要素と同じものである。
【００３８】図６において、ノッチフィルタ係数選択部
（ｃ）１７は電力スペクトルより推定したクラッタレベ
ル、帯域幅と、予め設定した複数の減衰量、阻止帯域幅
とを比較し、クラッタレベル、帯域幅に相当する減衰
量、阻止帯域幅に対応したフィルタ係数をノッチフィル
タ係数発生部（ｂ）１２へ出力する。減衰量、阻止帯域
幅に対応したフィルタ係数はノッチフィルタ係数選択部
（ｃ）１７に予め記憶しておく。図６に示す装置では、
クラッタのレベル相当の減衰量と、クラッタの帯域幅相
当の阻止帯域幅となるノッチフィルタを形成するので、
クラッタのレベルと帯域幅が変動しても常に追従し、ク
ラッタの消え残りや過剰抑圧を押さえることができる。
（ｃ）１７は電力スペクトルより推定したクラッタレベ
ル、帯域幅と、予め設定した複数の減衰量、阻止帯域幅
とを比較し、クラッタレベル、帯域幅に相当する減衰
量、阻止帯域幅に対応したフィルタ係数をノッチフィル
タ係数発生部（ｂ）１２へ出力する。減衰量、阻止帯域
幅に対応したフィルタ係数はノッチフィルタ係数選択部
（ｃ）１７に予め記憶しておく。図６に示す装置では、
クラッタのレベル相当の減衰量と、クラッタの帯域幅相
当の阻止帯域幅となるノッチフィルタを形成するので、
クラッタのレベルと帯域幅が変動しても常に追従し、ク
ラッタの消え残りや過剰抑圧を押さえることができる。
【００３９】実施の形態４．上記実施の形態では、つね
にクラッタ抑圧処理を行っている。このため、クラッタ
抑圧後の目標検出処理において誤警報を誘発する恐れの
ない小クラッタも抑圧する。一方でクラッタ抑圧処理に
より信号も減衰する。そこでクラッタ抑圧判定処理で抑
圧した方がよいかどうかを判断する回路を加えること
で、抑圧の必要がない小クラッタの過剰抑圧を押さえ、
クラッタの近傍または、クラッタ内に存在する目標信号
の減衰を防ぐことができる。図６に示す装置にクラッタ
抑圧判定処理を加えた場合について説明する。なお、図
１、３〜５に示す装置にクラッタ抑圧判定処理を加えた
場合も同様な効果が得られる。図７はこの発明の本実施
の形態におけるクラッタ抑圧装置のブロック図である。
新規な要素として、１８は位相検波ビデオ信号１の電力
スペクトルより雑音レベルを算出する雑音レベル推定
部、１９は算出されたクラッタレベルと雑音レベルとか
らＣ／Ｎ比を算出するＣ／Ｎ比算出部、２０はクラッタ
抑圧処理の入・切を判定する抑圧処理判定部、２１は抑
圧処理判定部２０の判定結果よりクラッタ抑圧の入り切
りを行うスイッチ、２２はスイッチ２１により出力され
る位相検波ビデオ信号１またはフィルタ処理部８により
位相検波ビデオ信号１からクラッタが除去された出力信
号である。その他の要素は先の実施の形態における同番
号の要素と同じものである。
にクラッタ抑圧処理を行っている。このため、クラッタ
抑圧後の目標検出処理において誤警報を誘発する恐れの
ない小クラッタも抑圧する。一方でクラッタ抑圧処理に
より信号も減衰する。そこでクラッタ抑圧判定処理で抑
圧した方がよいかどうかを判断する回路を加えること
で、抑圧の必要がない小クラッタの過剰抑圧を押さえ、
クラッタの近傍または、クラッタ内に存在する目標信号
の減衰を防ぐことができる。図６に示す装置にクラッタ
抑圧判定処理を加えた場合について説明する。なお、図
１、３〜５に示す装置にクラッタ抑圧判定処理を加えた
場合も同様な効果が得られる。図７はこの発明の本実施
の形態におけるクラッタ抑圧装置のブロック図である。
新規な要素として、１８は位相検波ビデオ信号１の電力
スペクトルより雑音レベルを算出する雑音レベル推定
部、１９は算出されたクラッタレベルと雑音レベルとか
らＣ／Ｎ比を算出するＣ／Ｎ比算出部、２０はクラッタ
抑圧処理の入・切を判定する抑圧処理判定部、２１は抑
圧処理判定部２０の判定結果よりクラッタ抑圧の入り切
りを行うスイッチ、２２はスイッチ２１により出力され
る位相検波ビデオ信号１またはフィルタ処理部８により
位相検波ビデオ信号１からクラッタが除去された出力信
号である。その他の要素は先の実施の形態における同番
号の要素と同じものである。
【００４０】図７において、抑圧処理判定部２０は算出
されたＣ／Ｎ比がユーザにより設定されたしきい値より
高い場合にはクラッタ抑圧”入”の判定を行い、低い場
合にクラッタ抑圧”切”の判定を行う。スイッチ２１は
抑圧処理判定部２０の結果より、クラッタ抑圧入り判定
の場合は位相検波ビデオ信号１からクラッタが除去され
た信号を出力し、クラッタ抑圧切り判定の場合は位相検
波ビデオ信号１を出力する。図７に示す装置では、Ｃ／
Ｎ比が一定値以下の場合はクラッタを抑圧せず、従って
目標検出に影響の少ない小クラッタは抑圧しないので、
過剰抑圧を押さえ、小クラッタの近傍または、小クラッ
タ内の目標はそのまま検出される。
されたＣ／Ｎ比がユーザにより設定されたしきい値より
高い場合にはクラッタ抑圧”入”の判定を行い、低い場
合にクラッタ抑圧”切”の判定を行う。スイッチ２１は
抑圧処理判定部２０の結果より、クラッタ抑圧入り判定
の場合は位相検波ビデオ信号１からクラッタが除去され
た信号を出力し、クラッタ抑圧切り判定の場合は位相検
波ビデオ信号１を出力する。図７に示す装置では、Ｃ／
Ｎ比が一定値以下の場合はクラッタを抑圧せず、従って
目標検出に影響の少ない小クラッタは抑圧しないので、
過剰抑圧を押さえ、小クラッタの近傍または、小クラッ
タ内の目標はそのまま検出される。
【００４１】
【発明の効果】以上のようにこの発明によれば、クラッ
タ中心周波数算出手段と、フィルタの係数を算出するノ
ッチフィルタ係数発生手段とを備えて受信した入力信号
にノッチフィルタ処理をするので、クラッタの近くの周
波数の目標からの信号を減衰させないで識別できる効果
がある。
タ中心周波数算出手段と、フィルタの係数を算出するノ
ッチフィルタ係数発生手段とを備えて受信した入力信号
にノッチフィルタ処理をするので、クラッタの近くの周
波数の目標からの信号を減衰させないで識別できる効果
がある。
【００４２】また更に、クラッタレベル推定手段と、推
定クラッタ・レベル相当の減衰量のフィルタ係数を選択
するノッチフィルタ係数選択手段を付加したので、クラ
ッタ・レベルが変動する場合でも必要な目標信号が得ら
れる効果がある。
定クラッタ・レベル相当の減衰量のフィルタ係数を選択
するノッチフィルタ係数選択手段を付加したので、クラ
ッタ・レベルが変動する場合でも必要な目標信号が得ら
れる効果がある。
【００４３】また更に、レーダ送信周波数でクラッタ中
心周波数を推定することに換え、入力信号電力スペクト
ルとレーダ送信周波数を用いて推定するクラッタ電力中
心周波数算出手段とし、クラッタレベル推定手段で推定
されたクラッタ・レベルに相当する減衰量のフィルタ係
数を選択するノッチフィルタ係数選択手段を付加したの
で、クラッタの中心周波数が変動する場合でも必要な目
標信号が得られる効果がある。
心周波数を推定することに換え、入力信号電力スペクト
ルとレーダ送信周波数を用いて推定するクラッタ電力中
心周波数算出手段とし、クラッタレベル推定手段で推定
されたクラッタ・レベルに相当する減衰量のフィルタ係
数を選択するノッチフィルタ係数選択手段を付加したの
で、クラッタの中心周波数が変動する場合でも必要な目
標信号が得られる効果がある。
【００４４】また更に、クラッタ・レベル相当の減衰量
を設定したフィルタ係数に換えてクラッタ帯域幅相当の
阻止帯域幅を設定したフィルタ係数として選択するノッ
チフィルタ係数選択手段とし、クラッタレベル推定手段
に換えて、クラッタの帯域幅を推定するクラッタ帯域幅
推定手段とし、推定クラッタの帯域幅対応のノッチフィ
ルタ係数を発生させたので、クラッタに周波数と帯域変
動があっても必要な目標信号が得られる効果がある。
を設定したフィルタ係数に換えてクラッタ帯域幅相当の
阻止帯域幅を設定したフィルタ係数として選択するノッ
チフィルタ係数選択手段とし、クラッタレベル推定手段
に換えて、クラッタの帯域幅を推定するクラッタ帯域幅
推定手段とし、推定クラッタの帯域幅対応のノッチフィ
ルタ係数を発生させたので、クラッタに周波数と帯域変
動があっても必要な目標信号が得られる効果がある。
【００４５】また更に、クラッタ・レベル相当の減衰量
を設定したフィルタ係数に換えてクラッタ・レベル相当
の減衰量とクラッタ帯域幅相当の阻止帯域幅を設定した
フィルタ係数として選択するレベル帯域幅ノッチフィル
タ係数選択手段とし、クラッタレベル推定手段に換え
て、クラッタ・レベル帯域幅推定手段とし、推定クラッ
タ・レベルとクラッタ帯域幅対応のノッチフィルタ係数
を発生させたので、クラッタに周波数と帯域変動があっ
ても必要な目標信号が得られる効果がある。
を設定したフィルタ係数に換えてクラッタ・レベル相当
の減衰量とクラッタ帯域幅相当の阻止帯域幅を設定した
フィルタ係数として選択するレベル帯域幅ノッチフィル
タ係数選択手段とし、クラッタレベル推定手段に換え
て、クラッタ・レベル帯域幅推定手段とし、推定クラッ
タ・レベルとクラッタ帯域幅対応のノッチフィルタ係数
を発生させたので、クラッタに周波数と帯域変動があっ
ても必要な目標信号が得られる効果がある。
【００４６】また更に、入力信号の電力スペクトルから
クラッタ対雑音比を算出するＣ／Ｎ比算出手段と、この
Ｃ／Ｎ比でフィルタ処理の実行の有無を切り換えるスイ
ッチを付加したので、影響の少ないクラッタの除去によ
る目標信号の減衰を防ぐ効果がある。
クラッタ対雑音比を算出するＣ／Ｎ比算出手段と、この
Ｃ／Ｎ比でフィルタ処理の実行の有無を切り換えるスイ
ッチを付加したので、影響の少ないクラッタの除去によ
る目標信号の減衰を防ぐ効果がある。
【図１】 この発明の実施の形態１のクラッタ抑圧装置
を示す構成ブロック図である。
を示す構成ブロック図である。
【図２】 海面からの反射信号に含まれるシークラッタ
のドップラスペクトルを説明する図である。
のドップラスペクトルを説明する図である。
【図３】 この発明の実施の形態２のクラッタ抑圧装置
を示す構成ブロック図である。
を示す構成ブロック図である。
【図４】 この発明の実施の形態３のクラッタ抑圧装置
を示す構成ブロック図である。
を示す構成ブロック図である。
【図５】 この発明の実施の形態３の他のクラッタ抑圧
装置を示す構成ブロック図である。
装置を示す構成ブロック図である。
【図６】 この発明の実施の形態３の他のクラッタ抑圧
装置を示す構成ブロック図である。
装置を示す構成ブロック図である。
【図７】 この発明の実施の形態４のクラッタ抑圧装置
を示す構成ブロック図である。
を示す構成ブロック図である。
【図８】 従来技術のＭＴＩの機能系統図である。
【図９】 従来技術の単一消去ＭＴＩの特性を表す図で
ある。
ある。
１ 位相検波ビデオ信号、２ 遅延回路、３ 演算回
路、４ クラッタ除去後の位相検波ビデオ、５ レーダ
送信周波数、６ クラッタ中心周波数算出部（ａ）、７
ノッチフィルタ係数発生部（ａ）、８ フィルタ処理
部、９ 電力スペクトル算出部、１０ クラッタレベル
推定部、１１ ノッチフィルタ係数選択部（ａ）、１２
ノッチフィルタ係数発生部（ｂ）、１３ クラッタ中
心周波数算出部（ｂ）、１４ クラッタ帯域幅推定部、
１５ ノッチフィルタ係数選択部（ｂ）、１６ クラッ
タレベル・帯域幅推定部、１７ ノッチフィルタ係数選
択部（ｃ）、１８ 雑音レベル推定部、１９ Ｃ／Ｎ比
算出部、２０ 抑圧処理判定部、２１ スイッチ、２２
出力信号。
路、４ クラッタ除去後の位相検波ビデオ、５ レーダ
送信周波数、６ クラッタ中心周波数算出部（ａ）、７
ノッチフィルタ係数発生部（ａ）、８ フィルタ処理
部、９ 電力スペクトル算出部、１０ クラッタレベル
推定部、１１ ノッチフィルタ係数選択部（ａ）、１２
ノッチフィルタ係数発生部（ｂ）、１３ クラッタ中
心周波数算出部（ｂ）、１４ クラッタ帯域幅推定部、
１５ ノッチフィルタ係数選択部（ｂ）、１６ クラッ
タレベル・帯域幅推定部、１７ ノッチフィルタ係数選
択部（ｃ）、１８ 雑音レベル推定部、１９ Ｃ／Ｎ比
算出部、２０ 抑圧処理判定部、２１ スイッチ、２２
出力信号。
───────────────────────────────────────────────────── フロントページの続き (72)発明者 高橋 和之 東京都千代田区丸の内二丁目２番３号 三 菱電機株式会社内 (72)発明者 水谷 明義 東京都千代田区丸の内二丁目２番３号 三 菱電機株式会社内
Claims (6)
- 【請求項１】 レーダ送信周波数よりクラッタ中心周波
数を推定するクラッタ中心周波数算出手段と、 上記推定されたクラッタ中心周波数と、予め設定したク
ラッタレベルとクラッタ帯域幅に対応したノッチフィル
タ係数とからフィルタの係数を算出するノッチフィルタ
係数発生手段とを備え、 受信した入力信号に上記フィルタ係数を用いたノッチフ
ィルタ処理によりクラッタを除去するクラッタ抑圧装
置。 - 【請求項２】 入力信号の電力スペクトルと推定された
クラッタ中心周波数からクラッタのレベルを算出するク
ラッタレベル推定手段と、予め各種のクラッタ・レベル
に相当する減衰量を設定したフィルタ係数から上記推定
されたクラッタ・レベルに相当する減衰量のフィルタ係
数を選択するノッチフィルタ係数選択手段を付加し、 上記推定されたクラッタ・レベル対応のノッチフィルタ
係数を発生させてノッチフィルタ処理によりクラッタを
除去することを特徴とする請求項１記載のクラッタ抑圧
装置。 - 【請求項３】 レーダ送信周波数よりクラッタ中心周波
数を推定するクラッタ中心周波数算出手段に換えて、入
力信号の電力スペクトルとレーダ送信周波数からクラッ
タ中心周波数を推定するクラッタ電力中心周波数算出手
段とし、入力信号の電力スペクトルと上記推定されたク
ラッタ電力中心周波数からクラッタのレベルを算出する
クラッタレベル推定手段と、予め各種のクラッタ・レベ
ルに相当する減衰量を設定したフィルタ係数から上記推
定されたクラッタ・レベルに相当する減衰量のフィルタ
係数を選択するノッチフィルタ係数選択手段を付加し、
上記推定されたクラッタ・レベル対応のノッチフィルタ
係数を発生させてノッチフィルタ処理によりクラッタを
除去することを特徴とする請求項１記載のクラッタ抑圧
装置。 - 【請求項４】 各種のクラッタ・レベルに相当する減衰
量を設定したフィルタ係数から対応フィルタ係数を選択
するノッチフィルタ係数選択手段に換えて、各種のクラ
ッタ帯域幅に相当する阻止帯域幅を設定したフィルタ係
数から対応帯域幅のフィルタ係数を選択するノッチフィ
ルタ係数選択手段とし、クラッタレベル推定手段に換え
て、入力信号の電力スペクトルと上記推定されたクラッ
タ電力中心周波数からクラッタの帯域幅を推定するクラ
ッタ帯域幅推定手段とし、上記推定されたクラッタの帯
域幅対応のノッチフィルタ係数を発生させてノッチフィ
ルタ処理によりクラッタを除去することを特徴とする請
求項３記載のクラッタ抑圧装置。 - 【請求項５】 各種のクラッタ・レベルに相当する減衰
量を設定したフィルタ係数から対応フィルタ係数を選択
するノッチフィルタ係数選択手段に換えて、各種のクラ
ッタ・レベルに相当する減衰量とクラッタ帯域幅に相当
する阻止帯域幅を設定したフィルタ係数から対応レベル
と帯域幅のフィルタ係数を選択するノッチフィルタ係数
選択手段とし、クラッタレベル推定手段に換えて、入力
信号の電力スペクトルと上記推定されたクラッタ電力中
心周波数からクラッタ・レベルとクラッタ帯域幅を推定
するクラッタ・レベル帯域幅推定手段とし、 上記推定されたクラッタ・レベルとクラッタ帯域幅対応
のノッチフィルタ係数を発生させてノッチフィルタ処理
によりクラッタを除去することを特徴とする請求項３記
載のクラッタ抑圧装置。 - 【請求項６】 入力信号の電力スペクトルからクラッタ
対雑音比を算出するＣ／Ｎ比算出手段と、上記算出され
たＣ／Ｎ比によりフィルタ処理を行うかまたは入力信号
を出力するかを切り換えるスイッチを付加し、 上記算出されたＣ／Ｎ比によりノッチフィルタ処理によ
るクラッタ除去の採否をすることを特徴とする請求項２
ないし請求項５いずれか記載のクラッタ抑圧装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP03392697A JP3145943B2 (ja) | 1997-02-18 | 1997-02-18 | クラッタ抑圧装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP03392697A JP3145943B2 (ja) | 1997-02-18 | 1997-02-18 | クラッタ抑圧装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH10227851A true JPH10227851A (ja) | 1998-08-25 |
| JP3145943B2 JP3145943B2 (ja) | 2001-03-12 |
Family
ID=12400133
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP03392697A Expired - Lifetime JP3145943B2 (ja) | 1997-02-18 | 1997-02-18 | クラッタ抑圧装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3145943B2 (ja) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006349477A (ja) * | 2005-06-15 | 2006-12-28 | Mitsubishi Electric Corp | レーダ装置 |
| KR100851888B1 (ko) | 2006-10-26 | 2008-08-13 | 국방과학연구소 | 클러터와 송신 누설 전력 제거가 가능한 레이더 |
| WO2012170135A1 (en) * | 2011-06-08 | 2012-12-13 | Raytheon Company | Method and apparatus for mitigating an effect of user movement in motion detecting radar |
| JP2013124971A (ja) * | 2011-12-15 | 2013-06-24 | Mitsubishi Electric Corp | クラッタ抑圧装置 |
| RU174360U1 (ru) * | 2017-04-10 | 2017-10-11 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Рязанский государственный радиотехнический университет" | Вычислительное устройство режекции пассивных помех |
| RU2646330C1 (ru) * | 2017-01-24 | 2018-03-02 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Рязанский государственный радиотехнический университет" | Вычислитель для режекторной фильтрации помех |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU182317U1 (ru) * | 2018-05-07 | 2018-08-14 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Рязанский государственный радиотехнический университет" | Вычислитель для адаптивной фильтрации помех |
| RU183016U1 (ru) * | 2018-05-29 | 2018-09-07 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Рязанский государственный радиотехнический университет" | Вычислительное устройство подавления помех |
| RU182703U1 (ru) * | 2018-05-29 | 2018-08-28 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Рязанский государственный радиотехнический университет" | Вычислительное устройство режекции помех |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH05142339A (ja) * | 1991-11-19 | 1993-06-08 | Mitsubishi Electric Corp | クラツタ抑圧装置 |
| JPH05203733A (ja) * | 1992-01-27 | 1993-08-10 | Mitsubishi Electric Corp | レーダ装置 |
-
1997
- 1997-02-18 JP JP03392697A patent/JP3145943B2/ja not_active Expired - Lifetime
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH05142339A (ja) * | 1991-11-19 | 1993-06-08 | Mitsubishi Electric Corp | クラツタ抑圧装置 |
| JPH05203733A (ja) * | 1992-01-27 | 1993-08-10 | Mitsubishi Electric Corp | レーダ装置 |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006349477A (ja) * | 2005-06-15 | 2006-12-28 | Mitsubishi Electric Corp | レーダ装置 |
| JP4727311B2 (ja) * | 2005-06-15 | 2011-07-20 | 三菱電機株式会社 | レーダ装置 |
| KR100851888B1 (ko) | 2006-10-26 | 2008-08-13 | 국방과학연구소 | 클러터와 송신 누설 전력 제거가 가능한 레이더 |
| WO2012170135A1 (en) * | 2011-06-08 | 2012-12-13 | Raytheon Company | Method and apparatus for mitigating an effect of user movement in motion detecting radar |
| JP2013124971A (ja) * | 2011-12-15 | 2013-06-24 | Mitsubishi Electric Corp | クラッタ抑圧装置 |
| RU2646330C1 (ru) * | 2017-01-24 | 2018-03-02 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Рязанский государственный радиотехнический университет" | Вычислитель для режекторной фильтрации помех |
| RU174360U1 (ru) * | 2017-04-10 | 2017-10-11 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Рязанский государственный радиотехнический университет" | Вычислительное устройство режекции пассивных помех |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3145943B2 (ja) | 2001-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6031267B2 (ja) | 干渉検出装置、干渉除去器、レーダ装置、干渉検出方法および干渉検出用プログラム | |
| JP2005530164A (ja) | 物体検出用システムにおける妨害を抑制するための方法 | |
| JP3145943B2 (ja) | クラッタ抑圧装置 | |
| US20200371236A1 (en) | Distance measurement apparatus | |
| JP2008256626A (ja) | 自動利得制御装置 | |
| JPS6111380B2 (ja) | ||
| JPH0529875B2 (ja) | ||
| JP2021067461A (ja) | レーダ装置及びレーダ信号処理方法 | |
| KR100522607B1 (ko) | 노이즈 상태를 고려한 적응적 비디오 신호 처리 장치 및방법 | |
| JP5858755B2 (ja) | レーダ装置 | |
| US4453223A (en) | Device for computing the pseudo-variance of the instantaneous periods of an alternating-current signal | |
| EP0284498A1 (en) | Frequency counter in particular for counting partly varying frequencies in real time | |
| JP3296234B2 (ja) | クラッタ抑圧装置およびクラッタ抑圧方法 | |
| JP3060012B2 (ja) | レーダ装置 | |
| JP2889226B1 (ja) | レーダ装置 | |
| Szajnowski et al. | Implicit averaging and delay determination of random binary waveforms | |
| JP2002141885A (ja) | Ｏｆｄｍ受信機 | |
| JPH05142339A (ja) | クラツタ抑圧装置 | |
| JP2005017143A (ja) | 気象レーダ信号処理装置 | |
| JPH0990025A (ja) | レーダクラッタ抑圧装置 | |
| JP3389127B2 (ja) | レーダー信号処理装置およびレーダー信号処理装置の信号処理方法 | |
| JP2837614B2 (ja) | 不要信号抑圧装置 | |
| JP2685041B2 (ja) | 移動目標の自動検出方式 | |
| JP3607846B2 (ja) | クラッタ抑圧装置 | |
| Agarwal et al. | VLSI architecture for IEEE single precision floating point moving average calculator |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090105 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100105 Year of fee payment: 9 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110105 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110105 Year of fee payment: 10 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |