Next Article in Journal
Three-Dimensional Blood Vessel Model with Temperature-Indicating Function for Evaluation of Thermal Damage during Surgery
Previous Article in Journal
Dual-Mode Gas Sensor Composed of a Silicon Nanoribbon Field Effect Transistor and a Bulk Acoustic Wave Resonator: A Case Study in Freons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Noncircular MUSIC Algorithm Based on the Concept of the Difference and Sum Coarray

School of Information and Electronics, Beijing Institute of Technology, 5 South Zhongguancun Street, Haidian District, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(2), 344; https://doi.org/10.3390/s18020344
Submission received: 15 December 2017 / Revised: 17 January 2018 / Accepted: 22 January 2018 / Published: 25 January 2018
(This article belongs to the Section Physical Sensors)

Abstract

:
In this paper, we propose a vectorized noncircular MUSIC (VNCM) algorithm based on the concept of the coarray, which can construct the difference and sum (diff–sum) coarray, for direction finding of the noncircular (NC) quasi-stationary sources. Utilizing both the NC property and the concept of the Khatri–Rao product, the proposed method can be applied to not only the ULA but also sparse arrays. In addition, we utilize the quasi-stationary characteristic instead of the spatial smoothing method to solve the coherent issue generated by the Khatri–Rao product operation so that the available degree of freedom (DOF) of the constructed virtual array will not be reduced by half. Compared with the traditional NC virtual array obtained in the NC MUSIC method, the diff–sum coarray achieves a higher number of DOFs as it comprises both the difference set and the sum set. Due to the complementarity between the difference set and the sum set for the coprime array, we choose the coprime array with multiperiod subarrays (CAMpS) as the array model and summarize the properties of the corresponding diff–sum coarray. Furthermore, we develop a diff–sum coprime array with multiperiod subarrays (DsCAMpS) whose diff–sum coarray has a higher DOF. Simulation results validate the effectiveness of the proposed method and the high DOF of the diff–sum coarray.

1. Introduction

Noncircular (NC) signals, such as amplitude-modulated (AM) signals and binary phase-shift keying (BPSK)-modulated signals, have been widely applied in various communication systems [1,2,3,4,5,6,7,8]. Different from circular signals, which can only use the information in the covariance matrix for direction finding, NC signals can use the information in both the covariance matrix and the elliptic covariance matrix for direction finding. This NC property can be utilized to increase the degrees of freedom (DOFs) and improve the estimation performance. A lot of DOA estimation algorithms for NC sources have been developed, such as the NC MUSIC method [9], NC Root-MUSIC method [1], NC ESPRIT method [10] and NC Unitary ESPRIT method [11]. These traditional NC DOA estimation algorithms utilize the complex conjugate counterpart of the received signals to obtain the NC covariance matrix, which corresponds to a virtual array consisting of the physical array and its flipped array [12]. These algorithms mostly utilize the uniform linear array (ULA) as the array model and can detect at most 2 ( N 1 ) sources with N physical sensors. In order to detect more sources, some NC high-order cumulant MUSIC methods based on the non-Gaussian characteristic of many NC sources, such as the NC 2q-MUSIC method [13] and ROOT NC 4-MUSIC method [2], have been proposed. However, the array model in these methods is still the ULA, and the computation complexity of the cumulant-based methods are large. As we know, the array aperture is a fundamental parameter that affects the performance of the DOA estimation. For the ULA, its aperture is usually less than that of the sparse array constructed based on the concept of the coarray [14,15,16,17,18,19]. Thus, an estimator, which jointly utilizes the NC property, the sparse array and the concept of the coarray, is likely to detect more sources than the traditional NC DOA estimators.
In recent years, sparse arrays have been attracting more and more attention due to the high number of DOFs [20,21,22,23]. The concept of the coarray based on the Khatri–Rao (KR) product [21,24] is applied to design various of sparse arrays, such as the coprime array [22]. The coprime array consists of two uniform linear subarrays with M and N sensors, respectively, where M and N are a coprime pair of positive integers. The M element subarray with spacing N units and the N element subarray with spacing M units share the first sensor. By applying the KR product to the covariance matrix of the received signals, the difference coarray [21,22,23,24,25,26,27] of the coprime array is constructed. With only M + N 1 physical sensors, the virtual array can detect O ( M N ) sources. In order to achieve a higher number of DOFs, many novel methods and optimized configurations based on the difference coarray for the coprime array have been proposed, such as the sparsity enforced recovery technique for the coprime array [18], the coarray interpolation method for the coprime array [28], the extended coprime array [29], the generalized coprime arrays [24] and the coprime array with multiperiod subarrays (CAMpS) [30,31]. Since these methods and improved configurations are developed based on the properties of the difference coarray for the coprime array, their DOFs cannot be more than twice the physical aperture. Constructing a novel virtual array with larger array aperture than the difference coarray is another useful way to increase the DOFs. In [16], we utilized both the temporal information and the spatial information of the received signals to propose the vectorized conjugate augmented MUSIC (VCAM) algorithm, which can construct the difference and sum (diff–sum) coarray. The diff–sum coarray comprises both the difference set and the sum set so that it can achieve a higher number of DOFs than the above difference coarrays. In addition, the aperture of the diff–sum coarray can be more than twice the physical aperture, which could help to decrease the array size. However, the source signals applied in both the VCAM algorithm and the methods constructing the difference coarray are circular signals.
In this paper, we propose an improved NC MUSIC algorithm based on the concept of the coarray to perform the DOA estimation of NC quasi-stationary sources. We name the novel method as the vectorized NC MUSIC (VNCM) method. By applying the KR product operation to the NC covariance matrix, we can obtain an equivalent received signal, which seems to be received from a diff–sum coarray. The diff–sum coarray, which is symmetrical with the zero point as the center, can be divided into three parts: the difference coarray, the nonpositive sum coarray and the non-negative sum coarray. Thus, the diff–sum coarray is likely to achieve a higher number of DOFs than the traditional NC virtual array consisting of the physical array and its flipped array. Also, we utilize the characteristic of the quasi-stationary sources instead of the spatial smoothing method to solve the single snapshot issue (the coherent issue) of the equivalent received signal so that the available DOFs will not be reduced by half. Due to the complementarity between the difference set and the sum set for the coprime array [16], we utilize the CAMpS, which contains the prototype coprime array as a special case, as the basic array model in this paper. In particular, we summarize and prove the properties and DOFs of the diff–sum coarray for the CAMpS. Furthermore, we improve the CAMpS to propose a diff–sum coprime array with multiperiod subarrays (DsCAMpS) which can achieve a higher number of DOFs than the CAMpS. Extensive simulations are provided to verify the performance of the diff–sum coarray and the effectiveness of the VNCM method.
Notations: Throughout the paper, we utilize lowercase bold italic letters to denote vectors, for example, a . We utilize capital bold italic letters to denote matrices, for example, A . ( . ) T , ( . ) * and ( . ) H respectively represent transpose, conjugation and conjugate transpose. E is used to denote the expectation operation and v e c ( . ) represents the vectorizing operation. ⊗ and ⊙ respectively denote the left Kronecker product and the Khatri–Rao product.

2. Data Model

Consider a linear antenna array consisting of R elements. The set of the sensor locations is denoted as d = { d 1 , , d R } . The first sensor is selected as the reference, that is, d 1 = 0 . Q far-field, narrowband, uncorrelated and NC sources impinge on this array from directions { θ 1 , , θ Q } . The NC sources in the paper are assumed as wide-sense quasi-stationary [32,33] with the NC rate as ρ = 1 [13,34]. Then, the received signal can be modeled as
x ( t ) = q = 1 Q a ( θ q ) s q ( t ) + n ( t ) = As ( t ) + n ( t ) ,
where s ( t ) = [ s 1 ( t ) , s 2 ( t ) , , s Q ( t ) ] T is the observed signal vector, n ( t ) = [ n 1 ( t ) , n 2 ( t ) , , n R ( t ) ] T represents the white circular complex Gaussian noise vector with zero mean and variance σ n 2 , and A = [ a ( θ 1 ) , a ( θ 2 ) , , a ( θ Q ) ] is the steering matrix. The qth ( q = 1 , 2 , , Q ) column vector of A is a ( θ q ) = [ 1 , e j 2 π d 2 sin ( θ q ) / λ , , e j 2 π d R sin ( θ q ) / λ ] T with λ being the signal wavelength.
According to the noncircularity of the signals, the following relationships exist:
s q t = e j φ q / 2 s q t , q = 1 , 2 , , Q
and
E [ s q 2 ( t ) ] = e j φ q E [ s q ( t ) s q * ( t ) ] ,
where φ q is the NC phase of s q ( t ) , s q t is real-valued with zero-phase, E [ s q ( t ) s q * ( t ) ] is the covariance and E [ s q 2 ( t ) ] is the eliptic covariance. As the signals are quasi-stationary, the following assumptions [33] hold as well.
(A1) Each signal s q ( t ) is wide-sense quasi-stationary with the frame length being L and the total number of frames being L ˜ , that is,
E s q i 2 t = ( i 1 ) L i L 1 s q ( t ) s q * ( t ) / L σ q 2 [ i ] , i = 1 , 2 , , L ˜ ,
where s q i = s q ( t ) , t = ( i 1 ) L , ( i 1 ) L + 1 , i L 1 represents the snapshots of the qth signal in the ith frame.
(A2) The variance sequence σ q 2 [ i ] , i = 1 , , L ˜ is wide-sense stationary and uncorrelated with each other. That is to say, E σ q 2 [ i ] = m ¯ q , E σ q 2 [ i ] m ¯ q 2 = σ ˜ q 2 and E σ q 1 2 [ i ] m ¯ q 1 σ q 2 2 [ i ] m ¯ q 2 = 0 , q 1 q 2 .
For the ith ( i = 1 , , L ˜ ) frame, combining the received signal and its conjugate version together results in the following NC observation vector,
r t = x ( t ) x * ( t ) = A Γ A Γ * s ( t ) + n ( t ) n * ( t ) = A s t + n ( t ) ,
where Γ = d i a g ( e j φ 1 / 2 , e j φ 2 / 2 , , e j φ Q / 2 ) , s ( t ) = s 1 ( t ) , s 2 ( t ) , , s Q ( t ) T , A = A Γ T , A Γ H T and n ( t ) = n T ( t ) , n H ( t ) T . The qth column vector of A is a θ q = e j φ q / 2 a T θ q , e j φ q / 2 a H θ q T , which can be considered to correspond to a traditional NC virtual array consisting of 2 R 1 different virtual sensors with the zero point as the center. From (5), the covariance matrix of r ( t ) can be obtained as
R r r i = E r ( t ) r H ( t ) = A N C R s s A N C H + σ n 2 I ,
where R s s = d i a g σ 1 2 i , σ 2 2 i , , σ Q 2 i and A N C = A T , A ϕ H T = a N C θ 1 , a N C θ 2 , , a N C θ Q with a N C θ q = a T θ q , e j φ q a H θ q T and ϕ = d i a g ( e j φ 1 , e j φ 2 , , e j φ Q ) . The traditional NC DOA estimators, of which the array model is the ULA, directly perform eigenvalue decomposition on R r r i to obtain the noise subspace U n . According to [1,9], the maximum number of detectable signals (MNDS) is determined by the degree of det { a b l k H U n U n H a b l k } with a b l k = b l k d i a g ( a ( θ q ) , a * ( θ q ) ) . Therefore, at most 2 R 2 sources can be detected. Since the DOF of the traditional NC virtual array with R sensors is 2 R 1 , we have MNDS = DOF−1.

3. Vectorized NC MUSIC Algorithm

