A BDS Interference Suppression Technique Based on Linear Phase Adaptive IIR Notch Filters

The continuous wave interferences (CWIs) and the narrow-band interferences (NBIs) have significantly impacted the acquisition, tracking and positioning accuracy of Beidou Navigation Satellite System (BDS) receivers. As an interference suppression technology with a simple structure and low hardware cost, the adaptive infinite-duration impulse response (IIR) notch filter has been widely employed in the receivers to mitigate the CWIs and the NBIs. However, the nonlinear phase characteristics introduced by the IIR notch filters into the navigation receivers, may cause the distortion of navigation signals. It also leads to amplitude and phase distortion of the correlation peak in acquisition and loop tracking, which consequently brings on positioning errors in the measurement domain. This problem, however, has been ignored by many previous papers. Meanwhile, some other researchers came up with the idea of equalizers and all-pass networks compensating the distortion, which is also highly infeasible. Therefore, we propose a new method of an adaptive linear phase IIR notch filter with low hardware cost which is composed of three parts—the complex IIR notch filters, the IIR linear phase structure, and the adaptive and variable step-size algorithms. Applying this method, interference suppression and the correlation peak distortion compensation can be achieved with a modest increase hardware cost. This paper compares the performance of the new method with other IIR filters in both CWI and NBI scene and presents the effects of its parameters on the anti-jamming performance. Simulation results show that the proposed module has better anti-jamming performance in NBI scene and can compensate for the correlation peak distortion in the meantime.


Introduction
China has successfully launched the third-generation satellites of its BDS in the first step of establishing a vast positioning and navigation network with global coverage. However, the radio frequency interference severely influences BDS system on integrity, accuracy, continuity, and availability. Although the spread-spectrum processing gains vary from 30 to 40 dB, high interference signal ratio (ISR) interference with still cause significant positioning errors or paralysis on BDS receivers. Generally, there are several approaches to mitigate the interferences: time domain [1][2][3][4], frequency domain [5][6][7] and space domain [8][9][10].
As a time-domain interference suppression technology with a simple structure and low hardware cost, IIR notch filters have been widely employed in the navigation satellite receivers to mitigate the CWIs and the NBIs. Currently, many authors have paid considerable attention to such method. In [1], the authors used two-pole IIR notch filters to estimate the frequency of interferences and suppress them with another notch-depth adjustable IIR notch filter. Punchalard [3] implemented a modified gradient algorithm to improve the convergence speed of a first-order complex adaptive IIR notch filter. The complex IIR notch filter is mainly employed on Quadrature Phase Shift Keying (QPSK) systems. The coefficients are complex-valued, and the filter can simultaneously mitigate the signals on the in-phase (I) branch and the quadrature (Q) branch. In [2], a single order complex IIR notch filter was also proposed based on a least mean square (LMS) algorithm and a cascade structure to mitigate multi-CWIs. In [11], a new adaptive all-pass based notch filter was employed for narrowband/FM interference suppression and frequency estimation in GPS receivers. However, in its anti-jamming modules, no impact of nonlinear phase characteristic introduced by the IIR notch filters was considered, although such nonlinear phase characteristics may cause the amplitude and phase distortions of the correlation peak in the acquisition and tracking loops, and consequently lead to positioning errors in the measurement domain.
Some researchers have also analyzed the reasons for the distortion caused by the IIR notch filter. In [12], it was concluded that IIR notch filter would bring in a time delay error when it measures the pseudo-range. In [13], the authors used an adaptive linear interpolation filter and linear prediction filter in a direct sequence spread spectrum system and discussed their effects on the correlation peaks. They also indicated that a high side lobe level caused by the nonlinear characteristic could severely impact the acquisition loops. Chien [14] proposed a lattice IIR notch filter and analyzed its effects on the acquisition and tracking loop. Simulation results showed that such problem could be effectively omitted if the −3 dB bandwidth of the adaptive notch filters is adequately controlled.
To eliminate the nonlinear phase, Powell and Chau [15] proposed a linear phase IIR filter structure with last in first out modules (LIFO) and over-lap add modules, but they didn't apply an adaptive algorithm for interference suppression. In [16], a digital approximate linear phase notch filter was proposed by using two parallel all-pass filters. But the adaptive algorithm still cannot be used in such linear filter. Besides, others brought on some all-pass networks and equalizers to realize the linear phase of the IIR notch filter. In [17], the authors put forward an IIR notch filter with approximately linear phase with two all-pass sub-filters, whereas the computing complexity of coefficients was not low enough.
Considering that BDS B1 signal employs the QPSK modulation, in this paper, we combine the single order complex IIR notch filters and linear phase filter structure to accomplish the linear phase of the IIR notch filter and designed the variable step-size LMS algorithm for adaptively suppressing the interferences. This filter can be employed on BDS receivers to mitigate the CWI and NBI and to make a compensation for the correlation peak distortion.
The following parts of the paper are organized as follows: Section 2 describes the primary signal and system models, in which reasons of correlation peak distortion caused by traditional IIR notch filter are analyzed. Section 3 proposes the linear phase IIR notch filters structure, presenting the design of corresponding adaptive algorithm and variable step-size algorithm are also presented. Section 4 demonstrates the feasibility and availability of proposed method with a comparison between the filter performance of this new filter and others. Finally, the conclusions are drawn in Section 5.

