Least-squares design of IIR filters with prescribed magnitude andphase responses and a pole radius constraint

This paper presents a method for the frequency domain design of infinite impulse response (IIR) digital filters. The proposed method designs filters approximating prescribed magnitude and phase responses. IIR filters of this kind can have approximately linear-phase responses in their passbands, or they can equalize magnitude and phase responses of given systems. In many cases, these filters can be implemented with less memory and with fewer computations per output sample than equivalent finite impulse response (FIR) digital filters. An important feature of the proposed method is the possibility to specify a maximum radius for the poles of the designed rational transfer function. Consequently, stability can be guaranteed, and undesired effects of implementations using fixed-point arithmetic can be alleviated by restricting the poles to keep a prescribed distance from the unit circle. This is achieved by applying Rouche's theorem in the proposed design algorithm. We motivate the use of IIR filters with an unequal number of poles and zeros outside the origin of the complex plane. In order to satisfy simultaneous specifications on magnitude and phase responses, it is advantageous to use IIR filters with only a few poles outside the origin of the z-plane and an arbitrary number of zeros. Filters of this type are a compromise between IIR filters with optimum magnitude responses and phase-approximating FIR filters. We use design examples to compare filters designed by the proposed method to those obtained by other methods. In addition, we compare the proposed general IIR filters with other popular more specialized structures such as FIR filters and cascaded systems consisting of frequency-selective IIR filters and phase-equalizing allpass filters     

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     

A technique for realizing linear phase IIR filters

A real-time IIR filter structure is presented that possesses exact phase linearity with 10~1000 times fewer general multiplies than conventional FIR filters of similar performance and better magnitude characteristics than equiripple or maximally flat group delay IIR filters. This structure is based on a technique using local time reversal and single pass sectioned convolution methods to realized a real-time recursive implementation of the noncausal transfer function H(z^-1). The time reversed section technique used to realize exactly linear phase IIR filters is described. The effects of finite section length on the sectional convolution are analyzed. A simulation methodology is developed to address the special requirements of simulating a time reversed section filter. A design example is presented, with computer simulation to illustrate performance, in terms of overall magnitude response and phase linearity, as a function of finite section length. Nine example filter specifications are used to compare the performance and complexity of the time reversed section technique to those of a direct FIR implementation     

Interpolated LIR Mth-band filter design with allpass subfilters

We propose a new allpass-based structure for the IIR Mth-and 2Mth-band filters. These filters consist of M allpass filters and an interpolation filter (sum of two allpasses). Consequently, the proposed structure is very efficient in implementation. By choosing the allpass phase appropriately, the resulting phase response of the IIR Mth-band filter becomes approximately linear. An example is designed and compared with FIR Mth-band filters     

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     

基于FIR与相位补偿IIR滤波器的雷达回波信号处理

FIR与IIR频率选择滤波器的设计,被广泛应用于数字信号处理领域之中。文章以雷达回波信号的数字处理为例,首先分别设计FIR,IIR滤波器完成了对信号特定频率分量的滤除。进而,针对IIR滤波器的非线性相位,基于最优化设计全通系统实现了相位补偿,并对FIR,IIR滤波器进行了综合比较。     

A direct approach to H2 optimal deconvolution ofperiodic digital channels

This paper studies the H₂ optimal deconvolution problem for periodic finite impulse response (FIR) and infinite impulse response (IIR) channels. It shows that the H₂ norm of a periodic filter can be directly quantified in terms of periodic system matrices and linear matrix inequalities (LMIs) without resorting to the commonly used lifting technique. The optimal signal reconstruction problem is then formulated as an optimization problem subject to a set of matrix inequality constraints. Under this framework, the optimization of both the FIR and IIR periodic deconvolution filters can be made convex, solved using the interior point method, and computed by using the Matlab LMI Toolbox. The robust deconvolution problem for periodic FIR and IIR channels with polytopic uncertainties are further formulated and solved, also by convex optimization and the LMIs. Compared with the lifting approach to the design of periodic filters, the proposed approach is simpler yet more powerful in dealing with multiobjective deconvolution problems and channel uncertainties, especially for IIR deconvolution filter design. The obtained solutions are applied to the design of an optimal filterbank yielding satisfactory performance     

Broadband Beamforming Using TDL-Form IIR Filters

Conventional broadband beamforming structures make use of finite-impulse-response (FIR) filters in each channel. Large numbers of coefficients are required to retain the desired signal-to-interference-plus-noise-ratio (SINR) performance as the operating bandwidth increases. It has been proven that the optimal frequency-dependent array weighting of broadband beamformers could be better approximated by infinite-impulse-response (IIR) filters. However, some potential problems, such as stability monitoring and sensitivity to quantization errors, of the IIR filters make the implementation of the IIR beamformers difficult. In this paper, new broadband IIR beamformers are proposed to solve these problems. The main contributions of this paper include 1) the Frost-based and generalized sidelobe canceller (GSC)-based broadband beamformers utilizing a kind of tapped-delay-line-form (TDL-form) IIR filters are proposed; 2) the combined recursive Gauss-Newton (RGN) algorithm is designed to compute the feedforward and feedback weights in the Frost-based implementation; and 3) in the GSC-based structure, the unconstrained RGN algorithm is customized for the TDL-form IIR filters in the adaptive beamforming part. Compared with the beamformer using direct-form IIR filters, the new IIR beamformers offer much easier stability monitoring and less sensitivity to the coefficient quantization, while comparable SINR improvement over the conventional FIR beamformer is achieved     

Design of low-order linear-phase IIR filters via orthogonalprojection

This paper presents an indirect linear-phase IIR filter design technique based on a reduction of linear-phase FIR filters. The desired filter is obtained by minimizing the L₂ norm of the difference between the original FIR filter and the lower order IIR filter. We first establish a relationship between the Hankel singular values of the discarded part of the FIR filter and the L₂ norm of the corresponding filter approximation error based on model truncation. This result motivates us to propose a simple finite search method that will achieve better approximation results than commonly used truncation methods such as the balanced truncation (BT) and the impulse response gramian (IRG) methods. We then develop an iterative algorithm for finding an optimal IIR filter based on a matrix projection of the original FIR filter. The convergence of the proposed algorithm is established. Filters designed using the proposed algorithm are compared with those obtained by other techniques with respect to the amplitude response and group delay characteristics in the passband. Numerical examples show that the proposed algorithm offers the best performance     

DESIGN OF POLYPHASE WAVE DIGITAL FILTERS WITH APPROXIMATELY LINEAR PHASE

It is well known that IIR digital filters require quite fewer computations,comparedwith FIR filters,in order to meet stringent magnitude specifications when the phase distortioncan be tolerated.An approximately linear phase,however,can be also obtained with the IIRfilter by making use of a technique without increasing the complexity.Based on a certain numberof attenuation zeros in the pass band,a new approach is developed for the design of polyphasewave digital filters with exact magnitude responses and Chebyshev approximation of the desiredphase responses.The minimum number of attenuation zeros is estimated,and some examples areincluded.     

Frequency domain synthesis of IIR digital filters by means of associated FIR filters

This correspondence describes a synthesis technique for IIR digital filters which allows the use of approximation methods already developed for designing FIR filters. The technique is based on the definition, by means of a suitable transformation, of a FIR filter associated to the desired IIR filter. An example of application is given and the related results are briefly discussed.     

Multidimensional IIR filters and robust rational interpolation

It is well-known that IIR filters can have a much lower order than FIR filters with the same performance. On the downside is that the implementation of an IIR filter is an iterative procedure while that of an FIR filter is a one-shot computation. But in higher dimensions IIR filters are definitely more attractive. We offer a technique where the filter’s performance specifications, stability constraints, its convergence speed and a protection against possible adverse effects of perturbations are all included in the design from the start. The technique only needs an off-the-shelf LP solver because the filter is obtained as a Chebyshev center of a convex polytope. The method deals with general non-causal non-separable filters.     

Revisiting the doppler filter of LEO satellite GNSS receivers for precise velocity estimation

The theoretical aspects of the precise velocity determination of Low Earth Orbit (LEO) satellites’on board Global Navigation Satellite Systems (GNSS) receivers are derived. It shows that the receiver’s Phase Lock Loop (PLL) is required to feature extremely small group delay within its low frequency band, which is in contrast to existing work that proposed wide band linear phase filters. Following this theory, a Finite Impulse Response (FIR) filter is proposed. To corroborate, the proposed FIR filter and an Infinite Impulse Response (IIR) filter lately proposed in literals are implemented in a LEO satellite onboard GNSS receiver. Tests are conducted using a third party commercial GPS signal generator. The results show that the GNSS receiver with the proposed FIR achieves 11 mm/s R.M.S precision, while the GNSS receiver with the IIR filter has a filter-caused velocity error that can not be ignored for space borne GNSS receivers.     

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.     

基于Laguerre滤波器等价设计的IIR宽带波束形成

常规IIR宽带波束形成器可以获得比FIR宽带波束形成器更好的性能,但其需要多极点的自适应调整过程,存在稳定性无法保证,计算复杂度较高等问题.本文提出一种新的基于IIR滤波器的宽带波束形成算法.该算法基于高阶Laguerre宽带波束形成器,利用双线性变换和函数束方法设计相应的低阶等价IIR宽带波束形成器.仿真实验及理论分析表明,该方法无需常规IIR宽带波束形成器的多极点自适应调整过程,在保证算法稳定性的同时,减少了计算复杂度,并提高了输出信干噪比(SINR).     

具有近似线性相位特性的多相波数字滤波器的设计

不考虑相位失真,而要求满足比较苛刻的振幅特性的时候,和FIR滤波器比较起来,采用 IIR 滤波器,计算量要小很多.但对某些IIR滤波器,采取适当的方法,也可以使其相位失真变得很小.本文提出一种新的设计方法,用来设计振幅响应完全满足设计要求,而相位特性用Chebyshev 近似来实现的多相波数字滤波器.这种方法的关键在于确定通带中的若干衰减零点.文中对需要的最小衰减零点数作了估计.实例表明,用这种方法设计的滤波器,时延小,相位特性非常接近线性.     

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.     

Reduced order IIR approximation to FIR digital filters

An algorithm for designing an infinite-impulse-response (IIR) stable filter using a finite-impulse-response (FIR) given filter, with the objective of reducing the delay and order, is described. The design is in the time domain using the least-squares-inverse algorithm, which is briefly described. In this method, the numerator of the approximated filter is part of the FIR filter itself and no calculations and minimization are needed to find the numerator coefficients (except finding the FIR roots). An error analysis between the given FIR and approximated IIR filters is provided. This error analysis enables the designer to fix a design parameter, often unnoted, keeping the energies of the approximated and original filters equal. Results and two illustrative examples are presented     

Optimum FIR and IIR multistage multirate filter design

A theorem is introduced which is useful in deriving equivalent multirate filter structures. Frequency responses of multistage multirate filters are derived and defined by deriving their equivalent one-stage filters. A design principle is proposed to reduce filtering requirements at each stage and move the filter operations to low-sampling-rate stages and thus result in a lower arithmetic rate. Optimum FIR and IIR multistage multirate filter designs are developed based on this principle. The new design has a one-point passband specification for each decimator and/or interpolator stage resulting in a wider transition region and lower filter order. Examples are given to explain the design procedure, and comparisons are made to show the superiority of the new filters.     

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.