def test02(self):
# Test QPSK sync
M = 4
theta = 0
loop_bw = cmath.pi/100.0
fmin = -0.5
fmax = 0.5
mu = 0.5
gain_mu = 0.01
omega = 2
gain_omega = 0.001
omega_rel = 0.001
self.test = digital.mpsk_receiver_cc(M, theta, loop_bw,
fmin, fmax, mu, gain_mu,
omega, gain_omega,
omega_rel)
data = 10000*[complex( 0.707, 0.707),
complex(-0.707, 0.707),
complex(-0.707, -0.707),
complex( 0.707, -0.707)]
data = [0.5*d for d in data]
self.src = gr.vector_source_c(data, False)
self.snk = gr.vector_sink_c()
# pulse shaping interpolation filter
nfilts = 32
excess_bw = 0.35
ntaps = 11 * int(omega*nfilts)
rrc_taps0 = filter.firdes.root_raised_cosine(
nfilts, nfilts, 1.0, excess_bw, ntaps)
rrc_taps1 = filter.firdes.root_raised_cosine(
1, omega, 1.0, excess_bw, 11*omega)
self.rrc0 = filter.pfb_arb_resampler_ccf(omega, rrc_taps0)
self.rrc1 = filter.fir_filter_ccf(1, rrc_taps1)
self.tb.connect(self.src, self.rrc0, self.rrc1, self.test, self.snk)
self.tb.run()
expected_result = 10000*[complex(-0.5, +0.0), complex(+0.0, -0.5),
complex(+0.5, +0.0), complex(+0.0, +0.5)]
# get data after a settling period
dst_data = self.snk.data()[200:]
# Only compare last Ncmp samples
Ncmp = 1000
len_e = len(expected_result)
len_d = len(dst_data)
expected_result = expected_result[len_e - Ncmp - 1:-1]
dst_data = dst_data[len_d - Ncmp:]
#for e,d in zip(expected_result, dst_data):
# print "{0:+.02f} {1:+.02f}".format(e, d)
self.assertComplexTuplesAlmostEqual(expected_result, dst_data, 1)
def test01(self):
# Test BPSK sync
M = 2
theta = 0
loop_bw = cmath.pi / 100.0
fmin = -0.5
fmax = 0.5
mu = 0.5
gain_mu = 0.01
omega = 2
gain_omega = 0.001
omega_rel = 0.001
self.test = digital.mpsk_receiver_cc(M, theta, loop_bw, fmin, fmax, mu,
gain_mu, omega, gain_omega,
omega_rel)
data = 10000 * [complex(1, 0), complex(-1, 0)]
#data = [2*random.randint(0,1)-1 for x in xrange(10000)]
self.src = blocks.vector_source_c(data, False)
self.snk = blocks.vector_sink_c()
# pulse shaping interpolation filter
nfilts = 32
excess_bw = 0.35
ntaps = 11 * int(omega * nfilts)
rrc_taps0 = filter.firdes.root_raised_cosine(nfilts, nfilts, 1.0,
excess_bw, ntaps)
rrc_taps1 = filter.firdes.root_raised_cosine(1, omega, 1.0, excess_bw,
11 * omega)
self.rrc0 = filter.pfb_arb_resampler_ccf(omega, rrc_taps0)
self.rrc1 = filter.fir_filter_ccf(1, rrc_taps1)
self.tb.connect(self.src, self.rrc0, self.rrc1, self.test, self.snk)
self.tb.run()
expected_result = [0.5 * d for d in data]
dst_data = self.snk.data()
# Only compare last Ncmp samples
Ncmp = 1000
len_e = len(expected_result)
len_d = len(dst_data)
expected_result = expected_result[len_e - Ncmp - 1:-1]
dst_data = dst_data[len_d - Ncmp:]
#for e,d in zip(expected_result, dst_data):
# print "{0:+.02f} {1:+.02f}".format(e, d)
self.assertComplexTuplesAlmostEqual(expected_result, dst_data, 1)
def test_fir_filter_ccf_002(self):
src_data = 40*[1+1j, 2+2j, 3+3j, 4+4j]
expected_data = ((0.5+0.5j), (5.5+5.5j), (10.5+10.5j), (15.5+15.5j),
(20.5+20.5j), (25.5+25.5j), (30.5+30.5j), (35.5+35.5j),
(40.5+40.5j), (45.5+45.5j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j))
src = gr.vector_source_c(src_data)
op = filter.fir_filter_ccf(4, 20*[0.5, 0.5])
dst = gr.vector_sink_c()
self.tb.connect(src, op, dst)
self.tb.run()
result_data = dst.data()
self.assertComplexTuplesAlmostEqual(expected_data, result_data, 5)
def test_fir_filter_ccf_002(self):
src_data = 40*[1+1j, 2+2j, 3+3j, 4+4j]
expected_data = ((0.5+0.5j), (5.5+5.5j), (10.5+10.5j), (15.5+15.5j),
(20.5+20.5j), (25.5+25.5j), (30.5+30.5j), (35.5+35.5j),
(40.5+40.5j), (45.5+45.5j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j))
src = blocks.vector_source_c(src_data)
op = filter.fir_filter_ccf(4, 20*[0.5, 0.5])
dst = blocks.vector_sink_c()
self.tb.connect(src, op, dst)
self.tb.run()
result_data = dst.data()
self.assertComplexTuplesAlmostEqual(expected_data, result_data, 5)
def test_fir_filter_ccf_001(self):
src_data = 40*[1+1j, 2+2j, 3+3j, 4+4j]
expected_data = ((0.5+0.5j), (1.5+1.5j), (3+3j), (5+5j), (5.5+5.5j),
(6.5+6.5j), (8+8j), (10+10j), (10.5+10.5j), (11.5+11.5j),
(13+13j), (15+15j), (15.5+15.5j), (16.5+16.5j), (18+18j),
(20+20j), (20.5+20.5j), (21.5+21.5j), (23+23j), (25+25j),
(25.5+25.5j), (26.5+26.5j), (28+28j), (30+30j), (30.5+30.5j),
(31.5+31.5j), (33+33j), (35+35j), (35.5+35.5j), (36.5+36.5j),
(38+38j), (40+40j), (40.5+40.5j), (41.5+41.5j), (43+43j),
(45+45j), (45.5+45.5j), (46.5+46.5j), (48+48j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j))
src = gr.vector_source_c(src_data)
op = filter.fir_filter_ccf(1, 20*[0.5, 0.5])
dst = gr.vector_sink_c()
self.tb.connect(src, op, dst)
self.tb.run()
result_data = dst.data()
self.assertComplexTuplesAlmostEqual(expected_data, result_data, 5)
def test_fir_filter_ccf_001(self):
src_data = 40*[1+1j, 2+2j, 3+3j, 4+4j]
expected_data = ((0.5+0.5j), (1.5+1.5j), (3+3j), (5+5j), (5.5+5.5j),
(6.5+6.5j), (8+8j), (10+10j), (10.5+10.5j), (11.5+11.5j),
(13+13j), (15+15j), (15.5+15.5j), (16.5+16.5j), (18+18j),
(20+20j), (20.5+20.5j), (21.5+21.5j), (23+23j), (25+25j),
(25.5+25.5j), (26.5+26.5j), (28+28j), (30+30j), (30.5+30.5j),
(31.5+31.5j), (33+33j), (35+35j), (35.5+35.5j), (36.5+36.5j),
(38+38j), (40+40j), (40.5+40.5j), (41.5+41.5j), (43+43j),
(45+45j), (45.5+45.5j), (46.5+46.5j), (48+48j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j),
(50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j))
src = blocks.vector_source_c(src_data)
op = filter.fir_filter_ccf(1, 20*[0.5, 0.5])
dst = blocks.vector_sink_c()
self.tb.connect(src, op, dst)
self.tb.run()
result_data = dst.data()
self.assertComplexTuplesAlmostEqual(expected_data, result_data, 5)
