[go: up one dir, main page]

CN102098258A - Method for removing narrow-band interference and self-adapting filter - Google Patents

Method for removing narrow-band interference and self-adapting filter Download PDF

Info

Publication number
CN102098258A
CN102098258A CN2009102254497A CN200910225449A CN102098258A CN 102098258 A CN102098258 A CN 102098258A CN 2009102254497 A CN2009102254497 A CN 2009102254497A CN 200910225449 A CN200910225449 A CN 200910225449A CN 102098258 A CN102098258 A CN 102098258A
Authority
CN
China
Prior art keywords
current sampling
signal
weight
sampling moment
reference 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
Application number
CN2009102254497A
Other languages
Chinese (zh)
Other versions
CN102098258B (en
Inventor
王守义
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
ZTE Corp
Original Assignee
ZTE Corp
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 ZTE Corp filed Critical ZTE Corp
Priority to CN 200910225449 priority Critical patent/CN102098258B/en
Priority to PCT/CN2010/072409 priority patent/WO2010148785A1/en
Publication of CN102098258A publication Critical patent/CN102098258A/en
Application granted granted Critical
Publication of CN102098258B publication Critical patent/CN102098258B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H21/00Adaptive networks
    • H03H21/0012Digital adaptive filters

Landscapes

  • Noise Elimination (AREA)
  • Soundproofing, Sound Blocking, And Sound Damping (AREA)

Abstract

本发明实施例涉及信号处理领域,特别涉及一种去除窄带干扰的方法和自适应滤波器,用以解决现有技术中存在的通过固定系数滤波器除去干扰会引入处理时延,实时性比较差的问题。本发明实施例的方法包括:自适应滤波器确定当前采样时刻对应的权值;所述自适应滤波器根据确定的权值和在当前采样时刻采样的参考信号,确定当前采样时刻对应的抵消信号,其中所述参考信号是根据窄带干扰的中心频率生成的;所述自适应滤波器根据确定的当前采样时刻对应的抵消信号,对在当前采样时刻采样的基带信号进行滤波,得到除去窄带干扰后的基带信号。采用本发明实施例能够减少处理时延,提高实时性,改善干扰抵消效果。

Figure 200910225449

The embodiment of the present invention relates to the field of signal processing, in particular to a method for removing narrow-band interference and an adaptive filter, which are used to solve the problem in the prior art that removing interference through a fixed-coefficient filter will introduce processing delay and poor real-time performance The problem. The method in the embodiment of the present invention includes: the adaptive filter determines the weight value corresponding to the current sampling moment; the adaptive filter determines the cancellation signal corresponding to the current sampling moment according to the determined weight value and the reference signal sampled at the current sampling moment , wherein the reference signal is generated according to the center frequency of the narrowband interference; the adaptive filter filters the baseband signal sampled at the current sampling moment according to the determined canceling signal corresponding to the current sampling moment, and obtains after removing the narrowband interference baseband signal. By adopting the embodiment of the present invention, the processing delay can be reduced, the real-time performance can be improved, and the effect of interference cancellation can be improved.

Figure 200910225449

Description

一种去除窄带干扰的方法和自适应滤波器A Method and Adaptive Filter for Removing Narrowband Interference

技术领域technical field

本发明涉及信号处理领域,特别涉及一种去除窄带干扰的方法和自适应滤波器。The invention relates to the field of signal processing, in particular to a method for removing narrowband interference and an adaptive filter.

背景技术Background technique

正交频分复用(Orthogonal Frequency Division Multiplexing,OFDM)技术是一种多载波数字调制技术,也是一种频率复用技术。相比较其他调制技术和复用技术,OFDM技术具有高效的频谱利用效率和良好的抗多径干扰能力,因此已经被广泛应用于广播式的音频和视频领域,以及军用和民用通信系统中。OFDM技术主要的应用包括:高比特率数字用户线系统(High-speed DigitalSubscriber Line,HDSL),非对称的数字用户环路(Asymmetric Digital SubscriberLine,ADSL)、ETSI(European Telecommunications Standards Institute,欧洲电信标准化协会)标准的数字音频广播(Digital Audio Broadcasting,DAB)、数字视频广播(Digital Video Broadcasting,DVB)、高清晰度电视(High DefinitionTelevision,HDTV)、无线城域网(WMAN)、无线局域网(WLAN),等等。特别是在无线城域网和无线局域网领域,以OFDM技术作为主要物理层技术的802.16系列和802.11系列标准已经得到了大量应用。Orthogonal Frequency Division Multiplexing (OFDM) technology is a multi-carrier digital modulation technology and a frequency multiplexing technology. Compared with other modulation techniques and multiplexing techniques, OFDM technology has efficient spectrum utilization efficiency and good anti-multipath interference ability, so it has been widely used in broadcast audio and video fields, as well as military and civilian communication systems. The main applications of OFDM technology include: High-speed Digital Subscriber Line (HDSL), Asymmetric Digital Subscriber Line (ADSL), ETSI (European Telecommunications Standards Institute, European Telecommunications Standards Institute) ) standard digital audio broadcasting (Digital Audio Broadcasting, DAB), digital video broadcasting (Digital Video Broadcasting, DVB), high definition television (High DefinitionTelevision, HDTV), wireless metropolitan area network (WMAN), wireless local area network (WLAN), etc. Especially in the field of wireless metropolitan area network and wireless local area network, the 802.16 series and 802.11 series standards with OFDM technology as the main physical layer technology have been widely used.

OFDM技术是以大量正交子载波来承载需要传输的信息。因此,OFDM系统极易受到窄带干扰的不利影响。此处的窄带干扰即包括通常意义上的窄带干扰,而且还包括单频干扰。因为经过OFDM系统中的FFT(快速傅立叶变换)变换后,单频信号的频谱不再是线谱而是占有一定带宽的窄带谱。当窄带干扰持续存在时,将会对系统造成长时间的干扰,严重影响OFDM系统的性能。而当窄带干扰在频谱上与用于信道估计的导频子载波重叠时,这种不利影响将会通过信道估计误差得到扩散,使得整个系统的性能恶化程度更加严重。OFDM technology uses a large number of orthogonal sub-carriers to carry information that needs to be transmitted. Therefore, OFDM systems are extremely vulnerable to adverse effects of narrowband interference. The narrow-band interference here includes narrow-band interference in the usual sense, and also includes single-frequency interference. Because after the FFT (Fast Fourier Transform) transformation in the OFDM system, the spectrum of the single-frequency signal is no longer a line spectrum but a narrow-band spectrum occupying a certain bandwidth. When the narrowband interference persists, it will cause long-term interference to the system and seriously affect the performance of the OFDM system. However, when the narrowband interference overlaps with the pilot subcarriers used for channel estimation on the spectrum, this adverse effect will be diffused through channel estimation errors, making the performance of the entire system worse.

目前用于干扰抵消的方法是通过构造一个复数系数的固定系数滤波器来抵消窄带干扰。但是在构造滤波器系数时,要求已知窄带干扰的中心频率或者其较高精度的估计值,否则不能得到满意的干扰抵消效果。固定系数的滤波器主要的缺点是系数不能根据环境变化实时调整,因此其性能不能满足时变的环境要求;并且在构造某些固定系数的滤波器的系数时,需要关于信号统计特性的先验知识,这些先验知识往往是无法预知的,或者是需要实时估计的,特别是干扰的中心频率,不可能得到其真实值或者很高精度的估计值。因此,通过固定系数滤波器除去干扰会引入处理时延,实时性比较差。The current method for interference cancellation is to cancel narrowband interference by constructing a fixed coefficient filter with complex coefficients. However, when constructing the filter coefficients, it is required to know the center frequency of the narrowband interference or its high-precision estimation value, otherwise a satisfactory interference cancellation effect cannot be obtained. The main disadvantage of the filter with fixed coefficients is that the coefficients cannot be adjusted in real time according to environmental changes, so its performance cannot meet the time-varying environmental requirements; and when constructing the coefficients of certain fixed-coefficient filters, prior knowledge about the statistical characteristics of the signal is required Knowledge, these prior knowledge are often unpredictable, or need to be estimated in real time, especially the center frequency of interference, it is impossible to get its true value or a very high-precision estimated value. Therefore, removing interference through a fixed coefficient filter will introduce processing delay, and the real-time performance is relatively poor.

