JPH04361186A - レーダ用クラッタ信号抑圧装置 - Google Patents
レーダ用クラッタ信号抑圧装置Info
- Publication number
- JPH04361186A JPH04361186A JP3136286A JP13628691A JPH04361186A JP H04361186 A JPH04361186 A JP H04361186A JP 3136286 A JP3136286 A JP 3136286A JP 13628691 A JP13628691 A JP 13628691A JP H04361186 A JPH04361186 A JP H04361186A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- circuit
- subtraction
- log
- cfar
- 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
- 238000004364 calculation method Methods 0.000 claims abstract description 47
- 230000001629 suppression Effects 0.000 claims description 3
- 230000003321 amplification Effects 0.000 claims 1
- 238000003199 nucleic acid amplification method Methods 0.000 claims 1
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 230000003111 delayed effect Effects 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 101000582320 Homo sapiens Neurogenic differentiation factor 6 Proteins 0.000 description 1
- 102100030589 Neurogenic differentiation factor 6 Human genes 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000034 method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
め要約のデータは記録されません。
Description
らの不必要な信号としてレーダ信号に混在するクラッタ
信号を抑圧するレーダ用クラッタ信号抑圧装置に関する
。
る船舶や航空機等の物標からの必要な信号と、クラッタ
と呼ばれる海面、雨、雲、大地等からの不必要な信号と
が混在している。レーダの信号処理においては、クラッ
タをできるだけ低く抑え、ターゲットを検出することを
目的としている。更に、クラッタが誤ってターゲットと
判定される確率を示す誤警報確率を一定値以下に抑え、
一定誤警報確率CFAR(Constant Fals
e Alarm Rate)を得る処理が非常に重要で
ある。
を抑圧する信号処理方式としてlog−CFAR処理方
式が知られている。従来のlog−CFAR処理方式を
図4によって説明する。図4において、400は対数増
幅器であり入力信号を対数変換する。401は遅延手段
であり、対数変換されたレーダービデオ信号をN段分遅
延させる。402は加算回路であり、遅延回路402の
N段の各遅延タップから出力される個々出力信号の総和
を演算する。403は平均値演算回路であり、加算回路
402からの総和信号を1/Nして平均値信号を算出す
る。204は減算回路であり、遅延回路401を通過し
たレーダビデオ信号から平均値演算回路403による平
均値信号を減算してlog−CFAR出力信号Wを算出
する。
信号が(1)式に示すレーリー分布信号にしたがう場合
には、その出力信号の確率密度関数p(t)が(2)式
に示すごとくレーリー分布信号の分散値bによらない信
号となる。尚、(2)式においてγはオイラー定数であ
る。
がある一定閾値T以上になる場合の確率を誤警報確率と
考えると、誤警報確率の値は(2)式をTから∞まで積
分することによって得られる。計算の結果、誤警報確率
は(3)式となりレーリー分布信号の分散値bによらず
一定の値となり、log−CFAR出力信号はCFAR
化されたことになる。
信号に於いてクラッタ信号はレーリー分布に従うばかり
ではなく、特に海面反射信号などはログノーマル分布(
log−normal分布)に従う場合も多く知られて
いる。このようにログノーマル分布に従うクラッタ信号
に対して従来のlog−CFAR処理を適用しても、そ
の出力信号はCFAR化されずにクラッタ信号の消え残
りが生じることになり、クラッタ中のターゲットを検出
することが困難となるという問題があった。
てなされたもので、従来のlog−CFAR処理におい
てはCFAR化が困難であったログノーマル分布による
クラッタ信号に対して誤警報確率を一定にしてCFAR
化できるレーダ用クラッタ信号抑圧装置を提供すること
を目的とする。
め本発明のレーダ用クラッタ信号抑圧装置にあっては次
のように構成する。尚、実施例図面中の対応する符号を
併せて示す。即ち、本発明のレーダ用クラッタ信号抑圧
装置は、レーダ信号(X)を対数変換した対数信号(Y
)を出力する対数増幅手段(1)と、対数信号(Y)の
移動平均値(<Y>)を算出する平均値演算手段(2)
と、対数信号(Y)から移動平均値(<Y>)を減算し
て第1の減算信号(V)を算出する第1の減算手段(3
)と、平均値演算手段(2)で演算対象となる範囲の対
数信号の最大値を算出する最大値算出手段(4)と、平
均値演算手段(2)で演算対象となる範囲の対数信号の
最小値を算出する最小値算出手段(5)と、前記最大値
から最小値を減算して第2の減算信号を算出する第2の
減算手段(6)と、第2の減算信号と、目標とする誤警
報確率と平均値演算手段(2)で演算対象となる範囲に
相当するサンプルデータ数によって定められる定数との
積を算出し閾値信号(T1)を出力する掛算手段(7)
と、第1の減算信号(V)から閾値信号(T1)を減算
する第3の減算手段とを備えたことを特徴とする。
ッタ信号抑圧装置によれば、従来のlog−CFAR処
理ではCFAR化が困難であったログノーマル分布信号
に対しても、ログノーマル分布信号の対数変換信号が正
規分布信号となることに着目し、正規分布信号の標準偏
差値を不変推定値を算出するのではなく、正規分布信号
の最大値と最小値から標準偏差値を推定し、この標準偏
差値から誤警報率を一定に保つ閾値を求め、最終的にl
og−CFAR出力信号から差し引くことで、ログノー
マル分布に従う信号からCFAR化された信号を得るこ
とができる。
ず確率密度関数p(x)が(4)式で示される確率変数
xをログノーマル分布という。
出力信号Yは(5)式で表され、その平均値<Y>は(
6)式となる。
FAR出力Vは(7)式で表される。
数P(V)は(8)式で表されるので、閾値をTと定め
たときのlog−CFAR出力Vの誤警報確率FVは(
9)式となる。
og−CFAR出力における誤警報確率Fv(T)は、
ログノーマル分布信号xの分散σと対数増幅器の定数a
の関数となる。いま、対数増幅器の特性を一定とすると
、ログノーマル分布信号xのlog−CFAR出力誤警
報確率は分散σの値によって変動し、CFAR化出力が
得られないことが理解される。
たときの閾値T1の値を求める式に変形すると(10)
式となる。これにより閾値T1は目標とする誤警報確率
Fvとlog−CFAR出力の標準偏差(VRY)1/
2 によって定められることが解る。尚、対数増幅器通
過出力信号Yの分散値は(11)式によって与えられる
。
偏差信号を求めて(10)式に基づき閾値T1を設定す
ればログノーマル分布信号xに対するCFAR出力が求
められることになる。一方、対数増幅器通過出力信号Y
はログノーマル分布信号xを対数増幅した信号であり、
その確率分布は正規分布に従う。従って、その標準偏差
は、推定区間のデータの最大値、最小値から推定可能と
なる。
差値<VRY>1/2 は(12)式に示すごとく、サ
ンプル集合の最大値Ymaxと最小値Yminiの差に
係数Aを掛けることで推定できる。サンプル集合のデー
タ個数による定数Aの値はティペット(Tippet)
により計算されており、統計学者のためのバイオメトリ
カ表(Biometrika Table for
Statisticians)に示されている。
誤警報確率を一定に保つための閾値T1は(13)式の
ごとく、対数増幅器通過出力信号Yの最大値、最小値の
差に定数Bを掛けることで求められる。
式で定められる閾値T1を減算すると一定の誤警報確率
Fv以上の信号のみ出力されることになり、CFAR化
出力が得られることになる。尚、ここでの最大値、最小
値を対数増幅器通過出力信号Yについて説明したが、l
og−CFAR出力信号V(減算信号V)について最大
値、最小値を求めても同様の結果が得られることは言う
までもない。
成図である。図1において、1は対数増幅器であり、2
は対数変換された対数信号Yの移動平均値<Y>を算出
する平均値演算回路である。3は対数信号Yからその平
均値<Y>を減算し従来のlog−CFAR化出力に相
当する第1の減算信号Vを出力する第1の減算回路であ
る。
われる演算対象となる信号の範囲でその最大値を算出す
る最大値算出回路であり、5は平均値演算回路2で移動
平均演算が行われる演算対象となる信号の範囲でその最
小値を算出する最小値算出回路である。6は最大値算出
回路4による最大値から最小値算出回路5による最小値
を減算し第2の減算信号を算出する第2の減算回路であ
る。7は第2の減算信号と目標とする誤警報確率値と移
動平均演算に用いられるデータサンプル数に従って定め
られる定数の積を算出し閾値信号T1を出力する掛算回
路である。更に、8は第1の減算信号Vから閾値信号T
1を減算し最終的にCFAR化出力信号を算出する第3
の減算回路である。
について説明する。ログノーマル分布信号Xは対数増幅
器1によって正規分布信号Yに変換される。正規分布信
号Yの移動平均値<Y>は平均値演算回路2によって算
出される。ここで平均値算出回路2は従来のlog−C
FAR処理で用いられていた図4に示した遅延回路40
1、加算回路402及び平均値演算回路403によって
実現できる。
その移動平均値<Y>の減算が行われ、従来のlog−
CFAR化出力に相当する第1の減算信号Vが算出され
る。最大値算出回路4では平均値演算回路2で演算対象
となる区間に対応する信号Yの最大値を算出して出力す
る。例えば、平均値演算回路2における図4に示した遅
延回路401の遅延段数が8段で構成されているとする
と、それに対応する最大値算出回路4は図2の構成によ
って実現できる。
それぞれ入力信号V0,V1,V2を連続して遅延させ
て出力する1段、2段、4段の遅延回路、202、20
4、206は2つの入力信号の振幅の大きい方を選択し
て出力する比較回路である。図3は図2の最大値算出回
路における各部の信号を示す。
1段遅延回路201で遅らせた信号を比較し、大きい方
を信号V1として出力する。これを図3の最初の2つの
V1信号では(1,2)で示しており、V0信号1とV
0信号2のうちの大きい方を出力することを意味する。 以下、残りのV0信号2〜12についても同様である。
段遅延回路203で遅らせた信号とを比較し、大きい方
を信号V2として出力する。これを図3の最初のV2信
号では(1〜4)で示しており、4つのV0信号1〜4
のうちの最も大きいものを出力することを意味する。以
下、残りのV2信号(2〜5)、(3〜6)、・・・(
9〜12)についても同様である。
段遅延回路205で遅らせた信号とを比較し、大きい方
を信号V3として出力する。これを図3の最初のV3信
号では(1〜8)で示しており、8つのV0信号1〜8
のうちの最も大きいものを出力することを意味する。以
下、残りのV3信号(2〜9)、(3〜10)・・・(
5〜12)についても同様である。
に相当する信号区間の中の最大値がV3信号として算出
されることになる。再び図1を参照するに、最小値算出
回路5では平均値演算回路2で演算対象となる区間に対
応する信号Yの最小値を算出し出力する。最小値算出回
路5の構成としては、図2に示した最大値算出回路にお
ける比較回路202、204、206の機能を2つの入
力信号の内振幅の小さい信号を選択し出力するようにす
れば実現できる。
4からの最大値から最小値算出回路5からの最小値との
差、即ち第2の減算信号が算出される。この第2の減算
信号は掛算回路7に与えられ、第2の減算信号に目標と
する誤警報確率値と移動平均演算に用いられるデータサ
ンプル数に従って定められる定数との積が演算され閾値
信号T1として第3の減算回路8に出力される。
用いられるデータサンプル数が十分に大きければ、前記
(13)式における定数Bとほぼ一致した値となる。デ
ータサンプル数が小さい場合にはコンピュータシミュレ
ーションによって目的とする誤警報確率が得られる定数
値を求めることが容易に可能である。第3の減算回路8
ではlog−CFAR化出力信号としての第1の減算信
号から閾値信号T1が減算され、最終結果としてログノ
ーマル分布信号xについてCFAR化された出力信号が
求められる。
来のlog−CFAR処理ではCFAR化が困難であっ
たログノーマル分布信号に対しても、CFAR化された
出力信号を算出することができる。またCFAR化の演
算過程においてログノーマル分布信号の対数変換信号が
正規分布信号となることに着目し、その標準偏差値を不
変推定値を算出することなく対数変換信号の最大値と最
小値から推定する構成とすることで、被推定信号の自乗
値を算出することなく求めることが可能となり、装置構
成を安価に実現できる。
図1の最大値選択回路の実施例を示した回路ブロック図
図
示した回路ブロック図
Claims (1)
- 【請求項1】レーダ信号を対数変換した対数信号を出力
する対数増幅手段と、該対数信号の移動平均値を算出す
る平均値演算手段と、前記対数信号から前記移動平均値
を減算して第1の減算信号を算出する第1の減算手段と
、前記平均値演算手段で演算対象となる範囲の対数信号
の最大値を算出する最大値算出手段と、前記平均値演算
手段で演算対象となる範囲の対数信号の最小値を算出す
る最小値算出手段と、前記最大値から最小値を減算して
第2の減算信号を算出する第2の減算手段と、該第2の
減算信号と、目標とする誤警報確率と前記平均値演算手
段で演算対象となる範囲に相当するサンプルデータ数に
よって定められる定数との積を算出し閾値信号を出力す
る掛算手段と、前記第1の減算信号から前記閾値信号を
減算する第3の減算手段と、を備えたことを特徴とする
レーダ用クラッタ信号抑圧装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP3136286A JP2902155B2 (ja) | 1991-06-07 | 1991-06-07 | レーダ用クラッタ信号抑圧装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP3136286A JP2902155B2 (ja) | 1991-06-07 | 1991-06-07 | レーダ用クラッタ信号抑圧装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH04361186A true JPH04361186A (ja) | 1992-12-14 |
JP2902155B2 JP2902155B2 (ja) | 1999-06-07 |
Family
ID=15171626
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP3136286A Expired - Fee Related JP2902155B2 (ja) | 1991-06-07 | 1991-06-07 | レーダ用クラッタ信号抑圧装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2902155B2 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008309606A (ja) * | 2007-06-14 | 2008-12-25 | Japan Radio Co Ltd | レーダ装置 |
CN107621635A (zh) * | 2017-08-21 | 2018-01-23 | 电子科技大学 | 一种前视海面目标角超分辨方法 |
CN107728111A (zh) * | 2017-09-22 | 2018-02-23 | 合肥工业大学 | 基于空间相关特性的sar图像联合cfar检测方法 |
-
1991
- 1991-06-07 JP JP3136286A patent/JP2902155B2/ja not_active Expired - Fee Related
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008309606A (ja) * | 2007-06-14 | 2008-12-25 | Japan Radio Co Ltd | レーダ装置 |
CN107621635A (zh) * | 2017-08-21 | 2018-01-23 | 电子科技大学 | 一种前视海面目标角超分辨方法 |
CN107621635B (zh) * | 2017-08-21 | 2020-09-01 | 电子科技大学 | 一种前视海面目标角超分辨方法 |
CN107728111A (zh) * | 2017-09-22 | 2018-02-23 | 合肥工业大学 | 基于空间相关特性的sar图像联合cfar检测方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2902155B2 (ja) | 1999-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5465405A (en) | Apparatus and method for detecting signals | |
JPS6111380B2 (ja) | ||
JPH0320714B2 (ja) | ||
JPH04361186A (ja) | レーダ用クラッタ信号抑圧装置 | |
JP2801904B2 (ja) | Ca log/cfar装置 | |
JPH023472B2 (ja) | ||
JPH04216484A (ja) | レーダ用クラッタ信号抑圧装置 | |
JPH0125342Y2 (ja) | ||
JPH035909Y2 (ja) | ||
JPS6216694Y2 (ja) | ||
JPH07311256A (ja) | クラッタ抑圧回路 | |
JP2596614B2 (ja) | 信号処理方法および信号処理装置 | |
JP2614157B2 (ja) | Stc制御信号発生回路 | |
JPS61205878A (ja) | レ−ダ−信号処理装置 | |
JP2678072B2 (ja) | レーダ信号処理方法及び装置 | |
JPS6035283A (ja) | 雑音抑圧装置 | |
JP2003227871A (ja) | レーダ装置のクラッタ抑圧回路 | |
JPH06148312A (ja) | 目標検出装置 | |
EP0083213A2 (en) | Digital background normalizer using two dimensional integration techniques | |
JPH062260U (ja) | 目標検出装置 | |
JPS61164172A (ja) | クラツタ除去装置 | |
JP2590118B2 (ja) | レーダ信号処理装置 | |
JPH05249227A (ja) | レーダ用クラッタ信号抑圧装置 | |
JPS5941148B2 (ja) | 信号処理装置 | |
JPH0464078A (ja) | レーダ信号処理装置 |
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 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080319 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090319 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090319 Year of fee payment: 10 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090319 Year of fee payment: 10 |
|
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: 20090319 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100319 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110319 Year of fee payment: 12 |
|
LAPS | Cancellation because of no payment of annual fees |