PLS Algorithms for Equality-Constrained Linear-Phase FIR Filter Design

The least-squares design of linear-phase finite impulse response (FIR) filters with equality constraints in the time and frequency domains can be formulated as an equality-constrained quadratic programming (QP) problem. This paper presents an effective and robust novel algorithm, the projected least-squares (PLS) algorithm, for equality-constrained QP problems. The algorithm eventually projects an unconstrained minimization solution successively onto the feasible hyperplane of the problem. An additional term is added into the Hessian matrix of the cost function, thus ensuring the positive definiteness of the Hessian matrices during the iterative procedure for the design of constrained minimax filters. To demonstrate its effectiveness and robustness, the PLS algorithm is applied to design optimal linear-phase Nyquist filters that may satisfy some frequency domain constraints besides the zero-crossing impulse response.     

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.     

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     

A fast design algorithm for elliptic-error and phase-error constrained LS 2-D FIR filters

Two-dimensional (2-D) nonlinear-phase finite impulse response (FIR) filters have found many applications in signal processing and communication systems. This paper considers the elliptic-error and phase-error constrained least-squares design of 2-D nonlinear-phase FIR filters, and develops a matrix-based algorithm to solve the design problem directly for the filter’s coefficient matrix rather than vectorizing it first as in the conventional methods. The matrix-based algorithm makes the design to consume much less design time than existing algorithms. Design examples and comparisons with existing methods demonstrate the effectiveness and high efficiency of the proposed design method.     

等式约束FIR滤波器设计的投影最小二乘算法

本文考虑具有频域和时域等式约束的FIR滤波器设计问题,提出一个非常有效的新算法——投影最小二乘算法.该算法由两部分组成,前一部分产生一个解析的最小二乘解,后一部分将此解逐次投影到每个等式约束上.该算法有两个显著特点:一是目标函数的Hessian矩阵不要求正定;二是由于采用平方根因子分解来计算增广Hessian矩阵及投影算子矩阵,算法具有很好的数字稳定性.以此算法为核心构成了一个迭代算法,用于实现FIR Nyquist滤波器的minimax设计.设计例子表明了所提算法的有效性和数字稳定性.     

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.     

Design of two-dimensional FIR digital filters by McClellantransform and quadratic programming

A quadratic programming (QP) approach for determining the coefficients of the McClellan transform is presented for the design of 2-D FIR digital filters. Three features of the proposed method are as follows. First, the transform parameters are determined by minimising the integration of the squared errors along the desired contour. Second, a set of linear constraints are incorporated into the QP formulation such that the conventional scaling problem of the transform can be avoided. Third, the optimal cutoff frequencies of a 1-D prototype filter are obtained directly from the QP solution. Several design examples, including fan filters, elliptic filters, diamond filters and bandpass filters, are illustrated to demonstrate the effectiveness of the QP method     

Design of linear phase FIR filters using fractional derivative constraints

In this paper, the designs of linear phase FIR filters using fractional derivative constraints are investigated. First, the definition of fractional derivative is reviewed briefly. Then, the linear phase FIR filters are designed by minimizing integral squares error under the constraint that the ideal response and actual response have several same fractional derivatives at the prescribed frequency point. Next, the fractional maximally flat FIR filters are designed by letting the number of fractional derivative constraints be equal to the number of filter coefficients. Finally, numerical examples are demonstrated to show that the proposed method has larger design flexibility than the conventional integer derivative constrained methods.     

Constrained eigenfilter design without specified transition bands

The design of a finite-impulse response (FIR) filter with constraints in the frequency domain and/or time domain is considered. We further consider the design specification without explicitly specified transition band bandedges. An iterative algorithm without transition band specification is proposed to design FIR filters with various design constraints. We suggest the possible design tradeoff between transition band bandwidth and the ripple size of the filter. The proposed algorithm can be used to design filters with an optimal tradeoff from the design specification. The eigenfilter formulation further allows the propose algorithm to incorporate time-domain constraints simultaneously. Various design examples are presented to illustrate the versatility of the digital filter obtained by the propose algorithm. Although we have not proven the convergence of the proposed algorithm, it is found to converge efficiently in all the simulations.     