In this paper, we focus on combining the NC property and the concept of the coarray. NC DOA estimators based on the concept of the coarray have two problems to be resolved: (1) how to solve the single snapshot issue (the coherent issue) of the equivalent received signal obtained by the KR product operation; and (2) how to generate a coarray with a large array aperture. In order to resolve the two problems, we propose a vectorized NC MUSIC (VNCM) algorithm.
Within the ith ( i = 1 , , L ˜ ) frame, vectorizing R r r i
z i = v e c R rr i = A N C * A N C p i + σ n 2 v ,
where p i = σ 1 2 i , σ 2 2 i , , σ Q 2 i T , v = v e c I and the qth ( q = 1 , 2 , , Q ) column vector of A N C * A N C has the form
a N C * θ q a N C θ q = a θ q e j φ q a * θ q * a θ q e j φ q a * θ q   = a * θ q a θ q a * θ q e j φ q a * θ q e j φ q a θ q a θ q e j φ q a θ q e j φ q a * θ q = N C d i f f 1 N C s u m 1 N C s u m 2 N C d i f f 2 .
Having the similar form as the received signal x ( t ) in (1), z [ i ] in (7) can be seen as the equivalent received signal at a virtual array whose steering matrix is given by A N C * A N C . According to (8), we can find that the virtual sensor locations can be represented as d r 1 d r 2 ± d r 1 + d r 2 , 1 r 1 , r 2 R . Thus, the novel NC virtual array is a diff–sum coarray consisting of both the difference and sum results. Specially, in (8), N C d i f f 1 and N C d i f f 2 correspond to the difference coarray, N C s u m 1 corresponds to the nonpositive sum coarray and N C s u m 2 corresponds to the non-negative sum coarray. In the following, we would show that the MNDS of the VNCM method is determined by the DOFs of the diff–sum coarray, regardless of the NC phase φ q . Assume that the consecutive range of the diff–sum coarray is [ l c , l c ] . As the performance of the MUSIC-class methods are mainly determined by the ULA part of the virtual array, we remove the repeated and discrete lags in (7) to obtain
z i = B p i + σ n 2 v ,
where v is a ( 2 l c + 1 ) × 1 vector extracted from v and B = [ b ( θ 1 ) , b ( θ 2 ) , , b ( θ Q ) ] is a ( 2 l c + 1 ) × Q matrix with the qth column being
b θ q = e j φ q a ˜ 1 θ q a ˜ 2 θ q e j φ q a ˜ 3 θ q = A ˜ ( θ q ) φ ˜ ( θ q ) ,
where a ˜ 1 θ q , a ˜ 2 θ q and a ˜ 3 θ q are respectively extracted from a * ( θ q ) a * ( θ q ) , a * ( θ q ) a ( θ q ) and a ( θ q ) a ( θ q ) , A ˜ ( θ q ) = b l k d i a g a ˜ 1 θ q , a ˜ 2 θ q , a ˜ 3 θ q is a block diagonal matrix and φ ˜ ( θ q ) = e j φ q , 1 , e j φ q T .
Combining all the NC virtual received signals z i , i = 1 , 2 , , L ˜ , we can obtain the NC virtual frame-data matrix over all the L ˜ frames as
Z = z 1 , z 2 , , z L ˜ = B P + σ n 2 V ,
where P = p 1 , p 2 , , p L ˜ and V = v , v , , v ( 2 l c + 1 ) × L ˜ . According to A2) in Section 2, it is obvious that each row in P is a wide-sense stationary process with the expectation being m ¯ q . Thus, the expectation vector of z i can be expressed as
z ¯ = E z i 1 L ˜ i = 1 L ˜ z i = B 1 L ˜ i = 1 L ˜ p i + σ n 2 v = B p ¯ + σ n 2 v ,
where p ¯ = m ¯ 1 , m ¯ 2 , , m ¯ Q T . Subtracting z ¯ from each column vector of the NC virtual frame-data matrix Z , we have
Z ˜ = Z z ¯ 1 , 1 , , 1 1 × L ˜   = B P p ¯ 1 , 1 , , 1 1 × L ˜ + σ n 2 V σ n 2 V   = B P ˜ ,
where Z ˜ = z ˜ 1 , z ˜ 2 , , z ˜ L ˜ with z ˜ i = z i z ¯ , i = 1 , 2 , , L ˜ , and P ˜ = p ˜ 1 , p ˜ 2 , , p ˜ L ˜ with p ˜ i = p i p ¯ . It is obvious that each row in P ˜ is a zero-mean wide-sense stationary process. Besides, according to (A2) in Section 2, each row sequence in P ˜ is uncorrelated with the other row sequences. Thus, similar to (6), we can obtain the correlation matrix of z ˜ i , i = 1 , 2 , , L ˜ as
R z ˜ z ˜ = E z ˜ i z ˜ H i 1 L ˜ Z ˜ Z ˜ H = B R p ˜ p ˜ B H ,
where R p ˜ p ˜ = d i a g σ ˜ 1 2 , σ ˜ 2 2 , , σ ˜ Q 2 . As R p ˜ p ˜ is a full-rank matrix, the eigenvalue decomposition of R z ˜ z ˜ can be obtained as
R z ˜ z ˜ = U p ˜ Σ p ˜ U p ˜ H + U n ˜ Σ n ˜ U n ˜ H ,
where U p ˜ is the signal subspace whose columns represent the signal subspace eigenvectors of R z ˜ z ˜ , U n ˜ is the noise subspace whose columns represent the noise subspace eigenvectors of R z ˜ z ˜ , Σ p ˜ = d i a g λ 1 , λ 2 , , λ Q with λ 1 , λ 2 , , λ Q representing the Q largest eigenvalues of R z ˜ z ˜ , and Σ n ˜ = 0 2 l c + 1 Q × 2 l c + 1 Q . Due to the orthogonality between the signal subspace and the noise subspace, any direction θ q , q = 1 , 2 , , Q from θ 1 , , θ Q satisfies the following equation
b H θ q U n ˜ U n ˜ H b θ q = 0 .
Associated with (10), (16) can be rewritten as
φ ˜ H ( θ q ) A ˜ H ( θ q ) U n ˜ U n ˜ H A ˜ ( θ q ) φ ˜ ( θ q ) = 0 .
According to [1,35,36], A ˜ H ( θ ) U n ˜ U n ˜ H A ˜ ( θ ) is rank deficient at θ = θ q , q = 1 , 2 , , Q . Thus, the DOAs can be estimated by the following estimator,
f θ = 1 d e t A ˜ H ( θ ) U n ˜ U n ˜ H A ˜ ( θ ) .
Searching the direction θ over the range π π 2 2 , π π 2 2 , the DOAs can be obtained from the peaks in f θ .
Remarks. 
When any element e j 2 π l sin θ / λ ( l c l l c ) in A ˜ ( θ ) is replaced by z l ( z = e j 2 π sin θ / λ ), d e t A ˜ H ( θ ) U n ˜ U n ˜ H A ˜ ( θ ) can be seen as a polynomial of degree 4 ( l c 1 ) , whose roots appear in reciprocal conjugate pairs. Thus, the VNCM algorithm can detect up to 2 ( l c 1 ) signals. Since the DOF of the diff–sum coarray is 2 l c + 1 , it can be concluded that M N D S = D O F 3 for the VNCM method. From the discussions above, it is obvious that the MNDS of the VNCM method is determined by the DOF of the diff–sum coarray, regardless of the NC phase. With R physical sensors, the MNDS of the traditional NC DOA estimators is 2 ( R 1 ) . As l c is obtained by combining the difference and sum results of the physical sensor locations, we can conclude l c > R . Therefore, the VNCM algorithm can detect more signals than the traditional NC DOA estimators. Furthermore, designing a sparse array, of which the diff–sum coarray achieves a high number of DOFs, can help further improve the performance of the proposed method. It is noted that in theory, the high-order cumulant-based MUSIC method can also be applied to the sparse array to detect NC sources. The virtual array generated in this kind of method should be able to achieve the same number of DOFs as that generated in the VNCM method. However, this kind of method should first solve one issue, which is how to separate the NC phases from the cumulant matrix to perform the eigenvalue decomposition. The method proposed in this paper can help to solve this issue, which would be a future work for us to do. The difference between the NC high-order cumulant-based method and the VNCM method is that the NC high-order cumulant-based method is restricted to non-Gaussian signal sources, but the signal model in the VNCM method is not necessarily non-Gaussian.

4. The Diff–Sum Coprime Array with Multiperiod Subarrays Based on the Concept of the Diff–Sum Coarray

Now, since the difference set and the sum set for the coprime array are complementary [16,37], we choose the CAMpS [30,31] as the basic array model. In this section, we would summarize and prove the properties and the DOF of the diff–sum coarray for the CAMpS. Then, based on these properties, we improve the CAMpS to propose a diff–sum coprime array with multiperiod subarrays (DsCAMpS) of which the diff–sum coarray achieves a higher number of DOFs.

4.1. The CAMpS and the Concept of the Diff–Sum Coarray

Firstly, we have a quick review of the CAMpS. As shown in Figure 1, the CAMpS, which is the multiperiod extension of the prototype coprime array, consists of two uniform linear subarrays. Subarray 1 contains P 1 M sensors with the intersensor spacing of N units, and Subarray 2 contains P 2 N sensors with the intersensor spacing of M units. We use d to denote the unit interelement spacing. Then, the sensors of the CAMpS are located at S C A M p S = S 1 S 2 , where S 1 = m N d 0 m P 1 M 1 and S 2 = n M d 0 n P 2 N 1 . Due to M and N being coprime, there are min { P 1 , P 2 } common elements between the two subarrays. Hence, the number of the elements in the CAMpS is P 1 M + P 2 N min { P 1 , P 2 } . For convenience, in the following sections of this paper, we normalize all the locations by the unit interelement spacing d.
Applying the VNCM method to the CAMpS, the resulting virtual array can be represented as
S d s = S d i f f S s u m + S s u m ,
where S d s is the diff–sum set, S d i f f = d r 1 d r 2 1 r 1 , r 2 R is the difference set, S s u m + = d r 1 + d r 2 1 r 1 , r 2 R is the non-negative sum set and S s u m = ( d r 1 + d r 2 ) 1 r 1 , r 2 R is the nonpositive sum set.
When choosing the CAMpS in Figure 1 as the basic array model, S d i f f can be expressed as the following union set,
S d i f f = S c d S c d S s d S s d ,
where S c d = d S 2 d S 1 d S 2 S 2 , d S 1 S 1 is the cross-difference set between S 2 and S 1 , S c d = d S 1 d S 2 is the mirrored set of S c d , S s d = d S 1 d S 2 is the self-difference set of S 1 and S 2 and S s d = d S 1 d S 2 is the mirrored set of S s d . Similarly, the total sum set S s u m = S s u m + S s u m can also be expressed as a union set
S s u m = S c s + S s s + S s u m + S c s S s s S s u m ,
where S s s + = d S 1 + d S 1 d S 1 , d S 1 S 1 d S 2 + d S 2 d S 2 , d S 2 S 2 is the self-sum set of S 1 and S 2 , S s s = d S 1 + d S 1 d S 2 + d S 2 is the mirrored set of S s s + , S c s + = d S 1 + d S 2 is the cross-sum set between S 1 and S 2 , S c s = d S 1 + d S 2 is the mirrored set of S c s + .

