Numerical differentiation of experimental data: local versus global methods

DOI10.1016/J.CPC.2007.03.009MaRDI QIDQ710123zbMATHOpenAlexWikidataFDO

Authors Karsten Ahnert, Markus Abel

Publication date 18 October 2010

Published in Computer Physics Communications (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/physics/0510176

zbMATH Keywords

nonparametric regression data analysis numerical differentiation filtering methods smoothing splines

Mathematics Subject Classification ID

Numerical differentiation (65D25)

Abstract: In the context of the analysis of measured data, one is often faced with the task to differentiate data numerically. Typically, this occurs when measured data are concerned or data are evaluated numerically during the evolution of partial or ordinary differential equations. Usually, one does not take care for accuracy of the resulting estimates of derivatives because modern computers are assumed to be accurate to many digits. But measurements yield intrinsic errors, which are often much less accurate than the limit of the machine used, and there exists the effect of ``loss of significance, well known in numerical mathematics and computational physics. The problem occurs primarily in numerical subtraction, and clearly, the estimation of derivatives involves the approximation of differences. In this article, we discuss several techniques for the estimation of derivatives. As a novel aspect, we divide into local and global methods, and explain the respective shortcomings. We have developed a general scheme for global methods, and illustrate our ideas by spline smoothing and spectral smoothing. The results from these less known techniques are confronted with the ones from local methods. As typical for the latter, we chose Savitzky-Golay filtering and finite differences. Two basic quantities are used for characterization of results: The variance of the difference of the true derivative and its estimate, and as important new characteristic, the smoothness of the estimate. We apply the different techniques to numerically produced data and demonstrate the application to data from an aeroacoustic experiment. As a result, we find that global methods are generally preferable if a smooth process is considered. For rough estimates local methods work acceptably well.

Recommendations

Cites work

Cited in

(13)

This page was built for publication: Numerical differentiation of experimental data: local versus global methods

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q710123)