point_process 0.2.3

A crate for simulating random point processes.
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
44
45
46
47
48
49
50
51
52
53
54
55

extern crate pointprocesses;
extern crate gnuplot;
extern crate serde_json;

use gnuplot::{Figure,Caption,Color,PointSymbol,PointSize};
use gnuplot::AxesCommon;

use pointprocesses::variable_poisson;

fn main() {
    
    let tmax = 60.0;
    let f: fn(f64) -> f64 = |t| {
        1.0 + (0.5*t).sin()*(-0.05*t).exp()
    };
    let events = variable_poisson(tmax, f, 2.0);

    println!("{}", serde_json::to_string_pretty(&events).unwrap());
    
    // Plotting
    let num_points = 100;
    let times: Vec<f64> = (0..num_points).map(|x| {
        tmax*x as f64/num_points as f64
    }).collect();

    let lambda_values: Vec<f64> = times.iter()
        .map(|&x| f(x)).collect();

    let mut event_times: Vec<f64> = vec!();
    let mut event_intens: Vec<f64> = vec!();
    for i in 0..events.len() {
        let event = &events[i];
        event_times.push(event.timestamp());
        event_intens.push(event.intensity());
    }

    let mut fg = Figure::new();

    fg.axes2d()
        .lines(&times, &lambda_values,
            &[
                Caption("λ(t)"),
                Color("#1E90FF"),
            ])
        .points(&event_times, &event_intens, 
            &[
                Caption("Events"),
                Color("black"),
                PointSymbol('O'),
                PointSize(0.8)])
        .set_x_label("Temps t", &[])
        .set_y_label("Intensité", &[]);
    fg.echo_to_file("test.gnuplot");
    fg.show();
}