CN112487887B - Method for improving signal-to-noise ratio of active target detection - Google Patents

Method for improving signal-to-noise ratio of active target detection Download PDF

Info

Publication number
CN112487887B
CN112487887B CN202011278977.1A CN202011278977A CN112487887B CN 112487887 B CN112487887 B CN 112487887B CN 202011278977 A CN202011278977 A CN 202011278977A CN 112487887 B CN112487887 B CN 112487887B
Authority
CN
China
Prior art keywords
signal
sequence
active target
active
copy
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.)
Active
Application number
CN202011278977.1A
Other languages
Chinese (zh)
Other versions
CN112487887A (en
Inventor
韩宁
温子旭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN202011278977.1A priority Critical patent/CN112487887B/en
Publication of CN112487887A publication Critical patent/CN112487887A/en
Application granted granted Critical
Publication of CN112487887B publication Critical patent/CN112487887B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02Preprocessing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/523Details of pulse systems
    • G01S7/526Receivers
    • G01S7/527Extracting wanted echo signals
    • G01S7/5273Extracting wanted echo signals using digital techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/42Global feature extraction by analysis of the whole pattern, e.g. using frequency domain transformations or autocorrelation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12Classification; Matching
    • G06F2218/14Classification; Matching by matching peak patterns

Landscapes

  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Multimedia (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a method for improving the signal-to-noise ratio of active target detection, which is suitable for the detection of an active target. The method comprises the steps of designing a plurality of copy signals by setting the scale and the speed range of an active target and performing matched filtering operation with received signals one by one; searching the maximum value in the matched filtering results to obtain the scale and speed information which are most matched with the echo signal of the active target, and designing a copy signal according to the scale and speed information; and performing matched filtering operation again by using the optimized copy signal and the echo signal, so that the signal-to-noise ratio of the active target detection can be improved. By the optimization design of the copy signal, the signal-to-noise ratio of the matched filtering result is effectively improved, and therefore the detection capability of the active target is improved.

Description

Method for improving signal-to-noise ratio of active target detection
Technical Field
The invention relates to a method for improving the detection signal-to-noise ratio of an active target, in particular to a method for improving the detection signal-to-noise ratio of the active target, which is suitable for the detection of the active target.
Background
Signal detection is an important task of detection equipment such as sonar and radar. The active sonar or the radar transmits sound waves or electromagnetic waves in a certain form, and performs signal detection and feature extraction on the reflected echoes of the target, so that the active target is positioned, classified and tracked. Classical detection theory states that the best receiver to detect a known signal in a white gaussian noise background is a matched filter, which is also the most basic receiver commonly used in many detection systems such as sonar. The response function of the matched filter is a delayed conjugate mirror waveform of the matched signal u (t), which is also referred to as the replica signal or reference signal of the matched filter. For a typical input signal R (t), the output of the matched filter is y (t) = R vu (t 0 -t), wherein R vu Is the cross-correlation value of signal r (t) against u (t). t = t 0 When, y (t) 0 )=R vu (0) Thus, the matched filter acts as a cross-correlator which can calculate the cross-correlation function.
In practical applications, the copy signal u (t) is generally a transmission signal in an active sonar or radar system. However, due to the reasons of complex channel influence, serious background interference and the like, and the sound stealth capability of the active target has been greatly improved in recent years, the functional requirements of an active detection system cannot be met by applying a conventional matched filtering method to active target detection. How to improve the detection capability of the active target, especially the performance of the matched filtering process, is one of the hot problems of the current research.
Disclosure of Invention
The purpose of the invention is as follows: in order to overcome the defects of the prior art, the method for improving the signal-to-noise ratio of the active target detection is provided, and the signal-to-noise ratio of the matched filtering result is effectively improved through the optimization design of the copy signal, so that the detection capability of the active target is improved.
In order to achieve the technical purpose, the method for improving the signal-to-noise ratio of active target detection of the invention uses an active sonar to emit an acoustic signal into water, the signal meets an active target to generate an echo signal when being transmitted in water, the emitted signal is selected from a linear frequency modulation signal or a hyperbolic frequency modulation signal, the pulse width of the emitted signal is T seconds, and the lower limit frequency is f 1 Hertz at an upper limit frequency f 2 Hertz, the sampling frequency of the active sonar system is f s Hertz, the propagation speed of sound waves in water is c m/s; the radial speed of the active target is v meters per second, and the dimension is l meters; the signal-to-noise ratio of the received signal is A decibels;
the method comprises the following steps:
firstly, setting a scale range of an active target, and extracting M scale values in the scale range of the active target; then setting a speed range of the active target, and extracting N speed values from the range; time series X constructed sequence X formed by active sonar emission signal i (ii) a Will sequence X i Performing variable sampling processing to obtain a copy signal sequence X corresponding to each speed ij (ii) a Then, the active sonar transmits the echo signal sequence R and copy signal sequence X received after the emission signal ij Performing matched filtering operation and taking the maximum value to obtain a matrix with the size of M multiplied by N; the mth row and the nth column of the maximum position in the matrix are the positions where the matching peak values are maximum; using sequence X i Constructing a sequence Y, and carrying out variable sampling processing on the sequence Y to obtain a copy signal sequence U; and matching and filtering the echo signal sequence R and the copy signal sequence U received by the active target detection system, so that the signal-to-noise ratio of the active target detection can be improved.
The scale range of the active target is Len 1 ~Len 2 And M scale values extracted from the range are respectively L i =Len 1 +i*(Len 2 -Len 1 ) (M-1), i =0, 1.., M-1, the dimension of the active target to be measured does not exceed the lower boundary Len 1 And an upper bound Len 2 A range of (a); master and slaveThe speed range of the moving target is Vel 1 ~Vel 2 Extracting N speed values, each V j =Vel 1 +j*(Vel 2 -Vel 1 ) V (N-1), j =0, 1.., N-1, with a speed range of the lower bound Vel 1 To the upper boundary Vel 2 In the meantime.
Using the formula: x i =x 1i +x 2i Wherein
Figure BDA0002780078150000021
Time series X construction sequence X formed by active sonar emission signal i And c is the speed of wave propagation.
Sequence X i I =0, 1.., M-1 is subjected to a variable sampling process to obtain a signal corresponding to each velocity V j Copy signal sequence X of ij : sampling rate of f s Is changed into f s c/(c-2V j ),j=0,1,...,N-1;
An echo signal sequence R and a copy signal sequence X received by the active sonar ij Performing matched filtering operation and taking the maximum value Max ij I =0,1, ·, M-1, j =0,1, ·, N-1; this results in a matrix of size mxn:
Figure BDA0002780078150000022
sequence Y = Y 1 +Y 2 Wherein, in the step (A),
Figure BDA0002780078150000023
the sequence Y is subjected to variable sampling processing, namely the sampling rate is f s Is changed into f s c/(c-2V n ) And obtaining a copy signal sequence U.
All copy signals in the previous period can be stored and directly selected to be U = X mn As a copy signal, or the sequence Y is subjected to a variable sampling process to regenerate the copy signal and applied to avoid the need to apply X ij I =0, 1., M-1, j =0, 1., N-1 are all stored, thereby reducing the occupation of storage resources.
Has the advantages that: the invention searches the scale information and the speed information of the active target in a certain range, and optimizes the copy signal, thereby occupying less storage resources and effectively improving the signal-to-noise ratio of the active target detection.
Drawings
FIG. 1 is a flow chart of a method of improving active target detection signal-to-noise ratio in accordance with the present invention;
fig. 2 is a graph comparing the matched filtering result of the transmit signal as the copy signal with the matched filtering result of the optimized copy signal according to the present invention.
Detailed Description
Embodiments of the invention are further described below with reference to the accompanying drawings:
as shown in figure 1, the method for improving the signal-to-noise ratio of active target detection of the invention uses an active sonar to emit an acoustic signal into water, the signal meets the active target to generate an echo signal when propagating in water, the emitted signal is selected from a linear frequency modulation signal or a hyperbolic frequency modulation signal, the pulse width of the emitted signal is T seconds, and the lower limit frequency is f 1 Hertz at an upper frequency limit of f 2 Hertz, the sampling frequency of the active sonar system is f s Hertz, the propagation speed of sound waves in water is c meters per second; the radial speed of the active target is v meters per second, and the dimension is l meters; the signal-to-noise ratio of the received signal is A decibels;
the method comprises the following steps:
firstly, setting the scale range of an active target, wherein the scale range of the active target is Len 1 ~Len 2 And M scale values extracted from the range are respectively L i =Len 1 +i*(Len 2 -Len 1 ) (M-1), i =0, 1.., M-1, the dimension of the active target to be measured does not exceed the lower boundary Len 1 And an upper bound Len 2 A range of (a); the speed range of the active target is Vel 1 ~Vel 2 Extracting N speed values, each V j =Vel 1 +j*(Vel 2 -Vel 1 ) V (N-1), j =0, 1.., N-1, with a speed range of the lower bound Vel 1 To the upper boundary Vel 2 To (c) to (d);
then setting the active targetAnd extracting N speed values from the range; time series X construction sequence X formed by active sonar emission signals i Using the formula: x i =x 1i +x 2i In which
Figure BDA0002780078150000031
Figure BDA0002780078150000041
Time series X construction sequence X formed by active sonar emission signal i And c is the speed of wave propagation;
will sequence X i I =0, 1., M-1 is subjected to a variable sampling process to obtain a signal corresponding to each velocity V j Copy signal sequence X of ij : sampling rate of f s Is changed into f s c/(c-2V j ) J =0, 1.., N-1. Copy signal generation by simulation software, such as variable sampling by using sample in matlab, internal algorithm involving FIR filter design and rate conversion, and other methods
Figure BDA0002780078150000042
The sampling interval is from 0 second to
Figure BDA0002780078150000043
Sampling is carried out within seconds, namely T is set as the pulse width of the signal, then
Figure BDA0002780078150000044
At a frequency f 0 Taking the single-frequency signal of (1) as an example, the sequence x = sin (2 π f) can be obtained 0 t) wherein F s =f s When x is the original copy signal; f s =f s c/(c-2V j ) X is the copy signal after variable sampling.
Then, an echo signal sequence R and a copy signal sequence X received by the active sonar are obtained ij Performing matched filtering operation and taking the maximum value Max ij I =0,1, · M-1, j =0,1,. · N-1; this results in a matrix of size mxn:
Figure BDA0002780078150000045
the mth row and the nth column of the maximum position in the matrix are the positions where the matching peak values are maximum;
using sequence X i Constructing a sequence Y, and performing variable sampling processing on the sequence Y to obtain a copy signal sequence U; sequence Y = Y 1 +Y 2 Wherein, in the step (A),
Figure BDA0002780078150000046
the sequence Y is subjected to variable sampling processing, namely the sampling rate is changed from f s Is changed to f s c/(c-2V n ) Obtaining a copy signal sequence U;
the echo signal sequence R received by the active target detection system and the copy signal sequence U are subjected to matched filtering operation, so that the signal-to-noise ratio of active target detection can be improved; all the copy signals in the previous stage can be directly stored, and U = X is selected mn As a copy signal, or by regenerating the copy signal and applying it in the previous step, to avoid copying X ij I =0, 1., M-1, j =0, 1., N-1 are all stored, thereby reducing the occupation of storage resources.
The first embodiment is as follows:
assuming that the emission signal of the active sonar is a linear frequency modulation signal, the pulse width is 1 second, the lower limit frequency is 950Hz, the upper limit frequency is 1050Hz, and the sampling frequency f of the active sonar system s 20kHz, and the propagation speed c of sound waves in water is 1500 m/s; the radial speed of the active target is 2 m/s, and the dimension is 100 m; the signal-to-noise ratio of the received signal is-20 dB. The signal-to-noise ratio of active target detection is improved according to the following steps:
(1) Setting the scale range of the active target to be 20-200 m, and sampling at equal intervals in the range to obtain 91 scale values L i = (20 + i × 2) meters, i =0,1, ·,90;
(2) Construction sequence X i =x 1i +x 2i Wherein
Figure BDA0002780078150000051
i =0,1.., 90,x is the sequence of transmitted signals in an active target detection system;
(3) Setting the speed range of the active target to be-15 m/s, and sampling at equal intervals in the range to obtain 61 speed values which are V respectively j = (-15 + j 0.5) m/s, j =0,1, ·,60;
(4) Will sequence X i I =0, 1.., 90 is subjected to a variable sampling process, i.e., the sampling rate is represented by f s Is changed into f s c/(c-2V j ) J =0, 1.., 60, obtained for each speed V j Copy signal X of ij
(5) The echo signal sequence R and the copy signal sequence X are combined ij Performing matched filtering operation and taking the maximum value Max ij I =0, 1., 90, j =0, 1., 60, resulting in a matrix size of 91 × 61
Figure BDA0002780078150000052
(6) Searching a maximum value MAX in a matrix Max to obtain an m-th row and an n-th column of the matrix at the position of the matrix, wherein m and n can change due to the fact that noise in each experiment is not identical;
(7) Construction sequence Y = Y 1 +Y 2 Wherein
Figure BDA0002780078150000053
And the sequence Y is subjected to variable sampling processing, namely the sampling rate is changed from f s Is changed into f s c/(c-2V n ) Obtaining a copy signal sequence U;
(8) And (4) carrying out matched filtering operation on the echo signal sequence R and the copy signal sequence U in (7). Fig. 2 is a graph comparing the matched filtering results of a transmit signal as a copy signal and the optimized copy signal according to the present invention. As indicated by the circle marks in the figure, after the copy signals are optimally designed, whether echo signals exist or not can be more easily distinguished from the matched filtering results, and the arrival time of the echo signals is obtained; compared with the traditional matched filtering result which takes the emission signal as the copy signal, the detection effect of the active target can be effectively improved.

