CN108680786B - Adaptive filtering envelope extraction method for pulse signal frequency domain - Google Patents

Adaptive filtering envelope extraction method for pulse signal frequency domain Download PDF

Info

Publication number
CN108680786B
CN108680786B CN201810325134.9A CN201810325134A CN108680786B CN 108680786 B CN108680786 B CN 108680786B CN 201810325134 A CN201810325134 A CN 201810325134A CN 108680786 B CN108680786 B CN 108680786B
Authority
CN
China
Prior art keywords
power spectrum
smoothing
pulse signal
signal
envelope
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
CN201810325134.9A
Other languages
Chinese (zh)
Other versions
CN108680786A (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 CN201810325134.9A priority Critical patent/CN108680786B/en
Publication of CN108680786A publication Critical patent/CN108680786A/en
Application granted granted Critical
Publication of CN108680786B publication Critical patent/CN108680786B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a method for extracting a pulse signal frequency domain self-adaptive filtering envelope, which comprises the following steps: the first step is as follows: acquiring a sampling data sequence of a pulse signal to be processed; the second step is that: calculating a power spectrum of the sampling data sequence; the third step: iterative smoothing and filter bandwidth parameter initialization; the fourth step: carrying out iterative smoothing processing on the power spectrum to obtain a smooth power spectrum; the fifth step: extracting power spectrum main lobe characteristic parameters; and a sixth step: generating adaptive filtering parameters; the seventh step: self-adaptive filtering in a frequency domain; eighth step: and extracting the envelope of the pulse signal. The extraction method generates the self-adaptive filtering parameters by extracting the main lobe characteristics of the pulse signals in real time, carries out self-adaptive filtering on the pulse signals, can realize an envelope extraction method matched with the type and the frequency of the pulse signals under the non-cooperative condition, has small operand and strong practicability, and is suitable for carrying out real-time processing on the signals.

Description

Adaptive filtering envelope extraction method for pulse signal frequency domain
Technical Field
The invention belongs to the technical field of signal processing, and particularly relates to a pulse signal frequency domain adaptive filtering envelope extraction method.
Background
The signal envelope is an important parameter for representing the time domain characteristics of the pulse signals, the accurate extraction of the pulse signal envelope can provide a basis for further accurate estimation of the arrival time and the pulse width of the pulse signals, the time delay difference of the two paths of pulse signal envelopes can be used for estimating the incoming wave direction of the signals, and the method has important theoretical and application values in the fields of sonar, radar, electronic warfare and the like, and is particularly prominent in the aspect of radar and sonar signal processing.
The conventional signal envelope extraction methods mainly comprise four types, namely (1) a Hilbert transform method, (2) a wavelet transform method, (3) an improved Hilbert transform method, such as an envelope extraction method based on Hilbert transform (or called Hilbert-yellow transform) of empirical mode decomposition and a Hilbert transform method based on blind source extraction, and (4) an envelope extraction method based on mathematical morphology processing, wherein the Hilbert transform method firstly converts an original signal into a complex analysis signal, and then takes a modulus value of the complex analysis signal as an envelope of the signal, and the method is based on FFT, simple in implementation, small in operand and good in envelope extraction effect under high signal-to-noise ratio, but is performed in the whole frequency interval and poor in noise suppression capability, and is very useful for envelope extraction of a narrow-band carrier signal, but has two inevitable defects, namely ① extraction of high-frequency components in the signal, so that time-varying adaptive analysis is not achieved ②.
The wavelet transformation method firstly carries out wavelet transformation on signals to obtain analytic signals, wherein the wavelet transformation is equivalent to multi-passband filtering on the signals by utilizing a series of bandpass filters, so that the obtained analytic signals are subjected to the bandpass filtering, and the real parts and the imaginary parts of the obtained analytic signals are orthogonal, and then the modulus values of the analytic signals are taken to obtain the envelopes of the signals. Compared with the Hilbert transform method, the envelope extracted by the method is more ideal, but the method has the difficulty that the selection of the wavelet transform scale is proper, the scale selection is proper, the frequency band covered by the band-pass filter can be matched with the signal frequency, and the envelope extraction effect is good; the scale selection is improper, so that the frequency band covered by the band-pass filter is mismatched with the signal frequency, the envelope extraction effect is rapidly deteriorated, the selection of the scale depends on the experience of a user to a great extent, and the universality is poor. The improved Hilbert transform method firstly utilizes a certain specific preprocessing method to carry out filtering and denoising processing on signals, and then further utilizes Hilbert transform to extract envelopes, such as the empirical mode decomposition and blind source extraction, and the like. The processing method based on mathematical morphology comprises the steps of preprocessing signals by utilizing morphological filtering, extracting signal envelopes by utilizing morphological closed operation, and eliminating noise envelopes by utilizing morphological open operation.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the problems and the defects in the prior art, the invention provides a pulse signal frequency domain adaptive filtering envelope extraction method to meet the requirements of radar and sonar signal processing. In addition, the method has the same operation amount as the traditional Hilbert transform envelope extraction method, but has better noise suppression capability and self-adaption capability compared with the Hilbert transform method, clear principle, simple realization and strong engineering applicability.
The technical scheme is as follows: in order to achieve the purpose, the invention adopts the following technical scheme: a method for extracting a pulse signal frequency domain adaptive filtering envelope comprises the following steps:
(1) acquiring a pulse signal sampling data sequence to be processed: receiving real-time acquisition data of N sampling points from a sensor as a data sequence x (N) to be processed, wherein N is 0,1, …, N-1, or extracting N sampling point data containing the whole pulse signal from a memory as the data sequence x (N) to be processed, wherein N is 0,1, …, N-1, the number of the sampling points containing the whole pulse signal is N, and the value is an integer power of 2.
(2) Calculating a pulse signal power spectrum P (k) according to the data sequence x (n).
(3) Initialization of parameters related to the bandwidth of the iterative smoothing and filter: setting a maximum iteration number threshold of power spectrum iterative smoothing, wherein the precision control threshold and the iterative smoothing window length are respectively I11And M1The maximum iteration number threshold, the precision control threshold and the iteration smoothing window length of the initial envelope iteration smoothing are respectively I22And M2Bandwidth coefficient omega of narrow-band signal1And bandwidth factor omega of broadband signal2And judging characteristic parameter thresholds η of main lobes of the narrow-band pulse and the wide-band pulse.
(4) And performing iterative smoothing processing on the P (k) to obtain a smoothed power spectrum S (k).
(5) Extracting the original power spectrum and the main lobe characteristic parameter gamma of the smoothed power spectrumpAnd gammaS
(6) Adaptive filtering parameters h (k) are generated.
(7) And carrying out adaptive filtering on the pulse signal according to the filtering parameters H (k).
(8) The pulse signal envelope y (n) is extracted.
Preferably, in the step (2), fast fourier transform is performed on the data sequence x (n), and discrete fourier transform x (l) and power spectrum p (k) of the data sequence are obtained through calculation, specifically including the following steps:
step 2-1: calculating the discrete Fourier transform of x (n):
Figure GDA0002462721590000031
where l is the discrete frequency index of X (l), and j represents an imaginary unit, i.e.
Figure GDA0002462721590000032
This equation can be implemented by fast fourier transform.
Step 2-2: calculating the power spectrum of x (n) according to x (l):
Figure GDA0002462721590000033
where k is the discrete frequency index of the power spectrum p (k), | represents the modulo operation, and | x (k) | is the signal amplitude spectrum.
Preferably, in the step (2), the discrete fourier transform x (l) of the step 2-1 calculating x (n) can be realized by fast fourier transform;
preferably, in step (3), a maximum iteration number threshold of power spectrum iterative smoothing is set, and the precision control threshold and the iterative smoothing window length are respectively I11And M1The maximum iteration number threshold, the precision control threshold and the iteration smoothing window length of the initial envelope iteration smoothing are respectively I22And M2Bandwidth coefficient omega of narrow-band signal1And bandwidth factor omega of broadband signal2Judging the characteristic parameter threshold η of the main lobe of the narrow-band pulse and the wide-band pulse, wherein I1And I2A value is a positive integer greater than 2,1and2the value is positive number less than 1, the length of the iterative smoothing window of the power spectrum satisfies 3 and is less than or equal to M1An odd number less than or equal to N/2-2, and an envelope iteration smoothing window length of 5 less than or equal to M2Odd number, omega, not more than 491To satisfy 0 ≤ ω1Any number, omega, not more than 0.052To satisfy 0.05 < omega2Any number not more than 0.1, η is any number satisfying 1.5 not more than η not more than 2.51The value is 10, and the number is,1value of 0.01, M1A value of 5, I2The value is 5, and the number of the grooves,2a value of 0.02, M2Value of 11, omega1Has a value of 0.05, omega2Is 0.1 and η is 2.0.
Preferably, in the step (4), the iterative smoothing processing is performed on the power spectrum p (k) to obtain a smoothed power spectrum s (k), and the method specifically includes the following steps:
step 4-1: the iterative smoothing times i of the initialized power spectrum is 0, and the iterative smoothing result S of the power spectrum0(k) Comprises the following steps:
S0(k)=P(k),k=0,1,2…,N/2-1
step 4-2: let the iteration number i be i +1, and smooth the last iteration of the power spectrum to obtain a result Si-1(k) Performing smoothing to obtain current smoothing result Si(k):
Figure GDA0002462721590000041
Step 4-3: judging whether the maximum iterative smoothing times of the power spectrum is reached, namely judging that I is less than or equal to I1If yes, entering the step 4-4, otherwise jumping to the step 4-6;
step 4-4: respectively calculating last smoothing results S of power spectrumi-1(k) And current smoothing result Si(k) Sum of squares J of residuals from original power spectrum p (k)i-1And JiAre respectively as
Figure GDA0002462721590000042
And 4-5: judgment of | Ji-1-Ji|≤1JiIf yes, entering the step 4-6, otherwise returning to the step 4-2; wherein1A threshold is controlled for power spectrum iterative smoothing precision;
and 4-6: let S (k) become Si(k) And k is 0,1,2 …, N/2-1, and a smoothed power spectrum s (k) is obtained.
Preferably, in the step (5), the original power spectrum and the smoothed power spectrum main lobe characteristic parameter γ are extractedpAnd gammaSThe method specifically comprises the following steps:
step 5-1: searching discrete frequency index k corresponding to maximum value of smoothed power spectrum S (k)M
Figure GDA0002462721590000043
Wherein
Figure GDA0002462721590000044
Representing that the discrete frequency index corresponding to the maximum value of S (k) is searched within the range of 1 ≦ k ≦ N/2-1;
step 5-2: subtracting half of the maximum value from s (k) to obtain a half-power difference power spectrum z (k):
Z(k)=S(k)-S(kM)/2,k=0,1,2…,N/2-1
step 5-3: searching for the discrete frequency index k of the first zero crossing on the left of the Z (k) peakLThe searching process comprises the following steps:
step 5-3-1: initializing a frequency start search index k of a first zero crossing on the left side of a Z (k) peak valuel=kM-1;
Step 5-3-2: judgment of k l0 or Z (k)l) If the result is less than 0, jumping to the step 5-3-4 if the result is true, and otherwise, entering the step 5-3-3;
step 5-3-3: let klk l1 and returning to step 5-3-2;
step 5-3-4: let kL=klObtaining the discrete frequency index k of the first zero-crossing point on the left side of the Z (k) peak valueL
Step 5-4: searching for the first zero-crossing discrete frequency index k to the right of the Z (k) peakRThe searching process comprises the following steps:
step 5-4-1: initializing a frequency start search index k of a first zero crossing on the right side of a Z (k) peak valuer=kM+1;
Step 5-4-2: judgment of krN/2-1 or Z (k)r) If the value is less than 0, jumping to the step 5-4-4 if the value is less than 0, otherwise, entering the step 5-4-3;
step 5-4-3: let kr=kr+1 and returning to step 5-4-2;
step 5-4-4: let kR=krObtaining the discrete frequency index k of the first zero-crossing point on the right side of the Z (k) peak valueR
Step (ii) of5-5: calculating the primary power spectrum P (k) main lobe characteristic parameter gammap
Figure GDA0002462721590000051
And 5-6: calculating a main lobe characteristic parameter gamma of the smoothed power spectrum S (k)S
Figure GDA0002462721590000052
Preferably, in the step (6), generating the adaptive filtering parameter h (k) specifically includes the following steps:
step 6-1: initializing adaptive filter parameters
H(k)=0,k=0,1,2…,N-1
Step 6-2: respectively calculating the discrete frequency indexes T of the start and the end of the transition band of the narrow-band pulse signal filterN1And TN2And a passband start and end discrete frequency index PN1And PN2
TN1=max(round[(1-0.5ω1)kM],0),TN2=min(round[(1+0.5ω1)kM],N/2-1)
PN1=max(kM-1,TN1),PN2=min(kM+1,TN2)
Where round [ ] represents a rounding operation, and max (,) and min (,) represent a large and small value operation, respectively.
Step 6-3: respectively calculating the start and end discrete frequency indexes T of the transition band of the broadband pulse signal filterW1And TW2And a passband start and end discrete frequency index PW1And PW2
TW1=max(round[(1-0.5ω2)kM],0),TW2=min(round[(1+0.5ω2)kM],N/2-1)
PW1=max(kL-1,TW1),PW2=max(kR+1,TW2)
Step 6-4: obtaining discrete frequency index T from adaptive filter transition band start and end according to signal type1And T2And a passband start and end discrete frequency index P1And P2
Namely, max (gamma)ps)/min(γps) If > η is true, if true, the signal is judged to be a narrow-band pulse signal and the following steps are carried out:
T1=TN1,T2=TN2,P1=PN1,P2=PN2
otherwise, judging the signal to be a broadband pulse signal and ordering:
T1=TW1,T2=TW2,P1=PW1,P2=PW2
step 6-5: the transition band part of the adaptive filter coefficients H (k) is modified so that k satisfies T1≤k≤T2The part of (A) is:
Figure GDA0002462721590000061
step 6-6: the pass band part of the adaptive filter coefficients H (k) is modified such that k satisfies P1≤k≤P2The part of (A) is: h (k) ═ 1.
Preferably, in step (7), the result of discrete fourier transform x (k) of the pulse signal is adaptively filtered to obtain a filtered discrete fourier transform y (k):
Y(k)=X(k)H(k),k=0,1,2…,N-1
preferably, in the step (8), the extracting of the pulse signal envelope y (N), where N is 0,1 …, N-1 specifically includes the following steps:
step 8-1: performing inverse discrete Fourier transform on the filtered discrete Fourier transform Y (k) to obtain a pulse signal filtered time-domain complex signal yc(n):
Figure GDA0002462721590000062
Step 8-2: taking a complex number yc(n) obtaining the initial envelope data sequence y of the pulse signal by the modulus value0(n):
y0(n)=|yc(n)|,n=0,1…,N-1
Step 8-3: initializing the iterative smoothing times q of the envelope data to be 0;
step 8-4: let the iteration number q be q +1 and smooth the result y for the last iteration of the envelope data sequenceq-1(n) smoothing to obtain current smoothing result yq(n):
Figure GDA0002462721590000071
And 8-5: judging whether the maximum iterative smoothing times of the envelope is reached, namely judging that q is less than or equal to I2If yes, entering the step 8-6, otherwise jumping to the step 8-8;
and 8-6: respectively calculating last smoothing result y of envelope data sequenceq-1(n) and current smoothing result yq(n) and the initial envelope data sequence y0Residual sum of squares Ψ of (n)q-1And ΨqRespectively is as follows:
Figure GDA0002462721590000072
and 8-7: determine | Ψq-1q|≤2ΨqIf yes, entering the step 8-8, otherwise returning to the step 8-4;
and 8-8: let y (n) be yq(N), N is 0,1,2 …, N-1, resulting in a sequence of signal envelope data:
y(n),n=0,1,2…,N-1。
preferably, in the step (8), the step 8-1 performs inverse discrete fourier transform on the filtered discrete fourier transform y (k) to obtain the pulse signal filtered time-domain complex signal yc(n) may be implemented by a fast fourier transform.
Has the advantages that: compared with the prior art, the method has the following beneficial effects:
(1) the envelope extraction method of the invention automatically judges the type of the pulse signal under the non-cooperative condition by extracting the main lobe characteristic of the pulse signal in real time, estimates the bandwidth parameter of the signal in real time, generates the adaptive filter parameter matched with the type and the frequency of the pulse signal, and carries out adaptive filtering on the pulse signal, thereby obtaining the signal processing gain of approximate matched filtering.
(2) According to the envelope extraction method, after the discrete Fourier transform result of the filtered signal is subjected to the inverse discrete Fourier transform and a modulus value is taken to obtain a signal initial envelope data sequence, the initial envelope is further subjected to iterative smoothing processing, burrs existing in the initial envelope are reduced, and the problem that the burrs in the envelope are more due to the existence of high-frequency components in the conventional Hilbert transform method is effectively solved.
Drawings
FIG. 1 is a schematic flow diagram of the process of the present invention;
FIG. 2 is the original power spectrum of the simulated pulse signal of example 1;
FIG. 3 is the power spectrum of the simulation pulse signal after smoothing in example 1;
FIG. 4 is the adaptive filter coefficients of example 1;
FIG. 5 is a simulated signal envelope extracted in example 1;
FIG. 6 is the original power spectrum of the simulated pulse signal of example 2;
FIG. 7 is the power spectrum of the simulation pulse signal after smoothing in example 2;
FIG. 8 is the adaptive filter coefficients of example 2;
fig. 9 is a simulated signal envelope extracted in example 2.
Detailed Description
The invention is further described with reference to the following figures and examples:
as shown in fig. 1, a method for extracting a frequency domain adaptive filtering envelope of an impulse signal includes the following steps:
(1) acquiring a pulse signal sampling data sequence to be processed: receiving real-time acquisition data of N sampling points from a sensor as a data sequence x (N) to be processed, wherein N is 0,1, …, N-1, or extracting N sampling point data containing the whole pulse signal from a memory as the data sequence x (N) to be processed, wherein N is 0,1, …, N-1, the number of the sampling points containing the whole pulse signal is N, and the value is an integer power of 2.
(2) Calculating a pulse signal power spectrum P (k) from the data sequence x (n): performing fast Fourier transform on the data sequence x (n), and calculating to obtain discrete Fourier transform X (l) and a power spectrum P (k) of the data sequence, wherein the method specifically comprises the following steps:
step 2-1: calculating the discrete Fourier transform of x (n):
Figure GDA0002462721590000081
where l is the discrete frequency index of X (l), and j represents an imaginary unit, i.e.
Figure GDA0002462721590000082
This equation can be implemented by fast fourier transform.
Step 2-2: calculating the power spectrum of x (n) according to x (l):
Figure GDA0002462721590000091
where k is the discrete frequency index of the power spectrum p (k), | represents the modulo operation, and | x (k) | is the signal amplitude spectrum.
In the step (2), the discrete fourier transform x (l) of the step 2-1 for calculating x (n) can be realized by fast fourier transform, so that the calculation efficiency of the algorithm is improved.
(3) Initialization of parameters related to the bandwidth of the iterative smoothing and filter: setting a maximum iteration number threshold of power spectrum iterative smoothing, wherein the precision control threshold and the iterative smoothing window length are respectively I11And M1Threshold of maximum iteration number of initial envelope iteration smoothing, precision controlThe system threshold and the iterative smoothing window are respectively I22And M2Bandwidth coefficient omega of narrow-band signal1And bandwidth factor omega of broadband signal2Judging characteristic parameter threshold η of main lobe of narrow-band pulse and wide-band pulse, where I1And I2A value is a positive integer greater than 2,1and2the value is positive number less than 1, the length of the iterative smoothing window of the power spectrum satisfies 3 and is less than or equal to M1An odd number less than or equal to N/2-2, and an envelope iteration smoothing window length of 5 less than or equal to M2An odd number of not more than 49, η is an arbitrary number satisfying 1.5 not more than η not more than 2.5, omega1To satisfy 0 ≤ ω1Any number, omega, not more than 0.052To satisfy 0.05 < omega2Any number less than or equal to 0.1. Preferred is I1The value is 10, and the number is,1value of 0.01, M1A value of 5, I2The value is 5, and the number of the grooves,2a value of 0.02, M2Value of 11, omega1Has a value of 0.05, omega2Is 0.1 and η is 2.0.
(4) Performing iterative smoothing processing on P (k) to obtain a smoothed power spectrum S (k): performing iterative smoothing processing on the power spectrum P (k) to obtain a smoothed power spectrum S (k), and specifically comprising the following steps:
step 4-1: the iterative smoothing times i of the initialized power spectrum is 0, and the iterative smoothing result S of the power spectrum0(k) Comprises the following steps:
S0(k)=P(k),k=0,1,2…,N2-1
step 4-2: let the iteration number i be i +1, and smooth the last iteration of the power spectrum to obtain a result Si-1(k) Performing smoothing to obtain current smoothing result Si(k):
Figure GDA0002462721590000101
Step 4-3: judging whether the maximum iterative smoothing times of the power spectrum is reached, namely judging that I is less than or equal to I1If yes, entering the step 4-4, otherwise jumping to the step 4-6;
step 4-4: respectively calculating last smoothing results S of power spectrumi-1(k) And current smoothing result Si(k) Residue of the original power spectrum P (k)Sum of squared differences Ji-1And JiRespectively is as follows:
Figure GDA0002462721590000102
and 4-5: judgment of | Ji-1-Ji|≤1JiIf yes, entering the step 4-6, otherwise returning to the step 4-2; wherein1A threshold is controlled for power spectrum iterative smoothing precision;
and 4-6: let S (k) become Si(k) And k is 0,1,2 … and N2-1, and a smoothed power spectrum s (k) is obtained.
(5) Extracting the original power spectrum and the main lobe characteristic parameter gamma of the smoothed power spectrumpAnd gammaSThe method specifically comprises the following steps:
step 5-1: searching discrete frequency index k corresponding to maximum value of smoothed power spectrum S (k)M
Figure GDA0002462721590000103
Wherein
Figure GDA0002462721590000104
Representing that the discrete frequency index corresponding to the maximum value of S (k) is searched within the range of 1 ≦ k ≦ N/2-1;
step 5-2: subtracting half of the maximum value from s (k) to obtain a half-power difference power spectrum z (k):
Z(k)=S(k)-S(kM)/2,k=0,1,2…,N/2-1
step 5-3: searching for the discrete frequency index k of the first zero crossing on the left of the Z (k) peakLThe searching process comprises the following steps:
step 5-3-1: initializing a frequency start search index k of a first zero crossing on the left side of a Z (k) peak valuel=kM-1;
Step 5-3-2: judgment of k l0 or Z (k)l) If the result is less than 0, jumping to the step 5-3-4 if the result is true, and otherwise, entering the step 5-3-3;
step 5-3-3: let klk l1 and returning to step 5-3-2;
step 5-3-4: let kL=klObtaining the discrete frequency index k of the first zero-crossing point on the left side of the Z (k) peak valueL
Step 5-4: searching for the first zero-crossing discrete frequency index k to the right of the Z (k) peakRThe searching process comprises the following steps:
step 5-4-1: initializing a frequency start search index k of a first zero crossing on the right side of a Z (k) peak valuer=kM+1;
Step 5-4-2: judgment of krN/2-1 or Z (k)r) If the value is less than 0, jumping to the step 5-4-4 if the value is less than 0, otherwise, entering the step 5-4-3;
step 5-4-3: let kr=kr+1 and returning to step 5-4-2;
step 5-4-4: let kR=krObtaining the discrete frequency index k of the first zero-crossing point on the right side of the Z (k) peak valueR
Step 5-5: calculating the primary power spectrum P (k) main lobe characteristic parameter gammap
Figure GDA0002462721590000111
And 5-6: calculating a main lobe characteristic parameter gamma of the smoothed power spectrum S (k)S
Figure GDA0002462721590000112
(6) Generating adaptive filtering parameters h (k), specifically including the following steps:
step 6-1: initializing adaptive filter parameters
H(k)=0,k=0,1,2…,N-1
Step 6-2: respectively calculating the discrete frequency indexes T of the start and the end of the transition band of the narrow-band pulse signal filterN1And TN2And a passband start and end discrete frequency index PN1And PN2
TN1=max(round[(1-0.5ω1)kM],0),TN2=min(round[(1+0.5ω1)kM],N/2-1)
PN1=max(kM-1,TN1),PN2=min(kM+1,TN2)
Where round [ ] represents a rounding operation, and max (,) and min (,) represent a large and small value operation, respectively.
Step 6-3: respectively calculating the start and end discrete frequency indexes T of the transition band of the broadband pulse signal filterW1And TW2And a passband start and end discrete frequency index PW1And PW2
TW1=max(round[(1-0.5ω2)kM],0),TW2=min(round[(1+0.5ω2)kM],N/2-1)
PW1=max(kL-1,TW1),PW2=max(kR+1,TW2)
Step 6-4: obtaining discrete frequency index T from adaptive filter transition band start and end according to signal type1And T2And a passband start and end discrete frequency index P1And P2
Namely, max (gamma)ps)/min(γps) If > η is true, if true, the signal is judged to be a narrow-band pulse signal and the following steps are carried out:
T1=TN1,T2=TN2,P1=PN1,P2=PN2
otherwise, judging the signal to be a broadband pulse signal and ordering:
T1=TW1,T2=TW2,P1=PW1,P2=PW2
step 6-5: the transition band part of the adaptive filter coefficients H (k) is modified so that k satisfies T1≤k≤T2The part of (A) is:
Figure GDA0002462721590000121
step 6-6: the pass band part of the adaptive filter coefficients H (k) is modified such that k satisfies P1≤k≤P2The part of (A) is: h (k) ═ 1;
(7) adaptively filtering the pulse signal according to the filter parameter h (k), that is, adaptively filtering the result x (k) of the discrete fourier transform of the pulse signal to obtain a filtered discrete fourier transform y (k):
Y(k)=X(k)H(k),k=0,1,2…,N-1
(8) the method for extracting the pulse signal envelope y (n) specifically comprises the following steps:
step 8-1: performing inverse discrete Fourier transform on the filtered discrete Fourier transform Y (k) to obtain a pulse signal filtered time-domain complex signal yc(n):
Figure GDA0002462721590000122
Step 8-2: taking a complex number yc(n) obtaining the initial envelope data sequence y of the pulse signal by the modulus value0(n):
y0(n)=|yc(n)|,n=0,1…,N-1
Step 8-3: initializing the iterative smoothing times q of the envelope data to be 0;
step 8-4: let the iteration number q be q +1 and smooth the result y for the last iteration of the envelope data sequenceq-1(n) smoothing to obtain current smoothing result yq(n)
Figure GDA0002462721590000131
And 8-5: judging whether the maximum number of iterative smoothing times of the envelope is reached, namely judging that q is less than or equal to I2If yes, entering the step 8-6, otherwise jumping to the step 8-8;
and 8-6: respectively calculating last smoothing result y of envelope data sequenceq-1(n) and current smoothing result yq(n) and the initial envelope data sequence y0Residual sum of squares Ψ of (n)q-1And ΨqRespectively is as follows:
Figure GDA0002462721590000132
and 8-7: determine | Ψq-1q|≤2ΨqIf yes, entering the step 8-8, otherwise returning to the step 8-4;
and 8-8: let y (n) be yq(N), N is 0,1,2 …, N-1, resulting in a sequence of signal envelope data:
y(n),n=0,1,2…,N-1。
wherein, the step 8-1 can be realized by fast Fourier transform, and the calculation efficiency of the algorithm is improved.
In the embodiment of the invention, the simulation rectangular envelope pulse signal model is as follows:
Figure GDA0002462721590000133
where a is the amplitude of the signal and,
Figure GDA0002462721590000134
to an initial phase, τ0For the pulse signal arrival time, τ is the pulse width, T is the received signal duration, f1Is the signal starting frequency, f2For signal termination frequency, when f2=f1Time is a narrow band signal, otherwise, is a wide band signal, w (t) is a mean value of 0, and variance is σ2White gaussian noise, variance σ2Is determined by the signal-to-noise ratio SNR: SNR of 10log10[A2/(2σ2)]。
At a sampling frequency fsThe pulse signal is subjected to discrete sampling to obtain a pulse signal sampling data sequence:
Figure GDA0002462721590000141
wherein n is0=int(τ0fs),M=int(τfs),N=int(Tfs)。
Example 1:
the simulation signal parameters are respectively set as: signal amplitude a 1, initial phase
Figure GDA0002462721590000142
Pulse width tau is 0.512s, received signal duration T is 1.024s, and pulse signal arrival time tau00.256s, signal start frequency f1400Hz, signal termination frequency f2400Hz, i.e. the artificial signal is a narrow-band pulse signal, the sampling frequency fs2000Hz, 2048 observation data sequence points N, and 3dB SNR.
According to the step (2), calculating the power spectrum P (k) of the data sequence x (n), as shown in FIG. 2.
In the step (3), setting a power spectrum iteration smoothing maximum iteration number threshold I110, precision control threshold10.01, iterative smoothing window length M1Envelope iteration smoothing maximum iteration number threshold I of 525, the precision control threshold20.02, iterative smoothing window length M2Narrow band signal bandwidth factor ω 1110.05, bandwidth factor ω of the wideband signal2And (4) judging that the main lobe characteristic parameter threshold η of the narrow-band pulse and the wide-band pulse is 2 when the main lobe characteristic parameter threshold is 0.1.
According to the step (4), performing iterative smoothing processing on p (k) to obtain a smoothed power spectrum s (k), as shown in fig. 3, it can be seen from comparing fig. 2 and fig. 3 that: after iterative smoothing processing, the main lobe of the power spectrum of the narrow-band pulse signal is obviously widened, and the whole frequency band spectrum is smoother.
According to the step (5), extracting the original power spectrum and the main lobe characteristic parameter gamma of the smoothed power spectrumpAnd gammaSAre respectively as
γp=45.6772,γS=1.4108
The adaptive filter parameters h (k) generated according to step (6) are shown in fig. 4.
According to the steps (7) and (8), the extracted pulse signal envelope is shown in fig. 5, and as can be seen from fig. 5, the front and back edge features of the extracted narrow-band pulse signal envelope are obvious, so that a basis can be provided for further extracting the arrival time and the end time of the narrow-band pulse signal.
Example 2:
the simulation signal parameters are respectively set as: signal amplitude a 2, initial phase
Figure GDA0002462721590000151
Pulse width τ is 0.5s, received signal duration T is 1s, and pulse signal arrival time τ is00.25s, signal start frequency f1350Hz, signal termination frequency f2450Hz, i.e. the simulated signal is a broadband pulse signal, the sampling frequency fs2048Hz, 2048 points N of the observation data sequence, and 6dB SNR of the signal-to-noise ratio.
According to the step (2), calculating the power spectrum P (k) of the data sequence x (n), as shown in FIG. 6.
In the step (3), setting a power spectrum iteration smoothing maximum iteration number threshold I110, precision control threshold10.01, iterative smoothing window length M17, envelope iteration smoothing maximum iteration number threshold I210, precision control threshold20.02, iterative smoothing window length M213, narrow-band signal bandwidth factor ω10.05, bandwidth factor ω of the wideband signal2And (4) judging that the main lobe characteristic parameter threshold η of the narrow-band pulse and the wide-band pulse is 2 when the main lobe characteristic parameter threshold is 0.1.
According to the step (4), performing iterative smoothing processing on p (k) to obtain a smoothed power spectrum s (k), as shown in fig. 7, comparing fig. 6 and fig. 7, it can be seen that: after iterative smoothing processing, the main lobe change of the power spectrum of the broadband pulse signal is not obvious, and the whole frequency band spectrum is smoother.
According to the step (5), extracting the original power spectrum and the main lobe characteristic parameter gamma of the smoothed power spectrumpAnd gammaSRespectively as follows:
γp=1.1081,γS=1.1974
the adaptive filter parameters h (k) generated according to step (6) are shown in fig. 8.
According to the steps (7) and (8), the extracted pulse signal envelope is shown in fig. 9, and as can be seen from fig. 9, the front and rear edge features of the extracted wideband pulse signal envelope are obvious, so that a basis can be provided for further extracting the arrival time and the end time of the wideband pulse signal.

