The Scientific World Journal

Volume 2014 (2014), Article ID 950860, 9 pages

http://dx.doi.org/10.1155/2014/950860

Research Article

Optimal Sharpening of Compensated Comb Decimation Filters: Analysis and Design

David Ernesto Troncoso Romero,¹ Massimiliano Laddomada,² and Gordana Jovanovic Dolecek¹

¹Electronics Department, INAOE Institute, 72840 Tonantzintla, PUE, Mexico

²Electrical Engineering Department, Texas A&M University, Texarkana, TX 75505, USA

Received 31 August 2013; Accepted 20 October 2013; Published 22 January 2014

Academic Editors: Z. Hou and J. C. Nieto-Borge

Copyright © 2014 David Ernesto Troncoso Romero 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.

Abstract.

Comb filters are a class of low-complexity filters especially useful for multistage decimation processes. However, the magnitude response of comb filters presents a droop in the passband region and low stopband attenuation, which is undesirable in many applications. In this work, it is shown that, for stringent magnitude specifications, sharpening compensated comb filters requires a lower-degree sharpening polynomial compared to sharpening comb filters without compensation, resulting in a solution with lower computational complexity. Using a simple three-addition compensator and an optimization-based derivation of sharpening polynomials, we introduce an effective low-complexity filtering scheme. Design examples are presented in order to show the performance improvement in terms of passband distortion and selectivity compared to other methods based on the traditional Kaiser-Hamming sharpening and the Chebyshev sharpening techniques recently introduced in the literature.

1. Introduction

Efficient decimation filtering for oversampled discrete-time signals is key in the development of low-power hardware platforms for reconfigurable communication transceivers [1–26]. From a practical point of view, decimation is usually accomplished using a cascade of two (or more) stages. The filter in the first stage is a comb filter of order decimating by a factor , with -transfer function and zero-phase frequency response, respectively, given as

Comb filters are used in the first stage of the decimation chain because their system function is simple and it does not require any multiplier. However, their magnitude response exhibits a considerable passband droop in the passband ,

where is the total decimation factor and is the residual decimation factor of the remaining stages in the multistage architecture. Furthermore, comb filters have low attenuation in the folding bands defined as

where stands for the integer part of .

Owing to their reduced computational complexity, research on comb filters to date has been focused on (1) improving the magnitude characteristic, (2) preserving linearity of phase, and (3) having the least possible increase of computational complexity [2–24]. With this background, let us review the literature in these three categories.

From the representative sample of works improving the magnitude characteristics of comb filters, we observe that the rotated-comb-based schemes [2–7] have the disadvantage of being susceptible to imperfect pole-zero cancelation. An effective way to prevent this problem consists in designing nonrecursive filters [3, 4, 7] with filtering implemented in polyphase form for ensuring power savings. However, this may result in higher demand for chip area. Other approaches improve the passband with low-order compensators and stopband attenuation by either increasing the order of the comb filter [8–14] or exploiting additional filtering at high rate [15–17]. These approaches provide low-complexity solutions, but the passband improvement cannot be completely controlled. Therefore, these methods are convenient when the desired magnitude characteristics are not too stringent, and when the bandwidth of interest is narrow.

On the other hand, the techniques relying on sharpening of comb-based filters in [18–24] are effective because they can take advantage of the structure proposed in [18], which harnesses the recursive form of comb filters resulting in Cascaded Integrator Comb- (CIC-) like architectures that move part of the filtering at lower rate. Additionally, this structure has all the sharpening coefficients at lower rate and, when integer coefficients scaled by a power of two are used, an effective overall structure is obtained, which does not suffer from finite-precision effects as rotated-comb-based methods.

Moreover, two-stage comb-based decimation schemes have gained great popularity because the comb decimation filter in the first stage, designed in nonrecursive form, can be implemented at lower rate by polyphase decomposition, thus resulting in lower power consumption. The second-stage filtering operates at lower rate as well, but it can take advantage of CIC-like architectures for area reduction. By doing so, the overall comb-based decimation scheme achieves power and area savings. This approach has been applied to traditional comb filters [25, 26] and to magnitude-improved comb filters [5, 15, 17, 20, 21, 24].

2. Problem Motivation, Contributions, and Paper Organization

The reasons at the very basis of this work stem from the following observations.

(a)Sharpened compensated comb filters [23, 24] based on the simplest polynomial of the traditional Kaiser-Hamming sharpening from [27] provide a good passband improvement over conventional comb filters. However, in methods [23, 24] the filter designer does not have control on the exact passband deviation and stopband attenuation achieved by the designed filter.(b)In two-stage comb-based decimation schemes, magnitude response improvements over the passband and the first folding band can be achieved by improving only the second-stage comb filter. However, in these cases, the filter in the first stage introduces a passband droop that cannot be corrected neither by resorting to traditional Kaiser-Hamming sharpening [27] nor by using the recent Chebyshev sharpening [22] applied to the comb filter placed in the second stage. Thus, a different sharpening approach has to be pursued.

In the light of the previous observations, the contributions of this work are the following.

(a)We show that, for similar magnitude characteristics, sharpened compensated comb filters guarantee lower complexity than sharpened comb filters without compensation, especially when stringent specifications must be met.(b)We introduce a low-complexity structure in which the simple multiplierless compensator can be embedded into the cascaded chain of comb filters working at lower rate.(c)We detail the optimization framework to design sharpened comb-based filters to attain given specifications on the acceptable maximum passband distortion and selectivity. The optimized sharpening coefficients are finite-precision values resulting in multiplierless structures, which are important for low-power applications. The optimization problem can be straightforwardly solved with a simple routine of the MATLAB Optimization Toolbox (available online at [28]).

The rest of this paper is organized as follows. Section 3 presents a summary of the generalized perspective of sharpening comb-based filters for decimation. The proposed filtering structure and the corresponding guidelines to decide when to use sharpened compensated filters instead of sharpened comb filters without compensation are introduced in Section 4. The optimization framework to design sharpened compensated comb filters along with the key design steps is provided in Section 5. Section 6 highlights the characteristics to be considered for the sharpening of the second-stage filter in a two-stage comb-based architecture. Design examples are presented in Section 7 where the goal is to contrast the magnitude responses of the proposed method against existing techniques, namely, the ones based on traditional Kaiser and Hamming sharpening and the Chebyshev sharpening methods recently proposed in the literature. Comparisons in terms of computational complexity quantified in Additions Per Output Sample (APOS) are also included in that section. Finally, concluding remarks are presented in Section 8.

3. Generalized Perspective of Sharpening Comb-Based Filters for Decimation

Let and be, respectively, the transfer function and the zero-phase frequency response of an arbitrary comb-based filter to be sharpened (referred to hereafter as subfilter). Any arbitrary th degree sharpening polynomial,

allows mapping the amplitude values to new amplitude values . The new values must approximate the desired values for , where and , with tolerances (for the approximation over the region ) and (for the approximation over the region ). In this way, the zero-phase frequency response of the sharpened filter achieves the desired values with a maximum absolute passband deviation over the range of where and a maximum absolute stopband deviation over the range of where . Usually, for and for . Thus, the sharpening polynomial must meet the following simultaneous conditions:

For the comb-based decimation filter, the range limits for and are

where and are, respectively, given in (3) and (4).

4. Proposed Sharpened Compensated Comb Filters

Let us consider as subfilter the simplest compensated comb filter, which has the following transfer function [11]:

The zero-phase frequency response is

The work [23] has shown that sharpening the subfilter with the polynomial (a polynomial obtained using method [27]) results in significant improvement of the passband characteristic. In that method, the magnitude in the stopband regions can be arbitrarily improved with the order of the comb filter, , and the parameter must be adjusted accordingly. The polynomial has been chosen in [23] because this is the simplest sharpening polynomial from [27] that can improve the passband. Due to that simplicity, it is inferred in [23] that the resulting sharpened compensated comb filter will have a low computational complexity.

In this paper, we propose to use the general sharpening polynomial from (5), finding the coefficients through optimization. The transfer function of the sharpened comb-based filter and its zero-phase frequency response are, respectively, given as

Figure 1 presents the proposed structure to efficiently implement a decimation filter in a CIC-like form. The structure is straightforwardly derived from the combination of both the structure from [18] and the structure introduced in [23] for the special case . Note that must be an even value to avoid fractional delays.

Figure 1: Proposed CIC-like structure for decimation filtering with sharpened compensated comb filters.

The computational complexity of this structure measured in Additions Per Output Sample (APOS) is given by

where indicates the number of adders required to implement the sharpening coefficient and is the number of nonzero sharpening coefficients. For comparison purposes, we present the computational complexity of a sharpening structure for comb filters without compensation (structure from [18]), which is given as

with being the degree of the sharpening polynomial used in that structure.

From (11) and (12), we can see that the highest impact on the APOS complexity metric depends on the products and , respectively. Hence, and (or ) should be chosen as smaller values as possible for any arbitrary decimation factor . Both and (or ) have the same impact on the APOS metric. However, can only take even values, whereas (or ) can also be odd. Since decreasing is therefore more convenient, we can set in advance. This choice leads us to use in agreement with [23].

Let us discard the computational complexity introduced by the sharpening coefficients in both (11) and (12) and assume for the reason discussed above. We will compare the terms in (11) and in (12), assuming that , with being an integer. With this setup, we have

Note that and differ in the terms versus . Clearly, the proposed structure can have a lower computational complexity (i.e., ) when ; that is, when . In that case, sharpened compensated comb filters become convenient when

At this point, it is important to mention that we can take advantage of the frequency transformation approach [29] to estimate the minimum degree of the sharpening polynomial with any of the formulas from [30–32]. These formulas are expressible as , where and are, respectively, the desired maximum absolute passband and stopband deviation of the sharpened filter, whereas and are given as [29]

with , , , and given in (7). Obviously, this is a preliminary estimation that depends on the accuracy of the formula being used. Substituting (7) in (15) and using the Kaiser formula [31], we can estimate (or ) as

Upon noticing that the shape of the magnitude response of comb filters changes very little with [33], we set in advance and we estimate the degrees and using the Kaiser formula for some typical values of , , and ; namely, = 2, 4, 6, and 8; (≈0.01 dB) and (60 dB), 0.0001 (80 dB) and 0.00001 (100 dB) (to estimate , the values , , , and must be obtained from (7) but first replacing from (2) instead of in these equations. To estimate , we use from (9) in these equations, first replacing and in (9)). These cases are shown in Figure 2. Note that, for these specifications, sharpened compensated comb filters are convenient when the residual decimation factor is equal to 2 or 4, that is, for small values. However, generally speaking, sharpened compensated comb filters become effective as the passband and stopband specifications become more stringent. From the previous analysis, we derive the following two important observations.

(a)In sharpened compensated comb filters, a lower computational complexity is obtained if . This is because both and have the same impact on the APOS metric. However, can only take even values, whereas can also be odd. Therefore, preserving a simple sharpening polynomial and improving the stopbands with the increase of , as suggested in [23], do not guarantee a result with low computational complexity.(b)Upon comparing the sharpened comb and sharpened compensated comb filters using , the former requires a polynomial with higher degree. As a consequence, its complexity is higher, despite the use of compensators in the latter. The reason is that the increased complexity in the sharpened compensated comb structures amounts to only 3 extra additions per polynomial degree (when the compensator from [11] is used), and these additions work at lower rate.

(a)

(b)

(c)

Figure 2: Estimated degree of the sharpening polynomial to sharpen comb filters (dashed line) and compensated comb filters (solid line), with and (a), (b), and (c).

5. The Optimization Framework to Design Sharpened Comb Filters

Now, we introduce the optimization framework to obtain the discrete coefficients of for which the maximum deviation of with respect to , denoted by , is minimized. Note that this polynomial will attain the desired passband and stopband deviations with a proper polynomial degree. To find the sharpening polynomial coefficients, we evaluate the conditions (6) over a dense grid of points covering the ranges and .

Let us consider the following notation in order to formalize the optimization problem.

(i) and are the desired passband and stopband deviations after sharpening.(ii) is the number of points partitioning the frequency sets and ; that is, the overall number of grid points in the region is .(iii) and are the th points belonging to the sets and , respectively. To find these points, we divide the range of frequencies into equally spaced points and the range of frequencies into equally spaced points , with . Then, we set and .(iv) and are the desired amplitudes of the polynomial at the points and , respectively. Usually, and for all .(v) denotes the entry in the th row and th column of the underlined matrix .(vi) denotes the th element of the underlined vector .(vii) is an arbitrary word-length for the fractional part in a fixed-point representation of the sharpening coefficients. In other words, every sharpening coefficient is an integer scaled by :

By this setup, the optimization problem proposed in this paper can be written as

where and are vectors of size , is a matrix of size and is a vector of size . In addition, we have

The optimization problem in (18) is a constrained mixed integer linear programming (MILP) problem whose solution can be obtained with generic MILP solvers. As Coleman pointed out in [22], these optimization resources could be inaccessible to many designers. However, the size of this problem is generally small and the simple MATLAB code available online [28] can be used straightforwardly. Such routine is based on the linprog function belonging to the MATLAB Optimization Toolbox. Once the vector has been obtained, the sharpening coefficients can be found as follows:

We notice in passing that a somewhat similar optimization approach was derived by Candan and made available online at [34], along with an extensive MATLAB code that, in general terms, finds the infinite-precision coefficients using the linprog function. However, the work [34] does not provide any method to find optimal discrete coefficients and simple rounding has been applied to the infinite precision solution, making pointless the infinite-precision optimization. Moreover, method [34] is focused on sharpening traditional comb filters without compensation.

5.1. Design Steps of the Proposed Method

Given the desired passband and stopband deviations and , the design steps of the proposed method can be summarized as follows

(1)Find , , , and using (7). Then, find the values and with (15) and estimate the degree of the sharpening polynomial using (16).(2)Obtain equally spaced points and , for , over the regions and , respectively, assigning , , , and , with given in (3). Then, set and , with given in (9).(3)Choose the desired word-length and the desired values and , for .(4)Create , and using (20)–(23). Then, solve the problem (18) for . A straightforward way is using the MATLAB routine available online at [28].(5)Obtain the sharpening coefficients using (24).

6. Sharpening the Second-Stage Filter in a Two-Stage Architecture

Earlier in this paper we pointed out that the two-stage comb-based structure, which can be formed when the decimation factor can be expressed as , is effective to balance area and power consumptions. When this structure is chosen, the second-stage comb filter must be carefully designed since this is the filter where the worst-case magnitude characteristic of the overall cascade does occur. Moreover, the first-stage comb filter introduces a passband droop that should be corrected as well. It is interesting to note that, with the proposed sharpening approach, we can obtain an overall magnitude response attaining desired passband and stopband deviations by improving only the second-stage filter. However, we must have monotonic magnitude characteristic over the passband region of the filter to be sharpened.

The transfer function of the proposed two-stage filter is

where , with given in (1) but replacing by , and , with given in (8) but replacing by , , and . The zero-phase frequency response is

where , with given in (2) but replacing by , and , with given in (9) but replacing by , , and .

Let us identify

To correct the passband droop of the first-stage comb filter, whose zero-phase frequency response is , the zero-phase frequency response of the second-stage sharpened filter, , must be designed to follow an amplitude given by over the frequency interval . Therefore, the desired values with must be chosen as

where is the th point of the equally spaced partition of (see Section 5). Since is monotonically decreasing over , all the desired values are different from each other. Thus, can approximate every value , with ), if for all , . To meet this condition, must be monotonic over . Finally, the entries and in (21) must be multiplied by . Similarly, the entries and must be multiplied by , where is the th point of the equally spaced partition of (see Section 5). This is done in order to achieve an equiripple passband deviation in the overall filter .

7. Examples and Discussion of Results

The following examples are discussed to show the improvement of magnitude characteristics of comb filters achieved with the proposed method in comparison to other sharpening-based schemes recently introduced in the literature.

Example 1 (see the example in Section 4 of [22]). Consider , , and . The goal is to attain at least −100 dB gain in the folding bands, with an additional passband improvement without any specific constraint.

Let us consider the following solutions.

(a)A 1st-order comb filter (), presharpened by the polynomial and then sharpened with the 5th-order degree, first kind, Chebyshev polynomial (solution using the presharpening approach introduced in Section 6 (a) of [22]). This filter is identified by .(b)A compensated comb filter with and , sharpened with the polynomial (solution using the method of [23]). Let us call this filter .(c), a compensated comb filter with and compensation parameter , sharpened with the polynomial (solution using the scheme proposed in this work with and ).

Figure 3 shows the magnitude response of these filters, along with detail in passband and the first folding band. Notice that the three filters attain the −100 dB requirement in the folding bands. In the passband, the behaviors of filters proposed in [22, 23] are similar. The proposed filter, on the other hand, achieves better passband droop correction, which meets the 0.01 dB ripple () specification.

Figure 3: Magnitude responses of the filters (method [22]), (method [23]), and (proposed method), presented in Example 1 with , , and .

When it comes to the complexities in terms of APOS, the proposed solution achieves better results too. The filter , implemented with a CIC-like structure, requires 20 integrators working at high rate, due to its double-sharpening scheme. Therefore, its APOS metric would be higher than . On the other hand, the APOS of is 211 (see [23] for calculation of APOS in such structure). In the proposed method, we substitute , , , , , , and in (11) obtaining an APOS of 154. These results are summarized in Table 1.

Table 1: Comparisons in terms of APOS for Examples 1 and 2.

Method
Example 1
Example 2

[22]
APOS > 360
—

[23]
APOS = 211
APOS = 223

