Open access

# Fractional Delay Digital Filters

Written By

Javier Diaz-Carmona and Gordana Jovanovic Dolecek

Submitted: 25 March 2011 Published: 09 September 2011

DOI: 10.5772/31051

From the Edited Volume

## Applications of MATLAB in Science and Engineering

Chapter metrics overview

View Full Metrics

## 1. Introduction

The chapter goal is focused to introduce the concept of fractional delay filters (FDF), as well as a concise description of most of the existing design techniques. For this purpose, several illustrative examples are presented, where each design method is implemented by MATLAB programs.

A fractional delay filter is a filter of digital type having as main function to delay the processed input signal a fractional of the sampling period time. There are several applications where such signal delay value is required, examples of such systems are: timing adjustment in all-digital receivers (symbol synchronization), conversion between arbitrary sampling frequencies, echo cancellation, speech coding and synthesis, musical instruments modelling etc. (Laakson et al., 1996).

In order to achieve the fractional delay filter function, two main frequency-domain specifications must be met by the filter. The filter magnitude frequency response must have an all-pass behaviour in a wide frequency range, as well as its phase frequency response must be linear with a fixed fractional slope through the bandwidth.

Several FIR design methods have been reported during the last two decades. There are two main design approaches: time-domain and frequency-domain design methods. In first one, the fractional delay filter coefficients are easily obtained through classical mathematical interpolation formulas, but there is a small flexibility to meet frequency-domain specifications. On the other hand, the frequency-domain methods are based on frequency optimization process, and a more frequency specification control is available. One important result of frequency-domain methods is a highly efficient implementation structure called Farrow structure, which allows online fractional value update.

The chapter is organized as follows. Next section gives the formal definition of fractional delay filter. In the third section, some design methods are briefly described. Two efficient implementation structures for wideband fractional delay filter, as well as description of recently reported design methods for such structures, are illustrated in fourth section. MATLAB designed examples and concluding remarks are presented in fifth and sixth sections, respectively.

## 2. Fractional delay filter definition

The continuous-time output signal ya(t) of a general signal delay system is defined by:

ya(t)=x(ttl)E1

where x(t) is the continuous-time input signal and tl the obtained time delay value. In a discrete-time system, the input-output relationship of a signal delay system is expressed as:

y(lT)=x(nTDT)E2

where the delay value is given by DT, y(lT) and x(nT) are the discrete-time versions of output and input signals, respectively, and T is the sampling period time.

A signal delay value equal to a multiple of the sampling period, D as an integer N, can be easily implemented in a discrete-time system by memory elements storing the signal value for a time of NT:

y(lT)=x(nTNT)E3

In this case, the signal delay value is limited to be only N time the sampling period, tl=NT. For instance in telephone quality signals, with a sampling frequency of 8 KHz, only delays values multiple of 125μseconds are allowed.

Let us introduce the FDF function using time-domain signals sketched in Fig 1. The FDF output y(lT), squared samples, is obtained a delay time tl after input x(nl), with a delay value μlT given as a fraction of the sampling period time, 0l<1. As shown in Fig. 1, the fractional delay value μl may be variable; this way, it can be changed at any desired time.

The fundamental design problem of a FDF is to obtain the FDF unit impulse response hFD(n,μ), in such a way that the obtained output value y(lT) = ya(DT) be as close as possible to ya(tl) for 0l <1. The simplified block diagram for a FDF is shown in Fig. 2, which output for a no causal FIR FDF filter is given by the discrete-time convolution:

y(lT)=k=NFD/2NFD/21x(nlk)hFD(k,μl)E4

where NFD is the even length of the FDF. The system function H(z) of the FDF can be expressed as:

H(z)=zDE5

where the delay value is given as: D = Dfix+μl, Dfix is a fixed delay value and μl is the desired fractional delay value. As a consequence, the ideal frequency response of a FDF Hid(ω,μl) is:

Hid(ω,μl)=ej(Dfix+μl)ωE6

Hence the ideal FDF frequency response has an all-band unity magnitude response:

|Hid(ω,μl)|=1,|ω|πE7

and a linear frequency phase response with a constant phase delay given, respectively, by:

θid(ω,μl)=DωE8
τpid(ω,μl)=DE9

The main goal of all existing FDF design methods, based on a frequency design approach, is to obtain the FDF filter coefficients through approximating this ideal frequency performance.

Applying inverse discrete Fourier transform to the ideal FDF frequency response, the ideal FDF filter unit impulse response hid(n,μ) is obtained as:

hid(n,μ)=sin[π(nD)]π(nD)=sinc(nD)E10

Given a desired factional delay value, the FDF coefficients can be easily obtained with this infinite length delayed sinc function. Due to this infinite length, it is evident that an FIR FDF will be always an approximation to the ideal case.

