An algorithm for discrete fractional Hadamard transform

We present a novel algorithm for calculating the discrete fractional Hadamard transform for data vectors whose size N is a power of two. A direct method for calculation of the discrete fractional Hadamard transform requires $N^2$ multiplications, while in proposed algorithm the number of real multiplications is reduced to $N$log$_2N$.Comment: 22 pages, 4 figure

SU(2) and SU(1,1) Approaches to Phase Operators and Temporally Stable Phase States: Applications to Mutually Unbiased Bases and Discrete Fourier Transforms

We propose a group-theoretical approach to the generalized oscillator algebra Ak recently investigated in J. Phys. A: Math. Theor. 43 (2010) 115303. The case k > or 0 corresponds to the noncompact group SU(1,1) (as for the harmonic oscillator and the Poeschl-Teller systems) while the case k < 0 is described by the compact group SU(2) (as for the Morse system). We construct the phase operators and the corresponding temporally stable phase eigenstates for Ak in this group-theoretical context. The SU(2) case is exploited for deriving families of mutually unbiased bases used in quantum information. Along this vein, we examine some characteristics of a quadratic discrete Fourier transform in connection with generalized quadratic Gauss sums and generalized Hadamard matrices

Efficient multistep methods for tempered fractional calculus: Algorithms and Simulations

In this work, we extend the fractional linear multistep methods in [C. Lubich, SIAM J. Math. Anal., 17 (1986), pp.704--719] to the tempered fractional integral and derivative operators in the sense that the tempered fractional derivative operator is interpreted in terms of the Hadamard finite-part integral. We develop two fast methods, Fast Method I and Fast Method II, with linear complexity to calculate the discrete convolution for the approximation of the (tempered) fractional operator. Fast Method I is based on a local approximation for the contour integral that represents the convolution weight. Fast Method II is based on a globally uniform approximation of the trapezoidal rule for the integral on the real line. Both methods are efficient, but numerical experimentation reveals that Fast Method II outperforms Fast Method I in terms of accuracy, efficiency, and coding simplicity. The memory requirement and computational cost of Fast Method II are $O(Q)$ and $O(Qn_T)$, respectively, where $n_T$ is the number of the final time steps and $Q$ is the number of quadrature points used in the trapezoidal rule. The effectiveness of the fast methods is verified through a series of numerical examples for long-time integration, including a numerical study of a fractional reaction-diffusion model

Numerical Approximations to Fractional Problems of the Calculus of Variations and Optimal Control

This chapter presents some numerical methods to solve problems in the fractional calculus of variations and fractional optimal control. Although there are plenty of methods available in the literature, we concentrate mainly on approximating the fractional problem either by discretizing the fractional term or expanding the fractional derivatives as a series involving integer order derivatives. The former method, as a subclass of direct methods in the theory of calculus of variations, uses finite differences, Grunwald-Letnikov definition in this case, to discretize the fractional term. Any quadrature rule for integration, regarding the desired accuracy, is then used to discretize the whole problem including constraints. The final task in this method is to solve a static optimization problem to reach approximated values of the unknown functions on some mesh points. The latter method, however, approximates fractional problems by classical ones in which only derivatives of integer order are present. Precisely, two continuous approximations for fractional derivatives by series involving ordinary derivatives are introduced. Local upper bounds for truncation errors are provided and, through some test functions, the accuracy of the approximations are justified. Then we substitute the fractional term in the original problem with these series and transform the fractional problem to an ordinary one. Hereafter, we use indirect methods of classical theory, e.g. Euler-Lagrange equations, to solve the approximated problem. The methods are mainly developed through some concrete examples which either have obvious solutions or the solution is computed using the fractional Euler-Lagrange equation.Comment: This is a preprint of a paper whose final and definite form appeared in: Chapter V, Fractional Calculus in Analysis, Dynamics and Optimal Control (Editor: Jacky Cresson), Series: Mathematics Research Developments, Nova Science Publishers, New York, 2014. (See http://www.novapublishers.com/catalog/product_info.php?products_id=46851). Consists of 39 page

Some fundamental properties on the sampling free nabla Laplace transform

Discrete fractional order systems have attracted more and more attention in recent years. Nabla Laplace transform is an important tool to deal with the problem of nabla discrete fractional order systems, but there is still much room for its development. In this paper, 14 lemmas are listed to conclude the existing properties and 14 theorems are developed to describe the innovative features. On one hand, these properties make the N-transform more effective and efficient. On the other hand, they enrich the discrete fractional order system theor