Proposed
APOS = 154
APOS = 164

Example 2. Consider a two-stage decimation filter with , , , and . The goal is to attain an attenuation of 60 dB in the folding bands with an additional passband improvement (without any given constraint). Assuming that the first-stage filter is a comb filter with , let us consider the following solutions for the second-stage filter:

(a), a compensated comb filter with and , sharpened with the polynomial (solution using the method in [23]);(b), a compensated comb filter with and compensation parameter , sharpened with the polynomial (solution using the proposed scheme with and ).

Figure 4 shows the magnitude response characteristics of these filters along with passband and first folding band details. Clearly, the filter designed with the proposed method presents both improvements: (1) better magnitude characteristic and (2) lower complexity, as summarized in Table 1. For this example, the APOS in Table 1 corresponds to the second-stage filter (the first-stage filtering is the same in both solutions and therefore it is omitted).

Figure 4: Magnitude responses of the filters (method [23]) and (proposed method), presented in Example 2 with , , and .

8. Conclusion

This paper proposed an optimization framework to design sharpening polynomials specifically suited to comb-based decimation filtering. The goal of the optimization problem was to minimize the min-max error over the frequency bands of interest of the sharpened filter. The optimization problem can be solved straightforwardly using the MATLAB Optimization Toolbox. The sharpening coefficients are guaranteed to be integers scaled by power-of-2 terms, thus resulting in low-complexity structures. Moreover, it was shown that the use of compensated comb filters, instead of combs only as basic building blocks in the sharpened filter, results in lower complexity structures (in terms of Additions Per Output Sample) for the same magnitude characteristics. Finally, it was shown that the proposed method provides better magnitude characteristic than other sharpening-based approaches for two-stage comb-based structures since it is able to correct the passband droop introduced by the first-stage comb filter.

Conflict of Interests

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

Acknowledgments

This work was supported by CONACYT Mexico (Grant no. 179587) and by US National Science Foundation (Grant no. 0925080).

References

E. B. Hogenauer, “An economical class of digital filters for decimation and interpolation,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 29, no. 2, pp. 155–162, 1981.

L. L. Presti, “Efficient modified-sinc filters for sigma-delta A/D converters,” IEEE Transactions on Circuits and Systems II, vol. 47, no. 11, pp. 1204–1213, 2000.

M. Laddomada, “Generalized comb decimation filters for ∑Δ A/D converters: analysis and design,” IEEE Transactions on Circuits and Systems I, vol. 54, no. 5, pp. 994–1005, 2007.

M. Laddomada, “On the polyphase decomposition for design of generalized comb decimation filters,” IEEE Transactions on Circuits and Systems I, vol. 55, no. 8, pp. 2287–2299, 2008.

