Next Article in Journal
Mapping and Tracking Forest Burnt Areas in the Indio Maiz Biological Reserve Using Sentinel-3 SLSTR and VIIRS-DNB Imagery
Next Article in Special Issue
SPICE-ML Algorithm for Direction-of-Arrival Estimation
Previous Article in Journal
Individual Tree Crown Delineation Using Multispectral LiDAR Data
Previous Article in Special Issue
Directional Modulation Technique Using a Polarization Sensitive Array for Physical Layer Security Enhancement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Joint Angle and Frequency Estimation Using One-Bit Measurements

1
National Demonstration Center for Experimental Electrical and Electronic Education, Yangtze University, Jingzhou 434023, China
2
State Key Laboratory of Marine Resource Utilization in South China Sea, Hainan University, Haikou 570228, China
3
National University of Defense Technology, Hefei 230037, China
4
Nanjing Marine Radar Institute, Nanjing 211153, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(24), 5422; https://doi.org/10.3390/s19245422
Submission received: 28 October 2019 / Revised: 2 December 2019 / Accepted: 3 December 2019 / Published: 9 December 2019

Abstract

:
Joint angle and frequency estimation is an important branch in array signal processing with numerous applications in radar, sonar, wireless communications, etc. Extensive attention has been paid and numerous algorithms have been developed. However, existing algorithms rely on accurately quantified measurements. In this paper, we stress the problem of angle and frequency estimation for sensor arrays using one-bit measurements. The relationship between the covariance matrices of one-bit measurement and that of the accurately quantified measurement is extended to the tensor domain. Moreover, a one-bit parallel factor analysis (PARAFAC) estimator is proposed. The simulation results show that the angle and frequency estimation can be quickly achieved and correctly paired.

1. Introduction

