Design of digital FIR variable fractional order integrator and differentiator

This paper presents an easy and simple method to design variable fractional order digital FIR integrators and differentiators based on fractional order systems. First, closed-form digital IIR fractional order integrators and differentiators have been obtained from the analog rational functions approximations, in a given frequency band, of the fractional order integrator s ^?m and differentiator s ^m (0?<?m?<?1) through the Tustin generating function. Then, closed-form digital FIR fractional order integrators and differentiators by truncation of the digital IIR ones have been derived. Next, polynomial interpolation has been used to design digital FIR variable fractional order integrators and differentiators that can be implemented by the Farrow structure. The main feature of variable fractional order operator is that its order can be changed without re-designing a new fractional order operator. Some examples have been presented through the paper to illustrate the performance and the effectiveness of the proposed design method. The results obtained have been discussed and have been compared to one of the most recent works in the literature using the same design parameters.     

分数阶数字FIR微分器的快速WSLD设计算法

一阶逼近格林瓦尔-莱特尼科夫(G-L)加权系数的计算具有准确快速的递推公式,而高阶逼近鲁比希加权系数的求解则复杂度高,计算消耗时间长。本文通过傅里叶变换证明了鲁比希算子的逼近阶,并基于移位鲁比希算子提出一类四阶逼近的加权移位鲁比希差分(WSLD)算子。从数字信号处理角度分析WSLD算子滤波特性,设计基于WSLD算子的分数阶数字FIR微分滤波器并进行数值仿真验证。对比Al-Alaoui、鲁比希2种典型分数阶算子,结果表明,利用WSLD算子求解分数阶数字FIR滤波器滤波系数的算法简单、高效,且相对其他算子能有效减小吉布斯效应影响。     

On the practical realization of higher-order filters with fractional stepping

We propose the use of a compact integer-order transfer function approximation of the fractional-order Laplacian operator s^α to realize fractional-step filters. Lowpass and bandpass filters of orders (n+α) and 2(n+α), where n is an integer and 0<α<1, can, respectively, be designed. A 5th-order lowpass filter with fractional steps from 0.1 to 0.9 (i.e. 5.1→5.9) is given as an example with its characteristics compared to 5th- and 6th-order Butterworth filters. Spice simulations and experimental results are shown.     

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.     

Thiele’s continued fractions in digital implementation of noninteger differintegrators

A rational approximation is the preliminary step of all the indirect methods for implementing digital fractional differintegrators s ^??, with ${\nu \in \mathbb{R}, 0<|\nu| <1 }$ , and where ${s \in \mathbb{C}}$ . This paper employs the convergents of two Thiele??s continued fractions as rational approximations of s ^??. In a second step, it uses known s-to-z transformation rules to obtain a rational, stable, and minimum-phase z-transfer function, with zeros interlacing poles. The paper concludes with a comparative analysis of the quality of the proposed approximations in dependence of the used s-to-z transformations and of the sampling period.     

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.     

The power spectrum of digital FM as produced by digital circuits

A simple exact method is given for calculating the power spectrum of a frequency modulated filtered M-ary data signal. The power spectrum and its discrete components are given in closed form for any pulse shape. The pulses are those generated by a digital transversal filter and the FM modulator is that composed of an ideal digital integrator followed by a phase modulator (look-up table).     

Chebyshev approximation problem of multidimensional IIR digital filters in the frequency domain

This paper investigates the problem of simultaneous approximation of a prescribed multidimensional frequency response. The frequency responses of multidimensional IIR digital filters are used as nonlinear approximating functions. Chebyshev approximation theory and the notion of line homotopy are used to reveal the approximation properties of this set of IIR functions. A sign condition is derived to characterize a convex stable domain in this set. This sign condition can be incorporated into the optimization of the Chebyshev simultaneous approximation. The generally sufficient global Kolmogorov criterion is shown to be a necessary condition, for the characterization of best approximation, in the considered set of approximating functions. Thus, it can be incorporated, as a stopping constraint, in the design of the optimal filter. Moreover, the local Kolmogorov criterion is shown to be also necessary for the set of approximating functions. Finally, it is proved that the best approximation is a global minimum.

     

Comparison between two approximation methods of state space fractional systems

A large number of real dynamic systems are better characterized using a non-integer order dynamic model based on fractional order differential equations. Nevertheless, their analytical solutions are rather complicated. In this paper we present two alternative methods of simulation of the fractional order dynamical system in state space. We particularly present a comparison between two approximation methods. The first one is based on the differentiation operator approximation and thus uses the usual form of the state space representation D^(α)(x)=Ax+Bu. The second method is based on the approximation of the integration operator and uses, in this case, the new representation w=AI^(α)(w)+Bu. The comparison between these two methods is made on the basis of initial and steady state approximation state and output errors which are characterized in this present paper. Numerical examples are also given to support this comparison.     

FPGA-based Novel Adaptive Scheme Using PN Sequences for Self-Calibration and Self-Testing of MEMS-based Inertial Sensors

We propose a novel adaptive technique based on pseudo-random (PN) sequences for self-calibration and self-testing of MEMS-based inertial sensors (accelerometers and gyroscopes). The method relies on using a parameterized behavioral model implemented on FPGA, whose parameters values are adaptively tuned, based on the response to test pseudo-random actuation of the physical structure. Dedicated comb drives actuate the movable mass with binary maximum length pseudo-random sequences of small amplitude, to keep the device within the linear operating regime. The frequency of the stimulus is chosen within the mechanical spectral operating range of the micro-device, such that the induced response leads to the identification of the mechanical transfer function, and to the tuning of the associated digital behavioral model. In case of a micro-gyroscope, experimental results demonstrate the adaptive tracking of the damping coefficient from 5.57?×?10^?5? Kg/s to 7.12?×?10^?5? Kg/s and of the stiffness coefficient from 132?N/m to 137.7?N/m. In the case of a MEMS accelerometer, the damping and stiffness coefficients are correctly tracked from 3.4?×?10^?3? Kg/s and 49.56?N/m to 4.57?×?10^?3? Kg/s and 51.48?N/m, respectively—the former values are designer-specified target values, while the latter are experimentally measured parameters for fabricated devices operating in a real environment. Hardware resources estimation confirms the small area the proposed algorithm occupies on the targeted FPGA device.     

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     

基于FPGA的IIR数字滤波器的实现

数字信号处理在科学和工程技术许多领域中得到广泛的应用,与FIR数字滤波器相比,IIR数字滤波器可以用较低的阶数获得较高的选择性,本文采用一种基于FPGA的IIR数字滤波器的设计方案,首先分析了IIR数字滤波器的原理及设计方法,然后通过MAX＋PLUSⅡ的设计平台,采用自顶向下的模块化设计思想将整个IIR数字滤波器分为：时序控制、延时、补码乘加和累加4个功能模块。分别对各模块采用VHDL进行描述后,进行了仿真和综合。仿真结果表明,本课题所设计的IIR数字滤波器运算速度较快,系数改变灵活,有较好的参考价值。     

基于免疫算法的IIR数字滤波器优化与实现

由于IIR数字滤波器设计实质上是一个非线性高维复杂函数优化问题,文中提出基于具有全局搜索能力强,收敛速度快特点的免疫算法实现IIR数字滤波器优化设计的新方法,给出了IIR滤波器优化设计的数学模型,描述了应用免疫算法优化设计IIR数字滤波器的具体实现步骤。通过低通和高通IIR数字滤波器设计的仿真结果表明方法的有效性和高效性。     

A fast algorithm for FIR digital filtering with a sum-of-triangles weighting function