G. J. Dolecek, “A new modified comb-rotated sinc (RS) decimator with improved magnitude response,” in Proceedings of the 14th IEEE International Conference on Electronics, Circuits and Systems (ICECS '07), pp. 250–253, Marrakech, Morocco, December 2007.

G. J. Dolecek and M. Laddomada, “An economical class of droop-compensated generalized comb filters: analysis and design,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 57, no. 4, pp. 275–279, 2010.

G. J. Dolecek and M. Laddomada, “A novel two-stage nonrecursive architecture for the design of generalized comb filters,” Digital Signal Processing, vol. 22, pp. 859–868, 2012.

H. J. Oh, S. Kim, G. Choi, and Y. H. Lee, “On the use of interpolated second-order polynomials for efficient filter design in programmable downconversion,” IEEE Journal on Selected Areas in Communications, vol. 17, no. 4, pp. 551–560, 1999.

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.

S. Kim, W.-C. Lee, S. Ahn, and S. Choi, “Design of CIC roll-off compensation filter in a W-CDMA digital IF receiver,” Digital Signal Processing, vol. 16, no. 6, pp. 846–854, 2006.

G. Jovanovic Dolecek and S. K. Mitra, “Simple method for compensation of CIC decimation filter,” Electronics Letters, vol. 44, no. 19, pp. 1162–1163, 2008.

G. Jovanovic Dolecek, “Simple wideband CIC compensator,” Electronics Letters, vol. 45, no. 24, pp. 1270–1272, 2009.

G. Molnar and M. Vucic, “Closed-form design of cic compensators based on maximally flat error criterion,” IEEE Transactions on Circuits and Systems II, vol. 58, no. 12, pp. 926–930, 2011.

D. E. T. Romero and G. J. Dolecek, “Application of amplitude transformation for compensation of comb decimation filters,” Electronics Letters, vol. 49, no. 16, pp. 985–987, 2013.

G. J. Dolecek and S. K. Mitra, “Stepped triangular CIC-cosine decimation filter,” in Proceedings of the 7th Nordic Signal Processing Symposium (NORSIG '06), pp. 26–29, June 2006.

G. J. Dolecek, “Simplified rotated sinc (RS) filter for sigma-delta A/D converters,” in Proceedings of the 1st International Conference on Green Circuits and Systems (ICGCS '10), pp. 283–288, Shanghai, China, June 2010.

G. J. Dolecek and S. K. Mitra, “Two-stage CIC-based decimator with improved characteristics,” IET Signal Process, vol. 4, no. 1, pp. 22–29, 2010.

T. Saramaki and T. Ritoniemi, “Modified comb filter structure for decimation,” in Proceedings of the IEEE International Symposium on Circuits and Systems (ISCAS '97), pp. 2353–2356, June 1997.

A. Y. Kwentus, Z. Jiang, and A. N. Willson Jr., “Application of filter sharpening to cascaded integrator-comb decimation filters,” IEEE Transactions on Signal Processing, vol. 45, no. 2, pp. 457–467, 1997.

G. Stephen and R. W. Stewart, “High-speed sharpening of decimating CIC filter,” Electronics Letters, vol. 40, no. 21, pp. 1383–1384, 2004.

G. Jovanovic-Dolecek and S. K. Mitra, “A new two-stage sharpened comb decimator,” IEEE Transactions on Circuits and Systems I, vol. 52, no. 7, pp. 1414–1420, 2005.

J. O. Coleman, “Chebyshev stopbands for CIC decimation filters and CIC-implemented array tapers in 1D and 2D,” IEEE Transactions on Circuits and Systems I, vol. 59, no. 12, pp. 2956–2968, 2012.

G. J. Dolecek and F. Harris, “Design of wideband CIC compensator filter for a digital IF receiver,” Digital Signal Processing, vol. 19, no. 5, pp. 827–837, 2009.

G. J. Dolecek and F. Harris, “On design of two-stage CIC compensation filter,” in Proceedings of the IEEE International Symposium on Industrial Electronics (ISIE '09), pp. 903–908, July 2009.

L. Xiong, “A high speed digital decimation filter with parallel cascaded integrator-comb pre-filters,” in Proceedings of the 2nd International Congress on Image and Signal Processing (CISP '09), pp. 1–4, October 2009.

X. Liu and A. N. Willson Jr., “A 1.2 Gb/s recursive polyphase cascaded integrator-comb prefilter for high speed digital decimation filters in 0.18-μm CMOS,” in Proceedings of the IEEE International Symposium on Circuits and Systems (ISCAS '10), pp. 2115–2118, Los Angeles, Calif, USA, June 2010.

J. F. Kaiser and R. W. Hamming, “Sharpening the response of a symmetric non-recursive filter by multiple use of the same filter,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 25, no. 5, pp. 415–422, 1977.

S. Tawfik, “MATLAB function to solve MILP problems,” 2013, http://www.mathworks.com/matlabcentral/fileexchange/6990-mixed-integer-lp/content/IP.m

T. Saramäki, “Design of FIR filters as a tapped cascaded interconnection of identical subfilters,” IEEE Transactions on Circuits and Systems I, vol. 34, no. 9, pp. 1011–1029, 1987.

O. Herrmann, L. R. Rabiner, and D. S. K. Chan, “Practical design rules for optimum finite impulse response low pass digital filters,” Bell System Technical Journal, vol. 52, no. 6, pp. 769–799, 1973.

J. F. Kaiser, “Non-recursive digital filter design using I ₀-sinh window function,” in Proceedings of the IEEE International Symposium on Circuits and Systems, pp. 20–23, April 1974.

K. Ichige, M. Iwaki, and R. Ishii, “Accurate estimation of minimum filter length for optimum FIR digital filters,” IEEE Transactions on Circuits and Systems II, vol. 47, no. 10, pp. 1008–1016, 2000.

R. Lyons, “Sample rate conversion,” in Understanding Digital Signal Processing, Prentice Hall, Upper Saddle River, NJ, USA, 2nd edition, 2004.

C. Candan, “Optimal Sharpening of CIC filters and an efficient implementation through Saramaki-Ritoniemi decimation filter structure,” 2013, http://www.eee.metu.edu.tr/~ccandan/pub_dir/opt_sharpened_CIC_filt__extended__new.pdf