Source localization is one of the most important branches of array signal processing [1,2]. It has been actively conducted in the fields of communication, radar, sonar, seismic exploration, and cognitive radio [3,4,5,6,7]. Source localization using a sensor array always involves spectrum estimation, such as direction-of-arrival (DOA) estimation, delay estimation, frequency estimation, polarization estimation, or a combination of them. Among the enormous research topics, joint angle and frequency estimation are particularly prominent, since the two parameters are very important in various fields, and they can improve the detection ability and anti-interference ability of the spatial source signals. For instance, the two parameters can be adopted in passive radar systems for target locating and tracking; in space division multiple access-based wireless communications systems, these two parameters can be utilized to locate the user and allocate pilot tones; also, these two parameters are useful for channel estimation and beamforming. Moreover, the algorithms used for joint DOA and frequency estimation can be easily extended for angle-delay estimation, delay-frequency estimation, angle-delay-frequency estimation, etc., since problems of multiple parameter estimation using a sensor array are very similar to each other. Due to the above reasons, we focus on the problem of joint DOA and frequency estimation in this paper.
In the past decades, various spectrum estimation algorithms have been proposed. Typical algorithms including multiple signal classification (MUSIC) [8,9,10], estimating signal parameters via rotational invariance technique (ESPRIT) [11,12], propagator method [13], maximum likelihood (ML) [14,15], tensor-based approaches [16,17,18,19,20], and optimization-aware algorithms [21,22,23,24,25,26,27]. Generally speaking, MUSIC is computationally inefficient as it requires multiple peak search. Also, an ML estimator has high complexity due to exhaustive iteration. Unlike MUSIC and ML, ESPRIT can obtain closed-form solution, at the expense of decreased array aperture. Both MUSIC and ESPRIT need eigen decomposition to obtain the signal subspace or the noise subspace, and the complexity of eigen decomposition is on the third order of the matrix dimension. To avoid eigen decomposition, the propagator method has been introduced, which can obtain the subspaces via the least squares (LS) method. Optimization-aware algorithms always make sure accurate spectrum estimation; however, they are often too complex to be applied. Usually, the above-mentioned algorithm rely on matrix decomposition—approaches based on tensor decomposition are often superior than the matrix-based methods as they have better de-noising performance than the latter. A tensor can be viewed as a multi-dimensional (more than three) vector (matrix can be interpreted as a two-dimensional vector). Two tensor models are frequently used, namely parallel factor analysis (PARAFAC) and Tucker tensor. The former factorize a high-dimensional low-rank tensor into sums of rank-one tensors, the latter is highly analogy to multi-dimensional eigen decomposition.
To perform high-resolution spectrum estimation, some optional techniques are helpful, e.g., large-scaler sensor array, wideband signals, high-precision, analog-to-digital converter (ADC). However, these techniques would bring massive measurements and challenge the sampling system. On the one hand, high-precision ADCs are expensive and are of high energy consumption, which leads to an increased cost and complexity. On the other hand, massive measurements require more storage and processing resource; thus the real-time performance of the system is difficult to guarantee. To overcome these disadvantages, the concept of compressed sensing (CS) has been proposed [28,29], which may provide new inspirations for signal acquisition and processing. Several CS-based sampling frameworks have been developed, e.g., random sampler, random demodulation, modulated wideband converter, time encoding machine. Nevertheless, these architectures rely on high-precision quantification. Recently, the concept of one-bit CS [30], i.e., quantization with only one bit, has been proposed. With one-bit quantization, only one data bit need to be stored and processed, the system complexity can be reduced accordingly. Owing to its potential prospect, one-bit quantization has aroused much attentions in massive multiple-input multiple-output (MIMO) communications and DOA estimation [31]. Generally, spectrum estimation with one-bit measurements is linked to a sparse inverse problem, which is resolved via the optimization method. More recently, a one-bit MUSIC framework was driven in [32]. It is dproven that the covariance matrix with one-bit measurements can be approximated by a scaled unquantized covariance matrix, thus the traditional subspace algorithms can be directly applied. Besides, many efforts have been devoted to the sparse recovery problem from one-bit measurement [33,34]. In addition, some works have been done to the waveform design problems in MIMO radar with the problem of one-bit DAC [35].
It should be pointed out that, as mentioned previously, one-bit quantization is a lossy compression method; thus, the performance of spectrum estimation algorithms with one-bit quantization are suffering from degradation. For performance enhancement, this paper tries to integrate the tensor approach with joint DOA and frequency estimation in the presence of one-bit measurement. To this end, a new spatial-time sampling framework is presented, in which a one-bit ADCs are adopted. The relationship between the covariance matrices of one-bit measurement and that of the unquantized covariance matrix is extended to the tensor domain. A one-bit PARAFAC algorithm is proposed for joint DOA and frequency estimation, in which closed-form and automatically paired parameters are achieved. Compared with the one-bit ESPRIT algorithm, the proposed algorithm offers a more accurate estimation performance. Numerical simulations verify the effectiveness of the proposed algorithm.
The rest of the paper is organized as follows. In Section 2, we present the signal model and analyze the impact of measurement process on signal noise. Section 3 provides the proposed ASCS scheme. The simulation results are given in Section 4. Finally, conclusions are given in Section 5.
Notation: Lower case and capital letters in bold denote, respectively, vectors and matrices. The superscript ( ) T , ( ) H , ( ) 1 , and ( ) represent the operators of transpose, Hermitian transpose, inverse, and pseudo-inverse, respectively; F denotes the Frobenius norm. angle ( · ) returns the phase of a vector in radian. E [ ] is to get the mathematical expectation of a variable. D m ( A ) returns a diagonal matrix with the diagonal entities are the m-th row of A .

2. Signal Model and One-Bit Quantization