Signal Models and Distortion Analysis
Before presenting the proposed method, firstly we need to introduce the signal model used in this paper. It is started with the BDS signal and interference model followed by the complex adaptive IIR notch filter. After that, we analyze the reason for the correlation peak distortion in theoretical derivation. Figure 1 shows the BDS anti-jamming receiver structures based on the complex IIR notch filters. The BDS B1 signal is composed of I-signal and Q-signal which are in quadrature with each other. It may suffer from CWIs and NBIs intentionally or unintentionally. The signal is received by the receiver antenna, after which it passes the band-pass filter, down converter, automatic gain control (AGC) and analog to digital converter (ADC). Finally, the sampled navigation intermediate-frequency signals r B (n) are obtained. Its complex expressions are as follows:

Signal Model
where B is the quantization bits. N sat represents the number of the visible navigation satellites. K is the number of interferences. s i (n) denotes the i-th navigation satellites signals. J m (n) is the m-th interference. N(n) represents the white Gaussian noise with zero mean and one-sided power spectral density N 0 . Concretely, the models of the navigation signals are as follows: where the index I represents the in-phase channel. Q represents the quadrature channel. P i I denotes the signal power. D i I (nT s ) and c i I (nT s ) represent navigation data and pseudo-random sequence (PRN), respectively. f L and f i dI are the transmitting frequency and Doppler frequency shift, respectively. θ i I is the carrier phase with uniform distribution. And in the practical conditions, P i I equals to P i Q . f i dI equals to f i dQ . The jamming signals can be modeled as follows: where Equation (3) is the CWIs model. P m I , f m I and θ m I denote the jam power, frequency, and carrier phase of the CWIs, respectively. Equation (4) represents the Auto-Regression (AR) model. In essence, the output of AR model can be considered as the narrowband interferences. Figure 2a shows the single order complex adaptive IIR notch filter in the satellite navigation receivers. r B (n), y BF (n) and g(n) are the input signals, output signals, and adaptive signals, respectively. D represents the delay module. The complex IIR filter's system function and its difference equation can be represented as follows: where X(z) and Y(z) are the z-transform of the input and output signals. r and α are the notch bandwidth and frequency parameters. They control the −3 dB bandwidth and the notch frequency. The IIR notch filter is stable if the r is smaller than 1. They can be written as: where BW denotes the −3 dB bandwidth. ω is the notch frequency. And the system function of IIR can be given by: where h i is the coefficient of the unit impulse response. δ(n) is the Dirac function.
To adjust the notch frequency parameter according to the interference frequency, we employ adaptive algorithms. In the adaptive filter theory [18], Wiener filter is a linear optimal filter based on the minimum mean square error (MSE). It extracts useful signals from stationary noise. In the navigation signal anti-jamming fields, the navigation signals, and white Gauss noise are considered to be the stationary noise. And the narrowband interferences are regarded as the useful signals. Figure 2b shows the equivalent Wiener Filter structure of Figure 2a. We can get the predicted interferenceĴ m (n) from the input signal r B (n). We also regard the received signal r B (n) as the desired signal in Figure 2b, so the output signal y BF (n) = r B (n) −Ĵ m (n) only contains the navigation signals and white Gauss noise. The adaptive algorithm adjusts its frequency parameters α to make the mean square minimum. The mean square can be express as follows: The solution of the frequency parameter is the solution of the Wiener-Hopf equation. However, it needs to calculate the inverse of the matrix which requires a tremendous hardware cost. So, the least mean square (LMS) algorithm has been adopted to simplify the computation. The adaptive algorithm is introduced in Section 3.3.
Without loss of generality, we assume that the interference has been entirely notched after the signals pass the adaptive IIR notch filter. Therefore, we obtain the output signals of filter y BF (n):

Distortion Analysis
The correlation peak will distort in acquisition and tracking modules. We analysis the distortion reason in this section.
After multiplying y BF (n) with the local carrier q(n) = exp{−j2π( f IF + f d )nT s } and the local PRN sequence c i (mT s −τ), the signals enter the integrate and dump modules. The output of the prompt correlator I p can be expressed as follows [14]: where G is the function of the AGC gain and signal-received power.τ is the time delay estimation. R(τ) is the correlation function. ∆τ is the delay time offset. N coh is the sample number of the integrate and dump filter for each integration. W(n) is the white Gaussian noise after processing. Equation (10) indicates that the prompt correlation result I p is the sum of the delayed correlation peaks weighed by the filter coefficients. Figure 3 shows the distorted correlation peaks with different −3 dB notch bandwidth parameters. Equation (10) shows that the correlation peaks will distort on amplitude and phase compared with the ideal one. There are two reasons for the distortion:

•
After the interference suppression, the partial power in the notch bandwidth of navigation signals is lost, which causes the decline of the correlation peak value, yet without affecting the symmetry. • From Equation (10), the nonlinear phase characteristic of IIR notch filters adds the ideal correlation peaks with the delayed correlation peaks together. And we can see the amplitude and phase distortion on the correlation results from Equation (5), where, h D = (r − 1)r D−1 (α * ) D (D > 1). Now we consider the absolute value of h D . That is due to the effect of absolute values of h D on the correlation amplitude. And from the equation above, the |h D | increases when r decreases, which means the delayed correlation peaks value increases and causes serious distortion. The simulation results in Figure 3 prove that when the parameter r decreases, the nonlinear phase characteristic causes more severe distortions.
The paper [4] shows that the for C/A code receivers, NBIs at frequencies away from band center degrade code accuracy more than that at band center. To illustrate this, we also compare the correlation peaks with different notch center frequencies-the simulation results shown in Figure 4 indicate that with the different jam offset frequency the distortion degrees also differ from each other.
To describe the distortion degree of the correlation peaks, we bring in the ∆-ratio test [19], which represents the symmetry of the correlation peaks. It can be expressed as Equation (12): where I e , I l and I p represent the coherent integration or non-coherent integration output of early, late, and prompt branches, respectively. ε is a constant value related to the correlator spacing. Generally, the correlator spacing between early and late correlator is half of the code chip; and the ε equals to 0.5. Under ideal circumstances, the ∆-ratio is zero. However, when the correlation peaks distort, the difference between I e and I l increases; and the value of ∆ also becomes larger. Assuming that the tracking loop has adopted the coherent dot power method in DLL, the formula is as follows [20]: where, the ρ cp is the code phase error. Then, the relationship between code phase error and ∆-ratio test is ρ cp = ∆ε/4. We can, therefore, conclude that in this case, the smaller ∆ is, the smaller the code error is.

