Minimax design of two-channel nonuniform-division FIR filter banks

The paper deals with the minimax design of two-channel nonuniform-division filter (NDF) banks. Based on a linearisation scheme, the design problem is formulated as an optimisation problem with linear constraints. The authors present a method to design a two-channel NDF bank using a modified dual-affine scaling variant of Karmarkar's (1984) algorithm. This method provides the optimal results that the linear-phase FIR analysis and synthesis filters have equiripple stopband response and the resulting NDF bank also shows equiripple reconstruction error behaviour. The effectiveness of the proposed design technique is demonstrated by several simulation examples     

Design of IIR digital filters in the complex domain by transforming the desired response

In this paper, we present a new design method of infinite impulse response (IIR) digital filters with quasi-equiripple absolute error in the complex domain. This method is based on solving a least squares solution iteratively. At each iteration, the desired response for the least squares approximation is transformed to have equiripple error. This algorithm is efficient because there is no need for any initial value or complex optimization algorithm. By this method, a quasi-equiripple solution is obtained very quickly with less computational complexity. Moreover, by multiplying an arbitrary weighting function on the desired responses of passband and stopband, respectively, the error at the passband and stopband can be controlled. Finally, we show some examples to validate the proposed method.     

A weighted least squares algorithm for quasi-equiripple FIR and IIRdigital filter design

It has been demonstrated by several authors that if a suitable frequency response weighting function is used in the design of a finite impulse response (FIR) filter, the weighted least squares solution is equiripple. The crux of the problem lies in the determination of the necessary least squares frequency response weighting function. A novel iterative algorithm for deriving the least squares frequency response weighting function which will produce a quasi-equiripple design is presented. The algorithm converges very rapidly. It typically produces a design which is only about 1 dB away from the minimax optimum solution in two iterations and converges to within 0.1 dB in six iterations. Convergence speed is independent of the order of the filter. It can be used to design filters with arbitrarily prescribed phase and amplitude response     

Weighted least squares near-equiripple approximation of variable fractional delay FIR filters

The weighted least squares (WLS) method is a well-known method for designing a finite impulse response (FIR) filter. And some authors have reported that if a suitable frequency response weighting function is used to design the filter, the WLS method can produce an equiripple result. However, the weighting function for minimax optimality of WLS design is hard to derive analytically. By an iterative method with an adjustable elaborately constructed weighting function, this idea is extended to design a near-equiripple variable fractional delay FIR filter. The proposed method is superior to the fixed-weighting WLS design in the peak absolute error by about 6.6874 dB. The algorithm converges very rapidly. From the simulation, it typically produces a design which is only about 1 dB away from the truly equiripple solution in two iterations and converges to within 0.0056 dB in eight iterations.     

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     

Quasi-equiripple approximation of minimum phase FIR filters by updating desired response

The authors present a numerical method for the Chebyshev approximation of minimum phase FIR digital filters. This method is based on solving a least squares (LS) problem iteratively. At each iteration, the desired response is transformed so as to have an equiripple magnitude error. This method makes it possible to design minimum phase FIR filters whose magnitude error is quasi-equiripple. Using this method, a quasi-equiripple solution is obtained very quickly. Since the proposed methods do not require any time-consuming optimisation procedure, they require less computational complexity than conventional methods. Finally, some examples to illustrate the advantage of the proposed methods are shown.     

Eigenfilter approach for the design of allpass filtersapproximating a given phase response

We apply the eigenfilter method to design an allpass filter that approximates a given phase response in the least-squares (LS) sense. As it is not possible to express the exact LS phase error as a quadratic form suitable for eigenfilter formulation, alternative error measures that approximate the ideal LS error are proposed. For each of these new formulations, the allpass coefficients are obtained as the elements of the eigenvector corresponding to the minimum eigenvalue of a real, symmetric, and positive definite matrix. We propose a fast-converging iterative technique to approximate the ideal LS phase error solution. By employing an iterative weighting technique, the phase error can he made approximately equiripple. The design methods are illustrated with various practical examples and the results are compared to allpass filters designs reported in the literature     

Constrained least squares design of 2-D FIR filters

We consider the design of 2-D linear phase finite impulse response (FIR) filters according to the least squares (LS) error criterion subject to equality and/or inequality constraints. Since we use a frequency domain formulation, these constraints can be used to explicitly prescribe (frequency-dependent) error tolerances, the maximum, minimum, or fixed values of the frequency response at certain points and/or regions. Our method combines Lagrange multiplier and Kuhn-Tucker theory to solve a much wider class of problems than do standard methods. It allows arbitrary compromises between the LS and the equiripple design     