综上所述,目前通过固定系数滤波器除去干扰会引入处理时延,实时性比较差。To sum up, at present, removing interference through a fixed coefficient filter will introduce processing delay, and the real-time performance is relatively poor.

发明内容Contents of the invention

本发明实施例提供一种去除窄带干扰的方法和自适应滤波器,用以解决现有技术中存在的通过固定系数滤波器除去干扰会引入处理时延,实时性比较差的问题。Embodiments of the present invention provide a method for removing narrowband interference and an adaptive filter to solve the problem in the prior art that removing interference through a fixed-coefficient filter introduces processing delay and poor real-time performance.

本发明实施例提供的一种去除窄带干扰的方法,该方法包括:A method for removing narrowband interference provided by an embodiment of the present invention includes:

自适应滤波器确定当前采样时刻对应的权值;The adaptive filter determines the weight corresponding to the current sampling moment;

所述自适应滤波器根据确定的权值和在当前采样时刻采样的参考信号,确定当前采样时刻对应的抵消信号,其中所述参考信号是根据窄带干扰的中心频率生成的;The adaptive filter determines the cancellation signal corresponding to the current sampling time according to the determined weight and the reference signal sampled at the current sampling time, wherein the reference signal is generated according to the center frequency of the narrowband interference;

所述自适应滤波器根据确定的当前采样时刻对应的抵消信号,对在当前采样时刻采样的基带信号进行滤波,得到除去窄带干扰后的基带信号。The adaptive filter filters the baseband signal sampled at the current sampling moment according to the determined cancellation signal corresponding to the current sampling moment, to obtain the baseband signal after the narrowband interference is removed.

本发明实施例提供的一种自适应滤波器,该自适应滤波器包括:An adaptive filter provided by an embodiment of the present invention includes:

权值确定模块,用于确定当前采样时刻对应的权值;A weight determination module, configured to determine the weight corresponding to the current sampling moment;

信号确定模块,用于根据确定的权值和在当前采样时刻采样的参考信号,确定当前采样时刻对应的抵消信号,其中所述参考信号是根据窄带干扰的中心频率生成的;The signal determination module is used to determine the corresponding cancellation signal at the current sampling time according to the determined weight and the reference signal sampled at the current sampling time, wherein the reference signal is generated according to the center frequency of the narrowband interference;

滤波模块,用于根据确定的当前采样时刻对应的抵消信号,对在当前采样时刻采样的基带信号进行滤波,得到除去窄带干扰后的基带信号。The filtering module is configured to filter the baseband signal sampled at the current sampling moment according to the determined canceling signal corresponding to the current sampling moment, so as to obtain the baseband signal after the narrowband interference is removed.

本发明实施例自适应滤波器根据当前采样时刻对应的权值和在当前采样时刻采样的参考信号,确定当前采样时刻对应的抵消信号,以及根据确定的当前采样时刻对应的抵消信号,对在当前采样时刻采样的基带信号进行滤波,得到除去窄带干扰后的基带信号。According to the weight value corresponding to the current sampling moment and the reference signal sampled at the current sampling moment, the adaptive filter in the embodiment of the present invention determines the cancellation signal corresponding to the current sampling moment, and according to the determined cancellation signal corresponding to the current sampling moment, the The baseband signal sampled at the sampling time is filtered to obtain the baseband signal after the narrowband interference is removed.

由于本发明实施例通过自适应滤波器除去窄带干扰,从而减少了处理时延,提高了实时性,并且采用自适应滤波器能够实时调整权值,满足时变的环境要求改善干扰抵消效果。Since the embodiment of the present invention removes narrow-band interference through an adaptive filter, the processing delay is reduced and real-time performance is improved, and the adaptive filter can be used to adjust the weight in real time to meet the time-varying environmental requirements and improve the interference cancellation effect.

进一步的,本发明实施例的实现简单,复杂度低,反应速度快,更加适合正交频分复用系统。Furthermore, the implementation of the embodiment of the present invention is simple, the complexity is low, and the response speed is fast, which is more suitable for the OFDM system.

附图说明Description of drawings

图1为本发明实施例自适应滤波器的结构示意图;FIG. 1 is a schematic structural diagram of an adaptive filter according to an embodiment of the present invention;

图2为本发明实施例去除窄带干扰的方法流程示意图;FIG. 2 is a schematic flowchart of a method for removing narrowband interference according to an embodiment of the present invention;

图3A为本发明实施例利用Notch(凹口)滤波器进行滤波原理示意图中;3A is a schematic diagram of the principle of filtering using a Notch (notch) filter according to an embodiment of the present invention;

图3B为本发明实施例Notch滤波器的处理原理示意图。FIG. 3B is a schematic diagram of a processing principle of a Notch filter according to an embodiment of the present invention.

具体实施方式Detailed ways

本发明实施例自适应滤波器根据当前采样时刻对应的权值和在当前采样时刻采样的参考信号,确定当前采样时刻对应的抵消信号,以及根据确定的当前采样时刻对应的抵消信号,对在当前采样时刻采样的基带信号进行滤波,得到除去窄带干扰后的基带信号。由于本发明实施例通过自适应滤波器除去窄带干扰,从而减少了处理时延,提高了实时性,并且采用自适应滤波器能够实时调整权值,满足时变的环境要求改善干扰抵消效果。According to the weight value corresponding to the current sampling moment and the reference signal sampled at the current sampling moment, the adaptive filter in the embodiment of the present invention determines the cancellation signal corresponding to the current sampling moment, and according to the determined cancellation signal corresponding to the current sampling moment, the The baseband signal sampled at the sampling time is filtered to obtain the baseband signal after the narrowband interference is removed. Since the embodiment of the present invention removes narrow-band interference through an adaptive filter, the processing delay is reduced and real-time performance is improved, and the adaptive filter can be used to adjust the weight in real time to meet the time-varying environmental requirements and improve the interference cancellation effect.

其中,参考信号是根据窄带干扰的中心频率生成的。Wherein, the reference signal is generated according to the center frequency of the narrowband interference.

具体的,对有无窄带干扰进行检测和判定,在有窄带干扰时对窄带干扰的中心频率进行估计,得到窄带干扰中心频率估计值fr,此频率将会作为自适应滤波器参考信号的中心频率。Specifically, detect and determine whether there is narrow-band interference, estimate the center frequency of narrow-band interference when there is narrow-band interference, and obtain the estimated value f r of the narrow-band interference center frequency, which will be used as the center of the adaptive filter reference signal frequency.

在具体实施过程中,可以通过窄带干扰检测器对有无窄带干扰进行检测和判定,并且对窄带干扰的中心频率进行估计,得到窄带干扰中心频率估计值。In a specific implementation process, the presence or absence of narrowband interference may be detected and determined by a narrowband interference detector, and the center frequency of the narrowband interference may be estimated to obtain an estimated value of the center frequency of the narrowband interference.

需要说明的是,本发明实施例并不局限于窄带干扰检测器,其他能够对有无窄带干扰进行检测和判定,以及对窄带干扰的中心频率进行估计得到窄带干扰中心频率估计值的方式都适用本发明实施例。It should be noted that the embodiment of the present invention is not limited to the narrowband interference detector, and other methods that can detect and determine the presence or absence of narrowband interference, and estimate the center frequency of narrowband interference to obtain the estimated value of the narrowband interference center frequency are applicable Embodiment of the present invention.

其中,根据窄带干扰的中心频率就可以生成参考信号。Wherein, the reference signal can be generated according to the center frequency of the narrowband interference.

比如需要两个参考信号,则可以根据公式一和公式二生成两路正交单频参考信号rs和rc,二者分别为以fr为中心频率的正弦参考信号和余弦参考信号:For example, if two reference signals are needed, two orthogonal single-frequency reference signals rs and r c can be generated according to formula 1 and formula 2, which are sine reference signals and cosine reference signals with f r as the center frequency respectively:

rs(k)=Csin(2πfrk)..............公式一;r s (k)=Csin(2πf r k)...Formula 1;

rc(k)=Ccos(2πfrk).............公式二。r c (k)=Ccos(2πf r k)..........Formula 2.

其中,rs(k)表示第k个采样时刻采样的正弦参考信号;rc(k)表示第k个采样时刻采样的余弦参考信号;C是参考信号的幅值,可以根据需要进行设定。Among them, rs (k) represents the sine reference signal sampled at the kth sampling moment; r c (k) represents the cosine reference signal sampled at the kth sampling moment; C is the amplitude of the reference signal, which can be set as required .

