Abstract

In this paper, we investigate the problem of direction of arrival (DOA) and polarization estimation for non-Gaussian signal in polarization-sensitive augmented coprime array. Instead of the second-order cumulant statistics, the fourth-order cumulant statistics of the received signals are used for parameter estimation because they can detect more information of the non-Gaussian signal. First, polarization-sensitive augmented coprime array is designed, where each sensor element is equipped with a pair of orthogonal electric dipoles. Furthermore, a low-complexity reduced-dimensional Capon algorithm which uses the fourth-order cumulant of the received array signal is proposed for DOA and polarization estimation. Only one-dimensional peak search is required by reconstructing peak search function. Theoretical analysis has proven the effectiveness of the algorithm, and simulation results demonstrate that the proposed fourth-order cumulant reduced-dimensional Capon algorithm outperforms the other algorithms.

1. Introduction

The research on directional of arrival (DOA), as a fundamental project in array signal processing, has attracted great attention from relevant scholars all around the world [1, 2]. In recent years, it is proved to be feasible in many fields, such as sonar, radar, and communication system. With the proposal and improvement of related algorithms for different scenarios, the estimate accuracy of DOA is gradually improved, the complexity is also gradually reduced, and more sources could be detected [37]. However, existing DOA estimation techniques are usually based on scalar sensor array structures with equal element spacing not greater than half-of-wavelength, which is also called uniform linear array [8]. This array structure can obtain unambiguous DOA estimates while suffering from heavy mutual coupling effects and low system resolution [9]. Simultaneously, only relying on DOA information received by scalar sensor array, it is hard to distinguish different signals impinging on the array from similar directions of arrival. Meanwhile, polarization-sensitive array (PSA) [10] is proposed to obtain another two-dimensional information combined with DOA information, which can distinguish different signals in higher dimensions, especially signals with similar incident angles. Many kinds of polarization-sensitive sensors have been designed for receiving different electromagnetic vector signals, like triad [11], cocentered orthogonal loop and dipole (COLD) [12], and crossed dipoles [13]. Relative algorithms are also proposed for measuring joint DOA and polarization estimates [1416].

In addition, coprime array (CA) [17, 18] is proposed by extending the array aperture and reducing mutual coupling effect sparse array, in order to improve the array performance. It solves the problem of compact array structure in uniform linear array, whose phase ambiguity can be eliminated by vectorizing the covariance matrix of the received signal to construct a virtual array. Based on CA, some high-performance DOA estimation algorithms have been proposed in the field of PSA [19, 20]. In [19], the CA is split up into two uniform linear arrays, and a minimum distance parallel factor algorithm is proposed to obtain an accurate estimate, though the advantage of high degrees of freedom (DOF) in CA is not exploited. In order to take full advantage of the large virtual aperture of the CA, a coarray interpolation method is proposed in [20]. However, it suffers from high computational complexity when solving the large stacked recovered covariance matrix.

Augmented CA (ACA) [21], as an improved array, is equipped with more virtual elements than other CAs after vectorization and a closed-form formula in its virtual array’s continuous part can be derived, which is more suitable for covariance matrix vectorization. Furthermore, compared with the covariance matrix which only uses second-order statistics of the received signals, fourth-order cumulants can demonstrate more information of non-Gaussian signals [22]. Besides, by constructing the fourth-order cumulants of the received signals, the array aperture of the equivalent virtual array is significantly extended [23], while the noise is suppressed in the fourth-order cumulant due to its Gaussian property [24]. As a result, fourth-order cumulant and its corresponding algorithms have shown great advantages in parameter estimation. However, the CA research on PSA and the corresponding algorithms has just begun, and more and more updated array structure designs and high-precision estimation algorithms are required in this field.

