automatica 1.0.0

Automatic control systems library
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
56
57
58
59
60
61

use approx::{assert_abs_diff_eq, assert_relative_eq};

use automatica::{signals::continuous, Seconds, Ss, Tf};
use polynomen::{poly, Poly};

/// TC3.1
#[test]
fn stability() {
    let stable_poles = [-1.0_f64, -1.3, -15.];
    let den = Poly::new_from_roots(&stable_poles);
    assert!(Tf::new([0.5, 1.5], den).is_stable());

    let unstable_poles = [-1.0_f32, 0.3, -5.];
    let den2 = Poly::new_from_roots(&unstable_poles);
    assert!(!Tf::new([0.5, 1.5], den2).is_stable());
}

/// TC3.2
#[test]
fn equilibrium() {
    // Es 3.7
    let a = [0.0_f64, 1., -1., -1.];
    let b = [0., 1.];
    let c = [0., 1.];
    let d = [0.];
    let sys = Ss::new_from_slice(2, 1, 1, &a, &b, &c, &d);
    let eq = sys.equilibrium(&[1.]).unwrap();
    assert_eq!(2, eq.x().len());
    assert_relative_eq!(1., eq.x()[0]);
    assert_relative_eq!(0., eq.x()[1]);
    assert_eq!(1, eq.y().len());
    assert_relative_eq!(0., eq.y()[0]);
}

/// TC3.3
#[test]
fn no_equilibrium() {
    // Es 3.9
    let a = [0.0_f32, 0., 1., 2.];
    let b = [1., 2., 3., 4.];
    let c = [5., 6.];
    let d = [0., 0.];
    let sys = Ss::new_from_slice(2, 2, 1, &a, &b, &c, &d);
    let no_eq = sys.equilibrium(&[1., 1.]);
    assert!(no_eq.is_none());
}

/// TC3.4
#[test]
fn static_gain() {
    // 5.4.4
    let tf = Tf::new([1., 5.], poly!(1., 2.) * poly!(1., 1.));
    let sys = Ss::new_observability_realization(&tf).unwrap();

    let step = continuous::step(1., 1);
    let evo = sys.rk2(&step, &[0., 0.], Seconds(0.1), 150);
    let last = evo.last().unwrap();

    let gain = tf.static_gain();
    assert_abs_diff_eq!(gain, last.output()[0], epsilon = 1e-2);
}