4.2. The Properties of the Diff–Sum Set for the CAMpS

Without loss of generality, we assume the period P 1 of Subarray 1 is no larger than the period P 2 of Subarray 2, that is, P 2 P 1 1 . Denote S 1 i = N m i 1 M m i M 1 and S 2 j = M n j 1 N n j N 1 as the location set of Period i ( 1 i P 1 ) in Subarray 1 and the location set of Period j ( 1 j P 2 ) in Subarray 2.
According to [30,31], some properties of the difference set for the CAMpS have been summarised as follows:
(1)
S c d contains all the consecutive lags in the range P 1 1 M N M + 1 , P 2 1 M N + N 1 .
(2)
When P 2 > P 1 , the difference set S d i f f contains all the consecutive lags in the range P 2 1 M N N + 1 , P 2 1 M N + N 1 .
(3)
When P 2 = P 1 , the result becomes P 2 1 M N M N + 1 , P 2 1 M N + M + N 1 .
Combining these properties, we can find that when M and N are fixed, the maximum value s c d m a x in the consecutive range of S c d is determined by the period P 2 of Subarray 2, and the minimum value is determined by the period P 1 of Subarray 1. Considering the consecutive range of S d i f f , when P 2 > P 1 , both the minimum and the maximum have the same absolute value as s c d m a x and the two values are only related to P 2 . When P 2 = P 1 , the two values have greater absolute value than s c d m a x since some holes in S d i f f with P 1 = P 2 1 can be filled by the cross-difference results between S 1 P 2 and S 2 . Thus, we can obtain the following two conclusions: (1) When P 2 > P 1 1 , the difference set for the CAMpS has the same consecutive range as that for a simplified CAMpS of which the two subarrays are respectively S 11 and S 2 . (2) When P 2 = P 1 1 , the equivalent array is a simplified CAMpS of which the two subarrays are respectively S 11 S 1 P 1 and S 2 .
In Figure 2, we assume ( M , N ) = ( 3 , 4 ) and show two examples of the difference sets for the CAMpS with ( P 1 , P 2 ) = ( 2 , 3 ) and ( P 1 , P 2 ) = ( 3 , 3 ) as the illustrative examples of the above properties. Figure 2 shows that the difference coarray is symmetrical with the zero point as the center. When P 2 > P 1 , the number of the consecutive elements in the difference coarray is 2 s c d m a x + 1 . When P 1 increases to P 2 , the result becomes 2 s c d m a x + 2 M + 1 .
The following proposition reveals the properties of the non-negative sum set for the CAMpS. The properties of the nonpositive sum set can be deduced by reversing the results of the following proposition.
Proposition 1.
The following facts hold for S c s + and S s u m + :
(a) 
The consecutive range of S c s + is M 1 N 1 , P 1 + P 2 1 M N 1 .
(b) 
When P 2 > P 1 , the non-negative sum set S s u m + contains all the consecutive lags in the range M 1 N 1 , P 1 + P 2 1 M N + N 1 .
(c) 
When P 2 = P 1 : S s u m + with M = 2 contains all the consecutive lags in the range M 1 N 1 , P 1 + P 2 1 M N + N 1 ; when N = 2 , the result becomes M 1 N 1 , P 1 + P 2 1 M N + M 1 ; when M , N > 2 , the result is M 1 N 1 , P 1 + P 2 1 M N + M + N 1 .
(d) 
The non-negative sum set for the CAMpS has the same consecutive range as that for a simplified CAMpS whose two subarrays are respectively S 11 S 1 P 1 and S 2 .
Proof. 
See Appendix A. ☐
In the proof of Proposition 1, we have proved that some holes in the cross-sum set can be aligned with the elements in the self-sum set. Thus, S s u m + has a wider consecutive range than S c s + . In addition, when M and N are fixed, the maximums in the consecutive range of S c s + and S s u m + are determined by the periods of the two subarrays and the minimum is fixed. In Figure 3, we assume ( M , N ) = ( 3 , 4 ) and show two examples of the non-negative sum coarray for the CAMpS with ( P 1 , P 2 ) = ( 2 , 3 ) and ( P 1 , P 2 ) = ( 3 , 3 ) as the illustrative examples of the above properties. It is obvious that when the period of one subarray adds 1, the maximums in the consecutive range of S c s + and S s u m + would add O ( M N ) . Since S s u m is the mirrored set of S s u m + , the consecutive range of the total sum set S s u m = S s u m + S s u m is symmetrical with the zero point as the center. Then, combining Proposition 1 and the properties of the difference set, we give the properties of the diff–sum set S d s for the CAMpS in the following proposition. Defining s d s max 1 = P 1 + P 2 1 M N + N 1 , s d s max 2 = P 1 + P 2 1 M N + M 1 and s d s max 3 = P 1 + P 2 1 M N + M + N 1 , we conclude:
Proposition 2.
When P 2 > P 1 , the consecutive range of S d s for the CAMpS is s d s max 1 , s d s max 1 . When P 2 = P 1 , the consecutive range of S d s is summarised based on the following three cases: (1) when M = 2 , the consecutive range is s d s max 1 , s d s max 1 ; (2) when N = 2 , the consecutive range is s d s max 2 , s d s max 2 ; (3) when M , N > 2 , the consecutive range is s d s max 3 , s d s max 3 .
According to Proposition 2, in Table 1, we summarize the DOF (the number of the consecutive lags) of S d s . The corresponding MNDS of the VNCM method can also be obtained by the equation MNDS = DOF−3.

4.3. The Diff–Sum Coprime Array with Multiperiod Subarrays

