Intelligent sensing decision-based binary array space-time-frequency code communication anti-interference method

Abstract

The application provides a binary array space-time-frequency code communication anti-interference method based on intelligent sensing decision, adopts a spectrum sensing technology to listen to and distinguish enemy interference signals, obtains relevant attribute parameters, and intelligently decides a communication mode of a self side according to different interference signal differentiation characteristics, and adopts a binary array space-time-frequency code joint anti-interference technology to improve communication anti-interference capability in various interference scenes including a full-band ultra-wideband interference. The application improves survival capability of a communication system in a complex electromagnetic interference environment. A control strategy of interference sensing and intelligent decision is introduced, a space-time-frequency-code joint anti-interference method is used according to interference signal differentiation characteristics, the system is more intelligent and robust in signal receiving and processing, and reliability of the system is enhanced.

Description

Technical Field

[0001] This invention relates to the field of frequency hopping communication anti-interference, and specifically to a communication anti-interference method combining space-time-frequency codes. Background Technology

[0002] In current data link communication technologies, the basic anti-jamming method relies on high-speed frequency hopping to evade enemy communication interference and information interception. Due to the complexity and secrecy of frequency hopping patterns, it is difficult for the enemy to directly crack them. However, once the enemy obtains the specific range of one's own communication frequency band through signal detection and other means, they can use ultra-wideband jamming signals covering the entire communication band to suppress interference with high power, causing one's own communication system to malfunction. At the same time, the current SJL communication system lacks the ability to detect, sense, and analyze the characteristics of jamming signals, and cannot intelligently adjust its own communication parameters to target the weaknesses of the enemy's jamming methods, thus limiting the improvement of the communication system's anti-jamming capability.

[0003] Based on the existing problems of full-band ultra-wideband interference and lack of intelligent sensing and decision-making, this invention provides a communication anti-interference method based on the joint use of binary array space-time-frequency codes with intelligent sensing and decision-making. Summary of the Invention

[0004] To overcome the shortcomings of existing technologies, this invention provides a binary array space-time-frequency code communication anti-interference method based on intelligent perception and decision-making. It employs spectrum sensing technology to detect and identify enemy interference signals, obtain their relevant attribute parameters, and intelligently decides on the user's communication mode based on the differentiated characteristics of different interference signals. Furthermore, it utilizes a binary array space-time-frequency code joint anti-interference technology to improve communication anti-interference capabilities under various interference scenarios, including full-band ultra-wideband interference.

[0005] The technical solution adopted by this invention to solve its technical problem includes the following steps:

[0006] Step 1: Divide the received signal into channels, such as... Figure 1 As shown, let the instantaneous bandwidth of the entire system be BW, and the number of frequency bands be K = (f1 - f0) / BW, where f1 represents the termination frequency and f0 represents the starting frequency. The entire frequency band is divided into K sub-bands, and the center frequency of each sub-band is w. k = k*BW, where k takes the value 0, 1, 2, ..., K;

[0007] Step 2: Perform intermediate frequency (IF) to baseband signal conversion on the divided channel. According to the Nyquist sampling theorem, time-domain sampling of the signal is equivalent to periodically expanding the signal's frequency domain by N*Fs, such as... Figure 2 As shown, after spectrum shifting, the effective frequency in the range (0, Fs / 2) is Fc, and the bandwidth is BW;

[0008] Step 3: After receiving the signal, the communication terminal uses orthogonal digital down-conversion to convert the intermediate frequency signal to baseband, such as... Figure 3 As shown, after the data is acquired by the ADC, it first undergoes digital down-conversion. The DDS generates a pair of orthogonal carrier signals COS and SIN with a frequency of Fc. Then, the ADC data is mixed to mix the Fc intermediate frequency signal to the baseband. Since a harmonic signal is generated during the mixing process, a low-pass filter (LPF) is used to filter out the harmonic signal. Finally, before the spectrum calculation, since the bandwidth of the interference signal to the main signal is unknown, a second frequency conversion is required. The signal rate is decimated using an integral-comb cascaded filter (CIC filter) to further reduce the signal sampling rate.

