Step 6a: Make FIR

This tutorial describes how a time-domain filter can be created based on theoretical Z, P and K values. In particular, a process of generating a combined response of the inverse analog and digital anti-aliasing (AA and dAA) filters is shown. <ul>

# inverse of the analog Anti-Aliasing (AA^-1)
dsgn_AA_z = [7680.37]
dsgn_AA_p = [8190, 8190]

def myZPKhz(z, p, k):
    z2pi = [-2 * np.pi * x for x in z]
    p2pi = [-2 * np.pi * x for x in p]
    g = k * np.prod(p2pi) / np.prod(z2pi)
    return sg.lti(z2pi, p2pi, g)

dsgn_AA = myZPKhz(dsgn_AA_z, dsgn_AA_p, 1.0)

Output

KAGRA/Subgroups/CAL/GstLAL/tutorials/step6a (last edited 2018-08-30 10:17:41 by darkhan.tuyenbayev)