According to the properties of the difference set, Propositions 1 and 2, the consecutive range of S d s for the CAMpS is the same as that of S d s for a simplified CAMpS, of which the location sets of the two subarrays can be expressed as S 11 S 1 P 1 and S 2 . We name the equivalent array as the diff–sum coprime array with multiperiod subarrays (DsCAMpS). Since the DOF of the DsCAMpS with P 1 = P 2 is greater than that with P 1 < P 2 , we define P 1 = P 2 in the DsCAMpS. Then, the structure of the DsCAMpS can be shown in Figure 4.
When P 2 2 , Subarray 1 in the DsCAMpS contains 2 periods, of which each period contains M sensors with the intersensor spacing of N units. The displacement between the two periods is ( P 2 2 ) M N + N . Subarray 2 in the DsCAMpS contains P 2 N sensors with the intersensor spacing of M units. The location set of the DsCAMpS can be expressed as
S D s C A M p S = m N 0 m M 1 , ( P 2 1 ) M m P 2 M 1   n M 0 n P 2 N 1 .
Between the two subarrays, there are two common elements which locate at 0 and ( P 2 1 ) M N . Thus, the number of the elements in the DsCAMpS is 2 M + P 2 N 2 . When P 2 = 1 , the DsCAMpS becomes the prototype coprime array with M + N 1 sensors, which means the prototype coprime array is a special kind of the DsCAMpS. Defining s d s max a = 2 P 2 1 M N + N 1 , s d s max b = 2 P 2 1 M N + M 1 and s d s max c = 2 P 2 1 M N + M + N 1 , we can show the DOF of the diff–sum coarray for the DsCAMpS in Table 2. Compared with the CAMpS, the DsCAMpS can achieve a higher number of DOFs when P 2 > 2 . When P 1 = P 2 2 , the two arrays have the same structure and DOF.
In Figure 5, we depict the consisting sets of the diff–sum coarray for the DsCAMpS with M , N , P 2 = 4 , 3 , 3 . Figure 5a shows that the consecutive range of the difference set is [ 30 , 30 ] , which is symmetrical with the zero point as the center. Figure 5b depicts the total sum set S s u m = S s u m + S s u m , which contains all the consecutive lags in the range 66 , 6 6 , 66 . The consecutive range of the difference set for the DsCAMpS with P 2 > 2 overlaps with part of the consecutive range of the total sum set. Thus, as shown in Figure 5c, the diff–sum set, which is the union set of S d i f f and S s u m , contains all the consecutive lags in the range [ 66 , 66 ] . Figure 5 verifies the properties of the diff–sum coarray for the DsCAMpS and the complementarity between the corresponding difference set and sum set.

5. Simulation Results

In this section, we consider the sensor number of all configurations as R = 15 . The unit interelement spacing is d = λ / 2 . Since the DsCAMpS and the CAMpS are the same when P 1 = P 2 2 , we consider the period P 2 in Subarray 2 of both the CAMpS and the DsCAMpS satisfying P 2 > 2 . The configurations utilized in this section are respectively the CAMpS with M , N , P 1 , P 2 = 4 , 3 , 2 , 3 , the DsCAMpS with M , N , P 2 = 4 , 3 , 3 and the ULA = { 0 , 1 , 14 } .

5.1. DOF Comparison

Figure 6 depicts four virtual configurations, which are respectively the traditional NC virtual array for the ULA, the diff–sum coarrays for the ULA, the CAMpS and the DsCAMpS. The traditional NC virtual array is constructed in the NC MUSIC algorithm and the three diff–sum coarrays are obtained by using the VNCM algorithm. As shown in Figure 6a, the traditional NC virtual array for the ULA, which consists of the physical array and its flipped array, contains all the consecutive lags in the range [ 14 , 14 ] . In Figure 6b–d, the consecutive ranges of the three diff–sum coarrays for the ULA, the CAMpS and the DsCAMpS are respectively [ 28 , 28 ] , [ 50 , 50 ] and [ 66 , 66 ] . Comparing Figure 6a,b, it is obvious that the diff–sum coarray achieves a higher number of DOFs than the traditional NC virtual array. From Figure 6b–d, we can find that the diff–sum coarray constructed by using the sparse array has a larger consecutive range than that constructed by using the ULA. In addition, the diff–sum coarray for the DsCAMpS contains more consecutive lags than that for the CAMpS with the same number of sensors.

5.2. MUSIC Spectra

Figure 7 presents the normalized MUSIC spectra of the uniform distributed signals detected by the four virtual configurations in Figure 6. Here, the frame length L, the number L ˜ of the frames in the VNCM method and the snapshots L s in the NC MUSIC satisfy L = L ˜ = L s = 800 . We consider the input SNR = 10 dB, and we suppose Q = 101 sources, which are uniformly distributed between 60 and 60 . It is noted that the frame length is actually the snapshots in one frame.
Figure 7d shows that the DsCAMpS with the VNCM used can detect all the 101 sources since the DOF of the diff–sum coarray for the DsCAMpS is 133. The corresponding MNDS is 130. According to Figure 7a–c, the other three virtual configurations fail to obtain the correct DOA estimations. This is because the numbers of the consecutive lags in the traditional NC virtual array for the ULA, the diff–sum coarrays for the ULA and the CAMpS are respectively 29, 57 and 101. Then, the corresponding MNDSs are respectively 28, 54 and 98. Thus, with the same number of sensors, the DsCAMpS with the VNCM used achieves a better performance than the ULA with the NC MUSIC used, the ULA with the VNCM used and the CAMpS with the VNCM used.
In order to demonstrate the estimated DOAs when the signal source distribution changes, in Figure 8, we further simulate the normalized MUSIC spectra of the non-uniformly distributed signals detected by the four virtual configurations. Here, Q = 101 sources are non-uniformly distributed between 60 and 60 , and the other parameters are the same as those in Figure 7. The non-uniform distribution of sources could make some DOAs get very close, which may deteriorate the DOA estimation performance. However, Figure 8d shows that the DsCAMpS with the VNCM used can still detect all the sources. In contrast, the other three virtual configurations fail to obtain the correct DOA estimations, which is shown in Figure 8a–c. Thus, regardless of whether the signal sources are uniformly distributed or non-uniformly distributed, the DsCAMpS with the VNCM used achieves a better performance than the other three virtual arrays.

5.3. Root Mean Square Error (RMSE)

We further conduct Monte Carlo simulations to compare the DOA estimation performance of the four virtual configurations in Figure 6. Here, we use the RMSE of the estimated DOAs as the performance metric. The RMSE is defined as
R M S E = 1 J Q j = 1 J q = 1 Q θ ^ q ( j ) θ q 2 ,
where J is the number of Monte Carlo simulations, θ q denotes the real DOA of the qth signal source and θ ^ q ( j ) denotes the estimate of θ q for the jth trial, j = 1 , , J . In all the simulations, we consider the signal source number Q = 25 and utilize 500 independent Monte Carlo simulations. Figure 9a depicts the RMSE performance as a function of the input SNR. In this simulation, we suppose L = L ˜ = L s = 500 . It is clear that all the arrays with the VNCM used outperform the ULA with the NC MUSIC used due to the diff–sum coarray containing more consecutive elements than the traditional NC virtual array. Among the three arrays with the VNCM used, the DsCAMpS achieves the best performance. In Figure 9b, we suppose SNR = 10 dB to compare the RMSE performance as a function of the snapshots. Here, we consider L = L ˜ = L s . From the results, all the arrays with the VNCM used still perform much better than the ULA with the NC MUSIC used. Compared with the RMSE of the ULA and the CAMpS with the VNCM used, the RMSE of the DsCAMpS with the VNCM used is smaller. The results of the two simulations suggest that the diff–sum coarray achieves a higher number of DOF than the traditional NC virtual array and the DsCAMpS is a novel array of which the diff–sum coarray has higher DOF than that of the ULA and the CAMpS.

6. Conclusions