[0009] Step 4: Perform spectrum scanning according to the channel division, and perform multiple scans with a fixed bandwidth to obtain the signal spectrum information of the entire operating frequency band;

[0010] Step 5: Perform spectrum analysis based on the spectrum scan results to obtain the number of interference signal frequencies N, the power of the interference signal Pg, and the period of the interference signal Tg;

[0011] Step 6: Make intelligent decisions based on spectrum analysis and adjust the communication system mode. Let the actual number of operating frequencies be fn, the interference signal frequency fg ∈ the communication signal frequency fs, and the communication signal power be Ps. The intelligent decision-making process is as follows:

[0012] If the interference signal is a sudden signal:

[0013] (1) When N < 30% * fn, a waveform with a code rate less than 2 / 3 is used for communication. At this time, the communication system can work normally without loss of key performance.

[0014] (2) When 30%*fn≤N<60%*fn, a waveform with a code rate of less than 1 / 3 is used for communication. At this time, the communication system can work normally and the key performance is not lost.

[0015] (3) When 60%*fn≤N<90%*fn, a waveform with a code rate less than 1 / 10 is used for communication. At this time, the communication system can work normally and the key performance is not lost.

[0016] (4) When 60%*fn≤N<90%*fn, the spatial anti-interference method is used to filter out interference;

[0017] (5) When N≥100%*fn, the spatial domain filtering anti-interference method is used to filter out interference;

[0018] (6) When the interference signal power \(P_g > P_s\), and \(10\times\lg(P_g / P_s)<15\ dB\), \(f_n < f_s\), change the frequency hopping pattern to avoid interference or use spread - spectrum communication for data transmission;

[0019] (7) When the interference signal power \(P_g > P_s\), and \(10\times\lg(P_g / P_s)>15\ dB\), use the spatial domain filtering anti - interference method to filter out interference;

[0020] If the interference signal is a periodic interference signal, assume the interference signal period is \(T_g\), and take the communication signal with a 50% duty - cycle signal as an example, the period is denoted as \(T_s\);

[0021] (1) When \(T_g < 1 / 2\times T_s\), use burst communication to avoid the periodic interference signal by adjusting the start and end times of the signal;

[0022] (2) When \(1 / 2\times T_s\leq T_g\leq T_s\), communicate by increasing the code rate, compress the message communication time, and reduce the signal exposure time in space.

[0023] The spatial domain filtering method in step 6 adopts the spatial anti - interference method of a binary array, and the steps are as follows:

[0024] Step 1: For the frequency - hopping communication signals of two channels, use the local oscillator frequency \(f\) with a rapidly variable frequency response nco , and down - mix the narrow - band frequency - hopping pulse signal with the current center frequency \(f_0\) and bandwidth \(f\) band to the analog intermediate frequency according to the frequency - hopping pattern, and obtain two analog intermediate - frequency signals with a frequency of \(f\) IF ;

[0025] The center frequency \(f_0\) of the narrow - band frequency - hopping pulse is a discrete point, and the value range is The minimum frequency interval is \(f\) band , and the maximum frequency interval is \(band-(num - 1)\cdot f\) band ; where \(f\) L , \(f\) R , and \(band\) are the left boundary, right boundary of the actual working frequency band of \(f_{qTX}\) and the frequency band width respectively, \(num\) is the number of effective working frequency bands within the working frequency range, Here, \(fix()\) represents rounding down; the frequency \(f\) of the rapidly variable local oscillator nco =\(f_0 - f\) IF ;

[0026] Step 2: Sample the two analog intermediate - frequency signals using an AD with a sampling rate of \(f\) sample , \(f\) sampleThe signal frequency band bandwidth is no less than twice that of the signal band to obtain a digital intermediate frequency (IF) signal; a digital frequency synthesizer (DDS) is used to generate a frequency of f using a frequency control word. IF The sine and cosine signals are used to perform frequency mixing on the digital intermediate frequency signal, transforming the intermediate frequency signal into a baseband complex signal;