As an illustrative example, the ideal FDF unit impulse responses for two delay values D= 3.0 (Dfix=3.0 and μ = 0) and D=3.65 (Dfix=3.0 and μ = 0.65) are shown in Fig. 3 and 4, respectively. The unit impulse responses were obtained using MATLAB function sinc. The FDF unit impulse responses are shown as solid lines, and the delayed sinc function as dot line. In the first case, only one three-sample delay is needed, which can be easily implemented with memory components as described above. However, the FDF unit impulse response for the second case has an infinite number of nonzero coefficients (IIR) and it is a no causal sequence, which makes it impractical for implementing in real-time applications.

## 3. FDF Design methods

The existing design methods for FIR FDF use a large range of strategies to approximate as close as possible the ideal FDF unit impulse response hid(n,μ). It is possible to highlight three main strategies:

• Magnitude frequency response approximation: The FDF unit impulse response is obtained such that its frequency magnitude response is as close as possible to the ideal FDF one, accordingly to some defined error criterion.

• Interpolation design method: The design approach is based on computing FDF coefficients through classical mathematical interpolation methods, such as Lagrange or B-spline. The design is a completely time-domain approach.

• Hybrid analogue-digital model approach: The FDF design is accomplished through the use of an analogue-digital model. The design methods using this strategy are based on a frequency-domain approach.

A concise description of each one of these strategies is presented in the following.

### 3.1. Magnitude frequency response approximation

The design method goal is to obtain the FDF unit impulse response hFD(n,μ) based on comparing its magnitude frequency response with the ideal one. The frequency response of the designed FDF with even-length NFD is given by:

HFD(ω,μ)=k=NFD/2+1NFD/2hFD(k,μ)ejωkE11

One of the criterions used for the magnitude frequency response comparison is the least squares magnitude error defined as:

e2(ω)=1π0ωp|HDF(ω,μ)Hid(ω,μ)|2dωE12

The error function e2(ω) is minimized by truncating the ideal unit impulse response to NFD samples, which can be interpreted as applying a delayed M-length window w(n) to the ideal IIR FDF unit impulse response:

hFD(n,μ)=hid(nD,μ)w(nD)E13

where ω(n) is equal to unity in the interval 0≤nNFD-1 and zero otherwise.

The windowing process on the ideal unit impulse response causes not-desired effects on the FDF frequency response, in particular the Gibbs phenomenon for rectangular window (Proakis & Manolakis, 1995).

In general, the performance of a FDF obtained by truncating the sinc function is usually not acceptable in practice. As a design example, the FDF frequency magnitude and phase responses for D=3.65, using a rectangular window with NFD=50, are shown in Fig 5. We can see that the obtained FDF bandwidth is less than 0.9π and although the IIR sinc function has been truncated up to 50 taps, neither its frequency magnitude nor its phase response are constant.

The windowed unit impulse response hFD(n,μ) has a low-pass frequency response, in this way it can be modified to approximate only a desired pass-band interval (0,απ) as follows:

hFD(n,μ)={αsinc[α(nD)]for0nNFD10otherwiseE14

The magnitude and phase responses of a FDF with NFD= 8 and α=0.5 are shown in Fig. 6, which were obtained using MATLAB. The phase delay range is from D=3.0 to 3.5 samples with an increment of 0.1. More constant phase delay responses and narrower bandwidth is achieved.

In principle, window-based design is fast and easy. However, in practical applications it is difficult to meet a desired magnitude and phase specifications by adjusting window parameters. In order to meet a variable fractional delay specification, a real-time coefficient update method is required. This can be achieved storing the window values in memory and computing the values of the sinc function on line, but this would require large memory size for fine fractional delay resolution (Vesma, 1999).

The smallest least squares error can be achieved by defining its response only in a desired frequency band and by leaving the rest as a “don’t care” band. This can be done using a frequency-domain weighting as follows (Laakson et al., 1996):

e3(ω)=1π0ωpW(ω)|HDF(ω,μ)Hid(ω,μ)|2E15

where ωp is the desired pass-band frequency and W(ω) represents the weighting frequency function, which defines the corresponding weight to each band. In this way, the error is defined only in the FDF pass-band, hence the optimization process is applied in this particular frequency range.

In Fig. 7 are shown the FDF frequency responses designed with this method using W(ω)=1, ΝFD = 8 and α =0.5. We can see a notable improvement in the resulting FDF bandwidth compared with the one obtained using the least square method, Fig. 6.

There is another design method based on the magnitude frequency response approach, which computes the FDF coefficients by minimizing the error function:

e4(ω)=max0ωωp|HFD(ω,μ)Hid(ω,μ)|E16

The solution to this optimization problem is given by the minimax method proposed by (Oetken, 1979). The obtained FDF has an equiripple pass-band magnitude response. As an illustrative example, the frequency response of an FDF designed through this minimax method is shown in Fig. 8, where NFD=20 and ωp=0.9π.

### 3.2. Interpolation design approach

Instead of minimizing an error function, the FDF coefficients are computed from making the error function maximally-flat at ω=0. This means that the derivatives of an error function are equal to zero at this frequency point:

nec(ω)ωn|ω=0=0,n=0,1,2,....NFD1E17

the complex error function is defined as:

ec(ω)=HFD(ω,μl)Hid(ω,μl)E18

where HFD(ω,μl) is the designed FDF frequency response, and Hid(ω,μl) is the ideal FDF frequency response, given by equation (Eq. 6). The solution of this approximation is the classical Lagrange interpolation formula, where the FDF coefficients are computed with the closed form equation:

hL(n)=k=0knNFDDknkn=0,1,2,....NFDE19

where NFD is the FDF length and the desired delayD=NFD/2+μl. We can note that the filter length is the unique design parameter for this method.

The FDF frequency responses, designed with Lagrange interpolation, with a length of 10 are shown in Fig. 9. As expected, a flat magnitude response at low frequencies is presented; a narrow bandwidth is also obtained.

The use of this design method has three main advantages (Laakson et al., 1994): 1) the ease to compute the FDF coefficients from one closed form equation, 2) the FDF magnitude frequency response at low frequencies is completely flat, 3) a FDF with polynomial-defined coefficients allows the use of an efficient implementation structure called Farrow structure, which will be described in section 3.3.

