CN115657086A - Satellite navigation terminal interference and anti-interference simulation verification platform - Google Patents

Satellite navigation terminal interference and anti-interference simulation verification platform Download PDF

Info

Publication number
CN115657086A
CN115657086A CN202211341909.4A CN202211341909A CN115657086A CN 115657086 A CN115657086 A CN 115657086A CN 202211341909 A CN202211341909 A CN 202211341909A CN 115657086 A CN115657086 A CN 115657086A
Authority
CN
China
Prior art keywords
interference
signal
satellite navigation
frequency
data stream
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211341909.4A
Other languages
Chinese (zh)
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.)
Chinese People's Liberation Army 63883
Original Assignee
Chinese People's Liberation Army 63883
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 Chinese People's Liberation Army 63883 filed Critical Chinese People's Liberation Army 63883
Priority to CN202211341909.4A priority Critical patent/CN115657086A/en
Publication of CN115657086A publication Critical patent/CN115657086A/en
Pending legal-status Critical Current

Links

Images

Abstract

A satellite navigation terminal interference and anti-interference simulation verification platform comprises a satellite navigation signal generation module, an interference signal generation module, a radio frequency front end simulation module, a signal capture process simulation module, a signal tracking process simulation module and an anti-interference processing simulation module, can generate satellite navigation signals and interference signal data streams with various frequency points, performs power spectrum analysis, can simulate the processing processes of band-pass filtering, down-conversion and analog-to-digital conversion of the radio frequency front end to the radio frequency signals, can visualize the signal capture and tracking demodulation processes, is internally provided with a space-domain, space-time and space-frequency self-adaptive zero-adaptive anti-interference processing algorithm, and can output processed array antenna directional diagrams, output signal power, signal-interference-noise ratio and other information.

Description

Satellite navigation terminal interference and anti-interference simulation verification platform
Technical Field
The invention relates to a satellite navigation anti-interference terminal and an array antenna anti-interference technology, in particular to a satellite navigation interference and anti-interference simulation verification platform.
Background
The satellite navigation terminal is used as the core and the foundation of a navigation, positioning and time service system, plays an extremely important role in the battle mission carried out by military weaponry, but the defects of weak signals and poor anti-interference capability are gradually shown, various intentional and unintentional interferences easily cause that a satellite navigation receiver cannot normally lock satellite signals, so that the performance of the satellite navigation receiver is reduced or the satellite navigation receiver cannot be used within a period of time, which undoubtedly forms a potential safety hazard for the application of the satellite navigation terminal, and the anti-interference performance of the satellite navigation terminal becomes one of the most important performance requirements of the receiver. Navigation wars have penetrated into each link of a war, and play an important role in each stage of the progress of the war. In order to improve the anti-interference capability of the weapon equipment, the vehicle-mounted, airborne, missile-borne and other weapon navigation terminals mostly adopt a multi-array element antenna self-adaptive zeroing technology, and the radio-frequency signals received by different array elements can be weighted to achieve the effect of inhibiting interference signals and improve the navigation and positioning capability of the weapon equipment in a complex electromagnetic environment.
The satellite navigation anti-interference terminal generally comprises a radio frequency front end module, an anti-interference processing module, a signal capturing module and a tracking demodulation module, wherein the radio frequency front end mainly completes the processes of filtering, down-conversion, automatic gain control, analog-to-digital conversion and the like of an antenna receiving signal, the anti-interference processing module utilizes adaptive filtering algorithms such as airspace, space-time or space-frequency to obtain a weight value to complete the suppression of an interference signal, and the signal capturing and tracking demodulation module completes the capturing of a digital baseband signal and the demodulation of a navigation message.
In the prior art, the effect of an interference signal on a satellite navigation terminal is researched mainly based on an outdoor experiment, although the method is a powerful means for verifying the interference effect, the method needs repeatability verification, is based on a large amount of experimental data, needs to comprehensively consider various test conditions such as an interference signal pattern and an interference direction, is easily influenced by factors such as weather and starry sky conditions, releases the interference signal for a long time to generate certain influence on peripheral navigation equipment, and has a short time window for research when a flight experiment for testing the anti-interference performance of missile-loading and airborne navigation equipment is carried out. On the other hand, outdoor testing is mainly carried out in a complete machine mode, the specific influence of interference signals on the radio frequency front-end module, the anti-interference processing module, the signal capturing module and the tracking demodulation module is difficult to qualitatively and quantitatively evaluate, the action mechanism of the interference signals cannot be clearly shown, and theoretical support is difficult to find.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a simulation analysis and verification platform capable of analyzing the action effect and the influence rule of interference signals on the radio frequency front end of the satellite navigation anti-interference terminal, the anti-interference processing, capturing and tracking demodulation process.
The above object of the present invention can be achieved by the following measures, a satellite navigation interference and anti-interference simulation verification platform, comprising:
the satellite navigation signal generation module: the Beidou satellite navigation signal data stream of the modulated custom data codes and the range finding codes can be generated, the satellite navigation signal data stream of the Beidou B1I, B3I, B1C, B2a, B2B and GPS L1C/A satellite navigation signals of the modulated custom data codes and the custom range finding codes can also be generated, the satellite navigation signal data stream of the Beidou B1Q, B2I/Q, B3Q, GPS P, GPS M and any BOC (M, n) system can also be generated, the code phase and the carrier Doppler frequency shift value of the satellite navigation signals can be set, and the power spectrum of the generated satellite navigation signal data stream can be displayed.
An interference signal generation module: the interference signal data stream of single tone, white noise, band-limited white noise, natural noise, noise amplitude modulation, noise frequency modulation, linear frequency sweep, BPSK, QPSK, 2FSK, 2ASK, MSK, BOC (m, n) may be generated, and the power spectrum of the interference signal data stream may be displayed.
The radio frequency front end simulation module: the band-pass filtering, frequency mixing, band-pass sampling and A/D conversion of the radio frequency signal data stream are realized, a digital intermediate frequency signal with specified frequency can be generated in a configurable manner, the power spectrum of a processing signal at each stage can be displayed, and the digital intermediate frequency signal data stream is output;
a signal capture process simulation module: acquisition algorithm processing can be performed on the digital intermediate frequency signal data stream to determine the PRN number, code phase and carrier doppler shift of the navigation satellite.
A signal tracking process simulation module: the acquisition result of the signal acquisition process simulation module can be analyzed, the frame bit synchronization, the de-spread and the demodulation of the digital intermediate frequency signal data stream are realized, and finally, the navigation message is output.
The anti-interference processing simulation module: the spatial domain, space-time and space-frequency adaptive filtering algorithm is built in, a four-array-element anti-interference antenna model can be simulated, different satellite navigation signals and interference signal incoming wave directions can be simulated, and an antenna directional diagram, an output signal-to-interference-and-noise ratio and output power after adaptive filtering can be displayed.
The satellite navigation interference and anti-interference simulation verification platform is operated and used based on an MATLAB environment.
Compared with the prior art, the invention has the following beneficial effects: the effect and the mechanism of action of the interference signal on the satellite navigation terminal can be analyzed in a simulation mode, the specific influence of the interference signal on the radio frequency front end of the terminal, the anti-interference processing, capturing and tracking demodulation process is shown, and a research platform is provided for navigation signal characteristic analysis, novel interference signal design, digital signal power spectrum analysis, interference countermeasure scheme research and anti-interference performance test.
Drawings
FIG. 1 is a flow chart of the operation of a satellite navigation signal simulation module according to the present invention;
FIG. 2 is a flow chart of the operation of the interference signal simulation module of the present invention;
FIG. 3 is a flow chart of the RF front end simulation module according to the present invention;
FIG. 4 is a flow chart of the signal acquisition process simulation module of the present invention;
FIG. 5 is a flow chart of the signal tracking process simulation module of the present invention;
FIG. 6 is a flowchart of the operation of the anti-interference processing simulation module of the present invention;
FIG. 7 is a flow chart of space-domain, space-time algorithm of the anti-interference processing simulation module of the present invention;
fig. 8 is a flow chart of the space-frequency algorithm of the anti-interference processing simulation module of the present invention.
Detailed Description
The invention is described in detail below with reference to the drawings and examples so that the advantages and features of the invention may be more readily understood by those skilled in the art, and thus the scope of the invention will be more clearly and clearly defined.
Referring to fig. 1, the satellite navigation signal generating module is implemented as follows:
firstly, selecting a satellite navigation frequency point and a satellite number, setting a carrier Doppler frequency shift value, a code phase value, a self-defined data code, signal power and data frequency parameter information, then carrying out carrier modulation according to a central frequency, a modulation mode, a ranging code constitution and a code rate which correspond to the navigation frequency point to generate a corresponding satellite navigation signal data stream, and finally carrying out spectrum analysis on the data stream.
(1) Satellite navigation frequency point and satellite number setting
The satellite navigation frequency points supporting selection comprise Beidou No. two B1I/Q, B2I/Q and B3I/Q, beidou No. three B1C, B2a and B2B and GPS L1C/A, P and M, and the frequency points supporting selection of satellite numbers comprise B1I, B3I, B1C, B2a, B2B and GPS L1C/A.
(2) Satellite navigation signal generation
(1) Beidou second satellite navigation signal generation
The Beidou I, B1I/Q, B2I/Q and B3I/Q frequency point navigation signals are modulated by QPSK, and the jth signal reaching a satellite navigation terminal receiving antenna can be expressed as follows:
Figure BDA0003914858830000041
wherein
Figure BDA0003914858830000042
For ranging code, according to Beidou space signalThe interface file is required to be acquired,
Figure BDA0003914858830000043
Figure BDA0003914858830000044
to define the data code, A BI 、A BQ Is the signal amplitude,
Figure BDA0003914858830000045
Is the initial phase, f, of the carrier wave 0 Is the center frequency, f d Is the carrier doppler shift.
(2) Beidou third satellite navigation signal generation
The design of the Beidou third frequency point signal comprises a pilot frequency component and a data component, wherein the data component modulates a ranging code and a navigation message, and the pilot frequency component modulates only the ranging code. The Beidou B1C data component is modulated by BOC (1, 1), the pilot frequency component is modulated by QMBOC (6, 1/11), and the data component and the pilot frequency component signal of the Beidou B1C frequency point reaching the receiving antenna at the jth point can be expressed as follows:
Figure BDA0003914858830000051
Figure BDA0003914858830000052
in the formula f B1C_a At 1.023MHz B1C_b Is at a frequency of 6.138MHz,
Figure BDA0003914858830000053
in order to be the amplitude of the signal,
Figure BDA0003914858830000054
for the data component of the ranging code,
Figure BDA0003914858830000055
for the pilot component ranging codes, the pilot component,
Figure BDA0003914858830000056
as initial phase of carrier wave, f j Is the superposition value of the carrier center frequency and the Doppler frequency shift.
The data component and the pilot frequency component of the Beidou B2a frequency point signal are modulated by BPSK (10), only the phases are different, and the data component and the pilot frequency component signal of the Beidou B2a frequency point signal of the jth navigation satellite reaching the receiving antenna can be expressed as follows:
Figure BDA0003914858830000057
in the formula
Figure BDA0003914858830000058
In order to be the amplitude of the signal,
Figure BDA0003914858830000059
for the data component of the ranging code,
Figure BDA00039148588300000510
for the pilot component of the ranging code to be,
Figure BDA00039148588300000511
in order to be the data code,
Figure BDA00039148588300000512
for the initial phase of the signal carrier, f j Is the superposition value of the carrier center frequency and the Doppler frequency shift.
The B2B signal also adopts a BPSK (10) modulation mode, only has an I branch, and is the same as the Beidou B2a data component expression.
(3) GPS satellite navigation signal generation
The C/A code of the L1 frequency point and the carrier signal modulated by the P code are orthogonal with each other, and the signal expression mode is as follows:
Figure BDA00039148588300000513
in the formula A C/A 、A P In order to be the amplitude of the signal,
Figure BDA0003914858830000061
in order to be a signal ranging code,
Figure BDA0003914858830000062
a data code on the ranging code and a ranging code,
Figure BDA0003914858830000063
for the initial phase of the signal carrier, f j Is the superposition value of the carrier center frequency and the Doppler frequency shift.
The GPS M code signal is a satellite navigation signal modulated by BOC (10, 5), and the signal is expressed as follows:
Figure BDA0003914858830000064
wherein f is Ma Is at a frequency of 10.23MHz and,
Figure BDA0003914858830000065
in order to be the amplitude of the signal,
Figure BDA0003914858830000066
for the data component of the ranging code to be,
Figure BDA0003914858830000067
in order to be the data code,
Figure BDA0003914858830000068
for the initial phase of the signal carrier, f j Is the superposition value of the carrier center frequency and the Doppler frequency shift.
(3) Spectral analysis
And setting the center frequency and the display bandwidth of the spectrum analysis, namely displaying the spectrum condition of the generated satellite navigation signal data stream through a program.
Referring to fig. 2, the implementation process of the interference signal generation module is described as follows: firstly, an interference signal pattern is selected, interference signal center frequency and bandwidth parameters are set, then, a corresponding interference signal data stream is generated according to a various interference signal generation mode, and finally, frequency spectrum analysis can be carried out on the data stream.
(1) Interference signal pattern selection
The examples support selected interference signal patterns including single tone, natural noise, band-limited white noise, noise amplitude modulation, noise frequency modulation, linear frequency sweep, BPSK, QPSK, 2FSK, 2ASK, MSK, BOC (m, n), pseudo code dependent interference.
(2) Interference signal parameter setting
The method comprises the steps of setting the center frequency and the bandwidth of an interference signal, setting parameters for a noise frequency modulation signal, further comprising a modulation slope, further comprising a frequency sweep rate for a linear frequency sweep signal, and further comprising a multiple of a subcarrier frequency relative to a basic code rate and a multiple of a pseudo code rate relative to the basic code rate for BOC modulation.
(3) Interference signal generation
(1) Single-tone interference:
i.e. interference signals with fixed signal frequency, the expression formula is:
J=Acos(2πf 0 t) (7)
wherein A is the carrier amplitude, f 0 Is the center frequency.
(2) Frequency sweep interference
I.e., an interference signal whose signal frequency varies with time, the expression formula is:
Figure BDA0003914858830000071
wherein A is the carrier amplitude, T is the sweep period, f 1 As the starting frequency, f 1 To terminate the frequency, f (t) is the instantaneous frequency of the signal during the sweep.
(3) Natural noise
I.e. a noise signal with a power density of-174 dB/Hz.
(4) Band limited white noise
Generated by passing a white gaussian noise signal through a band pass filter.
(5) Amplitude modulation of noise
The noise is taken as a modulation signal, the AM modulation is carried out on the carrier signal, so that the amplitude of the carrier signal is randomly changed along with the baseband noise, and the expression formula is as follows:
J=(A+U n (t))cos(2πf 0 t) (9)
wherein A is the carrier amplitude, U n And (t) is baseband noise, and the bandwidth of the noise amplitude modulation signal is 2 times of the bandwidth of the baseband noise.
(6) Noise frequency modulation
Taking noise as a modulation signal, carrying out FM modulation on a carrier signal to ensure that the frequency of the carrier signal randomly changes along with baseband noise, and the expression formula is as follows:
Figure BDA0003914858830000072
wherein A is the carrier amplitude, U n (t) is the baseband noise, K FM In order to modulate the slope of the light,
Figure BDA0003914858830000073
is the carrier signal phase.
(7) BPSK modulation (binary phase shift keying)
A set of 1-1 symbol sequences is randomly generated at a code rate fv and is continued with a time period, and then a carrier is modulated according to a phase relation of 0 and pi, wherein the signal bandwidth is 2fv.
(8) QPSK modulation (Quadrature phase Shift keying)
Two groups of code element sequences of {1, -1} are randomly generated according to a code rate fv, and are continued along with a time period, and then carriers are modulated according to a phase relation of { pi/4, 3 pi/4, 5 pi/4 and 7 pi/4 } and the signal bandwidth is 2fv.
(9) 2FSK modulation (binary frequency Shift keying)
A set of 1-1 symbol sequences is randomly generated at a code rate fv and is spread over a time period, and then modulated at different carrier frequencies.
(10) 2ASK modulation (binary amplitude keying)
A set of 1-1 symbol sequences is randomly generated at a code rate fv, and is spread over a period of time, and then modulated at different amplitudes.
(11) MSK modulation (minimum phase shift keying)
First according to code rate f v Randomly generating a group of {1, -1} code element sequences, extending with time period, and modulating as follows, the bandwidth of the modulation signal is 1.5f v
Figure BDA0003914858830000081
In the formula a n Is the symbol value.
(12) BOC (m, n) modulation (binary offset carrier)
Firstly, randomly generating a group of {1, -1} code element sequences according to the code rate nfv, continuing along with a time period, simultaneously generating mfv subcarrier signals with frequency along with time, carrying out exclusive-or multiplication on the code element time sequences and subcarriers, and carrying out BPSK modulation to obtain BOC modulation signals.
(13) Pseudo code dependent interference
Namely, the generated {1, -1} code element time sequence has certain correlation with the satellite navigation signal, continues along with the time period, and then is modulated according to the modulation method corresponding to the navigation frequency point.
(4) Spectral analysis
And setting the center frequency and the display bandwidth of the spectrum analysis, namely displaying the spectrum condition of the generated satellite navigation signal data stream through a program.
Referring to fig. 3, the implementation process of the rf front-end simulation module is described as follows: firstly, the front-end filtering bandwidth, the analog intermediate frequency value, the A/D sampling frequency and the AGC (automatic gain control) output power are set, then the superposed satellite navigation signal data stream and the superposed interference signal data stream can be subjected to radio frequency filtering, LNA (low noise amplifier), image filtering, frequency mixing, intermediate frequency filtering, AGC and ADC (analog-to-digital converter), and finally the output of the digital intermediate data stream is realized.
(1) A radio frequency filter: filtering out-of-band interference, extracting useful signals in a working frequency band and improving the signal-to-noise ratio.
(2) LNA: and amplifying the signal processed by the radio frequency filter.
(3) An image filter: the image frequency components are filtered out to prevent the image signal from being fed into the mixer.
(4) A mixer: and mixing the image-filtered signal with a local oscillation generation signal is realized.
(5) Local carrier generator: a carrier signal of a specified frequency is provided to the mixer.
(6) An intermediate frequency filter: and filtering high-frequency signal components in the mixing signals to obtain analog intermediate-frequency signals.
(7) AGC: the power of the analog intermediate frequency signal is automatically adjusted to tend to a certain constant value, so that ADC sampling is facilitated.
(8) ADC: and sampling the analog intermediate frequency signal after AGC to obtain digital intermediate frequency data.
Referring to fig. 4, the implementation of the simulation module for the signal acquisition process is described as follows: the process performs acquisition algorithm processing on the digital intermediate frequency signal data stream to determine the PRN number, code phase and carrier doppler shift of the navigation satellite.
(1) Performing a Capture Algorithm
The known jth beidou B3I signal arriving at the receiving antenna can be expressed as follows:
Figure BDA0003914858830000101
setting the working frequency of the local carrier generator to be
Figure BDA0003914858830000102
Figure BDA0003914858830000103
Is an estimate of the carrier doppler shift. In the figure, INCO is a carrier cosine signal, QNCO is a carrier sine signal, an intermediate frequency signal is multiplied by INCO and QNCO generated by a local carrier generator respectively, and the multiplied intermediate frequency signal passes through a low-pass filter to obtain a signal:
Figure BDA0003914858830000104
the local code signal is processed by FFT, conjugate, and multiplied by the result after FFT, and then IFFT is processed to obtain the following result:
Figure BDA0003914858830000105
above formula T s The coherent integration duration of single acquisition is generally an integral multiple of the ranging code period, and is obtained by further performing modulus on Y:
Figure BDA0003914858830000106
(2) peak determination
Outputting Y value distribution chart and obtaining Y max . If Y- max If the signal is greater than the acquisition threshold, the satellite navigation signal is acquired.
(3) Outputting the captured result
After successful capture, | Y- max The corresponding navigation signal number, code phase and carrier Doppler frequency shift are the acquisition results.
Referring to fig. 5, the implementation of the simulation module for the signal tracking process is described as follows: the process realizes the frame bit synchronization, the de-spread and the demodulation of the digital intermediate frequency signal data flow according to the capture result, and finally outputs the navigation message.
(1) Frame bit synchronization
The bit synchronization is used for locking each data bit boundary of the navigation message, and the frame synchronization is used for locking the boundary between navigation message frames so as to determine the initial position of the frames;
(2) despreading (removing ranging codes)
If the dynamically adjusted local code phase is τ, the signal despreading process is as follows:
Figure BDA0003914858830000111
C τ (t) is the generated local ranging code signal with phase τ, since the local code remains aligned with the received code, i.e., C (t) C τ (t) =1, so there are:
Figure BDA0003914858830000112
(3) demodulation (removing carrier)
The tracking loop makes the frequency and phase of the local carrier wave and the receiving carrier wave synchronous by dynamically adjusting the frequency difference and phase difference of the local carrier wave and the receiving carrier wave to obtain the Doppler frequency shift of the receiving carrier wave
Figure BDA0003914858830000113
And
Figure BDA0003914858830000114
the signal demodulation process is as follows:
the satellite navigation receiver generates a local carrier:
Figure BDA0003914858830000115
then with s 1 (t) are multiplied and passed through a low pass filter to obtain:
Figure BDA0003914858830000116
if it is
Figure BDA0003914858830000117
I branch output is
Figure BDA0003914858830000118
The Q branch output is 0. Because A is amplitude and is always greater than 0, and D (t) value is +/-1, namely the I branch output is B3I data code, after integration, the I branch can be obviously identified, the local ranging code of B3Q is generated in the same way, the B3Q data code can be demodulated by the above process, and the Q branch output。
Referring to fig. 6, the implementation process of the anti-interference processing simulation module is described as follows: firstly, setting an interference signal data stream incoming wave direction, a satellite navigation signal data stream incoming wave direction and a self-adaptive filtering algorithm, then generating corresponding array element guide vectors, calculating self-adaptive weights by a self-adaptive filter, and finally outputting a power directional diagram, a signal-to-interference-and-noise ratio and output power.
(1) Incoming wave direction and adaptive algorithm setting
The incoming wave directions are the incidence azimuth angle and the pitch angle of the signal, and the self-adaptive algorithm comprises an airspace algorithm, a space-time algorithm and a space-frequency algorithm.
(2) Steering vector generation
The design of the satellite navigation anti-interference antenna mostly adopts a uniform circular array, and if M array elements are uniformly distributed on a circle with the radius of R. Wherein θ and
Figure BDA0003914858830000121
the pitch angle and azimuth angle of the incident signal S, and gamma m is the included angle of the array element relative to the circle center. The coordinates of the mth array element can be expressed as:
Figure BDA0003914858830000122
the incident signal vector is represented as:
r=(sinθcosφ,sinθsinφ,cosθ) (19)
then the phase lag of the mth array relative to the center is:
Figure BDA0003914858830000123
from this, the steering vector of the incident signal can be derived:
Figure BDA0003914858830000124
(3) adaptive filtering
(1) Spatial and spatial filtering
The filter structure is shown in fig. 7, assuming that the antenna array has M array elements, and L time-domain taps are added behind each array element, a space-time filter with L taps can be obtained (L =1 is pure spatial filtering)
The resulting input signal of the filter can be expressed as:
Figure BDA0003914858830000125
the weights corresponding to the input signals are set as:
w=[w 11 ,…,w 1L ,…,w M1 ,…,w ML ] T (23)
the final output of the space-time filtering can be obtained as:
y(n)=w H x(n) (24)
the way of solving the weight value according to the SMI algorithm is as follows:
Figure BDA0003914858830000131
rxx is the autocorrelation matrix of the input signal, and if the limited sampling data of N snapshots is x (N) (N =1,2, \8230;, N), the autocorrelation matrix of the input signal is:
Figure BDA0003914858830000132
(2) Space frequency filtering
The structure of the filter is shown in fig. 8, and the space-frequency adaptive filtering is to perform anti-interference processing on a radio frequency signal in a joint space domain and a frequency domain, and the principle is to divide a broadband signal into a plurality of frequency spectrum components after an input signal is subjected to FFT, and then perform space-domain filtering on each frequency spectrum component. Supposing that the antenna array is composed of M array elements, performing J-point FFT on J data received by each array element, and obtaining input signals of J spectrum components after processing, and supposing that the jth spectrum value of the mth array element is Xm (fj), the input data for performing spatial filtering anti-interference processing on the jth spectrum component is:
X(f j )=[X 1 (f j ),X 1 (f j ),…,X M (f j )] T (27)
because J spectral components of an array element input signal are required to be obtained, FFT transformation operation is performed on J data received by the array element, and an input signal received by an array element from 1 st to J th time is called a data block, that is, an nth data block of an mth array element is:
X m (n)=[x m (n),x m (n+1),…,x M (n+J-1)] T (28)
the value of the jth spectral component obtained by performing FFT on the mth array element by the nth data block is Xm (fj, n), and the input data of the jth spectral component at this time may be represented as:
X(f j ,n)=[X 1 (f j ,n),X 1 (f j ,n),…,X M (f j ,n)] T (29)
the space-frequency SMI algorithm is to perform space-domain SMI filtering anti-interference processing on each frequency spectrum component, and to input signals X (fj) of the jth frequency spectrum component, limited sampling data X (f) of N frequency-domain snapshots are also passed through j, N), N =1,2, \ 8230;, N, where a spatial SMI filter input signal correlation matrix is obtained, denoted as:
Figure BDA0003914858830000141
after the autocorrelation matrix of the input signal for each spectral component is found, the weights w (fj) for the corresponding spectral components can be found according to the following equation.
Figure BDA0003914858830000142
Therefore, the space-time filtering algorithm only needs to invert the M × M dimensional matrix, and then according to the following formula:
Y(f j ,n)=w(f j ) H X(f j ,n) (32)
obtaining the anti-interference output value Y (fJ, n) of J frequency spectrum components, and finally obtaining the anti-interference values Y (f 1, n), Y (f 2, n), \8230, and the anti-interference time domain data Y (n), Y (n + 1), \8230, and Y (n + J-1) of the J frequency spectrum components through IFFT transformation.
(4) Outputting anti-interference processing result
After the signal is processed by the adaptive filter, the adaptive weight is output, and a power directional diagram, a signal-to-interference-and-noise ratio and output power can be further calculated.
The above description is only an embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes performed by the present specification and drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.

