CN114488208B - Beidou signal anti-interference method combining empirical wavelet and SPWVD conversion - Google Patents
Beidou signal anti-interference method combining empirical wavelet and SPWVD conversion Download PDFInfo
- Publication number
- CN114488208B CN114488208B CN202210145552.6A CN202210145552A CN114488208B CN 114488208 B CN114488208 B CN 114488208B CN 202210145552 A CN202210145552 A CN 202210145552A CN 114488208 B CN114488208 B CN 114488208B
- Authority
- CN
- China
- Prior art keywords
- signal
- interference
- omega
- function
- frequency
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000006243 chemical reaction Methods 0.000 title description 6
- 230000009466 transformation Effects 0.000 claims abstract description 15
- 238000001228 spectrum Methods 0.000 claims description 17
- 238000001514 detection method Methods 0.000 claims description 14
- 230000003595 spectral effect Effects 0.000 claims description 12
- 230000001131 transforming effect Effects 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 7
- 125000003275 alpha amino acid group Chemical group 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 3
- 238000005314 correlation function Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 5
- 238000000354 decomposition reaction Methods 0.000 abstract description 3
- 238000012545 processing Methods 0.000 abstract description 3
- 238000012216 screening Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 4
- 238000009499 grossing Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000002776 aggregation Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
-
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
- G01S19/215—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service issues related to spoofing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/148—Wavelet transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
-
- 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
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Computer Networks & Wireless Communication (AREA)
- Discrete Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention discloses a Beidou signal anti-interference method combining empirical wavelet and SPWVD transformation, which comprises the following steps: 1. the complex interference signals are represented by analysis; 2. the analysis signals are decomposed into empirical wavelet functions which are arranged according to the frequency sequence, so that the defect that a conventional wavelet base is difficult to construct is overcome, and the mode aliasing phenomenon during the decomposition of complex interference signals is effectively avoided; 3. separating and screening useful information and interference signals for the signals after the empirical wavelet transformation; 4. for the cross terms between existing identical frequency signal components, a smooth pseudo Wigner-Ville distribution (SPWVD) transform is employed for further cancellation. The invention can effectively inhibit the influence of interference and noise of the cross item, keeps good time-frequency positioning characteristics, can be widely applied to GPS, GLONASS, galileo and Beidou navigation satellite signal processing, and has wide application prospect.
Description
Technical Field
The invention belongs to the field of satellite navigation positioning, and particularly relates to a Beidou receiver interference suppression method combining Empirical Wavelet Transform (EWT) and smooth pseudo Wigner-Ville transform.
Background
Beidou satellite navigation technology has become an increasingly important information technology in civil and military applications, and the characteristic of easy interference is exposed in the application field of complex interference environments nowadays. Therefore, it is important to study the anti-interference technology of the Beidou satellite navigation system receiver.
The main three sources causing the interference of the Beidou receiver are: firstly, transmission path interference is caused by interference experienced by Beidou navigation signals propagating from space to ground, and mainly comprises tropospheric interference, ionospheric interference and multipath interference; secondly, non-malicious electromagnetic wave interference is interference generated by a communication system with similar Beidou navigation frequency bands, such as a mobile communication system and a radar system, and the interference is generally weak, but the electromagnetic interference in a specific scene is strong, such as the electromagnetic interference in a wind power plant and a transmission line environment is a main interference source; thirdly, artificial malicious interference, which refers to hold-down interference and deception interference released by enemies in a battlefield environment. The invention is applied to civil scenes, so that the third artificial malicious interference is not considered.
The Direct Sequence Spread Spectrum (DSSS) scheme is a scheme adopted by most satellite navigation systems at present, and can spread the received navigation signal power to a wider bandwidth, but the broadcasting distance of the navigation satellite from the ground signal is too far, so that the signal power reaching the ground end is very weak, and even if the navigation signal adopts the direct spread spectrum technology, the power intensity is weaker than the intensity of the background normal noise.
The methods commonly used in the navigation receiving interference detection mainly include an Automatic Gain Control (AGC) method, a time domain method, a frequency domain method and a wavelet transformation method. The AGC adjusts the navigation input signal level to the range of an analog-to-digital converter (ADC), and performs interference detection by monitoring the AGC level; however, in an environment where the interference is weak, the detection performance of the AGC method may be greatly reduced. For the time domain approach, this approach may be implemented at the digital Intermediate Frequency (IF) level after the ADC at the front end of the navigation receiver. This time domain approach is only effective for narrowband RFI sources, as wideband interference is difficult to distinguish from thermal noise. For frequency domain methods, due to their spectral characteristics, narrowband carrier interference is typically detected. With respect to navigation reception interference detection using time domain methods or frequency domain methods, they cannot fully describe the nature and characteristics of the time-varying interference present in the received navigation signal. The wavelet transform method is a time-frequency domain transform method, which can increase the time scale while analyzing the signal frequency, and push the time-frequency analysis to research climax, however, the wavelet basis of the common wavelet transform is difficult to construct and has the defects of time shift and frequency aliasing. Thus, interference detection has become a critical role in navigation applications.
Disclosure of Invention
The invention aims to solve the defects of the prior art, and provides a Beidou signal anti-interference method combining empirical wavelet and SPWVD conversion, so that the influence of cross term interference and noise can be effectively restrained, and good time-frequency positioning characteristics are reserved, thereby solving the problems of time-frequency compromise, cross term interference, energy disaggregation and the like of navigation receiving signals in scanning interference detection.
The invention solves the technical problems by adopting the following technical scheme:
the invention relates to a Beidou signal anti-interference method combining empirical wavelet and SPWVD conversion, which is characterized by comprising the following steps:
step 1: analyzing and transforming the original continuous time signal received by the Beidou receiver by utilizing the (1) to obtain an analysis signal y a (t):
In the formula (1), j is a virtual root unit, and y (t) represents an original continuous time signal;representing y (t) a Hilbert transformed signal;
step 2: resolving the signal y a (t) performing empirical wavelet transformation to obtain modal function components with different frequencies and mutually independent;
step 2.1: for the analytic signal y a (t) Fourier transforming to obtain corresponding frequency spectrum y a (ω);
Step 2.2: in the frequency spectrum y a Finding M local maximum points on (omega), and arranging the M local maximum points in a descending order to obtain ordered local maximum points;
according to the ordered local maximum value points, the local maximum value points are in [0, pi ]]Up-to the spectrum y a (omega) dividing into N consecutive parts to obtain N divided spectral bandsAnd forming n+1 endpoints;
if M is greater than N, reserving the first N-1 local maximum points after sequencing;
if M < N, resetting the parameter N such that n=m;
let omega n Representing boundary points of adjacent n-1 th and n-th spectrum segments, letting the initial boundary point omega 0 =0, end boundary point ω N =pi; n=1, …, N; then the nth spectral band Λ n =[ω n-1 ,ω n ]Thereby obtainingWherein omega n-1 Representing boundary points of adjacent n-2 th and n-1 th spectral bins;
step 2.3: an expression of an empirical scale function and a wavelet function expression of an empirical wavelet are respectively constructed using the expression (2) and the expression (3):
in the formula (2) and the formula (3), phi n (omega) represents the nth empirical scale function, ψ n (ω) represents the nth empirical wavelet function, γ represents a constant between 0 and 1, and β (·) is a value of [0,1 ]]A function of (a); omega n+1 Representing boundary points of adjacent nth and n+1th spectral bins;
step 2.4: defining detail correlation coefficients by means of equations (4) and (5), respectivelyApproximate correlation coefficient
In the formula (4), "-" represents a conjugate operation, and when n=1, φ 1 (t) is defined by phi n (ω) inverse fourier transforming to obtain a function value; psi phi type n (t) consists of psi n (omega) function obtained by inverse Fourier transform, y a (t) is defined by y a (ω) inverse fourier transforming the resulting function; f (F) -1 (. Cndot.) represents the inverse Fourier transform;
step 2.5: by using (6) the analysis signal y a (t) reconstructing to obtain an analysis signal y a The empirical mode components of (t) are shown in the formulas (7) and (8):
in equation (6), the convolution operation,is->Fourier transform form of>Represents n=0 +.>Is a value of (2);
in formula (7), y 0 (t) represents an amplitude modulation-frequency modulation single component when n=0,representation->Is an inverse fourier transform form of (a);
in formula (8), y n (t) represents amplitude-modulated-frequency components of different frequencies, n=1, 2,3, N;
step 3: the analysis signal y is compared with the reference signal y by the method (9) a Each modal component y of (t) n (t) performing smooth pseudo-Wigner-Ville transformation to obtain a smooth pseudo-Wigner-Ville time-frequency distribution expression
In the formula (9), the amino acid sequence of the compound,and->Representing y n Transient correlation function of (t),>representation->Conjugation of (2); g (s-t) is a smoothing window function in the frequency axis direction, h (τ) is a smoothing window function in the time domain direction, and g (0) =h (0) =1; s and τ are integral variables;
step 4: using (10) to divide each mode component y n Linear superposition of the smoothed pseudo Wigner-Ville analysis results of (t) to obtain an analysis signal y a (t) smooth pseudo-Wigner-Ville distribution
Step 5: for the analytic signal y a Sliding pseudo-Wigner-Ville distribution of (t)And performing interference detection, and filtering the obtained interference signal by adopting a notch IIR filter, so as to obtain the sliding pseudo Wigner-Ville distribution after filtering.
Compared with the prior art, the invention has the beneficial effects that:
1. the navigation receiving signals subjected to the frequency sweep interference are represented by analysis and then subjected to empirical wavelet transformation, so that the defect that a wavelet basis is difficult to construct in conventional wavelet transformation is overcome, the modal aliasing phenomenon during the decomposition of complex interference signals is effectively avoided, the spectral characteristics of the interference signals and the useful spectral characteristics of the navigation signals can be clearly distinguished, and the method can be widely applied to the processing of GPS, GLONASS, galileo and Beidou navigation satellite signals and has wide application prospect.
2. The invention adopts smooth pseudo Wigner-Ville distribution (SPWVD) time-frequency transformation, can further eliminate the cross terms between the existing signal components with the same frequency while enhancing the self terms, and effectively suppresses the interference and noise influence of the cross terms.
3. The invention adopts the combined empirical wavelet and SPWVD conversion method, improves the energy aggregation performance in the navigation signal interference detection, overcomes the time-frequency compromise problem, and the influence of cross term interference and noise, and obviously improves the interference detection performance.
Drawings
FIG. 1 is a flow chart of a Beidou signal anti-interference method combining empirical wavelet and SPWVD conversion;
FIG. 2a is a WVD scale plot of LFM signals in the presence of AWGN in the prior art;
FIG. 2b is a WVD profile of the LFM signal at AWGN in the prior art;
FIG. 3 is an experimental test protocol of the present invention;
FIG. 4a is a WVD scale plot of a GPS L1-C/A signal with frequency sweep interference in the prior art;
FIG. 4b is a WVD profile of a GPS L1-C/A signal with frequency sweep disturbance in the prior art;
FIG. 5a is a graph of the joint empirical wavelet and SPWVD transform scales of a swept-disturbed GPS L1-C/A signal of the present invention;
FIG. 5b is a graph of the combined empirical wavelet and SPWVD transform profile of a swept-disturbed GPS L1-C/A signal of the present invention.
Detailed Description
In this embodiment, a method for suppressing interference of a beidou receiver combining empirical wavelet transform and smooth pseudo-Wigner-Ville transform, as shown in fig. 1, includes the following steps:
step 1: resolving and transforming the original continuous time signal received by the Beidou receiver by utilizing the (1),obtain the analysis signal y a (t):
In the formula (1), j is a virtual root unit, and y (t) represents an original continuous time signal;representing y (t) a Hilbert transformed signal;
step 2: will resolve the signal y a (t) performing empirical wavelet transformation to obtain modal function components with different frequencies and mutually independent;
step 2.1: for the analysis signal y a (t) Fourier transforming to obtain corresponding frequency spectrum y a (ω);
Step 2.2: in the frequency spectrum y a Finding M local maximum points on (omega), and arranging the M local maximum points in a descending order to obtain ordered local maximum points;
according to the ordered local maximum value points, the local maximum value points are in [0, pi ]]Upper pair spectrum y a (omega) dividing into N consecutive parts to obtain N divided spectral bandsAnd forming n+1 endpoints;
if M is greater than N, reserving the first N-1 local maximum points after sequencing;
if M < N, resetting the parameter N such that n=m;
let omega n Representing the boundary points of the n-1 th and n th adjacent spectrum segments, wherein the selection of the boundary points is determined by the median between two adjacent local maximum points, so that the initial boundary point omega 0 =0, end boundary point ω N =pi; n=1, …, N; then the nth spectral band Λ n =[ω n-1 ,ω n ]Thereby obtainingWherein omega n-1 Represents adjacent n-2 th and n-1 boundary point of spectrum segment;
step 2.3: an expression of an empirical scale function and a wavelet function expression of an empirical wavelet are respectively constructed using the expression (2) and the expression (3):
in the formula (2) and the formula (3), phi n (omega) represents the nth empirical scale function, ψ n (ω) represents the nth empirical wavelet function, γ represents a constant between 0 and 1, and β (·) is a value of [0,1 ]]Any function above; omega n+1 Representing boundary points of the n-th and n+1-th adjacent spectrum segments, the equations (2) and (3) correspond to the functions of a low-pass filter and a band-pass filter, and each signal component can be obtained after filtering the analysis signal.
Step 2.4: defining detail correlation coefficients by means of equations (4) and (5), respectivelyApproximate correlation coefficient
In the formula (4), "-" represents a conjugate operation, and when n=1, φ 1 (t) is defined by phi n (ω) inverse fourier transforming to obtain a function value; psi phi type n (t) consists of psi n (omega) function obtained by inverse Fourier transform, y a (t) is defined by y a (omega) Fourier inverseTransforming the obtained function; f (F) -1 (. Cndot.) represents the inverse Fourier transform;
step 2.5: by using (6) the analysis signal y a (t) reconstructing to obtain an analysis signal y a The empirical mode components of (t) are shown in the formulas (7) and (8):
signal y a (t) decomposing to obtain amplitude modulation-frequency modulation single component components with frequencies from low to high; in equation (6), the convolution operation,is->Fourier transform form of>Represents n=0 +.>Is a value of (2);
in formula (7), y 0 (t) represents an amplitude modulation-frequency modulation single component when n=0,representation->Is an inverse fourier transform form of (a);
in formula (8), y n (t) represents amplitude-modulated-frequency components of different frequencies, n=1, 2,3, N;
step 3: if a signal contains multiple components in the TF plane, its Wigner-Ville transform (WVD) can be affected by spurious features containing cross terms that occur between each pair of automatic terms. For example, a noisy chirp (Linearly Frequency Modulated, LFM) signal in the presence of AWGN is analyzed, with the SNR set to 3dB. The WVD of the noise LFM signal is provided in fig. 2a, and correspondingly the calculated profile of the WVD is given in fig. 2 b. From the figure, the TF peak representing the noise LFM signal can be observed, but there is a severe cross term in the TF plane, and no energy at all is expected. The presence of cross terms of the analysis signal in the TF plane does not have any physical significance, which makes correct signal interpretation very difficult, which is also a major drawback of the WVD method. The invention thus employs SPWVD methods.
Step 3.1: the analysis signal y is compared with the reference signal y by the method (9) a Each modal component y of (t) n (t) performing smooth pseudo-Wigner-Ville transformation to obtain a smooth pseudo-Wigner-Ville time-frequency distribution expression
In the formula (9), the amino acid sequence of the compound,and->Representing y n Transient correlation function of (t),>representation->Conjugation of (2); g (s-t) is the frequency axisA smoothing window function of direction, h (τ) is a smoothing window function of time domain direction, and g (0) =h (0) =1; s and τ are both integral variables.
Step 4: using (10) to divide each mode component y n Linear superposition of the smoothed pseudo Wigner-Ville analysis results of (t) to obtain an analysis signal y a (t) smooth pseudo-Wigner-Ville distribution
Step 5: analysis signal y obtained by linear superposition a (t)And performing interference detection, and filtering interference signals contained in the interference detection by adopting a notch IIR filter so as to obtain the sliding pseudo Wigner-Ville distribution after filtering.
The experimental test scheme of the invention is shown in fig. 3, a software jammer is adopted to generate a sweep frequency interference signal for GNSS application, and the generated sweep frequency interference is added into GPS samples collected at the front end of a GNSS receiver, so that a GPS signal subjected to sweep frequency interference is generated.
To verify the effectiveness of the method, it was compared with the existing classical method WVD. FIG. 4a depicts the WVD of the GPS L1-C/A signal in the presence of scanning disturbances, with the profile of the WVD correspondingly shown in FIG. 4 b. Although a straight line can be seen on the TF plane representing the frequency modulation law of the scanning disturbance, due to the bilinear nature of WVD, severe cross terms can also be observed, which inevitably leads to errors in the instantaneous frequency estimation. Scanning the interference signal brings difficulty to the detection of the interference, so that the GNSS interference characteristic parameters cannot be extracted correctly by using WVD.
FIG. 5a is C/N 0 =40 dB, jnr=2 dB, the joint empirical wavelet and SPWVD time-frequency transformed scale plot of the swept-interfered GPS L1-C/a signal, fig. 5b is a profile plot. It can be seen that the cross terms are in the TF planeThe method has the advantages that effective inhibition is obtained, TF energy peaks are linearly distributed on a straight line, and the TF energy peaks represent that sweep frequency interference exists in the received GPS L1-C/A signals. The energy aggregation performance in GNSS interference detection is improved, so that the instantaneous frequency estimation of an interference signal is improved.
In summary, the invention provides a Beidou receiver interference suppression method combining empirical wavelet transformation and smooth pseudo Wigner-Ville transformation, which decomposes complex interference signals into empirical wavelet functions arranged according to frequency sequence through analysis representation, so that the defect that a conventional wavelet base is difficult to construct is overcome, and the modal aliasing phenomenon during the decomposition of the complex interference signals is effectively avoided. The signal after the empirical wavelet transformation separates and screens useful information from interfering signals, and for the cross terms between existing signal components of the same frequency, a smooth pseudo Wigner-Ville distribution (SPWVD) transformation is used for further cancellation. The method effectively suppresses the influence of interference and noise of the cross item, maintains good time-frequency positioning characteristics, can be widely applied to GPS, GLONASS, galileo and Beidou navigation satellite signal processing, and has wide application prospect.
Claims (1)
1. A Beidou signal anti-interference method combining empirical wavelet and SPWVD transformation is characterized by comprising the following steps:
step 1: analyzing and transforming the original continuous time signal received by the Beidou receiver by utilizing the (1) to obtain an analysis signal y a (t):
In the formula (1), j is a virtual root unit, and y (t) represents an original continuous time signal;representing y (t) a Hilbert transformed signal;
step 2: resolving the signal y a (t) performing empirical wavelet transform to obtainEach modal function component with different frequencies and mutually independent;
step 2.1: for the analytic signal y a (t) Fourier transforming to obtain corresponding frequency spectrum y a (ω);
Step 2.2: in the frequency spectrum y a Finding M local maximum points on (omega), and arranging the M local maximum points in a descending order to obtain ordered local maximum points;
according to the ordered local maximum value points, the local maximum value points are in [0, pi ]]Up-to the spectrum y a (omega) dividing into N consecutive parts to obtain N divided spectral bandsAnd forming n+1 endpoints;
if M is greater than N, reserving the first N-1 local maximum points after sequencing;
if M < N, resetting the parameter N such that n=m;
let omega n Representing boundary points of adjacent n-1 th and n-th spectrum segments, letting the initial boundary point omega 0 =0, end boundary point ω N =pi; n=1, …, N; then the nth spectral band Λ n =[ω n-1 ,ω n ]Thereby obtainingWherein omega n-1 Representing boundary points of adjacent n-2 th and n-1 th spectral bins;
step 2.3: an expression of an empirical scale function and a wavelet function expression of an empirical wavelet are respectively constructed using the expression (2) and the expression (3):
in the formula (2) and the formula (3), phi n (omega) represents the nth empirical scale function, ψ n (ω) represents the nth empirical wavelet function, γ represents a constant between 0 and 1, and β (·) is a value of [0,1 ]]A function of (a); omega n+1 Represents adjacent nth and nth +
1 boundary point of spectrum segment; omega is the angular frequency;
step 2.4: defining detail correlation coefficients by means of equations (4) and (5), respectivelyApproximate correlation coefficient W f ε (0,t):
In the formula (4), "-" represents a conjugate operation, and when n=1, φ 1 (t) is defined by phi n (ω) inverse fourier transforming to obtain a function value; psi phi type n (t) consists of psi n (omega) function obtained by inverse Fourier transform, y a (t) is defined by y a (ω) inverse fourier transforming the resulting function; f (F) -1 (. Cndot.) represents the inverse Fourier transform; τ is an integral variable;
step 2.5: by using (6) the analysis signal y a (t) reconstructing to obtain an analysis signal y a The empirical mode components of (t) are shown in the formulas (7) and (8):
in equation (6), the convolution operation,is->Fourier transform form of>When n=0 is expressedIs a value of (2);
in formula (7), y 0 (t) represents an amplitude modulation-frequency modulation single component when n=0,representation->Is an inverse fourier transform form of (a);
in formula (8), y n (t) represents amplitude-modulated-frequency components of different frequencies, n=1, 2,3, N;
step 3: the analysis signal y is compared with the reference signal y by the method (9) a Each modal component y of (t) n (t) performing smooth pseudo-Wigner-Ville transformation to obtain a smooth pseudo-Wigner-Ville time-frequency distribution expression
In the formula (9), the amino acid sequence of the compound,and->Representing y n Transient correlation function of (t),>representation->Conjugation of (2); g (s-t) is a time domain window function, h (τ) is a frequency domain window function, and g (0) =h (0) =1, s and τ are both integral variables; omega is the angular frequency;
step 4: using (10) to divide each mode component y n Linear superposition of the smoothed pseudo Wigner-Ville analysis results of (t) to obtain an analysis signal y a (t) smooth pseudo-Wigner-Ville distribution
Step 5: for the analytic signal y a Sliding pseudo-Wigner-Ville distribution of (t)And performing interference detection, and filtering the obtained interference signal by adopting a notch IIR filter, so as to obtain the sliding pseudo Wigner-Ville distribution after filtering.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210145552.6A CN114488208B (en) | 2022-02-17 | 2022-02-17 | Beidou signal anti-interference method combining empirical wavelet and SPWVD conversion |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210145552.6A CN114488208B (en) | 2022-02-17 | 2022-02-17 | Beidou signal anti-interference method combining empirical wavelet and SPWVD conversion |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114488208A CN114488208A (en) | 2022-05-13 |
CN114488208B true CN114488208B (en) | 2024-04-05 |
Family
ID=81482207
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210145552.6A Active CN114488208B (en) | 2022-02-17 | 2022-02-17 | Beidou signal anti-interference method combining empirical wavelet and SPWVD conversion |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114488208B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106597408A (en) * | 2016-12-16 | 2017-04-26 | 重庆邮电大学 | Method for estimating high-order PPS signal parameter based on time-frequency analysis and instantaneous frequency curve-fitting |
WO2017205382A1 (en) * | 2016-05-23 | 2017-11-30 | The University Of New Hampshire | Techniques for empirical mode decomposition (emd)-based signal de-noising using statistical properties of intrinsic mode functions (imfs) |
US9838151B1 (en) * | 2015-03-03 | 2017-12-05 | Rockwell Collins, Inc. | Wavelet transformation based anti-jam processing techniques |
CN107783200A (en) * | 2017-11-21 | 2018-03-09 | 吉林大学 | Joint EMD and TFPF algorithms a kind of all-wave magnetic resonance signal random noise method for reducing |
CN108362942A (en) * | 2018-02-11 | 2018-08-03 | 中国铁道科学研究院 | The time-frequency spectrum acquisition methods and device of multicomponent signal |
CN111767811A (en) * | 2020-06-18 | 2020-10-13 | 东南大学 | Non-stationary signal analysis method based on bandwidth Fourier decomposition |
-
2022
- 2022-02-17 CN CN202210145552.6A patent/CN114488208B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9838151B1 (en) * | 2015-03-03 | 2017-12-05 | Rockwell Collins, Inc. | Wavelet transformation based anti-jam processing techniques |
WO2017205382A1 (en) * | 2016-05-23 | 2017-11-30 | The University Of New Hampshire | Techniques for empirical mode decomposition (emd)-based signal de-noising using statistical properties of intrinsic mode functions (imfs) |
CN106597408A (en) * | 2016-12-16 | 2017-04-26 | 重庆邮电大学 | Method for estimating high-order PPS signal parameter based on time-frequency analysis and instantaneous frequency curve-fitting |
CN107783200A (en) * | 2017-11-21 | 2018-03-09 | 吉林大学 | Joint EMD and TFPF algorithms a kind of all-wave magnetic resonance signal random noise method for reducing |
CN108362942A (en) * | 2018-02-11 | 2018-08-03 | 中国铁道科学研究院 | The time-frequency spectrum acquisition methods and device of multicomponent signal |
CN111767811A (en) * | 2020-06-18 | 2020-10-13 | 东南大学 | Non-stationary signal analysis method based on bandwidth Fourier decomposition |
Non-Patent Citations (3)
Title |
---|
A Novel GNSS Interference Detection Method Based on Smoothed Pseudo-Wigner–Hough Transform;Kewen Sun 等;sensors;20210624;第21卷(第4306期);1-22 * |
地震信号线性与非线性时频分析方法对比;黄昱丞;郑晓东;栾奕;杨廷强;;石油地球物理勘探;20181015(第05期);101-115 * |
基于重排小波变换的 GNSS 接收机干扰检测;赵慧子 等;科技创新与应用;20171231(第11期);13-15 * |
Also Published As
Publication number | Publication date |
---|---|
CN114488208A (en) | 2022-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0932837B1 (en) | Procedure for the elimination of interference in a radar unit of the fmcw type | |
EP3816665B1 (en) | Interference suppression in a fmcw radar system | |
Ford et al. | RFI mitigation techniques in radio astronomy | |
Turley | Impulsive noise rejection in HF radar using a linear prediction technique | |
Ellingson et al. | Mitigation of radar interference in L-band radio astronomy | |
CN115549709B (en) | Satellite communication system and method for inhibiting multi-channel mutual interference | |
WO2018009476A1 (en) | System and method for mitigating narrow-band interference | |
CN112684419A (en) | Anti-intermittent sampling forwarding type interference processing method and system based on double LFM cancellation | |
CN116148802A (en) | Bistatic SAR phase synchronization anti-interference method | |
Lamont-Smith et al. | Filtering approaches for interference suppression in low‐frequency SAR | |
CN114488208B (en) | Beidou signal anti-interference method combining empirical wavelet and SPWVD conversion | |
US6731233B2 (en) | Method of suppressing jammer signals | |
CN101639530B (en) | SAR echo signal de-noising preprocessing method based on two-dimensional mixed transformation | |
CN113238194B (en) | Broadband phased array radar anti-decoy interference method based on fractional domain-frequency domain processing | |
Miller et al. | Near-least-squares radio frequency interference suppression | |
Ferrell | Interference suppression in UHF synthetic aperture radar | |
JP2022120352A (en) | Radar signal processing device, radar device, radar signal processing method, and radar signal processing program | |
Lu et al. | Impulsive noise excision and performance analysis | |
서보석 et al. | Performance of Interference Mitigation with Different Wavelets in Global Positioning Systems | |
CN114675307B (en) | Satellite navigation receiver interference detection method based on fractional order pseudo Wigner-Ville distribution | |
CN114184848B (en) | Goertzel algorithm-based satellite-borne VHF transient signal point-by-point scanning real-time processing method | |
CN117665810B (en) | Ionosphere electron density detection method, system and equipment for linear frequency modulation signals | |
Luo et al. | Adaptive interference suppression based on Lv distribution for DSSS communications | |
Bartlett et al. | Extraction, analysis and interpretation of digital ionograms | |
Qiu et al. | A system design of real-time narrowband RFI detection and mitigation for GNSS-R receiver |
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 |