899eeb608b
this makes it easier for callers to avoid excessive slope calculations Thanks to Jon Challinger for the idea
54 lines
1.6 KiB
C++
54 lines
1.6 KiB
C++
// -*- tab-width: 4; Mode: C++; c-basic-offset: 4; indent-tabs-mode: nil -*-
|
|
//
|
|
// This is free software; you can redistribute it and/or modify it under
|
|
// the terms of the GNU Lesser General Public License as published by the
|
|
// Free Software Foundation; either version 2.1 of the License, or (at
|
|
// your option) any later version.
|
|
//
|
|
|
|
/// @file Derivative.h
|
|
/// @brief A class to implement a derivative (slope) filter
|
|
/// See http://www.holoborodko.com/pavel/numerical-methods/numerical-derivative/smooth-low-noise-differentiators/
|
|
|
|
#ifndef Derivative_h
|
|
#define Derivative_h
|
|
|
|
#include <inttypes.h>
|
|
#include <Filter.h>
|
|
#include <FilterWithBuffer.h>
|
|
|
|
// 1st parameter <T> is the type of data being filtered.
|
|
// 2nd parameter <FILTER_SIZE> is the number of elements in the filter
|
|
template <class T, uint8_t FILTER_SIZE>
|
|
class DerivativeFilter : public FilterWithBuffer<T,FILTER_SIZE>
|
|
{
|
|
public:
|
|
// constructor
|
|
DerivativeFilter() : FilterWithBuffer<T,FILTER_SIZE>() {};
|
|
|
|
// update - Add a new raw value to the filter, but don't recalculate
|
|
virtual void update(T sample, uint32_t timestamp);
|
|
|
|
// return the derivative value
|
|
virtual float slope(void);
|
|
|
|
// reset - clear the filter
|
|
virtual void reset();
|
|
|
|
private:
|
|
bool _new_data;
|
|
float _last_slope;
|
|
|
|
// microsecond timestamps for samples. This is needed
|
|
// to cope with non-uniform time spacing of the data
|
|
uint32_t _timestamps[FILTER_SIZE];
|
|
};
|
|
|
|
typedef DerivativeFilter<float,5> DerivativeFilterFloat_Size5;
|
|
typedef DerivativeFilter<float,7> DerivativeFilterFloat_Size7;
|
|
typedef DerivativeFilter<float,9> DerivativeFilterFloat_Size9;
|
|
|
|
|
|
#endif // Derivative_h
|
|
|