JPS6189575A - 雑音抑圧装置 - Google Patents
雑音抑圧装置Info
- Publication number
- JPS6189575A JPS6189575A JP59211118A JP21111884A JPS6189575A JP S6189575 A JPS6189575 A JP S6189575A JP 59211118 A JP59211118 A JP 59211118A JP 21111884 A JP21111884 A JP 21111884A JP S6189575 A JPS6189575 A JP S6189575A
- Authority
- JP
- Japan
- Prior art keywords
- output
- circuit
- outputs
- moving average
- input signal
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔産業上の利用分野〕
この発明はレーダ信号処理装置において。
レーダ受信信号の振幅分布が任意のワイプル分布をする
クラッタに対して、クラッタを抑圧し一定誤警報率(0
onstan−t False Alarm&te)を
実現する雑音抑圧装置に関するものである。
クラッタに対して、クラッタを抑圧し一定誤警報率(0
onstan−t False Alarm&te)を
実現する雑音抑圧装置に関するものである。
第4図は例えば、電子通信学会論文誌゛79/ l v
ol 、 J 62−B ff1l pp 45−4
9 「ワイプル分布をするレーダクラッタの抑圧」の4
8P図6に示された従来の雑音抑圧装置を示す図でめり
2図において(ハ)は入力信号を2乗し出力する2乗回
路、 (3a)はこの2乗回路出力の移動平均を出力
する第1の移動平均回路、α2はこの第1の移動平均回
路出力の平方根を出力する平方根回路、 (3b)は
前記入力信号の移動平均を出力する第2の移動平均回路
、(5)はこの第2の移動平均回路出力を前記平方根回
路α2で割った商Qを出力する除算器、 (6b)は
この除算器出力Qを導入し第1の出力としてYl、第2
の出力としてY2.ここに。
ol 、 J 62−B ff1l pp 45−4
9 「ワイプル分布をするレーダクラッタの抑圧」の4
8P図6に示された従来の雑音抑圧装置を示す図でめり
2図において(ハ)は入力信号を2乗し出力する2乗回
路、 (3a)はこの2乗回路出力の移動平均を出力
する第1の移動平均回路、α2はこの第1の移動平均回
路出力の平方根を出力する平方根回路、 (3b)は
前記入力信号の移動平均を出力する第2の移動平均回路
、(5)はこの第2の移動平均回路出力を前記平方根回
路α2で割った商Qを出力する除算器、 (6b)は
この除算器出力Qを導入し第1の出力としてYl、第2
の出力としてY2.ここに。
Yl= −・・・・・・・・・・・・・・・・・・・・
・・・・・・・・・・・・・・・・・・・・・ (1+
Y、 = −・・・・・・・・・・・・・・・・・・・
・・・・・・・・・・・・・・・・・・・・ +21と
なるようなr Y1+ Y2 を出力する演算器。
・・・・・・・・・・・・・・・・・・・・・ (1+
Y、 = −・・・・・・・・・・・・・・・・・・・
・・・・・・・・・・・・・・・・・・・・ +21と
なるようなr Y1+ Y2 を出力する演算器。
(4)は前記入力信号を対数変換し出力する対数変換回
路、 (3c)はこの対数変換回路の出力の移動平均
を第1の出力とし、その移動平均幅の2分の1遅延した
前記対数変換回路(4)出力信号を第2の出力として出
力する第3の移動平均回路、(7)は前記演算回路(6
)の第1の出力と前記第3の移動平均回路(3C)の第
1の出力を加算し出力する加算器、(8)は前記第3の
移動平均回路(3りの第2の出力から前記加算器(7)
の出力を減算し出力する減算器、(9)はこの減算器出
力と前記演算器(6b)の第2の出力の積を出力する乗
算器、αQはこの乗算器出力の逆対数変換を出力する逆
対数変換回路であるここで、入力信号Xがワイブル分布
に従う場合、ワイブル分布の確率照度関数P (x)は
次式で表せる。
路、 (3c)はこの対数変換回路の出力の移動平均
を第1の出力とし、その移動平均幅の2分の1遅延した
前記対数変換回路(4)出力信号を第2の出力として出
力する第3の移動平均回路、(7)は前記演算回路(6
)の第1の出力と前記第3の移動平均回路(3C)の第
1の出力を加算し出力する加算器、(8)は前記第3の
移動平均回路(3りの第2の出力から前記加算器(7)
の出力を減算し出力する減算器、(9)はこの減算器出
力と前記演算器(6b)の第2の出力の積を出力する乗
算器、αQはこの乗算器出力の逆対数変換を出力する逆
対数変換回路であるここで、入力信号Xがワイブル分布
に従う場合、ワイブル分布の確率照度関数P (x)は
次式で表せる。
ここで、aは形状パラメータ、bはスケール・パラメー
タである。
タである。
Xの平均値E(X’lと2乗平均値ECX”lは次式%
式% E(x’lとJ:[i[x’]の比をとると。
式% E(x’lとJ:[i[x’]の比をとると。
とlる。E[Xlを入力Xの移動平均で推定し。
E[X′〕をXの2乗の移動平均の平方根で推定すれば
。
。
となり、形状パラメータaのみの関数となる。したがっ
てQを求めれば式(8)から2が求まる。Qは除算器(
5)の出力でめる。
てQを求めれば式(8)から2が求まる。Qは除算器(
5)の出力でめる。
式(4)のワイブル分布に従う入力信号Xを対数変換回
路(4)に通し、その出力をyとおくとy吐lnx
・・・・・・・・・・・・・・・
(91このyの平均値ECy)は EIJ1=J。 lnx −P(x)d x” ln
b−−・・・・・・・・・・・・・・α1となる。ここ
でrはオイラーの定数で、γ=0.57721 でる
る。平均値E Cy〕を移動平均回路(3C)で推定す
れば。
路(4)に通し、その出力をyとおくとy吐lnx
・・・・・・・・・・・・・・・
(91このyの平均値ECy)は EIJ1=J。 lnx −P(x)d x” ln
b−−・・・・・・・・・・・・・・α1となる。ここ
でrはオイラーの定数で、γ=0.57721 でる
る。平均値E Cy〕を移動平均回路(3C)で推定す
れば。
ln b = E [y]+−
となる。これは加算器(7)の出力である。
ここで、対数変換回路(4)の出力信号yと1jnbの
差は次式のようになる。
差は次式のようになる。
Y /?nb = lnX −enb =
lln (、) −−・ αtこれが減算器
(8)の出力になる。
lln (、) −−・ αtこれが減算器
(8)の出力になる。
これにa/2を乗すると。
X !!−
一・/!!n () = 6n () 2 ・・
・・・・・・・・αJ、2 bb となり1乗算器(9)の出力になる。
・・・・・・・・αJ、2 bb となり1乗算器(9)の出力になる。
逆対数変換回路Qlの出力Zは次式で与えられる。
したがって、Zの確率密度関数Pz(Z)は次式で与え
られる。
られる。
Pz(句=P−(x ) a z ・・・・・
・・・・・・・・・・aSx =2Ze” このZが雑音抑圧装置の出力となる。式0Sは、第4図
に示した従来の雑音抑圧装置の出力の分布が、入力信号
の分布の形状パラメータaやスケール・パラメータbに
依存せずに一定でるることを示している。したがって。
・・・・・・・・・・aSx =2Ze” このZが雑音抑圧装置の出力となる。式0Sは、第4図
に示した従来の雑音抑圧装置の出力の分布が、入力信号
の分布の形状パラメータaやスケール・パラメータbに
依存せずに一定でるることを示している。したがって。
固定スレッショルドTに対する誤警報確率Pfaは。
■ −が
Pea = fo2Ze aZ=(−e JT=8
イ“ 11109101090.20.aeの
ように、一定になる。
イ“ 11109101090.20.aeの
ように、一定になる。
上記のような従来の雑音抑圧装置は、多くの非線形変換
すなわち、2乗回路aυ、平方根回路aZ、除算器(5
)、演算器(6b)、乗算器(9)。
すなわち、2乗回路aυ、平方根回路aZ、除算器(5
)、演算器(6b)、乗算器(9)。
対数変換回路(4)、逆対数変換回路a1かめるため9
ハードウエアで実現する場合、特にディジタル信号処理
を行う場合、ハードウェア量が多くなるという問題点が
めった。
ハードウエアで実現する場合、特にディジタル信号処理
を行う場合、ハードウェア量が多くなるという問題点が
めった。
この発明は、かかる問題点を解決するためになされたも
ので、2乗回路αυと平方根回路Q3が不要と永る雑音
抑圧装置を得ることを目的とする。
ので、2乗回路αυと平方根回路Q3が不要と永る雑音
抑圧装置を得ることを目的とする。
この発明に係る雑音抑圧装置は、ワイプル分布の形状パ
ラメータaを求める原理を変えたものである。すなわち
、従来の雑音抑圧装置が入力の2乗平均平方根と平均値
の比から形状パラメータaを求めていたのに対して。
ラメータaを求める原理を変えたものである。すなわち
、従来の雑音抑圧装置が入力の2乗平均平方根と平均値
の比から形状パラメータaを求めていたのに対して。
この発明は入力の区間最小値の平均と入力の平均値の比
から形状パラメータaを求めるようにしたものである。
から形状パラメータaを求めるようにしたものである。
この発明においては、入力の区間最小値の平均と入力の
平均値の比から形状パラメータaを求めているため、2
乗回路と平方根回路が不要となり、ハードウェアを減ら
すことができる。
平均値の比から形状パラメータaを求めているため、2
乗回路と平方根回路が不要となり、ハードウェアを減ら
すことができる。
また、この発明においては、演算器(6)の入出力関係
を解析的に解けるため、演算器(6)の設計が容易にな
るという作用もるる。
を解析的に解けるため、演算器(6)の設計が容易にな
るという作用もるる。
第1図はこの発明の一実施例を示す図でろり9演算器(
6a)を除り(3)〜tttu家上記従来の装置と全く
同一のものでるる。(1)は入力信号の相関時間以上の
時間だけ入力信号を遅延きせる遅延回路、(2)はこの
遅延回路出力と前記入力信号の振幅を比軟し、振幅の小
さい方を出力する最小値検出回路、 (6a)は除算器
(5)の出力Qを導入し、第1の出力としてYl、第2
の出力としてYl、ここに。
6a)を除り(3)〜tttu家上記従来の装置と全く
同一のものでるる。(1)は入力信号の相関時間以上の
時間だけ入力信号を遅延きせる遅延回路、(2)はこの
遅延回路出力と前記入力信号の振幅を比軟し、振幅の小
さい方を出力する最小値検出回路、 (6a)は除算器
(5)の出力Qを導入し、第1の出力としてYl、第2
の出力としてYl、ここに。
Y、 =γlog2 Q ・・・・・・
・・・・・・・・−αn(rはオイラーの定数) となるようなYl、Yl を出力する演算器でるる。
・・・・・・・・−αn(rはオイラーの定数) となるようなYl、Yl を出力する演算器でるる。
第2図(a)は移動平均回路(3a)(3b)の構成例
を示す図でめり9図中、 (la)は入力を連続して
遅延させ出力するN段(N=8)の遅延回路、 (7a
)はこの遅延回路のM個の出力を加算し出力する加算器
である。
を示す図でめり9図中、 (la)は入力を連続して
遅延させ出力するN段(N=8)の遅延回路、 (7a
)はこの遅延回路のM個の出力を加算し出力する加算器
である。
第2図(b)は移動平均回路(3C)の構成例を示す図
であり9図中、 (Ib)は入力を連続して遅延させ出
力するN段(N=8)の出力と遅延量の%の時点の出力
を持つ遅延回路、 (7a)はこの遅延回路のN個の
出力を加算し出力する加算器でろる。
であり9図中、 (Ib)は入力を連続して遅延させ出
力するN段(N=8)の出力と遅延量の%の時点の出力
を持つ遅延回路、 (7a)はこの遅延回路のN個の
出力を加算し出力する加算器でろる。
第3図(a)は最小値検出回路(2)の構成例を示す図
であり1図中、αJは第1及び第2の入力の振幅を比較
し、第1の入力が第2の入力よりも太きいときに論理”
l”、そうでないときに論理”o”を出力する比較器、
α心はこの比較器出力を制御信号とし、前記第1及び第
2の入力を入力し、制御信号が論理“0”のとき第1の
入力を、制御信号が論理“0″のとき第2の入力を選択
し出力するセレクタである。第3図(b)はこのセレク
タの入出力論理を示す図である。
であり1図中、αJは第1及び第2の入力の振幅を比較
し、第1の入力が第2の入力よりも太きいときに論理”
l”、そうでないときに論理”o”を出力する比較器、
α心はこの比較器出力を制御信号とし、前記第1及び第
2の入力を入力し、制御信号が論理“0”のとき第1の
入力を、制御信号が論理“0″のとき第2の入力を選択
し出力するセレクタである。第3図(b)はこのセレク
タの入出力論理を示す図である。
上記のように構成された雑音抑圧装置では、ワイブル分
布の形状パラメータaを求める方法のみ従来の装置と異
る。
布の形状パラメータaを求める方法のみ従来の装置と異
る。
ここで、入力信号Xがワイブル分布に従う場合、ワイブ
ル分布の確率密度関数Px(x)は式(4)で与えられ
るので、最小値検出回路の出力Uの確率密度関数Pu(
u)は次式で与えられる。
ル分布の確率密度関数Px(x)は式(4)で与えられ
るので、最小値検出回路の出力Uの確率密度関数Pu(
u)は次式で与えられる。
Pu(u) = 2Px(u) ・f:PK(x)d
x式09のごと<、 Pu(u)は形状パラメータ
がa、スケール・パラメータがb/21Aのワイブル分
布となる。したがってUの平均値EI:u〕は。
x式09のごと<、 Pu(u)は形状パラメータ
がa、スケール・パラメータがb/21Aのワイブル分
布となる。したがってUの平均値EI:u〕は。
となる。これは、移動平均回路(3a)出力の推定値と
なる。
なる。
したがって、除算器(5)の出力Qは。
となり、aについて解くと。
1 、、、、、、、、、、、、、、
、 。
、 。
og2 Q
となるので、演算器(6a)の第1の出力Yl及び第2
の出力Y2は次式・で与えられる。
の出力Y2は次式・で与えられる。
Yl= −−γlog2 Q ・・・・・・・
・・・・・・・・■Y2=−=□ ・・・・・・・
・・・・・・・・(至)2 21og2Q 〔発明の効果〕 この発明は以上説明したとおり、入力信号の区間最小値
の平均を用いることによって。
・・・・・・・・■Y2=−=□ ・・・・・・・
・・・・・・・・(至)2 21og2Q 〔発明の効果〕 この発明は以上説明したとおり、入力信号の区間最小値
の平均を用いることによって。
入力の2乗平均平方根を用いる従来の方法よりもハード
ウェア量か減るという効果かめるまた。演算器(6)の
入出力関係を解析的に解り−るため、演算器(6)の設
計か容易になるという効果tろる。
ウェア量か減るという効果かめるまた。演算器(6)の
入出力関係を解析的に解り−るため、演算器(6)の設
計か容易になるという効果tろる。
第1図はこの発明の一実施例を示す図、第2図は移動平
均回路の一構成例を示す図、第3図は最小値検出回路の
一構成例を示1図、第4図は従来の雑音抑圧装置を示す
図である。 図において、(1)は遅延回路、(2)は最小値検出回
路、(3]は移動平均回路、(4)は約数変換回路。 (5)は除算器、(6;は演算器、(7)は加算器、(
8)は減算器、(9)は乗算器、αlは逆対数変換回路
、αDは2乗回路、 aaは平方根回路、09は比較器
、α4はセレクタでるる。 なお2図中同一行号は同一または相当部分を示す。
均回路の一構成例を示す図、第3図は最小値検出回路の
一構成例を示1図、第4図は従来の雑音抑圧装置を示す
図である。 図において、(1)は遅延回路、(2)は最小値検出回
路、(3]は移動平均回路、(4)は約数変換回路。 (5)は除算器、(6;は演算器、(7)は加算器、(
8)は減算器、(9)は乗算器、αlは逆対数変換回路
、αDは2乗回路、 aaは平方根回路、09は比較器
、α4はセレクタでるる。 なお2図中同一行号は同一または相当部分を示す。
Claims (1)
- 入力信号を導入し遅延させる遅延回路と、前記入力信号
と前記遅延回路出力の振幅の小さい方を出力する最小値
検出回路と、前記最小値検出回路出力の移動平均を出力
する第1の移動平均回路と、前記入力信号の移動平均を
出力する第2の移動平均回路と、前記入力信号の対数を
出力する対数変換回路と、前記対数変換回路出力の移動
平均を第1の出力とし、その移動平均幅の2分の1遅延
した前記対数変換回路出力信号を第2の出力として出力
する第3の移動平均回路と、前記第2の移動平均回路出
力を前記第1の移動平均回路出力で除算した商Qを出力
する除算器と、前記除算器出力を導入し第1の出力とし
てγlog_2Q(γはオイラーの定数)を出力し、第
2の出力として1/(2log_2Q)を出力する演算
器と、前記第3の移動平均回路の第1の出力と前記演算
器の第1の出力を加算し出力する加算器と、前記第3の
移動平均回路の第2の出力から前記加算器出力を減算し
出力する減算器と、前記減算器出力と前記演算器の第2
の出力との積を出力する乗算器と、前記乗算器出力の逆
対数を出力する逆対数変換回路を備えたことを特徴とす
る雑音抑圧装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59211118A JPS6189575A (ja) | 1984-10-08 | 1984-10-08 | 雑音抑圧装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59211118A JPS6189575A (ja) | 1984-10-08 | 1984-10-08 | 雑音抑圧装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPS6189575A true JPS6189575A (ja) | 1986-05-07 |
Family
ID=16600697
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP59211118A Pending JPS6189575A (ja) | 1984-10-08 | 1984-10-08 | 雑音抑圧装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS6189575A (ja) |
-
1984
- 1984-10-08 JP JP59211118A patent/JPS6189575A/ja active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5465405A (en) | Apparatus and method for detecting signals | |
| JPH0158466B2 (ja) | ||
| US6065031A (en) | Log converter utilizing offset and method of use thereof | |
| JPS6111380B2 (ja) | ||
| JPS6189575A (ja) | 雑音抑圧装置 | |
| EP0637797A1 (en) | Calculation apparatus | |
| JPS6189576A (ja) | 雑音抑圧装置 | |
| JPS61180165A (ja) | 雑音抑圧装置 | |
| JPS61180166A (ja) | 雑音抑圧装置 | |
| KR910008382B1 (ko) | 제 산 기 | |
| JP2902155B2 (ja) | レーダ用クラッタ信号抑圧装置 | |
| US3548325A (en) | Digital transmission of television | |
| JPH0662284A (ja) | 雑音除去装置 | |
| JP2826678B2 (ja) | ディジタル信号の圧伸方法及びその装置 | |
| JPS61180164A (ja) | 雑音抑圧装置 | |
| US3686489A (en) | Digital slope limiter | |
| RU2013806C1 (ru) | Устройство для определения оптимального периода технического обслуживания изделия | |
| JPS61205878A (ja) | レ−ダ−信号処理装置 | |
| JP2701501B2 (ja) | 平滑化装置 | |
| JPS61219988A (ja) | 映像信号処理回路 | |
| KR900006007B1 (ko) | 디지탈신호 처리회로 | |
| JPH0215835B2 (ja) | ||
| SU1013987A1 (ru) | Устройство дл распознавани случайных сигналов | |
| JP2678072B2 (ja) | レーダ信号処理方法及び装置 | |
| JPS62221212A (ja) | デイジタル信号処理装置 |