JPS5870182A - Processor of signal for supressing clutter - Google Patents
Processor of signal for supressing clutterInfo
- Publication number
- JPS5870182A JPS5870182A JP56168156A JP16815681A JPS5870182A JP S5870182 A JPS5870182 A JP S5870182A JP 56168156 A JP56168156 A JP 56168156A JP 16815681 A JP16815681 A JP 16815681A JP S5870182 A JPS5870182 A JP S5870182A
- Authority
- JP
- Japan
- Prior art keywords
- clutter
- signal
- transmission
- operating circuit
- digital filter
- 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
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/2806—Employing storage or delay devices which preserve the pulse form of the echo signal, e.g. for comparing and combining echoes received during different periods
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
Description
【発明の詳細な説明】
この発明は、レーダで受信される信号のうち海面あるい
は陸で散乱されて受信される信号。DETAILED DESCRIPTION OF THE INVENTION The present invention relates to signals received by radar that are scattered on the sea surface or land.
即ち、クラッタ(CLUテTIR)を抑圧する信号処理
装置に関するものである。That is, the present invention relates to a signal processing device that suppresses clutter (CLU and TIR).
この種の信号処理装置はクラッタKllれた目標信号を
検出することを目的とするもので、処理後の目標信号ピ
ーク電力対クラッタ平均電力(以下これを出力−比と呼
ぶ)を最大にする一種のディジタルフィルタを構成する
ものであるO
このような出力810比を最大にするディジタルフィル
タは、伝達関数!!(f)が次式で今えられるフィルタ
によって得られることは周知の事実となっている。This type of signal processing device is intended to detect a target signal with clutter, and is a type of signal processing device that maximizes the peak power of the target signal after processing to the average power of the clutter (hereinafter referred to as the output-ratio). The digital filter that maximizes the output 810 ratio is the transfer function! ! It is a well-known fact that (f) can be obtained by the filter given by the following equation.
H(f) −R8*(す/P0(リ (1)ここ
に8(つけ上記ディジタルフィルタに入力する目標信号
のスペクトルであ!7 I P、(f)は同ディジタル
フィルタに入力するクラッタの電カスベクトルでToi
RFiOでない任意定数でアシ。H(f) -R8*(S/P0(ri) (1) Here is the spectrum of the target signal input to the above digital filter!7 I P, (f) is the spectrum of the clutter input to the same digital filter. Toi with electric waste vector
Use any constant that is not RFiO.
本は複素共役である。A book is a complex conjugate.
8(f)はレーダから送信される送信信号の数学的表現
から計算される量で、既知量である。即ち、送信信号−
、!(1)が第(2)式で与えられるとき第(萄式で計
算される
ここに。8(f) is a quantity calculated from a mathematical expression of a transmission signal transmitted from the radar, and is a known quantity. That is, the transmitted signal -
,! When (1) is given by equation (2), the equation (2) is calculated here.
bn=送信パルスの複素振幅 M:送信パルスのm& T、:送信パルス間隔 fo:送信周波数 である。bn=complex amplitude of transmitted pulse M: m& of transmitted pulse T: Transmission pulse interval fo: Transmission frequency It is.
一方、クラッタの電カスベクトルp、<rtは、クラッ
タの性質によって変化するものであ〕、一般には未知量
であシ、tた。Po(つを正確に知るために鉱1通常、
多くの観測データを必要とする。そのために従来では式
(重)のようなディジタルフィルタは、クラッタの性質
がよ(知られておシ、かつ、その変化の小さい場合に適
用されてお〕、上述の条件が満えされないクラッタに対
して紘全く適用されていなかった。On the other hand, the electric sludge vector p,<rt of the clutter changes depending on the properties of the clutter, and is generally an unknown quantity. Po
Requires a lot of observation data. For this reason, conventionally, digital filters such as the formula (heavy) have been applied to cases where the properties of clutter are well known and their changes are small, and where the above conditions are not met. On the other hand, Hiro was not applied at all.
この発明鉱、クラッタの電カスベクトルを受信信号から
正確に推定できる周知のスペクトル推定法であるMIM
(Maximum Intropy M@thod)
の演算アルゴリズムを実施する演算回路を用い。This invention is based on MIM, a well-known spectrum estimation method that can accurately estimate the electric scum vector of clutter from the received signal.
(Maximum Intropy M@thod)
using an arithmetic circuit that implements the arithmetic algorithm.
第(1)式に示したディジタルフィルタが、予め性質の
全くわからないクラッタに対しても効果があるようにす
ることを目的とする。以下1図面について詳細に説明す
る。It is an object of the present invention to make the digital filter shown in equation (1) effective even for clutter whose properties are completely unknown in advance. One drawing will be explained in detail below.
第1図は、この発明の一!施例であって1図中、(1)
はクラッタの電カスベクトルを計算するM、KM演算回
路、(2)はディジタルフィルタ、(3)はメモリ、(
4)はフィルタ係数演算回路、(団はバッファである。Figure 1 is one of the inventions! As an example, (1) in Figure 1
(2) is a digital filter, (3) is a memory, (
4) is a filter coefficient calculation circuit (the group is a buffer).
レーダ受信信号y(nT)(””#1#2t・・・M−
1)は、tず9M!!輩演算回路(1)およびバッファ
(鴎に入力される。バッファ(5)に入力された受信信
号y (nT。)は後に述べるフィルタ係数0・の計算
が完了するまで一時蓄えられる。一方。Radar received signal y (nT) (""#1#2t...M-
1) is tzu9M! ! The received signal y (nT.) input to the buffer (5) is input to the filter calculation circuit (1) and the buffer (Oro) and is temporarily stored until the calculation of the filter coefficient 0, which will be described later, is completed.
■演算回゛路(1)は受信信号y(n’r−の電カスベ
クトルP、(イ)を次式の形で推定する。(2) The calculation circuit (1) estimates the electric waste vector P, (a) of the received signal y(n'r-) in the form of the following equation.
MIMはパラメータM 、 PM 、 ai (L ■
1 g 2 # −” IM)をr(my−(nm(1
# 1 # 2 e−tw−1)から計算し、P−を推
定する会知のアルゴリズムである。M]:Mは、不規則
変動する受信信号r (nT、)に第(6)式に示す自
己回帰式を適用しe Ms ’Msai(1−1,2s
”・、M)を計算するアルゴリズムであって。MIM has parameters M , PM , ai (L ■
1 g 2 # −”IM) as r(my−(nm(1
#1 #2 e-tw-1) and estimates P-. M]: M applies the autoregressive equation shown in Equation (6) to the irregularly fluctuating received signal r (nT,).
”・,M).
r (”s)−JT ((”−1)テ、) −a、r(
(n−2)丁、)■−−−−aMr((n−M)? )
+e(nT、)・(nT、) :白色ガラス雑音
第(6)式がsr(”)の予測値y(nT、)を*ff
)式の回帰予測式で計算し。r (”s)−JT ((”−1)te,) −a, r(
(n-2) Ding,)■----aMr((n-M)? )
+e(nT,)・(nT,): White glass noise Equation (6) sets the predicted value y(nT,) of sr('') as *ff
) is calculated using the regression prediction formula.
r(nT、) −−J a、g(n−m)?、)
(7)m■1
その予測誤差・Ca?、)−y(nT、)−y(nT、
)を白色ガラス雑音にするための一種の数値フィルタと
解釈できるためeMti数値フィルタのオーダ。r(nT,) --J a,g(n-m)? ,)
(7) m■1 The prediction error・Ca? , )-y(nT,)-y(nT,
) can be interpreted as a kind of numerical filter to make white glass noise, so the order of eMti numerical filter.
al(t−t H2# ”・s M)は数値フィルタ係
数。al(t-t H2#”・s M) is a numerical filter coefficient.
PMは予測誤差・(nT s )の分散と呼ばれる。PM is called the variance of prediction error/(nTs).
MIM演算回路(1)はこのような公知のアルゴリズム
を実施する演算回路である。一般にy(nT、)にはク
ラッタ0(nT )と目標信号−(11?、)が會まれ
でいるが、第(8)式の関係が成立するから。The MIM arithmetic circuit (1) is an arithmetic circuit that implements such a known algorithm. Generally, clutter 0(nT) and target signal -(11?,) are combined in y(nT,), and the relationship in equation (8) holds true.
MIM演算回路(11によってクラッタの電力スペク)
/L/ P、(f)が精度よく計算される。MIM calculation circuit (clutter power spec by 11)
/L/P, (f) is calculated with high accuracy.
l C(nT、)l>> I I(nT、)1
1g)さて、クラッタの電力スペクトルア。(
f′)が計Xできれば、ディジタルフィルタ(2)の伝
達関数H(f)蝶、第(3)式、第(5)式を第(1)
式に代入して次式で与えられる。l C(nT,)l>> I I(nT,)1
1g) Now, the power spectrum of clutter. (
If f') can be summed up to
Substituting into the equation, it is given by the following equation.
第(9)式に示すような伝達関数H(f)をもつディジ
タルフィルタは第2図に示すようにg−zm+1段のト
ランスバーサルディジタルフィルタ(テransver
sal digital Filt@r )によって実
現されることは公知の事実である。第2図において(6
)は時間Tの遅延素子、 ff)は複素乗算器、(a)
は複素加算器である。フィルタ係数0.Ct−o。A digital filter with a transfer function H(f) as shown in equation (9) is a g-zm+1 stage transversal digital filter (transversal digital filter) as shown in FIG.
It is a well-known fact that this is realized by digital filter@r). In Figure 2 (6
) is a delay element of time T, ff) is a complex multiplier, (a)
is a complex adder. Filter coefficient 0. Ct-o.
1、・・・、L−1)はMIM演算回路(1)から転送
されるパラメータMおよびai(1ml、−・・t M
)と送信パルス複素重みbn(nm(1、−、M−1)
とを入力とするフィルタ係数演算回路(4によ)計算さ
れ、その結果はメモリ(3)に転送され、記憶される。1, ..., L-1) are the parameters M and ai (1ml, ..., t M
) and the transmitted pulse complex weight bn (nm (1, -, M-1)
is calculated by the filter coefficient calculation circuit (4) which receives as input, and the result is transferred to and stored in the memory (3).
フィルタ係数演算回路(4では次式に示すような演算が
行なわれてフィルタ係数at(Z−・。In the filter coefficient calculation circuit (4), calculations as shown in the following equation are performed and the filter coefficient at(Z-.
1.2.・・・、L−1)が計算される@N−4−M−
$43
但し、ao−1
’p” @ # P >M or p ) 0bp
−O12>N Orp<。1.2. ..., L-1) is calculated @N-4-M-
$43 However, ao-1 'p" @ # P > M or p ) 0bp
-O12>N Orp<.
である。It is.
以上のように2本発明によるクラッタ抑圧用信号処理装
置ではクラッタの性質に応じて、アダプティブにその伝
達関数をかえることができるから、クラッタの性質およ
びその変化にかかわらず常に出力84比を最大にするこ
と゛ができるO
以上述べたように、この発明によれば、クラッタの性質
に関係なくクヂツタの抑圧が可能でIh!>、クラッタ
に埋れた目標を検出する装置に用いて、その効果拡大き
い。As described above, the signal processing device for clutter suppression according to the present invention can adaptively change its transfer function according to the characteristics of clutter, so the output 84 ratio can always be maximized regardless of the characteristics of clutter and its changes. As described above, according to the present invention, it is possible to suppress ivy regardless of the nature of the clutter. >It can be used in a device to detect targets buried in clutter, increasing its effectiveness.
第1図はこの発明の一実施例を示す図、!2図はトラン
スバーサルディジタルフィルタの詳細図である。
図中、(1)はM1eM演算回路、(りdディジタルフ
ィルタ、(3)はメモリ、(菊はフィルタ係数演算回路
、(51Fiバツフア、(6)は遅延素子、(7)は複
素乗算器、(8)は複素加算器である〇
代理人 弁理士葛野信−FIG. 1 is a diagram showing an embodiment of this invention. FIG. 2 is a detailed diagram of the transversal digital filter. In the figure, (1) is an M1eM calculation circuit, (digital filter), (3) is a memory, (chrysanthemum is a filter coefficient calculation circuit, (51Fi buffer), (6) is a delay element, (7) is a complex multiplier, (8) is a complex adder 〇 Agent Patent attorney Makoto Kuzuno -
Claims (1)
れて受信されるクラッタ(ahaテテU)を抑圧するク
ラッタ抑圧用信号処理装置において* bB(”” s
1 a 2 m”” 1 )を送信パルスの複素振
幅1Mを送信パルス数、!、をパルス送信間隔、foを
送信周波数、tを送信パルス幅、 rect(t)を で定義される関終として。 の数学的表現式で表わされる夏個のパルス電液倉ら構成
される送信信号を送信した結果得られるN個のクラッタ
のうちレーダから同一距離にある成分を抽出して得られ
るN個の複素ディジタル信号の電カスベクトル1’、(
f)t−で与えられる数学的表現式で表わしたときに。 輩、 PM、 ai(im*1 、2、−1M)を計算
−j ルM]liM(Maximum ]!:ntro
py Method )アルゴリズムを実行する演算装
置と、上記N個の複素ディジタル信号を入力信号とし、
伝達−数B(f)がRを任意定数、*を複素共役として H(f)−11+ X aiexp(−コIgifT
)121−1′ であるトランスバーサルディジタルフィルタ(τran
−マersal Digital Filter )と
を具備したことを特徴とするクラッタ抑圧用信号処理装
置。[Claims] In a signal processing device for clutter suppression that suppresses clutter (aha tete U) that is scattered and received on the sea surface or land out of the signal O received by a radar, *bB("" s
1 a 2 m”” 1 ) is the complex amplitude of the transmitted pulse 1M is the number of transmitted pulses,! , where is the pulse transmission interval, fo is the transmission frequency, t is the transmission pulse width, and rect(t) is the endpoint defined by. N complex clutter obtained by extracting components at the same distance from the radar from among N clutter obtained as a result of transmitting a transmission signal consisting of a number of pulsed electric liquids represented by the mathematical expression Electric waste vector 1' of digital signal, (
f) When expressed in the mathematical expression given by t-. Calculate PM, ai (im*1, 2, -1M) -j leM]liM (Maximum]!: ntro
py Method) an arithmetic unit that executes the algorithm, and the above N complex digital signals as input signals,
The transmission number B(f) is H(f)-11+X aiexp(-koIgifT
)121-1' is a transversal digital filter (τran
1. A signal processing device for suppressing clutter, comprising: - a digital filter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP56168156A JPS5870182A (en) | 1981-10-21 | 1981-10-21 | Processor of signal for supressing clutter |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP56168156A JPS5870182A (en) | 1981-10-21 | 1981-10-21 | Processor of signal for supressing clutter |
Publications (2)
Publication Number | Publication Date |
---|---|
JPS5870182A true JPS5870182A (en) | 1983-04-26 |
JPH0318156B2 JPH0318156B2 (en) | 1991-03-11 |
Family
ID=15862846
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP56168156A Granted JPS5870182A (en) | 1981-10-21 | 1981-10-21 | Processor of signal for supressing clutter |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPS5870182A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4680589A (en) * | 1984-10-02 | 1987-07-14 | The United States Of America As Represented By The Secretary Of The Air Force | Adaptive fast Fourier transform weighting technique to increase small target sensitivity |
JPS6443780A (en) * | 1987-08-11 | 1989-02-16 | Mitsubishi Electric Corp | Clutch suppressing device |
JPH07146350A (en) * | 1994-08-02 | 1995-06-06 | Nec Corp | Method for orienting azimuth of sound source |
WO1998021602A1 (en) * | 1996-11-08 | 1998-05-22 | Bayerische Motoren Werke Aktiengesellschaft | Process for determining the distance between two objects |
JP2011242288A (en) * | 2010-05-19 | 2011-12-01 | Honda Elesys Co Ltd | Electronic scanning radar device, reception wave direction estimation method, and reception wave direction estimation program |
CN105425218A (en) * | 2015-12-31 | 2016-03-23 | 电子科技大学 | Radar-communication integrated realizing method |
RU2819294C1 (en) * | 2023-07-07 | 2024-05-16 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Рязанский государственный радиотехнический университет имени В.Ф. Уткина" | Interference suppression computer |
-
1981
- 1981-10-21 JP JP56168156A patent/JPS5870182A/en active Granted
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4680589A (en) * | 1984-10-02 | 1987-07-14 | The United States Of America As Represented By The Secretary Of The Air Force | Adaptive fast Fourier transform weighting technique to increase small target sensitivity |
JPS6443780A (en) * | 1987-08-11 | 1989-02-16 | Mitsubishi Electric Corp | Clutch suppressing device |
JPH07146350A (en) * | 1994-08-02 | 1995-06-06 | Nec Corp | Method for orienting azimuth of sound source |
WO1998021602A1 (en) * | 1996-11-08 | 1998-05-22 | Bayerische Motoren Werke Aktiengesellschaft | Process for determining the distance between two objects |
JP2011242288A (en) * | 2010-05-19 | 2011-12-01 | Honda Elesys Co Ltd | Electronic scanning radar device, reception wave direction estimation method, and reception wave direction estimation program |
CN105425218A (en) * | 2015-12-31 | 2016-03-23 | 电子科技大学 | Radar-communication integrated realizing method |
RU2819294C1 (en) * | 2023-07-07 | 2024-05-16 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Рязанский государственный радиотехнический университет имени В.Ф. Уткина" | Interference suppression computer |
Also Published As
Publication number | Publication date |
---|---|
JPH0318156B2 (en) | 1991-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5900835A (en) | Coherent hidden markov model | |
US4783660A (en) | Signal source distortion compensator | |
Frost et al. | An adaptive filter for smoothing noisy radar images | |
Kaveh | High resolution spectral estimation for noisy signals | |
Zedel et al. | Resolving velocity ambiguity in multifrequency, pulse-to-pulse coherent Doppler sonar | |
EP0149981A2 (en) | Adaptive radar signal processor for the detection of the useful echo and the cancellation of clutter | |
JP2005530164A (en) | Method for suppressing interference in an object detection system | |
JPH0527287B2 (en) | ||
JPS5870182A (en) | Processor of signal for supressing clutter | |
Papoulis et al. | Improvement of range resolution by spectral extrapolation | |
JPH032433B2 (en) | ||
Burkhardt et al. | A modified Newton-Raphson search for the model-adaptive identification of delays | |
Lee et al. | Evaluation of reduced-rank, adaptive matched field processing algorithms for passive sonar detection in a shallow-water environment | |
CN109143217B (en) | Distance folding echo suppression method based on MISO system | |
JPS62132186A (en) | Radar signal processor | |
US4916453A (en) | Spatial filtering system | |
JP4109918B2 (en) | Clutter suppression apparatus and method | |
JP2513784B2 (en) | Radar device | |
JPS59191608A (en) | Improvement in processing through put speed | |
SU693342A1 (en) | Device for identification of linear dynamic objects | |
CA1283200C (en) | Correlator circuit and device for ultrasound flow mapping comprising such as circuit | |
RU2191405C1 (en) | Procedure determining radial velocity of object | |
SU1753601A1 (en) | Device for compensation of noises from closely placed transmitter | |
Gawrys et al. | Divergence and the fading memory filter | |
SU839024A2 (en) | Random pulse generator |