mrrs 0.1.1

Multi rate filtering tools
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

from scipy.signal import remez
import numpy as np
np.set_printoptions(precision=20)

N = 30
assert N % 2 == 0
assert N % 4 != 0
flt_2 = remez(N//2+1,
              [0, 2*0.2125, 0.5, 0.5],
              [1, 0],
              weight=[1, 1])
# Insert zeros
flt = [item for sublist in zip(flt_2, np.zeros(N)) for item in sublist][:-1]
# Half coefficients
flt = np.array(flt)/2
# Set center tap
flt[N//2] = 0.5
print(len(flt))
print(flt)
print(flt_2[:len(flt_2)//2]/2)