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 , , , and represent the operators of transpose, Hermitian transpose, inverse, and pseudo-inverse, respectively; denotes the Frobenius norm. returns the phase of a vector in radian. is to get the mathematical expectation of a variable. returns a diagonal matrix with the diagonal entities are the m-th row of .
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]
where
denotes the response entity of the
m-th sensor with respect to the
k-th
signal,
represents the speed of light,
represents the
k-th incident far-field narrowband signal,
and
are the DOA and carrier frequency of the
k-th signal.
denotes the noise signal of the
m-th sensor. Then the received signal is sampled using a one-bit ADC, i.e.,
where
is the complex value quantization processing of the received signal, such as
where if z > 0,
returns 1, otherwise returns −1.
and
represent the real and imaginary parts of the complex value, respectively.
The relationship between the covariance matrix
of
and the covariance matrix
of the one-bit measurement
is discussed as follows. First of all, the
m-th diagonal element of
is
where
and
represent the power of the
k-th signal and noise, the element correlation coefficient of the (
)-th position of
can be expressed as follows
where
is defined as
Since all the sensors are identical,
. Therefore, Equation (4) can be simplified to
Assuming that the powers of the signals are the same, i.e.,
, then we have
which implies
. More importantly, we find the value of both the real and imaginary parts of the correlation coefficient decrease as
decreases.
Next, we focus on
. It is obviously that
is a zero-mean, unit variance, i.e.,
,
. The (
m,
n)-th entity of
can be expressed as
It is obvious that if
m =
n,
. According to the inverse sine law [
37], we have
Hence the relationship between
and
is
It can be seen from Equation (11) that
where
,
is an unknown scaling factor. Furthermore, it is proven that [
32]
where
is an identity matrix. Although
is unknown, we can see from Equation (13) that
and
share the same eigenvectors. This is why the traditional subspace-based algorithm can be directly applied to
for spectrum estimation. Now we consider L one-bit snapshots are available, and the quantified data is arranged into matrix format as
where
denotes the direction matrix,
accounts for the steering vector,
is an unquantized signal measurement matrix,
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
where
and
are scalers,
is a unitary matrix,
,
. Generally,
. It should be pointed out that
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
are follows the sensor array. If
, the array signal from the
p-th delay units can be expressed as
where
Define the delay matrix as
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
where
is the (
)-th element of matrix
,
is the (
)-th element of matrix
and
is the ()-th element of matrix
.
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
is to solve
where
is the
p-th signal slice of
from the ‘source’ direction. The LS update for
(denoted by
) is then given by
Similarly, the
l-th slice of
from the ‘DOA’ direction can be expressed as
,
,
is the corresponding noise. The LS update for
(denoted by
) is given by
In addition, the slice of
from the ‘frequency’ direction can be formulated as
,
is the corresponding noise. The LS update to
(denoted by
) is
Before the first calculation of TLAS,
and
must be initialized. Usually,
and
are randomly initialized or initialized with ESPRIT or PM. Firstly, Equation (21) is computed to estimate
. Then, Equation (22) is calculated to update
(based on the initialized
and previously estimated
). 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,, andthat establish the PARAFAC model in (19). If the k-rank of,,and(denoted by,, and) satisfyThen the estimation of,
, andare unique up to permutation and scaling of columns, which can be expressed aswhereis a permutation matrix,,, andstand for the fitting errors, and,, andare diagonal scaling matrices with. The k-th column of the delay matrix
is
Thus, we can get the phase of
as
It is easy to find
where
where
is a scaler that we do not care about. Suppose the k-th column of
is
, and let the estimation of
is
. The LS solution of
is
From the second entity of
, we can get
. On the other hand, we define
Similarly, we have
with
where
is a uninteresting constant. Replace
with
, we can get
. Let the estimation of
be
, then the LS solution of
can be obtained via
and finally the DOA can be estimated via
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
,
and
, 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
s. The signal-to-noise ratio (SNR) is defined as
, where
and
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
where
is the estimated DOA or frequency of the
q-th Monte Carlo trial,
is the true value of the DOA or frequency,
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.