Next Article in Journal
A Novel Wearable Sensor-Based Human Activity Recognition Approach Using Artificial Hydrocarbon Networks
Previous Article in Journal
Exchange Bias Tuning for Magnetoresistive Sensors by Inclusion of Non-Magnetic Impurities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Channelization-Based DOA Estimation Method for Wideband Signals

Science and Technology on Automatic Target Recognition Laboratory (ATR), National University of Defense Technology, Changsha 410073, China
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(7), 1031; https://doi.org/10.3390/s16071031
Submission received: 4 May 2016 / Revised: 12 June 2016 / Accepted: 28 June 2016 / Published: 4 July 2016
(This article belongs to the Section Physical Sensors)

Abstract

:
In this paper, we propose a novel direction of arrival (DOA) estimation method for wideband signals with sensor arrays. The proposed method splits the wideband array output into multiple frequency sub-channels and estimates the signal parameters using a digital channelization receiver. Based on the output sub-channels, a channelization-based incoherent signal subspace method (Channelization-ISM) and a channelization-based test of orthogonality of projected subspaces method (Channelization-TOPS) are proposed. Channelization-ISM applies narrowband signal subspace methods on each sub-channel independently. Then the arithmetic mean or geometric mean of the estimated DOAs from each sub-channel gives the final result. Channelization-TOPS measures the orthogonality between the signal and the noise subspaces of the output sub-channels to estimate DOAs. The proposed channelization-based method isolates signals in different bandwidths reasonably and improves the output SNR. It outperforms the conventional ISM and TOPS methods on estimation accuracy and dynamic range, especially in real environments. Besides, the parallel processing architecture makes it easy to implement on hardware. A wideband digital array radar (DAR) using direct wideband radio frequency (RF) digitization is presented. Experiments carried out in a microwave anechoic chamber with the wideband DAR are presented to demonstrate the performance. The results verify the effectiveness of the proposed method.

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 K omnidirectional sensors receives P ( P K ) wideband signals s p ( t ) ( p = 1 , 2 , P ) from unknown directions { θ 1 , θ 2 , , θ P } . The kth sensor output x k ( t ) is:
x k ( t ) = p = 1 P s p ( t τ k ( θ p ) ) + n k ( t ) k = 1 , 2 , K
where n k ( t ) is the additive white Gaussian noise of kth sensor which is assumed to be statistically independent with signal sources. τ k ( θ p ) is the propagation delay of the kth sensor:
τ k ( θ p ) = ( k 1 ) d sin θ p / c
where c 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:
X k ( ω ) = p = 1 P S p ( ω ) e j ω τ k ( θ p ) + N k ( ω )
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 J frequencies as:
X ( ω i ) = A ( ω i ) S ( ω i ) + N ( ω i ) i = 0 , , J 1
where:
X ( ω i ) = [ X 1 ( ω i ) , X 2 ( ω i ) , , X K ( ω i ) ] T
S ( ω i ) = [ S 1 ( ω i ) , S 2 ( ω i ) , , S P ( ω i ) ] T
[·]T denotes the transpose operator. A ( ω i ) is the K × P steering matrix. Its columns are the K × 1 array manifolds a ( ω i , θ p ) :
A ( ω i ) = [ a ( ω i , θ 1 ) , a ( ω i , θ 2 ) , , a ( ω i , θ P ) ]
a ( ω i , θ p ) = [ e j ω i τ 1 ( θ p ) , e j ω i τ 2 ( θ p ) , , e j ω i τ K ( θ p ) ] T

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. x k ( n ) represents the digitization of kth sensor output x k ( t ) with the sampling rate f s , where f s is greater than twice of bandwidth of the input wideband signal. means decimation. The array output is channelized into J frequency channels. The carrier frequency of the i-th channel is at:
ω i = 2 π i J
h L P ( n ) is a low-pass filter.
The output of i-th sub-channel can be expressed as:
y k i ( m ) = [ x k ( n ) e j ω i n ] h L P ( n ) | n = m J = r x k ( n r ) e j ω i ( n r ) h L P ( r ) | n = m J = r x k ( m J r ) e j ω i ( m J r ) h L P ( r )
The order of prototype filter is D and D = L × J 1 , where L is a positive integer. By polyphase decomposition with the prototype filter, Equation (10) can be rewritten as:
y k i ( m ) = r x k ( m J r ) e j ω i ( m J r ) h L P ( r ) = p = 0 J 1 l = 0 L 1 x k ( m J l J p ) e j ω i ( m J l J p ) h L P ( l J + p )
Define x k p ( m ) = x k ( m J p ) and h p ( m ) = h L P ( m J + p ) . Then:
y k i ( m ) = p = 0 J 1 l = 0 L 1 x k ( m J l J p ) e j ω i ( m J l J ) h L P ( l J + p ) e j ω i p = p = 0 J 1 [ l = 0 L 1 x k p ( m l ) e j ω i ( m l ) J h p ( l ) ] e j ω i p
Substituting Equation (9) into Equation (12), we can rewrite y k i ( m ) as:
y k i ( m ) = p = 0 J 1 [ l = 0 L 1 x k p ( m l ) e j ω i ( m l ) J h p ( l ) ] e j ω i p = p = 0 J 1 [ l = 0 L 1 x k p ( m l ) e j 2 π i J ( m l ) J h p ( l ) ] e j 2 π i J p = p = 0 J 1 [ l = 0 L 1 x k p ( m l ) h p ( l ) ] e j 2 π i J p = p = 0 J 1 [ x k p ( m ) h p ( m ) ] e j 2 π i J p
where represents convolution operation. Replacing x k p ( m ) h p ( m ) as x ^ k p ( m ) , then:
y k i ( m ) = J I D F T ( x ^ k p ( m ) )
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 J is moved to the left of the filtering operation. h 0 ( m ) to h J 1 ( m ) represent the polyphase components of the low-pass filter h L P ( n ) . This requires fewer components and is a more realistic architecture.
The output of digital channelization receiver y k i ( m ) can be regarded as a narrowband signal sampled with a sampling rate of f s / J . Assume that N is the number of snapshots of x k ( n ) and M that of y k i ( m ) , then we have the conclusion that M = N / J . 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 x k ( n ) consists of a wideband signal s k ( n ) and white Gaussian noise n k ( n ) N ( 0 , σ n 2 ) . s k ( n ) has a power of A 2 and the power spreads on S ( S < J ) sub-channels equally. The input SNR is:
SNR in = A 2 σ n 2
The output y k i ( m ) is composed of the response due to the noise y k n i ( m ) and the signal y k s i ( m ) , respectively. The noise response is a complex filtered baseband noise with an autocorrelation function γ y k n [29]:
γ y k n = E [ y k n y k n * ] = σ n 2 n = 0 N | h k ( m ) | 2 σ n 2 J
Only the S sub-channels within the wideband signal bandwidth are expected in further processing. The SNR of these sub-channels are:
SNR out = A 2 J σ n 2 S
The SNR is improved approximately J / S ( J / S > 1 ) 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 S 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:
R ^ ( ω i ) = 1 M m = 1 M y ( ω i , m ) y H ( ω i , m ) , i = 0 , 1 , , S 1
where [ ] H denotes the complex conjugate transpose y ( ω i , m ) = [ y 1 i ( m ) , y 2 i ( m ) , , y K i ( m ) ] T . Performing an eigendecomposition on R ^ ( ω i ) , the eigenvalues λ 1 , λ 2 , , λ K in descending order and their associated eigenvectors u 1 , u 2 , , u K can be obtained. Define matrices F ( ω i ) = [ u 1 , , u P ] and U ( ω i ) = [ u P + 1 , , u K ] , the columns of F ( ω i ) and U ( ω i ) span the signal subspace and the noise subspace respectively. DOAs of every sub-channels can be estimated via the spatial spectrum function:
P S U B ( ω i , θ ) = 1 U T ( ω i ) a ( ω i , θ ) 2
Then we combine the estimated DOA θ i 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:
P C h a n n e l i z e d I S M 1 ( ω i , θ ) = 1 1 S i = 0 S 1 1 P S U B ( ω i , θ )
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:
P C h a n n e l i z e d I S S M 2 ( ω i , θ ) = 1 ( i = 0 S 1 1 P S U B ( ω i , θ ) ) 1 / S
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 R ^ ( ω i ) and obtain the noise subspace U ( ω i ) . Select one output sub-channel as reference to obtain the reference signal subspace F ( ω 0 ) .
Define a diagonal unitary transformation matrix whose size is k × k and diagonal elements are [ T ( ω i , θ i ) ] k , k :
[ T ( ω i , θ i ) ] k , k = exp { j ω i ( k 1 ) d sin θ i / c }
[ T ( ω i , θ i ) ] k , k can transform the array manifold of one frequency into that of another frequency without changing the DOA.
Define a K × P matrix F i ( ϕ ) :
F i ( ϕ ) = T ( Δ ω i , ϕ ) F ( ω 0 ) , i = 1 , 2 , S 1
where Δ ω i = ω i ω 0 , ϕ is a hypothetic arrival angle of signal:
D ( ϕ ) = [ F 1 H U ( ω 1 ) | F 2 H U ( ω 2 ) | F S 1 H U ( ω S 1 ) ]
D ( ϕ ) is a full rank matrix if ϕ is not equal to the arrival angle of signal otherwise D ( ϕ ) 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.,
θ ^ = arg max ϕ 1 σ min ( ϕ )
where σ min ( ϕ ) is the smallest singular value of D ( ϕ ) .

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 J / S ( J / S > 1 ) 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 ( K 1 ). Suppose two wideband signals with the same power located at θ 1 = 22 ° and θ 2 = 28 ° impinging on the ULA. The two sources have the same center frequency f 0 = 2.65   GHz and the same bandwidth B = 400   MHz . The sampling frequency is chosen to be f s = 1.2   GHz . 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:
R M S E = 1 200 t = 1 200 ( θ w θ ) 2
where θ w 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 | θ ^ k θ k | is smaller than 0.5, where θ ^ k and θ k 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.

5. Experiments

In this section, we present experiments carried out in a microwave anechoic chamber to further demonstrate the performance of the proposed algorithm.

5.1. The Architecture of DAR Test-Bed

Figure 11 shows the architecture of the wideband DAR test-bed, which is composed of wideband antennas, RF front modules, digital receiver, digital signal processor, built-in calibration circuitry and data recorder.
The wideband antenna is realized using eight tapered slot antennas (TSAs) shown in Figure 12. Figure 13 shows the architecture of the RF front-end, which is composed of low noise amplifier (LNA), programmable attenuator, band pass filters, power splitter and RF switch. Signals from the RF front-end are directly digitized by the wideband digital receiver. As shown in Figure 14, the digital receiver consists of ADC, large-scale field programmable gate array (FPGA) and multichannel optical modules. It is capable of directly sampling four input signals up to 3 GHz at the sampling rate of 1.2 GHz utilizing band-pass sampling theories. The digital signal processor shown in Figure 15 is based on the latest generation of high-performance multicore fixed and floating-point DSP TMS320C6678 and Xilinx Virtex-6 family’s high-density, FPGA, which is applicable for high-speed data communication, processing. The data recorder is designed to store kinds of data for system test and post-processing.

5.2. Experiment Carried Out in Microwave Anechoic Chamber

