kshana 0.7.0

Open hybrid quantum/classical PNT performance simulator
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
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101

// SPDX-License-Identifier: Apache-2.0
use crate::scenario::GnssState;
use crate::types::Seconds;

/// Disciplined to truth while GNSS is nominal; open-loop during outage,
/// predicting phase with calibrated frequency AND aging (quadratic):
///   predicted = sync_phase + f_s*dt + 0.5*drift*dt^2
/// Deterministic offset+aging are removed, so the reported residual is the
/// stochastic (white-FM + random-walk-FM) holdover error — the fundamental limit.
pub struct HoldoverEstimator {
    f_s: f64,
    drift_est: f64,
    last_sync_t: Seconds,
    last_sync_phase: Seconds,
    synced: bool,
}

impl HoldoverEstimator {
    pub fn new() -> Self {
        Self {
            f_s: 0.0,
            drift_est: 0.0,
            last_sync_t: 0.0,
            last_sync_phase: 0.0,
            synced: false,
        }
    }

    pub fn timing_error(
        &mut self,
        t: Seconds,
        true_phase: Seconds,
        det_freq: f64,
        drift: f64,
        gnss: GnssState,
    ) -> Seconds {
        match gnss {
            GnssState::Nominal => {
                self.f_s = det_freq;
                self.drift_est = drift;
                self.last_sync_t = t;
                self.last_sync_phase = true_phase;
                self.synced = true;
                0.0
            }
            _ => {
                if !self.synced {
                    return 0.0;
                }
                let dt = t - self.last_sync_t;
                let predicted =
                    self.last_sync_phase + self.f_s * dt + 0.5 * self.drift_est * dt * dt;
                true_phase - predicted
            }
        }
    }
}

impl Default for HoldoverEstimator {
    fn default() -> Self {
        Self::new()
    }
}

#[cfg(test)]
mod tests {
    use super::*;
    use crate::scenario::GnssState::*;

    #[test]
    fn synced_then_perfect_prediction_is_zero_error() {
        let mut e = HoldoverEstimator::new();
        assert_eq!(e.timing_error(0.0, 0.0, 1e-9, 0.0, Nominal), 0.0);
        let true_phase = 1e-9 * 10.0;
        assert!((e.timing_error(10.0, true_phase, 1e-9, 0.0, Denied)).abs() < 1e-18);
    }

    #[test]
    fn quadratic_aging_is_removed() {
        // sync t=0: phase 0, f_s=1e-9, drift=2e-13.
        // true deterministic phase at t=10: 1e-9*10 + 2e-13*10*10/2.
        let mut e = HoldoverEstimator::new();
        e.timing_error(0.0, 0.0, 1e-9, 2e-13, Nominal);
        let true_phase = 1e-9 * 10.0 + 2e-13 * 10.0 * 10.0 / 2.0;
        let err = e.timing_error(10.0, true_phase, 1.002e-9, 2e-13, Denied);
        assert!(err.abs() < 1e-22, "err={err}");
    }

    #[test]
    fn residual_drift_during_outage_is_reported() {
        let mut e = HoldoverEstimator::new();
        e.timing_error(0.0, 0.0, 0.0, 0.0, Nominal);
        assert!((e.timing_error(5.0, 2e-7, 0.0, 0.0, Denied) - 2e-7).abs() < 1e-18);
    }

    #[test]
    fn never_synced_returns_zero() {
        let mut e = HoldoverEstimator::new();
        assert_eq!(e.timing_error(3.0, 1.0, 0.0, 0.0, Denied), 0.0);
    }
}