Self-Calibration for Sparse Uniform Linear Arrays with Unknown Direction-Dependent Sensor Phase by Deploying an Individual Standard Sensor

: Calibration of the unknown direction-dependent (DD) sensor phase and aliasing-free directions of arrival (DOA) estimation for sparse linear arrays are difﬁcult tasks. In this work, we deploy an individual standard sensor with a known sensor phase response along the axis of the uncalibrated sparse linear array, a self-calibration method is proposed, in which the unknown DD sensor phase and the aliasing-free DOAs are both estimated. The proposed method is realized with a two-step scheme. In the ﬁrst step, the sensor phase is eliminated by the Kronecker product of the covariance matrices in two different frequency bins, and the frequency difference satisﬁes the spatial Nyquist sampling theorem. Then, the DOAs can be robustly estimated without the inﬂuences of grating lobes and unknown sensor phase parameters. In the second step, the steering matrix is estimated with the known phase parameters of the deployed standard sensor. Then, the DD sensor phase is extracted from the steering matrix using the DOAs obtained in the ﬁrst step. Hence, the disadvantages of iteration-based strategies in conventional calibration algorithms (e.g., local minima convergence) can be avoided. The performance of the proposed method is evaluated using simulation data and compared with that of Cramer–Rao bounds.


Introduction
Array signal processing techniques, such as the high-resolution direction-of-arrival (DOA) estimation [1,2] and adaptive beamforming [3], are sensitive to the unknown sensor responses. Calibration of sensor response is an important step to improve not only the DOA estimation but also the following adaptive beamformers: hydrophones ( [4], Chapter 4), the fiber optic strain sensor [5] in an acoustic system, or antennas [6] in an electromagnetic system. Self-calibration is a type of calibration algorithm that can simultaneously estimate source DOA and array response vector parameters, such as coupling errors [7,8] and unknown sensor phase responses [9]. In applications such as speech or underwater acoustics, the bandwidth of the signals generated by speakers or sources typically exceeds 1 oct [10,11]. Accordingly, the half wavelength of the high frequency components is smaller than the interelement spacing, and the sensor array can be regarded as a sparse one. The coupling error between sensors can be ignored because of the large interelement spacing, whereas the unknown sensor phase occupies the array response mismatches.
In contrast, it is difficult to distinguish whether the phase mismatch of the array response is caused by the DOA estimation errors or the unknown sensor phase under the DD model. Accordingly, the DOAs must first be obtained on the basis of the imperfect array response [24] or estimated with the partly calibrated sensors [23]. A DOA and steering vector estimation algorithm in [24] is a classical calibration algorithm for the DD sensor phase models. This algorithm alternately estimates the DOAs and the steering vectors with an iterative strategy. By contrast, the number of precalibrated sensors should be larger than that of sources to obtain the initial value of the DOAs, and the convergence heavily depends on the initializations. A weighted alternating least squares algorithm in [25] is also an iteration-based calibration algorithm for DD sensor parameters. The iteration easily converges to local minima due to the poor initialization. The instrumental sensor method in [26] is a noniterative algorithm that can estimate the DD sensor phase with instrumental sensors (i.e., precalibrated sensors). However, the number of precalibrated sensors should be also larger than that of the sources to make the unknown DOA and sensor responses solvable. In our previous work [27], the DD sensor phase parameters are extracted from the steering matrix without iterations. However, the DOAs should be first estimated with two precalibrated sensors.
These existing algorithms are not suitable for sparse uniform linear arrays (ULAs) because of the ambiguities in the spatial spectrum [28,29]. One disambiguation scheme for addressing this issue is to apply dual-size spatial invariance arrays; the ambiguous directions are eliminated with half-wavelength-spaced subarrays [30,31]. Another disambiguation scheme is the multistage DOA estimation scheme [32], in which the outputs of sensors are decomposed into multiple subbands and the influence of aliasing components is minimized due to the different structures of the aliasing components in each subband. However, the robustness against the unknown sensor response is not considered. To the authors' knowledge, only a few studies have been conducted on self-calibration methods for sparse linear arrays.
This study proposes a self-calibration method for sparse ULAs with an unknown DD sensor phase. An individual standard sensor is first deployed along the axis of the sparse ULA. Then, a Kronecker product of covariance matrices in two different frequency bins is conducted. In this manner, the unknown and frequency invariant DD sensor phase is eliminated. Additionally, the frequency difference of the two different frequency bins satisfies the spatial Nyquist sampling theorem and the spatial aliasing is also avoided. The sensor phase is then calibrated with the constant modulus (CM) algorithm by using the estimated DOAs and the known phase response of the deployed standard sensor. In the proposed method, iterations are not required, and the disadvantages of iteration-based strategies in conventional calibration algorithms (e.g., local minima convergence) can be avoided. Moreover, the deployment of an additional standard sensor is easy to implement in real tasks.

Sensor Phase Models
In this work, a ULA consisting of M sensors with an unknown DD sensor phase is considered. We assume that the sensor phase response in the array does not have a uniform increment across the array but is frequency invariant [25] when processing low frequency signals, especially for hydrophones in sonar arrays ( [4], Chapter 4). Thus, the direction-dependent sensor phase Φ can be modeled as follows: where ϕ m,k is the phase response of the mth sensor in the direction of the kth signal. A standard acoustic sensor with an accurate phase response is deployed along the axis of the uncalibrated ULA to calibrate the sensor phase parameters in the ULA. Thus, a new ULA with M = M + 1 sensors is formed, and the deployed sensor becomes the first sensor in the new ULA, as shown in Figure 1. Without losing generality, we assume that the phase parameter of the deployed standard sensor is ϕ 1,k = 0, k = 1, · · · , K, where K is the number of signals.
are the additive Gaussian noise of the deployed standard sensor and the original uncalibrated ULA, respectively.
θ is the unknown direction of the th k incident signal; and  is the Hadamard product. Considering that 1, 0, 1, , On the basis of (2), the data covariance matrix in the th q frequency bin is represented as follows:

Signal Model
K broadband and uncorrelated signals impinge on the ULA from the far field. The half wavelength of the incident signals is smaller than the interelement spacing, i.e., λ/2 <d, where d is the interelement spacing. Thus, the sensor array is regarded as the sparse ULA.
The sparse ULA and the deployed standard sensor synchronously acquire observation data. The observed time interval is partitioned into L sections without overlapping. Then, a Q-point FFT is applied to each section. The frequency domain observation in the lth section and the qth frequency bin can be expressed as follows: where x 1,l ( f q ) and x ' l ( f q ) are the observations of the deployed standard sensor and the original uncalibrated ULA, respectively; s l ( f q ) = [s 1,l ( f q ), · · · , s K,l ( f q )] T contains the signal spectra; and n 1,l ( f q ) and n ' l ( f q ) = [n 2,l ( f q ), · · · , n M ,l ( f q )] T are the additive Gaussian noise of the deployed standard sensor and the original uncalibrated ULA, respectively.
The steering matrix A q in (2) is A q = P • A q , where P = e jΦ ; A q is the ideal steering matrix in the qth frequency bin; A q m,k = exp j2π f q (m − 1)d cos θ k /c , m = 1, 2, · · · , M , c is the wave speed; θ k is the unknown direction of the kth incident signal; and • is the Hadamard product. Considering that ϕ 1,k = 0, k = 1, · · · , K, we have A q On the basis of (2), the data covariance matrix in the qth frequency bin is represented as follows: where Π q = diag σ q , and σ q = [σ 1,q , · · · , σ K,q ] T denote the powers of the sources in the qth frequency bin, η q is the noise power in the qth frequency bin; I is an M × M identity matrix, and (•) H stands for the conjugate transpose. In practice, the data covariance matrix R q is calculated from the sample covariance matrix as follows:R whereR q represents the sample covariance matrix.

Proposed Method
A robust DOA is first estimated without the influence of spatial aliasing or sensor phase uncertainties. Then, the sensor phase is calibrated using the estimated DOA values.

Robust DOA Estimation without Spatial Aliasing
Spatial aliasing is generated when the interelement spacing is larger than the half wavelength of the signal of interest (SOI). To overcome this issue, a frequency difference operation is conducted and the frequency of the SOI is pulled to a reference frequency f r , which satisfies the spatial Nyquist sampling theorem (i.e., f r = c/(2d)).
After the Q-point FFT operation, the index of the frequency bin for f r can be calculated as q r = round{Q f r / f s }, where round{•} is the rounding operation and f s is the sample rate.
Considering the properties of the Kronecker product, i.e., (AB) ' , we can obtain the following expression: where ⊗ is the Kronecker product, * is the conjugate operation, and B q,q+q r = A q ⊗ A * q+q r . In (5), the noise power η q and η q+q r can be estimated by the average of the smallest M − K eigenvalues of the sample covariance matricesR q andR q+q r , respectively, because the noise subspace estimation from the eigendecomposition of the sample covariance matrix is not affected by the unknown DD sensor phase.
The ((m 1 − 1)M + m 2 )th row and the ((k 1 − 1)K + k 2 )th column of matrix B q,q+q r can be expressed as follows: where 1 ≤ k 1 ≤ K and 1 ≤ k 2 ≤ K. Given the inconsistency of the order of the DOAs in A q and A q+q r , the k 1 th column of the matrix A q corresponds to the direction of the k 1 th signal and the k 2 th column of the matrix A q+q r is the steering vector of the k 2 th signal.
where θ k 1 ,k 2 = arccos f q cos θ k 1 − f q+q r cos θ k 2 /(− f r ) + 2k , and k is an integer that makes the inequality f q cos θ k 1 − f q+q r cos θ k 2 /(− f r ) + 2k ≤ 1 hold. In this case, we have θ k 1 = θ k 2 = θ k . Thus, θ k 1 ,k 2 is also equal to θ k , which is the direction of the kth signal. Given that exp j(ϕ m,k 1 − ϕ m,k 2 ) = 1, the influence of the sensor phase is canceled. At this time, the steering matrix in (7) can be regarded as the equivalent steering matrix at the reference frequency in the direction of θ k , k = 1, · · · , K.
Considering that the reference frequency f r = c/(2d), the negative effect of the spatial aliasing is avoided. CASE II: k 1 = k 2 In this case, the direction of θ k 1 ,k 2 is regarded as the cross term between the incident signals in θ k 1 and θ k 2 . A total of K 2 − K cross terms are generated and θ k 1 ,k 2 , k 1 = k 2 varies with the frequency.
Then, we collect the to form an equivalent covariance matrix at the reference frequency as follows: where J is an M × M 2 section matrix with [J] m,(m−1)M +m = 1, and the other entries are equal to zero. The calculation of the equivalent covariance matrix at f r is shown in Figure 2.
the inequality ( ) ( )  (7) can be regarded as the equivalent steering matrix at the reference frequency in the direction of , 1, , k k K θ =  . Considering that the reference frequency (2 ) r f c d = , the negative effect of the spatial aliasing is avoided.
In this case, the direction of Then, we collect the to form an equivalent covariance matrix at the reference frequency as follows: where J is an , and the other entries are equal to zero. The calculation of the equivalent covariance matrix at r f is shown in Figure 2.  Given the Kronecker product between σ q and σ q+q r in (5), the number of signals (i.e., the incident signals and the cross terms) in the equivalent covariance matrix is K 2 . The condition that M > K 2 should be satisfied to identify all the signals.
After the eigendecomposition of R q,q+q r , we obtain unitary matrices U q,q+q r ∈ C M ×K 2 and V q,q+q r ∈ C M ×(M −K 2 ) , which are composed of signal and noise subspace eigenvectors, respectively. Similar to the traditional multiple signal classification algorithm (MUSIC), the spatial spectrum for the frequency pairs of f q and f q+q r is estimated as follows: In this spatial spectrum, the negative effects of the grating lobes and sensor phase uncertainties are avoided. However, the cross terms in Case II are the new interferences that degrade the DOA estimation performance of the incident signals. Given the sensor phase term (i.e., e j(ϕ m,k 1 −ϕ m,k 2 ) ) in (7), the vector b r in (9) does not match B q,q+q r in Case II. Consequently, the peaks of the cross terms are weaker than those of the incident signals. Additionally, the cross term θ k 1 ,k 2 , k 1 = k 2 varies with the frequency. These peaks are further weakened by the summation of P q (θ) in the frequency domain: where q L and q H are the lowest and highest frequency bins, respectively. According to Parseval's theorem, the spatial spectrum P(θ) is also the power of the broadband signals.

DD Sensor Phase Calibration
We first apply the CM algorithm [27,28,33] to estimate the steering matrix in each frequency bin. Then, the DD sensor phase responses are extracted from the estimated steering matrix and the permutation problems are finally solved.

Steering Matrix Estimation
We assume that U q ∈ C M ×K is the signal subspace of the covariance matrix R q through eigendecomposition. Thus, we have The norm of A q m,k can be written as follows: where |•| is the modulus of a complex number, u H q(m) is the mth row of U q , and t q(k) is the kth column of T q .
Considering that the modulus of A q m,k is always equal to one for any indices of m and k, we have where , and vec(•) is the matrix vectorization. The M − 1 equations in (12) are written in a matrix form as follows: where 0 M−1 is (M − 1) × 1 vector with all elements equal to zero, and K q is an (M − 1) × K 2 matrix whose (m − 1)th row is vec T (K q,m ), m = 2, · · · , M . The vector t q(k) ⊗ t * q(k) belongs to the null space of K q (denoted as N K q ). The null space N K q can be calculated from the singular value decomposition of matrix K q . If the right singular matrix of K q is Z q = [z 1 , · · · , z K 2 ] and z k is the kth right singular vectors of K q , we find that In accordance with [27,33], t q(k) ⊗ t * q(k) , k = 1, . . . , K is also a group of linearly independent vectors of N K q . Thus, any vector w q ∈ N K q can be expressed as follows: where α q t and β q t are the coefficients under the different basis of null space N K q , vec −1 is the inverse operation of matrix vectorization, and Σ q = diag β q 1 , · · · , β q K . If Σ q is a nonsingular matrix, then we obtain where (15) is similar to the eigendecomposition of the matrix W −1 q W ' q . Additionally, we can calculate the steering matrix A q as follows by considering that A q = U q T q and A q 1,k : , E q is the unitary matrix composed of the eigenvectors of W −1 q W ' q , J ' q = [e n 1 , e n 2 , · · · , e n K ] is a K × K selection matrix, n 1 = n 2 = · · · = n K , n 1 , · · · , n K ∈ {1, 2, · · · , K}, and e n k is a vector with the n k th element equal to one and with the other equal to zero. The arrangement of the columns of the steering matrix at different frequency bins is adjusted by the selection matrix J ' q , which is calculated later.

Sensor Phase Estimation
In the qth frequency bin, we define where angle{•} is the phase of the bracketed variable. The relationship of the phase of the steering matrix can be expressed as follows: where amodb is the remainder of a divided by b and a selection matrix J " q = [e n 1 , e n 2 , · · · , e n K ]; n 1 = n 2 = · · · = n K , n 1 , · · · , n K ∈ {1, 2, · · · , K} is introduced in (18) to align the sequence of the sensor phase responses between ψ q and ψ q . Then, we have The sensor phase is frequency invariant and can be estimated by averaging (19) along the frequency domain:

Permutation Problem Solving
Two permutation problems, namely, J ' q in (16) and J " q in (19), need to be addressed.
To align the columns of the estimated steering matrix in (16) at different frequency bins, we calculate the selection matrix J ' q by minimizing the difference of the steering matrix at neighboring frequency bins as [27] which means that^ A q , q = q L + 1, · · · , q H aligns with^ A q L . In (19), the selection matrix J " q is used to align the sequence of the sensor phase responses between ψ q and ψ q . Given that the sensor phase response is frequency invariant, the variance of ψ q − ψ q J " q mod(2π) along the frequency domain approaches zero for an appropriate J " q . Hence, J " q can be calculated as follows: (22) where var{•} is the variance of the bracketed vector.

