## 1 Introduction  
The Global Navigation Satellite System (GNSS) is now widely used to provide position, velocity, and timing information in a wide variety of applications [1]. Unfortunately, the GNSS signal is easily disturbed by both intentional and unintentional interference [2]. Therefore, it is of great importance to keep GNSS receivers operating properly in environments containing interference. The adaptive antenna array is a state-of-the-art anti-jamming approach that can maintain the functionality of GNSS receivers in a harsh interference environment [1].  
Blind anti-jamming algorithms, such as the power inversion (PI) algorithm [3], can suppress interference without any prior knowledge, making them popular in practical engineering [4]. However, apart from the problem of satellite signal attenuation [5], one of the main problems with these algorithms is that they cause phase distortions and introduce bias errors into the carrier phase [6]. The bias errors cannot be precalibrated since they vary in response to the operating environment [7]. Moreover, the biases vary from satellite to satellite [1], and can change significantly during interference suppression [8]. This measurement bias cannot be ignored, since it consequently results in a position deviation [9]. Therefore, maintaining high-precision measurement accuracy while suppressing interference is urgently needed.  
To solve this problem, some researchers have attempted to optimize the array configuration. A special centrosymmetric array configuration that does not cause any phase distortion has been proposed [10; 11]. However, this type of approach is not suitable when the nonideal factors of the array elements are taken into account. Moreover, sometimes, it is not easy to modify the array configuration in practice. Recently, a phase compensation technique was studied in the literature. The phase bias term is estimated and compensated for in this type of approach. The estimation of the phase bias can be divided into two approaches. The first approach is performed in the adaptive array processing. The phase-rotated steering vector of the satellite signal is estimated and used to obtain the phase bias estimate [12; 13]. However, these methods significantly increase the computational complexity. Furthermore, the steering vector is estimated from its orthogonal component in the interference subspace. The accuracy of the estimates cannot be guaranteed and will therefore introduce biases. The second approach is accomplished in the receiver signal processing. The main advantage of this method is that it involves only the addition logic in the receiver, which does not require significant modifications to existing systems. The phase discriminator is used to obtain the phase bias estimate [14; 15]. However, the estimation accuracy is easily affected by the phase discriminator noise. To address the above issues, the phase bias detection and estimation methods are investigated, and a robust phase compensation technique is proposed.  
The remainder of this paper is organized as follows. Section 2 gives the signal model and formulates the problem. Section 3 presents the proposed phase compensation technique. Section 4 performs simulation experiments and analyzes the results. Finally, the conclusions are drawn in Section 5.