由于参考信号为单频信号,因此可以只存储一个周期内的采样值,在滤波时循环调用,以节省存储器资源;如果窄带干扰中心频率为缓慢时变的,可以通过查表计算参考信号采样值。Since the reference signal is a single-frequency signal, it can only store the sampling value in one cycle, and call it cyclically during filtering to save memory resources; if the center frequency of narrow-band interference is slowly time-varying, the sampling value of the reference signal can be calculated by looking up the table .

需要说明的是,本发明实施例并不局限于上述采用公式一和公式二生成参考信号的方式,其他能够根据窄带干扰的中心频率生成参考信号的方式都适用本发明实施例。It should be noted that the embodiment of the present invention is not limited to the above-mentioned method of generating the reference signal using Formula 1 and Formula 2, and other methods that can generate the reference signal according to the center frequency of the narrowband interference are applicable to the embodiment of the present invention.

在具体实施过程中,可以通过设置自适应滤波器参数得到一个合理的自适应滤波器的滤波带宽,而且自适应滤波器可以自动调整自身的滤波频带中心频率,即可以跟踪窄带干扰的中心频率,因此不需要已知窄带干扰中心频率的精确值,只需要一个有一定精度的大概值。In the specific implementation process, a reasonable filter bandwidth of the adaptive filter can be obtained by setting the parameters of the adaptive filter, and the adaptive filter can automatically adjust the center frequency of its own filter frequency band, that is, it can track the center frequency of narrow-band interference, Therefore, it is not necessary to know the exact value of the central frequency of the narrow-band interference, but only an approximate value with certain precision.

其中,本发明实施例可以应用于OFDM系统,还可以应用于以大量正交子载波来承载需要传输的信息的系统。Wherein, the embodiment of the present invention can be applied to an OFDM system, and can also be applied to a system in which information to be transmitted is carried by a large number of orthogonal subcarriers.

下面结合说明书附图对本发明实施例作进一步详细描述。The embodiments of the present invention will be further described in detail below in conjunction with the accompanying drawings.

如图1所示,本发明实施例自适应滤波器包括:权值确定模块10、信号确定模块20和滤波模块30。As shown in FIG. 1 , the adaptive filter according to the embodiment of the present invention includes: a weight determination module 10 , a signal determination module 20 and a filtering module 30 .

权值确定模块10,用于确定当前采样时刻对应的权值。The weight determination module 10 is configured to determine the weight corresponding to the current sampling moment.

信号确定模块20,用于根据权值确定模块10确定的权值和在当前采样时刻采样的参考信号,确定当前采样时刻对应的抵消信号,其中参考信号是根据窄带干扰的中心频率生成的。The signal determination module 20 is configured to determine the corresponding cancellation signal at the current sampling time according to the weight determined by the weight determination module 10 and the reference signal sampled at the current sampling time, wherein the reference signal is generated according to the center frequency of the narrowband interference.

滤波模块30,用于根据信号确定模块20确定的当前采样时刻对应的抵消信号,对在当前采样时刻采样的基带信号进行滤波,得到除去窄带干扰后的基带信号。The filtering module 30 is configured to filter the baseband signal sampled at the current sampling moment according to the cancellation signal corresponding to the current sampling moment determined by the signal determining module 20 to obtain the baseband signal after narrowband interference is removed.

其中,权值确定模块10根据上一个采样时刻对应的权值、上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的参考信号,确定当前采样时刻对应的权值。Wherein, the weight determination module 10 determines the weight corresponding to the current sampling time according to the weight corresponding to the previous sampling time, the baseband signal after removing the narrowband interference at the previous sampling time, and the reference signal sampled at the previous sampling time.

具体的权值确定模块10可以采用自行应算法确定当前采样时刻对应的权值。自适应算法有很多种,比如最小均方误差(Least Mean Square,LMS)算法、递归最小二乘(Recursive Least Squares,RLS)算法,归一化LMS算法,变步长LMS算法等。The specific weight determination module 10 may determine the weight corresponding to the current sampling moment by adopting an adaptive algorithm. There are many kinds of adaptive algorithms, such as Least Mean Square (LMS) algorithm, Recursive Least Squares (RLS) algorithm, normalized LMS algorithm, variable step size LMS algorithm, etc.

下面以最小均方误差算法为例进行说明,其他算法与最小均方误差算法类似,只是算法不同,在此不再赘述。The following uses the minimum mean square error algorithm as an example for illustration. Other algorithms are similar to the minimum mean square error algorithm, but the algorithm is different, and will not be repeated here.

权值确定模块10如果采用最小均方误差算法可以根据公式三确定当前采样时刻对应的权值:If the weight value determination module 10 adopts the minimum mean square error algorithm, it can determine the weight value corresponding to the current sampling moment according to formula three:

w ( k ) = w ( k - 1 ) + μ [ y ^ ( k - 1 ) + n ^ ( k - 1 ) ] r ( k - 1 ) ...........公式三。 w ( k ) = w ( k - 1 ) + μ [ the y ^ ( k - 1 ) + no ^ ( k - 1 ) ] r ( k - 1 ) ..........Formula three.

其中,w(k)是当前采样时刻对应的权值;w(k-1)是上一个采样时刻对应的权值;μ是自适应算法的迭代计算步长;

Figure G2009102254497D00061
是上一个采样时刻除去窄带干扰后的基带信号中的有用信号,是上一个采样时刻除去窄带干扰后的基带信号中除有用信号外的其他信号(比如噪声干扰等),
Figure G2009102254497D00063
Figure G2009102254497D00064
构成了上一个采样时刻除去窄带干扰后的基带信号;r(k-1)是在上一个采样时刻采样的参考信号。Among them, w(k) is the weight corresponding to the current sampling moment; w(k-1) is the weight corresponding to the previous sampling moment; μ is the iterative calculation step of the adaptive algorithm;
Figure G2009102254497D00061
is the useful signal in the baseband signal after removing the narrowband interference at the last sampling moment, It is the signal other than the useful signal (such as noise interference, etc.) in the baseband signal after removing the narrowband interference at the last sampling moment,
Figure G2009102254497D00063
and
Figure G2009102254497D00064
It constitutes the baseband signal after removing the narrowband interference at the last sampling moment; r(k-1) is the reference signal sampled at the last sampling moment.

如果当前采样时刻是第一个采样时刻,则上一个采样时刻对应的权值就是预先设定的权值,上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的参考信号都是0。If the current sampling moment is the first sampling moment, the weight corresponding to the last sampling moment is the preset weight, and the baseband signal after removing the narrowband interference at the last sampling moment and the reference signal sampled at the last sampling moment are both is 0.

其中,不同的自适应滤波器需要的权值数量也不相同,以Notch滤波器为例,需要两个权值,则权值包括第一权值和第二权值;相应的,参考信号包括第一参考信号和第二参考信号。Among them, the number of weights required by different adaptive filters is also different. Taking the Notch filter as an example, two weights are required, and the weights include the first weight and the second weight; correspondingly, the reference signal includes A first reference signal and a second reference signal.

权值确定模块10需要根据上一个采样时刻对应的第一权值、上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的第一参考信号,确定当前采样时刻对应的第一权值,以及根据上一个采样时刻对应的第二权值、上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的第二参考信号,确定当前采样时刻对应的第二权值。The weight determination module 10 needs to determine the first weight value corresponding to the current sampling time according to the first weight value corresponding to the previous sampling time, the baseband signal after removing the narrowband interference at the last sampling time, and the first reference signal sampled at the last sampling time. Weight, and according to the second weight corresponding to the last sampling moment, the baseband signal after removing the narrowband interference at the last sampling moment, and the second reference signal sampled at the last sampling moment, determine the second weight corresponding to the current sampling moment .

以最小均方误差算法为例,则根据公式三可以衍生出公式四和公式五,权值确定模块10根据公式四确定第一权值,根据公式五确定第二权值:Taking the minimum mean square error algorithm as an example, formula 4 and formula 5 can be derived according to formula 3, the weight determination module 10 determines the first weight according to formula 4, and determines the second weight according to formula 5:

w s ( k ) = w s ( k - 1 ) + μ [ y ^ ( k - 1 ) + n ^ ( k - 1 ) ] r s ( k - 1 ) ...........公式四; w the s ( k ) = w the s ( k - 1 ) + μ [ the y ^ ( k - 1 ) + no ^ ( k - 1 ) ] r the s ( k - 1 ) .....Formula 4;

其中,ws(k)是当前采样时刻对应的第一权值;ws(k-1)是上一个采样时刻对应的第一权值;μ是自适应算法的迭代计算步长;

Figure G2009102254497D00066
是上一个采样时刻除去窄带干扰后的基带信号中的有用信号,
Figure G2009102254497D00067
是上一个采样时刻除去窄带干扰后的基带信号中处有用信号外的其他信号(比如噪声干扰等),
Figure G2009102254497D00068
Figure G2009102254497D00071
构成了上一个采样时刻除去窄带干扰后的基带信号;rs(k-1)是在上一个采样时刻采样的第一参考信号。Among them, w s (k) is the first weight corresponding to the current sampling moment; w s (k-1) is the first weight corresponding to the previous sampling moment; μ is the iterative calculation step size of the adaptive algorithm;
Figure G2009102254497D00066
is the useful signal in the baseband signal after removing the narrowband interference at the last sampling moment,
Figure G2009102254497D00067
It is the signal other than the useful signal (such as noise interference, etc.) in the baseband signal after removing the narrowband interference at the last sampling moment,
Figure G2009102254497D00068
and
Figure G2009102254497D00071
constitutes the baseband signal after removing the narrowband interference at the last sampling moment; r s (k-1) is the first reference signal sampled at the last sampling moment.

w c ( k ) = w c ( k - 1 ) + μ [ y ^ ( k - 1 ) + n ^ ( k - 1 ) ] r c ( k - 1 ) ...........公式五; w c ( k ) = w c ( k - 1 ) + μ [ the y ^ ( k - 1 ) + no ^ ( k - 1 ) ] r c ( k - 1 ) ......... Formula five;

其中,wc(k)是当前采样时刻对应的第一权值;wc(k-1)是上一个采样时刻对应的第二权值;μ是自适应算法的迭代计算步长;

Figure G2009102254497D00073
是上一个采样时刻除去窄带干扰后的基带信号中的有用信号,
Figure G2009102254497D00074
是上一个采样时刻除去窄带干扰后的基带信号中处有用信号外的其他信号(比如噪声干扰等),
Figure G2009102254497D00075
Figure G2009102254497D00076
构成了上一个采样时刻除去窄带干扰后的基带信号;rc(k-1)是在上一个采样时刻采样的第二参考信号。Among them, w c (k) is the first weight value corresponding to the current sampling moment; w c (k-1) is the second weight value corresponding to the previous sampling moment; μ is the iterative calculation step size of the adaptive algorithm;
Figure G2009102254497D00073
is the useful signal in the baseband signal after removing the narrowband interference at the last sampling moment,
Figure G2009102254497D00074
It is the signal other than the useful signal (such as noise interference, etc.) in the baseband signal after removing the narrowband interference at the last sampling moment,
Figure G2009102254497D00075
and
Figure G2009102254497D00076
constitutes the baseband signal after the narrowband interference has been removed at the last sampling moment; r c (k-1) is the second reference signal sampled at the last sampling moment.

需要说明的是,如果第一权值是根据第一参考信号确定的,则后续采样时刻对应的第一权值都必须根据第一参考信号;相应的,如果第二权值是根据第一参考信号确定的,则后续采样时刻对应的第二权值都必须根据第二参考信号。It should be noted that if the first weight value is determined according to the first reference signal, the first weight value corresponding to the subsequent sampling time must be based on the first reference signal; correspondingly, if the second weight value is determined according to the first reference signal If the signal is determined, the second weights corresponding to subsequent sampling moments must be based on the second reference signal.

比如第一参考信号是正弦参考信号,第二参考信号是余弦参考信号。如果第一个采样时刻对应的第一权值是根据正弦参考信号确定的,第一个采样时刻对应的第二权值是根据余弦参考信号确定的,则后续采样时刻对应的第一权值必须根据正弦参考信号确定,后续采样时刻对应的第二权值必须根据余弦参考信号确定。For example, the first reference signal is a sine reference signal, and the second reference signal is a cosine reference signal. If the first weight corresponding to the first sampling moment is determined according to the sine reference signal, and the second weight corresponding to the first sampling moment is determined according to the cosine reference signal, then the first weight corresponding to the subsequent sampling moment must Determined according to the sine reference signal, the second weight corresponding to the subsequent sampling moment must be determined according to the cosine reference signal.

对于信号确定模块20可以根据公式六确定当前采样时刻对应的抵消信号:For the signal determination module 20, the offset signal corresponding to the current sampling moment can be determined according to formula six:

n ^ I ( k ) = w ( k ) r ( k ) ...............公式六; no ^ I ( k ) = w ( k ) r ( k ) ..........Formula 6;

其中,

Figure G2009102254497D00078
是抵消信号;w(k)是当前采样时刻对应的权值;r(k)是在当前采样时刻采样的参考信号。in,
Figure G2009102254497D00078
is the offset signal; w(k) is the weight corresponding to the current sampling moment; r(k) is the reference signal sampled at the current sampling moment.

如果本发明实施例的自适应滤波器需要两个权值,则信号确定模块20将确定的当前采样时刻对应的第一权值和当前采样时刻采样的第一参考信号相乘,得到当前采样时刻对应的第一处理信号,将确定的当前采样时刻对应的第二权值和当前采样时刻采样的第二参考信号相乘,得到当前采样时刻对应的第二处理信号,将第一处理信号和第二处理信号相加,得到当前采样时刻对应的抵消信号。If the adaptive filter of the embodiment of the present invention needs two weights, the signal determination module 20 multiplies the determined first weight corresponding to the current sampling time with the first reference signal sampled at the current sampling time to obtain the current sampling time The corresponding first processed signal is multiplied by the determined second weight value corresponding to the current sampling time and the second reference signal sampled at the current sampling time to obtain the second processed signal corresponding to the current sampling time, and the first processed signal and the second The two processed signals are added to obtain the cancellation signal corresponding to the current sampling moment.

具体的,根据公式六可以衍生出公式七,信号确定模块20可以根据公式七得到当前采样时刻对应的抵消信号:Specifically, formula 7 can be derived according to formula 6, and the signal determination module 20 can obtain the offset signal corresponding to the current sampling moment according to formula 7:

n ^ I ( k ) = w s ( k ) r s ( k ) + w c ( k ) r c ( k ) .................公式七; no ^ I ( k ) = w the s ( k ) r the s ( k ) + w c ( k ) r c ( k ) ..........Formula VII;

其中,

Figure G2009102254497D00082
是当前采样时刻对应的抵消信号;ws(k)是当前采样时刻对应的第一权值;rs(k)是在当前采样时刻采样的第一参考信号;ws(k)rs(k)是当前采样时刻对应的第一处理信号;wc(k)是当前采样时刻对应的第二权值;rc(k)是在当前采样时刻采样的第二参考信号;wc(k)rc(k)是当前采样时刻对应的第二处理信号。in,
Figure G2009102254497D00082
is the offset signal corresponding to the current sampling moment; w s (k) is the first weight corresponding to the current sampling moment; rs (k) is the first reference signal sampled at the current sampling moment; w s (k) rs ( k) is the first processed signal corresponding to the current sampling moment; w c (k) is the second weight value corresponding to the current sampling moment; r c (k) is the second reference signal sampled at the current sampling moment; w c (k ) rc (k) is the second processed signal corresponding to the current sampling moment.

其中,滤波模块30可以将在当前采样时刻采样的基带信号和确定的当前采样时刻对应的抵消信号相减,得到的信号是除去窄带干扰后的基带信号。Wherein, the filtering module 30 may subtract the baseband signal sampled at the current sampling moment from the determined cancellation signal corresponding to the current sampling moment, and the obtained signal is the baseband signal after narrowband interference is removed.

具体的,滤波模块30可以根据公式八得到的信号是除去窄带干扰后的基带信号:Specifically, the signal obtained by the filtering module 30 according to Formula 8 is the baseband signal after the narrowband interference is removed:

y ^ ( k ) + n ^ ( k ) = s ( k ) - n ^ I ( k ) ............公式八; the y ^ ( k ) + no ^ ( k ) = the s ( k ) - no ^ I ( k ) .........Formula eight;

其中,

Figure G2009102254497D00084
是当前采样时刻除去窄带干扰后的基带信号中的有用信号,
Figure G2009102254497D00085
是当前采样时刻除去窄带干扰后的基带信号中除有用信号外的其他信号,
Figure G2009102254497D00086
构成了当前采样时刻除去窄带干扰后的基带信号;s(k)是在当前采样时刻采样的基带信号;
Figure G2009102254497D00087
是当前采样时刻确定的当前采样时刻对应的抵消信号。in,
Figure G2009102254497D00084
is the useful signal in the baseband signal after removing the narrowband interference at the current sampling moment,
Figure G2009102254497D00085
It is the signal other than the useful signal in the baseband signal after removing the narrowband interference at the current sampling time,
Figure G2009102254497D00086
Constitute the baseband signal after the narrowband interference is removed at the current sampling moment; s(k) is the baseband signal sampled at the current sampling moment;
Figure G2009102254497D00087
is the cancellation signal corresponding to the current sampling time determined by the current sampling time.

如图2所示,本发明实施例去除窄带干扰的方法包括下列步骤:As shown in Figure 2, the method for removing narrowband interference in the embodiment of the present invention includes the following steps:

步骤201、自适应滤波器确定当前采样时刻对应的权值。Step 201, the adaptive filter determines the weight corresponding to the current sampling moment.

步骤202、自适应滤波器根据确定的权值和在当前采样时刻采样的参考信号,确定当前采样时刻对应的抵消信号,其中参考信号是根据窄带干扰的中心频率生成的。Step 202, the adaptive filter determines the corresponding cancellation signal at the current sampling time according to the determined weight and the reference signal sampled at the current sampling time, wherein the reference signal is generated according to the center frequency of the narrowband interference.

步骤203、自适应滤波器根据确定的当前采样时刻对应的抵消信号,对在当前采样时刻采样的基带信号进行滤波,得到除去窄带干扰后的基带信号。Step 203 , the adaptive filter filters the baseband signal sampled at the current sampling moment according to the determined cancellation signal corresponding to the current sampling moment, to obtain the baseband signal after narrowband interference is removed.

步骤201中,自适应滤波器根据上一个采样时刻对应的权值、上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的参考信号,确定当前采样时刻对应的权值。In step 201, the adaptive filter determines the weight corresponding to the current sampling time according to the weight corresponding to the previous sampling time, the baseband signal after removing the narrowband interference at the previous sampling time, and the reference signal sampled at the previous sampling time.

具体的自适应滤波器可以采用自行应算法确定当前采样时刻对应的权值。自适应算法有很多种,比如最小均方误差算法、递归最小二乘算法,归一化LMS算法,变步长LMS算法等。A specific adaptive filter may use a self-response algorithm to determine the weight corresponding to the current sampling moment. There are many kinds of adaptive algorithms, such as minimum mean square error algorithm, recursive least squares algorithm, normalized LMS algorithm, variable step size LMS algorithm, etc.

下面以最小均方误差算法为例进行说明,其他算法与最小均方误差算法类似,只是算法不同,在此不再赘述。The following uses the minimum mean square error algorithm as an example for illustration. Other algorithms are similar to the minimum mean square error algorithm, but the algorithm is different, and will not be repeated here.

自适应滤波器如果采用最小均方误差算法可以根据公式三确定当前采样时刻对应的权值。If the adaptive filter adopts the minimum mean square error algorithm, the weight corresponding to the current sampling moment can be determined according to formula three.

如果当前采样时刻是第一个采样时刻,则上一个采样时刻对应的权值就是预先设定的权值,上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的参考信号都是0。If the current sampling moment is the first sampling moment, the weight corresponding to the last sampling moment is the preset weight, and the baseband signal after removing the narrowband interference at the last sampling moment and the reference signal sampled at the last sampling moment are both is 0.

其中,不同的自适应滤波器需要的权值数量也不相同,以Notch滤波器为例,需要两个权值,则权值包括第一权值和第二权值;相应的,参考信号包括第一参考信号和第二参考信号。Among them, the number of weights required by different adaptive filters is also different. Taking the Notch filter as an example, two weights are required, and the weights include the first weight and the second weight; correspondingly, the reference signal includes A first reference signal and a second reference signal.

步骤201中,自适应滤波器需要根据上一个采样时刻对应的第一权值、上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的第一参考信号,确定当前采样时刻对应的第一权值,以及根据上一个采样时刻对应的第二权值、上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的第二参考信号,确定当前采样时刻对应的第二权值。In step 201, the adaptive filter needs to determine the weight corresponding to the current sampling time according to the first weight value corresponding to the previous sampling time, the baseband signal after removing narrowband interference at the previous sampling time, and the first reference signal sampled at the previous sampling time. , and according to the second weight value corresponding to the last sampling moment, the baseband signal after removing the narrowband interference at the last sampling moment, and the second reference signal sampled at the last sampling moment, determine the first weight corresponding to the current sampling moment Two weights.

以最小均方误差算法为例,则根据公式三可以衍生出公式四和公式五,自适应滤波器根据公式四确定第一权值,根据公式五确定第二权值。Taking the minimum mean square error algorithm as an example, formula 4 and formula 5 can be derived according to formula 3. The adaptive filter determines the first weight value according to formula 4, and determines the second weight value according to formula 5.

需要说明的是,如果第一权值是根据第一参考信号确定的,则后续采样时刻对应的第一权值都必须根据第一参考信号;相应的,如果第二权值是根据第一参考信号确定的,则后续采样时刻对应的第二权值都必须根据第二参考信号。It should be noted that if the first weight value is determined according to the first reference signal, the first weight value corresponding to the subsequent sampling time must be based on the first reference signal; correspondingly, if the second weight value is determined according to the first reference signal If the signal is determined, the second weights corresponding to subsequent sampling moments must be based on the second reference signal.

比如第一参考信号是正弦参考信号,第二参考信号是余弦参考信号。如果第一个采样时刻对应的第一权值是根据正弦参考信号确定的,第一个采样时刻对应的第二权值是根据余弦参考信号确定的,则后续采样时刻对应的第一权值必须根据正弦参考信号确定,后续采样时刻对应的第二权值必须根据余弦参考信号确定。For example, the first reference signal is a sine reference signal, and the second reference signal is a cosine reference signal. If the first weight corresponding to the first sampling moment is determined according to the sine reference signal, and the second weight corresponding to the first sampling moment is determined according to the cosine reference signal, then the first weight corresponding to the subsequent sampling moment must Determined according to the sine reference signal, the second weight corresponding to the subsequent sampling moment must be determined according to the cosine reference signal.

步骤202中,对于自适应滤波器可以根据公式六确定当前采样时刻对应的抵消信号。In step 202, for the adaptive filter, the cancellation signal corresponding to the current sampling moment can be determined according to Formula 6.

如果本发明实施例的自适应滤波器需要两个权值,则步骤202中,自适应滤波器将确定的当前采样时刻对应的第一权值和当前采样时刻采样的第一参考信号相乘,得到当前采样时刻对应的第一处理信号,将确定的当前采样时刻对应的第二权值和当前采样时刻采样的第二参考信号相乘,得到当前采样时刻对应的第二处理信号,将第一处理信号和第二处理信号相加,得到当前采样时刻对应的抵消信号。If the adaptive filter of the embodiment of the present invention needs two weights, then in step 202, the adaptive filter multiplies the determined first weight corresponding to the current sampling moment by the first reference signal sampled at the current sampling moment, Obtain the first processed signal corresponding to the current sampling moment, multiply the determined second weight value corresponding to the current sampling moment by the second reference signal sampled at the current sampling moment, obtain the second processed signal corresponding to the current sampling moment, and multiply the first The processed signal is added to the second processed signal to obtain a cancellation signal corresponding to the current sampling moment.

具体的,根据公式六可以衍生出公式七,自适应滤波器可以根据公式七得到当前采样时刻对应的抵消信号。Specifically, Formula 7 can be derived according to Formula 6, and the adaptive filter can obtain the cancellation signal corresponding to the current sampling moment according to Formula 7.

