Table of Contents
Fetching ...

On reconstructing high derivatives of noisy time-series with confidence intervals

Mazen Alamir

TL;DR

Reconstructing high-order derivatives from noisy time-series is challenging; the authors introduce a bandwidth- and noise-aware dictionary of linear operators $A_{(j, ell,d)}$ learned by cross-validation, and apply moving-window averaging to produce $oldsymbol{ ilde s}^{[d]}$ with a corresponding confidence interval $oldsymbol{ ilde oldsymbol ext{sigma}}^{[d]}$. The method demonstrates higher accuracy than Kalman-filtering, spectral derivation, Savitzky–Golay, and implicit sliding-mode approaches across derivative orders up to $d_{max}=4$, and remains effective for extensions to higher derivatives. Core components include a bandwidth-limited basis $B$, projection operators $oldsymbol{ ilde s}(oldsymbol{ar ext{Ω}})$ and $oldsymbol{ ilde ext{P}}$, and a design pulsation grid $ar{ ext{Ω}}_{design}$ used to populate the dictionary. The framework is particularly valuable for offline continuous-time identification and normality characterization, while maintaining potential for online use on shorter windows due to favorable computation times.

Abstract

Reconstructing high derivatives of noisy measurements is an important step in many control, identification and diagnosis problems. In this paper, a heuristic is proposed to address this challenging issue. The framework is based on a dictionary of identified models indexed by the bandwidth, the noise level and the required degrees of derivation. Each model in the dictionary is identified via cross-validation using tailored learning data. It is also shown that the proposed approach provides heuristically defined confidence intervals on the resulting estimation. The performance of the framework is compared to the state-of-the-art available algorithms showing noticeably higher accuracy. Although the results are shown for up to the 4-th derivative, higher derivation orders can be used with comparable results.

On reconstructing high derivatives of noisy time-series with confidence intervals

TL;DR

Reconstructing high-order derivatives from noisy time-series is challenging; the authors introduce a bandwidth- and noise-aware dictionary of linear operators learned by cross-validation, and apply moving-window averaging to produce with a corresponding confidence interval . The method demonstrates higher accuracy than Kalman-filtering, spectral derivation, Savitzky–Golay, and implicit sliding-mode approaches across derivative orders up to , and remains effective for extensions to higher derivatives. Core components include a bandwidth-limited basis , projection operators and , and a design pulsation grid used to populate the dictionary. The framework is particularly valuable for offline continuous-time identification and normality characterization, while maintaining potential for online use on shorter windows due to favorable computation times.

Abstract

Reconstructing high derivatives of noisy measurements is an important step in many control, identification and diagnosis problems. In this paper, a heuristic is proposed to address this challenging issue. The framework is based on a dictionary of identified models indexed by the bandwidth, the noise level and the required degrees of derivation. Each model in the dictionary is identified via cross-validation using tailored learning data. It is also shown that the proposed approach provides heuristically defined confidence intervals on the resulting estimation. The performance of the framework is compared to the state-of-the-art available algorithms showing noticeably higher accuracy. Although the results are shown for up to the 4-th derivative, higher derivation orders can be used with comparable results.

Paper Structure

This paper contains 11 sections, 35 equations, 3 figures, 1 table.

Figures (3)

  • Figure 1: Typical results of derivative reconstruction up to order 4 under noisy measurements with standard deviation of 5%. Raw signal of bandwidth 1.94 Hz. Yellow bands show the 3$\sigma$-confidence intervals. Notice that both exact and estimated values are plotted.
  • Figure 2: Percentiles of reconstruction errors on the first four derivatives for the different algorithms (logarithmic scale). These statistics are computed over all the scenarios including the ones where a standard deviation of 10% is used for the noise. The total number of samples is equal to 192000. The $y$-scale is lower bounded by $10^{-2}$ for a better readability of the comparison. The best parameters is used based on the ground truth for each reconstruction of the alternative solutions.
  • Figure 3: Evolution of the statistics of the error as the measurement noise increases for the different algorithms. The best parameters is used based on the ground truth for each reconstruction of the alternative solutions.