Fractional order differentiation by integration with Jacobi polynomials

The differentiation by integration method with Jacobi polynomials was originally introduced by Mboup, Join and Fliess. This paper generalizes this method from the integer order to the fractional order for estimating the fractional order derivatives of noisy signals. The proposed fractional order differentiator is deduced from the Jacobi orthogonal polynomial filter and the Riemann-Liouville fractional order derivative definition. Exact and simple formula for this differentiator is given where an integral formula involving Jacobi polynomials and the noisy signal is used without complex mathematical deduction. Hence, it can be used both for continuous-time and discrete-time models. The comparison between our differentiator and the recently introduced digital fractional order Savitzky-Golay differentiator is given in numerical simulations so as to show its accuracy and robustness with respect to corrupting noises

The fractional orthogonal derivative

This paper builds on the notion of the so-called orthogonal derivative, where an n-th order derivative is approximated by an integral involving an orthogonal polynomial of degree n. This notion was reviewed in great detail in a paper in J. Approx. Theory (2012) by the author and Koornwinder. Here an approximation of the Weyl or Riemann-Liouville fractional derivative is considered by replacing the n-th derivative by its approximation in the formula for the fractional derivative. In the case of, for instance, Jacobi polynomials an explicit formula for the kernel of this approximate fractional derivative can be given. Next we consider the fractional derivative as a filter and compute the transfer function in the continuous case for the Jacobi polynomials and in the discrete case for the Hahn polynomials. The transfer function in the Jacobi case is a confluent hypergeometric function. A different approach is discussed which starts with this explicit transfer function and then obtains the approximate fractional derivative by taking the inverse Fourier transform. The theory is finally illustrated with an application of a fractional differentiating filter. In particular, graphs are presented of the absolute value of the modulus of the transfer function. These make clear that for a good insight in the behavior of a fractional differentiating filter one has to look for the modulus of its transfer function in a log-log plot, rather than for plots in the time domain.Comment: 32 pages, 7 figures. The section between formula (4.15) and (4.20) is correcte

Non-asymptotic fractional order differentiators via an algebraic parametric method

Recently, Mboup, Join and Fliess [27], [28] introduced non-asymptotic integer order differentiators by using an algebraic parametric estimation method [7], [8]. In this paper, in order to obtain non-asymptotic fractional order differentiators we apply this algebraic parametric method to truncated expansions of fractional Taylor series based on the Jumarie's modified Riemann-Liouville derivative [14]. Exact and simple formulae for these differentiators are given where a sliding integration window of a noisy signal involving Jacobi polynomials is used without complex mathematical deduction. The efficiency and the stability with respect to corrupting noises of the proposed fractional order differentiators are shown in numerical simulations

Chebfun and numerical quadrature

Chebfun is a Matlab-based software system that overloads Matlab’s discrete operations for vectors and matrices to analogous continuous operations for functions and operators. We begin by describing Chebfun’s fast capabilities for Clenshaw–Curtis and also Gauss–Legendre, –Jacobi, –Hermite, and –Laguerre quadrature, based on algorithms of Waldvogel and Glaser, Liu, and Rokhlin. Then we consider how such methods can be applied to quadrature problems including 2D integrals over rectangles, fractional derivatives and integrals, functions defined on unbounded intervals, and the fast computation of weights for barycentric interpolation