reda-spice 0.1.0

Spice simulate and parse 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
62
63
64
65
66
67
68

use core::f64;
use derive_builder::Builder;
use reda_unit::{num, Frequency, Number, Time, Voltage, u};

#[derive(Debug, Clone, Builder)]
#[builder(setter(into))]
pub struct SineVoltage {
    pub vo: Voltage,
    pub va: Voltage,
    pub freq_hz: Frequency,
    #[builder(default = "u!(0 s)")]
    pub delay: Time,
    #[builder(default = "u!(0 hz)")]
    pub damping: Frequency,
    #[builder(default = "num!(0)")]
    pub phase_deg: Number,
}

impl SineVoltage {
    // F(t) = A sin(2pi f t)
    pub fn sin(amplitude: Voltage, frequency: Frequency) -> Self {
        SineVoltageBuilder::default()
            .vo(u!(0 V))
            .va(amplitude)
            .freq_hz(frequency)
            .build().unwrap()
    }

    // F(t) = A cos(2pi f t)
    pub fn cos(amplitude: Voltage, frequency: Frequency) -> Self {
        let period = frequency.to_period();
        SineVoltageBuilder::default()
            .vo(u!(0 V))
            .va(amplitude)
            .freq_hz(frequency)
            .delay(-0.25 * period)
            .build().unwrap()
    }

    pub fn period(&self) -> Time {
        self.freq_hz.to_period()
    }

    pub fn voltage_at(&self, time: Time) -> Voltage {
        if time < self.delay {
            return self.vo; 
        }

        let td = time - self.delay;
        let envelope = (-self.damping * td).exp();
        let omega = (2.0 * f64::consts::PI) * self.freq_hz;
        let sine = (omega * td + self.phase_deg / 360.).sin();

        self.vo + self.va * envelope * sine
    }

    pub fn to_spice(&self) -> String {
        format!(
            "SIN({} {} {} {} {} {})",
            self.vo,
            self.va,
            self.freq_hz,
            self.delay,
            self.damping,
            self.phase_deg
        )
    }
}