On the other hand, there are some disadvantages to be taken into account when a Lagrange interpolation is used in FDF design: 1) the achieved bandwidth is narrow, 2) the design is made in time-domain and then any frequency information of the processed signal is not taken into account; this is a big problem because the time-domain characteristics of the signals are not usually known, and what is known is their frequency band, 3) if the polynomial order is NFD; then the FDF length will be NFD, 4) since only one design parameter is used, the design control of FDF specifications in frequency-domain is limited.

The use of Lagrange interpolation for FDF design is proposed in (Ging-Shing & Che-Ho, 1990, 1992), where closed form equations are presented for coefficients computing of the desired FDF filter. A combination of a multirate structure and a Lagrange-designed FDF is described in (Murphy et al., 1994), where an improved bandwidth is achieved.

The interpolation design approach is not limited only to Lagrange interpolation; some design methods using spline and parabolic interpolations were reported in (Vesma, 1995) and (Erup et al., 1993), respectively.

### 3.3. Hybrid analogue-digital model approach

In this approach, the FDF design methods are based on the hybrid analogue-digital model proposed by (Ramstad, 1984), which is shown in Fig. 10. The fractional delay of the digital signal x(n) is made in the analogue domain through a re-sampling process at the desired time delay tl. Hence a digital to analogue converter is taken into account in the model, where a reconstruction analog filter ha(t) is used.

An important result of this modelling is the relationship between the analogue reconstruction filer ha(t) and the discrete-time FDF unit impulse response hFD(n,μ), which is given by:

hFD(n,μ)=ha((n+μl)T)E20

where n=-NFD/2,-NFD/2+1,…., NFD/2-1, and T is the signal sampling frequency. The model output is obtained by the convolution expression:

y(l)=k=0NFD1x(nlk+NFD/2)ha((k+μlNFD/2)T)E21

This means that for a given desired fractional value, the FDF coefficients can be obtained from a designed continuous-time filter.

The design methods using this approach approximate the reconstruction filter ha(t) in each interval of length T by means of a polynomial-based interpolation as follows:

ha((n+μl)T)=m=0Mcm(n)μlmE22

for k=-NFD/2,-NFD/2+1,…., NFD/2-1. The cm(k)’s are the unknown polynomial coefficients and M is the polynomials order.

If equation (Eq. 22) is substituted in equation (Eq. 21), the resulted output signal can be expressed as:

y(l)=m=0Mvm(nl)μlmE23

where:

vm(nl)=k=0NFD1x(nlk+NFD/2)cm(kNFD/2)E24

are the output samples of the M+1 FIR filters with a system function:

Cm(z)=k=0NFD1cm(kNFD/2)zkE25

The implementation of such polynomial-based approach results in the Farrow structure, (Farrow, 1988), sketched in Fig. 11. This implementation is a highly efficient structure composed of a parallel connection of M+1 fixed filters, having online fractional delay value update capability. This structure allows that the FDF design problem be focused to obtain each one of the fixed branch filters cm(k) and the FDF structure output is computed from the desired fractional delay given online μl.

The coefficients of each branch filter Cm(z) are determined from the polynomial coefficients of the reconstruction filter impulse response ha(t). Two mainly polynomial-based interpolation filters are used: 1) conventional time-domain design such as Lagrange interpolation, 2) frequency-domain design such as minimax and least mean squares optimization.

As were pointed out previously, Lagrange interpolation has several disadvantages. A better polynomial approximation of the reconstruction filter is using a frequency-domain approach, which is achieved by optimizing the polynomial coefficients of the impulse response ha(t) directly in the frequency-domain. Some of the design methods are based on the optimization of the discrete-time filter hFD(n,μl)) and others on making the optimization of the reconstruction filter ha(t). Once that this filter is obtained, the Farrow structure branch filters cm(k) are related to hFD(n,ml) using equations (Eq. 20) and (Eq. 22). One of main advantages of frequency-domain design methods is that they have at least three design parameters: filter length NFD, interpolation order M, and pass-band frequency ωp.

