CN109639303A - A kind of Interference Detection and suppressing method based on windowing process - Google Patents
A kind of Interference Detection and suppressing method based on windowing process Download PDFInfo
- Publication number
- CN109639303A CN109639303A CN201811638158.6A CN201811638158A CN109639303A CN 109639303 A CN109639303 A CN 109639303A CN 201811638158 A CN201811638158 A CN 201811638158A CN 109639303 A CN109639303 A CN 109639303A
- Authority
- CN
- China
- Prior art keywords
- windowing
- interference
- digital signal
- sequence
- signal
- 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.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 69
- 230000008569 process Effects 0.000 title abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 38
- 238000004891 communication Methods 0.000 claims abstract description 36
- 230000009466 transformation Effects 0.000 claims abstract description 21
- 230000001629 suppression Effects 0.000 claims description 51
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 230000001131 transforming effect Effects 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 abstract description 12
- 230000002829 reductive effect Effects 0.000 abstract description 10
- 230000003595 spectral effect Effects 0.000 description 24
- 238000010586 diagram Methods 0.000 description 7
- 238000010187 selection method Methods 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 230000000875 corresponding effect Effects 0.000 description 5
- 230000008030 elimination Effects 0.000 description 5
- 238000003379 elimination reaction Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 108010076504 Protein Sorting Signals Proteins 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005562 fading Methods 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000006735 deficit Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
- H04B1/1027—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
- H04B1/1027—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal
- H04B1/1036—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal with automatic suppression of narrow band noise or interference, e.g. by using tuneable notch filters
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/16—Circuits
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Noise Elimination (AREA)
Abstract
The invention discloses a kind of Interference Detection and suppressing method based on windowing process.The method comprising the steps of frequency-domain transform, Interference Detection, AF panel and time domain transformation, overlapping windowing process has wherein been carried out after frequency-domain transform, it has then carried out going overlap processing after time domain transformation, in addition it can including carrying out ambient noise estimation to transformed frequency domain digital signal and adaptive determining interference detection threshold in Interference Detection, the accuracy to interferer signal detection is improved, the noise specific loss of windowing process bring is thus also reduced.This method is versatile, especially suitable for the interferer signal detection in satellite communication, microwave communication.
Description
Technical Field
The present invention relates to the field of communication signal processing technologies, and in particular, to an interference detection and suppression method based on windowing for an interference signal.
Background
In radio communications such as satellite communications and mobile communications, various interference signals are often mixed into a communication signal when the communication signal is received by a receiver, and detection and cancellation of these interference signals are advantageous for correct reception of the communication signal.
In the prior art, the signal-to-noise ratio loss of communication signals occurs in the process of detecting and eliminating interference signals, and the threshold setting of interference detection is difficult to achieve optimization due to the fact that background noise cannot be accurately estimated.
Disclosure of Invention
The invention mainly solves the technical problem of providing an interference detection and suppression method based on windowing processing, and solves the problems of inaccurate background noise detection, signal to noise ratio reduction and the like existing in the detection and elimination of interference signals in the prior art.
In order to solve the above technical problem, one technical solution adopted by the present invention is to provide an interference detection and suppression method based on windowing processing, including the following steps: the method comprises the steps of frequency domain transformation, namely transforming a communication signal mixed with an interference signal into an input time domain digital signal through an AD sampler, windowing the input time domain digital signal, and then transforming the input time domain digital signal into a frequency domain digital signal through FFT; interference detection, namely performing background noise estimation on the frequency domain digital signal and then determining an interference detection threshold; interference suppression, wherein the component of the frequency domain digital signal which is greater than the interference detection threshold is regarded as an interference signal component, and the interference signal component is suppressed; and time domain transformation, namely, after the interference signal component in the frequency domain digital signal is suppressed, carrying out IFFT transformation on the frequency domain digital signal to restore the frequency domain digital signal into an output time domain digital signal.
In another embodiment of the interference detection and suppression method based on windowing processing, the windowing processing uses a window function including windowing processing of the input time-domain digital signal using a buttlet window, a hanning window, a Hamming window, a Blackman window, or a Blackman-Karris window.
In another embodiment of the interference detection and suppression method based on windowing, the input time-domain digital signal is a sequence x (k):
x(k)＝Ap(k)+n(k)
wherein, p (k) is a PN sequence with equal probability value of +/-l, and the length of the PN sequence is N; n (k) is zero mean and varianceA is the amplitude of the communication signal;
performing windowing on x (k), wherein the window function is w (k), and the sequence after windowing is as follows:
xw(k)＝Ap(k)w(k)+n(k)w(k)。
in another embodiment of the interference detection and suppression method based on windowing processing, the windowing processing is overlap windowing processing, and the time domain digital signal output by IFFT is subjected to overlap removing processing in the time domain transformation.
In another embodiment of the interference detection and suppression method based on windowing processing, the data sequence of the input time domain digital signal is equally segmented, the length of each data segment is N, the overlapping length of the front and rear adjacent data segments is Nr during overlapping windowing processing, and r (0 ≦ r < 1) is an overlapping factor.
In another embodiment of the interference detection and suppression method based on windowing, the overlap factor r is 1/2, and the input time-domain digital signal is a sequence x (k) ═ ap (k) + N (k), where p (k) is a PN sequence with equal probability value ± l, and its length is N; n (k) is zero mean and varianceA is the amplitude of the communication signal;
the output sequence after the overlapping windowing is
Wherein,w (k) is a window function.
In another embodiment of the interference detection and suppression method based on windowing, the overlap factor r is 1/2, and the input time-domain digital signal is a sequence x (k) ═ ap (k) + N (k), where p (k) is a PN sequence with equal probability value ± l, and its length is N; n (k) is zero mean and varianceA is the amplitude of the communication signal;
the output sequence after the overlapping windowing is
Wherein,w (k) is a window function.
In another embodiment of the interference detection and suppression method based on windowing, the input time domain digital signal is divided into two input sequences, the first input sequence is subjected to windowing, and then subjected to FFT, interference detection, interference suppression and IFFT to obtain a first output sequence, the second input sequence is subjected to N/2 delay and then windowing, and then subjected to FFT, interference detection, interference suppression and IFFT to obtain a second output sequence, and the first output sequence and the second output sequence are added to obtain the output time domain digital signal.
The invention has the beneficial effects that: the invention discloses an interference detection and suppression method based on windowing processing. The method comprises the steps of frequency domain transformation, interference detection, interference suppression and time domain transformation, wherein overlapping windowing is carried out after the frequency domain transformation, overlap removing processing is carried out after the time domain transformation, in addition, background noise estimation and self-adaption determination of an interference detection threshold can be carried out on a transformed frequency domain digital signal in the interference detection, the accuracy of interference signal detection is improved, and the signal-to-noise ratio loss brought by the windowing processing is reduced. The method has strong universality and is particularly suitable for detecting the interference signals in satellite communication and microwave communication.
Drawings
FIG. 1 is a schematic block diagram illustrating an embodiment of a method for interference detection and suppression based on windowing in accordance with the present invention;
fig. 2 is a flow diagram of another embodiment of an interference detection and suppression method based on windowing in accordance with the present invention;
FIG. 3 is a block diagram of an overlay windowing process in accordance with another embodiment of the interference detection and suppression method based on windowing;
FIG. 4 is an illustration of an overlap-windowing process in accordance with another embodiment of the interference detection and suppression method based on windowing in accordance with the present invention;
FIG. 5 is a schematic diagram of an overlap-windowing process according to another embodiment of the interference detection and suppression method based on windowing of the present invention;
fig. 6 is a flow diagram of another embodiment of an interference detection and suppression method based on windowing in accordance with the present invention;
FIG. 7 is a graph of interference rejection threshold versus bit error rate for another embodiment of a method for interference detection and rejection based on windowing in accordance with the present invention;
FIG. 8 is a graph of trap magnitude versus bit error rate for another embodiment of a method for interference detection and suppression based on windowing in accordance with the present invention;
FIG. 9 is a graph of windowing overlap versus bit error rate according to another embodiment of the interference detection and suppression method based on windowing in accordance with the present invention;
fig. 10 is a schematic block diagram of another embodiment of an interference detection and suppression method based on windowing according to the present invention.
Detailed Description
In order to facilitate an understanding of the invention, the invention is described in more detail below with reference to the accompanying drawings and specific examples. Preferred embodiments of the present invention are shown in the drawings. This invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete.
It is to be noted that, unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
With reference to the schematic diagram of fig. 1 and the flowchart of fig. 2, a flowchart of an embodiment of an interference detection and suppression method based on windowing according to the present invention is disclosed. The method comprises the following steps:
first step S101: the method comprises the steps of frequency domain transformation, wherein a communication signal mixed with an interference signal is transformed into an input time domain digital signal through an AD sampler, and the input time domain digital signal is transformed into a frequency domain digital signal after FFT transformation;
second step S102: interference detection, namely performing background noise estimation on the frequency domain digital signal and then determining an interference detection threshold;
the third step S103: and interference suppression, wherein the component of the frequency domain digital signal which is greater than the interference detection threshold is regarded as an interference signal component, and the interference signal component is suppressed to eliminate interference signal energy.
Here, for the detected interference signal component, the amplitude of the interference signal component needs to be modified, which is also called threshold processing, and this process inevitably brings loss to the useful communication signal while processing the interference, resulting in a reduction in the output signal-to-noise ratio.
The fourth step S104: and time domain transformation, namely, after the interference signal component in the frequency domain digital signal is suppressed, carrying out IFFT transformation on the frequency domain digital signal to restore the frequency domain digital signal into an output time domain digital signal.
Preferably, the windowing process is further performed on the input time-domain digital signal before the frequency-domain transform in the first step S101.
Here, it is an important improvement to perform windowing, and if the windowing is not performed, performing an N-point FFT operation is equivalent to a rectangular window with an N-point added, and its first side lobe is only 13dB lower than the main lobe, i.e. the degree of side lobe suppression is only-13 dB, and for an interference signal several tens of dB greater than the communication signal, its side lobe is also much greater than the communication signal. Thus, from the frequency domain, the entire communication signal frequency domain is contaminated by interference. Therefore, windowing is required to focus more energy on the main lobe of the communication signal and reduce the amplitude of the side lobe. In time domain, the windowing process is to weight the input data, and the window function coefficients are gradually attenuated from the center to the two ends, so that the two ends of the data segment are smooth, and the frequency spectrum leakage is reduced. However, the window function is attenuated towards both ends, so that the input signal is distorted, and additional signal-to-noise ratio loss is brought.
Preferably, x (k) is the input time domain digital signal sequence:
x(k)＝Ap(k)+n(k)
wherein, p (k) is a PN sequence with equal probability value of +/-l, and the length of the PN sequence is N, which indicates that the sequence is a direct sequence spread spectrum signal; n (k) is zero mean and varianceAnd a is the signal amplitude.
Windowing x (k), the window function being w (k), the sequence after windowing:
xw(k)＝Ap(k)w(k)+n(k)w(k)
for xw(k) And performing correlated despreading with the length of N and integrating to obtain:
mean and variance of the sequence z:
the signal-to-noise ratio of the windowed sequence is:
the correlated output signal-to-noise ratio without windowing is:
the signal-to-noise ratio loss due to windowing is:
it can be seen that the loss of the snr is related to the coefficients of the window functions, and that different window functions bring different snr losses.
Table 1 shows the characteristics of several window functions and the resulting snr loss, and also compares the relative sidelobe peak amplitude and main lobe bandwidth of each window function. At a given length, the rectangular window has its main lobe narrowest, but its relative side lobe amplitude is greatest. From top to bottom, the main lobe of the various windows becomes wider and the relative side lobe amplitudes become smaller.
TABLE 1 fundamental characteristics of the windowing function and resulting loss of SNR
Most energy of the narrow-band interference signal can be limited in a plurality of limited spectral lines by selecting a window function with lower side lobes, so that the number of spectral lines needing to be suppressed is reduced, and the influence of interference on a useful signal is reduced to the greatest extent; meanwhile, the distortion of the useful signal can be reduced by selecting a smaller main lobe broadband.
The choice of the window function in practical applications requires consideration of the trade-off between sidelobe attenuation and main lobe bandwidth. Since the lower the side lobe is, the wider the main lobe is, and the more the impairment to the useful signal is at the same time as suppressing the narrowband interference, when selecting the window function, an appropriate window function is selected in accordance with the dynamic range of the received signal of the receiver and the intensity of the narrowband interference that needs to be suppressed.
Further preferably, in order to compensate for the snr loss caused by windowing, the frequency domain digital signal needs to be subjected to overlap windowing. The overlap-and-window process has the advantage of reducing the loss of signal-to-noise ratio due to windowing, at the cost of increased computational complexity. And, the main factors affecting the snr loss are the window function type, the window length, and the windowing overlap.
From the selection of the type of the window function, when the relative amplitude of the side lobe of the window function is low, the main lobe bandwidth is wider, more frequency points need to be processed, and when the main lobe bandwidth is narrow, the relative amplitude of the side lobe is high, the frequency spectrum leakage is serious, and more frequency points are polluted, so when the window function is selected, the dynamic range of the received signal and the intensity of the narrow-band interference signal needing to be suppressed need to be considered, a proper window function is selected, and the loss of the signal-to-noise ratio is reduced. For the windowing overlapping degree, the larger the windowing overlapping proportion is, the smaller the signal-to-noise ratio loss is, but the implementation complexity is high, the hardware resource consumption is high, and compromise consideration is required.
In further correspondence with the overlap windowing process, in conjunction with fig. 3, it is preferable that the output time-domain digital signal is also subjected to a de-overlap process after the time-domain transform at the fourth step S104.
The windowing causes distortion of the desired signal and the input data is overlap windowed in order to compensate for the loss of signal to noise ratio caused by the windowing. Assuming that the data length into which the signal sequence is segmented is N, the overlap factor r (0 ≦ r < 1), the overlap windowing principle is shown in FIG. 4. Here, the data sequence inputted after windowing is segmented, the data length of each segment is N, and when overlapping is performed, the data length of each segment is also N, but the two segments of data overlap each other, and the overlap length is Nr. Therefore, another problem to be considered when performing overlap windowing on signals is how to combine two paths of data into one path of data as a final output result.
Preferably, there are two methods for the overlap windowing process: selection methods and addition methods. For example, 1/2 overlap windowing is performed, and the selection method is to take 1/2 blocks of data at the center of the window of each signal, discard 1/4 blocks of data on the left and right sides, and combine 1/2 blocks of data on the upper and lower sides into a complete data sequence. The addition rule is to add the overlapping portions of two pieces of data as the final output signal.
The loss in signal-to-noise ratio of the overlap output of the selection and addition method with the overlap factor of 1/2 is analyzed below. Preferably, in the selection method output mode, the output sequence after the overlapping windowing is as follows:
wherein
w1(k) The signal weighted value under the output of the selection method is known from the derivation result of the signal-to-noise ratio loss, and the signal-to-noise ratio loss under the output mode of the selection method is as follows:
preferably, in the additive output mode, the output sequence after the overlapping windowing is:
wherein
w2(k) Is the weighted value of the signal under the output of the summation, the loss of the signal-to-noise ratio of the summation is as follows:
table 2 shows the snr loss due to the selection and summation outputs at different overlap factors.
TABLE 2 SNR loss for different overlap factors for two data synthesis modes
As can be seen from table 2, the same overlap factor results in a lower signal-to-noise ratio for the summation method than for the selection method, and this difference can be intuitively explained from the time domain. Considering hardware implementation, the addition output is more than the selection output by Nxr times of addition operation; for a non-rectangular window, theoretically, the larger the overlap ratio of two consecutive segments of data is, the smaller the introduced windowing loss is, the larger the corresponding computation amount is, and the implementation is not easy, and the selection of the overlap ratio in practical application depends on hardware conditions and system design performance requirements.
Preferably, the selected addition and subtraction factor is 1/2, and another embodiment of the interference detection and suppression method based on overlap windowing is shown in fig. 5, where the input time domain digital signal is divided into two input sequences, a first input sequence x (N) is subjected to windowing, and then to FFT conversion, interference detection, interference suppression and IFFT conversion to obtain a first output sequence, a second input sequence is subjected to N/2 delay x (N + N/2) and then to windowing, and then is subjected to FFT conversion, interference detection, interference suppression and IFFT conversion to obtain a second output sequence, and the first output sequence and the second output sequence are added to obtain the output time domain digital signal. Here, the first path is a waveform obtained by windowing a single-path signal, and the second path is a waveform obtained by delaying and windowing a signal x (N) by N/2. The data fidelity is highest in the center of the window, and the closer to the ends of the window, the more severe the distortion. However, the two signals are complementary, and the place with the most serious distortion of one signal can be compensated by the other signal. By time domain overlap windowing, useful signals can be better recovered when IFFT conversion is performed after interference cancellation. Preferably, Blackman overlap windowing with a length N of 1024 is chosen to verify the effect of overlap windowing on signal distortion reduction, with an overlap factor of 1/2, the number of overlap windowing being 2.
Further preferably, in the interference detection, the background noise estimation is the most basic condition for the interference detection, and the background noise power is set to be highIs small asAnd performing FFT after windowing, so that the components of the interference frequency are concentrated in the main lobe. The frequency bands of several interfering signals do not occupy more than half of the signal bandwidth and the useful communication signal is small and buried in noise, so that the frequency components of white gaussian noise are dominant in the frequency band without interference. Mathematical expectation is zero and variance isThe one-dimensional distribution of the envelope of the stationary gaussian band-limited noise follows rayleigh distribution. Since the frequency component of Gaussian white noise follows a normal distribution with a variance ofTherefore, the amplitude of the frequency component obeys Rayleigh distribution, the mean value of the amplitude of the whole frequency component can be obtained according to the statistical characteristics of the Rayleigh distribution, the amplitude variance of the frequency domain component of the background noise is estimated according to the relation between the mean value of the Rayleigh distribution and the normal distribution variance, and the time domain variance, namely the background noise power, is obtained according to the property of Fourier transform.
Here, assuming that the amplitude x of the frequency component of gaussian white noise follows rayleigh distribution, its probability density function:
probability distribution function of amplitude x:
the mean of the amplitudes x is:
the frequency component data with 1/2 ratio is taken from small to large, F (A) is 1/2, if the number of points obtained in the frequency domain is large, half of the minimum points are selected, and the frequency component data falls in the interval [0, A]In (1). It can be calculated that the variable x is in the interval [0, A]The mean value of (d) is recorded as μ1I.e. the mean of the minimum points of half of the frequency components, comparing mu and mu1The relationship between the two can be obtained, and then the average value of the whole frequency component can be deduced from the undisturbed frequency component in practical application. When F (A) is 1/2, obtainμ1Comprises the following steps:
relationship u/u between half minimum point mean of frequency component and mean of ensemble13.4328, i.e. obtaining mu1Then multiplied by 3.4328 to calculate mu. Therefore, the average of the entire frequency component is 3.4328 times the average of the minimum point of half of the frequency component, which is derived from the statistical properties of the rayleigh distribution.
Then according to the mean value of Rayleigh distributionσFIs the mean square error of the normal distribution of the frequency components, based onThe background noise power is obtained.
Thus, μ can be determined first and then1Then find σ againFFinally, find outAt this point, the estimation of the background noise is completed.
Further preferably, the interference detection threshold design directly affects the anti-interference performance, and the design of the interference detection threshold is the key of the whole interference detection and elimination. On the one hand, the actual interference and signal are time-varying, so the interference detection threshold should not be fixed, and the threshold should be selected based on the statistical characteristics of the received signal. On the other hand, in order to overcome the disadvantage that the threshold detection algorithm is not strong in adaptability to the environment, it is necessary to be able to dynamically and adaptively detect the interference signal by multi-threshold interference.
Preferably, the mean value is obtained by counting x (k) after FFT of the input time domain digital signalSum mean square errorN is the number of spectral lines, ζnThe threshold of the first moment algorithm is defined as Th mu + β sigma, and β is a threshold optimization factor, which can be selected from a preset set of weighting factors according to different channel environments (such as fading, multipath delay, etc.).
When an interference signal with power far larger than that of a spread spectrum signal and background noise exists in an input time domain digital signal, through FFT (fast Fourier transform), the output of a sub-band with interference does not meet Gaussian distribution, the output power of the sub-band polluted by strong interference is far larger than that of the sub-band polluted by weak interference or not polluted by interference, and the missed detection of part of the interfered sub-band is caused, so that the signal-to-noise ratio of an output decision variable is reduced, and the error rate of a system is increased. In addition, the common interference detection algorithm has low adaptability to the environment, and the characteristics of the channel and the interference and the conditions of multipath delay, fading and the like of the signal in the transmission process need to be estimated in advance when the threshold value is selected each time.
The self-adaptive threshold needs to meet two conditions, namely when no interference signal exists, the self-adaptive threshold value is higher than spectral lines of most communication signals, and useful spectral lines of the communication signals cannot be used as interference processing, so that the false alarm probability is minimum; secondly, when there is an interference signal, the adaptive threshold value should be lower than the spectral lines of all interference signals and higher than the spectral lines of most communication signals, so that the detection probability is maximum.
To satisfy these two conditions, the distribution of the spectral lines is first analyzed. Let x (n) ═ s (n) + n (n), s (n) be a signal, and n (n) be noise. The adaptive multi-threshold interference detection method here considers that the desired signal s (n) in the received signal is swamped by the channel noise n (n), for example, for a direct sequence spread spectrum signal. The DFT of the N-point sequence x (N) is defined as:
according to the property of fourier transform, the output spectral line is x (k) ═ s (k) + N (k), s (k) is DFT of N-point sequence s (N), and N (k) is DFT of N-point sequence N (N). The DFT transform may also be implemented by an FFT transform. Without interference, | X (k) non-woven2＝|S(k)+N(k)|2Subject to an exponential distribution with a parameter λ, havingSetting the adaptive threshold value as Th, then | X (k) & gtY2The probability p lower than the threshold Th is:
it can be calculated that when Th is 1/λ, p is 0.6321; when Th is 2/lambda, p is 0.8647; when Th is 3/lambda, p is 0.9502; when Th is 4/lambda, p is 0.9817; when Th is 5/λ, p is 0.9933.
Assuming that the output signal is a superposition of the useful communication signal and noise, the square of the amplitudes of the N spectral lines after the discrete fourier transform should follow an exponential distribution with a parameter λ, according to the above analysis. In the case of a non-interfering signal, the probability of the square of the N spectral line amplitudes being less than 5/λ is 0.9933, i.e. the probability of being greater than 5/λ is only 0.0067. It can be considered that the spectral lines of the N spectral lines whose square of the amplitude is greater than 5/λ are almost nonexistent. Then if there is interference, the spectral line whose square of the amplitude of the spectral line after the discrete fourier transform is greater than 5/λ can be considered as the spectral line containing the interference, and the interference is processed by using a zeroing or clamping algorithm.
However, in practical tests it has been found that when the gate is limited to the theoretical optimum value of 5/λ, the probability of treating the communication signal as interference is almost zero, i.e. the false alarm probability is lowest. However, after the strong interference is subjected to the clamping processing, the interference suppression effect is not thorough, and the detection probability of the weak interference is reduced. Therefore, from comprehensive factors such as false alarm probability, strong interference suppression effect, weak interference detection probability and the like, limiting the gate to 3/lambda is a compromise choice. Therefore, if the discrete fourier transform of the N-point input sequence x (N) is known as x (k), then in practical applications, when the number of points N of the DFT transform is large, the estimate of the statistical mean 1/λ may be replaced by the mean of the sum of squared spectral line amplitudes, and the adaptive threshold may be:
in addition, there is a problem that the first interference threshold determination in the dynamic adaptive multi-threshold interference detection technology cannot be performed by itself, because the above analysis assumes that there is no interference signal or at least needs to ensure that there is no strong interference, but actually the signals are all superimposed with interference signals, so other means are needed for the first interference threshold determination.
In this embodiment, a simplified approximate background noise estimation algorithm is used to assist in determining the first interference threshold. Due to the desired signal s (n) in the received signal being channeledThe noise n (n) is drowned, i.e. S (k) < N (k) in the frequency domain, μ is the mean of the entire undisturbed frequency components obtained by the background noise estimation algorithm, so the first adaptive threshold can be Th05 mu. Thus, the method at least ensures that most interference is suppressed, and further determines the threshold by using a dynamic self-adaptive multi-threshold interference detection algorithm. Fig. 6 shows a flow of an adaptive multi-threshold interference detection algorithm based on background noise estimation.
As shown in fig. 6, when the input time domain digital signal is FFT transformed to be represented as a sequence y (k), the background noise estimation is performed on y (k), and then the first detection threshold Th is determined0And 5 μ, then comparing the amplitude of each frequency component in the sequence y (k) with a threshold value of 5 μ, and if the amplitude of a frequency component is greater than 5 μ, performing interference suppression processing on the frequency component, for example, performing amplitude trapping processing on the amplitude of the frequency component, and if the amplitude of the frequency component is less than or equal to 5 μ, performing no interference suppression processing. Therefore, each frequency component in Y (K) is subjected to threshold detection and interference suppression for the first time to obtain a sequence X (K), and then adaptive interference threshold detection is carried out, wherein the detection threshold at this time is as follows:then comparing the amplitude of each frequency component in the sequence X (K) with a threshold value Th, when the amplitude of the frequency component is larger than Th, performing interference suppression processing on the frequency component, then re-encoding the frequency component subjected to the interference suppression processing into the sequence X (K), continuing to compare the frequency component with the current adaptive detection threshold Th (since the frequency component in the sequence X (K) is transformed, the threshold is correspondingly and adaptively adjusted), until the amplitude of all the frequency components is smaller than or equal to the current adaptive detection threshold Th, then performing IFFT conversion on the sequence X (K) and outputting the result.
Therefore, the interference detection problem is converted into a hypothesis test problem by the dynamic self-adaptive multi-threshold interference detection, a narrow-band interference component does not exist in a hypothesis receiving signal, then whether the interference component exists in the receiving signal is tested by setting a threshold value for the distribution obeyed by the transformed frequency spectrum, and the first interference detection threshold is determined by simplified background noise estimation. If the hypothesis is found to be established through detection, if no narrow-band interference exists in the received signal, the IFFT can be directly carried out on the spectral line, the spectral line is converted into a time domain, and then the data is sent to the next stage for processing. And if the detection finds that the hypothesis is not true, the interference exists in the received signal, the interference elimination processing is carried out on the spectral line of which the spectral line value is greater than the threshold value, then the average value of the processed spectral line is recalculated, and a new threshold is set to detect the spectral line until the spectral line of which the hypothesis is not true does not exist. Therefore, the interference frequency point is detected by setting the threshold in a multi-time self-adaptive manner, and the interference is removed by an interference elimination algorithm, so that the interference suppression is realized.
Further, the embodiments shown in fig. 6 and fig. 5 can be combined, so that the adaptive threshold detection can be performed on both paths of sequences in the process of the superimposition window according to the method shown in fig. 6.
Preferably, the effect of the above processing method can be checked by simulation. Fig. 7 shows the effect of different interference detection thresholds on the bit error rate. Taking the interference detection threshold as a variable, and analyzing the interference detection threshold to be 0-10 sigmaFThe corresponding bit error rate. As can be seen from fig. 7, the interference detection threshold takes 4 σFThe time error code rate is minimum. From the "3 σ rule" of normal distribution, it is known that Gaussian random variables are distributed in the interval (-3 σ) with a probability of 95%F,+3σF) In this case, the frequency component of the background noise is mainly distributed in the region (-3 σ)F,+3σF) When the threshold value of interference detection is less than 4 sigmaFWhen the frequency component of the partial communication signal is judged as interference and eliminated, the error rate is increased, and when the threshold value of the interference detection is 4 sigmaFTo 10 sigmaFWhen the interference detection threshold value is within the range of (3), the interference detection threshold value is larger than the amplitude value of the frequency component of the communication signal, the frequency component of the useful communication signal which is not interfered is reserved, and the range is smaller than the amplitude of the interference signal, and the frequency component of the interference signal can be effectively eliminated, so the interference elimination threshold value is 4 sigmaFTo 10σFWhen the bit error rate is within the range of (2), the bit error rate is within a normal range.
Fig. 8 shows the error rate of the clipping process for interference in interference suppression. The amplitude of the interference frequency component is analyzed to be 0 to 30 sigma by using the trap amplitude value as a variableFThe corresponding error rate is shown in fig. 8, where the trap amplitude is equal to the background noise amplitude, the error rate is within the normal range, and when the trap amplitude is 2 σFTo 6 sigmaFWhen the trap amplitude is larger than 10 sigma, the error rate is minimumFIn the process, the interference power is eliminated too little, and the error rate rapidly becomes larger and exceeds the normal range along with the increase of the trap amplitude value.
Fig. 9 shows the effect of the windowing overlap degree on the error rate, and the error rate is analyzed when the windowing overlap degree is 10% to 100% by taking the windowing overlap degree as a variable, and as can be seen from fig. 9, the error rates corresponding to the windowing overlap degree in the range of (0, 50%) and the range of (50%, 100%) are symmetrical, and the error rate is the smallest when the windowing overlap degree is 1/2. Due to the influence of windowing, the signal time domain waveform presents the shape of a window, the output signal subjected to interference suppression still presents the current situation of the window, when the windowing overlapping degree is 1/2, two paths of signals are orthogonal, and the signal time domain waveform after addition becomes gentle. Corresponding to the previously described embodiment shown in fig. 5.
Therefore, based on the simulation analysis of fig. 7 to 9, when the interference cancellation threshold value takes 4 σFThe trap amplitude is 2 sigmaFTo 6 sigmaFAnd when the windowing overlapping degree is 1/2, the algorithm can obtain the best anti-interference performance.
Further, after the interference detection and interference suppression processing, the amplitude (or power) of the output time-domain digital signal output from the time-domain transform and the amplitude (or power) of the input time-domain digital signal are changed, and therefore, it is desirable that the original amplitude (or power) of the input time-domain digital signal can be maintained after the processing, or that the amplitude (or power) of the output time-domain digital signal can be controlled to a stable power value for output. As shown in fig. 10, on the basis of the embodiment shown in fig. 3, a function of detecting the amplitude of the input time-domain digital signal and the amplitude of the output time-domain digital signal, that is, an amplitude detection function, is added, so that a gain control signal is further generated and output by comparing the amplitudes of the input and output signals, and the gain control signal is output to an amplifier that performs gain control on the output time-domain digital signal, so that the amplitude of the output time-domain digital signal can be controlled, and therefore, fig. 10 also adds a function of performing gain control on the output time-domain digital signal as compared with fig. 3, so that the amplitude of the output time-domain digital signal can be gain controlled as needed.
Therefore, the invention discloses an interference detection and suppression method based on windowing processing. The method comprises the steps of frequency domain transformation, interference detection, interference suppression and time domain transformation, wherein overlapping windowing is carried out after the frequency domain transformation, overlap removing processing is carried out after the time domain transformation, in addition, background noise estimation and self-adaption determination of an interference detection threshold can be carried out on a transformed frequency domain digital signal in the interference detection, the accuracy of interference signal detection is improved, and the signal-to-noise ratio loss brought by the windowing processing is reduced. The method has strong universality and is particularly suitable for detecting the interference signals in satellite communication and microwave communication.
The above description is only an embodiment of the present invention, and not intended to limit the scope of the present invention, and all equivalent structural changes made by using the contents of the present specification and the drawings, or applied directly or indirectly to other related technical fields, are included in the scope of the present invention.
Claims (8)
1. An interference detection and suppression method based on windowing processing is characterized by comprising the following steps:
the method comprises the steps of frequency domain transformation, namely transforming a communication signal mixed with an interference signal into an input time domain digital signal through an AD sampler, windowing the input time domain digital signal, and then transforming the input time domain digital signal into a frequency domain digital signal through FFT;
interference detection, namely performing background noise estimation on the frequency domain digital signal and then determining an interference detection threshold;
interference suppression, wherein the component of the frequency domain digital signal which is greater than the interference detection threshold is regarded as an interference signal component, and the interference signal component is suppressed;
and time domain transformation, namely, after the interference signal component in the frequency domain digital signal is suppressed, carrying out IFFT transformation on the frequency domain digital signal to restore the frequency domain digital signal into an output time domain digital signal.
2. The windowing based interference detection and suppression method according to claim 1, wherein said windowing uses a window function comprising windowing said input time domain digital signal using a buttlett window, a hanning window, a Hamming window, a Blackman window, or a Blackman-Karris window.
3. The windowing based interference detection and suppression method according to claim 1, wherein said input time domain digital signal is the sequence x (k):
x(k)＝Ap(k)+n(k)
wherein, p (k) is a PN sequence with equal probability value of +/-l, and the length of the PN sequence is N; n (k) is zero mean and varianceA is the amplitude of the communication signal;
performing windowing on x (k), wherein the window function is w (k), and the sequence after windowing is as follows:
xw(k)＝Ap(k)w(k)+n(k)w(k)。
4. the method of claim 1, wherein the windowing is an overlap windowing corresponding to a de-overlap processing of the time domain digital signal output by the IFFT in the time domain transform.
5. The method of claim 4, wherein the data sequence of the input time-domain digital signal is equally segmented, each data segment has a length N, and when performing the overlap windowing, the overlapping length of the adjacent data segments is Nr, and r (0 ≦ r < 1) is an overlap factor.
6. The method of claim 5, wherein the overlap factor r is 1/2, and the input time-domain digital signal is a sequence x (k) ap (k) + N (k), where p (k) is a PN sequence with equal probability ± l, and the length of p (k) is N; n (k) is zero mean and varianceA is the amplitude of the communication signal;
the output sequence after the overlapping windowing is
Wherein,w (k) is a window function.
7. The method of claim 6, wherein the overlap factor r is 1/2, and the input time-domain digital signal is a sequence x (k) ap (k) + N (k), where p (k) is a PN sequence with equal probability ± l, and the length is N; n (k) is zero mean and varianceA is the amplitude of the communication signal;
the output sequence after the overlapping windowing is
Wherein,w (k) is a window function.
8. The interference detection and suppression method based on windowing processing according to claim 7, wherein the input time domain digital signal is divided into two input sequences, the first input sequence is subjected to windowing processing and then subjected to FFT conversion, interference detection, interference suppression and IFFT conversion to obtain a first output sequence, the second input sequence is subjected to windowing processing after N/2 delay, and then subjected to FFT conversion, interference detection, interference suppression and IFFT conversion to obtain a second output sequence, and the first output sequence and the second output sequence are added to obtain the output time domain digital signal.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811638158.6A CN109639303B (en) | 2018-12-29 | 2018-12-29 | Interference detection and suppression method based on windowing processing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811638158.6A CN109639303B (en) | 2018-12-29 | 2018-12-29 | Interference detection and suppression method based on windowing processing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109639303A true CN109639303A (en) | 2019-04-16 |
| CN109639303B CN109639303B (en) | 2021-06-18 |
Family
ID=66054690
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811638158.6A Active CN109639303B (en) | 2018-12-29 | 2018-12-29 | Interference detection and suppression method based on windowing processing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109639303B (en) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110391820A (en) * | 2019-06-11 | 2019-10-29 | 东南大学 | A kind of Novel Communication method of reseptance for evading co-channel interference based on DFT |
| CN110441600A (en) * | 2019-07-17 | 2019-11-12 | 纳瓦电子（上海）有限公司 | One kind being based on the jamproof adding window method of anti-leak |
| CN110445733A (en) * | 2019-06-27 | 2019-11-12 | 熊军 | Iteration self-adapting channel denoising method and iteration self-adapting channel denoise device |
| CN110933007A (en) * | 2019-08-20 | 2020-03-27 | 熊军 | Device and method for eliminating interference aiming at OFDM broadband signals |
| CN112003632A (en) * | 2020-07-06 | 2020-11-27 | 南京天际砺剑科技有限公司 | Direct sequence spread spectrum communication multi-address interference suppression method |
| CN112649678A (en) * | 2020-12-24 | 2021-04-13 | 广州山锋测控技术有限公司 | Antenna feeder measuring method and device, antenna feeder measuring device and tester |
| CN112995084A (en) * | 2021-02-07 | 2021-06-18 | 比科奇微电子(杭州)有限公司 | Signal processing method and processing device |
| CN113302511A (en) * | 2021-03-31 | 2021-08-24 | 华为技术有限公司 | Interference processing method and device |
| CN113391122A (en) * | 2021-06-09 | 2021-09-14 | 中电科思仪科技股份有限公司 | Method for improving selectivity of frequency spectrum monitoring channel |
| CN113541706A (en) * | 2021-06-23 | 2021-10-22 | 中国电子科技集团公司第三十八研究所 | Narrow-band interference suppression method based on transform domain processing |
| CN114779282A (en) * | 2022-03-30 | 2022-07-22 | 中国科学院国家授时中心 | Continuous wave interference detection method in Loran-C timing and positioning terminal |
| CN116232376A (en) * | 2023-02-20 | 2023-06-06 | 西安宇飞电子技术有限公司 | Frequency domain notch method, device and equipment for resisting narrow-band interference and readable storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101895370A (en) * | 2010-04-01 | 2010-11-24 | 电子科技大学 | Method for detecting interference of OFDM communication system |
| CN105634543A (en) * | 2015-12-30 | 2016-06-01 | 航天恒星科技有限公司 | Narrow-band interference prevention method and system |
-
2018
- 2018-12-29 CN CN201811638158.6A patent/CN109639303B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101895370A (en) * | 2010-04-01 | 2010-11-24 | 电子科技大学 | Method for detecting interference of OFDM communication system |
| CN105634543A (en) * | 2015-12-30 | 2016-06-01 | 航天恒星科技有限公司 | Narrow-band interference prevention method and system |
Non-Patent Citations (5)
| Title |
|---|
| 张春海等: "基于自适应多门限算法的变换域窄带干扰抑制", 《电子与信息学报》 * |
| 曾祥华等: "扩频系统频域窄带干扰抑制算法加窗损耗研究", 《电子与信息学报》 * |
| 李健伟等: "DSSS系统频域干扰抑制算法的信噪比损耗分析", 《通信技术》 * |
| 李平博等: "改进的频域窄带干扰抑制方法", 《空军工程大学学报（自然科学版）》 * |
| 邹宁等: "一种重叠加窗频域抑制窄带干扰算法及研究", 《现代防御技术》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110391820B (en) * | 2019-06-11 | 2021-06-11 | 东南大学 | Communication receiving method for avoiding same frequency interference based on DFT |
| CN110391820A (en) * | 2019-06-11 | 2019-10-29 | 东南大学 | A kind of Novel Communication method of reseptance for evading co-channel interference based on DFT |
| CN110445733A (en) * | 2019-06-27 | 2019-11-12 | 熊军 | Iteration self-adapting channel denoising method and iteration self-adapting channel denoise device |
| CN110441600A (en) * | 2019-07-17 | 2019-11-12 | 纳瓦电子（上海）有限公司 | One kind being based on the jamproof adding window method of anti-leak |
| CN110933007A (en) * | 2019-08-20 | 2020-03-27 | 熊军 | Device and method for eliminating interference aiming at OFDM broadband signals |
| CN112003632A (en) * | 2020-07-06 | 2020-11-27 | 南京天际砺剑科技有限公司 | Direct sequence spread spectrum communication multi-address interference suppression method |
| CN112649678B (en) * | 2020-12-24 | 2024-05-28 | 广州山锋测控技术有限公司 | Antenna feeder measuring method and device, antenna feeder measuring device and tester |
| CN112649678A (en) * | 2020-12-24 | 2021-04-13 | 广州山锋测控技术有限公司 | Antenna feeder measuring method and device, antenna feeder measuring device and tester |
| CN112995084A (en) * | 2021-02-07 | 2021-06-18 | 比科奇微电子(杭州)有限公司 | Signal processing method and processing device |
| CN113302511A (en) * | 2021-03-31 | 2021-08-24 | 华为技术有限公司 | Interference processing method and device |
| CN113302511B (en) * | 2021-03-31 | 2022-05-24 | 华为技术有限公司 | Interference processing method and device |
| CN113391122A (en) * | 2021-06-09 | 2021-09-14 | 中电科思仪科技股份有限公司 | Method for improving selectivity of frequency spectrum monitoring channel |
| CN113541706A (en) * | 2021-06-23 | 2021-10-22 | 中国电子科技集团公司第三十八研究所 | Narrow-band interference suppression method based on transform domain processing |
| CN113541706B (en) * | 2021-06-23 | 2022-04-08 | 中国电子科技集团公司第三十八研究所 | Narrow-band interference suppression method based on transform domain processing |
| CN114779282A (en) * | 2022-03-30 | 2022-07-22 | 中国科学院国家授时中心 | Continuous wave interference detection method in Loran-C timing and positioning terminal |
| CN114779282B (en) * | 2022-03-30 | 2024-05-10 | 中国科学院国家授时中心 | Continuous wave interference detection method in Loran-C timing and positioning terminal |
| CN116232376A (en) * | 2023-02-20 | 2023-06-06 | 西安宇飞电子技术有限公司 | Frequency domain notch method, device and equipment for resisting narrow-band interference and readable storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109639303B (en) | 2021-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109639303B (en) | Interference detection and suppression method based on windowing processing | |
| US5612978A (en) | Method and apparatus for real-time adaptive interference cancellation in dynamic environments | |
| US9118401B1 (en) | Method of adaptive interference mitigation in wide band spectrum | |
| CN102904604B (en) | Narrow-band interference suppression method and device | |
| US7283956B2 (en) | Noise suppression | |
| CN106772457B (en) | A kind of satellite navigation anti-interference method based on empty time-frequency framework | |
| US7324616B2 (en) | Low cost and high performance narrowband interference cancellation system | |
| CN110034834A (en) | A kind of Interference Detection and inhibit optimization method | |
| CN110445733B (en) | Self-adaptive channel denoising method and self-adaptive channel denoising device | |
| CN109474550A (en) | Interference detection suppressor | |
| EP2752848A1 (en) | Method and apparatus for generating a noise reduced audio signal using a microphone array | |
| Akay et al. | Broadband interference excision in spread spectrum communication systems via fractional Fourier transform | |
| CN113824488B (en) | Satellite communication non-malicious interference suppression method based on decision feedback adaptive cancellation | |
| CN105911565B (en) | Method and device for suppressing narrow-band interference | |
| Yu et al. | Cognitive radio wideband spectrum sensing using multitap windowing and power detection with threshold adaptation | |
| US8509341B2 (en) | Communications device including a filter for notching wideband transmit signals and associated methods | |
| CN101154983B (en) | Implementing method for anti-interference technology of satellite communication based on single-time sampling | |
| Huang et al. | Study on interference suppression based on joint fractional Fourier domain and time domain | |
| US9118370B2 (en) | Method and apparatus for impulsive noise mitigation using adaptive blanker based on BPSK modulation system | |
| CN116318246A (en) | Direct sequence spread spectrum anti-interference method based on narrow-band interference suppression | |
| CN113541706B (en) | Narrow-band interference suppression method based on transform domain processing | |
| KR102338762B1 (en) | Apparatus for receiving wireless signal with anti-jamming | |
| Chen et al. | Research on anti-jamming algorithm of BDS signal in WiFi6 environment | |
| Kobayashi et al. | Improved MB Cognitive Radio Spectrum Sensing Using Wavelet Spectrum Filtering | |
| JP2019062385A (en) | Radio communication device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |