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
53
54
55




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

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

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

impl<T, L> Idht<T> for L
where
    T: ComplexFloat,
    L: Lists<T>,
    Self: Idft<T>,
    Self::Mapped<Complex<T::Real>>: Lists<Complex<T::Real>, Mapped<T::Real> = Self::Mapped<T::Real>>,
{
    fn idht(self) -> Self::Mapped<T::Real>
    {
        let y = self.idft();
        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, Idht}};

    #[test]
    fn test()
    {
        const N: usize = 1024;
        const T: f64 = 0.1;
        const F: f64 = 220.0;
        
        let t: [_; N] = (0.0..T).linspace_array();
        let x = t.map(|t| (TAU*F*t).sin());

        let xf = x.dht();
        let y = xf.idht();

        plot::plot_curves("x(t)", "plots/x_t_idht.png", [&t.zip(y.map(|y| y))])
            .unwrap()
    }
}