IIR Notch Filter Based on the Linear Phase Structure
In [14], Chien concluded that if the −3 dB bandwidth of the IIR notch filter is less than about 30 KHz, the filters brings almost zero offset on both the acquisition and tracking loops. However, the power amplifiers, up-converters, and other nonlinear components in the jammer widen the CWIs' bandwidth. And the signal processing, such as signal truncation, may also result in a wider CWIs bandwidth. What's more, for improving the application range, the −3 dB bandwidth of the IIR filter should not be limited under 30 KHz. Therefore, when the −3 dB bandwidth gets wider, the effects of the nonlinear phase characteristic on the acquisition and tracking loop cannot be neglected.
In this section, we combine the single order complex IIR notch filter with the linear phase IIR structure in [15] to compensate the nonlinear effects. The filters can apply to the BDS receivers to reject the CWIs and NBIs and we also design the adaptive algorithm and variable step-size algorithm to the corresponding linear structure.

Linear Phase IIR Notch Filter Structure
In this section, we propose the complex linear phase IIR notch filter, which is based on the linear phase IIR structure, single order complex IIR notch filters, and the corresponding adaptive algorithm and variable-step algorithm. The structure presents as follows ( Figure 5):  Figure 5 shows that this module is composed of the LIFO modules, an overlap-add module, and an adaptive module. The LIFO modules are composed of the L registers, which can make the L signal samples reversed in the time domain. The overlapped modules are including a delay module, an adder, two complex IIR filters and three switches. Its primary task is to realize the sectioned convolution. The input of the two IIR filters is controlled by the switch-K, which change its direction every L samples. And the output signals are composed of a superposition of the two single order complex IIR notch filters output. Figure 6 shows the signal sequence chart to explain how the overlap-add module works to achieve the sectioned convolution. Its theoretical analysis is demonstrated as follows. The input of the sectioned convolution can be expressed as follows: The output signal v(n) can be expressed as: and the v k (n) can be shown by: The k-th part and the (k + 1)-th part of the convolution result can be expressed as follows: We can see that the convolution results v(n) in Equation (16) are composed of two parts. The first part is from the second item of the Equation (17), and the second part is from the Equation (18)'s first item. In summary, the convolution results of the whole sequence can be calculated by the sum of the sectioned convolution results.
And according to the adaptive algorithm mentioned below, the input signal of the adaptive modules are the gradient signal of the output signal of the IIR filter, which can be represented as Equation (24). And the output of the adaptive modules adjusts the notch frequency parameter r and the step parameter µ of the three notch filters simultaneously. In summary, the signal processing procedure is as follows:

•
Step 1: The navigation signals enter the LIFO modules to make the signal time reversed, which indicates u(n) = x(−n); • Step 2: u(n) passes the overlap-add modules for sectioned convolution; • Step 3: The signal v(n) passes the second LIFO modules. After that, the input signals x(n) are equivalent to pass a non-causal filter, whose unit impulse response is expressed as h(−n); • Step 4: The IIR3 outputs the final signals which is also the output of the linear phase IIR notch filter.

