FIR, Allpass, and IIR Variable Fractional Delay Digital Filter Design

This paper presents two-step design methodologies and performance analyses of finite-impulse response (FIR), allpass, and infinite-impulse response (IIR) variable fractional delay (VFD) digital filters. In the first step, a set of fractional delay (FD) filters are designed. In the second step, these FD filter coefficients are approximated by polynomial functions of FD. The FIR FD filter design problem is formulated in the peak-constrained weighted least-squares (PCWLS) sense and solved by the projected least-squares (PLS) algorithm. For the allpass and IIR FD filters, the design problem is nonconvex and a global solution is difficult to obtain. The allpass FD filters are directly designed as a linearly constrained quadratic programming problem and solved using the PLS algorithm. For IIR FD filters, the fixed denominator is obtained by model reduction of a time-domain average FIR filter. The remaining numerators of the IIR FD filters are designed by solving linear equations derived from the orthogonality principle. Analyses on the relative performances indicate that the IIR VFD filter with a low-order fixed denominator offers a combination of the following desirable properties including small number of denominator coefficients, lowest group delay, easily achievable stable design, avoidance of transients due to nonvariable denominator coefficients, and good overall magnitude and group delay performances especially for high passband cutoff frequency ( ges 0.9pi) . Filter examples covering three adjacent ranges of wideband cutoff frequencies [0.95, 0.925, 0.9], [0.875, 0.85, 0.825], and [0.8, 0.775, 0.75] are given to illustrate the design methodologies and the relative performances of the proposed methods.     

Tunable FIR and IIR Fractional-Delay Filter Design and Structure Based on Complex Cepstrum

A cepstrum-based approach is proposed to design finite- and infinite-impulse-response (IIR) fractional-delay (FD) filters. The maximal-flatness criteria on frequency responses are formulated as a system of linear equations to solve the truncated complex cepstrum. The closed-form solutions to cepstrum sequences can be derived. Moreover, it is very attractive that the resultant cepstrum coefficients are directly proportional to the desired FD. Under a fixed filter order, the set of normalized complex cepstra needs to be computed once and stored, and the specific set for an arbitrary FD is obtained by simply multiplying the stored set with the delay value. According to this observation, we also design two kinds of tunable filter structures consisting of several linear-phase filters, in which it is more flexible to obtain better performance by adding the extra substructure without modifying the present one. Moreover, the tunable FD is simply controlled by a single parameter, and the usage of linear-phase filters saves half of the multipliers, largely reducing the cost of hardware implementation. In addition, we obtain an IIR all-pass filter with a wider useful band than that based on Thiran's design.     

Design of stable IIR digital filter based on least P-power error criterion

In this paper, the least p-power error criterion is presented to design digital infinite impulse response (IIR) filters to have an arbitrarily prescribed frequency response. First, an iterative quadratic programming (QP) method is used to design a stable unconstrained one-dimensional IIR filter whose optimal filter coefficients are obtained by solving the QP problem in each iteration. Then, the proposed method is extended to design constrained IIR filters and two-dimensional IIR filters with a separable denominator polynomial. Finally, design examples of the low-pass filter are demonstrated to illustrate the effectiveness of the proposed iterative QP method.     

Closed-Form Design of Generalized Maxflat$R$-Regular FIR$M$th-Band Filters Using Waveform Moments

$M$th-band filters have found numerous applications in multirate signal processing systems, filter banks, and wavelets. In this paper, the design problem of generalized maxflat$R$-regular finite impulse response (FIR)$M$th-band filters with a specified integer group delay at$ omega =0 $is considered, and the closed-form expression for its impulse response is presented. The filter coefficients are directly derived by solving a linear system of Vandermonde equations that are obtained from the regularity condition of the maxflat$R$-regular FIR$M$th-band filters via the blockwise waveform moments. Differing from the conventional FIR$M$th-band filters with exactly linear phase responses, the generalized FIR$M$th-band filters proposed in this paper have an arbitrarily specified integer group delay at$ omega =0 $. Moreover, a new efficient implementation of the generalized maxflat$R$-regular FIR$M$th-band filters is proposed by making use of the relationship between the filter coefficients in the closed-form solution. Finally, several design examples are presented to demonstrate the effectiveness of the proposed FIR$M$th-band filters.     