[0027] The frequency control word of the DDS has a bit width of n bits, where 8 ≤ n ≤ 64, and the value is calculated as follows: Among them, round() means rounding to the nearest integer, and dec2hex() means converting decimal to hexadecimal;

[0028] Step 3: Perform CIC filtering and R-fold speed reduction decimation on the baseband complex signal to obtain a low-speed baseband complex signal, and use an FIR low-pass filter to filter out useless out-of-band signals. At this time, two low-speed, spectrally clean baseband complex signals are obtained, denoted by din(1) and din(2) respectively.

[0029] The CIC filter design, based on a trade-off between performance and resources, employs a 4-stage cascaded design. The position of the first zero point on the frequency response curve of the CIC filter is chosen to be no greater than [missing value]. The maximum value within the range, the range of the deceleration factor R is: The passband flatness of the FIR filter should not exceed 1dB, and the out-of-band rejection should not be less than 40dB.

[0030] Step 4: Using the processed baseband signals din(1) and din(2) as input signals for spatial anti-interference processing, perform N-order data taps on each input signal, namely dex1(i) and dex2(i), i = 1, 2, ..., N, to obtain a total of 2N output signals, such as... Figure 5 As shown.

[0031] Step 5: In the absence of communication signal transmission and reception, set the local oscillator's operating frequency to the first operating frequency point;

[0032] Step 6: Combine the 2N tapped data points into a column vector x. Calculate the statistical average of the instantaneous covariance matrix of vector x at point p to obtain the statistical mean of the covariance matrix. Where x = [dex1(1), dex1(2), ..., dex1(N), dex2(1), dex2(2), ..., dex2(N)] T , here,[] T Represents the transpose of a vector;

[0033] Step 7: Using the formula Solve for w; where s∈C 2N×1The first element of s is 1, and the rest are 0. for The inverse matrix, w is the spatial anti-interference weight, w∈C 2N×1 ;

[0034] Step 8: Adjust to the working frequency, repeat steps 6 and 7 at each working frequency, traverse all working frequencies, and obtain the spatial anti-interference weights of all working frequencies;

[0035] Step 9: Under normal communication signal transmission and reception conditions, perform a weighted operation on the data tap vector x of the current working frequency and the anti-interference weight w corresponding to the current working frequency, thereby obtaining the anti-interference output of the current frequency. Thus, steps 5-9 complete the spatial filtering anti-interference operation.

[0036] An electronic device includes one or more processors; a memory; and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs are configured to perform the methods described above.

[0037] A computer-readable storage medium storing program code that can be invoked by a processor to perform the method described above.

[0038] The beneficial effects of this invention are as follows:

[0039] (1) This invention achieves rapid suppression of large bandwidth interference in high-speed frequency hopping communication, solves the problem that conventional communication systems cannot resist high-power suppression interference across the entire frequency band, and improves the survivability of communication systems in complex electromagnetic interference environments.

[0040] (2) This invention introduces a control strategy of interference perception and intelligent decision-making. Based on the differentiated characteristics of interference signals, it adopts a targeted anti-interference method that combines space, time, frequency and code, making the system more intelligent and robust in signal reception and processing, and enhancing the reliability of the system. Attached Figure Description

[0041] Figure 1 This is a channel partitioning diagram of the real signal in this invention.

[0042] Figure 2 This is a schematic diagram of the spectrum shifting process of the present invention.

[0043] Figure 3 This is a diagram of the orthogonal downconversion process of the present invention.

[0044] Figure 4 This is the spectrum clipping diagram of the present invention.

[0045] Figure 5 This is a diagram of the output signal of the present invention. Detailed Implementation

[0046] The present invention will be further described below with reference to the accompanying drawings and embodiments.

[0047] The technical solution adopted in this invention includes the following steps:

[0048] Step 1: Divide the received signal into channels, such as... Figure 1 Let the instantaneous bandwidth of the entire system be BW, and the number of frequency bands be K = (f1 - f0) / BW, where f1 represents the termination frequency and f0 represents the starting frequency. The entire frequency band is divided into k sub-bands, and the center frequency of each sub-band is w. k = k*BW, where k takes the value (0, 1, 2, ..., K).

[0049] Step 2: Perform intermediate frequency (IF) to baseband signal conversion on the divided channel. According to the Nyquist sampling theorem, time-domain sampling of the signal is equivalent to periodically expanding the signal's frequency domain by N*Fs. Figure 2 As shown, after spectrum shifting, the effective frequency in the range (0, Fs / 2) is Fc, and the bandwidth is BW.

[0050] Step 3: After receiving the signal, the communication terminal uses orthogonal digital down-conversion to convert the intermediate frequency signal to baseband, such as... Figure 3 As shown, after data is acquired by the ADC, it first undergoes digital down-conversion. The DDS generates a pair of orthogonal carrier signals, COS and SIN, with a frequency of Fc. Then, the ADC data is mixed to transfer the Fc intermediate frequency signal to the baseband. Since harmonic signals are generated during the mixing process, they are filtered out by a low-pass filter. It is important to note that before spectrum calculation, a second frequency conversion is required because the bandwidth of interference signals to the main signal is unknown. CIC is used to perform signal rate decimation, further reducing the signal sampling rate.

[0051] Step 4: Perform a spectrum scan based on the channel allocation, such as... Figure 4 As shown, multiple scans are performed each time with a fixed signal bandwidth. In order to complete the splicing of the entire frequency band spectrum, the actual calculated spectrum needs to be cropped to avoid spectral aliasing. Since the actual signal bandwidth is only BW during the FFT calculation process, and the signal sampling rate is Fs, according to the spectrum shifting law, when the signal is a real signal, the spectrum of the signal on the left half axis of the baseband (0 frequency) is exactly the same as the spectrum shifted to the left side of the signal at the (Fs) sampling rate, as shown in the shaded part of the figure. The two spectra are completely equivalent.

[0052] Therefore, the signal spectrum information of the entire operating frequency band can be obtained by using spectrum splicing.

[0053] Step 5: Perform spectrum analysis based on the spectrum scanning results to obtain the number N of interference signal frequencies, the power magnitude Pg of the interference signal, and the period Tg of the interference signal.

[0054] Step 6: Make an intelligent decision based on the spectrum analysis and adjust the communication system mode. Let the actual number of working frequencies be fn, the interference signal frequency fg ∈ the communication signal frequency fs, and the communication signal power be Ps. The decision-making process is as follows:

[0055] If the interference signal is a burst signal:

[0056] (1) When N < 30% * fn, use a waveform with a code rate less than 2 / 3 for communication. At this time, the communication system can work normally without loss of key performance.

[0057] (2) When 30% * fn ≤ N < 60% * fn, use a waveform with a code rate less than 1 / 3 for communication. At this time, the communication system can work normally without loss of key performance.

[0058] (3) When 60% * fn ≤ N < 90% * fn, use a waveform with a code rate less than 1 / 10 for communication. At this time, the communication system can work normally without loss of key performance.

[0059] (4) When 60% * fn ≤ N < 90% * fn, use a spatial domain anti-jamming method to filter out interference.

[0060] (5) When N ≥ 100% * fn, use a spatial domain filtering anti-jamming method to filter out interference.

[0061] (6) When the interference signal power Pg > Ps and 10 * lg(Pg / Ps) < 15 dB, and fn < fs, change the frequency hopping pattern to avoid interference or use spread spectrum communication for data transmission.

[0062] (7) When the interference signal power Pg > Ps and 10 * lg(Pg / Ps) > 15 dB, use a spatial domain filtering anti-jamming method to filter out interference.

[0063] If the interference signal is a periodic interference signal, let the period of the interference signal be Tg, and take the communication signal with a 50% duty cycle as an example, and the period is recorded as Ts;