A doubly recursive algorithm for time domain convolution with a piecewise linear weighting function is presented that combines the speed of a recursive (IIR) digital filter with the flexibility and ease of design of a nonrecursive (FIR) digital filter. The approach approximates the desired FIR weighting function by a sum-of-triangles weighting function. ForL triangles (or triangle pairs for a linear phase filter) the algorithm is of orderLN. The approximation improves with the number of triangles. A significant advantage of the algorithm compared to FFT filtering or direct convolution is that there is no necessity of a tradeoff between frequency response accuracy and computation time per output point as the data spacing decreases in the filtered signal. The computational complexity is dependent on the number of triangles chosen, not the width of the weighting function, so the algorithm is especially effective for filters with an inherently wide FIR weighting function.     

The transmission and reflection coefficients of an electromagnetic wave for a gradient dielectric layer

A gradient dielectric layer is considered. In the layer, permittivity ?(z) depends on coordinate z as follows: ?(z) = ?_L a ⁴/(z ? a)⁴ (a being the gradientness parameter). For this layer, reflection and transmission coefficients r _r and r _t are determined. It is found that, in a wide frequency range, the considered gradient layer has amplitude- and phase-constant reflection and transmission coefficients and can serve as a low-pass filter in radio devices involved in the transmission of wideband signals.     

Hubbard energy of two-electron tin centers in PbS1−zTez solid solutions

It is shown by Mössbauer spectroscopy of the isotope ¹¹⁹Sn that an isovalent tin impurity in PbS_1?z Te_z solid solutions is a two-electron donor with a negative correlation energy, where the energy levels associated with tin centers are situated against the background of the valence band continuum. The Hubbard energy U is estimated for impurity tin atoms in PbS_1?z Te_z (|U|>0.2 eV), which is found to be substantially higher than for analogous tin centers in PbS_1?z Se_z solid solutions (|U|=0.058 eV).     

PCLS IIR digital filters with simultaneous frequency responsemagnitude and group delay specifications

This paper presents the peak-constrained least-squares (PCLS) approach to designing IIR digital filters. PCLS IIR digital filters that meet simultaneous specifications on the frequency response magnitude and the group delay are introduced. As a point of reference, we consider the IIR digital filter design problem that appears in Deczky's (1972) classic paper and in the popular textbook by Oppenheim and Schafer (1989). In addition, the same design problem appears in the IIR filter design chapter by Higgins and Munson (1993) in the Handbook for Digital Signal Processing. By using our new algorithm with simultaneous optimization of the frequency response magnitude and the group delay, we obtain a dramatic improvement in the solution of this classic IIR digital filter design problem. Starting from the same filter structure and the same specifications for the frequency response magnitude as in the works of Deczky, Oppenheim and Schafer, and Higgins and Munson, we are able to reduce the group delay ripple by a factor of 35. In another design problem that originated in Deczky's work, we use PCLS optimization to reduce the group delay ripple by a factor of 40 at the same time we reduce the stopband energy by 6 dB, without sacrificing any other performance measure. The group delay ripple in this IIR digital filter example is reduced to only ±0.002 samples     

Suppression of transients in variable recursive digital filterswith a novel and efficient cancellation method

A new method for suppressing transients in recursive infinite impulse response (IIR) digital filters is proposed. The technique is based on modifying the state (delay) variables of the filter when coefficients are changed so that the filter enters a new state smoothly without transient attacks, as originally proposed by Zetterberg and Zhang (1988). In this correspondence, we modify the Zetterberg-Zhang algorithm to render it feasible for efficient implementation. We define a mean square error (MSE) measure for transients and determine the optimal transient suppressor to cancel the transients down to a desired level at the minimum complexity of implementation. The application of the method to all-pole and direct-form II (DF II) IIR filter sections is studied in detail. Time-varying recursive filtering with transient elimination is illustrated for tunable fractional delay filters and variable-bandwidth lowpass filters     

Analytical impulse response of a fractional second order filter and its impulse response invariant discretization

In this paper, we derive the impulse response of a fractional second order filter of the form (s²+as+b)^−γ, where a,b?0 and γ>0. The asymptotic properties of the impulse responses are obtained. Moreover, based on the derived analytical impulse response, we show how to perform the discretization of the above fractional second order filter. Finally, a number of illustrated examples in time and frequency domains are provided as proofs of concepts.