Herein, we consider a uniform linear array scenario with M-element sensors, the inter-element interval is d. Suppose that there are K uncorrelated narrow-band sources appearing on the far-field of the array. Besides, the signal and noise are uncorrelated, and both of them are modeled as independent, zero-mean, circular, complex Gaussian random processes, then the received signal of the m-th sensor can be expressed as [36]
y m ( t ) = k = 1 K a m ( θ k , f k ) s k ( t ) + n m ( t ) ( m = 1 , , M )
where a m ( θ k , f k ) = e j 2 π ( m 1 ) d f k sin ( θ k ) / c denotes the response entity of the m-th sensor with respect to the k-th ( 1 k K ) signal, c represents the speed of light, s k ( t ) represents the k-th incident far-field narrowband signal, θ k and f k are the DOA and carrier frequency of the k-th signal. n m ( t ) denotes the noise signal of the m-th sensor. Then the received signal is sampled using a one-bit ADC, i.e.,
x m ( t ) = O ( y m ( t ) )
where O ( · ) is the complex value quantization processing of the received signal, such as
O ( z ) = 1 2 ( s i g n ( r e a l ( z ) ) + j s i g n ( i m a g ( z ) ) )
where if z > 0, s i g n ( z ) returns 1, otherwise returns −1. r e a l ( z ) and i m a g ( z ) represent the real and imaginary parts of the complex value, respectively.
The relationship between the covariance matrix R y of y m ( t ) and the covariance matrix R x of the one-bit measurement x m ( t ) is discussed as follows. First of all, the m-th diagonal element of R y is δ y m 2
δ y m 2 = [ R y ] m m = k = 1 K | a m ( θ k , f k ) | 2 δ k 2 + δ n 2
where δ k 2 and δ n 2 represent the power of the k-th signal and noise, the element correlation coefficient of the ( m , n )-th position of R y can be expressed as follows
ρ y m y n = E [ y m ( t ) y n * ( t ) ] δ y m δ y n = [ R y ] m n [ R y ] m m [ R y ] n n = k = 1 K a m ( θ k , f k ) a n * ( θ k , f k ) δ k 2 k = 1 K | a m ( θ k , f k ) | 2 δ k 2 + δ n 2 k = 1 K | a n ( θ k , f k ) | 2 δ k 2 + δ n 2 = k = 1 K a m ( θ k , f k ) a n * ( θ k , f k ) ξ k k = 1 K | a m ( θ k , f k ) | 2 ξ k + 1 k = 1 K | a m ( θ k , f k ) | 2 ξ k + 1
where ξ k is defined as
ξ k = δ k 2 δ n 2
Since all the sensors are identical, | a m ( θ k , f k ) | = 1 . Therefore, Equation (4) can be simplified to
ρ = [ R y ] m m = k = 1 K δ k 2 + δ n 2
Assuming that the powers of the signals are the same, i.e., ξ k = ξ , then we have
ρ y m y n = k = 1 K e j 2 π f ( τ m ( θ i ) τ n ( θ i ) ) K + ξ 1
which implies | ρ y m y n | < 1 . More importantly, we find the value of both the real and imaginary parts of the correlation coefficient decrease as ξ decreases.
Next, we focus on R x . It is obviously that x m ( t ) is a zero-mean, unit variance, i.e., E [ x m ( t ) ] = 0 , δ 2 x m = 1 . The (m,n)-th entity of R x can be expressed as
ρ x m x n = E [ x m ( t ) x n * ( t ) ] δ x m δ x n = [ R x ] m n
It is obvious that if m = n, ρ x m x n = [ R y ] m m = 1 . According to the inverse sine law [37], we have
ρ x m x n = 2 arcsin e ( ρ y m y n ) π 2 π ( arcsin ( r e a l { ρ y m y n } ) + j arcsin ( i m a g { ρ y m y n } ) )
Hence the relationship between R x and R y is
R x = 2 π arcsin e ( 1 a R y )
It can be seen from Equation (11) that
R y = a sin e ( π 2 R x )
where sin e ( z ) sin ( r e a l { z } ) + j sin ( i m a g { z } ) , a is an unknown scaling factor. Furthermore, it is proven that [32]
R x 2 R y a π + ( 1 2 π ) I
where I is an identity matrix. Although a is unknown, we can see from Equation (13) that R x and R y share the same eigenvectors. This is why the traditional subspace-based algorithm can be directly applied to R x for spectrum estimation. Now we consider L one-bit snapshots are available, and the quantified data is arranged into matrix format as
X = O ( A S ˜ T + N )
where A = [ a ( θ 1 , f 1 ) , a ( θ 2 , f 2 ) , , a ( θ K , f K ) ] M × K denotes the direction matrix, a ( θ k , f k ) = [ 1 , a 2 ( θ k , f k ) , , a M ( θ k , f k ) ] T accounts for the steering vector, S ˜ L × K is an unquantized signal measurement matrix, N is the array noise sample matrix. We know it that the signal subspace obtained from eigenvalue decomposition of the covariance matrix span the same subspace of the signal subspace achieved from singular value decomposition. Combined with the result of in Equations (13) and (14) can be rewritten as
X 0 ε A S ˜ T U + β N U = ε A S T + β E 0
where ε and β are scalers, U L × L is a unitary matrix, S = U T S ˜ , E 0 = N U . Generally, K < min { M , P , L } . It should be pointed out that E 0 is a Gaussian white noise matrix, as proven in [32].

3. The Proposed Framework

