Variable digital filter design using the outer product expansion

Digital filters with adjustable frequency domain characteristics are referred to as variable digital filters. Variable filters are useful in the applications where the filter characteristics are required to be changeable during the course of signal processing. Especially in real time applications, variable filters are needed to change their coefficients instantaneously such that the real time signal processing can be performed. The present paper proposes a very efficient technique for variable 1D digital filter design. Generally speaking, the variable coefficients of variable digital filters are multidimensional functions of a set of spectral parameters which define the desired frequency domain characteristics. The authors first sample the given variable 1D magnitude specification and use the samples to construct a multidimensional array, then propose an outer product expansion method for expanding the multidimensional array as the sum of outer products of 1D arrays (vectors). Based on the outer product expansion, one can reduce the difficult problem of designing a variable 1D digital filter to the easy one that only needs constant 1D filter designs and 1D polynomial approximations. The technique can obtain variable 1D filters having arbitrary desired magnitude characteristics with a high design accuracy     

Design of arbitrary-phase variable digital filters using SVD-based vector-array decomposition

This paper proposes a straightforward method for designing variable digital filters with arbitrary variable magnitude as well as arbitrary fixed-phase or variable fractional delay (VFD) responses. The basic idea is to avoid the complicated direct design of one-dimensional (1-D) variable digital filters by decomposing the original variable filter design problem into easier subproblems that only require constant 1-D filter designs and multidimensional polynomial approximations. Through constant 1-D filter designs and multidimensional polynomial fits, we can easily obtain a variable digital filter satisfying the given variable design specifications. To decompose the original variable filter design into constant 1-D filter designs and multidimensional polynomial fits, a new multidimensional complex array decomposition called vector array decomposition (VAD) is proposed, which is based on two new theorems using the singular value decomposition (SVD). Once the VAD is obtained, the subproblems can be easily solved. Furthermore, we show that the VAD can also be generalized to the weighted least squares (WLS) case (WLS-VAD) for the WLS variable filter design. Three design examples are given to illustrate that the WLS-VAD and VAD-based techniques are considerably efficient for designing variable digital filters with arbitrary variable magnitude and arbitrary fixed-phase or VFD responses.     

二维零相位FIR数字滤波器设计的闭式最小二乘解

本文二维零相位ＦＩＲ数字滤波器的解析最小二乘设计技术。通过建立频域误差差函数的矩阵形式，并运用与设计问题有关的矩阵的一些性质，得到了滤波器系数的闭式解，使得由给定的频响指标可直接计算滤波器系数，而不必对矩阵进行数值示逆，也不需要基于迭代运算的优化过程。文中给出了滤波器实例，其结果证实了该设计方法的简便性与有效性。     

A novel nonnegative decomposition method and its application to 2-D digital filter design

In designing two-dimensional (2-D) digital filters in the frequency domain, an efficient technique is to first decompose the given 2-D frequency domain design specifications into one-dimensional (1-D) ones, and then approximate the resulting 1-D magnitude specifications using the well-developed 1-D filter design techniques. Finally, by interconnecting the designed 1-D filters one can obtain a 2-D digital filter. However, since the magnitude responses of digital filters must be nonnegative, it is required that the decomposition of 2-D magnitude specifications result in nonnegative 1-D magnitude specifications. We call such a decomposition the nonnegative decomposition. This paper proposes a nonnegative decomposition method for decomposing the given 2-D magnitude specifications into 1-D ones, and then transforms the problem of designing a 2-D digital filter into that of designing 1-D filters. Consequently, the original problem of designing a 2-D filter is significantly simplified.     

具有任意幅度频响的二维线性相位FIR数字滤波器的设计

本文提出了设计任意幅度频响的二维线性相位ＦＩＲ数字滤波器的解析最小二乘方法，通过最小化频域平方误差函数得到了滤波器系数的闭式解，运用导出的闭式式，可根据给定的任意幅度频响指标直接计算滤波器的系数，从而简化了滤波器设计程序，并大大降低了运算量。     

Closed-form design and efficient implementation of variable digital filters with simultaneously tunable magnitude and fractional delay

