ukalman 0.1.0

Super tiny zero allocation filters for embedded
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43

use std::io::Read;
use plotters::prelude::*;
use ukalman::Kalman1D;

fn main() {
    println!("Reading data...");
    let mut f = std::fs::File::open("test_assets/cap_full_throttle.txt").unwrap();
    let mut str = String::new();
    f.read_to_string(&mut str).unwrap();

    println!("Creating plot...");
    let root_drawing_area = BitMapBackend::new("../2.4.png", (1000, 1000))
        .into_drawing_area();

    root_drawing_area.fill(&WHITE).unwrap();

    let mut ctx = ChartBuilder::on(&root_drawing_area)
        .caption("Filtering noisy signal", ("Arial", 30))
        .set_label_area_size(LabelAreaPosition::Left, 40)
        .set_label_area_size(LabelAreaPosition::Bottom, 40)
        .build_cartesian_2d(0..20000isize, 0..3000isize)
        .unwrap();

    ctx.configure_mesh().draw().unwrap();

    println!("starting drawing...");

    ctx.draw_series(LineSeries::new(str.lines().flat_map(|s| s.parse()).scan(0isize, |x, y: isize| {
        *x += 32;
        Some((*x, y as isize))
    }), &BLUE)).unwrap();

    println!("Done drawing normal...");

    let mut kal = Kalman1D::new(0.0, 10.0);

    ctx.draw_series(LineSeries::new(str.lines().flat_map(|s| s.parse()).scan(0isize, |x, y: isize| {
        *x += 32;
        Some((*x, kal.filter(y as f32, 0.005, |x| x, |s| s + 0.0001) as isize))
    }), &RED)).unwrap();

    println!("Done drawing kalman...");
}