Pade table, continued fraction expansion, and perfectreconstruction filter banks

We investigate the relationships among the Pade table, continued fraction expansions and perfect reconstruction (PR) filter banks. We show how the Pade table can be utilized to develop a new lattice structure for general two-channel bi-orthogonal perfect reconstruction (PR) filter banks. This is achieved through characterization of all two-channel biorthogonal PR filter banks. The parameterization found using this method is unique for each filter bank. Similar to any other lattice structure, the PR property is achieved structurally and the quantization of the parameters of the lattice does not effect this property. Furthermore, we demonstrate that for a given filter, the set of all complementary filters can be uniquely specified by two parameters, namely, the end-to-end delay of the system and a scalar quantity. Finally, we investigate the convergence of the successive filters found through the proposed lattice structure and develop a sufficient condition for this convergence     

Design of FIR Hilbert transformers using prescribed subfilters and nested FT technique

In this article, a transformation (subfilter) for designing finite impulse response (FIR) Hilbert transformers (HTs) is proposed. With our approach, except simple search procedures, neither optimisation nor any filter design algorithm is needed to obtain the transformation. The proposed subfilter requires only two multipliers regardless of the subfilter order. For the frequency transformation design, the purpose of the subfilter is to provide a “rough” shape of the desired HT; the two coefficients of the subfilter can be implemented as the form of sum of power of two (SOPOT) with only a few bits, thus leading to a multiplierless realisation. Moreover, by applying the transformation on the subfilter again, a technique named as nested frequency transformation (nested FT) is introduced. This technique can further reduce the number of multipliers needed in the overall HT.     

Two-channel IIR QMF banks with approximately linear-phase analysisand synthesis filters

Perfect linear-phase two-channel QMF banks require the use of finite impulse response (FIR) analysis and synthesis filters. Although they are less expensive and yield superior stopband characteristics, perfect linear phase cannot be achieved with stable infinite impulse response (IIR) filters. Thus, IIR designs usually incorporate a postprocessing equalizer that is optimized to reduce the phase distortion of the entire filter bank. However, the analysis and synthesis filters of such an IIR filter bank are not linear phase. In this paper, a computationally simple method to obtain IIR analysis and synthesis filters that possess negligible phase distortion is presented. The method is based on first applying the balanced reduction procedure to obtain nearly allpass IIR polyphase components and then approximating these with perfect allpass IIR polyphase components. The resulting IIR designs already have only negligible phase distortion. However, if required, further improvement may be achieved through optimization of the filter parameters. For this purpose, a suitable objective function is presented. Bounds for the magnitude and phase errors of the designs are also derived. Design examples indicate that the derived IIR filter banks are more efficient in terms of computational complexity than the FIR prototypes and perfect reconstruction FIR filter banks. Although the PR FIR filter banks when implemented with the one-multiplier lattice structure and IIR filter banks are comparable in terms of computational complexity, the former is very sensitive to coefficient quantization effects     

Design of Biorthogonal Filter Banks Using a Multi-objective Genetic Algorithm for an Image Coding Scheme

In this paper, we present an optimization method based on a multi-objective Genetic Algorithm (GA) for the design of linear phase filter banks for an image coding scheme. To be effective, the filter banks should satisfy a number of desirable criteria related to such a scheme. Instead of imposing the entire PR condition as in conventional designs, we introduce flexibility in the design by relaxing the Perfect Reconstruction (PR) condition and defining a PR violation measure as an objective criterion to maintain near perfect reconstruction (N-PR) filter banks. Particularly in this work, the designed filter banks are near-orthogonal. This has been made possible by minimizing the deviation from the orthogonality in the optimization process. The optimization problem is formulated as a constrained multi-objective, and a modified Nondominated Sorting Genetic Algorithm NSGAII is proposed in this work to find the Pareto optimal solutions that achieve the best compromise between the different objective criteria. The experimental results show that the filter banks designed with the proposed method outperform significantly the 9/7 filter bank of JPEG2000 in most cases. Furthermore, the filter banks are near orthogonal. This is very helpful, especially where embedded coding is required.     

Design of multiplierless fir filters by multiple use of the same filter

