orion-sdr 0.0.12

DSP/SDR block library targeting HF-to-UHF, satellites, and Python bindings. Roadmap inside.
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
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121

// src/demodulate/am.rs
use num_complex::Complex32 as C32;
use crate::core::{Block, WorkReport};
use crate::dsp::LpDcCascade;

#[derive(Debug, Clone, Copy)]
pub enum Envelope {
    /// env = sqrt(I^2 + Q^2) after IIR LP (highest fidelity)
    PowerSqrt,
    /// env ≈ a*|I| + b*|Q| (very fast, small amplitude error)
    AbsApprox { k1: f32, k2: f32 },
}

#[derive(Debug, Clone)]
pub struct AmEnvelopeDemod {
    filt: LpDcCascade,
    method: Envelope,
}

impl AmEnvelopeDemod {
    pub fn new(fs: f32, audio_bw_hz: f32) -> Self {
        // 4th-order LR at ~0.9*BW gives clean audio passband.
        Self { filt: LpDcCascade::design(fs, audio_bw_hz * 0.9, 2.0), method: Envelope::PowerSqrt }
    }

    /// Optional: switch to abs-approx envelope (faster, small error).
    pub fn with_abs_approx(mut self, k1: f32, k2: f32) -> Self {
        self.method = Envelope::AbsApprox { k1, k2 };
        self
    }
}

impl Block for AmEnvelopeDemod {
    type In  = C32; // complex baseband AM
    type Out = f32; // audio

    #[inline(always)]
    fn process(&mut self, input: &[Self::In], output: &mut [Self::Out]) -> WorkReport {
        let n = input.len().min(output.len());
        let mut i = 0;
        let nn = n & !3;

        match self.method {
            Envelope::PowerSqrt => {
                while i < nn {
                    // 0
                    let z0 = input[i];
                    let p0 = z0.re.mul_add(z0.re, z0.im * z0.im);
                    output[i] = self.filt.process_mapped(p0, f32::sqrt);

                    // 1
                    let z1 = input[i+1];
                    let p1 = z1.re.mul_add(z1.re, z1.im * z1.im);
                    output[i+1] = self.filt.process_mapped(p1, f32::sqrt);

                    // 2
                    let z2 = input[i+2];
                    let p2 = z2.re.mul_add(z2.re, z2.im * z2.im);
                    output[i+2] = self.filt.process_mapped(p2, f32::sqrt);

                    // 3
                    let z3 = input[i+3];
                    let p3 = z3.re.mul_add(z3.re, z3.im * z3.im);
                    output[i+3] = self.filt.process_mapped(p3, f32::sqrt);

                    i += 4;
                }
                while i < n {
                    let z = input[i];
                    let p = z.re.mul_add(z.re, z.im * z.im);
                    output[i] = self.filt.process_mapped(p, f32::sqrt);
                    i += 1;
                }
            }
            Envelope::AbsApprox { k1, k2 } => {
                while i < nn {
                    // 0
                    let z0 = input[i];
                    let a0 = z0.re.abs();
                    let b0 = z0.im.abs();
                    // e0 = k1*a0 + k2*b0, fused as e0 = k1.mul_add(a0, k2 * b0)
                    let e0 = k1.mul_add(a0, k2 * b0);
                    output[i] = self.filt.process(e0);

                    // 1
                    let z1 = input[i+1];
                    let a1 = z1.re.abs();
                    let b1 = z1.im.abs();
                    let e1 = k1.mul_add(a1, k2 * b1);
                    output[i+1] = self.filt.process(e1);

                    // 2
                    let z2 = input[i+2];
                    let a2 = z2.re.abs();
                    let b2 = z2.im.abs();
                    let e2 = k1.mul_add(a2, k2 * b2);
                    output[i+2] = self.filt.process(e2);

                    // 3
                    let z3 = input[i+3];
                    let a3 = z3.re.abs();
                    let b3 = z3.im.abs();
                    let e3 = k1.mul_add(a3, k2 * b3);
                    output[i+3] = self.filt.process(e3);

                    i += 4;
                }
                while i < n {
                    let z = input[i];
                    let a = z.re.abs();
                    let b = z.im.abs();
                    let e = k1.mul_add(a, k2 * b);
                    output[i] = self.filt.process(e);
                    i += 1;
                }
            }
        }

        WorkReport { in_read: n, out_written: n }
    }
}