CN116047488A - Active sonar interference suppression method for pressed interference source - Google Patents
Active sonar interference suppression method for pressed interference source Download PDFInfo
- Publication number
- CN116047488A CN116047488A CN202211473519.2A CN202211473519A CN116047488A CN 116047488 A CN116047488 A CN 116047488A CN 202211473519 A CN202211473519 A CN 202211473519A CN 116047488 A CN116047488 A CN 116047488A
- Authority
- CN
- China
- Prior art keywords
- interference
- theta
- azimuth
- signal
- pressed
- 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
- 230000001629 suppression Effects 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000000737 periodic effect Effects 0.000 claims abstract description 20
- 238000012545 processing Methods 0.000 claims abstract description 20
- 238000001514 detection method Methods 0.000 claims abstract description 19
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 238000000605 extraction Methods 0.000 claims description 10
- 230000006835 compression Effects 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 9
- 230000003044 adaptive effect Effects 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000002452 interceptive effect Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000012216 screening Methods 0.000 claims description 6
- 230000021615 conjugation Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
Images
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/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/537—Counter-measures or counter-counter-measures, e.g. jamming, anti-jamming
-
- 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/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses an active sonar interference suppression method for a pressed interference source, and relates to the field of acoustic signal processing target identification. The method utilizes the repeated periodicity of the suppressed interference to extract the periodic interference signals in the interference azimuth, and uses the extracted periodic interference signals as reference signals to carry out self-adaptive filtering on the beam signals in all the azimuth in the interference range, thereby realizing the interference suppression of the full-space multi-beam signals. The method can effectively inhibit periodic suppression type interference, greatly improves the active target detection rate under the interference background, and can support the improvement of the anti-interference capability of sonar equipment.
Description
Technical Field
The invention relates to the field of acoustic signal processing target identification, in particular to an active sonar interference suppression method aiming at a pressed interference source.
Background
The underwater sound active target detection technology is one of the main technical means for underwater target detection in various countries at present, and various countries successively develop underwater countermeasure equipment with rich types so as to weaken or even inhibit the underwater sound target detection efficiency, so that the underwater countermeasure situation is increasingly complex. The method is characterized in that a high-power interference signal is emitted by a pressing type interference device to press the azimuth of the target, which is one of typical underwater countermeasure technical means, so that the azimuth detection background of the target can be greatly improved, and the effect of 'submerging' the target is achieved. Aiming at the difficult problem of sonar target detection under the pressed interference, the existing interference suppression methods such as self-adaptive beam forming, spatial filtering and the like are adopted, and the time-frequency characteristics of the interference are not considered by utilizing the space azimuth difference of the interference and the target.
The common suppression type interference mode is aiming frequency interference and noise interference, wherein the aiming frequency interference is mainly used for suppressing an active target echo signal, and the active sonar mostly utilizes pulse signals to detect, so that the aiming frequency interference signal has the property of a periodic pulse signal. The noise interference adopts broadband white noise as an interference signal, and can be used for suppressing an active target echo signal and a passive target radiation noise signal at the same time, and the current part of model noise interference equipment is used for generating a continuous suppressing interference effect by utilizing a cyclic playing mode by recording a section of white noise signal in advance when suppressing is needed. The repetitive period characteristic of the suppression interference is determined by the design mode of the suppression interference equipment, the periodic characteristic of the target under the active sonar detection is based on the active sonar emission period, and the two are not consistent, so that the suppression interference can be suppressed by utilizing the characteristic difference. Different from the aforementioned adaptive beam forming, spatial filtering and other spatial interference suppression methods, the interference suppression method based on the time-frequency difference characteristics of the interference and the target signal can suppress the interference when the interference and the target azimuth overlap each other, and acquire the high signal-to-noise ratio active target echo signal under the main lobe interference.
Disclosure of Invention
The invention provides an active sonar interference suppression method for a pressed interference source aiming at the problem of difficult detection of an active sonar target under the pressed interference, and the method can effectively suppress periodic pressed interference, greatly improve the detection rate of the active target under the interference background and improve the anti-interference capability of supportable sonar equipment.
In order to achieve the above purpose, the present invention provides the following technical solutions: an active sonar interference suppression method for a pressed interference source comprises the following steps:
s1: processing the matrix data to obtain multi-beam data covering all spatial orientations;
s2: determining an estimated value B (p, t) of the suppressed interference main lobe beam signal;
s3: intercepting at least one full period of suppressed interfering signal
S4: for a pair of
Performing autocorrelation processing to obtain an autocorrelation envelope X (tau) with a plurality of correlation peaks;
s5: peak value extraction is carried out on the correlation peak in the S4, and finally the suppression type interference cycle period T is obtained 2 ；
S6: determining the impact range theta of the piezoelectric interference I ；
S7: extracting periodic interference for all beams within the suppressed interference impact range
S8: to be used for
Performing adaptive filtering on each wave beam in the current active transmission period in the interference influence range for the reference signal to realize interference suppression;
s9: and taking the transmitted signal as a reference, and carrying out matching processing on the result after the error signal is filtered. The invention is further configured that the S1 is specifically as follows: processing the matrix data by conventional beam forming method to obtain multi-beam data covering all spatial orientations as follows
B(θ，t)＝S 0 (θ，t)+S I (θ，t)+N(θ，t)
Wherein B (θ, t) represents time-domain beam data at time of θ azimuth t, S 0 (θ, t) is a leakage component of the target signal at the time of θ azimuth t, S I (theta, t) is the time of the interference signal at the theta azimuth tN (θ, t) is the background noise at time θ azimuth t.
The invention further provides that the step S2 is specifically as follows: determining the azimuth of the pressed interference main lobe beam, and detecting the energy of the whole azimuth of B (theta, t)
Wherein E (θ) is the θ azimuth beam energy; searching for the azimuth p based on the suppressed interference characteristics to satisfy
E(p)＝max(E(θ))
The beam signal of the p azimuth is
B(p，t)＝S 0 (p，t)+S I (p，t)+N(p，t)≈S I (t)
Wherein S is I (t) is time domain data of the original interference signal at the time t; due to S I The energy intensity of (p, t) is much greater than S 0 (p, t) and N (p, t), and S I (p, t) is an interference signal S I An approximate estimate of (t), so B (p, t) can be taken as S I An estimate of (t).
The invention is further configured that the S3 specifically includes: estimating periodicity of the interference by B (p, T), and intercepting the time length T in B (p, T) 1 A section of arbitrary data B T1 (p，t)，T 1 It is required to be long enough to be consistent with the length of B (p, t) to ensure that at least one complete cycle of interference data is contained.
The invention further provides that the step S4 is specifically as follows: for a pair of
Performing autocorrelation processing to obtain
Wherein X (τ) is the autocorrelation envelope under the condition of delay τ; due to
There is at least one period of interfering signal and thus there are a plurality of correlation peaks in X (τ).
The invention is further configured that the S5 specifically includes: peak extraction of correlation peaks due to
The dry-to-noise ratio of each point X (tau) is extremely high, so that the peak value extraction is carried out by utilizing a method of threshold value screening and maximum value adding, and the dry-to-noise ratio of each point X (tau) is calculated
INR(τ)＝20lg|X(τ)/X b |
Where INR (τ) represents the dry-to-noise ratio of each point of X (τ); x is X b Is an autocorrelation envelope background, which can be obtained by carrying out a sliding average on X (tau), and carrying out threshold value screening and maximum value searching on INR (tau) to obtain a sequence Q= [ Q ] of the positions of a group of correlation peaks 1 ，q 2 ，...，q k ]Then
T 2 ＝q i+1 -q i
Wherein T is 2 The compression type interference cycle period is represented, i is more than or equal to 1 and less than k.
The invention is further configured that the step S6 is specifically as follows: determining the influence range of the compression type interference, and calculating the omnibearing dry-to-noise ratio by combining E (theta) obtained in the step 2
I(θ)＝20lg(E(θ)/E 0 (θ))
Wherein, I (theta) is the theta azimuth dry-to-noise ratio, E 0 (theta) is the theta azimuth background noise energy before the interference appears, and the influence range of the suppressed interference is theta I 0 is then I Needs to meet the requirements of
I(θ I )≥δ
Where δ is a threshold of the interference-to-noise ratio, and δ is generally set to be more than 10dB in order to prevent the interference suppression effect from being affected by too small interference-to-noise ratio.
The invention is further configured that the S7 is specifically as follows: extracting periodic interference from all wave beams in the influence range of the suppressed interference; to avoid the influence of active emission factors such as direct wave, reverberation and the like, T before active detection pulse emission is intercepted 2 Time period beam data as periodic interference, expressed as
The invention is further configured that the S8 specifically includes: to be used for
Performing adaptive filtering on each wave beam in the current active transmission period in the interference influence range for the reference signal to realize interference suppression; a least mean square error (LMS) algorithm is used here, a transversal filter structure is used;
let theta in the nth iteration process I The directional filter weight coefficient vector is w (θ) I N), then n+1th iteration
w(θ I ，n+1)＝w(θ I ，n)+μ(θ I )u(θ I ，n)(d * (n)-u H (θ I ，n)w(θ I ，n))
Wherein μ (θ) I ) For theta I Direction step, u (θ) I N) is theta I Azimuth beam input sampling, u H (θ I N) is u (theta) I Transpose of n), d * (n) is
Conjugation of the samples; the error signal can be expressed as
e(θ I ，n)＝d(n)-w T (θ I ，n)u(θ I ，n)
Wherein e (θ) I N) is theta I The error of the nth iteration of the direction comprises filtering the target echo signal after the compression type interference.
The invention is further configured that the S9 specifically includes: with reference to the transmitted signal, pair e (θ I N) performing matching processing
Wherein Z (θ) I τ) is θ I The matching envelope samples of the directional delay tau,
to sample the transmitted signal S L (n) performing delay conjugation; for each azimuth Z (theta) I τ), and can obtain azimuth distance information about the target; since the interference is suppressed at this time, the target echo detection performance is better than before the interference suppression.
In general, the above technical solutions conceived by the present invention, compared with the prior art, enable the following beneficial effects to be obtained:
the method can effectively inhibit periodic suppression type interference, greatly improve the active target detection rate under the interference background, and improve the anti-interference capability of the supportable sonar equipment.
Drawings
Fig. 1 is a schematic diagram of the interference suppression principle of periodic suppressed interference;
FIG. 2 is a block diagram of the technical scheme of the invention;
FIG. 3 is a schematic diagram of periodic disturbance extraction;
FIG. 4 is a graph of the active target detection effect prior to periodic hold-down interference suppression;
FIG. 5 is a graph of the active target detection effect after periodic hold-down interference suppression;
Detailed Description
The following detailed description of the invention, taken in conjunction with fig. 1-5, will provide those skilled in the art with a more clear understanding of how to practice the invention. While the present invention has been described in connection with the preferred embodiments thereof, these embodiments are set forth only and are not intended to limit the scope of the invention.
As shown in fig. 1 to 3, an active sonar interference suppression method for a pressed interference source includes the following steps:
s1: processing the matrix data to obtain multi-beam data covering all spatial orientations;
s2: determining an estimated value B (p, t) of the suppressed interference main lobe beam signal;
s3: intercepting at least one full period of suppressed interfering signal
S4: for a pair of
Performing autocorrelation processing to obtain an autocorrelation envelope X (tau) with a plurality of correlation peaks;
s5: peak value extraction is carried out on the correlation peak in the S4, and finally the suppression type interference cycle period T is obtained 2 ；
S6: determining the impact range theta of the piezoelectric interference I ；
S7: extracting periodic interference for all beams within the suppressed interference impact range
S8: to be used for
Performing adaptive filtering on each wave beam in the current active transmission period in the interference influence range for the reference signal to realize interference suppression;
s9: and taking the transmitted signal as a reference, and carrying out matching processing on the result after the error signal is filtered. The S1 specifically comprises the following steps: processing the matrix data by conventional beam forming method to obtain multi-beam data covering all spatial orientations as follows
B(θ，t)＝S 0 (θ，t)+S I (θ，t)+N(θ，t)
Wherein B (θ, t) represents time-domain beam data at time of θ azimuth t, S 0 (θ, t) is a leakage component of the target signal at the time of θ azimuth t, S I (θ, t) is a leakage component of the interference signal at the θ azimuth t, and N (θ, t) is background noise at the θ azimuth t.
The step S2 is specifically as follows: determining the azimuth of the pressed interference main lobe beam, and detecting the energy of the whole azimuth of B (theta, t)
Wherein E (θ) is the θ azimuth beam energy; searching for the azimuth p based on the suppressed interference characteristics to satisfy
E(p)＝max (E(θ))
The beam signal of the p azimuth is
B(p，t)＝S 0 (p，t)+S I (p，t)+N(p，t)≈S I (t)
Wherein S is I (t) is time domain data of the original interference signal at the time t; due to S I The energy intensity of (p, t) is much greater than S 0 (p, t) and N (p, t), and S I (p, t) is an interference signal S I An approximate estimate of (t), so B (p, t) can be taken as S I An estimate of (t).
The step S3 is specifically as follows: estimating periodicity of the interference by B (p, T), and intercepting the time length T in B (p, T) 1 Is a piece of arbitrary data of (1)
T 1 It is required to be long enough to be consistent with the length of B (p, t) to ensure that at least one complete cycle of interference data is contained. The step S4 is specifically as follows: for->
Performing autocorrelation processing to obtain
Wherein X (τ) is the autocorrelation envelope under the condition of delay τ; due to
There is at least one period of interfering signal and thus there are a plurality of correlation peaks in X (τ).
The S5 toolThe body is as follows: peak extraction of correlation peaks due to
The dry-to-noise ratio of each point X (tau) is extremely high, so that the peak value extraction is carried out by utilizing a method of threshold value screening and maximum value adding, and the dry-to-noise ratio of each point X (tau) is calculated
INR(τ)＝20lg|X(τ)/X b |
Where INR (τ) represents the dry-to-noise ratio of each point of X (τ); x is X b Is an autocorrelation envelope background, which can be obtained by carrying out a sliding average on X (tau), and carrying out threshold value screening and maximum value searching on INR (tau) to obtain a sequence Q= [ Q ] of the positions of a group of correlation peaks 1 ，q 2 ，...，q k ]Then
T 2 ＝q i+1 -q i
Wherein T is 2 The compression type interference cycle period is represented, i is more than or equal to 1 and less than k.
The step S6 is specifically as follows: determining the influence range of the compression type interference, and calculating the omnibearing dry-to-noise ratio by combining E (theta) obtained in the step 2
I(θ)＝20lg(E(θ)/E 0 (θ))
Wherein, I (theta) is the theta azimuth dry-to-noise ratio, E 0 (theta) is the theta azimuth background noise energy before the interference appears, and the influence range of the suppressed interference is theta I Theta is then I Needs to meet the requirements of
I(θ I )≥δ
Where δ is a threshold of the interference-to-noise ratio, and δ is generally set to be more than 10dB in order to prevent the interference suppression effect from being affected by too small interference-to-noise ratio.
The step S7 is specifically as follows: extracting periodic interference from all wave beams in the influence range of the suppressed interference; to avoid the influence of active emission factors such as direct wave, reverberation and the like, T before active detection pulse emission is intercepted 2 Time period beam data as periodic interference, expressed as
The step S8 is specifically as follows: to be used for
Performing adaptive filtering on each wave beam in the current active transmission period in the interference influence range for the reference signal to realize interference suppression; a least mean square error (LMS) algorithm is used here, a transversal filter structure is used;
let theta in the nth iteration process I The directional filter weight coefficient vector is w (θ) I N), then n+1th iteration
w(θ I ，n+1)＝w(θ I ，n)+μ(θ I )u(θ I ，n)(d * (n)-u H (θ I ，n)w(θ I ，n))
Wherein μ (θ) I ) For theta I Direction step, u (θ) I N) is theta I Azimuth beam input sampling, u H (θ I N) is u (theta) I Transpose of n), d * (n) is
Conjugation of the samples; the error signal can be expressed as
e(θ I ，n)＝d(n)-w T (θ I ，n)u(θ I ，n)
Wherein e (θ) I N) is theta I The error of the nth iteration of the direction comprises filtering the target echo signal after the compression type interference.
The step S9 is specifically as follows: with reference to the transmitted signal, pair e (θ I N) performing matching processing
Wherein Z (θ) I τ) is θ I The matching envelope samples of the directional delay tau,
to sample the transmitted signal S L (n) performing delay conjugation; for each azimuth Z (theta) I τ) to perform target detection,azimuth distance information about the target can be obtained; since the interference is suppressed at this time, the target echo detection performance is better than before the interference suppression.
The effect is shown in fig. 4 and 5, in which the interference signal before interference suppression completely suppresses the target echo signal, the target cannot be detected, the interference is basically eliminated after interference suppression, and the target echo is clearly visible.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the present invention, and various modifications and variations may be made to the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (10)
1. An active sonar interference suppression method for a pressed interference source is characterized by comprising the following steps:
s1: processing the matrix data to obtain multi-beam data covering all spatial orientations;
s2: determining an estimated value B (p, t) of the suppressed interference main lobe beam signal;
s3: intercepting at least one full period of suppressed interfering signal
S4: for a pair of
Performing autocorrelation processing to obtain an autocorrelation envelope X (tau) with a plurality of correlation peaks;
s5: peak value extraction is carried out on the correlation peak in the S4, and finally the suppression type interference cycle period T is obtained 2 ；
S6: determining the impact range theta of the piezoelectric interference I ；
S7: extracting periodic interference for all beams within the suppressed interference impact range
S8: to be used for
Performing adaptive filtering on each wave beam in the current active transmission period in the interference influence range for the reference signal to realize interference suppression;
s9: and taking the transmitted signal as a reference, and carrying out matching processing on the result after the error signal is filtered.
2. The active sonar interference rejection method for a pressed interference source of claim 1 wherein S1 is specifically as follows: processing the matrix data by conventional beam forming method to obtain multi-beam data covering all spatial orientations as follows
B(θ,t)＝S 0 (θ,t)+S I (θ,t)+N(θ,t)
Wherein B (θ, t) represents time-domain beam data at time of θ azimuth t, S 0 (θ, t) is a leakage component of the target signal at the time of θ azimuth t, S I (θ, t) is a leakage component of the interference signal at the θ azimuth t, and N (θ, t) is background noise at the θ azimuth t.
3. The active sonar interference rejection method for a pressed interference source of claim 1 wherein S2 is specifically as follows: determining the azimuth of the pressed interference main lobe beam, and detecting the energy of the whole azimuth of B (theta, t)
Wherein E (θ) is the θ azimuth beam energy; searching for the azimuth p based on the suppressed interference characteristics to satisfy
E(p)＝max(E(θ))
The beam signal of the p azimuth is
B(p,t)＝S 0 (p,t)+S I (p,t)+N(p,t)≈S I (t)
Wherein S is I (t) is time domain data of the original interference signal at the time t; due to S I The energy intensity of (p, t) is much greater than S 0 (p, t) and N (p, t), and S I (p, t) is an interference signal S I An approximate estimate of (t), so B (p, t) can be taken as S I An estimate of (t).
4. The active sonar interference rejection method for a pressed interference source of claim 1 wherein S3 is specifically as follows: estimating periodicity of the interference by B (p, T), and intercepting the time length T in B (p, T) 1 Is a piece of arbitrary data of (1)
T 1 It is required to be long enough to be consistent with the length of B (p, t) to ensure that at least one complete cycle of interference data is contained.
5. The active sonar interference rejection method for a pressed interference source of claim 1 wherein S4 is specifically as follows: for a pair of
Performing autocorrelation processing to obtain
Wherein X (τ) is the autocorrelation envelope under the condition of delay τ; due to
There is at least one period of interfering signal and thus there are a plurality of correlation peaks in X (τ).
6. The active sonar interference rejection method for a pressed interference source of claim 1 wherein S5 is specifically as follows:peak extraction of correlation peaks due to
The dry-to-noise ratio of each point X (tau) is extremely high, so that the peak value extraction is carried out by utilizing a method of threshold value screening and maximum value adding, and the dry-to-noise ratio of each point X (tau) is calculated
INR(τ)＝20lg|X(τ)/X b |
Where INR (τ) represents the dry-to-noise ratio of each point of X (τ); x is X b Is an autocorrelation envelope background, which can be obtained by carrying out a sliding average on X (tau), and carrying out threshold value screening and maximum value searching on INR (tau) to obtain a sequence Q= [ Q ] of the positions of a group of correlation peaks 1 ,q 2 ,…,q k ]Then
T 2 ＝q i+1 -q i
Wherein T is 2 Representing the compression type interference cycle period, i is more than or equal to 1<k。
7. The active sonar interference rejection method for a pressed interference source of claim 1 wherein S6 is specifically as follows: determining the influence range of the compression type interference, and calculating the omnibearing dry-to-noise ratio by combining E (theta) obtained in the step 2
I(θ)＝20lg(E(θ)/E 0 (θ))
Wherein, I (theta) is the theta azimuth dry-to-noise ratio, E 0 (theta) is the theta azimuth background noise energy before the interference appears, and the influence range of the suppressed interference is theta I Theta is then I Needs to meet the requirements of
I(θ I )≥δ
Where δ is a threshold of the interference-to-noise ratio, and δ is generally set to be more than 10dB in order to prevent the interference suppression effect from being affected by too small interference-to-noise ratio.
8. The active sonar interference rejection method for a pressed interference source of claim 1 wherein S7 is specifically as follows: extracting periodic interference from all wave beams in the influence range of the suppressed interference; to avoid the influence of active emission factors such as direct wave, reverberation and the like, T before active detection pulse emission is intercepted 2 Time period beam data as periodic interference, expressed as
9. The active sonar interference rejection method for a pressed interference source of claim 1 wherein S8 is specifically as follows: to be used for
Performing adaptive filtering on each wave beam in the current active transmission period in the interference influence range for the reference signal to realize interference suppression; a least mean square error (LMS) algorithm is used here, a transversal filter structure is used;
let theta in the nth iteration process I The directional filter weight coefficient vector is w (θ) I N), then n+1th iteration
w(θ I ,n+1)＝w(θ I ,n)+μ(θ I )u(θ I ,n)(d * (n)-u H (θ I ,n)w(θ I ,n))
Wherein μ (θ) I ) For theta I Direction step, u (θ) I N) is theta I Azimuth beam input sampling, u H (θ I N) is u (theta) I Transpose of n), d * (n) is
Conjugation of the samples; the error signal can be expressed as
e(θ I ,n)＝d(n)-w T (θ I ,n)u(θ I ,n)
Wherein e (θ) I N) is theta I The error of the nth iteration of the direction comprises filtering the target echo signal after the compression type interference.
10. The active sonar interference rejection method for a pressed interference source of claim 1 wherein S9 is specifically as follows: hair-to-hairThe shot signal is used as a reference, and the correlation e (theta I N) performing matching processing
Wherein Z (θ) I τ) is θ I The matching envelope samples of the directional delay tau,
to sample the transmitted signal S L (n) performing delay conjugation; for each azimuth Z (theta) I τ), and can obtain azimuth distance information about the target; since the interference is suppressed at this time, the target echo detection performance is better than before the interference suppression. />
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211473519.2A CN116047488A (en) | 2022-11-22 | 2022-11-22 | Active sonar interference suppression method for pressed interference source |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211473519.2A CN116047488A (en) | 2022-11-22 | 2022-11-22 | Active sonar interference suppression method for pressed interference source |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116047488A true CN116047488A (en) | 2023-05-02 |
Family
ID=86128454
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211473519.2A Pending CN116047488A (en) | 2022-11-22 | 2022-11-22 | Active sonar interference suppression method for pressed interference source |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116047488A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116819430A (en) * | 2023-06-30 | 2023-09-29 | 中国人民解放军军事科学院系统工程研究院 | Direction finding method for same-frequency signal under strong radiation source background |
-
2022
- 2022-11-22 CN CN202211473519.2A patent/CN116047488A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116819430A (en) * | 2023-06-30 | 2023-09-29 | 中国人民解放军军事科学院系统工程研究院 | Direction finding method for same-frequency signal under strong radiation source background |
| CN116819430B (en) * | 2023-06-30 | 2024-03-15 | 中国人民解放军军事科学院系统工程研究院 | Direction finding method for same-frequency signal under strong radiation source background |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109521430B (en) | Distance extension target detection method for inhibiting narrow-band interference | |
| JP2011158471A (en) | Method for detecting target in time-space adaptive processing system | |
| CN111024209B (en) | Line spectrum detection method suitable for vector hydrophone | |
| CN104835503A (en) | Improved GSC self-adaptive speech enhancement method | |
| CN101876715A (en) | Method for suppressing noises from towing ships by using towed acoustic array | |
| CN110673116B (en) | Same frequency interference suppression method | |
| CN110085247A (en) | A kind of dual microphone noise-reduction method for complicated noise | |
| CN111693971B (en) | Wide-beam interference suppression method for weak target detection | |
| CN116047488A (en) | Active sonar interference suppression method for pressed interference source | |
| CN116520333A (en) | Multi-base buoy fusion detection method and system based on complex interference | |
| CN111090089B (en) | Space-time adaptive detection method based on two types of auxiliary data | |
| He et al. | Forward scattering detection of a submerged moving target based on adaptive filtering technique | |
| CN112505665B (en) | Space-time self-adaptive detection method and system suitable for partial uniform reverberation environment | |
| CN107728132B (en) | A method of improving interference blocking algorithm output signal-to-noise ratio | |
| CN116106879A (en) | Linear array line spectrum coherent accumulation detection method in multi-path environment | |
| CN115993591A (en) | Self-adaptive interference suppression method based on characteristic subspace judgment | |
| CN115656940A (en) | Space-frequency domain narrow-band composite interference resisting method and system | |
| CN115542304A (en) | Self-interference signal suppression method based on transfer function fitting and phase matching | |
| Saric et al. | Adaptive microphone array based on pause detection | |
| CN110632579B (en) | Iterative beam forming method using subarray beam domain characteristics | |
| CN109407078B (en) | Double-adaptive copy correlator method | |
| CN110007294A (en) | A kind of time domain interference cancellation method based on energy compensating | |
| CN110007296A (en) | A kind of time domain interference cancellation method based on guidance signal correction | |
| Shang et al. | Research on Feature Extraction and Classification Recognition of Underwater Small Target Active Sonar | |
| CN111257859B (en) | Wave beam domain self-correlation underwater target identification method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |