def test_001_simple_receiver(self):
# make sure advanced receiver works like simple receiver
# in case no IC iterations are applied!
reps = 5
alpha = 0.5
M = 127
K = 16
L = 2
taps = filters.get_frequency_domain_filter("rrc", alpha, M, K, L)
data = np.array([], dtype=complex)
ref = np.array([], dtype=complex)
for _11 in range(reps):
d = utils.get_random_qpsk(M * K)
ref = np.append(ref, gfdm_demodulate_block(d, taps, K, M, L))
data = np.append(data, d)
# print data
# print ref
# print "MAXIMUM ref value: ", np.max(abs(ref))
src = blocks.vector_source_c(data)
est_data = np.ones(len(data), dtype=complex)
est_src = blocks.vector_source_c(est_data)
gfdm_constellation = digital.constellation_qpsk().base()
mod = gfdm.advanced_receiver_sb_cc(M, K, L, 0, taps,
gfdm_constellation, np.arange(K), 0)
dst = blocks.vector_sink_c()
self.tb.connect(src, (mod, 0), dst)
self.tb.connect(est_src, (mod, 1))
# set up fg
self.tb.run()
# check data
res = np.array(dst.data())
self.assertComplexTuplesAlmostEqual(ref, res, 4)
def test_003_setIC(self):
ic = 2
timeslots = 9
subcarriers = 32
active_subcarriers = 20
overlap = 2
f_taps = filters.get_frequency_domain_filter('rrc', .5, timeslots,
subcarriers, overlap)
gfdm_constellation = digital.constellation_qpsk().base()
subcarrier_map = get_subcarrier_map(subcarriers, active_subcarriers)
demod = gfdm.advanced_receiver_sb_cc(timeslots, subcarriers, overlap,
64, f_taps, gfdm_constellation,
subcarrier_map, 0)
demod.set_ic(ic)
self.assertEqual(ic, demod.get_ic())
def test_002_t(self):
n_frames = 1
gfdm_var = struct({
"subcarriers": 64,
"timeslots": 9,
"alpha": 0.5,
"overlap": 2,
})
gfdm_constellation = digital.constellation_qpsk().base()
self.f_taps = f_taps = filters.get_frequency_domain_filter(
"rrc", 1.0, gfdm_var.timeslots, gfdm_var.subcarriers,
gfdm_var.overlap)
source_bits = np.random.randint(
0,
len(gfdm_constellation.points()),
n_frames * gfdm_var.timeslots * gfdm_var.subcarriers,
).astype(np.uint8)
self.random_bits = blocks.vector_source_b(source_bits, False)
self.bits_to_symbols = digital.chunks_to_symbols_bc(
(gfdm_constellation.points()), 1)
self.mod = gfdm.simple_modulator_cc(gfdm_var.timeslots,
gfdm_var.subcarriers,
gfdm_var.overlap, f_taps)
self.demod = gfdm.advanced_receiver_sb_cc(
gfdm_var.timeslots,
gfdm_var.subcarriers,
gfdm_var.overlap,
64,
f_taps,
gfdm_constellation,
np.arange(gfdm_var.subcarriers),
0,
)
self.tx_symbols = blocks.vector_sink_c()
self.rx_symbols = blocks.vector_sink_c()
self.tb.connect((self.random_bits, 0), (self.bits_to_symbols, 0))
self.tb.connect((self.bits_to_symbols, 0), (self.tx_symbols, 0))
self.tb.connect((self.bits_to_symbols, 0), (self.mod, 0))
self.tb.connect((self.mod, 0), (self.demod, 0))
self.tb.connect((self.demod, 0), (self.rx_symbols, 0))
self.tb.run()
ref = np.array(self.tx_symbols.data())
res = np.array(self.rx_symbols.data())
# more or less make sure all symbols have their correct sign.
self.assertComplexTuplesAlmostEqual(ref, res, 2)
def test_004_active_subcarriers(self):
n_frames = 1
timeslots = 9
subcarriers = 32
active_subcarriers = 20
overlap = 2
ic_iterations = 64
f_taps = filters.get_frequency_domain_filter("rrc", 0.5, timeslots,
subcarriers, overlap)
gfdm_constellation = digital.constellation_qpsk().base()
print(gfdm_constellation.points())
subcarrier_map = get_subcarrier_map(subcarriers, active_subcarriers)
data = get_random_qpsk(n_frames * timeslots * active_subcarriers)
data *= 2.0
src = blocks.vector_source_c(data)
mapper = gfdm.resource_mapper_cc(timeslots, subcarriers,
active_subcarriers, subcarrier_map,
True)
mod = gfdm.simple_modulator_cc(timeslots, subcarriers, overlap, f_taps)
demod = gfdm.advanced_receiver_sb_cc(
timeslots,
subcarriers,
overlap,
ic_iterations,
f_taps,
gfdm_constellation,
subcarrier_map,
0,
)
demod.set_ic(64)
demapper = gfdm.resource_demapper_cc(timeslots, subcarriers,
active_subcarriers,
subcarrier_map, True)
snk = blocks.vector_sink_c()
self.tb.connect(src, mapper, mod, demod, demapper, snk)
self.tb.run()
res = np.array(snk.data())
print(data[0:10])
print(res[0:10])
self.assertComplexTuplesAlmostEqual(data, res, 1)