There are several methods using the frequency design method (Vesma, 1999). In (Farrow, 1988) a least-mean-squares optimization is proposed in such a way that the squared error between HFD(ω,μl) and the ideal response Hid(ω,μl) is minimized for 0≤ωωp and for 0≤μl<1. The design method reported in (Laakson et al., 1995) is based on optimizing cm(k) to minimize the squared error between ha(t) and the hFD(n,μl) filters, which is designed through the magnitude frequency response approximation approach, see section 3.1. The design method introduced in (Vesma et al., 1998) is based on approximating the Farrow structure output samples vm(nl) as an mth order differentiator; this is a Taylor series approximation of the input signal. In this sense, Cm(ω) approximates in a minimax or L2 sense the ideal response of the mth order differentiator, denoted as Dm(ω), in the desired pass-band frequencies. In (Vesma & Saramaki, 1997) the designed FDF phase delay approximates the ideal phase delay value μl in a minimax sense for 0≤ωωp and for 0≤μl<1 with the restriction that the maximum pass-band amplitude deviation from unity be smaller than the worst-case amplitude deviation, occurring when μ=0.5.

## 4. FDF Implementation structures

As were described in section 3.3, one of the most important results of the analogue-digital model in designing FDF filters is the highly efficient Farrow structure implementation, which was deduced from a piecewise approximation of the reconstruction filter through a polynomial based interpolation. The interpolation process is made as a frequency-domain optimization in most of the existing design methods.

An important design parameter is the FDF bandwidth. A wideband specification, meaning a pass-band frequency of 0.9π or wider, imposes a high polynomial order M as well as high branch filters length NFD. The resulting number of products in the Farrow structure is given by NFD(M+1)+M, hence in order to reduce the number of arithmetic operations per output sample in the Farrow structure, a reduction either in the polynomial order or in the FDF length is required.

Some design approaches for efficient implementation structures have been proposed to reduce the number of arithmetic operations in a wideband FDF. A modified Farrow structure, reported in (Vesma & Samaraki, 1996), is an extension of the polynomial based interpolation method. In (Johansson & Lowerborg, 2003), a frequency optimization technique is used a modified Farrow structure achieving a lower arithmetic complexity with different branch filters lengths. In (Yli-Kaakinen & Saramaki, 2006a, 2006a, 2007), multiplierless techniques were proposed for minimizing the number of arithmetic operations in the branch filters of the modified Farrow structure. A combination of a two-rate factor multirate structure and a time-domain designed FDF (Lagrange) was reported in (Murphy et al., 1994). The same approach is reported in (Hermanowicz, 2004), where symmetric Farrow structure branch filters are computed in time-domain with a symbolic approach. A combination of the two-rate factor multirate structure with a frequency-domain optimization process was firstly proposed in (Jovanovic-Docelek & Diaz-Carmona, 2002). In subsequence methods (Hermanowicz & Johansson, 2005) and (Johansson & Hermanowicz &, 2006), different optimization processes were applied to the same multirate structure. In (Hermanowicz & Johansson, 2005), a two stage FDF jointly optimized technique is applied. In (Johansson & Hermanowicz, 2006) a complexity reduction is achieved by using an approximately linear phase IIR filter instead of a linear phase FIR in the interpolation process.

Most of the recently reported FDF design methods are based on the modified Farrow structure as well as on the multirate Farrow structure. Such implementation structures are briefly described in the following.

### 4.1. Modified Farrow structure

The modified Farrow structure is obtained by approximating the reconstruction filter with the interpolation variable 2μl -1 instead of μl in equation (Eq. 22):

ha((n+μl)T)=m=0Mcm'(k)(2μl1)mE26

for k=-NFD/2,-NFD/2+1,…., NFD/2-1. The first four basis polynomials are shown in Fig. 12. The symmetry property ha(-t)= ha(t) is achieved by:

cm'(n)=(1)mcm'(n1)E27

for m= 0, 1, 2,…,M and n=0, 1,….,NFD/2. Using this condition, the number of unknowns is reduced to half.

The reconstruction filter ha(t) can be now approximated as follows:

ha(t)=n=0NFD/2m=0Mcm'(n)g(n,m,t)E28

where cm(n) are the unknown coefficients and g(n,m,t)’s are basis functions reported in (Vesma & Samaraki, 1996).

The modified Farrow structure has the following properties: 1) polynomial coefficients cm(n) are symmetrical, according to equation (Eq. 27); 2) The factional value μl is substituted by 2μl -1, the resulting implementation of the modified Farrow structure is shown in Fig. 13; 3) the number of products per output sample is reduced from NFD(M+1)+M to NFD(M+1)/2+M.

The frequency design method in (Vesma et al., 1998) is based on the following properties of the branch digital filters Cm(z):

• The FIR filter Cm(z), 0≤mM, in the original Farrow structure is the mth order Taylor approximation to the continuous-time interpolated input signal.

• In the modified Farrow structure, the FIR filters C’m(z) are linear phase type II filters when m is even and type IV when m is odd.

Each filter Cm(z) approximates in magnitude the function Kmwm, where Km is a constant. The ideal frequency response of an mth order differentiator is ()m, hence the ideal response of each Cm(z) filter in the Farrow structure is an mth order differentiator.

In same way, it is possible to approximate the input signal through Taylor series in a modified Farrow structure for each C’m(z), (Vesma et al., 1998). The mth order differential approximation to the continuous-time interpolated input signal is done through the branch filter C’m(z), with a frequency response given as:

C'm(ω)=ejω(NFD1)/2(jω)m2mm!E29

The input design parameters are: the filter length NFD, the polynomial order M, and the desired pass-band frequency ωp.

The NFD coefficients of the M+1 C’m(z) FIR filters are computed in such a way that the following error function is minimized in a least square sense through the frequency range [0,ωp]:

em(ω)=|n=oNFD/21[cm(NFD/2n)]γ(m,n,ω)D(m,ω)|E30

where:

D(m,ω)=(ω)m2mm!,γ(m,n,ω)=2cos[(n+1/2)ω],mevenγ(m,n,ω)=2sin[(n+1/2)ω],moddE31

Hence the objective function is given as:

E1=0ωp[n=0NFD/21cm(NFD/21)γ(m,n,ω)D(m,ω)]dωE32

From this equation it can be observed that the design of a wide bandwidth FDF requires an extensive computing workload. For high fractional delay resolution FDF, high precise differentiator approximations are required; this imply high branch filters length, NFD, and high polynomial order, M. Hence a FDF structure with high number of arithmetic operations per output sample is obtained.

### 4.2. Multirate Farrow structure

A two-rate-factor structure in (Murphy et al., 1994), is proposed for designing FDF in time-domain. The input signal bandwidth is reduced by increasing to a double sampling frequency value. In this way Lagrange interpolation is used in the filter coefficients computing, resulting in a wideband FDF.

The multirate structure, shown in Fig. 14, is composed of three stages. The first one is an upsampler and a half-band image suppressor HHB(z) for incrementing twice the input sampling frequency. Second stage is the FDF HDF(z), which is designed in time-domain through Lagrange interpolation. Since the signal processing frequency of HDF(z) is twice the input sampling frequency, such filter can be designed to meet only half of the required bandwidth. Last stage deals with a downsampler for decreasing the sampling frequency to its original value. Notice that the fractional delay is doubled because the sampling frequency is twice. Such multirate structure can be implemented as the single-sampling-frequency structure shown in Fig. 15, where H0(z) and H1(z) are the first and second polyphase components of the half-band filter HHB(z), respectively. In the same way HFD0(z) and HFD1(z) are the polyphase components of the FDF HFD(z) (Murphy et al, 1994).

The resulting implementation structure for HDF(z) designed as a modified Farrow structure and after some structure reductions (Jovanovic-Dolecek & Diaz-Carmona, 2002) is shown in Fig. 16. The filters Cm,0(z) and Cm,1(z) are the first and second polyphase components of the branch filter Cm(z), respectively.

The use of the obtained structure in combination with a frequency optimization method for computing the branch filters Cm(z) coefficients was exploited in (Jovanovic-Dolecek & Diaz-Carmona, 2002). The approach is a least mean square approximation of each one of the mth differentiator of input signal, which is applied through the half of the desired pass-band. The resulting objective function, obtained this way from equation (Eq. 32), is:

E2=0ωp2[n=0NFD/21cm(NFD/21)γ(m,n,ω)D(m,ω)]dωE33

The decrease in the optimization frequency range allows an abrupt reduction in the coefficient computation time for wideband FDF, and this less severe condition allows a resulting structure with smaller length of filters Cm(z).

The half-band HHB(z) filter plays a key role in the bandwidth and fractional delay resolution of the FDF filter. The higher stop-band attenuation of filter HHB(z), the higher resulting fractional delay resolution. Similarly, the narrower transition band of HHB(z) provides the wider resulting bandwidth.

In (Ramirez-Conejo, 2010) and (Ramirez-Conejo et al., 2010a), the branch filters coefficients cm(n) are obtained approximating each mth differentiator with the use of another frequency optimization method. The magnitude and phase frequency response errors are defined, for 0≤wwp and 0≤μl≤1, respectively as:

emag(ω)=|HFD(ω)|1,E34
epha(ω)=ϕ(ω)ω(Dfix+μl),E35

where HFD(ω) and ϕ(ω) are, respectively, the frequency and phase responses of the FDF filter to be designed. In the same way, this method can also be extended for designing FDF with complex specifications, where the complex error used is given by equation (Eq. 18).

The coefficients computing of the resulting FDF structure, shown in Fig. 16, is done through frequency optimization for global magnitude approximation to the ideal frequency response in a minimax sense. The objective function is defined as:

Δm=max0μl1[max0ωωp|em(ω)|]E36

The objective function is minimized until a magnitude error specification δm is met. In order to meet both magnitude and phase errors, the global phase delay error is constrained to meet the phase delay restriction:

Δp=max0μl1[max0ωωp|ep(ω)|]δpE37

where δp is the FDF phase delay error specification. The minimax optimization can de performed using the function fminmax available in the MATLAB Optimization Toolbox.

As is well known, the initial solution plays a key role in a minimax optimization process, (Johansson & Lowenborg, 2003), the proposed initial solution is the individual branch filters approximations to the mth differentiator in a least mean squares sense, accordingly to (Jovanovic-Delecek & Diaz-Carmona, 2002):

Em=0ωp2[em(ω)]2dωE38

The initial half-band filter HHB(z) to the frequency optimization process can be designed as a Doph-Chebyshev window or as an equirriple filter. The final hafband coefficients are obtained as a result of the optimization.

The fact of using the proposed optimization process allows the design of a wideband FDF structure with small arithmetic complexity. Examples of such designing are presented in section 5.

An implementation of this FDF design method is reported in (Ramirez-Conejo et al., 2010b), where the resulting structure, as one shown in Fig. 16, is implemented in a reconfigurable hardware platform.

## 5. FDF Design examples

The results obtained with FDF design methods described in (Diaz-Carmona et al., 2010) and (Ramirez-Conejo et al., 2010) are shown through three design examples, that were implemented in MATLAB.

Example 1:

The design example is based on the method described in (Diaz-Carmona et al., 2010). The desired FDF bandwidth is 0.9π, and a fractional delay resolution of 1/10000.

A half-band filter HHB(z) with 241 coefficients was used, which was designed with a Dolph-Chebyshev window, with a stop-band attenuation of 140 dBs. The design parameters are: M=12 and NFD=10 with a resulting structure arithmetic of 202 products per output sample.

The frequency optimization is applied up to only ωp=0.45π, causing a notably computing workload reduction, compared with an optimization on the whole desired bandwidth (Vesma et al., 1998). As a matter of comparison, the MATLAB computing time in a PC running at 2GHz for the optimization on half of the desired pass-band is 1.94 seconds and 110 seconds for the optimization on the whole pass-band. The first seven differentiator approximations for both cases are shown in Fig. 17 and Fig. 18.

The frequency responses of the resulted FDF from μ=0.008 to 0.01 samples for the half pass-band and for the whole pass-band optimization process, are shown in Fig. 19 and Fig. 20, respectively.

The use of the optimization process (Vesma et al., 1998) with design parameters of M=12 and NFD=104 results in a total number of 688 products per output sample. Accordingly to the described example in (Zhao & Yu, 2006), using a weighted least squares design method, an implementation structure with NFD=67 and M=7 is required to meet ωp=0.9π, which results in arithmetic complexity of 543 products per output sample.

In order to compare the frequency-domain approximation achieved by the described method with existing design methods results, the frequency-domain absolute error e(ω,μ), the maximum absolute error emax, and the root mean square error erms are defined, like in (Zhao & Yu, 2006), by:

e(ω,μ)=|HFD(ω,μ)Hid(ω,μ)|E39
emax=max{e(ω,μ)},0ωωp,0μ1E40
erms=[0ωp01e2(ω,μ)dμdω]1/2E41

The maximum absolute magnitude error and the root mean square error obtained are shown in Table 1, reported in (Diaz-Carmona et al., 2010), as well as the results reported by some design methods.

 Method emax(dBs) erms (Tarczynski et al., 1997) -100.0088 2.9107x10-6 (Wu-Sheng, & Tian-Bo, 1999) -100.7215 2.7706x10-6 (Tian-Bo, 2001) -99.9208 4.931x10-4 (Zhao & Yu, 2006) -99.3669 2.8119x10-6 (Vesma et al., 1998) -93.69 4.81x10-4 (Diaz-Carmona et al., 2010) -86.17 2.78x10-4

### Table 1.

Magnitude frequency response error comparison.

Example 2:

The FDF is designed using the explained minimax optimization approach applied on the single-sampling-frequency structure, Fig. 16, according to (Ramirez et al., 2010a). The FDF specifications are: ωp = 0.9π, δm = 0.01 and δp =0.001, the same ones as in the design example of (Yli-Kaakinen & Saramaki, 2006a). The given criterion is met with NFD = 7 and M = 4 and a half-band filter length of 55. The overall structure requires Prod = 32 multipliers, Add = 47 adders, resulting in a Δm = 0.0094448 and Δp = 0.00096649. The magnitude and phase delay responses obtained for μl = 0 to 0.5 with 0.1 delay increment are depicted in Fig. 21. The results obtained, and compared with those reported by other design methods, are shown in Table 2. The design described requires less multipliers and adders than (Vesma & Saramaki, 1997), (Johansson & Lowenborg, 2003), the same number of multipliers and nine less adders than (Yli-Kaakinen & Saramaki, 2006a), one more multiplier and three less adders than (Yli-Kaakinen & Saramaki, 2006b), and two more multipliers than (Yli-Kaakinen, & Saramaki, 2007).

 Method Arithmetic complexity NFD M Prod Add Δm Δp (Vesma & Saramaki, 1997) 26 4 69 91 0.006571 0.0006571 (Johansson, & Lowenborg, 2003) 28 5 57 72 0.005608 0.0005608 (Yli-Kaakinen & Saramaki, 2006a) 28 4 32 56 0.009069 0.0009069 (Yli-Kaakinen & Saramaki, 2006b) 28 4 31 50 0.009742 0.0009742 (Yli-Kaakinen & Saramaki, 2007) 28 4 30 - 0.009501 0.0009501 (Ramirez-Conejo et al.,2010) 7 4 32 47 0.0094448 0.0009664

### Table 2.

Arithmetic complexity results for example 2.- Not reported

Example 3:

This example shows that the same minimax optimization approach can be extended for approximating a global complex error. For this purpose, the filter design example described in (Johansson & Lowenborg 2003) is used, which specifications are ωp = 0.9π, and maximum global complex error of δc = 0.0042. Such specifications are met with NFD = 7 and M = 4 and a half-band filter length of 69. The overall structure requires Prod = 35 multipliers with a resulting maximum complex error Δc = 0.0036195. The results obtained are compared in Table 3 with the reported ones in existing methods. The described method requires less multipliers than (Johansson & Lowenborg 2003), (Hermanowicz, 2004) and case A of (Hermanowicz & Johansson, 2005). Reported multipliers of (Johansson & Hermanowicz, 2006) and case B of (Hermanowicz & Johansson, 2005) are less than the obtained with the presented design method. It should be pointed out that in (Johansson & Hermanowicz, 2006) an IIR half-band filter is used and in case B of (Hermanowicz & Johansson, 2005) and (Johansson & Hermanowicz, 2006) a switching technique between two multirate structures must be implemented. The resulted complex error magnitude is shown in Fig. 23 for fractional delay values from D =17.5 to 18.0 with 0.1 increment, magnitude response of the designed FDF is shown in Fig. 24 and errors of magnitude and phase frequency responses are presented in Fig 25.

 Method Arithmetic complexity NFD M Prod (Johansson & Lowenborg 2003) 39 6 73 (Johansson & Lowenborg 2003)a 31 5 50 (Hermanowicz, 2004) 11 6 60(54) (Hermanowicz & Johansson, 2005) 7 5 36 (Hermanowicz & Johansson, 2005)b 7 3 26 (Johansson & Hermanowicz, 2006) - 6 32 (Johansson & Hermanowicz, 2006)b - 3 22 (Ramirez-Conejo et al., 2010) 7 4 35

### Table 3.

Arithmetic complexity results for example 3.a. Minimax design with subfilters jointly optimized.

## 6. Conclusion

The concept of fractional delay filter is introduced, as well as a general description of most of the existing design methods for FIR fractional delay filters is presented. Accordingly to the explained concepts and to the results of recently reported design methods, one of the most challenging approaches for designing fractional delay filters is the use of frequency-domain optimization methods. The use of MATLAB as a design and simulation platform is a very useful tool to achieve a fractional delay filter that meets best the required frequency specifications dictated by a particular application.

## Acknowledgments

Authors would like to thank to the Technological Institute of Celaya at Guanajuato State, Mexico for the facilities in the project development, and CONACYT for the support.

## References

1. 1. Diaz-CarmonaJ.Jovanovic-DolecekG.Ramirez-AgundisA.2010Frequency-based optimization design for fractional delay FIR filters. International Journal of Digital Multimedia Broadcasting, 2010January 2010), 161687-7578
2. 2. ErupL.GardnerF.HarrisF.1993Interpolation in digital modems-part II: implementation and performance. IEEE Trans. on Communications, 41June 1993), 9981008
3. 3. FarrowC.1988A continuously variable digital delay element, Proceedings of IEEE Int. Symp. Circuits and Systems, 26412645Espoo, Finland, June, 1988.
4. 4. GardnerF.1993Interpolation in digital modems-part I: fundamentals. IEEE Trans. on Communications, 41March 1993), 501507
5. 5. Ging-ShingL.Che-HoW.1992A new variable fractional sample delay filter with nonlinear interpolation. IEEE Trans. on Circuits Syst., 39February 1992), 123126
6. 6. Ging-ShingL.Che-HoW.1990Programmable fractional sample delay filter with Lagrange interpolator. Electronics Letters, 2619September 1990), 16081610
7. 7. HermanowiczE.2004On designing a wideband fractional delay filter using the Farrow approach, Proceedings of XII European Signal Processing Conference, 961964Vienna, Austria, September 6-10, 2004.
8. 8. HermanowiczE.JohanssonH.2005On designing minimax adjustable wideband fractional delay FIR filters using two-rate approach, Proceedings of European Conference on Circuit Theory and Design, 473440Cork, Ireland, August 29- September 1, 2005.
9. 9. JohanssonH.LowenborgP.2003On the design of adjustable fractional delay FIR filters. IEEE Trans. on Circuits and Syst.-II, Analog and Digital Signal Processing, 50April 2003), 164169
10. 10. JohanssonH.HermanowiczE.2006Adjustable fractional delay filters utilizing the Farrow structure and multirate techniques, Proceedings Sixth Int. Workshop Spectral Methods Multirate Signal Processing,Florence, Italy, September 12
11. 11. Jovanovic-DolecekG.Diaz-CarmonaJ.2002One structure for fractional delay filter with small number of multipliers, Electronics Letters, 1819September 2002), 10831084
12. 12. LaaksonT.ValimakiV.KarjalainenM.LaineU.1996Splitting the unit delay. IEEE Signal Processing Magazine, 131January 1996), 3060
13. 13. MurphyN.KrukowskiA.KaleI.1994Implementation of wideband integer and fractional delay element. Electronics Letters, 3020September 1994), 16541659
14. 14. OetkenG.1979A new approach for the design of digital interpolation filters. IEEE Trans. on Acoust., Speech, Signal Process., ASSP-27December 1979), 637643
15. 15. ProakisJ.ManolakisD.1995Digital Signal Processing: Principles, Algorithms and Applications, Prentice Hall, 978-0-13373-762-2USA.
16. 16. Ramirez-ConejoG.2010Diseño e implementación de filtros de retraso fraccionario, Master in Science thesis, Technological Institute of Celaya, Celaya Mexico, ( June, 2010).
17. 17. Ramirez-ConejoG.Diaz-CarmonaJ.Delgado-FriasJ.Jovanovic-DolecekG.Prado-OlivarezJ.2010aAdjustable fractional delay FIR filters design using multirate and frequency optimization techniques, Proceedings of Trigésima Convención de Centroamérica y Panáma del IEEE, CONCAPAN XXX, San José, Costa Rica, November 1718
18. 18. Ramirez-ConejoG.Diaz-CarmonaJ.Delgado-FriasJ.Padilla-MedinaA.Ramirez-AgundisA.2010bFPGA implementation of adjustable wideband fractional delay FIR filters, Proceedings International Conference On Reconfigurable Computing and FPGAs, December 1315
19. 19. RamstadT.1984Digital methods for conversion between arbitrary sampling frequencies. IEEE Trans. on Acoust. Speech, Signal Processing, ASSP-32June 1984), 577591
20. 20. TarczynskiA.CainG.HermanoviczE.RojewskiM.1997WLS design of variable frequency response FIR response, Proceedings IEEE International Symp. On Circuits and Systems, 22442247Hong Kong, June 9-12, 1997.
21. 21. Tian-BoD.2001Discretization-free design of variable fractional-delay FIR filters. IEEE Trans. On Circuits and Systems-II: Analog and Digital Signal Processing, 48June 2001), 637644
22. 22. VesmaJ.1999Optimization Applications of Polynomial-Based Interpolation Filters, PhD thesis, University of Technology, Tampere Finland, (May 1999), 9-52150-206-1
23. 23. VesmaJ.1995Timing adjustment in digital receivers using interpolation, Master in Science thesis, University of Technology, Tampere Finland, (November 1995).
24. 24. VesmaJ.SamarakiT.1996Interpolation filters with arbitrary frequency response for all digital receivers, Proceedings of IEEE Int. Sust. Circuits Syst., 568571Georgia, USA, May 12-15, 1996.
25. 25. VesmaJ.HamilaR.SaramakiT.RenforsM.1998Design of polynomial-based interpolation filters based on Taylor series, Proceedings of IX European signal Processing Conference, 283286Rodhes, Greece, September 8-11, 1998.
26. 26. VesmaJ.SaramakiT.1997Optimization and efficient implementation of FIR filters with adjustable fractional delay, Proceedings IEEE Int. Symp. Circuits and Systems, 22562259Hong Kong, June 9-12, 1997.
27. 27. Wu-ShengL.Tian-BoD.1999An improved weighted least-squares design for variable fractional delay FIR filters. IEEE Trans. On Circuits and Systems-II: Analog and Digital Signal Processing, 46August 1999), 10351049
28. 28. Yli-KaakinenJ.SaramakiT.2006aMultiplication-free polynomial based FIR filters with an adjustable fractional delay. Springer Circuits, syst., Signal Processing, 25April 2006), 265294
29. 29. Yli-KaakinenJ.SaramakiT.2006bAn efficient structure for FIR filters with an adjustable fractional delay, Proceedings of Digital Signal Processing Applications, 617623Moscow, Russia, March 29-31, 2006.
30. 30. Yli-KaakinenJ.SaramakiT.2007A simplified structure for FIR filters with an adjustable fractional delay, Proceedings of IEEE Int. Symp. Circuits and Systems, 34393442New Orleans, USA, May 27-30, 2007.
31. 31. ZhaoH.YuJ.2006A simple and efficient design of variable fractional delay FIR filters. IEEE Trans. On Circuits and Systems-II: Express Brief, 53February 2006), 157160

Written By

Javier Diaz-Carmona and Gordana Jovanovic Dolecek

Submitted: 25 March 2011 Published: 09 September 2011