signal_processing 0.3.0

A signal processing 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

use num::{complex::ComplexFloat, Complex};

use crate::{quantities::{Lists, ContainerOrSingle}, transforms::fourier::Dft};

pub trait Dht<T>: Lists<T>
where
    T: ComplexFloat
{
    fn dht(self) -> Self::Mapped<T::Real>;
}

impl<T, L> Dht<T> for L
where
    T: ComplexFloat,
    L: Lists<T>,
    Self: Dft<T>,
    Self::Mapped<Complex<T::Real>>: Lists<Complex<T::Real>, Mapped<T::Real> = Self::Mapped<T::Real>>,
{
    fn dht(self) -> Self::Mapped<T::Real>
    {
        let y = self.dft();
        y.map_into_owned(|y| y.re - y.im)
    }
}

#[cfg(test)]
mod test
{
    use core::f64::consts::TAU;

    use array_math::ArrayOps;
    use linspace::LinspaceArray;

    use crate::{plot, transforms::fourier::Dht};

    #[test]
    fn test()
    {
        const N: usize = 1024;
        const T: f64 = 1.0;
        const F: f64 = 220.0;
        
        let x: [_; N] = ArrayOps::fill(|i| (TAU*F*i as f64/N as f64*T).sin());

        let xf = x.dht();

        let w = (0.0..TAU).linspace_array();

        plot::plot_curves("X(e^jw)", "plots/x_z_dht.png", [&w.zip(xf)])
            .unwrap()
    }
}