A novel approach to the design of multiplierless filters, based on the ACF (amplitude change function), is discussed. The prototype filter chosen is a CCOS (the cascade of the cosine functions) which requires no multipliers and only some adders. The required filter specifications are met by multiple use of the same CCOS filter. Effects due to coefficient quantisation do not arise when using the new approach. No multipliers are required to implement this filter.<>     

Design of linear phase M-channel perfect reconstruction FIR filterbanks

We propose two approaches to design M channel nonparaunitary filter banks that satisfy perfect reconstruction (PR) and linear phase (LP) properties. In the first approach, the PR condition is imposed on only a high-pass filter. Although this method does not require nonlinear optimization, it has a demerit in that the order of a high-pass filter becomes high. In the second approach, two filters are optimized simultaneously using a Lagrange-Newton method. We can design PR filter banks that have the same length. The PR constraint is also formulated as a linear and nonlinear equation of the analysis filter coefficients. Finally, some design examples are included     

Orthogonal and Biorthogonal FIR Hexagonal Filter Banks With Sixfold Symmetry

Recently, hexagonal image processing has attracted attention. The hexagonal lattice has several advantages in comparison with the rectangular lattice, the conventionally used lattice for image sampling and processing. For example, a hexagonal lattice needs less sampling points; it has better consistent connectivity; it has higher symmetry; and its structure is plausible to human vision systems. The multiresolution analysis method has been used for hexagonal image processing. Since the hexagonal lattice has high degree of symmetry, it is desirable that the hexagonal filter banks designed for multiresolution hexagonal image processing also have high order of symmetry, which is pertinent to the symmetry structure of the hexagonal lattice. The orthogonal or prefect reconstruction (PR) hexagonal filter banks that are available in the literature have only threefold symmetry. In this paper, we investigate the construction of orthogonal and PR finite impulse response (FIR) hexagonal filter banks with sixfold symmetry. We obtain block structures of 7-size refinement (seven-channel two-dimensional) orthogonal and PR FIR hexagonal filter banks with sixfold rotational symmetry. $sqrt{7}$-refinement orthogonal and biorthogonal wavelets based on these block structures are constructed. In this paper, we also consider FIR hexagonal filter banks with axial (line) symmetry, and we present a block structure of FIR hexagonal filter banks with pseudo-sixfold axial symmetry.      

Analytical design of 3-D wavelet filter banks using themultivariate Bernstein polynomial

The design of 3-D multirate filter banks where the downsampling/upsampling is on the FCO (face centred orthorhombic) lattice is addressed. With such a sampling lattice, the ideal 3-D sub-band of the low-pass filter is of the TRO (truncated octahedron) shape. The transformation of variables has been shown previously to be an effective technique for designing M-D (multidimensional) filter banks. A design technique is presented for the transformation function using the multivariate Bernstein polynomial which provides a good approximation to the TRO sub-band shape. The method is analytically based and does not require any optimisation procedure. Closed form expressions are obtained for the filters of any order. Another advantage of this technique is that it yields filters with a flat frequency response at the aliasing frequency (ω₁, ω₂ , ω₃)=(π, π, π). This flatness is important for giving regular discrete wavelet transform systems     

Two-stage CIC-based decimator with improved characteristics

A simple two-stage multiplierless cascaded-integrator-comb (CIC)-based decimator is presented. The first stage is a cascaded CIC filter whereas the second stage is a cascaded CIC filter and a second-order multiplierless compensator. The proposed decimator can be realised without filtering at high input rate by making use of the polyphase decomposition of the comb filter in the first stage. The proposed filter exhibits high aliasing attenuation and a low passband droop. The design parameters are the decimation factors, M ₁ and M ₂, numbers of cascaded CIC filters L and K, and parameter b of the compensator.     

Design of Multiplierless Lattice QMF: Structure and Algorithm Development

In this brief, efficient multiplierless design of lattice quadrature mirror filter bank is presented. Previous work by the authors has shown that splitting each lattice stage into cascade of subrotations results in larger stopband attenuation of filter than the conventional direct quantization. This brief extends the work further by exploiting the subrotations which yield more flexible sum of signed powers-of-two quantization. This enables us to find more possible discrete representations, and hence to reduce the quantization error. Also, an algorithm for the efficient gathering of candidate discrete coefficients is developed, based on the trellis-based searching approach. It substantially alleviates the overheads of optimization program, especially when the wordlengths and the number of nonzero digits are large. Several design examples are provided to show that the proposed structure with the candidate gathering algorithm provides improved frequency response.     