[0064] (1) When Tg < 1 / 2 * Ts, use burst communication to avoid the periodic interference signal by adjusting the start and end times of the signal.

[0065] (2) When Tg ≥ 1 / 2 * Ts and is less than Ts, communicate by increasing the code rate, compress the message communication time, and reduce the signal exposure time in space.

[0066] For the spatial filtering method in step 6, this invention adopts a spatial anti-interference method using a binary array, the steps of which are as follows:

[0067] Step 1: Simultaneously use a fast frequency response variable local oscillator (local oscillator frequency is f) to convert the fq frequency hopping TX signals of both channels. nco According to TPTA, the current center frequency is f0 and the bandwidth is f. band The narrowband frequency-hopping pulse signal is down-mixed to the analog intermediate frequency to obtain two channels with frequency f. IF The analog intermediate frequency signal.

[0068] Here, the center frequency f0 of the frequency-hopping pulse is a discrete point, with a range of values ​​of... The minimum frequency interval is f band The maximum frequency interval is band-(num-1)·f band ; where f L f R `band` and `band` represent the left and right frequency boundaries and bandwidth of the actual operating frequency band of `fqTX`, respectively, and `num` represents the number of effective operating frequency bands within the operating frequency range. Here, `fix()` represents rounding down; the frequency f of the fast variable local oscillator. nco =f0-f IF ;

[0069] Step 2: Using a sampling rate of f sample The AD converter samples two analog intermediate frequency signals, f sample The signal bandwidth must be at least twice that of the signal band to obtain a digital intermediate frequency (IF) signal; a frequency control word is used to control a digital frequency synthesizer to generate a frequency of f. IF The sine and cosine signals are used to perform mixing processing on the digital intermediate frequency signal, transforming the intermediate frequency signal into a baseband complex signal.

[0070] Here, the bit width of the DDS frequency control word is n bits, where 8 ≤ n ≤ 64, and its value is calculated as follows: Among them, round() means rounding to the nearest integer, and dec2hex() means converting decimal to hexadecimal;

[0071] Step 3: Perform CIC filtering and R-fold speed reduction decimation on the baseband complex signal to obtain a low-speed baseband complex signal. Then, use an FIR low-pass filter to filter out useless out-of-band signals. At this point, two low-speed, spectrally clean baseband complex signals are obtained, denoted by din(1) and din(2) respectively.

[0072] The CIC filter design, based on a trade-off between performance and resources, employs a 4-stage cascaded design. The position of the first zero point on the frequency response curve of the CIC filter is chosen to be no greater than [missing value]. The maximum value within the range, the range of the deceleration factor R is: The passband flatness of the FIR filter should not exceed 1dB, and the out-of-band rejection should not be less than 40dB.

[0073] Step 4: Using the processed baseband signals din(1) and din(2) as input signals for spatial anti-interference processing, perform N-order data taps on each input signal, namely dex1(i) and dex2(i), i = 1, 2, ..., N, to obtain a total of 2N output signals, such as... Figure 5 As shown.

[0074] Step 5: In the absence of TX signal transmission and reception, preset the operating frequency of the hardware receiver section's control RF variable local oscillator to the first operating frequency point.

[0075] Step 6: Combine the 2N tapped data points into a column vector x. Calculate the statistical average of the instantaneous covariance matrix of vector x at point p to obtain the statistical mean of the covariance matrix. Where x = [dex1(1), dex1(2), ..., dex1(N), dex2(1), dex2(2), ..., dex2(N)] T , here,[] T This represents the transpose of a vector.

[0076] Step 7: Using the formula Solve for w. Where s∈C 2N×1 The first element of s is 1, and the rest are 0. for The inverse matrix, w is the spatial anti-interference weight, w∈C 2N×1 .

[0077] Step 8: Adjust to the working frequency. Repeat steps 6 and 7 at each working frequency. Adjust and traverse to all working frequencies to obtain the spatial anti-interference weights of all working frequencies.

[0078] Step 9: Under the normal transmission and reception state of TX signal, according to the current working frequency, perform a weighting operation on the data tap vector x of the current frequency and its corresponding anti-interference weight w to obtain the anti-interference output of the current frequency. Thus, steps 5-9 have completed the spatial anti-interference operation.

Claims

1. A binary array space-time-frequency code communication anti-interference method based on intelligent perception decision-making, characterized in that... It includes the following steps: Step 1: Divide the received signal into channels. Let the instantaneous bandwidth of the entire system be BW, and the number of frequency bands be determined. f1 represents the termination frequency, f0 represents the starting frequency, and the entire frequency band is divided into K sub-bands, with the center frequency of each sub-band being... Where k takes the values ​​0, 1, 2, ..., K; Step 2: Transform the divided channel from intermediate frequency signal to baseband signal. According to the Nyquist sampling theorem, the time-domain sampling of a signal is equivalent to the periodic extension of the signal's frequency domain by N*Fs. After spectrum shifting, the effective frequency within the range (0, Fs / 2) is Fc, and the bandwidth is BW; Step 3: After the communication terminal receives the signal, it uses the method of quadrature digital down-conversion to complete the baseband conversion of the intermediate frequency signal. When the data is collected by the ADC, first perform digital down-conversion. Generate a pair of quadrature carrier signals COS and SIN with a frequency of Fc by the DDS, and then mix the ADC data. Through mixing, the Fc intermediate frequency signal is mixed to the baseband. Since frequency doubling signals are generated during the mixing process, use a low-pass filter to filter out the frequency doubling signals; Finally, since the bandwidth of the interference signal to the main signal is unknown before spectrum calculation, secondary conversion is required. Use an integrator-comb filter cascade to complete signal rate extraction and further reduce the signal sampling rate; Step 4: Perform spectrum scanning according to the channel division, and perform multiple scans with a fixed bandwidth to obtain the signal spectrum information of the entire working frequency band; Step 5: Perform spectrum analysis according to the spectrum scanning results to obtain the number N of interference signal frequencies, the power magnitude Pg of the interference signal, and the period Tg of the interference signal; Step 6: Make intelligent decisions based on spectrum analysis and adjust the communication system mode. Let the actual number of operating frequencies be fn and the interference signal frequency be fg. The communication signal frequency is fs, and the communication signal power is Ps; The intelligent decision-making process in Step 6 is as follows: If the interference signal is a burst signal: (1) When N < 30%*fn, use a waveform with a code rate less than 2 / 3 for communication. At this time, the communication system can work normally without loss of key performance; (2) When 30%*fn ≤ N < 60%*fn, use a waveform with a code rate less than 1 / 3 for communication. At this time, the communication system can work normally without loss of key performance; (3) When 60%*fn ≤ N < 90%*fn, use a waveform with a code rate less than 1 / 10 for communication. At this time, the communication system can work normally without loss of key performance; (4) When 60%*fn ≤ N < 90%*fn, use a spatial anti-jamming method to filter out interference; (5) When N ≥ 100%*fn, use a spatial filtering anti-jamming method to filter out interference; (6) When the interference signal power Pg > Ps, and 10*lg(Pg / Ps) < 15dB, fn < fs, change the frequency hopping pattern to avoid interference or use spread-spectrum communication for data transmission; (7) When the interference signal power Pg > Ps, and 10*lg(Pg / Ps) > 15dB, use a spatial filtering anti-jamming method to filter out interference; If the interference signal is a periodic interference signal, let the interference signal period be Tg, and take the communication signal with a 50% duty cycle signal as an example, and the period is denoted as Ts; (1) When Tg < 1 / 2*Ts, use burst communication to avoid the periodic interference signal by adjusting the start and end times of the signal; (2) When 1 / 2*Ts ≤ Tg ≤ Ts, communicate by increasing the code rate, compress the message communication time, and reduce the signal exposure time in space.

