Table of Contents Author Guidelines Submit a Manuscript
International Journal of Digital Multimedia Broadcasting
Volume 2010 (2010), Article ID 194306, 6 pages
Research Article

Frequency-Based Optimization Design for Fractional Delay FIR Filters with Software-Defined Radio Applications

1Electronics Department, Technological Institute of Celaya, Celaya, GTO 38010, Mexico
2Electronics Department, National Institute of Astrophysics Optics and Electronics, Tonantzintla, PUE 72840, Mexico

Received 8 November 2009; Accepted 1 January 2010

Academic Editor: Massimiliano Laddomada

Copyright © 2010 Javier Díaz-Carmona et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


A frequency-designed fractional delay FIR structure, which is suitable for software radio applications, is presented. The design method is based on frequency optimization of a combination of modified Farrow and mutirate structures. As a result the optimization frequency range is made only in half of desired total bandwidth. According to the obtained results the proposed fractional delay structure allows online desired fractional delay update, with a high fractional delay value resolution.

1. Introduction

Software-defined radio represents a major change in the design paradigm for radios, [1], in which most of the functionality is made through programmable signal processing devices, giving the radio the ability to change its operating parameters to accommodate new features and capabilities [2].

A software-defined radio platform is designed to make mobile systems more flexible with respect to the bandwidth requirements of different mobile standards. This flexibility is achieved by performing channel selection in the digital domain through sample rate conversion (SRC) with programmable digital filters. Fractional Delay (FD) filters are key components used to perform nonrational SRC [35].

Additionally, Software radio systems employ direct conversion receivers with asynchronous sampling such that the actual sampling instants are not synchronized with the incoming symbols. In order to evaluate the received symbols a digital symbol synchronization is implemented through FD filter structures [6, 7].

One of the key requirements for FD in software radio applications is to have the flexibility to change among different communication protocols and to be able to perform fractional delay value update on line, known as variable fractional delay (VFD) filters [8]. Other important FD characteristics are wide bandwidth, high fractional delay resolution value, and a small number of arithmetic operations per output sample.

There are several FD design methods [9]; among them the use of a polynomial approach allows online desired fractional delay value update using a Farrow structure [10] or a modified Farrow structure [11, 12]. Both structures are composed of parallel FIR filters , each one with length , where is the chosen polynomial order, as it is shown in Figure 1. In a modified Farrow structure , where is the required fractional delay value, , and are linear phase filters (symmetrical coefficient values). In the original Farrow structure and parallel filters are not symmetrical.

Figure 1: Farrow structure.

There are two main FD design approaches based on the polynomial approach. The first one is completely time domain design based on either Lagrange interpolation [13] or B spline functions [14, 15]. The implementation of this design approach is made through original Farrow structure having as main advantage that filter coefficients are obtained by closed form expressions. The disadvantage is the small flexibility available for this approach to meet FD filter frequency domain specifications. This is because there is only one design parameter, that is, the polynomial order . Most of the recently reported FD design methods belong to this time domain approach such as [16, 17].

The second FD design approach is made in frequency domain using optimization techniques for coefficients computing. The main advantage of this design approach is an improved control on frequency specifications. This is because three design parameters are available: polynomial order , filter length , and desired frequency passband . The disadvantage of this approach is the need of an optimization method for filter coefficient computing. Several design methods have been proposed such as, for example, [18], where the FD filter is implemented in a modified Farrow structure and a Taylor approximation is achieved. Similarly in [1921] the implementation is made using an original Farrow structure and a weighted least square optimization is accomplished.

The use of frequency domain design methods for FD filters with a wide bandwidth requires that an optimization method be applied over a large frequency range. On the other hand large filter length and polynomial order are obtained when high fractional delay resolution is required. Hence the design process requires high coefficients computing time and a high number of arithmetic operations per output sample in the resulting FD filter implementation.

This paper describes the use of a multirate structure in a frequency design approach in order to reduce the optimization workload in coefficients computing for FD filters with a wide bandwidth, high fractional delay resolution, and online fractional delay value update capability. In this way a flexible frequency design method with a reduced optimization workload as well as a resulting structure with a reduced number of arithmetic operations per output sample is obtained.

The used frequency design method is the modified Farrow structure [18], where each parallel filter is designed as a minimum least square approximation of an order differentiator. In the same way it is possible to extend such proposal through other optimization frequency design methods.

Section 2 describes in a general way the frequency design method basis. The multirate structure is given in Section 3. In Section 4 the proposed design method is shown, which is illustrated through one design example. Conclusions are presented in Section 5.

2. Frequency Design

The frequency design method in [18] is based on the following properties of the parallel digital filters .

(1)FIR filters , , in original Farrow structure are an order Taylor approximation to the continuous-time interpolated input signal.(2)In the modified Farrow structure the FIR filters are linear phase type II for even and type IV for odd.

Each filter approximates in magnitude the function , where is a constant. The ideal frequency response of an order differentiator is ; hence the ideal response of each filter in the Farrow structure is an order differentiator.

In the same way it is possible to approximate the input signal through Taylor series in a modified Farrow structure. The order differential approximation to the continuous-time interpolated input signal is done through the branch filter , with a frequency response given as The input design parameters are the filter length , the polynomial order , and the desired passband frequency .

The coefficients of the FIR filters are computed in such a way that the following error function is minimized in a least square sense in the frequency range [0, ]: where Hence the objective function is given as From this equation it can be observed that the design of a wide bandwidth FD filter requires an extensive computing workload. For high fractional delay resolution FD filters high precise differentiator approximations are required; this implies high branch filters length, , and high polynomial order, . Hence an FD filter structure with high number of arithmetic operations per output sample is obtained.

The arithmetic complexity of the resulting implemented structure is an important factor to be considered. The comparative parameters are the following:

(1)number of multipliers per output sample (MPS),(2)number of additions per output sample (APS).

In the modified Farrow structure the MPS1 and APS1 are given as

3. Multirate Structure

The multirate structure in [22] is proposed for designing FD filters in time domain. The input signal bandwidth is reduced by incrementing the sampling frequency. In this way Lagrange interpolation is used in filter coefficients computing for an FD filter with a wide bandwidth.

The multirate structure shown in Figure 2 is composed of three sections. The first one is an upsampler and a half band image suppressor filter for incrementing twice the input sampling frequency. The second section is the FD filter , which is designed in time domain through Lagrange interpolation [11]. Since the signal processing frequency of filter is two times more than the input sampling frequency, such filter can be designed to meet only half of the required bandwidth. Last section deals with a downsampler for decreasing the sampling frequency to its original value. The upsamplig process is made through insertion of one zero between every two input samples. Hence for each output sample only half of the FIR filter coefficients are used. This means that in one time instant the input samples are processed through the even coefficients and next time instant through the odd coefficients of the filters and . According to this technique and using multirate processing noble identities [23], such processing can be represented as shown in Figure 3, where filters and are the first and second polyphase components of the half band filter . In the same way and are the first and second polyphase components of fractional delay filter . As can be seen the input sampling frequency is the same for all filters in the resulting structure.

Figure 2: Multirate structure for FD filter.
Figure 3: Resulting structure for FD filter.

4. Proposed Design Method

The proposed method for FD filter design with a wide bandwidth and high fractional delay resolution is based on a frequency domain optimization approach, described in the second section, applied to the FD multirate structure, described in last the section.

As mentioned before the maximum frequency of the FD filter in the multirate structure is half of the desired bandwidth. In this way the frequency optimization is made only on the half of the required passband. That means that the upper frequency limit in (4) is . This optimization frequency range decrease allows an abrupt coefficient computing time reduction for the wide bandwidth FD filters and the resulting structure requires filters with smaller length .

In Figure 4 is shown the initial structure, where a double fractional delay value is considered in the update parameter , as a result of the doubled processing sampling frequency. The resulting structure after applying noble identities is shown in Figure 5, where the filters and are the first and second polyphase components of , respectively.

Figure 4: Initial structure of the proposed method.
Figure 5: Resulting proposed structure.

The filter plays a key role in resulting bandwidth and fractional delay resolution of FD filter. The higher the stopband attenuation of the filter , the higher the resulting fractional delay resolution. Similarly the smaller the transition band of , the higher the resulting bandwidth of the FD filter. Both conditions imply the use of a high-order filter.