A general formulation of modulated filter banks

This paper presents a general framework for maximally decimated modulated filter banks. The theory covers the known classes of cosine modulation and relates them to complex-modulated filter banks. The prototype filters have arbitrary lengths, and the overall delay of the filter bank is arbitrary, within fundamental limits. Necessary and sufficient conditions for perfect reconstruction (PR) are derived using the polyphase representation. It is shown that these PR conditions are identical for all types of modulation-modulation based on the discrete cosine transform (DCT), both DCT-III/DCT-IV and DCT-I/DCT-II, and modulation based on the modified discrete Fourier transform (MDFT). A quadratic-constrained design method for prototype filters yielding PR with arbitrary length and system delay is derived, and design examples are presented to illustrate the tradeoff between overall system delay and stopband attenuation (subchannelization)     

Digital Filter Design Using Computer Algebra Systems

Digital filter design can be performed very efficiently using modern computer tools. The drawback of the numeric-based tools is that they usually generate a tremendous amount of numeric data, and the user might easily lose insight into the phenomenon being investigated. The computer algebra systems successfully overcome some problems encountered in the traditional numeric-only approach. In this paper, we introduce an original approach to algorithm development and digital filter design using a computer algebra system. The main result of the paper is the development of an algorithm for an infinite impulse response (IIR) filter design that, theoretically, is impossible to be implemented using the traditional approach. We present a step-by-step procedure which includes derivations of closed-form expressions for (1) the transfer functions of the implemented digital filter which contains the algebraic loop; (2) the closed-form expression for computing the number of requested iteration steps; and (3) the error function representing the difference of the output sample values of the new filter and that of the conventional filter. We demonstrate how one can use some already-known multiplierless digital filter as a start-up filter to design a new digital filter whose passband edge frequency can be simply moved by using a single parameter. As a result, we obtain a multiplierless IIR filter, which belongs to the family of low-power digital filters where multipliers are replaced with a small number of adders and shifters.     

过采样精确重构余弦调制滤波器组的设计

本文首先推导出过采样滤波器组精确重构的条件,由于此时所需的约束条件数比临界采样时少,因而可以设计出频域衰减特性更好的滤波器用,然后提出了精确重构约束条件下原型概通江波器的一种新的设计方法,采用矢的二次型约束优化算法,该算法优化方便,收剑速度快,与其它方法相比,滤波器的阻带衰减大。     

Theory and design of two-parallelogram filter banks

It is well known that the analysis and synthesis filters of orthonormal DFT filter banks can not have good frequency selectivity. The reason for this is that each of the analysis and synthesis filters have only one passband. Such frequency stacking (or configuration) in general does not allow alias cancellation when the individual filters have good stopband attenuation. A frequency stacking of this nature is called nonpermissible and should be avoided if good filters are desired. In a usual M-channel filter bank with real-coefficient filters, the analysis and synthesis filters have two passbands. It can be shown that the configuration is permissible in this case. Many designs proposed in the past demonstrate that filter banks with such configurations can have perfect reconstruction and be good filters at the same time. We develop the two-parallelogram filter banks, which is the class of 2-D filter banks in which the supports of the analysis and synthesis filters consist of two parallelograms. The two-parallelogram filter banks are analyzed from a pictorial viewpoint by exploiting the concept of permissibility. Based on this analysis, we construct and design a special type of two-parallelogram filter banks, namely, cosine-modulated filter banks (CMFB). In two-parallelogram CMFB, the analysis and synthesis filters are cosine-modulated versions of a prototype that has a parallelogram support. Necessary and sufficient conditions for perfect reconstruction of two-parallelogram CMFB are derived     

一类冗余比为2的二维DFT调制滤波器组的设计