基于加权总体最小二乘正则化方法的混合滤波器组最优化设计

模拟分析滤波器组的实现欠理想、系统噪声以及数字综合滤波器有效阶数实现所带来的系统误差均有可能造成混合滤波器组的设计出现解不稳定、无唯一解等病态问题,影响混合滤波器组的准确重构效果。本文首先给出了满足准确重构条件下,以综合滤波器组频域响应为求解变量的混合滤波器组线性求解模型。针对线性方程中系数矩阵以及目标向量受扰动误差影响特点,提出一种新的基于加权总体最小二乘正则化算法的IIR形式综合滤波器设计方法。算法以系统扰动误差最小化为目标函数,根据随机误差变量的二阶统计特性,采用加权总体最小二乘算法抑制滤波器实现误差以及随机噪声等扰动因素影响,使得到的综合滤波器组频域响应解的加权误差平方和最小化,并通过Tikhonov正则化方法优化病态情况下方程组解的稳定性。提出一种IIR类型的综合滤波器系数的求解算法,并利用正则化方法优化滤波器系数,提高系统稳定性。该方法可应用于过采样混合滤波器组的设计。仿真结果表明该算法的有效提高系统鲁棒性和改善重构性能。      

A new method for the design of IIR filters with flat magnitude response

This paper presents a new direct design of infinite-impulse response (IIR) filters with a flat magnitude response in both passband and stopband (Butterworth filters). The design specifications are passband and stopband frequencies and passband droop and stopband attenuation. The approach is based on an allpass filter with flatness at frequency points /spl omega/=0 and /spl omega/=/spl pi/. Depending on the parity of the IIR filter order, the allpass filter is either real or complex. However, in both cases, the resulting IIR filter is real.     

Minimax design of two-channel IIR QMF banks with arbitrary groupdelay

The paper deals with the minimax design of two-channel infinite impulse response (IIR) QMF banks with arbitrary group delay, for which the IIR analysis filters and the resulting filter bank possess the frequency response optimal in the minimax (L_∞) sense. Utilising a lattice structure for the denominators of the IIR analysis filters, a design technique is presented based on an approximation scheme and a weighted least-squares (WLS) algorithm, previously developed by one of the authors for solving the resulting design problem that is basically a nonlinear optimisation problem. During the design process, this technique finds the tap coefficients for the numerator and the reflection coefficients for the denominator of the prototype IIR analysis filter simultaneously. The stability of the designed prototype IIR analysis filter is ensured by incorporating an efficient stabilisation procedure to make all of the reflection coefficient values fall between -1 and +1. Computer simulations show the effectiveness of the proposed design technique     

Noise reduction in recursive digital filters using high-order errorfeedback

The problem of solving the optimal (minimum-noise) error feedback coefficients for recursive digital filters is addressed in the general high-order case. It is shown that when minimum noise variance at the filter output is required, the optimization problem leads to set of familiar Wiener-Hopf or Yule-Walker equations, demonstrating that the optimal error feedback can be interpreted as a special case of Wiener filtering. As an alternative to the optimal solution, the formulas for suboptimal error feedback with symmetric or antisymmetric coefficients are derived. In addition, the design of error feedback using power-of-two coefficients is discussed. The efficiency of high order error feedback is examined by test implementations of the set of standard filters. It is concluded that error feedback is a very powerful and versatile method for cutting down the quantization noise in any classical infinite impulse response (IIR) filter implemented as a cascade of second-order direct form sections. The high-order schemes are attractive for use with high-order direct form sections     

Design of 1-D Stable Variable Fractional Delay IIR Filters

In this brief, a two-stage approach for the design of 1-D stable variable fractional delay infinite-impulse response (IIR) digital filters is proposed. In the first stage, a set of fixed delay stable IIR filters are designed by minimizing a quadratic objective function, which is defined by integrating error criterion with IIR filter stability constraint condition. Then, the final design is determined by fitting each of the fixed delay filter coefficients as a 1-D polynomial. Two design examples are given to show the effectiveness of the proposed design method     

A new class of two-channel biorthogonal filter banks and waveletbases

Proposes a novel framework for a new class of two-channel biorthogonal filter banks. The framework covers two useful subclasses: i) causal stable IIR filter banks. ii) linear phase FIR filter banks. There exists a very efficient structurally perfect reconstruction implementation for such a class. Filter banks of high frequency selectivity can be achieved by using the proposed framework with low complexity. The properties of such a class are discussed in detail. The design of the analysis/synthesis systems reduces to the design of a single transfer function. Very simple design methods are given both for FIR and IIR cases. Zeros of arbitrary multiplicity at aliasing frequency can be easily imposed, for the purpose of generating wavelets with regularity property. In the IIR case, two new classes of IIR maximally flat filters different from Butterworth filters are introduced. The filter coefficients are given in closed form. The wavelet bases corresponding to the biorthogonal systems are generated. the authors also provide a novel mapping of the proposed 1-D framework into 2-D. The mapping preserves the following: i) perfect reconstruction; ii) stability in the IIR case; iii) linear phase in the FIR case; iv) zeros at aliasing frequency; v) frequency characteristic of the filters     

Design of two-channel linear-phase QMF banks based on real IIR all-pass filters

The design of two-channel linear-phase quadrature mirror filter (QMF) banks constructed by real infinite impulse response (IIR) digital all-pass filters is considered. The design problem is appropriately formulated to result in a simple optimisation problem. Using a variant of Karmarkar's algorithm, the optimisation problem can be efficiently solved through a frequency sampling and iterative approximation method to find the real coefficients for the IIR digital all-pass filters. The resulting two-channel QMF banks possess an approximately linear phase response without magnitude distortion. The effectiveness of the proposed technique is achieved by forming an appropriate Chebyshev approximation of the desired phase response and then finding its solution from a linear subspace in a few iterations. Finally, several simulation examples are presented for illustration and comparison     

Optimization-based design and implementation of multidimensional zero-phase IIR filters

This paper considers multidimensional infinite-impulse response (IIR) filters that are iteratively implemented. The focus is on zero-phase filters with symmetric polynomials in the numerator and denominator of the multivariable transfer function. A rigorous optimization-based design of the filter is considered. Transfer function magnitude specifications, convergence speed requirements for the iterative implementation, and spatial decay of the filter impulse response (which defines the boundary condition influence in the spatial domain of the filtered signal) are all formulated as optimization constraints. When the denominator of the zero-phase IIR filter is strictly positive, these frequency domain specifications can be cast as a linear program and then efficiently solved. The method is illustrated with two two-dimensional IIR filter design examples.     

Design and Factorization of Two-Channel Perfect Reconstruction Filter Banks with Causal-Stable IIR Filters

In this paper, new design and factorization methods of two-channel perfect reconstruction (PR) filter banks (FBs) with casual-stable IIR filters are introduced. The polyphase components of the analysis filters are assumed to have an identical denominator in order to simplify the PR condition. A modified model reduction is employed to derive a nearly PR causal-stable IIR FB as the initial guess to obtain a PR IIR FB from a PR FIR FB. To obtain high quality PR FIR FBs for carrying out model reduction, cosine-rolloff FIR filters are used as the initial guess to a nonlinear optimization software for solving to the PR solution. A factorization based on the lifting scheme is proposed to convert the IIR FB so obtained to a structurally PR system. The arithmetic complexity of this FB, after factorization, can be reduced asymptotically by a factor of two. Multiplier-less IIR FB can be obtained by replacing the lifting coefficients with the canonical signal digitals (CSD) or sum of powers of two (SOPOT) coefficients.     

General IIR optical filter design for WDM applications usingall-pass filters

