Robust Wideband Interference Suppression Method for GNSS Array Antenna Receiver via Hybrid Beamforming Technique

: Global navigation satellite system (GNSS) array antenna receivers are widely used to suppress wideband interference in navigation countermeasures. However, existing array antenna receivers all adopt a digital array structure and digital beamforming technique, and they have limited analog-front-end (AFE) dynamic range. In strong interference scenarios, AFE saturation will occur, which limits the maximum interference suppression ability of the array receiver. Aiming at this issue, this paper proposes a robust wideband interference suppression method for GNSS array antenna receivers based on a hybrid beamforming technique. Firstly, a novel, fully connected hybrid array receiver structure is proposed. Secondly, the corresponding hybrid beamforming method is proposed at the same time, and it realizes the complete elimination of the strong wideband interference by joint suppression in the analog domain and digital domain. After mathematical simulations, it is veri ﬁ ed that, compared to the digital beamforming-based anti-jamming technique, the proposed method can e ﬀ ectively suppress strong wideband interference, and the maximum interference suppression ability is improved by 36 dB.


Introduction
The global navigation satellite system (GNSS) has the characteristics of all-weather and all-time, and it is widely used in various fields of society to obtain position, velocity, and time (PVT) information [1].Due to the low power of GNSS signals, they are highly susceptible to various intentional or unintentional interferences [2].In order to ensure the effective performance of GNSS receivers in harsh civilian and military environments, various anti-jamming techniques have been developed.They eliminate narrowband and wideband interference in the time domain and space domain [3,4].
However, various anti-jamming techniques can only be effectively implemented when the receiver analog-front-end (AFE, including the antenna, radio frequency (RF), and analog-to-digital converter (ADC) in a broad sense) responds linearly to the input combination of signal, noise, and interference [5].As a weak link, the AFE of GNSS receivers is susceptible to nonlinear effects such as saturation that distort the input waveforms due to the strong interference received [6,7].AFE saturation is an extreme case of the interference effect in strong interference scenarios.The linear dynamic range of the AFE sets the maximum interference suppression ability of the anti-jamming techniques [8].Therefore, how to relax the linear dynamic range of the receiver AFE to prevent saturation in strong interference scenarios is of great significance to improving the anti-jamming performance of the receiver.
In recent years, some techniques have been developed to deal with AFE saturation issues in GNSS receivers.The paper [9] proposes stopband filtering techniques to alleviate saturation caused by out-of-band and near-band interference.Paper [10] proposes two antenna modules composed of dynamically configured component topology networks.These two configurable antenna modules can detect, identify, and isolate potential interference.They switch modules with different gains based on the detected interference power to prevent AFE saturation.In the literature [11], an ideal and very slow automatic gain control (AGC) is placed at the AFE to prevent AFE saturation when continuous wave interference occurs.Papers [10,11] both ensure that the AFE does not experience saturation in strong interference scenarios by adaptively adjusting the AFE gain.However, all the methods mentioned above for dealing with AFE saturation focus on narrowband interference, and they are mainly applied to single antenna receivers.The methods for dealing with saturation caused by the strong wideband interference are rarely mentioned, and they are mainly applied to array antenna receivers.Firstly, wideband interference belongs to in-band interference, and the stopband filtering technology in the literature [9] cannot be used to prevent AFE saturation.Secondly, wideband interference is mostly intentional interference, and it has a greater interference power than narrowband interference, which mostly belongs to unintentional interference.Excessive AFE gain adjustment in the literature [10,11] may lead to insufficient amplification of weak GNSS signals, and the ADC sampled signal cannot obtain a sufficient signal-to-noise ratio (SNR).Therefore, the existing AFE saturation processing methods cannot be used in a strong wideband interference scenario.
There is some literature on the AFE saturation processing methods caused by the strong wideband interference in radar and communication fields.Paper [12] proposes a hybrid beamforming method to filter out the strong interference.A phase shifter array is added to the AFE of the antenna subarray, and the strong interference is jointly suppressed in the analog domain and the digital domain.Papers [13,14] propose a fully connected hybrid array chip design method for a massive multiple-input multiple-output (MIMO) system, and it prevents AFE saturation by suppressing interference in the analog domain.Papers [15,16] propose a hybrid array form with a fully connected structure and a hybrid beamforming method.The high-power wideband self-interference signal is suppressed in both the analog and digital domains.However, the aforementioned hybrid beamforming method still cannot be used in GNSS array antennas.Firstly, the GNSS array is a small aperture array, and it has a limited number of antennas [17].Thus, the method proposed for antenna subarray in the literature [12][13][14] is not applicable to GNSS array antennas.Secondly, the receiver only receives but does not transmit signals, and the prior information of the interference signal waveform cannot be obtained, so the method proposed in the literature [15,16] cannot be adopted either.
Aiming at the AFE saturation caused by the intentional strong wideband interference in the GNSS array antenna, a hybrid beamforming technique is proposed that simultaneously performs interference suppression in the analog domain and digital domain.Firstly, a novel, fully connected hybrid array anti-jamming receiver structure is proposed.Secondly, the corresponding hybrid beamforming technique is proposed at the same time.In the analog domain, multiple different 2D spatial analog notch filters are used to partially suppress interference and prevent AFE saturation.In the digital domain, the existing digital beamforming-based anti-jamming technique is calibrated and used to filter out residual interference.
The remainder of this paper is organized as follows: The second section briefly introduces the signal model and digital beamforming-based anti-jamming technique; the third section proposes the hybrid beamforming technique, which includes analog beamforming technique and digital beamforming technique; the fourth section gives the mathematical simulation results; and finally, the fifth section gives some conclusions about this paper.