Minimax design of two-channel low-delay perfect-reconstruction FIRfilter banks

The authors deal with the design problem of low-delay perfect-reconstruction filter banks for which the FIR analysis and synthesis filters have equiripple magnitude response. Based on the minimax error criterion, the design problem is formulated in such a manner that the coefficients for the FIR analysis filters can be found by minimising the weighted peak error of the designed analysis filters, subject to the perfect-reconstruction constraints. A design technique based on a modified dual-affine scaling variant of Karmarkar's (1989) algorithm, in conjunction with approximation schemes, is then developed for solving the resulting nonlinear optimisation problem. The effectiveness of the proposed design technique is demonstrated by several simulation examples     

数字滤波器的近似等波纹设计法

本文证明了数字滤波器的自适应设计法等价于MMSE设计法，并提出用幅度误差函数对代价函数进行自适应迭代的算法达到近似等滤纹数字滤波器的自适应设计法．实验证明该设计方法简单有效，适合用于滤波器的工程设计。     

Perfect-reconstruction filter banks having linear-phase FIR filterswith equiripple response

The design of equiripple linear-phase analysis and synthesis FIR filters of two-channel perfect-reconstruction (PR) filter banks is formulated as the minimization of a weighted peak-error under both linear inequality (arising from the desired responses of the analysis filters) and nonlinear equality (PR) constraints. The effectiveness of a proposed method to solve the design problem (a modified dual-affine scaling variant of Karmarkar's (1989) algorithm and an approximation scheme) is illustrated through several design examples     

Analytical design method for optimal equiripple comb FIR filters

A novel analytical design method for highly selective digital optimal equiripple comb finite-impulse response (FIR) filters is presented. The equiripple comb FIR filters are optimal in the Chebyshev sense. The number of notch bands, the width of the notch bands and the attenuation in the passbands can be independently specified. The degree formula and the differential equation for the generating polynomial of the filter is presented. Based on the differential equation, a fast simple algebraic recursive procedure for the evaluation of the impulse response of the filter is described. Its arithmetic robustness outperforms, by far, the known analytical design method. Highly selective equiripple comb FIR filters with thousands of coefficients can be designed. One example demonstrates the efficiency of the filter design.     

Optimal Design of All-Pass Variable Fractional-Delay Digital Filters

This paper presents a computational method for the optimal design of all-pass variable fractional-delay (VFD) filters aiming to minimize the squared error of the fractional group delay subject to a low level of squared error in the phase response. The constrained optimization problem thus formulated is converted to an unconstrained least-squares (LS) optimization problem which is highly nonlinear. However, it can be approximated by a linear LS optimization problem which in turn simply requires the solution of a linear system. The proposed method can efficiently minimize the total error energy of the fractional group delay while maintaining constraints on the level of the error energy of the phase response. To make the error distribution as flat as possible, a weighted LS (WLS) design method is also developed. An error weighting function is obtained according to the solution of the previous constrained LS design. The maximum peak error is then further reduced by an iterative updating of the error weighting function. Numerical examples are included in order to compare the performance of the filters designed using the proposed methods with those designed by several existing methods.     

一种设计对数FIR数字滤波器的方法

本文提出了一种设计具有等波纹对数幅度响应的线性相位ＦＩＲ数字滤波器的方法,该设计方法以多次交换算法为基础。在给定通带与阻带误差比、通带误差和阻带误差三种情况下讨论该设计方法。介绍几个低通对数ＦＩＲ滤波器的设计例子,来说明该设计方法的效率。     

Linear periodic systems and multirate filter design

We present a systematic procedure for the design of filters intended for multirate systems. This procedure Is motivated by viewing the equiripple design of filters in linear time-invariant systems as a process of obtaining optimum minimax filters for a class of bounded energy input signals. The philosophy of designing optimum minimax filters for classes of input signals is extended to multirate systems, which are not time-invariant. We develop a generalized Fourier analysis appropriate for linear periodic systems and use it to derive new error criteria for multirate filter design. Using such criteria yields optimum minimax multirate filters for the input signal class. The utility of our method is demonstrated by using it to analyze several multirate systems. We give numerical results on the design of a multirate implementation of a narrowband filter and compare our work to previous work on multirate filter design. Our numerical analysis is based upon a new formulation of the design as a semi-infinite linear programming problem     

最小平方误差受限的FIR滤波器设计

在设计FIR滤波器时,往往会指定过渡带大小,但过渡带的引入只是为了便于滤波器的设计,而并不是物理上的需要,所以在设计中仅需指定截至频率。这是第一个设计理念。此外,在FIR滤波器的设计中存在两种准则:一是等波纹设计准则 (即最大误差最小化或者Chebyshev准则 ),另一种是平方误差最小准则。但在现实中两种准则往往要同时兼顾,所以仅基于其中一种准则来设计不能得到最佳结果。这是第二个设计理念。基于上述两种设计理念,提出了一种新的FIR滤波器设计算法。该算法采用最陡梯度下降法来对平方误差最小化下的最佳滤波器系数进行迭代修正,得到最佳结果。     

Generalized WLS Method for Designing All-Pass Variable Fractional-Delay Digital Filters

This paper presents an optimal weighted least squares (WLS) method for designing low-complexity all-pass variable fractional-delay (VFD) digital filters. Instead of using a fixed range for the VFD parameter p and same-order constant-coefficient filters (subfilters), both the VFD parameter range p isin [p _Min,p _Max] and subfilter orders are optimized such that a low-complexity all-pass VFD filter can be achieved for the LS design. To suppress the peak errors of variable frequency response, weighting functions are adopted and optimized such that the boundary peak errors can be further reduced but without noticeably increasing the total error energy (integral of squared error) of variable frequency response. After optimizing the variable range of the VFD parameter, weighting functions, and subfilter orders, an all-pass VFD filter can be designed by using a generalized noniterative WLS method, which yields a closed-form solution. Design examples are given to illustrate that utilizing different-order subfilters, along with the optimal range and optimal weighting functions, can yield an all-pass VFD filter with significantly reduced complexity and design errors as compared with existing ones.     

FIR数字滤波器幅频响应约束最大加权相位误差最小化设计

多数信号滤波应用,对滤波器幅频响应的要求高于相频响应.本文研究了满足幅频响应约束的有限脉冲响应（Infinite Impulse Response,FIR）数字滤波器设计,提出了最大加权相位误差最小化方法.用凸的椭圆误差约束代替非凸的幅值误差约束,将设计问题转化为凸问题;通过与二分技术结合,提出了给定权函数的幅值误差约束最大加权相位误差最小化设计的求解算法.以此算法为核心,构建了迭代重加权最大加权相位误差最小化算法,其中的权函数不再固定,而是基于修改的群延迟误差包络线在迭代中不断更新.权函数收敛后,所得滤波器具有近似等纹波的群延迟误差,最大群延迟误差得到了有效减小.仿真实验表明,与现有相位误差约束最大幅值误差最小化方法相比,得到的FIR滤波器具有更小的最大相位误差和最大群延迟误差.     

Microwave Networks with Equiripple Delay Characteristics

The problem of finding a characteristic polynomial to provide equiripple delay characteristics in commensurate TEM microwave networks is solved numerically. These polynomials enable the synthesis of such filters in the form of cascaded transmission lines or in the form of stub filters to be undertaken. Results for the orders 2 through 6 are presented. The synthesis of constant delay all-pass networks can also be performed and the method for doing so is presented. The attenuation characteristic corresponding to the equiripple delay filter does not, of course, possess a sharp cutoff, since these filters are of the minimum phase type.     

Equiripple FIR filter design by the FFT algorithm

The fast Fourier transform (FFT) algorithm has been used in a variety of applications in signal and image processing. In this article, a simple procedure for designing finite-extent impulse response (FIR) discrete-time filters using the FFT algorithm is described. The zero-phase (or linear phase) FIR filter design problem is formulated to alternately satisfy the frequency domain constraints on the magnitude response bounds and time domain constraints on the impulse response support. The design scheme is iterative in which each iteration requires two FFT computations. The resultant filter is an equiripple approximation to the desired frequency response. The main advantage of the FFT-based design method is its implementational simplicity and versatility. Furthermore, the way the algorithm works is intuitive and any additional constraint can be incorporated in the iterations, as long as the convexity property of the overall operations is preserved. In one-dimensional cases, the most widely used equiripple FIR filter design algorithm is the Parks-McClellan algorithm (1972). This algorithm is based on linear programming, and it is computationally efficient. However, it cannot be generalized to higher dimensions. Extension of our design method to higher dimensions is straightforward. In this case two multidimensional FFT computations are needed in each iteration