二维线性相位FIR滤波器设计的投影最小二乘算法

考虑二维线性相位矩形对称FIR滤波器的约束最小二乘设计问题,即在通带和阻带逼近误差不超过给定值的约束下使逼近误差平方和最小.提出一个投影最小二乘算法,它是一个交替地更新有效约束集及将二次误差无约束极小点(最小二乘解)逐次投影到有效约束边界的迭代过程.通过二维FIR低通圆形滤波器和方形滤波器的设计例子,对算法的性能进行了仿真,并与基于内点算法和有效集方法的设计程序进行了比较,结果表明本文算法具有很高的效率.     

多通带低群延迟误差FIR滤波器设计的迭代算法

常数低群延迟有限冲击响应（FIR）滤波器在通信等领域得到了广泛应用,尤其是要求无波形失真、信号延迟小的场合。而低群延迟的FIR滤波器,其相位响应只能做到近似线性相位,其群延迟只能做到近似常数。为了减小与期望常数群延迟之间的误差,最近提出的通过迭代更新相位误差上界函数来逐步减小群延迟误差的方法,只考虑了单通带滤波器的minimax设计。本文将把该方法推广至多通带FIR滤波器的minimax设计和约束最小二乘设计,先对各通带单独处理使每个通带的最大群延迟误差有效降低后,再考虑各通带之间平衡,对各子带的最大群延迟误差进行折中,进而使整个通带上的最大群延迟误差继续减小。对约束最小二乘设计还特别考虑了通过修改收敛参数来解决相位误差约束过紧时设计问题无解的问题。仿真实例表明,该方法能有效减小多通带滤波器的最大群延迟误差。      

Constrained Chebyshev design of FIR filters

In many filter-design problems, additional constraints are often imposed on the optimal filter in the sense of, say, minimal Chebyshev error norm. Based on the characteristic properties of the optimal filter for the Chebyshev design with frequency equation constraints, a modified Remez (MRemez) algorithm is proposed in this paper. The central problem of this paper is the constrained Chebyshev design of finite-impulse response filters with equation and inequality constraints in the frequency domain. By converting the problem into a series of Chebyshev design problems with equation constraints, an iterative MRemez algorithm which uses the MRemez algorithm as the computational core of the iteration is proposed, and the convergence of the algorithm is obtained. Design examples demonstrate the effectiveness and the fast convergence of the MRemez algorithm and the iterative MRemez algorithm.     

Optimal Design of Nonlinear-Phase FIR Filters With Prescribed Phase Error

Constrained least-squares design and constrained Chebyshev design of one- and two-dimensional nonlinear-phase FIR filters with prescribed phase error are considered in this paper by a unified semi-infinite positive-definite quadratic programming approach. In order to obtain unique optimal solutions, we propose to impose constraints on the complex approximation error and the phase error. By introducing a sigmoid phase-error constraint bound function, the group-delay error can be greatly reduced. A Goldfarb–Idnani based algorithm is presented to solve the semi-infinite positive-definite quadratic program resulting from the constrained least-squares design problem, and then applied after some modifications to the constrained Chebyshev design problem, which is proved in this paper to be equivalent also to a semi-infinite positive-definite quadratic program. Through design examples, the proposed method is compared with several existing methods. Simulation results demonstrate the effectiveness and efficiency of the proposed method.      

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     

Optical half-band filters