This paper proposes a closed-form solution for designing variable one-dimensional (1-D) finite-impulse-response (FIR) digital filters with simultaneously tunable magnitude and tunable fractional phase-delay responses. First, each coefficient of a variable FIR filter is expressed as a two-dimensional (2-D) polynomial of a pair of parameters called spectral parameters; one is for independently tuning the cutoff frequency of the magnitude response, and the other is for independently tuning fractional phase-delay. Then, the closed-form error function between the desired and actual variable frequency responses is derived without discretizing any design parameters such as the frequency and the two spectral parameters. Finally, the optimal solution for the 2-D polynomial coefficients can be easily determined through minimizing the closed-form error function. We also show that the resulting variable FIR filter can be efficiently implemented by generalizing Farrow structure to our two-parameter case. The generalized Farrow structure requires only a small number of multiplications and additions for obtaining any new frequency characteristic, which is particularly suitable for high-speed tuning.     

基于神经网络的FRM滤波器优化设计研究

以频率响应屏蔽(FRM)技术为基础,提出了一种基于神经网络的窄过渡带FIR数字滤波器的优化设计新方法.该算法主要通过使频率响应平方误差函数最小化来获得FRM滤波器系数.文中详细介绍了基于神经网络的基本FRM滤波器和多层FRM滤波器的设计算法及设计步骤,证明了该算法的稳定性定理,给出了仿真实例,并与已有的设计方法进行了比较,设计结果表明用该方法设计的窄过渡带FIR数字滤波器性能更为优越.     

Two-Dimensional Farrow Structure and the Design of Variable Fractional-Delay 2-D FIR Digital Filters

In this paper, a 2-D Farrow structure is proposed and used to implement variable fractional-delay (VFD) 2-D FIR digital filters. Compared with the existing literature, the desired response of a VFD 2-D digital filter is analyzed in detail, and it is found that there are four types of 2-D symmetric/antisymmetric sequences that need to be used for the design of VFD 2-D FIR digital filters. Moreover, due to the orthogonality among the approach functions, the four types of 2-D sequences can be determined independently, such that the dimension for each computation can be reduced drastically. For simplicity, only the designs of even–even- and odd–odd-order VFD 2-D filters are presented in this paper, and the other cases can be achieved in the same manner. To reveal the coefficient characteristics, the symmetric/antisymmetric properties of filter coefficients and the relationships between coefficients are all tabulated. Also, design examples such as nonseparable circularly symmetric low-pass VFD filters are presented to demonstrate the effectiveness of the proposed method.      

Efficient design of digital filters for 2-D and 3-D depth migration

Two- and three-dimensional (2-D and 3-D) depth migration can be performed using 1-D and 2-D extrapolation digital filters, respectively. The depth extrapolation is done, one frequency at a time, by convolving the seismic wavefield with a complex-valued, frequency- and velocity-dependent, digital filter. This process requires the design of a complete set of extrapolation filters: one filter for each possible frequency-velocity pair. Instead of independently designing the frequency- and velocity-dependent filters, an efficient procedure is introduced for designing a complete set of 1-D and 2-D extrapolation filters using transformations. The problem of designing a desired set of migration filters is thus reduced to the design of a single 1-D filter, which is then mapped to produce all the desired 1-D or 2-D migration filters. The new design procedure has the additional advantage that both the 1-D and 2-D migration filters can be realized efficiently and need not have their coefficients precomputed or tabulated     

QPSO算法在FIR数字滤波器设计中的应用

FIR数字滤波器的设计实际上是一个多维变量寻优问题,滤波器的设计可转化为滤波器参数优化的问题。文章介绍了粒子群优化算法、量子粒子群优化算法,然后利用算法对参数空间进行搜索以获得参数的最优化,根据预期频率特性的设计要求,通过仿真实验表明设计方法的有效性。     

Design of Variable Two-Dimensional FIR Digital Filters by McClellan Transformation

In this paper, the technique of McClellan transformation is applied to design variable 2-D FIR digital filters. Compared with the conventional transformation, the 2-D transformation subfilter and the 1-D prototype filter are designed such that their frequency characteristics are adjustable. Moreover, they are tunable by the same variable parameter, so the variable characteristics of 1-D prototype filters are coincident with those of 2-D subfilters. Several examples, including variable fan filters, variable circularly symmetric filters, and variable elliptically symmetric filters with arbitrary orientation, are presented to demonstrate the effectiveness and the flexibility of the presented method.      

A Very Efficient Approach for the Synthesis of 2-D Recursive Fan Filters

ＡＶｅｒｙＥｆｆｉｃｉｅｎｔＡｐｐｒｏａｃｈｆｏｒｔｈｅＳｙｎｔｈｅｓｉｓｏｆ２－ＤＲｅｃｕｒｓｉｖｅＦａｎＦｉｌｔｅｒｓＺｈｕＷｅｉｐｉｎｇ（ＤｅｐａｒｔｍｅｎｔｏｆＴｅｌｅｃｏｍｍｕｎｉｃａｔｉｏｎＥｎｇｉｎｅｅｒｉｎｇ，ＮａｎｊｉｎｇＵｎｉｖｅ...     

