A Channelization-Based DOA Estimation Method for Wideband Signals

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.


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. 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.

The Wideband Signal Model
Assuming a uniform linear array (ULA) of K omnidirectional sensors receives PpP ď Kq wideband signals s p ptqpp " 1, 2,¨¨¨Pq from unknown directions tθ 1 , θ 2 ,¨¨¨, θ P u. The kth sensor output x k ptq is: x k ptq " P ÿ p"1 s p pt´τ k pθ p qq`n k ptq k " 1, 2,¨¨¨K (1) where n k ptq is the additive white Gaussian noise of kth sensor which is assumed to be statistically independent with signal sources. τ k pθ p q is the propagation delay of the kth sensor: 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: 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: where: Xpω i q " rX 1 pω i q, X 2 pω i q,¨¨¨, X K pω i qs T Spω i q " rS 1 pω i q, S 2 pω i q,¨¨¨, S P pω i qs T [¨] T denotes the transpose operator. Apω i q is the KˆP steering matrix. Its columns are the Kˆ1 array manifolds apω i , θ p q: Apω i q " rapω i , θ 1 q, apω i , θ 2 q,¨¨¨, apω i , θ P qs (7) apω i , θ p q " " e´j ω i τ 1 pθ p q , e´j ω i τ 2 pθ p q ,¨¨¨, e´j ω i τ K pθ p q ı T (8)

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.

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 1. The array output is channelized into J frequency channels. The carrier frequency of the i-th channel is at: h n is a low-pass filter. The output of i-th sub-channel can be expressed as: The order of prototype filter is D and 1 D L J    , where L is a positive integer. By polyphase decomposition with the prototype filter, Equation (10) can be rewritten as:

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 pnq represents the digitization of kth sensor output x k ptq 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: h LP pnq is a low-pass filter.

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 1. The array output is channelized into J frequency channels. The carrier frequency of the i-th channel is at: The output of i-th sub-channel can be expressed as: The order of prototype filter is D and 1 D L J    , where L is a positive integer. By polyphase decomposition with the prototype filter, Equation (10) can be rewritten as:  The output of i-th sub-channel can be expressed as: 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: Substituting Equation (9) into Equation (12), we can rewrite y i k pmq as: where˚represents convolution operation. Replacing x p k pmq˚h p pmq asx p k pmq, 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 J is moved to the left of the filtering operation. h 0 pmq to h J´1 pmq represent the polyphase components of the low-pass filter h LP pnq. This requires fewer components and is a more realistic architecture. . Then: Substituting Equation (9) into Equation (12), we can rewrite ( ) i k y m as: where  represents convolution operation. Replacing x m , then: 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]. The output of digital channelization receiver y i k pmq 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 pnq and M that of y i k pmq, 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].

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 pnq consists of a wideband signal s k pnq and white Gaussian noise n k pnq " Np0, σ 2 n q. s k pnq has a power of A 2 and the power spreads on SpS ă Jq sub-channels equally. The input SNR is: The output y i k pmq is composed of the response due to the noise y i kn pmq and the signal y i ks pmq, respectively. The noise response is a complex filtered baseband noise with an autocorrelation function γ y kn [29]: Only the S sub-channels within the wideband signal bandwidth are expected in further processing. The SNR of these sub-channels are: The SNR is improved approximately J{SpJ{S ą 1q 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.

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: where r¨s H denotes the complex conjugate transpose ypω i , mq " ry i 1 pmq, y i 2 pmq,¨¨¨, y i K pmqs T .
Performing an eigendecomposition onRpω i q, the eigenvalues λ 1 , λ 2 ,¨¨¨, λ K in descending order and their associated eigenvectors u 1 , u 2 ,¨¨¨, u K can be obtained. Define matrices Fpω i q " ru 1 ,¨¨¨, u P s and Upω i q " ru P`1 ,¨¨¨, u K s, the columns of Fpω i q and Upω i q 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 θ 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: 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.

Channelized Based TOPS Method
Following the standard TOPS approach introduced in [7], we first perform an eigendecomposition onRpω i q and obtain the noise subspace Upω i q. Select one output sub-channel as reference to obtain the reference signal subspace Fpω 0 q.
Define a diagonal unitary transformation matrix whose size is kˆk and diagonal elements are rTpω i , θ i qs k,k : rTpω i , θ i qs 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 pφq: where ∆ω i " ω i´ω0 , φ is a hypothetic arrival angle of signal: Dpφq is a full rank matrix if φ is not equal to the arrival angle of signal otherwise Dpφq 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 σ min pφq is the smallest singular value of Dpφq.

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{SpJ{S ą 1q 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.

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: RMSE " 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 Figures 4 and 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, Figures 6 and 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.  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 7 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   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 7 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 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. 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 8 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 9 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  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 8 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 9 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.   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.   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.

Experiments
In this section, we present experiments carried out in a microwave anechoic chamber to further demonstrate the performance of the proposed algorithm. Figure 10 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 11. Figure 12 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 13, 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 14 is based on the latest generation of high-performance multicore fixed and floating-point

Experiments
In this section, we present experiments carried out in a microwave anechoic chamber to further demonstrate the performance of the proposed algorithm. 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.

Experiments
In this section, we present experiments carried out in a microwave anechoic chamber to further demonstrate the performance of the proposed algorithm. Figure 10 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 11. Figure 12 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 13, 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 14 is based on the latest generation of high-performance multicore fixed and floating-point 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. 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.    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.    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.

Experiment Carried Out in Microwave Anechoic Chamber
The experiment is carried out in a microwave anechoic chamber. As shown in Figure 15, 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 16. 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.

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 f 0 = 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.

Experiment Carried Out in Microwave Anechoic Chamber
The experiment is carried out in a microwave anechoic chamber. As shown in Figure 15, 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 16. 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 17 shows the amplitude of received signals in time domain. Eight curves with different colors represent the different sensor outputs of the ULA. Figure 18 depicts the received signals in the frequency domain. The spatial spectrum estimations of the received signals with different methods are plotted in Figure 19 and Figure 20. 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.     Figure 17 shows the amplitude of received signals in time domain. Eight curves with different colors represent the different sensor outputs of the ULA. Figure 18 depicts the received signals in the frequency domain. The spatial spectrum estimations of the received signals with different methods are plotted in Figure 19 and Figure 20. 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.

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

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

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.