Signal Model and Conventional Anti-Jamming Technique
In this section, the signal model of the anti-jamming receiver is first introduced, and then two conventional anti-jamming methods based on the digital beamforming technique are introduced.

Signal Model
The GNSS array antenna receiver generally adopts a small aperture circular array, as shown in Figure 1.The array is composed of M isotropic elements.Without loss of generality, assuming there is one GNSS signal and Q interference signals received by the array antenna, the array form of the received signal can be represented as follows: where t is the time, ( ) is the array snapshot with M-element array antennas.( ) s t is the GNSS signal, ( ) are the steering vectors of GNSS signal ( ) s t and interference signal ( ) q j t , respectively, and is the additive white noise vector.The steering vector of the array in Equation ( 1) can be represented as follows: ( ) ( ) where Here, 1 0 = ψ , λ represents the wavelength corresponding to GNSS signal center fre- quency f, i.e., λ= f c ; the ( ) denotes the pattern radiated by each element.Here, to simplify the derivation process, ( )

Digital Beamforming-Based Anti-Jamming Processing Technique
In order to suppress interference in Equation (1), a weight vector w is applied to the received signal vector, and the output signal of the array ( ) y n is shown below.
( ) ( ) where ( ) is the conjugate transpose operation, and ( ) n x denotes the signal ( ) In the field of conventional digital beamforming-based anti-jamming processing techniques, there are two techniques for obtaining the weight w: the minimum variance distortionless response (MVDR) technique and the power inversion (PI) technique [18,19].
The MVDR technique is widely used in various fields, and its objective function is as follows: where represents the covariance matrix of the received signal, and N denotes the number of snapshots.The solution of Equation ( 5) is given by The objective function of the PI technique is as follows: Different anti-jamming techniques are adopted according to different application scenarios.However, both of the above techniques belong to the digital beamforming technique.If there is signal distortion caused by AFE saturation in strong interference scenarios, both techniques cannot work.The interference suppression ultimately cannot be completed, and then the GNSS PVT solution cannot be performed.

Proposed Technique
The proposed hybrid beamforming technique follows in this section.The hybrid array antenna receiver structure is shown in Figure 2.The receiver adopts a fully connected hybrid array structure containing multiple R-modules at the RFE.Each R-module contains an attenuator and a phase shifter, as shown in Figure 3.
Hybrid beamforming techniques include analog beamforming techniques and digital beamforming techniques, as shown in Figure 2.
(1) Analog beamforming technique.Taking the direction of arrival (DOA) of interference as a priori information, the multiple-notch filter (MNF) composed of M different 2D spatial-notch filters (SNF) is generated simultaneously in the baseband, as shown in Figure 2. The signal is filtered in the analog domain to partially suppress the strong interference and retain the original DOA distribution of the signal (including the GNSS signal and interference signal).
(2) Digital beamforming technique.After the data are filtered by the analog beamforming technique, they are sampled into the baseband through the RF channel and ADC.
The original digital beamforming-based anti-jamming processing technique, that is, Equations ( 6) and ( 8), needs to be calibrated and then adopted to suppress interference in the digital domain.In the end, the final suppression of interference is completed.
In this paper, the analog beamforming weights are represented by h, and the digital beamforming weights are represented by w.

Analog Beamforming Technique
Firstly, the analog beamforming technique is introduced.As shown in Figure 2, the GNSS hybrid array receiver includes a total of M SNFs, represented by Each SNF cannot be generated independently.The MNF not only needs to achieve partial suppression of interference in the analog domain but also strives to retain the original signal's DOA information as much as possible.After GNSS and interference signals are filtered by the MNF, the signal can continue to retain the original DOA distribution in the baseband.Therefore, effective implementation of subsequent baseband anti-jamming algorithms can be ensured.The bampattern synthesis method is used to generate the MNF [20].
In the synthesis process of the MNF, different steering vector forms are used as follows: Here,

( )
m Ω a represents the steering vector corresponding to the m-th SNF; that is, the corresponding antenna is used as the reference antenna for each The SNF and MNF constraints are simultaneously applied to the synthesized pattern.The first SNF 1 h is generated by the SNF constraint.Then, using 1 h as a reference mask, the other notch filters m m M , 2,3, , = h  are generated by the MNF constraints.
(1) SNF Constraints The SNF constraints include constraints on both non-interference and interference areas.The non-interference area and interference area are constrained separately.The 3 dB main lobe width of the interference spatial spectrum is used as a basis to distinguish between the non-interference area and the interference area.The area outside the 3 dB main lobe width of the interference spatial spectrum is the non-interference area, and the other area is the interference area.
The constraint for the non-interference area is as follows: ( ) where the lower and upper magnitude bounds are applied [20].In Equation ( 10), the noninterference area is divided into P uniformly spaced angles, Θ represents the set of P angles in the non-interference area, i.e.,

{ }
, and δ is the ripple term.
The constraint for interference area is as In Equation ( 11), Q represents the number of interferences, Ψ represents the set of Q , and ς is the nulling term.
As shown in Equations ( 10) and ( 11), P angles are used to constrain the non-interference area, but only one single angle is used to constrain each interference.This situation will lead to the inability of the objective function to converge and generate a deep enough null that meets the requirements.To solve this problem, the method of virtual interference is adopted, and the K interferences at the same interference direction are assumed.In practice, the value of K should be greater than or equal to the number of non-interference area constraint angles P, that is, K ≥ P.And then, Equation ( 11) is written as follows: ( ) where ( ( ) ( ) ( ) Here, ( ) x t represents the GNSS signal and interference signal.Substitute Equation (9) into Equation ( 14), the array input signal for each SNF in MNF can be represented as follows: The output of the signal ( ) G represents the beampattern of the m-th SNF.
Therefore, the vector form of the output signal, after being filtered by the MNF, is In the MNF beampattern synthesis, the below constraints are applied.
( ) ( ) Then, Equation ( 17) can be written as follows: Through Equation ( 19), it can be seen that the output data of the MNF can retain the original DOA distribution.The Equation ( 18) for the MNF constraints can be written as follows: Here, Equations ( 20) and ( 21) represent the consistency constraints for non-interference and interference areas between the m-th SNF and the first SNF, respectively.η and ξ represent the response error terms.The 2 l norm of 1 h is adopted as the objective function, which is used to receive useful signals with high gain.The objective function is shown below.
( ) ( ) ( ) ( ) Due to Equation (10), the objective function ( 22) is a nonconvex optimization problem, and the alternating directions method of multipliers (ADMM) optimization method is used to solve the objective function.It decomposes the problem into multiple sub-problems, iteratively solves the sub-problems in sequence, and obtains the final optimization result.The specific derivation process can be referred to in the literature [20,21]  is obtained through Equation ( 22), and it requires less than 100 iterations.After quantization, the analog beamforming coefficient

Digital Beamforming Technique
Secondly, a digital beamforming technique is introduced.The MVDR technique and PI technique in Equations ( 6) and ( 8) are still used in the digital beamforming part.However, through Equations ( 20) and ( 21), it can be seen that although we have retained the DOA distribution of the original GNSS signal as much as possible, a small amount of amplitude/phase error is inevitably introduced.Therefore, in order to maximize the interference suppression ability of the digital domain anti-jamming technique, two processes should be added: amplitude/phase error calculation and calibration.
Then, the calculation process of amplitude/phase error is as follows: Finally, the updated MVDR algorithm is obtained as follows: Where ( ) ( ) ( ) , and represents the operation of transforming a vector into a diagonal matrix.
Due to the fact that a small amount of amplitude/phase error has no impact on the PI algorithm, Equation ( 8) is not updated in this section.

Implementation of Proposed Technique
Figure 2 illustrates the data processing flowchart of the robust wideband interference suppression method for the GNSS array antenna receiver based on a hybrid beamforming technique.The steps of the method are as follows: Step 1: Confirm whether the interference's DOA information has been obtained.If the interference's DOA information is obtained, move to Step 5.If not, move to Step 2.
Step 2: Confirm whether ADC saturation occurs.If ADC saturation occurs, move to Step 3. If not, move to Step 4.
Step 3: The RF gain is adjusted from 30 dB to −10 dB to ensure that ADC does not experience saturation.The interference DOA information of the first epoch can be obtained.
Step Step 5: Confirm whether the RF gain is set to 30 dB and set the RF gain to 30 dB.The prior information on interference direction is first obtained from DOA estimation, and then the interference angle set q  is written down to the corresponding phase shifters and attenuators in the hybrid array anti-jamming receiver.
Step 10: After filtering through the analog MNF, the signal is amplified, filtered, and sampled to the digital baseband signal Step 11: DOA estimation results can be obtained using a multiple signal classification (MUSIC) algorithm [22] on baseband data, and then DOA estimation results are sent to Step 5 in real-time.
Step 12: Using Equation (25), a digital beamforming technique is performed on the baseband data The robust wideband interference suppression process based on the hybrid beamforming technique is completed.
This section presents a hybrid beamforming-based GNSS anti-jamming method.Compared with the conventional digital beamforming-based anti-jamming method, there will be an obvious increase in hardware structure and implementation complexity.In terms of hardware structure, M 2 phase shifters and M 2 attenuations need to be added to the AFE, as shown in Figure 2. In terms of implementation, DOA estimation, analog MNF calculation, amplitude/phase error calculation, and calibration need to be added to the implementation, as shown in Figure 4.

Simulations and Results
In this section, the performance of the proposed method in this paper is evaluated through simulations, which are performed using the MATLAB R2018b software platform.In Section 4.1, the simulation settings are presented.In Section 4.2, Simulation 1 verifies the analog beamforming performance of the proposed method in two types of interference scenarios.In Section 4.3, Simulation 2 verifies the strong interference suppression performance.In Section 4.4, Simulation 3 verifies the maximum interference suppression ability of the proposed method.In Simulations 2 and 3, the performance of the proposed method is compared with the digital beamforming-based anti-jamming technique.In order to fully verify the interference suppression ability of the proposed technique, the MVDR technique is used in hybrid beamforming-based and digital beamforming-based anti-jamming processing, respectively.

Simulation Settings
In the mathematical simulation verification, we use a 7-element array antenna, as shown in Figure 1.The GNSS receiver simulation parameters are shown in Table 1.The  2. The simulation includes six interference scenarios.Scenarios 1-3 are the single-interference scenarios, and Scenarios 4-6 are the twointerference scenarios.The GNSS constellation adopts the real constellation, as shown in Table 3.The parameter settings of the analog beamforming technique in Section 3.1 are shown in Table 4.

Simulation 1: Analog Beamforming Performance Verification
The analog beamforming technique is the key part of the proposed method.In this section, the performance of the analog MNF is verified, including the performance verification of each SNF and the consistency analysis between each SNF in the MNF.The performance of each SNF is calculated by Equation ( 23).The gain G in the non-interference area and the nulling depth D in the interference area are calculated.The calculation formulas of the amplitude error m A and phase error m P between each SNF and the first SNF in the non-interference area are as follows: where ⋅ represents the absolute value operation and ∠ represents the operation of taking complex phase.Figures 5 and 6 show the amplitude/phase response of the analog MNF.The performance of analog MNF is calculated, where max G , min G , mean A , max A , mean P , and max P represent the maximum gain, minimum gain, average amplitude error, maximum amplitude error, average phase error, and maximum phase error in the non-interference area, respectively, and D represents the nulling depth in the interference direction.The detailed performance results are shown in Table 5.
From Table 5, it can be seen that in the analog MNF, the maximum gain of all SNFs max G is 2.687 dB, the minimum gain min G is −4.683 dB, and the minimum nulling depth of each SNF in the analog MNF D is −35.07 dB.The amplitude/phase error between each SNF and reference SNF-1 in the analog MNF is that the maximum average amplitude error mean A is 1.812 dB, the maximum amplitude error max A is 4.508 dB, the maximum average phase error mean P is 3.96°, and the maximum phase error max P is 17.47°.This amplitude/phase error could meet the requirements of the array antenna.

Simulation 2: Strong Interference Suppression Performance Verification
In this section, the strong interference suppression performance of the proposed method is verified and compared with the digital beamforming-based technique.The interference power is set to 0 dBm; that is, the input interference-to-noise ratio (INR) is 100 dB.
Firstly, the interference suppression ability of the analog MNF in the analog domain is verified.The nulling depth values for two scenarios (Scenarios 1 and 4) can be seen in Tables 5 and 6  it can be seen that, due to the partial suppression of interference in the analog domain, the technique proposed in this paper can effectively relax the linear dynamic range of AFE to prevent saturation in strong interference scenarios.However, due to the presence of strong wideband interference, digital beamforming-based antijamming array antenna receivers will experience AFE saturation, which results in signal distortion and the inability to operate baseband processing.After the interference suppression ability of the analog MNF is verified, the interference suppression ability of a hybrid beamforming-based array antenna receiver in strong interference scenarios will be verified.The interference suppression ability of a hybrid beamforming-based array antenna receiver is verified by comparing the output signal-to-interference-noise ratio (SINR) of the two array forms under strong interference conditions.Due to AFE saturation, the sampled signal is distorted, and the digital beamforming in two array forms is ineffective.Therefore, the statistics of SINR after AFE saturation are no longer meaningful.The SINR is uniformly set to −50 dB when AFE experiences saturation.

Simulation 3: Maximum Interference Suppression Ability Performance Verification
In this section, the maximum interference suppression ability of the proposed method is verified and compared with the beamforming-based technique.
To verify the improvement of the proposed technique in its maximum interference suppression ability, the output SINR of all GNSS signals in Table 2 is calculated under different input INR conditions in six scenarios.The interference power is from −40 dBm to 30 dBm, and the step length of the interference power is 2 dBm; that is, the INR is from 60 dB to 130 dB.
Figures 13 and 14, respectively, show the output SINR curves of all GNSS signals under different INR conditions for the two beamforming techniques in Scenarios 1-3.From Figure 13, it can be seen that when the input INR is ≤ 80 dB, the digital beamformingbased array antenna receiver can complete anti-jamming processing normally.The number of satellites that can be captured in three scenarios is 10, 9, and 10, respectively.And the output SINR of the captured satellites is greater than −25 dB.Therefore, the GNSS PVT solution can be completed.When the input INR is greater than 80 dB, the digital beamforming-based array antenna receiver experiences AFE saturation.All GNSS signals' output SINRs are −50 dB, and the number of available satellites is 0. The receiver is unable to complete the GNSS PVT solution.Therefore, the maximum interference suppression ability of the digital beamforming-based array antenna receiver is INR = 80 dB.From Figure 14, it can be seen that the AFE saturation of the hybrid beamforming-based array antenna receiver occurs when the input INR > 116 dB.When the input INR is ≤ 116 dB, the number of satellites that can be captured in three scenarios is 10, 9, and 10, respectively.Moreover, the output SINR of the captured satellites is greater than −25 dB.The receiver can complete the GNSS PVT solution normally.Therefore, the maximum interference suppression ability of a hybrid beamforming-based array antenna receiver is INR = 116 dB.Therefore, in Scenarios 1-3, compared with the digital beamforming-based array antenna receiver, the maximum interference suppression ability of the hybrid beamforming-based anti-jamming receiver is improved by 36 dB.

Discussion
In the present study, the method proposed in this paper takes the interference's DOA as the prior information, and it requires continuous tracking of the interference direction from the weak interference scenario to the strong interference scenario.Figure 17

Conclusions
Aiming at the AFE saturation problem caused by the strong wideband interference in GNSS array antenna receivers.The hybrid beamforming-based anti-jamming technique is proposed.The main contributions are as follows: (1) A novel, fully connected hybrid array anti-jamming receiver structure is proposed.It is used for hybrid beamforming by adding the attenuators and phase shifters to the AFE.
(2) A hybrid beamforming method for the GNSS array antenna receiver is proposed.In the analog domain, the interference is partially suppressed (better than 30 dB) by analog MNF generated in the baseband, and it relaxes the linear dynamic range of AFE to prevent saturation under the strong interference scenario.In the digital domain, the calibrated MVDR algorithm is used for the final suppression of residual interference.
(3) Through six interference scenarios, compared with the digital beamforming-based array antenna receiver, it is verified that the maximum interference suppression ability of the hybrid beamforming-based array antenna receiver is improved by 36 dB under the conditions of using 7-bit attenuators and 8-bit phase shifters.
There are several issues that should be considered in future work: (1) Space-time hybrid array structures should be considered.Firstly, it can provide higher interference suppression ability in the analog domain.Secondly, it can increase the degree of freedom to ensure that the signal has better amplitude consistency, especially phase consistency, after it is filtered by the analog MNF.(2) The algorithm in this paper assumes that the array antenna is ideal and ignores the influence of amplitude/phase error in practical applications.In analog beamforming, the influence of non-ideal factors will be further considered, and the analog beamforming algorithm will be modified accordingly.(3) This article only verifies the improvement of anti-jamming ability under static conditions.The antijamming performance under dynamic conditions needs to be verified through field experiments, and further consideration needs to be given to the width parameter setting of nulling in analog beamforming.

Figure 2 .
Figure 2. The hybrid array receiver structure.

Figure 3 .
Figure 3.The components meaning of the R-module.

4 :
Figure2illustrates the data processing flowchart of the robust wideband interference suppression method for the GNSS array antenna receiver based on a hybrid beamforming technique.The steps of the method are as follows:Step 1: Confirm whether the interference's DOA information has been obtained.If the interference's DOA information is obtained, move to Step 5.If not, move to Step 2.Step 2: Confirm whether ADC saturation occurs.If ADC saturation occurs, move to Step 3. If not, move to Step 4.Step 3: The RF gain is adjusted from 30 dB to −10 dB to ensure that ADC does not experience saturation.The interference DOA information of the first epoch can be obtained.Step 4: Set the weight of the 2D spatial notch filter to m m M T [1, 0,..., 0] , 1, 2, , = = h  , and then move to Step 7.Step 5: Confirm whether the RF gain is set to 30 dB and set the RF gain to 30 dB.The prior information on interference direction is first obtained from DOA estimation, and then the interference angle set q Ω q Q , 1, 2, , ∈ = Ψ  and the non-interference angle set

Figure 4 .
Figure 4.The flowchart of the proposed method.

N
Obtain interference direction Ω j , interference area and non interference angle set Θ and ΨStartAccording to Equation (22), calculate M SNFs hm, m = 1, 2, ..., M Whether the Equation (22) meets the convergence requirements Quantify the weight hm, m = 1, 2, ..., M，the analog beamforming weight hQuan-m, m = 1, ..., M is obtained Write quantified weights into the corresponding phase shifters and attenuators for analog interference suppression The signal that completes analog interference suppression is amplified, filtered, and sampled by ADC DOA estimation of interference Using Equation (25) to perform digital beamforming anti-jamming processing on the baseband signal End N Y Calculate the amplitude/phase error according to Equation (24) and obtain the matrix Γ Whether ADC saturation occurs under 30dB RF gain?The RF gain is adjusted from 30dB to -10dBhm=[1,0,…,0] T , m = 1, 2, ..., MWhether the interference's DOA information is obtained?
, but they do not accurately reflect the interference suppression ability of the MNF.The interference suppression ability of the MNF is accurately demonstrated by comparing the input INR and INR before ADC sampling in the receiver.The INR of the two scenarios after their respective MNFs is shown in Figures 9 and 10.In the figures, red represents the INR after the analog MNF, the blue color indicates the input INR, and both are also the INR before ADC sampling in two array forms.In Figure 9, it can be seen that after the interference is filtered by the analog MNF, the input INR of the hybrid array decays from 100 dB to 61.6 dB, 59.95 dB, 58.06 dB, 61.77 dB, 59.98 dB, 62.82 dB, and 64.8 dB, respectively.The wideband interference suppression ability of each SNF is 38.4 dB, 40.05 dB, 41.94 dB, 38.23 dB, 40.02 dB, 37.18 dB, and 35.2 dB, respectively.In Scenario 4, after two interferences enter the channel, the interference's power should be the sum of the two interferences, so the input INR increases to 106 dB.In Figure 10, the input INR of the hybrid array decays from 106 dB to 64.91 dB, 60.5 dB, 60.49 dB, 67.1 dB, 63.31 dB, 64.61 dB, and 67.03 dB, respectively.The wideband interference suppression ability of each SNF is 41.09 dB, 45.95 dB, 45.51 dB, 38.9 dB, 42.69 dB, 41.39 dB, and 38.97 dB, respectively.The interference suppression ability of Scenario 4 is achieved through the combination of two nulls.From Figures 9 and 10 ,

Figure 9 .
Figure 9.The INR of each RF channel before ADC sampling (Scenario 1).

Figure 10 .
Figure 10.The INR of each RF channel before ADC sampling (Scenario 4).

Figures 11 and 12
, respectively, show the output SINR of the two scenarios (Scenarios 1 and 4) in Table 2 under input INR = 100 dB.Due to AFE saturation in the digital array at this INR, the output SINR of all GNSS signals is −50 dB.However, in the hybrid array, due to the partial suppression of interference by analog MNF in the analog domain, the AFE does not experience saturation, and the digital beamforming of the hybrid array can work normally.The output SINR of each GNSS signal is calculated after digital beamforming in the hybrid array.As shown in Figure 11, the output SINR of GNSS signal in Scenario 1 are −12.33 dB, −17.49dB, −14.69 dB, −10.43 dB, −9.49dB, −11.46 dB, −9.56 dB, −9.01 dB, −9.78 dB, and −9.33 dB, respectively; As shown in Figure 12, the output SINR of GNSS signal in Scenario 4 are −22.39dB, −11.96 dB, −15.41 dB, −10.32 dB, −15.87 dB, −22.92 dB, −10.02 dB, −9.21 dB, −9.85 dB, and −38.22 dB, respectively.In the two scenarios, after the signal is processed by the hybrid beamforming technique, the number of GNSS signals with an output SINR greater than −25 dB is 10 and 9, respectively, which can be captured.The hybrid beamforming-based array antenna receiver can effectively complete the GNSS PVT solution under input INR = 100 dB.However, due to AFE saturation, the digital beamforming-based array antenna receiver cannot work normally.

Figures 15
Figures 15 and 16, respectively, show the output SINR curves of all GNSS signals under different INR conditions of the two beamforming techniques in Scenarios 4-6.The AFE saturation position of the digital beamforming-based array antenna receiver is input INR = 74 dB.However, the AFE saturation position of the hybrid beamforming-based array antenna receiver is input INR = 110 dB.Therefore, in Scenarios 4-6, compared with the digital beamforming-based array antenna receiver, the maximum interference

Table 1 .
Parameters of the GNSS receiver.

Table 2 .
The interference distribution.

Table 3 .
The GNSS constellation distribution.

Table 4 .
Parameters of the analog beamforming.

Table 5 .
Performance statistics of analog MNF * (Scenario 1).Figures7 and 8show the amplitude/phase response of the analog MNF in Scenario 4, and Table6presents the detailed performance of the analog MNF.In this scenario, the minimum nulling depth of each SNF in the analog MNF D is −34.22 dB, the maximum gain
min G max G 1 D 2 D mean A max A mean P max P are degree.