Next Article in Journal
Deforestation as the Prominent Driver of the Intensifying Wildfire in Cambodia, Revealed through Geospatial Analysis
Next Article in Special Issue
Anti-Jamming GNSS Antenna Array Receiver with Reduced Phase Distortions Using a Robust Phase Compensation Technique
Previous Article in Journal
Spatiotemporal Features of the Surface Urban Heat Island of Bacău City (Romania) during the Warm Season and Local Trends of LST Imposed by Land Use Changes during the Last 20 Years
Previous Article in Special Issue
Evaluation and Analysis of the Accuracy of Open-Source Software and Online Services for PPP Processing in Static Mode
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

Blind Adaptive Beamforming for a Global Navigation Satellite System Array Receiver Based on Direction Lock Loop

College of Electronic Science and Technology, National University of Defense Technology, Changsha 410073, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(13), 3387; https://doi.org/10.3390/rs15133387
Submission received: 14 May 2023 / Revised: 27 June 2023 / Accepted: 28 June 2023 / Published: 3 July 2023
(This article belongs to the Special Issue Advancement of GNSS Signal Processing and Navigation)

Abstract

:
The adaptive beamforming algorithm can realize interference suppression and navigation signal enhancement, and has been widely used. However, achieving high-precision real-time estimation of the direction of arrival (DOA) parameters of navigation signals in strong-interference scenarios with low complexity is still a challenge. In this paper, a blind adaptive beamforming algorithm for a Global Navigation Satellite System (GNSS) array receiver based on direction lock loop is proposed without using the prior information of the DOA parameter. The direction lock loop is used for DOA tracking and estimation after interference suppression, which uses the spatial correlation of the array beam pattern to construct a closed direction-tracking loop. The DOA estimation value is adjusted in real time based on the loop errors. A blind beamformer is constructed using the DOA estimation results to provide gain by forming a beam in the satellite direction. This method improves the accuracy and dynamic adaptability of DOA estimation while significantly reducing the computational complexity. The theoretical analysis and simulation results verify the effectiveness of the proposed algorithm.

1. Introduction

The Global Navigation Satellite System (GNSS) can provide high-precision positioning, navigation and timing services [1]. However, the satellite navigation signal is very weak and buried in noise, which is easily affected by intentional and unintentional interference [2,3,4]. The array anti-jamming method can improve the interference suppression capabilities by using an adaptive antenna array, and has become one of the most useful anti-jamming methods [5,6].
The blind interference suppression algorithms, such as the power inversion (PI) algorithm [7,8], can suppress the interference without using the prior knowledge of the direction of arrival (DOA) parameter of the satellite signal. However, these methods are unable to enhance the navigation signal energy, and may even reduce the signal-to-noise ratio [9].
The beamforming method [10] forms a null in the direction of the interference and points the beam in the direction of the navigation signal, thereby achieving interference suppression while increasing the energy of the useful signal, which can effectively improve the weak signal reception capability. The Minimum Variance Distortionless Response (MVDR) method [11] is the most commonly used beamforming method, which uses the prior DOA information of the incident signal to construct the beamformer. This method generally needs to obtain the real-time attitude information of the array antenna by means of inertial navigation; however, most array antennas cannot provide this.
Blind beamforming algorithms can achieve similar effects to MVDR by using the DOA estimation method without prior information [12,13]. The DOA parameters of the incident signal are estimated in real time by the array signal processing method. The attitude information of the array antenna is not required as a priori information. The commonly used DOA parameter estimation methods include subspace-based methods such as MUSIC [14] and ESPRIT [15], the maximum likelihood estimator (MLE) method [16], the parse recovery method [17], etc. However, the DOA estimation performance of these methods in dynamic scenes will decrease because the target motion will degrade the spectral resolution of the array [18]. At the same time, the computational complexity of these methods is high, which is not applicable for the actual navigation receiver.
Direction lock loop (DiLL) [19,20,21] can achieve high-precision DOA parameter tracking, which is very suitable for a direct sequence code division multiple access (DS-CDMA) system. DiLL uses the left-shifted and right-shifted spatial array response vectors to construct a closed loop for direction tracking, which can reduce the computational complexity significantly [22].
In Ref. [23], DiLL is used for the DOA estimation of mobile target, which can achieve a high-resolution DOA estimate signal. In Ref. [24], a DOA tracking method based on DiLL with tracking bias compensation for GNSS receivers was proposed, in which the tracking bias was effectively eliminated by estimating and compensating the error in real time. In Ref. [22], the rotate-to-zero DiLL (RZ-DiLL) method was proposed, which could effectively improve the DOA parameter estimation accuracy of large-angle incident signals. In Ref. [25], a DiLL method for a three-element L-shaped array was proposed to achieve the fast acquisition and tracking of two-dimensional DOA parameters, which significantly simplified the antenna structure and eliminated the dead-zone effect. However, most of the above methods are only suitable for non-interference scenarios. The realization of DOA tracking in strong interference scenarios is in urgent need of research.
In this paper, a blind beamforming method for a GNSS array receiver based on DiLL is proposed. Firstly, the subspace projection method is used to suppress the interference signal. Secondly, the DOA parameter of the navigation signal is tracked and estimated by using the DiLL method after interference suppression. Thirdly, the estimated DOA parameters are used to construct the beamformer, which can achieve interference suppression and signal enhancement. The contribution and novelty are reflected in the following aspects:
  • The DOA parameter of the GNSS signal can be estimated accurately in strong interference scenarios. It solves the problem that it is difficult to achieve both interference suppression and DOA estimation with traditional methods.
  • The proposed method can realize DOA parameter estimation with low computational complexity, and thus stable DOA tracking in dynamic scenes. It is very suitable for the GNSS receiver, which needs real-time processing.
  • The blind adaptive beamformer is constructed using the estimated DOA parameter, which can realize interference suppression and navigation signal enhancement. It can achieve almost the same performance as the MVDR method.
The remainder of this paper is organized as follows: Section 2 presents the signal model. Section 3 presents the proposed blind adaptive beamforming method based on DiLL. Section 4 gives the simulation results, and Section 5 concludes this paper.

2. Signal Model

In order to simplify the analysis without loss of generality, under the far-field condition, it is assumed that a single navigation signal and Q interference signal occur on the uniform line array (ULA) simultaneously. The ULA has N antenna elements. The array configuration of the ULA is shown in Figure 1.
The coupled signal received by the GNSS receiver can be expressed as
x ( t ) = a ( θ 0 ) s ( t τ 0 ) + q = 1 Q a ( θ q ) j q ( t ) + n ( t )
where x ( t ) = [ x 1 ( t ) , x 2 ( t ) , , x N ( t ) ] T stands for array data vector. s ( t ) represents the baseband GNSS signal, τ 0 and θ 0 are the time delay and DOA parameter of GNSS signal, respectively. j q denotes the q th interference signal component with DOA parameter θ q . n ( t ) represents the additive white Gaussian noise vector, and a ( θ ) is the array guide vector, which can be expressed as
a ( θ ) = 1 , e j 2 π d sin ( θ ) / λ , , e j 2 π d ( N 1 ) sin ( θ ) / λ T
where θ is the signal incident angle, λ is the wavelength, and d = λ / 2 is the array spacing.
The baseband GNSS signal s ( t - τ 0 ) can be expressed as
s ( t τ 0 ) = A D ( t τ 0 ) c ( t τ 0 ) e j 2 π f d t + φ
where D ( ) indicates the message bit. In this paper, the correlation accumulation time is all within a bit range, so the navigation message is assumed to be constant, i.e., 1. c ( ) represents a pseudo-random noise (PRN) code. A , f d , and φ represent navigation signal amplitude, Doppler frequency, and carrier phase, respectively.
By constructing an array beamformer, the array can be weighted to achieve an effective suppression of interference. The output after array weighting is
y ( t ) = w H x = w H a ( θ 0 ) s ( t τ 0 ) + q = 1 Q w H a ( θ q ) j q ( t ) + w H n ( t ) w H a ( θ 0 ) s ( t τ 0 ) + w H n ( t )
where w stands for array weight vector.

3. Blind Adaptive Beamforming Based on Direction Lock Loop

3.1. Interference Suppression