We have proposed the VNCM method, which utilizes both the NC property and the concept of the coarray, to obtain a novel NC virtual array named as the diff–sum coarray. Due to comprising both the difference set and the sum set, the diff–sum coarray has a higher DOF than the traditional NC virtual array. Also, we utilize the quasi-stationary characteristic instead of the spatial smoothing method to solve the coherent issue generated by the KR product operation. Thus, the available DOFs would not be reduced by half. Taking the CAMpS as the array model, we have summarized the properties of the corresponding diff–sum coarray and then further proposed the DsCAMpS to achieve a higher number of DOFs. The high DOF of the diff–sum coarray and the performance of the novel method were numerically studied and evaluated.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant No. 61271113).

Author Contributions

The main idea was proposed by Zhenhong Chen. Zhenhong Chen and Zhiming Chen conceived and designed the experiments. Zhenhong Chen and Shiwei Ren performed the experiments and analyzed the data. Zhenhong Chen and Shiwei Ren wrote the paper and all the authors participated in amending the original manuscript. In the stage of the major revision, the design and simulation of the new experiments together with the revision of the manuscript were mainly completed by Yingtao Ding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proof of Proposition 1

(a) Define S c s l + = M n + N m 0 n P 2 N 1 , 0 m M 1 as the cross-sum set between S 2 and S 11 . Firstly, we would show that S c s l + contains all the consecutive lags in the range M 1 N 1 , P 2 M N 1 . Given an arbitrary integer x in the above range, we need to demonstrate there exist m 2 0 , M 1 and n 2 0 , P 2 N 1 such that x = M n 2 + N m 2 . According to [22], we know that there exist n and m , which can vary over all the integers, such that x = M n + N m . Then, we can further have x = M n + y N + N m y M , where y is an arbitrary integer. Denoting n 2 = n + y N and m 2 = m y M , we can choose the value of y to make m 2 satisfy 0 m 2 M 1 . Thus, we obtain 1 M M n 2 = x N m 2 P 2 M N 1 , that is, 0 n 2 P 2 N 1 . This demonstrates S c s l + contains all the consecutive lags in the range M 1 N 1 , P 2 M N 1 . Then, it is obvious that the consecutive range of the cross-sum set between S 2 and S 1 i is M 1 N 1 + i 1 M N , P 2 M N 1 + i 1 M N . When P 2 2 , the length of the consecutive range is greater than M N . Thus, the cross-sum set between S 2 and S 1 contains all the consecutive lags in the range M 1 N 1 , P 1 + P 2 1 M N 1 .
(b) In [16], we have proved that in the range M N , M N + M + N 1 , the holes in the cross-sum result M n + N m 0 n N 1 , 0 m M 1 between the two subarrays of the coprime array can be expressed as S p h 1 S p h 2 , where
S p h 1 = M n N n M N + M + N 1 / M
and
S p h 2 = N m M m M N + M + N 1 / N .
In the CAMpS, Period 1 of Subarray 1 is actually the first subarray of the coprime array and Period 1 of Subarray 2 is the second subarray of the coprime array. In addition, the elements in S 1 P 1 are equivalent to the elements in S 11 adding ( P 1 1 ) M N , and the elements in S 2 P 2 are equivalent to the elements in S 21 adding ( P 2 1 ) M N . Thus, in the range P 1 + P 2 1 M N , P 1 + P 2 1 M N + M + N 1 , the holes in the cross-sum set between S 1 P 1 and S 2 P 2 can be expressed as S m h 1 S m h 2 , where
S m h 1 = M n P 1 + P 2 1 N n P 1 + P 2 1 M N + M + N 1 / M
and
S m h 2 = N m P 1 + P 2 1 M m P 1 + P 2 1 M N + M + N 1 / N .
When P 2 > P 1 , the self-sum set of Subarray 2 contains all the elements in S m h 1 . However, any element in S m h 2 is larger than the maximum in the self-sum set of Subarray 1. Therefore, only the holes in the range P 1 + P 2 1 M N , P 1 + P 2 1 M N + N 1 can be filled by the self-sum set of the CAMpS. Then, the consecutive range of the non-negative sum set S s u m + is M 1 N 1 , P 1 + P 2 1 M N + N 1 when P 2 > P 1 .
(c) Denote
S p s u m + = M n + N m M n + M n N m + N m , 0 n N 1 , 0 m M 1
as the sum result for the prototype coprime array. In [16], we have proved that when M = 2 or N = 2 , the consecutive range of S p s u m + is M 1 N 1 , M N + N 1 or M 1 N 1 , M N + M 1 . Thus, when P 2 = P 1 , we can know that with M = 2 , the cross-sum result between S 1 P 1 and S 2 P 2 contains all the consecutive lags in the range ( P 1 + P 2 1 ) M N M N + 1 , P 1 + P 2 1 M N + N 1 . According to the property (a) of Proposition 1, it is obvious that the consecutive range of S s u m + is M 1 N 1 , P 1 + P 2 1 M N + N 1 . When N = 2 , the result becomes M 1 N 1 , P 1 + P 2 1 M N + M 1 . When M , N 2 , the set S p s u m + contains all the consecutive lags in the range M 1 N 1 , M N + M + N 1 . Then, we can utilize the similar method as the above procedure to demonstrate the consecutive range of S s u m + is M 1 N 1 , P 1 + P 2 1 M N + M + N 1 .
(d) In the proof of the property (a), we have shown that the consecutive range of the cross-sum set between S 2 and S 1 i is M 1 N 1 + i 1 M N , P 2 M N 1 + i 1 M N . It is easy to find that the minimum in the consecutive range of the cross-sum set between S 2 and S 1 P 1 is less than the maximum in the consecutive range of the cross-sum set between S 2 and S 11 . Thus, the consecutive range of the non-negative sum set for the CAMpS is the same as that of the non-negative sum set for a simplified CAMpS, of which the location sets of the two subarrays can be expressed as S 11 S 1 P 1 and S 2 .

