Abstract

In a broadband adaptive array antenna with a space-time filter, a delay filter is required before digital beamforming when the Frost algorithm is used to obtain the weight vector. In this paper, we propose a Farrow structure instead of a direct form FIR structure to implement the time delay filter since it can satisfy the demand for real-time update and is very suited for the FPGA platform. Furthermore, a new off-line algorithm to calculate the Farrow filter coefficient is presented if the filter coefficient is symmetric. Finally, simulations are presented to illustrate that the design methodology for the Farrow filter is correct. The simulations also prove that the Frost broadband adaptive antenna effectively mitigates the interferences with a Farrow filter.

1. Introduction

In a broadband adaptive array antenna with space-time beamforming, each beamformer channel is a transversal filter. The Frost algorithm is a classic method to calculate the weight vector in space-time beamforming. In contrast to a narrowband system [1], Frost’s broadband adaptive array demands that the desired signal must impinge normally on the array face [2]. A delay filter is needed before beamforming if the desired signal does not satisfy the normal incident requirement. There is a large body of work focusing on how to obtain the weight vector with the Frost algorithm ([35] and the references therein). However, very limited research has been done on the design of the delay filter in a broadband Frost array.

The traditional delay filter takes a direct form structure. With this structure, the methodology for delay filter design can be divided into two categories: (a) time-domain methodology [6] and (b) frequency-domain methodology [7]. Both methodologies are online methodologies. Typically, the Frost broadband array delay filter [8] is designed using the time-domain methodology [6]. One downside of using the time-domain methodology is that different delay times result in different filter coefficients, which implies that we have to redesign the filter coefficients if the direction of arrival (DOA) of the desired signal changes because DOA decides the delay time. Often, the DOA of the desired signal changes very fast, especially for high-speed vehicles such as missiles. So, updating the filter in real time is almost impossible when the delay filter is implemented on the FPGA [9]. To overcome this disadvantage, we propose the Farrow structure to realize the delay filter for a Frost broadband adaptive antenna. The Farrow filter structure for a delay filter was first proposed by Farrow [10]. The design for a Farrow filter is done in an off-line manner; that is, the coefficients of the Farrow filter need to be designed only once and kept constant even if the delay requirement is changed. So, the Farrow filter is often used in multirate signal processing [11, 12], but a broadband adaptive array employing a Farrow filter has not been researched before.

The Farrow filter can be designed in an off-line manner. Convex optimization [13] and neural networks [14] can be used to compute the coefficients of each subfilter in the Farrow filter, but these methods have large implementation complexity. The weighted least squares (WLS) algorithm presented in [15] is another approach for Farrow filter design, but [15] did not show how to choose the weight value. An improved WLS algorithm is proposed in [16] that improves the performance under the assumption that the weight function is separable and piecewise constant. In [17], it is pointed out that the coefficients of the Farrow filter are symmetric, which allows the modification of the WLS algorithm to reduce the computational cost further.

This paper is organized as follows. Background information is given in this section. The Frost broadband adaptive array with space-time beamforming and the Farrow delay filter are discussed in Section 2. Farrow filter design and performance analysis are detailed in Section 3. The simulation work is presented in Section 4. Finally, the paper is concluded in Section 5.

2. Broadband Space-Time Adaptive Filtering

2.1. Space-Time Digital Beamforming

Different from the narrow band adaptive array, the space-time beamformer in a broadband adaptive array uses a transversal filter while the narrow band array uses only a complex multiplier at each channel. The space-time beamformer is shown in Figure 1 [2], where the array antenna has elements and the order of each transversal filter is .

Supposing that the received signal by the element is , the signal at each node could be described as

Clearly, is the sampling period. Then, the whole signal in the space-time filter is an dimension column vector; that is, where

Obviously, the corresponding weight for the whole beamformer is

And we see

So the output of the array is

The minimum variance distortionless response (MVDR) beamformer assumes the desired signal induced into the array from the normal direction, so that all the elements receive the same desired signal without any phase difference [2, 8]. Therefore, the optimal problem for space-time beamforming is as follows: where the covariance matrix is

