Binary array space-time-frequency code communication anti-interference method based on intelligent perception decision

Abstract

The invention provides a binary array space-time-frequency code communication anti-interference method based on intelligent sensing decision, which comprises the following steps of: monitoring and judging enemy interference signals by adopting a frequency spectrum sensing technology to obtain related attribute parameters of the enemy interference signals, and intelligently deciding an own communication mode according to differential characteristics of different interference signals, so that the anti-interference performance of the enemy interference signals is improved. And a binary array space-time-frequency code combined anti-interference technology is adopted, so that the communication anti-interference capability under various interference scenes including full-band ultra-wideband interference is improved. According to the invention, the survivability of the communication system in a complex electromagnetic interference environment is improved. A control strategy of interference perception and intelligent decision is introduced, and a space, time, frequency and code combined anti-interference method is adopted in a targeted manner according to differential characteristics of interference signals, so that the system is more intelligent and robust in signal receiving and processing, and the reliability of the system is enhanced.

Description

Technical Field

[0001] The invention relates to the field of frequency hopping communication anti-interference, and particularly to a space-time-frequency code combined communication anti-interference method. Background Art

[0002] In current data link communication technology, the basic anti-interference means is to rely on high-speed frequency hopping to avoid enemy communication interference and information interception. Due to the complexity and confidentiality of the frequency hopping pattern, it is difficult for the enemy to crack it directly. However, once the enemy obtains the specific range of its own communication frequency band through signal detection and other means, it can use ultra-wideband interference signals covering the entire communication frequency band to suppress high-power interference, causing its own communication system to fail to work normally. At the same time, the current SJL communication system does not have the ability to detect, perceive and analyze the characteristics of interference signals, and cannot intelligently adjust its own communication parameters based on the weaknesses of the enemy's interference methods, thereby limiting the improvement of the communication system's anti-interference ability.

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

[0004] In order to overcome the shortcomings of the prior art, the present invention provides a binary array space-time-frequency code communication anti-interference method based on intelligent perception decision-making. Spectrum sensing technology is used to listen and identify enemy interference signals, obtain their relevant attribute parameters, and intelligently decide on one's own communication mode based on the differentiated characteristics of different interference signals, and binary array space-time-frequency code joint anti-interference technology is used to improve the communication anti-interference capability in various interference scenarios including full-band ultra-wideband interference.

[0005] The technical solution adopted by the present invention to solve the technical problem comprises the following steps:

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

