<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JCC</journal-id><journal-title-group><journal-title>Journal of Computer and Communications</journal-title></journal-title-group><issn pub-type="epub">2327-5219</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jcc.2021.95005</article-id><article-id pub-id-type="publisher-id">JCC-109204</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject></subj-group></article-categories><title-group><article-title>
 
 
  Multi-Frequency Interference Detection and Mitigation Using Multiple Adaptive IIR Notch Filter with Lattice Structure
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Abdelrahman</surname><given-names>El Gebali</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>René</surname><given-names>Jr Landry</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Lassena Laboratory, Department of Electrical Engineering, école de Technologie Supérieure (éTS), Montreal, Canada</addr-line></aff><pub-date pub-type="epub"><day>12</day><month>05</month><year>2021</year></pub-date><volume>09</volume><issue>05</issue><fpage>58</fpage><lpage>77</lpage><history><date date-type="received"><day>20,</day>	<month>April</month>	<year>2021</year></date><date date-type="rev-recd"><day>17,</day>	<month>May</month>	<year>2021</year>	</date><date date-type="accepted"><day>20,</day>	<month>May</month>	<year>2021</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Radio Frequency Interferences (RFI), such as strong Continuous Wave Interferences (CWI), can influence the Quality of Service (QoS) of communications, increasing the Bit Error Rate (BER) and decreasing the Signal-to-Noise Ratio (SNR) in any wireless transmission, including in a Digital Video Broadcasting (DVB-S2) receiver. Therefore, this paper presents an algorithm for detecting and mitigating a Multi-tone Continuous Wave Interference (MCWI) using a Multiple Adaptive Notch Filter (MANF), based on the lattice form structure. The Adaptive Notch Filter (ANF) is constructed using the second-order IIR NF. The approach consists in developing a robust low-complexity algorithm for removing unknown MCWI. The MANF model is a multistage model, with each stage consisting of two ANFs: the adaptive IIR notch filter 
  H
  <sub>l</sub>(
  z) and the adaptive IIR notch filter 
  H
  <sub>N</sub>(
  z), which can detect and mitigate CWI. In this model, the ANF is used for estimating the Jamming-to-Signal Ratio (JSR) and the frequency of the interference (
  w(0)) by using an LMS-based algorithm. The depth of the notch is then adjusted based on the estimation of the JSR. In contrast, the ANF 
  H
  <sub>N</sub>(
  z) is used to mitigate the CW interference. Simulation results show that the proposed ANF is an effective method for eliminating/reducing the effects of MCWI, and yields better system performance than full suppression (
  k<sub>N</sub>=1) for low JSR values, and mostly the same performance for high JSR values. Moreover, the proposed can detect low and high JSR and track hopping frequency interference and provides better Bit error ratio (BER) performance compared to the case without an IIR notch filter.
 
</p></abstract><kwd-group><kwd>Component</kwd><kwd> Radio Frequency Interference (RFI)</kwd><kwd> Multiple Adaptive Notch Filter (MANF)</kwd><kwd> Jamming-to-Signal Ratio (JSR)</kwd><kwd> Quality of Service (QoS)</kwd><kwd> BER</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Many communications systems, such as wireless and satellite communication systems, suffer from Radio Frequency Interference (RFI) thanks to the ever-increasing number of wireless applications in use, as well as human-made RFI. RFI may be unintentional or could be an intentional act, such as jamming. It can also reduce the Signal of Interest (SoI) and result in a degraded QoS, a decrease in SNR, and an increase in BER or a loss of the receiver’s communication signal. Various jamming mitigation techniques have been proposed in the literature, and can be classified as adaptive antenna-based techniques (space domain) [<xref ref-type="bibr" rid="scirp.109204-ref1">1</xref>], Time-Frequency Filtering (TFF)-based methods [<xref ref-type="bibr" rid="scirp.109204-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref3">3</xref>], and adaptive filtering-based techniques [<xref ref-type="bibr" rid="scirp.109204-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref7">7</xref>]. The adaptive antenna methods (spatial filtering) are appropriate for Narrow-Band (NB) and wideband (WB) jamming mitigations. The space domain (antenna array) can also achieve a high jammer suppression ratio but is not suitable for mobile devices [<xref ref-type="bibr" rid="scirp.109204-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref9">9</xref>]. The TFF methods remove signal features in both the Time Domain (TD) and the Frequency Domain (FD) simultaneously, but the major concern with the FD method is computational complexity, which arises from FFT and IFFT in hardware. [<xref ref-type="bibr" rid="scirp.109204-ref10">10</xref>]. The methods are appropriate for NB jamming mitigation, such as Short Time Fourier transforms (STFT) [<xref ref-type="bibr" rid="scirp.109204-ref11">11</xref>], bilinear signal distribution [<xref ref-type="bibr" rid="scirp.109204-ref12">12</xref>], filter banks [<xref ref-type="bibr" rid="scirp.109204-ref13">13</xref>], and Wavelet Transforms (WT) [<xref ref-type="bibr" rid="scirp.109204-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref15">15</xref>]. The STFT and the filter bank techniques employ fixed windows, which are useful for analyzing nonstationary signals. The adaptive filtering-based methods are commonly used for narrowband CWI jamming rejection and can be classified into two approaches, namely, TD and FD. For the FD approach, [<xref ref-type="bibr" rid="scirp.109204-ref16">16</xref>] proposed an N-sigma method to reject CWI. Also, [<xref ref-type="bibr" rid="scirp.109204-ref17">17</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref18">18</xref>] proposed methods to detect and mitigate CWI and NB interference. Also, [<xref ref-type="bibr" rid="scirp.109204-ref19">19</xref>] proposed a method to detect and mitigate MCWI using a first-order IIR notch filter based on a simplified Welch algorithm and notch filter. However, the main concerns with the FD approach are computational complexity and higher costs. Using the Finite Impulse Response (FIR) and the Infinite Impulse Response (IIR) filters to detect and mitigate interference for TD approaches. The papers by [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref21">21</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref22">22</xref>] present methods for suppressing the CWI in the TD using an ANF based on a lattice form structure. The authors in [<xref ref-type="bibr" rid="scirp.109204-ref23">23</xref>] proposed a method for detecting and mitigating CW interference using two-pole and multiple NF. Also, [<xref ref-type="bibr" rid="scirp.109204-ref9">9</xref>] presented a technique using ANF to detect and mitigate MCWI in GPS. In [<xref ref-type="bibr" rid="scirp.109204-ref10">10</xref>], the ANF is used to mitigate both CWI and NB interference for GNSS. Moreover, the papers by [<xref ref-type="bibr" rid="scirp.109204-ref24">24</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref25">25</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref26">26</xref>] present linear prediction and nonlinear prediction methods. Both methods require longer tap filters to effectively mitigate narrowband interference, which increases the system’s complexity [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>]. In the literature, NF is used to reduce and remove interference in different applications, such as biomedical applications and [<xref ref-type="bibr" rid="scirp.109204-ref27">27</xref>] and GNSS applications (e.g., DSSS systems and GPS receivers) [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref23">23</xref>]. The lattice IIR ANF performs better than the direct IIR ANF form [<xref ref-type="bibr" rid="scirp.109204-ref21">21</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref22">22</xref>]. The IIR NF in lattice form is less computationally complex than are direct form implementations [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref21">21</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref23">23</xref>].</p><p>A typical IIR notch digital filter has constrained zeros on the unit circle, which makes the notch more depth “results into infinite” leading to complete removal of the interference. Nevertheless, this method creates self-noise due to the distortion of the data “desired signals” at the notch frequency as the interference is removed [<xref ref-type="bibr" rid="scirp.109204-ref22">22</xref>], called self-noise, as briefly described in [<xref ref-type="bibr" rid="scirp.109204-ref28">28</xref>]. In DSSS systems, Choi [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>] proposed a new CWI excision method that uses an ANF to reduce self-noise. The authors in [<xref ref-type="bibr" rid="scirp.109204-ref29">29</xref>] proposed a novel low-complexity anti-jamming filter to reduce the CWI for general wireless communications. With the exception of the works in [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>] and [<xref ref-type="bibr" rid="scirp.109204-ref29">29</xref>], these previous works concentrated on detecting and removing CWI and MCWI without focusing on controlling the notch filter’s depth. The two divergent papers [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref29">29</xref>] focused only on the excision of the CWI. The present paper proposes a low-complexity algorithm operating in TD using multiple ANF to detect and mitigate MCWI with considering controlling the notch depth based on an estimation of the interference power, called an anti-jamming receiver system model. This model is constructed using the second-order IIR NF based on a lattice form structure described in Section 3.</p><p>The rest of the paper is organized as follows. The signal model of the transmitted and received signals in the presence of MCWI is described in Section 2. Section 3 reviews and describes the adaptive IIR NF structure and the proposed adaptation algorithm for mitigating MCWI using MANF. Section 4 describes and defines the optimal notch depth that maximizes the output SNR, while Section 5 presents the simulation performed and discusses the results obtained. Finally, Section 6 concludes the paper.</p></sec><sec id="s2"><title>2. Signal Model</title><p>In this section, <xref ref-type="fig" rid="fig1">Figure 1</xref> presents the transmitted and received signal model, where a Quadrature Phase Shift Keying (QPSK) signal and an interference signal with the AWGN are generated using a MATLAB simulation model.</p><sec id="s2_1"><title>2.1. Interference Signal Model</title><p>Continuous Wave Interference (CWI) is one of the most common sources of interference. CWI can be modeled as a sinusoidal wave in time, such as a single or multi-CWI. A multitoned CWIs is considered in this paper, and can be modeled as:</p><p>J i ( t ) = ∑ i = 1 m A i cos ( 2 π f i t + θ i ) (1)</p><p>where:</p><p>• m is the number of CW interferences;</p><p>• A i is the amplitude of the i<sup>th</sup> interference signal;</p><p>• f i is the frequency of the i<sup>th</sup> interference signal;</p><p>• θ i is the phase delay of the i<sup>th</sup> interference signal, which is assumed to be a random variable uniformly distributed in the range [ − π , π ] .</p></sec><sec id="s2_2"><title>2.2. Received Signal Model</title><p>The received signal r ( t ) , which is the input signal to the IIR notch filter, is the sum of the QPSK modulated baseband signal S ( t ) , the AWGN, w ( t ) , and the interference signal, J i ( t ) , which can be modeled as:</p><p>r ( t ) = S ( t ) + w ( t ) + J i ( t ) (2)</p><p>In Section 3, the case of a multitoned CWI is considered, and m = 2 . The input baseband signal to the IIR notch filter r ( t ) given by Equation (2) is then sampled at chip rate to convert it into discrete-time samples for further processing as shown here:</p><p>r ( n T s ) = S ( n T s ) + w ( n T s ) + ∑ i = 1 m A i cos ( 2 π f i n T s + θ i ) (3)</p><p>For ease of notation, the sampling interval T s will be ignored in the rest of this paper, and Equation (3) will be adopted as:</p><p>r ( n ) = S ( n ) + w ( n ) + ∑ i = 1 m A i cos ( 2 π f i n + θ i ) (4)</p><p>where:</p><p>• S ( n ) is the desired QPSK modulated signal;</p><p>• w ( n ) is the AWGN.</p></sec></sec><sec id="s3"><title>3. Proposed Algorithm and Structure of the Filter</title><p>This paper proposes using a MANF based on a second-order IIR lattice form structure as an anti-jamming receiver system model, as shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>. The model presents as a low-complexity algorithm operating in the TD approach. The reason the approach operates in TD is to avoid the costs associated with the transfer domain approach. The MANF system model is a multistage model, with each of its stages consisting of two ANFs: the adaptive IIR notch filter H l ( z ) and the adaptive IIR notch filter H N ( z ) . The transfer functions of the two ANFs [ H l ( z ) and H N ( z ) ] are described in Section 3.1. The ANF H l ( z ) transfer function has zeros on the unit circle, resulting in the notch depth going</p><p>into infinite, and leading to the complete removal of the interference and exclusion of some desired signals. To avoid this, the ANF H N ( z ) is added to adjust and control the depth of the notch ( k N ) according to the estimated JSR<sub>l</sub> in order to reduce the elimination of some useful signals while removing interferences. The adaptive filter can be performed in various structures or realizations. There are two major implementations of an adaptive digital notch filter: the FIR and the IIR. FIR filters are also called moving average (MA) filters, which implement an all-zero transfer function. While the IIR filter has two parts, the auto-regression (AR) and the moving average (MA) (or an ARMA), the filter implements all-pole and all-zero transfer functions. In this work, the IIR filter is used since the IIR NF obtains a frequency response closer to that of the ideal NF than do FIR notch filters of the same length. <xref ref-type="fig" rid="fig3">Figure 3</xref> describes the adaptive IIR notch filter algorithm structure H l ( z ) [<xref ref-type="bibr" rid="scirp.109204-ref21">21</xref>] that attempts to adjust the k o l of the notch filter, which consists of two parts, the AR block and the MA block. The MA block is responsible for removing the interference signal, but it also removes the signal near the interference frequency. The ANF H l ( z ) is used to detect and estimate JSR<sub>l</sub> and the frequency of the interference ( ω o l ), while the ANF ( H N ( z ) ) is used for interference cancellation. Note that in this work, the number of states is chosen as two (k = 2).</p><p><xref ref-type="fig" rid="fig4">Figure 4</xref> introduces the multiple adaptive IIR NF implementation, similar to the cascade filter, but with an individual adaptive algorithm. The multiple adaptive IIR NF is implemented by having two different notch frequencies.</p><sec id="s3_1"><title>3.1. Transfer Function of IIR Notch Filter</title><p>This section introduces the second-order IIR NF transfer function H l ( z ) based on the lattice form structure. A decision was made to use the IIR NF with an LMS algorithm to design and implement the Adaptive Notch Filter (ANF). By indicating the adaptive coefficient of l<sup>th</sup> ANF by k o l , the second-order IIR NF transfer function H l ( z ) [<xref ref-type="bibr" rid="scirp.109204-ref16">16</xref>] is defined as:</p><p>H l ( z ) = N l ( z ) D l ( Z ) = 1 + 2 k o l Z − 1 + Z − 2 1 + k o l ( 1 + β ) Z − 1 + β Z − 2 (5)</p><p>where β is the pole radius factor that controls the filter’s bandwidth ( 0 &lt; β &lt; 1 ); k o l is a coefficient of lattice IIR NF parameter that rejects an unknown CWI, which is defined as k o l = − cos ( ω o l ) . If the absolute values of both β and k o l are less than 1, the H l ( z ) is stable. Also, in this work, the parameter β is fixed and only adapt k o l to avoid notch too much QPSK signal.</p></sec><sec id="s3_2"><title>3.2. Adaptation of Parameter k<sub>ol</sub> and LMS Adaptive Algorithm</title><p>In this subsection, the derivation of the LMS adaptive algorithm is utilized in transversal tapped (IIR) filters with the structure shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>. It is assumed that the MCWI and the AWGN are added to the QPSK signal described in Equation (4). y ( n ) and u ( n ) are the outputs of the MA and the AR block, respectively, as shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>, which can be defined in terms of the signal that was received r ( n ) as given by:</p><p>y ( n ) = u ( n ) + 2 k o l u ( n − 1 ) + u ( n − 2 ) (6)</p><p>u ( n ) = r ( n ) − k o l ( 1 + β ) u ( n − 1 ) − β u ( n − 2 ) (7)</p><p>The cost function of the filter can be given as [<xref ref-type="bibr" rid="scirp.109204-ref30">30</xref>]:</p><p>J ( k o l ) = E { | y ( n ) | 2 } (8)</p><p>where E { ∗ } is an expectation operator. Replacing the expectation value in Equation (8) with the ensemble-averaged value as shown in Equation (9), we get:</p><p>J ˙ ( k o l ) = 1 L ∑ n = 0 L − 1 y 2 ( n ) (9)</p><p>where L is the length of data. The goal of the filter is to adjust the notch parameter k o for minimizing the cost function, and the gradient-based method of steepest descent is used to achieve the goal:</p><p>k o l ( n + 1 ) = k o l ( n ) − μ ( d J ˙ ( k o l ) d k o l ( n ) ) (10)</p><p>The instantaneous value of J ˙ ( k o ) is used to replace the ensemble-averaged value as:</p><p>k o l ( n + 1 ) = k o l ( n ) − μ ( d y 2 ( n ) d k o l ( n ) ) (11)</p><p>To simplify, Equation (11) thus becomes:</p><p>k o l ( n + 1 ) = k o l ( n ) − μ y ( n ) ( g ( J ( k o l ) ) ) (12)</p><p>where g ( J ( k o l ) ) is the gradient signal for the adaptation of k o l , and is defined as:</p><p>g ( J ( k o l ) ) = − ( 1 + β ) x ( n − 1 ) (13)</p><p>Substituting into Equation (13):</p><p>k o l ( n + 1 ) = k o l ( n ) − μ y ( n ) ( ( 1 − β ) x ( n − 1 ) ) (14)</p><p>where &#181; is the step-size factor that controls the convergence speed and stability, and Equation (14) is a new adaptive algorithm for our proposed ANF that allows and controls the update filter coefficient in this work.</p><p>Since the normalized notch frequencies, f N and k o l , are related by k o l = − cos ( ω o l ) , where ω o l = 2 π f N , at sample n, the estimated frequency is given by:</p><p>f ^ N ( n ) = 1 2 π arccos ( k ^ o l ( n ) ) (15)</p><p>The estimated JSR<sub>l</sub> can be calculated by subtracting the output of the NF y ( n ) from the received signal, r ( n ) , and is given by:</p><p>JSR l dB = 10 log 10 { | r [ n ] − y [ n ] | 2 } (16)</p><p>The transfer function of the ANF H N ( z ) [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>] is given by:</p><p>H N ( z ) = N N ( z ) D N ( Z ) = 1 + K 0 l ( 1 + K N ) Z − 1 + K N Z − 2 1 + K 0 l ( 1 + β K N ) Z − 1 + β K N Z − 2 (17)</p><p>where k N is the depth of the NF and is a function of the estimated JSR<sub>l</sub>. β is the pole radius that controls the bandwidth of the NF. Note that N = l . While many implementation schemes can obtain H N ( z ) , this work, however, implements H N ( z ) by cascading the lattice IIR NF model [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>], as shown in <xref ref-type="fig" rid="fig5">Figure 5</xref>. If k N is set equal to 1 ( k N = 1 ), the result of H N ( z ) will be the same as that</p><p>of H l ( z ) , which has zero on the unit circle, resulting in an infinite depth of the notch that completely removes the interference [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>]. As stated previously, with the NF’s infinite depth, the interference is removed completely, but nonetheless creates self-noise at the NF, resulting from the exclusion of some useful signals. Hence, the depth of the notch k N should thus be adjusted with the interference power to reduce the distortion of the signal. Therefore, the same structure can be used for ANF H l ( z ) and ANF H N ( z ) , and the notch of Equation (17) can be placed on the IF by copying the k o l of ANF H l ( z ) that is tuned to IF by Equation (14) at any stage. These parameters, JSR<sub>l</sub>, k o l and k N , are needed to design the IIR notch filter H N ( z ) to maximize the output SNR at H N ( z ) .</p><p>The IIR NF H N ( z ) output can be defined as:</p><p>y ( n ) = H N ( z ) r ( n ) ≜ S o ( n ) + w o ( n ) + J o ( n ) (18)</p><p>where S o ( n ) , w o ( n ) and J o ( n ) are the output components of the desired QPSK modulated signal, the white Gaussian noise, and the MCWI.</p></sec></sec><sec id="s4"><title>4. Optimal Depth of the Notch that Maximizes SNR<sub>out</sub></title><p>This section described how the notch’s optimal depth maximizes the SNR output of the IIR NF as a function of its parameters [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref31">31</xref>] [<xref ref-type="bibr" rid="scirp.109204-ref32">32</xref>]. SNR<sub>out</sub> can be expressed as:</p><p>SNR out = E [ S 2 ( n ) ] E [ ( y ( n ) − S ( n ) ) 2 ] (19)</p><p>where y ( n ) is the output of the ANF. Equation (20) describes SNR<sub>out</sub> in terms of the filter parameters as given [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>]:</p><p>SNR out = 1 ( 1 + σ 2 ) ( ( 1 + k N 2 − 2 β k N 2 ) 1 − β 2 k N 2 ) + JSR l ( ( 1 − k N ) 2 ( 1 − β k N ) 2 ) − 1 (20)</p><p>where σ 2 is the variance of AWGN, and JSR l = A i 2 2 .</p><p>Let us assume that the first part of the Equation (20) denominator is D, called self-noise, as shown in Equation (21)</p><p>D = ( 1 + σ 2 ) ( 1 + k N 2 − 2 β k N 2 1 − β 2 k N 2 ) (21)</p><p>From Equation (20), the SNR output is affected by JSR<sub>l</sub> and D in Equation (21). D shows a large value since k N increases as the notch becomes deeper, whereas JSR<sub>l</sub> effect indicated by the second part of Equation (20). To maximize SNR<sub>out</sub>, the optimal value of k N needs to be found in Equation (20). The denominator in Equation (20) can be rewritten after some modification as a function of k N [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>]:</p><p>f ( k N ) = ( 1 + k N 2 − 2 β k N 2 1 − β 2 k N 2 ) + G l ( ( 1 − k N ) 2 ( 1 − β k N ) 2 ) − 1 1 + σ 2 (22)</p><p>where G l = JSR l 1 + σ 2 .</p><p>To find the optimal k N , we differentiate f ( k N ) and solve f ′ ( k N ) = 0 for all possible roots of k N , as shown in Equation (23):</p><p>f ′ ( k N ) = β 2 G l k N 3 + [ 2 β G − β 2 G l + β 2 − β ] k N 2     + [ 1 − β + G l − 2 β G ] k N − G l = 0 (23)</p><p>Equation (23) has at least one real root in the [0 1] range that gives the optimal k N as a function of JSR<sub>l</sub>. <xref ref-type="fig" rid="fig6">Figure 6</xref> describes the flow chart for the system model. After a signal is received, the ANF H l ( z ) estimates JSR<sub>l</sub> and the frequency of the interference ( ω o l ) using Equation (16) and Equation (15), respectively. Next, and based on the estimation of JSR<sub>l</sub> and f ^ N ( n ) , the estimated frequency f ^ N ( n ) is used to place the notch by calculating k o l = − cos ( 2 π f ^ N ) , and then using Equation (23) to calculate the optimal value of k N to maximize SNR<sub>out</sub>. Then, these parameters are used to design the IIR notch filter H N ( z ) as described by Equation (17).</p></sec><sec id="s5"><title>5. Simulation and Discussion of Results</title><p>The performance of the proposed low-complexity anti-jamming receiver system model of a MANF with MCWI is demonstrated in this section. The system performance is investigated in terms of the BER and the SNR<sub>out</sub> of the IIR NF for varying JSR and SOI power ( E b / N o ). In this work, for the i<sup>th</sup> jamming signal, the number of CWIs is considered, and thus, m = 2 with center frequencies ( f 1 = 0.03 and f 2 = 0.08 ) in Equation (2) and with equal magnitude in the input. However, frequencies could be any value in the [0 - 1] range. The simulation was run for 1e7 data bits for the MCWI. β is chosen to be 0.98 to give more satisfactory results for various situations. As β is chosen, there is a trade-off between some properties of ANF H l ( z ) and ANF H N ( z ) . A large β value gave more accurate JSR<sub>l</sub> and frequency estimates, but slower convergence and tracking [<xref ref-type="bibr" rid="scirp.109204-ref21">21</xref>] in ANF H l ( z ) . However, for ANF H N ( z ) , a large β gave a larger NB noise at the output, but a smaller broadband noise [<xref ref-type="bibr" rid="scirp.109204-ref20">20</xref>]. The BER is measured as JSR changes.</p><p>The estimated optimal value of k N vs. JSR is obtained by calculating the roots as JSR is changed, with a different parameter of β and E b / N o , as shown in <xref ref-type="fig" rid="fig7">Figure 7</xref>. As a result, k N approaches 1 as JSR<sub>l</sub> increases to a large value. Also, this approach depends on the parameters of ( E b / N o ) and β. Therefore, as the parameter of E b / N o decrease, JSR<sub>l</sub> increases as well, as β → 1 .</p><p><xref ref-type="fig" rid="fig8">Figure 8</xref>, part (a), describes the frequency response of the lattice IIR notch filter for various frequencies and a fixed notch width “ β = 0.98 ”, where part (b) describes the frequency response of the lattice IIR notch filter for various values of β. As β → 1 , the bandwidth of the notch is decreased as well as the depth of the magnitude responses of the notch, but in some cases, the depth of the notch is infinity.</p><p><xref ref-type="fig" rid="fig9">Figure 9</xref> and <xref ref-type="fig" rid="fig1">Figure 1</xref>0 show that the notch depth depends on the estimated JSR and interference frequency. Therefore, as JSR → increases, the notch magnitude response depth decrease, as does the notch bandwidth. On the other hand, the notch magnitude response depth depends on the power of the CWI; the higher power of the CWI, the deeper the notch magnitude response.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>1 and <xref ref-type="fig" rid="fig1">Figure 1</xref>2 depict the BER’s performance as E b / N o changes. The simulations were run with different values of E b / N o and with a fixed value of JSR. The results show that the simulation with “the proposed IIR NF” filter achieved a better BER in the presence of interference, especially for a low JSR and a large β. <xref ref-type="fig" rid="fig1">Figure 1</xref>3 and <xref ref-type="fig" rid="fig1">Figure 1</xref>4 show the performance of the BER vs. E b / N o in the presence of CWI and MCWI, with fixed values of JSR and different β parameters. The results show that the filter can be able to mitigate MCWI and achieve a better BER in the presence of CWI and MCWI, with differentβs.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>5 and <xref ref-type="fig" rid="fig1">Figure 1</xref>6 show the BER’s performance as JSR changes. Simulations were run with different values of JSR and a fixed E b / N o value, with filtering and without filtering. Therefore, the results show a comparison between the simulations with filtering and without filtering. Hence, the proposed ANF algorithm with filtering provides a better result than that without filtering. <xref ref-type="fig" rid="fig1">Figure 1</xref>7 shows a comparison between the proposed ANF algorithm and a full suppression ( k N = 1 ). It can be seen that the proposed algorithm effectively controls the notch depth and obtains better results than in the case of full suppression for a low JSR, and mostly the same for a high JSR.</p><p>The performance of the SNR vs. JSR is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>8 and <xref ref-type="fig" rid="fig1">Figure 1</xref>9. The proposed IIR NF algorithm effectively controls the notch depth, to be deeper for high JSR in each stage. The results show that full suppression achieves better results than does the case without a notch filter for a high JSR. Hence, full suppression is better for a high JSR, and without a notch filter is better for a low</p><p>JSR. The results also show that the proposed IIR NF at any stage approaches the case of full suppression when the JSR is increasing “high” because the notch depth becomes deeper in this case. So the notch depth becomes smaller for a lower JSR value and deeper for a higher JSR value.</p><p><xref ref-type="fig" rid="fig2">Figure 2</xref>0 shows the received signal of the Power Spectrum (PS) before and after processing by the IIR notch filter. The results show that the proposed IIR NF algorithm is able to detect and mitigate the interference. They also show that with the adjusted depth of the IIR NF, the interference is removed according to the estimated JSR<sub>l</sub> to avoid the reduction of the Signal of Interest (SoI).</p><p>These results will be compared with [<xref ref-type="bibr" rid="scirp.109204-ref19">19</xref>]. [<xref ref-type="bibr" rid="scirp.109204-ref19">19</xref>] proposed a method to detect and mitigate MCWI using a simplified Welch algorithm and notch filter operating in a frequency-domain approach. The approach used a first-order IIR notch filter. Although the CWI has been removed, a portion of the useful signal has been lost, compromising overall signal integrity. Also, the main concerns with the FD approach are computational complexity and higher costs. On the other hand, the proposed method proposes a low-complexity algorithm operating in a time-domain approach using multiple ANF to detect and mitigate MCWI. The approach used the second-order IIR notch filter based on a lattice form structure. The IIR notch depth is adjusted based on interference power estimation. <xref ref-type="table" rid="table1">Table 1</xref> shows a comparison between the algorithm proposed in [<xref ref-type="bibr" rid="scirp.109204-ref19">19</xref>] with our proposed algorithm.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Comparison between the algorithm proposed in [<xref ref-type="bibr" rid="scirp.109204-ref19">19</xref>] with our proposed algorithm</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >No.</th><th align="center" valign="middle" >Ref. [<xref ref-type="bibr" rid="scirp.109204-ref19">19</xref>]</th><th align="center" valign="middle" >Proposed algorithm</th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >They used QPSK signal</td><td align="center" valign="middle" >They used QPSK signal</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >They proposed a method to detect and mitigate MCWI using a simplified Welch algorithm and notch filter</td><td align="center" valign="middle" >We proposed a low-complexity algorithm to detect and mitigate MCWI using a multiple adaptive notch filter based on the lattice form structure</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >They used a frequency-domain approach</td><td align="center" valign="middle" >We used a time-domain approach</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >They used a first-order IIR adaptive notch filter</td><td align="center" valign="middle" >we used a second-order IIR adaptive notch filter</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >They do not adjust the depth of the notch</td><td align="center" valign="middle" >We adjusted the depth of the notch based on the estimation of the JSR</td></tr></tbody></table></table-wrap></sec><sec id="s6"><title>6. Conclusion</title><p>This paper proposes a novel low-complexity algorithm for mitigating multi-tone continuous wave interference using multiple adaptive notch filters based on a second-order IIR lattice form structure. This structure detects, estimates, and removes the MCWI. The algorithm operates in TD, which reduces hardware costs. Adaptive frequency estimators are used to adjust the notch filter’s zeros, placing them on the interference frequency. However, if zeros on the unit circle and has a notch of infinite depth ( k N = 1 ), and leading to the complete removal of the interference and exclusion of some desired signals This causes a reduction of the Signal of Interest (SoI) and results in a degraded QoS, a decrease in SNR, and an increase in BER. The adaptive notch filter algorithm for the given notch filter is thus proposed to adjust the notch depth and frequency. Therefore, the proposed method can effectively detect and mitigate the MCWI and adjust the notch depth for any given value of JSR, and can therefore effectively control the notch depth. It also provides better results than does a full suppression for low JSR values, and mostly the same performance for high JSR values, and provides a better BER performance. Also, the resulting SNR<sub>out</sub> of the NF is maximized for lower and higher JSR values with different values of ( E b / N o ) power. Therefore, this technique can be applied in a DVB-S2 receiver or any other communication and navigation receivers.</p></sec><sec id="s7"><title>Acknowledgements</title><p>The work reported in this paper was done under the AVIO 601-Interference Mitigation in Satellite Communication Project of the LASSENA Lab, &#201;cole de Technologie Sup&#233;rieure (&#201;TS). This research was also supported by the Natural Sciences and Engineering Research Council of Canada (NSERC), Thales, Telesat, Vigilant Global, CRIAQ, and Atem Canada.</p></sec><sec id="s8"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s9"><title>Cite this paper</title><p>El Gebali, A. and Landry, R.J. (2021) Multi-Frequency Interference Detection and Mitigation Using Multiple Adaptive IIR Notch Filter with Lattice Structure. Journal of Computer and Communications, 9, 58-77. https://doi.org/10.4236/jcc.2021.95005</p></sec></body><back><ref-list><title>References</title><ref id="scirp.109204-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Lu, D., Wu, R. and Liu, H. (2013) Global Positioning System Anti-Jamming Algorithm Based on Period Repetitive CLEAN. IET Radar, Sonar &amp; Navigation, 7, 164-169. https://doi.org/10.1049/iet-rsn.2010.0353</mixed-citation></ref><ref id="scirp.109204-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Lin, T., Abdizadeh, M., Broumandan, A., Wang, D., O’Keefe, K. and Lachapelle, G. (2011) Interference Suppression for High Precision Navigation Using Vector-Based GNSS Software Receivers. Proceedings of the 24th International Technical Meeting of The Satellite Division of the Institute of Navigation, Portland, 19-23 September 2011, 372-383.</mixed-citation></ref><ref id="scirp.109204-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Glennon, E.P. and Dempster, A.G. (2011) Delayed PIC for Postcorrelation Mitigation of Continuous Wave and Multiple Access Interference in GPS Receivers. IEEE Transactions on Aerospace and Electronic Systems, 47, 2544-2557. https://doi.org/10.1109/TAES.2011.6034650</mixed-citation></ref><ref id="scirp.109204-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Choi, J.W. and Cho, N.I. (2001) Narrow-Band Interference Suppression in Direct Sequence Spread Spectrum Systems Using a Lattice IIR Notch Filter. 2001 IEEE International Conference on Acoustics, Speech, and Signal Processing, Salt Lake City, 7-11 May 2001, 2237-2240. https://doi.org/10.1109/ICASSP.2001.940443</mixed-citation></ref><ref id="scirp.109204-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Barbarossa, S. and Scaglione, A. (1999) Adaptive Time-Varying Cancellation of Wideband Interferences in Spread-Spectrum Communications Based on Time-Frequency Distributions. IEEE Transactions on Signal Processing, 47, 957-965. https://doi.org/10.1109/78.752594</mixed-citation></ref><ref id="scirp.109204-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Amin, M.G., Wang, C. and Lindsey, A.R. (1999) Optimum Interference Excision in Spread Spectrum Communications Using Open-Loop Adaptive Filters. IEEE Transactions on Signal Processing, 47, 1966-1976. https://doi.org/10.1109/78.771045</mixed-citation></ref><ref id="scirp.109204-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Ma, W.J., Mao, W.-L. and Chang, F.R. (2005) Design of Adaptive All-Pass Based Notch Filter for Narrowband Anti-Jamming GPS System. 2005 International Symposium on Intelligent Signal Processing and Communication Systems, Hong Kong, 13-16 December 2005, 305-308. https://doi.org/10.1109/ISPACS.2005.1595407</mixed-citation></ref><ref id="scirp.109204-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Lv, Q. and Qin, H. (2018) A Novel Algorithm for Adaptive Notch Filter to Detect and Mitigate the CWI for GNSS Receivers. 2018 IEEE 3rd International Conference on Signal and Image Processing, Shenzhen, 13-15 July 2018, 444-451. https://doi.org/10.1109/SIPROCESS.2018.8600453</mixed-citation></ref><ref id="scirp.109204-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Chien, Y.R. (2013) Design of GPS Anti-Jamming Systems Using Adaptive Notch Filters. IEEE Systems Journal, 9, 451-460. https://doi.org/10.1109/JSYST.2013.2283753</mixed-citation></ref><ref id="scirp.109204-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Lv, Q. and Qin, H. (2020) General Method to Mitigate the Continuous Wave Interference and Narrowband Interference for GNSS Receivers. IET Radar, Sonar &amp; Navigation, 14, 1430-1435. https://doi.org/10.1049/iet-rsn.2020.0115</mixed-citation></ref><ref id="scirp.109204-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Ouyang, X. and Amin, M.G. (2001) Short-Time Fourier Transform Receiver for Nonstationary Interference Excision in Direct Sequence Spread Spectrum Communications. IEEE Transactions on Signal Processing, 49, 851-863. https://doi.org/10.1109/78.912929</mixed-citation></ref><ref id="scirp.109204-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Zhang, Y., Amin, M.G. and Lindsey, A.R. (2001) Anti-Jamming GPS Receivers Based on Bilinear Signal Distributions. 2001 MILCOM Proceedings Communications for Network-Centric Operations: Creating the Information Force, McLean, 28-31 October 2001, 1070-1074. https://doi.org/10.1109/MILCOM.2001.986006</mixed-citation></ref><ref id="scirp.109204-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Jiang, F.W., Xiao, Z., Yi, K.C. and Sun, Y.J. (2008) Adaptive Two-Stage Filter Bank Based Narrowband Interference Suppression in DSSS Systems. 2008 9th International Conference on Signal Processing, Beijing, 26-29 October 2008, 88-91. https://doi.org/10.1109/ICOSP.2008.4697075</mixed-citation></ref><ref id="scirp.109204-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Dovis, F. and Musumeci, L. (2011) Use of Wavelet Transforms for Interference Mitigation. 2011 International Conference on Localization and GNSS, Tampere, 29-30 June 2011, 116-121. https://doi.org/10.1109/ICL-GNSS.2011.5955275</mixed-citation></ref><ref id="scirp.109204-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Landry, R.J., Mouyon, P. and Leka&amp;iuml;m, D. (1998) Interference Mitigation in Spread Spectrum Systems by Wavelet Coefficients Thresholding. European Transactions on Telecommunications, 9, 191-202. https://doi.org/10.1002/ett.4460090209</mixed-citation></ref><ref id="scirp.109204-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Capozza, P.T., Holland, B.J., Hopkinson, T.M. and Landrau, R.L. (2000) A Single-Chip Narrow-Band Frequency-Domain Excisor for a Global Positioning System (GPS) Receiver. IEEE Journal of Solid-State Circuits, 35, 401-411.https://doi.org/10.1109/4.826823</mixed-citation></ref><ref id="scirp.109204-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Lotz, T. (2008) Adaptive Analog-to-Digital Conversion and Pre-Correlation Interference Mitigation Techniques in a GNSS Receiver. Diploma Thesis, Technical University of Kaiserslautern, Kaiserslautern.</mixed-citation></ref><ref id="scirp.109204-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Balaei, A.T. and Dempster, A.G. (2009) A Statistical Inference Technique for GPS Interference Detection. IEEE Transactions on Aerospace and Electronic Systems, 45, 1499-1511. https://doi.org/10.1109/TAES.2009.5310313</mixed-citation></ref><ref id="scirp.109204-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Same, M.H., Gleeton, G., Gandubert, G., Ivanov, P. and Landry, R. (2021) Multiple Narrowband Interferences Characterization, Detection and Mitigation Using Simplified Welch Algorithm and Notch Filtering. Applied Sciences, 11, Article No。 1331. https://doi.org/10.3390/app11031331</mixed-citation></ref><ref id="scirp.109204-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Choi, J.W. and Cho, N.I. (2002) Suppression of Narrow-Band Interference in DS-Spread Spectrum Systems Using Adaptive IIR Notch Filter. Signal Processing, 82, 2003-2013. https://doi.org/10.1016/S0165-1684(02)00385-7</mixed-citation></ref><ref id="scirp.109204-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Cho, N.I. and Lee, S.U. (1993) On the Adaptive Lattice Notch Filter for the Detection of Sinusoids. IEEE Transactions on Circuits and Systems II: Analog and Digital Signal Processing, 40, 405-416. https://doi.org/10.1109/82.238368</mixed-citation></ref><ref id="scirp.109204-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Cho, N.I., Choi, C.-H. and Lee, S.U. (1989) Adaptive Line Enhancement by Using an IIR Lattice Notch Filter. IEEE Transactions on Acoustics, Speech, and Signal Processing, 37, 585-589. https://doi.org/10.1109/29.17543</mixed-citation></ref><ref id="scirp.109204-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Borio, D., Camoriano, L. and Presti, L.L. (2008) Two-Pole and Multi-Pole Notch Filters: A Computationally Effective Solution for GNSS Interference Detection and Mitigation. IEEE Systems Journal, 2, 38-47. https://doi.org/10.1109/JSYST.2007.914780</mixed-citation></ref><ref id="scirp.109204-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Vijayan, R. and Poor, H.V. (1990) Nonlinear Techniques for Interference Suppression in Spread-Spectrum Systems. IEEE Transactions on Communications, 38, 1060-1065. https://doi.org/10.1109/26.57504</mixed-citation></ref><ref id="scirp.109204-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Wang, K.-J. and Yao, Y. (1999) New Nonlinear Algorithms for Narrowband Interference Suppression in CDMA Spread-Spectrum Systems. IEEE Journal on Selected Areas in Communications, 17, 2148-2153. https://doi.org/10.1109/49.814831</mixed-citation></ref><ref id="scirp.109204-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Wu, W.-R. and Yu, F.-F. (1996) New Nonlinear Algorithms for Estimating and Suppressing Narrowband Interference in DS Spread Spectrum Systems. IEEE Transactions on Communications, 44, 508-515. https://doi.org/10.1109/26.489097</mixed-citation></ref><ref id="scirp.109204-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Ferdjallah, M. and Barr, R.E. (1994) Adaptive Digital Notch Filter Design on the Unit Circle for the Removal of Powerline Noise from Biomedical Signals. IEEE Transactions on Biomedical Engineering, 41, 529-536. https://doi.org/10.1109/10.293240</mixed-citation></ref><ref id="scirp.109204-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Borio, D., Camoriano, L., Savasta, S. and Presti, L.L. (2008) Time-Frequency Excision for GNSS Applications. IEEE Systems Journal, 2, 27-37. https://doi.org/10.1109/JSYST.2007.914914</mixed-citation></ref><ref id="scirp.109204-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">El Gebali, A. and Landry Jr., R. (2020) Mitigation of Continuous Wave Narrow-Band Interference in QPSK Demodulation Using Adaptive IIR Notch Filter. American Journal of Signal Processing, Scientific &amp; Academic Publishing, 10, 10-18.</mixed-citation></ref><ref id="scirp.109204-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Punchalard, R. (2002) On Adaptive IIR Lattice Notch Filter Using a Robust Variable Step-Size for the Detection of Sinusoid. The 2002 8th International Conference on Communication Systems, Singapore, 28 November 2002, 800-804. https://doi.org/10.1109/ICCS.2002.1183240</mixed-citation></ref><ref id="scirp.109204-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Batra, A. and Zeidler, J.R. (2008) Narrowband Interference Mitigation in OFDM Systems. 2008 IEEE Military Communications Conference, San Diego, 16-19 November 2008, 1-7. https://doi.org/10.1109/MILCOM.2008.4753296</mixed-citation></ref><ref id="scirp.109204-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Ahn, J., Kim, J. and Kim, K. (2009) Adaptive Interference Suppression Methods Using Transform Domain Approach in an On-Board Filter Bank for Satellite Communications. TENCON 2009-2009 IEEE Region 10 Conference, Singapore, 23-26 January 2009, 1-5. https://doi.org/10.1109/TENCON.2009.5395990</mixed-citation></ref></ref-list></back></article>