Fitting polynomials for computing time derivatives

Christophe LETELLIER

- Motivations

As soon as derivatives estimated from experimental data are involved, comes the question whether the contaminating noise is a problem to safely compute derivatives from such data. The answer depends on the way derivatives are estimated. Using Euler scheme or finite difference scheme leads to a negative answer, mainly because the noise contamination is very quickly amplified from one derivative to the other. Using Legendre polynomials as discussed in [1] provides slightly better results but, we observed a significant improvement in the use of derivative coordinates when locally fitted polynomials were analytically derived. Such a technique presents the great advantage to spread the noise contribution in a homogeneous way in all the coordinates. Moreover, successive derivatives share a coherence that other methods do not provide. For instance, in the case of Legendre polynomials, each derivative is estimated independently, a property not recommended for global modelling.

Zip - 3.7 kb
Algorithm for computing time derivatives

- The algorithm

The central idea of our technique to compute derivatives is to use local interpolating polynomials to compute the derivatives instead of computing the derivatives directly from the data. To compute the derivatives throughout the whole time series, a moving window method is employed. The procedure can be summarized as follows:

  • Choose a window w_{L}(1) of lenght L in the beginning of the time series, that is, w_{L}(1)=(x(1),x(2),\ldots,x(L));
  • Find an interpolating polynomial P_{L}(t) of degree L-1 that passes through all the points of W_{L}(1);
  • Compute the derivatives at the center of the window t_{*} by analytically derivating the polynomial at this point;
  • Move the window one sample forward in time and repeat all the steps until all the time series have been processed.
  • There is a low-pass filter in a frequency domain (function smooth). The filtering parameter (vois) is expressed in terms of number of data points. It should be taken as small as possible.

This method was never published before this recent paper about Frequently Asked Questions about global modelling [2].

- An example

The space reconstructed spanned by the successive time derivatives computed from the anode current I(t) recorded in the copper electrodissolution experiments is here shown. It was produced with the Fortran code you can download from this page with the parameters left in the code. In fact the model we got from this experiments was obtained with derivative using Legendre polynomials [3]. This is only later that Loïs Le Sceller developed this procedure to compute derivatives [4]. Then, this was applied to a velocity component of the flow in a mixing tank [5], a Belousov-Zhabotinsky reaction  [6] and the Lynx population [7].

PNG - 7.7 kb
Differential space from the copper electrodissolution

[1] J. F. Gibson, J. D. Farmer, M. Casdagli & S. Eubank, An analytic approach to practical state space reconstruction, Physica D, 57, 1-30, 1992.

[2] C. Letellier, L. A. Aguirre & U. S. Freitas, Frequently Asked Questions about global modelling, submitted.

[3] C. Letellier, L. Le Sceller, E. Maréchal, P. Dutertre, B. Maheu, G. Gouesbet, Z. Fei & J. L. Hudson, Global vector field reconstruction from a chaotic experimental signal in copper electrodissolution, Physical Review E, 51 (5), 4262-4266, 1995.

[4] L. Le Sceller, Reconstruction globale de champ de vecteurs et Applications, Ph’D Thesis defended on February 28, 1997 at CORIA UMR 6614, Université de Rouen.

[5] C. Letellier, L. Le Sceller, G. Gouesbet, F. Lusseyran, A. Kemoun & B. Izrar, Recovering deterministic behavior from experimental time series in a standard mixing reactor, AIChE Journal, 43 (9), 2194-2202, 1997.

[6] C. Letellier, J. Maquet, H. Labro, L. Le Sceller, G. Gouesbet, F. Argoul & A. Arn\’eodo, Analyzing chaotic behaviour in a Belousov-Zhabotinskii reaction by using a global vector field reconstruction, Journal of Physical Chemistry A, 102, 10265-10273, 1998.

[7] J. Maquet, C. Letellier & L. A. Aguirre, Global models from the Canadian Lynx cycles as a first evidence for chaos in real ecosystems, Journal of Mathematical Biology, 55 (1), 21-39, 2007.


Algorithm for computing time derivatives
Zip · 3.7 kb
2044 - 22/03/17

ATOMOSYD © 2007-2017 |  Follow-up of the site's activity  |  SPIP  |  scoty  |  MàJ . 06/02/2017