Claims (7)

1. A method for improving the signal-to-noise ratio of active target detection uses an active sonar to emit an acoustic signal into water, the signal meets an active target to generate an echo signal when being transmitted in the water, and the emitted signal is selected to be a linear frequency modulation signal or a hyperbolic frequency modulation signal, which is characterized in that: the pulse width of the transmitted signal is T seconds, and the lower limit frequency is f 1 Hertz at an upper frequency limit of f 2 Hertz, the sampling frequency of the active sonar system is f s Hertz, the propagation speed of sound waves in water is c m/s; the radial speed of the active target is v m/s, and the scale is l m; the signal-to-noise ratio of the received signal is A decibels;
the method comprises the following steps:
firstly, setting a scale range of an active target, and extracting M scale values in the scale range of the active target; then setting a speed range of the active target, and extracting N speed values from the range; time series X constructed sequence X formed by active sonar emission signal i (ii) a Will sequence X i Performing variable sampling processing to obtain a copy signal sequence X corresponding to each speed ij (ii) a Then, the active sonar transmits the echo signal sequence R and copy signal sequence X received after the emission signal ij Performing matched filtering operation and taking the maximum value to obtain a matrix with the size of M multiplied by N; the mth row and the nth column of the maximum value position in the matrix are the positions where the matching peak values are maximum; using the sequence X i Constructing a sequence Y, and carrying out variable sampling processing on the sequence Y to obtain a copy signal sequence U; and matching and filtering the echo signal sequence R and the copy signal sequence U received by the active target detection system, so that the signal-to-noise ratio of the active target detection can be improved.
2. The method of claim 1 for improving signal-to-noise ratio of active target detection, wherein: the scale range of the active target is Len 1 ~Len 2 And M scale values extracted from the range are respectively L i =Len 1 +i*(Len 2 -Len 1 ) (M-1), i =0, 1.., M-1, the dimension of the active target to be measured does not exceed the lower boundary Len 1 And an upper bound Len 2 A range of (d); the speed range of the active target is Vel 1 ~Vel 2 Extracting N speed values, each V j =Vel 1 +j*(Vel 2 -Vel 1 ) V (N-1), j =0, 1.., N-1, with a speed range of the lower bound Vel 1 To the upper boundary Vel 2 In between.
3. The method of claim 1, wherein the formula is used to: x i =x 1i +x 2i Wherein
Figure FDA0002780078140000011
Time series X construction sequence X formed by active sonar emission signal i And c is the speed of wave propagation.
4. The method of claim 1 for improving signal-to-noise ratio of active target detection, wherein: sequence X i I =0, 1.., M-1 is subjected to a variable sampling process to obtain a signal corresponding to each velocity V j Copy signal sequence X of ij : sampling rate of f s Is changed into f s c/(c-2V j ),j=0,1,...,N-1。
5. The method of claim 1 for improving the signal-to-noise ratio of active target detection, wherein: an echo signal sequence R and a copy signal sequence X received by the active sonar ij Performing matched filtering operation and taking the maximum value Max ij I =0,1, · M-1, j =0,1,. · N-1; this results in a matrix of size mxn:
Figure FDA0002780078140000021
6. according to the rightThe method for improving the signal-to-noise ratio of active target detection according to claim 1, characterized in that: sequence Y = Y 1 +Y 2 Wherein, in the step (A),
Figure FDA0002780078140000022
the sequence Y is subjected to variable sampling processing, namely the sampling rate is changed from f s Is changed into f s c/(c-2V n ) And obtaining a copy signal sequence U.
7. The method of claim 5, wherein the step of improving the signal-to-noise ratio of active target detection comprises: all copy signals in the previous period can be stored and directly selected to be U = X mn As a copy signal, or the sequence Y is subjected to a variable sampling process to regenerate the copy signal and applied to avoid the need to apply X ij I =0, 1., M-1, j =0, 1., N-1 are all stored, thereby reducing the occupation of storage resources.
CN202011278977.1A 2020-11-16 2020-11-16 Method for improving signal-to-noise ratio of active target detection Active CN112487887B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011278977.1A CN112487887B (en) 2020-11-16 2020-11-16 Method for improving signal-to-noise ratio of active target detection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011278977.1A CN112487887B (en) 2020-11-16 2020-11-16 Method for improving signal-to-noise ratio of active target detection