In order to reduce the total number of arithmetic operations per output sample the filter is designed as a half hand FIR filter. In this way the total number of MPS and APS is given as where is the number of coefficients of .

5. Obtained Results

The design method was implemented in MATLAB. An illustrative design example is presented with an FD filter bandwidth of 0.9 and a fractional delay resolution of 1/10000.

The FD filter design using WLS method [20] results in an implementation processing arithmetic of and with design parameters and and weighting functions given by

In the same way the resulting processing arithmetic for the implementation of the FD filter example based on variable FD design method of [21] is and with design parameters and and next weighting function:

The direct use of the frequency domain method [18] with design parameters of and results in a total number of and .

For the proposed design method an interpolator filter with 241 coefficients was used, designed with a Dolph-Chevishev window having a stopband attenuation of 140 dBs. The design parameters are and with a resulting processing arithmetic of and .

In the proposed method the frequency optimization is applied up to and in the direct method to , as is depicted in Figure 6, where the first seven differentiator approximations are shown for both methods. This half frequency range optimization implies a notable computing workload reduction. The MATLAB computing time in a PC running at 2 GHz for the proposed method is 6.95 seconds and for direct frequency optimization of [18] is 322 seconds.

Figure 6: Differentiators ideal responses (dot line) and approximations (solid lines). (a) Proposed method , , and . (b) Direct all band optimization , , and .

The all band magnitude responses and group delays for fractional delay values range from 0.0080 to 0.0090 using the direct frequency FD design method and the proposed method results are shown in Figures 7 and 8, respectively.

Figure 7: All band magnitude responses for (a) proposed method and (b) direct all band optimization.
Figure 8: All band fractional delay responses for (a) proposed method and (b) direct all band optimization.

According to obtained results the proposed method has smaller number of operations per output sample. In order to compare the achieved proposed method approximation with the one obtained with existing methods, the frequency domain error , the maximum absolute error , and the root mean square error are defined, like in [21], as where and are the frequency response of the designed and ideal FD filters, respectively, and is passband frequency of the FD filter.

The obtained maximum absolute error and the root mean square error are presented in Table 1; for comparison purpose the results obtained by using the approaches in [1821, 24] are also presented.

Table 1: Comparison of approximation errors for several methods.

6. Conclusions

A frequency optimization design approach for wide bandwidth and high fractional value resolution FD filters has been proposed. These specifications coupled with the capability of updating the fractional delay value in real-time make the resulting structure suitable to perform important physical layer functions for software-defined radio applications, such as digital symbol synchronization and sample rate conversion.

