Abstract

In this study, a robust adaptive beamformer based on constant modulus (CM) criteria is developed to improve the robustness of the array beamforming, which is a reconstructing minimal optimization for solving the mismatch problem of weight vector caused by steering vector mismatch. In the global positioning system (GPS) L1 band, firstly, a GPS array signal is modelled by designing a dual-polarized antenna array. Secondly, the distortion problem of beamforming is formulated in the traditional minimum variance distortionless response (MVDR) beamformer. For repairing the weight vector mismatch problem, a novel beamformer based on the CM envelope response is proposed to reconstruct MVDR beamforming in the array processing. Besides, min-max penalty criteria are used to enable the beamformer to allocate more degrees of freedom (DOFs) when penalizing the MVDR beamformer responses. Finally, an auxiliary two-element real variable is designed to plan the proposed beamformer. But it is still a nonconvex quadratic programming problem, so an alternating direction method of multipliers (ADMM) is employed to transform the objective function into several subproblems. Illustrative numerical simulation results are provided for validating the effectiveness of the proposed beamformer by comparing it with other existing approaches.

1. Introduction

With the development of a robust adaptive beamformer (RAB), the concerns over its reliability also heighten in the global positioning system (GPS). But, because of multiple interferences (i.e., jamming signals, calibration errors, and distorted antenna shape), the antijamming robustness of RAB severely degrades in real applications, wherein the mismatching error between the presumed and actual response is the most commonly affected by the interferences [1]. To solve this problem, a special case of null-steering technique-based dual-polarized array response is proposed to provide double degrees of freedom (DOFs) [2, 3]. Here, it cannot only help the beamformer obtain the desired signal from the mismatch polarization [2] but also prove the more 2 N-1 DOFs in comparison with the N-element single-polarized array.

Beyond that, some diagonal loading (DL) techniques with various uncertainty sets of the steering vector (SV) are studied to obtain a good antijamming RAB. Such as in [4], the sample covariance matrix of RAB is regularized by using a symmetric tridiagonal loading factor. Undeniably, a suitable factor of the projection subspace is difficult to achieve, affecting the robustness of antijamming beamforming. Then, a time-variant factor via the uncertain set constraints is proposed in the paper [5]. However, it is still faced with a crucial problem with the selection of the loading level. Then, a fixed DL factor-based eigenvalue is selected to enhance the antijamming robustness of array processing in [6]; however, this approach must be built on a fixing transformation of the covariance matrix. In [7], a DL beamformer based on the sample matrix inversion is presented to adaptively steer the beamforming, but it is still sensitive to the accuracy of the sample covariance matrix of the array response. Additionally, in [8], adaptive beamforming is obtained by using the DL technique, which is able to steer the beamforming at a small number of samples. Unfortunately, this DL beamforming is subject to the limitations of the small sampling snapshots so that exists a mismatch error between the likelihood estimation and the reality in the signal processing response. In particular, if the desired signal is contained in the small sampling snapshots, the imprecise SV estimation of these DL beamforming severely degrades the antijamming robustness.

In another class of RAB techniques, eigenspace-based approaches are exploited to construct the subspace of the input signal. For example, in [9], the antijamming weight vector of the subaperture is obtained to generate some nulls from the jamming direction in the eigenspace-based beamformer; however, some prior conditions must be known when the array beamforming is achieved, otherwise not. In [10], the sampling covariance matrix is reconstructed to obtain the antijamming RAB by using the orthogonality between the SV subspace and the interference subspace, but there still exists a calibration error in steering beamforming. Then, an eigenspace-based beamformer is reconstructed by solving the linear problem of SV subspace projection in the array response, i.e., the Capon spatial spectrum in [11]. But this eigenspace is still sensitive to the sample snapshots and the input signal-to-noise ratio (SNR). Additionally, there are some antijamming beamformers that combat the small snapshots, enhancing the robustness of beamforming. For example, in [12], a Lagrange multiplier based on loading sample matrix inversion (LSMI) is used to obtain the RAB in each iteration. But due to the array calibration error, the no-vanishing SV error still severely degrades the robustness of beamforming. In another example from [13, 14], there are two different beamformers being designed by using the modified NCCB (norm-constrained Capon beamformer) approach and the WC (worst case) method, but the computational time can cost higher than that of LSMI for the same number of antennas. In [15], an antijamming problem of RAB is solved by using the spatial filtering techniques in the radar beamforming, but the precision of beamforming also suffers from the SV error in the iteration; meanwhile, the computational cost of that RAB is higher than the smaller interval needed. Then, a RAB based on minimum variance distortionless response (MVDR) is considered by solving the output power maximization problem and SV mismatch error in [16]. Although the MVDR beamformer can obtain the signal-of-interest steering vector through a double-sided norm perturbation constraint and a similarity quadratic constraint, it is still a quadratically constrained quadratic programming (QCQP) problem with inhomogeneous constraints. When facing the QCQP problem of RAB, the mismatch error of the array response can seriously reduce the robustness of the beamformer. For solving the mentioned QCQP problem, this paper employs a variable splitting technique such as the alternating direction method of multipliers (ADMM), where the ADMM technique is an efficient dual-decomposed optimisation approach that combines the dual-ascent method and the convergence bounded multiplier method. In particular, the dual-ascent performance by adding the convex penalty terms can bring numerical robustness [17].

In this paper, a novel RAB based on the constant modulus (CM) penalty criterion is developed to reconstruct the traditional MVDR beamformer [16]. By contrast, due to CM equalization with constant envelope features, this paper presents the advantage of higher accuracy than the MVDR beamformer [16] under the optimizing weight of RAB per iteration. And the proposed method is computationally superior to the state-of-the-art QCQP problem by using the ADMM variable splitting technique; this is because ADMM has the characteristic of simple update rules and stable-fast convergence to solve the corresponding control variables at each iteration. In addition, our method highlights, compared with the beamformer [16], the following points: (a)The sample covariance subspace is reconstructed by using the envelope response of the CM array processing in the collected snapshots. This allows the beamformer to enhance the robust adaptation of weight from the spatial blind beamforming(b)A min-max penalized optimization criterion can manage the CM spatial responses from the sampling interference subspace, which can minimize the maximal penalty term to allocate more DOFs to suppress the interference as well as improve the robust adaptation of the potential SV mismatch constraints in blind beamforming(c)By designing two auxiliary real-valued conditions, the objective function is converted to solve the rank uncertainty problem, followed by the CVX tool. Then, the beamformer can obtain better robustness than other approaches at a low cost

This paper is organized as follows. In Section 2, a GPS signal is modelled by using a dual-polarized GPS antenna array. In Section 3, the antijamming distortion problem is formulated from a traditional MVDR beamforming. Next, a novel RAB based on the constant modulus (CM) penalty criterion is proposed in Section 4. Finally, the numerical simulations and conclusions are drawn in Section 5 and Section 6, respectively.

Notation 1. Table 1 lists this article’s operators.

2. GPS Signal Model

This paper considers a uniform square antenna array (USAA), as shown in Figure 1, where dual-polarized GPS antenna elements are designed on a dual-layer substrate of “Rogers RO3010 (tm)” substrate with a relative permittivity of 10.2 and a dielectric loss tangent of 0.0035, and the distance between antenna elements is half λ-wavelength (i.e., about  mm in the GPS L1 band). In USAA, the 8-ports are matched by a 50 Ohm impedance, and then the radiation gains are motivated to obtain the desired GPS signal by orthogonal linear polarization, where is an column vector at the time instant as follows: where is the amplitude of the minimum guaranteed signal power from the GPS receiver, the subscripts and stand for the horizontal and vertical polarization (HP and VP), respectively, and is the maximum gain of HP or VP of the -th antenna element. The jamming signal has a similar form to the desired signal , which is only considered as deception jamming. Then, in this paper, the GPS transmission signal is designed in the GPS L1 band, which is a composite model including a desired GPS signal , jamming signal , and noise , and all given signals are uncorrelated, as explained in the following equation: where and denote the SV at the azimuth-pitch angle (, ) and (, ) direction, obeying to as a Gaussian noise is used in , the symbol Γ is defined as an manifold error matrix to repair the desired signal “” and jamming signal “”, as follows: where represents the -th gain of USAA, is the -th phase shift error corresponding to the -th independent random variables, i.e., the -th steering element as the desired -signal.

Then, the corresponding expressions are shown as follows: where is the relative ratio between the element distance and the wavelength and is a pair of the azimuth angle and pitch angle. Let represent the total SV with array error. Note that RAB is implemented in the expected environment, removing pseudorandom code and carrier error. Namely, the uncorrelated characteristics of a GPS signal are assumed in the stage after despreading and before carrier synchronization.

3. Problem Formulation

In traditional MVDR beamformer [16], the directions of arrival (DOAs) of transmission signal and SV are usually known in advance. And let as the complex weight variable of beamforming, then the output signal from snapshots is as follows:

Thus, the output energy is

When , is rewritten as where is the covariance vector of the array response. This beamformer based on the MVDR criterion has the minimum output energy, while the target directional gain keep unchanged, i.e., , then the model of the MVDR beamformer is obtained as follows: where each variable of is explained as weight vector, covariance matrix, steering vector, and azimuth-pitch angle about the beamforming. By using the Lagrange multiplier method, the cost function of the above problem is produced as follows:

When the partial derivative of about is zero, the final solution to (10) is obtained as follows: where . In MVDR beamformer, due to the SV mismatch error including of random error and calibration error [18], the array response can cause weight mismatch problems as follows:

Then, in the following section, to reduce the weight mismatch error of beamforming, a novel RAB on based the CM penalty criterion is developed to improve the robustness of beamforming.

4. Proposed Method

4.1. Formulation

Firstly, for limiting the target SV distortion, a nonlinear inequality relax constraint of is used to repair the constraint condition of the MVDR beamformer. Secondly, a CM envelope response is employed to reconstruct the target output energy of the MVDR beamformer in the collected snapshots. This CM blind channel equalization can obviously improve the robust adaptation of the potential constraints from weight mismatch processing. Finally, a novel optimization problem is formulated as follows: where and are defined as the discrete angular set of azimuth and pitch direction, the parameter is assumed to be the MVDR constraint supremum, is a given constant parameter, and is the received signal amplitude of the first snapshot, and the superscript represents the number of constant modulus. In (13), the first constraint states the estimated weight modulus of the allocation variables , and the second constraint represents a nonlinear inequality relaxation to restrict the sidelobe level of beamforming for repairing the SV, i.e., (-30 dB). Next, to consolidate the beamformer, this work is mainly explored from two perspectives:

(1) A min-max recursive criterion-based inequality relaxed constraints is proposed to penalize the spatial responses of interferences as follows where the penalty term as possesses the minimum output energy as from the MVDR beamformer. Namely, the maximum penalty parameter is regarded as the suppression preference; interferences with the larger have a higher priority to be suppressed. This allows to minimize the energy function to obtain the robust response in beamforming. Note that the inequality constraints in (14) are always feasible, and the penalty term can automatically allocate DOFs given by the dual-polarized USAA.

(2) The antijamming weight vector is obtained by reconstructing the sample covariance subspace of the MVDR beamformer while solving the minimum CM problem of the output beamforming , as follows where .

By combining (14) and (15), the estimated problem in (13) can be recast to a new object function as follows: where is used for the trade-off between antijamming and reducing noise. In practice, the optimization problem in (16) is still nonconvex programming. In general, semidefinite relaxation (SDR) could be used to solve it. But if the rank of the solution is not 1 when SDR is used to eliminate the equation, , a nonoptimal solution is possibly obtained [19]. Then, a variable splitting technique based on two-element real-valued conditions is employed to avoid this situation in the following subsection, such as ADMM.

4.2. Solution Procedure via ADMM

In this subsection, the ADMM method is used to transform the cost function (12) into several subproblems by variable splitting technique, and for bypassing the difficulty of this work, a penalty parameter as well as setting , are used to rewrite the object functions (16) as follows:

Motivated by [20], there are two auxiliary real-valued variables being used as follows:

where .

Then

Based on the above equations with , the multiple constraints with only the real-valued variables can be tackled separately as

After that, a novel equation based on the augmented Lagrangian method is defined as: where are the Lagrangian multipliers vectors and is the step size. Each computational process is shown as follows: (1)Solving

By omitting the irrelevant terms, the solution for in equation (21) is where and are explained as follows: (2)Solving and

Noted where is independent with . For and , then (21) can be reformulated to two subproblems: where and .

Solutions of (23c) are (3)Maximum penalty parameter

The Karush-Kuhn-Tucker conditions [21] of primal feasibility is the necessary condition to minimize a feasible point of constrained nonlinear programming. Thus, equation (21) with respect to is equivalent as follows:

A Bisection method (The bisection method is a successive nonlinear approximation method that keeps cut the interval (sign of of ) into 2 halves and check which half interval contains a root of the nonlinear function.) [22] is employed to find the root with respect to in (23e). That result is shown as: where . (4)Updating the Lagrangian multipliers ,

Finally, the proposed RAB based on the CM penalty criterion is summarized in Algorithm 1. Here, the computational burden is to find the root in (23f), but the finding root needs to run only once by using the bisection method. And the running time of the designed beamformer mainly depends on the iterative process of solving the weight vector, where the computational complexity is at each iteration.

Initialize
  Compute Ŝ using ((1)-(2)); Determine ã(θ,φ) using (5);
Repeat for t=1,…,L;
  Obtain m(t+1) and n(t+1) using (23d);
  Determine by solving (23f) using Bisection method;
  Obtain using (23a);
  Update using (23g) and (23h).
Until end for t=L;
  Output
  Final anti-jamming weight

5. Simulation Results

5.1. The Signal Model Based GPS Antenna Array

In Ansoft HFSS software, a four-element dual-polarized USAA with ports is designed to obtain the radiation gains of the antenna array at the GPS L1 band, which can model the GPS signal . In this work, the results are shown in Figure 2, which shows the normalized gain of the E/H-plane at 1.57 GHz when the first port is active and other ports are terminated in match loads. It is clearly found that the E-plane gain pattern is more symmetrical and has similar “apple-shape” radiation than the H-plane. Although the measured results are lower than the simulation due to the scattering structure errors in the testing environment, the maximal radiation efficiency is almost up to 77.26%. Then, we employ the maximal gains of the E/H-plane to obtain the signal . In addition, the deception jamming has a similar form with except for a different power, and the noise adopts a Gaussian noise. In short, the performance of dual-polarized USAA meets the requirement of a GPS antenna terminal, and the results are used to receive a GPS signal model .

5.2. Antijamming RAB Performance Comparison

To show the antijamming robustness of that novel beamforming, we list some exciting methods as the references, where there are the optimal beamforming, DL [6], LSMI [12], NCCB [13], WC [14], MVDR [16], and sample matrix inversion-based reconstruction (SMI) [23]. The interference-to-noise ratio (INR) is equal to -15 dB, and the noise obeys a Gaussian distribution . The snapshot is set to be , and the computational results are averaged by 200 Monte Carlo trials. The input SNR is set from -20 dB to 20 dB, and the simulation is running with the iteration from 10 to 100.

Firstly, the estimated directions of input signals are obtained as the a priori condition of the simulation experiment. MUSIC (multiple signal classification) [24], where the signal source both interferences and noise are shown in Figure 3. Here, it especially explains that the spectrum data are obtained by the covariance of different signals. In Figure 3, a solid line indicates the spectrum of interference from three peaks at . While the spectrum peak coincides with the peak of noise, then the DOAs are chosen as the directions of interferences. Then, the obtained DOAs and maximum radiation gains from HP and VP are introduced into formula (2) to generate the signal model.

Next step, the antijamming comparisons of beamforming based on the a priori DOAs condition are shown in Figure 4. Here, the  dB and  dB are chosen as other a priori conditions in the simulation experiment. It is clearly found that all valley points can be well located in the DOAs from the spectrum peaks shown in Figure 4. Note that the number of interference sources is while there are three interference sources [i.e., ] in Figure 4; this is because the azimuth angle is only chosen as the half number of the whole angle from -180° to 180°. By comparing the proposed method with LSMI and RAB, it can be observed that the antijamming null attenuates stronger than 38 dB in the jamming direction and the beamforming has a lower side lobe in all nonjamming directions. The conclusion of this experiment is that the proposed approach has good robustness to antenna array perturbation due to its constant modulus smooth characteristic in the beamforming.

Besides, for showing the comparison of beamforming power, the interference suppression (IS) levels are used as an important index of before and after antijamming work. Figure 5 plots IS levels by using the existing beamformers, where the black and blue graphs are pointed to the left black label and the right blue label, respectively. It can be found that the proposed method obtains the best antijamming level from  dB to 35 dB in all approaches. This is because the penalty variable can automatically enable the reconstructed beamformer to distribute more DOFs to suppress interferences in the spatial signal responses. Although its IS level is slightly worse than RAB, from 35 dB to 50 dB, the IS level gradually enhances along with the INR increase.

Finally, we list four scenarios as presented in Table 2, where the different example errors are shown to explore the robustness of antijamming beamforming by using the signal-to-interference plus noise ratio (SINR) comparison versus snapshots and SNR.

Scenario 1 considers the SINR comparison versus input SNR and the number of snapshots in random look direction errors. As shown in Figure 6(a), all existing beamformers can perform good SINR performance along with SNR increasing, but when SNR exceeds zero dB, RAB shows worse performance due to SINR decline. Additionally in Figure 6(b), we can clearly find that the proposed method performs best among all methods along with the growth of snapshots. This is because the used CM criteria can enable the beamformer to collect the max-stable power from look direction errors when SNR and snapshots gradually increase. But the SINR of LSMI and NCCB degrade along with the snapshots increasing. In scenario 2, we discuss the effects of antijamming beamforming in gain and phase perturbations. In Figure 6(c), along with the input SNR increasing, we can observe that the antijamming performance of all existing beamformers can show good robustness except RAB. In Figure 6(d), due to the gain and phase smooth perturbations, the proposed method performs the best SINR level versus the number of snapshots. Scenario 3 considers that the antenna element error is caused by the position mismatch. Figures 6(e) and 6(f) depict the robustness against antenna position error. We can observe that the SINR level of all beamformers can obtain the enhancement at  dB. Especially in the field of navigation, the proposed method performs best when the SNR is generally -30 dB. And when the number of snapshots exceeds 60, the proposed method outperforms other methods. In scenario 4, the coherent local scattering can cause the SV distortion in (12). The final results are shown in Figures 6(g) and 6(h), where the proposed method provides the best performance among all beamformers when the number of snapshots exceeds 50, while the SINR of NCCB, WC, LSMI, and RAB show a downward trend versus the input SNR increasing.

From the above analysis, it is illustrated that the proposed RAB can obtain good robustness of beamforming in all existing beamformers. However, the SINR of RAM suffers from performance degradation when the input SNR increases in all scenarios. Although the weight vectors of WC and NCCB can be solved by the CVX tool, the objective solution sometimes appears in the local optimization when the nonconvex nonlinear optimization is solved by the long running time. LSMI can give an efficient solution to find the optimal loading level, but an inaccurate loading sample leads to an imprecise Lagrangian multiplier. In contrast, the proposed RAB based on CM envelope response can obtain the robust adaptation of weight vector mismatch constraints in the spatial beamforming; further, the antijamming robustness of beamforming is also improved by the simple update rules and stable-fast convergence.

6. Conclusion

In this paper, we propose a novel beamformer via the constant modulus (CM) penalization criteria to improve the antijamming robustness of beamforming. Wherein, the CM envelope response is used to obtain the robust adaptation of weight vector mismatch by reconstructing the MVDR beamformer in collected snapshots, and a min-max penalized function is proposed to penalize the spatial responses of the array signal model. This can enable the reconstructed beamformer to allocate more DOFs to suppress the interferences. However, a nonconvex quadratic program problem is formulated in this beamformer, and then we use the ADMM optimization method to solve this quadratic problem. According to the SINR, some simulations show that the proposed beamformer can achieve better antijamming robustness in beamforming by competing with the exciting beamformers.

Data Availability

For example, “Data availability statement: The data that support the findings of this study are available on request from the corresponding author.” The (gains of antenna array) data used to support the findings of this study were supplied by (ORCID ID: 0000-0002-9623-8296) under license and so cannot be made freely available. Requests for access to these data should be made to (ORCID ID: 0000-0002-9623-8296).

Disclosure

This was performed as part of the employment of Xubao Sun because Xubao Sun as a funder was involved in the manuscript writing.

Conflicts of Interest

The authors declare that they have no conflicts of interest.