Claims (5)

1. A satellite navigation terminal interference and anti-interference simulation verification platform is characterized in that: the method comprises the following steps:
the satellite navigation signal generation module is used for generating a satellite navigation signal data stream, setting a code phase and a carrier Doppler frequency shift value of a satellite navigation signal and displaying a power spectrum of the generated satellite navigation signal data stream;
the interference signal generating module is used for generating an interference signal data stream and displaying a power spectrum of the interference signal data stream;
the radio frequency front end analog module is used for performing band-pass filtering, frequency mixing, band-pass sampling and A/D conversion on the satellite navigation signal data stream and the interference signal data stream, displaying the power spectrum of the processing signal at each stage and outputting a digital intermediate frequency signal data stream;
the signal acquisition process simulation module is used for carrying out acquisition algorithm processing on the digital intermediate frequency signal data stream so as to determine a PRN (pseudo random number), a code phase and a carrier Doppler frequency shift of a navigation satellite;
the signal tracking process simulation module is used for analyzing the capturing result of the signal capturing process simulation module, realizing frame bit synchronization, de-spreading and demodulation of digital intermediate frequency signal data streams and finally outputting navigation messages;
the anti-interference processing simulation module is internally provided with an airspace, space-time and space-frequency self-adaptive filtering algorithm and is used for simulating a four-array-element anti-interference antenna model, simulating different satellite navigation signals and interference signal incoming wave directions, and displaying an antenna directional diagram, an output signal-to-interference-plus-noise ratio and output power after self-adaptive filtering.
2. The satellite navigation terminal interference and interference rejection simulation verification platform according to claim 1, wherein: the implementation steps of the satellite navigation signal generation module are as follows: selecting a satellite navigation frequency point and a satellite number;
setting carrier Doppler frequency shift values, code phase values, user-defined data codes, signal power and data frequency parameter information;
then, carrying out carrier modulation according to the center frequency, the modulation mode, the ranging code constitution and the code rate corresponding to the navigation frequency point;
generating a corresponding satellite navigation signal data stream;
the data stream may be spectrally analyzed.
3. The satellite navigation terminal interference and interference rejection simulation verification platform according to claim 1, wherein: the implementation steps of the interference signal generation module are as follows:
selecting an interference signal pattern;
setting interference signal center frequency and bandwidth parameters;
generating a corresponding interference signal data stream according to the various interference signal generation mode;
a spectral analysis is performed on the data stream.
4. The satellite navigation terminal interference and interference rejection simulation verification platform according to claim 1, wherein: the implementation steps of the radio frequency front end simulation module are as follows:
setting a front-end filtering bandwidth, an analog intermediate frequency value, an A/D sampling frequency and AGC output power; performing radio frequency filtering, LNA, image filtering, frequency mixing, intermediate frequency filtering, AGC and ADC on the superposed satellite navigation signal data stream and interference signal data stream;
and finally, outputting the data stream in the digital.
5. The satellite navigation terminal interference and interference rejection simulation verification platform according to claim 1, wherein: the anti-interference processing simulation module comprises the following implementation steps: setting the incoming wave direction of an interference signal data stream, the incoming wave direction of a satellite navigation signal data stream and a self-adaptive filtering algorithm;
generating corresponding array element guide vectors, and calculating self-adaptive weight values by a self-adaptive filter;
an output power pattern, a signal to interference and noise ratio, and an output power.
CN202211341909.4A 2022-10-28 2022-10-28 Satellite navigation terminal interference and anti-interference simulation verification platform Pending CN115657086A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211341909.4A CN115657086A (en) 2022-10-28 2022-10-28 Satellite navigation terminal interference and anti-interference simulation verification platform

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211341909.4A CN115657086A (en) 2022-10-28 2022-10-28 Satellite navigation terminal interference and anti-interference simulation verification platform

Publications (1)

Publication Number Publication Date
CN115657086A true CN115657086A (en) 2023-01-31

Family

ID=84993104

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211341909.4A Pending CN115657086A (en) 2022-10-28 2022-10-28 Satellite navigation terminal interference and anti-interference simulation verification platform

Country Status (1)

Country Link
CN (1) CN115657086A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117452446A (en) * 2023-12-26 2024-01-26 中国人民解放军国防科技大学 Anti-interference processing architecture based on satellite navigation signal bi-component fusion application
CN117538905A (en) * 2023-11-22 2024-02-09 湖南跨线桥航天科技有限公司 Satellite navigation system with strong anti-interference capability and design method of terminal
CN117590430A (en) * 2023-11-27 2024-02-23 湖南跨线桥航天科技有限公司 M code signal smart interference method for GPS captured by punching code

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117538905A (en) * 2023-11-22 2024-02-09 湖南跨线桥航天科技有限公司 Satellite navigation system with strong anti-interference capability and design method of terminal
CN117590430A (en) * 2023-11-27 2024-02-23 湖南跨线桥航天科技有限公司 M code signal smart interference method for GPS captured by punching code
CN117590430B (en) * 2023-11-27 2024-04-19 湖南跨线桥航天科技有限公司 M code signal smart interference method for GPS captured by punching code
CN117452446A (en) * 2023-12-26 2024-01-26 中国人民解放军国防科技大学 Anti-interference processing architecture based on satellite navigation signal bi-component fusion application
CN117452446B (en) * 2023-12-26 2024-03-12 中国人民解放军国防科技大学 Anti-interference processing architecture based on satellite navigation signal bi-component fusion application

Similar Documents

Publication Publication Date Title
CN115657086A (en) Satellite navigation terminal interference and anti-interference simulation verification platform
US8688759B2 (en) Efficient detection algorithm system for a broad class of signals using higher-order statistics in time as well as frequency domains
US7885688B2 (en) Methods and systems for signal selection
Benedetto et al. A fast unambiguous acquisition algorithm for BOC-modulated signals
US9057772B1 (en) Referenced-based radio frequency feature estimation
TW201003110A (en) Method and apparatus for mitigating the effects of narrowband interfering signals in a GPS receiver
CN104158777B (en) United orthogonal projects the rangefinder pulse interference suppression method with circulation beam forming
CN102055512A (en) Circulating self-correlation-based signal phase difference estimation device and method for antenna array
CN106646529A (en) GNSS antenna array anti-interference method based on multi-beam optimization
CN102707270B (en) Automatic estimation method for antenna pattern relative to high-frequency ground wave radar
Nishio et al. A high-speed adaptive antenna array with simultaneous multibeam-forming capability
Elghamrawy et al. High-resolution spectral estimation for continuous wave jamming mitigation of GNSS signals in autonomous vehicles
CN101515912B (en) Detecting method for various interference signals of DS-CDMA system based on converting identifying schemes
KR102295191B1 (en) HILS(Hardware In Loop Simulation) system for Radar deception apparatus and simulation method thereof
Broumandan et al. Demonstration of a multi-layer spoofing detection implemented in a high precision gnss receiver
CN106772277A (en) BDS direct wave interference suppression methods based on signal reconstruction
Hornbostel et al. Simulation of multi-element antenna systems for navigation applications
US11670848B2 (en) Anti-jamming system
Lombardo et al. PBR activity at INFOCOM: adaptive processing techniques and experimental results
CN114839652A (en) Robust deception jamming suppression method under array error
Shen et al. Multiple-access interference and multipath influence mitigation for multicarrier code-division multiple-access signals
CN110441798B (en) Beidou RDSS weak signal capturing method based on multiplication accumulation integration and satellite selection assistance
CN113938221A (en) Method for testing anti-interference capability of blind source separation module
GB2418536A (en) Calibration of antenna array elements
US8743854B1 (en) Signal separation and SINR enhancement

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