Abstract

We prove that the Chebyshev sharpening technique, recently introduced in literature, provides filters with a Minimum Phase (MP) characteristic when it is applied to cosine filters. Additionally, we demonstrate that cascaded expanded Chebyshev-Sharpened Cosine Filters (CSCFs) are also MP filters, and we show that they achieve a lower group delay for similar magnitude characteristics in comparison with traditional cascaded expanded cosine filters. The importance of the characteristics of cascaded expanded CSCFs is also elaborated. The developed examples show improvements in the group delay ranged from 23% to 47% at the cost of a slight increase of usage of hardware resources. For an application of a low-delay decimation filter, the proposed scheme exhibits a 24% lower group delay, with 35% less computational complexity (estimated in Additions per Output Sample) and slightly less usage of hardware elements.

1. Introduction

A Minimum Phase (MP) digital filter has all zeros on or inside the unit circle [1]. We consider in this paper MP Finite Impulse Response (FIR) filters, which find applications in cases where a long delay, usually introduced by Linear Phase (LP) FIR filters, is not allowed. Examples of these cases include data communication systems or speech and audio processing systems [24].

The basic building block analyzed in this paper, the cosine filter, is a simple FIR filter whose transfer function and frequency response are, respectively, given by This filter is of special interest because of the following main reasons: (a)It has MP property because its zero lies on the unit circle.(b)It has a low computational complexity because it does not require multipliers, which are the most costly and power-consuming elements in a digital filter [5].(c)It has a low usage of hardware elements, which can be translated into a low demand of chip area for implementation.

The fact that a cosine filter has the Minimum Phase characteristic becomes significant because these basic building blocks can be used to design filters with a low delay and simultaneously a low computational complexity and a low usage of hardware resources.

Due to the aforementioned characteristics, cascaded expanded cosine filters were investigated in [6]. A cascaded expanded cosine filter is defined as a filter with transfer function and frequency response, respectively, given bywith being the cosine filter given in (1), whereas and are arbitrary integers. In method [6], Rouche’s theorem was employed to demonstrate that when cascaded expanded cosine filters are sharpened with a modified version of the technique from [7] (originally devised for LP FIR filters), the result is an overall FIR filter that has all its zeros on or inside the unit circle; that is, it satisfies the MP characteristic. Nevertheless, the resulting filter still has a high group delay because a large number of cascaded expanded cosine filters are needed to meet a given attenuation specification, as we see in the examples from [6].

On the other hand, the Chebyshev sharpening technique was recently introduced in [8] to design LP FIR filters based on comb subfilters for decimation applications. In that method the sharpening is performed with an th degree Chebyshev polynomial of first kind, defined as where are integers [9]. When applied to comb filters, Chebyshev sharpening provides solutions with advantages like a simple and elegant design method, a low-complexity resulting LP FIR filter, and improved attenuation characteristics in the resulting filter. However, filters from [8] are not guaranteed to have MP characteristic.

From the aforementioned literature we can extract the following observations:(a)In MP FIR filters the reduction of the group delay is a priority.(b)The use of cosine filters results in low-complexity multiplierless MP FIR filters.(c)The recent Chebyshev sharpening method from [8] can improve the attenuation of cosine filters and is a potentially useful approach to preserve a simple multiplierless solution with a lower group delay in comparison with simple cascaded expanded cosine filters.

Motivated by the remarks listed above, we present in this paper the following contributions:(1)The mathematical proof that the use of Chebyshev sharpening in cosine filters, which produces Chebyshev-Sharpened Cosine Filters (CSCFs), guarantees resulting multiplierless filters with all of their zeros placed on the unit circle, that is, with MP property: this demonstration hinges upon the factorization of the transfer function of the CSCF into second-order sections, taking advantage of the antisymmetry of the roots of the Chebyshev polynomial.(2)The mathematical proof that cascaded expanded CSCFs have also all of their zeros placed on the unit circle: this demonstration is a simple extension of aforementioned proof for CSCFs.(3)The explanation of how cascaded expanded CSCFs can be efficiently employed in the design of MP FIR filters.(4)Examples where it is shown that cascaded expanded CSCFs are useful to obtain the same stopband attenuation as cascaded expanded cosine filters but with a lower group delay: from these examples we see an improvement from 23% to 47% in the reduction of the group delay, at the cost of a slight increase of the usage of hardware resources. For an application in a decimation filter embedded into a low-delay oversampled Analog-to-Digital Converter (ADC), the proposed scheme exhibits a 24% lower group delay referred to high rate, with 35% less computational complexity (estimated in Additions per Output Sample) and slightly less usage of hardware elements.

Following this introduction, Section 2 presents the definition of CSCFs and cascaded expanded CSCFs. The proofs of MP characteristic in CSCFs and cascaded expanded CSCFs are given in Sections 3 and 4, respectively. In Section 5 we provide details on the characteristics and applications of the cascaded expanded CSCFs. Then, Section 6 presents examples and discussion of results. Finally, concluding remarks are given in Section 7.

2. Definition of Chebyshev-Sharpened Cosine Filter (CSCF) and Cascaded Expanded CSCF

We define the transfer function and the frequency response of an th order Chebyshev-Sharpened Cosine Filter (CSCF), respectively, aswithwhere are the coefficients of the Chebyshev polynomial of first kind, represented in (5), and is given in (1). To obtain a low-complexity multiplierless implementation, the constant must be expressible as a Sum of Powers of Two (SOPOT). To this end, we set where denotes “the closest value less than or equal to that can be realized with at most adders” and denotes rounding to the closest integer less than or equal to . To provide an improved attenuation around the zero of the cosine filter, must be as close as possible to its upper limit [8]. This is achieved by increasing the integer . The value in (8)-(9) is usually set as an integer equal to or greater than 2 for applications in decimation processes [8]. However, we will allow for more flexibility to the parameter in this paper, as will be explained in the next section.

The transfer function and frequency response of a cascaded expanded CSCF are, respectively, defined aswhere the integer indicates the number of cascaded CSCF blocks, each of them repeated times, with . Every value of is a distinct factor that expands a different CSCF whose corresponding order is . These CSCFs have different factors , which can be obtained using (9), just replacing by and by , where and are integer parameters that correspond to the th CSCF in the cascade. Figure 1(a) shows the structure of the CSCF, where we have that , with and with if is odd or if is even. Dashed blocks in Figure 1(a) appear only if is odd. Figure 1(b) presents the structure of the cascaded expanded CSCF whose transfer function is given in (10).

3. Proof of Minimum Phase Property in CSCFs

The proof starts with the expression of the Chebyshev polynomial from (5) in the form of a product of first-order terms as [9] On the other hand, we rewrite the transfer function of the CSCF from (6) as Using (12), and after simple rearrangement of terms, we express as follows:which can be rewritten aswhere denotes rounding to the closest integer greater than or equal to .

At this point, it is worth highlighting that the antisymmetry relationshold [9]. Thus, replacing (17) in (16), and after simple manipulation of terms, we haveFrom (18) we have that consists of a product of either several terms if is even or several terms and a term if is odd, with . Thus, to prove the MP property of the CSCF, it is only necessary to ensure that and have MP characteristic for all values .

Using (1), it is easy to see that the term has a root on the unit circle and thus it corresponds to a MP filter. On the other hand, replacing (1) into (19) and after simple rearrangement of terms, we get From (20) it is easy to show that the roots of () are placed on the unit circle; that is, if the argument in (22) is preserved into the range . From (13) we have that holds. Additionally, by settingin (8)-(9), we ensure . Under this condition for , we have that holds. In this case, has its roots on the unit circle for all the valid values and, as a consequence, the filter has a MP characteristic.

Figure 2 shows the pole-zero plots for the filters , , , and . For all of these filters, we have , which is implemented with just one subtraction.

4. Proof of Minimum Phase Property in Cascaded Expanded CSCFs

The proof starts with the expression of every CSCF of the cascaded expanded CSCFs from (10) in the form of a product of second-order expanded transfer functions using (18) and (20); that is,where and . Since the transfer function of the cascaded expanded CSCF from (10) consists of a product of several terms with different values , it is only necessary to ensure that has a MP characteristic for all values . Moreover, from (24) we see that is expressed as a product of either several terms if is even or several terms and if is odd. Thus, to prove the MP property in cascaded expanded CSCFs we only need to ensure that and have MP characteristic for all values and .

By replacing (1) in the term and then making the resulting expression equal to zero, we can find the roots of . These roots turn out to be the complex roots of , which have unitary magnitude. Thus, has MP characteristic, since its roots are placed on the unit circle. On the other hand, using (21) we can express (25) as follows: To preserve the argument in (27) into the range , we setUnder this condition for , we have that holds. In this case, the respective roots of factors () and () in (25) are the roots of the complex numbers and , which have unitary magnitude for all the valid values . Therefore, has MP characteristic, since its roots are placed on the unit circle. Finally, since and have MP characteristic, the overall cascaded expanded CSCF from (9), , also has MP characteristic.

Figure 3 shows the pole-zero plots for the filters , , , and . For all these filters, we have , which is implemented with just one subtraction.

5. Characteristics and Applications of Cascaded Expanded CSCFs

A cascaded expanded CSCF has both MP and LP characteristics. The former was proven in Section 4, whereas the latter is easily seen from the frequency response given in (11). A consequence of this is that the cascaded expanded CSCF has a passband droop in its magnitude response. Due to this passband droop, the cascaded expanded CSCF should be employed only to provide a given attenuation requirement of an overall MP FIR filter over a prescribed stopband region (depending on the application). Thus, the resulting structure to design an overall MP FIR filter can be associated with the prefilter-equalizer scheme of [10], shown in Figure 4, where the prefilter provides the required attenuation whereas the equalizer corrects the passband droop of the prefilter. The cascaded expanded CSCF, with transfer function defined in (10), can be used as prefilter. Note that since a cascaded expanded cosine filter (whose transfer function is defined in (3)) also has both LP and MP properties, it is used as prefilter in [6].

Even though this paper is not focused on the design of the equalizer, it is worthwhile to spend some words on how this filter could be designed. An Infinite Impulse Response (IIR) filter with optimally located poles based, for example, in the Least Squares criterion as shown in [11] can form a proper equalizer. However, FIR filters are usually preferred over their IIR counterparts because they have guaranteed stability, they are free of limit-cycle oscillations, and their polyphase decomposition in multirate schemes allows them to reduce the computational load, among other characteristics [12]. Thus, we are more concerned here with FIR equalizers. Since a FIR equalizer with LP characteristic has its zeros placed in quadruplets around the unit circle [1], it does not accomplish the MP characteristic. Therefore, a MP FIR equalizer (i.e., that filter whose zeros appear inside the unit circle) does not have a Linear Phase.

When a LP FIR filter is designed by sharpening cascaded expanded cosine filters with the traditional sharpening polynomial from [7], the resulting filter has a prefilter given by and a LP equalizer given by , where is the group delay of used to preserve the Linear Phase characteristic. In method [6] the delay has been removed to obtain a MP FIR equalizer. Thus, a first option would be to use the same approach of [6] to design a FIR equalizer. However, it is worth highlighting that the recent LP droop compensators proposed in literature (see, e.g., [13, 14]) are novel low-complexity alternatives to the aforementioned LP equalizer based on the traditional sharpening. Inspired by these alternatives, a more convenient approach would be to design MP droop compensators as counterparts of the MP equalizer based on sharpening, proposed in [6].

Besides method [6], other design methods for MP FIR filters have been introduced, for example, in [1522]. They can be classified in general terms as methods based on cepstrum [1517] and methods based on the design of a LP FIR filter [1822]. However, in general, these methods have the inconvenience of producing filtering solutions that require multipliers, which are the most costly elements in a digital filter [5]. To solve this problem, the cascaded expanded CSCF can be used as a prefilter to implement an overall MP FIR filter using several multiplierless CSCFs. A similar approach can be followed with the use of a cascaded expanded cosine filter from [6]. Nevertheless, as we mentioned earlier, the problem with method [6] is that the resulting filter requires a large cascade of expanded cosine filters, increasing the group delay of the resulting filter.

Finally, it is worth highlighting that, in comparison to the filter , the filter has many more parameters, namely, , , , and , with , to be tuned in order to find a desired attenuation. This characteristic provides more flexibility for the design of MP FIR filters in comparison to . Moreover, by setting , = 1, , and for all in (10), we obtain the same expression as (3). Thus, the cascaded expanded CSCF from (10) is a generalized case of (3).

6. Examples and Discussion

This section presents a couple of examples (Examples 1 and 2) that compare the cascaded expanded CSCFs with cascaded expanded cosine filters from [6]. This comparison is made in terms of(a)group delay, measured in samples and defined as follows [1]:where is the frequency response of the corresponding filter;(b)implementation complexity, measured in the required number of adders and delays for a given attenuation over a prescribed stopband region.

Additionally, an engineering application is provided in Example 3, namely, the antialiasing filtering process used in the first stage of a two-stage decimation structure applied in a low-delay Sigma-Delta ADC for audio systems, detailed in [3]. In this case, comparisons are made in terms of group delay referred to high rate, computational complexity counted in Additions per Output Sample (APOS), and number of hardware elements assuming that both filters, the one used in method [3] and the proposed filter, are implemented in recursive form.

Example 1. Design a MP FIR filter with minimum attenuation equal to 60 dB over the range from to (see Figure 1 of [6]).

In [6], the filter employed to accomplish such characteristic is obtained from (3) using and . The group delay is obtained by replacing these values in (4) and then using (4) in (29). This filter requires 15 adders and 45 delays, but it has a group delay of 22.5 samples.

If we use , , , , , , and , with and for all in (10), we get a filter whose group delay, obtained by replacing the aforementioned parameters in (11) and then using (11) in (29), is 16 samples, that is, nearly 30% less delay than that of [6]. Since this filter uses 30 adders and 44 delays, the price to pay is of additional implementation complexity. Figure 5 shows the magnitude responses and group delays of both filters. Moreover, Tables 1 and 2 present, respectively, the first half of the symmetric impulse response of the filter designed with method [6] and the proposed filter.

Example 2. Design a MP FIR filter with minimum attenuation equal to 100 dB over the range from ω = 0.15π to ω = π (see Example 1 of [6]).

In [6], the filter employed to accomplish such characteristic is obtained from (3) using and . The group delay is obtained by replacing these values in (4) and then using (4) in (29). This filter requires 24 adders and 120 delays. However, its group delay is 56 samples.

By using , = 4, = 6, = 4, = 8, = 6, = 2, = 0.8, = 1.2, = 2, = = 3, = 4, and = 1 for all in (10), we get a filter whose group delay, obtained by replacing the aforementioned parameters in (11) and then using (11) in (29), is 30 samples, that is, approximately 47% less delay than that of [6]. This filter uses 58 adders and 90 delays; thus the price to pay is a of additional implementation complexity. Figure 6 shows the magnitude responses and group delays of both filters. Tables 3 and 4, respectively, show the first half of the symmetric impulse response of the filter of method [6] and the proposed filter in Example 2.

Table 5 summarizes the results from the previous examples. From them we observe that the cascaded expanded CSCFs achieve a lower group delay in comparison to the cascaded expanded cosine filters from [6]. This characteristic, desirable for MP filters, occurs because CSCFs take advantage of the zero-rotation effect (see Figures 2 and 3), provided by Chebyshev sharpening, to achieve a given attenuation using less cascaded filters. Nevertheless, since every CSCF needs in general a few more hardware resources (adders and delays) than their cascaded cosine counterparts, the price to pay is an increase in the implementation complexity. Note, however, that the resulting filters still are low-complexity multiplierless solutions.

In the following example we show the design of the antialiasing filter used in the first stage of a two-stage decimation structure applied in a low-delay Sigma-Delta ADC for audio systems, where the second stage is a FIR filter with droop compensation characteristic followed by a downsampler by 2 [3].

Example 3. Design a MP FIR filter for a decimation factor and residual decimation factor , with minimum attenuation equal to 95 dB over the range of frequencies from to , where these frequencies are given by

In [3], the filter employed to accomplish such characteristic is obtained from method [23]. Its transfer function and frequency response are, respectively, given as This filter has a group delay of samples. Moreover, its Cascaded Integrator-Comb (CIC) structure, shown in Figure 7, performs 10 32 + 10 = 330 Additions per Output Sample (APOS) and uses 20 adders and 20 delays.

Now consider the proposed filter, whose transfer function is given by where is an expanded CSCF given byThe frequency response is given byThis filter has a group delay of samples, which is approximately 24% less delay than that of [3]. The filter can be moved after a downsampling by 16 because it is actually a CSCF expanded by 16. Thus, the resulting CIC-based structure, shown in Figure 8, performs (6 × 32 + 6) + (3 × 2 + 5 + 6) = 215 Additions per Output Sample (APOS), that is, nearly 35% less computational complexity with regard to the filter used in [3]. Moreover, this filter uses 20 adders and 16 delays, which represents 10% lower usage of hardware resources compared with method [3].

Figure 9 shows the magnitude responses and group delays referred to high rate of both the filter used in method [3] and the proposed filter. Tables 6 and 7 show the first half of the symmetric impulse response referred to high rate of the filters obtained with method [3] and the proposed method, respectively, whereas the summary of results is given in Table 8.

It is worth highlighting that the implementation of the comb decimation filter in a CIC structure has been employed in method [3] due to its regularity and simplicity, which has a low usage of hardware resources (see Figure 7). However, the price to pay for such desirable characteristics is a high computational complexity. Our proposed solution has taken advantage of the possibility of factorize the decimation factor = 32 as = 16 × 2. With this we have used an expanded-by-16 CSCF as an additional filter that contributes to improving the attenuation in the first stopband, where the comb filter has the worst attenuation. The first advantage of doing so is observed in the reduction of the number of Integrator-Comb stages from 10 to 6. Moreover, since the CSCF can operate at a sampling rate reduced by 16, the computational complexity of the decimation process is reduced and the number of hardware elements is also reduced. Of course one can resort to other types of architectures, such as the nonrecursive form of the comb filter and its subsequent polyphase decomposition. However, this decreases the computational complexity at the cost of a considerable increase of the number of hardware elements.

7. Conclusion

In this paper we have presented the mathematical demonstration that the application of Chebyshev sharpening to cosine filters results in filters with zeros on the unit circle, that is, with Minimum Phase (MP) characteristic. From this, we have proven that filters composed by a cascade of Chebyshev-Sharpened Cosine Filters (CSCFs) expanded by different factors, called cascaded expanded CSCFs, also have MP property. The cascaded expanded CSCFs are useful prefilters that provide the attenuation in an overall MP FIR filter. Moreover, these filters are a general case where the cascaded expanded cosine filters are a subset. The CSCFs are low-complexity filters, since they do not need multipliers.

It has been shown with three examples that, for a desired attenuation in the magnitude response, cascaded expanded CSCFs achieve a lower group delay in comparison to cascaded expanded cosine filters. This lower group delay is desirable in the design of MP FIR filters. Since the purpose of this paper is to prove the suitability of cascaded expanded CSCFs as MP prefilters, the CSCF-based solutions provided in the examples of this work are suboptimal.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.