WO2014019331A1 - Self-adaptive quick narrowband interference suppression device - Google Patents

Self-adaptive quick narrowband interference suppression device Download PDF

Info

Publication number
WO2014019331A1
WO2014019331A1 PCT/CN2013/000117 CN2013000117W WO2014019331A1 WO 2014019331 A1 WO2014019331 A1 WO 2014019331A1 CN 2013000117 W CN2013000117 W CN 2013000117W WO 2014019331 A1 WO2014019331 A1 WO 2014019331A1
Authority
WO
WIPO (PCT)
Prior art keywords
interference
module
frequency
dft
detection
Prior art date
Application number
PCT/CN2013/000117
Other languages
French (fr)
Chinese (zh)
Inventor
蔺晓龙
翟昆朋
何文涛
金亮良
李尤
徐建华
殷明
Original Assignee
江苏物联网研究发展中心
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 江苏物联网研究发展中心 filed Critical 江苏物联网研究发展中心
Publication of WO2014019331A1 publication Critical patent/WO2014019331A1/en

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details 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/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7097Interference-related aspects
    • H04B1/71Interference-related aspects the interference being narrowband interference
    • H04B1/7102Interference-related aspects the interference being narrowband interference with transform to frequency domain

Definitions

  • the invention belongs to the field of spread spectrum communication, and in particular relates to an adaptive fast narrowband interference suppression device. Background technique
  • GNSS Global Navigation Satellite System
  • GPS Global Positioning System
  • GL0NASS Russia's Gruss system
  • Galileo European Union's Galileo system
  • BD Beidou system
  • the signal intensity of the GNSS received by the ground is around -130dBm, far below the noise power. Interference from other communication equipment leaking into the GNSS signal band can severely degrade the performance of the receiver and even render the receiver inoperable.
  • the most common type of interference is the Continuous Wave Interference, which is a narrowband interference.
  • Frequency domain based narrowband interference suppression technology It is necessary to use Fast Fourier Transform (FFT) and its inverse transform (IFFT) to eliminate narrowband interference in the frequency domain.
  • FFT Fast Fourier Transform
  • IFFT inverse transform
  • Chinese Patent Application No. 201010610524.4, "Frequency-Frequency Receiver Front-End Domain Filtering Anti-Jamming Device and Its Implementation Method” the disadvantage is due to the use of FFT It requires a lot of computing resources, which is complicated to implement and costly.
  • Narrowband interference suppression technology for adaptive FIR filters FIR filter + adaptive algorithm, that is, adaptive FIR filter.
  • FIR filter + adaptive algorithm that is, adaptive FIR filter.
  • the Chinese patent entitled “Blind Adaptive Filtering Device and Its Application for Eliminating Narrowband Interference" with the application number of 201010199697.1 has the disadvantage that the number of filters is larger in order to suppress multiple interferences, and the adaptive algorithm is computationally intensive, and the design is completed. It is fixed afterwards and cannot be automatically configured according to the number of interferences.
  • a narrowband interference suppression technique that cooperates with the cancellation device. Firstly, the characteristics of the frequency, phase and amplitude of the narrowband interference are detected, and then the interference is filtered out by using a parameter-adjustable filter (notch filter); or the reference local signal is subtracted from the input signal to complete the narrowband interference suppression.
  • notch filter parameter-adjustable filter
  • the Chinese patent “Application No. 00803299.8” maintains the performance quality of broadband systems with narrowband interference.
  • the disadvantage is that the implementation is complex, costly, and has no adaptive tracking ability to change the frequency of interference.
  • the object of the present invention is to overcome the deficiencies of the prior art and to provide an adaptive fast narrowband interference suppression apparatus for achieving fast interference detection, suppression and adaptive tracking with less resource consumption.
  • the adaptive fast narrowband interference suppression apparatus includes: one or more digital Fourier transform (DFT) modules, performing digital Fourier transform calculation on the input signal;
  • DFT digital Fourier transform
  • Two interference suppression modules suppress the detected narrowband interference, and the interference suppression module includes a multi-stage cascaded IIR trap;
  • An interference capture control module that controls the DFT module to detect a narrowband interference frequency in the input signal;
  • the input signal is connected to the interference suppression module input, and is output by the interference suppression module; the input signal is also connected to the signal input of the DFT module, the output of the DFT module is connected to the input of the interference capture control module, and the interference capture control module output is connected to the DFT module. Control input and coefficient input of the interference suppression module;
  • the input signal is connected to the interference suppression module input, and is output by the interference suppression module;
  • the input signal is also connected to the signal input of the DFT module, the output of the DFT module is connected to the input of the interference tracking control module, and the interference tracking control module output is connected to the DFT module.
  • the control input and the coefficient input of the interference suppression module, and the coefficient output of the interference suppression module is connected to the input of the interference tracking control module;
  • the detection and acquisition process of the narrowband interference signal is: the input signal received by the narrowband interference is used as the input of the DFT module, and the interference acquisition control module uses the multi-level DFT frequency detection algorithm to control the DFT module to calculate the frequency spectrum of the input signal using the digital Fourier transform; After completion, the interference capture control module detects the narrowband interference frequency in the frequency spectrum, and converts the frequency information into an interference suppression module in which the IIR notch coefficient is placed in the interference suppression module; the frequency of the narrowband interference of the interference suppression module on the input signal Filtering is performed; the tracking process of the narrowband interference signal that has been detected is: After the detection capture is completed and the interference suppression module is started, the interference tracking control module selects a running IIR trap in the interference suppression module to obtain the IIR.
  • the center frequency and bandwidth of the trap, the central frequency and bandwidth are used as configuration parameters to control the DFT module to perform DFT transform calculation on the input signal; the interference tracking control module determines the maximum amplitude of the calculation result: if it is greater than the threshold, the maximum is based on the amplitude Frequency update of the value The coefficient of the IIR notch; if it is less than the threshold, the interference of the frequency has disappeared, and the IIR trap is turned off.
  • the interference capture control module uses a multi-stage DFT frequency detection algorithm as follows: assuming that the frequency range to be detected is S, a total of M levels are used for detection; first, level 1 detection, the range of ⁇ f is performed at a frequency interval of /; After the DFT detection, after the result is obtained, the result of the frequency amplitude exceeding the threshold is recorded in the first table Table-1; then the second level detection is started, and the frequency components recorded in the first table Table-1 are sequentially performed at a frequency interval of / 2 .
  • /i, / 2 ,..., /M respectively represent the frequency interval of the first, second, ..., M level detection, ie the frequency resolution, fu represents the final frequency resolution, determined by the design target;
  • indicates the number of frequencies recorded in each level, which is also the number of narrowband interference suppression defined in the design;
  • the initial value ⁇ 0 is obtained for the series, and then the minimum value is obtained for the equation (5), and the frequency interval of each level detection is determined, ..., / ⁇ , substituted into the equation (5), and the minimum time of the multi-stage DFT frequency detection time is obtained.
  • the interference tracking control module uses the DFT module to the IIR notch Frequency detection is performed in the frequency range of 3 ⁇ 3 ⁇ 4 near the center frequency; if the amplitude of the result is greater than the threshold, the narrowband interference still exists, and the coefficient of the IIR trap is updated according to the frequency at which the amplitude is maximum; If the maximum amplitude of the result is less than the threshold, indicating that the narrowband interference of the frequency has disappeared, the IIR trap is turned off; the interference tracking control module can continuously detect the presence or absence of interference and continuously track the frequency variation of the interference.
  • the design target determines the value, firstly determine J and SPT 3 (3 ⁇ 4 value, according to the design goal and use formula (8) to determine the value of ⁇ ; compare ⁇ and ⁇ size, check whether it meets ⁇ ⁇ 3 / ⁇ , if not satisfied , you need to reduce V, or increase S 3i3 ⁇ 4 , repeat the above steps until you get the required 3dB bandwidth as Sff 3i3 ⁇ 4 and the detected frequency resolution ⁇ .
  • the IIR trap of the present invention uses a 1st order or a 2nd order structure.
  • the present invention can use multiple DFT modules to simultaneously track multiple IIR traps in the interference suppression module, or can use a DFT module to track multiple IIR traps in a time division multiplexed sequence.
  • the results calculated by the DFT module are outputted in order of magnitude values from largest to smallest.
  • Figure 1 is a schematic diagram of a navigation satellite positioning receiver with narrowband interference suppression.
  • Figure 2 (a) is a structural diagram of the acquisition mode of the adaptive fast narrowband interference suppression device.
  • Figure 2 (b) is a structural diagram of the tracking mode of the adaptive fast narrowband interference suppression device.
  • Figure 3 is a block diagram of the DFT module.
  • Figure 4 is a block diagram of the interference suppression module.
  • Figure 5 is a flow chart of interference capture control (multi-level DFT frequency detection).
  • Figure 6 is a flow chart of interference tracking control.
  • Figure 1 shows the structure of a navigation satellite positioning receiver 107 with narrowband interference suppression.
  • the satellite 101 and narrowband interfering signals 102 of a global navigation satellite system (such as a GPS system, a GLONASS system, a Galileo system, a Beidou system, or other global navigation satellite system) are received by the receiver antenna 103 and transmitted to the receiver RF 104 for mixing and filtering.
  • the IF digital signal is generated after the operation such as quantization.
  • a receiver that does not have narrowband interference suppression directly transmits the intermediate frequency digital signal to the receiver baseband 106 for processing. When the interference signal 102 is present, this will cause the receiver baseband 106 to fail to process the satellite signal normally, and the positioning navigation function cannot be realized.
  • the adaptive fast narrowband interference suppression apparatus 105 of the present invention performs frequency detection on the intermediate frequency digital signal output from the radio frequency 104, identifies and suppresses a plurality of narrowband interferences, and simultaneously tracks the frequencies of the interferences, thereby continuously adaptively suppressing narrowband interference. Finally, the output signal after the interference suppression is transmitted to the receiver baseband 106 for measurement and calculation, and the positioning and navigation functions are realized.
  • the receiver 107 is made to have the ability to suppress narrowband interference.
  • the adaptive fast narrowband interference suppression apparatus 105 includes: one or more DFT modules 201, an interference suppression module 202, an interference capture control module 203, and a
  • the interference tracking control module 204 has the following modules connected.
  • the capture phase the input signal is connected to the interference suppression module 202 input, and is output by the interference suppression module 202; the input signal is also connected to the signal input of the DFT module 201, the output of the DFT module 201 is connected to the input of the interference capture control module 203, and the interference capture control module 203 outputs
  • the control input of the DFT module 201 and the coefficient input of the interference suppression module 202 are connected.
  • Tracking phase the input signal is connected to the interference suppression module 202 input, and is output by the interference suppression module 202; the input signal is also connected to the signal input of the DFT module 201, the output of the DFT module 201 is connected to the input of the interference tracking control module 204, and the interference tracking control module 204 outputs
  • the control input of the DFT module 201 and the coefficient input of the interference suppression module 202 are connected, and the coefficient output of the interference suppression module 202 is connected to the input of the interference tracking control module 204.
  • the detection acquisition process for the narrowband interference signal is as follows:
  • the input signal subjected to narrowband interference is input to the DFT module 201, and the interference capture control module 203 controls the DFT module 201 to perform calculation using DFT (Digital Fourier Transform) using a multi-stage DFT frequency detection algorithm.
  • the frequency spectrum of the input signal After the calculation is completed, the interference capture control module 203 detects the narrowband interference frequency in the frequency spectrum, and converts the frequency information into the interference suppression module 202 to place the IIR notch coefficient into the interference suppression module 202.
  • the interference suppression module 202 filters out the frequency of narrowband interference in the input signal.
  • the tracking process for the narrowband interference signal that has been detected is as follows: After the detection capture is completed and the interference suppression module is activated, the interference tracking control module 204 selects an operating IIR trap in the interference suppression module 202 to obtain the IIR trap. The center frequency and bandwidth of the wave, as the center frequency and bandwidth The configuration parameter controls the DFT module 201 to perform a DFT transform calculation on the input signal. The interference tracking control module 204 determines the maximum value of the calculation result: if it is greater than the threshold, the coefficient of the IIR notch is updated according to the frequency of the maximum amplitude; if it is less than the threshold, the interference of the frequency has disappeared, and the IIR trap.
  • the DFT module 201 performs a digital Fourier transform (DFT) on the input signal to calculate a plurality of frequency amplitudes, and then outputs the results in amplitude.
  • the configurable parameter input 301 can control the frequency of the local COS, the SIN carrier 303, the integration time length of the accumulator 304, and the number of DFT calculation results 306.
  • the input signal is multiplied by the local carrier 303 via the mixer 302 and sent to the accumulator 304 for accumulation.
  • the amplitude value of the frequency component is obtained using the amplitude equalizer 305 when the integration is completed, and the value is stored in the register set 306. This is repeated until the number of settlement results input by the input control configuration parameter 301 is completed.
  • the DFT calculation results are arranged as output according to the amplitude from large to small.
  • the DFT module uses Digital Fourier Transform (DFT) technology, whose expression:
  • / s is the sampling frequency and ⁇ is the cumulative time of the L sampling points.
  • the frequency bandwidth to be detected 2 ⁇ 10 ⁇ ⁇ , to reach 100 Hz
  • the detection time ⁇ 1 : ⁇ : 2 is required .
  • Traditional DFT detection time is very long.
  • the interference suppression module 202 implements a narrowband interference suppression function by a plurality of cascaded IIR traps 401. Where -al, bl, -a2, b2 are all configurable coefficients of the IIR notch 401. By inputting the appropriate coefficients, the IIR notch 401 can filter out a specified frequency. Thus, by cascading the IIR filter 401, filtering of a plurality of specified frequencies can be achieved.
  • Each stage of the IIR trap can use a 1st or 2nd order structure.
  • the most common second-order IIR notch is used to illustrate its working principle.
  • is the frequency at which the interference is located;
  • k ⁇ l which is a parameter related to the bandwidth of the filter. The closer the value is to 1, the narrower the bandwidth of the notch, and the wider the bandwidth.
  • Figure 5 is a flow chart of interference capture control (multi-level DFT frequency detection).
  • the interference capture control module uses a multi-level DFT frequency detection algorithm that significantly reduces detection time compared to traditional DFT detection algorithms.
  • the principle is as follows: When each level is detected, only a small part of the narrowband interference frequency component will have a larger amplitude than the threshold, so it is recorded and continues to be detected by the next stage. Most of the frequency range in which no narrowband interference exists is Excluded, reducing overall inspection time.
  • /i , / 2 ,..., /M indicates the frequency interval (frequency resolution) of the 1, 2, ..., M level detection.
  • Fu represents the final frequency resolution and is determined by the design goals.
  • N represents the number of frequencies recorded for each level, and is also the narrowband interference suppression defined in the design.
  • the frequency bandwidth to be detected 2 ⁇ 10 6 ⁇ .
  • the final frequency resolution is 100Hz, but the detection time is only 1% of the traditional detection time (200s), which greatly shortens the detection time.
  • the specific implementation manner is as follows: Assuming that the frequency range to be detected is, a total of M levels are used for detection. First, the first level detection, DFT detection of the range of S by the frequency interval (frequency resolution), after the result is obtained, the result of the frequency amplitude exceeding the threshold is recorded in the first table Table_l; then the second level detection is started, The frequency interval of / 2 is sequentially subjected to DFT detection for the frequency components recorded in Table-1, and after the result is obtained, in the second table Table_2, the frequency amplitude exceeds the threshold result; then the third level detection, method and method Level 2 detection is similar. And so on, until the M-level detection is completed, the final DFT with resolution is obtained. Taking the 2-level DFT frequency detection algorithm as an example, the control flow is as follows:
  • First step 501 Perform level 1 detection.
  • the detection frequency range and the detection frequency resolution are /;, and the DFT module 201 is controlled to operate.
  • the output result 306 of the DFT module 201 is judged: the result whose amplitude is greater than the threshold is recorded in the first table Table-1, which is discarded less than the threshold result.
  • Step 502 Perform a level 2 detection.
  • the frequency components in the first table Table-1 are sequentially used for level 2 detection: detection frequency range /;, detection frequency resolution / 2 , and control DFT module 201 to operate.
  • the output result 306 of the DFT module 201 is judged: the result whose amplitude is greater than the threshold ⁇ 2 is recorded in the second table Table_2, and the result smaller than the threshold ⁇ ⁇ 2 is discarded.
  • Step 503 Determine whether the result of the first 3 ⁇ 4 Table-1 has completed the second level detection, if not, return to the second step 502, and perform the second level detection on the result of not performing the second level detection; otherwise, the description Level 2 testing has been completed and the fourth step continues.
  • Step 505 Convert the detected frequency result into coefficient values -al, bl, -a2, and b2 accepted by the IIR trap in the interference suppression module 202, and sequentially put them into the interference suppression module 202, so that the detected Multiple narrowband interference frequencies are suppressed and the capture control flow is completed.
  • Figure 6 is a flow chart of interference tracking control.
  • the 3dB bandwidth of the filter is S 3DS .
  • Design target decision values first determine the W and 3T3 ⁇ 4 values, and determine the tracking parameters based on the design goals and using equation (8). . Compare ⁇ with 3 3 3 ⁇ 4 size, check if BW 3dB ⁇ 3/ 0 is satisfied, if not, then reduce V, or increase, repeat the above steps until the required 3dB bandwidth is S 3£3 ⁇ 4 and detected. Frequency resolution /. .
  • the interference tracking control module uses DFT The module performs frequency detection on the frequency range of 3rf near the center frequency w where the IIR trap is located. If the maximum amplitude of the result is greater than the threshold, it indicates that the narrowband interference still exists, and the coefficient of the IIR notch is updated according to the frequency at which the maximum value is located; if the maximum amplitude of the result is less than the threshold, the narrowband interference of the frequency has disappeared. Turn off the IIR trap.
  • the interference tracking control module can continuously detect the presence or absence of interference and continuously track the frequency change of the interference.
  • the specific tracking control process is as follows:
  • the first step 601 Check the interference suppression module 202, and select a running IIR trap 401 to perform the ⁇ ⁇ ⁇ .
  • Step 602 Using parameters: detecting frequency range 3 (bandwidth of IIR trap 401), detecting frequency resolution ⁇ , turning on DFT module 201 to perform DFT frequency detection on IIR trap 401.
  • the third step 603 processing the result of the interference detection.
  • These processes include signal processing operations such as fitting, culling, filtering, etc., to make the results more accurate.
  • Step 604 Determine the result: If the maximum amplitude of the DFT detection result is greater than the threshold, the interference still exists, and operation 605 is performed to update the IIR notch coefficient; otherwise, the interference of the frequency has disappeared, and the operation is performed. 606 - Turn off the IIR trap 401.

Abstract

Disclosed is a self-adaptive quick narrowband interference suppression device. The device comprises one or more DFT modules, an interference suppression module, an interference acquisition control module, and an interference tracking control module. The device completes the steps of: firstly, acquiring a narrowband interference signal, wherein the interference acquisition control module uses a multistage DFT frequency detection algorithm to control the DFT module to subject an input signal to narrowband interference frequency detection; after detecting the narrowband interference, initializing and starting up the interference suppression module to filter a plurality of narrowband interferences; and then tracking the narrowband interference signal by using the interference tracking control module to control the DFT module to subject the detected narrowband interference frequency to detection tracking and updating the coefficient of an IIR wave trap in the interference suppression module, thereby conducting self-adaptive tracking and suppression of the narrowband interference. The present invention can quickly detect a plurality of narrowband interferences and continuously track the interference frequency variation, and has the advantages that the implementation cost is low, the detection is quick, and the self-adaptive tracking of the frequency variation of the narrowband interference can be conducted.

Description

自适应快速窄带干扰抑制装置 技术领域  Adaptive fast narrowband interference suppression device
本发明属于扩频通信领域, 具体涉及一种自适应快速窄带干扰抑制装置。 背景技术  The invention belongs to the field of spread spectrum communication, and in particular relates to an adaptive fast narrowband interference suppression device. Background technique
全球导航卫星系统(GNSS)为用户提供定位导航服务, 已广泛应用于消费、 交通、 电力、 国防等领域。 GNSS包括美国的全球定位系统(GPS)、俄罗斯的格 鲁纳斯系统(GL0NASS)、 欧盟的伽利略系统(Galileo)、 中国的北斗系统(BD) 以及其他一些增强系统。  The Global Navigation Satellite System (GNSS) provides users with positioning and navigation services, and has been widely used in the fields of consumption, transportation, power, and defense. GNSS includes the United States' Global Positioning System (GPS), Russia's Gruss system (GL0NASS), the European Union's Galileo system (Galileo), China's Beidou system (BD) and other enhancement systems.
地面接收到的全球导航卫星系统的信号强度在 -130dBm左右, 远在噪声功率 之下。 其他通信设备泄露到全球导航卫星系统信号频段的干扰会严重降低接收 机的性能, 甚至导致接收机无法工作。 其中最常见的干扰类型是连续波干扰 ( Continue Wave Interference ), 它属于窄带干扰。  The signal intensity of the GNSS received by the ground is around -130dBm, far below the noise power. Interference from other communication equipment leaking into the GNSS signal band can severely degrade the performance of the receiver and even render the receiver inoperable. The most common type of interference is the Continuous Wave Interference, which is a narrowband interference.
目前的窄带干扰抑制技术包括有以下这些。  Current narrowband interference suppression techniques include the following.
基于频率域的窄带干扰抑制技术。 需要使用快速傅里叶变换(FFT)及其逆 变换(IFFT),在频率域内消除窄带干扰。如申请号为 200580013652.8的中国专 利《用于窄带干扰消除的接收机》, 申请号为 201010610524.4的中国专利《扩频 接收机前端频域滤波抗干扰装置及其实现方法》,其缺点是由于使用 FFT需要大 量计算资源, 实现复杂, 成本过大。  Frequency domain based narrowband interference suppression technology. It is necessary to use Fast Fourier Transform (FFT) and its inverse transform (IFFT) to eliminate narrowband interference in the frequency domain. For example, Chinese Patent Application No. 200580013652.8, "Receiver for Narrowband Interference Cancellation", Chinese Patent Application No. 201010610524.4, "Frequency-Frequency Receiver Front-End Domain Filtering Anti-Jamming Device and Its Implementation Method", the disadvantage is due to the use of FFT It requires a lot of computing resources, which is complicated to implement and costly.
自适应 FIR滤波器的窄带干扰抑制技术。 FIR滤波器 +自适应算法, 即自适应 FIR滤波器。 如申请号为 201010199697.1的中国专利 《能消除窄带干扰的盲自适 应滤波装置及其应用》, 其缺点是为了抑制多个干扰要较多的滤波器阶数, 自适 应算法计算量大, 设计完成后即固定, 不可根据干扰的个数自动配置。  Narrowband interference suppression technology for adaptive FIR filters. FIR filter + adaptive algorithm, that is, adaptive FIR filter. For example, the Chinese patent entitled "Blind Adaptive Filtering Device and Its Application for Eliminating Narrowband Interference" with the application number of 201010199697.1 has the disadvantage that the number of filters is larger in order to suppress multiple interferences, and the adaptive algorithm is computationally intensive, and the design is completed. It is fixed afterwards and cannot be automatically configured according to the number of interferences.
检测与消除装置配合的窄带干扰抑制技术。首先检测窄带干扰的频率、相位、 幅度等特点, 然后使用参数可调的滤波器(陷波器)滤除干扰; 或者产生参考 本地信号与输入信号相减完成窄带干扰抑制。如申请号为 00803299.8的中国专利 《保持存在窄带干扰的宽带系统的性能质量》。 缺点是实现复杂, 代价大, 无自 适应跟踪干扰频率变化能力。  A narrowband interference suppression technique that cooperates with the cancellation device. Firstly, the characteristics of the frequency, phase and amplitude of the narrowband interference are detected, and then the interference is filtered out by using a parameter-adjustable filter (notch filter); or the reference local signal is subtracted from the input signal to complete the narrowband interference suppression. For example, the Chinese patent “Application No. 00803299.8” maintains the performance quality of broadband systems with narrowband interference. The disadvantage is that the implementation is complex, costly, and has no adaptive tracking ability to change the frequency of interference.
发明内容 Summary of the invention
本发明的目的是克服现有技术的不足, 提供一种自适应快速窄带干扰抑制装 置, 以较少的资源消耗实现快速的干扰检测、 抑制和自适应跟踪。  SUMMARY OF THE INVENTION The object of the present invention is to overcome the deficiencies of the prior art and to provide an adaptive fast narrowband interference suppression apparatus for achieving fast interference detection, suppression and adaptive tracking with less resource consumption.
按照本发明提供的技术方案, 所述的自适应快速窄带干扰抑制装置包括: 一个或多个数字傅里叶变换 (DFT)模块, 对输入信号进行数字傅里叶变换 计算;  According to the technical solution provided by the present invention, the adaptive fast narrowband interference suppression apparatus includes: one or more digital Fourier transform (DFT) modules, performing digital Fourier transform calculation on the input signal;
二个干扰抑制模块, 抑制检测到的窄带干扰, 干扰抑制模块包括多级级联的 IIR陷波器;  Two interference suppression modules suppress the detected narrowband interference, and the interference suppression module includes a multi-stage cascaded IIR trap;
一个干扰捕获控制模块, 控制 DFT模块检测输入信号中的窄带干扰频率; 一个或多个干扰跟踪控制模块,控制 DFT模块和干扰抑制模块对窄带干扰进 行跟踪; An interference capture control module that controls the DFT module to detect a narrowband interference frequency in the input signal; One or more interference tracking control modules, controlling the DFT module and the interference suppression module to track narrowband interference;
在捕获阶段, 输入信号连接干扰抑制模块输入, 由干扰抑制模块输出; 所述 输入信号还连接 DFT模块的信号输入, DFT模块的输出连接干扰捕获控制模块 的输入, 干扰捕获控制模块输出连接 DFT模块的控制输入和干扰抑制模块的系 数输入;  In the capture phase, the input signal is connected to the interference suppression module input, and is output by the interference suppression module; the input signal is also connected to the signal input of the DFT module, the output of the DFT module is connected to the input of the interference capture control module, and the interference capture control module output is connected to the DFT module. Control input and coefficient input of the interference suppression module;
在跟踪阶段, 输入信号连接干扰抑制模块输入, 由干扰抑制模块输出; 所述 输入信号还连接 DFT模块的信号输入, DFT模块的输出连接干扰跟踪控制模块 的输入, 干扰跟踪控制模块输出连接 DFT模块的控制输入和干扰抑制模块的系 数输入, 干扰抑制模块的系数输出连接干扰跟踪控制模块的输入;  In the tracking phase, the input signal is connected to the interference suppression module input, and is output by the interference suppression module; the input signal is also connected to the signal input of the DFT module, the output of the DFT module is connected to the input of the interference tracking control module, and the interference tracking control module output is connected to the DFT module. The control input and the coefficient input of the interference suppression module, and the coefficient output of the interference suppression module is connected to the input of the interference tracking control module;
对窄带干扰信号的检测捕获过程为:受到窄带干扰的输入信号作为 DFT模块 的输入,干扰捕获控制模块使用多级 DFT频率检测算法控制 DFT模块使用数字 傅里叶变换计算输入信号的频率谱; 计算完成后, 干扰捕获控制模块检测频率 谱中的窄带干扰频率,并将该频率信息转化成干扰抑制模块中 IIR陷波器系数置 入干扰抑制模块中; 干扰抑制模块对输入信号中窄带干扰的频率进行滤除; 对已经检测捕获到的窄带干扰信号的跟踪过程为: 在完成了检测捕获并启动 干扰抑制模块后, 干扰跟踪控制模块选取干扰抑制模块中一个运行的 IIR陷波 器,得到该 IIR陷波器的中心频率和带宽,将该中心频率和带宽作为配置参数控 制 DFT模块对输入信号进行 DFT变换计算;干扰跟踪控制模块判断计算结果的 幅度最大值:如果大于门限,则根据该幅度最大值所在的频率更新该 IIR陷波器 的系数; 如果小于门限, 则说明该频率的干扰已经消失, 关闭该 IIR陷波器。  The detection and acquisition process of the narrowband interference signal is: the input signal received by the narrowband interference is used as the input of the DFT module, and the interference acquisition control module uses the multi-level DFT frequency detection algorithm to control the DFT module to calculate the frequency spectrum of the input signal using the digital Fourier transform; After completion, the interference capture control module detects the narrowband interference frequency in the frequency spectrum, and converts the frequency information into an interference suppression module in which the IIR notch coefficient is placed in the interference suppression module; the frequency of the narrowband interference of the interference suppression module on the input signal Filtering is performed; the tracking process of the narrowband interference signal that has been detected is: After the detection capture is completed and the interference suppression module is started, the interference tracking control module selects a running IIR trap in the interference suppression module to obtain the IIR. The center frequency and bandwidth of the trap, the central frequency and bandwidth are used as configuration parameters to control the DFT module to perform DFT transform calculation on the input signal; the interference tracking control module determines the maximum amplitude of the calculation result: if it is greater than the threshold, the maximum is based on the amplitude Frequency update of the value The coefficient of the IIR notch; if it is less than the threshold, the interference of the frequency has disappeared, and the IIR trap is turned off.
所述干扰捕获控制模块使用多级 DFT频率检测算法如下:假设需要检测的频 率范围为 S , 共使用 M个等级进行检测; 首先第 1级检测, 以 /;的频率间隔 对^ f的范围进行 DFT检测,得到结果后,在第一表 Table— 1中记录频率幅度超 过门限的结果; 然后开始第 2级检测, 以 /2的频率间隔对—第一表 Table— 1中记 录的频率分量依次进行 DFT检测, 得到结果后, 在第二表 Table— 2中 录频率 幅度超过门限的结果; 以此类推进行下一级检测; 直到第 M等级—检测完成, 即 得到分辨率为 fM的最终 DFT检测结果。 The interference capture control module uses a multi-stage DFT frequency detection algorithm as follows: assuming that the frequency range to be detected is S, a total of M levels are used for detection; first, level 1 detection, the range of ^f is performed at a frequency interval of /; After the DFT detection, after the result is obtained, the result of the frequency amplitude exceeding the threshold is recorded in the first table Table-1; then the second level detection is started, and the frequency components recorded in the first table Table-1 are sequentially performed at a frequency interval of / 2 . After the DFT test is performed, after the result is obtained, the result of the frequency range exceeding the threshold is recorded in the second table Table-2; the next level detection is performed by analogy; until the Mth level - the detection is completed, the final resolution is f M is obtained. DFT test results.
所述多级 DFT频率检测算法中, 所需要耗费的总时间建模如下:
Figure imgf000004_0001
In the multi-level DFT frequency detection algorithm, the total time required to consume is modeled as follows:
Figure imgf000004_0001
其中: /i,/2,...,/M分别表示第 1,2,...,M等级检测的频率间隔, 即频率分辨率, fu表示最终频率分辨率, ώ设计目标所确定; Ν表示每个等级记录的频率个数, 也是设计中所限定的窄带干扰抑制的个数; Where: /i, / 2 ,..., /M respectively represent the frequency interval of the first, second, ..., M level detection, ie the frequency resolution, fu represents the final frequency resolution, determined by the design target; Ν indicates the number of frequencies recorded in each level, which is also the number of narrowband interference suppression defined in the design;
首先对级数 Μ取初值 Μ0, 然后对式 (5 ) 求极小值, 确定每级检测的频率 间隔 ,…,/^, 代入式(5), 得到多级 DFT频率检测时间最小时间的 Γ值, 然后检査该 Γ值是否满足要求,如果不满足则增加级数 Μ的值,重复上述步骤, 直至得到满足要求的检测时间 Γ和对应的级数 M、 检测的频率间隔 假设 IIR陷波器工作中心频率为 , 其 3dB带宽为 ^3^, 所述干扰跟踪控 制模块使用 DFT模块对 IIR陷波器所处的中心频率 附近 3<¾的频率范围进 行频率检测; 如果结果的幅度最大值大于门限, 说明窄带千扰依然存在, 根据 该幅度最大值所处的频率更新 IIR陷波器的系数;如果结果的幅度最大值小于门 限, 说明该频率的窄带干扰已经消失, 关闭该 IIR陷波器; 千扰跟踪控制模块能 够连续检测干扰存在与否并连续跟踪干扰的频率变化。 First, the initial value Μ0 is obtained for the series, and then the minimum value is obtained for the equation (5), and the frequency interval of each level detection is determined, ..., /^, substituted into the equation (5), and the minimum time of the multi-stage DFT frequency detection time is obtained. Depreciate, then check if the Γ value meets the requirements, if not, increase the value of the Μ ,, repeat the above steps, Until the required detection time Γ and the corresponding number of stages M, the detected frequency interval assumes that the IIR trap operating center frequency is 3 dB bandwidth is ^3^, the interference tracking control module uses the DFT module to the IIR notch Frequency detection is performed in the frequency range of 3<3⁄4 near the center frequency; if the amplitude of the result is greater than the threshold, the narrowband interference still exists, and the coefficient of the IIR trap is updated according to the frequency at which the amplitude is maximum; If the maximum amplitude of the result is less than the threshold, indicating that the narrowband interference of the frequency has disappeared, the IIR trap is turned off; the interference tracking control module can continuously detect the presence or absence of interference and continuously track the frequency variation of the interference.
以 Λ为频率分辨率, 使用 DFT检测 3ί/β的频率范围, 其跟踪频率变化的 能力为-
Figure imgf000005_0001
Using Λ as the frequency resolution, using DFT to detect the frequency range of 3ί/β , its ability to track frequency changes is -
Figure imgf000005_0001
当只有一个 DFT模块, 需要使用时分复用技术, 即一个 DFT模块轮流分配 给多个 IIR陷波器使用时, 则其跟踪频率变化的能力变为:  When there is only one DFT module, it is necessary to use time division multiplexing technology, that is, when one DFT module is allocated to multiple IIR traps in turn, its ability to track frequency changes becomes:
Sweep = ^ Hz I s) (8) 其中: ^为复用的 IIR陷波器个数; Sweep = ^ Hz I s) ( 8 ) where: ^ is the number of multiplexed IIR notches;
设计目标决定 值,首先初步确定 J 和 SPT3(¾值,根据设计目标和利用 公式 (8) 确定 Λ的值; 比较 Λ与 ^的大小, 检査是否满足 ^≥3/ο , 如果不满足, 则需要减小 V, 或者增大 S 3i¾, 重复以上步骤, 直到得到满足 要求的 3dB带宽为 Sff3i¾和检测的频率分辨率 Λ。 The design target determines the value, firstly determine J and SPT 3 (3⁄4 value, according to the design goal and use formula (8) to determine the value of Λ; compare Λ and ^ size, check whether it meets ^ ≥ 3 / ο, if not satisfied , you need to reduce V, or increase S 3i3⁄4 , repeat the above steps until you get the required 3dB bandwidth as Sff 3i3⁄4 and the detected frequency resolution Λ.
本发明所述的 IIR陷波器使用 1阶或 2阶结构。  The IIR trap of the present invention uses a 1st order or a 2nd order structure.
本发明可以使用多个 DFT模块同时跟踪干扰抑制模块中的多个 IIR陷波器, 也可以使用一个 DFT模块时分复用的顺序跟踪多个 IIR陷波器。  The present invention can use multiple DFT modules to simultaneously track multiple IIR traps in the interference suppression module, or can use a DFT module to track multiple IIR traps in a time division multiplexed sequence.
所述 DFT模块计算的结果按幅度值从大到小排列输出。  The results calculated by the DFT module are outputted in order of magnitude values from largest to smallest.
本发明的优点是:  The advantages of the invention are:
1、使用数字傅里叶变换(DFT)模块计算输入信号的频谱, 减少硬件计算资 源消耗;  1. Calculate the spectrum of the input signal using a digital Fourier transform (DFT) module to reduce hardware computational resource consumption;
2、使用多级 DFT频率检测算法进行设计和控制 DFT模块, 大幅度减少检测 时间;  2. Design and control the DFT module using a multi-level DFT frequency detection algorithm to greatly reduce the detection time;
3、 使用级联的 IIR陷波器抑制多个窄带干扰, 减少硬件计算资源消耗; 3. Use cascaded IIR notch to suppress multiple narrowband interferences and reduce hardware computing resource consumption;
4、 使用干扰跟踪算法设计和控制 DFT模块和干扰抑制模块, 实现自适应跟 踪抑制窄带干扰。 4. Use the interference tracking algorithm to design and control the DFT module and the interference suppression module to achieve adaptive tracking to suppress narrowband interference.
附图说明 DRAWINGS
图 1是具有窄带干扰抑制的导航卫星定位接收机示意图。  Figure 1 is a schematic diagram of a navigation satellite positioning receiver with narrowband interference suppression.
图 2 (a) 是自适应快速窄带干扰抑制装置捕获模式的结构示意图。 图 2 (b) 是自适应快速窄带干扰抑制装置跟踪模式的结构示意图。 图 3是 DFT模块结构图。 Figure 2 (a) is a structural diagram of the acquisition mode of the adaptive fast narrowband interference suppression device. Figure 2 (b) is a structural diagram of the tracking mode of the adaptive fast narrowband interference suppression device. Figure 3 is a block diagram of the DFT module.
图 4是干扰抑制模块结构图。  Figure 4 is a block diagram of the interference suppression module.
图 5是干扰捕获控制 (多级 DFT频率检测) 流程图。  Figure 5 is a flow chart of interference capture control (multi-level DFT frequency detection).
图 6是干扰跟踪控制流程图。  Figure 6 is a flow chart of interference tracking control.
具体实施方式 detailed description
下面结合说明书附图和实施例对本发明做进一步描述。  The invention is further described below in conjunction with the drawings and embodiments of the specification.
图 1所示是具有窄带干扰抑制的导航卫星定位接收机 107的结构。全球导航 卫星系统(如 GPS系统、 GLONASS系统、 Galileo系统、 北斗系统或其他全球 导航卫星系统) 的卫星 101与窄带干扰信号 102被接收机天线 103接收并传输 给接收机射频 104进行混频、 滤波、 量化等操作后产生中频数字信号。 不具备 窄带干扰抑制能力的接收机直接将中频数字信号传输给接收机基带 106进行处 理。 在干扰信号 102存在时, 这将会导致接收机基带 106无法正常处理卫星信 号, 无法实现定位导航功能。 本发明的自适应快速窄带干扰抑制装置 105对射 频 104输出的中频数字信号进行频率检测, 识别出多个窄带干扰并进行抑制, 同时跟踪这些干扰的频率, 从而连续自适应抑制窄带干扰。 最后将抑制干扰后 的输出信号传输给接收机基带 106进行测量、 解算, 实现定位与导航功能。 使 接收机 107具有抑制窄带干扰的能力。  Figure 1 shows the structure of a navigation satellite positioning receiver 107 with narrowband interference suppression. The satellite 101 and narrowband interfering signals 102 of a global navigation satellite system (such as a GPS system, a GLONASS system, a Galileo system, a Beidou system, or other global navigation satellite system) are received by the receiver antenna 103 and transmitted to the receiver RF 104 for mixing and filtering. The IF digital signal is generated after the operation such as quantization. A receiver that does not have narrowband interference suppression directly transmits the intermediate frequency digital signal to the receiver baseband 106 for processing. When the interference signal 102 is present, this will cause the receiver baseband 106 to fail to process the satellite signal normally, and the positioning navigation function cannot be realized. The adaptive fast narrowband interference suppression apparatus 105 of the present invention performs frequency detection on the intermediate frequency digital signal output from the radio frequency 104, identifies and suppresses a plurality of narrowband interferences, and simultaneously tracks the frequencies of the interferences, thereby continuously adaptively suppressing narrowband interference. Finally, the output signal after the interference suppression is transmitted to the receiver baseband 106 for measurement and calculation, and the positioning and navigation functions are realized. The receiver 107 is made to have the ability to suppress narrowband interference.
如图 2(a), 图 2(b)所示, 所述的自适应快速窄带干扰抑制装置 105包含: 一 个或多个 DFT模块 201、 一个干扰抑制模块 202、 一个干扰捕获控制模块 203、 一个干扰跟踪控制模块 204, 模块连接方式如下。  As shown in FIG. 2(a) and FIG. 2(b), the adaptive fast narrowband interference suppression apparatus 105 includes: one or more DFT modules 201, an interference suppression module 202, an interference capture control module 203, and a The interference tracking control module 204 has the following modules connected.
捕获阶段:输入信号连接干扰抑制模块 202输入,由干扰抑制模块 202输出; 输入信号还连接 DFT模块 201的信号输入, DFT模块 201的输出连接干扰捕获 控制模块 203的输入, 干扰捕获控制模块 203输出连接 DFT模块 201的控制输 入和干扰抑制模块 202的系数输入。  The capture phase: the input signal is connected to the interference suppression module 202 input, and is output by the interference suppression module 202; the input signal is also connected to the signal input of the DFT module 201, the output of the DFT module 201 is connected to the input of the interference capture control module 203, and the interference capture control module 203 outputs The control input of the DFT module 201 and the coefficient input of the interference suppression module 202 are connected.
跟踪阶段:输入信号连接干扰抑制模块 202输入,由干扰抑制模块 202输出; 输入信号还连接 DFT模块 201的信号输入, DFT模块 201的输出连接干扰跟踪 控制模块 204的输入, 干扰跟踪控制模块 204输出连接 DFT模块 201的控制输 入和干扰抑制模块 202的系数输入, 干扰抑制模块 202的系数输出连接干扰跟 踪控制模块 204的输入。  Tracking phase: the input signal is connected to the interference suppression module 202 input, and is output by the interference suppression module 202; the input signal is also connected to the signal input of the DFT module 201, the output of the DFT module 201 is connected to the input of the interference tracking control module 204, and the interference tracking control module 204 outputs The control input of the DFT module 201 and the coefficient input of the interference suppression module 202 are connected, and the coefficient output of the interference suppression module 202 is connected to the input of the interference tracking control module 204.
对窄带干扰信号的检测捕获过程如下: 受到窄带干扰的输入信号作为 DFT 模块 201的输入, 干扰捕获控制模块 203使用多级 DFT频率检测算法控制 DFT 模块 201使用 DFT (数字傅里叶变换)变换计算输入信号的频率谱。 计算完成 后, 干扰捕获控制模块 203检测频率谱中的窄带干扰频率, 并将该频率信息转 化成干扰抑制模块 202中 IIR陷波器系数置入干扰抑制模块 202中。干扰抑制模 块 202对输入信号中窄带干扰的频率进行滤除。  The detection acquisition process for the narrowband interference signal is as follows: The input signal subjected to narrowband interference is input to the DFT module 201, and the interference capture control module 203 controls the DFT module 201 to perform calculation using DFT (Digital Fourier Transform) using a multi-stage DFT frequency detection algorithm. The frequency spectrum of the input signal. After the calculation is completed, the interference capture control module 203 detects the narrowband interference frequency in the frequency spectrum, and converts the frequency information into the interference suppression module 202 to place the IIR notch coefficient into the interference suppression module 202. The interference suppression module 202 filters out the frequency of narrowband interference in the input signal.
对已经检测捕获到的窄带干扰信号的跟踪过程如下: 在完成了检测捕获并启 动干扰抑制模块后, 干扰跟踪控制模块 204选取干扰抑制模块 202中某个运行 的 IIR陷波器, 得到该 IIR陷波器的中心频率和带宽, 将该中心频率和带宽作为 配置参数控制 DFT模块 201对输入信号进行 DFT变换计算。干扰跟踪控制模块 204判断计算结果的幅度最大值: 如果大于门限, 则根据该幅度最大值所在的频 率更新该 IIR陷波器的系数; 如果小于门限, 则说明该频率的干扰已经消失, 关 闭该 IIR陷波器。 The tracking process for the narrowband interference signal that has been detected is as follows: After the detection capture is completed and the interference suppression module is activated, the interference tracking control module 204 selects an operating IIR trap in the interference suppression module 202 to obtain the IIR trap. The center frequency and bandwidth of the wave, as the center frequency and bandwidth The configuration parameter controls the DFT module 201 to perform a DFT transform calculation on the input signal. The interference tracking control module 204 determines the maximum value of the calculation result: if it is greater than the threshold, the coefficient of the IIR notch is updated according to the frequency of the maximum amplitude; if it is less than the threshold, the interference of the frequency has disappeared, and the IIR trap.
图 3是本发明实施例的 DFT模块结构图。 DFT模块 201对输入信号进行数 字傅里叶变换(DFT)计算多个频率幅度, 然后将结果按幅度大小排列输出。其 中可配置参数输入 301可以控制本地 COS、 SIN载波 303的频率、 累加器 304 的积分时间长度、 DFT计算结果 306结果的个数。 输入信号与本地载波 303经 过混频器 302相乘后送入累加器 304进行累加。 当积分累加完成后使用求幅度 器 305得到频率分量 的幅度值, 将该值存储进寄存器组 306中。 如此反复直 到完成输入控制配置参数 301所输入的结算结果个数。 最后将 DFT计算结果按 照幅度从大到小排列作为输出。  3 is a structural diagram of a DFT module according to an embodiment of the present invention. The DFT module 201 performs a digital Fourier transform (DFT) on the input signal to calculate a plurality of frequency amplitudes, and then outputs the results in amplitude. The configurable parameter input 301 can control the frequency of the local COS, the SIN carrier 303, the integration time length of the accumulator 304, and the number of DFT calculation results 306. The input signal is multiplied by the local carrier 303 via the mixer 302 and sent to the accumulator 304 for accumulation. The amplitude value of the frequency component is obtained using the amplitude equalizer 305 when the integration is completed, and the value is stored in the register set 306. This is repeated until the number of settlement results input by the input control configuration parameter 301 is completed. Finally, the DFT calculation results are arranged as output according to the amplitude from large to small.
DFT模块使用数字傅里叶变换 (DFT) 技术, 其表达式:  The DFT module uses Digital Fourier Transform (DFT) technology, whose expression:
X{ ) = j x{n)e-j<an 其中: L是累加的采样点 DFT结果的频率分辨率:
Figure imgf000007_0001
X{ ) = j x{n)e- j<an where: L is the frequency resolution of the accumulated sample point DFT result:
Figure imgf000007_0001
其中: /s是采样频率, ^是 L个采样点的累加时间。 Where: / s is the sampling frequency and ^ is the cumulative time of the L sampling points.
假设分辨率 /。, 检测带 S 的频率范围所需要的时间为:
Figure imgf000007_0002
Assume resolution /. , the time required to detect the frequency range of the band S is:
Figure imgf000007_0002
以 GPS信号为例, 需要检测的频率带宽 = 2 χ 10ό ίζ, 要达到 100Hz的 Taking the GPS signal as an example, the frequency bandwidth to be detected = 2 χ 10 ό ίζ, to reach 100 Hz
2 χ ΐθ6 2 χ ΐθ 6
频率分辨率, 则需要检测时间^1 : ^^ : 2。。 传统 DFT检测时间很长。 For the frequency resolution, the detection time ^ 1 : ^^ : 2 is required . . Traditional DFT detection time is very long.
图 4是本发明实施例的干扰抑制模块结构图。 干扰抑制模块 202是由多个级 联的 IIR陷波器 401来实现窄带干扰抑制功能的。其中 -al、 bl、 -a2、 b2都是 IIR 陷波器 401的可配置系数。 输入适当的系数, IIR陷波器 401就能实现对一个指 定频率的滤除。如此对级联的 IIR滤波器 401操作,就能实现对多个指定频率的 滤除。  4 is a structural diagram of an interference suppression module according to an embodiment of the present invention. The interference suppression module 202 implements a narrowband interference suppression function by a plurality of cascaded IIR traps 401. Where -al, bl, -a2, b2 are all configurable coefficients of the IIR notch 401. By inputting the appropriate coefficients, the IIR notch 401 can filter out a specified frequency. Thus, by cascading the IIR filter 401, filtering of a plurality of specified frequencies can be achieved.
每一级的 IIR陷波器可以使用 1阶或 2阶的结构。 以最常用的 2阶 IIR陷波 器来说明其工作原理, 2阶 IIR陷波器系统传递函数:  Each stage of the IIR trap can use a 1st or 2nd order structure. The most common second-order IIR notch is used to illustrate its working principle. The second-order IIR notch system transfer function:
1 - a z~ - a2z 1 - az~ - a 2 z
,- 2 cos(iy)、 1 - 2 (4 ) 1 z - z  , - 2 cos(iy), 1 - 2 (4 ) 1 z - z
1 - 2k cos( )z - k z 其中: ^是千扰所在处的频率; k≤l, 是与滤波器带宽有关的参数, 的值越 接近 1, 陷波器的带宽越窄, 反之带宽越宽。 1 - 2k cos( )z - kz Where: ^ is the frequency at which the interference is located; k ≤ l, which is a parameter related to the bandwidth of the filter. The closer the value is to 1, the narrower the bandwidth of the notch, and the wider the bandwidth.
图 5是干扰捕获控制 (多级 DFT频率检测) 流程图。  Figure 5 is a flow chart of interference capture control (multi-level DFT frequency detection).
干扰捕获控制模块使用多级 DFT频率检测算法, 比传统 DFT检测算法大幅 减少检测时间。 其原理为: 每个等级检测时, 仅有小部分窄带干扰所在的频率 分量会有较大的幅度超过门限, 从而被记录下来由下一级继续检测, 大部分不 存在窄带干扰的频率范围被排除掉, 从而减少总体检测时间。  The interference capture control module uses a multi-level DFT frequency detection algorithm that significantly reduces detection time compared to traditional DFT detection algorithms. The principle is as follows: When each level is detected, only a small part of the narrowband interference frequency component will have a larger amplitude than the threshold, so it is recorded and continues to be detected by the next stage. Most of the frequency range in which no narrowband interference exists is Excluded, reducing overall inspection time.
首先需要设计多级 DFT频率检测算法的具体参数。这些参数包括: 级数 M、 每级检测频率范围 /;,/2,/3……、 整体检测时间 Γ。 First, we need to design specific parameters of the multi-level DFT frequency detection algorithm. These parameters include: the number of stages M, the range of detection frequencies per level /;, / 2 , / 3 ..., the overall detection time Γ.
假设使用 Μ等级检测 的频率范围, 所需要耗费的总时间建模如下:
Figure imgf000008_0001
Assuming the frequency range of the Μ level detection is used, the total time required to spend is modeled as follows:
Figure imgf000008_0001
其中: among them:
/i ,/2,...,/M表示第 1,2,...,M等级检测的频率间隔 (频率分辨率)。 /i , / 2 ,..., /M indicates the frequency interval (frequency resolution) of the 1, 2, ..., M level detection.
fu 表示最终频率分辨率, 是由设计目标所确定的。  Fu represents the final frequency resolution and is determined by the design goals.
N表示每个等级记录的频率个数, 也是设计中所限定的窄带干扰抑制的个 N represents the number of frequencies recorded for each level, and is also the narrowband interference suppression defined in the design.
、H- 数。 , H-number.
首先对级数 M取初值 M0, 然后对式 (5 ) 求极小值, 确定每级检测的频率 间隔/;,/2,...,^, 代入表达式 (5 ), 得到多级 DFT频率检测时间最小时间7\ 然后检査 Γ值是否满足要求, 如果不满足则增加级数 M的值, 重复上述步骤, 直至得到满足要求的检测时间 Γ和对应的级数 M、 检测的频率间隔First, take the initial value M0 for the series M, and then find the minimum value for the equation (5), determine the frequency interval /;, / 2 ,..., ^, for each level of detection, and substitute the expression (5) to obtain multiple levels. DFT frequency detection time minimum time 7\ Then check whether the threshold value meets the requirements. If it is not satisfied, increase the value of the level M. Repeat the above steps until the required detection time Γ and the corresponding number of stages M and the frequency of detection are obtained. interval
/ /*2 "··,/Λ/。 / /*2 "··, /Λ/.
同样以 GPS信号为例,需要检测的频率带宽 = 2χ 106 Ηζ,要达到 100Hz 的频率分辨率, 使用多级' DFT频率检测算法, 假设级数 M = 2, 抗干扰个数 N = 8 , /, = 1710Hz, f2 = l00Hz , 检测时间为: Similarly, taking the GPS signal as an example, the frequency bandwidth to be detected = 2 χ 10 6 Ηζ. To achieve a frequency resolution of 100 Hz, a multi-stage 'DFT frequency detection algorithm is used, assuming that the number of stages is M = 2 and the number of anti-interferences is N = 8, /, = 1710Hz, f 2 = l00Hz, detection time is:
T = ^2 + 8 * ΛΛ2 = 2.05* (6) T = ^2 + 8 * ΛΛ2 = 2.05* (6)
100; 100 ;
与传统频率检测相比: 最终频率分辨率都是 100Hz, 然而检测时间仅仅是传 统检测时间 (200s) 的 1%, 大幅度缩短了检测时间。  Compared with the traditional frequency detection: the final frequency resolution is 100Hz, but the detection time is only 1% of the traditional detection time (200s), which greatly shortens the detection time.
确定了多级 DFT频率检测算法的参数后,具体实现方式如下:假设需要检测 的频率范围为 , 共使用 M个等级进行检测。 首先第 1级检测, 以 ;的频率 间隔 (频率分辨率)对 S 的范围进行 DFT检测, 得到结果后, 在第一表 Table_l 中记录频率幅度超过门限的结果; 然后开始第 2级检测, 以/ 2的频率间隔对 一表 Table— 1 中记录的频率分量依次进行 DFT检测, 得到结果后, 在第二表 Table_2中 ^录频率幅度超过门限的结果; 之后是第 3级检测, 方法与第 2级检 测类似。 以此类推, 直到第 M等级检测完成, 即得到分辨率为 的最终 DFT 以 2级 DFT频率检测算法为例, 其控制流程如下: After determining the parameters of the multi-level DFT frequency detection algorithm, the specific implementation manner is as follows: Assuming that the frequency range to be detected is, a total of M levels are used for detection. First, the first level detection, DFT detection of the range of S by the frequency interval (frequency resolution), after the result is obtained, the result of the frequency amplitude exceeding the threshold is recorded in the first table Table_l; then the second level detection is started, The frequency interval of / 2 is sequentially subjected to DFT detection for the frequency components recorded in Table-1, and after the result is obtained, in the second table Table_2, the frequency amplitude exceeds the threshold result; then the third level detection, method and method Level 2 detection is similar. And so on, until the M-level detection is completed, the final DFT with resolution is obtained. Taking the 2-level DFT frequency detection algorithm as an example, the control flow is as follows:
第一步 501 : 进行第 1 .级检测。 检测频率范围 , 检测频率分辨率为 /;, 控制 DFT模块 201运行。计算完成后,对 DFT模块 201的输出结果 306进行判 断: 幅度大于门限 的结果记录在第一表 Table—l中, 小于门限 结果丢 弃。  First step 501: Perform level 1 detection. The detection frequency range and the detection frequency resolution are /;, and the DFT module 201 is controlled to operate. After the calculation is completed, the output result 306 of the DFT module 201 is judged: the result whose amplitude is greater than the threshold is recorded in the first table Table-1, which is discarded less than the threshold result.
第二步 502: 进行第 2级检测。 对第 1级检测第一表 Table— 1中的频率分量 依次使用 2级检测: 检测频率范围 /;, 检测频率分辨率 /2 , 控制 DFT模块 201 运行。计算完成后,对 DFT模块 201的输出结果 306进行判断:幅度大于门限 ΊΉ2 的结果记录在第二表 Table_2中, 小于门限 Γ ί2的结果丢弃。 Step 502: Perform a level 2 detection. For the first-level detection, the frequency components in the first table Table-1 are sequentially used for level 2 detection: detection frequency range /;, detection frequency resolution / 2 , and control DFT module 201 to operate. After the calculation is completed, the output result 306 of the DFT module 201 is judged: the result whose amplitude is greater than the threshold ΊΉ 2 is recorded in the second table Table_2, and the result smaller than the threshold ί ί 2 is discarded.
第三步 503 : 判断第一 ¾ Table— 1中结果是否都已完成第 2级检测, 如果没 有则返回第二步 502, 对未进行第 2级检测的结果进行第 2级检测; 否则说明第 2级检测已全部完成, 继续第四步。  Step 503: Determine whether the result of the first 3⁄4 Table-1 has completed the second level detection, if not, return to the second step 502, and perform the second level detection on the result of not performing the second level detection; otherwise, the description Level 2 testing has been completed and the fourth step continues.
第四步 504: 对第 2级频率检测的结果进行处理, 这些处理包括拟合、剔除、 滤波等信号处理操作, 使结果更加准确。  Step 4: 504: Process the results of the second-level frequency detection, including signal processing operations such as fitting, culling, filtering, etc., to make the result more accurate.
第五步 505 : 将检测到的频率结果转换成干扰抑制模块 202中 IIR陷波器所 接受的系数值 -al、 bl、 -a2、 b2, 并依次置入干扰抑制模块 202中, 使检测到的 多个窄带干扰频率被抑制, 捕获控制流程完成。  Step 505: Convert the detected frequency result into coefficient values -al, bl, -a2, and b2 accepted by the IIR trap in the interference suppression module 202, and sequentially put them into the interference suppression module 202, so that the detected Multiple narrowband interference frequencies are suppressed and the capture control flow is completed.
图 6是干扰跟踪控制流程图。  Figure 6 is a flow chart of interference tracking control.
首先需要设计自适应跟踪算法的参数, 包括: 滤波器的 3dB带宽为 S 3DSFirst, you need to design the parameters of the adaptive tracking algorithm, including: The 3dB bandwidth of the filter is S 3DS ,
DFT检测的频率分辨率/ o。 Frequency resolution of the DFT detection / o.
假设陷波器工作中心频率为《, 其 3dB带宽为 ^^β, 以 Λ为频率分辨率, 使用 DFT检测 ^3ίίΰ的频率范围, 其跟踪频率变化的能力为- Sweep = (Hz/ iAssume that the trap center frequency is ", its 3dB bandwidth is ^^ β , with Λ as the frequency resolution, using DFT to detect the frequency range of ^ 3ίίΰ , its ability to track the frequency change is - Sweep = (Hz/i
Figure imgf000009_0001
Figure imgf000009_0001
当只有一个 DFT模块, 需要使用时分复用技术, 即一个 DFT模块轮流分配 给多个 IIR陷波器使用时, 则其跟踪频率变化的能力变为:  When there is only one DFT module, it is necessary to use time division multiplexing technology, that is, when one DFT module is allocated to multiple IIR traps in turn, its ability to track frequency changes becomes:
Sweep = ^{Hz I s) ( 8) 其中: V为复用的 IIR陷波器个数。 Sweep = ^{Hz I s) ( 8 ) where: V is the number of multiplexed IIR notches.
设计目标决定 值,首先初步确定 W和 3T¾值,根据设计目标和利用 公式 (8 ) 确定跟踪的参数/。。 比较 Λ与 3ί¾的大小, 检查是否满足 BW3dB≥ 3/0 , 如果不满足, 则需要减小 V, 或者增大 , 重复以上步骤, 直到得到满足要求的 3dB带宽为 S 3£¾和检测的频率分辨率 /。。 Design target decision values, first determine the W and 3T3⁄4 values, and determine the tracking parameters based on the design goals and using equation (8). . Compare Λ with 3 3 3⁄4 size, check if BW 3dB ≥ 3/ 0 is satisfied, if not, then reduce V, or increase, repeat the above steps until the required 3dB bandwidth is S 3£3⁄4 and detected. Frequency resolution /. .
确定了干扰跟踪参数后, 具体实现方式如下: 干扰跟踪控制模块使用 DFT 模块对 IIR陷波器所处的中心频率 w附近 3rf 的频率范围进行频率检测。如果 结果的幅度最大值大于门限, 说明窄带干扰依然存在, 根据该幅度最大值所处 的频率更新 IIR陷波器的系数;如果结果的幅度最大值小于门限,说明该频率的 窄带干扰已经消失,关闭该 IIR陷波器。干扰跟踪控制模块能够连续检测干扰存 在与否并连续跟踪干扰的频率变化。 After the interference tracking parameters are determined, the specific implementation is as follows: The interference tracking control module uses DFT The module performs frequency detection on the frequency range of 3rf near the center frequency w where the IIR trap is located. If the maximum amplitude of the result is greater than the threshold, it indicates that the narrowband interference still exists, and the coefficient of the IIR notch is updated according to the frequency at which the maximum value is located; if the maximum amplitude of the result is less than the threshold, the narrowband interference of the frequency has disappeared. Turn off the IIR trap. The interference tracking control module can continuously detect the presence or absence of interference and continuously track the frequency change of the interference.
具体跟踪控制流程如下:  The specific tracking control process is as follows:
第一步 601 :检查干扰抑制模块 202,选择一个运行中的 IIR陷波器 401进行 ί艮足宗。  The first step 601: Check the interference suppression module 202, and select a running IIR trap 401 to perform the 艮 艮 。.
第二步 602: 使用参数: 检测频率范围 3 (IIR陷波器 401的带宽), 检 测频率分辨率 Λ, 开启 DFT模块 201对 IIR陷波器 401进行 DFT频率检测。 Step 602: Using parameters: detecting frequency range 3 (bandwidth of IIR trap 401), detecting frequency resolution Λ, turning on DFT module 201 to perform DFT frequency detection on IIR trap 401.
第三步 603 : 处理干扰检测的结果。 这些处理包括拟合、 剔除、 滤波等信号 处理操作, 使结果更加准确。  The third step 603: processing the result of the interference detection. These processes include signal processing operations such as fitting, culling, filtering, etc., to make the results more accurate.
第四步 604: 对结果进行判断: 如果 DFT检测结果的幅度最大值大于门限, 则说明干扰依然存在, 运行操作 605——更新 IIR陷波器系数; 否则, 该频率的 干扰已经消失, 运行操作 606——关闭 IIR陷波器 401。  Step 604: Determine the result: If the maximum amplitude of the DFT detection result is greater than the threshold, the interference still exists, and operation 605 is performed to update the IIR notch coefficient; otherwise, the interference of the frequency has disappeared, and the operation is performed. 606 - Turn off the IIR trap 401.

Claims

1、 自适应快速窄带干扰抑制装置, 其特征在于, 包括: 1. Adaptive fast narrowband interference suppression device, characterized by: including:
一个或多个 DFT模块(201 ), 对输入信号进行数字傅里叶变换计算; 一个干扰抑制模块 (202), 抑制检测到的窄带干扰, 干扰抑制模块(202) 包括多级级联的 IIR陷波器(401 ); One or more DFT modules (201) perform digital Fourier transform calculations on the input signal; an interference suppression module (202) suppresses the detected narrowband interference, and the interference suppression module (202) includes a multi-stage cascaded IIR trap. wave device(401);
一个干扰捕获控制模块(203), 控制 DFT模块 (201 )检测输入信号中的窄 带干扰频率; An interference capture control module (203) controls the DFT module (201) to detect the narrow-band interference frequency in the input signal;
一个或多个干扰跟踪控制模块(204), 控制 DFT模块(201 )和干扰抑制模 块(202)对窄带干扰进行跟踪; One or more interference tracking control modules (204) control the DFT module (201) and the interference suppression module (202) to track narrowband interference;
在捕获阶段,输入信号连接干扰抑制模块(202)输入,由干扰抑制模块(202) 输出; 所述输入信号还连接 DFT模块(201 ) 的信号输入, DFT模块(201 ) 的 输出连接干扰捕获控制模块(203 ) 的输入, 干扰捕获控制模块(203)输出连 接 DFT模块(201 ) 的控制输入和干扰抑制模块(202) 的系数输入; In the capture phase, the input signal is connected to the input of the interference suppression module (202) and output by the interference suppression module (202); the input signal is also connected to the signal input of the DFT module (201), and the output of the DFT module (201) is connected to the interference capture control The input of the module (203), the output of the interference capture control module (203) is connected to the control input of the DFT module (201) and the coefficient input of the interference suppression module (202);
在跟踪阶段,输入信号连接干扰抑制模块(202)输入,由干扰抑制模块(202) 输出; 所述输入信号还连接 DFT模块(201 ) 的信号输入, DFT模块(201 ) 的 输出连接干扰跟踪控制模块(204) 的输入, 干扰跟踪控制模块(204)输出连 接 DFT模块(201 ) 的控制输入和干扰抑制模块(202) 的系数输入, 干扰抑制 模块(202) 的系数输出连接干扰跟踪控制模块 (204) 的输入; In the tracking phase, the input signal is connected to the input of the interference suppression module (202) and output by the interference suppression module (202); the input signal is also connected to the signal input of the DFT module (201), and the output of the DFT module (201) is connected to the interference tracking control The input of the module (204), the output of the interference tracking control module (204) is connected to the control input of the DFT module (201) and the coefficient input of the interference suppression module (202), and the coefficient output of the interference suppression module (202) is connected to the interference tracking control module (202). 204) input;
对窄带干扰信号的检测捕获过程为:受到窄带干扰的输入信号作为 DFT模块 (201 )的输入,干扰捕获控制模块(203 )使用多级 DFT频率检测算法控制 DFT 模块(201 )使用数字傅里叶变换计算输入信号的频率谱; 计算完成后, 干扰捕 获控制模块(203 )检测频率谱中的窄带干扰频率, 并将该频率信息转化成干扰 抑制模块(202) 中 IIR陷波器(401 )系数置入干扰抑制模块(202) 中; 干扰 抑制模块(202)对输入信号中窄带干扰的频率进行滤除; The detection and capture process of narrowband interference signals is: the input signal subject to narrowband interference is used as the input of the DFT module (201), and the interference capture control module (203) uses a multi-level DFT frequency detection algorithm to control the DFT module (201) using digital Fourier Transform and calculate the frequency spectrum of the input signal; after the calculation is completed, the interference capture control module (203) detects the narrow-band interference frequency in the frequency spectrum, and converts the frequency information into the IIR notch (401) coefficient in the interference suppression module (202) Placed in the interference suppression module (202); The interference suppression module (202) filters out the frequency of narrowband interference in the input signal;
对已经检测捕获到的窄带干扰信号的跟踪过程为: 在完成了检测捕获并启动 干扰抑制模块后, 干扰跟踪控制模块 (204)选取干扰抑制模块 (202) 中一个 运行的 IIR陷波器 (401 ), 得到该 IIR陷波器 (401 ) 的中心频率和带宽, 将该 中心频率和带宽作为配置参数控制 DFT模块(201 )对输入信号进行 DFT变换 计算; 干扰跟踪控制模块(204)判断计算结果的幅度最大值: 如果大于门限, 则根据该幅度最大值所在的频率更新该 IIR陷波器(401 ) 的系数; 如果小于门 限, 则说明该频率的干扰已经消失, 关闭该 IIR陷波器 (401 )。 The tracking process of the narrowband interference signal that has been detected and captured is: After completing the detection and capture and starting the interference suppression module, the interference tracking control module (204) selects an operating IIR trap (401) in the interference suppression module (202) ), obtain the center frequency and bandwidth of the IIR trap (401), and use the center frequency and bandwidth as configuration parameters to control the DFT module (201) to perform DFT transformation calculations on the input signal; the interference tracking control module (204) judges the calculation results The maximum amplitude of 401).
2、 根据权利要求 1 所述自适应快速窄带干扰抑制装置, 其特征在于, 所述 干扰捕获控制模块(203)使用多级 DFT频率检测算法如下:假设需要检测的频 率范围为^^, 共使用 M个等级进行检测; 首先第 1级检测, 以/;的频率间隔 对 S 的范围进行 DFT检测,得到结果后,在第一表 Table一 1中记录频率幅度超 过门限的结果; 然后开始第 2级检测, 以 /2的频率间隔对一第一表 Table— 1中记 录的频率分量依次进行 DFT检测, 得到结果后, 在第二表 Table_2中 录频率 幅度超过门限的结果; 以此类推进行下一级检测; 直到第 M等级检测完成, 即 得到分辨率为 fM的最终 DFT检测结果。 2. The adaptive fast narrowband interference suppression device according to claim 1, characterized in that the interference capture control module (203) uses a multi-level DFT frequency detection algorithm as follows: assuming that the frequency range to be detected is ^^, a total of M levels are tested; first, the first level is tested, and the range of S is tested by DFT at a frequency interval of Level detection, perform DFT detection on the frequency components recorded in the first table Table-1 in sequence at a / 2 frequency interval. After obtaining the results, record the frequency in the second table Table_2 The amplitude exceeds the threshold; and so on for the next level of detection; until the M-th level detection is completed, the final DFT detection result with a resolution of f M is obtained.
3、 根据权利要求 2所述自适应快速窄带干扰抑制装置, 其特征在于, 所述 多级 DFT频率检测算法中, 所需要耗费的总时间建模如下-
Figure imgf000012_0001
3. The adaptive fast narrowband interference suppression device according to claim 2, characterized in that, in the multi-level DFT frequency detection algorithm, the total time required is modeled as follows -
Figure imgf000012_0001
其中: /1 ,/2,...,/Α分别表示第 1,2,...,M等级检测的频率间隔, 即频率分辨率,Among them: /1,/2,...,/Α respectively represent the frequency intervals of the 1st, 2nd ,..., M level detection, that is, the frequency resolution,
/M表示最终频率分辨率,由设计目标所确定; 表示每个等级记录的频率个数, 也是设计中所限定的窄带干扰抑制的个数; /M represents the final frequency resolution, which is determined by the design goal; represents the number of frequencies recorded at each level, which is also the number of narrowband interference suppression limited in the design;
首先对级数 M取初值 M0 , 然后对式 (5 ) 求极小值, 确定每级检测的频率 间隔 , ,…,^, 代入式 (5), 得到多级 DFT频率检测时间最小时间的 Γ值, 然后检査该 Γ值是否满足要求,如果不满足则增加级数 M的值, 重复上述步骤, 直至得到满足要求的检测时间 Γ和对应的级数 M、 检测的频率间隔 i,fi,""f i ° First, take the initial value M0 for the series M, and then find the minimum value of equation (5) to determine the frequency interval of each level of detection, ,..., ^. Substitute into equation (5) to get the minimum time of multi-level DFT frequency detection time Γ value, and then check whether the Γ value meets the requirements. If not, increase the value of the series M, and repeat the above steps until the detection time Γ that meets the requirements, the corresponding series M, and the detection frequency interval i, fi are obtained. ,""f i °
4、 根据权利要求 1 所述自适应快速窄带干扰抑制装置, 其特征在于, 假设 IIR陷波器 (401 ) 工作中心频率为 ^, 其 3dB带宽为 3ί/β, 所述干扰跟踪控 制模块(204)使用 DFT模块(201 )对 IIR陷波器(401 )所处的中心频率《附 近 S 3rf 的频率范围进行频率检测; 如果结果的幅度最大值大于门限, 说明窄 带干扰依然存在, 根据该幅度最大值所处的频率更新 IIR陷波器(401 )的系数; 如果结果的幅度最大值小于门限, 说明该频率的窄带干扰已经消失, 关闭该 IIR 陷波器(401 ); 干扰跟踪控制模块(204) 能够连续检测干扰存在与否并连续跟 踪干扰的频率变化。 4. The adaptive fast narrowband interference suppression device according to claim 1, characterized in that, assuming that the IIR trap (401) has a working center frequency of ^ and its 3dB bandwidth is 3ί/β , the interference tracking control module (204 ) Use the DFT module (201) to perform frequency detection in the frequency range S 3rf near the center frequency of the IIR trap (401); if the maximum amplitude of the result is greater than the threshold, it means that narrowband interference still exists, according to the maximum amplitude Update the coefficient of the IIR notch filter (401) at the frequency where the value is located; if the maximum amplitude of the result is less than the threshold, it means that the narrowband interference at this frequency has disappeared, and the IIR notch filter (401) is turned off; the interference tracking control module (204 ) can continuously detect the presence or absence of interference and continuously track the frequency changes of interference.
5、根据权利要求 4所述自适应快速窄带干扰抑制装置,其特征在于,以 Λ为 频率分辨率, 使用 DFT检测 3ί¾的频率范围, 其跟踪频率变化的能力为- 5. The adaptive fast narrowband interference suppression device according to claim 4, characterized in that, taking Λ as the frequency resolution, using DFT to detect a frequency range of 3ί¾ , its ability to track frequency changes is -
Sweep = ( 7 )
Figure imgf000012_0002
Sweep = (7)
Figure imgf000012_0002
当只有一个 DFT模块 (201 ),需要使用时分复用技术,即一个 DFT模块 (201 ) 轮流分配给多个 IIR陷波器 (401 ) 使用时, 则其跟踪频率变化的能力变为- When there is only one DFT module (201) and time division multiplexing technology needs to be used, that is, when one DFT module (201) is assigned to multiple IIR notches (401) in turn, its ability to track frequency changes becomes -
Sweep = ^ (Hz / s) ( 8) 其中: 为复用的 IIR陷波器 (401 ) 个数; Sweep = ^ (Hz / s) ( 8 ) where: is the number of multiplexed IIR traps (401);
设计目标决定^ eep值,首先初步确定 N和 β 3ί¾值,根据设计目标和利用 公式 (8) 确定 /。的值; 比较 /。与^^ ^的大小, 检査是否满足 3ί¾≥3/。, 如果不满足, 则需要减小 V, 或者增大 BW3dB, 重复以上步骤, 直到得到满足 要求的 3dB带宽为 3 和检测的频率分辨率 Λ。 The design goal determines the ^ eep value. First, initially determine the N and β 3ί¾ values, and determine / according to the design goal and the use of formula (8). The value of ; compares /. With the size of ^^ ^, check whether it satisfies 3ί¾ ≥ 3/. , if not satisfied, you need to reduce V, or increase BW 3dB , repeat the above steps until satisfied The required 3dB bandwidth is 3 and the detection frequency resolution Λ.
6、 根据权利要求 1 所述自适应快速窄带干扰抑制装置, 其特征在于, 所述 IIR陷波器(401 ) 使用 1阶或 2阶结构。 6. The adaptive fast narrowband interference suppression device according to claim 1, characterized in that the IIR trap (401) uses a 1st-order or 2nd-order structure.
7、 根据权利要求 1 所述自适应快速窄带干扰抑制装置, 其特征在于, 使用 多个 DFT模块(201 )同时跟踪干扰抑制模块(202 )中的多个 IIR陷波器(401 ), 或使用一个 DFT模块 (201 ) 时分复用的顺序跟踪多个 IIR陷波器(401 )。 7. The adaptive fast narrowband interference suppression device according to claim 1, characterized by using multiple DFT modules (201) to simultaneously track multiple IIR traps (401) in the interference suppression module (202), or using A DFT module (201) time-division multiplexes sequential tracking of multiple IIR notches (401).
8、 根据权利要求 1 所述自适应快速窄带干扰抑制装置, 其特征在于, 所述 DFT模块(201 ) 计算的结果按幅度值从大到小排列输出。 8. The adaptive fast narrowband interference suppression device according to claim 1, characterized in that the results calculated by the DFT module (201) are arranged and output in order of amplitude values from large to small.
PCT/CN2013/000117 2012-08-02 2013-02-04 Self-adaptive quick narrowband interference suppression device WO2014019331A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201210273228.9 2012-08-02
CN201210273228.9A CN102769478B (en) 2012-08-02 2012-08-02 Self-adaption rapid narrow-band interference suppression device

Publications (1)

Publication Number Publication Date
WO2014019331A1 true WO2014019331A1 (en) 2014-02-06

Family

ID=47096744

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2013/000117 WO2014019331A1 (en) 2012-08-02 2013-02-04 Self-adaptive quick narrowband interference suppression device

Country Status (2)

Country Link
CN (1) CN102769478B (en)
WO (1) WO2014019331A1 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2933928A1 (en) * 2014-04-16 2015-10-21 Digi International Inc. Low complexity narrowband interference suppression
EP3671283A1 (en) * 2018-12-21 2020-06-24 Accord Ideation Private Limited Tones processing system in a global navigation satellite system receiver
US11671133B2 (en) 2020-10-16 2023-06-06 Deere & Company Adaptive narrowband and wideband interference rejection for satellite navigation receiver
US11742883B2 (en) 2020-10-16 2023-08-29 Deere & Company Adaptive narrowband interference rejection for satellite navigation receiver
US11750274B2 (en) 2020-10-16 2023-09-05 Deere & Company Adaptive narrowband interference rejection for satellite navigation receiver

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102769478B (en) * 2012-08-02 2014-05-21 江苏物联网研究发展中心 Self-adaption rapid narrow-band interference suppression device
CN104215976A (en) * 2014-09-12 2014-12-17 中国电子科技集团公司第二十研究所 Narrow-band interference resisting system based on cascaded wave traps
CN104319988B (en) * 2014-09-28 2017-09-22 安徽中科医药成像技术科技有限公司 The system for eliminating the method for interfering frequency and eliminating interfering frequency
CN105549035B (en) * 2015-12-22 2018-05-29 武汉梦芯科技有限公司 A kind of baseband signal frequency domain narrowband Interference Detection cancellation element and method
CN106501819B (en) * 2016-11-16 2018-11-16 西安电子科技大学 Interfere implementation method in the anti-narrowband of satellite receiver based on FPGA
CN106656226A (en) * 2016-12-06 2017-05-10 大唐半导体设计有限公司 Narrowband interference processing method and device
CN107290758B (en) * 2017-05-02 2020-07-07 南京航空航天大学 GNSS interference signal multi-stage identification and detection system and method
CN110493860B (en) * 2018-05-14 2022-09-16 海能达通信股份有限公司 Multimode communication method, mobile terminal and storage medium
CN111443362B (en) * 2018-12-29 2022-04-01 上海司南卫星导航技术股份有限公司 Receiving device, terminal device and computer readable storage medium
CN110244127B (en) * 2019-07-30 2021-03-12 广东电网有限责任公司 Device and method for testing frequency sweeping impedance of transformer
CN112589794A (en) * 2020-12-02 2021-04-02 法奥意威(苏州)机器人系统有限公司 Method for suppressing vibration of robot
CN114916071A (en) * 2021-02-10 2022-08-16 Oppo广东移动通信有限公司 Near field communication module, electronic equipment and near field communication method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102291160A (en) * 2011-05-18 2011-12-21 西安电子科技大学 TDM (Time Division Multiplexing)-based high-order adaptive wave trap of parallel structure and adaptive wave trap method
CN102571139A (en) * 2011-12-31 2012-07-11 北京华信创通科技有限公司 Adaptive anti-interference system and method
CN102590829A (en) * 2012-03-14 2012-07-18 西安电子科技大学 Complete self-adaptive notch filter for satellite navigation system and notch filtering method of same
CN102611672A (en) * 2012-04-23 2012-07-25 哈尔滨工业大学 Narrowband interference suppression method based on weighted score Fourier transform domain
CN102769478A (en) * 2012-08-02 2012-11-07 江苏物联网研究发展中心 Self-adaption rapid narrow-band interference suppression device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101571585B (en) * 2009-03-18 2011-08-24 杭州中科微电子有限公司 Method for eliminating interference signal of GPS receiver and system thereof

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102291160A (en) * 2011-05-18 2011-12-21 西安电子科技大学 TDM (Time Division Multiplexing)-based high-order adaptive wave trap of parallel structure and adaptive wave trap method
CN102571139A (en) * 2011-12-31 2012-07-11 北京华信创通科技有限公司 Adaptive anti-interference system and method
CN102590829A (en) * 2012-03-14 2012-07-18 西安电子科技大学 Complete self-adaptive notch filter for satellite navigation system and notch filtering method of same
CN102611672A (en) * 2012-04-23 2012-07-25 哈尔滨工业大学 Narrowband interference suppression method based on weighted score Fourier transform domain
CN102769478A (en) * 2012-08-02 2012-11-07 江苏物联网研究发展中心 Self-adaption rapid narrow-band interference suppression device

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2933928A1 (en) * 2014-04-16 2015-10-21 Digi International Inc. Low complexity narrowband interference suppression
US9276797B2 (en) 2014-04-16 2016-03-01 Digi International Inc. Low complexity narrowband interference suppression
EP3671283A1 (en) * 2018-12-21 2020-06-24 Accord Ideation Private Limited Tones processing system in a global navigation satellite system receiver
US11671133B2 (en) 2020-10-16 2023-06-06 Deere & Company Adaptive narrowband and wideband interference rejection for satellite navigation receiver
US11742883B2 (en) 2020-10-16 2023-08-29 Deere & Company Adaptive narrowband interference rejection for satellite navigation receiver
US11750274B2 (en) 2020-10-16 2023-09-05 Deere & Company Adaptive narrowband interference rejection for satellite navigation receiver
US11764862B2 (en) 2020-10-16 2023-09-19 Deere & Company Adaptive narrowband interference rejection for satellite navigation receiver

Also Published As

Publication number Publication date
CN102769478A (en) 2012-11-07
CN102769478B (en) 2014-05-21

Similar Documents

Publication Publication Date Title
WO2014019331A1 (en) Self-adaptive quick narrowband interference suppression device
Borio et al. Two-pole and multi-pole notch filters: A computationally effective solution for GNSS interference detection and mitigation
US9065686B2 (en) Spur detection, cancellation and tracking in a wireless signal receiver
CN102904604B (en) Narrow-band interference suppression method and device
WO2006092093A1 (en) Method and device for suppressing narrowband interference
RU2513028C2 (en) Device for suppressing narrow-band interference in satellite navigation receiver
RU2012117904A (en) FILTER DEVICE
CN105429918B (en) A kind of quick capturing method for deep space TT&C signal
CN109741762B (en) Voice activity detection method and device and computer readable storage medium
CN106797517B (en) Multi-ear MMSE analysis techniques for cleaning audio signals
RU2007136612A (en) METHOD AND DEVICE FOR REMOVING NARROW-BAND RADIO INTERFERENCE IN A BROADBAND SYSTEM
Chen et al. Chirp-like jamming mitigation for GPS receivers using wavelet-packet-transform-assisted adaptive filters
KR100996430B1 (en) Apparatus and method for collecting frequency hopping signal using direction of arrival
RU2351005C1 (en) Method of evaluating signal parameters and device to this end (versions)
Yaghoobi et al. A low-complexity sub-Nyquist sampling system for wideband radar ESM receivers
CN107728167A (en) A kind of anti-drying method of GPS satellite navigation receiver analyzed based on cyclostationary characteristic with circulating filtering
US10880029B2 (en) System and method for improved data decoding, tracking, and other receiver functions in the presence of interference
KR101501483B1 (en) Method of blind demodulation using hopping phase estimator in hopped FM/BFSK Systems
Qin et al. Effects of interference mitigation methods on scintillation detection
RU2395158C1 (en) Digital signal filtration method
FR2974443A1 (en) METHOD AND SYSTEM FOR REDUCING NOISE
Abdizadeh GNSS signal acquisition in the presence of narrowband interference
Lu et al. Impulsive noise excision and performance analysis
JP2010028530A (en) Signal separation device and signal separation method
JPH11118906A (en) Receiver

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13825994

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 13825994

Country of ref document: EP

Kind code of ref document: A1