Claims (8)

1. A method for extracting a self-adaptive filtering envelope of a pulse signal frequency domain is characterized by comprising the following steps:
(1) acquiring a pulse signal sampling data sequence x (N) to be processed, wherein N is 0,1, … and N-1, and N is the number of sampling points corresponding to the pulse width length of a detected pulse signal and is an integer power of 2;
(2) calculating an original power spectrum P (k) of the pulse signal according to the data sequence x (n);
(3) carrying out parameter initialization on iterative smoothing and filter bandwidth;
(4) carrying out iterative smoothing processing on the original power spectrum P (k) to obtain a smoothed power spectrum S (k);
(5) extracting main lobe characteristic parameters of the original power spectrum P (k) and the smooth power spectrum S (k) respectively;
(6) generating adaptive filtering parameters h (k);
(7) performing adaptive filtering on the discrete Fourier transform result X (k) of the pulse signal according to the adaptive filtering parameter H (k);
(8) extracting a pulse signal envelope y (n);
in the step (5), the main lobe characteristic parameters of the original power spectrum p (k) and the smoothed power spectrum s (k) are respectively extracted, which specifically includes the following steps:
step 5-1: searching discrete frequency index k corresponding to maximum value of smoothed power spectrum S (k)M
Figure FDA0002462721580000011
Wherein the content of the first and second substances,
Figure FDA0002462721580000012
representing that the discrete frequency index corresponding to the maximum value of S (k) is searched within the range of 1 ≦ k ≦ N/2-1;
step 5-2: subtracting half of the maximum value from s (k) to obtain a half-power difference power spectrum z (k):
Z(k)=S(k)-S(kM)/2,k=0,1,2…,N/2-1
step 5-3: searching for the discrete frequency index k of the first zero crossing on the left of the Z (k) peakLThe searching process comprises the following steps:
step 5-3-1: initializing a frequency start search index k of a first zero crossing on the left side of a Z (k) peak valuel=kM-1;
Step 5-3-2: judgment of kl0 or Z (k)l) If the result is less than 0, jumping to the step 5-3-4 if the result is true, and otherwise, entering the step 5-3-3;
step 5-3-3: let kl=kl1 and returning to step 5-3-2;
step 5-3-4: let kL=klObtaining the discrete frequency index k of the first zero-crossing point on the left side of the Z (k) peak valueL
Step 5-4: searching for the first zero-crossing discrete frequency index k to the right of the Z (k) peakRThe searching process comprises the following steps:
step 5-4-1: initializing a frequency start search index k of a first zero crossing on the right side of a Z (k) peak valuer=kM+1;
Step 5-4-2: judgment of krN/2-1 or Z (k)r) If the value is less than 0, jumping to the step 5-4-4 if the value is less than 0, otherwise, entering the step 5-4-3;
step 5-4-3: let kr=kr+1 and returning to step 5-4-2;
step 5-4-4: let kR=krObtaining the discrete frequency index k of the first zero-crossing point on the right side of the Z (k) peak valueR
Step 5-5: calculating the primary power spectrum P (k) main lobe characteristic parameter gammap
Figure FDA0002462721580000021
And 5-6: calculating a smoothed power spectrum S (k) main lobe characteristic parameter gammaS
Figure FDA0002462721580000022
2. The method according to claim 1, wherein in step (1), the processed pulse signal sample data sequence is obtained by the following method: receiving real-time acquisition data of N sampling points from a sensor as a data sequence x (N) to be processed; or extracting the N sampling point data containing the whole pulse signal from the memory as the data sequence x (N) to be processed.
3. The method for extracting adaptive filtering envelopes in frequency domain of impulse signal as claimed in claim 1, wherein in the step (2), the original power spectrum p (k) of the impulse signal is calculated for the data sequence x (n) by the following method: performing fast Fourier transform on the data sequence x (n), and calculating to obtain discrete Fourier transform X (k) and an original power spectrum P (k) of the data sequence, wherein the method specifically comprises the following steps:
step 2-1: calculating the discrete Fourier transform of x (n):
Figure FDA0002462721580000023
where k is the discrete frequency index of X (k) and j represents the unit of an imaginary number, i.e.
Figure FDA0002462721580000031
This formula is implemented by fast fourier transform;
step 2-2: calculating a power spectrum for x (n) according to x (k):
Figure FDA0002462721580000032
wherein k is the discrete frequency index of the original power spectrum p (k), | represents the modulo value operation, and | x (k) | is the signal amplitude spectrum.
4. The method for extracting the adaptive filtering envelope in the frequency domain of the impulse signal as claimed in claim 1, wherein in the step (3), the parameters of the iterative smoothing and the filter bandwidth are initialized by the following method: setting a maximum iteration number threshold of power spectrum iterative smoothing, wherein the precision control threshold and the iterative smoothing window length are respectively I11And M1The maximum iteration number threshold, the precision control threshold and the iteration smoothing window length of the initial envelope iteration smoothing are respectively I22And M2Bandwidth coefficient omega of narrow-band signal1And bandwidth factor omega of broadband signal2Judging the characteristic parameter threshold η of the main lobe of the narrow-band pulse and the wide-band pulse, wherein I1And I2A value is a positive integer greater than 2,1and2the value is positive number less than 1, the length of the iterative smoothing window of the power spectrum satisfies 3 and is less than or equal to M1An odd number less than or equal to N/2-2, and an envelope iteration smoothing window length of 5 less than or equal to M2Odd number, omega, not more than 491To satisfy 0 ≤ ω1Any number, omega, not more than 0.052To satisfy 0.05 < omega2Any number not less than 0, η is any number satisfying 1.5 not less than η not more than 2.5.
5. The method for extracting adaptive filtering envelopes of a pulse signal frequency domain according to claim 1, wherein in the step (4), an iterative smoothing process is performed on an original power spectrum p (k) by using a method as follows to obtain a smoothed power spectrum s (k), and the method specifically includes the following steps:
step 4-1: the iterative smoothing times i of the initialized power spectrum is 0, and the iterative smoothing result S of the power spectrum0(k) Comprises the following steps:
S0(k)=P(k),k=0,1,2…,N/2-1
step 4-2: let the iteration number i be i +1, and smooth the last iteration of the power spectrum to obtain a result Si-1(k) Performing smoothing to obtain current smoothing result Si(k):
Figure FDA0002462721580000041
Wherein the content of the first and second substances,M1iteratively smoothing the window length for the power spectrum;
step 4-3: judging whether the maximum iterative smoothing times of the power spectrum is reached, namely judging that I is less than or equal to I1If yes, entering the step 4-4, otherwise jumping to the step 4-6; wherein I1A maximum iteration number threshold for power spectrum iterative smoothing;
step 4-4: respectively calculating last smoothing results S of power spectrumi-1(k) And current smoothing result Si(k) Sum of squares J of residuals from original power spectrum p (k)i-1And JiRespectively is as follows:
Figure FDA0002462721580000042
and 4-5: judgment of | Ji-1-Ji|≤1JiIf yes, entering the step 4-6, otherwise returning to the step 4-2; wherein1A threshold is controlled for power spectrum iterative smoothing precision;
and 4-6: let S (k) become Si(k) And k is 0,1,2 …, N/2-1, and the smoothed power spectrum s (k) is obtained.
6. The method according to claim 1, wherein in the step (6), the adaptive filtering parameters h (k) are generated by using the following method, k being 0,1,2 …, N-1, and specifically comprising the following steps:
step 6-1: initializing adaptive filtering parameters:
H(k)=0,k=0,1,2…,N-1
step 6-2: respectively calculating the discrete frequency indexes T of the start and the end of the transition band of the narrow-band pulse signal filterN1And TN2And a passband start and end discrete frequency index PN1And PN2
TN1=max(round[(1-0.5ω1)kM],0),TN2=min(round[(1+0.5ω1)kM],N/2-1)
PN1=max(kM-1,TN1),PN2=min(kM+1,TN2)
Wherein, ω is1Is a narrow band signal bandwidth coefficient, round [ ]]Representing a rounding operation, max (,) and min (,) representing a large and small value operation, respectively, kMRepresenting the discrete frequency index corresponding to the maximum value of the smoothed power spectrum S (k);
step 6-3: respectively calculating the start and end discrete frequency indexes T of the transition band of the broadband pulse signal filterW1And TW2And a passband start and end discrete frequency index PW1And PW2
TW1=max(round[(1-0.5ω2)kM],0),TW2=min(round[(1+0.5ω2)kM],N/2-1)
PW1=max(kL-1,TW1),PW2=max(kR+1,TW2)
Wherein, ω is2For the bandwidth factor, k, of the broadband signalLDiscrete frequency index for the first zero crossing point on the left side of peak value of half-power difference power spectrum Z (k)RIndexing the first zero-crossing discrete frequency on the right side of a half-power difference power spectrum Z (k) peak value;
step 6-4: obtaining discrete frequency index T from adaptive filter transition band start and end according to signal type1And T2And a passband start and end discrete frequency index P1And P2
Namely, max (gamma)ps)/min(γps) If > η is true, if true, the signal is judged to be a narrow-band pulse signal and the following steps are carried out:
T1=TN1,T2=TN2,P1=PN1,P2=PN2
otherwise, judging the signal to be a broadband pulse signal and ordering:
T1=TW1,T2=TW2,P1=PW1,P2=PW2
wherein, γpIs the main lobe characteristic parameter, gamma, of the original power spectrum P (k)STo smooth the power spectrum S: (k) η is a main lobe characteristic parameter threshold for judging the narrow-band pulse and the wide-band pulse;
step 6-5: modifying the transition band part of the adaptive filter parameters H (k) so that k satisfies T1≤k≤T2The part of (A) is:
Figure FDA0002462721580000051
step 6-6: modifying the passband portion of the adaptive filter parameters H (k) such that k satisfies P1≤k≤P2The part of (A) is: h (k) ═ 1.
7. The method for extracting adaptive filtering envelopes in frequency domain of impulse signals as claimed in claim 1, wherein in the step (7), the discrete fourier transform result x (k) of the impulse signals is adaptively filtered by the following method to obtain the filtered discrete fourier transform y (k):
Y(k)=X(k)H(k),k=0,1,2…,N-1
where h (k) is an adaptive filtering parameter.
8. The method for extracting the envelope of the pulse signal by frequency domain adaptive filtering according to claim 7, wherein in the step (8), the envelope y (N) of the pulse signal is extracted by the following method, where N is 0,1 …, and N-1, and the method specifically includes the following steps:
step 8-1: performing inverse discrete Fourier transform on the filtered discrete Fourier transform Y (k) to obtain a pulse signal filtered time-domain complex signal yc(n):
Figure FDA0002462721580000061
Step 8-2: taking a complex number yc(n) obtaining the initial envelope data sequence y of the pulse signal by the modulus value0(n):
y0(n)=|yc(n)|,n=0,1…,N-1
Step 8-3: initializing the iterative smoothing times q of the envelope data to be 0;
step 8-4: let the iteration number q be q +1 and smooth the result y for the last iteration of the envelope data sequenceq-1(n) smoothing to obtain current smoothing result yq(n):
Figure FDA0002462721580000062
Wherein M is2Iterating the smoothing window length for the initial envelope;
and 8-5: judging whether the maximum iterative smoothing times of the envelope is reached, namely judging that q is less than or equal to I2If yes, entering the step 8-6, otherwise jumping to the step 8-8; wherein, I2A maximum iteration number threshold for initial envelope iteration smoothing;
and 8-6: respectively calculating last smoothing result y of envelope data sequenceq-1(n) and current smoothing result yq(n) and the initial envelope data sequence y0Residual sum of squares Ψ of (n)q-1And ΨqRespectively is as follows:
Figure FDA0002462721580000071
and 8-7: determine | Ψq-1q|≤2ΨqIf yes, entering the step 8-8, otherwise returning to the step 8-4;
wherein the content of the first and second substances,2a threshold is controlled for the initial envelope iteration smoothing precision;
and 8-8: let y (n) be yq(N), N is 0,1,2 …, N-1, resulting in a sequence of signal envelope data:
y(n),n=0,1,2…,N-1。
CN201810325134.9A 2018-04-11 2018-04-11 Adaptive filtering envelope extraction method for pulse signal frequency domain Active CN108680786B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810325134.9A CN108680786B (en) 2018-04-11 2018-04-11 Adaptive filtering envelope extraction method for pulse signal frequency domain

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810325134.9A CN108680786B (en) 2018-04-11 2018-04-11 Adaptive filtering envelope extraction method for pulse signal frequency domain