The experiment is carried out in a microwave anechoic chamber. As shown in Figure 16, the antenna array is placed on a rotation platform. Wideband signals impinging on the array are transformed into parallel waves through a reflector which is shown in Figure 17. There center frequency f0 = 2.7 GHz and bandwidth B = 500 MHz. The gain-phase and mutual coupling errors are calibrated by the method given in [34]. The wideband signals are directly sampled by the 10 bits ADCs with a sampling rate of 1.2 Gsps utilizing the band-pass sampling theory. The complex signal data speed is 600 MHz after digital quadrature demodulation. Then we split the signal into 32 sub-channels with the FPGA on digital signal platform. The output data speed of each sub-channel is 18.75 MHz, which is lower than the highest frequency of the Virtex-6 FPGA. Besides, the hard resources consumption is calculated as shown in Table 1. The results illustrate that the hardware platform adequately meets the resources constraint.
Figure 18 shows the amplitude of received signals in time domain. Eight curves with different colors represent the different sensor outputs of the ULA. Figure 19 depicts the received signals in the frequency domain. The spatial spectrum estimations of the received signals with different methods are plotted in Figure 20 and Figure 21. The true DOA could be measured according to the records of the rotation platform with a precision of 0.05°. It is clear that the DOA estimations of channelization-based methods are more approximate to the true DOA. The performances of the proposed channelization-based methods are better than those of the conventional ISM and TOPS methods. Rotating the array from −50° to 50° at the step of 1°, 101 groups of original data are obtained. DOA estimations with different incident angles are calculated using different methods. Table 2 presents the RMSEs of 101 estimates with different methods. The proposed channelization methods outperform the conventional ISM and TOPS methods. The experimental results are in good agreement with the simulation results.

6. Conclusions

In this paper, we propose a hardware-efficient, channelization-based DOA estimation method for wideband signals. It splits the wideband array output into multiple frequency sub-channels and estimate signal parameters of each sub-channel. The sub-channels within the wideband signal bandwidth are selected on the basis of the power of every sub-channel. The channelization-ISM method applies narrow techniques to the selected channels and gets the final DOA estimate by averaging the results from all the selected channels. The Channelization-TOPS method measures the orthogonality between the signal and the noise subspaces of the output sub-channels to estimate DOAs. The proposed channelization methods outperform the conventional ISM method and TOPS method, which results from the proposed channelization method being able to estimate the signal parameters more accurately and improve the output SNR more. Besides, the parallel processing architecture makes it easy to implement on hardware. Simulations and experiments in a microwave anechoic chamber are carried out as well. The results illustrate that the proposed method can achieve better estimation accuracy and dynamic range performance than the conventional ISM method and TOPS method.

Acknowledgments

The authors are thankful to the anonymous reviewers and the editors for their help in improving the manuscript and the ATR Lab for collecting the data used in this research.

Author Contributions

All authors contributed extensively to this paper. Rui Guo and Yue Zhang provided the main idea of this work. Qianqiang Lin and Yue Zhang conceived and designed the experiments. Rui Guo performed the experiments and wrote the paper. Zengping Chen modified the manuscript and provided many valuable suggestions.

Conflicts of Interest

The authors declare no conflict of interests.