Algorithmic Steps
In the proposed method, the DOAs are first estimated without the influences of the unknown sensor phase and spatial aliasing. Then, the DD sensor phase is calibrated with the estimated DOAs. The algorithmic steps for the sparse ULA are described in Table 1. Table 1. Algorithmic steps in the proposed method. Algorithm: Self-calibration for sparse ULAs with unknown DD sensor phase Input: Phase parameters of the standard sensor, observations of the sparse ULA and the individual standard sensor Output: Estimated DOAs and the DD sensor phase responses of the sparse ULA 1: x l ( f q ) is obtained by the Q-point FFTs on the L sections of the array observations. 2: The sample covariance matrixR q is calculated in (4).

3:
Noise power η q is calculated by averaging the small M − K eigenvalues of the sample covariance matrixR q . 4: // Step 1 Estimate the DOAs

6:
The DOA for the frequency pairs of f q 1 and f q 1 +q r is estimated in (9). 7: The broadband spatial spectrum is calculated in (10).

7:
The broadband spatial spectrum is calculated in (10). 8: // Step 2 Estimate the DD sensor phase 9: The steering matrix ˆq A is estimated in (16) and the selection matrix q ′ J is calculated in (21). 10: The sensor phase is estimated in (20) and the selection matrix q ′′ J is calculated in (22).

Discussion on the Deployed Standard Sensors
Although an individual standard sensor is deployed in the previous derivation, we still need to discuss why the individual standard sensor is required in the proposed method.
Lemma. At least one individual standard sensor is required to calibrate the DD sensor phase.
is estimated in (16) and the selection matrix J ' q is calculated in (21). 10: The sensor phase is estimated in (20) and the selection matrix J " q is calculated in (22).

Discussion on the Deployed Standard Sensors
Although an individual standard sensor is deployed in the previous derivation, we still need to discuss why the individual standard sensor is required in the proposed method.

Lemma 1. At least one individual standard sensor is required to calibrate the DD sensor phase.
Proof. After deploying an individual standard sensor along the axis of the original uncalibrated sparse ULAs, a new ULA is formed and the deployed standard sensor becomes the first sensor of the new ULA. In accordance with (16), the estimated steering matrix is normalized by the first row of the matrixÃ q . If the phase responses of the first sensor (i.e., the deployed standard sensor) are symbolized as ϕ ' 1,• = [ϕ 1,1 , ϕ 1,2 , · · · , ϕ 1,K ] T , then the estimated steering matrix can be expressed as follows: whereÃ norm q =Ã q diag −1 e T 1Ã q . If ϕ ' 1,• can be estimated, then the deployed sensor phase parameter is not necessary. In other words, specifically, no deployed sensor is required. However, if ϕ ' 1,• cannot be identified, at least one deployed sensor is needed.
The covariance matrix reconstructed by the steering matrix^ A ' q is represented as In consideration of the commutativity of the Hadamard product, Equation (24) can be reformed as where^ A q is the steering matrix defined in (16).
According to (25), the deployed sensor phase response ϕ ' 1,• is canceled, indicating that the variable ϕ ' 1,• cannot be identified from the covariance matrix. Hence, at least one deployed standard sensor with known phase response is required to estimate the unknown sensor phase responses of the original uncalibrated ULA.

Comparison with the Current Self-Calibration Algorithms
In the current self-calibration algorithms, an iteration-based strategy is commonly used, in which the DOAs and sensor phase are alternatively estimated during the iterations [12][13][14][15][16]. However, the local minima are easily converged due to the poor initialization. The initializations of the DOAs are difficult to choose, especially for the sparse linear arrays, because of the spatial aliasing.
To avoid the disadvantages of the iteration strategy, the self-calibration method in [27] and the proposed method both work without iterations. In [27], the DOAs are first obtained from the estimated steering matrix with two precalibrated sensors and the DD sensor phase is then calculated. The DOA estimation accuracy is significantly influenced by the number of the precalibrated sensors. The estimation of the DD sensor phase is also limited. Moreover, the assumption of the multiple precalibrated sensors that belong to the uncalibrated array is unrealistic in practice.
By contrast, the DOA estimation accuracy of the proposed method is improved owing to the use of the whole aperture of the sensor array and elimination of the unknown sensor phase. Additionally, spatial aliasing is no longer a problem when the proposed method is applied to the sparse ULAs since it is based on the frequency difference methodology. Notably, only an additional standard sensor is exploited for the extraction of the sensor phase, which is easier to implement in real tasks since the sensor is not a part of the original uncalibrated array.