本文构建了一类冗余比为2的二维线性相位的双原型离散傅立叶变换（DFT）调制滤波器组。利用原型滤波器的多相位分解,推导出了该滤波器组的完全重构(PR)条件。基于该PR条件,我们将滤波器组的设计归结为一个关于原型滤波器的多相位分量的无约束优化问题。由于原型滤波器是线性相位的,多相位分量之间具有一定的关系,因此我们可以简化该优化问题。仿真结果验证了滤波器组PR条件的正确性。同时,仿真表明了优化算法的有效性,设计所得的滤波器组重构误差很小、频率特性较好,基本满足实际应用的需要。      

Eigenfilter Approach to the Design of One-Dimensional and Multidimensional Two-Channel Linear-Phase FIR Perfect Reconstruction Filter Banks

We present an eigenfilter-based approach for the design of two-channel linear-phase FIR perfect-reconstruction (PR) filter banks. This approach can be used to design 1-D two-channel filter banks, as well as multidimensional nonseparable two-channel filter banks. Our method consists of first designing the low-pass analysis filter. Given the low-pass analysis filter, the PR conditions can be expressed as a set of linear constraints on the complementary-synthesis low-pass filter. We design the complementary-synthesis filter by using the eigenfilter design method with linear constraints. We show that, by an appropriate choice of the length of the filters, we can ensure the existence of a solution to the constrained eigenfilter design problem for the complementary-synthesis filter. Thus, our approach gives an eigenfilter-based method of designing the complementary filter, given a “predesigned” analysis filter, with the filter lengths satisfying certain conditions. We present several design examples to demonstrate the effectiveness of the method.      

用余弦调制PR－FIR方法构造正交小波基

离散小波变换将离散时间信号分解为一系列不同分辨率下的离散近似信号和离散细节，紧支的正交规范小波与完全重构正交镜象滤波器（ＰＲ－ＱＭＦ）相对应。本文在“二带”正交小波基的构造条件下，利用余弦调制完全重构滤波器组的方法，实现了正交小波基的构造，计算模拟表明该方法非常简单、有效。     

Genetic approach to design of multiplierless FIR filters

The design of multiplierless FIR filters using genetic algorithms is presented. The proposed algorithm uses simple operators (reproduction, crossover, and mutation) to search through the discrete coefficient space of predefined power-of-two coefficients. This approach has proved to be highly effective and outperformed existing multiplierless FIR design techniques.<>     

The theory and design of arbitrary-length cosine-modulated filterbanks and wavelets, satisfying perfect reconstruction

It is well known that FIR filter banks that satisfy the perfect-reconstruction (PR) property can be obtained by cosine modulation of a linear-phase prototype filter of length N=2mM, where M is the number of channels. In this paper, we present a PR cosine-modulated filter bank where the length of the prototype filter is arbitrary. The design is formulated as a quadratic-constrained least-squares optimization problem, where the optimized parameters are the prototype filter coefficients. Additional regularity conditions are imposed on the filter bank to obtain the cosine-modulated orthonormal bases of compactly supported wavelets. Design examples are given     

Design of multiplierless elliptic IIR filters with a smallquantization error

We present a new technique for the design of multiplierless IIR elliptic filters. The multiplierless filter has all multiplication constants implemented with a small number of shifters and adders. The proposed technique is based on sensitivity analysis. An analytical expression for amplitude response sensitivity is derived for the filter structures consisting of two allpass subfilters in parallel. It is shown that the amplitude response sensitivity to some constant x can be expressed as a product of the filter reflectance function and the phase sensitivity of the allpass section that implements the constant. The closed-form expressions for the phase sensitivities of the first- and second-order allpass sections are also developed. It is shown in the paper that the (n+1)/2 most sensitive constants can be directly controlled by the transfer function parameters if the transfer function is derived by the bilinear transformation from an elliptic minimal Q-factors (EMQF) analog prototype. This way, (n+1)/2 multiplication constants can be implemented without quantization, leaving the filter characteristic strictly elliptic. This is achieved for a class of low-noise allpass sections and for the wave lattice digital filter as well. The quantization of the remaining (n-1)/2 less-sensitive constants is performed using the phase-tolerance scheme and phase-sensitivity functions. The proposed design technique is straight-forward and, consequently, very fast. The application is demonstrated on the examples of narrowband, wideband, and halfband filters