Publications (2)

Publication Number Publication Date
CN108680786A CN108680786A (en) 2018-10-19
CN108680786B true CN108680786B (en) 2020-07-31

Family

ID=63800312

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810325134.9A Active CN108680786B (en) 2018-04-11 2018-04-11 Adaptive filtering envelope extraction method for pulse signal frequency domain

Country Status (1)

Country Link
CN (1) CN108680786B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110897621B (en) * 2019-10-21 2022-10-28 思澜科技(成都)有限公司 Physiological signal extraction method and system
CN110852201B (en) * 2019-10-28 2023-08-11 东南大学 Pulse signal detection method based on multi-pulse envelope spectrum matching
CN111812404B (en) * 2020-09-14 2021-01-01 湖南国科雷电子科技有限公司 Signal processing method and processing device
CN113191317B (en) * 2021-05-21 2022-09-27 江西理工大学 Signal envelope extraction method and device based on pole construction low-pass filter
CN113848391B (en) * 2021-09-26 2024-04-09 成都华日通讯技术股份有限公司 Pulse signal detection and extraction method
CN113791391B (en) * 2021-10-15 2024-03-26 西安黄河机电有限公司 Interference characteristic parameter identification method, device, storage and equipment
CN114499702B (en) * 2022-03-28 2022-07-12 成都锢德科技有限公司 Portable real-time signal acquisition, analysis and recognition system
CN114826541B (en) * 2022-04-08 2023-12-26 西南石油大学 Low-sample-number signal center frequency estimation method and system
CN116032310B (en) * 2023-02-20 2023-07-04 西安瀚博电子科技有限公司 Signal self-adaptive detection reconstruction method based on channelized filtering

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100525713C (en) * 2005-06-10 2009-08-12 深圳迈瑞生物医疗电子股份有限公司 Method for extracting envelope curve of sound spectrogram
CN100496409C (en) * 2005-08-02 2009-06-10 深圳迈瑞生物医疗电子股份有限公司 Automatic detection method of frequency spectrum Doppler blood flow velocity
US7877252B2 (en) * 2007-05-18 2011-01-25 Stmicroelectronics S.R.L. Automatic speech recognition method and apparatus, using non-linear envelope detection of signal power spectra
CN103393435B (en) * 2013-07-31 2015-01-07 深圳市理邦精密仪器股份有限公司 Method and device for acquiring fetal heart sound signal envelope
KR20170045990A (en) * 2015-10-20 2017-04-28 울산대학교 산학협력단 Method for machine condition monitoring
CN106073822B (en) * 2016-05-31 2018-11-09 青岛惠尔医疗科技有限公司 A kind of superonic spectrum adaptive noise estimation envelope extraction method
CN106153333B (en) * 2016-06-29 2018-06-22 潍坊学院 A kind of envelope Analysis Method based on wavelet decomposition filtering
CN106208974A (en) * 2016-08-30 2016-12-07 厦门信同信息技术有限公司 A kind of design method of digital power amplifier based on adaptive digital envelop following technology
CN107064629B (en) * 2017-06-07 2019-07-12 东南大学 A kind of comprehensive simple signal frequency estimating methods of the segmentation estimated based on frequency relative bias