A general design algorithm is presented for infinite impulse response (IIR) bandpass and arbitrary magnitude response filters that use optical all-pass filters as building blocks. Examples are given for an IIR multichannel frequency selector, an amplifier gain equalizer, a linear square-magnitude response, and a multi-level response. Major advantages are the efficiency of the IIR filter compared to finite impulse response (FIR) filters, the simplicity of the optical architecture, and its tolerance for loss. A reduced set of unique operating states is discussed for implementing a reconfigurable multichannel selection filter     

Design of Constrained Causal Stable IIR Filters Using a New Second- Order-Cone-Programming-Based Model-Reduction Technique

This brief proposes a new method for designing infinite-impulse response (IIR) filter with peak error and prescribed flatness constraints. It is based on the model reduction of a finite-impulse response function that satisfies the specification by extending a method previously proposed by Brandenstein. The proposed model-reduction method retains the denominator of the conventional techniques and formulates the optimal design of the numerator as a second-order cone programming problem. Therefore, linear and convex quadratic inequalities such as peak error constraints and prescribed number of zeros at the stopband for IIR filters can be imposed and solved optimally. Moreover, a method is proposed to express the denominator of the model-reduced IIR filter as a polynomial in integer power of z, which efficiently facilitates its polyphase implementation in multirate applications. Design examples show that the proposed method gives better performance, and more flexibility in incorporating a wide variety of constraints than conventional methods     

Blind adaptation of stable discrete-time IIR filters in state-space form

Blind deconvolution consists of extracting a source sequence and impulse response of a linear system from their convolution. In the presence of system zeros close to the unit circle, which give rise to very long impulse responses, infinite-impulse-response (IIR) adaptive structures are of use, whose adaptation should be carefully designed in order to guarantee stability. In this paper, we propose a blind-type discrete-time IIR adaptive filter structure realized in state-space form that, with a suitable parameterization of its coefficients, remains stable. The theory is first developed for a two-pole filter, whose numerical behavior is investigated via numerical experiments. The proposed structure/adaptation theory is then extended to a multipole structure realized as a cascade of two-pole filters. Computer-based experiments are proposed and discussed, which aim at illustrating the behavior of the filter cascade on several cases of study. The numerical results obtained show the proposed filters remain stable during adaptation and provide satisfactory deconvolution results.     

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     

Boundary implications for frequency response of interval FIR andIIR filters

It is shown that vertex implication results in parameter space apply to interval trigonometric polynomials. Subsequently, it is shown that the frequency responses of both interval FIR and IIR filters are bounded by the frequency responses of certain extreme filters. The results apply directly in the evaluation of properties of designed filters, especially because it is more realistic to bound the filter coefficients from above and below instead of determining those with infinite precision because of finite arithmetic effects. Illustrative examples are provided to show how the extreme filters might be easily derived in any specific interval FIR or IIR filter design problem     

Optimal design of two-channel nonuniform-division FIR filter bankswith -1, 0, and +1 coefficients

This paper deals with the optimal design of two-channel nonuniform-division filter (NDF) banks whose linear-phase FIR analysis and synthesis filters have coefficients constrained to -1, 0, and +1 only. Utilizing an approximation scheme and a weighted least squares algorithm, we present a method to design a two-channel NDF bank with continuous coefficients under each of two design criteria, namely, least-squares reconstruction error and stopband response for analysis filters and equiripple reconstruction error and least-squares stopband response for analysis filters. It is shown that the optimal filter coefficients can be obtained by solving only linear equations. In conjunction with the proposed filter structure, a method is then presented to obtain the desired design result with filter coefficients constrained to -1, 0, and +1 only. The effectiveness of the proposed design technique is demonstrated by several simulation examples     

Coefficients of maximally flat nonrecursive digital filters

A general method is presented for determining the coefficients of maximally flat nonrecursive digital filters. Recursive relations with integer coefficients are derived for these coefficients by explicit calculations imitating Crout's method for solving the linear equations obtained by setting the derivatives of the symmetric, linear phase transfer function to zero at w = 0 and w = π. For a given order of the filter, the amount of computation required in this method is uniformly the same for any combination of permissible degrees of flatness at w = 0 and w = π.