References

  1. Chargé, P.; Wang, Y.; Saillard, J. A non-circular sources direction finding method using polynomial rooting. Signal Process. 2001, 81, 1765–1770. [Google Scholar] [CrossRef]
  2. Shen, L.; Liu, Z.; Gou, X.; Xu, Y. Polynomial-rooting based fourth-order MUSIC for direction-of-arrival estimation of noncircular signals. Syst. Eng. Electron. 2014, 25, 942–948. [Google Scholar] [CrossRef]
  3. Wang, Z.; Xiaofei, Z.; Huapu, S.; Renzheng, C. Non-circular generalised-ESPRIT algorithm for direction of arrival estimation. IET Radar Sonar Navig. 2017, 11, 736–744. [Google Scholar] [CrossRef]
  4. Steinwandt, J.; Roemer, F.; Haardt, M. Analytical performance assessment of esprit-type algorithms for coexisting circular and strictly non-circular signals. In Proceedings of the 2016 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Shanghai, China, 20–25 March 2016; pp. 2931–2935. [Google Scholar]
  5. Yin, J.; Wang, D.; Wu, Y.; Yao, X. ML-based single-step estimation of the locations of strictly noncircular sources. Digit. Signal Process. 2017, 69, 224–236. [Google Scholar] [CrossRef]
  6. Shi, H.; Leng, W.; Guan, Z.; Jin, T. Two novel two-stage direction of arrival estimation algorithms for two-dimensional mixed noncircular and circular sources. Sensors 2017, 17, 1433. [Google Scholar] [CrossRef] [PubMed]
  7. Wang, Q.; Zhu, X.; Chen, H.; Wang, L.; Yan, W.; Fang, H. Computationally efficient direction finding for a mixture of circular and strictly noncircular sources with uniform rectangular arrays. Sensors 2017, 17, 1269. [Google Scholar] [CrossRef] [PubMed]
  8. Wang, X.; Wang, W.; Li, X.; Liu, Q.; Liu, J. Sparsity-aware DOA estimation scheme for noncircular source in MIMO radar. Sensors 2016, 16, 539. [Google Scholar] [CrossRef] [PubMed]
  9. Abeida, H.; Delmas, J. MUSIC-like estimation of direction of arrival for noncircular sources. IEEE Trans. Signal Process. 2006, 54, 2678–2690. [Google Scholar] [CrossRef]
  10. Zoubir, A.; Chargé, P.; Wang, Y. Non circular sources localization with ESPRIT. In Proceedings of the European Conference on Wireless Technology (ECWT 2003), Munich, Germany, 8–10 October 2003. [Google Scholar]
  11. Haardt, M.; Romer, F. Enhancements of unitary ESPRIT for non-circular sources. In Proceedings of the 2004 IEEE International Conference on Acoustics, Speech, and Signal Processing, Montreal, QC, Canada, 17–21 May 2004. [Google Scholar]
  12. Xu, Y.; Liu, Z. Modified virtual spatial smoothing algorithm. Acta Electron. Sin. 2009, 37, 2646–2650. [Google Scholar]
  13. Liu, J.; Huang, Z.; Zhou, Y. Extended 2q-MUSIC algorithm for noncircular signals. Signal Process. 2008, 88, 1327–1339. [Google Scholar] [CrossRef]
  14. Krim, H.; Viberg, M. Two decades of array signal processing research: the parametric approach. IEEE Signal Process. Mag. 1996, 13, 67–94. [Google Scholar] [CrossRef]
  15. Liu, J.; Zhang, Y.; Lu, Y.; Wang, W. DOA estimation based on multi-resolution difference co-array perspective. Digit. Signal Process. 2017, 62, 187–196. [Google Scholar] [CrossRef]
  16. Wang, X.; Chen, Z.; Ren, S.; Cao, S. DOA estimation based on the difference and sum coarray for coprime arrays. Digit. Signal Process. 2017, 69, 22–31. [Google Scholar] [CrossRef]
  17. Qin, S.; Zhang, Y.D.; Amin, M.G. DOA estimation of mixed coherent and uncorrelated targets exploiting coprime MIMO radar. Digit. Signal Process. 2017, 61, 26–34. [Google Scholar] [CrossRef]
  18. Sun, F.; Wu, Q.; Sun, Y.; Ding, G.; Lan, P. An iterative approach for sparse direction-of-arrival estimation in co-prime arrays with off-grid targets. Digit. Signal Process. 2017, 61, 35–42. [Google Scholar] [CrossRef]
  19. Harry, L.; Trees, V. Optimum Array Processing: Part IV of Detection, Estimation, and Modulation Theory; John Wiler & Sons: New York, NY, USA, 2002. [Google Scholar]
  20. Moffet, A. Minimum-redundancy linear arrays. IEEE Trans. Antennas Propag. 1968, 16, 172–175. [Google Scholar] [CrossRef]
  21. Pal, P.; Vaidyanathan, P. Nested arrays: A novel approach to array processing with enhanced degrees of freedom. IEEE Trans. Signal Process. 2010, 58, 4167–4181. [Google Scholar] [CrossRef]
  22. Vaidyanathan, P.; Pal, P. Sparse sensing with co-prime samplers and arrays. IEEE Trans. Signal Process. 2011, 59, 573–586. [Google Scholar] [CrossRef]
  23. Liu, C.L.; Vaidyanathan, P.P. Super nested arrays: Linear sparse arrays with reduced mutual coupling–Part I: Fundamentals. IEEE Trans. Signal Process. 2016, 64, 3997–4012. [Google Scholar] [CrossRef]
  24. Qin, S.; Zhang, Y.; Amin, M. Generalized coprime array configurations for direction-of-arrival estimation. IEEE Trans. Signal Process. 2015, 63, 1377–1390. [Google Scholar] [CrossRef]
  25. Abramovich, Y.I.; Gray, D.A.; Gorokhov, A.Y.; Spencer, N.K. Positive-definite toeplitz completion in DOA estimation for nonuniform linear antenna arrays. I. Fully augmentable arrays. IEEE Trans. Signal Process. 1998, 46, 2458–2471. [Google Scholar] [CrossRef]
  26. Abramovich, Y.I.; Spencer, N.K.; Gorokhov, A.Y. Positive-definite toeplitz completion in DOA estimation for nonuniform linear antenna arrays. II. Partially augmentable arrays. IEEE Trans. Signal Process. 1999, 47, 1502–1521. [Google Scholar] [CrossRef]
  27. Pillai, S.; Haber, F. Statistical analysis of a high resolution spatial spectrum estimator utilizing an augmented covariance matrix. IEEE Trans. Acoust. Speech Signal Process. 1987, 35, 1517–1523. [Google Scholar] [CrossRef]
  28. Liu, C.L.; Vaidyanathan, P.P.; Pal, P. Coprime coarray interpolation for DOA estimation via nuclear norm minimization. In Proceedings of the 2016 IEEE International Symposium on Circuits and Systems (ISCAS), Montreal, QC, Canada, 22–25 May 2016; pp. 2639–2642. [Google Scholar]
  29. Pal, P.; Vaidyanathan, P. Coprime sampling and the MUSIC algorithm. In Proceedings of the 2011 IEEE Digital Signal Processing Workshop and IEEE Signal Processing Education Workshop (DSP/SPE), Sedona, AZ, USA, 4–7 January 2011; pp. 289–294. [Google Scholar]
  30. Ren, S.; Wang, W.; Chen, Z. DOA estimation exploiting unified coprime array with multiperiod subarrays. In Proceedings of the 2016 CIE International Conference on Radar (RADAR), Guangzhou, China, 10–13 October 2016. [Google Scholar]
  31. Wang, W.; Ren, S.; Chen, Z. Unified coprime array with multiperiod subarrays for direction-of-arrival estimation. Digit. Signal Process. 2018, 74, 30–42. [Google Scholar] [CrossRef]
  32. Ma, W.K.; Hsieh, T.H.; Chi, C.Y. DOA estimation of quasi-stationary signals with less sensors than sources and unknown spatial noise covariance: a Khatri–Rao subspace approach. IEEE Trans. Signal Process. 2010, 58, 2168–2180. [Google Scholar] [CrossRef]
  33. Shen, Q.; Liu, W.; Cui, W.; Wu, S. Extension of co-prime arrays based on the fourth-order difference co-array concept. IEEE Signal Process. Lett. 2016, 23, 615–619. [Google Scholar] [CrossRef]
  34. Wan, L.; Xie, L. An improved DOA estimation algorithm for circular and non-circular signals with high resolution. In Proceedings of the 2016 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Shanghai, China, 20–25 March 2016; pp. 3051–3055. [Google Scholar]
  35. Gao, F.; Nallanathan, A.; Wang, Y. Improved MUSIC under the coexistence of both circular and noncircular sources. IEEE Trans. Signal Process. 2008, 56, 3033–3038. [Google Scholar] [CrossRef]
  36. Chen, H.; Hou, C.; Liu, W.; Zhu, W.P.; Swamy, M.N.S. Efficient two-dimensional direction-of-arrival estimation for a mixture of circular and noncircular sources. IEEE Sens. J. 2016, 16, 2527–2536. [Google Scholar] [CrossRef]
  37. Chen, Z.; Ren, S.; Wang, W. DOA estimation exploiting extended co-array of coprime array. In Proceedings of the 2016 CIE International Conference on Radar (RADAR), Guangzhou, China, 10–13 October 2016. [Google Scholar]