In this paper, we introduce the augmented CA to polarization-sensitive arrays in electromagnetic environment and design a kind of new sparse array called polarization-sensitive ACA (PSACA). Each sensor element is equipped with a pair of orthogonal dipoles, which are orthogonal to each other while the array element location is the same as ACA. On the other hand, a fourth-order cumulant reduced-dimensional Capon (FOC-RD-Capon) algorithm is proposed for non-Gaussian signals. First, the fourth-order cumulant is constructed by received signals, replacing the covariance matrix in traditional Capon algorithm, which can obtain a large amount of consecutive virtual elements. Benefited from the fourth-order cumulant, DOF are increased and virtual array aperture is extended in contrast to the virtual array by vectorization of covariance matrix. Second, considering the high calculation burden of three-dimensional spectral peak search for DOA and polarization estimates, a reduced-dimensional Capon algorithm is proposed, which only requires one-dimensional spectral peak search. The accuracy of the algorithm can be improved by shortening the search interval. The proposed algorithm is proved to have much lower complexity than traditional Capon algorithm. Numerous simulations illustrate the effectiveness of the ACA in polarization-sensitive scenarios and the proposed FOC-RD-Capon algorithm.

To summarize, the contributions of this paper are as follows:(1)We design the ACA in polarization-sensitive scenario, where each sensor element is a pair of orthogonal dipoles to receive electric field strength vector. Compared with traditional CA, ACA can achieve higher DOF.(2)We use the fourth-order cumulant signal instead of second-order cumulant signal for non-Gaussian signal, which can demonstrate more information to obtain high estimation performance. Moreover, the equivalent virtual array for the fourth-order cumulant signal has more consecutive elements.(3)We propose a low-complexity reduced-dimensional algorithm for polarization-sensitive augmented CA, avoiding three-dimensional peak search. Meanwhile, the proposed algorithm enjoys high estimation accuracy with multi-parameter autopairing.

The rest of the paper is arranged as follows.

Section 2 demonstrates the data model of the signal impinging on the PSA. Section 3 introduces the array structure of PSACA and its virtual array element model after employing the fourth-order cumulant. Section 4 elaborates the proposed reduced-dimensional Capon algorithm. Section 5 depicts the simulation results, and Section 6 concludes the paper.

Notations. We use lower-case (upper-case) bold character to imply vector (matrix). , , and denote the conjugate, transpose, and the conjugate transpose of a matrix or vector, respectively. denotes the Kronecker product, and represents Khatri–Rao product. represents matrix inverse. means phase operator. means diagonalization operator. means the integer set from to , and means the mathematical expectation of a vector or matrix.

2. Data Model

Consider that K far-field non-Gaussian narrowband signals impinge on the linear array with sensors where each sensor element consists of a pair of orthogonal electric dipoles to ensure vector receiving capability. The two dipoles with different polarization modes in an element are supposed to be distributed in the direction along with -axis and -axis, respectively, which are able to receive the electric component in the corresponding direction. The elevation angles of the signals are . The corresponding polarization information in the signal can be modeled as two polarization parameters: polarization auxiliary angle and polarization phase difference . As a result, the single snapshot data model of the received signal on the array is written as [16]where denotes the directional vector and is the location of sensor element. represents the non-Gaussian signal vector while is the white Gaussian noise vector. symbolizes the polarization-space steering vector of signal for orthogonal electric dipoles, which can be represented as [13]where denotes the coordinate matrix only containing DOA information, which is determined by the polarization mode of the electric dipoles, and is the polarization vector. Assume that snapshots are received during a period of time for parameter estimation, and the overall received signal is modeled as [25]where , are the directional matrix and electric matrix, respectively. denotes the non-Gaussian signal vector while symbolizes the additive Gaussian noise matrix.

Note that the number of the signals to be estimated in this paper is known. If the number is uncertain, methods such as matrix decomposition [26], information theory [27], or Geist's circle [28] can be used for estimation.

3. ACA and Virtual Elements

Based on the array model in [21], Figure 1 depicts the array structure of the ACA, which is composed of two subarrays. The two subarrays are both uniform linear arrays. Subarray 1, which is marked with white rectangles, is equipped with array elements with adjacent interval . Meanwhile, subarray 2 is marked with black circles, whose total array elements are with adjacent interval . The total element number of the ACA is because the first element of subarray 1 and subarray 2 is located at the same place. ACA is considered as the improvement of traditional CA because it keeps subarray 1 unchangeable while extending the array aperture of subarray 2 twice as that of CA by equipping sensor elements.

It can be obviously concluded from Figure 1 that the sensor element location of each subarray is thatwhere denotes the integers from to . Therefore, the overall position set of the sensor elements is described as .