•
Step 5: And consider the output signal y(n) and the adaptive signals g(n) of the last level notch filter as the adaptive modules input. It adjusts the notch frequency and −3 dB bandwidth parameters.
We prove that the proposed filter has a linear phase characteristic from the time domain and frequency domain. In the time domain, the proposed modules, in essence, are composed of two cascaded single order notch filters. There unit impulse responses are h(−n) and h(n), respectively. The cascade of the two filters constitutes a second order linear phase IIR notch filter. Its unit impulse response is an even symmetry sequence. The system function can be expressed as: h (n) = h(−n) * h(n). Its coefficients are even-symmetric. The symbol * represents the convolution. The sufficient and necessary conditions for the filter's linear phase characteristic is that the coefficients of the time-domain system function are even or odd sequence. As a conclusion, the proposed structure in Figure 6 has the linear phase characteristic.
From the frequency domain, the proofs are expressed as follows: where U(z), V(z) and W(z) are the Z-transform of the u(n), v(n) and w(n), respectively in Figure 5. Equation (22) shows that the output signals are relevant to the magnitude characteristic of the IIR filter without non-linear phase characteristic. Therefore, this structure can be employed to mitigate the correlation peaks distortion caused by the nonlinear phase characteristic. In essence, the proposed filter is a second order IIR filter. It has almost the equal anti-jamming performance comparing with the other traditional notch filters. Generally, in some cases which need strict linear phase, such as high-accuracy survey and mapping, the all-pass networks and equalizers will be added to compensate the nonlinear phase effect. But the design of all-pass network and equalizer is more challenging to design an IIR notch filter. The structure of the linear phase IIR notch filter only needs three single order complex IIR filters to meet the linear phase characteristic needs.
The system complexity and the hardware cost of the proposed filter are discussed as follows. To simplify the calculation, we adopt the following numbers to calculate for every L samples. Concerning system complexity, samples in IIR notch filter need 2L multiplication and 2L addition. L additions are needed in over-add modules. 4L addition, 4L + 3 multiplication, and 2 subtractions are also needed in adaptive algorithm modules, so the proposed filter's system complexity is O(14L + 5), that is O(L). Regarding the hardware cost, only 6L + 3 registers, 2L + 7 summing units, 2L + 10 multiplying units and two dividers are needed to meet the filter's need. Consequently, both system complexity and hardware cost are low enough for use in low-cost receivers.

Corresponding Adaptive and Variable
Step-Size Algorithm Figure 5 shows that the linear phase structure adapted a block processing method to realize the time-reversal and sectioned-convolution. To maintain the system coherence, we rewrite the adaptive algorithm [14] to the block form. The adaptive algorithm adjusts the α every K samples. K is the block length. The adaptive algorithm formulations are inferred as follows, it's the N-th renewal process: where Y BF (N) represents the output signals column vectors with length K. w(NK) is the input signal of the IIR3 filter in Figure 5. Y BF (N)'s derivative functions for α are expressed as follows. It is regarded as the adaptive signals: According to the LMS algorithm, the instantaneous gradient can be expressed as follows: the iteration formula of the parameter α is as follows: where the superscript H represents the complex conjugate transpose. µ represents the update step-size. Considering G(N) and Y BF (N) as the input of the adaptive modules, the filters updated α every K samples. Figure 5 shows that the three IIR notch filters need the same system function when handling the same block signals. This linear phase structure cannot make the system function ultimately the same. But the simulation results indicate that if the step parameter µ is small enough, the effects on the filter performance can be ignored.
To accelerate the convergence rate, we employ the variable step-size algorithm. The Equation (26) is rewritten as follows: where s is the constant step-size for adjusting the convergence rate. ϕ(N) is an instantaneous power estimation of G(N), whose initial value affects the initial value of the first step-size. And ρ (0 < ρ < 1) is the forgetting factor which presents the correlation between power estimation of G(N) and G(N − 1). After the iteration, the α(N) convergence to a fixed value. That means the filter has entered the steady state. What we are concerned about in steady state is the mean square error (the steady state error) shown in Equation (29). It affects the convergence performance and the stability of the filter:

Analysis of Filter Performance
As mentioned above, the block length L, the notch bandwidth parameter r and the block length K of the adaptive module are the main factors which affect the performance of the IIR notch filter. In this section, therefore, they are examined in detail.
The block length L the number of the LIFO registers, which means the input signals are divided into L samples per block. It affects the over-added performance on approximating to the ideal transfer function. The more samples per block, the better each section's transfer function approximates H(z). We assume that L is larger than impulse response length of the physical realization. Then when |i| > L, the notch filter coefficient h i is small enough to be neglected.
The −3 dB notch bandwidth parameter r can significantly affect the correlation peaks' symmetry in the traditional IIR filters. But it does little influence on it in the proposed IIR structure. It mainly changes the value of the correlation peaks. Two reasons mainly cause the distortion of the correlation peaks: the power loss of the normal signals and the linear characteristic. As a result, when we need to increase the −3 dB bandwidth of the filter, the power loss also increases, leading to the decrease of the peaks value.
The block length K of the adaptive modules also affect the filter performance in the aspects of convergence rate and steady-state error. With the block length increasing and without the variable step-size algorithm, the filter gives a slower convergence rate but less steady-state errors. However, while applying the variable step-size algorithm, the filter converges to smaller steady-state errors than the fixed step-size. Therefore, to guarantee no conflict between the update of the filter coefficients and the sectioned convolution, the relationship between block length K and LIFO length L should be as follows: but the simulation shows that K = L is the best choice. The larger K is, the more complex the IIR notch filter is.

Simulation and Verification
This section divides into two parts: one is comparing the proposed filter with other filters under CWIs and NBIs scene, respectively; the other is comparing their performance with different parameters. First, the simulation conditions present as follows. Figure 7 demonstrates the BDS anti-jamming simulation and verification platform, which is composed of jammers, BDS signal acquisition system, anti-jamming modules, and software receiver.  Table 1 lists the verification conditions on the system parameters.

Comparison with Other Filters
By applying different filters, including adaptive IIR notch filter with LMS algorithm [1], frequency excision method [21], Chien's filter [14] and the proposed filter, we compare their anti-jamming capabilities and the correlation peak symmetry and set up two application scenarios. One is that the BDS receiver suffers from the CWI, where we set the −3 dB bandwidth of the filters to the same and narrow level. The other is suffering from the NBI, where we also set the −3 dB bandwidth to the same value but a wider level. Figure 8 shows the acquisition chart results of the four filters in the CWI scenario. The correlation peak values and ∆ ratio test results are listed in Table 2: From Table 2 and Figure 8, the results show that in the CWI scene, the correlation peak values are almost the same except for the frequency excision. The ∆ ratio tests also have low values. The proposed filter has the same performance in this scenario. However, in the NBI scenario, we can find that all the peak values decrease. Such a phenomenon occurs because the filter mitigates more navigation signal with wider −3 dB bandwidth. Compared with other two IIR filters, the proposed filter has higher correlation peak and lower ∆ ratio value. In the given verification conditions, from the ρ cp results, the proposed method reduces the code phase error from 22.11 m to 0.16 m without any other errors. In conclusion, the proposed filter is found to have better interference suppression performance and effectively mitigates the correlation peak distortion. Figures 9 and 10 illustrate the correlation peaks of the four filters with CWI and NBI, respectively.
The frequency excision has almost equal values on the correlation peaks and ∆ ratio value in different scenarios. Figure 9 shows that the proposed filter has the same anti-jamming performance as other filters. Figure 10 shows that in NBI scene, the adaptive IIR notch filter and Chien's filter could have high side-lobes alongside the main peak without linear characteristics. The nonlinear phase characteristic also causes the asymmetry of the correlation peak. However, the proposed filter has little effects on the correlation peaks.
In conclusion, in the CWI scenario, the IIR notch filters have almost the same anti-jamming performance. However, in the NBI scene, the proposed filter performs better because of its distortionless correlation peak and lower side-lobes.

