def test_8_4_3(self):
order = 8
osr = 4
h_inf = 3
z, p, k = synthesizeChebyshevNTF(order, osr, 1, h_inf)
e_k = 1
e_z = [
716.675171843707e-003 + 697.407125044473e-003j,
716.675171843707e-003 - 697.407125044473e-003j,
987.024523733349e-003 + 160.569578528934e-003j,
987.024523733349e-003 - 160.569578528934e-003j,
787.924470008031e-003 + 615.771897347195e-003j,
787.924470008031e-003 - 615.771897347195e-003j,
899.412094416429e-003 + 437.101686587289e-003j,
899.412094416429e-003 - 437.101686587289e-003j
]
e_p = [
626.054825241768e-003 + 122.563529859845e-003j,
626.054825241768e-003 - 122.563529859845e-003j,
603.812988962951e-003 + 354.087262468194e-003j,
603.812988962951e-003 - 354.087262468194e-003j,
583.053683660931e-003 + 552.205015613472e-003j,
583.053683660931e-003 - 552.205015613472e-003j,
599.974045233308e-003 + 710.003702019636e-003j,
599.974045233308e-003 - 710.003702019636e-003j
]
e_z = cplxpair(e_z)
e_p = cplxpair(e_p)
np.testing.assert_almost_equal(k, e_k, 4)
np.testing.assert_almost_equal(z, e_z, 4)
np.testing.assert_almost_equal(p, e_p, 4)
def test_opt3(self):
z, p, k = synthesizeNTF(opt=3)
e_k = 1
e_z = [1.0000, 0.9993 - 0.0382j, 0.9993 + 0.0382j]
e_z = cplxpair(e_z)
e_p = [0.6692, 0.7652 - 0.2795j, 0.7652 + 0.2795j]
e_p = cplxpair(e_p)
np.testing.assert_almost_equal(k, e_k, 6)
np.testing.assert_almost_equal(z, e_z, 4)
np.testing.assert_almost_equal(p, e_p, 4)
def test_uc_zeros(self):
order = 7
OSR = 32
zeros1 = spread_fir_uc_zeros(order, OSR,
quantization_noise_gain,
cf_kwargs={'bounds': (0, 0.5/OSR)})
zeros1 = cplxpair(zeros1)
zeros2 = ds_optzeros(order)
zeros2 = np.exp(1j*np.pi*zeros2/OSR)
zeros2 = cplxpair(zeros2)
np.testing.assert_almost_equal(zeros1, zeros2, 4)
def test_uc_zeros(self):
order = 7
OSR = 32
zeros1 = spread_fir_uc_zeros(order,
OSR,
quantization_noise_gain,
cf_kwargs={'bounds': (0, 0.5 / OSR)})
zeros1 = cplxpair(zeros1)
zeros2 = ds_optzeros(order)
zeros2 = np.exp(1j * np.pi * zeros2 / OSR)
zeros2 = cplxpair(zeros2)
np.testing.assert_almost_equal(zeros1, zeros2, 4)
def test_BP6(self):
z, p, k = synthesizeNTF(order=6, f0=0.3)
e_k = 1
e_z = np.concatenate(((-0.3090 + 0.9511j) * np.ones(3),
(-0.3090 - 0.9511j) * np.ones(3)))
e_z = cplxpair(e_z)
e_p = [
-0.4238 - 0.7906j, -0.4238 + 0.7906j, -0.2480 - 0.7802j,
-0.2480 + 0.7802j, -0.1290 - 0.8985j, -0.1290 + 0.8985j
]
e_p = cplxpair(e_p)
np.testing.assert_almost_equal(k, e_k, 6)
np.testing.assert_almost_equal(z, e_z, 4)
np.testing.assert_almost_equal(p, e_p, 4)
def test_default(self):
z, p, k = synthesizeNTF()
e_k = 1
e_z = np.ones(3)
e_p = [0.7654 - 0.2793j, 0.7654 + 0.2793j, 0.6694]
e_p = cplxpair(e_p)
np.testing.assert_almost_equal(k, e_k, 6)
np.testing.assert_almost_equal(z, e_z, 6)
np.testing.assert_almost_equal(p, e_p, 4)
def test_clans(self):
"""Test function for clans"""
e_z = [1,
0.99604531 + 0.08884669j,
0.99604531 - 0.08884669j,
0.99860302 + 0.05283948j,
0.99860302 - 0.05283948j]
e_p = [0.41835234,
0.48922229 + 0.1709716j,
0.48922229 - 0.1709716j,
0.65244884 + 0.38172239j,
0.65244884 - 0.38172239j]
e_k = 1
e_z = cplxpair(e_z)
e_p = cplxpair(e_p)
z, p, k = clans(5, 32, 5, .95, 1)
z = cplxpair(z)
p = cplxpair(p)
np.testing.assert_almost_equal(k, e_k, 6)
np.testing.assert_almost_equal(z, e_z, 6)
np.testing.assert_almost_equal(p, e_p, 6)
def test_5_32_5_095_1(self):
z, p, k = clans(5, 32, 5, 0.95, 1)
e_k = 1
e_z = [
1.00000000000000e+000,
998.603018798634e-003 + 52.8394818885978e-003j,
998.603018798634e-003 - 52.8394818885978e-003j,
996.045312843931e-003 + 88.8466924631155e-003j,
996.045312843931e-003 - 88.8466924631155e-003j
]
e_p = [
418.352260015271e-003,
489.222137893601e-003 + 170.971685170035e-003j,
489.222137893601e-003 - 170.971685170035e-003j,
652.448968088646e-003 + 381.722332008184e-003j,
652.448968088646e-003 - 381.722332008184e-003j
]
e_z = cplxpair(e_z)
e_p = cplxpair(e_p)
z = cplxpair(z)
p = cplxpair(p)
np.testing.assert_almost_equal(k, e_k, 4)
np.testing.assert_almost_equal(z, e_z, 4)
np.testing.assert_almost_equal(p, e_p, 4)