This paper proposes two kinds of novel 2×2 circuit configuration for finite-impulse response (FIR) half-band filters. These configurations can be transformed into each other by a symmetric transformation and their power transmittance is identical. The configurations have only about half the elements of conventional FIR lattice-form filters. We derive a design algorithm for achieving desired power transmittance spectra. We also describe 2×2 circuit configurations for infinite-impulse response (IIR) half-band filters. These configurations are designed to realize arbitrary-order IIR half-band filter characteristics by extending the conventional half-band circuit configuration used in millimeter-wave devices. We discuss their filter characteristics and confirm that they have a power half-band property. We demonstrate design examples including FIR maximally flat half-band filters, an FIR Chebyshev half-band filter, and an IIR elliptic half-band filter     

Weighted least-squares design and characterization of complex FIRfilters

The paper presents two novel weighted least-squares methods for the design of complex coefficient finite impulse response (FIR) filters to attain specified arbitrary multiband magnitude and linear or arbitrary phase responses. These methods are computationally efficient, requiring only the solution of a Toeplitz system of N linear equations for an N-length filter that can be obtained in o(N²) operations. Illustrative filter design examples are presented     

Digital Laguerre filter design with maximum passband-to-stopband energy ratio subject to peak and Group delay constraints

In this paper, we formulate a general design of transversal filter structures with maximum relative passband-to-stopband energy ratio subject to complex frequency response constraints in the passband and the stopband as well as additional constraints such as constraints. These constraints are important for applications where the suppression of noise at certain frequencies are important. Additional constraints are introduced allowing approximately linear phase and constant group delay in the passband. For a given set of basis functions, the design problem can be formulated as a semi-infinite quadratic optimization problem in the filter coefficients, which are the decision variables to be optimized. In this paper, we focus on the design of digital Laguerre filter and digital finite impulse response (FIR) filter structures. A modified bridging algorithm is developed for searching for the optimum pole of the Laguerre filters. Design examples are given to demonstrate the effectiveness of the proposed algorithm.     

Design of FIR filters with extra constraints using modified LMSalgorithm

An adaptive approach to the design of linear phase low-pass FIR filters with extra constraints on filter coefficients is presented. In this approach, the procedures using the LMS adaptive algorithm are modified to include the constraints on the filter coefficients. Numerical examples are presented and compared to the results obtained using least-square design in the frequency domain in which the filter design problem is transformed into an equivalent nonlinear optimisation problem     

A modified algorithm for constrained least square design ofmultiband FIR filters without specified transition bands

In a previous paper, we described a constrained least square approach to FIR filter design that does not use “don't care” regions. In that paper, we described a simple algorithm for the design of lowpass filters according to that approach. In this paper, we describe a modification of that algorithm that makes it converge for many multiband filter designs. Although no proof of convergence is given, the modified algorithm remains simple and converges rapidly in many cases. In this approach, the user supplies a lower and upper bound constraint that is exactly satisfied by the local minima and maxima of the frequency response amplitude. Yet, the constraints can be made as tight as desired-the transition band automatically adjusts (widens) to accommodate the constraints     

Theory and design of signal-adapted FIR paraunitary filter banks

We study the design of signal-adapted FIR paraunitary filter banks, using energy compaction as the adaptation criterion. We present some important properties that globally optimal solutions to this optimization problem satisfy. In particular, we show that the optimal filters in the first channel of the filter bank are spectral factors of the solution to a linear semi-infinite programming (SIP) problem. The remaining filters are related to the first through a matrix eigenvector decomposition. We discuss uniqueness and sensitivity issues. The SIP problem is solved using a discretization method and a standard simplex algorithm. We also show how regularity constraints may be incorporated into the design problem to obtain globally optimal (in the energy compaction sense) filter banks with specified regularity. We also consider a problem in which the polyphase matrix implementation of the filter bank is constrained to be DCT based. Such constraints may also be incorporated into our optimization algorithm; therefore, we are able to obtain globally optimal filter banks subject to regularity and/or computational complexity constraints. Numerous experiments are presented to illustrate the main features that distinguish adapted and nonadapted filters, as well as the effects of the various constraints. The conjecture that energy compaction and coding gain optimization are equivalent design criteria is shown not to hold for FIR filter banks