The obtained results show that the design method notably reduces the coefficients computing workload. The resulting structure allows fractional delay values of 1/10000 of sample and a bandwidth of , with a reduced number of arithmetic operations per output sample. The described method is based on a least mean square frequency optimization for coefficients computation. In a future work we will consider the use of other optimization methods in the same proposed approach.


  1. J. Mitola, “The software radio architecture,” IEEE Communications Magazine, vol. 33, no. 5, pp. 26–38, 1995. View at Google Scholar
  2. J. Vankka, Digital Synthesizers and Transmitters for Software Radio, Springer, Dordrecht, The Netherlands, 2005.
  3. T. Hentschel and G. Fettweis, “Sample rate conversion for software radio,” IEEE Communications Magazine, vol. 38, no. 8, pp. 142–150, 2000. View at Google Scholar
  4. K. S. Yeung and S. C. Chan, “The design and multiplier-less realization of software radio receivers with reduced system delay,” IEEE Transactions on Circuits and Systems I, vol. 51, no. 12, pp. 2444–2459, 2004. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  5. V. Lehtinen, D. Babic, and M. Renfors, “On impulse response symmetry of farrow interpolators in rational sample rate conversion,” in Proceedings of the International Symposium on Control, Communications and Signal Processing, (ISCCSP '04), pp. 693–696, 2004.
  6. M. Makundi, T. Laakso, and A. Hjorungnes, “Iterative symbol synchronization using variable IIR and FIR fractional-delay filters with arbitrary oversampling ratios,” in Proceedings of the 1st International Symposium on Wireless Communication Systems, pp. 105–109, September 2004.
  7. S. Pietro, G. Paolo, and F. Lorenzo, “An iterative algorithm for joint symbol timming recovery and equalization of short bursts,” Journal of Communications, vol. 3, no. 4, pp. 34–40, 2008. View at Google Scholar
  8. D. Tian-Bo, “Design and parallel implementation of FIR digital filters with simultaneous variable magnitude and non-integer phase-delay,” IEEE Transactions on Circuits and Systems II, vol. 50, no. 5, pp. 243–250, 2003. View at Google Scholar
  9. T. I. Laakson, V. Valimaki, M. Karjalainen, and U. K. Laine, “Splitting the unit delay: tools for fractional filter design,” IEEE Signal Processing Magazine, vol. 13, no. 1, pp. 30–60, 1996. View at Google Scholar
  10. C. W. Farrow, “A continuously variable digital delay element,” in Proceedings of the IEEE International Symposium on Circuits and Systems, vol. 3, pp. 2641–2645, June 1988.
  11. J. Vesma and T. Saramaki, “Optimization and efficient implementation of FIR filters with adjustable fractional delay,” in Proceedings of the IEEE International Symposium on Circuits and Systems, pp. 2256–2259, June 1997.
  12. J. Vesma and T. Saramaki, “Design and properties of polynomial-based fractional delay filters,” in Proceedings of the IEEE International Symposium on Circuits and Systems, pp. 104–107, May 2000.
  13. G.-S. Liu, C.-H. Wei et al., “A new variable fractional sample delay filter with nonlinear interpolation,” IEEE Transactions on Circuits and Systems II, vol. 39, no. 2, pp. 123–126, 1992. View at Google Scholar
  14. J. Vesma, Timing adjustment in digital receivers using interpolation, M.S. thesis, Department of Information Technology, Tampere University of Technology, Tampere, Finland, 1995.
  15. S. Samadi, M. O. Ahmad, and M. Swamy, “Characterization of B-spline digital filters,” IEEE Transactions on Circuits and Systems I, vol. 51, no. 4, pp. 808–816, 2004. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  16. C. Cagatay, “An efficient filter structure for lagrange interpolation,” IEEE Signal Processing Letters, vol. 14, no. 1, pp. 17–19, 2007. View at Google Scholar
  17. A. Haghparast and V. Välimäki, “A computationally efficient coefficient update technique for lagrange fractional delay filters,” in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '08), pp. 3737–3740, March 2008. View at Publisher · View at Google Scholar
  18. J. Vesma, R. Hamila, T. Saramaki, and M. Renfors, “Design of polynomial-based interpolation filters based on Taylor series,” in Proceedings of the 10th European Signal Processing Conference, pp. 283–286, September 1998.
  19. A. Tarczynski, G. D. Cain, E. Hermanovicz, and M. Rojewski, “WLS design of variable frequency response FIR response,” in Proceedings of the IEEE International Symposium on Circuits and Systems, pp. 2244–2247, June 1997.
  20. T.-B. Deng, “Discretization-free design of variable fractional-delay FIR digital filters,” IEEE Transactions on Circuits and Systems II, vol. 48, no. 6, pp. 637–644, 2001. View at Publisher · View at Google Scholar · View at Scopus
  21. H. Zhao and J. Yu, “A simple and efficient design of variable fractional delay FIR filters,” IEEE Transactions on Circuits and Systems II, vol. 53, no. 2, pp. 157–160, 2006. View at Publisher · View at Google Scholar · View at Scopus
  22. N. Murphy, A. Krukowski, and I. Kale, “Implementation of a wide-band integer and fractional delay element,” Electronics Letters, vol. 30, no. 20, pp. 1658–1659, 1994. View at Google Scholar
  23. P. Vaidyanathan, Multirate Systems and Filter Banks, Signal Processing Series, Prentice Hall, Upper Saddle River, NJ, USA, 1993.
  24. W.-S. Lu and T.-B. Deng, “An improved weighted least-squares design for variable fractional delay fir filters,” IEEE Transactions on Circuits and Systems II, vol. 46, no. 8, pp. 1035–1040, 1999. View at Google Scholar · View at Scopus