def test_qpsk_channel(self):
upper_bound = tuple(50.0 * numpy.ones((self.num_data, )))
lower_bound = tuple(0.0 * numpy.zeros((self.num_data, )))
self.cons = cons = digital.constellation_qpsk().base()
self.data = data = [
random.randrange(len(cons.points())) for x in range(self.num_data)
]
self.symbols = symbols = numpy.squeeze(
[cons.map_to_points_v(i) for i in data])
chan = channels.channel_model(noise_voltage=0.1,
frequency_offset=0.0,
epsilon=1.0,
taps=[1.0 + 0.0j],
noise_seed=0,
block_tags=False)
evm = digital.meas_evm_cc(cons, digital.evm_measurement_t_EVM_PERCENT)
vso = blocks.vector_source_c(symbols, False, 1, [])
mc = blocks.multiply_const_cc(3.0 + 2.0j)
vsi = blocks.vector_sink_f()
self.tb.connect(vso, chan, evm, vsi)
self.tb.run()
# check data
output_data = vsi.data()
self.assertLess(output_data, upper_bound)
self.assertGreater(output_data, lower_bound)
def __init__(self, constellation=[-1, 1]):
gr.top_block.__init__(self, "Channel Model Gui")
##################################################
# Parameters
##################################################
self.constellation = constellation
##################################################
# Variables
##################################################
self.const_points = const_points = constellation
self.taps = taps = [1.0, 0.25-0.25j, 0.50 + 0.10j, 0.3 + 0.2j]
self.samp_rate = samp_rate = 8e6
self.noise_level = noise_level = 0.3
self.frequency = frequency = 2000000
self.const_type = const_type = 1
self.const_dist = const_dist = digital.constellation_calcdist(const_points,[],0,1)
self.const = const = digital.constellation_qpsk()
##################################################
# Blocks
##################################################
self.digital_constellation_decoder_cb_0_0 = digital.constellation_decoder_cb(const_dist.base())
self.digital_constellation_decoder_cb_0 = digital.constellation_decoder_cb(const_dist.base())
self.digital_chunks_to_symbols_xx_0 = digital.chunks_to_symbols_bc((const_dist.points()), 1)
self.custom_ber_0 = custom_ber(
n_bits=10000,
bits_per_symbol=3,
)
self.channels_channel_model_0_0 = channels.channel_model(
noise_voltage=0,
frequency_offset=0,
epsilon=1,
taps=([1,0,0,0]),
noise_seed=0,
block_tags=False
)
self.channels_channel_model_0 = channels.channel_model(
noise_voltage=noise_level,
frequency_offset=0,
epsilon=1,
taps=(taps),
noise_seed=0,
block_tags=False
)
self.blocks_vector_sink_x_1 = blocks.vector_sink_f(1)
self.analog_random_source_x_0 = blocks.vector_source_b(map(int, numpy.random.randint(0, len(const_points), 10000)), False)
##################################################
# Connections
##################################################
self.connect((self.custom_ber_0, 0), (self.blocks_vector_sink_x_1, 0))
self.connect((self.channels_channel_model_0, 0), (self.digital_constellation_decoder_cb_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0, 0), (self.channels_channel_model_0, 0))
self.connect((self.analog_random_source_x_0, 0), (self.digital_chunks_to_symbols_xx_0, 0))
self.connect((self.digital_constellation_decoder_cb_0_0, 0), (self.custom_ber_0, 0))
self.connect((self.digital_constellation_decoder_cb_0, 0), (self.custom_ber_0, 1))
self.connect((self.digital_chunks_to_symbols_xx_0, 0), (self.channels_channel_model_0_0, 0))
self.connect((self.channels_channel_model_0_0, 0), (self.digital_constellation_decoder_cb_0_0, 0))
def test_002_static_wo_tags (self):
fft_len = 8
# 4 5 6 7 0 1 2 3
tx_data = [-1, -1, 1, 2, -1, 3, 0, -1, # 0
-1, -1, 0, 2, -1, 2, 0, -1, # 8
-1, -1, 3, 0, -1, 1, 0, -1, # 16 (Pilot symbols)
-1, -1, 1, 1, -1, 0, 2, -1] # 24
cnst = digital.constellation_qpsk()
tx_signal = [cnst.map_to_points_v(x)[0] if x != -1 else 0 for x in tx_data]
occupied_carriers = ((1, 2, 6, 7),)
pilot_carriers = ((), (), (1, 2, 6, 7), ())
pilot_symbols = (
[], [], [cnst.map_to_points_v(x)[0] for x in (1, 0, 3, 0)], []
)
equalizer = digital.ofdm_equalizer_static(fft_len, occupied_carriers, pilot_carriers, pilot_symbols)
channel = [
0, 0, 1, 1, 0, 1, 1, 0,
0, 0, 1, 1, 0, 1, 1, 0, # These coefficients will be rotated slightly...
0, 0, 1j, 1j, 0, 1j, 1j, 0, # Go crazy here!
0, 0, 1j, 1j, 0, 1j, 1j, 0 # ...and again here.
]
for idx in range(fft_len, 2*fft_len):
channel[idx] = channel[idx-fft_len] * numpy.exp(1j * .1 * numpy.pi * (numpy.random.rand()-.5))
idx2 = idx+2*fft_len
channel[idx2] = channel[idx2] * numpy.exp(1j * 0 * numpy.pi * (numpy.random.rand()-.5))
src = blocks.vector_source_c(numpy.multiply(tx_signal, channel), False, fft_len)
sink = blocks.vector_sink_c(fft_len)
eq = digital.ofdm_frame_equalizer_vcvc(equalizer.base(), 0, "", False, 4)
self.tb.connect(src, eq, sink)
self.tb.run ()
rx_data = [cnst.decision_maker_v((x,)) if x != 0 else -1 for x in sink.data()]
self.assertEqual(tx_data, rx_data)
def test_001_simple_receiver(self):
# make sure advanced receiver works like simple receiver in case no IC iterations are applied!
reps = 5
alpha = .5
M = 127
K = 16
L = 2
taps = filters.get_frequency_domain_filter('rrc', alpha, M, K, L)
data = np.array([], dtype=np.complex)
ref = np.array([], dtype=np.complex)
for i in xrange(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=np.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_001_simple_receiver(self):
# make sure advanced receiver works like simple receiver in case no IC iterations are applied!
reps = 5
alpha = .5
M = 127
K = 16
L = 2
taps = filters.get_frequency_domain_filter('rrc', alpha, M, K, L)
data = np.array([], dtype=np.complex)
ref = np.array([], dtype=np.complex)
for i in xrange(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=np.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))
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_active_subcarriers(self):
n_frames = 1
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)
data = get_random_qpsk(n_frames * timeslots * active_subcarriers)
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,
64, f_taps, gfdm_constellation,
subcarrier_map, 0)
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())
self.assertComplexTuplesAlmostEqual(data, res, 2)
def get_constellation(bps):
map = {
1: digital.constellation_bpsk(),
2: digital.constellation_qpsk(),
3: digital.constellation_8psk()
}
return map[bps]
def get_constellation_from_string(self, const_string):
self.const = {
'bpsk': digital.constellation_bpsk().base(),
'qpsk': digital.constellation_qpsk().base(),
'8psk': digital.constellation_8psk().base(),
'16qam': digital.constellation_16qam().base()
}.get(const_string, digital.constellation_bpsk().base())
def __init__(self,
constellation_points=_def_constellation_points,
gray_coded=_def_gray_coded,
*args,
**kwargs):
"""
Hierarchical block for RRC-filtered QPSK modulation.
The input is a byte stream (unsigned char) and the
output is the complex modulated signal at baseband.
See generic_mod block for list of parameters.
"""
constellation_points = _def_constellation_points
constellation = digital_swig.constellation_qpsk()
if constellation_points != 4:
raise ValueError("QPSK can only have 4 constellation points.")
if not gray_coded:
raise ValueError(
"This QPSK mod/demod works only for gray-coded constellations."
)
super(qpsk_mod, self).__init__(constellation=constellation,
gray_coded=gray_coded,
*args,
**kwargs)
def test_002_t(self):
n_frames = 1
self.gfdm_var = gfdm_var = struct({'subcarriers': 64, 'timeslots': 9, 'alpha': 0.5, 'overlap': 2,})
self.gfdm_constellation = 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)
self.random_bits = blocks.vector_source_b(map(int, np.random.randint(0, len(gfdm_constellation.points()),
n_frames * gfdm_var.timeslots * gfdm_var.subcarriers)),
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))
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 __init__(self, ipp1="127.0.0.1", ipp2="127.0.0.1", ipp3="127.0.0.1", ipp4="127.0.0.1", iptx="127.0.0.1", samp_rate=10000):
gr.top_block.__init__(self, "OFDM Rx")
##################################################
# Parameters
##################################################
self.ipp1 = ipp1
self.ipp2 = ipp2
self.ipp3 = ipp3
self.ipp4 = ipp4
self.iptx = iptx
self.samp_rate = samp_rate
##################################################
# Variables
##################################################
self.pilot_symbols = pilot_symbols = ((1, 1, 1, -1,),)
self.pilot_carriers = pilot_carriers = ((-21, -7, 7, 21,),)
self.payload_mod = payload_mod = digital.constellation_qpsk()
self.packet_length_tag_key = packet_length_tag_key = "packet_len"
self.occupied_carriers = occupied_carriers = (range(-26, -21) + range(-20, -7) + range(-6, 0) + range(1, 7) + range(8, 21) + range(22, 27),)
self.length_tag_key = length_tag_key = "frame_len"
self.header_mod = header_mod = digital.constellation_bpsk()
self.fft_len = fft_len = 64
self.sync_word2 = sync_word2 = [0j, 0j, 0j, 0j, 0j, 0j, (-1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1 +0j), (1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), 0j, (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), 0j, 0j, 0j, 0j, 0j]
self.sync_word1 = sync_word1 = [0., 0., 0., 0., 0., 0., 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 0., 0., 0., 0., 0.]
self.payload_equalizer = payload_equalizer = digital.ofdm_equalizer_simpledfe(fft_len, payload_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols, 1)
self.packet_len = packet_len = 96
self.header_formatter = header_formatter = digital.packet_header_ofdm(occupied_carriers, n_syms=1, len_tag_key=packet_length_tag_key, frame_len_tag_key=length_tag_key, bits_per_header_sym=header_mod.bits_per_symbol(), bits_per_payload_sym=payload_mod.bits_per_symbol(), scramble_header=False)
self.header_equalizer = header_equalizer = digital.ofdm_equalizer_simpledfe(fft_len, header_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols)
##################################################
# Blocks
##################################################
self.zeromq_push_sink_0_0_1_0 = zeromq.push_sink(gr.sizeof_gr_complex, 64, "tcp://"+ ipp2 + ":55521", 100, True)
self.zeromq_push_sink_0_0_1 = zeromq.push_sink(gr.sizeof_gr_complex, 64, "tcp://"+ ipp2 + ":55520", 100, True)
self.zeromq_pull_source_0_0_0 = zeromq.pull_source(gr.sizeof_char, 1, "tcp://"+ ipp1 + ":55511", 100, True)
self.zeromq_pull_source_0_0 = zeromq.pull_source(gr.sizeof_gr_complex, 1, "tcp://"+ ipp1 + ":55510", 100, True)
self.zeromq_pull_msg_source_0 = zeromq.pull_msg_source("tcp://"+ ipp3 + ":55530", 100)
self.digital_header_payload_demux_0 = digital.header_payload_demux(
3,
fft_len,
fft_len/4,
length_tag_key,
"",
True,
gr.sizeof_gr_complex,
"rx_time",
samp_rate,
(),
)
##################################################
# Connections
##################################################
self.msg_connect((self.zeromq_pull_msg_source_0, 'out'), (self.digital_header_payload_demux_0, 'header_data'))
self.connect((self.digital_header_payload_demux_0, 0), (self.zeromq_push_sink_0_0_1, 0))
self.connect((self.digital_header_payload_demux_0, 1), (self.zeromq_push_sink_0_0_1_0, 0))
self.connect((self.zeromq_pull_source_0_0, 0), (self.digital_header_payload_demux_0, 0))
self.connect((self.zeromq_pull_source_0_0_0, 0), (self.digital_header_payload_demux_0, 1))
def test_001_t (self):
"""AGC on random noisy QPSK symbols"""
# Parameters
snr_db = 0.0
const_mult = 1.5
agc_rate = 1e-4
agc_ref = 1.0
agc_gain = 1.0
N_last = 1000 # analyze the last symbols
# Constants
n_symbols = int(10*(1/agc_rate))
noise_v = 1/sqrt((10**(float(snr_db)/10)))
rndm = random.Random()
# Input data
in_vec = tuple([rndm.randint(0,1) for i in range(0, n_symbols)])
# Flowgraph
src = blocks.vector_source_b(in_vec)
pack = blocks.repack_bits_bb(1, 2, "", False, gr.GR_MSB_FIRST)
const = digital.constellation_qpsk().base()
cmap = digital.chunks_to_symbols_bc(const.points())
mult = blocks.multiply_const_cc(const_mult)
nadder = blocks.add_cc()
noise = analog.noise_source_c(analog.GR_GAUSSIAN, noise_v, 0)
agc = blocksat.agc_cc(agc_rate, agc_ref, agc_gain)
snk = blocks.vector_sink_c()
snk2 = blocks.vector_sink_c()
# Reference AGC approach
rms_cf = blocks.rms_cf(0.0001)
f2c = blocks.float_to_complex()
div = blocks.divide_cc()
snk3 = blocks.vector_sink_c()
self.tb.connect(src, pack, cmap, mult)
self.tb.connect(mult, (nadder, 0))
self.tb.connect(noise, (nadder, 1))
self.tb.connect(nadder, agc, snk)
self.tb.connect(nadder, snk2)
self.tb.connect(nadder, (div, 0))
self.tb.connect(nadder, rms_cf, (f2c, 0), (div, 1))
self.tb.connect(div, snk3)
self.tb.run()
# Collect results
agc_syms = snk.data()
pre_agc_syms = snk2.data()
ref_agc_syms = snk3.data()
rms_agc = self.rms(agc_syms, N_last)
rms_pre_agc = self.rms(pre_agc_syms, N_last)
rms_agc_ref = self.rms(ref_agc_syms, N_last)
print('RMS before AGC: %f' %(rms_pre_agc))
print('RMS after AGC: %f' %(rms_agc))
print('RMS after reference AGC: %f' %(rms_agc_ref))
# Check results
self.assertAlmostEqual(rms_agc, 1.0, places=1)
def test_003_t(self):
"""Encoder to Decoder - noisy QPSK w/ soft constellation de-mapper"""
# Parameters
N = 18444
K = 6144
pct_en = False
n_ite = 6
M = 4
snr_db = SNR_DB
max_len = 10 * K
flip_llrs = False
sym_rate = 2 * 156250
# NOTE: just like in the previous example, flipping of the LLRs is
# necessary here.
# Constants
noise_v = 1 / math.sqrt((10**(float(snr_db) / 10)))
rndm = random.Random()
# Input data
in_vec = tuple([rndm.randint(0, 1) for i in range(0, max_len)])
# Flowgraph
src = blocks.vector_source_b(in_vec)
enc = blocksattx.turbo_encoder(K, pct_en)
pack = blocks.repack_bits_bb(1, 2, "", False, gr.GR_MSB_FIRST)
const = digital.constellation_qpsk().base()
cmap = digital.chunks_to_symbols_bc(const.points())
nadder = blocks.add_cc()
noise = analog.noise_source_c(analog.GR_GAUSSIAN, noise_v, 0)
cdemap = blocksat.soft_decoder_cf(M, noise_v)
dec = blocksat.turbo_decoder(K, pct_en, n_ite, flip_llrs)
snk = blocks.vector_sink_b()
snk_sym = blocks.vector_sink_c()
snk_llr = blocks.vector_sink_f()
self.tb.connect(src, enc, pack, cmap)
self.tb.connect(cmap, (nadder, 0))
self.tb.connect(noise, (nadder, 1))
self.tb.connect(nadder, cdemap, dec, snk)
self.tb.connect(cmap, snk_sym)
self.tb.connect(cdemap, snk_llr)
self.tb.run()
# Collect results
out_vec = snk.data()
llrs = snk_llr.data()
syms = snk_sym.data()
diff = array(in_vec) - array(out_vec)
err = [0 if i == 0 else 1 for i in diff]
self.dump(const.points(), in_vec, syms, llrs, out_vec)
print('Number of errors: %d' % (sum(err)))
# Check results
self.assertEqual(sum(err), 0)
def __init__(self, snr_db_ae = 15, signal_len = 1024, samp_rate = 100000, samples = 1000, const_index = 3):
gr.top_block.__init__(self, "Eve Sim")
##################################################
# Variables
##################################################
self.const_qpsk = const_qpsk = digital.constellation_qpsk().base()
self.const_bpsk = const_bpsk = digital.constellation_bpsk().base()
self.const_8psk = const_8psk = digital.constellation_8psk().base()
self.const_16qam = const_16qam = digital.constellation_16qam().base()
self.snr_db_ae = snr_db_ae # = 15
self.signal_len = signal_len # = 1024
self.samp_rate = samp_rate # = 100000
self.constellations = constellations = [const_bpsk, const_qpsk, const_8psk, const_16qam]
self.const_index = const_index
##################################################
# Blocks
##################################################
self.digital_chunks_to_symbols_xx_0 = digital.chunks_to_symbols_bc((constellations[const_index].points()), 1)
self.classify_trained_model_classifier_vc_0 = classify.trained_model_classifier_vc(64, '/home/gvanhoy/gr-classify/apps/cumulant_classifier.pkl')
self.channels_channel_model_0_0 = channels.channel_model(
noise_voltage=numpy.sqrt(10.0**(-snr_db_ae/10.0)/2),
frequency_offset=0.0,
epsilon=1.0,
taps=(1.0, ),
noise_seed=0,
block_tags=False
)
self.blocks_throttle_0_0 = blocks.throttle(gr.sizeof_gr_complex*1, samp_rate,True)
self.blocks_stream_to_vector_0 = blocks.stream_to_vector(gr.sizeof_gr_complex*1, 64)
self.blocks_repack_bits_bb_0 = blocks.repack_bits_bb(8, int(np.log2(constellations[const_index].arity())), "", False, gr.GR_MSB_FIRST)
self.blocks_head_1 = blocks.head(gr.sizeof_gr_complex*64, samples)
#message block to pull messages out of
self.blocks_message_debug_0 = blocks.message_debug()
self.analog_random_source_x_0 = blocks.vector_source_b(map(int, numpy.random.randint(0, 256, 10000)), True)
##################################################
# Connections
##################################################
#self.msg_connect((self.classify_trained_model_classifier_vc_0, 'classification_info'), (self.blocks_message_debug_0, 'print'))
self.msg_connect((self.classify_trained_model_classifier_vc_0, 'classification_info'), (self.blocks_message_debug_0, 'store'))
self.connect((self.analog_random_source_x_0, 0), (self.blocks_repack_bits_bb_0, 0))
self.connect((self.blocks_repack_bits_bb_0, 0), (self.digital_chunks_to_symbols_xx_0, 0))
#self.connect((self.blocks_stream_to_vector_0, 0), (self.classify_trained_model_classifier_vc_0, 0))
self.connect((self.blocks_stream_to_vector_0, 0), (self.blocks_head_1, 0))
self.connect((self.blocks_head_1, 0), (self.classify_trained_model_classifier_vc_0, 0))
self.connect((self.blocks_throttle_0_0, 0), (self.blocks_stream_to_vector_0, 0))
self.connect((self.channels_channel_model_0_0, 0), (self.blocks_throttle_0_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0, 0), (self.channels_channel_model_0_0, 0))
def test_001_identity(self):
# Constant modulus signal so no adjustments
const = digital.constellation_qpsk()
src_data = const.points()*1000
N = 100 # settling time
expected_data = src_data[N:]
result = self.transform(src_data, 0.1, const)[N:]
N = -500
self.assertComplexTuplesAlmostEqual(expected_data[N:], result[N:], 5)
def test_001_identity(self):
# Constant modulus signal so no adjustments
const = digital.constellation_qpsk()
src_data = const.points() * 1000
N = 100 # settling time
expected_data = src_data[N:]
result = self.transform(src_data, 0.1, const)[N:]
N = -500
self.assertComplexTuplesAlmostEqual(expected_data[N:], result[N:], 5)
def test_001_identity_lms(self):
# Constant modulus signal so no adjustments
const = digital.constellation_qpsk()
src_data = const.points() * 1000
alg = digital.adaptive_algorithm_lms(const, .1).base()
N = 100 # settling time
expected_data = src_data[N:]
result = self.transform(src_data, const, alg)[N:]
N = -500
self.assertComplexTuplesAlmostEqual(expected_data[N:], result[N:], 5)
def test_004_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)
demod.set_ic(ic)
self.assertEqual(ic, demod.get_ic())
def __init__(self):
ofdm_source.__init__(self, mod_name="ofdm_64_qpsk", fft_len=128)
self.const = digital.constellation_qpsk().base()
self.pack = blocks.packed_to_unpacked_bb(self.const.bits_per_symbol(),
gr.GR_MSB_FIRST)
self.map = digital.chunks_to_symbols_bc((self.const.points()), 1)
self.connect(self.random_source, self.pack, self.map)
self.connect(self.null, (self.mux, 0))
self.connect(self.map, (self.mux, 1))
self.connect(self.null, (self.mux, 2))
self.connect(self.mux, self.s2v, self.fft, self.cp, self.mult, self)
def _get_constellation(bps):
""" Returns a modulator block for a given number of bits per symbol """
constellation = {
1: digital.constellation_bpsk(),
2: digital.constellation_qpsk(),
3: digital.constellation_8psk()
}
try:
return constellation[bps]
except KeyError:
print 'Modulation not supported.'
exit(1)
def __init__(self, ipp1="127.0.0.1", ipp2="127.0.0.1", ipp3="127.0.0.1", ipp4="127.0.0.1", iptx="127.0.0.1", samp_rate=100):
gr.top_block.__init__(self, "OFDM Rx")
##################################################
# Parameters
##################################################
self.ipp1 = ipp1
self.ipp2 = ipp2
self.ipp3 = ipp3
self.ipp4 = ipp4
self.iptx = iptx
self.samp_rate = samp_rate
##################################################
# Variables
##################################################
self.pilot_symbols = pilot_symbols = ((1, 1, 1, -1,),)
self.pilot_carriers = pilot_carriers = ((-21, -7, 7, 21,),)
self.payload_mod = payload_mod = digital.constellation_qpsk()
self.packet_length_tag_key = packet_length_tag_key = "packet_len"
self.occupied_carriers = occupied_carriers = (range(-26, -21) + range(-20, -7) + range(-6, 0) + range(1, 7) + range(8, 21) + range(22, 27),)
self.length_tag_key = length_tag_key = "frame_len"
self.header_mod = header_mod = digital.constellation_bpsk()
self.fft_len = fft_len = 64
self.sync_word2 = sync_word2 = [0j, 0j, 0j, 0j, 0j, 0j, (-1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1 +0j), (1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), 0j, (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), 0j, 0j, 0j, 0j, 0j]
self.sync_word1 = sync_word1 = [0., 0., 0., 0., 0., 0., 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 0., 0., 0., 0., 0.]
self.payload_equalizer = payload_equalizer = digital.ofdm_equalizer_simpledfe(fft_len, payload_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols, 1)
self.packet_len = packet_len = 96
self.header_formatter = header_formatter = digital.packet_header_ofdm(occupied_carriers, n_syms=1, len_tag_key=packet_length_tag_key, frame_len_tag_key=length_tag_key, bits_per_header_sym=header_mod.bits_per_symbol(), bits_per_payload_sym=payload_mod.bits_per_symbol(), scramble_header=False)
self.header_equalizer = header_equalizer = digital.ofdm_equalizer_simpledfe(fft_len, header_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols)
##################################################
# Blocks
##################################################
self.zeromq_pull_source_0_0_0_0_0_0 = zeromq.pull_source(gr.sizeof_gr_complex, 64, "tcp://" + ipp2 + ":55521", 100, True)
self.my_number_sync_timestamp_0 = my.number_sync_timestamp()
self.fft_vxx_1 = fft.fft_vcc(fft_len, True, (), True, 1)
self.digital_ofdm_serializer_vcc_payload = digital.ofdm_serializer_vcc(fft_len, occupied_carriers, length_tag_key, packet_length_tag_key, 1, "", True)
self.digital_ofdm_frame_equalizer_vcvc_1 = digital.ofdm_frame_equalizer_vcvc(payload_equalizer.base(), fft_len/4, length_tag_key, True, 0)
self.digital_crc32_bb_0 = digital.crc32_bb(True, packet_length_tag_key, True)
self.digital_constellation_decoder_cb_1 = digital.constellation_decoder_cb(payload_mod.base())
self.blocks_repack_bits_bb_0 = blocks.repack_bits_bb(payload_mod.bits_per_symbol(), 8, packet_length_tag_key, True, gr.GR_LSB_FIRST)
##################################################
# Connections
##################################################
self.connect((self.blocks_repack_bits_bb_0, 0), (self.digital_crc32_bb_0, 0))
self.connect((self.digital_constellation_decoder_cb_1, 0), (self.blocks_repack_bits_bb_0, 0))
self.connect((self.digital_crc32_bb_0, 0), (self.my_number_sync_timestamp_0, 0))
self.connect((self.digital_ofdm_frame_equalizer_vcvc_1, 0), (self.digital_ofdm_serializer_vcc_payload, 0))
self.connect((self.digital_ofdm_serializer_vcc_payload, 0), (self.digital_constellation_decoder_cb_1, 0))
self.connect((self.fft_vxx_1, 0), (self.digital_ofdm_frame_equalizer_vcvc_1, 0))
self.connect((self.zeromq_pull_source_0_0_0_0_0_0, 0), (self.fft_vxx_1, 0))
def test_qpsk_3tap_lms_training(self):
# set up fg
gain = 0.01 # LMS gain
num_taps = 16
num_samp = 2000
num_test = 500
cons = digital.constellation_qpsk().base()
rxmod = digital.generic_mod(cons, False, self.sps, True, self.eb,
False, False)
modulated_sync_word_pre = digital.modulate_vector_bc(
rxmod.to_basic_block(), self.preamble + self.preamble, [1])
modulated_sync_word = modulated_sync_word_pre[86:(
512 + 86)] # compensate for the RRC filter delay
corr_max = numpy.abs(
numpy.dot(modulated_sync_word, numpy.conj(modulated_sync_word)))
corr_calc = self.corr_thresh / (corr_max * corr_max)
preamble_symbols = self.map_symbols_to_constellation(
self.unpack_values(self.preamble, 8, 2), cons)
alg = digital.adaptive_algorithm_lms(cons, gain).base()
evm = digital.meas_evm_cc(cons, digital.evm_measurement_t.EVM_PERCENT)
leq = digital.linear_equalizer(num_taps, self.sps, alg, False,
preamble_symbols, 'corr_est')
correst = digital.corr_est_cc(modulated_sync_word, self.sps, 12,
corr_calc, digital.THRESHOLD_ABSOLUTE)
constmod = digital.generic_mod(constellation=cons,
differential=False,
samples_per_symbol=4,
pre_diff_code=True,
excess_bw=0.35,
verbose=False,
log=False)
chan = channels.channel_model(noise_voltage=0.0,
frequency_offset=0.0,
epsilon=1.0,
taps=(1.0 + 1.0j, 0.63 - .22j,
-.1 + .07j),
noise_seed=0,
block_tags=False)
vso = blocks.vector_source_b(self.preamble + self.data, True, 1, [])
head = blocks.head(gr.sizeof_float * 1, num_samp)
vsi = blocks.vector_sink_f()
self.tb.connect(vso, constmod, chan, correst, leq, evm, head, vsi)
self.tb.run()
# look at the last 1000 samples, should converge quickly, below 5% EVM
upper_bound = list(20.0 * numpy.ones((num_test, )))
lower_bound = list(0.0 * numpy.zeros((num_test, )))
output_data = vsi.data()
output_data = output_data[-num_test:]
self.assertLess(output_data, upper_bound)
self.assertGreater(output_data, lower_bound)
def _get_constellation(bps):
""" Returns a modulator block for a given number of bits per symbol """
constellation = {
1: digital.constellation_bpsk(),
2: digital.constellation_qpsk(),
3: digital.constellation_8psk()
}
try:
return constellation[bps]
except KeyError:
print 'Modulation not supported.'
exit(1)
def __init__(self, ipp1="127.0.0.1", ipp2="127.0.0.1", ipp3="127.0.0.1", ipp4="127.0.0.1", iptx="127.0.0.1", samp_rate=10000):
gr.top_block.__init__(self, "OFDM Rx")
##################################################
# Parameters
##################################################
self.ipp1 = ipp1
self.ipp2 = ipp2
self.ipp3 = ipp3
self.ipp4 = ipp4
self.iptx = iptx
self.samp_rate = samp_rate
##################################################
# Variables
##################################################
self.pilot_symbols = pilot_symbols = ((1, 1, 1, -1,),)
self.pilot_carriers = pilot_carriers = ((-21, -7, 7, 21,),)
self.payload_mod = payload_mod = digital.constellation_qpsk()
self.packet_length_tag_key = packet_length_tag_key = "packet_len"
self.occupied_carriers = occupied_carriers = (range(-26, -21) + range(-20, -7) + range(-6, 0) + range(1, 7) + range(8, 21) + range(22, 27),)
self.length_tag_key = length_tag_key = "frame_len"
self.header_mod = header_mod = digital.constellation_bpsk()
self.fft_len = fft_len = 64
self.sync_word2 = sync_word2 = [0j, 0j, 0j, 0j, 0j, 0j, (-1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1 +0j), (1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), 0j, (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), 0j, 0j, 0j, 0j, 0j]
self.sync_word1 = sync_word1 = [0., 0., 0., 0., 0., 0., 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 0., 0., 0., 0., 0.]
self.packet_len = packet_len = 96
self.header_formatter = header_formatter = digital.packet_header_ofdm(occupied_carriers, n_syms=1, len_tag_key=packet_length_tag_key, frame_len_tag_key=length_tag_key, bits_per_header_sym=header_mod.bits_per_symbol(), bits_per_payload_sym=payload_mod.bits_per_symbol(), scramble_header=False)
self.header_equalizer = header_equalizer = digital.ofdm_equalizer_simpledfe(fft_len, header_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols)
##################################################
# Blocks
##################################################
self.zeromq_push_sink_0_0_0 = zeromq.push_sink(gr.sizeof_char, 1, "tcp://"+ ipp1 + ":55511", 100, True)
self.zeromq_push_sink_0_0 = zeromq.push_sink(gr.sizeof_gr_complex, 1, "tcp://"+ ipp1 + ":55510", 100, True)
self.zeromq_pull_source_0 = zeromq.pull_source(gr.sizeof_gr_complex, 1, "tcp://"+ iptx + ":55500", 100, True)
self.digital_ofdm_sync_sc_cfb_0 = digital.ofdm_sync_sc_cfb(fft_len, fft_len/4, False)
self.blocks_multiply_xx_1 = blocks.multiply_vcc(1)
self.blocks_delay_0 = blocks.delay(gr.sizeof_gr_complex*1, fft_len+fft_len/4)
self.analog_frequency_modulator_fc_0 = analog.frequency_modulator_fc(-2.0/fft_len)
##################################################
# Connections
##################################################
self.connect((self.analog_frequency_modulator_fc_0, 0), (self.blocks_multiply_xx_1, 0))
self.connect((self.blocks_delay_0, 0), (self.blocks_multiply_xx_1, 1))
self.connect((self.blocks_multiply_xx_1, 0), (self.zeromq_push_sink_0_0, 0))
self.connect((self.digital_ofdm_sync_sc_cfb_0, 0), (self.analog_frequency_modulator_fc_0, 0))
self.connect((self.digital_ofdm_sync_sc_cfb_0, 1), (self.zeromq_push_sink_0_0_0, 0))
self.connect((self.zeromq_pull_source_0, 0), (self.blocks_delay_0, 0))
self.connect((self.zeromq_pull_source_0, 0), (self.digital_ofdm_sync_sc_cfb_0, 0))
def test_002_simpledfe (self):
""" Use the simple DFE equalizer. """
fft_len = 8
# 4 5 6 7 0 1 2 3
tx_data = [-1, -1, 1, 2, -1, 3, 0, -1, # 0
-1, -1, 0, 2, -1, 2, 0, -1, # 8
-1, -1, 3, 0, -1, 1, 0, -1, # 16 (Pilot symbols)
-1, -1, 1, 1, -1, 0, 2, -1] # 24
cnst = digital.constellation_qpsk()
tx_signal = [cnst.map_to_points_v(x)[0] if x != -1 else 0 for x in tx_data]
occupied_carriers = ((1, 2, 6, 7),)
pilot_carriers = ((), (), (1, 2, 6, 7), ())
pilot_symbols = (
[], [], [cnst.map_to_points_v(x)[0] for x in (1, 0, 3, 0)], []
)
equalizer = digital.ofdm_equalizer_simpledfe(
fft_len, cnst.base(), occupied_carriers, pilot_carriers, pilot_symbols, 0, 0.01
)
channel = [
0, 0, 1, 1, 0, 1, 1, 0,
0, 0, 1, 1, 0, 1, 1, 0, # These coefficients will be rotated slightly...
0, 0, 1j, 1j, 0, 1j, 1j, 0, # Go crazy here!
0, 0, 1j, 1j, 0, 1j, 1j, 0 # ...and again here.
]
for idx in range(fft_len, 2*fft_len):
channel[idx] = channel[idx-fft_len] * numpy.exp(1j * .1 * numpy.pi * (numpy.random.rand()-.5))
idx2 = idx+2*fft_len
channel[idx2] = channel[idx2] * numpy.exp(1j * 0 * numpy.pi * (numpy.random.rand()-.5))
len_tag_key = "frame_len"
len_tag = gr.tag_t()
len_tag.offset = 0
len_tag.key = pmt.string_to_symbol(len_tag_key)
len_tag.value = pmt.from_long(4)
chan_tag = gr.tag_t()
chan_tag.offset = 0
chan_tag.key = pmt.string_to_symbol("ofdm_sync_chan_taps")
chan_tag.value = pmt.init_c32vector(fft_len, channel[:fft_len])
src = blocks.vector_source_c(numpy.multiply(tx_signal, channel), False, fft_len, (len_tag, chan_tag))
eq = digital.ofdm_frame_equalizer_vcvc(equalizer.base(), 0, len_tag_key, True)
sink = blocks.vector_sink_c(fft_len)
self.tb.connect(src, eq, sink)
self.tb.run ()
rx_data = [cnst.decision_maker_v((x,)) if x != 0 else -1 for x in sink.data()]
self.assertEqual(tx_data, rx_data)
for tag in sink.tags():
if pmt.symbol_to_string(tag.key) == len_tag_key:
self.assertEqual(pmt.to_long(tag.value), 4)
if pmt.symbol_to_string(tag.key) == "ofdm_sync_chan_taps":
self.assertComplexTuplesAlmostEqual(list(pmt.c32vector_elements(tag.value)), channel[-fft_len:], places=1)
def test_004_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 get_constellation_from_string(self, const_string):
self.const = {
'ook': constellations.constellation_ook(),
'bpsk': digital.constellation_bpsk().base(),
'qpsk': digital.constellation_qpsk().base(),
'4ask': constellations.constellation_4_ask(),
'8psk': digital.constellation_8psk().base(),
'8qam_cross': constellations.constellation_8qam_cross(),
'8qam_circular': constellations.constellation_8qam_circular(),
'16qam': digital.constellation_16qam().base(),
'16psk': constellations.constellation_16_psk(),
'32qam_cross': constellations.constellation_32qam_cross(),
'64qam': constellations.constellation_64qam(),
}.get(const_string,
digital.constellation_bpsk().base())
def test_002_simpledfe (self):
""" Use the simple DFE equalizer. """
fft_len = 8
# 4 5 6 7 0 1 2 3
tx_data = [-1, -1, 1, 2, -1, 3, 0, -1, # 0
-1, -1, 0, 2, -1, 2, 0, -1, # 8
-1, -1, 3, 0, -1, 1, 0, -1, # 16 (Pilot symbols)
-1, -1, 1, 1, -1, 0, 2, -1] # 24
cnst = digital.constellation_qpsk()
tx_signal = [cnst.map_to_points_v(x)[0] if x != -1 else 0 for x in tx_data]
occupied_carriers = ((1, 2, 6, 7),)
pilot_carriers = ((), (), (1, 2, 6, 7), ())
pilot_symbols = (
[], [], [cnst.map_to_points_v(x)[0] for x in (1, 0, 3, 0)], []
)
equalizer = digital.ofdm_equalizer_simpledfe(
fft_len, cnst.base(), occupied_carriers, pilot_carriers, pilot_symbols, 0, 0.01
)
channel = [
0, 0, 1, 1, 0, 1, 1, 0,
0, 0, 1, 1, 0, 1, 1, 0, # These coefficients will be rotated slightly...
0, 0, 1j, 1j, 0, 1j, 1j, 0, # Go crazy here!
0, 0, 1j, 1j, 0, 1j, 1j, 0 # ...and again here.
]
for idx in range(fft_len, 2*fft_len):
channel[idx] = channel[idx-fft_len] * numpy.exp(1j * .1 * numpy.pi * (numpy.random.rand()-.5))
idx2 = idx+2*fft_len
channel[idx2] = channel[idx2] * numpy.exp(1j * 0 * numpy.pi * (numpy.random.rand()-.5))
chan_tag = gr.tag_t()
chan_tag.offset = 0
chan_tag.key = pmt.string_to_symbol("ofdm_sync_chan_taps")
chan_tag.value = pmt.init_c32vector(fft_len, channel[:fft_len])
src = blocks.vector_source_c(numpy.multiply(tx_signal, channel), False, fft_len, (chan_tag,))
eq = digital.ofdm_frame_equalizer_vcvc(equalizer.base(), 0, self.tsb_key, True)
sink = blocks.tsb_vector_sink_c(fft_len, tsb_key=self.tsb_key)
self.tb.connect(
src,
blocks.stream_to_tagged_stream(gr.sizeof_gr_complex, fft_len, len(tx_data) // fft_len, self.tsb_key),
eq,
sink
)
self.tb.run ()
rx_data = [cnst.decision_maker_v((x,)) if x != 0 else -1 for x in sink.data()[0]]
self.assertEqual(tx_data, rx_data)
self.assertEqual(len(sink.tags()), 1)
tag = sink.tags()[0]
self.assertEqual(pmt.symbol_to_string(tag.key), "ofdm_sync_chan_taps")
self.assertComplexTuplesAlmostEqual(list(pmt.c32vector_elements(tag.value)), channel[-fft_len:], places=1)
def test_constellation_encoder_bc_qpsk_random(self):
cnst = digital.constellation_qpsk()
src_data = np.random.randint(0, 4, size=20000)
expected_result = [cnst.points()[x] for x in src_data]
src = blocks.vector_source_b(src_data)
op = digital.constellation_encoder_bc(cnst.base())
dst = blocks.vector_sink_c()
self.tb.connect(src, op)
self.tb.connect(op, dst)
self.tb.run() # run the graph and wait for it to finish
actual_result = dst.data() # fetch the contents of the sink
# print "actual result", actual_result
# print "expected result", expected_result
self.assertFloatTuplesAlmostEqual(expected_result, actual_result)
def test_002_static_wo_tags (self):
""" Same as before, but the input stream has no tag.
We specify the frame size in the constructor.
We also specify a tag key, so the output stream *should* have
a length tag.
"""
fft_len = 8
n_syms = 4
# 4 5 6 7 0 1 2 3
tx_data = [-1, -1, 1, 2, -1, 3, 0, -1, # 0
-1, -1, 0, 2, -1, 2, 0, -1, # 8
-1, -1, 3, 0, -1, 1, 0, -1, # 16 (Pilot symbols)
-1, -1, 1, 1, -1, 0, 2, -1] # 24
cnst = digital.constellation_qpsk()
tx_signal = [cnst.map_to_points_v(x)[0] if x != -1 else 0 for x in tx_data]
occupied_carriers = ((1, 2, 6, 7),)
pilot_carriers = ((), (), (1, 2, 6, 7), ())
pilot_symbols = (
[], [], [cnst.map_to_points_v(x)[0] for x in (1, 0, 3, 0)], []
)
equalizer = digital.ofdm_equalizer_static(fft_len, occupied_carriers, pilot_carriers, pilot_symbols)
channel = [
0, 0, 1, 1, 0, 1, 1, 0,
0, 0, 1, 1, 0, 1, 1, 0, # These coefficients will be rotated slightly (below)...
0, 0, 1j, 1j, 0, 1j, 1j, 0, # Go crazy here!
0, 0, 1j, 1j, 0, 1j, 1j, 0 # ...and again here.
]
for idx in range(fft_len, 2*fft_len):
channel[idx] = channel[idx-fft_len] * numpy.exp(1j * .1 * numpy.pi * (numpy.random.rand()-.5))
idx2 = idx+2*fft_len
channel[idx2] = channel[idx2] * numpy.exp(1j * 0 * numpy.pi * (numpy.random.rand()-.5))
src = gr.vector_source_c(numpy.multiply(tx_signal, channel), False, fft_len)
# We do specify a length tag, it should then appear at the output
eq = digital.ofdm_frame_equalizer_vcvc(equalizer.base(), 0, "frame_len", False, n_syms)
sink = blocks.vector_sink_c(fft_len)
self.tb.connect(src, eq, sink)
self.tb.run ()
rx_data = [cnst.decision_maker_v((x,)) if x != 0 else -1 for x in sink.data()]
self.assertEqual(tx_data, rx_data)
# Check len tag
tags = sink.tags()
len_tag = dict()
for tag in tags:
ptag = gr.tag_to_python(tag)
if ptag.key == 'frame_len':
len_tag[ptag.key] = ptag.value
self.assertEqual(len_tag, {'frame_len': 4})
def test_002_static_wo_tags (self):
""" Same as before, but the input stream has no tag.
We specify the frame size in the constructor.
We also specify a tag key, so the output stream *should* have
a TSB tag.
"""
fft_len = 8
n_syms = 4
# 4 5 6 7 0 1 2 3
tx_data = [-1, -1, 1, 2, -1, 3, 0, -1, # 0
-1, -1, 0, 2, -1, 2, 0, -1, # 8
-1, -1, 3, 0, -1, 1, 0, -1, # 16 (Pilot symbols)
-1, -1, 1, 1, -1, 0, 2, -1] # 24
cnst = digital.constellation_qpsk()
tx_signal = [cnst.map_to_points_v(x)[0] if x != -1 else 0 for x in tx_data]
occupied_carriers = ((1, 2, 6, 7),)
pilot_carriers = ((), (), (1, 2, 6, 7), ())
pilot_symbols = (
[], [], [cnst.map_to_points_v(x)[0] for x in (1, 0, 3, 0)], []
)
equalizer = digital.ofdm_equalizer_static(fft_len, occupied_carriers, pilot_carriers, pilot_symbols)
channel = [
0, 0, 1, 1, 0, 1, 1, 0,
0, 0, 1, 1, 0, 1, 1, 0, # These coefficients will be rotated slightly (below)...
0, 0, 1j, 1j, 0, 1j, 1j, 0, # Go crazy here!
0, 0, 1j, 1j, 0, 1j, 1j, 0 # ...and again here.
]
for idx in range(fft_len, 2*fft_len):
channel[idx] = channel[idx-fft_len] * numpy.exp(1j * .1 * numpy.pi * (numpy.random.rand()-.5))
idx2 = idx+2*fft_len
channel[idx2] = channel[idx2] * numpy.exp(1j * 0 * numpy.pi * (numpy.random.rand()-.5))
src = blocks.vector_source_c(numpy.multiply(tx_signal, channel), False, fft_len)
eq = digital.ofdm_frame_equalizer_vcvc(equalizer.base(), 0, self.tsb_key, False, n_syms)
sink = blocks.tsb_vector_sink_c(vlen=fft_len, tsb_key=self.tsb_key)
self.tb.connect(
src,
blocks.stream_to_tagged_stream(gr.sizeof_gr_complex, fft_len, len(tx_data) // fft_len, self.tsb_key),
eq,
sink
)
self.tb.run ()
rx_data = [cnst.decision_maker_v((x,)) if x != 0 else -1 for x in sink.data()[0]]
self.assertEqual(tx_data, rx_data)
# Check TSB Functionality
packets = sink.data()
self.assertEqual(len(packets), 1)
self.assertEqual(len(packets[0]), len(tx_data))
def _test_constellation_decoder_cb_qpsk(self):
cnst = digital.constellation_qpsk()
src_data = (0.5 + 0.5j, 0.1 - 1.2j, -0.8 - 0.1j, -0.45 + 0.8j,
0.8 + 1.0j, -0.5 + 0.1j, 0.1 - 1.2j)
expected_result = (3, 1, 0, 2, 3, 2, 1)
src = blocks.vector_source_c(src_data)
op = digital.constellation_decoder_cb(cnst.base())
dst = blocks.vector_sink_b()
self.tb.connect(src, op)
self.tb.connect(op, dst)
self.tb.run() # run the graph and wait for it to finish
actual_result = dst.data() # fetch the contents of the sink
#print "actual result", actual_result
#print "expected result", expected_result
self.assertFloatTuplesAlmostEqual(expected_result, actual_result)
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 __init__(self, constellation_points=_def_constellation_points,
*args, **kwargs):
"""
Hierarchical block for RRC-filtered QPSK modulation.
The input is a byte stream (unsigned char) and the
output is the complex modulated signal at baseband.
See generic_demod block for list of parameters.
"""
constellation_points = _def_constellation_points
constellation = digital_swig.constellation_qpsk()
if constellation_points != 4:
raise ValueError('Number of constellation points must be 4 for QPSK.')
super(qpsk_demod, self).__init__(constellation=constellation,
*args, **kwargs)
def _test_constellation_decoder_cb_qpsk(self):
cnst = digital.constellation_qpsk()
src_data = (0.5 + 0.5j, 0.1 - 1.2j, -0.8 - 0.1j, -0.45 + 0.8j,
0.8 + 1.0j, -0.5 + 0.1j, 0.1 - 1.2j)
expected_result = ( 3, 1, 0, 2,
3, 2, 1)
src = blocks.vector_source_c(src_data)
op = digital.constellation_decoder_cb(cnst.base())
dst = blocks.vector_sink_b()
self.tb.connect(src, op)
self.tb.connect(op, dst)
self.tb.run() # run the graph and wait for it to finish
actual_result = dst.data() # fetch the contents of the sink
#print "actual result", actual_result
#print "expected result", expected_result
self.assertFloatTuplesAlmostEqual(expected_result, actual_result)
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)
def test_001_t (self):
# Trellislength
K = 3072 # Bit in one packet
n = 2 # Bit in codeword
k = 5 # Constraint length
start_state = 0
end_state = -1 # Endstate not defined
fsm = fsm_args["awgn1o2_16"]
## setup dummy data
os = numpy.array(fsm[4], dtype=int) # Encoder output matrix
data = numpy.random.randint(0,2,K) # Create 0s and 1s
sym_table = digital.constellation_qpsk()
# Setup blocks
data_src = blocks.vector_source_s(map(int, data))
src_head = blocks.head(gr.sizeof_short*1, K)
# TX Sim
encoder = trellis.encoder_ss(trellis.fsm(*fsm), 0)
modulator = digital.chunks_to_symbols_sc(sym_table.points(), 1)
# Decoder
viterbi_cel = celec.gen_viterbi_fi(n, k, K, start_state, end_state,
sym_table.points(), os)
# Sinks
rx_sink = blocks.vector_sink_b(1)
# Connections
self.tb.connect(data_src, src_head, encoder)
self.tb.connect(encoder, modulator)
self.tb.connect(modulator, viterbi_cel)
self.tb.connect(viterbi_cel, rx_sink)
self.tb.run ()
# # Check data
rx_output = numpy.array(rx_sink.data())
for k in range(0, K):
self.assertEqual(int(rx_output[k]), int(data[k]))
def test_001_t(self):
# Trellislength
K = 3072 # Bit in one packet
n = 2 # Bit in codeword
k = 5 # Constraint length
start_state = 0
end_state = -1 # Endstate not defined
fsm = fsm_args["awgn1o2_16"]
## setup dummy data
os = numpy.array(fsm[4], dtype=int) # Encoder output matrix
data = numpy.random.randint(0, 2, K) # Create 0s and 1s
sym_table = digital.constellation_qpsk()
# Setup blocks
data_src = blocks.vector_source_s(map(int, data))
src_head = blocks.head(gr.sizeof_short * 1, K)
# TX Sim
encoder = trellis.encoder_ss(trellis.fsm(*fsm), 0)
modulator = digital.chunks_to_symbols_sc(sym_table.points(), 1)
# Decoder
viterbi_cel = celec.gen_viterbi_fi(n, k, K, start_state, end_state,
sym_table.points(), os)
# Sinks
rx_sink = blocks.vector_sink_b(1)
# Connections
self.tb.connect(data_src, src_head, encoder)
self.tb.connect(encoder, modulator)
self.tb.connect(modulator, viterbi_cel)
self.tb.connect(viterbi_cel, rx_sink)
self.tb.run()
# # Check data
rx_output = numpy.array(rx_sink.data())
for k in range(0, K):
self.assertEqual(int(rx_output[k]), int(data[k]))
def __init__(self,
constellation_points=_def_constellation_points,
*args,
**kwargs):
"""
Hierarchical block for RRC-filtered QPSK modulation.
The input is a byte stream (unsigned char) and the
output is the complex modulated signal at baseband.
See generic_demod block for list of parameters.
"""
constellation_points = _def_constellation_points
constellation = digital_swig.constellation_qpsk()
if constellation_points != 4:
raise ValueError(
'Number of constellation points must be 4 for QPSK.')
super(qpsk_demod, self).__init__(constellation=constellation,
*args,
**kwargs)
def test_002_t(self):
"""Decoding of noisy QPSK symbols"""
# Parameters
M = 4
N0 = 1.0
max_len = 50
snr_db = 20
# Constants
rndm = random.Random()
noise_v = 1 / math.sqrt((10**(float(snr_db) / 10)))
in_vec = tuple([rndm.randint(0, 1) for i in range(0, max_len)])
# Flowgraph
src = blocks.vector_source_b(in_vec)
pack = blocks.repack_bits_bb(1, 2, "", False, gr.GR_MSB_FIRST)
const = digital.constellation_qpsk().base()
cmap = digital.chunks_to_symbols_bc(const.points())
nadder = blocks.add_cc()
noise = analog.noise_source_c(analog.GR_GAUSSIAN, noise_v, 0)
soft_decoder = blocksat.soft_decoder_cf(M, N0)
slicer = digital.binary_slicer_fb()
inverter = digital.map_bb([1, 0])
snk = blocks.vector_sink_b()
snk_llr = blocks.vector_sink_f()
self.tb.connect(src, pack, cmap)
self.tb.connect(cmap, (nadder, 0))
self.tb.connect(noise, (nadder, 1))
self.tb.connect(nadder, soft_decoder, slicer, inverter, snk)
self.tb.connect(soft_decoder, snk_llr)
self.tb.run()
# check data
dec_vec = snk.data()
llrs = snk_llr.data()
print(llrs)
print(in_vec)
print(dec_vec)
self.assertFloatTuplesAlmostEqual(in_vec, dec_vec, places=4)
def __init__(self, constellation_points=_def_constellation_points,
gray_coded=_def_gray_coded,
*args, **kwargs):
"""
Hierarchical block for RRC-filtered QPSK modulation.
The input is a byte stream (unsigned char) and the
output is the complex modulated signal at baseband.
See generic_mod block for list of parameters.
"""
constellation_points = _def_constellation_points
constellation = digital_swig.constellation_qpsk()
if constellation_points != 4:
raise ValueError("QPSK can only have 4 constellation points.")
if not gray_coded:
raise ValueError("This QPSK mod/demod works only for gray-coded constellations.")
super(qpsk_mod, self).__init__(constellation=constellation,
gray_coded=gray_coded,
*args, **kwargs)
def test_003_active_subcarriers(self):
n_frames = 1
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)
data = get_random_qpsk(n_frames * timeslots * active_subcarriers)
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, 64, f_taps, gfdm_constellation, subcarrier_map)
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())
self.assertComplexTuplesAlmostEqual(data, res, 2)
def test_qpsk_nonzeroevm(self):
# set up fg
expected_result = tuple(numpy.zeros((self.num_data, )))
self.cons = cons = digital.constellation_qpsk().base()
self.data = data = [
random.randrange(len(cons.points())) for x in range(self.num_data)
]
self.symbols = symbols = numpy.squeeze(
[cons.map_to_points_v(i) for i in data])
evm = digital.meas_evm_cc(cons, digital.evm_measurement_t_EVM_PERCENT)
vso = blocks.vector_source_c(symbols, False, 1, [])
mc = blocks.multiply_const_cc(3.0 + 2.0j)
vsi = blocks.vector_sink_f()
self.tb.connect(vso, mc, evm, vsi)
self.tb.run()
# check data
output_data = vsi.data()
self.assertNotEqual(expected_result, output_data)
# FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL # THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER # LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, # OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN # THE SOFTWARE. #!/usr/bin/env python from gnuradio import gr, digital from math import pi import dvb_swig M = 4 # QPSK is of order 4 nfilts = 32 mod_constellation = digital.constellation_qpsk() # Receiver parameters (fmin, fmax) = (-0.5, 0.5) # Allow Costas frequency swing of +/- half of the sample rate freq_alpha = 0.01 # FLL alpha gain freq_bw = 0.035 # Will be overriden by set_alpha(), set_beta() timing_alpha = 0.10 # Costas loop gain for frequency adjustments timing_beta = 0.01 # Costas loop gain for phase adjustments phase_alpha = 0.01 phase_bw = 0.035 # Will be overriden by set_alpha(), set_beta() class s_modulator_bc(gr.hier_block2): """ DVB-S modulator Input: unsigned char (unpacked bits)
def __init__(self):
gr.top_block.__init__(self, "Audio Tx")
Qt.QWidget.__init__(self)
self.setWindowTitle("Audio Tx")
qtgui.util.check_set_qss()
try:
self.setWindowIcon(Qt.QIcon.fromTheme('gnuradio-grc'))
except:
pass
self.top_scroll_layout = Qt.QVBoxLayout()
self.setLayout(self.top_scroll_layout)
self.top_scroll = Qt.QScrollArea()
self.top_scroll.setFrameStyle(Qt.QFrame.NoFrame)
self.top_scroll_layout.addWidget(self.top_scroll)
self.top_scroll.setWidgetResizable(True)
self.top_widget = Qt.QWidget()
self.top_scroll.setWidget(self.top_widget)
self.top_layout = Qt.QVBoxLayout(self.top_widget)
self.top_grid_layout = Qt.QGridLayout()
self.top_layout.addLayout(self.top_grid_layout)
self.settings = Qt.QSettings("GNU Radio", "audio_tx")
if StrictVersion(Qt.qVersion()) < StrictVersion("5.0.0"):
self.restoreGeometry(self.settings.value("geometry").toByteArray())
else:
self.restoreGeometry(
self.settings.value("geometry", type=QtCore.QByteArray))
##################################################
# Variables
##################################################
self.sps = sps = 1
self.nfilts = nfilts = 25
self.samp_rate = samp_rate = 32E3
self.rrc_taps = rrc_taps = firdes.root_raised_cosine(
nfilts, nfilts, 1.0 / float(sps), 0.35, 45 * nfilts)
self.qpsk = qpsk = digital.constellation_rect(([
0.707 + 0.707j, -0.707 + 0.707j, -0.707 - 0.707j, 0.707 - 0.707j
]), ([0, 1, 2, 3]), 4, 2, 2, 1, 1).base()
self.phase_bw = phase_bw = 6.28 / 100.0
self.excess_bw = excess_bw = 0.35
self.BPSK = BPSK = digital.constellation_qpsk().base()
##################################################
# Blocks
##################################################
self._phase_bw_range = Range(0.0, 1.0, 0.01, 6.28 / 100.0, 200)
self._phase_bw_win = RangeWidget(self._phase_bw_range,
self.set_phase_bw, 'Phase: Bandwidth',
"slider", float)
self.top_layout.addWidget(self._phase_bw_win)
self.qtgui_time_sink_x_1_1 = qtgui.time_sink_f(
1024, #size
samp_rate, #samp_rate
"", #name
1 #number of inputs
)
self.qtgui_time_sink_x_1_1.set_update_time(0.10)
self.qtgui_time_sink_x_1_1.set_y_axis(-1, 1)
self.qtgui_time_sink_x_1_1.set_y_label('Amplitude', "")
self.qtgui_time_sink_x_1_1.enable_tags(-1, True)
self.qtgui_time_sink_x_1_1.set_trigger_mode(qtgui.TRIG_MODE_FREE,
qtgui.TRIG_SLOPE_POS, 0.0,
0, 0, "")
self.qtgui_time_sink_x_1_1.enable_autoscale(False)
self.qtgui_time_sink_x_1_1.enable_grid(False)
self.qtgui_time_sink_x_1_1.enable_axis_labels(True)
self.qtgui_time_sink_x_1_1.enable_control_panel(False)
self.qtgui_time_sink_x_1_1.enable_stem_plot(False)
if not True:
self.qtgui_time_sink_x_1_1.disable_legend()
labels = ['', '', '', '', '', '', '', '', '', '']
widths = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
colors = [
"blue", "red", "green", "black", "cyan", "magenta", "yellow",
"dark red", "dark green", "blue"
]
styles = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1, -1, -1, -1, -1, -1]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_time_sink_x_1_1.set_line_label(
i, "Data {0}".format(i))
else:
self.qtgui_time_sink_x_1_1.set_line_label(i, labels[i])
self.qtgui_time_sink_x_1_1.set_line_width(i, widths[i])
self.qtgui_time_sink_x_1_1.set_line_color(i, colors[i])
self.qtgui_time_sink_x_1_1.set_line_style(i, styles[i])
self.qtgui_time_sink_x_1_1.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_1_1.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_1_1_win = sip.wrapinstance(
self.qtgui_time_sink_x_1_1.pyqwidget(), Qt.QWidget)
self.top_layout.addWidget(self._qtgui_time_sink_x_1_1_win)
self.qtgui_time_sink_x_1_0 = qtgui.time_sink_f(
1024, #size
samp_rate, #samp_rate
"", #name
1 #number of inputs
)
self.qtgui_time_sink_x_1_0.set_update_time(0.10)
self.qtgui_time_sink_x_1_0.set_y_axis(-1, 1)
self.qtgui_time_sink_x_1_0.set_y_label('Amplitude', "")
self.qtgui_time_sink_x_1_0.enable_tags(-1, True)
self.qtgui_time_sink_x_1_0.set_trigger_mode(qtgui.TRIG_MODE_FREE,
qtgui.TRIG_SLOPE_POS, 0.0,
0, 0, "")
self.qtgui_time_sink_x_1_0.enable_autoscale(False)
self.qtgui_time_sink_x_1_0.enable_grid(False)
self.qtgui_time_sink_x_1_0.enable_axis_labels(True)
self.qtgui_time_sink_x_1_0.enable_control_panel(False)
self.qtgui_time_sink_x_1_0.enable_stem_plot(False)
if not True:
self.qtgui_time_sink_x_1_0.disable_legend()
labels = ['', '', '', '', '', '', '', '', '', '']
widths = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
colors = [
"blue", "red", "green", "black", "cyan", "magenta", "yellow",
"dark red", "dark green", "blue"
]
styles = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1, -1, -1, -1, -1, -1]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_time_sink_x_1_0.set_line_label(
i, "Data {0}".format(i))
else:
self.qtgui_time_sink_x_1_0.set_line_label(i, labels[i])
self.qtgui_time_sink_x_1_0.set_line_width(i, widths[i])
self.qtgui_time_sink_x_1_0.set_line_color(i, colors[i])
self.qtgui_time_sink_x_1_0.set_line_style(i, styles[i])
self.qtgui_time_sink_x_1_0.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_1_0.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_1_0_win = sip.wrapinstance(
self.qtgui_time_sink_x_1_0.pyqwidget(), Qt.QWidget)
self.top_layout.addWidget(self._qtgui_time_sink_x_1_0_win)
self.qtgui_time_sink_x_1 = qtgui.time_sink_f(
1024, #size
samp_rate, #samp_rate
"", #name
1 #number of inputs
)
self.qtgui_time_sink_x_1.set_update_time(0.10)
self.qtgui_time_sink_x_1.set_y_axis(-1, 1)
self.qtgui_time_sink_x_1.set_y_label('Amplitude', "")
self.qtgui_time_sink_x_1.enable_tags(-1, True)
self.qtgui_time_sink_x_1.set_trigger_mode(qtgui.TRIG_MODE_FREE,
qtgui.TRIG_SLOPE_POS, 0.0, 0,
0, "")
self.qtgui_time_sink_x_1.enable_autoscale(False)
self.qtgui_time_sink_x_1.enable_grid(False)
self.qtgui_time_sink_x_1.enable_axis_labels(True)
self.qtgui_time_sink_x_1.enable_control_panel(False)
self.qtgui_time_sink_x_1.enable_stem_plot(False)
if not True:
self.qtgui_time_sink_x_1.disable_legend()
labels = ['', '', '', '', '', '', '', '', '', '']
widths = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
colors = [
"blue", "red", "green", "black", "cyan", "magenta", "yellow",
"dark red", "dark green", "blue"
]
styles = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1, -1, -1, -1, -1, -1]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_time_sink_x_1.set_line_label(
i, "Data {0}".format(i))
else:
self.qtgui_time_sink_x_1.set_line_label(i, labels[i])
self.qtgui_time_sink_x_1.set_line_width(i, widths[i])
self.qtgui_time_sink_x_1.set_line_color(i, colors[i])
self.qtgui_time_sink_x_1.set_line_style(i, styles[i])
self.qtgui_time_sink_x_1.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_1.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_1_win = sip.wrapinstance(
self.qtgui_time_sink_x_1.pyqwidget(), Qt.QWidget)
self.top_layout.addWidget(self._qtgui_time_sink_x_1_win)
self.qtgui_time_sink_x_0 = qtgui.time_sink_f(
1024, #size
samp_rate, #samp_rate
"", #name
1 #number of inputs
)
self.qtgui_time_sink_x_0.set_update_time(0.10)
self.qtgui_time_sink_x_0.set_y_axis(-1, 1)
self.qtgui_time_sink_x_0.set_y_label('Amplitude', "")
self.qtgui_time_sink_x_0.enable_tags(-1, True)
self.qtgui_time_sink_x_0.set_trigger_mode(qtgui.TRIG_MODE_FREE,
qtgui.TRIG_SLOPE_POS, 0.0, 0,
0, "")
self.qtgui_time_sink_x_0.enable_autoscale(False)
self.qtgui_time_sink_x_0.enable_grid(False)
self.qtgui_time_sink_x_0.enable_axis_labels(True)
self.qtgui_time_sink_x_0.enable_control_panel(False)
self.qtgui_time_sink_x_0.enable_stem_plot(False)
if not True:
self.qtgui_time_sink_x_0.disable_legend()
labels = ['', '', '', '', '', '', '', '', '', '']
widths = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
colors = [
"blue", "red", "green", "black", "cyan", "magenta", "yellow",
"dark red", "dark green", "blue"
]
styles = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1, -1, -1, -1, -1, -1]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_time_sink_x_0.set_line_label(
i, "Data {0}".format(i))
else:
self.qtgui_time_sink_x_0.set_line_label(i, labels[i])
self.qtgui_time_sink_x_0.set_line_width(i, widths[i])
self.qtgui_time_sink_x_0.set_line_color(i, colors[i])
self.qtgui_time_sink_x_0.set_line_style(i, styles[i])
self.qtgui_time_sink_x_0.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_0.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_0_win = sip.wrapinstance(
self.qtgui_time_sink_x_0.pyqwidget(), Qt.QWidget)
self.top_layout.addWidget(self._qtgui_time_sink_x_0_win)
self.qtgui_const_sink_x_0_0 = qtgui.const_sink_c(
1024, #size
"demod", #name
1 #number of inputs
)
self.qtgui_const_sink_x_0_0.set_update_time(0.10)
self.qtgui_const_sink_x_0_0.set_y_axis(-2, 2)
self.qtgui_const_sink_x_0_0.set_x_axis(-2, 2)
self.qtgui_const_sink_x_0_0.set_trigger_mode(qtgui.TRIG_MODE_FREE,
qtgui.TRIG_SLOPE_POS, 0.0,
0, "")
self.qtgui_const_sink_x_0_0.enable_autoscale(False)
self.qtgui_const_sink_x_0_0.enable_grid(True)
self.qtgui_const_sink_x_0_0.enable_axis_labels(True)
if not True:
self.qtgui_const_sink_x_0_0.disable_legend()
labels = ['', '', '', '', '', '', '', '', '', '']
widths = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
colors = [
"blue", "red", "red", "red", "red", "red", "red", "red", "red",
"red"
]
styles = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
markers = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_const_sink_x_0_0.set_line_label(
i, "Data {0}".format(i))
else:
self.qtgui_const_sink_x_0_0.set_line_label(i, labels[i])
self.qtgui_const_sink_x_0_0.set_line_width(i, widths[i])
self.qtgui_const_sink_x_0_0.set_line_color(i, colors[i])
self.qtgui_const_sink_x_0_0.set_line_style(i, styles[i])
self.qtgui_const_sink_x_0_0.set_line_marker(i, markers[i])
self.qtgui_const_sink_x_0_0.set_line_alpha(i, alphas[i])
self._qtgui_const_sink_x_0_0_win = sip.wrapinstance(
self.qtgui_const_sink_x_0_0.pyqwidget(), Qt.QWidget)
self.top_layout.addWidget(self._qtgui_const_sink_x_0_0_win)
self.qtgui_const_sink_x_0 = qtgui.const_sink_c(
1024, #size
"Rx", #name
1 #number of inputs
)
self.qtgui_const_sink_x_0.set_update_time(0.10)
self.qtgui_const_sink_x_0.set_y_axis(-2, 2)
self.qtgui_const_sink_x_0.set_x_axis(-2, 2)
self.qtgui_const_sink_x_0.set_trigger_mode(qtgui.TRIG_MODE_FREE,
qtgui.TRIG_SLOPE_POS, 0.0,
0, "")
self.qtgui_const_sink_x_0.enable_autoscale(False)
self.qtgui_const_sink_x_0.enable_grid(True)
self.qtgui_const_sink_x_0.enable_axis_labels(True)
if not True:
self.qtgui_const_sink_x_0.disable_legend()
labels = ['', '', '', '', '', '', '', '', '', '']
widths = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
colors = [
"blue", "red", "red", "red", "red", "red", "red", "red", "red",
"red"
]
styles = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
markers = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_const_sink_x_0.set_line_label(
i, "Data {0}".format(i))
else:
self.qtgui_const_sink_x_0.set_line_label(i, labels[i])
self.qtgui_const_sink_x_0.set_line_width(i, widths[i])
self.qtgui_const_sink_x_0.set_line_color(i, colors[i])
self.qtgui_const_sink_x_0.set_line_style(i, styles[i])
self.qtgui_const_sink_x_0.set_line_marker(i, markers[i])
self.qtgui_const_sink_x_0.set_line_alpha(i, alphas[i])
self._qtgui_const_sink_x_0_win = sip.wrapinstance(
self.qtgui_const_sink_x_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_const_sink_x_0_win, 0, 1, 1,
1)
[self.top_grid_layout.setRowStretch(r, 1) for r in range(0, 1)]
[self.top_grid_layout.setColumnStretch(c, 1) for c in range(1, 2)]
self.low_pass_filter_0 = filter.fir_filter_ccf(
1,
firdes.low_pass(1, samp_rate, 1.6E3, .6E3, firdes.WIN_HAMMING,
6.76))
self.hilbert_fc_0 = filter.hilbert_fc(65, firdes.WIN_HAMMING, 6.76)
self.digital_pfb_clock_sync_xxx_0 = digital.pfb_clock_sync_ccf(
sps, .063, (rrc_taps), nfilts, nfilts / 2, 1.5, 1)
self.digital_lms_dd_equalizer_cc_0 = digital.lms_dd_equalizer_cc(
8, .01, 1, BPSK)
self.digital_diff_decoder_bb_0 = digital.diff_decoder_bb(4)
self.digital_costas_loop_cc_0 = digital.costas_loop_cc(
phase_bw, 4, False)
self.digital_constellation_receiver_cb_0 = digital.constellation_receiver_cb(
BPSK, phase_bw, 0, 2 * 3.14)
self.digital_constellation_modulator_0 = digital.generic_mod(
constellation=BPSK,
differential=True,
samples_per_symbol=sps,
pre_diff_code=True,
excess_bw=2,
verbose=False,
log=False,
)
self.blocks_multiply_xx_3_0 = blocks.multiply_vff(1)
self.blocks_multiply_xx_3 = blocks.multiply_vff(1)
self.blocks_multiply_xx_2_1 = blocks.multiply_vff(1)
self.blocks_multiply_xx_2_0_0 = blocks.multiply_vcc(1)
self.blocks_multiply_xx_2_0 = blocks.multiply_vcc(1)
self.blocks_multiply_xx_2 = blocks.multiply_vff(1)
self.blocks_multiply_xx_1 = blocks.multiply_vff(1)
self.blocks_multiply_xx_0 = blocks.multiply_vff(1)
self.blocks_float_to_complex_0 = blocks.float_to_complex(1)
self.blocks_file_source_0 = blocks.file_source(
gr.sizeof_char * 1,
'/home/peter/Desktop/acoustic_radio/test_input.txt', True)
self.blocks_file_sink_0 = blocks.file_sink(
gr.sizeof_char * 1, '/home/peter/Desktop/output.txt', False)
self.blocks_file_sink_0.set_unbuffered(False)
self.blocks_complex_to_imag_1 = blocks.complex_to_imag(1)
self.blocks_complex_to_imag_0 = blocks.complex_to_imag(1)
self.blocks_complex_to_float_0 = blocks.complex_to_float(1)
self.blocks_add_xx_3 = blocks.add_vff(1)
self.blocks_add_xx_1 = blocks.add_vff(1)
self.blocks_add_xx_0 = blocks.add_vff(1)
self.blks2_packet_encoder_0 = grc_blks2.packet_mod_b(
grc_blks2.packet_encoder(
samples_per_symbol=sps,
bits_per_symbol=1,
preamble='',
access_code='',
pad_for_usrp=False,
),
payload_length=1,
)
self.blks2_packet_decoder_0 = grc_blks2.packet_demod_b(
grc_blks2.packet_decoder(
access_code='',
threshold=-1,
callback=lambda ok, payload: self.blks2_packet_decoder_0.
recv_pkt(ok, payload),
), )
self.analog_sig_source_x_0_1_0_0 = analog.sig_source_c(
samp_rate, analog.GR_SIN_WAVE, 4E3, 1, 0)
self.analog_sig_source_x_0_1_0 = analog.sig_source_c(
samp_rate, analog.GR_COS_WAVE, 4E3, 1, 0)
self.analog_sig_source_x_0_1 = analog.sig_source_f(
samp_rate, analog.GR_COS_WAVE, 4E3, 1, 0)
self.analog_sig_source_x_0_0_0 = analog.sig_source_f(
samp_rate, analog.GR_SIN_WAVE, 4E3, 1, 0)
self.analog_sig_source_x_0_0 = analog.sig_source_f(
samp_rate, analog.GR_SIN_WAVE, 4E3, 1, 0)
self.analog_sig_source_x_0 = analog.sig_source_f(
samp_rate, analog.GR_COS_WAVE, 4E3, 1, 0)
self.analog_feedforward_agc_cc_0 = analog.feedforward_agc_cc(
1024, 1.55)
self.analog_const_source_x_0_0 = analog.sig_source_f(
0, analog.GR_CONST_WAVE, 0, 0, -1)
self.analog_const_source_x_0 = analog.sig_source_f(
0, analog.GR_CONST_WAVE, 0, 0, -1)
##################################################
# Connections
##################################################
self.connect((self.analog_const_source_x_0, 0),
(self.blocks_multiply_xx_3, 1))
self.connect((self.analog_const_source_x_0_0, 0),
(self.blocks_multiply_xx_3_0, 1))
self.connect((self.analog_feedforward_agc_cc_0, 0),
(self.digital_pfb_clock_sync_xxx_0, 0))
self.connect((self.analog_sig_source_x_0, 0),
(self.blocks_multiply_xx_0, 0))
self.connect((self.analog_sig_source_x_0_0, 0),
(self.blocks_multiply_xx_3_0, 0))
self.connect((self.analog_sig_source_x_0_0_0, 0),
(self.blocks_multiply_xx_3, 0))
self.connect((self.analog_sig_source_x_0_1, 0),
(self.blocks_multiply_xx_2, 1))
self.connect((self.analog_sig_source_x_0_1_0, 0),
(self.blocks_multiply_xx_2_0, 1))
self.connect((self.analog_sig_source_x_0_1_0_0, 0),
(self.blocks_multiply_xx_2_0_0, 1))
self.connect((self.blks2_packet_decoder_0, 0),
(self.blocks_file_sink_0, 0))
self.connect((self.blks2_packet_encoder_0, 0),
(self.digital_constellation_modulator_0, 0))
self.connect((self.blocks_add_xx_0, 0), (self.blocks_multiply_xx_2, 0))
self.connect((self.blocks_add_xx_0, 0),
(self.blocks_multiply_xx_2_1, 0))
self.connect((self.blocks_add_xx_0, 0), (self.hilbert_fc_0, 0))
self.connect((self.blocks_add_xx_0, 0), (self.qtgui_time_sink_x_0, 0))
self.connect((self.blocks_add_xx_1, 0),
(self.blocks_float_to_complex_0, 1))
self.connect((self.blocks_add_xx_3, 0),
(self.blocks_float_to_complex_0, 0))
self.connect((self.blocks_complex_to_float_0, 0),
(self.blocks_multiply_xx_0, 1))
self.connect((self.blocks_complex_to_float_0, 1),
(self.blocks_multiply_xx_1, 0))
self.connect((self.blocks_complex_to_imag_0, 0),
(self.blocks_add_xx_1, 0))
self.connect((self.blocks_complex_to_imag_1, 0),
(self.blocks_add_xx_3, 1))
self.connect((self.blocks_file_source_0, 0),
(self.blks2_packet_encoder_0, 0))
self.connect((self.blocks_float_to_complex_0, 0),
(self.low_pass_filter_0, 0))
self.connect((self.blocks_multiply_xx_0, 0), (self.blocks_add_xx_0, 0))
self.connect((self.blocks_multiply_xx_1, 0), (self.blocks_add_xx_0, 1))
self.connect((self.blocks_multiply_xx_2, 0), (self.blocks_add_xx_1, 1))
self.connect((self.blocks_multiply_xx_2_0, 0),
(self.blocks_complex_to_imag_0, 0))
self.connect((self.blocks_multiply_xx_2_0_0, 0),
(self.blocks_complex_to_imag_1, 0))
self.connect((self.blocks_multiply_xx_2_1, 0),
(self.blocks_add_xx_3, 0))
self.connect((self.blocks_multiply_xx_3, 0),
(self.blocks_multiply_xx_2_1, 1))
self.connect((self.blocks_multiply_xx_3_0, 0),
(self.blocks_multiply_xx_1, 1))
self.connect((self.digital_constellation_modulator_0, 0),
(self.blocks_complex_to_float_0, 0))
self.connect((self.digital_constellation_receiver_cb_0, 0),
(self.digital_diff_decoder_bb_0, 0))
self.connect((self.digital_constellation_receiver_cb_0, 4),
(self.qtgui_const_sink_x_0_0, 0))
self.connect((self.digital_constellation_receiver_cb_0, 3),
(self.qtgui_time_sink_x_1, 0))
self.connect((self.digital_constellation_receiver_cb_0, 1),
(self.qtgui_time_sink_x_1_0, 0))
self.connect((self.digital_constellation_receiver_cb_0, 2),
(self.qtgui_time_sink_x_1_1, 0))
self.connect((self.digital_costas_loop_cc_0, 0),
(self.digital_lms_dd_equalizer_cc_0, 0))
self.connect((self.digital_diff_decoder_bb_0, 0),
(self.blks2_packet_decoder_0, 0))
self.connect((self.digital_lms_dd_equalizer_cc_0, 0),
(self.digital_constellation_receiver_cb_0, 0))
self.connect((self.digital_lms_dd_equalizer_cc_0, 0),
(self.qtgui_const_sink_x_0, 0))
self.connect((self.digital_pfb_clock_sync_xxx_0, 0),
(self.digital_costas_loop_cc_0, 0))
self.connect((self.hilbert_fc_0, 0), (self.blocks_multiply_xx_2_0, 0))
self.connect((self.hilbert_fc_0, 0),
(self.blocks_multiply_xx_2_0_0, 0))
self.connect((self.low_pass_filter_0, 0),
(self.analog_feedforward_agc_cc_0, 0))
bits_per_header = (symbols_per_header + 1) * header_mod.bits_per_symbol()
symbols_per_header = bits_per_header / header_mod.bits_per_symbol()
#header_formatter = cdma.packet_header(bits_per_header,length_tag_name,num_tag_name,header_mod.bits_per_symbol());
#header_formatter = digital.packet_header_default(bits_per_header, length_tag_name,num_tag_name,header_mod.bits_per_symbol());
#tcm_indicator_symbols_per_frame=4; #Zhe added, 4 bits are used as tcm mode indicator, it is used as a part of header.
# Achilles' comment: this may change later when filler bits are introduced...
print "bits_per_header=", bits_per_header
print "symbols_per_header=", symbols_per_header
#print "tcm_indicator_symbols_per_frame=",tcm_indicator_symbols_per_frame
print "\n"
#trellis coding and modulation info
payload_mod = [digital.constellation_qpsk(), digital.constellation_8psk()]
pdir = prefix + "/python/fsm_files/"
fsm = [pdir + "awgn2o2_1.fsm", pdir + "awgn2o3_8ungerboecka.fsm"]
uncoded_fsm = [
trellis.fsm(2, 2, [1, 0, 0, 1]),
trellis.fsm(3, 3, [1, 0, 0, 0, 1, 0, 0, 0, 1])
]
bits_per_coded_symbol = [
int(math.log(trellis.fsm(fsm[i]).O(), 2)) for i in range(len(payload_mod))
]
#coding_rate=[Fraction(int(math.log(trellis.fsm(fsm[i]).I(),2)), int(math.log(trellis.fsm(fsm[i]).O(),2))) for i in range(len(fsm))]
if bits_per_coded_symbol != [
def __init__(self):
gr.top_block.__init__(self, "Transmitter NC-OFDM")
filepath = os.getcwd()
global expduration
##################################################
# Read parameter for configuration
##################################################
file_handle = open('/root/gr-ncofdm/examples/config_ncofdm_tx', 'r')
for line in file_handle: # read rest of lines
linefromfile = ([x for x in line.split()])
linesize = len(linefromfile)-1
if linesize > 1:
linedata = map(float, linefromfile[1:linesize+1])
else:
linedata = float(linefromfile[1])
if linefromfile[0] == "fft_len":
fft_len = int(linedata)
elif linefromfile[0] == "cp_len":
cp_len = int(linedata)
elif linefromfile[0] == "shseq_len":
shseq_len = int(linedata)
elif linefromfile[0] == "lgseq_len":
lgseq_len = int(linedata)
elif linefromfile[0] == "dataseq_len":
dataseq_len = int(linedata)
elif linefromfile[0] == "shseq_rep":
shseq_rep = int(linedata)
elif linefromfile[0] == "samp_rate":
samp_rate = int(linedata)
elif linefromfile[0] == "pnseq_offset":
pnseq_offset = int(linedata)
elif linefromfile[0] == "shseq":
shseq = [int(z) for z in linedata]
elif linefromfile[0] == "lgseq":
lgseq = [int(z) for z in linedata]
elif linefromfile[0] == "dataseq":
dataseq = [int(z) for z in linedata]
elif linefromfile[0] == "occupied_carriers":
oc = []
for z in linedata:
oc.append(int(z))
occupied_carriers = ((oc),)#[int(z) for z in linedata]
elif linefromfile[0] == "pilot_carriers":
pc = []
for z in linedata:
pc.append(int(z))
pilot_carriers = ((pc),)#[int(z) for z in linedata]
elif linefromfile[0] == "pilot_symbols":
ps = []
for z in linedata:
ps.append(int(z))
pilot_symbols = ((ps),)#[int(z) for z in linedata]
elif linefromfile[0] == "cen_freq":
cen_freq = linedata
elif linefromfile[0] == "sig_ul_ratio_db":
sig_ul_ratio_db = linedata
elif linefromfile[0] == "usrp_gain":
usrp_gain = linedata
elif linefromfile[0] == "final_gain":
final_gain = linedata
elif linefromfile[0] == "expduration":
expduration = linedata
##################################################
# Variables
##################################################
self.pn_symbols = pn_symbols = (1, -1, -1, 1, 1, -1, -1, 1, 1, 1, -1, 1, -1, 1, 1, 1, -1, 1, 1, 1, 1, -1, -1, -1, 1, -1, -1, -1, -1, -1, 1, 1, -1, 1, -1, 1, 1, 1, 1, -1, -1, -1, 1, 1, -1, 1, 1, -1, -1, -1, 1, -1, -1, 1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, 1, -1, 1, -1, -1, -1, 1, -1, -1, -1, -1, 1, -1, 1, 1, 1, 1, -1, -1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, 1, 1, 1, -1, 1, -1, -1, 1, -1, 1, -1, 1, -1, -1, -1, -1, -1, 1, -1, -1, 1, 1, 1, -1, 1, -1, 1, 1, 1, 1, 1, 1, -1, 1, 1, 1, 1, -1, 1, -1, -1, 1, 1, -1, -1, -1, -1, -1, 1, -1, 1, 1, -1, -1, 1, 1, -1, 1, 1, 1, -1, 1, -1, 1, 1, -1, -1, 1, -1, 1, -1, -1, -1, 1, -1, 1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, 1, -1, -1, -1, 1, 1, 1, -1, 1, 1, -1, 1, 1, -1, -1, 1, 1, -1, -1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, 1, -1, 1, -1, 1, 1, -1, -1, 1, -1, -1, -1, 1, -1, -1, -1, 1, -1, -1, -1, -1, -1, -1, 1, -1, -1, -1, 1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, 1, -1, 1, -1, -1, 1, 1, -1, 1, -1, -1, -1, 1, -1, 1, 1, -1, 1, -1, -1, -1, 1, -1, -1, -1, 1, 1, -1, -1, 1, -1, 1, -1, 1, -1, -1, 1, -1, -1, 1, 1, 1, -1, -1, 1, 1, 1, 1, 1, -1, 1, -1, -1, 1, 1, 1, 1, 1, 1, -1, -1, -1, -1, 1, -1, -1, 1, -1, -1, -1, -1, 1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, 1, -1, 1, 1, -1, -1, 1, -1, -1, 1, -1, -1, -1, 1, -1, 1, -1, -1, 1, -1, -1, -1, -1, 1, -1, 1, -1, 1, -1, -1, -1, -1, 1, 1, 1, -1, -1, 1, 1, -1, -1, -1, -1, -1, 1, 1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, -1, 1, 1, -1, -1, -1, 1, -1, -1, -1, -1, -1, -1, 1, -1, 1, 1, -1, 1, 1, -1, 1, 1, -1, -1, 1, 1, 1, 1, 1, -1, 1, -1, -1, -1, -1, 1, -1, -1, 1, -1, -1, -1, 1, -1, -1, -1, 1, -1, 1, -1, -1, 1, 1, 1, -1, -1, 1, -1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, 1, 1, -1, -1, 1, -1, -1, -1, 1, -1, 1, -1, -1, -1, 1, -1, -1, -1, -1, -1, -1, 1, -1, -1, -1, 1, -1, -1, 1, -1, -1, 1, -1, 1, -1, -1, 1, 1, -1, 1, -1, -1, -1, -1, -1, 1, -1, -1, 1, 1, -1, -1, 1, -1, 1, -1, 1, 1, 1, 1, 1, 1, -1, -1, -1, 1, -1, 1, -1, -1, 1, -1, -1, 1, -1, -1, -1, -1, -1, 1, 1, 1, -1, 1, 1, -1, -1, 1, -1, 1, -1, -1, 1, 1, 1, -1, 1, -1, -1, 1, -1, 1, 1, 1, 1, 1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, 1, 1, 1, 1, -1, 1, -1, 1, 1, 1, 1, 1, -1, -1, 1, 1, -1, -1, -1, -1, 1, 1, 1, -1, 1, 1, 1, 1, 1, 1, -1, -1, -1, 1, -1, -1, -1, -1, 1, -1, 1, -1, -1, -1, 1, 1, 1, -1, -1, 1, -1, 1, -1, 1, 1, 1, 1, 1, -1, -1, -1, -1, 1, 1, -1, 1, -1, -1, -1, 1, -1, 1, -1, 1, 1, -1, -1, 1, -1, -1, 1, -1, 1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, 1, 1, 1, 1, -1, 1, 1, 1, 1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1, 1, -1, -1, 1, 1, 1, 1, 1, -1, -1, -1, 1, -1, 1, 1, 1, -1, 1, -1, -1, -1, -1, 1, 1, 1, -1, -1, -1, 1, -1, 1, 1, 1, 1, -1, 1, -1, -1, 1, 1, -1, 1, -1, 1, 1, 1, 1, 1, 1, 1, 1, -1, 1, 1, -1, 1, -1, -1, -1, -1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 1, -1, 1, 1, -1, 1, 1, -1, -1, 1, 1, 1, 1, 1, 1, 1, 1, -1, 1, -1, 1, -1, 1, 1, -1, -1, -1, 1, 1, 1, 1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, -1, -1, 1, -1, -1, -1, -1, 1, -1, 1, -1, 1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, -1, 1, 1, -1, 1, 1, 1, 1, -1, -1, 1, 1, 1, 1, -1, -1, 1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1, 1, -1, 1, 1, -1, -1, -1, -1, 1, 1, -1, -1, 1, -1, -1, -1, -1, 1, -1, -1, 1, -1, 1, -1, 1, -1, -1, -1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, -1, 1, -1, 1, 1, -1, 1, 1, 1, 1, 1, -1, 1, 1, -1, -1, -1, -1, -1, -1, -1, -1, 1, -1, -1, -1, -1, -1, -1, 1, 1, 1, -1, 1, -1, 1, 1, 1, 1, 1, -1, 1, -1, 1, 1, 1, 1, -1, -1, -1, -1, -1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, 1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, -1, 1, -1, -1, -1, -1, 1, -1, 1, 1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1, 1, 1, -1, 1, 1, -1, -1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, 1, -1, -1, -1, 1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, 1, 1, -1, -1, -1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, -1, -1, -1, 1, -1, 1, -1, -1, 1, -1, 1, 1, -1, 1, 1, 1, 1, 1, -1, 1, 1, -1, -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, 1, -1, 1, -1, 1, 1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, -1, 1, -1, -1, -1, 1, -1, 1, -1, -1, 1, -1, 1, -1, -1, 1, 1, -1, 1, -1, 1, -1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, -1, 1, 1, -1, 1, 1, -1, 1, 1, 1, 1, 1, 1, -1, -1, 1, -1, 1, 1, 1, -1, -1, 1, 1, 1, -1, -1, -1, 1, 1, 1, -1, 1, -1, -1, -1, 1, -1, 1, 1, -1, 1, -1, 1, 1, -1, -1, 1, -1, -1, 1, 1, 1, -1, -1, -1, -1, 1, -1, -1, 1, 1, 1, 1, 1, -1, 1, 1, 1, 1, 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, -1, 1, 1, 1, -1, 1, 1, -1, -1, 1, -1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, -1, -1, 1, -1, -1, -1, 1, -1, -1, 1, 1, 1, -1, -1, 1, -1, 1, 1, -1, -1, -1, 1, 1, -1, -1, -1, 1, 1, -1, -1, 1, -1, 1, -1, 1, 1, 1, -1, -1, 1, -1, -1, -1, -1, 1, -1, -1, 1, -1, -1, -1, -1, -1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, -1, 1, 1, -1, 1, -1, 1, -1, 1, -1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, 1, 1, -1, -1, 1, 1, 1, -1, -1, 1, 1, -1, -1, -1, -1, -1, 1, 1, -1, -1, 1, 1, -1, 1, 1, 1, 1, -1, 1, -1, -1, 1, -1, -1, -1, 1, -1, -1, -1, 1, -1, -1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, -1, 1, 1, -1, 1, -1, -1, -1, -1, -1, 1, 1, -1, 1, -1, -1, 1, -1, -1, -1, 1, -1, -1, 1, -1, 1, 1, 1, -1, -1, 1, -1, -1, 1, 1, 1, -1, 1, 1, -1, -1, -1, -1, 1, -1, -1, 1, -1, -1, -1, 1, 1, -1, 1, -1, -1, -1, 1, -1, -1, -1, 1, -1, 1, -1, 1, 1, 1, 1, 1, -1, -1, -1, -1, 1, -1, 1, 1, 1, -1, -1, 1, 1, -1, 1, 1, -1, -1, 1, -1, 1, -1, 1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, -1, -1, -1, 1, -1, 1, 1, 1, 1, -1, 1, -1, -1, -1, 1, -1, -1, -1, 1, 1, 1, -1, 1, 1, -1, -1, 1, -1, 1, 1, 1, -1, -1, 1, 1, 1, 1, 1, -1, 1, -1, -1, -1, 1, 1, -1, -1, -1, -1, 1, 1, -1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, -1, -1, -1, -1, -1, -1, 1, -1, 1, -1, 1, -1, 1, 1, 1, -1, 1, -1, -1, -1, -1, -1, -1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1, 1, -1, 1, 1, 1, 1, 1, 1, -1, 1, 1, -1, -1, 1, -1, -1, -1, -1, -1, -1, -1, -1, 1, -1, 1, 1, 1, 1, 1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, -1, 1, -1, 1, 1, -1, 1, -1, -1, -1, -1, 1, -1, -1, 1, -1, 1, 1, -1, 1, 1, 1, 1, 1, -1, 1, -1, 1, -1, 1, 1, 1, -1, -1, -1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1, 1, -1, 1, -1, -1, 1, -1, 1, 1, -1, -1, -1, 1, 1, 1, -1, -1, 1, 1, 1, -1, -1, -1, 1, -1, 1, 1, 1, 1, -1, 1, -1, -1, -1, -1, -1, -1, -1, 1, -1, -1, -1, 1, 1, -1, 1, 1, 1, 1, 1, 1, -1, 1, -1, -1, -1, -1, -1, -1, 1, 1, 1, -1, -1, -1, -1, 1, -1, 1, -1, 1, -1, -1, -1, -1, -1, -1, -1, 1, 1, -1, -1, 1, 1, 1, -1, 1, 1, 1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, 1, 1, -1, 1, 1, -1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, -1, 1, -1, 1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, -1, 1, -1, -1, -1, -1, 1, 1, 1, -1, -1, -1, -1, -1, -1, -1, -1, 1, -1, -1, -1, 1, -1, -1, -1, -1, -1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, 1, -1, 1, -1, 1, 1, -1, 1, 1, -1, 1, 1, 1, -1, -1, -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, -1, 1, -1, -1, 1, 1, 1, -1, 1, 1, -1, 1, 1, 1, 1, -1, -1, 1, 1, 1, -1, -1, -1, -1, -1, -1, 1, -1, -1, 1, 1, 1, 1, 1, -1, 1, 1, -1, -1, 1, -1, -1, 1, 1, 1, 1, 1, 1, 1, 1, 1, -1, -1, 1, -1, 1, 1, 1, 1, 1, 1, 1, -1, 1, -1, 1, 1, 1, 1, -1, -1, -1, -1, 1, 1, -1, -1, -1, 1, -1, -1, -1, 1, 1, -1, -1, -1, 1, -1, 1, -1, 1, -1, 1, 1, 1, -1, -1, -1, 1, 1, 1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, 1, -1, -1, 1, -1, 1, -1, -1, -1, -1, -1, -1, 1, -1, -1, 1, 1, 1, -1, 1, 1, 1, 1, -1, 1, -1, 1, 1, 1, 1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, -1, 1, -1, -1, -1, -1, 1, 1, 1, 1, 1, 1, -1, -1, 1, -1, 1, -1, -1, -1, -1, -1, -1, 1, -1, 1, 1, 1, 1, -1, 1, 1, -1, -1, 1, 1, 1, 1, 1, 1, -1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, -1, 1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1, -1, 1, 1, -1, 1, -1, 1, -1, 1, -1, -1, 1, -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, 1, -1, -1, -1, 1, 1, 1, 1, 1, -1, -1, -1, -1, -1, 1, -1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, 1, 1, 1, 1, -1, 1, 1, -1, 1, -1, 1, -1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, -1, -1, 1, 1, -1, 1, 1, -1, -1, 1, 1, 1, -1, 1, 1, -1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, -1, -1, -1, 1, -1, -1, -1, 1, -1, 1, -1, 1, -1, -1, 1, -1, -1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, 1, 1, -1, -1, 1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, -1, 1, 1, 1, 1, 1, 1, 1, -1, -1, -1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, -1, 1, -1, 1, -1, -1, 1, 1, -1, 1, -1, 1, -1, -1, -1, -1, 1, -1, -1, -1, 1, -1, -1, 1, 1, 1, -1, -1, 1, 1, -1, -1, -1, 1, -1, 1, -1, 1, 1, -1, -1, -1, -1, 1, 1, -1, 1, 1, -1, 1, 1, 1, 1, 1, -1, -1, -1, 1, -1, 1, 1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, -1, -1, -1, -1, 1, -1, -1, -1, 1, -1, 1, 1, 1, -1, -1, 1, 1, 1, -1, -1, 1, -1, -1, -1, -1, -1, -1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, 1, 1, -1, 1, 1, -1, 1, 1, -1, -1, 1, 1, 1, -1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, 1, -1, 1, 1, 1, -1, -1, -1, 1, -1, 1, -1, -1, -1, 1, -1, 1, -1, -1, 1, -1, -1, 1, -1, -1, 1, 1, 1, 1, 1, 1, 1, -1, 1, -1, -1, -1, -1, -1, 1, -1, 1, -1, -1, 1, 1, 1, -1, 1, -1, -1, -1, -1, 1, 1, -1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, -1, 1, 1, 1, 1, 1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, 1, -1, 1, 1, 1, -1, -1, -1, 1, -1, 1, 1, -1, 1, -1, 1, -1, -1, 1, -1, -1, -1, 1, -1, 1, -1, -1, -1, -1, -1, 1, 1, -1, 1, 1, -1, 1, 1, 1, 1, 1, 1, -1, 1, -1, -1, 1, 1, 1, 1, 1, -1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, 1, -1, -1, 1, 1, 1, -1, 1, -1, 1, 1, -1, 1, -1, -1, -1, -1, 1, 1, -1, 1, 1, -1, -1, 1, -1, 1, -1, 1, 1, 1, -1, -1, 1, 1, 1, -1, -1, 1, 1, -1, -1, -1, 1, 1, -1, 1, -1, 1, 1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1, 1, 1, -1, 1, -1, 1, 1, -1, 1, 1, -1, 1, 1, -1, 1, -1, -1, -1, 1, 1, 1, 1, 1, -1, 1, -1, 1, 1, -1, -1, -1, 1, -1, 1, 1, -1, -1, 1, -1, -1, -1, -1, 1, -1, 1, 1, -1, 1, -1, 1, 1, 1, 1, -1, -1, 1, -1, 1, -1, 1, -1, -1, -1, 1, -1, -1, -1, -1, -1, 1, -1, 1, 1, -1, -1, -1, 1, 1, 1, -1, 1, -1, -1, -1, -1, 1, -1, -1, -1, 1, 1, 1, 1, 1, 1, -1, -1, -1, 1, -1, -1, 1, -1, 1, 1, 1, -1, 1, 1, -1, -1, -1, 1, -1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, -1, 1, -1, 1, 1, 1, -1, 1, 1, -1, -1, 1, 1, 1, 1, -1, 1, 1, 1, -1, 1, 1, -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, 1, 1, 1, 1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, -1, -1, 1, 1, 1, -1, -1, 1, 1, 1, 1, -1, 1, 1, 1, 1, 1, -1, 1, 1, 1, 1, -1, 1, 1, 1, 1, 1, -1, 1, 1, -1, 1, 1, 1, 1, 1, -1, -1, 1, -1, 1, 1, -1, -1, -1, 1, -1, 1, -1, 1, -1, 1, 1, 1, 1, 1, -1, -1, 1, 1, -1, 1, -1, 1, 1, 1, 1, -1, -1, 1, 1, 1, -1, -1, 1, -1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, 1, 1, 1, 1, 1, 1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, 1, 1, -1, 1, 1, -1, 1, 1, 1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, 1, 1, -1, -1, 1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1, -1, -1, 1, -1, -1, -1, -1, 1, -1, -1, -1, -1, 1, 1, -1, -1, -1, 1, -1, -1, 1, 1, 1, -1, 1, 1, 1, 1, -1, 1, 1, -1, 1, 1, 1, 1, 1, -1, 1, 1, -1, 1, -1, -1, -1, 1, -1, 1, 1, -1, -1, 1, -1, 1, 1, 1, 1, -1, -1, 1, 1, 1, -1, -1, -1, -1, -1, 1, -1, 1, -1, 1, -1, 1, 1, -1, 1, -1, -1, -1, -1, -1, 1, -1, -1, -1, 1, 1, -1, 1, 1, 1, 1, 1, 1, -1, -1, 1, 1, 1, -1, -1, 1, 1, -1, 1, 1, -1, 1, -1, 1, 1, 1, 1, 1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1, 1, -1, -1, -1, -1, 1, -1, 1, 1, 1, 1, -1, -1, -1, -1, 1, -1, -1, -1, -1, 1, 1, -1, 1, 1, -1, 1, -1, 1, 1, 1, 1, -1, -1, -1, 1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, 1, -1, -1, 1, -1, -1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, 1, 1, 1, 1, -1, -1, 1, -1, -1, -1, -1, -1, 1, -1, 1, -1, 1, -1, 1, 1, -1, 1, -1, 1, -1, 1, 1, 1, 1, 1, 1, 1, -1, -1, -1, 1, -1, -1, -1, -1, -1, 1, -1, 1, 1, -1, -1, -1, 1, -1, -1, -1, -1, -1, 1, -1, 1, -1, -1, -1, 1, 1, 1, 1, 1, -1, -1, -1, 1, -1, -1, 1, 1, -1, -1, -1, 1, 1, 1, -1, 1, 1, 1, -1, -1, -1, 1, -1, 1, 1, 1, 1, 1, 1, -1, 1, 1, -1, 1, 1, 1, 1, 1, -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, -1, 1, 1, -1, -1, -1, 1, 1, 1, -1, -1, -1, 1, -1, -1, -1, -1, 1, 1, -1, 1, 1, 1, -1, 1, 1, 1, -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, -1, -1, -1, 1, -1, -1, 1, 1, -1, 1, 1, 1, 1, -1, 1, -1, 1, 1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, -1, -1, -1, 1, 1, -1, -1, 1, 1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, -1, 1, 1, 1, 1, -1, 1, 1, 1, 1, 1, -1, -1, 1, -1, -1, -1, -1, -1, 1, -1, 1, -1, -1, -1, -1, -1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, -1, 1, 1, 1, 1, 1, 1, -1, -1, 1, -1, -1, -1, 1, -1, 1, 1, -1, -1, 1, 1, 1, -1, -1, 1, -1, -1, 1, 1, -1, -1, -1, -1, 1, 1, 1, -1, 1, -1, 1, 1, -1, 1, 1, -1, 1, -1, 1, -1, 1, 1, 1, -1, 1, 1, -1, -1, -1, 1, 1, 1, -1, -1, 1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1, 1, -1, 1, 1, 1, -1, 1, -1, -1, -1, 1, 1, 1, -1, 1, 1, -1, -1, -1, -1, 1, 1, -1, 1, 1, -1, -1, -1, -1, 1, -1, 1, 1, -1, -1, -1, -1, 1, -1, 1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, 1, -1, 1, 1, -1, -1, -1, 1, 1, -1, 1, -1, 1, -1, 1, 1, -1, 1, -1, -1, -1, -1, 1, -1, -1, -1, -1, -1, -1, 1, -1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, 1, -1, 1, -1, 1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, -1, -1, 1, -1, 1, -1, 1, -1, 1, -1, -1, -1, -1, 1, -1, 1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, -1, -1, 1, 1, 1, 1, 1, -1, -1, 1, -1, -1, -1, -1, 1, -1, 1, -1, 1, -1, 1, 1, -1, -1, 1, -1, -1, -1, 1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, -1, 1, -1, 1, 1, -1, -1, -1, -1, -1, -1, -1, 1, -1, -1, 1, -1, -1, -1, -1, 1, 1, -1, 1, 1, 1, 1, 1, 1, -1, -1, 1, -1, 1, -1, -1, 1, -1, 1, 1, -1, -1, 1, 1, 1, -1, 1, 1, -1, 1, -1, -1, 1, 1, -1, -1, -1, -1, -1, -1, 1, -1, -1, -1, 1, 1, 1, -1, 1, 1, 1, -1, -1, -1, -1, -1, 1, -1, 1, -1, -1, 1, -1, -1, -1, 1, -1, 1, -1, 1, 1, 1, -1, -1, 1, -1, -1, 1, -1, 1, 1, -1, 1, -1, 1, 1, 1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, -1, -1, 1, 1, 1, 1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, -1, -1, 1, 1, -1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1, 1, 1, 1, -1, 1, 1, -1, 1, -1, 1, -1, -1, -1, -1, -1, 1, 1, 1, -1, -1, -1, 1, 1, 1, -1, -1, -1, -1, 1, 1, -1, -1, -1, 1, 1, -1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1, 1, 1, 1, -1, -1, -1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, -1, 1, 1, -1, 1, -1, 1, 1, -1, 1, 1, 1, -1, 1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, 1, -1, 1, 1, -1, 1, -1, 1, 1, -1, 1, 1, -1, 1, -1, -1, -1, 1, -1, -1, -1, 1, 1, -1, 1, -1, -1, -1, 1, -1, 1, 1, -1, 1, -1, -1, -1, -1, -1, -1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, -1, 1, 1, 1, 1, -1, 1, 1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, -1, -1, 1, 1, 1, 1, 1, -1, 1, -1, -1, 1, -1, -1, -1, 1, 1, -1, -1, 1, -1, 1, -1, -1, 1, 1, -1, 1, 1, 1, 1, -1, 1, -1, -1, 1, 1, -1, 1, -1, -1, -1, 1, 1, 1, -1, -1, 1, 1, 1, -1, -1, -1, 1, -1, -1, 1, -1, -1, -1, 1, -1, 1, -1, 1, -1, 1, -1, 1, 1, 1, 1, -1, -1, 1, 1, -1, 1, -1, -1, -1, -1, -1, -1, 1, 1, 1, -1, -1, -1, -1, -1, 1, 1, 1, -1, -1, -1, 1, -1, -1, -1, 1, -1, 1, -1, -1, -1, 1, -1, -1, -1, -1, 1, -1, 1, -1, 1, -1, -1, -1, -1, -1, 1, 1, -1, 1, 1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1, 1, -1, 1, 1, 1, -1, 1, -1, -1, 1, -1, 1, 1, -1, -1, 1, -1, 1, -1, -1, 1, -1, 1, -1, -1, -1, -1, 1, -1, -1, 1, -1, -1, 1, 1, 1, -1, 1, 1, -1, -1, -1, -1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, 1, -1, 1, 1, 1, 1, 1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, 1, 1, -1, -1, 1, 1, -1, -1, -1, -1, 1, -1, -1, 1, 1, -1, 1, 1, -1, -1, -1, -1, 1, -1, -1, -1, -1, -1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, -1, -1, -1, -1, 1, 1, -1, 1, -1, 1, -1, 1, 1, -1, -1, -1, -1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, -1, -1, -1, 1, -1, -1, -1, 1, -1, -1, -1, -1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, 1, -1, -1, 1, 1, 1, 1, 1, 1, 1, -1, -1, -1, 1, 1, 1, 1, -1, -1, 1, 1, 1, 1, 1, 1, -1, -1, -1, 1, 1, 1, 1, 1, -1, 1, 1, 1, 1, 1, 1, -1, 1, -1, -1, -1, 1, 1, 1, -1, 1, 1, 1, -1, -1, 1, 1, 1, -1, -1, 1, 1, 1, 1, 1, -1, 1, -1, -1, 1, -1, 1, 1, 1, -1, 1, 1, 1, 1, -1, 1, 1, -1, -1, 1, -1, -1, -1, 1, -1, -1, 1, 1, -1, -1, -1, 1, -1, 1, 1, -1, 1, 1, -1, -1, 1, -1, -1, 1, -1, 1, -1, 1, -1, 1, -1, -1, -1, -1, 1, -1, -1, 1, -1, -1, 1, -1, 1, 1, 1, 1, -1, -1, 1, -1, 1, 1, -1, 1, -1, 1, -1, -1, -1, -1, 1, -1, -1, 1, 1, 1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, 1, 1, 1, 1, -1, -1, 1, -1, 1, 1, -1, 1, 1, 1, 1, -1, -1, -1, 1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, -1, -1, -1, 1, 1, 1, 1, -1, -1, 1, -1, -1, -1, -1, -1, -1, 1, -1, 1, 1, -1, -1, -1, 1, 1, 1, -1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, 1, 1, -1, -1, 1, -1, -1, -1, -1, 1, 1, 1, -1, -1, -1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, -1, -1, 1, -1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1, 1, 1, 1, 1, 1, -1, 1, -1, -1, -1, 1, -1, 1, -1, -1, -1, -1, -1, 1, -1, -1, -1, 1, 1, 1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1, 1, -1, -1, -1, 1, -1, -1, 1, 1, 1, 1, 1, -1, -1, 1, 1, 1, -1, -1, 1, -1, 1, 1, -1, -1, -1, 1, -1, -1, -1, 1, -1, 1, -1, 1, 1, -1, -1, -1, 1, 1, 1, -1, 1, 1, 1, 1, -1, -1, -1, -1, 1, -1, 1, -1, -1, 1, 1, 1, -1, -1, -1, 1, 1, -1, -1, 1, 1, 1, 1, 1, -1, -1, 1, 1, -1, -1, -1, 1, -1, -1, 1, 1, 1, -1, 1, 1, -1, -1, -1, -1, -1, 1, -1, 1, 1, -1, 1, 1, -1, -1, -1, -1, -1, 1, 1, -1, -1, 1, 1, 1, -1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, 1, -1, 1, -1, 1, 1, 1, -1, 1, -1, 1, 1, 1, -1, -1, -1, 1, 1, -1, -1, 1, 1, 1, -1, 1, 1, 1, -1, 1, 1, 1, 1, 1, 1, -1, -1, -1, -1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, 1, 1, -1, 1, 1, -1, 1, -1, 1, -1, 1, -1, -1, 1, 1, -1, 1, -1, 1, -1, -1, 1, -1, 1, -1, 1, -1, 1, -1, -1, 1, 1, -1, 1, 1, -1, -1, -1, 1, 1, -1, -1, -1, 1, 1, -1, -1, -1, 1, -1, -1, 1, -1, 1, 1, 1, -1, 1, 1, -1, -1, -1, -1, 1, 1, -1, -1, -1, 1, -1, -1, 1, 1, 1, 1, -1, 1, -1, -1, -1, -1, -1, 1, -1, 1, -1, -1, 1, 1, -1, -1, 1, -1, -1, -1, 1, -1, -1, 1, -1, 1, -1, 1, -1, -1, -1, 1, 1, 1, -1, 1, -1, 1, 1, -1, 1, 1, 1, -1, -1, 1, 1, -1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, 1, 1, -1, 1, -1, -1, -1, 1, 1, 1, -1, -1, -1, -1, -1, 1, -1, 1, 1, -1, 1, -1, -1, 1, 1, 1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, 1, 1, -1, 1, 1, -1, 1, 1, 1, 1, 1, 1, 1, -1, -1, -1, 1, -1, -1, -1, -1, -1, 1, -1, -1, 1, -1, 1, 1, -1, -1, -1, 1, 1, 1, 1, -1, -1, 1, -1, -1, -1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, 1, -1, -1, 1, -1, 1, 1, -1, 1, 1, 1, 1, 1, -1, 1, -1, -1, 1, 1, 1)
self.shseq_len = shseq_len
self.shseq_rep = shseq_rep
self.occupied_carriers = occupied_carriers
self.lgseq_len = lgseq_len
self.length_tag_name = length_tag_name = "packet_len"
self.fft_len = fft_len
self.dataseq_len = dataseq_len
self.sync_word2 = sync_word2 = (0, 0, 0, 0, 0, 1, 1, -1.0, -1, 1.0, 1, 1.0, -1, -1.0, -1, 1.0, 1, -1.0, 1, 1.0, 1, -1.0, -1, -1.0, -1, 1.0, -1, 1.0, -1, 1.0, 1, -1.0, 0, 1.0, 1, -1.0, 1, 1.0, -1, -1.0, 1, -1.0, -1, -1.0, 1, 1.0, 1, -1.0, 1, 1.0, -1, 1.0, -1, -1.0, -1, 1.0, 1, -1.0, 0, 0, 0, 0, 0, 0)
self.sync_word1 = sync_word1 = (0, 0, 0, 0, 0, 0, 0, -1.0, 0, 1.0, 0, 1.0, 0, -1.0, 0, 1.0, 0, -1.0, 0, 1.0, 0, -1.0, 0, -1.0, 0, 1.0, 0, 1.0, 0, 1.0, 0, -1.0, 0, 1.0, 0, -1.0, 0, 1.0, 0, -1.0, 0, -1.0, 0, -1.0, 0, 1.0, 0, -1.0, 0, 1.0, 0, 1.0, 0, -1.0, 0, 1.0, 0, -1.0, 0, 0, 0, 0, 0, 0)
self.sync_len = sync_len = 20
self.sig_ul_ratio_db = sig_ul_ratio_db
self.shseq = shseq = pn_symbols[pnseq_offset:pnseq_offset+shseq_len]
self.samp_rate = samp_rate
self.pilot_symbols = pilot_symbols
self.pilot_carriers = pilot_carriers
self.payload_mod = payload_mod = digital.constellation_qpsk()
self.packet_len = packet_len = len(occupied_carriers[0])
self.lgseq = lgseq = pn_symbols[pnseq_offset+shseq_len:pnseq_offset+shseq_len+lgseq_len]
self.header_formatter = header_formatter = digital.packet_header_ofdm(occupied_carriers, 1, length_tag_name)
self.final_gain = final_gain
self.usrp_gain = usrp_gain
self.dataseq = dataseq = pn_symbols[pnseq_offset+shseq_len+lgseq_len:pnseq_offset+shseq_len+lgseq_len+dataseq_len]
self.data_len = data_len = dataseq_len*20
self.cp_len = cp_len
self.cen_freq = cen_freq
##################################################
# Blocks
##################################################
self.uhd_usrp_sink = uhd.usrp_sink(
",".join(("", "")),
uhd.stream_args(
cpu_format="fc32",
channels=range(1),
),
)
self.uhd_usrp_sink.set_samp_rate(samp_rate)
self.uhd_usrp_sink.set_center_freq(cen_freq, 0)
self.uhd_usrp_sink.set_gain(usrp_gain, 0)
self.uhd_usrp_sink.set_antenna("TX/RX", 0)
self.ncofdm_ncofdm_carrier_allocator = ncofdm.ncofdm_carrier_allocator(fft_len, occupied_carriers, (), (), "packet_len")
self.ncofdm_add_cp_underlay = ncofdm.add_cp_underlay(fft_len, cp_len, sig_ul_ratio_db, data_len, shseq_len, shseq_rep, shseq, dataseq_len, dataseq, lgseq_len, lgseq, "packet_len")
self.main_source = blocks.vector_source_b(map(int, numpy.random.randint(0, 200, 10000)), True)
self.fft_vxx = fft.fft_vcc(fft_len, False, (()), True, 1)
self.digital_chunks_to_symbols = digital.chunks_to_symbols_bc((payload_mod.points()), 1)
self.blocks_stream_to_tagged_stream = blocks.stream_to_tagged_stream(gr.sizeof_char, 1, packet_len, "packet_len")
self.blocks_repack_bits = blocks.repack_bits_bb(8, payload_mod.bits_per_symbol(), length_tag_name, False, gr.GR_LSB_FIRST)
self.blocks_overlaygain = blocks.multiply_const_vcc(([1/(math.sqrt(44))]*fft_len))
self.blocks_final_gain = blocks.multiply_const_vcc((final_gain, ))
##################################################
# Connections
##################################################
self.connect((self.blocks_final_gain, 0), (self.uhd_usrp_sink, 0))
self.connect((self.blocks_overlaygain, 0), (self.ncofdm_add_cp_underlay, 0))
self.connect((self.blocks_repack_bits, 0), (self.digital_chunks_to_symbols, 0))
self.connect((self.blocks_stream_to_tagged_stream, 0), (self.blocks_repack_bits, 0))
self.connect((self.digital_chunks_to_symbols, 0), (self.ncofdm_ncofdm_carrier_allocator, 0))
self.connect((self.fft_vxx, 0), (self.blocks_overlaygain, 0))
self.connect((self.main_source, 0), (self.blocks_stream_to_tagged_stream, 0))
self.connect((self.ncofdm_add_cp_underlay, 0), (self.blocks_final_gain, 0))
self.connect((self.ncofdm_ncofdm_carrier_allocator, 0), (self.fft_vxx, 0))
def main():
#self.data = ( [[random.randint(0,255) for x in xrange(self.data_len)]
# +50*[0,],] )
#self.data = [[0 for x in xrange(self.data_len)],] \
# + self.nofdm_frames*[[random.randint(0,255) \
# for x in xrange(self.data_len)],]
#self.data = self.nofdm_frames*[[x for x in xrange(self.data_len)],]\
# + [50*[0,],]
args = get_arguments()
constellation = {
1:digital.constellation_bpsk(),
2:digital.constellation_qpsk(),
3:digital.constellation_8psk(),
}
packet_len_tag = "packet_length"
#fft_len = 64
#cp_len = 16
#occupied_carriers = (range(-26, -21) + range(-20, -7) +
# range(-6, 0) + range(1, 7) +
# range(8, 21) + range(22, 27),)
#pilot_carriers = ((-21, -7, 7, 21),)
#pilot_symbols = tuple([(1, -1, 1, -1),])
fft_len = 16
cp_len = 4
occupied_carriers = ((-5, -4, -2, -1, 1, 2, 4, 5),)
pilot_carriers = ((-3, 3),)
pilot_symbols = tuple([(1, -1),])
data_len = utils.ofdm_get_data_len(
nofdm_symbols=args.nsymbols,
noccupied_carriers=len(occupied_carriers[0]),
constellation=constellation[args.bits]
)
data = args.nframes*[[39 for x in xrange(data_len)],]
# sometimes Schmidl-Cox-Correlator ignores first frame
# so a dummy-frame is a good idea
if args.dummy_frame_start == True:
data.insert(0, data_len*[0])
# if file-output GNURadio needs extra frame at the and to loop over all
# OFDM-Frames before, last OFDM-Frame is ignored
# not so in case of using UHD-devices because there exists incoming
# input-data all the time
if args.dummy_frame_end == True:
data.append(data_len*[0])
tb = gr.top_block()
(data_tosend, tags) = packet_utils.packets_to_vectors(
data,
packet_len_tag
)
data_source = blocks.vector_source_b(
data=data_tosend,
vlen=1,
tags=tags,
repeat=False
)
ofdm_mapper = mimoots.ofdm_symbol_mapper_bc(
constellation=constellation[args.bits],
packet_len_tag=packet_len_tag,
verbose=args.verbose
)
ofdm_framer = mimoots.ofdm_symbols_to_frame_cvc(
fft_len=fft_len,
cp_len=cp_len,
occupied_carriers=occupied_carriers,
pilot_carriers=pilot_carriers,
pilot_symbols=pilot_symbols,
packet_len_tag=packet_len_tag,
verbose=args.verbose
)
ofdm_basebander = mimoots.ofdm_frames_to_basebandsignal_vcc(
fft_len=fft_len,
cp_len=cp_len,
packet_len_tag=packet_len_tag,
verbose=args.verbose
)
if args.freq == None:
data_sink = blocks.file_sink(
itemsize=gr.sizeof_gr_complex,
filename=args.to_file
)
else:
data_sink = mimoots.uhd_sink(freq=args.freq, gain=args.gain)
tb.connect(data_source, ofdm_mapper, ofdm_framer, ofdm_basebander,
data_sink)
tb.run()
# need to wait until the GNURadio-graph is finished
time.sleep(5)
def __init__(self):
gr.top_block.__init__(self, "OFDM Transceiver")
Qt.QWidget.__init__(self)
self.setWindowTitle("OFDM Transceiver")
try:
self.setWindowIcon(Qt.QIcon.fromTheme('gnuradio-grc'))
except:
pass
self.top_scroll_layout = Qt.QVBoxLayout()
self.setLayout(self.top_scroll_layout)
self.top_scroll = Qt.QScrollArea()
self.top_scroll.setFrameStyle(Qt.QFrame.NoFrame)
self.top_scroll_layout.addWidget(self.top_scroll)
self.top_scroll.setWidgetResizable(True)
self.top_widget = Qt.QWidget()
self.top_scroll.setWidget(self.top_widget)
self.top_layout = Qt.QVBoxLayout(self.top_widget)
self.top_grid_layout = Qt.QGridLayout()
self.top_layout.addLayout(self.top_grid_layout)
self.settings = Qt.QSettings("GNU Radio", "ofdm_transceiver")
self.restoreGeometry(self.settings.value("geometry").toByteArray())
self._lock = threading.RLock()
##################################################
# Variables
##################################################
self.pilot_symbols = pilot_symbols = ((1, 1, 1, -1,),)
self.pilot_carriers = pilot_carriers = ((-21, -7, 7, 21,),)
self.payload_mod = payload_mod = digital.constellation_qpsk()
self.packet_length_tag_key = packet_length_tag_key = "packet_len"
self.occupied_carriers = occupied_carriers = (range(-26, -21) + range(-20, -7) + range(-6, 0) + range(1, 7) + range(8, 21) + range(22, 27),)
self.length_tag_key = length_tag_key = "frame_len"
self.header_mod = header_mod = digital.constellation_bpsk()
self.fft_len = fft_len = 64
self.waterfall_min = waterfall_min = -140
self.waterfall_max = waterfall_max = -20
self.tx_gain = tx_gain = 0.01
self.th = th = 10**(-40/10)
self.sync_word2 = sync_word2 = [0j, 0j, 0j, 0j, 0j, 0j, (-1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1 +0j), (1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), 0j, (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), 0j, 0j, 0j, 0j, 0j]
self.sync_word1 = sync_word1 = [0., 0., 0., 0., 0., 0., 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 0., 0., 0., 0., 0.]
self.size = size = 256
self.samp_rate = samp_rate = 2e6
self.period = period = 10
self.payload_equalizer = payload_equalizer = digital.ofdm_equalizer_simpledfe(fft_len, payload_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols, 1)
self.nsamples = nsamples = 1000
self.header_formatter = header_formatter = digital.packet_header_ofdm(occupied_carriers, n_syms=1, len_tag_key=packet_length_tag_key, frame_len_tag_key=length_tag_key, bits_per_header_sym=header_mod.bits_per_symbol(), bits_per_payload_sym=payload_mod.bits_per_symbol(), scramble_header=False)
self.header_equalizer = header_equalizer = digital.ofdm_equalizer_simpledfe(fft_len, header_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols)
self.center_freq = center_freq = 2.412e9
##################################################
# Blocks
##################################################
self.tab = Qt.QTabWidget()
self.tab_widget_0 = Qt.QWidget()
self.tab_layout_0 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.tab_widget_0)
self.tab_grid_layout_0 = Qt.QGridLayout()
self.tab_layout_0.addLayout(self.tab_grid_layout_0)
self.tab.addTab(self.tab_widget_0, "Spectrum")
self.tab_widget_1 = Qt.QWidget()
self.tab_layout_1 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.tab_widget_1)
self.tab_grid_layout_1 = Qt.QGridLayout()
self.tab_layout_1.addLayout(self.tab_grid_layout_1)
self.tab.addTab(self.tab_widget_1, "Settings")
self.top_layout.addWidget(self.tab)
self._waterfall_min_range = Range(-200, -50, 1, -140, 200)
self._waterfall_min_win = RangeWidget(self._waterfall_min_range, self.set_waterfall_min, "Waterfall Min", "counter_slider", float)
self.tab_layout_1.addWidget(self._waterfall_min_win)
self._waterfall_max_range = Range(-50, 0, 1, -20, 200)
self._waterfall_max_win = RangeWidget(self._waterfall_max_range, self.set_waterfall_max, "Waterfall Max", "counter_slider", float)
self.tab_layout_1.addWidget(self._waterfall_max_win)
self._samp_rate_options = [1.5e6, 2e6, 4e6, 5e6, 10e6, 20e6 ]
self._samp_rate_labels = ["1.5 MHz", "2 MHz", "4 MHz", "5 MHz", "10 MHz", "20 MHz"]
self._samp_rate_tool_bar = Qt.QToolBar(self)
self._samp_rate_tool_bar.addWidget(Qt.QLabel("Sample Rate [MHz]"+": "))
self._samp_rate_combo_box = Qt.QComboBox()
self._samp_rate_tool_bar.addWidget(self._samp_rate_combo_box)
for label in self._samp_rate_labels: self._samp_rate_combo_box.addItem(label)
self._samp_rate_callback = lambda i: Qt.QMetaObject.invokeMethod(self._samp_rate_combo_box, "setCurrentIndex", Qt.Q_ARG("int", self._samp_rate_options.index(i)))
self._samp_rate_callback(self.samp_rate)
self._samp_rate_combo_box.currentIndexChanged.connect(
lambda i: self.set_samp_rate(self._samp_rate_options[i]))
self.tab_layout_0.addWidget(self._samp_rate_tool_bar)
self._center_freq_options = [868e6, 2.412e9, 2.417e9, 2.422e9, 2.427e9, 2.432e9, 2.437e9, 2.442e9, 2.447e9, 2.452e9, 2.457e9, 2.462e9, 2.467e9, 2.472e9]
self._center_freq_labels = ["868M", "1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13"]
self._center_freq_tool_bar = Qt.QToolBar(self)
self._center_freq_tool_bar.addWidget(Qt.QLabel("Channel"+": "))
self._center_freq_combo_box = Qt.QComboBox()
self._center_freq_tool_bar.addWidget(self._center_freq_combo_box)
for label in self._center_freq_labels: self._center_freq_combo_box.addItem(label)
self._center_freq_callback = lambda i: Qt.QMetaObject.invokeMethod(self._center_freq_combo_box, "setCurrentIndex", Qt.Q_ARG("int", self._center_freq_options.index(i)))
self._center_freq_callback(self.center_freq)
self._center_freq_combo_box.currentIndexChanged.connect(
lambda i: self.set_center_freq(self._center_freq_options[i]))
self.top_layout.addWidget(self._center_freq_tool_bar)
self._tx_gain_range = Range(0, 1, 0.001, 0.01, 200)
self._tx_gain_win = RangeWidget(self._tx_gain_range, self.set_tx_gain, "TX Gain", "counter_slider", float)
self.top_layout.addWidget(self._tx_gain_win)
self._th_range = Range(10**(-80/10), 10**(-10/10), 0.0001, 10**(-40/10), 200)
self._th_win = RangeWidget(self._th_range, self.set_th, "Threshold", "counter_slider", float)
self.tab_layout_0.addWidget(self._th_win)
self.qtgui_waterfall_sink_x_0_0 = qtgui.waterfall_sink_c(
1024, #size
firdes.WIN_BLACKMAN_hARRIS, #wintype
center_freq, #fc
samp_rate, #bw
"", #name
1 #number of inputs
)
self.qtgui_waterfall_sink_x_0_0.set_update_time(0.01)
self.qtgui_waterfall_sink_x_0_0.enable_grid(False)
if not True:
self.qtgui_waterfall_sink_x_0_0.disable_legend()
if complex == type(float()):
self.qtgui_waterfall_sink_x_0_0.set_plot_pos_half(not True)
labels = ["", "", "", "", "",
"", "", "", "", ""]
colors = [5, 0, 0, 0, 0,
0, 0, 0, 0, 0]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_waterfall_sink_x_0_0.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_waterfall_sink_x_0_0.set_line_label(i, labels[i])
self.qtgui_waterfall_sink_x_0_0.set_color_map(i, colors[i])
self.qtgui_waterfall_sink_x_0_0.set_line_alpha(i, alphas[i])
self.qtgui_waterfall_sink_x_0_0.set_intensity_range(waterfall_min, waterfall_max)
self._qtgui_waterfall_sink_x_0_0_win = sip.wrapinstance(self.qtgui_waterfall_sink_x_0_0.pyqwidget(), Qt.QWidget)
self.tab_layout_0.addWidget(self._qtgui_waterfall_sink_x_0_0_win)
self._period_range = Range(1, 10000, 1, 10, 200)
self._period_win = RangeWidget(self._period_range, self.set_period, "Period", "counter_slider", float)
self.tab_layout_1.addWidget(self._period_win)
self.osmosdr_source_0 = osmosdr.source( args="numchan=" + str(1) + " " + "bladerf=0" )
self.osmosdr_source_0.set_sample_rate(samp_rate)
self.osmosdr_source_0.set_center_freq(center_freq, 0)
self.osmosdr_source_0.set_freq_corr(0, 0)
self.osmosdr_source_0.set_dc_offset_mode(1, 0)
self.osmosdr_source_0.set_iq_balance_mode(1, 0)
self.osmosdr_source_0.set_gain_mode(False, 0)
self.osmosdr_source_0.set_gain(10, 0)
self.osmosdr_source_0.set_if_gain(20, 0)
self.osmosdr_source_0.set_bb_gain(20, 0)
self.osmosdr_source_0.set_antenna("", 0)
self.osmosdr_source_0.set_bandwidth(samp_rate, 0)
self.digital_ofdm_rx_0 = digital.ofdm_rx(
fft_len=64, cp_len=fft_len/4,
frame_length_tag_key='frame_'+packet_length_tag_key,
packet_length_tag_key=packet_length_tag_key,
occupied_carriers=occupied_carriers,
pilot_carriers=pilot_carriers,
pilot_symbols=pilot_symbols,
sync_word1=sync_word1,
sync_word2=sync_word2,
bps_header=1,
bps_payload=1,
debug_log=False,
scramble_bits=False
)
self.dc_blocker_xx_0 = filter.dc_blocker_cc(128, True)
self.blocks_tagged_stream_to_pdu_0 = blocks.tagged_stream_to_pdu(blocks.byte_t, packet_length_tag_key)
self.blocks_message_debug_0 = blocks.message_debug()
##################################################
# Connections
##################################################
self.msg_connect((self.blocks_tagged_stream_to_pdu_0, 'pdus'), (self.blocks_message_debug_0, 'print_pdu'))
self.connect((self.dc_blocker_xx_0, 0), (self.digital_ofdm_rx_0, 0))
self.connect((self.dc_blocker_xx_0, 0), (self.qtgui_waterfall_sink_x_0_0, 0))
self.connect((self.digital_ofdm_rx_0, 0), (self.blocks_tagged_stream_to_pdu_0, 0))
self.connect((self.osmosdr_source_0, 0), (self.dc_blocker_xx_0, 0))
def __init__(self): grc_wxgui.top_block_gui.__init__(self, title="OFDM Rx") _icon_path = "/usr/share/icons/hicolor/32x32/apps/gnuradio-grc.png" self.SetIcon(wx.Icon(_icon_path, wx.BITMAP_TYPE_ANY)) ################################################## # Variables ################################################## self.occupied_carriers = occupied_carriers = (range(-26, -21) + range(-20, -7) + range(-6, 0) + range(1, 7) + range(8, 21) + range(22, 27),) self.length_tag_name = length_tag_name = "frame_len" self.sync_word2 = sync_word2 = (0, 0, 0, 0, 0, 1, 1, -1.0, -1, 1.0, 1, 1.0, -1, -1.0, -1, 1.0, 1, -1.0, 1, 1.0, 1, -1.0, -1, -1.0, -1, 1.0, -1, 1.0, -1, 1.0, 1, -1.0, 0, 1.0, 1, -1.0, 1, 1.0, -1, -1.0, 1, -1.0, -1, -1.0, 1, 1.0, 1, -1.0, 1, 1.0, -1, 1.0, -1, -1.0, -1, 1.0, 1, -1.0, 0, 0, 0, 0, 0, 0) self.sync_word1 = sync_word1 = (0, 0, 0, 0, 0, 0, 0, -1.0, 0, 1.0, 0, 1.0, 0, -1.0, 0, 1.0, 0, -1.0, 0, 1.0, 0, -1.0, 0, -1.0, 0, 1.0, 0, 1.0, 0, 1.0, 0, -1.0, 0, 1.0, 0, -1.0, 0, 1.0, 0, -1.0, 0, -1.0, 0, -1.0, 0, 1.0, 0, -1.0, 0, 1.0, 0, 1.0, 0, -1.0, 0, 1.0, 0, -1.0, 0, 0, 0, 0, 0, 0) self.samp_rate = samp_rate = 3200000 self.pilot_symbols = pilot_symbols = ((1, 1, 1, -1,),) self.pilot_carriers = pilot_carriers = ((-21, -7, 7, 21,),) self.payload_mod = payload_mod = digital.constellation_qpsk() self.header_mod = header_mod = digital.constellation_bpsk() self.header_formatter = header_formatter = digital.packet_header_ofdm(occupied_carriers, 1, length_tag_name) self.fft_len = fft_len = 64 ################################################## # Blocks ################################################## self.gr_delay_0 = gr.delay(gr.sizeof_gr_complex*1, fft_len+fft_len/4) self.fft_vxx_0_0 = fft.fft_vcc(fft_len, True, (), True, 1) self.fft_vxx_0 = fft.fft_vcc(fft_len, True, (), True, 1) self.digital_packet_headerparser_b_0 = digital.packet_headerparser_b(header_formatter.formatter()) self.digital_ofdm_sync_sc_cfb_0 = digital.ofdm_sync_sc_cfb(fft_len, fft_len/4, False) self.digital_ofdm_serializer_vcc_1 = digital.ofdm_serializer_vcc(fft_len, occupied_carriers, "length_tag_key", "", 1, "", True) self.digital_ofdm_serializer_vcc_0 = digital.ofdm_serializer_vcc(fft_len, occupied_carriers, length_tag_name, "", 0, "", True) self.digital_ofdm_frame_equalizer_vcvc_0_0 = digital.ofdm_frame_equalizer_vcvc(digital.ofdm_equalizer_simpledfe(fft_len, header_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols).base(), fft_len/4, length_tag_name, True, 0) self.digital_ofdm_frame_equalizer_vcvc_0 = digital.ofdm_frame_equalizer_vcvc(digital.ofdm_equalizer_simpledfe(fft_len, header_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols, 2).base(), fft_len/4, "length_tag_key", False, 0) self.digital_ofdm_chanest_vcvc_0 = digital.ofdm_chanest_vcvc((sync_word1), (sync_word2), 2, 0, -1, False) self.digital_header_payload_demux_0 = digital.header_payload_demux(3, fft_len, fft_len/4, length_tag_name, "", True, gr.sizeof_gr_complex) self.digital_constellation_decoder_cb_0_0 = digital.constellation_decoder_cb(header_mod.base()) self.digital_constellation_decoder_cb_0 = digital.constellation_decoder_cb(payload_mod.base()) self.blocks_throttle_0 = blocks.throttle(gr.sizeof_gr_complex*1, samp_rate) self.blocks_tag_debug_0 = blocks.tag_debug(gr.sizeof_char*1, "Rx Packets") self.blocks_multiply_xx_0 = blocks.multiply_vcc(1) self.analog_noise_source_x_0 = analog.noise_source_c(analog.GR_GAUSSIAN, 1, 0) self.analog_frequency_modulator_fc_0 = analog.frequency_modulator_fc(-2.0/fft_len) ################################################## # Connections ################################################## self.connect((self.digital_ofdm_frame_equalizer_vcvc_0_0, 0), (self.digital_ofdm_serializer_vcc_0, 0)) self.connect((self.digital_header_payload_demux_0, 0), (self.fft_vxx_0, 0)) self.connect((self.fft_vxx_0, 0), (self.digital_ofdm_chanest_vcvc_0, 0)) self.connect((self.digital_ofdm_chanest_vcvc_0, 0), (self.digital_ofdm_frame_equalizer_vcvc_0_0, 0)) self.connect((self.digital_constellation_decoder_cb_0_0, 0), (self.digital_packet_headerparser_b_0, 0)) self.connect((self.digital_ofdm_serializer_vcc_0, 0), (self.digital_constellation_decoder_cb_0_0, 0)) self.connect((self.digital_constellation_decoder_cb_0, 0), (self.blocks_tag_debug_0, 0)) self.connect((self.fft_vxx_0_0, 0), (self.digital_ofdm_frame_equalizer_vcvc_0, 0)) self.connect((self.digital_ofdm_frame_equalizer_vcvc_0, 0), (self.digital_ofdm_serializer_vcc_1, 0)) self.connect((self.digital_header_payload_demux_0, 1), (self.fft_vxx_0_0, 0)) self.connect((self.digital_ofdm_serializer_vcc_1, 0), (self.digital_constellation_decoder_cb_0, 0)) self.connect((self.blocks_throttle_0, 0), (self.digital_ofdm_sync_sc_cfb_0, 0)) self.connect((self.blocks_throttle_0, 0), (self.gr_delay_0, 0)) self.connect((self.analog_noise_source_x_0, 0), (self.blocks_throttle_0, 0)) self.connect((self.analog_frequency_modulator_fc_0, 0), (self.blocks_multiply_xx_0, 0)) self.connect((self.digital_ofdm_sync_sc_cfb_0, 0), (self.analog_frequency_modulator_fc_0, 0)) self.connect((self.gr_delay_0, 0), (self.blocks_multiply_xx_0, 1)) self.connect((self.digital_ofdm_sync_sc_cfb_0, 1), (self.digital_header_payload_demux_0, 1)) self.connect((self.blocks_multiply_xx_0, 0), (self.digital_header_payload_demux_0, 0)) ################################################## # Asynch Message Connections ################################################## self.msg_connect(self.digital_packet_headerparser_b_0, "header_data", self.digital_header_payload_demux_0, "header_data")
def __init__(self):
gr.top_block.__init__(self, "Uhd Snr Receiver")
Qt.QWidget.__init__(self)
self.setWindowTitle("Uhd Snr Receiver")
self.setWindowIcon(Qt.QIcon.fromTheme('gnuradio-grc'))
self.top_scroll_layout = Qt.QVBoxLayout()
self.setLayout(self.top_scroll_layout)
self.top_scroll = Qt.QScrollArea()
self.top_scroll.setFrameStyle(Qt.QFrame.NoFrame)
self.top_scroll_layout.addWidget(self.top_scroll)
self.top_scroll.setWidgetResizable(True)
self.top_widget = Qt.QWidget()
self.top_scroll.setWidget(self.top_widget)
self.top_layout = Qt.QVBoxLayout(self.top_widget)
self.top_grid_layout = Qt.QGridLayout()
self.top_layout.addLayout(self.top_grid_layout)
self.settings = Qt.QSettings("GNU Radio", "uhd_snr_receiver")
self.restoreGeometry(self.settings.value("geometry").toByteArray())
##################################################
# Variables
##################################################
self.sps = sps = 4
self.nfilts = nfilts = 32
self.samp_rate = samp_rate = 1e6
self.rrc_taps = rrc_taps = filter.firdes.root_raised_cosine(nfilts, nfilts, 1.0/float(sps), 0.35, 11*sps*nfilts)
self.gain = gain = 15
self.freq = freq = 520e6
self.fine_freq = fine_freq = -28400
##################################################
# Blocks
##################################################
self._gain_layout = Qt.QVBoxLayout()
self._gain_tool_bar = Qt.QToolBar(self)
self._gain_layout.addWidget(self._gain_tool_bar)
self._gain_tool_bar.addWidget(Qt.QLabel("RX Gain"+": "))
self._gain_counter = Qwt.QwtCounter()
self._gain_counter.setRange(0, 31.5, 0.5)
self._gain_counter.setNumButtons(2)
self._gain_counter.setValue(self.gain)
self._gain_tool_bar.addWidget(self._gain_counter)
self._gain_counter.valueChanged.connect(self.set_gain)
self._gain_slider = Qwt.QwtSlider(None, Qt.Qt.Horizontal, Qwt.QwtSlider.BottomScale, Qwt.QwtSlider.BgSlot)
self._gain_slider.setRange(0, 31.5, 0.5)
self._gain_slider.setValue(self.gain)
self._gain_slider.setMinimumWidth(200)
self._gain_slider.valueChanged.connect(self.set_gain)
self._gain_layout.addWidget(self._gain_slider)
self.top_layout.addLayout(self._gain_layout)
self._freq_layout = Qt.QVBoxLayout()
self._freq_tool_bar = Qt.QToolBar(self)
self._freq_layout.addWidget(self._freq_tool_bar)
self._freq_tool_bar.addWidget(Qt.QLabel("Frequency"+": "))
self._freq_counter = Qwt.QwtCounter()
self._freq_counter.setRange(514e6, 526e6, 1e6)
self._freq_counter.setNumButtons(2)
self._freq_counter.setValue(self.freq)
self._freq_tool_bar.addWidget(self._freq_counter)
self._freq_counter.valueChanged.connect(self.set_freq)
self._freq_slider = Qwt.QwtSlider(None, Qt.Qt.Horizontal, Qwt.QwtSlider.BottomScale, Qwt.QwtSlider.BgSlot)
self._freq_slider.setRange(514e6, 526e6, 1e6)
self._freq_slider.setValue(self.freq)
self._freq_slider.setMinimumWidth(200)
self._freq_slider.valueChanged.connect(self.set_freq)
self._freq_layout.addWidget(self._freq_slider)
self.top_grid_layout.addLayout(self._freq_layout, 2,0,1,1)
self._fine_freq_layout = Qt.QVBoxLayout()
self._fine_freq_tool_bar = Qt.QToolBar(self)
self._fine_freq_layout.addWidget(self._fine_freq_tool_bar)
self._fine_freq_tool_bar.addWidget(Qt.QLabel("Fine Frequency"+": "))
self._fine_freq_counter = Qwt.QwtCounter()
self._fine_freq_counter.setRange(-50e3, 50e3, 100)
self._fine_freq_counter.setNumButtons(2)
self._fine_freq_counter.setValue(self.fine_freq)
self._fine_freq_tool_bar.addWidget(self._fine_freq_counter)
self._fine_freq_counter.valueChanged.connect(self.set_fine_freq)
self._fine_freq_slider = Qwt.QwtSlider(None, Qt.Qt.Horizontal, Qwt.QwtSlider.BottomScale, Qwt.QwtSlider.BgSlot)
self._fine_freq_slider.setRange(-50e3, 50e3, 100)
self._fine_freq_slider.setValue(self.fine_freq)
self._fine_freq_slider.setMinimumWidth(200)
self._fine_freq_slider.valueChanged.connect(self.set_fine_freq)
self._fine_freq_layout.addWidget(self._fine_freq_slider)
self.top_grid_layout.addLayout(self._fine_freq_layout, 2,1,1,1)
self.uhd_usrp_source_0 = uhd.usrp_source(
device_addr="",
stream_args=uhd.stream_args(
cpu_format="fc32",
channels=range(1),
),
)
self.uhd_usrp_source_0.set_samp_rate(samp_rate)
self.uhd_usrp_source_0.set_center_freq(freq + fine_freq, 0)
self.uhd_usrp_source_0.set_gain(gain, 0)
self.uhd_usrp_source_0.set_antenna("TX/RX", 0)
self.qtgui_time_sink_x_0_0 = qtgui.time_sink_c(
500, #size
samp_rate, #bw
"QT GUI Plot", #name
3 #number of inputs
)
self._qtgui_time_sink_x_0_0_win = sip.wrapinstance(self.qtgui_time_sink_x_0_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_time_sink_x_0_0_win, 0,0,1,1)
self.qtgui_freq_sink_x_0 = qtgui.freq_sink_c(
1024, #size
firdes.WIN_BLACKMAN_hARRIS, #wintype
0, #fc
samp_rate, #bw
"QT GUI Plot", #name
1 #number of inputs
)
self._qtgui_freq_sink_x_0_win = sip.wrapinstance(self.qtgui_freq_sink_x_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_freq_sink_x_0_win, 1,0,1,2)
self.qtgui_const_sink_x_0 = qtgui.const_sink_c(
1024, #size
"QT GUI Plot", #name
2 #number of inputs
)
self._qtgui_const_sink_x_0_win = sip.wrapinstance(self.qtgui_const_sink_x_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_const_sink_x_0_win, 0,1,1,1)
self.gr_multiply_xx_0 = gr.multiply_vcc(1)
self.gr_agc2_xx_0 = gr.agc2_cc(1e-1, 1e-2, 1.0, 1.0, 0.0)
self.digital_pfb_clock_sync_xxx_0 = digital.pfb_clock_sync_ccf(sps, 2*3.14/100.0, (rrc_taps), nfilts, nfilts/2, 1.5, 1)
self.digital_mpsk_snr_est_cc_0_1 = digital.mpsk_snr_est_cc(3, 10000, 0.001)
self.digital_mpsk_snr_est_cc_0_0 = digital.mpsk_snr_est_cc(2, 10000, 0.001)
self.digital_mpsk_snr_est_cc_0 = digital.mpsk_snr_est_cc(0, 10000, 0.001)
self.digital_lms_dd_equalizer_cc_0 = digital.lms_dd_equalizer_cc(15, 0.010, 1, digital.constellation_qpsk().base())
self.digital_costas_loop_cc_0 = digital.costas_loop_cc(2*3.14/100.0, 4)
##################################################
# Connections
##################################################
self.connect((self.digital_pfb_clock_sync_xxx_0, 0), (self.digital_mpsk_snr_est_cc_0, 0))
self.connect((self.gr_multiply_xx_0, 0), (self.qtgui_freq_sink_x_0, 0))
self.connect((self.digital_lms_dd_equalizer_cc_0, 0), (self.qtgui_const_sink_x_0, 1))
self.connect((self.gr_agc2_xx_0, 0), (self.gr_multiply_xx_0, 3))
self.connect((self.gr_agc2_xx_0, 0), (self.gr_multiply_xx_0, 2))
self.connect((self.gr_agc2_xx_0, 0), (self.gr_multiply_xx_0, 1))
self.connect((self.gr_agc2_xx_0, 0), (self.gr_multiply_xx_0, 0))
self.connect((self.gr_agc2_xx_0, 0), (self.digital_pfb_clock_sync_xxx_0, 0))
self.connect((self.digital_pfb_clock_sync_xxx_0, 0), (self.digital_costas_loop_cc_0, 0))
self.connect((self.digital_costas_loop_cc_0, 0), (self.digital_lms_dd_equalizer_cc_0, 0))
self.connect((self.digital_costas_loop_cc_0, 0), (self.qtgui_const_sink_x_0, 0))
self.connect((self.uhd_usrp_source_0, 0), (self.gr_agc2_xx_0, 0))
self.connect((self.digital_pfb_clock_sync_xxx_0, 0), (self.digital_mpsk_snr_est_cc_0_0, 0))
self.connect((self.digital_pfb_clock_sync_xxx_0, 0), (self.digital_mpsk_snr_est_cc_0_1, 0))
self.connect((self.digital_mpsk_snr_est_cc_0, 0), (self.qtgui_time_sink_x_0_0, 0))
self.connect((self.digital_mpsk_snr_est_cc_0_0, 0), (self.qtgui_time_sink_x_0_0, 1))
self.connect((self.digital_mpsk_snr_est_cc_0_1, 0), (self.qtgui_time_sink_x_0_0, 2))
def select_adaptive_mod(self):
# TODO: Do real adaptive modulation
self.header_const = digital.constellation_bpsk()
self.payload_const = digital.constellation_qpsk()
def __init__(
self, ipp1="127.0.0.1", ipp2="127.0.0.1", ipp3="127.0.0.1", ipp4="127.0.0.1", iptx="127.0.0.1", samp_rate=1000
):
gr.top_block.__init__(self, "OFDM Rx")
##################################################
# Parameters
##################################################
self.ipp1 = ipp1
self.ipp2 = ipp2
self.ipp3 = ipp3
self.ipp4 = ipp4
self.iptx = iptx
self.samp_rate = samp_rate
##################################################
# Variables
##################################################
self.pilot_symbols = pilot_symbols = ((1, 1, 1, -1),)
self.pilot_carriers = pilot_carriers = ((-21, -7, 7, 21),)
self.payload_mod = payload_mod = digital.constellation_qpsk()
self.packet_length_tag_key = packet_length_tag_key = "packet_len"
self.occupied_carriers = occupied_carriers = (
range(-26, -21) + range(-20, -7) + range(-6, 0) + range(1, 7) + range(8, 21) + range(22, 27),
)
self.length_tag_key = length_tag_key = "frame_len"
self.header_mod = header_mod = digital.constellation_bpsk()
self.fft_len = fft_len = 64
self.sync_word2 = sync_word2 = [
0j,
0j,
0j,
0j,
0j,
0j,
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(1 + 0j),
(1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(1 + 0j),
(-1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(1 + 0j),
(-1 + 0j),
(-1 + 0j),
(1 + 0j),
(-1 + 0j),
0j,
(1 + 0j),
(-1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(-1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(-1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(-1 + 0j),
(1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(1 + 0j),
(-1 + 0j),
(1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
0j,
0j,
0j,
0j,
0j,
]
self.sync_word1 = sync_word1 = [
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
-1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
-1.41421356,
0.0,
-1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
1.41421356,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
1.41421356,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
1.41421356,
0.0,
1.41421356,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
]
self.payload_equalizer = payload_equalizer = digital.ofdm_equalizer_simpledfe(
fft_len, payload_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols, 1
)
self.packet_len = packet_len = 96
self.header_formatter = header_formatter = digital.packet_header_ofdm(
occupied_carriers,
n_syms=1,
len_tag_key=packet_length_tag_key,
frame_len_tag_key=length_tag_key,
bits_per_header_sym=header_mod.bits_per_symbol(),
bits_per_payload_sym=payload_mod.bits_per_symbol(),
scramble_header=False,
)
self.header_equalizer = header_equalizer = digital.ofdm_equalizer_simpledfe(
fft_len, header_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols
)
##################################################
# Blocks
##################################################
self.zeromq_push_sink_0_0_1_0_0 = zeromq.push_sink(gr.sizeof_char, 1, "tcp://" + ipp3 + ":55530", 100, True)
self.zeromq_push_sink_0_0_1_0 = zeromq.push_sink(
gr.sizeof_gr_complex, 64, "tcp://" + ipp2 + ":55521", 100, True
)
self.zeromq_push_sink_0_0_1 = zeromq.push_sink(gr.sizeof_gr_complex, 64, "tcp://" + ipp2 + ":55520", 100, True)
self.zeromq_push_sink_0_0_0 = zeromq.push_sink(gr.sizeof_char, 1, "tcp://" + ipp1 + ":55511", 100, True)
self.zeromq_push_sink_0_0 = zeromq.push_sink(gr.sizeof_gr_complex, 1, "tcp://" + ipp1 + ":55510", 100, True)
self.zeromq_push_sink_0 = zeromq.push_sink(gr.sizeof_gr_complex, 1, "tcp://" + iptx + ":55500", 100, True)
self.zeromq_pull_source_0_0_0_0_0_0 = zeromq.pull_source(
gr.sizeof_gr_complex, 64, "tcp://" + ipp2 + ":55521", 100, True
)
self.zeromq_pull_source_0_0_0_0_0 = zeromq.pull_source(
gr.sizeof_gr_complex, 64, "tcp://" + ipp2 + ":55520", 100, True
)
self.zeromq_pull_source_0_0_0_0 = zeromq.pull_source(gr.sizeof_char, 1, "tcp://" + ipp3 + ":55530", 100, True)
self.zeromq_pull_source_0_0_0 = zeromq.pull_source(gr.sizeof_char, 1, "tcp://" + ipp1 + ":55511", 100, True)
self.zeromq_pull_source_0_0 = zeromq.pull_source(gr.sizeof_gr_complex, 1, "tcp://" + ipp1 + ":55510", 100, True)
self.zeromq_pull_source_0 = zeromq.pull_source(gr.sizeof_gr_complex, 1, "tcp://" + iptx + ":55500", 100, True)
self.my_random_source_limit_rate_0 = my.random_source_limit_rate(1000)
self.my_number_sync_timestamp_0 = my.number_sync_timestamp()
self.fft_vxx_1 = fft.fft_vcc(fft_len, True, (), True, 1)
self.fft_vxx_0 = fft.fft_vcc(fft_len, True, (()), True, 1)
self.digital_packet_headerparser_b_0 = digital.packet_headerparser_b(header_formatter.base())
self.digital_ofdm_tx_0 = digital.ofdm_tx(
fft_len=fft_len,
cp_len=fft_len / 4,
packet_length_tag_key=packet_length_tag_key,
occupied_carriers=occupied_carriers,
pilot_carriers=pilot_carriers,
pilot_symbols=pilot_symbols,
sync_word1=sync_word1,
sync_word2=sync_word2,
bps_header=1,
bps_payload=2,
rolloff=0,
debug_log=True,
scramble_bits=False,
)
self.digital_ofdm_sync_sc_cfb_0 = digital.ofdm_sync_sc_cfb(fft_len, fft_len / 4, False)
self.digital_ofdm_serializer_vcc_payload = digital.ofdm_serializer_vcc(
fft_len, occupied_carriers, length_tag_key, packet_length_tag_key, 1, "", True
)
self.digital_ofdm_serializer_vcc_header = digital.ofdm_serializer_vcc(
fft_len, occupied_carriers, length_tag_key, "", 0, "", True
)
self.digital_ofdm_frame_equalizer_vcvc_1 = digital.ofdm_frame_equalizer_vcvc(
payload_equalizer.base(), fft_len / 4, length_tag_key, True, 0
)
self.digital_ofdm_frame_equalizer_vcvc_0 = digital.ofdm_frame_equalizer_vcvc(
header_equalizer.base(), fft_len / 4, length_tag_key, True, 1
)
self.digital_ofdm_chanest_vcvc_0 = digital.ofdm_chanest_vcvc((sync_word1), (sync_word2), 1, 0, 3, False)
self.digital_header_payload_demux_0 = digital.header_payload_demux(
3, fft_len, fft_len / 4, length_tag_key, "", True, gr.sizeof_gr_complex, "rx_time", samp_rate, ()
)
self.digital_crc32_bb_0 = digital.crc32_bb(True, packet_length_tag_key, True)
self.digital_constellation_decoder_cb_1 = digital.constellation_decoder_cb(payload_mod.base())
self.digital_constellation_decoder_cb_0 = digital.constellation_decoder_cb(header_mod.base())
self.channels_channel_model_0 = channels.channel_model(
noise_voltage=0.1,
frequency_offset=0 * 1.0 / fft_len,
epsilon=1.0,
taps=(1.0,),
noise_seed=0,
block_tags=True,
)
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_gr_complex * 1, samp_rate, True)
self.blocks_tag_debug_1 = blocks.tag_debug(gr.sizeof_char * 1, "Rx Bytes", "")
self.blocks_tag_debug_1.set_display(True)
self.blocks_stream_to_tagged_stream_0 = blocks.stream_to_tagged_stream(
gr.sizeof_char, 1, packet_len, packet_length_tag_key
)
self.blocks_repack_bits_bb_0 = blocks.repack_bits_bb(
payload_mod.bits_per_symbol(), 8, packet_length_tag_key, True, gr.GR_LSB_FIRST
)
self.blocks_multiply_xx_0 = blocks.multiply_vcc(1)
self.blocks_delay_0 = blocks.delay(gr.sizeof_gr_complex * 1, fft_len + fft_len / 4)
self.analog_frequency_modulator_fc_0 = analog.frequency_modulator_fc(-2.0 / fft_len)
##################################################
# Connections
##################################################
self.msg_connect(
(self.digital_packet_headerparser_b_0, "header_data"), (self.digital_header_payload_demux_0, "header_data")
)
self.connect((self.analog_frequency_modulator_fc_0, 0), (self.blocks_multiply_xx_0, 0))
self.connect((self.blocks_delay_0, 0), (self.blocks_multiply_xx_0, 1))
self.connect((self.blocks_multiply_xx_0, 0), (self.zeromq_push_sink_0_0, 0))
self.connect((self.blocks_repack_bits_bb_0, 0), (self.digital_crc32_bb_0, 0))
self.connect((self.blocks_stream_to_tagged_stream_0, 0), (self.digital_ofdm_tx_0, 0))
self.connect((self.blocks_throttle_0, 0), (self.zeromq_push_sink_0, 0))
self.connect((self.channels_channel_model_0, 0), (self.blocks_throttle_0, 0))
self.connect((self.digital_constellation_decoder_cb_0, 0), (self.zeromq_push_sink_0_0_1_0_0, 0))
self.connect((self.digital_constellation_decoder_cb_1, 0), (self.blocks_repack_bits_bb_0, 0))
self.connect((self.digital_crc32_bb_0, 0), (self.blocks_tag_debug_1, 0))
self.connect((self.digital_crc32_bb_0, 0), (self.my_number_sync_timestamp_0, 0))
self.connect((self.digital_header_payload_demux_0, 0), (self.zeromq_push_sink_0_0_1, 0))
self.connect((self.digital_header_payload_demux_0, 1), (self.zeromq_push_sink_0_0_1_0, 0))
self.connect((self.digital_ofdm_chanest_vcvc_0, 0), (self.digital_ofdm_frame_equalizer_vcvc_0, 0))
self.connect((self.digital_ofdm_frame_equalizer_vcvc_0, 0), (self.digital_ofdm_serializer_vcc_header, 0))
self.connect((self.digital_ofdm_frame_equalizer_vcvc_1, 0), (self.digital_ofdm_serializer_vcc_payload, 0))
self.connect((self.digital_ofdm_serializer_vcc_header, 0), (self.digital_constellation_decoder_cb_0, 0))
self.connect((self.digital_ofdm_serializer_vcc_payload, 0), (self.digital_constellation_decoder_cb_1, 0))
self.connect((self.digital_ofdm_sync_sc_cfb_0, 0), (self.analog_frequency_modulator_fc_0, 0))
self.connect((self.digital_ofdm_sync_sc_cfb_0, 1), (self.zeromq_push_sink_0_0_0, 0))
self.connect((self.digital_ofdm_tx_0, 0), (self.channels_channel_model_0, 0))
self.connect((self.fft_vxx_0, 0), (self.digital_ofdm_chanest_vcvc_0, 0))
self.connect((self.fft_vxx_1, 0), (self.digital_ofdm_frame_equalizer_vcvc_1, 0))
self.connect((self.my_random_source_limit_rate_0, 0), (self.blocks_stream_to_tagged_stream_0, 0))
self.connect((self.zeromq_pull_source_0, 0), (self.blocks_delay_0, 0))
self.connect((self.zeromq_pull_source_0, 0), (self.digital_ofdm_sync_sc_cfb_0, 0))
self.connect((self.zeromq_pull_source_0_0, 0), (self.digital_header_payload_demux_0, 0))
self.connect((self.zeromq_pull_source_0_0_0, 0), (self.digital_header_payload_demux_0, 1))
self.connect((self.zeromq_pull_source_0_0_0_0, 0), (self.digital_packet_headerparser_b_0, 0))
self.connect((self.zeromq_pull_source_0_0_0_0_0, 0), (self.fft_vxx_0, 0))
self.connect((self.zeromq_pull_source_0_0_0_0_0_0, 0), (self.fft_vxx_1, 0))
def __init__(self, freq=0, gain=40, loopbw=100, loopbw_0=100, fllbw=0.002):
gr.top_block.__init__(self, "Rx Gui")
Qt.QWidget.__init__(self)
self.setWindowTitle("Rx Gui")
try:
self.setWindowIcon(Qt.QIcon.fromTheme('gnuradio-grc'))
except:
pass
self.top_scroll_layout = Qt.QVBoxLayout()
self.setLayout(self.top_scroll_layout)
self.top_scroll = Qt.QScrollArea()
self.top_scroll.setFrameStyle(Qt.QFrame.NoFrame)
self.top_scroll_layout.addWidget(self.top_scroll)
self.top_scroll.setWidgetResizable(True)
self.top_widget = Qt.QWidget()
self.top_scroll.setWidget(self.top_widget)
self.top_layout = Qt.QVBoxLayout(self.top_widget)
self.top_grid_layout = Qt.QGridLayout()
self.top_layout.addLayout(self.top_grid_layout)
self.settings = Qt.QSettings("GNU Radio", "rx_gui")
self.restoreGeometry(self.settings.value("geometry").toByteArray())
##################################################
# Parameters
##################################################
self.freq = freq
self.gain = gain
self.loopbw = loopbw
self.loopbw_0 = loopbw_0
self.fllbw = fllbw
##################################################
# Variables
##################################################
self.sps = sps = 8
self.excess_bw = excess_bw = 0.25
self.target_samp_rate = target_samp_rate = sps*(200e3/(1 + excess_bw))
self.qpsk_const = qpsk_const = digital.constellation_qpsk().base()
self.dsp_rate = dsp_rate = 100e6
self.const_choice = const_choice = "qpsk"
self.bpsk_const = bpsk_const = digital.constellation_bpsk().base()
self.barker_code_two_dim = barker_code_two_dim = [-1.0000 - 1.0000j, -1.0000 - 1.0000j, -1.0000 - 1.0000j, -1.0000 - 1.0000j, -1.0000 - 1.0000j, 1.0000 + 1.0000j, 1.0000 + 1.0000j, -1.0000 - 1.0000j, -1.0000 - 1.0000j, 1.0000 + 1.0000j, -1.0000 - 1.0000j, 1.0000 + 1.0000j, -1.0000 - 1.0000j]
self.barker_code_one_dim = barker_code_one_dim = sqrt(2)*numpy.real([-1.0000 - 1.0000j, -1.0000 - 1.0000j, -1.0000 - 1.0000j, -1.0000 - 1.0000j, -1.0000 - 1.0000j, 1.0000 + 1.0000j, 1.0000 + 1.0000j, -1.0000 - 1.0000j, -1.0000 - 1.0000j, 1.0000 + 1.0000j, -1.0000 - 1.0000j, 1.0000 + 1.0000j, -1.0000 - 1.0000j])
self.rrc_delay = rrc_delay = int(round(-44*excess_bw + 33))
self.nfilts = nfilts = 32
self.n_barker_rep = n_barker_rep = 10
self.dec_factor = dec_factor = ceil(dsp_rate/target_samp_rate)
self.constellation = constellation = qpsk_const if (const_choice=="qpsk") else bpsk_const
self.barker_code = barker_code = barker_code_two_dim if (const_choice == "qpsk") else barker_code_one_dim
self.preamble_syms = preamble_syms = numpy.matlib.repmat(barker_code, 1, n_barker_rep)[0]
self.n_rrc_taps = n_rrc_taps = rrc_delay * int(sps*nfilts)
self.n_codewords = n_codewords = 1
self.even_dec_factor = even_dec_factor = dec_factor if (dec_factor % 1 == 1) else (dec_factor+1)
self.const_order = const_order = pow(2,constellation.bits_per_symbol())
self.codeword_len = codeword_len = 18444
self.usrp_rx_addr = usrp_rx_addr = "192.168.10.2"
self.samp_rate = samp_rate = dsp_rate/even_dec_factor
self.rrc_taps = rrc_taps = firdes.root_raised_cosine(nfilts, nfilts*sps, 1.0, excess_bw, n_rrc_taps)
self.rf_center_freq = rf_center_freq = 1428.4309e6
self.preamble_size = preamble_size = len(preamble_syms)
self.pmf_peak_threshold = pmf_peak_threshold = 0.6
self.payload_size = payload_size = codeword_len*n_codewords/int(numpy.log2(const_order))
self.dataword_len = dataword_len = 6144
self.barker_len = barker_len = 13
##################################################
# Blocks
##################################################
self.tabs = Qt.QTabWidget()
self.tabs_widget_0 = Qt.QWidget()
self.tabs_layout_0 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.tabs_widget_0)
self.tabs_grid_layout_0 = Qt.QGridLayout()
self.tabs_layout_0.addLayout(self.tabs_grid_layout_0)
self.tabs.addTab(self.tabs_widget_0, 'PMF Out')
self.tabs_widget_1 = Qt.QWidget()
self.tabs_layout_1 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.tabs_widget_1)
self.tabs_grid_layout_1 = Qt.QGridLayout()
self.tabs_layout_1.addLayout(self.tabs_grid_layout_1)
self.tabs.addTab(self.tabs_widget_1, 'Abs PMF Out')
self.tabs_widget_2 = Qt.QWidget()
self.tabs_layout_2 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.tabs_widget_2)
self.tabs_grid_layout_2 = Qt.QGridLayout()
self.tabs_layout_2.addLayout(self.tabs_grid_layout_2)
self.tabs.addTab(self.tabs_widget_2, 'FLL In')
self.tabs_widget_3 = Qt.QWidget()
self.tabs_layout_3 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.tabs_widget_3)
self.tabs_grid_layout_3 = Qt.QGridLayout()
self.tabs_layout_3.addLayout(self.tabs_grid_layout_3)
self.tabs.addTab(self.tabs_widget_3, 'FLL Out')
self.tabs_widget_4 = Qt.QWidget()
self.tabs_layout_4 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.tabs_widget_4)
self.tabs_grid_layout_4 = Qt.QGridLayout()
self.tabs_layout_4.addLayout(self.tabs_grid_layout_4)
self.tabs.addTab(self.tabs_widget_4, 'FLL State')
self.tabs_widget_5 = Qt.QWidget()
self.tabs_layout_5 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.tabs_widget_5)
self.tabs_grid_layout_5 = Qt.QGridLayout()
self.tabs_layout_5.addLayout(self.tabs_grid_layout_5)
self.tabs.addTab(self.tabs_widget_5, 'PFB Sync Out')
self.tabs_widget_6 = Qt.QWidget()
self.tabs_layout_6 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.tabs_widget_6)
self.tabs_grid_layout_6 = Qt.QGridLayout()
self.tabs_layout_6.addLayout(self.tabs_grid_layout_6)
self.tabs.addTab(self.tabs_widget_6, 'Costas State')
self.tabs_widget_7 = Qt.QWidget()
self.tabs_layout_7 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.tabs_widget_7)
self.tabs_grid_layout_7 = Qt.QGridLayout()
self.tabs_layout_7.addLayout(self.tabs_grid_layout_7)
self.tabs.addTab(self.tabs_widget_7, 'Demod Bits')
self.tabs_widget_8 = Qt.QWidget()
self.tabs_layout_8 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.tabs_widget_8)
self.tabs_grid_layout_8 = Qt.QGridLayout()
self.tabs_layout_8.addLayout(self.tabs_grid_layout_8)
self.tabs.addTab(self.tabs_widget_8, 'Costas Sym Out')
self.tabs_widget_9 = Qt.QWidget()
self.tabs_layout_9 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.tabs_widget_9)
self.tabs_grid_layout_9 = Qt.QGridLayout()
self.tabs_layout_9.addLayout(self.tabs_grid_layout_9)
self.tabs.addTab(self.tabs_widget_9, 'Payload Symbols')
self.top_layout.addWidget(self.tabs)
self.rtlsdr_source_0 = osmosdr.source( args="numchan=" + str(1) + " " + '' )
self.rtlsdr_source_0.set_sample_rate(samp_rate)
self.rtlsdr_source_0.set_center_freq(freq, 0)
self.rtlsdr_source_0.set_freq_corr(0, 0)
self.rtlsdr_source_0.set_dc_offset_mode(0, 0)
self.rtlsdr_source_0.set_iq_balance_mode(0, 0)
self.rtlsdr_source_0.set_gain_mode(False, 0)
self.rtlsdr_source_0.set_gain(gain, 0)
self.rtlsdr_source_0.set_if_gain(20, 0)
self.rtlsdr_source_0.set_bb_gain(20, 0)
self.rtlsdr_source_0.set_antenna('', 0)
self.rtlsdr_source_0.set_bandwidth(0, 0)
self.qtgui_time_sink_x_2 = qtgui.time_sink_c(
preamble_size + payload_size, #size
samp_rate, #samp_rate
"", #name
1 #number of inputs
)
self.qtgui_time_sink_x_2.set_update_time(0.10)
self.qtgui_time_sink_x_2.set_y_axis(-1, 1)
self.qtgui_time_sink_x_2.set_y_label('Amplitude', "")
self.qtgui_time_sink_x_2.enable_tags(-1, True)
self.qtgui_time_sink_x_2.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, 0, "")
self.qtgui_time_sink_x_2.enable_autoscale(True)
self.qtgui_time_sink_x_2.enable_grid(False)
self.qtgui_time_sink_x_2.enable_axis_labels(True)
self.qtgui_time_sink_x_2.enable_control_panel(False)
if not True:
self.qtgui_time_sink_x_2.disable_legend()
labels = ['', '', '', '', '',
'', '', '', '', '']
widths = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
colors = ["blue", "red", "green", "black", "cyan",
"magenta", "yellow", "dark red", "dark green", "blue"]
styles = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1,
-1, -1, -1, -1, -1]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(2*1):
if len(labels[i]) == 0:
if(i % 2 == 0):
self.qtgui_time_sink_x_2.set_line_label(i, "Re{{Data {0}}}".format(i/2))
else:
self.qtgui_time_sink_x_2.set_line_label(i, "Im{{Data {0}}}".format(i/2))
else:
self.qtgui_time_sink_x_2.set_line_label(i, labels[i])
self.qtgui_time_sink_x_2.set_line_width(i, widths[i])
self.qtgui_time_sink_x_2.set_line_color(i, colors[i])
self.qtgui_time_sink_x_2.set_line_style(i, styles[i])
self.qtgui_time_sink_x_2.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_2.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_2_win = sip.wrapinstance(self.qtgui_time_sink_x_2.pyqwidget(), Qt.QWidget)
self.tabs_layout_0.addWidget(self._qtgui_time_sink_x_2_win)
self.qtgui_time_sink_x_1_0_0 = qtgui.time_sink_f(
1024, #size
samp_rate, #samp_rate
"", #name
1 #number of inputs
)
self.qtgui_time_sink_x_1_0_0.set_update_time(0.10)
self.qtgui_time_sink_x_1_0_0.set_y_axis(-128, 128)
self.qtgui_time_sink_x_1_0_0.set_y_label('Amplitude', "")
self.qtgui_time_sink_x_1_0_0.enable_tags(-1, False)
self.qtgui_time_sink_x_1_0_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, 0, "")
self.qtgui_time_sink_x_1_0_0.enable_autoscale(False)
self.qtgui_time_sink_x_1_0_0.enable_grid(False)
self.qtgui_time_sink_x_1_0_0.enable_axis_labels(True)
self.qtgui_time_sink_x_1_0_0.enable_control_panel(False)
if not True:
self.qtgui_time_sink_x_1_0_0.disable_legend()
labels = ['', '', '', '', '',
'', '', '', '', '']
widths = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
colors = ["blue", "red", "green", "black", "cyan",
"magenta", "yellow", "dark red", "dark green", "blue"]
styles = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1,
-1, -1, -1, -1, -1]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_time_sink_x_1_0_0.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_time_sink_x_1_0_0.set_line_label(i, labels[i])
self.qtgui_time_sink_x_1_0_0.set_line_width(i, widths[i])
self.qtgui_time_sink_x_1_0_0.set_line_color(i, colors[i])
self.qtgui_time_sink_x_1_0_0.set_line_style(i, styles[i])
self.qtgui_time_sink_x_1_0_0.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_1_0_0.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_1_0_0_win = sip.wrapinstance(self.qtgui_time_sink_x_1_0_0.pyqwidget(), Qt.QWidget)
self.tabs_layout_7.addWidget(self._qtgui_time_sink_x_1_0_0_win)
self.qtgui_time_sink_x_1_0 = qtgui.time_sink_f(
1024, #size
samp_rate, #samp_rate
"", #name
1 #number of inputs
)
self.qtgui_time_sink_x_1_0.set_update_time(0.10)
self.qtgui_time_sink_x_1_0.set_y_axis(-1, 1)
self.qtgui_time_sink_x_1_0.set_y_label('Amplitude', "")
self.qtgui_time_sink_x_1_0.enable_tags(-1, True)
self.qtgui_time_sink_x_1_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, 0, "")
self.qtgui_time_sink_x_1_0.enable_autoscale(False)
self.qtgui_time_sink_x_1_0.enable_grid(False)
self.qtgui_time_sink_x_1_0.enable_axis_labels(True)
self.qtgui_time_sink_x_1_0.enable_control_panel(False)
if not True:
self.qtgui_time_sink_x_1_0.disable_legend()
labels = ['', '', '', '', '',
'', '', '', '', '']
widths = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
colors = ["blue", "red", "green", "black", "cyan",
"magenta", "yellow", "dark red", "dark green", "blue"]
styles = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1,
-1, -1, -1, -1, -1]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_time_sink_x_1_0.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_time_sink_x_1_0.set_line_label(i, labels[i])
self.qtgui_time_sink_x_1_0.set_line_width(i, widths[i])
self.qtgui_time_sink_x_1_0.set_line_color(i, colors[i])
self.qtgui_time_sink_x_1_0.set_line_style(i, styles[i])
self.qtgui_time_sink_x_1_0.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_1_0.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_1_0_win = sip.wrapinstance(self.qtgui_time_sink_x_1_0.pyqwidget(), Qt.QWidget)
self.tabs_layout_7.addWidget(self._qtgui_time_sink_x_1_0_win)
self.qtgui_time_sink_x_1 = qtgui.time_sink_f(
8192, #size
samp_rate, #samp_rate
"", #name
3 #number of inputs
)
self.qtgui_time_sink_x_1.set_update_time(0.10)
self.qtgui_time_sink_x_1.set_y_axis(-1, 1)
self.qtgui_time_sink_x_1.set_y_label('Amplitude', "")
self.qtgui_time_sink_x_1.enable_tags(-1, True)
self.qtgui_time_sink_x_1.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, 0, "")
self.qtgui_time_sink_x_1.enable_autoscale(False)
self.qtgui_time_sink_x_1.enable_grid(False)
self.qtgui_time_sink_x_1.enable_axis_labels(True)
self.qtgui_time_sink_x_1.enable_control_panel(False)
if not True:
self.qtgui_time_sink_x_1.disable_legend()
labels = ['FLL Freq (PI Output)', 'FLL Phase Accum', 'FLL Error', '', '',
'', '', '', '', '']
widths = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
colors = ["blue", "red", "green", "black", "cyan",
"magenta", "yellow", "dark red", "dark green", "blue"]
styles = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1,
-1, -1, -1, -1, -1]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(3):
if len(labels[i]) == 0:
self.qtgui_time_sink_x_1.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_time_sink_x_1.set_line_label(i, labels[i])
self.qtgui_time_sink_x_1.set_line_width(i, widths[i])
self.qtgui_time_sink_x_1.set_line_color(i, colors[i])
self.qtgui_time_sink_x_1.set_line_style(i, styles[i])
self.qtgui_time_sink_x_1.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_1.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_1_win = sip.wrapinstance(self.qtgui_time_sink_x_1.pyqwidget(), Qt.QWidget)
self.tabs_layout_4.addWidget(self._qtgui_time_sink_x_1_win)
self.qtgui_time_sink_x_0_3 = qtgui.time_sink_f(
1024*4, #size
samp_rate, #samp_rate
"Error", #name
1 #number of inputs
)
self.qtgui_time_sink_x_0_3.set_update_time(0.10)
self.qtgui_time_sink_x_0_3.set_y_axis(-1, 1)
self.qtgui_time_sink_x_0_3.set_y_label('Amplitude', "")
self.qtgui_time_sink_x_0_3.enable_tags(-1, True)
self.qtgui_time_sink_x_0_3.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, 0, "")
self.qtgui_time_sink_x_0_3.enable_autoscale(False)
self.qtgui_time_sink_x_0_3.enable_grid(False)
self.qtgui_time_sink_x_0_3.enable_axis_labels(True)
self.qtgui_time_sink_x_0_3.enable_control_panel(False)
if not True:
self.qtgui_time_sink_x_0_3.disable_legend()
labels = ['', '', '', '', '',
'', '', '', '', '']
widths = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
colors = ["blue", "red", "green", "black", "cyan",
"magenta", "yellow", "dark red", "dark green", "blue"]
styles = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1,
-1, -1, -1, -1, -1]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_time_sink_x_0_3.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_time_sink_x_0_3.set_line_label(i, labels[i])
self.qtgui_time_sink_x_0_3.set_line_width(i, widths[i])
self.qtgui_time_sink_x_0_3.set_line_color(i, colors[i])
self.qtgui_time_sink_x_0_3.set_line_style(i, styles[i])
self.qtgui_time_sink_x_0_3.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_0_3.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_0_3_win = sip.wrapinstance(self.qtgui_time_sink_x_0_3.pyqwidget(), Qt.QWidget)
self.tabs_layout_6.addWidget(self._qtgui_time_sink_x_0_3_win)
self.qtgui_time_sink_x_0_0 = qtgui.time_sink_f(
preamble_size + payload_size, #size
samp_rate, #samp_rate
"", #name
1 #number of inputs
)
self.qtgui_time_sink_x_0_0.set_update_time(0.10)
self.qtgui_time_sink_x_0_0.set_y_axis(-1, 1)
self.qtgui_time_sink_x_0_0.set_y_label('Amplitude', "")
self.qtgui_time_sink_x_0_0.enable_tags(-1, True)
self.qtgui_time_sink_x_0_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, 0, "")
self.qtgui_time_sink_x_0_0.enable_autoscale(True)
self.qtgui_time_sink_x_0_0.enable_grid(False)
self.qtgui_time_sink_x_0_0.enable_axis_labels(True)
self.qtgui_time_sink_x_0_0.enable_control_panel(False)
if not True:
self.qtgui_time_sink_x_0_0.disable_legend()
labels = ['Mag Sq', 'Mag', '', '', '',
'', '', '', '', '']
widths = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
colors = ["blue", "red", "green", "black", "cyan",
"magenta", "yellow", "dark red", "dark green", "blue"]
styles = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1,
-1, -1, -1, -1, -1]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_time_sink_x_0_0.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_time_sink_x_0_0.set_line_label(i, labels[i])
self.qtgui_time_sink_x_0_0.set_line_width(i, widths[i])
self.qtgui_time_sink_x_0_0.set_line_color(i, colors[i])
self.qtgui_time_sink_x_0_0.set_line_style(i, styles[i])
self.qtgui_time_sink_x_0_0.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_0_0.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_0_0_win = sip.wrapinstance(self.qtgui_time_sink_x_0_0.pyqwidget(), Qt.QWidget)
self.tabs_layout_1.addWidget(self._qtgui_time_sink_x_0_0_win)
self.qtgui_sink_x_5 = qtgui.sink_c(
1024, #fftsize
firdes.WIN_BLACKMAN_hARRIS, #wintype
0, #fc
samp_rate, #bw
"", #name
True, #plotfreq
False, #plotwaterfall
False, #plottime
True, #plotconst
)
self.qtgui_sink_x_5.set_update_time(1.0/10)
self._qtgui_sink_x_5_win = sip.wrapinstance(self.qtgui_sink_x_5.pyqwidget(), Qt.QWidget)
self.tabs_layout_3.addWidget(self._qtgui_sink_x_5_win)
self.qtgui_sink_x_5.enable_rf_freq(False)
self.qtgui_sink_x_1 = qtgui.sink_c(
1024, #fftsize
firdes.WIN_BLACKMAN_hARRIS, #wintype
0, #fc
samp_rate, #bw
"", #name
True, #plotfreq
False, #plotwaterfall
False, #plottime
True, #plotconst
)
self.qtgui_sink_x_1.set_update_time(1.0/10)
self._qtgui_sink_x_1_win = sip.wrapinstance(self.qtgui_sink_x_1.pyqwidget(), Qt.QWidget)
self.tabs_layout_2.addWidget(self._qtgui_sink_x_1_win)
self.qtgui_sink_x_1.enable_rf_freq(False)
self.qtgui_sink_x_0 = qtgui.sink_c(
1024, #fftsize
firdes.WIN_BLACKMAN_hARRIS, #wintype
0, #fc
samp_rate, #bw
"", #name
False, #plotfreq
False, #plotwaterfall
False, #plottime
True, #plotconst
)
self.qtgui_sink_x_0.set_update_time(1.0/10)
self._qtgui_sink_x_0_win = sip.wrapinstance(self.qtgui_sink_x_0.pyqwidget(), Qt.QWidget)
self.tabs_layout_5.addWidget(self._qtgui_sink_x_0_win)
self.qtgui_sink_x_0.enable_rf_freq(False)
self.qtgui_freq_sink_x_0 = qtgui.freq_sink_c(
1024, #size
firdes.WIN_BLACKMAN_hARRIS, #wintype
0, #fc
samp_rate, #bw
"", #name
1 #number of inputs
)
self.qtgui_freq_sink_x_0.set_update_time(0.10)
self.qtgui_freq_sink_x_0.set_y_axis(-140, 10)
self.qtgui_freq_sink_x_0.set_y_label('Relative Gain', 'dB')
self.qtgui_freq_sink_x_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, 0.0, 0, "")
self.qtgui_freq_sink_x_0.enable_autoscale(False)
self.qtgui_freq_sink_x_0.enable_grid(False)
self.qtgui_freq_sink_x_0.set_fft_average(1.0)
self.qtgui_freq_sink_x_0.enable_axis_labels(True)
self.qtgui_freq_sink_x_0.enable_control_panel(False)
if not True:
self.qtgui_freq_sink_x_0.disable_legend()
if "complex" == "float" or "complex" == "msg_float":
self.qtgui_freq_sink_x_0.set_plot_pos_half(not True)
labels = ['', '', '', '', '',
'', '', '', '', '']
widths = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
colors = ["blue", "red", "green", "black", "cyan",
"magenta", "yellow", "dark red", "dark green", "dark blue"]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_freq_sink_x_0.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_freq_sink_x_0.set_line_label(i, labels[i])
self.qtgui_freq_sink_x_0.set_line_width(i, widths[i])
self.qtgui_freq_sink_x_0.set_line_color(i, colors[i])
self.qtgui_freq_sink_x_0.set_line_alpha(i, alphas[i])
self._qtgui_freq_sink_x_0_win = sip.wrapinstance(self.qtgui_freq_sink_x_0.pyqwidget(), Qt.QWidget)
self.tabs_layout_2.addWidget(self._qtgui_freq_sink_x_0_win)
self.qtgui_const_sink_x_1 = qtgui.const_sink_c(
1024, #size
"", #name
1 #number of inputs
)
self.qtgui_const_sink_x_1.set_update_time(0.10)
self.qtgui_const_sink_x_1.set_y_axis(-2, 2)
self.qtgui_const_sink_x_1.set_x_axis(-2, 2)
self.qtgui_const_sink_x_1.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, "")
self.qtgui_const_sink_x_1.enable_autoscale(False)
self.qtgui_const_sink_x_1.enable_grid(False)
self.qtgui_const_sink_x_1.enable_axis_labels(True)
if not True:
self.qtgui_const_sink_x_1.disable_legend()
labels = ['', '', '', '', '',
'', '', '', '', '']
widths = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
colors = ["blue", "red", "red", "red", "red",
"red", "red", "red", "red", "red"]
styles = [0, 0, 0, 0, 0,
0, 0, 0, 0, 0]
markers = [0, 0, 0, 0, 0,
0, 0, 0, 0, 0]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_const_sink_x_1.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_const_sink_x_1.set_line_label(i, labels[i])
self.qtgui_const_sink_x_1.set_line_width(i, widths[i])
self.qtgui_const_sink_x_1.set_line_color(i, colors[i])
self.qtgui_const_sink_x_1.set_line_style(i, styles[i])
self.qtgui_const_sink_x_1.set_line_marker(i, markers[i])
self.qtgui_const_sink_x_1.set_line_alpha(i, alphas[i])
self._qtgui_const_sink_x_1_win = sip.wrapinstance(self.qtgui_const_sink_x_1.pyqwidget(), Qt.QWidget)
self.tabs_layout_9.addWidget(self._qtgui_const_sink_x_1_win)
self.qtgui_const_sink_x_0 = qtgui.const_sink_c(
1024, #size
"", #name
1 #number of inputs
)
self.qtgui_const_sink_x_0.set_update_time(0.10)
self.qtgui_const_sink_x_0.set_y_axis(-2, 2)
self.qtgui_const_sink_x_0.set_x_axis(-2, 2)
self.qtgui_const_sink_x_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, "")
self.qtgui_const_sink_x_0.enable_autoscale(False)
self.qtgui_const_sink_x_0.enable_grid(False)
self.qtgui_const_sink_x_0.enable_axis_labels(True)
if not True:
self.qtgui_const_sink_x_0.disable_legend()
labels = ['', '', '', '', '',
'', '', '', '', '']
widths = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
colors = ["blue", "red", "red", "red", "red",
"red", "red", "red", "red", "red"]
styles = [0, 0, 0, 0, 0,
0, 0, 0, 0, 0]
markers = [0, 0, 0, 0, 0,
0, 0, 0, 0, 0]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_const_sink_x_0.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_const_sink_x_0.set_line_label(i, labels[i])
self.qtgui_const_sink_x_0.set_line_width(i, widths[i])
self.qtgui_const_sink_x_0.set_line_color(i, colors[i])
self.qtgui_const_sink_x_0.set_line_style(i, styles[i])
self.qtgui_const_sink_x_0.set_line_marker(i, markers[i])
self.qtgui_const_sink_x_0.set_line_alpha(i, alphas[i])
self._qtgui_const_sink_x_0_win = sip.wrapinstance(self.qtgui_const_sink_x_0.pyqwidget(), Qt.QWidget)
self.tabs_layout_8.addWidget(self._qtgui_const_sink_x_0_win)
self.mods_turbo_decoder_0 = mods.turbo_decoder(codeword_len, dataword_len)
self.mods_frame_sync_fast_0 = mods.frame_sync_fast(pmf_peak_threshold, preamble_size, payload_size, 0, 1, 1, int(const_order))
self.mods_fifo_async_sink_0 = mods.fifo_async_sink('/tmp/async_rx')
self.interp_fir_filter_xxx_0_0 = filter.interp_fir_filter_fff(1, ( numpy.ones(n_barker_rep*barker_len)))
self.interp_fir_filter_xxx_0_0.declare_sample_delay(0)
self.interp_fir_filter_xxx_0 = filter.interp_fir_filter_ccc(1, ( numpy.flipud(numpy.conj(preamble_syms))))
self.interp_fir_filter_xxx_0.declare_sample_delay(0)
self.framers_gr_hdlc_deframer_b_0 = framers.gr_hdlc_deframer_b(0)
self.digital_pfb_clock_sync_xxx_0 = digital.pfb_clock_sync_ccf(sps, 2*pi/50, (rrc_taps), nfilts, nfilts/2, pi/8, 1)
self.digital_map_bb_0_0_0 = digital.map_bb(([1,- 1]))
self.digital_fll_band_edge_cc_1 = digital.fll_band_edge_cc(sps, excess_bw, rrc_delay * int(sps) + 1, fllbw)
self.digital_descrambler_bb_0 = digital.descrambler_bb(0x21, 0x7F, 16)
self.digital_costas_loop_cc_0 = digital.costas_loop_cc(2*pi/loopbw, 2**constellation.bits_per_symbol(), False)
self.digital_constellation_decoder_cb_0 = digital.constellation_decoder_cb(constellation.base())
self.blocks_unpack_k_bits_bb_0 = blocks.unpack_k_bits_bb(constellation.bits_per_symbol())
self.blocks_rms_xx_1 = blocks.rms_cf(0.0001)
self.blocks_pack_k_bits_bb_1 = blocks.pack_k_bits_bb(8)
self.blocks_multiply_xx_0 = blocks.multiply_vff(1)
self.blocks_multiply_const_vxx_1_1 = blocks.multiply_const_vcc((1.0/sqrt(2), ))
self.blocks_multiply_const_vxx_1 = blocks.multiply_const_vcc((1.0/(preamble_size*sqrt(2)), ))
self.blocks_float_to_complex_0 = blocks.float_to_complex(1)
self.blocks_divide_xx_1 = blocks.divide_ff(1)
self.blocks_divide_xx_0 = blocks.divide_cc(1)
self.blocks_complex_to_mag_squared_0 = blocks.complex_to_mag_squared(1)
self.blocks_complex_to_mag_1 = blocks.complex_to_mag(1)
self.blocks_char_to_float_0_1 = blocks.char_to_float(1, 1)
self.blocks_char_to_float_0_0 = blocks.char_to_float(1, 1)
##################################################
# Connections
##################################################
self.msg_connect((self.framers_gr_hdlc_deframer_b_0, 'pdu'), (self.mods_fifo_async_sink_0, 'async_pdu'))
self.connect((self.blocks_char_to_float_0_0, 0), (self.qtgui_time_sink_x_1_0, 0))
self.connect((self.blocks_char_to_float_0_1, 0), (self.qtgui_time_sink_x_1_0_0, 0))
self.connect((self.blocks_complex_to_mag_1, 0), (self.blocks_divide_xx_1, 0))
self.connect((self.blocks_complex_to_mag_squared_0, 0), (self.interp_fir_filter_xxx_0_0, 0))
self.connect((self.blocks_divide_xx_0, 0), (self.digital_fll_band_edge_cc_1, 0))
self.connect((self.blocks_divide_xx_0, 0), (self.qtgui_sink_x_1, 0))
self.connect((self.blocks_divide_xx_1, 0), (self.blocks_multiply_xx_0, 0))
self.connect((self.blocks_divide_xx_1, 0), (self.blocks_multiply_xx_0, 1))
self.connect((self.blocks_float_to_complex_0, 0), (self.blocks_divide_xx_0, 1))
self.connect((self.blocks_multiply_const_vxx_1, 0), (self.mods_frame_sync_fast_0, 2))
self.connect((self.blocks_multiply_const_vxx_1, 0), (self.qtgui_time_sink_x_2, 0))
self.connect((self.blocks_multiply_const_vxx_1_1, 0), (self.blocks_complex_to_mag_1, 0))
self.connect((self.blocks_multiply_xx_0, 0), (self.mods_frame_sync_fast_0, 1))
self.connect((self.blocks_multiply_xx_0, 0), (self.qtgui_time_sink_x_0_0, 0))
self.connect((self.blocks_pack_k_bits_bb_1, 0), (self.blocks_char_to_float_0_1, 0))
self.connect((self.blocks_rms_xx_1, 0), (self.blocks_float_to_complex_0, 0))
self.connect((self.blocks_unpack_k_bits_bb_0, 0), (self.digital_map_bb_0_0_0, 0))
self.connect((self.digital_constellation_decoder_cb_0, 0), (self.blocks_unpack_k_bits_bb_0, 0))
self.connect((self.digital_costas_loop_cc_0, 0), (self.blocks_complex_to_mag_squared_0, 0))
self.connect((self.digital_costas_loop_cc_0, 0), (self.interp_fir_filter_xxx_0, 0))
self.connect((self.digital_costas_loop_cc_0, 0), (self.mods_frame_sync_fast_0, 0))
self.connect((self.digital_costas_loop_cc_0, 0), (self.qtgui_const_sink_x_0, 0))
self.connect((self.digital_costas_loop_cc_0, 1), (self.qtgui_time_sink_x_0_3, 0))
self.connect((self.digital_descrambler_bb_0, 0), (self.blocks_char_to_float_0_0, 0))
self.connect((self.digital_descrambler_bb_0, 0), (self.blocks_pack_k_bits_bb_1, 0))
self.connect((self.digital_descrambler_bb_0, 0), (self.framers_gr_hdlc_deframer_b_0, 0))
self.connect((self.digital_fll_band_edge_cc_1, 0), (self.digital_pfb_clock_sync_xxx_0, 0))
self.connect((self.digital_fll_band_edge_cc_1, 0), (self.qtgui_sink_x_5, 0))
self.connect((self.digital_fll_band_edge_cc_1, 3), (self.qtgui_time_sink_x_1, 2))
self.connect((self.digital_fll_band_edge_cc_1, 1), (self.qtgui_time_sink_x_1, 0))
self.connect((self.digital_fll_band_edge_cc_1, 2), (self.qtgui_time_sink_x_1, 1))
self.connect((self.digital_map_bb_0_0_0, 0), (self.mods_turbo_decoder_0, 0))
self.connect((self.digital_pfb_clock_sync_xxx_0, 0), (self.digital_costas_loop_cc_0, 0))
self.connect((self.digital_pfb_clock_sync_xxx_0, 0), (self.qtgui_sink_x_0, 0))
self.connect((self.interp_fir_filter_xxx_0, 0), (self.blocks_multiply_const_vxx_1, 0))
self.connect((self.interp_fir_filter_xxx_0, 0), (self.blocks_multiply_const_vxx_1_1, 0))
self.connect((self.interp_fir_filter_xxx_0_0, 0), (self.blocks_divide_xx_1, 1))
self.connect((self.mods_frame_sync_fast_0, 0), (self.digital_constellation_decoder_cb_0, 0))
self.connect((self.mods_frame_sync_fast_0, 0), (self.qtgui_const_sink_x_1, 0))
self.connect((self.mods_turbo_decoder_0, 0), (self.digital_descrambler_bb_0, 0))
self.connect((self.rtlsdr_source_0, 0), (self.blocks_divide_xx_0, 0))
self.connect((self.rtlsdr_source_0, 0), (self.blocks_rms_xx_1, 0))
self.connect((self.rtlsdr_source_0, 0), (self.qtgui_freq_sink_x_0, 0))
symbols_per_header = bits_per_header/header_mod.bits_per_symbol() if (1.0*bits_per_header)/header_mod.bits_per_symbol() != symbols_per_header: print "Error in evaluating symbols per header; adjusting bits per header" bits_per_header=(symbols_per_header+1)*header_mod.bits_per_symbol() symbols_per_header = bits_per_header/header_mod.bits_per_symbol() header_formatter = digital.packet_header_default(bits_per_header, length_tag_name,num_tag_name,header_mod.bits_per_symbol()); print "symbols_per_header=",symbols_per_header #payload info payload_bytes_per_frame = 50; crc_bytes=4; coded_payload_bytes_per_frame = payload_bytes_per_frame+crc_bytes payload_mod = digital.constellation_qpsk() coded_payload_symbols_per_frame = (coded_payload_bytes_per_frame * 8)/payload_mod.bits_per_symbol() if (coded_payload_bytes_per_frame * 8.0)/payload_mod.bits_per_symbol() != coded_payload_symbols_per_frame: print "Error in evaluating payload symbols per frame; adjusting payload bytes per frame" # add code to fix payload_bytes_per_frame print "coded_payload_symbols_per_frame=", coded_payload_symbols_per_frame symbols_per_frame = symbols_per_header + coded_payload_symbols_per_frame print "symbols_per_frame=", symbols_per_frame numpy.random.seed(666) training_long=[]; for i in range(symbols_per_frame): x=random.randint(0,1)*2-1
def __init__(self):
gr.top_block.__init__(self, "Send Test Rx")
##################################################
# Variables
##################################################
self.pilot_symbols = pilot_symbols = ((1, 1, 1, -1,),)
self.pilot_carriers = pilot_carriers = ((-21, -7, 7, 21,),)
self.payload_mod = payload_mod = digital.constellation_qpsk()
self.packet_length_tag_key = packet_length_tag_key = "packet_len"
self.occupied_carriers = occupied_carriers = (range(-26, -21) + range(-20, -7) + range(-6, 0) + range(1, 7) + range(8, 21) + range(22, 27),)
self.length_tag_key = length_tag_key = "frame_len"
self.header_mod = header_mod = digital.constellation_bpsk()
self.fft_len = fft_len = 64
self.sync_word2 = sync_word2 = [0j, 0j, 0j, 0j, 0j, 0j, (-1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1 +0j), (1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), 0j, (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (1+0j), (1+0j), (1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (1+0j), (-1+0j), (1+0j), (-1+0j), (-1+0j), (-1+0j), (-1+0j), 0j, 0j, 0j, 0j, 0j]
self.sync_word1 = sync_word1 = [0., 0., 0., 0., 0., 0., 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., -1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., -1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 1.41421356, 0., 0., 0., 0., 0., 0.]
self.samp_rate = samp_rate = 5000000
self.rolloff = rolloff = 0
self.payload_equalizer = payload_equalizer = digital.ofdm_equalizer_simpledfe(fft_len, payload_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols, 1)
self.packet_len = packet_len = 96
self.header_formatter = header_formatter = digital.packet_header_ofdm(occupied_carriers, n_syms=1, len_tag_key=packet_length_tag_key, frame_len_tag_key=length_tag_key, bits_per_header_sym=header_mod.bits_per_symbol(), bits_per_payload_sym=payload_mod.bits_per_symbol(), scramble_header=False)
self.header_equalizer = header_equalizer = digital.ofdm_equalizer_simpledfe(fft_len, header_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols)
##################################################
# Blocks
##################################################
self.uhd_usrp_source_0 = uhd.usrp_source(
",".join(("serial=30CEECB", "")),
uhd.stream_args(
cpu_format="fc32",
channels=range(1),
),
)
self.uhd_usrp_source_0.set_samp_rate(samp_rate)
self.uhd_usrp_source_0.set_center_freq(2.492e9, 0)
self.uhd_usrp_source_0.set_gain(40, 0)
self.uhd_usrp_source_0.set_antenna("TX/RX", 0)
self.pir_get_timestamp_delta_cm_0 = pir.get_timestamp_delta_cm()
self.fft_vxx_0 = fft.fft_vcc(fft_len, True, (()), True, 1)
self.digital_packet_headerparser_b_0 = digital.packet_headerparser_b(header_formatter.base())
self.digital_ofdm_sync_sc_cfb_0 = digital.ofdm_sync_sc_cfb(fft_len, fft_len/4, False)
self.digital_ofdm_serializer_vcc_header_0 = digital.ofdm_serializer_vcc(fft_len, occupied_carriers, length_tag_key, "", 0, "", True)
self.digital_ofdm_serializer_vcc_header = digital.ofdm_serializer_vcc(fft_len, occupied_carriers, length_tag_key, "", 0, "", True)
self.digital_ofdm_frame_equalizer_vcvc_0 = digital.ofdm_frame_equalizer_vcvc(header_equalizer.base(), fft_len/4, length_tag_key, True, 1)
self.digital_ofdm_chanest_vcvc_0 = digital.ofdm_chanest_vcvc((sync_word1), (sync_word2), 1, 0, 3, False)
self.digital_header_payload_demux_0 = digital.header_payload_demux(
3,
fft_len,
fft_len/4,
length_tag_key,
"",
True,
gr.sizeof_gr_complex,
"rx_time",
samp_rate,
(),
)
self.digital_corr_est_cc_0 = digital.corr_est_cc(((0.00000 - 0.00000j, -0.23756 + 0.06817j, 0.02455 - 0.06790j, -0.02923 + 0.05414j, 0.00701 - 0.10558j, 0.07080 - 0.00739j, 0.04334 + 0.00660j, -0.11691 + 0.08759j, 0.09375 + 0.06250j, 0.04313 - 0.10143j, 0.00862 - 0.03317j, -0.05702 + 0.07170j, -0.09857 - 0.09157j, 0.02654 + 0.03204j, -0.03675 - 0.04537j, 0.06443 + 0.17533j, 0.04419 + 0.00000j, 0.06443 - 0.17533j, -0.03675 + 0.04537j, 0.02654 - 0.03204j, -0.09857 + 0.09157j, -0.05702 - 0.07170j, 0.00862 + 0.03317j, 0.04313 + 0.10143j, 0.09375 - 0.06250j, -0.11691 - 0.08759j, 0.04334 - 0.00660j, 0.07080 + 0.00739j, 0.00701 + 0.10558j, -0.02923 - 0.05414j, 0.02455 + 0.06790j, -0.23756 - 0.06817j, 0.00000 + 0.00000j, 0.23756 - 0.06817j, -0.02455 + 0.06790j, 0.02923 - 0.05414j, -0.00701 + 0.10558j, -0.07080 + 0.00739j, -0.04334 - 0.00660j, 0.11691 - 0.08759j, -0.09375 - 0.06250j, -0.04313 + 0.10143j, -0.00862 + 0.03317j, 0.05702 - 0.07170j, 0.09857 + 0.09157j, -0.02654 - 0.03204j, 0.03675 + 0.04537j, -0.06443 - 0.17533j, -0.04419 + 0.00000j, -0.06443 + 0.17533j, 0.03675 - 0.04537j, -0.02654 + 0.03204j, 0.09857 - 0.09157j, 0.05702 + 0.07170j, -0.00862 - 0.03317j, -0.04313 - 0.10143j, -0.09375 + 0.06250j, 0.11691 + 0.08759j, -0.04334 + 0.00660j, -0.07080 - 0.00739j, -0.00701 - 0.10558j, 0.02923 + 0.05414j, -0.02455 - 0.06790j, 0.23756 + 0.06817j, 0.00000 - 0.00000j, -0.23756 + 0.06817j, 0.02455 - 0.06790j, -0.02923 + 0.05414j, 0.00701 - 0.10558j, 0.07080 - 0.00739j, 0.04334 + 0.00660j, -0.11691 + 0.08759j, 0.09375 + 0.06250j, 0.04313 - 0.10143j, 0.00862 - 0.03317j, -0.05702 + 0.07170j, -0.09857 - 0.09157j, 0.02654 + 0.03204j, -0.03675 - 0.04537j, 0.06443 + 0.17533j, 0.06250 + 0.03125j, -0.10790 - 0.08255j, 0.00182 + 0.03337j, 0.09586 - 0.07711j, 0.06872 + 0.05573j, 0.00373 + 0.04257j, -0.03565 - 0.03853j, -0.14193 + 0.01821j, -0.15089 - 0.03504j, 0.14492 - 0.01154j, 0.07984 + 0.05798j, -0.06544 - 0.02403j, 0.10927 - 0.05236j, 0.14303 + 0.12879j, -0.04601 - 0.13003j, -0.11594 + 0.05444j, -0.06250 + 0.00000j, -0.11594 - 0.05444j, -0.04601 + 0.13003j, 0.14303 - 0.12879j, 0.10927 + 0.05236j, -0.06544 + 0.02403j, 0.07984 - 0.05798j, 0.14492 + 0.01154j, -0.15089 + 0.03504j, -0.14193 - 0.01821j, -0.03565 + 0.03853j, 0.00373 - 0.04257j, 0.06872 - 0.05573j, 0.09586 + 0.07711j, 0.00182 - 0.03337j, -0.10790 + 0.08255j, 0.06250 - 0.03125j, -0.00527 - 0.11799j, 0.00182 + 0.05304j, 0.05923 + 0.02687j, -0.09460 + 0.11347j, 0.06227 - 0.14667j, -0.03565 + 0.03377j, -0.10488 - 0.05448j, 0.02589 + 0.09754j, 0.06134 - 0.07480j, 0.07984 + 0.06226j, -0.02766 + 0.02419j, 0.04162 - 0.02844j, -0.09425 - 0.14799j, -0.04601 + 0.09145j, 0.09289 + 0.12172j, -0.06250 + 0.00000j, 0.09289 - 0.12172j, -0.04601 - 0.09145j, -0.09425 + 0.14799j, 0.04162 + 0.02844j, -0.02766 - 0.02419j, 0.07984 - 0.06226j, 0.06134 + 0.07480j, 0.02589 - 0.09754j, -0.10488 + 0.05448j, -0.03565 - 0.03377j, 0.06227 + 0.14667j, -0.09460 - 0.11347j, 0.05923 - 0.02687j, 0.00182 - 0.05304j, -0.00527 + 0.11799j, 0.06250 + 0.03125j, -0.10790 - 0.08255j, 0.00182 + 0.03337j, 0.09586 - 0.07711j, 0.06872 + 0.05573j, 0.00373 + 0.04257j, -0.03565 - 0.03853j, -0.14193 + 0.01821j, -0.15089 - 0.03504j, 0.14492 - 0.01154j, 0.07984 + 0.05798j, -0.06544 - 0.02403j, 0.10927 - 0.05236j, 0.14303 + 0.12879j, -0.04601 - 0.13003j, -0.11594 + 0.05444j)), 1, 0, 0.999999999)
self.digital_constellation_decoder_cb_0 = digital.constellation_decoder_cb(header_mod.base())
self.blocks_tag_gate_0 = blocks.tag_gate(gr.sizeof_gr_complex * 1, False)
self.blocks_null_sink_0 = blocks.null_sink(gr.sizeof_gr_complex*1)
self.blocks_multiply_xx_0 = blocks.multiply_vcc(1)
self.blocks_message_debug_0 = blocks.message_debug()
self.blocks_delay_0 = blocks.delay(gr.sizeof_gr_complex*1, fft_len+fft_len/4)
self.analog_frequency_modulator_fc_0 = analog.frequency_modulator_fc(-2.0/fft_len)
##################################################
# Connections
##################################################
self.msg_connect((self.digital_packet_headerparser_b_0, 'header_data'), (self.digital_header_payload_demux_0, 'header_data'))
self.msg_connect((self.pir_get_timestamp_delta_cm_0, 'out'), (self.blocks_message_debug_0, 'print'))
self.connect((self.analog_frequency_modulator_fc_0, 0), (self.blocks_multiply_xx_0, 0))
self.connect((self.blocks_delay_0, 0), (self.blocks_multiply_xx_0, 1))
self.connect((self.blocks_multiply_xx_0, 0), (self.digital_corr_est_cc_0, 0))
self.connect((self.blocks_multiply_xx_0, 0), (self.digital_header_payload_demux_0, 0))
self.connect((self.blocks_tag_gate_0, 0), (self.blocks_delay_0, 0))
self.connect((self.blocks_tag_gate_0, 0), (self.digital_ofdm_sync_sc_cfb_0, 0))
self.connect((self.digital_constellation_decoder_cb_0, 0), (self.digital_packet_headerparser_b_0, 0))
self.connect((self.digital_corr_est_cc_0, 0), (self.pir_get_timestamp_delta_cm_0, 0))
self.connect((self.digital_header_payload_demux_0, 1), (self.digital_ofdm_serializer_vcc_header_0, 0))
self.connect((self.digital_header_payload_demux_0, 0), (self.fft_vxx_0, 0))
self.connect((self.digital_ofdm_chanest_vcvc_0, 0), (self.digital_ofdm_frame_equalizer_vcvc_0, 0))
self.connect((self.digital_ofdm_frame_equalizer_vcvc_0, 0), (self.digital_ofdm_serializer_vcc_header, 0))
self.connect((self.digital_ofdm_serializer_vcc_header, 0), (self.digital_constellation_decoder_cb_0, 0))
self.connect((self.digital_ofdm_serializer_vcc_header_0, 0), (self.blocks_null_sink_0, 0))
self.connect((self.digital_ofdm_sync_sc_cfb_0, 0), (self.analog_frequency_modulator_fc_0, 0))
self.connect((self.digital_ofdm_sync_sc_cfb_0, 1), (self.digital_header_payload_demux_0, 1))
self.connect((self.fft_vxx_0, 0), (self.digital_ofdm_chanest_vcvc_0, 0))
self.connect((self.uhd_usrp_source_0, 0), (self.blocks_tag_gate_0, 0))
def __init__(self):
gr.top_block.__init__(self, "Bpsk Sender1")
##################################################
# Variables
##################################################
self.sps = sps = 8
self.nfilts = nfilts = 32
self.transistion = transistion = 100
self.sideband_rx = sideband_rx = 500
self.sideband = sideband = 500
self.samp_rate = samp_rate = 48000
self.rrc_taps = rrc_taps = firdes.root_raised_cosine(nfilts, nfilts, 1.0/float(sps), 0.35, 11*sps*nfilts)
self.qpsk = qpsk = digital.constellation_rect(([0.707+0.707j, -0.707+0.707j, -0.707-0.707j, 0.707-0.707j]), ([0, 1, 2, 3]), 4, 2, 2, 1, 1).base()
self.preamble = preamble = [1,-1,1,-1,1,1,-1,-1,1,1,-1,1,1,1,-1,1,1,-1,1,-1,-1,1,-1,-1,1,1,1,-1,-1,-1,1,-1,1,1,1,1,-1,-1,1,-1,1,-1,-1,-1,1,1,-1,-1,-1,-1,1,-1,-1,-1,-1,-1,1,1,1,1,1,1,-1,-1]
self.payload_size = payload_size = 10
self.interpolation = interpolation = 2000
self.gap = gap = 0
self.eb = eb = 0.35
self.delay = delay = 0
self.decimation = decimation = 1
self.constel = constel = digital.constellation_calcdist(([1,- 1]), ([0,1]), 2, 1).base()
self.const_type = const_type = 1
self.const = const = (digital.constellation_bpsk(), digital.constellation_qpsk(), digital.constellation_8psk())
self.carrier = carrier = 10000
self.arity = arity = 2
##################################################
# Blocks
##################################################
self.rational_resampler_xxx_0 = filter.rational_resampler_ccc(
interpolation=interpolation,
decimation=decimation,
taps=None,
fractional_bw=None,
)
self.freq_xlating_fir_filter_xxx_0 = filter.freq_xlating_fir_filter_ccc(1, (firdes.band_pass (0.5,samp_rate,10000-sideband,10000+sideband,transistion)), -carrier, samp_rate)
self.digital_constellation_modulator_0 = digital.generic_mod(
constellation=constel,
differential=True,
samples_per_symbol=sps,
pre_diff_code=True,
excess_bw=0.35,
verbose=False,
log=False,
)
self.blocks_wavfile_sink_1 = blocks.wavfile_sink("BPSK_output.wav", 1, 48000, 16)
self.blocks_unpack_k_bits_bb_0_0 = blocks.unpack_k_bits_bb(8)
self.blocks_multiply_const_vxx_0 = blocks.multiply_const_vff((0.1, ))
self.blocks_file_source_0 = blocks.file_source(gr.sizeof_char*1, "TestData2", False)
self.blocks_file_sink_0_0_0 = blocks.file_sink(gr.sizeof_char*1, "inputBinary", False)
self.blocks_file_sink_0_0_0.set_unbuffered(False)
self.blocks_file_sink_0_0 = blocks.file_sink(gr.sizeof_char*1, "inputBinary2", False)
self.blocks_file_sink_0_0.set_unbuffered(False)
script, sdelay= argv
self.blocks_delay_0_0 = blocks.delay(gr.sizeof_char*1, int(sdelay))
self.blocks_complex_to_real_0 = blocks.complex_to_real(1)
##################################################
# Connections
##################################################
self.connect((self.blocks_complex_to_real_0, 0), (self.blocks_multiply_const_vxx_0, 0))
self.connect((self.blocks_delay_0_0, 0), (self.blocks_file_sink_0_0, 0))
self.connect((self.blocks_delay_0_0, 0), (self.blocks_file_sink_0_0_0, 0))
self.connect((self.blocks_file_source_0, 0), (self.blocks_unpack_k_bits_bb_0_0, 0))
self.connect((self.blocks_file_source_0, 0), (self.digital_constellation_modulator_0, 0))
self.connect((self.blocks_multiply_const_vxx_0, 0), (self.blocks_wavfile_sink_1, 0))
self.connect((self.blocks_unpack_k_bits_bb_0_0, 0), (self.blocks_delay_0_0, 0))
self.connect((self.digital_constellation_modulator_0, 0), (self.rational_resampler_xxx_0, 0))
self.connect((self.freq_xlating_fir_filter_xxx_0, 0), (self.blocks_complex_to_real_0, 0))
self.connect((self.rational_resampler_xxx_0, 0), (self.freq_xlating_fir_filter_xxx_0, 0))
def __init__(
self, ipp1="127.0.0.1", ipp2="127.0.0.1", ipp3="127.0.0.1", ipp4="127.0.0.1", iptx="127.0.0.1", samp_rate=10000
):
gr.top_block.__init__(self, "OFDM Rx")
##################################################
# Parameters
##################################################
self.ipp1 = ipp1
self.ipp2 = ipp2
self.ipp3 = ipp3
self.ipp4 = ipp4
self.iptx = iptx
self.samp_rate = samp_rate
##################################################
# Variables
##################################################
self.pilot_symbols = pilot_symbols = ((1, 1, 1, -1),)
self.pilot_carriers = pilot_carriers = ((-21, -7, 7, 21),)
self.payload_mod = payload_mod = digital.constellation_qpsk()
self.packet_length_tag_key = packet_length_tag_key = "packet_len"
self.occupied_carriers = occupied_carriers = (
range(-26, -21) + range(-20, -7) + range(-6, 0) + range(1, 7) + range(8, 21) + range(22, 27),
)
self.length_tag_key = length_tag_key = "frame_len"
self.header_mod = header_mod = digital.constellation_bpsk()
self.fft_len = fft_len = 64
self.sync_word2 = sync_word2 = [
0j,
0j,
0j,
0j,
0j,
0j,
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(1 + 0j),
(1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(1 + 0j),
(-1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(1 + 0j),
(-1 + 0j),
(-1 + 0j),
(1 + 0j),
(-1 + 0j),
0j,
(1 + 0j),
(-1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(-1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(-1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(1 + 0j),
(-1 + 0j),
(1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(1 + 0j),
(-1 + 0j),
(1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
(-1 + 0j),
0j,
0j,
0j,
0j,
0j,
]
self.sync_word1 = sync_word1 = [
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
-1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
-1.41421356,
0.0,
-1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
1.41421356,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
1.41421356,
0.0,
1.41421356,
0.0,
-1.41421356,
0.0,
1.41421356,
0.0,
1.41421356,
0.0,
1.41421356,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
]
self.payload_equalizer = payload_equalizer = digital.ofdm_equalizer_simpledfe(
fft_len, payload_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols, 1
)
self.packet_len = packet_len = 96
self.header_formatter = header_formatter = digital.packet_header_ofdm(
occupied_carriers,
n_syms=1,
len_tag_key=packet_length_tag_key,
frame_len_tag_key=length_tag_key,
bits_per_header_sym=header_mod.bits_per_symbol(),
bits_per_payload_sym=payload_mod.bits_per_symbol(),
scramble_header=False,
)
self.header_equalizer = header_equalizer = digital.ofdm_equalizer_simpledfe(
fft_len, header_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols
)
##################################################
# Blocks
##################################################
self.zeromq_push_sink_0 = zeromq.push_sink(gr.sizeof_gr_complex, 1, "tcp://" + iptx + ":55500", 100, True)
self.my_random_source_limit_rate_0 = my.random_source_limit_rate(10000000)
self.my_number_sync_timestamp_0 = my.number_sync_timestamp()
self.digital_ofdm_tx_0 = digital.ofdm_tx(
fft_len=fft_len,
cp_len=fft_len / 4,
packet_length_tag_key=packet_length_tag_key,
occupied_carriers=occupied_carriers,
pilot_carriers=pilot_carriers,
pilot_symbols=pilot_symbols,
sync_word1=sync_word1,
sync_word2=sync_word2,
bps_header=1,
bps_payload=2,
rolloff=0,
debug_log=True,
scramble_bits=False,
)
self.channels_channel_model_0 = channels.channel_model(
noise_voltage=0.0, frequency_offset=0.0, epsilon=1.0, taps=(1.0,), noise_seed=0, block_tags=False
)
self.blocks_stream_to_tagged_stream_0 = blocks.stream_to_tagged_stream(
gr.sizeof_char, 1, packet_len, packet_length_tag_key
)
(self.blocks_stream_to_tagged_stream_0).set_max_output_buffer(1)
##################################################
# Connections
##################################################
self.connect((self.blocks_stream_to_tagged_stream_0, 0), (self.digital_ofdm_tx_0, 0))
self.connect((self.blocks_stream_to_tagged_stream_0, 0), (self.my_number_sync_timestamp_0, 0))
self.connect((self.channels_channel_model_0, 0), (self.zeromq_push_sink_0, 0))
self.connect((self.digital_ofdm_tx_0, 0), (self.channels_channel_model_0, 0))
self.connect((self.my_random_source_limit_rate_0, 0), (self.blocks_stream_to_tagged_stream_0, 0))
