A novel adaptive beamforming scheme for array signal data processing

When the desired signal data exists in the array received data or the steering vector has a mismatch problem, the current traditional adaptive beamformers will suffer from the effect of the desired signal cancelation phenomenon, resulting in a sharp decline in performance. To address the occurrence of desired signal cancelation, an improved matrix projection-based efficient beamforming method is proposed. Firstly, based on spatial partitioning (SP) technology, a significant projection matrix for interference-plus-noise space (INS) is constructed. Secondly, using the constructed key projection matrix, the sample data covariance matrix is projected into the INS to achieve the goal of suppressing the desired signal data information. Finally, the weight data vector is calculated by Capon beamformer. The proposed algorithm does not require an iterative search for the optimal solution, which has the advantage of a small amount of calculation. Simulation experiments have verified that the proposed method has significant advantages in suppressing the desired data signals. Especially when the desired data signal has large power, the signal-to-interference-plus-noise ratio (SINR) of the proposed algorithm is better than that of the compared algorithms under the conditions of random directionality errors or local scattering errors between the desired signal and interference.

In literature [13], a robust beamforming algorithm based on the worst-case optimization (WCO) idea is proposed, which is robust to various multiple mismatch situations.For the steering vector error caused by the steering deviation, the steering vector cyclic iterative search beamforming method, which is defined as sequential quadratic programming (SQP) [14,15], iteratively searches for the optimal steering vector from a preset initial value in a characteristic subspace composed of a predefined observation sector, overcoming the steering vector error through optimizing the selection of orthogonal vectors of steering vectors to achieve the maximum output SINR.The literature [16] is based on the literature [14] and studies robust beamforming with as little as possible prior (LP) information.Specifically, based on steering vector estimation, a robust beamforming algorithm is proposed, making the proposed algorithm unnecessary for other prior information except for the predefined observation sectors.In literature [17], an efficient beamforming method based on sparse multi-input multi-output (MIMO) array and spatial filter bank (SFB) is proposed.
In literature [18], an adaptation beamforming called continual learning-based beamforming neural network (CL-BNN) method is addressed, the advantage of which is its satisfying performance in a time varying environment without large computation.In literature [19], an efficient beamforming method called alternating direction method of multipliers (ADMM)-based is proposed for dealing with all kinds of complicated situation, which obtains better result with low complexity.An anti-array mutual coupling algorithm [20] for covariance matrix reconstruction is proposed, which only requires the information of the desired direction of the signals without requiring other prior information to achieve a good anti-mutual coupling effect.
In this paper, an improved efficient matrix projection-based beamforming algorithm is proposed, which is different from the traditional projection beamforming technology.Firstly, the projection matrix of the INS is established by adopting the method of spatial division aforementioned in literature [16] and combined with the range of known desired signals.Then, the data covariance matrix is projected into the INS to suppress the information of the desired signals in the data covariance matrix.Finally, the matrix obtained after the projection and the corrected desired steering vector is substituted into the Capon beamformer to obtain a weighted vector.Among these steps, the "matrix projection" mentioned in the second step refers to projecting the steering vectors of each signal constituting the new reconstruction matrix into the target space, respectively.Furthermore, it is theoretically proven that the matrix projection can be easily implemented through the data covariance matrix.Compared with the currently recognized robust adaptive beamforming algorithms, the proposed algorithm has advantages in performance, which is also able to suppress the influence of observation errors, and converge faster with the number of sampled snapshots.

Array data model
Assume a uniform linear array (ULA) consisting of M isotropic sensors, with an array element spacing of half a wavelength, and the received signals being a far-field narrowband source signal.The desired signal and interference are uncorrelated.At moment k , the data received by the array can be represented as Among them, s(k) = s(k)d(θ s ) , i(k) , and n(k) stands for statistically independent the expected signals, interference, and noise, respectively.s(k) is the source signal, and d(θ s ) is the steering vector.Therefore, the steering vector ULA-based is described as follows Furthermore, the received signal covariance matrix is derived as follows stands for the direction-of-arrival (DOA) of these source sig- nals, θ 0 represents the desired source signal, and n I is the noise covariance matrix, and I is a well-known identity matrix.In practical applications, the theoretical covariance matrix cannot be realized.Therefore, the data sampling covariance matrix is used for calculation in the simulation experiments.And define the data covariance matrix obtained from the N samples of data as bellow The adaptive beamforming can be further expressed as follows where w = [w 1 , ..., w M ] T ∈ C M denotes the complex weighted vector of the sensor array.
Therefore, the optimal output SINR of beamforming can be calculated by using the following formula where R i+n stands for the covariance matrix of the noise-plus-interference.
Further, the maximization problem of SINR can be transformed into the following optimization scheme The solution to the above problem is the MVDR beamformer, also known as the Capon beamformer.The specific formula is as follows (1) x As can be found from Eq. ( 8) that the essence of the Capon beamformer is to only pass the desired steering vector without distortion, while these signals in other directions are regarded as interference, and its output should be suppressed, resulting in deep null in the output directional diagram at the position of strong interference points.In practice, the desired source signals will be suppressed as interference when there are undesirable conditions such as sensor array error and observation error, which is known as desired signal cancelation [21].Research shows that the greater the steering vector mismatch or the stronger the desired signal, the weaker the interference suppression ability of the Capon beamformer, the more serious the cancelation of the desired signal, and the worse the beamforming performance.Therefore, if the energy of the desired signal direction in the data covariance matrix can be suppressed, the performance of the Capon beamformer will be effectively improved.

The proposed beamforming method
Based on the spatial partitioning method [16], combined with the known region of the incoming signal, the projection matrix that divides the signal space and the interferenceplus-noise space can be constructed.By further analyzing the structure of the theoretical covariance matrix, Eq. ( 3) can be repressed as follows Description the covariance matrix projection expression as follows Among them, ⌢ P U 1 is composed by the region of the desired source signals.And Eq. ( 10) completes the projection transformation from the covariance matrix to the INS.It is known that in the projected matrix, the desired signal information will be suppressed.In the following simulation experiments, the theoretical covariance matrix will be replaced by the data covariance matrix to complete the projection operation.
Theoretical and experimental data indicate that the matrix obtained in Eq. ( 10) still requires diagonal processing to achieve normal noise levels, with the loading amount being the average power of the noise.The specific expression is as follows Among them, σ 2 n is the estimated value of noise power, which can be equivalently replaced by 10 times the minimum eigenvalue of the data covariance matrix.
By substituting R into the Capon beamformer, the weighted vector of the sensor array element can be achieved (8) From the above analysis, it can be concluded that the calculated amount of the proposed beamforming algorithm is mainly concentrated on the inverse operation of the matrix in Eq. ( 12), so the complexity of the proposed algorithm is O M 3 , which is the same as that of the Capon beamforming algorithm.In addition, the proposed algorithm does not require iterative search for the optimal solution operation, which greatly saves computing resources and is easier to implement in practical computers and hardware.

Results and discussion
In this section, the simulation conditions are set to use a ULA composed of 30 array elements, with a spacing of half a wavelength between the elements.Let the noise be independent additive Gaussian white noise with zero mean and standard deviation.The directions of the two interference waves are, respectively, − 20° and − 50°, and signalto-noise ratio (SNR) is both 30 dB.Assume the desired source signals and the interference signals are not correlated with each other, and the detected direction of the desired signal is 5°.All incident signals conform to the model of plane wave.In all simulation experiments, each data result is averaged from 500 Monte Carlo simulations.
In each simulation experiment, the proposed matrix projection beamforming algorithm is compared with the worst-case beamformer [13], the SQP beamformer [14], and the LP beamformer [16], among which the uncertainty set parameter used in the worstcase beamformer is ε = 0.3M in the literature [13].The main eigenvalues of matrix C of K = M − 1 are used in the proposed algorithm and taking = [0 o , 10 o ] , and 10 times of the minimum eigenvalue of R are selected as the diagonal loading amount σ 2 n .In each figure, the optimal output SINR is given, which is calculated from covariance matrix Eq. ( 6) of the true desired signal arrival direction and interference noise.The change relationship between the output SINR of beamformer and input signal to SNR of the desired signal is shown in Fig. 1.The number of sample snapshots used in this simulation is set to 60.It can be found that the output SINR of the proposed algorithm is lower than that of the other three comparison algorithms when the input SNR is negative.However, the performance of the proposed algorithm is significantly higher than that of the comparison algorithm when the input SNR is above 10 dB, but the difference of the optimal SINR is about 6 dB.
Figure 2 shows the conditions of output SINR of the four types of beamformers varying with the number of sampled snapshots.In this simulation, the desired signal input SNR is 20 dB.As shown in the figure, the convergence rate of the proposed algorithm is obviously higher than that of the algorithms compared.(12)

Experiment 2: Non-correlated Local Scattering Jamming Existed
The distributed signal or non-correlated local scattering signal is caused by multipath scattering of local scattering sources, which are widely used in the applications such as radar, sonar, space radio, and wireless communication.In this simulation experiment, we assume that the desired signal is the non-correlated scattering signal with time-varying spatial characteristics, and the steering vector is expressed as follows where s 0 (k) and s p (k) (p = 1, 2, 3, 4) are independent and identically distributed zero-mean complex Gaussian random variables, and the variance is 1.The value of θ p (p = 1, 2, 3, 4) follows Gaussian distribution with the mean value is θ s = 5 • , and the variance is 4° in each Monte Carlo simulation, but the value of θ p (p = 1, 2, 3, 4) does not change in different sampling snapshots in the same Monte Carlo simulation, and the value of s 0 (k) and s p (k) (p = 1, 2, 3, 4) will be different in different sampling snapshots.Therefore, the norm a(k) of the steering vector of the desired signal is time varying.It can be found from non-correlated scattering signal model, the signal covariance matrix R is no longer a rank-one matrix, and the output SINR is no longer given by Eq. ( 6), which should be given by the definition formula Eq. ( 14).
The optimal weight of the output maximizing SINR [22] can be achieved where P{•} denotes the principal eigenvalue of the matrix.
Figure 3 shows the change relationship between the output SINR of the four types of beamformers and the input SNR of the desired signal.The number of sampling snapshots is set to 60.It can be seen that when the input SNR is above 15 dB, the performance of the proposed algorithm is obviously higher than that of the comparison algorithms, but the difference of the optimal SINR is about 6 dB.
Figure 4 shows the conditions of output SINR of the four types of beamformers varying with the number of sampled snapshots.In this simulation, the desired signal input SNR is 20 dB.It can be found that the convergence rate of the proposed algorithm is obviously higher than that of the algorithms compared.A novel fast and efficient beamforming technique that relies on matrix projection is addressed in this paper.Specifically, firstly, under the premise of knowing the area covered by the desired source signals, the spatial partitioning technique is used to construct a projection matrix composed of INS.Secondly, the data covariance matrix is projected into the INS to effectively restrain the desired source signals, greatly improving the problem of desired signal cancelation in the case of excessive input SNR and steering vector mismatch.One of the advantages of the proposed algorithm is that its computational complexity is O M 3 without iterative operations, which is the same as that of Capon beamformer.Therefore, it is able to facilitate computer deployment and hardware implementation.The simulation results prove that when the input SNR is high, the performance of the proposed beamforming algorithm is superior to that of the comparison algorithms in the presence of mismatch issues.

Experiment 1 :
The Random Directivity Error of the Desired Signal and the Interference In the first experiment, assume the random directivity error is uniformly distributed in [− 4°, 4°], i.e., the true arrival direction of the desired signal is the arbitrary value in [1°, 9°] in each Monte Carlo simulation, the arrival directions of two interferences are the arbitrary values in [− 24°, − 16°] and [− 54°, − 46°], respectively, but the arrival direction values are not changed in different sampling snapshots in each Monte Carlo simulation.