Simulation Results
A sparse ULA composed of 11 acoustic sensors is considered in the simulation. The wave speed in the water medium is 1500 m/s. The interelement spacing is equal to 3.75 m (i.e., the half wavelength of 200 Hz). Two uncorrelated broadband signals impinge on the array from the directions of 50 • and 80 • , respectively, and they cover the frequency band from 200 to 600 Hz. The test array is sparse relative to the frequency band of the incident signals.
A standard sensor is deployed along the axis of the sparse ULA and a new ULA with M = 12 sensors is formed. The deployed sensor now becomes the first sensor of the new ULA. The accurate sensor phase of the deployed sensor and the unknown sensor phase responses, which are frequency invariant and randomly generated from uniform distributions in the range −π ∼ π, are listed in Table 2. The sample frequency is set to 5 kHz.

Example for DOAs and DD Sensor Phase Estimation
The observed data are partitioned into L = 200 sections, with each section containing 500 samples. A Q-point FFT is applied to each section, and Q = 500. A total of 41 frequency bins are located in the frequency band. We set the reference frequency f r = 200 Hz to avoid spatial aliasing. SNR is set to 5 dB. In the proposed method, the first sensor is assumed to be precalibrated.
The spatial spectra of the proposed method with f q = 200, 250, 300, 350, and 400 Hz (i.e., Equation (9)) and the broadband (i.e., Equation (10)) are plotted in Figure 3a. The spectral values of the cross terms are much lower than the peaks of the incident signals due to the mismatch between b r in (9) and B q 1 ,q 1 +q r in (7). Hence, two sharp peaks can be observed in all the spatial spectra in Figure 3a. After summing the spatial spectra along the frequency domain, we observe that the broadband spectrum has relatively flat sidelobe sectors. Additionally, spatial aliasing is effectively avoided.
(i.e., Equation (9)) and the broadband (i.e., Equation (10)) are plotted in Figure 3a. The spectral values of the cross terms are much lower than the peaks of the incident signals due to the mismatch between r b in (9) and 1 1 , r q q q + B in (7). Hence, two sharp peaks can be observed in all the spatial spectra in Figure 3a. After summing the spatial spectra along the frequency domain, we observe that the broadband spectrum has relatively flat sidelobe sectors. Additionally, spatial aliasing is effectively avoided. The unknown DD sensor phase parameters are then obtained by applying the estimated DOAs. The true and the estimated sensor phase values are displayed in Figure 3b. We can observe that the proposed method shows a high estimation accuracy of the DD sensor phase.

Statistical Performance of the Proposed Algorithm
The statistical performance of the proposed algorithm is shown in the following simulation. Considering that studies on the self-calibration of the unknown DD sensor phase are limited, Weiss' method in [25], Yang's method in [27], and the Cramer-Rao bounds (CRBs) [25,27] are applied for comparison. In Weiss' method, the number of precalibrated sensors is set to three. Yang's method requires two precalibrated sensors. In addition, the Weiss' method is a narrowband processing algorithm that is applied in each frequency bin and the sensor phase parameters and the DOAs are averaged over the frequency band. The estimates of DOAs and sensor phase responses are measured by the root mean square error (RMSE). A total of 200 Monte Carlo experiments are conducted.
The statistical performance of the estimates of the DOAs and sensor phase responses versus the SNR are shown in Figure 4a,b, respectively. In Weiss' method, the low estimation accuracy of the DOAs and sensor phase parameters is caused by the poor initialization of the DOAs resulting from the spatial aliasing and the short aperture of the precalibrated subarray. Yang's method can calibrate the sensor phase responses whose The unknown DD sensor phase parameters are then obtained by applying the estimated DOAs. The true and the estimated sensor phase values are displayed in Figure 3b. We can observe that the proposed method shows a high estimation accuracy of the DD sensor phase.