Anti-jamming Performance with Different Filter Parameters
Besides comparing the correlation peak value and the convergence performance, respectively, we also analyze how the block length L and the −3 dB bandwidth parameter r affect the proposed filter performance itself. Figure 11 shows the influence the block length L has on the correlation peak values: Figure 11. Comparison of the acquisition result with different L.
The result shows the correlation peak values produce nearly the same results with different block length L. When L is 60, the correlation peak value is 1.052 × 10 12 , and when L is 20, the correlation peak is 1.041 × 10 12 . The reason (for what) is that when the length L rises to a particular value, the second term on the right side of the Equation (26), Y H (N)G(N)/L is equal to the same constant value. In each iteration, we can get α L=20 (n) = α L=30 (n) = . . . = α L=60 (n). That means the filters with different L converge to the same notch frequency in the steady state. With the same −3 dB notch bandwidth and the same notch frequency in the steady state, the block length L doesn't significantly affect the correlation values.
Whats' more, Equation (27) shows that the length K (K = L) affects the convergence rate of parameter α(n), due to which, we analyze the convergence rate with the same variable step-parameter µ(n). Figures 12 and 13 show the convergence performance under the same initial conditions of the parameter α(n) and µ(n). Figure 12 indicates that with the increasing of parameter L, it costs more time to make the notch frequency parameter convergence close to steady-state. The fastest convergence rate is 0.2 ms as the L equals to 20. It is three times faster than it is when L equals to 60. The block length 60 is also three times larger than L = 20. On the premise that the notch frequency parameter α(n) and step parameter µ(n) start from identical initial states. We, therefore, conclude that from each iteration, we can get that α L=20 (n) = α L=30 (n) = . . . = α L=60 (n). However, the iteration time has a linear relation with block length, so on the timescale, the convergence rate is slower with the increasing length L.  While Figure 12 indicates the relationship between convergence rate and the block length L. Figure 13 shows the stability of the filter. It can be noticed that the Normalized MSE (L = 100) has more unstable steady state errors because of its dependence on the variable step-size. Equation (27) and Equation (28) show that the µ(n) has the same value under the same G(n). Although the convergence rate becomes higher with the increase of L, the stability of the steady-state error decreases, which means a more stable state of the filter. As a result, when block length L is given, both the convergence rate and the stability of the filter should be taken into consideration.
Meanwhile, we also compare the convergence rate with different parameter r and present the results in Figure 14. The presented results indicate that the narrower the −3 dB bandwidth is, the faster the convergence rate is. This can be explained by the fact that when the filters do not notch the interferences, the ϕ(n) in Equation (28) remains the same-they have the same step-size to converge to steady-state. However, the interferences are firstly notched by the wider −3 dB bandwidth filter; and consequently, the MSE decreases, which makes the step-size smaller and the filter gets a slower convergence rate. Also, we can conclude that the convergence rate is sensitive to the parameter r.

Conclusions
This paper mainly focuses on the non-linear phase characteristics caused by the traditional IIR notch filter after interference suppression and proposes an adaptive linear phase second order IIR notch filter with a small hardware cost. Meanwhile, the corresponding variable step-size adaptive algorithm is designed to improve the convergence performance, which can help the BDS anti-jamming receivers to reject CWIs and NBIs. According to the simulation results, the proposed filter presents almost equal anti-jamming abilities in the CWI scenario compared with other methods, including the adaptive IIR notch filter, frequency excision, as well as Chien's filter. Contrastively, it performs better on anti-jamming in the NBI scenario, with the capability of mitigating the acquisition distortion and side-lobes at the same time. The block length L and the −3 dB bandwidth affect the performance as well. Besides, we also analyze their impacts on filter performance: As L decreases, the correlation peak values increase slightly; the convergence rate speeds up; the hardware cost reduces. Besides, the convergence rate also increases with narrower bandwidth. However, the simulation results also indicate that the length L should be larger than a certain value, or the filter cannot work normally. Therefore, when designing the linear phase IIR notch filters, researchers should consider the block length L and the −3 dB bandwidth as well.