1. Introduction
The problem of locating wideband emitters with sensor arrays is of growing interest in many applications such as sonars, radars, advanced satellites and cellular communication systems. Unfortunately, we cannot apply narrowband signal subspace methods such as multiple signal classification (MUSIC) [
1] and estimation of signal parameters via rotation invariance techniques ESPRIT [
2] on wideband sources because the phase difference between sensor outputs is no longer just dependent on the direction of arrival (DOA) alone, but also depends on the temporal frequency, which has a wide range.
The incoherent signal subspace method (ISM) [
3] decomposes wideband signals into many narrowband signals and uses narrowband methods on each decomposed narrowband signals. Then, the results from all the frequency bins are combined to create the final DOA estimation. The ISM method is simple and effective, but it cannot resolve coherent sources and suffers from low signal noise ratio (SNR). The coherent signal-subspace method (CSM) [
4] uses a transformation matrix to transform the correlation matrices from many frequency bins into one general correlation matrix at one focusing frequency. Then narrowband signal subspace methods can be applied to estimate the DOA. Although this method has been extended in [
5,
6] to decrease the resolution threshold and reduce the estimation bias, the method still requires initial DOA estimations and its estimation performance is very sensitive to these initial values. The test of orthogonality of projected subspaces (TOPS) [
7] and the test of orthogonality of frequency subspaces (TOFS) [
8] are newly introduced methods which estimate DOA by checking the orthogonality between the signal and noise subspaces of different frequency components of the sources. Both TOPS and TOFS are incoherent methods. The performance of TOPS is between that of the coherent and the incoherent methods in the whole SNR range. Recently, TOPS was extended in [
9,
10] to deal with coherent signals and achieve a higher accuracy. These algorithms perform satisfactorily when the received signals are simple, but this requirement is not guaranteed in real environments. Important parameters including density of signals, dynamic range, frequency range and bandwidth should be taken in consideration. That is why studies in which these algorithms are used in real environments are rare.
The DOA estimation is transformed into a sparse reconstruction problem in recent research [
11,
12,
13,
14]. These algorithms have a number of advantages, including increased resolution and improved robustness to noise. However, it is necessary to utilize iteration to achieve a prescribed accuracy in these methods. The intractable computational complexity for sparse recovery makes it difficult to implement on hardware platforms. To the best of our knowledge, there are no wideband DOA estimation methods based on sparse reconstruction that have been validated on hardware platforms. The main purpose of this paper is to present a hardware-efficient wideband DOA estimation method.
The digital channelization receiver shows several desirable characteristics, including high broadband instantaneous frequency coverage, good sensitivity and dynamic range, simultaneous signal detection, arbitration and parameter encoding, and fine frequency measurement [
15]. Digital channelization receivers play an important role in radar and electronic warfare applications. The channelization receiver partitions the received signal into a number of sub-channels and arbitrates which sub-channels the signal truly resides in. Afterwards, the parameters of the arbitrated channels including signal center frequency, pulse width and amplitude are then estimated, which is very important for DOA estimation. The channelization receiver architecture was first used to estimate DOA of ultra-wideband signal in [
16,
17,
18]. It splits the array output into multiple frequency channels and down-converts each channel into much lower frequency. Each channel can then be sampled at a lower sampling rate. Each branch of the proposed channelization receiver architecture consists of a mixer, Bessel band-pass filter and low sampling rate analog-to-digital converter (ADC). However, the frequency split and down-convert are implemented before digitization. The performance is highly sensitive to environmental factors, for example, the temperature.
With the development of high-speed ADCs, direct sample to wideband intermediate frequency (IF) signal conversion is possible. Recently, the Science and Technology for Automatic Target Recognition (ATR) Laboratory developed a digital array radar (DAR) test-bed using direct wideband radio frequency (RF) digitization [
19,
20]. The test-bed works at the IF of 0.6–3.0 GHz. It can achieve an instantaneous bandwidth of 500 MHz. Based on this test-bed, a new DOA estimate flow including IF digitization, digital channelization, channel selection and DOAs estimation is proposed. Firstly, the digital channelization receiver is used to split the wideband array output into multiple narrowband sub-channels. Secondly, the sub-channels within the bandwidth of wideband signals are selected out among the output sub-channels and signal parameters are estimated accurately. Finally, a channelization-based ISM method (Channelization-ISM) and a channelization-based TOPS (Channelization-TOPS) are proposed to estimate the DOA of wideband signals. Channelization-ISM first estimates the DOAs of each sub-channel independently by exploiting the narrowband signal subspace method. Then the estimated DOAs from each sub-channel are averaged to yield the final result. Channelization-TOPS measures the orthogonality between the signal and the noise subspaces of the output sub-channels to subsequently estimate DOAs. Based on the accurate signal parameters estimation, signals with different bandwidths are isolated reasonably. The output SNR is improved by removing the noises from the bandwidth of signals. The proposed channelization method outperforms the conventional ISM method and TOPS method, especially in real environments. Moreover, the output signals in a channelization receiver are in a parallel format, making it easy to implement on hardware.
2. The Wideband Signal Model
Assuming a uniform linear array (ULA) of
omnidirectional sensors receives
wideband signals
from unknown directions
. The
kth sensor output
is:
where
is the additive white Gaussian noise of
kth sensor which is assumed to be statistically independent with signal sources.
is the propagation delay of the
kth sensor:
where
is the speed of the signal propagation,
d is the spacing between adjacent sensors. The discrete-time Fourier transform (DTFT) of the
kth sensor output is:
Decomposing the sensor outputs into several narrowband signals using a filter bank or the discrete Fourier transform (DFT), the output can be written in vector format at
frequencies as:
where:
[·]
T denotes the transpose operator.
is the
K ×
P steering matrix. Its columns are the
K × 1 array manifolds
:
3. The Proposed Method
Inspired by the digital channelization receiver, we propose a new DOA estimate flow shown in
Figure 1, which is composed of IF digitization, digital channelization, channel selection and DOAs estimation. Detailed descriptions of each procedure are presented below. In the DOAs estimation, the channelization-based ISM method and channelization-based TOPS method are proposed.
3.1. IF Digitization and Digital Channelization
The wideband radar returns are directly digitized. A digital quadrature demodulation system [
21,
22,
23] is applied to get the amplitude and phase information of the received signals. A channelization receiver typically employs digital filter banks to extract multiple narrowband sub-channels from the received wideband signals [
24,
25,
26]. The architecture can be illustrated in
Figure 2.
represents the digitization of
kth sensor output
with the sampling rate
, where
is greater than twice of bandwidth of the input wideband signal.
means decimation. The array output is channelized into
frequency channels. The carrier frequency of the
i-th channel is at:
is a low-pass filter.
The output of
i-th sub-channel can be expressed as:
The order of prototype filter is
and
, where
is a positive integer. By polyphase decomposition with the prototype filter, Equation (10) can be rewritten as:
Define
and
. Then:
Substituting Equation (9) into Equation (12), we can rewrite
as:
where
represents convolution operation. Replacing
as
, then:
Figure 3 shows the optimized digital channelization receiver architecture, which is realized by implementing one low-pass filter and an inverse discrete Fourier transform (IDFT) module. The down sampling by
is moved to the left of the filtering operation.
to
represent the polyphase components of the low-pass filter
. This requires fewer components and is a more realistic architecture.
The output of digital channelization receiver
can be regarded as a narrowband signal sampled with a sampling rate of
. Assume that
N is the number of snapshots of
and
that of
, then we have the conclusion that
. The bandwidths of the output sub-channels are same. Different output bandwidths are required to deal with more complicated situation in real environments where signals with different unknown bandwidths coexist. The digital channelization receiver architecture can be extended to isolate channels with different bandwidths based on a multistage DFT filter [
27,
28].
3.2. Channelization Selection and Parameter Estimation
The frequency of received signals is not known in a real environment. Sometimes only part of the output sub-channels are within the wideband signal bandwidth. Before estimating the DOA, we need to find the sub-channels in which the wideband signals exist.
The output sub-channels within the bandwidth of wideband signals always result in higher power output. The sub-channel power is a criterion to determine whether there signals exist in the sub-channel. We can compare the power of all the sub-channels to select out the sub-channels within the wideband signal bandwidth.
Assume that
consists of a wideband signal
and white Gaussian noise
.
has a power of
and the power spreads on
sub-channels equally. The input SNR is:
The output
is composed of the response due to the noise
and the signal
, respectively. The noise response is a complex filtered baseband noise with an autocorrelation function
[
29]:
Only the
sub-channels within the wideband signal bandwidth are expected in further processing. The SNR of these sub-channels are:
The SNR is improved approximately times. The signal parameters of each selected sub-channel, including frequency, bandwidth and amplitude must be determined. The DOA of a narrow signal can be estimated with the only sub-channel which the signal resides in. However, wideband signals always span a number of sub-channels. We will introduce how to estimate the DOA of wideband signals in the next section.
3.3. Channelization-Based ISM Method
Assume that a total of
sub-channels within the wideband signal bandwidth are selected. Firstly, we use narrowband methods to estimate the DOAs of every selected sub-channel. In practice, the estimated
K ×
K covariance matrix can be expressed as:
where
denotes the complex conjugate transpose
. Performing an eigendecomposition on
, the eigenvalues
in descending order and their associated eigenvectors
can be obtained. Define matrices
and
, the columns of
and
span the signal subspace and the noise subspace respectively. DOAs of every sub-channels can be estimated via the spatial spectrum function:
Then we combine the estimated DOA
for the different sub-channels, and we can find the final DOA estimation. The first estimator exploits the arithmetic mean metric both for the individual sub-channels and the combination over the frequency range, which can be expressed as:
The second estimator uses the arithmetic mean metric for the individual sub-channels and the geometric mean metric for the combination over the frequency range:
The first estimator shows higher accuracy and lower resolution than the second one.
3.4. Channelized Based TOPS Method
Following the standard TOPS approach introduced in [
7], we first perform an eigendecomposition on
and obtain the noise subspace
. Select one output sub-channel as reference to obtain the reference signal subspace
.
Define a diagonal unitary transformation matrix whose size is
and diagonal elements are
:
can transform the array manifold of one frequency into that of another frequency without changing the DOA.
Define a
matrix
:
where
,
is a hypothetic arrival angle of signal:
is a full rank matrix if
is not equal to the arrival angle of signal otherwise
loses its rank. This theorem holds as long as the source signals are not fully correlated. The estimation of DOA can be obtained through a one-dimensional scan on
, i.e.,
where
is the smallest singular value of
.
3.5. Discussion
The channelization receiver shows several appealing characteristics such as wideband frequency coverage, high sensitivity and dynamic range, and accurate measurement of frequency resolution. The digital channelization receiver splits the wideband array output into multiple narrowband sub-channels and meanwhile estimates the parameters of the selected sub-channels. The estimation of the parameters including signal center frequency, pulse width and amplitude is very helpful to DOA estimation. Firstly, signals with different frequencies can be separated into different sub-channels. The DOA estimations of signals in different sub-channels can be obtained accurately without any influence of other signals. Secondly, to estimate the DOA of a signal, the knowledge of the signal parameters such as signal center frequency, pulse width and amplitude is necessary. The channelization receiver offers us a good way to obtain these parameters. Thirdly, the channelization receiver removes noises beyond the bandwidth of signals and selects the sub-channels within the bandwidth of signals with a high accuracy. The output SNR is improved approximately
times. Finally, the channelization receiver can be implemented in a parallel processing architecture, consequently it is able to tackle real-time processing on hardware platforms [
29,
30]. In virtue of these advantages, the proposed channelization-based method shows a better performance over the conventional ISM and TOPS methods. Compared with the method described in [
16], the ADC in the proposed method is closer to the antenna, which can significantly reduce the complexity and hardware cost.
The performances of the proposed method may degrade when the array manifold is affected by mutual coupling, sensor position errors and gain or phase uncertainties. Many calibration algorithms have been developed to compensate these uncalibration impairments [
31,
32,
33] In this paper we assume that the array manifold is deterministic and constant in time. Similar to TOPS, the proposed method works with arbitrary 1-D or 2-D arrays when the array manifolds are linearly independent. Only ULAs are used in the simulation for convenience.
In the DOA estimation procedure, we present the channelization-based ISM method and channelization-based TOPS method. The difference between these two methods is similar to that between the ISM method and TOPS, which is presented in [
31]. In this paper, we only focus on the improvement of the proposed channelization-based method instead of the conventional ISM method and TOPS method.
4. Simulation
In this section, some computer simulations are conducted to demonstrate the performance of proposed algorithm for DOA estimation of wideband signals. Consider a ULA composed of eight sensors. The maximum number of resolved sources is 7 (). Suppose two wideband signals with the same power located at and impinging on the ULA. The two sources have the same center frequency and the same bandwidth . The sampling frequency is chosen to be . The sensor noise is complex Gaussian noise.
The average root mean square error (RMSE) of the estimates from 200 Monte Carlo trials is defined as:
where
is the DOA estimation value of the incident signal at the
w-th Monte Carlo trial.
In the first simulation, the accuracy performances of the ISM method [
3] and the proposed Channelization-ISM method are compared. We also present the Cramér-Rao lower bound (CRLB) for the DOA estimations. The first estimator in Equation (17) is applied to yield the final results. The number of sub-channels in the Channelization-ISM method is the same as that of frequency bins in the ISM method, which is 64. The RMSEs against the SNR and the number of snapshots are plotted in
Figure 4 and
Figure 5, respectively. In
Figure 4, the SNR varies from 0 to 20 dB with a step size of 2 dB and the number of snapshots is 500. In
Figure 5, the number of snapshots varies from 50 to 1000 and the SNR is fixed at 10 dB. It is clear that with the increase of the SNR and the number of snapshots, the performances of both methods improves. The performance of the proposed method asymptotically approaches CRLB. The proposed Channelization-ISM method shows a better performance than the ISM method at an arbitrary level of SNR and an arbitrary number of snapshots. This is mainly because the channelization receiver has a higher frequency resolution. Signal parameters can be obtained accurately, which is very helpful for DOA estimation. The output SNR is improved by removing by removing the noises beyond the bandwidth of the signals. This result supports the discussions in
Section 3.5 and demonstrates the efficiency of our proposed method.
The second simulation considers the same scenario as the first simulation to compare the performances of the TOPS method [
7] and the proposed Channelization-TOPS method. Similar to the result of the first simulation,
Figure 6 and
Figure 7 show that the proposed Channelization-TOPS method always outperforms the TOPS method at any level of SNR and any number of snapshots. The proposed channelization-based method has a better DOA estimation accuracy than the conventional ISM and TOPS methods.
To study the detection capability of week signal with the proposed Channelization-ISM method, the resolution probabilities versus SNR for the proposed Channelization-ISM and ISM method are calculated. By definition, a signal is resolved successfully when
is smaller than 0.5, where
and
represent the estimated and true directions for the
kth incident signal, respectively. This simulation considers the scenario where there are two sources with one source (weak source) being weaker than the other one by 20 dB. The simulation parameters are the same with those of the second simulation.
Figure 8 shows the resolution probability of the two sources versus SNR that corresponds to the strong source with ISM and Channelization-ISM method. Source 1 is the weak source. The resolution probability for the weak source with ISM method is unsatisfactory, as it remains low even though the SNR increases to 6. The proposed Channelization-ISM can obtain similar performance with lower SNR. ISM suffers from bad performance degradation for the weak source while the proposed Channelization-ISM doesn’t. The results illustrate that the proposed channelization-based method performs better on the aspect of dynamic range.
Figure 9 which shows the comparison between the TOPS method and Channelization-TOPS also supports the conclusion.
In the last simulation, the performances of the proposed Channelization-ISM method and Channelization-TOPS method against the number of sub-channels are studied. The SNR is fixed at 10 dB and the number of snapshots is 256.
Figure 10 shows that with the increase of the number of sub-channels, the performances of the proposed methods improves. This is mainly because the sub-channel bandwidth becomes narrower with the increase of the number of sub-channels. Then, the influence over the phase difference between sensor outputs from the temporal frequency of sub-channels becomes weaker. These will lead to a better performance of the DOA estimation. However, the increase of the number of sub-channels may result an increase of the computational complexity. The number of channels used generally is limited when digital channelization is implemented on a hardware platform.