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 K decimating by a factor M, with z-transfer function and zero-phase frequency response, respectively, given as
(1)HK(z)=(1M1-z-M1-z-1)K,(2)HK(ω)=(1Msin(ωM/2)sin(ω/2))K.
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 Ωp,
(3)Ωp=[0,ωp],ωp=πD=πMυ,
where D=Mυ 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 Ωk defined as
(4)Ωk=[2πkM-ωp,2πkM+ωp],k=1,2,…,⌊M2⌋,
where ⌊x⌋ stands for the integer part of x.
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.
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.
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.
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.
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.
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 F(z) and F(ω) 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 Nth degree sharpening polynomial,
(5)P(x)=∑k=0Npkxk,
allows mapping the amplitude values x=F(ω) to new amplitude values y=P(x). The new values y=P(x) must approximate the desired values d=D(x) for x∈Xp∪Xs, where Xp=[xp,l,xp,u] and Xs=[xs,l,xs,u], with tolerances δp (for the approximation over the region Xp) and δs (for the approximation over the region Xs). In this way, the zero-phase frequency response of the sharpened filter achieves the desired values with a maximum absolute passband deviation δp over the range of ω where F(ω)∈Xp and a maximum absolute stopband deviation δs over the range of ω where F(ω)∈Xs. Usually, D(x)=1 for x∈Xp and D(x)=0 for x∈Xs. Thus, the sharpening polynomial must meet the following simultaneous conditions:
(6)D(x)-δp≤P(x)≤D(x)+δp,forx∈Xp=[xp,l,xp,u],D(x)-δs≤P(x)≤D(x)+δs,forx∈Xs=[xs,l,xs,u].
For the comb-based decimation filter, the range limits for Xp and Xs are
(7)xp,l=min{F(ω)}|ω∈Ωp,xp,u=max{F(ω)}|ω∈Ωp,xs,l=min{F(ω)}|ω∈Ω1,xs,u=max{F(ω)}|Ω1,
where Ωp and Ω1 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]:
(8)F(z)=HK(z)·{-2-(b+2)[1-(2b+2+2)z-M+z-2M]}.
The zero-phase frequency response is
(9)F(ω)=HK(ω)·(1+2-(b+1)-2-(b+1)cos(Mω)).
The work [23] has shown that sharpening the subfilter F(z) with the polynomial P(x)=2x-x2 (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, K, and the parameter b must be adjusted accordingly. The polynomial P(x)=2x-x2 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
(10)HSH(z)=∑k=0Nz-(N-k)[K(M-1)/2+M]pkFk(z),HSH(ω)=∑k=0NpkFk(ω).
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 P(x)=2x-x2. Note that K must be an even value to avoid fractional delays.
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
(11)Ac=N[K(M+1)+3]+Q-1+∑l=1QS(pkl),
where S(pkl) indicates the number of adders required to implement the sharpening coefficient pkl and Q 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
(12)A=N~K(M+1)+Q-1+∑l=1QS(pkl),
with N~ 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 NK[(M+1)+3] and N~K(M+1), respectively. Hence, K and N (or N~) should be chosen as smaller values as possible for any arbitrary decimation factor M. Both K and N (or N~) have the same impact on the APOS metric. However, K can only take even values, whereas N (or N~) can also be odd. Since decreasing K is therefore more convenient, we can set K=2 in advance. This choice leads us to use b=1 in agreement with [23].
Let us discard the computational complexity introduced by the sharpening coefficients in both (11) and (12) and assume K=2 for the reason discussed above. We will compare the terms NK[(M+1)+3] in (11) and N~K(M+1) in (12), assuming that N~=N+a, with a being an integer. With this setup, we have
(13)Ac≈A~c=N[K(M+1)+3]|K=2=2N(M+1)+3N,A≈A~=N~K(M+1)|K=2,N~=N+a=2N(M+1)+2a(M+1).
Note that A~c and A~ differ in the terms 3N versus 2a(M+1). Clearly, the proposed structure can have a lower computational complexity (i.e., A~c<A~) when a>0; that is, when N~>N. In that case, sharpened compensated comb filters become convenient when
(14)M>(3N2a)-1.
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 P(x) with any of the formulas from [30–32]. These formulas are expressible as f(δp,δs,φp,φs), where δp and δs are, respectively, the desired maximum absolute passband and stopband deviation of the sharpened filter, whereas φp and φs are given as [29]
(15)φp=cos-1(2xp,l-xp,u-xs,lxp,u-xs,l),φs=cos-1(2xs,u-xp,u-xs,lxp,u-xs,l),
with xp,l, xp,u, xs,l, and xs,u 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 N (or N~) as
(16)N≈fδp,δs,φp,φs2=[-20log10(δpδs)-13]π14.6(φs-φp).
Upon noticing that the shape of the magnitude response of comb filters changes very little with M [33], we set in advance M=16 and we estimate the degrees N and N~ using the Kaiser formula for some typical values of υ, δp, and δs; namely, υ = 2, 4, 6, and 8; δp=0.001 (≈0.01 dB) and δs=0.001 (60 dB), 0.0001 (80 dB) and 0.00001 (100 dB) (to estimate N~, the values xp,l, xp,u, xs,l, and xs,u must be obtained from (7) but first replacing H2(ω) from (2) instead of F(ω) in these equations. To estimate N, we use F(ω) from (9) in these equations, first replacing K=2 and b=1 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.
In sharpened compensated comb filters, a lower computational complexity is obtained if K=2. This is because both K and N have the same impact on the APOS metric. However, K can only take even values, whereas N can also be odd. Therefore, preserving a simple sharpening polynomial and improving the stopbands with the increase of K, as suggested in [23], do not guarantee a result with low computational complexity.
Upon comparing the sharpened comb and sharpened compensated comb filters using K=2, 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.
Estimated degree of the sharpening polynomial to sharpen comb filters (dashed line) and compensated comb filters (solid line), with δp=0.001 and δs=0.001 (a), δs=0.0001 (b), and δs=0.00001 (c).
5. The Optimization Framework to Design Sharpened Comb Filters
Now, we introduce the optimization framework to obtain the discrete coefficients of P(x) for which the maximum deviation of P(x) with respect to D(x), 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 x covering the ranges Xp and Xs.
Let us consider the following notation in order to formalize the optimization problem.
δp and δs are the desired passband and stopband deviations after sharpening.
m=10N is the number of points partitioning the frequency sets Xp and Xs; that is, the overall number of grid points in the region Xp∪Xs is 20N.
x~i,p and x~i,s are the ith points belonging to the sets Xp and Xs, respectively. To find these points, we divide the range of frequencies Ωp=[0,ωp] into m equally spaced points ω~i,p and the range of frequencies Ω1=[(2π/M)-ωp,(2π/M)+ωp] into m equally spaced points ω~i,s, with i=1,2,…,m. Then, we set x~i,p=F(ω~i,p) and x~i,s=F(ω~i,s).
di,p and di,s are the desired amplitudes of the polynomial P(x) at the points x=x~i,p and x=x~i,s, respectively. Usually, di,p=1 and di,s=0 for all i.
[M]i,j denotes the entry in the ith row and jth column of the underlined matrix M.
[v]i denotes the ith element of the underlined vector v.
B 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 2-B:
(17)pk=2-Bp~k,p~k∈{integers},0≤k≤N.
By this setup, the optimization problem proposed in this paper can be written as
(18)hhhlhhminsfTssubjecttoAs≤b,[s]i∈{integers},2≤i≤N+2,
where f and s are vectors of size (N+2)×1, A is a matrix of size 4m×(N+2) and b is a vector of size 4m×1. In addition, we have
(19)[s]i={δ;i=1,p~i-2;i=2,3,…,N+2,(20)[f]i={1;i=1,0;i=2,3,…,N+2,(21)[A]i,j={-1;1≤i≤2m,j=1,-δsδp;2m<i≤4m,j=1,(22)[A]i,j={2-B(x~i,p)j-2;1≤i≤m,2≤j≤N+2,-2-B(x~i-m,p)j-2;m<i≤2m,2≤j≤N+2,2-B(x~i-2m,s)j-2;2m<i≤3m,2≤j≤N+2,-2-B(x~i-3m,s)j-2;3m<i≤4m,2≤j≤N+2,(23)[b]i={di,p;1≤i≤m,-di-m,p;m<i≤2m,di-2m,s;2m<i≤3m,-di-3m,s;3m<i≤4m.
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 s has been obtained, the sharpening coefficients can be found as follows:
(24)pk=2-B[s]k+2,0≤k≤N.
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 δp and δs, the design steps of the proposed method can be summarized as follows
Find xp,l, xp,u, xs,l, and xs,u using (7). Then, find the values φp and φs with (15) and estimate the degree N of the sharpening polynomial using (16).
Obtain m=10N equally spaced points ω~i,p and ω~i,s, for i=1,2,…,m, over the regions Ωp=[0,ωp] and Ω1=[(2π/M)-ωp,(2π/M)+ωp], respectively, assigning ω~1,p=0, ω~m,p=ωp, ω~1,s=(2π/M)-ωp, and ω~m,s=(2π/M)+ωp, with ωp given in (3). Then, set x~i,p=F(ω~i,p) and x~i,s=F(ω~i,s), with F(ω) given in (9).
Choose the desired word-length B and the desired values di,p and di,s, for i=1,2,…,m.
Create f, A and b using (20)–(23). Then, solve the problem (18) for s. A straightforward way is using the MATLAB routine available online at [28].
Obtain the sharpening coefficients pk 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 M can be expressed as M=M1M2, 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
(25)HTS(z)=H1K1(z)·[∑k=0Nz-M1(N-k)pkF2k(zM1)],
where H1(z)=H(z), with H(z) given in (1) but replacing M by M1, and F2(z)=F(z), with F(z) given in (8) but replacing M by M2, K=2, and b=1. The zero-phase frequency response is
(26)HT(ω)=H1K1(ω)·∑k=0NpkF2k(M1ω),
where H1(ω)=H(ω), with H(ω) given in (2) but replacing M by M1, and F2(ω)=F(ω), with F(ω) given in (9) but replacing M by M2, K=2, and b=1.
Let us identify
(27)P(F2(M1ω))=∑k=0NpkF2k(M1ω).
To correct the passband droop of the first-stage comb filter, whose zero-phase frequency response is H1K1(ω), the zero-phase frequency response of the second-stage sharpened filter, P(F2(M1ω)), must be designed to follow an amplitude given by 1/H1K1(ω) over the frequency interval Ωp=[0,ωp]. Therefore, the desired values di,p with i=1,2,…,m must be chosen as
(28)di,p=1H1K1(ω~i,p),
where ω~i,p is the ith point of the equally spaced partition of Ωp (see Section 5). Since H1(ω) is monotonically decreasing over Ωp, all the desired values di,p are different from each other. Thus, P(x~i,p) can approximate every value di,p, with x~i,p=F2(M1ω~i,p), if x~i,p≠x~j,p for all i≠j, i,j=1,2,…,m. To meet this condition, F2(M1ω) must be monotonic over Ωp. Finally, the entries [A]i,1 and [A]i+m,1 in (21) must be multiplied by di,p. Similarly, the entries [A]i+2m,1 and [A]i+3m,1 must be multiplied by 1/H1K1(ω~i,s), where ω~i,s is the ith point of the equally spaced partition of Ω1=[(2π/M)-ωp,(2π/M)+ωp] (see Section 5). This is done in order to achieve an equiripple passband deviation in the overall filter HTS(z).
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 M=16, υ=4, and ωp=0.907π/(Mυ). 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 1st-order comb filter (K=1), presharpened by the polynomial P1(x)=2x2-x4 and then sharpened with the 5th-order degree, first kind, Chebyshev polynomial P(x)=5x-20x3+16x5 (solution using the presharpening approach introduced in Section 6 (a) of [22]). This filter is identified by Ha(z).
A compensated comb filter with K=6 and b=-1, sharpened with the polynomial P(x)=2x-x2 (solution using the method of [23]). Let us call this filter Hb(z).
Hc(z), a compensated comb filter with K=2 and compensation parameter b=1, sharpened with the polynomial P(x)=2-5(-3x2+131x3-96x4) (solution using the scheme proposed in this work with δp=0.0006 and δs=0.000032).
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 (δp=0.0006) specification.
Magnitude responses of the filters Ha(z) (method [22]), Hb(z) (method [23]), and Hc(z) (proposed method), presented in Example 1 with M=16, υ=4, and ωp=0.907π/(Mυ).
When it comes to the complexities in terms of APOS, the proposed solution achieves better results too. The filter Ha(z), 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 320=20×16. On the other hand, the APOS of Hb(z) is 211 (see [23] for calculation of APOS in such structure). In the proposed method, we substitute M=16, K=2, N=4, Q=3, S(p2=3)=1, S(p3=131)=2, and S(p4=96)=1 in (11) obtaining an APOS of 154. These results are summarized in 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 M1=8, M2=17, υ=2, and ωp=0.9π/(M1M2υ). 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 K1=4, let us consider the following solutions for the second-stage filter:
Ha(z), a compensated comb filter with K=6 and b=-1, sharpened with the polynomial P(x)=2x-x2 (solution using the method in [23]);
Hb(z), a compensated comb filter with K=2 and compensation parameter b=1, sharpened with the polynomial P(x)=2-5(-x+5x2+116x3-88x4) (solution using the proposed scheme with δp=0.0006 and δs=0.001).
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).
Magnitude responses of the filters Ha(z) (method [23]) and Hb(z) (proposed method), presented in Example 2 with M=M1M2=8*17=136, υ=2, and ωp=0.9π/(Mυ).
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).
HogenauerE. B.An economical class of digital filters for decimation and interpolation19812921551622-s2.0-0019558332PrestiL. L.Efficient modified-sinc filters for sigma-delta A/D converters20004711120412132-s2.0-003431813510.1109/82.885128LaddomadaM.Generalized comb decimation filters for ∑Δ A/D converters: analysis and design200754599410052-s2.0-3424864116710.1109/TCSI.2007.895528LaddomadaM.On the polyphase decomposition for design of generalized comb decimation filters2008558228722992-s2.0-5474915710410.1109/TCSI.2008.920136DolecekG. J.A new modified comb-rotated sinc (RS) decimator with improved magnitude responseProceedings of the 14th IEEE International Conference on Electronics, Circuits and Systems (ICECS '07)December 2007Marrakech, Morocco2502532-s2.0-5064909967710.1109/ICECS.2007.4510977DolecekG. J.LaddomadaM.An economical class of droop-compensated generalized comb filters: analysis and design20105742752792-s2.0-7795169864910.1109/TCSII.2010.2043469DolecekG. J.LaddomadaM.A novel two-stage nonrecursive architecture for the design of generalized comb filters2012228598682-s2.0-8485988492210.1016/j.dsp.2012.04.007OhH. J.KimS.ChoiG.LeeY. H.On the use of interpolated second-order polynomials for efficient filter design in programmable downconversion19991745515602-s2.0-003263285110.1109/49.761035YeungK. S.ChanS. C.The design and multiplier-less realization of software radio receivers with reduced system delay20045112244424592-s2.0-1094426267310.1109/TCSI.2004.838253KimS.LeeW.-C.AhnS.ChoiS.Design of CIC roll-off compensation filter in a W-CDMA digital IF receiver20061668468542-s2.0-3375116442210.1016/j.dsp.2006.06.003Jovanovic DolecekG.MitraS. K.Simple method for compensation of CIC decimation filter20084419116211632-s2.0-5154910711910.1049/el:20081603Jovanovic DolecekG.Simple wideband CIC compensator20094524127012722-s2.0-7054909139810.1049/el.2009.1860MolnarG.VucicM.Closed-form design of cic compensators based on maximally flat error criterion201158129269302-s2.0-8365518274310.1109/TCSII.2011.2172522RomeroD. E. T.DolecekG. J.Application of amplitude transformation for compensation of comb decimation filters20134916985987DolecekG. J.MitraS. K.Stepped triangular CIC-cosine decimation filterProceedings of the 7th Nordic Signal Processing Symposium (NORSIG '06)June 200626292-s2.0-3904914901710.1109/NORSIG.2006.275267DolecekG. J.Simplified rotated sinc (RS) filter for sigma-delta A/D convertersProceedings of the 1st International Conference on Green Circuits and Systems (ICGCS '10)June 2010Shanghai, China2832882-s2.0-7795657160310.1109/ICGCS.2010.5543051DolecekG. J.MitraS. K.Two-stage CIC-based decimator with improved characteristics2010412229SaramakiT.RitoniemiT.Modified comb filter structure for decimationProceedings of the IEEE International Symposium on Circuits and Systems (ISCAS '97)June 1997235323562-s2.0-0030652102KwentusA. Y.JiangZ.WillsonA. N.Jr.Application of filter sharpening to cascaded integrator-comb decimation filters19974524574672-s2.0-0031079365StephenG.StewartR. W.High-speed sharpening of decimating CIC filter20044021138313842-s2.0-754422293110.1049/el:20046387Jovanovic-DolecekG.MitraS. K.A new two-stage sharpened comb decimator2005527141414202-s2.0-2314443790610.1109/TCSI.2005.851390ColemanJ. O.Chebyshev stopbands for CIC decimation filters and CIC-implemented array tapers in 1D and 2D2012591229562968DolecekG. J.HarrisF.Design of wideband CIC compensator filter for a digital IF receiver20091958278372-s2.0-6734927745710.1016/j.dsp.2009.03.013DolecekG. J.HarrisF.On design of two-stage CIC compensation filterProceedings of the IEEE International Symposium on Industrial Electronics (ISIE '09)July 20099039082-s2.0-7795013441610.1109/ISIE.2009.5213597XiongL.A high speed digital decimation filter with parallel cascaded integrator-comb pre-filtersProceedings of the 2nd International Congress on Image and Signal Processing (CISP '09)October 2009142-s2.0-7384914046610.1109/CISP.2009.5301187LiuX.WillsonA. N.Jr.A 1.2 Gb/s recursive polyphase cascaded integrator-comb prefilter for high speed digital decimation filters in 0.18-μm CMOSProceedings of the IEEE International Symposium on Circuits and Systems (ISCAS '10)June 2010Los Angeles, Calif, USA211521182-s2.0-7795600135810.1109/ISCAS.2010.5537065KaiserJ. F.HammingR. W.Sharpening the response of a symmetric non-recursive filter by multiple use of the same filter19772554154222-s2.0-0017548514TawfikS.MATLAB function to solve MILP problems2013, http://www.mathworks.com/matlabcentral/fileexchange/6990-mixed-integer-lp/content/IP.mSaramäkiT.Design of FIR filters as a tapped cascaded interconnection of identical subfilters198734910111029HerrmannO.RabinerL. R.ChanD. S. K.Practical design rules for optimum finite impulse response low pass digital filters19735267697992-s2.0-0015648304KaiserJ. F.Non-recursive digital filter design using I0-sinh window functionProceedings of the IEEE International Symposium on Circuits and SystemsApril 197420232-s2.0-0015984922IchigeK.IwakiM.IshiiR.Accurate estimation of minimum filter length for optimum FIR digital filters20004710100810162-s2.0-003429621710.1109/82.877143LyonsR.Sample rate conversion20042ndUpper Saddle River, NJ, USAPrentice HallCandanC.Optimal Sharpening of CIC filters and an efficient implementation through Saramaki-Ritoniemi decimation filter structure2013, http://www.eee.metu.edu.tr/~ccandan/pub_dir/opt_sharpened_CIC_filt__extended__new.pdf