Numerical differentiation of noisy time series data in python¶
Numerical differentiation methods for noisy time series data in python includes:
from derivative import dxdt import numpy as np t = np.linspace(0,2*np.pi,50) x = np.sin(x) # 1. Finite differences with central differencing using 3 points. result1 = dxdt(x, t, kind="finite_difference", k=1) # 2. Savitzky-Golay using cubic polynomials to fit in a centered window of length 1 result2 = dxdt(x, t, kind="savitzky_golay", left=.5, right=.5, order=3) # 3. Spectral derivative result3 = dxdt(x, t, kind="spectral") # 4. Spline derivative with smoothing set to 0.01 result4 = dxdt(x, t, kind="spline", s=1e-2) # 5. Total variational derivative with regularization set to 0.01 result5 = dxdt(x, t, kind="trend_filtered", order=0, alpha=1e-2)
Symmetric finite difference schemes using arbitrary window size.
Savitzky-Galoy derivatives of any polynomial order with independent left and right window parameters.
Spectral derivatives with optional filter.
Spline derivatives of any order.
Polynomial-trend-filtered derivatives generalizing methods like total variational derivatives.
The goal of this package is to provide some common numerical differentiation techniques that showcase improvements that can be made on finite differences when data is noisy.
This package binds these common differentiation methods to a single easily implemented differentiation interface to encourage user adaptation.
 Numerical differentiation of experimental data: local versus global methods- K. Ahnert and M. Abel
 Numerical Differentiation of Noisy, Nonsmooth Data- Rick Chartrand
 The Solution Path of the Generalized LASSO- R.J. Tibshirani and J. Taylor