Publications (2)

Publication Number Publication Date
CN112487887A CN112487887A (en) 2021-03-12
CN112487887B true CN112487887B (en) 2022-11-08

Family

ID=74930523

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011278977.1A Active CN112487887B (en) 2020-11-16 2020-11-16 Method for improving signal-to-noise ratio of active target detection

Country Status (1)

Country Link
CN (1) CN112487887B (en)

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106019263B (en) * 2016-07-13 2018-03-20 东南大学 Target radial speed measuring method based on more bright spot echo models
CN109061625B (en) * 2018-06-22 2020-09-18 西安空间无线电技术研究所 Distance and speed matched moving target detection method
CN109507642A (en) * 2018-12-29 2019-03-22 广西科技大学 A kind of Nonlinear Parameter harmonic detecting method based on noise FM technology

Also Published As

Publication number Publication date
CN112487887A (en) 2021-03-12

Similar Documents

Publication Publication Date Title
CN107450055B (en) High-speed maneuvering target detection method based on discrete linear frequency modulation Fourier transform
CN109814073A (en) A kind of MTD radar ambiguity solution speed-measuring method
JP4899031B2 (en) Fast coherent processing of codes for periodic line spectra.
CN103018727A (en) Sample-training-based non-stationary clutter suppression method of vehicle-mounted radar
CN107356908B (en) Frequency agile signal coherent accumulation method
CN115219997A (en) Multi-intermittent sampling interference resisting method based on cognitive waveform and filter combined design
CN110308427A (en) LFM pulse radar frequency-domain impulse compression processing method based on FPGA
CN109407098B (en) Method for realizing parallel processing algorithm for processing orthogonal pulse train signals
CN112487887B (en) Method for improving signal-to-noise ratio of active target detection
CN100485413C (en) Frequency-scanning signal time-frequency correlation detection and time delay estimating method
CN111413702B (en) Efficient target segmentation method for broadband fish finder
CN110632605B (en) Wide-tolerance large-aperture towed linear array time domain single-beam processing method
CN104580055B (en) A kind of Doppler estimation based on construction target derived function determination section search
CN112014807A (en) Self-adaptive clutter suppression method for frequency agile radar
CN110673118A (en) Active sonar single-frequency pulse train waveform design and detection algorithm
CN114428244B (en) Super-resolution time delay estimation method and device based on sparse learning
CN104793197A (en) Direct-wave suppression method based on IFFT frequency spectrum division method and gradient adaptive lattice filter
CN113704685B (en) Deep sea blind deconvolution method based on vertical line array
CN113359095B (en) Coherent passive MIMO radar Clarithrome boundary calculation method
CN108459316B (en) Echo signal parameter estimation method under multi-path channel
CN114563777A (en) Time delay Doppler frequency joint estimation method based on two-dimensional off-grid sparse learning
Liu et al. Detection performance analysis of sub-band processing continuous active sonar
CN114325671B (en) Improved Doppler method based on underwater acoustic beacon signal multi-pulse accumulation
CN110632579A (en) Iterative beam forming method using subarray beam domain characteristics
Park et al. High-resolution range and velocity estimation method based on generalized sinusoidal frequency modulation for high-speed underwater vehicle detection

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
GR01 Patent grant
GR01 Patent grant