In this paper, the delay-based sampling framework is proposed for joint DOA and delayed estimation. As shown in Figure 1, the P delay units τ p ( p = 1 , 2 , , P ) are follows the sensor array. If 0 < τ 1 < τ 2 < < 1 / max ( f k ) , the array signal from the p-th delay units can be expressed as
X P = ε A Φ P S T + β E P , p = 1 , 2 , , P
where
Φ P = [ e j 2 π f 1 τ p e j 2 π f 2 τ p e j 2 π f K τ p ] K × K
Define the delay matrix as
Φ = [ 1 1 1 e j 2 π f 1 τ 1 e j 2 π f 2 τ 1 e j 2 π f K τ 1 e j 2 π f 1 τ P e j 2 π f 2 τ P e j 2 π f K τ P ] ( P + 1 ) × K
Then Equation (16) defines a PARAFAC slicing model of the array output. In addition, outputs from all the delay units can be expressed as a third-order tensor χ with the (m,l,p)-th entity given by
χ m , n , p = ε k = 1 K a m , k s ˜ l , k ϕ q , k + e m , n , q ( m = 1 , 2 , , M ; l = 1 , 2 , , L ; q = 1 , 2 , , P + 1 )
where a m , f is the ( m , k )-th element of matrix A , s ˜ l , k is the ( l , k )-th element of matrix S and ϕ q , k  is the ()-th element of matrix Φ . e m , n , q is the associate noise measurement.
The trilinear alternating least squares (TALS) algorithm is a very popular technique for trilinear models. The basic principle of TALS is to update one-factor matrix via least squares (LS) technique while treating other factor matrices as known parameters. Based on the previous estimation, TALS update the residual matrices successively. The above iterations will repeat until the converge conditions have been satisfied. The LS fitting for S is to solve
min A , S , Φ [ X 0 X 1 X P ] - [ A D 1 ( Φ ) A D 2 ( Φ ) A D P + 1 ( Φ ) ] S T F
where X p ( p = 1 , 2 , , P ) is the p-th signal slice of χ from the ‘source’ direction. The LS update for S (denoted by S ^ ) is then given by
S ^ T = [ A D 1 ( Φ ) A D 2 ( Φ ) A D P + 1 ( Φ ) ] [ X 0 X 1 X P ]
Similarly, the l-th slice of χ from the ‘DOA’ direction can be expressed as Y l = Φ D l ( S ) A T +   N ˜ l , l = 1 , 2 , , L , N ˜ l is the corresponding noise. The LS update for A (denoted by A ^ ) is given by
A ^ T = [ Φ D 1 ( S ) Φ D 2 ( S ) Φ D L ( S ) ] [ Y 1 Y 2 Y L ]
In addition, the slice of χ from the ‘frequency’ direction can be formulated as Z m ( t ) = S D m ( A ) Φ T + N ¯ m ( m = 1 , 2 , , M ) , N ¯ m is the corresponding noise. The LS update to Φ (denoted by Φ ^ ) is
Φ ^ T = [ S D 1 ( A ) S D 2 ( A ) S D M ( A ) ] [ Z 1 Z 2 Z M ]
Before the first calculation of TLAS, A and Φ must be initialized. Usually, A and Φ are randomly initialized or initialized with ESPRIT or PM. Firstly, Equation (21) is computed to estimate S . Then, Equation (22) is calculated to update A (based on the initialized Φ and previously estimated S ). Thereafter, Equation (23) is computed to update Φ . Finally, iterations in Equations (21)–(23) will repeat until convergence. In this paper, we adopted the COMFAC algorithm for PARAFAC decomposition [38], which can be quickly converge after only a few iteration.
It is well known to us that matrix decompositions are usually not unique unless some constrains are enforced. Unlike matrix decomposition, tensor decompositions are often unique under mild conditions. The following Theorem 1 gives the uniqueness of PARAFAC analysis.
Theorem 1.
Consider the matrices A , Φ , and S that establish the PARAFAC model in (19). If the k-rank of A , Φ ,and S (denoted by k r ( A ) , k r ( Φ ) , and k r ( S ) ) satisfy
k r ( A ) + k r ( Φ ) + k r ( S ) 2 K + 2
Then the estimation of A , Φ , and S are unique up to permutation and scaling of columns, which can be expressed as
{ A ^ = A Π Δ 1 + E 1 Φ ^ = Φ Π Δ 2 + E 2 S ^ = S Π Δ 3 + E 3
where Π is a permutation matrix, E 1 , E 2 , and E 3 stand for the fitting errors, and Δ 1 , Δ 2 , and Δ 3 are diagonal scaling matrices with Δ 1 Δ 2 Δ 3 = I k .
The k-th column of the delay matrix Φ is
g ( f k ) = [ 1 , e j 2 π f k τ 1 , , e j 2 π f k τ P ] T
Thus, we can get the phase of g ( f k ) as
h = - angle ( g ( f k ) ) = [ 0 , 2 π f k τ 1 , , 2 π f k τ P ] T
It is easy to find
P 1 b = h
where
P 1 = [ 1 0 1 2 π τ 1 1 2 π τ P ]
where b 0 is a scaler that we do not care about. Suppose the k-th column of Φ ^ is g ^ ( f k ) , and let the estimation of h is h ^ . The LS solution of b is
b ^ = P 1 h ^
From the second entity of b ^ , we can get f ^ k . On the other hand, we define
u = a n g l e ( a ( θ k , f k ) ) = [ 0 , 2 π d f k sin θ k / c , , 2 π d ( M 1 ) f k sin θ k / c ] T
Similarly, we have
P 2 c = u
with
P 2 = [ 1 0 1 2 π d f k / c 1 ( M 1 ) 2 π d f k / c ]
where c 1 is a uninteresting constant. Replace f k with f ^ k , we can get P 2 . Let the estimation of u be u ^ , then the LS solution of c can be obtained via
c ^ = P 2 u ^
and finally the DOA can be estimated via
θ ^ i = arcsin ( e ^ 1 )

4. Simulation Results and Discussions

In this section, numerical simulations are carried out to verify the effectiveness of the proposed framework. In the simulations, we consider there are K = 3 uncorrelated far-field narrow-band signals with DOA-frequency pairs are ( 10 , 0.5   MHz ) , ( 20 , 0.7   MHz ) and ( 30 , 0.9   MHz ) , and L snapshots are collected. An M-element ULA with half-wavelength spacing is adopted to receive the incoming signals. Assume there are 3 uniform delay units, with the delay interval is 10 7 s. The signal-to-noise ratio (SNR) is defined as 10 lg ( P s / P n ) ( dB ) , where P s and P n are the powers with respect to signal and noise counterparts, respectively. The root mean square error (RMSE) is utilized for performance assessment. Herein, RMSE is defined as
RMSE = 1 Q q = 1 Q ( ϑ q ϑ 0 ) 2
where ϑ q is the estimated DOA or frequency of the q-th Monte Carlo trial, ϑ 0 is the true value of the DOA or frequency, Q is the total number of Monte Carlo trials.
Firstly, we illustrate Q = 500 scatter results of the one-bit PARAFAC framework in Figure 2 and Figure 3, where M = 12 and L = 1000 are considered, and SNRs are set to 5 dB and 10 dB, respectively. It is shown that both DOA and frequency of the sources can be estimated and correctly paired from one-bit measurements. Moreover, it seems the estimated accuracy can be improved with increasing SNR, as the scatter results are more concentrated with larger SNR.
Secondly, the RMSE performance of the proposed algorithm with various SNR is depicted in Figure 4 and Figure 5, where M = 12 and L = 1000. For the performance comparison, the performance of traditional ESPRIT with one-bit measurement (marked with O-ESPRIT) and the traditional ESPRIT with un-quantified measurement (marked with U-ESPRIT) are added. RMSE with respect to DOA estimation and frequency estimation are shown in Figure 4 and Figure 5, respectively. It is seen that at low SNR regions (SNR < 0), there is no visible performance difference between the three algorithms. Notably, there is a performance gap between the algorithm with un-quantified measurements and that with one-bit measurement, this is because the one-bit measurement is lossy. As the SNR increases, the performance gap becomes larger. Another interesting finding is that when SNR is larger than a given threshold (SNR = 10), RMSE will not decrease with the increasing SNR. Besides, the proposed O-PARAFAC algorithm achieves better RMSE than O-ESPRIT, since the tensor structure has been exploited.
Thirdly, the RMSE curves with different snapshot number L are given in Figure 6 and Figure 7, where M = 12 and SNR = 0 dB are considered. Clearly, RMSE on DOA estimation and frequency estimation would improve with L increasing. Similarly, O-PARAFAC algorithm provides more accurate DOA and frequency estimation than O-ESPRIT. Also, both algorithms offer higher RMSE than the U-ESPRIT, since the one-bit measurement is lossy.
Fourthly, the RMSE performances versus sensor number M are depicted in Figure 8, Figure 9 and Figure 10, where L = 1000 and SNR = 0 dB. It is shown that estimation performance would improve with the increasing M. A similar observation can be seen that the performance corresponding to O-PARAFAC is between that of the O-ESPRIT and U-ESPRIT. In addition, one can observe that the proposed algorithm requires less calculation time then O-ESPRIT and U-ESPRIT when M is larger than 60, which implies that the proposed algorithms is much more efficient than the compared algorithms is the presence of massive antennas.
Finally, the RMSE curves with various delay P are depicted in Figure 11 and Figure 12, respectively, where M = 12 and SNR = 0 dB are considered. Similar to our previous findings, RMSE on DOA estimation and frequency estimation would improve with P increasing. In addition, the performance of O-PARAFAC is better than O-ESPRIT but worse than U-ESPRIT.

5. Conclusions

In this article, we stress the one-bit quantization problem in joint DOA and frequency estimation using a sensor array. A one-bit PARAFAC framework has been accomplished, in which the relationship of one-bit quantization and de-quantization measurement has been extended to tensor domain. Simulation results show the feasibility of the proposed framework. Since the one-bit quantization is lossy, the performance of the traditional subspace-based algorithms may degrade. However, the proposed one-bit PARAFAC can still offer good parameter estimation accuracy. As we all know, the one-bit quantization system is more flexible than the high-precision quantization and is less sensitive to storage capacity than the latter; the one-bit quantization system will have a bright prospect in further applications.

Author Contributions

Z.L. carried out the experiments and wrote the first draft. J.S. helped to revise the original manuscript. X.W. helps to revise the manuscript. F.W. participated in the design of this algorithm. All authors read and approved the final manuscript.

Funding

This research was funded by yhr National Natural Science Foundation of China, grant number 61701046, and 61871218.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wu, F.; Wu, T.; Yuce, M. An Internet-of-Things (IoT) Network System for Connected Safety and Health Monitoring Applications. Sensors 2018, 19, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Mamun, M.A.A.; Yuce, M.R. Sensors and Systems for Wearable Environmental Monitoring Toward IoT-Enabled Applications: A Review. IEEE Sens. J. 2019, 19, 7771–7788. [Google Scholar] [CrossRef]
  3. Kim, B.; Kim, S.; Lee, J. A novel DFT-based DOA estimation by a virtual array extension using simple multiplications for FMCW radar. Sensors 2018, 18, 1560. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Ahmed, A.; Arablouei, R.; de Hoog, F.; Kusy, B.; Jurdak, R.; Bergmann, N. Estimating angle-of-arrival and time-of-flight for multipath components using wifi channel state information. Sensors 2018, 18, 1753. [Google Scholar] [CrossRef] [Green Version]
  5. Hsu, K.C.; Kiang, J.F. Joint Estimation of DOA and Frequency of Multiple Sources with Orthogonal Coprime Arrays. Sensors 2019, 19, 335. [Google Scholar] [CrossRef] [Green Version]
  6. Suksiri, B.; Fukumoto, M. An Efficient Framework for Estimating the Direction of Multiple Sound Sources Using Higher-Order Generalized Singular Value Decomposition. Sensors 2019, 19, 2977. [Google Scholar] [CrossRef] [Green Version]
  7. Lee, S.; Kim, S.C. Logarithmic-Domain Array Interpolation for Improved Direction of Arrival Estimation in Automotive Radars. Sensors 2019, 19, 2410. [Google Scholar] [CrossRef] [Green Version]
  8. Zhou, C.; Shi, Z.; Gu, Y.; Shen, X. DECOM: DOA Estimation with Combined MUSIC for Coprime Array. In Proceedings of the International Conference on Wireless Communications and Signal Processing (WCSP), Hangzhou, China, 24–26 October 2013; pp. 1–5. [Google Scholar]
  9. Liao, B. Fast angle estimation for MIMO radar with nonorthogonal waveforms. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 2091–2096. [Google Scholar] [CrossRef]
  10. Wang, X.; Wan, L.; Huang, M.; Shen, C.; Zhang, K. Polarization channel estimation for circular and non-circular signals in massive MIMO systems. IEEE J. STSP 2019, 13, 1001–1016. [Google Scholar] [CrossRef]
  11. Zhou, C.; Zhou, J. Direction-of-Arrival Estimation with Coarray ESPRIT for Coprime Array. Sensors 2017, 17, 1779. [Google Scholar] [CrossRef] [Green Version]
  12. Wen, F.; Shi, J.; Zhang, Z. Direction finding for bistatic MIMO radar with unknown spatially colored noise. Circuits Syst. Signal Process. 2019. [Google Scholar] [CrossRef]
  13. Dong, Y.; Dong, C.; Xu, J.; Zhao, G. Computationally Efficient 2-D DOA Estimation for L-Shaped Array With Automatic Pairing. IEEE Antennas Wirel. Propag. Lett. 2016, 15, 1669–1672. [Google Scholar] [CrossRef]
  14. Meng, X.; Jiang, Z. A Generalized Sparse Bayesian Learning Algorithm for one-bit DOA Estimation. IEEE Commun. Lett. 2018, 22, 1414–1417. [Google Scholar] [CrossRef]
  15. Wang, H.; Wan, L.; Dong, M.; Ota, K.; Wang, X. Assistant Vehicle Localization Based on Three Collaborative Base Stations via SBL-based Robust DOA Estimation. IEEE Internet Things J. 2019, 6, 5766–5777. [Google Scholar] [CrossRef]
  16. Wang, X.; Wang, W.; Li, X.; Wang, J. A tensor-based subspace approach for bistatic MIMO radar in spatial colored noise. Sensors 2014, 14, 3897–3907. [Google Scholar] [CrossRef] [Green Version]
  17. Guo, Y.; Wang, X.; Wang, W.; Huang, M.; Shen, C.; Cao, C.; Bi, G. Tensor-Based Angle Estimation Approach for Strictly Noncircular Sources with Unknown Mutual Coupling in Bistatic MIMO Radar. Sensors 2018, 18, 2788. [Google Scholar] [CrossRef] [Green Version]
  18. Wen, F.; Zhang, Z.; Wang, K.; Sheng, G.; Zhang, G. Angle estimation and mutual coupling self-calibration for ULA-based bistatic MIMO radar. Signal Process. 2018, 144, 61–67. [Google Scholar] [CrossRef]
  19. Wen, F.; Mao, C.; Zhang, G. Direction finding in MIMO radar with large antenna arrays and nonorthogonal waveforms. Digit. Signal Process. 2019, 94, 75–83. [Google Scholar] [CrossRef]
  20. Wen, F.; Zhang, Z.; Zhang, G. Joint DOD and DOA estimation for bistatic MIMO radar: A covariance trilinear decomposition perspective. IEEE Access 2019, 7, 53273–53283. [Google Scholar] [CrossRef]
  21. Wang, J.; Wang, X.; Xu, D.; Bi, G. Robust angle estimation for MIMO radar with the coexistence of mutual coupling and colored noise. Sensors 2018, 18, 832. [Google Scholar] [CrossRef] [Green Version]
  22. Wang, X.; Huang, M.; Wu, X.; Bi, G. Direction of arrival estimation for MIMO radar via unitary nuclear norm minimization. Sensors 2017, 17, 939. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Shi, Z.; Zhou, C.; Gu, Y.; Goodman, N.A.; Qu, F. Source estimation using coprime array: A sparse reconstruction perspective. IEEE Sens. J. 2016, 17, 755–765. [Google Scholar] [CrossRef]
  24. Zhou, C.; Gu, Y.; He, S.; Shi, Z. A robust and efficient algorithm for coprime array adaptive beamforming. IEEE Trans. Veh. Technol. 2017, 67, 1099–1112. [Google Scholar] [CrossRef]
  25. Zhou, C.; Gu, Y.; Zhang, Y.D.; Shi, Z.; Jin, T.; Wu, X. Compressive sensing-based coprime array direction-of-arrival estimation. IET Commun. 2017, 11, 1719–1724. [Google Scholar] [CrossRef]
  26. Zhou, C.; Gu, Y.; Fan, X.; Shi, Z.; Mao, G.; Zhang, Y.D. Direction-of-arrival estimation for coprime array via virtual array interpolation. IEEE Trans. Signal Process. 2018, 66, 5956–5971. [Google Scholar] [CrossRef]
  27. Zhou, C.; Gu, Y.; Shi, Z.; Zhang, Y.D. Off-grid direction-of-arrival estimation using coprime array interpolation. IEEE Signal Process. Lett. 2018, 25, 1710–1714. [Google Scholar] [CrossRef]
  28. Candès, E.J.; Wakin, M.B. An introduction to compressive sampling. IEEE Signal Process. Mag. 2008, 25, 21–30. [Google Scholar] [CrossRef]
  29. Duarte, M.F.; Davenport, M.A.; Takhar, D.; Laska, J.N.; Sun, T.; Kelly, K.F.; Baraniuk, R.G. Single-pixel imaging via compressive sampling. IEEE Signal Process. Mag. 2008, 25, 83–91. [Google Scholar] [CrossRef] [Green Version]
  30. Yan, M.; Yang, Y.; Osher, S. Robust 1-bit compressive sensing using adaptive outlier pursuit. IEEE Trans. Signal Process. 2012, 60, 3868–3875. [Google Scholar]
  31. Liu, C.L.; Vaidyanathan, P.P. One-bit sparse array DOA estimation. In Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), New Orleans, LA, USA, 5–9 March 2017; pp. 3126–3130. [Google Scholar]
  32. Huang, X.; Liao, B. One-Bit MUSIC. IEEE Signal Process. Lett. 2019, 26, 961–965. [Google Scholar] [CrossRef] [Green Version]
  33. Xiao, P.; Liao, B.; Li, J. One-bit compressive sensing via Schur-concave function minimization. IEEE Trans. Signal Process. 2019, 67, 4139–4151. [Google Scholar] [CrossRef]
  34. Xiao, P.; Liao, B. Robust one-bit compressive sensing with weighted l1-norm minimization. Signal Process. 2019, 164, 380–385. [Google Scholar] [CrossRef]
  35. Cheng, Z.; Liao, B.; He, Z.; Li, J. Transmit signal design for large-scale MIMO system with 1-bit DACs. IEEE Trans. Wirel. Commun. 2019, 18, 4466–4478. [Google Scholar] [CrossRef]
  36. Zhang, X.; Wang, D.; Xu, D. Novel blind joint direction of arrival and frequency estimation for uniform linear array. Prog. Electromagn. Res. 2008, 86, 199–215. [Google Scholar] [CrossRef] [Green Version]
  37. Jacovitti, G.; Neri, A. Estimation of the autocorrelation function of complex Gaussian stationary processes by amplitude clipped signals. IEEE Trans. Inf. Theory 1994, 40, 239–245. [Google Scholar] [CrossRef]
  38. Wen, F.; Shi, J.; Zhang, Z. Joint 2D-DOD, 2D-DOA and polarization angles estimation for bistatic EMVS-MIMO radar via PARAFAC analysis. IEEE Trans. Veh. Technol. 2019. [Google Scholar] [CrossRef]
