Abstract
The increase in element spacing can increase the aperture of the array and improve its resolution performance. However, phase ambiguity will occur when the array element interval is larger than the minimum half wavelength of the incident signal. The three acoustic velocity components of the acoustic vector are ingeniously constructed into a new kind of quaternions because of the special structure of the acoustic vector sensor array, and the rough estimation of the direction of arrival (DOA) is obtained using the rotation relationship between the subarray steering vectors corresponding to quaternion data. The rough estimate is used to resolve the phase ambiguity of the spatial phase difference between the array elements, and the high-precision DOA estimation of the signal can be obtained. Simulation results show that the method is effective.
1. Introduction
The use of electromagnetic wave for communication or target location has become mainstream. However, the acoustic wave has more advantages than the electromagnetic wave in some special circumstances. For example, sound waves can travel around obstacles, such as mountains, rivers, and forests, and the sound source detection system can effectively detect the information of the target. The loss of acoustic wave propagation is small in water. Thus, the wave can spread for a long distance. Sonar is a key equipment in underwater communication. The acoustic vector sensor is the core part of new sonars. A single acoustic vector sensor not only can output the acoustic pressure information of the scalar sensor but also can measure the acoustic velocity information of each point directly and synchronously [1, 2]. The acoustic vector sensor consists of an acoustic pressure component sensor and a speed sensor in three-speed directions. The vibration speed sensor may be an acceleration, a displacement, or a sound pressure gradient sensor [3, 4]. In 2010, Wu et al. obtained the near-field source model of the acoustic vector sensor by using the wave equations of the acoustic and particle velocity fields [5]. In 2012, Song et al. proposed to place the components of the three speed directions at the coordinate origin and an acoustic pressure component on one of the coordinate axes; the former is used to estimate the angle of the sound source, and the latter is used to estimate the distance information of the sound source [6]. In the array direction finding, the resolution of the angle with large array spacing will improve correspondingly. However, when the array spacing increases to more than half the wavelength of the incident signal, the pattern of the array will show grating lobes, which will lead to the ambiguity problem of parameter estimation [7, 8]. Many ways can be used to solve the ambiguity problem. Wang et al. proposed a virtual baseline method of resolving ambiguity based on the virtual ring [9]. Wong and Zoltowski proposed a uniform plane array method for solving ambiguity [10]. A nonuniform array method for solving ambiguity problems was proposed in [11], and the problem of ambiguity resolution by virtual array expansion was studied in [12]. In fact, the general mode of the ambiguity-solving method is to obtain two different arrays through array arrangement, virtual reconstruction, or data processing. One nonsparse array without ambiguity can enable rough parameter estimation, and the other sparse array with phase ambiguity can realise accurate parameter estimation. The accurate and unambiguous estimations of the direction of arrival (DOA) are obtained by using the rough estimation of the nonsparse array to solve the phase ambiguity of the sparse array.
In recent years, array signal processing technology based on quaternion [13] has attracted increasing attention. In [14], the singular value decomposition of the quaternion matrix was studied, and a polarisation wave separation algorithm based on the real signal quaternion model was proposed. In [15, 16], the quaternion and biquaternion music algorithms for two- and three-component vector sensor arrays were proposed, respectively. In [17], the DOA estimation for the six-component electromagnetic vector sensor array was obtained using a quad-quaternion music operation. In [18], a quaternion ESPRIT algorithm based on the translation-invariant cross dipole array was proposed. These studies show that the parameter estimation algorithm based on the multivariate number has improved model error tolerance by using the constraint relationship between the different imaginary parts of the multivariate number and the strong constraint condition between the orthogonal vectors of multivariate number. In this work, a quaternion–ESPRIT algorithm for solving the direction-finding ambiguity based on the acoustic vector is proposed. This algorithm is different from the previous ambiguity-solving method. The DOA is roughly estimated using the rotation invariant relationship between the quaternion subarrays and is used to solve the phase ambiguity of the spatial steering vector. The phase ambiguity number vector is attained using the relationship between the estimated value of the phase difference between sparse elements with phase ambiguity and the rough value without ambiguity but with low accuracy. The high-accuracy estimated value without ambiguity is obtained in the end. This method has a small amount of computation and can enable high-precision DOA estimation without ambiguity.
2. Acoustic Vector Array Model
A uniform circular array is composed of elements; is the radius of the circular array, as shown in Figure 1. Phase ambiguity occurs when the interval between the adjacent elements is greater than the half wavelength, that is,, where is the minimum wavelength of the incident signal. The array elements are the acoustic vector sensors; each of these elements is composed of three orthogonal velocity sensors along the -axis, -axis, and -axis, which induce the velocity component of the incoming wave signal in the -axis, -axis, and -axis directions, respectively.
Assuming that narrow-band, zero-mean, and stationary signals are incident from the far-field onto the uniform circular array shown in Figure 1, the incident direction of the th signal source is , , and is the angular coordinates of the th element. The velocity components along the -, - and -axes of the th unit energy signal in the Cartesian coordinate system can be expressed as:
3. Quaternion Method
3.1. Mathematical Model of Quaternion Method
The -axis velocity component , -axis velocity component , and -axis velocity component on the same array element are expressed as quaternion models as : where and are the imaginary parts of a quaternion.
The spatial steering vector of the array consisting of the phase difference between the array elements and the reference origin on the circumference is as follows:
In the formula (3) above, represents the angular coordinates of the th element on the circumference, and is the wavelength of the th incident signal.
The received data of the array at time can be expressed as:
where is the th incident signal, represents the Kronecker product, is the steering vector matrix at time , is the incident signals, and is the additive white Gaussian noise at time .
Then, the received data of the array at time can be represented as: where is the interval between two snapshots, is the frequency of the th signal, is the steering vector matrices at time , is the additive white Gaussian noise at time , and
So the whole received data can be written as: where is the whole steering vector matrix, and is the whole additive white Gaussian noise.
3.2. Quaternion Ambiguity Resolution Algorithm
The quaternion correlation matrix of the received data is Rwhere denotes the mathematical expectation, represents the conjugate transpose, , , indicates the white noise power, and is the identity matrix.
The signal subspace and noise subspace can be obtained by quaternion eigendecomposition operations of the data correlation matrix , and are the first rows and the last rows of , respectively. According to the subspace theory, it can be got as follows: where is the nonsingular matrix; the matrices , , and are the respective estimates of , , and ; denotes the pseudoinverse operator; and represents the inverse operation. Let , then formula (10) can be expressed as: where the eigenvalue of constitutes the estimation of , which contains the frequency information of the source, and is the corresponding eigenvector. The matrix factorisation of the quaternion is different from that of the long vector, that is, the eigendecomposition operation of the quaternion matrices cannot be performed directly. After the quaternion matrix is obtained, it is usually used to extract one real part and the three imaginary parts of the quaternion matrix and reconstruct the complex term to form the complex adjoint matrix. The properties of quaternion matrices are obtained by studying the properties of the plural adjoint matrix. Therefore, the corresponding information should be extracted from the quaternion matrix . Accordingly, the estimate of the steering vector is obtained as follows: where is the estimation of .
3.2.1. Rough Estimate of DOA
The eigenvalue of can be expressed as follows:
where are the frequencies of the corresponding signals.
The frequency of the th signal source can be estimated as follows: where denotes the operation of the extracting phase and is the th row and th column element of the matrix .
According to formulas (2) and (12), the steering vector can be written as another form:
Among them, and satisfy the following relationships: where
According to formula (17), the rough estimate of DOA without ambiguity is obtained. where is the th row and th column element of the matrix .
3.2.2. Ambiguity Resolution Algorithm
According to the unambiguous rough estimate of and , the rough phase difference estimation can be obtained: where , , , , , represents the submatrix extracted from the first row to row of , and represents the submatrix extracted from the second row to row of , with
Let where and denote vectors composed of elements from row 1 to row and from row 2 to row in the column of the matrix , respectively.
From formula (21), the estimate of spatial steering vectors can be written as: where denotes the dot division which is the division corresponding elements of two vectors.
By taking the phase on both sides of equation (22), the estimated value containing the ambiguous phase difference can be expressed as:
The ambiguous vector can be solved by the following formula:
The expression represents the vector composed of the smallest integer satisfying the condition and is the possible vector of ambiguity numbers.
By substituting the obtained ambiguous vector into formula (25), the accurate unambiguous phase difference can be obtained as follows:
From (25), the accurate DOA estimate without ambiguity is estimated as follows:
4. Long Vector Method
4.1. Mathematical Model of the Long Vector Method
The received data of the acoustic velocity subarray in the -axis, -axis, and -axis directions can be expressed as: where is the incident signal; , , and are, respectively, the additive white Gaussian noise along the -axis, -axis, and -axis directions; and , , and are the subarray steering vector along the -axis, -axis, and -axis directions.
, , and can be further expressed as follows: where , , and , with the spatial steering vector of the array; its expression is the same as formula (3).
The relationship between the three subarray steering vectors is: where
The received data of the acoustic velocity subarray along the -axis, -axis, and -axis directions are arranged in sequence to form a long vector; then, the received data of the array at and time can be expressed as: where is an array steering vector at time , denotes the Kronecker product, is the th element of , is the receiving noise of the whole array at time , is the receiving noise of the whole array at time , is the noise received by the -axis vibration velocity subarray, and are the noises received by the -axis and -axis vibration velocity subarrays, respectively, and the expression of is the same as formula (1). is an array steering vector at time .
In summary, the whole received data can be written as:
where , .
4.2. Long Vector Ambiguity Resolution Algorithm
The plural received data correlation matrix is:
The signal and noise subspace and can be obtained by the complex domain eigendecomposition of the data correlation matrix . The signal subspace is divided into front and back halves, and the estimation of signal frequency and steering vector can be obtained by the ESPRIT method. The process is similar to Section 3.2.
where
4.2.1. Rough Estimate of DOA
According to formulas (29) and (34), the following estimates can be obtained:
The rough estimate value of DOA without ambiguity is given by: where and are the th row and th column element of . The ambiguity resolution algorithm is the same as in Section 3.2.2.
5. Performance Analysis of the Algorithm
5.1. Scope of Application Analysis
The quaternion ESPRIT parameter estimation method is convenient for estimating the roughly estimated angle of arrival of the signal. It ingeniously combines three velocity sensor subarrays into a quaternion subarray by using the rotation relationship between the quaternion subarrays. Comparing formulas (18) and (36) reveals that the quaternion method has evident advantages in calculating the rough estimation of the angle of arrival. The rotation invariant relationship between direct and delayed sampling data can be used to estimate the parameters for narrow-band far-field signals with different frequencies; in this scenario, no requirement is imposed for the array configuration. The translation invariant relationship between subarrays should be used for parameter estimation of the same-frequency signal; this approach is only suitable for uniform array and results in array aperture loss. For uniform L-array, the eigendecomposition should be performed using the translation-invariant relationship in the - and -axis directions to obtain the spatial steering vector in these directions.
5.2. Computational Complexity Analysis
The following is a comparison of the computational cost of the algorithms based on the two models. Given that the amount of algorithmic computation depends on hardware and software, we only focus on the computational complexity of the covariance matrix of the algorithm. Without loss of generality, an array of acoustic velocity vector sensors is considered. In data processing, snapshots are used to estimate the data covariance matrix. The data covariance matrices based on the quaternion and plural model are expressed as follows:
Each of the array elements is composed of three orthogonal velocity sensors along the -, -, and -axis directions. Quaternion is , contains quaternion elements, storage matrix requires machine memory units, and the multiplication of two quaternions is equivalent to 16 real multiplications and 12 real additions. The operation of the matrix requires real multiplications and real additions. Thus, real multiplications, real divisions, and real additions are required for the operation of the matrix . The plural is , contains complex units, storage matrix requires machine memory units, and the multiplication of two plurals is equivalent to 4 real multiplications and 2 real additions. The operation of the matrix requires real multiplications and real additions. Thus, the operation of the matrix requires real multiplications, real divisions, and real additions. Under the same array conditions, the storage units needed for the calculation of the quaternion covariance matrix and the real number addition, multiplication, and division operations are less than those for the long-vector method.
6. Simulation Experiment and Analysis
The proposed algorithm in this paper is applicable to uniform linear array, uniform circular array, and so on. Without loss of generality, the simulation experiment takes uniform circular array as an example. Two far-field, narrow-band, non-Gaussian stationary sound signals are incident into the uniform circular acoustic vector sensor array shown in Figure 1. The receiving array is composed of 9 array element s; the incident signal parameters are set as follows: , , digital frequencies , the interelement with ambiguity is , 200 Monte Carlo experiments, and the noise is zero-mean additive white Gaussian noise.
6.1. Simulation Experiment
Figures 2 and 3 are the scatter diagrams of the estimated results, the number of snapshots is 50, and the SNR is set to 20 dB. From the intuitive comparison of the two graphs, it can be seen that the elevation and azimuth angles of the quaternion deambiguity algorithm are more concentrated around the actual value, while the angle estimation of the long vector de-ambiguity algorithm is more scattered. It is shown that the quaternion de-ambiguity algorithm has a little improvement in parameter estimation accuracy compared with the long vector de-ambiguity algorithm, but the calculation amount and storage space of the quaternion de-ambiguity algorithm is much less than that of the long vector de-ambiguity algorithm.
Figures 4 and 5 plot the root mean square error (RMSE) of elevation angle and azimuth angle, respectively estimated by long vector and quaternion de-ambiguity method, at different snapshot number levels. The SNR is set as 10 dB. It can be seen from Figures 4 and 5 that the RMSE of the elevation and azimuth angles of the two algorithms decreases with the increase of the number of snapshots. The RMSE performance of the quaternion-based de-ambiguity algorithm is better than that of the long vector de-ambiguity algorithm in the snapshot number range from 25 to 75.
Figures 6 and 7 show the success probability of the elevation angle and azimuth angle, respectively, estimated by the long vector and quaternion de-ambiguity method, at different snapshot number levels. The SNR is set as 10 dB. It can be seen that the success probability of elevation and azimuth angles of the two algorithms increase with the increase of the number of snapshots. The success probability of the elevation angle and azimuth angle based on the quaternion-based de-ambiguity algorithm is better than that of the long vector de-ambiguity algorithm in the snapshot number range from 25 to 60 and from 25 to 55, respectively.
6.2. Performance Analysis
The proposed method needs a quaternion eigendecomposition and a complex eigendecomposition. The eigenvalue of the plural eigendecomposition can be used to calculate the rough estimate of the direction of arrival of the signal, and the eigenvectors constitute the estimation of the spatial steering vector. The parameters are automatically paired, and no additional pairing operations are required. However, the long-vector method needs to divide the signal subspace into three subarrays corresponding to the steering vectors along the , , and -axes. The complex eigendecomposition of the full data must be performed, and two additional complex eigendecompositions are needed. The eigenvectors of the two latter eigendecompositions are used to estimate the steering vector and the direction of arrival of the signal. The two eigendecompositions are performed separately, and additional pairing operations are required.
7. Conclusion
The proposed phase ambiguity resolution method of the sparse array can effectively improve the estimation accuracy of parameters and the flexibility of array placement. A set of rough and unambiguous estimates of DOA can be given very conveniently through the ingenious combination of three acoustic velocity sensors into quaternion data. Thus, the phase ambiguity of the sparse array steering vector can be solved, and high-precision DOA estimates are herein obtained. The vector orthogonality between the quaternion data components is fully utilised, and the automatic matching of the eigenvalues and eigenvectors does not require additional matching operations. The quaternion method is less computationally complex than the long-vector method. ESPRIT of the long vector can be applied to arbitrarily distributed arrays in space and has a wide application range.
Data Availability
All data included in this study are available upon request by contact with the corresponding author.
Conflicts of Interest
The author(s) declare(s) that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China under contract (61772398, 61972239), the Key Research and Development Program Projects of Shaanxi Province (2019SF-257, 2020GY-024), the Opening Foundation of Shaanxi University of Technology Shaanxi Key Laboratory of industrial Automation (SLGPT2019KF01-15), and the Science and Technology Program of Hantai District (2019KX-21). The authors would like to thank the anonymous reviewers and the associated editor for their valuable comments and suggestions that improved the clarity of this manuscript.