Since the power of the interference signal is much higher than that of the navigation signal, the covariance matrix of the array signal is determined mainly by the interference and noise terms.
R = E x x H = n = 1 N λ n u n u n H   = n = 1 Q ( λ n σ e 2 ) u n u n H + σ e 2 I N × N
where λ n and u n ( n = 1 , 2 , , N ) denote the eigenvalues and the corresponding eigenvectors of R . σ e 2 represents the power of the noise. The eigenvalues are arranged in descending order:
λ 1 λ 2 λ Q λ Q + 1 = = λ N = σ e 2
The inverse matrix R 1 can be calculated as
R 1 = n = 1 N 1 λ n u n u n H   = n = 1 Q 1 λ n u n u n H   + 1 σ e 2 n = Q + 1 N u n u n H = 1 σ e 2 I N × N n = 1 Q λ n σ e 2 λ n u n u n H
Assuming the interference power is much larger than the noise power, the orthogonal complementary space U I can be approximated by the inverse of the covariance matrix R [26], i.e.,
R - 1 1 σ e 2 U I
In Equation (8), the coefficient 1 / σ e 2 is a constant related to the noise power, which can be omitted. There is no need to perform eigenvalue decomposition, which can significantly reduce the amount of computation. By projecting the array signal into the orthogonal complementary space, the interference can be mitigated. The array signal after the interference suppression can be expressed as
x ^ ( t ) R 1 x  

3.2. DOA Tracking Based on DiLL

Since the navigation signal is very weak and buried under the noise, the correlation accumulation with the locally copied pseudo-code signal is required to improve the signal-to-noise ratio. The correlation results can be modeled as
r = E x ^ ( t ) c L * ( t τ ^ )   = R 1 a ( θ 0 ) r 0 ( Δ τ )
where r 0 ( Δ τ ) represents the correlation summation value between the received navigation signal and the local replica pseudo-code. c L ( t τ ^ ) is the local replica PRN code, and Δ τ = τ 0 τ ^ represents the estimation error of the time delay. The carrier is considered to be completely wiped off.
From Equation (10), it can be seen that the correlation accumulation vector is closely related to the incident angle of the received signal. The array direction diagram is constructed by using this vector as the weight vector. Then, the gain is maximum in the direction of the signal incident.
The normalized array pattern of the ULA is shown in Figure 2. The normalized spatial correlation function can be defined as
R ( θ , θ + Δ θ ) = 1 N a H ( θ + Δ θ ) a ( θ ) = 1 N n = 1 N e j π ( n 1 ) sin θ sin ( θ + Δ θ )
where ( ) H denotes the Hermitian transform and Δ θ is the direction bias.
The array pattern is symmetrical around the central direction, so
R ( θ , θ + Δ θ ) R ( θ , θ Δ θ )
The direction can be tracked as the same principle of the code DLL loop. Assuming that the estimated value of DOA at time k is θ ^ k , the array guide vectors a ( θ ^ k Δ θ ) and a ( θ ^ k + Δ θ ) of the left and right offset Δ θ of the current observed DOA are constructed. After correlating with the estimated signal vector, the spatial correlation functions can be obtained as
z R ( k ) = 1 N a H θ ^ k + Δ θ r = 1 N r 0 ( Δ τ ) a H θ ^ k + Δ θ R 1 a ( θ 0 ) = r 0 ( Δ τ ) R θ 0 , θ ^ k + Δ θ
z L ( k ) = 1 N a H θ ^ k Δ θ r = 1 N r 0 ( Δ τ ) a H θ ^ k Δ θ R 1 a ( θ 0 ) = r 0 ( Δ τ ) R θ 0 , θ ^ k Δ θ
where R ( θ 1 , θ 2 ) = a H ( θ 2 ) R 1 a ( θ 1 ) / N represents the normalized spatial correlation function after interference suppression.
The discrimination error can be expressed as
ε ( k ) = f z R ( k ) , z L ( k )
where f ( ) denotes the discriminator function. The coherent DiLL and non-coherent DiLL are the two mainly used direction discriminating functions, which can be expressed as
f c o h R 1 , R 1 = Re R 1 R 2
f n o n c o h R 1 , R 2 = R 1 2 R 2 2 R 1 2 + R 2 2
The S-curves of the coherent DiLL and non-coherent DiLL are shown in Figure 3, in which N = 4 ,   Δ θ = 5 and θ = 0 . As can be seen from Figure 3, the differential curve shows an approximate linear change near 0 degrees, so the DOA estimation error can be corrected by the discriminator. It also shown that the convergence of the non-coherent DiLL is larger than that of the coherent DiLL.
After loop filtering, the direction discrimination error is sent to the beamforming device to update the DOA parameter value at the next time. Then, the estimated DOA value at time k + 1 is
θ ^ k + 1 = θ ^ k + K 0 ε ( k ) h ( k )
where h ( k ) stands for loop filter, represents a convolution operation. If θ ^ k < θ , then ε ( k ) > 0 , and the estimated value of DOA at the next time will increase; otherwise, the DOA parameter at the next time will decrease. The accurate estimation of DOA parameters can be achieved through continuous iteration.

3.3. Beamformer Construction

The array vector a ( θ ^ k ) of the signal direction is constructed according to the DOA parameter estimated by DiLL in real time. The output power of the array is minimized, while the response of the signal direction is guaranteed with the unit array:
min w   { w H R w }   s . t .   w H a ( θ ^ k ) = 1
Using the linear constraint minimum variance (LCMV) method to solve the above equation, the optimal weight can be obtained as
w o p t = R 1 a ( θ ^ k ) a H ( θ ^ k ) R 1 a ( θ ^ k )
The array-weighted output can be calculated as
y o u t ( t ) = w o p t H x

3.4. The Procedure of the Proposed Method

The diagram of the proposed algorithm is shown in Figure 4 with the following steps.
(1)
The inverse matrix of the covariance matrix of the array signal was solved and multiplied by the array signal to achieve interference suppression.
(2)
The array signal after interference suppression was correlated with the local copied pseudo-code signal to obtain the correlation accumulation vector.
(3)
The DiLL method was used to realize the stable tracking processing of the DOA parameters of navigation signals. The spatial correlation values corresponding to the Δ θ bias besides the current DOA estimation were constructed, and the revised DOA estimate was obtained after processing and filtering by the discriminator.
(4)
The beamformer was constructed according to the estimated DOA value, the array was weighted, and the weighted output signal was used for follow-up tracking processing.

4. Simulation Results

4.1. DOA Tracking Performance Analysis

The tracking accuracy of DOA is generally characterized by mean square error (MSE), which can be defined as
MSE = 1 I i = 1 I θ ( i ) θ ^ ( i ) 2
where θ ( i ) and θ ^ ( i ) are the real DOA value and estimation value at the i th time, I is the number of processing times.
In this section, the performance of the DOA estimation is evaluated by Monte Carlo simulation. The Monte Carlo simulation was performed 1000 times; θ = 0 and Δ θ = 5 . Figure 5 shows the MSE results of the DOA estimation under different carrier-to-noise ratios ( C / N 0 ) and different array elements. The performances of the coherent DiLL method, non-coherent DiLL method, and MUSIC method are compared and analyzed. As can be seen from Figure 5, the accuracy of DOA estimation is higher with the increase in the C / N 0 . Under the same C / N 0 condition, the more array elements there are, the higher the estimation accuracy. When the number of array elements is small, the DiLL algorithm can achieve better performance than MUSIC, mainly because when the number of array elements is small, the resolution of the MUSIC algorithm is significantly reduced, while the loop tracking method can still achieve high-precision tracking results. With the increase in the number of array elements, the performance of the MUSIC algorithm is improved and the DiLL algorithm is also effectively improved. Therefore, when the number of array elements is small, the DiLL algorithm has better advantages. It can also be seen that the performance of the coherent DiLL method is better than that of the non-coherent DiLL algorithm.

4.2. Dynamic Adaptability Analysis

Through estimating the tracking error and the error compensation, DiLL realizes DOA estimation in real time, and thus stable DOA tracking in dynamic scenes. Typical DOA dynamic change scenarios include fixed DOA parameters, uniform acceleration change, sinusoidal change, and step change. Figure 6, Figure 7, Figure 8 and Figure 9 give the DiLL tracking performance in these scenarios, respectively. In the simulation, the angle interval is 5 degrees, the number of array elements is eight, and the carrier-to-noise ratio parameter is 40 dBHz.
Figure 6 shows the tracking results under fixed DOA conditions. In this scenario, the DOA parameter is a fixed value of 10°. The tracking results under the condition of uniform acceleration of DOA are given in Figure 7. The change rate of DOA is 20°/s2, and the angle range is −10–0°. Figure 8 shows the tracking results under the condition of sinusoidal variation in DOA. The variation period is 1 s and the variation range is −10–10°. Figure 9 shows the tracking results under the condition of DOA step change, and the DOA step is 10°. When there is a step change in the DOA, the tracking loop needs to converge to a new parameter. Therefore, the DOA error may increase in several loop update cycles. The DOA tracking errors in the above four scenarios are 0.053, 0.054, 0.06, and 0.11, respectively. The simulation results show that the DiLL method can adapt well to the DOA tracking problem in different dynamic scenes.

4.3. GNSS Signal Tracking Performance Analysis

The simulation module, which consists of three parts including the signal generation module, array signal processing module, and GNSS tracking module, was used to verify the efficacy of the algorithm. Figure 10 illustrates the simulation module. The signal generation module is used for simulating the array receiver signals. Specifically, the GNSS signal vector, interference signal vector, and noise signal vector were simulated individually, and then combined to form the received signal. On the other hand, the array signal processing module were mainly used for interference suppression, DOA estimation, and beamforming. After array weighting, the GNSS software-defined receiver [27] was used for GNSS signal tracking.
In the simulation, one satellite navigation signal and two interference signals were simulated, a uniform line array with eight elements was used, the sampling rate was 6 MHz, and the correlation integration time was 1 ms. The navigation signal was GPS L1C/A signal, the pseudo-code frequency was 1.023 MHz, C/N0 was 45 dBHz, and the DOA parameter was 30°. The interference signal was broadband Gaussian noise interference, the bandwidth was 2.046 MHz, the incidence direction was −40° and 55°, respectively, and the interference-to-noise-ratios were 50 dB and 55 dB, respectively. The simulation duration was 3 s, which was divided into three stages: the interference was turned off in the first second; the interference was turned on in the second second; and the C/N0 value was turned down to 30 dBHz in the third second.
In the software receiver, the performance of the proposed method was compared with the PI and MVDR methods. The PI algorithm is one of the most useful blind adaptive beamformers, and forms nulls in the direction of the jammers. The cost function of PI is as follows:
min   w P I w P I H R w P I   s . t .   w P I H δ = 1
where δ = [ 1 , 0 , , 0 ] T is an N × 1 vector.
The PI beamformer can be expressed as
w P I = R 1 δ δ H R 1 δ
The MVDR cost function is given as
min w M V D R w M V D R H R w M V D R   s . t .   w M V D R H a θ = 1
The solution of Equation (25) is given by
w M V D R = R 1 a θ a H θ R 1 a θ
The simulation results are shown in Figure 11, Figure 12, Figure 13 and Figure 14. Figure 11 gives the normalized array pattern results of the proposed method and the PI method with and without the influence of the interference. It shows that the PI method can achieve an effective suppression of interference, but cannot achieve signal direction enhancement, and may even lead to signal energy loss. In the proposed method, the DOA parameter can be estimated accurately using DiLL, so a blind adaptive beamformer can be constructed to point the beam in the direction of the navigation signals, which can enhance the signal energy, even in the interference scenario.
Figure 12 shows the punctual branch correlation value using the proposed method and the PI and MVDR methods. From the figure, it can be seen that the proposed method and MVDR can both achieve gain in the signal direction, and significantly increase signal energy compared with the PI algorithm. When the interference is turned on, the PI algorithm will cause signal energy loss while suppressing the interference as shown in Figure 12. When the C/N0 drops to 30 dBHz, the navigation signal energy after the PI algorithm is too low to cause the signal tracking loop to lose lock. However, both the proposed method and the MVDR method can still guarantee the stability tracking of the signal through the array gain, which significantly improves the tracking performance of the navigation receiver under weak signal conditions.
Figure 13 shows the pseudo-range error calculated by the software receiver. It can be clearly seen from the figure that the accuracy of the pseudo-range measurement processed by the proposed method and the MVDR method is better than that of the PI algorithm. Under the condition of a weak signal, the PI algorithm is unable to obtain effective pseudo-range measurement due to the lock-out of the signal tracking loop, which means that the observation output of the receiver is no longer available. Although the accuracy of the pseudo-range measurement obtained by the proposed method has decreased, it can still ensure the availability of the pseudo-range measurement. The above results show that the proposed method can achieve the same performance as MVDR.
Figure 14 shows the DOA tracking results. The DOA tracking accuracy is closely related to the C/N0 value. When C/N0 is 45 dBHz, the MSE results are 0.03 degree2 under no interference and 0.08 degree2 under interference; when the C/N0 is 30 dBHz, the MSE value is reduced to 0.48 degree2. The simulation results show that even under the condition of weak signal, high-precision DOA tracking performance can still be achieved.

5. Conclusions

In order to achieve interference suppression and DOA estimation, a blind adaptive beamforming algorithm for a GNSS array receiver based on DiLL is proposed. Firstly, the subspace projection principle is used to achieve interference suppression, and the direction lock loop is used to achieve high-precision tracking of the DOA of navigation signals in the correlation domain. Then, the beamformer is constructed with the estimated DOA estimation results, which can provide gain by forming a beam in the satellite direction. After the array signals are weighted, interference suppression and navigation signal enhancement can be realized. The proposed method can realize DOA estimation in the presence of interference with low computational complexity, and thus achieve stable DOA tracking in dynamic scenes. It deserves great attention in the field of GNSS array antenna receivers.

Author Contributions

J.W. performed the theoretical study, conducted the experiments, processed the data, and wrote the manuscript. X.T. and L.H. provided research suggestions and revised the manuscript. S.N. helped in performing the experiments. F.W. provided the experimental equipment and suggestions for the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported in part by the Natural Science Foundation of China (NSFC) grants U20A0193.

Data Availability Statement

No date supporting.

Acknowledgments

The authors would like to thank the editors and reviewers for their efforts towards the publication of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, J.; Gao, K.; Guo, W.; Cui, J.; Guo, C. Role, path, and vision of “5G + BDS/GNSS”. Satell. Navig. 2020, 1, 23. [Google Scholar] [CrossRef]
  2. Elghamrawy, H.; Karaim, M.; Korenberg, M.; Noureldin, A. High-Resolution Spectral Estimation for Continuous Wave Jamming Mitigation of GNSS Signals in Autonomous Vehicles. IEEE Trans. Intell. Transp. Syst. 2021, 23, 7881–7895. [Google Scholar] [CrossRef]
  3. Lu, Z.; Song, J.; Huang, L.; Ren, C.; Xiao, Z.; Li, B. Distortionless 1/2 Overlap Windowing in Frequency Domain Anti-Jamming of Satellite Navigation Receivers. Remote Sens. 2022, 14, 1801. [Google Scholar] [CrossRef]
  4. Chaparro, D.; Feldman, A.F.; Chaubell, M.J.; Yueh, S.H.; Entekhabi, D. Robustness of Vegetation Optical Depth Retrievals Based on L-Band Global Radiometry. IEEE Trans. Geosci. Remote Sens. 2022, 60, 4413417. [Google Scholar] [CrossRef]
  5. Sun, Y.; Chen, F.; Lu, Z.; Wang, F. Anti-Jamming Method and Implementation for GNSS Receiver Based on Array Antenna Rotation. Remote Sens. 2022, 14, 4774. [Google Scholar] [CrossRef]
  6. Wolfley, O.; Schvartzman, D.; Mansur, B.; Fusselman, J.D.; Zhang, Y.R.; Snelling, B. Three-dimensional anti-jamming array processing for GNSS-based navigational aid inspection. In Proceedings of the Radar Sensor Technology XXVI 2022, Virtual, 6–12 June 2022; The Society of Photo-Optical Instrumentation Engineers (SPIE): Bellingham, DC, USA, 2022. [Google Scholar]
  7. Compton, R.T., Jr. The Power-Inversion Adaptive Array: Concept and Performance. IEEE Trans. Aerosp. Electron. Syst. 1979, AES-15, 803–814. [Google Scholar] [CrossRef]
  8. Yang, X.; Liu, W.; Chen, F.; Lu, Z.; Wang, F. Analysis of the Effects Power-Inversion (PI) Adaptive Algorithm Have on GNSS Received Pseudorange Measurement. IEEE Access 2019, 10, 70242–70251. [Google Scholar] [CrossRef]
  9. Lu, Z.; Nie, J.; Wan, Y.; Ou, G. Optimal reference element for interference suppression in GNSS antenna arrays under channel mismatch. IET Radar Sonar Navig. 2017, 11, 1161–1169. [Google Scholar] [CrossRef]
  10. Fernandez-Prades, C.; Arribas, J.; Closas, P. Robust GNSS Receivers by Array Signal Processing: Theory and Implementation. Proc. IEEE 2016, 104, 1207–1220. [Google Scholar] [CrossRef]
  11. Dai, X.; Nie, J.; Chen, F.; Ou, G. Distortionless space-time adaptive processor based on MVDR beamformer for GNSS receiver. IET Radar Sonar Navig. 2017, 11, 1488–1494. [Google Scholar] [CrossRef]
  12. Jia, Q.; Wu, R.; Wang, W.; Lu, D.; Wang, L. Adaptive blind anti-jamming algorithm using acquisition information to reduce the carrier phase bias. GPS Solut. 2018, 22, 99. [Google Scholar] [CrossRef]
  13. Wu, J.; Huang, L.; Tang, X.; Li, B.; Wang, F. Distortionless Blind Beamformer for Interference Suppression in GNSS Antenna Array Receiver. In Proceedings of the 11th China Satellite Navigation Conference, CSNC 2020, Chengdu, China, 22–25 November 2020; Springer Science and Business Media Deutschland GmbH: Berlin, Germany, 2020; Volume 652, pp. 733–741. [Google Scholar]
  14. Tan, J.; Nie, Z.; Peng, S. Quadratic forward/backward spatial smoothing polarization smoothing MUSIC algorithm for low angle estimation. In Proceedings of the 2019 IEEE Radar Conference, RadarConf 2019, Boston, MA, USA, 22–26 April 2019; Institute of Electrical and Electronics Engineers Inc.: Piscateville, NJ, USA, 2019. [Google Scholar]
  15. Zhang, W.; Han, Y.; Jin, M.; Li, X.-S. An Improved ESPRIT-Like Algorithm for Coherent Signals DOA Estimation. IEEE Commun. Lett. 2019, 24, 339–343. [Google Scholar] [CrossRef]
  16. Chen, X.; Morton, Y. Iterative subspace alternating projection method for GNSS multipath DOA estimation. IET Radar Sonar Navig. 2016, 10, 1260–1269. [Google Scholar] [CrossRef]
  17. Weiland, L.; Wiese, T.; Utschick, W. Multipath mitigation using OMP and Newton’s method for multi-antenna GNSS receivers. In Proceedings of the 7th IEEE International Workshop on Computational Advances in Multi-Sensor Adaptive Processing, CAMSAP 2017, Curacao, The Netherlands, 10–13 December 2017; Institute of Electrical and Electronics Engineers Inc.: Piscateville, NJ, USA, 2017; pp. 1–5. [Google Scholar]
  18. Rao, C.; Sastry, C.; Zhou, B. Tracking the direction of arrival of multiple moving targets. IEEE Trans. Signal Process. 1994, 42, 1133–1144. [Google Scholar] [CrossRef]
  19. Denic, N.; Dimitrijevic, B.; Miloevic, N.; Nikolic, Z. Three dimensional direction lock loop for DS/CDMA systems. In Proceedings of the TELSIKS 2007-8th International Conference on Telecommunications in Modern Satellite, Cable and Broadcasting Services, Nis, Serbia, 26–28 September 2007; IEEE Computer Society: Washington, DC, USA, 2007; pp. 11–14. [Google Scholar]
  20. Min, S.; Seo, D.; Lee, K.B.; Kwon, H.; Lee, Y.-H. Direction-of-Arrival Tracking Scheme for DS/CDMA Systems: Direction Lock Loop. IEEE Trans. Wirel. Commun. 2004, 3, 191–202. [Google Scholar] [CrossRef]
  21. Walid, G.D.; Hassan, M.E. DOA tracking in multipath environment based on the Direction Lock Loop. In Proceedings of the 2007 IEEE International Conference on Signal Processing and Communications, ICSPC 2007, Dubai, United Arab Emirates, 24–27 November 27; IEEE Computer Society: Washington, DC, USA, 2007; pp. 1099–1102. [Google Scholar]
  22. Wang, B.; Shivaramaiah, N.C.; Akos, D.M.; Wei, J. GNSS direction of arrival tracking using the rotate-to-zero direction lock loop. GPS Solut. 2020, 24, 39. [Google Scholar] [CrossRef]
  23. Cheng, X.-Z.; Zhu, D.-R.; Zhang, S.; He, P. Tracking Positioning Algorithm for Direction of Arrival Based on Direction Lock Loop. Futur. Internet 2015, 7, 214–224. [Google Scholar] [CrossRef] [Green Version]
  24. Guan, G.; Xian, D.; Zhang, K.; Gao, Y.; Zhu, X. Direction-of-arrival (DOA) tracking using improved direction lock loop with real-time bias correction for antenna array based GNSS receivers. Adv. Space Res. 2017, 60, 2733–2741. [Google Scholar] [CrossRef]
  25. Yang, B.; Cui, W.; Liu, Y.; Gao, B.; Mei, F.; Xu, H. DOA tracking using an improved direction lock loop based on a three-element L-shaped array. GPS Solut. 2022, 26, 108. [Google Scholar] [CrossRef]
  26. Sgammini, M.; Antreich, F.; Kurz, L.; Michael, M.; Tobias, G. Blind adaptive beamformer based on orthogonal projections for GNSS. In Proceedings of the ION ITM 2012, Nashville, TN, USA, 17–21 September 2012; Institute of Navigation: Manassas, VA, USA, 2012; pp. 926–935. [Google Scholar]
  27. Kai, B.; Akos, D.M.; Bertelsen, N.; Rinder, P.; Jensen, S.H. A Softwaredefined GPS and Galileo Receiver: A Single-Frequency Approach; Springer: New York, NY, USA, 2007. [Google Scholar]
Figure 1. Array configuration of the ULA.
Figure 1. Array configuration of the ULA.
Remotesensing 15 03387 g001
Figure 2. Normalized array pattern of the ULA.
Figure 2. Normalized array pattern of the ULA.
Remotesensing 15 03387 g002
Figure 3. The S-curves of the coherent and non-coherent DiLL.
Figure 3. The S-curves of the coherent and non-coherent DiLL.
Remotesensing 15 03387 g003
Figure 4. Block diagram of the proposed method, ( ) denotes complex conjugation.
Figure 4. Block diagram of the proposed method, ( ) denotes complex conjugation.
Remotesensing 15 03387 g004
Figure 5. MSE estimation results of DOA under different carrier-to-noise ratios: (a) DOA estimation results under 6 array elements; (b) DOA estimation results under 10 array elements.
Figure 5. MSE estimation results of DOA under different carrier-to-noise ratios: (a) DOA estimation results under 6 array elements; (b) DOA estimation results under 10 array elements.
Remotesensing 15 03387 g005
Figure 6. Tracking results under the scenario of fixed DOA parameters.
Figure 6. Tracking results under the scenario of fixed DOA parameters.
Remotesensing 15 03387 g006
Figure 7. Tracking results of DOA parameters in uniformly accelerated changing scenarios.
Figure 7. Tracking results of DOA parameters in uniformly accelerated changing scenarios.
Remotesensing 15 03387 g007
Figure 8. Tracking results in the scenario of sinusoidal variation in DOA parameters.
Figure 8. Tracking results in the scenario of sinusoidal variation in DOA parameters.
Remotesensing 15 03387 g008
Figure 9. Tracking results in the scenario of DOA parameter step change: (a) DOA tracking results; (b) DOA tracking error.
Figure 9. Tracking results in the scenario of DOA parameter step change: (a) DOA tracking results; (b) DOA tracking error.
Remotesensing 15 03387 g009
Figure 10. The flowchart of the simulation model.
Figure 10. The flowchart of the simulation model.
Remotesensing 15 03387 g010
Figure 11. Normalized array patterns of the proposed method, the PI method, and the MVDR method: (a) no-interference scenario; (b) interference scenario.
Figure 11. Normalized array patterns of the proposed method, the PI method, and the MVDR method: (a) no-interference scenario; (b) interference scenario.
Remotesensing 15 03387 g011
Figure 12. The punctual branch correlation cumulative value.
Figure 12. The punctual branch correlation cumulative value.
Remotesensing 15 03387 g012
Figure 13. Pseudo-distance measurement value.
Figure 13. Pseudo-distance measurement value.
Remotesensing 15 03387 g013
Figure 14. DOA tracking results.
Figure 14. DOA tracking results.
Remotesensing 15 03387 g014
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wu, J.; Tang, X.; Huang, L.; Ni, S.; Wang, F. Blind Adaptive Beamforming for a Global Navigation Satellite System Array Receiver Based on Direction Lock Loop. Remote Sens. 2023, 15, 3387. https://doi.org/10.3390/rs15133387

AMA Style

Wu J, Tang X, Huang L, Ni S, Wang F. Blind Adaptive Beamforming for a Global Navigation Satellite System Array Receiver Based on Direction Lock Loop. Remote Sensing. 2023; 15(13):3387. https://doi.org/10.3390/rs15133387

Chicago/Turabian Style

Wu, Jian, Xiaomei Tang, Long Huang, Shaojie Ni, and Feixue Wang. 2023. "Blind Adaptive Beamforming for a Global Navigation Satellite System Array Receiver Based on Direction Lock Loop" Remote Sensing 15, no. 13: 3387. https://doi.org/10.3390/rs15133387

APA Style

Wu, J., Tang, X., Huang, L., Ni, S., & Wang, F. (2023). Blind Adaptive Beamforming for a Global Navigation Satellite System Array Receiver Based on Direction Lock Loop. Remote Sensing, 15(13), 3387. https://doi.org/10.3390/rs15133387

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop