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

use core::ops::{Sub, Mul};

use num::complex::ComplexFloat;
use option_trait::Maybe;

use crate::{quantities::{MaybeList, MaybeOwnedList}, systems::{Sos, Tf}, decompositions::SplitNumerDenom, transforms::system::{ToSos, ToTf}};

impl<T1, T2, T3, B1, B2, A1, A2, A3, S1, S2, S3, S4> Sub<Sos<T2, B2, A2, S2>> for Sos<T1, B1, A1, S1>
where
    T1: ComplexFloat,
    T2: ComplexFloat,
    T3: ComplexFloat,
    B1: Maybe<[T1; 3]> + MaybeOwnedList<T1> + Clone,
    B2: Maybe<[T2; 3]> + MaybeOwnedList<T2> + Clone,
    A1: Maybe<[T1; 3]> + MaybeOwnedList<T1> + Clone,
    A2: Maybe<[T2; 3]> + MaybeOwnedList<T2> + Clone,
    A3: Maybe<[T3; 3]> + MaybeOwnedList<T3>,
    S1: MaybeList<Tf<T1, B1, A1>>,
    S2: MaybeList<Tf<T2, B2, A2>>,
    S3: MaybeList<Tf<T3, (), A3>>,
    S4: MaybeList<Tf<T3, [T3; 3], A3>>,
    
    Self: SplitNumerDenom,
    Sos<T2, B2, A2, S2>: SplitNumerDenom,

    <Self as SplitNumerDenom>::OutputNum: ToTf<T1, Vec<T1>, (), (), ()>,
    <Sos<T2, B2, A2, S2> as SplitNumerDenom>::OutputNum: ToTf<T2, Vec<T2>, (), (), ()>,

    Tf<T1, Vec<T1>>: Sub<Tf<T2, Vec<T2>>, Output = Tf<T3, Vec<T3>>>,
    Tf<T3, Vec<T3>>: ToSos<T3, [T3; 3], (), Vec<Tf<T3, [T3; 3]>>, (), ()>,

    <Self as SplitNumerDenom>::OutputDen: Mul<<Sos<T2, B2, A2, S2> as SplitNumerDenom>::OutputDen, Output = Sos<T3, (), A3, S3>>,
    Sos<T3, [T3; 3], (), Vec<Tf<T3, [T3; 3]>>>: Mul<Sos<T3, (), A3, S3>, Output = Sos<T3, [T3; 3], A3, S4>>
{
    type Output = Sos<T3, [T3; 3], A3, S4>;

    fn sub(self, rhs: Sos<T2, B2, A2, S2>) -> Self::Output
    {
        let (n1, d1) = self.split_numer_denom();
        let (n2, d2) = rhs.split_numer_denom();

        let n1 = n1.to_tf((), ());
        let n2 = n2.to_tf((), ());

        let n = n1 - n2;
        let n = n.to_sos((), ());

        let d = d1*d2;
        n*d
    }
}