Two-dimensional digital filters with sparse coefficients

Is sparsity an issue in 2-D digital filter design problems to explore and why is it important? How a 2-D filter can be designed to retain a desired coefficient sparsity for efficient implementation while achieving best possible performance subject to that sparsity constraint? These are the focus of this paper in which we present a two-phase design method for 2-D FIR digital filters in two most common design settings, namely, the least squares and minimax designs. Simulation studies are presented to illustrate each phase of the proposed design method and to evaluate the performance of the filters designed.     

SVD-based design and new structures for variable fractional-delay digital filters

This paper shows that the problem of designing one-dimensional (1-D) variable fractional-delay (VFD) digital filter can be elegantly reduced to the easier subproblems that involve one-dimensional (1-D) constant filter (subfilter) designs and 1-D polynomial approximations. By utilizing the singular value decomposition (SVD) of the variable design specification, we prove that both 1-D constant filters and 1-D polynomials possess either symmetry or anti-symmetry simultaneously. Therefore, a VFD filter can be efficiently obtained by designing 1-D constant filters with symmetrical or antisymmetrical coefficients and performing 1-D symmetrical or antisymmetrical approximations. To perform the weighted-least-squares (WLS) VFD filter design, a new WLS-SVD method is also developed. Moreover, an objective criterion is proposed for selecting appropriate subfilter orders and polynomial degrees. Our computer simulations have shown that the SVD-based design and WLS-SVD design can achieve much higher design accuracy with significantly reduced filter, complexity than the existing WLS design method. Another important part of the paper proposes two new structures for efficiently implementing the resulting VFD filter, which require less computational complexity than the so-called Farrow structure.     

Three-dimensional rotated digital filters: Design,stability, and applications

Three-dimensional (3-D) signal processing offers many advantages over two-dimensional (2-D) processing, because it preserves 3-D correlations. In this paper the design and the stability of 3-D rotated filters are considered. These filters are designed by rotating a one-dimensional (1-D) digital filter in 3-D space. The rotated filters are valuable in the design of various 3-D filters which possess prescribed spectral specifications. An efficient algorithm for the design of 3-D lowpass (LP) digital filters, with approximately spherically symmetric magnitude responses, is introduced. To achieve the desirable spectral characteristics, a number of 3-D rotated filters is cascaded. The stability of the spherically symmetric filters designed is considered, and stable realizations are proposed. The relation between the cut-off isopotential sphere of the 3-D filter and the cut-off frequency of the 1-D filter employed in the design, is derived. Finally, configurations that result in highpass (HP) and bandpass (BP) filters are proposed. Examples of LP, HP, and BP filters, designed on the basis of the method proposed, are presented.This research was supported by the Public Benefit Foundation, Alexander S. Onassis.     

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

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

一维线性相位FIR数字滤波器的频域最小平方误差设计方法

本文提出了一种设计一维线性相位FIR数字滤波器的新方法。该方法采用频域的最小平方误差函数使所求的增益特性逼近所希望的增益特性 ,其计算公式非常简单 ,并可设计任意形状增益特性的FIR数字滤波器。设计实例表明提出的方法非常有效     

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.     

Realization of 2-D Variable IIR Digital Filter Structures with a Small Amount of Calculations for Coefficient Update

This paper proposes 2-D variable IIR digital filter structures with a small amount of calculations for coefficient update. The proposed realization method uses the 2-D parallel allpass structure derived from the separable denominator 2-D filter as the prototype structure for 2-D variable digital filters. In order to reduce the amount of calculations, all the redundant first-order complex allpass sections are combined by modularization of the variable structure. Furthermore, we can realize a very compact variable structure with a minimal number of first-order complex allpass sections by combining complex allpass sections with their complex conjugate allpass sections. Comparison of the calculation loads of the variable structures is presented to demonstrate that the amount of calculations for coefficient update of the proposed variable structure is far less than that of the original and the modular variable structure.     

基于FPGA的嵌入式程控数字滤波器的实现研究

提出了一种程控数字滤波器的FPGA实现方法.以FIR数字滤波器为例,通过在单片FPGA芯片中构建嵌入式微处理器,加入系数可配置FIR滤波器等功能模块实现了频谱可控的程控数字滤波器.通过Modelsim后仿真和在Xilinx公司XUPV5-LX110T开发板上进行的实际测试,结果表明本文提出的方法是可行的.