步骤203中,自适应滤波器可以将在当前采样时刻采样的基带信号和确定的当前采样时刻对应的抵消信号相减,得到的信号是除去窄带干扰后的基带信号。In step 203, the adaptive filter may subtract the baseband signal sampled at the current sampling moment from the determined cancellation signal corresponding to the current sampling moment, and the obtained signal is the baseband signal after narrowband interference is removed.

具体的,自适应滤波器可以根据公式八得到的信号是除去窄带干扰后的基带信号。Specifically, the signal obtained by the adaptive filter according to Formula 8 is the baseband signal after narrowband interference is removed.

本发明实施例的自适应滤波器可以是Notch滤波器,比如凹口滤波器、陷波器等,也可以是其他自适应滤波器。下面以本发明实施例的自适应滤波器是Notch滤波器为例进行说明。The adaptive filter in this embodiment of the present invention may be a Notch filter, such as a notch filter, a notch filter, etc., or other adaptive filters. In the following, description will be made by taking the adaptive filter of the embodiment of the present invention as a Notch filter as an example.

如图3A所示,本发明实施例利用Notch滤波器进行滤波原理示意图中:As shown in FIG. 3A, in the schematic diagram of the principle of filtering using the Notch filter in the embodiment of the present invention:

首先,窄带干扰检测器对有无窄带干扰进行检测和判定,在有窄带干扰时对窄带干扰的中心频率进行估计,得到窄带干扰中心频率估计值frFirst, the narrowband interference detector detects and judges whether there is narrowband interference, and estimates the center frequency of the narrowband interference when there is narrowband interference, and obtains an estimated value f r of the narrowband interference center frequency.

然后,根据窄带干扰的中心频率生成正弦参考信号rs和余弦参考信号rcThen, a sine reference signal rs and a cosine reference signal rc are generated according to the center frequency of the narrowband interference.

最后,Notch滤波器根据正弦参考信号rs和余弦参考信号rc,对输入的基带信号(s=y+n+nI)进行滤波,得到除去窄带干扰后的基带信号

Figure G2009102254497D00111
和抵消信号
Figure G2009102254497D00112
Finally, the Notch filter filters the input baseband signal (s=y+n+n I ) according to the sine reference signal rs and cosine reference signal rc to obtain the baseband signal after narrowband interference is removed
Figure G2009102254497D00111
and offset signal
Figure G2009102254497D00112

如图3B所示,本发明实施例Notch滤波器的处理原理示意图中,有三路信号,其中一路是基带信号s,包括有用信号y,窄带干扰信号nI和其他干扰信号n;另一路是正弦参考信号rs;最后一路是余弦参考信号rcAs shown in Figure 3B, in the schematic diagram of the processing principle of the Notch filter of the embodiment of the present invention, there are three signals, one of which is a baseband signal s, including a useful signal y, a narrowband interference signal n 1 and other interference signals n; the other is a sinusoidal signal Reference signal r s ; the last path is cosine reference signal r c .

其中,正弦参考信号rs和余弦参考信号rc存储到Notch滤波器中,比如存储一个周期内的采样值,则Notch滤波器会接收到一路信号,即基带信号S;也可以与基带信号s类似将正弦参考信号rs和余弦参考信号rc发送给Notch滤波器,则Notch滤波器会接收到三路信号,即基带信号S、正弦参考信号rs和余弦参考信号rcAmong them, the sine reference signal r s and the cosine reference signal r c are stored in the Notch filter, such as storing the sampling value in one cycle, then the Notch filter will receive a signal, that is, the baseband signal S; it can also be combined with the baseband signal s Similar to sending the sine reference signal rs and the cosine reference signal rc to the Notch filter, the Notch filter will receive three signals, namely the baseband signal S, the sine reference signal rs and the cosine reference signal rc .

Notch滤波器根据当前采样时刻采样的正弦参考信号rs(k)和当前采样时刻对应的第一权值ws(k),确定当前采样时刻对应的第一处理信号 The Notch filter determines the first processed signal corresponding to the current sampling time according to the sinusoidal reference signal r s (k) sampled at the current sampling time and the first weight w s (k) corresponding to the current sampling time

Notch滤波器根据根据当前采样时刻采样的余弦参考信号rc(k)和当前采样时刻对应的第二权值wc(k),确定当前采样时刻对应的第二处理信号

Figure G2009102254497D00114
The Notch filter determines the second processed signal corresponding to the current sampling time according to the cosine reference signal r c (k) sampled at the current sampling time and the second weight w c (k) corresponding to the current sampling time
Figure G2009102254497D00114

Notch滤波器根据当前采样时刻对应的第一处理信号

Figure G2009102254497D00115
和当前采样时刻对应的第二处理信号
Figure G2009102254497D00121
得到当前采样时刻对应的抵消信号
Figure G2009102254497D00122
Notch filter according to the first processing signal corresponding to the current sampling moment
Figure G2009102254497D00115
The second processed signal corresponding to the current sampling moment
Figure G2009102254497D00121
Get the offset signal corresponding to the current sampling moment
Figure G2009102254497D00122

Notch滤波器根据当前采样时刻对应的抵消信号

Figure G2009102254497D00123
对在当前采样时刻采样的基带信号s(k)进行滤波,得到当前采样时刻对应的除去窄带干扰后的基带信号
Figure G2009102254497D00124
The Notch filter cancels the signal according to the current sampling moment
Figure G2009102254497D00123
Filter the baseband signal s(k) sampled at the current sampling moment to obtain the baseband signal corresponding to the current sampling moment after removing narrowband interference
Figure G2009102254497D00124

在经过一段收敛过程后,Notch滤波器的两个权值就会收敛接近于最小均方误差准则下的最优线性滤波权值(假设采用最小均方误差算法),这时得到的就是抵消窄带干扰后的输出信号

Figure G2009102254497D00125
以及窄带干扰信号的估计波形
Figure G2009102254497D00126
After a period of convergence, the two weights of the Notch filter will converge close to the optimal linear filter weights under the minimum mean square error criterion (assuming that the minimum mean square error algorithm is used), then the result is to offset the narrowband Output signal after interference
Figure G2009102254497D00125
and the estimated waveform of the narrowband interfering signal
Figure G2009102254497D00126

在自适应收敛过程结束后,

Figure G2009102254497D00127
将是y+n最小均方误差估计,而
Figure G2009102254497D00128
将是nI最小均方误差估计。
Figure G2009102254497D00129
可以直接用下一步处理,例如信道估计、解调、译码等等,
Figure G2009102254497D001210
可以用于对原始信号的校准和修正等。After the adaptive convergence process ends,
Figure G2009102254497D00127
will be the y+n minimum mean square error estimate, while
Figure G2009102254497D00128
will be the n I minimum mean square error estimate.
Figure G2009102254497D00129
It can be processed directly in the next step, such as channel estimation, demodulation, decoding, etc.,
Figure G2009102254497D001210
It can be used for calibration and correction of the original signal, etc.

本发明实施例的自适应滤波可以采用多种硬件器件实现,不仅可以采用通用的处理器、单片机以及专用的数字信号处理器,还可以采用可编程逻辑器件等。The adaptive filtering in the embodiment of the present invention can be realized by using various hardware devices, not only general-purpose processors, single-chip microcomputers and dedicated digital signal processors, but also programmable logic devices and the like.

从上述实施例可以看出:本发明实施例自适应滤波器确定当前采样时刻对应的权值;自适应滤波器根据确定的权值和在当前采样时刻采样的参考信号,确定当前采样时刻对应的抵消信号,其中参考信号是根据窄带干扰的中心频率生成的;自适应滤波器根据确定的当前采样时刻对应的抵消信号,对在当前采样时刻采样的基带信号进行滤波,得到除去窄带干扰后的基带信号。It can be seen from the above embodiments that the adaptive filter in the embodiment of the present invention determines the weight corresponding to the current sampling moment; the adaptive filter determines the weight corresponding to the current sampling moment according to the determined weight and the reference signal sampled at the current sampling moment. Cancellation signal, wherein the reference signal is generated according to the center frequency of the narrowband interference; the adaptive filter filters the baseband signal sampled at the current sampling moment according to the determined cancellation signal corresponding to the current sampling moment, and obtains the baseband signal after removing the narrowband interference Signal.

本发明实施例充分利用了循环前缀时间作为自适应滤波器的收敛时间。The embodiment of the present invention fully utilizes the cyclic prefix time as the convergence time of the adaptive filter.