2. The anti-jamming method for binary array space-time-frequency code communication based on intelligent perception decision according to claim 1, wherein: The spatial filtering method in step 6 employs a spatial anti-interference method using a binary array, and its steps are as follows: Step 6.1: For frequency hopping communication signals in two channels, use a local oscillator frequency with a rapidly variable frequency response. Based on the frequency hopping pattern, the current center frequency is... Bandwidth is The narrowband frequency-hopping pulse signal is down-mixed to the analog intermediate frequency to obtain two channels with frequencies of The analog intermediate frequency signal; Step 6.2: Utilize the sampling rate The AD converter samples two analog intermediate frequency signals. Not less than the signal frequency band bandwidth This is twice the size of the digital intermediate frequency signal; The frequency control word is used to control the digital frequency synthesizer to generate a frequency of The sine and cosine signals are used to perform frequency mixing on the digital intermediate frequency signal, transforming the intermediate frequency signal into a baseband complex signal; Step 6.3: Perform CIC filtering on the baseband complex signal. The low-speed baseband complex signal is obtained by decimation by a factor of 1, and the useless out-of-band signal is filtered out by an FIR low-pass filter. At this time, two low-speed baseband complex signals with clean spectrum are obtained, which are represented by din(1) and din(2) respectively. The CIC filter design, based on a trade-off between performance and resources, employs a 4-stage cascaded design. The position of the first zero point on the frequency response curve of the CIC filter is chosen to be no greater than [missing value]. Maximum value within the range, deceleration factor The range is 2 The passband flatness of the FIR filter should not exceed 1 dB, and the out-of-band rejection should not be less than 40 dB. Step 6.4: Use the two processed baseband signals din(1) and din(2) as input signals for spatial anti-interference processing. Perform M-order data taps on each input signal, namely dex1(i) and dex2(i), i=1,2,…,N, to obtain a total of 2N output signals; Step 6.5: When there is no communication signal transmission or reception, set the local oscillator operating frequency to the first operating frequency point; Step 6.6: Combine the 2N tapped data points into a column vector. For vectors conduct The statistical mean of the covariance matrix is ​​obtained by statistically averaging the instantaneous covariance matrices of the points. ,in, , ,here, Represents the transpose of a vector; Step 6.7: Using the formula Solve ;in, , The first element is 1, and the rest are 0. for The inverse matrix, For airspace anti-interference weights, ; Step 6.8: Adjust to the working frequency, repeat steps 6 and 7 at each working frequency, traverse all working frequencies, and obtain the spatial anti-interference weights of all working frequencies; Step 6.9: Under normal communication signal transmission and reception conditions, adjust the data tap vector of the current operating frequency according to the current operating frequency. Anti-interference weights corresponding to the current operating frequency The weighting operation is completed to obtain the anti-interference output of the current frequency point. Thus, steps 5-9 complete the spatial filtering anti-interference operation.

3. The binary array space-time-frequency code communication anti-interference method based on intelligent perception decision-making according to claim 2, characterized in that: In step 6.1, the center frequency of the narrowband frequency-hopping pulse For discrete points, the range of values ​​is... The minimum frequency interval is The maximum frequency interval is ;in, , , These represent the left and right frequency boundaries and the bandwidth of the actual operating frequency band of the narrowband frequency-hopping pulse. The number of effective operating frequency bands within the operating frequency range. , here Indicates rounding down; frequency of the rapidly variable local oscillator. .

4. The binary array space-time-frequency code communication anti-interference method based on intelligent perception decision-making according to claim 2, characterized in that: The bit width of the frequency control word of the DDS is Bit, The numerical calculation method is as follows ,in, This indicates rounding to the nearest integer. This indicates converting decimal to hexadecimal.

5. An electronic device, characterized in that, include: One or more processors; Memory; One or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, the one or more programs being configured to perform the method as described in any one of claims 1-4.

6. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores program code that can be invoked by a processor to execute the method as described in any one of claims 1-4.

Images