Document Type

Article

Publication Date

12-8-2015

Publication Title

Academia

Abstract

In this work we develop a means to rapidly and accurately compute the Caputo fractional derivative of a function, using fast convolution. The key element to this approach is the compression of the fractional kernel into a sum of M decaying exponentials, where M is minimal. Specifically, after N time steps we find M= O (log N) leading to a scheme with O (N log N) complexity. We illustrate our method by solving the fractional differential equation representing the Kelvin-Voigt model of viscoelasticity, and the partial differential equations that model the propagation of electromagnetic pulses in the Cole-Cole model of induced dielectric polarization.

First Page

1

Last Page

21

ISSN

2077-9038

Rights

Academia

Included in

Mathematics Commons

Share

COinS