According to the data model for PSA presented in Section 2, the fourth-order cumulant matrix of received signal with snapshots can be expressed as [24]

Assume that the signals are independent with each other, and fourth-order cumulant matrix can be considered as the covariance matrix of virtual received signal corresponding to an equivalent virtual array whose element location set is the difference co-array of the physical sensor elements, which can be defined as

Figure 2 demonstrates the virtual array structure after fourth-order cumulant extension for ACA and CA with the same physical sensor elements. There are both 12 physical sensor elements with in ACA and in CA. As depicted in the figure, CA has 25 consecutive virtual elements while ACA has 47 consecutive virtual elements, which is obviously more than that of CA. Benefited from the merit, ACA can achieve high DOF to detect more signals [29].

4. FOC-RD-Capon Algorithm

4.1. Algorithm Process Introduction

Capon algorithm can obtain multi-parameter estimates simultaneously by a multi-dimensional peak search, which has been widely used in DOA estimation. However, the direct introduction of Capon algorithm from scalar sensor array to polarization-sensitive array will result in the dramatic increase in search dimension, which causes the inevitable high-dimensional search process. To reduce the calculation burden of the Capon algorithm in polarization-sensitive arrays while improving the estimation performance for non-Gaussian signals, a FOC-RD-Capon algorithm is proposed in this section. The introduction will be presented as follows.

Traditional Capon search function for DOA and polarization estimation is expressed as [30]where is the covariance matrix of the received signal and denotes the array manifold for Q parameters. In PSA scenarios, the array manifold includes three parameters, DOA , polarization auxiliary angle , and polarization phase difference , which means three-dimensional peak search is required, but it is unrealistic in engineering application. This part introduces a kind of parameter elimination method for polarization dimensional reduction.

Replacing the covariance matrix by the fourth-order cumulant matrix , define , where is the directional vector for virtual array structure which is the difference coarray of physical array. Then, the peak search function for PSA fourth-order cumulant scenario is given by

Note that and a characteristic of Kronecker product is defined as [25]. The array manifold can be reconstructed as ; therefore, (8) is written as

Classifying the terms in the function according to the parameters, the terms only about DOA is expressed as . Besides, there is an identity that holds. Then, the peak search function has been transformed into the minimum finding under constraints that

In order to find the minimum value of (10), we calculate the partial derivative of about that

Let (11) equal 0 to calculate the final estimates that . Instead of solving the equations, we assume it as the form of eigenvalue and eigenvector pairs, where is the eigenvalue and is the corresponding eigenvector. Based on the analysis above, the peak search function (9) is expressed as

Because is a constant, (12) is actually the search about the eigenvalue [32].

Every peak calculation can be considered as finding the minimum eigenvalue about which is only constructed by DOA information. Polarization information is eliminated in Capon search function to tremendously reduce the computational complexity.

After finding peaks about the DOA estimates , the corresponding is obtained as

Polarization matrix is measured by the eigen decomposition of . Note that has 2 eigenvalue and eigenvector pairs, and the eigenvector whose eigenvalue is the smallest is the estimate . Ultimately, polarization estimates are obtained aswhere and denotes the element in vector .

The main steps of the proposed FOC-RD-Capon algorithm are summarized as follows:Step 1. Compute the fourth-order cumulant matrix of received signal according to (5).Step 2. Obtain the peak search function with (8) and reconstruct it according to the characteristic of Kronecker product (9).Step 3. Establish the search dictionary to calculate with each dictionary element .Step 4. Perform eigen decomposition of and find its smallest eigenvalue .Step 5. Find peaks in (13) to estimate .Step 6. Calculate the eigenvector corresponding to the smallest eigenvalue in .

4.2. Discussion

The calculation burden of the FOC-RD-Capon algorithm mainly results from the following steps. Assume that there are sensor elements in the ACA. Computing the fourth-order cumulant of the received signal requires . The inverse of fourth-order cumulant matrix needs the complexity of . Every peak search consists of for -times peak search. Therefore, it can be summed up that the proposed FOC-RD-Capon algorithm is composed of the complexity of .