Figure 1. Architecture of the spatial-temporal sampling.
Figure 1. Architecture of the spatial-temporal sampling.
Sensors 19 05422 g001
Figure 2. Scatter results with SNR = 5 dB.
Figure 2. Scatter results with SNR = 5 dB.
Sensors 19 05422 g002
Figure 3. Scatter results with SNR = 10 dB.
Figure 3. Scatter results with SNR = 10 dB.
Sensors 19 05422 g003
Figure 4. RMSE vs. SNR for DOA estimation.
Figure 4. RMSE vs. SNR for DOA estimation.
Sensors 19 05422 g004
Figure 5. RMSE vs. SNR for frequency estimation.
Figure 5. RMSE vs. SNR for frequency estimation.
Sensors 19 05422 g005
Figure 6. RMSE vs. L for DOA estimation.
Figure 6. RMSE vs. L for DOA estimation.
Sensors 19 05422 g006
Figure 7. RMSE vs. L for frequency estimation.
Figure 7. RMSE vs. L for frequency estimation.
Sensors 19 05422 g007
Figure 8. RMSE vs. M for DOA estimation.
Figure 8. RMSE vs. M for DOA estimation.
Sensors 19 05422 g008
Figure 9. RMSE vs. M for frequency estimation.
Figure 9. RMSE vs. M for frequency estimation.
Sensors 19 05422 g009
Figure 10. Average running time vs. M.
Figure 10. Average running time vs. M.
Sensors 19 05422 g010
Figure 11. RMSE vs. P for DOA estimation.
Figure 11. RMSE vs. P for DOA estimation.
Sensors 19 05422 g011
Figure 12. RMSE vs. P for frequency estimation.
Figure 12. RMSE vs. P for frequency estimation.
Sensors 19 05422 g012

Share and Cite

MDPI and ACS Style

Li, Z.; Shi, J.; Wang, X.; Wen, F. Joint Angle and Frequency Estimation Using One-Bit Measurements. Sensors 2019, 19, 5422. https://doi.org/10.3390/s19245422

AMA Style

Li Z, Shi J, Wang X, Wen F. Joint Angle and Frequency Estimation Using One-Bit Measurements. Sensors. 2019; 19(24):5422. https://doi.org/10.3390/s19245422

Chicago/Turabian Style

Li, Zeyang, Junpeng Shi, Xinhai Wang, and Fangqing Wen. 2019. "Joint Angle and Frequency Estimation Using One-Bit Measurements" Sensors 19, no. 24: 5422. https://doi.org/10.3390/s19245422

APA Style

Li, Z., Shi, J., Wang, X., & Wen, F. (2019). Joint Angle and Frequency Estimation Using One-Bit Measurements. Sensors, 19(24), 5422. https://doi.org/10.3390/s19245422

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