[0007] Step 2: Convert the intermediate frequency signal to the baseband signal of the divided channel. According to the Nyquist sampling theorem, the time domain sampling of the signal is equivalent to the periodic expansion of the frequency domain of the signal according to N*Fs, such as Figure 2 As shown, after spectrum shifting, the effective frequency in the range of (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. Figure 3 As shown in the figure, after the data is collected by the ADC, it is first digitally down-converted, and a pair of orthogonal carrier signals COS and SIN with a frequency of Fc are generated by the DDS. Then, the ADC data is mixed, and the Fc intermediate frequency signal is mixed to the baseband through mixing. Since the multiplied signal is generated during the mixing process, a low-pass filter (LPF) is used to filter out the multiplied signal. Finally, before the spectrum calculation, since the bandwidth of the interference signal to the main signal is unknown, a secondary frequency conversion is required, and the integral-comb cascade filter (CIC filter) is used to complete the signal rate extraction, and further reduce the signal sampling rate.

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

[0010] Step 5: Perform spectrum analysis based on the spectrum scanning 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. Assume that the actual number of working frequencies is fn, the interference signal frequency fg∈communication signal frequency fs, and the communication signal power is Ps. The intelligent decision-making process is as follows:

[0012] If the interference signal is a burst 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 less than 1 / 3 is used for communication. At this time, the communication system can work normally without loss of key performance;

[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 without loss of key performance;

[0016] (4) When 60%*fn≤N<90%*fn, the spatial domain 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 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;

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

[0020] If the interference signal is a periodic interference signal, assume the interference signal period is Tg, and take the communication signal with a 50% duty - cycle signal as an example, and its period is denoted as Ts;

[0021] (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;

[0022] (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.

[0023] The spatial - domain filtering method in step 6 adopts the spatial - domain 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 nco with a rapidly variable frequency response. According to the frequency - hopping pattern, down - mix the narrow - band frequency - hopping pulse signal with the current center frequency of f0 and bandwidth of f band to the analog intermediate frequency, and obtain two analog intermediate - frequency signals with a frequency of f IF ;

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

[0026] Step 2: Sample the two analog intermediate - frequency signals using an AD with a sampling rate of f sample and f sampleThe frequency control word is used to control the digital frequency synthesizer (DDS) to generate a frequency of f IF The sine and cosine signals are mixed with the digital intermediate frequency signal to transform the intermediate frequency signal into a baseband complex signal;

[0027] The bit width of the DDS frequency control word is n bits, 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 processing 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 baseband complex signals with clean spectrum are obtained, which are represented by din(1) and din(2) respectively.

[0029] The design of the CIC filter adopts a 4-stage cascade design based on the compromise between performance and resources. The first zero position of the frequency response curve of the CIC filter is selected to be no greater than The maximum value in the range, the range of the speed reduction ratio R is The passband flatness fluctuation of the FIR filter is no more than 1dB, and the out-of-band suppression is no less than 40dB;

[0030] Step 4: The processed two-channel baseband signals din(1) and din(2) are used as input signals for spatial domain anti-interference processing. Each input signal is tapped in N steps, namely dex1(i) and dex2(i), i = 1, 2, ..., N, to obtain a total of 2N output signals, such as Figure 5 shown.

[0031] Step 5: When there is no communication signal being sent or received, set the operating frequency of the local oscillator to the first operating frequency point;

[0032] Step 6: The 2N tapped data are combined into a column vector x, and the instantaneous covariance matrix of vector x at point p is statistically averaged to obtain the statistical average of the covariance matrix. in, here,[] T Represents the transpose of a vector;

[0034] 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 of , w is the spatial anti-interference weight, w∈C2N×1 ;

[0035] 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;

[0036] Step 9: When the communication signal is transmitted and received normally, according to the current working frequency, the data tap vector x of the current working frequency and the anti-interference weight w corresponding to the current working frequency are weighted to obtain the anti-interference output of the current frequency. At this point, steps 5 to 9 complete the spatial domain filtering anti-interference operation.

[0037] An electronic device comprises 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 execute the method as described above.

[0038] A computer-readable storage medium stores program codes, and the program codes can be called by a processor to execute the method as described above.

[0039] The beneficial effects of the present invention are:

[0040] (1) The present invention realizes rapid suppression processing of large-bandwidth interference in high-speed frequency-hopping communication, solves the problem that conventional communication systems are unable to resist full-band high-power suppression interference, and improves the survivability of communication systems in complex electromagnetic interference environments.

[0041] (2) The present invention introduces a control strategy for interference perception and intelligent decision-making. According to the differentiated characteristics of the interference signal, a joint anti-interference method of space, time, frequency and code is used in a targeted manner, making the system more intelligent and robust in signal reception and processing, thereby enhancing the reliability of the system. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] Figure 1 It is a channel division diagram of the real signal of the present invention.

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

[0044] Figure 3 It is a diagram of the orthogonal down-conversion process of the present invention.

[0045] Figure 4 It is the spectrum clipping diagram of the present invention.

[0046] Figure 5 1 is an output signal diagram of the present invention. DETAILED DESCRIPTION

[0047] The present invention is further described below in conjunction with the accompanying drawings and embodiments.

[0048] The technical solution adopted by the present invention comprises the following steps:

[0049] Step 1: Divide the received signal into channels, such as Figure 1 , assuming that the instantaneous bandwidth of the entire system is BW, the number of divided frequency bands K = (f1-f0) / BW, where f1 represents the end frequency, f0 represents the start frequency, the entire frequency band is divided into k sub-bands, and the corresponding center frequency of each frequency band is w k =k*BW, where the value of k is (0,1,2..K).

[0050] Step 2: Convert the intermediate frequency signal to baseband signal for the divided channels. According to the Nyquist sampling theorem, the time domain sampling of the signal is equivalent to the periodic expansion of the frequency domain of the signal according to N*Fs. Figure 2 As shown, after spectrum shifting, it can be seen that the effective frequency in the range of (0, Fs / 2) is Fc and the bandwidth is BW.

[0051] Step 3: After receiving the signal, the communication terminal uses orthogonal digital down-conversion to convert the intermediate frequency signal to baseband. Figure 3 As shown in the figure, after the data is collected by the ADC, it is first digitally down-converted, and the DDS generates a pair of orthogonal carrier signals COS and SIN, whose frequency is Fc. Then the ADC data is mixed, and the Fc intermediate frequency signal is mixed to the baseband through mixing. Since the mixing process will generate a multiplied signal, the multiplied signal can be filtered out by a low-pass filter. It should be noted that before the spectrum calculation, since the bandwidth of the interference signal to the main signal is unknown, a secondary frequency conversion is required, and the signal rate extraction is completed using CIC to further reduce the signal sampling rate.

[0052] Step 4: Perform spectrum scanning according to channel division, 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 spectrum aliasing of the signal. Since the bandwidth of the actual signal is only BW during the FFT calculation process, and the signal sampling rate is Fs, it can be seen from the spectrum shifting rule that 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 on the left side of the signal shifted at the sampling rate of (Fs), as shown in the shaded part of the figure, the two spectra are completely equivalent.

[0053] Therefore, the signal spectrum information of the entire working frequency band is obtained by using spectrum splicing.

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

[0055] Step 6: Make an intelligent decision based on the 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 decision-making process is as follows:

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

[0057] (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 operate normally without loss of key performance.

[0058] (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 operate normally without loss of key performance.

[0059] (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 operate normally without loss of key performance.

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

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

[0062] (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.

[0063] (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 the interference.

[0064] 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 denoted as Ts;

[0065] (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.

[0066] (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.

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

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

[0069] Here, the center frequency f0 of the frequency hopping pulse is a discrete point with a value range of The minimum frequency interval is f band , the maximum frequency interval is band-(num-1)·f band ; Among them, f L 、f R , band are the left frequency boundary, right frequency boundary and frequency band bandwidth of the actual working frequency band of fqTX, num is the number of effective working frequency bands within the working frequency range, Here fix() means rounding down; the frequency f of the fast variable local oscillator nco =f0-f IF ;

[0070] Step 2: Using the sampling rate f sample The AD samples the two analog intermediate frequency signals. sample Not less than 2 times of the signal bandwidth band, thereby obtaining a digital intermediate frequency signal; using the frequency control word to control the digital frequency synthesizer to generate a frequency of f IF The sine and cosine signals are mixed with the digital intermediate frequency signal to transform the intermediate frequency signal into a baseband complex signal.

[0071] Here, the bit width of the DDS frequency control word is n bits, 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.

[0072] Step 3: Perform CIC filtering and R-fold speed reduction 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 baseband complex signals with clean spectra are obtained, represented by din(1) and din(2) respectively.

[0073] The design of CIC filter adopts a 4-stage cascade design based on the compromise between performance and resources. The first zero position of the frequency response curve of CIC filter is selected to be no greater than The maximum value in the range, the range of the speed reduction ratio R is The passband flatness fluctuation of the FIR filter is no more than 1dB, and the out-of-band suppression is no less than 40dB.

[0074] Step 4: The processed two-channel baseband signals din(1) and din(2) are used as input signals for spatial domain anti-interference processing. Each input signal is subjected to N-order data tapping, namely dex1(i) and dex2(i), i = 1, 2, ..., N, to obtain a total of 2N output signals, such as Figure 5 shown.

[0075] Step 5: When there is no TX signal transmission or reception, control the operating frequency of the RF variable local oscillator of the hardware receiving section and preset it to the first operating frequency point.

[0076] Step 6: The 2N tapped data are combined into a column vector x, and the instantaneous covariance matrix of vector x at point p is statistically averaged to obtain the statistical average of the covariance matrix. in, here,[] T Represents the transpose of a vector.

[0078] 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 of , w is the spatial anti-interference weight, w∈C 2N×1 .

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

[0080] Step 9: When the TX signal is normally transmitted and received, according to the current working frequency, the data tap vector x of the current frequency and its corresponding anti-interference weight w are weighted to obtain the anti-interference output of the current frequency. At this point, steps 5 to 9 complete the spatial domain 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. Assume that the instantaneous bandwidth of the entire system is BW, and the number of divided frequency bands is K = (f1-f0) / BW, where f1 represents the end frequency and f0 represents the start frequency. The entire frequency band is divided into K sub-bands, and the center frequency of each corresponding frequency band is w k =k*BW, where k is 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 of (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. After the data is collected by the ADC, digital down-conversion is first performed. A pair of quadrature carrier signals COS and SIN with a frequency of Fc are generated by the DDS, and then the ADC data is mixed. Through mixing, the Fc intermediate frequency signal is mixed to the baseband. Since double-frequency signals are generated during the mixing process, a low-pass filter is used to filter out the double-frequency signals. Finally, before spectrum calculation, since the bandwidth of the interference signal to the main signal is unknown, secondary frequency conversion is required. An integrator-comb cascaded filter is used 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 based on the spectrum scanning results to obtain the number of interference signal frequencies N, the power magnitude Pg of the interference signal, and the period Tg of the interference signal. Step 6: Perform intelligent decision-making 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.

2. The anti-interference method for binary array space-time-frequency code communication based on intelligent perception decision according to claim 1, characterized in that: 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 space-domain anti-interference method to filter out interference. (5) When N ≥ 100% * fn, use a space-domain filtering anti-interference method to filter out interference. (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. (7) When the interference signal power Pg > Ps and 10 * lg(Pg / Ps) > 15 dB, use a space-domain filtering anti-interference method to filter out interference.

3. The anti-interference method for binary array space-time-frequency code communication based on intelligent perception decision according to claim 1, characterized in that: The intelligent decision-making process in step 6 is as follows: If the interference signal is a periodic interference signal, let the interference signal period be Tg, and take the communication signal as an example of a 50% duty cycle signal, and the period is recorded as Ts; (1) When Tg < 1 / 2*Ts, burst communication is used to avoid periodic interference signals by adjusting the start and end time of the signal; (2) When 1 / 2*Ts≤Tg≤Ts, communication is performed by increasing the code rate, compressing the message communication time and reducing the signal exposure time in space.

4. The method for anti-interference of binary array space-time-frequency code communication based on intelligent perception decision according to claim 2 is characterized in that: The spatial domain filtering method in step 6 adopts a spatial domain anti-interference method of a binary array, and the steps are as follows: Step 6.1: For the two-channel frequency hopping communication signal, use a local oscillator frequency f with a fast variable frequency response. nco According to the frequency hopping pattern, the current center frequency is f0 and the bandwidth is f band The narrowband frequency-hopping pulse signal is mixed down to the analog intermediate frequency to obtain two channels with a frequency of f IF The analog intermediate frequency signal; Step 6.2: Using the sampling rate f sample The AD samples the two analog intermediate frequency signals. sample Not less than 2 times of the signal frequency band bandwidth, thereby obtaining a digital intermediate frequency signal; using the frequency control word to control the digital frequency synthesizer to generate a frequency of f IF The sine and cosine signals are mixed with the digital intermediate frequency signal to transform the intermediate frequency signal into a baseband complex signal; Step 6.3: Perform CIC filtering and R-fold speed reduction processing 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 baseband complex signals with clean spectrum are obtained, which are represented by din(1) and din(2) respectively; The design of the CIC filter adopts a 4-stage cascade design based on the compromise between performance and resources. The first zero position of the frequency response curve of the CIC filter is selected to be no greater than The maximum value in the range, the range of the speed reduction ratio R is The passband flatness fluctuation of the FIR filter is no more than 1dB, and the out-of-band suppression is no less than 40dB; Step 6.4: The processed two baseband signals din(1) and din(2) are used as input signals for spatial domain anti-interference processing, and each input signal is subjected to N-order data tapping, namely dexl(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 being sent or received, set the operating frequency of the local oscillator to the first operating frequency point; Step 6.6: The 2N tapped data are combined into a column vector x, and the instantaneous covariance matrix of vector x at p points is statistically averaged to obtain the statistical average 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; Step 6.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 of , w is the spatial anti-interference weight, w∈C 2N×1 ; 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: When the communication signal is transmitted and received normally, according to the current working frequency, the data tap vector x of the current working frequency and the anti-interference weight w corresponding to the current working frequency are weighted to obtain the anti-interference output of the current frequency. At this point, steps 5 to 9 complete the spatial domain filtering anti-interference operation.

5. The method for anti-interference of binary array space-time-frequency code communication based on intelligent perception decision according to claim 4 is characterized in that: In step 6.1, the center frequency f0 of the narrowband frequency hopping pulse is a discrete point with a value range of The minimum frequency interval is f band , the maximum frequency interval is band-(num-1)·f band ; Among them, f L 、f R , band are the left frequency boundary, right frequency boundary and frequency band bandwidth of the actual working frequency band of fqTX, num is the number of effective working frequency bands within the working frequency range, Here fix() means rounding down; the frequency f of the fast variable local oscillator nco =f0-f IF .

6. The method for anti-interference of binary array space-time-frequency code communication based on intelligent perception decision according to claim 4 is characterized in that: The bit width of the frequency control word of the DDS is n bits, 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.

7. 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, and the one or more programs are configured to execute the method according to any one of claims 1 to 6.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores program codes, and the program codes can be called by a processor to execute the method according to any one of claims 1 to 6.

Images