Meanwhile, the traditional Capon algorithm without dimensional reduction requires three-dimensional peak search. Computing the fourth-order cumulant and inverse of it requires totally , which is the same as the proposed algorithm. However, the three-dimensional peak search needs the complexity of for -times search in each dimension. To sum up, the approximate calculation burden of the traditional Capon algorithm is , which is obviously higher than that of the proposed algorithm. Figure 3 shows the comparison of the complexity of the two algorithms with the number of sensor elements.

In addition, ACA is employed for the proposed algorithm due to the following merits:(1)ACA is a kind of sparse array, where array aperture is extended and mutual coupling effect is eliminated. In practice, the coprime numbers and can be chosen according to the requirement to satisfy different scenarios. Benefited from the sparse structure, the DOA and polarization estimation results are tremendously improved.(2)The proposed algorithm uses the fourth-order cumulant statistic instead of the second-order cumulant, where the received signal can be regarded as the signal received by the equivalent virtual array. The virtual array is the difference coarray of physical array. Based on the conclusion and the analysis in Section 3, ACA can achieve more consecutive elements than CA, which means being able to detect more signals.

Cramér–Rao bound (CRB) is often used as the standard error which can be calculated. We derive the CRB formula to evaluate the RMSE performance of the algorithms on the designed ACA. Define , , and , and the formula of CRB for polarization-sensitive ACA is demonstrated as [32]where , , , , , . is the power of noise.

5. Simulation Results

Numerous simulations are performed to verify the effectiveness of the designed ACA and proposed FOC-RD-Capon algorithm, where root mean square error (RMSE) is employed as an evaluation standard to judge the performance. RMSE is defined aswhich is the RMSE of parameter . denotes the times of independent Monte Carlo simulations. is the actual parameter, and is its estimate in simulation. The two sinusoidal waveform signals are impinging on the ACA, and the DOA and polarization parameters are set as , . The array structure of ACA is , and the total number of physical array elements is 12, which means subarray 1 is equipped with sensor elements and subarray 2 has 8 sensor elements.

5.1. Scatter Plot of the Proposed Algorithm

Figures 4 and 5 depict the scatter plots of DOA and polarization estimates, where 100 independent Monte Carlo simulations are performed. Simulation environments are set as follows: the signal-to-noise ratio (SNR) is 15 dB and the number of snapshots is 300. As is exhibited in the figures, DOA estimates and two polarization estimates are accurately measured and different parameters are correctly paired. In addition, it can be noticed that the variances of polarization estimates are higher than DOA estimates.

5.2. RMSE Performance of Different Algorithms versus SNR

Figures 68 demonstrate the RMSE performance of the proposed FOC-RD-Capon algorithm versus SNR, where CRB curve is given as the standard. Snapshots are fixed as 200 while SNR varies from −10 dB to 20 dB. Compared with the propagator method (PM) [33], estimating signal parameter via rotational invariance techniques (ESPRIT) [34], and parallel factor (PARAFAC) [35] algorithms, the proposed algorithm has the lowest RMSE, and its curve is closest to the CRB. Meanwhile, the RMSE of all the four algorithms decreases as the SNR improves, which indicates its influence to the estimates.

5.3. RMSE Performance of Different Algorithms versus Snapshots

Figures 911 show the RMSE performance versus snapshots. Similar to Section 5.2, PM, ESPRIT, and PARAFAC algorithms are used for comparison and CRB curve is depicted as the standard. With the increase of snapshots, the estimation accuracy of all the algorithms is improved. Meanwhile, the proposed FOC-RD-Capon algorithm has the least RMSE whatever snapshots are, which shows its great performance.

6. Conclusion

In this paper, we make use of the inherent information in fourth-order cumulant of non-Gaussian signals, proposing a FOC-RD-Capon algorithm for PSACA. Compared with the scalar sensor arrays, each sensor element of PSASA is equipped with a pair of orthogonal electric dipoles to receive vector signals, which is able to obtain joint DOA and polarization estimates. The proposed algorithm reduces the search dimension of Capon method from three dimensions to one dimension, tremendously reducing the calculation burden. According to the analysis and the numerous simulations, the effectiveness of the proposed FOC-RD-Capon algorithm is confirmed.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by China NSF (61971217, 61971218, and 61631020), Jiangsu NSF (BK20200444), and National Key Research and Development Project (2020YFB1807602).