Statistical Performance of the Proposed Algorithm
The statistical performance of the proposed algorithm is shown in the following simulation. Considering that studies on the self-calibration of the unknown DD sensor phase are limited, Weiss' method in [25], Yang's method in [27], and the Cramer-Rao bounds (CRBs) [25,27] are applied for comparison. In Weiss' method, the number of precalibrated sensors is set to three. Yang's method requires two precalibrated sensors. In addition, the Weiss' method is a narrowband processing algorithm that is applied in each frequency bin and the sensor phase parameters and the DOAs are averaged over the frequency band. The estimates of DOAs and sensor phase responses are measured by the root mean square error (RMSE). A total of 200 Monte Carlo experiments are conducted.
The statistical performance of the estimates of the DOAs and sensor phase responses versus the SNR are shown in Figure 4a,b, respectively. In Weiss' method, the low estimation accuracy of the DOAs and sensor phase parameters is caused by the poor initialization of the DOAs resulting from the spatial aliasing and the short aperture of the precalibrated subarray. Yang's method can calibrate the sensor phase responses whose performance is also affected by the directions of the sources. In Yang's method, the DOAs are estimated with the first two precalibrated sensors at an equivalent frequency of ( f H − f L )/2 [27]. At this equivalent frequency, the spatial ambiguity can be avoided, whereas the estimation accuracy of the DOAs is still limited by the aperture of the first two precalibrated sensors. The proposed method has the best DOA estimation performance among the methods mainly because of the applications of the whole aperture of the sensor array and the antialiasing processing. The calibration performance of the sensor phase responses is also improved owing to the unambiguous and accurate DOAs. The CRBs of the DOA and sensor phase estimation in Figure 4a,b are also well approached by the proposed method.
whereas the estimation accuracy of the DOAs is still limited by the aperture of the first two precalibrated sensors. The proposed method has the best DOA estimation performance among the methods mainly because of the applications of the whole aperture of the sensor array and the anti-aliasing processing. The calibration performance of the sensor phase responses is also improved owing to the unambiguous and accurate DOAs. The CRBs of the DOA and sensor phase estimation in Figure 4a,b are also well approached by the proposed method.  We repeat the simulation by fixing SNR = 5 dB and varying the number of sections L from 20 to 300. Each section also contains 500 samples. The estimation performance of DOAs and sensor phase responses versus the number of sections are shown in Figure  5a,b, respectively. Similar to the results in Figure 4, the proposed method has the lowest RMSE values among the compared methods, and the CRB is approached. In the next simulation, the effect of the signal bandwidth on the performance of the mentioned self-calibration algorithms is discussed. We fix the lower limit frequency of the incident signals L f as 200 Hz, whose wavelength is equal to twice the interelement spacing. Additionally, the higher limit frequency H f ranges from 400 to 2400 Hz at a step of 200 Hz, indicating that the bandwidth of the incident signals increases from 1 oct We repeat the simulation by fixing SNR = 5 dB and varying the number of sections L from 20 to 300. Each section also contains 500 samples. The estimation performance of DOAs and sensor phase responses versus the number of sections are shown in Figure 5a,b, respectively. Similar to the results in Figure 4, the proposed method has the lowest RMSE values among the compared methods, and the CRB is approached.
whereas the estimation accuracy of the DOAs is still limited by the aperture of the first two precalibrated sensors. The proposed method has the best DOA estimation performance among the methods mainly because of the applications of the whole aperture of the sensor array and the anti-aliasing processing. The calibration performance of the sensor phase responses is also improved owing to the unambiguous and accurate DOAs. The CRBs of the DOA and sensor phase estimation in Figure 4a,b are also well approached by the proposed method.  We repeat the simulation by fixing SNR = 5 dB and varying the number of sections L from 20 to 300. Each section also contains 500 samples. The estimation performance of DOAs and sensor phase responses versus the number of sections are shown in Figure  5a,b, respectively. Similar to the results in Figure 4, the proposed method has the lowest RMSE values among the compared methods, and the CRB is approached. In the next simulation, the effect of the signal bandwidth on the performance of the mentioned self-calibration algorithms is discussed. We fix the lower limit frequency of the incident signals L f as 200 Hz, whose wavelength is equal to twice the interelement spacing. Additionally, the higher limit frequency H f ranges from 400 to 2400 Hz at a step of 200 Hz, indicating that the bandwidth of the incident signals increases from 1 oct In the next simulation, the effect of the signal bandwidth on the performance of the mentioned self-calibration algorithms is discussed. We fix the lower limit frequency of the incident signals f L as 200 Hz, whose wavelength is equal to twice the interelement spacing. Additionally, the higher limit frequency f H ranges from 400 to 2400 Hz at a step of 200 Hz, indicating that the bandwidth of the incident signals increases from 1 oct to approximately 3.5 oct. The larger the f H , the sparser of the sensor array. Additionally, SNR is set to 5 dB, and the number of sections L = 200.
The RMSEs of the DOAs and sensor phase responses versus the signal bandwidth are plotted in Figure 6a,b, respectively. Weiss' method fails to estimate the unknown parameters. In Yang's method, the equivalent frequency of ( f H − f L )/2 does not satisfy the spatial Nyquist sampling theorem when the higher limit frequency f H is larger than 600 Hz. Accordingly, the estimation performance is dramatically degraded due to the spatial ambiguity for f H > 600 Hz. However, the directions of the sources are estimated at the reference frequency, and the spatial Nyquist sampling theorem is always satisfied in the proposed method. Moreover, the DOA estimation performance can be improved with more frequency bins applied due to the increased signal bandwidth. The RMSE of the proposed method is also approached by the CRB. 600 Hz. Accordingly, the estimation performance is dramatically degraded due to the spatial ambiguity for H >600Hz f . However, the directions of the sources are estimated at the reference frequency, and the spatial Nyquist sampling theorem is always satisfied in the proposed method. Moreover, the DOA estimation performance can be improved with more frequency bins applied due to the increased signal bandwidth. The RMSE of the proposed method is also approached by the CRB.

