You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
54 lines
1.9 KiB
54 lines
1.9 KiB
// -*- tab-width: 4; Mode: C++; c-basic-offset: 4; indent-tabs-mode: nil -*- |
|
|
|
/* |
|
This program is free software: you can redistribute it and/or modify |
|
it under the terms of the GNU General Public License as published by |
|
the Free Software Foundation, either version 3 of the License, or |
|
(at your option) any later version. |
|
|
|
This program is distributed in the hope that it will be useful, |
|
but WITHOUT ANY WARRANTY; without even the implied warranty of |
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
|
GNU General Public License for more details. |
|
|
|
You should have received a copy of the GNU General Public License |
|
along with this program. If not, see <http://www.gnu.org/licenses/>. |
|
*/ |
|
|
|
/// @file LowPassFilter.cpp |
|
/// @brief A class to implement a second order low pass filter |
|
/// Author: Leonard Hall <LeonardTHall@gmail.com> |
|
|
|
#include <inttypes.h> |
|
#include <AP_Math.h> |
|
#include "LowPassFilter2p.h" |
|
|
|
void LowPassFilter2p::set_cutoff_frequency(float sample_freq, float cutoff_freq) |
|
{ |
|
_cutoff_freq = cutoff_freq; |
|
float fr = sample_freq/_cutoff_freq; |
|
float ohm = tanf(PI/fr); |
|
float c = 1.0f+2.0f*cosf(PI/4.0f)*ohm + ohm*ohm; |
|
_b0 = ohm*ohm/c; |
|
_b1 = 2.0f*_b0; |
|
_b2 = _b0; |
|
_a1 = 2.0f*(ohm*ohm-1.0f)/c; |
|
_a2 = (1.0f-2.0f*cosf(PI/4.0f)*ohm+ohm*ohm)/c; |
|
} |
|
|
|
float LowPassFilter2p::apply(float sample) |
|
{ |
|
// do the filtering |
|
float delay_element_0 = sample - _delay_element_1 * _a1 - _delay_element_2 * _a2; |
|
if (isnan(delay_element_0) || isinf(delay_element_0)) { |
|
// don't allow bad values to propogate via the filter |
|
delay_element_0 = sample; |
|
} |
|
float output = delay_element_0 * _b0 + _delay_element_1 * _b1 + _delay_element_2 * _b2; |
|
|
|
_delay_element_2 = _delay_element_1; |
|
_delay_element_1 = delay_element_0; |
|
|
|
// return the value. Should be no need to check limits |
|
return output; |
|
}
|
|
|