由于本发明实施例通过自适应滤波器除去窄带干扰,从而减少了处理时延,提高了实时性,并且采用自适应滤波器能够实时调整权值,满足时变的环境要求改善干扰抵消效果。Since the embodiment of the present invention removes narrow-band interference through an adaptive filter, the processing delay is reduced and real-time performance is improved, and the adaptive filter can be used to adjust the weight in real time to meet the time-varying environmental requirements and improve the interference cancellation effect.

进一步的,本发明实施例的实现简单,复杂度低,反应速度快,更加适合正交频分复用系统。Furthermore, the implementation of the embodiment of the present invention is simple, the complexity is low, and the response speed is fast, which is more suitable for the OFDM system.

显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and equivalent technologies thereof, the present invention also intends to include these modifications and variations.

Claims (10)

1.一种去除窄带干扰的方法,其特征在于,该方法包括:1. A method for removing narrowband interference, characterized in that the method comprises: 自适应滤波器确定当前采样时刻对应的权值;The adaptive filter determines the weight corresponding to the current sampling moment; 所述自适应滤波器根据确定的权值和在当前采样时刻采样的参考信号,确定当前采样时刻对应的抵消信号,其中所述参考信号是根据窄带干扰的中心频率生成的;The adaptive filter determines the cancellation signal corresponding to the current sampling time according to the determined weight and the reference signal sampled at the current sampling time, wherein the reference signal is generated according to the center frequency of the narrowband interference; 所述自适应滤波器根据确定的当前采样时刻对应的抵消信号,对在当前采样时刻采样的基带信号进行滤波,得到除去窄带干扰后的基带信号。The adaptive filter filters the baseband signal sampled at the current sampling moment according to the determined cancellation signal corresponding to the current sampling moment, to obtain the baseband signal after the narrowband interference is removed. 2.如权利要求1所述的方法,其特征在于,所述自适应滤波器确定当前采样点对应的权值包括:2. The method according to claim 1, wherein the adaptive filter determines the weight corresponding to the current sampling point comprising: 所述自适应滤波器根据上一个采样时刻对应的权值、上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的参考信号,确定当前采样时刻对应的权值。The adaptive filter determines the weight corresponding to the current sampling time according to the weight corresponding to the previous sampling time, the baseband signal after removing the narrowband interference at the previous sampling time, and the reference signal sampled at the previous sampling time. 3.如权利要求2所述的方法,其特征在于,所述权值包括第一权值和第二权值,所述参考信号包括第一参考信号和第二参考信号;3. The method according to claim 2, wherein the weight comprises a first weight and a second weight, and the reference signal comprises a first reference signal and a second reference signal; 所述自适应滤波器确定当前采样点对应的权值包括:The adaptive filter determining the weight corresponding to the current sampling point includes: 所述自适应滤波器根据上一个采样时刻对应的第一权值、上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的第一参考信号,确定当前采样时刻对应的第一权值;The adaptive filter determines the first weight corresponding to the current sampling time according to the first weight value corresponding to the previous sampling time, the baseband signal after removing narrowband interference at the previous sampling time, and the first reference signal sampled at the previous sampling time. Weight; 所述自适应滤波器根据上一个采样时刻对应的第二权值、上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的第二参考信号,确定当前采样时刻对应的第二权值。The adaptive filter determines the second weight corresponding to the current sampling time according to the second weight value corresponding to the previous sampling time, the baseband signal after removing the narrowband interference at the previous sampling time, and the second reference signal sampled at the previous sampling time. weight. 4.如权利要求3所述的方法,其特征在于,所述自适应滤波器确定当前采样时刻对应的抵消信号包括:4. The method according to claim 3, wherein the adaptive filter determines the corresponding cancellation signal at the current sampling moment comprising: 所述自适应滤波器将确定的当前采样时刻对应的第一权值和当前采样时刻采样的第一参考信号相乘,得到当前采样时刻对应的第一处理信号;The adaptive filter multiplies the determined first weight corresponding to the current sampling moment by the first reference signal sampled at the current sampling moment to obtain a first processed signal corresponding to the current sampling moment; 所述自适应滤波器将确定的当前采样时刻对应的第二权值和当前采样时刻采样的第二参考信号相乘,得到当前采样时刻对应的第二处理信号;The adaptive filter multiplies the determined second weight corresponding to the current sampling moment by the second reference signal sampled at the current sampling moment to obtain a second processed signal corresponding to the current sampling moment; 所述自适应滤波器将所述第一处理信号和所述第二处理信号相加,得到当前采样时刻对应的抵消信号。The adaptive filter adds the first processed signal and the second processed signal to obtain a cancellation signal corresponding to the current sampling moment. 5.如权利要求1~4任一权利要求所述的方法,其特征在于,所述自适应滤波器对在当前采样时刻采样的基带信号进行滤波包括:5. The method according to any one of claims 1 to 4, wherein the adaptive filter filtering the baseband signal sampled at the current sampling moment comprises: 所述自适应滤波器将在当前采样时刻采样的基带信号和确定的当前采样时刻对应的抵消信号相减,得到的信号是除去窄带干扰后的基带信号。The adaptive filter subtracts the baseband signal sampled at the current sampling moment from the determined cancellation signal corresponding to the current sampling moment, and the obtained signal is the baseband signal after narrowband interference is removed. 6.一种自适应滤波器,其特征在于,该自适应滤波器包括:6. An adaptive filter, characterized in that the adaptive filter comprises: 权值确定模块,用于确定当前采样时刻对应的权值;A weight determination module, configured to determine the weight corresponding to the current sampling moment; 信号确定模块,用于根据确定的权值和在当前采样时刻采样的参考信号,确定当前采样时刻对应的抵消信号,其中所述参考信号是根据窄带干扰的中心频率生成的;The signal determination module is used to determine the corresponding cancellation signal at the current sampling time according to the determined weight and the reference signal sampled at the current sampling time, wherein the reference signal is generated according to the center frequency of the narrowband interference; 滤波模块,用于根据确定的当前采样时刻对应的抵消信号,对在当前采样时刻采样的基带信号进行滤波,得到除去窄带干扰后的基带信号。The filtering module is configured to filter the baseband signal sampled at the current sampling moment according to the determined canceling signal corresponding to the current sampling moment, so as to obtain the baseband signal after the narrowband interference is removed. 7.如权利要求6所述的自适应滤波器,其特征在于,所述权值确定模块用于:7. adaptive filter as claimed in claim 6, is characterized in that, described weight determination module is used for: 根据上一个采样时刻对应的权值、上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的参考信号,确定当前采样时刻对应的权值。The weight corresponding to the current sampling time is determined according to the weight corresponding to the previous sampling time, the baseband signal after removing the narrowband interference at the previous sampling time, and the reference signal sampled at the previous sampling time. 8.如权利要求7所述的自适应滤波器,其特征在于,所述权值确定模块用于:8. adaptive filter as claimed in claim 7, is characterized in that, described weight determination module is used for: 在所述权值包括第一权值和第二权值,所述参考信号包括第一参考信号和第二参考信号时,根据上一个采样时刻对应的第一权值、上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的第一参考信号,确定当前采样时刻对应的第一权值,以及根据上一个采样时刻对应的第二权值、上一个采样时刻除去窄带干扰后的基带信号和在上一个采样时刻采样的第二参考信号,确定当前采样时刻对应的第二权值。When the weight includes the first weight and the second weight, and the reference signal includes the first reference signal and the second reference signal, the narrowband is removed according to the first weight corresponding to the last sampling moment and the last sampling moment Based on the interfered baseband signal and the first reference signal sampled at the last sampling moment, determine the first weight corresponding to the current sampling moment, and remove the narrowband interference based on the second weight corresponding to the last sampling moment and the last sampling moment The baseband signal and the second reference signal sampled at the last sampling time are used to determine the second weight corresponding to the current sampling time. 9.如权利要求8所述的自适应滤波器,其特征在于,所述信号确定模块用于:9. adaptive filter as claimed in claim 8, is characterized in that, described signal determination module is used for: 将确定的当前采样时刻对应的第一权值和当前采样时刻采样的第一参考信号相乘,得到当前采样时刻对应的第一处理信号,将确定的当前采样时刻对应的第二权值和当前采样时刻采样的第二参考信号相乘,得到当前采样时刻对应的第二处理信号,将所述第一处理信号和所述第二处理信号相加,得到当前采样时刻对应的抵消信号。Multiplying the determined first weight corresponding to the current sampling moment with the first reference signal sampled at the current sampling moment to obtain the first processed signal corresponding to the current sampling moment, and multiplying the determined second weight corresponding to the current sampling moment with the current The second reference signal sampled at the sampling time is multiplied to obtain a second processed signal corresponding to the current sampling time, and the first processed signal and the second processed signal are added to obtain a cancellation signal corresponding to the current sampling time. 10.如权利要求6~9任一权利要求所述的自适应滤波器,其特征在于,所述滤波模块用于:10. The adaptive filter according to any one of claims 6 to 9, wherein the filtering module is used for: 将在当前采样时刻采样的基带信号和确定的当前采样时刻对应的抵消信号相减,得到的信号是除去窄带干扰后的基带信号。The baseband signal sampled at the current sampling moment is subtracted from the determined cancellation signal corresponding to the current sampling moment, and the obtained signal is the baseband signal after the narrowband interference is removed.
CN 200910225449 2009-12-14 2009-12-14 Method for removing narrow-band interference and self-adapting filter Expired - Fee Related CN102098258B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN 200910225449 CN102098258B (en) 2009-12-14 2009-12-14 Method for removing narrow-band interference and self-adapting filter
PCT/CN2010/072409 WO2010148785A1 (en) 2009-12-14 2010-05-04 Method for removing narrowband interference and adaptive filter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200910225449 CN102098258B (en) 2009-12-14 2009-12-14 Method for removing narrow-band interference and self-adapting filter