Conclusions
A self-calibration method is proposed by deploying an individual standard sensor with a known sensor phase response along the axis of the uncalibrated sparse linear array, in which the unknown DD sensor phase and the aliasing-free DOAs are estimated. The DOA is first accurately estimated owing to the use of the whole aperture of the sensor array and the elimination of unknown sensor phase. Additionally, spatial aliasing is also solved based on the frequency difference methodology. Then, the steering matrix is estimated by the CM algorithm with the known phase parameters of the deployed standard sensor. The DD sensor phase is finally extracted from the steering matrix by using the estimated DOAs.
The advantages of the proposed method are summarized as follows: (1) The DOAs can be accurately estimated even if the sensor phase is unknown. (2) The spatial aliasing of the DOA estimation is avoided. (3) Only an additional standard sensor is exploited for the extraction of the DD sensor phase, which is easier to implement in real tasks.

Conclusions
A self-calibration method is proposed by deploying an individual standard sensor with a known sensor phase response along the axis of the uncalibrated sparse linear array, in which the unknown DD sensor phase and the aliasing-free DOAs are estimated. The DOA is first accurately estimated owing to the use of the whole aperture of the sensor array and the elimination of unknown sensor phase. Additionally, spatial aliasing is also solved based on the frequency difference methodology. Then, the steering matrix is estimated by the CM algorithm with the known phase parameters of the deployed standard sensor. The DD sensor phase is finally extracted from the steering matrix by using the estimated DOAs.
The advantages of the proposed method are summarized as follows: (1) The DOAs can be accurately estimated even if the sensor phase is unknown. (2) The spatial aliasing of the DOA estimation is avoided. (3) Only an additional standard sensor is exploited for the extraction of the DD sensor phase, which is easier to implement in real tasks. (4) The CRBs of the DOA and DD sensor phase estimation are well approached.