The symbol indicates the mean value. and

The constraint requires the signal to impinge normally on the array, but this is not true generally. Then, a delay filter is required in serial with beamforming to adjust the input signal, so that the signal phase for each element is equal. The scheme of the broadband adaptive array is illustrated in Figure 2.

2.2. Frost Algorithm for Beamforming

To solve the optimum problem of (7), the Frost algorithm is employed which is a kind of LMS adaptive filter. The weight update equation of the Frost algorithm is where is the step size to assure convergence of the iteration only if

In the above equation, is the maximum eigenvalue of the covariance matrix . The output term is

The constant matrix is

3. Farrow Filter Design

3.1. Farrow Filter Model

The digital delay filter in Figure 2 is an FIR filter. If the order of the filter is , the coefficient vector can be expressed as

The delay includes -time sampling periods and a term which is shorter than one sampling period ; namely,

In the FPGA hardware platform, realizing delay is easy through clock control. Hence, the role of the delay filter is to deal with the delay . This kind of delay filter is also called variable fractional delay filters [1117]. We define the delay parameter as

Both time-domain and frequency-domain methods can design the coefficient vector of the delay filter. Since the delay parameter is decided by the DOA of the desired signal, it is necessary to redesign the coefficient vector at each time when DOA changes. If FPGA is used to implement the delay filter, the real-time requirements cannot be guaranteed. Hence, a Farrow filter is proposed to solve this problem. Although the Farrow filter is used in software radio widely [1114], it has not been considered for a broadband array before.

As for the Farrow digital delay filter, its each coefficient is regarded as -order polynomial of the delay parameter ; namely,

So the Farrow filter can be derived into subfilters as shown in Figure 3. So long as we adjust the delay parameter , the delay will be manipulated without change .

The Farrow filter has symmetric coefficients to assure the linear phase [17]. For the sake of convenience, we redefine the coefficient vector of the Farrow filter as

Obviously, according to the definition of the Farrow structure, we have

Meanwhile, assuming that the received signal of the array is a low pass signal and the delay filter is a low pass filter, so the pass band is , in which

3.2. Coefficient-Symmetric Farrow Filter Design with WLS Methodology

The system function of the ideal delay filter is expressed as

On the other hand, the frequency characteristic of the Farrow filter is

Now, we define a weighted error function as

The error function in the equation above is

At the same time, is a nonnegative function, which has these two characteristics:

With regard to the definitions above, Farrow filter design is a process of looking for suitable coefficient set , namely, , to minimize according to the band width and delay parameter .

Since the coefficients of the Farrow Filter are symmetric, we have the conclusion that

By substituting the conclusion into (23), we get where

The element of the matrix is

Substituting the result above into error function and taking the conjugate property of the weighted function into consideration, we get the weighted error function.

Compared to (24), the integration of the delay parameter is defined in . Hence, the process to solve can be converted into the process to solve the matrix which minimizes .

The weighted error function can be expressed as in which while

Assuming that the weight is known, all the equations above can be solved. Considering that , , , and are symmetric and positively definite, they can be factorized into , , , and , respectively, via the Cholesky factorization.

Finally, according to the Lagrangian multiplier method, taking the derivative of and letting it to be zero, we get

According to the definition (31), (32), and (33), can be determined.

3.3. Weight Determination

Generally, the weight is unknown to us. In order to decide the correct , iterative design is proposed. When beginning to design, we often let

Following the methodology in the previous section, we can obtain the designed coefficient . Then, we define the relative error as

If the required relative error cannot be satisfied, we define the weight function as segmentation:

For some interval with greater error, we set a larger weight and redesign . Then, the error in this interval could be reduced.

4. Simulation Result

4.1. Proposed WLS Farrow Filter Design Methodology

In our simulations, we use a linear array with 5 elements. The radio frequency and intermediate frequency of the received signal are 2 GHz and 400 MHz, respectively. Assuming that the sampling frequency is 1.25 GHz and the bandwidth is 300 MHz, the signal frequency is located in the interval [] since the first Nyquist domain is 0~312 MHz. In addition, the delay interval is []. We set the order of the Farrow delay filter as and .

We prove the validity of the proposed methodology for the Farrow filter design by evaluating the error , as defined in (25). We set the initial weight function as in (39) and compute the coefficient of the Farrow filter by the steps listed in Section 3.2. After obtaining the coefficient vector , it is easy to get the amplitude error as shown in Figure 4.

We focus on the interval where the desired signal exits. We notice that the large error is located in the intervals [] and []. If this error meets the requirement, we finish the design and take as the final result. And if not, we enlarge the weight function such as and we redesign the Farrow filter again. The new error is shown in Figure 5, where the error is less than about 6 dB (Figure 4) due to the weight function. In fact, the error is also impacted by the order of the Farrow filter which will be discussed in the next section.

4.2. Influence of the Filter Order

At the same simulation condition as the above, we only change the order of the filter and compare the error between the proposed Farrow and direct form structure in [68]. A comparison of the relative error , defined in (40), for different filter structures is listed in Table 1. The result in Table 1 illustrates that the direct form structure filter outperforms the Farrow filter if the order is small. However, when the order of the filter is large, the Farrow structure outperforms the direct form filters. Note that the filter with the direct form structure has to be redesigned online when the time delay changes. And hence it is not suitable for the condition where the DOA changes fast. On the other hand, the Farrow filter needs to be designed in an off-line manner only once and adjusts the input to change the delay time as shown in Figure 3.

Figure 6 illustrates the relative error with different and . It shows that as the order of the Farrow filter increases, the relative error decreases rapidly. But once the order is large enough, the performance keeps constant.

4.3. Computational Complexity

The WLS algorithm to design the Farrow filter in [17] does not consider symmetric coefficient property. Under the same conditions, we compare the computational complexity between the proposed methodology and the one in [17]. The order of the Farrow delay filter is set with and . The running time of the MATLAB program and the error are listed in Table 2, where the CPU frequency is 2.4 GHz. We can see that the new approach could reduce the computational complexity about 50% with the similar relative error as compared to the reference methodology.

4.4. Beamforming Based on the Farrow Filter

To investigate the function of the Farrow filter in a space-time two-dimensional filter, we do the simulation about the array pattern. Suppose that the desired signal induces on the array with and the DOA of one strong interference is θ1 = 60°. The radio frequencies of the two signals are both 2 GHz and the bandwidths are both 400 MHz.

Without a Farrow filter, the pattern is shown in Figure 7 by a dashed red line, where there are two nulls at and , respectively. That means that the desired signal and interference are both suppressed. The other three lines in Figure 7 are the pattern when the delay filter is used through the three methodologies such as time domain [6, 8], frequency domain [7], and Farrow. They are very similar because the delay filters play the same functions. The interference is mitigated since the null still appears at . The main lobe points to where the desired signal is impinging due to the role of calibration by the delay filter. Clearly, the result means that the delay filter could keep the desired signal to satisfy the requirement of the Frost algorithm. And it also proves that our proposed methodology for Farrow filter design is right.

If we zoom in the null in Figure 8, it can be seen that the depth is different. This is because the error of the delay filter is different. With the same filter order 56, the result is consistent with the result shown in Table 1 and that less error means deeper null.

5. Conclusion

The Farrow delay filter is taken in the broadband adaptive array instead of the delay filter with a direct form structure. A new methodology is presented to design the Farrow filter according to the symmetric coefficient. The computational cost is lower than that of the classic WLS algorithm. Although the Farrow filter is suited to the FPGA platform, the coefficient computation is heavily complicated. How to choose the order of the Farrow filter and error requirement is an important work in the future.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no financial and personal relationships with other people or organizations that can inappropriately influence their work; there are no professional or other personal interests of any nature or kind in any product, service, and company that could be construed as influencing the position presented in the paper.

Acknowledgments

The authors would like to acknowledge the support of innovation funds from China Academy of Space Technology (Project no. CAST2016021) and funds from Chongqing University (Project no. 2018CDGFTX0014).