Publications (2)

Publication Number Publication Date
CN102098258A true CN102098258A (en) 2011-06-15
CN102098258B CN102098258B (en) 2013-11-06

Family

ID=43385933

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200910225449 Expired - Fee Related CN102098258B (en) 2009-12-14 2009-12-14 Method for removing narrow-band interference and self-adapting filter

Country Status (2)

Country Link
CN (1) CN102098258B (en)
WO (1) WO2010148785A1 (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103004086A (en) * 2012-09-27 2013-03-27 华为技术有限公司 Self-adaptive filtering method and system based on error sub-band
CN103248595A (en) * 2013-05-08 2013-08-14 桂林电子科技大学 Method and device for self-adaptive in-place interference offset
CN103609085A (en) * 2011-06-14 2014-02-26 三星电子株式会社 Apparatus and method for receiving signal in wireless communication system
CN104660217A (en) * 2015-01-09 2015-05-27 深圳市美的连电子科技有限公司 Signal analysis method and terminal
CN106027439A (en) * 2015-03-27 2016-10-12 意法半导体(鲁塞)公司 Method for processing an analog signal coming from a transmission channel, in particular a signal carried by power line communications
CN106789777A (en) * 2015-11-19 2017-05-31 北京科技大学 Wear digital interference in steel number energy simultaneous interpretation system eliminate, bypassing method and device
WO2017097161A1 (en) * 2015-12-11 2017-06-15 Huawei Technologies Co., Ltd. Time domain ingress noise detection and cancellation
CN111862925A (en) * 2020-07-03 2020-10-30 天津大学 An adaptive active noise control system and method based on lazy learning

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7720185B2 (en) * 2006-11-06 2010-05-18 Qualcomm Incorporated Narrow-band interference canceller
KR100906383B1 (en) * 2007-07-26 2009-07-06 인하대학교 산학협력단 Narrowband Interference Rejection Method for Ultra-Wideband Ranging Systems
CN101594166B (en) * 2008-05-31 2012-09-05 华为技术有限公司 Method and device for eliminating interference in signals in short distance wireless network

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9680593B2 (en) 2011-06-14 2017-06-13 Samsung Electronics Co., Ltd. Apparatus and method for receiving signal in wireless communication system
CN103609085A (en) * 2011-06-14 2014-02-26 三星电子株式会社 Apparatus and method for receiving signal in wireless communication system
WO2014047854A1 (en) * 2012-09-27 2014-04-03 华为技术有限公司 Adaptive filtering method and system based on error sub-band
CN103004086B (en) * 2012-09-27 2016-06-08 华为技术有限公司 Adaptive filter method and system based on error subband
US9419826B2 (en) 2012-09-27 2016-08-16 Huawei Technologies Co., Ltd. Adaptive filtering method and system based on error sub-band
CN103004086A (en) * 2012-09-27 2013-03-27 华为技术有限公司 Self-adaptive filtering method and system based on error sub-band
CN103248595A (en) * 2013-05-08 2013-08-14 桂林电子科技大学 Method and device for self-adaptive in-place interference offset
CN103248595B (en) * 2013-05-08 2016-01-20 桂林电子科技大学 A kind of self adaptation is with location interference cancellation method and device
CN104660217A (en) * 2015-01-09 2015-05-27 深圳市美的连电子科技有限公司 Signal analysis method and terminal
CN104660217B (en) * 2015-01-09 2018-01-23 深圳市美的连医疗电子股份有限公司 A kind of signal analysis method and terminal
CN106027439A (en) * 2015-03-27 2016-10-12 意法半导体(鲁塞)公司 Method for processing an analog signal coming from a transmission channel, in particular a signal carried by power line communications
CN106027439B (en) * 2015-03-27 2019-07-23 意法半导体(鲁塞)公司 To the method and receiver handled from the channel analog signal for sending channel
CN106789777A (en) * 2015-11-19 2017-05-31 北京科技大学 Wear digital interference in steel number energy simultaneous interpretation system eliminate, bypassing method and device
CN106789777B (en) * 2015-11-19 2019-11-29 北京科技大学 Wear digital interference elimination, bypassing method and the device in steel number energy simultaneous interpretation system
WO2017097161A1 (en) * 2015-12-11 2017-06-15 Huawei Technologies Co., Ltd. Time domain ingress noise detection and cancellation
US9860078B2 (en) 2015-12-11 2018-01-02 Futurewei Technologies, Inc. Time domain ingress noise detection and cancellation
CN111862925A (en) * 2020-07-03 2020-10-30 天津大学 An adaptive active noise control system and method based on lazy learning
CN111862925B (en) * 2020-07-03 2024-04-12 天津大学 Adaptive active noise control system and method based on inertia learning

Also Published As

Publication number Publication date
WO2010148785A1 (en) 2010-12-29
CN102098258B (en) 2013-11-06

Similar Documents

Publication Publication Date Title
CN102098258B (en) Method for removing narrow-band interference and self-adapting filter
KR100878430B1 (en) Staggered pilot transmission method and apparatus for channel estimation and time tracking
KR101061742B1 (en) Jammer Detection and Suppression for Wireless Communications
CN101888351B (en) Channel estimation device and method
JP5302887B2 (en) Transmission method and transmission apparatus for canceling inter-carrier interference
Yucek et al. Carrier frequency offset compensation with successive cancellation in uplink OFDMA systems
JP2007195218A (en) Digital radio frequency interference canceller
CN107659527B (en) Phase noise suppression system and method for simultaneous same-frequency full duplex cooperative communication
US20090207925A1 (en) Wireless communication system, OFDM communication apparatus and method thereof
TWI234950B (en) Mode detection for OFDM signals
WO2007020943A1 (en) Ofdm communication method
GB2548293A (en) Noise variance estimation and interference detection
US7340000B1 (en) Decision feedback equalizer in an OFDM WLAN receiver
CN105337911B (en) A kind of removing method and device of inter-sub-carrier interference
US7023929B2 (en) Digital pre-compensation filter for DMT type transmitter
Ambatali et al. Performance evaluation of the UFMC scheme under various transmission impairments
CN106797361B (en) System and method for circular convolution
CN102843326B (en) A receiver for OFDM system and its signal receiving method
TWI617146B (en) Device and method of handling effective path of channel impulse response
JP2001024619A (en) Ofdm signal receiver
CN109245839B (en) Method for detecting narrow-band interference, communication chip and communication device thereof
CN110661735A (en) A Delay Estimation Method for Sparse Channel Estimation Based on Orthogonal Matching Tracking
Sun et al. Optimal pilot based frequency-dependent I/Q imbalance compensation for wideband direct-conversion transmitters
Marelli et al. A subband approach to channel estimation and equalization for DMT and OFDM systems
CN109120560B (en) Self-adaptive ICI inter-carrier interference elimination method and device

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20131106