References

  1. Schmidt, R.O. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef]
  2. Roy, R.; Kailath, T. ESPRIT—Estimation of signal parameters via rotational invariance techniques. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 984–995. [Google Scholar] [CrossRef]
  3. Wax, M.; Shan, T.; Kailath, T. Spatio-temporal spectral analysis by eigenstructure methods. IEEE Trans. Acoust. Speech Signal Process. 1984, 32, 817–827. [Google Scholar] [CrossRef]
  4. Wang, H.; Kaveh, M. Coherent signal-subspace processing for the detection and estimation of angles of arrival of multiple wideband sources. IEEE Trans. Acoust. Speech Signal Process. 1985, 33, 823–831. [Google Scholar] [CrossRef]
  5. Hung, H.; Kaveh, M. Focussing matrices for coherent signal-subspace processing. IEEE Trans. Acoust. Speech Signal Process. 1988, 36, 1272–1281. [Google Scholar] [CrossRef]
  6. Valaee, S.; Kabal, P. Wideband array processing using a two-sided correlation transformation. IEEE Trans. Signal Process. 1988, 43, 160–172. [Google Scholar] [CrossRef]
  7. Yeo-sun, Y.; Kaplan, L.M.; McClellan, J.H. TOPS: New DOA estimator for wideband signals. IEEE Trans. Signal Process. 2006, 54, 1977–1989. [Google Scholar] [CrossRef]
  8. Yu, H.Q.; Huang, Z.T.; Liu, J.; Zhou, Y.Y. TOFS: A New DOA Estimator for Wideband Sources. J. Astronaut. 2007, 28, 1304–1309. [Google Scholar]
  9. Guo, R.; Li, W.; Zhang, Y.; Chen, Z. DOA Estimation of Coherent Wideband Signals Based on Extended TOPS Algorithm. Proc. SPIE 2015, 9812. [Google Scholar] [CrossRef]
  10. Zhang, J.; Dai, J.; Ye, Z. An extended TOPS algorithm based on incoherent signal subspace method. Signal Process. 2010, 90, 3317–3324. [Google Scholar] [CrossRef]
  11. Sha, Z.; Liu, Z.; Huang, Z.; Zhou, Y. Covariance-based direction-of-arrival estimation of wideband coherent chirp signals via sparse representation. Sensors 2013, 13, 11490–11497. [Google Scholar] [CrossRef] [PubMed]
  12. Chen, T.; Wu, H.; Guo, L.; Liu, L. A Modified rife algorithm for off-grid DOA estimation based on sparse representations. Sensors 2015, 15, 29721–29733. [Google Scholar] [CrossRef] [PubMed]
  13. Dai, J.; Hu, N.; Xu, W.; Chang, C. Sparse Bayesian Learning for DOA Estimation with Mutual Coupling. Sensors 2015, 15, 26276–26280. [Google Scholar] [CrossRef] [PubMed]
  14. Wang, L.; Zhao, L.; Bi, G.; Wan, C.; Zhang, L.; Zhang, H. Novel wideband DOA estimation based on sparse bayesian learning with dirichlet process priors. IEEE Trans. Signal Process. 2016, 64, 275–289. [Google Scholar] [CrossRef]
  15. Anderson, G.W.; Webb, D.C.; Spezio, A.E.; Lee, J.N. Advanced channelization technology for RF, microwave and millimeterwave applications. Proc. IEEE 1991, 7, 1296–1297. [Google Scholar]
  16. Lie, J.P.; See, C.M.; Ng, B.P. Ultra wideband direction finding using digital channelization receiver architecture. IEEE Commun. Lett. 2006, 10, 85–87. [Google Scholar]
  17. Lie, J.P.; Ng, B.P.; See, C.M.; Ma, C. Multiple UWB emitters DOA estimation employing TH-SS. In Proceedings of the 2005 International Conference on Wireless Communications, Networking and Mobile Computing, Wuhan, China, 23–26 September 2005; Volume 1, pp. 336–339.
  18. Lie, J.P.; Ng, B.P.; See, C.M. Multiple UWB emitters DOA estimation employing time hopping spread spectrum. Prog. Electromagn. Res. 2008, 78, 83–101. [Google Scholar] [CrossRef]
  19. Li, W.X.; Zhang, Y.; Lin, J.Z.; Chen, Z.P. IF digitization receiver of wideband digital array radar test-bed. In Proceedings of the SPIE International Society Option Engineering, Amsterdam, The Netherland, 24–26 September 2014.
  20. Zhang, Y.; Bao, Q.; Wu, J.; Li, S. Design and implementation of wideband all digital array radar test-bed. In Proceedings of the 2014 European Radar Conference (EuRAD), Rome, Italy, 8–10 October 2014; pp. 617–620.
  21. Ho, K.C.; Chan, Y.T.; Inkol, R. A digital quadrature demodulation system. IEEE Trans. Aerosp. Electron. Syst. 1996, 32, 1218–1227. [Google Scholar] [CrossRef]
  22. Liu, H.; Ghafoor, A.; Stockman, P.H. A new quadrature sampling and processing approach. IEEE Trans. Aerosp. Elecrrun. Syst. 1989, 25, 733–748. [Google Scholar] [CrossRef]
  23. Zhang, G.; I-Khaiili, D.A.; Inkol, R.; Saper, R. Novel approach to the design of I/Q demodulation filters. IEEE Proc. Vis. Image Signal Proc. 1994, 141, 154–160. [Google Scholar] [CrossRef]
  24. Namgoong, W. A channelized digital ultra wideband receiver. IEEE Trans. Wirel. Commun. 2003, 2, 502–510. [Google Scholar] [CrossRef]
  25. Mahesh, R.; Vinod, A.P.; Lai, E.M.-K.; Omondi, A. Filter bank channelizers for multi-standard software defined radio receivers. J. Signal Process. Syst. 2011, 62, 157–171. [Google Scholar] [CrossRef]
  26. Yung, W.H.; Jian, M.; Ho, Y.W. Polyphase decomposition channelizers for software radios. In Proceedings of the IEEE International Symposium on ISCAS, Geneva, Switzerland, 18–31 May 2000; Volume 2, pp. 353–356.
  27. Michael, N.; Vinod, A.P.; Moy, C.; Palicot, J. Design of multistandard channelization accelerators for software defined radio handsets. IEEE Trans. Signal Process. 2011, 59, 4767–4780. [Google Scholar] [CrossRef]
  28. Fung, C.Y.; Chan, S.C. A multistage filterbank-based channelizer and its multiplier-less realization. In Proceedings of the IEEE International Symposium on Circuits and Systems, Phoenix, AZ, USA, 26–29 May 2002; Volume 3, pp. 429–432.
  29. Zahirniak, D.R.; Sharpin, D.L.; Fields, T.W. A hardware-efficient, multirate, digital channelized receiver architecture. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 137–152. [Google Scholar] [CrossRef]
  30. Xu, S.; Liu, G.; Gao, M. Design and implementation of digital channelized receiver in multi-fpga. In Proceedings of the 1st International Conference on Information Science and Engineering, Nanjing, China, 26–28 December 2009.
  31. Inghelbrecht, V.; Verhaevert, J.; van Hecke, T.; Rogier, H. The influence of random element displacement on DOA estimates obtained with (Khatri-Rao-) Root-MUSIC. Sensors 2014, 14, 21258–21280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Wang, W.; Ren, S.; Ding, Y.; Wang, H. An efficient algorithm for direction finding against unknown mutual coupling. Sensors 2014, 14, 20064–20077. [Google Scholar] [CrossRef] [PubMed]
  33. Guo, R.; Li, W.; Zhang, Y.; Chen, Z. New direction of arrival estimation of coherent signals based on reconstructing matrix under unknown mutual coupling. J. Appl. Remote Sens. 2016, 10, 1–11. [Google Scholar] [CrossRef]
  34. See, C. Sensor array calibration in the presence of mutual coupling and unknown sensor gains and phases. Electron. Lett. 1994, 30, 373–374. [Google Scholar] [CrossRef]