Figure 1. The coprime array with multiperiod subarrays.
Figure 1. The coprime array with multiperiod subarrays.
Sensors 18 00344 g001
Figure 2. Two examples of the difference coarray for the CAMpS with ( M , N ) = ( 3 , 4 ) : (a) the set S c d with ( P 1 , P 2 ) = ( 2 , 3 ) ; (b) the set S d i f f with ( P 1 , P 2 ) = ( 2 , 3 ) ; (c) the set S c d with ( P 1 , P 2 ) = ( 3 , 3 ) ; and (d) the set S d i f f with ( P 1 , P 2 ) = ( 3 , 3 ) .
Figure 2. Two examples of the difference coarray for the CAMpS with ( M , N ) = ( 3 , 4 ) : (a) the set S c d with ( P 1 , P 2 ) = ( 2 , 3 ) ; (b) the set S d i f f with ( P 1 , P 2 ) = ( 2 , 3 ) ; (c) the set S c d with ( P 1 , P 2 ) = ( 3 , 3 ) ; and (d) the set S d i f f with ( P 1 , P 2 ) = ( 3 , 3 ) .
Sensors 18 00344 g002
Figure 3. Two examples of the non-negative sum coarray for the CAMpS with ( M , N ) = ( 3 , 4 ) : (a) the set S c s + with ( P 1 , P 2 ) = ( 2 , 3 ) ; (b) the set S s u m + with ( P 1 , P 2 ) = ( 2 , 3 ) ; (c) the set S c s + with ( P 1 , P 2 ) = ( 3 , 3 ) ; and (d) the set S s u m + with ( P 1 , P 2 ) = ( 3 , 3 ) .
Figure 3. Two examples of the non-negative sum coarray for the CAMpS with ( M , N ) = ( 3 , 4 ) : (a) the set S c s + with ( P 1 , P 2 ) = ( 2 , 3 ) ; (b) the set S s u m + with ( P 1 , P 2 ) = ( 2 , 3 ) ; (c) the set S c s + with ( P 1 , P 2 ) = ( 3 , 3 ) ; and (d) the set S s u m + with ( P 1 , P 2 ) = ( 3 , 3 ) .
Sensors 18 00344 g003
Figure 4. The diff–sum coprime array with multiperiod subarrays.
Figure 4. The diff–sum coprime array with multiperiod subarrays.
Sensors 18 00344 g004
Figure 5. The consisting sets of the diff–sum coarray for the DsCAMpS with M , N , P 2 = 4 , 3 , 3 : (a) the difference set S d i f f ; (b) the total sum set S s u m ; and (c) the diff–sum set S d s .
Figure 5. The consisting sets of the diff–sum coarray for the DsCAMpS with M , N , P 2 = 4 , 3 , 3 : (a) the difference set S d i f f ; (b) the total sum set S s u m ; and (c) the diff–sum set S d s .
Sensors 18 00344 g005
Figure 6. Four virtual configurations ( R = 15 ): (a) the traditional NC virtual array for the ULA; (b) the diff–sum coarray for the ULA; (c) the diff–sum coarray for the CAMpS with M , N , P 1 , P 2 = 4 , 3 , 2 , 3 ; and (d) the diff–sum coarray for the DsCAMpS with M , N , P 2 = 4 , 3 , 3 .
Figure 6. Four virtual configurations ( R = 15 ): (a) the traditional NC virtual array for the ULA; (b) the diff–sum coarray for the ULA; (c) the diff–sum coarray for the CAMpS with M , N , P 1 , P 2 = 4 , 3 , 2 , 3 ; and (d) the diff–sum coarray for the DsCAMpS with M , N , P 2 = 4 , 3 , 3 .
Sensors 18 00344 g006
Figure 7. Spatial spectra of the uniformly distributed signals detected by the four virtual configurations ( Q = 101 and SNR = 10 dB): (a) the ULA with the NC MUSIC used; (b) the ULA with the VNCM used; (c) the CAMpS with the VNCM used; and (d) the DsCAMpS with the VNCM used.
Figure 7. Spatial spectra of the uniformly distributed signals detected by the four virtual configurations ( Q = 101 and SNR = 10 dB): (a) the ULA with the NC MUSIC used; (b) the ULA with the VNCM used; (c) the CAMpS with the VNCM used; and (d) the DsCAMpS with the VNCM used.
Sensors 18 00344 g007
Figure 8. Spatial spectra of the non-uniformly distributed signals detected by the four virtual configurations ( Q = 101 and SNR = 10 dB): (a) the ULA with the NC MUSIC used; (b) the ULA with the VNCM used; (c) the CAMpS with the VNCM used; and (d) the DsCAMpS with the VNCM used.
Figure 8. Spatial spectra of the non-uniformly distributed signals detected by the four virtual configurations ( Q = 101 and SNR = 10 dB): (a) the ULA with the NC MUSIC used; (b) the ULA with the VNCM used; (c) the CAMpS with the VNCM used; and (d) the DsCAMpS with the VNCM used.
Sensors 18 00344 g008
Figure 9. Estimation precision ( Q = 25 ): (a) RMSE versus SNR; and (b) RMSE versus the number of snapshots.
Figure 9. Estimation precision ( Q = 25 ): (a) RMSE versus SNR; and (b) RMSE versus the number of snapshots.
Sensors 18 00344 g009
Table 1. The DOFs of the diff–sum coarray for the CAMpS.
Table 1. The DOFs of the diff–sum coarray for the CAMpS.
P 2 > P 1 1 P 2 = P 1 1
M , N 2 N > M = 2 M > N = 2 M , N > 2
DOF 2 s d s max 1 + 1 2 s d s max 1 + 1 2 s d s max 2 + 1 2 s d s max 3 + 1
Table 2. The DOFs of the diff–sum coarray for the DsCAMpS .
Table 2. The DOFs of the diff–sum coarray for the DsCAMpS .
N > M = 2 M > N = 2 M , N > 2
DOF 2 s d s max a + 1 2 s d s max b + 1 2 s d s max c + 1

Share and Cite

MDPI and ACS Style

Chen, Z.; Ding, Y.; Ren, S.; Chen, Z. A Novel Noncircular MUSIC Algorithm Based on the Concept of the Difference and Sum Coarray. Sensors 2018, 18, 344. https://doi.org/10.3390/s18020344

AMA Style

Chen Z, Ding Y, Ren S, Chen Z. A Novel Noncircular MUSIC Algorithm Based on the Concept of the Difference and Sum Coarray. Sensors. 2018; 18(2):344. https://doi.org/10.3390/s18020344

Chicago/Turabian Style

Chen, Zhenhong, Yingtao Ding, Shiwei Ren, and Zhiming Chen. 2018. "A Novel Noncircular MUSIC Algorithm Based on the Concept of the Difference and Sum Coarray" Sensors 18, no. 2: 344. https://doi.org/10.3390/s18020344

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop