Error analysis of numerical Weyl fractional derivatives in the case of certain Hölder continuous functions
Nissilä, Juhani (2018-08-02)
J. Nissilä, "Error analysis of numerical Weyl fractional derivatives in the case of certain Hölder continuous functions," 2018 14th IEEE/ASME International Conference on Mechatronic and Embedded Systems and Applications (MESA), Oulu, Finland, 2018, pp. 1-7. doi: 10.1109/MESA.2018.8449193
© 2018 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes,creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.
The calculation of fractional or integer order derivatives and integrals has been demonstrated to be simple and fast in the frequency domain. It is also the most sensible method if one wishes to calculate derivatives or integrals of periodic signals. In this paper, error analysis is carried out for the numerical algorithm for Weyl fractional derivatives. To derive an upper bound for the numerical error, some knowledge of the smoothness of the signal must be known in advance or it must be estimated. The derived error analysis is tested with sampled functions with known regularity and with real vibration measurements from rotating machines. Compared to previous publications which deal with error analysis of integer order numerical derivatives in the frequency domain using L 2 errors, the result of this paper is in terms of maximum absolute error and it is based on a novel result on the signal’s regularity. The general conclusion using either error estimates is the same: the error of numerical Weyl derivatives is bounded by some constant times the sequence length raised to a negative power. The exponent depends on the smoothness of the signal. This contrasts with using difference quotients in numerical differentiation, in which case the error is bounded by a constant times the sequence length raised to a some fixed negative power and the order of the method defines that exponent.
- Avoin saatavuus