Also Published As

Publication number Publication date
CN108680786A (en) 2018-10-19

Similar Documents

Publication Publication Date Title
CN108680786B (en) Adaptive filtering envelope extraction method for pulse signal frequency domain
CN109446928B (en) Signal noise reduction method based on variational modal decomposition and minimum mean square error adaptive filter
US7492814B1 (en) Method of removing noise and interference from signal using peak picking
CN106597408B (en) High-order PPS signal parameter estimation method based on time-frequency analysis and instantaneous frequency curve fitting
CN108469602B (en) Pulse signal type automatic discrimination method based on spectral feature extraction
US7676046B1 (en) Method of removing noise and interference from signal
CN108594177A (en) Based on radar signal modulation system analysis method, the signal processing system for improving HHT
WO2017147255A1 (en) Estimation and mitigation of swept-tone interferers in frequency-hopped systems
CN110808057A (en) Voice enhancement method for generating confrontation network based on constraint naive
CN108333568B (en) Broadband echo Doppler and time delay estimation method based on Sigmoid transformation in impact noise environment
Roenko et al. An overview of the adaptive robust DFT
CN109104215B (en) Frequency hopping signal code element rate blind estimation method based on wavelet transformation
CN103915102A (en) Method for noise abatement of LFM underwater sound multi-path signals
CN109460614B (en) Signal time-frequency decomposition method based on instantaneous bandwidth
CN111600665B (en) SAGE channel parameter estimation method based on ARIMA filter
Al-Khafaji et al. Frequency estimation of FM signals under non-Gaussian and colored noise
CN113030880A (en) ADMM-based signal estimation method for intermittent sampling forwarding type interference
CN106997766B (en) Homomorphic filtering speech enhancement method based on broadband noise
CN112230200A (en) Improved combined noise reduction method based on laser radar echo signals
CN113238206B (en) Signal detection method and system based on decision statistic design
Cang et al. Robust Deconvolution of Underwater Acoustic Channels Corrupted by Impulsive Noise
Sharma et al. Implementation of Sound Effects on Audio Signals
Zaw et al. Performance Comparison of Noise Detection and Elimination Methods for Audio Signals
CN111178325A (en) Signal extraction and processing method
CN114722345A (en) Method for estimating chirp rate of linear frequency modulation signal based on deskew processing and time delay correlation

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