Figure 1. Flow of the proposed method.
Figure 1. Flow of the proposed method.
Sensors 16 01031 g001
Figure 2. Digital channelization receiver architecture.
Figure 2. Digital channelization receiver architecture.
Sensors 16 01031 g002
Figure 3. The optimized digital channelization receiver architecture.
Figure 3. The optimized digital channelization receiver architecture.
Sensors 16 01031 g003
Figure 4. RMSE of the DOA estimates against the SNR.
Figure 4. RMSE of the DOA estimates against the SNR.
Sensors 16 01031 g004
Figure 5. RMSE of the DOA estimates against the number of snapshots.
Figure 5. RMSE of the DOA estimates against the number of snapshots.
Sensors 16 01031 g005
Figure 6. RMSE of the DOA estimates against the SNR.
Figure 6. RMSE of the DOA estimates against the SNR.
Sensors 16 01031 g006
Figure 7. RMSE of the DOA estimates against the number of snapshots.
Figure 7. RMSE of the DOA estimates against the number of snapshots.
Sensors 16 01031 g007
Figure 8. Resolution probability versus SNR with ISM and Channelization-ISM method.
Figure 8. Resolution probability versus SNR with ISM and Channelization-ISM method.
Sensors 16 01031 g008
Figure 9. Resolution probability versus SNR with TOPS and Channelization-TOPS method.
Figure 9. Resolution probability versus SNR with TOPS and Channelization-TOPS method.
Sensors 16 01031 g009
Figure 10. RMSE of the DOA estimates against the number of sub-channels.
Figure 10. RMSE of the DOA estimates against the number of sub-channels.
Sensors 16 01031 g010
Figure 11. The architecture of the wideband DAR test-bed.
Figure 11. The architecture of the wideband DAR test-bed.
Sensors 16 01031 g011
Figure 12. The wideband antenna.
Figure 12. The wideband antenna.
Sensors 16 01031 g012
Figure 13. The RF front-end.
Figure 13. The RF front-end.
Sensors 16 01031 g013
Figure 14. The digital receiver.
Figure 14. The digital receiver.
Sensors 16 01031 g014
Figure 15. The digital signal platform.
Figure 15. The digital signal platform.
Sensors 16 01031 g015
Figure 16. Experiment in a microwave anechoic chamber.
Figure 16. Experiment in a microwave anechoic chamber.
Sensors 16 01031 g016
Figure 17. The reflector.
Figure 17. The reflector.
Sensors 16 01031 g017
Figure 18. Received signals expressed in time domain.
Figure 18. Received signals expressed in time domain.
Sensors 16 01031 g018
Figure 19. Received signals expressed in frequency domain.
Figure 19. Received signals expressed in frequency domain.
Sensors 16 01031 g019
Figure 20. DOA estimations with ISM method and Channelization-ISM method.
Figure 20. DOA estimations with ISM method and Channelization-ISM method.
Sensors 16 01031 g020
Figure 21. DOA estimations with TOPS method and Channelization-TOPS method.
Figure 21. DOA estimations with TOPS method and Channelization-TOPS method.
Sensors 16 01031 g021
Table 1. Hardware resources consumption.
Table 1. Hardware resources consumption.
Hardware ResourcesTotalConsumptionPercentage
occupied Slices49,20024564.99%
Slice Registers393,60010,5232.67%
Slice LUTs196,80066183.36%
DSP48E1s134447735.49%
RAMB36E1/FIFO36E1s704567.95%
Table 2. Mean deviation of different methods.
Table 2. Mean deviation of different methods.
MethodRMSE
TOPS Method (degree)0.6792
Channelization-TOPS Method (degree)0.6121
ISM Method (degree)3.2509
Channelization-ISM Method (degree)0.9103

Share and Cite

MDPI and ACS Style

Guo, R.; Zhang, Y.; Lin, Q.; Chen, Z. A Channelization-Based DOA Estimation Method for Wideband Signals. Sensors 2016, 16, 1031. https://doi.org/10.3390/s16071031

AMA Style

Guo R, Zhang Y, Lin Q, Chen Z. A Channelization-Based DOA Estimation Method for Wideband Signals. Sensors. 2016; 16(7):1031. https://doi.org/10.3390/s16071031

Chicago/Turabian Style

Guo, Rui, Yue Zhang, Qianqiang Lin, and Zengping Chen. 2016. "A Channelization-Based DOA Estimation Method for Wideband Signals" Sensors 16, no. 7: 1031. https://doi.org/10.3390/s16071031

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