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 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 __init__(self, varicode_decode=True, differential_decode=True):
gr.hier_block2.__init__(
self,
"Coherent PSK31 Demodulator",
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
gr.io_signature(1, 1, gr.sizeof_char * 1),
)
our_blocks = [self]
constellation = digital.constellation_bpsk().base()
our_blocks.append(digital.constellation_decoder_cb(constellation))
if differential_decode:
# PSK31 defines 0 as a phase change, opposite the usual
# differential encoding which is: out = (next - prev) % 2
our_blocks.extend([
digital.diff_decoder_bb(2),
blocks.not_bb(),
blocks.and_const_bb(1),
])
if varicode_decode:
our_blocks.append(varicode_decode_bb())
our_blocks.append(self)
self.connect(*our_blocks)
def __init__(self):
gr.top_block.__init__(self, "Gotenna Rx Usrp")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 32000000
self.fsk_deviation_hz = fsk_deviation_hz = 12500
self.chan_spacing = chan_spacing = 500000
self.baud_rate = baud_rate = 24000
##################################################
# Blocks
##################################################
self.uhd_usrp_source_0 = uhd.usrp_source(
",".join(("", '')),
uhd.stream_args(
cpu_format="fc32",
args='',
channels=list(range(0,1)),
),
)
self.uhd_usrp_source_0.set_center_freq(915000000, 0)
self.uhd_usrp_source_0.set_gain(5, 0)
self.uhd_usrp_source_0.set_antenna('TX/RX', 0)
self.uhd_usrp_source_0.set_samp_rate(samp_rate)
self.uhd_usrp_source_0.set_time_unknown_pps(uhd.time_spec())
self.rational_resampler_xxx_0 = filter.rational_resampler_fff(
interpolation=1,
decimation=4,
taps=None,
fractional_bw=None)
self.gotenna_sink = gotenna_sink.blk()
self.digital_symbol_sync_xx_0 = digital.symbol_sync_ff(
digital.TED_DANDREA_AND_MENGALI_GEN_MSK,
float(chan_spacing) / baud_rate / 4,
0.05,
1.5,
1.0,
0.001 * float(chan_spacing) / baud_rate / 4,
1,
digital.constellation_bpsk().base(),
digital.IR_MMSE_8TAP,
128,
[])
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.blocks_keep_one_in_n_0 = blocks.keep_one_in_n(gr.sizeof_gr_complex*1, samp_rate // chan_spacing)
self.analog_quadrature_demod_cf_0 = analog.quadrature_demod_cf(chan_spacing/(2*math.pi*fsk_deviation_hz))
##################################################
# Connections
##################################################
self.connect((self.analog_quadrature_demod_cf_0, 0), (self.rational_resampler_xxx_0, 0))
self.connect((self.blocks_keep_one_in_n_0, 0), (self.analog_quadrature_demod_cf_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.gotenna_sink, 0))
self.connect((self.digital_symbol_sync_xx_0, 0), (self.digital_binary_slicer_fb_0, 0))
self.connect((self.rational_resampler_xxx_0, 0), (self.digital_symbol_sync_xx_0, 0))
self.connect((self.uhd_usrp_source_0, 0), (self.blocks_keep_one_in_n_0, 0))
def __init__(self, ebn0, nbits):
gr.hier_block2.__init__(self, 'BPSK',
gr.io_signature(1, 1, gr.sizeof_char),
gr.io_signature(1, 1, gr.sizeof_char))
samp_rate = 48e3
sps = 5
self.head = blocks.head(gr.sizeof_char, nbits)
self.pack = blocks.pack_k_bits_bb(8)
self.bpsk_constellation = digital.constellation_bpsk().base()
self.modulator = digital.generic_mod(
constellation = self.bpsk_constellation,
differential = False,
samples_per_symbol = sps,
pre_diff_code = True,
excess_bw = 0.35,
verbose = False,
log = False)
spb = sps
self.channel = channels.channel_model(np.sqrt(spb)/10**(ebn0/20), 0, 1.0, [1], RAND_SEED, False)
self.demod = bpsk_demodulator(samp_rate/sps, samp_rate, iq = True)
self.slice = digital.binary_slicer_fb()
self.connect(self, self.head, self.pack, self.modulator, self.channel,
self.demod, self.slice, self)
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 get_constellation(bps):
map = {
1: digital.constellation_bpsk(),
2: digital.constellation_qpsk(),
3: digital.constellation_8psk()
}
return map[bps]
def __init__(self):
gr.top_block.__init__(self, "Gotenna Rx Hackrf")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 500000
self.fsk_deviation_hz = fsk_deviation_hz = 12500
self.chan_spacing = chan_spacing = 500000
self.baud_rate = baud_rate = 24000
##################################################
# Blocks
##################################################
self.rational_resampler_xxx_0 = filter.rational_resampler_fff(
interpolation=1,
decimation=4,
taps=None,
fractional_bw=None)
self.osmosdr_source_0 = osmosdr.source(
args="numchan=" + str(1) + " " + ''
)
self.osmosdr_source_0.set_time_unknown_pps(osmosdr.time_spec_t())
self.osmosdr_source_0.set_sample_rate(samp_rate)
self.osmosdr_source_0.set_center_freq(915000000, 0)
self.osmosdr_source_0.set_freq_corr(0, 0)
self.osmosdr_source_0.set_gain(1, 0)
self.osmosdr_source_0.set_if_gain(16, 0)
self.osmosdr_source_0.set_bb_gain(1, 0)
self.osmosdr_source_0.set_antenna('', 0)
self.osmosdr_source_0.set_bandwidth(26000000, 0)
self.gotenna_sink = gotenna_sink.blk()
self.digital_symbol_sync_xx_0 = digital.symbol_sync_ff(
digital.TED_DANDREA_AND_MENGALI_GEN_MSK,
float(chan_spacing) / baud_rate / 4,
0.05,
1.5,
1.0,
0.001 * float(chan_spacing) / baud_rate / 4,
1,
digital.constellation_bpsk().base(),
digital.IR_MMSE_8TAP,
128,
[])
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.dc_blocker_xx_0 = filter.dc_blocker_ff(512, True)
self.analog_quadrature_demod_cf_0 = analog.quadrature_demod_cf(chan_spacing/(2*math.pi*fsk_deviation_hz))
##################################################
# Connections
##################################################
self.connect((self.analog_quadrature_demod_cf_0, 0), (self.dc_blocker_xx_0, 0))
self.connect((self.dc_blocker_xx_0, 0), (self.rational_resampler_xxx_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.gotenna_sink, 0))
self.connect((self.digital_symbol_sync_xx_0, 0), (self.digital_binary_slicer_fb_0, 0))
self.connect((self.osmosdr_source_0, 0), (self.analog_quadrature_demod_cf_0, 0))
self.connect((self.rational_resampler_xxx_0, 0), (self.digital_symbol_sync_xx_0, 0))
def test_002_t(self):
"""Encoder to Decoder - noisy BPSK w/ soft constellation de-mapper"""
# Parameters
N = 18444
K = 6144
pct_en = False
n_ite = 6
M = 2
snr_db = SNR_DB
max_len = 10 * K
flip_llrs = False
sym_rate = 2 * 156250
# NOTE: The soft decoder block outputs positive LLR if bit "1" is more
# likely, and negative LLR otherwise (bit = 0). That is, the soft
# de-mapper's convention is the oposite of the one adopted in the Turbo
# Decoder. To compensate that, we choose to flip the LLRs (multiply them
# by "-1") inside the decoder wrapper.
# 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)
const = digital.constellation_bpsk().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, 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 test_001_t(self):
# set up variables
sps = 4
samp_rate_rx = 3000000
divider = 2
delay = 6
v_max = 10
skip_data = 10
avg_length = 20
bpsk = digital.constellation_bpsk().base()
# generate blocks
lab_radar_simple_decimator = lab_radar.simple_decimator_cc(
decimation=divider)
lab_radar_signal_corr_estimator_cf_0 = lab_radar.signal_corr_estimator_cf(
samp_rate_rx, divider, 8, sps, avg_length, skip_data, v_max)
digital_constellation_modulator_0 = digital.generic_mod(
constellation=bpsk,
differential=False,
samples_per_symbol=sps * divider,
pre_diff_code=True,
excess_bw=0.35,
verbose=False,
log=False)
blocks_vector_source_x_0 = blocks.vector_source_b(
(154, 154, 154, 154, 154, 154, 154, 154, 154, 154, 154, 154, 154,
154, 154, 154), True, 1, [])
offsets = blocks.vector_sink_f(1, 1024)
peaks = blocks.vector_sink_f(1, 1024)
distances = blocks.vector_sink_f(1, 1024)
blocks_throttle_0 = blocks.throttle(gr.sizeof_char * 1,
samp_rate_rx / sps / 8, True)
blocks_delay_0 = blocks.delay(gr.sizeof_gr_complex * 1, delay)
# setup connections
self.tb.connect((blocks_delay_0, 0), (lab_radar_simple_decimator, 0))
self.tb.connect((blocks_throttle_0, 0),
(digital_constellation_modulator_0, 0))
self.tb.connect((blocks_vector_source_x_0, 0), (blocks_throttle_0, 0))
self.tb.connect((digital_constellation_modulator_0, 0),
(blocks_delay_0, 0))
self.tb.connect((digital_constellation_modulator_0, 0),
(lab_radar_signal_corr_estimator_cf_0, 1))
self.tb.connect((lab_radar_signal_corr_estimator_cf_0, 0),
(distances, 0))
self.tb.connect((lab_radar_signal_corr_estimator_cf_0, 1), (peaks, 0))
self.tb.connect((lab_radar_signal_corr_estimator_cf_0, 2),
(offsets, 0))
self.tb.connect((lab_radar_simple_decimator, 0),
(lab_radar_signal_corr_estimator_cf_0, 0))
# set up fg
self.tb.run()
print(distances.data())
print(offsets.data())
print(peaks.data())
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_t(self):
"""MER of noisy QPSK symbols"""
# Parameters
snr_db = 10.0
alpha = 0.01
M = 2
frame_len = 10000
# Constants
bits_per_sym = int(log(M, 2))
n_bits = int(bits_per_sym * frame_len)
rndm = random.Random()
# Iterate over SNR levels
for snr_db in reversed(range(3, 15)):
print("Trying SNR: %f dB" % (float(snr_db)))
# Random input data
in_vec = tuple([rndm.randint(0, 1) for i in range(0, n_bits)])
src = blocks.vector_source_b(in_vec)
# Random noise
noise_var = 1.0 / (10**(float(snr_db) / 10))
noise_std = sqrt(noise_var)
print("Noise amplitude: %f" % (noise_std))
noise = analog.noise_source_c(analog.GR_GAUSSIAN, noise_std, 0)
# Fixed flowgraph blocks
pack = blocks.repack_bits_bb(1, bits_per_sym, "", False,
gr.GR_MSB_FIRST)
const = digital.constellation_bpsk().base()
cmap = digital.chunks_to_symbols_bc(const.points())
nadder = blocks.add_cc()
mer = blocksat.mer_measurement(alpha, M, frame_len)
snk = blocks.vector_sink_f()
snk_n = blocks.vector_sink_c()
# Connect source into flowgraph
self.tb.connect(src, pack, cmap)
self.tb.connect(cmap, (nadder, 0))
self.tb.connect(noise, (nadder, 1))
self.tb.connect(nadder, mer, snk)
self.tb.connect(noise, snk_n)
self.tb.run()
# Collect results
mer_measurements = snk.data()
noise_samples = snk_n.data()
print(mer_measurements)
avg_mer = np.mean(mer_measurements)
# Check results
self.assertAlmostEqual(round(avg_mer), snr_db, places=1)
def __init__(self):
constellation_source.__init__(self,
mod_name="bpsk",
samp_per_sym=2,
excess_bw=.35)
self.const = digital.constellation_bpsk().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.rrc_filter,
self)
def __init__(self,
length_tag_key="package_size",
mod_encabezado=digital.constellation_bpsk(),
occupied_carriers=(range(-27, -21) + range(-20, -7) +
range(-6, 0) + range(1, 7) + range(8, 21) +
range(22, 28), ),
package_size=96):
gr.hier_block2.__init__(
self,
"HeaderModulator",
gr.io_signature(1, 1, gr.sizeof_char * 1),
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
)
##################################################
# Parameters
##################################################
self.length_tag_key = length_tag_key
self.mod_encabezado = mod_encabezado
self.occupied_carriers = occupied_carriers
self.package_size = package_size
##################################################
# Variables
##################################################
self.hdr_format = hdr_format = digital.header_format_ofdm(
occupied_carriers,
1,
length_tag_key,
)
##################################################
# Blocks
##################################################
self.digital_protocol_formatter_bb_0 = digital.protocol_formatter_bb(
hdr_format, length_tag_key)
self.digital_chunks_to_symbols_xx_0 = digital.chunks_to_symbols_bc(
(mod_encabezado.points()), 1)
self.blocks_stream_to_tagged_stream_0 = blocks.stream_to_tagged_stream(
gr.sizeof_char, 1, package_size, length_tag_key)
self.blocks_repack_bits_bb_0_0 = blocks.repack_bits_bb(
8, 1, length_tag_key, False, gr.GR_LSB_FIRST)
##################################################
# Connections
##################################################
self.connect((self.blocks_repack_bits_bb_0_0, 0),
(self.digital_chunks_to_symbols_xx_0, 0))
self.connect((self.blocks_stream_to_tagged_stream_0, 0),
(self.digital_protocol_formatter_bb_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0, 0), (self, 0))
self.connect((self.digital_protocol_formatter_bb_0, 0),
(self.blocks_repack_bits_bb_0_0, 0))
self.connect((self, 0), (self.blocks_stream_to_tagged_stream_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=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 __init__(self):
ofdm_source.__init__(self, mod_name="ofdm_64_bpsk", fft_len=128)
self.const = digital.constellation_bpsk().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 _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 __init__(self, file_name="", modulation="bpsk"):
gr.top_block.__init__(self, "File Demodulator")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 100000
self.file_name = file_name
self.const = const = digital.constellation_bpsk().base()
self.const = digital.constellation_bpsk().base()
self.get_constellation_from_string(modulation)
##################################################
# Blocks
##################################################
self.digital_constellation_decoder_cb_0 = digital.constellation_decoder_cb(
self.const)
self.blocks_vector_sink_x_0 = blocks.vector_sink_b(1)
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_char * 1,
self.samp_rate, True)
self.blocks_repack_bits_bb_0 = blocks.repack_bits_bb(
self.const.bits_per_symbol(), 1, "", False, gr.GR_LSB_FIRST)
self.blocks_file_source_0 = blocks.file_source(
gr.sizeof_gr_complex * 1, self.file_name, False)
self.blocks_file_source_0.set_begin_tag(pmt.PMT_NIL)
##################################################
# Connections
##################################################
self.connect((self.blocks_file_source_0, 0),
(self.digital_constellation_decoder_cb_0, 0))
self.connect((self.blocks_repack_bits_bb_0, 0),
(self.blocks_throttle_0, 0))
self.connect((self.blocks_throttle_0, 0),
(self.blocks_vector_sink_x_0, 0))
self.connect((self.digital_constellation_decoder_cb_0, 0),
(self.blocks_repack_bits_bb_0, 0))
def test_001_t (self):
packet_len_tag = "packet_length"
self.payload_constellation = digital.constellation_bpsk()
self.data = [[0,],]
self.create_frames = ofdm_create_frames_bc(
data=self.data,
constellation=digital.constellation_8psk())
self.data_sink = blocks.file_sink(gr.sizeof_gr_complex, 'testdata.dat')
self.tb.connect(self.create_frames, self.data_sink)
# set up fg
self.tb.run ()
def __init__(self, pkt_len, nb_pkt, EbN0dB, fsm):
gr.top_block.__init__(self,
"Transmitter, channel and metrics computation")
##################################################
# Variables
##################################################
self.pkt_len = pkt_len
Rc = 0.5
self.const = const = digital.constellation_bpsk().base()
var_c = 1
Nb = 1
Eb = var_c / (2.0 * Nb * Rc)
N0 = Eb * 10**(-EbN0dB / 10.0)
noisevar = 2 * N0
##################################################
# Blocks
##################################################
self.bits_src = blocks.vector_source_b(
map(int, numpy.random.randint(0, 2, nb_pkt * pkt_len)), False)
self.trellis_encoder = trellis.encoder_bb(trellis.fsm(fsm), 0, pkt_len)
self.unpack = blocks.unpack_k_bits_bb(2)
self.to_symbols = digital.chunks_to_symbols_bc((const.points()), 1)
self.noise_src = analog.noise_source_c(analog.GR_GAUSSIAN, noisevar, 0)
self.noise_adder = blocks.add_vcc(1)
self.metrics_computer = trellis.metrics_c(
4, 2, ([-1, -1, -1, 1, 1, -1, 1, 1]), digital.TRELLIS_EUCLIDEAN)
self.dst = blocks.vector_sink_f()
##################################################
# Connections
##################################################
self.connect((self.bits_src, 0), (self.trellis_encoder, 0))
self.connect((self.trellis_encoder, 0), (self.unpack, 0))
self.connect((self.unpack, 0), (self.to_symbols, 0))
self.connect((self.to_symbols, 0), (self.noise_adder, 0))
self.connect((self.noise_src, 0), (self.noise_adder, 1))
self.connect((self.noise_adder, 0), (self.metrics_computer, 0))
self.connect((self.metrics_computer, 0), (self.dst, 0))
def test_constellation_decoder_cb_bpsk(self):
cnst = digital.constellation_bpsk()
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 = (1, 1, 0, 0, 1, 0, 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 __init__(self):
gr.top_block.__init__(self, "Bpsk Tx")
##################################################
# Variables
##################################################
self.sync_word = sync_word = [complex(1, 0), complex(1, 0), complex(-1, 0), complex(-1, 0), complex(-1, 0), complex(-1, 0), complex(1, 0), complex(-1, 0), complex(1, 0), complex(1, 0), complex(-1, 0), complex(1, 0), complex(1, 0), complex(1, 0), complex(-1, 0), complex(-1, 0), complex(1, 0), complex(-1, 0), complex(1, 0), complex(-1, 0), complex(-1, 0), complex(1, 0), complex(1, 0), complex(-1, 0), complex(-1, 0), complex(1, 0), complex(1, 0), complex(-1, 0), complex(1, 0), complex(1, 0), complex(1, 0), complex(-1, 0), complex(1, 0), complex(-1, 0), complex(1, 0), complex(1, 0), complex(1, 0), complex(-1, 0), complex(-1, 0), complex(-1, 0), complex(-1, 0), complex(1, 0), complex(-1, 0), complex(-1, 0), complex(-1, 0), complex(-1, 0), complex(-1, 0), complex(1, 0), complex(1, 0), complex(-1, 0)]
self.samp_rate = samp_rate = 1000000
self.samp_per_symb = samp_per_symb = 8
self.payload_mod = payload_mod = digital.constellation_bpsk()
self.packet_len = packet_len = 96
self.length_tag_name = length_tag_name = "packet_len"
self.bandwidth_excess = bandwidth_excess = .9
##################################################
# Blocks
##################################################
self.root_raised_cosine_filter_0 = filter.interp_fir_filter_ccf(samp_per_symb, firdes.root_raised_cosine(
samp_per_symb, samp_per_symb, 1, bandwidth_excess, 11*samp_per_symb))
self.digital_chunks_to_symbols_xx_0 = digital.chunks_to_symbols_bc((payload_mod.points()), 1)
self.dctk_generic_tx_path_0 = dctk.generic_tx_path(
payload_packet_len=packet_len,
sync_packet_len=len(sync_word),
)
self.blocks_vector_source_x_0 = blocks.vector_source_c(sync_word, True, 1, [])
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_gr_complex*1, samp_rate)
self.blocks_file_sink_0_0_0 = blocks.file_sink(gr.sizeof_char*1, "/home/csmedina/WiComModules/gr-dctk/build/tx_bits.txt")
self.blocks_file_sink_0_0_0.set_unbuffered(False)
self.blocks_file_sink_0_0 = blocks.file_sink(gr.sizeof_gr_complex*1, "/home/csmedina/WiComModules/gr-dctk/build/tx_symbols.txt")
self.blocks_file_sink_0_0.set_unbuffered(False)
self.blocks_file_sink_0 = blocks.file_sink(gr.sizeof_gr_complex*1, "/home/csmedina/WiComModules/gr-dctk/build/tx_2usrp.txt")
self.blocks_file_sink_0.set_unbuffered(False)
self.analog_random_source_x_0 = blocks.vector_source_b(map(int, numpy.random.randint(0, 2, packet_len)), True)
##################################################
# Connections
##################################################
self.connect((self.blocks_vector_source_x_0, 0), (self.dctk_generic_tx_path_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0, 0), (self.dctk_generic_tx_path_0, 1))
self.connect((self.analog_random_source_x_0, 0), (self.digital_chunks_to_symbols_xx_0, 0))
self.connect((self.analog_random_source_x_0, 0), (self.blocks_file_sink_0_0_0, 0))
self.connect((self.root_raised_cosine_filter_0, 0), (self.blocks_file_sink_0, 0))
self.connect((self.dctk_generic_tx_path_0, 0), (self.blocks_throttle_0, 0))
self.connect((self.blocks_throttle_0, 0), (self.blocks_file_sink_0_0, 0))
self.connect((self.blocks_throttle_0, 0), (self.root_raised_cosine_filter_0, 0))
def __init__(self, samp_rate, symbol_rate=31.25):
super(psk31_receiver, self).__init__(
"psk31_receiver",
gr.io_signature(1, 1, gr.sizeof_gr_complex),
gr.io_signature(0, 0, 1))
self.symbol_rate = symbol_rate
self.costas = digital.costas_loop_cc(2*3.14/100, 4)
self.clock_recovery = digital.clock_recovery_mm_cc(
1.0*samp_rate/symbol_rate, 0.25 * 0.1*0.1, 0.05, 0.1, 0.001)
self.receiver = digital.constellation_receiver_cb(
digital.constellation_bpsk().base(), 2*3.14/100, -0.25, 0.25)
self.diff = gr.diff_decoder_bb(2)
self.decoder = ham.psk31_decode_bb(True)
self.msgq_out = gr.msg_queue()
self.snk = gr.message_sink(gr.sizeof_char, self.msgq_out, True)
self.connect(self, self.costas, self.clock_recovery,
self.receiver, self.diff, self.decoder, self.snk)
def test_constellation_decoder_cb_bpsk(self):
cnst = digital.constellation_bpsk()
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 = ( 1, 1, 0, 0,
1, 0, 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 __init__(self, constellation_points=_def_constellation_points,
differential=False, *args, **kwargs):
"""
Hierarchical block for RRC-filtered BPSK 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_bpsk()
if constellation_points != 2:
raise ValueError('Number of constellation points must be 2 for BPSK.')
super(bpsk_demod, self).__init__(constellation=constellation,
differential=differential, *args, **kwargs)
def test_constellation_encoder_bc_bpsk(self):
cnst = digital.constellation_bpsk()
src_data = (1, 1, 0, 0, 1, 0, 1)
const_map = [-1.0, 1.0]
expected_result = [const_map[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_001_t(self):
packet_len_tag = "packet_length"
self.payload_constellation = digital.constellation_bpsk()
self.data = [
[
0,
],
]
self.create_frames = ofdm_create_frames_bc(
data=self.data, constellation=digital.constellation_8psk())
self.data_sink = blocks.file_sink(gr.sizeof_gr_complex, 'testdata.dat')
self.tb.connect(self.create_frames, self.data_sink)
# set up fg
self.tb.run()
def test_001_t (self):
chaos_values = (1+0j, 0+1j, -1+0j, 0-1j)
expected_result = (0, 1, 0, 0, 1, 1)
zero = chaos_values + tuple(map(lambda v: -v, chaos_values))
one = 2 * chaos_values
src_data = zero + one + zero + zero + one + one
src = blocks.vector_source_c (src_data)
demod = chaos.dcsk_demod_cc (len(chaos_values), 0)
const = digital.constellation_bpsk()
const_decoder = digital.constellation_decoder_cb(const.base())
dst = blocks.vector_sink_b ()
self.tb.connect (src, demod, const_decoder, dst)
self.tb.run ()
result_data = dst.data ()
self.assertEqual (expected_result, result_data)
def __init__(self):
gr.top_block.__init__(self, "Top Block")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 32000
self.bpsk1 = bpsk1 = digital.constellation_bpsk().base()
##################################################
# Blocks
##################################################
self.digital_pfb_clock_sync_xxx_0 = digital.pfb_clock_sync_ccf(
4, .0628,
(firdes.root_raised_cosine(32, 32, 1.0 / float(1), 0.35, 45 * 32)),
32, 16, 1.5, 1)
self.digital_constellation_decoder_cb_0 = digital.constellation_decoder_cb(
bpsk1)
self.blocks_wavfile_source_0 = blocks.wavfile_source(
os.getenv("DIR", "/data") + os.sep + 'challenge.wav', False)
self.blocks_float_to_complex_0 = blocks.float_to_complex(1)
self.blocks_file_sink_0 = blocks.file_sink(
gr.sizeof_char * 1,
os.getenv("DIR", "/data") + os.sep + 'output.txt', False)
self.blocks_file_sink_0.set_unbuffered(True)
##################################################
# Connections
##################################################
self.connect((self.blocks_wavfile_source_0, 0),
(self.blocks_float_to_complex_0, 0))
self.connect((self.blocks_float_to_complex_0, 0),
(self.digital_pfb_clock_sync_xxx_0, 0))
self.connect((self.digital_pfb_clock_sync_xxx_0, 0),
(self.digital_constellation_decoder_cb_0, 0))
self.connect((self.digital_constellation_decoder_cb_0, 0),
(self.blocks_file_sink_0, 0))
def __init__(self,
length_tag_key="package_size",
package_size=96,
payload_mod=digital.constellation_bpsk()):
gr.hier_block2.__init__(
self,
"PayloadModulator",
gr.io_signature(1, 1, gr.sizeof_char * 1),
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
)
##################################################
# Parameters
##################################################
self.length_tag_key = length_tag_key
self.package_size = package_size
self.payload_mod = payload_mod
##################################################
# Blocks
##################################################
self.digital_chunks_to_symbols_xx_0_0 = digital.chunks_to_symbols_bc(
(payload_mod.points()), 1)
self.blocks_stream_to_tagged_stream_0_0 = blocks.stream_to_tagged_stream(
gr.sizeof_char, 1, package_size, length_tag_key)
self.blocks_repack_bits_bb_0 = blocks.repack_bits_bb(
8, payload_mod.bits_per_symbol(), length_tag_key, False,
gr.GR_LSB_FIRST)
##################################################
# Connections
##################################################
self.connect((self.blocks_repack_bits_bb_0, 0),
(self.digital_chunks_to_symbols_xx_0_0, 0))
self.connect((self.blocks_stream_to_tagged_stream_0_0, 0),
(self.blocks_repack_bits_bb_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0_0, 0), (self, 0))
self.connect((self, 0), (self.blocks_stream_to_tagged_stream_0_0, 0))
def __init__(self,
constellation_points=_def_constellation_points,
differential=False,
*args,
**kwargs):
"""
Hierarchical block for RRC-filtered BPSK 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_bpsk()
if constellation_points != 2:
raise ValueError(
'Number of constellation points must be 2 for BPSK.')
super(bpsk_demod, self).__init__(constellation=constellation,
differential=differential,
*args,
**kwargs)
def __init__(self):
gr.top_block.__init__(self, "Top Block")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 32000
self.bpsk1 = bpsk1 = digital.constellation_bpsk().base()
##################################################
# Blocks
##################################################
self.digital_constellation_modulator_0 = digital.generic_mod(
constellation=bpsk1,
differential=False,
samples_per_symbol=4,
pre_diff_code=False,
excess_bw=0.35,
verbose=False,
log=False,
)
self.blocks_wavfile_sink_0 = blocks.wavfile_sink(
'/tmp/' + filename, 1, samp_rate, 16)
self.blocks_file_source_0 = blocks.file_source(gr.sizeof_char * 1,
'/tmp/input.txt', False)
self.blocks_complex_to_float_0 = blocks.complex_to_float(1)
##################################################
# Connections
##################################################
self.connect((self.blocks_complex_to_float_0, 0),
(self.blocks_wavfile_sink_0, 0))
self.connect((self.blocks_file_source_0, 0),
(self.digital_constellation_modulator_0, 0))
self.connect((self.digital_constellation_modulator_0, 0),
(self.blocks_complex_to_float_0, 0))
#
import math
import os
from gnuradio import gr, gr_unittest, trellis, digital, analog, blocks
fsm_args = {"awgn1o2_4": (2, 4, 4,
(0, 2, 0, 2, 1, 3, 1, 3),
(0, 3, 3, 0, 1, 2, 2, 1),
),
"rep2": (2, 1, 4, (0, 0), (0, 3)),
"nothing": (2, 1, 2, (0, 0), (0, 1)),
}
constells = {2: digital.constellation_bpsk(),
4: digital.constellation_qpsk(),
}
class test_trellis (gr_unittest.TestCase):
def test_001_fsm (self):
f = trellis.fsm(*fsm_args["awgn1o2_4"])
self.assertEqual(fsm_args["awgn1o2_4"],(f.I(),f.S(),f.O(),f.NS(),f.OS()))
def test_002_fsm (self):
f = trellis.fsm(*fsm_args["awgn1o2_4"])
g = trellis.fsm(f)
self.assertEqual((g.I(),g.S(),g.O(),g.NS(),g.OS()),(f.I(),f.S(),f.O(),f.NS(),f.OS()))
def test_003_fsm (self):
def main():
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),])
data_len = utils.ofdm_get_data_len(
nofdm_symbols=args.nsymbols,
noccupied_carriers=len(occupied_carriers[0]),
constellation=constellation[args.bits]
)
data = args.nframes*[[1 for x in xrange(data_len)],]
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
)
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
)
ofdm_basebander = mimoots.ofdm_frames_to_basebandsignal_vcc(
fft_len=fft_len,
cp_len=cp_len,
packet_len_tag=packet_len_tag
)
shifter = blocks.delay(
itemsize=gr.sizeof_gr_complex,
delay=(fft_len+cp_len)*(3)
)
xor_source = blocks.vector_source_c(
data=5*[0]+(3*(fft_len+cp_len)-10)*[1]+5*[0]+3*(fft_len+cp_len)*[0],
repeat=True
)
xor1 = blocks.multiply_cc()
xor2 = blocks.multiply_cc()
if args.freq == None:
data_sink = mimoots.file_sink2(
itemsize=(gr.sizeof_gr_complex, gr.sizeof_gr_complex),
filename=(
'1.'.join(args.to_file.rsplit('.',1)),
'2.'.join(args.to_file.rsplit('.',1))
)
)
else:
data_sink = mimoots.uhd_sink2(freq=args.freq, gain=args.gain)
tb.connect(xor_source, (xor1, 1))
tb.connect(xor_source, (xor2, 1))
tb.connect(data_source, ofdm_mapper, ofdm_framer, ofdm_basebander)
tb.connect(ofdm_basebander, (xor1, 0), (data_sink, 0))
tb.connect(ofdm_basebander, (xor2, 0), shifter, (data_sink, 1))
tb.run()
time.sleep(5)
def __init__(
self,
pilot_carriers=((-40, -14, 13, 39),),
pilot_symbols=((1, 1, 1, -1),),
occupied_carriers=(
[
-54,
-53,
-52,
-51,
-50,
-49,
-48,
-47,
-46,
-45,
-44,
-43,
-42,
-41,
-39,
-38,
-37,
-36,
-35,
-34,
-33,
-32,
-31,
-30,
-29,
-28,
-27,
-26,
-25,
-24,
-23,
-22,
-21,
-20,
-19,
-18,
-17,
-16,
-15,
-13,
-12,
-11,
-10,
-9,
-8,
-7,
-6,
-5,
-4,
-3,
-2,
-1,
1,
2,
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26,
27,
28,
29,
30,
31,
32,
33,
34,
35,
36,
37,
38,
40,
41,
42,
43,
44,
45,
46,
47,
48,
49,
50,
51,
52,
53,
],
),
samp_rate=10000,
payload_mod="qpsk",
sync_word1=[
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
-1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
1.42,
0.0,
1.42,
0.0,
1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
1.42,
0.0,
1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
-1.42,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
0.0,
],
sync_word2=[
0j,
0j,
0j,
0j,
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),
(-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),
(-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),
(-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),
(1 + 0j),
0j,
0j,
0j,
0j,
0j,
0j,
0j,
0j,
0j,
0j,
],
scramble_mode=0,
crc_mode=0,
clipper_mode=0,
filter_mode=1,
clipping_factor=10,
):
gr.hier_block2.__init__(
self,
"Ofdm Radio Hier",
gr.io_signaturev(2, 2, [gr.sizeof_char * 1, gr.sizeof_gr_complex * 1]),
gr.io_signaturev(2, 2, [gr.sizeof_char * 1, gr.sizeof_gr_complex * 1]),
)
##################################################
# Parameters
##################################################
self.pilot_carriers = pilot_carriers
self.pilot_symbols = pilot_symbols
self.occupied_carriers = occupied_carriers
self.samp_rate = samp_rate
self.sync_word1 = sync_word1
self.sync_word2 = sync_word2
self.scramble_mode = scramble_mode
self.crc_mode = crc_mode
self.clipping_factor = clipping_factor
self.clipper_mode = clipper_mode
self.filter_mode = filter_mode
if payload_mod == "qpsk":
self.payload_mod = payload_mod = digital.constellation_qpsk()
elif payload_mod == "qam16":
self.payload_mod = payload_mod = digital.qam.qam_constellation(16, True, "none", False)
elif payload_mod == "bpsk":
self.payload_mod = payload_mod = digital.constellation_bpsk()
##################################################
# Variables
##################################################
self.packet_length_tag_key = packet_length_tag_key = "packet_len"
self.length_tag_key = length_tag_key = "frame_len"
self.header_mod = header_mod = digital.constellation_bpsk()
self.fft_len = fft_len = (len(sync_word1) + len(sync_word2)) / 2
self.scramble_seed = scramble_seed = 0x7F
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.len_ocup_carr = len_ocup_carr = len(occupied_carriers[0])
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=True,
)
self.header_equalizer = header_equalizer = digital.ofdm_equalizer_simpledfe(
fft_len, header_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols
)
self.forward_OOB = forward_OOB = [
0.40789374966665903,
3.2351160543115207,
11.253435139165413,
22.423991613997735,
27.99555756436666,
22.423991613997735,
11.253435139165425,
3.235116054311531,
0.40789374966666014,
]
self.feedback_OOB = feedback_OOB = [
1.0,
6.170110168740749,
16.888669609673336,
26.73762881119027,
26.75444043101795,
17.322358010203928,
7.091659316015212,
1.682084643429639,
0.17795354282083842,
]
self.cp_len_0 = cp_len_0 = fft_len / 4
self.cp_len = cp_len = fft_len / 4
self.active_carriers = active_carriers = len(occupied_carriers[0]) + 4
##################################################
# Blocks
##################################################
self.ofdm_tools_clipper_0 = ofdm_tools.clipper_cc(clipping_factor)
self.iir_filter_xxx_1 = filter.iir_filter_ccd((forward_OOB), (feedback_OOB), False)
self.fft_vxx_txpath = fft.fft_vcc(fft_len, False, (()), True, 1)
self.fft_vxx_2_rxpath = fft.fft_vcc(fft_len, True, (), True, 1)
self.fft_vxx_1_rxpath = fft.fft_vcc(fft_len, True, (()), True, 1)
self.digital_packet_headerparser_b_rxpath = digital.packet_headerparser_b(header_formatter.base())
self.digital_packet_headergenerator_bb_txpath = digital.packet_headergenerator_bb(
header_formatter.formatter(), "packet_len"
)
self.digital_ofdm_sync_sc_cfb_rxpath = digital.ofdm_sync_sc_cfb(fft_len, fft_len / 4, False)
self.digital_ofdm_serializer_vcc_payload_rxpath = digital.ofdm_serializer_vcc(
fft_len, occupied_carriers, length_tag_key, packet_length_tag_key, 1, "", True
)
self.digital_ofdm_serializer_vcc_header_rxpath = digital.ofdm_serializer_vcc(
fft_len, occupied_carriers, length_tag_key, "", 0, "", True
)
self.digital_ofdm_frame_equalizer_vcvc_2_rxpath = digital.ofdm_frame_equalizer_vcvc(
payload_equalizer.base(), cp_len, length_tag_key, True, 0
)
self.digital_ofdm_frame_equalizer_vcvc_1_rxpath = digital.ofdm_frame_equalizer_vcvc(
header_equalizer.base(), cp_len, length_tag_key, True, 1
)
self.digital_ofdm_cyclic_prefixer_txpath = digital.ofdm_cyclic_prefixer(
fft_len, fft_len + cp_len, rolloff, packet_length_tag_key
)
(self.digital_ofdm_cyclic_prefixer_txpath).set_min_output_buffer(24000)
self.digital_ofdm_chanest_vcvc_rxpath = digital.ofdm_chanest_vcvc((sync_word1), (sync_word2), 1, 0, 3, False)
self.digital_ofdm_carrier_allocator_cvc_txpath = digital.ofdm_carrier_allocator_cvc(
fft_len, occupied_carriers, pilot_carriers, pilot_symbols, (sync_word1, sync_word2), packet_length_tag_key
)
(self.digital_ofdm_carrier_allocator_cvc_txpath).set_min_output_buffer(16000)
self.digital_header_payload_demux_rxpath = digital.header_payload_demux(
3, fft_len, cp_len, length_tag_key, "", True, gr.sizeof_gr_complex, "rx_time", samp_rate, ()
)
self.digital_crc32_bb_txpath = digital.crc32_bb(False, packet_length_tag_key)
self.digital_crc32_bb_rxpath = digital.crc32_bb(True, packet_length_tag_key)
self.digital_constellation_decoder_cb_1_rxpath = digital.constellation_decoder_cb(payload_mod.base())
self.digital_constellation_decoder_cb_0 = digital.constellation_decoder_cb(header_mod.base())
self.digital_chunks_to_symbols_x_txpath = digital.chunks_to_symbols_bc((payload_mod.points()), 1)
self.digital_chunks_to_symbols_txpath = digital.chunks_to_symbols_bc((header_mod.points()), 1)
self.digital_additive_scrambler_bb_txpath_0 = digital.additive_scrambler_bb(
0x8A, 0x7F, 7, 0, bits_per_byte=8, reset_tag_key=self.packet_length_tag_key
)
self.digital_additive_scrambler_bb_rxpath_0 = digital.additive_scrambler_bb(
0x8A, 0x7F, 7, 0, bits_per_byte=8, reset_tag_key=self.packet_length_tag_key
)
self.blocks_tagged_stream_mux_txpath = blocks.tagged_stream_mux(
gr.sizeof_gr_complex * 1, packet_length_tag_key, 0
)
(self.blocks_tagged_stream_mux_txpath).set_min_output_buffer(16000)
self.blocks_tag_gate_txpath = blocks.tag_gate(gr.sizeof_gr_complex * 1, False)
self.blocks_repack_bits_bb_txpath = blocks.repack_bits_bb(
8, payload_mod.bits_per_symbol(), packet_length_tag_key, False
)
self.blocks_repack_bits_bb_rxpath = blocks.repack_bits_bb(
payload_mod.bits_per_symbol(), 8, packet_length_tag_key, True
)
self.blocks_multiply_xx_rxpath = blocks.multiply_vcc(1)
self.blocks_multiply_const_vxx_0 = blocks.multiply_const_vcc((0.01,))
self.blocks_delay_rxpath = blocks.delay(gr.sizeof_gr_complex * 1, fft_len + fft_len / 4)
self.blks2_selector_0_2 = grc_blks2.selector(
item_size=gr.sizeof_gr_complex * 1, num_inputs=2, num_outputs=1, input_index=clipper_mode, output_index=0
)
self.blks2_selector_0_1 = grc_blks2.selector(
item_size=gr.sizeof_char * 1, num_inputs=2, num_outputs=1, input_index=scramble_mode, output_index=0
)
self.blks2_selector_0_0_0 = grc_blks2.selector(
item_size=gr.sizeof_char * 1, num_inputs=2, num_outputs=1, input_index=scramble_mode, output_index=0
)
self.blks2_selector_0_0 = grc_blks2.selector(
item_size=gr.sizeof_char * 1, num_inputs=2, num_outputs=1, input_index=crc_mode, output_index=0
)
self.blks2_selector_0 = grc_blks2.selector(
item_size=gr.sizeof_char * 1, num_inputs=2, num_outputs=1, input_index=crc_mode, output_index=0
)
self.blks2_selector_0_2_0 = grc_blks2.selector(
item_size=gr.sizeof_gr_complex * 1, num_inputs=2, num_outputs=1, input_index=filter_mode, output_index=0
)
self.analog_frequency_modulator_fc_rxpath = analog.frequency_modulator_fc(-2.0 / fft_len)
self.analog_agc2_xx_0 = analog.agc2_cc(1e-1, 1e-2, 1.0, 1.0)
self.analog_agc2_xx_0.set_max_gain(65536)
##################################################
# Connections
##################################################
self.connect((self.digital_ofdm_sync_sc_cfb_rxpath, 0), (self.analog_frequency_modulator_fc_rxpath, 0))
self.connect((self.analog_agc2_xx_0, 0), (self.digital_ofdm_sync_sc_cfb_rxpath, 0))
self.connect((self.analog_agc2_xx_0, 0), (self.blocks_delay_rxpath, 0))
self.connect((self.digital_packet_headergenerator_bb_txpath, 0), (self.digital_chunks_to_symbols_txpath, 0))
self.connect((self.blocks_repack_bits_bb_txpath, 0), (self.digital_chunks_to_symbols_x_txpath, 0))
self.connect((self.digital_chunks_to_symbols_txpath, 0), (self.blocks_tagged_stream_mux_txpath, 0))
self.connect((self.digital_chunks_to_symbols_x_txpath, 0), (self.blocks_tagged_stream_mux_txpath, 1))
self.connect((self.digital_ofdm_cyclic_prefixer_txpath, 0), (self.blocks_tag_gate_txpath, 0))
self.connect((self.fft_vxx_txpath, 0), (self.digital_ofdm_cyclic_prefixer_txpath, 0))
self.connect((self.blocks_tagged_stream_mux_txpath, 0), (self.digital_ofdm_carrier_allocator_cvc_txpath, 0))
self.connect((self.digital_ofdm_carrier_allocator_cvc_txpath, 0), (self.fft_vxx_txpath, 0))
self.connect((self.digital_header_payload_demux_rxpath, 0), (self.fft_vxx_1_rxpath, 0))
self.connect((self.fft_vxx_1_rxpath, 0), (self.digital_ofdm_chanest_vcvc_rxpath, 0))
self.connect(
(self.digital_ofdm_frame_equalizer_vcvc_1_rxpath, 0), (self.digital_ofdm_serializer_vcc_header_rxpath, 0)
)
self.connect((self.digital_ofdm_chanest_vcvc_rxpath, 0), (self.digital_ofdm_frame_equalizer_vcvc_1_rxpath, 0))
self.connect((self.digital_header_payload_demux_rxpath, 1), (self.fft_vxx_2_rxpath, 0))
self.connect(
(self.digital_ofdm_frame_equalizer_vcvc_2_rxpath, 0), (self.digital_ofdm_serializer_vcc_payload_rxpath, 0)
)
self.connect((self.fft_vxx_2_rxpath, 0), (self.digital_ofdm_frame_equalizer_vcvc_2_rxpath, 0))
self.connect(
(self.digital_ofdm_serializer_vcc_payload_rxpath, 0), (self.digital_constellation_decoder_cb_1_rxpath, 0)
)
self.connect((self.digital_constellation_decoder_cb_1_rxpath, 0), (self.blocks_repack_bits_bb_rxpath, 0))
self.connect((self.digital_ofdm_serializer_vcc_header_rxpath, 0), (self.digital_constellation_decoder_cb_0, 0))
self.connect((self.analog_frequency_modulator_fc_rxpath, 0), (self.blocks_multiply_xx_rxpath, 0))
self.connect((self.digital_ofdm_sync_sc_cfb_rxpath, 1), (self.digital_header_payload_demux_rxpath, 1))
self.connect((self.blocks_multiply_xx_rxpath, 0), (self.digital_header_payload_demux_rxpath, 0))
self.connect((self.blocks_delay_rxpath, 0), (self.blocks_multiply_xx_rxpath, 1))
self.connect((self.digital_constellation_decoder_cb_0, 0), (self.digital_packet_headerparser_b_rxpath, 0))
self.connect((self, 0), (self.digital_crc32_bb_txpath, 0))
self.connect((self, 1), (self.analog_agc2_xx_0, 0))
self.connect((self.blocks_tag_gate_txpath, 0), (self.blocks_multiply_const_vxx_0, 0))
self.connect((self, 0), (self.blks2_selector_0, 0))
self.connect((self.digital_crc32_bb_txpath, 0), (self.blks2_selector_0, 1))
self.connect((self.blks2_selector_0_1, 0), (self.blocks_repack_bits_bb_txpath, 0))
self.connect((self.blks2_selector_0, 0), (self.digital_packet_headergenerator_bb_txpath, 0))
self.connect((self.blks2_selector_0, 0), (self.digital_additive_scrambler_bb_txpath_0, 0))
self.connect((self.digital_additive_scrambler_bb_txpath_0, 0), (self.blks2_selector_0_1, 1))
self.connect((self.blks2_selector_0, 0), (self.blks2_selector_0_1, 0))
self.connect((self.digital_crc32_bb_rxpath, 0), (self.blks2_selector_0_0, 1))
self.connect((self.digital_additive_scrambler_bb_rxpath_0, 0), (self.blks2_selector_0_0_0, 1))
self.connect((self.blocks_repack_bits_bb_rxpath, 0), (self.digital_additive_scrambler_bb_rxpath_0, 0))
self.connect((self.blocks_repack_bits_bb_rxpath, 0), (self.blks2_selector_0_0_0, 0))
self.connect((self.blks2_selector_0_0_0, 0), (self.digital_crc32_bb_rxpath, 0))
self.connect((self.blks2_selector_0_0_0, 0), (self.blks2_selector_0_0, 0))
self.connect((self.blks2_selector_0_0, 0), (self, 0))
self.connect((self.blocks_multiply_const_vxx_0, 0), (self.ofdm_tools_clipper_0, 0))
self.connect((self.blocks_multiply_const_vxx_0, 0), (self.blks2_selector_0_2, 0))
self.connect(self.blks2_selector_0_2, self.iir_filter_xxx_1)
self.connect(self.blks2_selector_0_2, (self.blks2_selector_0_2_0, 0))
self.connect(self.iir_filter_xxx_1, (self.blks2_selector_0_2_0, 1))
self.connect(self.blks2_selector_0_2_0, (self, 1))
self.connect((self.ofdm_tools_clipper_0, 0), (self.blks2_selector_0_2, 1))
##################################################
# Asynch Message Connections
##################################################
self.msg_connect(
self.digital_packet_headerparser_b_rxpath,
"header_data",
self.digital_header_payload_demux_rxpath,
"header_data",
)
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))
def main():
args = get_arguments()
constellation = {
1:digital.constellation_bpsk(),
2:digital.constellation_qpsk(),
3:digital.constellation_8psk(),
}
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),])
packet_len_tag = "packet_length"
fft_len = 64
cp_len = 16
tb = gr.top_block()
if args.freq == None:
data_source = mimoots.file_source2(
itemsize=(gr.sizeof_gr_complex, gr.sizeof_gr_complex),
filename=(
'1.'.join(args.from_file.rsplit('.',1)),
'2.'.join(args.from_file.rsplit('.',1))
)
)
else:
data_source = mimoots.uhd_source(freq=args.freq, gain=args.gain)
skip0 = blocks.skiphead(
itemsize=gr.sizeof_gr_complex,
nitems_to_skip=args.skiphead
)
skip1 = blocks.skiphead(
itemsize=gr.sizeof_gr_complex,
nitems_to_skip=args.skiphead
)
#ofdm_frames1 = mimoots.ofdm_receive_frames_cb(
# nofdm_frames=args.nframes,
# nofdm_symbols=args.nsymbols,
# constellation=constellation[args.bits],
# occupied_carriers=occupied_carriers,
# pilot_carriers=pilot_carriers,
# pilot_symbols=pilot_symbols,
# debug=args.debug
#)
#ofdm_frames2 = mimoots.ofdm_receive_frames_cb(
# nofdm_frames=args.nframes,
# nofdm_symbols=args.nsymbols,
# constellation=constellation[args.bits],
# occupied_carriers=occupied_carriers,
# pilot_carriers=pilot_carriers,
# pilot_symbols=pilot_symbols,
# debug=args.debug
#)
ofdm_framer0 = mimoots.ofdm_basebandsignal_to_frames_cvc(
fft_len=fft_len,
cp_len=cp_len,
nofdm_symbols=args.nsymbols
)
ofdm_framer1 = mimoots.ofdm_basebandsignal_to_frames_cvc(
fft_len=fft_len,
cp_len=cp_len,
nofdm_symbols=args.nsymbols
)
ofdm_symboler0 = mimoots.ofdm_frame_to_symbols_vcc(
fft_len=fft_len,
cp_len=cp_len,
occupied_carriers=occupied_carriers,
pilot_carriers=pilot_carriers,
pilot_symbols=pilot_symbols,
constellation=constellation[args.bits],
nofdm_symbols=args.nsymbols,
packet_len_tag=packet_len_tag
)
ofdm_symboler1 = mimoots.ofdm_frame_to_symbols_vcc(
fft_len=fft_len,
cp_len=cp_len,
occupied_carriers=occupied_carriers,
pilot_carriers=pilot_carriers,
pilot_symbols=pilot_symbols,
constellation=constellation[args.bits],
nofdm_symbols=args.nsymbols,
packet_len_tag=packet_len_tag
)
ofdm_demapper0 = mimoots.ofdm_symbol_demapper_cb(
constellation=constellation[args.bits],
packet_len_tag=packet_len_tag
)
ofdm_demapper1 = mimoots.ofdm_symbol_demapper_cb(
constellation=constellation[args.bits],
packet_len_tag=packet_len_tag
)
data_sink0 = blocks.vector_sink_b()
data_sink1 = blocks.vector_sink_b()
tb.connect((data_source, 0), skip0, ofdm_framer0, ofdm_symboler0,
ofdm_demapper0, data_sink0)
tb.connect((data_source, 1), skip1, ofdm_framer1, ofdm_symboler1,
ofdm_demapper1, data_sink1)
tb.run()
random.seed(42)
#data_expected = tuple(args.nframes*[random.randint(0,255) for x in xrange(args.bits*60)])
data_len = utils.ofdm_get_data_len(
nofdm_symbols=args.nsymbols,
noccupied_carriers=len(occupied_carriers[0]),
constellation=constellation[args.bits]
)
if args.dummy_frame == True:
data_expected = [
tuple(data_len*[0] + args.nframes*[1 for x in xrange(data_len)]),
tuple(data_len*[0] + args.nframes*[2 for x in xrange(data_len)])
]
else:
data_expected = [
tuple(args.nframes*[1 for x in xrange(data_len)]),
tuple(args.nframes*[2 for x in xrange(data_len)])
]
data = [data_sink0.data(), data_sink1.data()]
len_data = []
len_data_expected = []
ignored_correlations = []
biterrors_complete = []
biterrorrate_complete = []
for data_index in range(2):
biterror = 0
index = 0
for i,j in itertools.izip_longest(data[data_index], data_expected[data_index]):
index += 1
if i != j:
if i is None or j is None:
biterror += 8
else:
#print("Index: {} data: {} expected: {}".format(index, i, j))
biterror += bin(i^j).count("1")
bits = args.bits
len_data = len(data[data_index])
len_data_expected = len(data_expected[data_index])
ignored_correlations = (len(data_expected[data_index])-len(data[data_index]))/(args.bits*60)
biterrors_complete = float(biterror)/float(8*len(data_expected[data_index]))
biterrorrate_complete = float(biterror)/float(8*len(data_expected[data_index]))
print '{}:len(data): {}'.format(data_index, len_data)
print '{}:len(data_expected): {}'.format(data_index, len_data_expected)
print '{}:Correlation didn\'t work: {} times'.format(data_index, ignored_correlations)
print '{}:biterrors: {} biterrorrate: {}'.format(data_index, biterror, biterrorrate_complete)
print "\n"
def main():
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),])
tb = gr.top_block()
if args.freq == None:
data_source = blocks.file_source(
itemsize=gr.sizeof_gr_complex,
filename=args.from_file
)
else:
data_source = mimoots.uhd_source(freq=args.freq, gain=args.gain)
skip = blocks.skiphead(
itemsize=gr.sizeof_gr_complex,
nitems_to_skip=args.skiphead
)
#ofdm_frames = mimoots.ofdm_receive_frames_cb(
# nofdm_frames=args.nframes,
# nofdm_symbols=args.nsymbols,
# constellation=constellation[args.bits],
# occupied_carriers=occupied_carriers,
# pilot_carriers=pilot_carriers,
# pilot_symbols=pilot_symbols,
# debug=args.debug
#)
ofdm_framer = mimoots.ofdm_basebandsignal_to_frames_cvc(
fft_len=fft_len,
cp_len=cp_len,
nofdm_symbols=args.nsymbols
)
ofdm_symboler = mimoots.ofdm_frame_to_symbols_vcc(
fft_len=fft_len,
cp_len=cp_len,
occupied_carriers=occupied_carriers,
pilot_carriers=pilot_carriers,
pilot_symbols=pilot_symbols,
constellation=constellation[args.bits],
nofdm_symbols=args.nsymbols,
packet_len_tag=packet_len_tag
);
ofdm_demapper = mimoots.ofdm_symbol_demapper_cb(
constellation=constellation[args.bits],
packet_len_tag=packet_len_tag
);
data_sink = blocks.vector_sink_b()
tb.connect(data_source, skip, ofdm_framer, ofdm_symboler, ofdm_demapper, data_sink)
tb.run()
random.seed(42)
#data_expected = tuple(args.nframes*[random.randint(0,255) for x in xrange(args.bits*60)])
data_len = utils.ofdm_get_data_len(
nofdm_symbols=args.nsymbols,
noccupied_carriers=len(occupied_carriers[0]),
constellation=constellation[args.bits]
)
if args.dummy_frame == True:
data_expected = tuple(data_len*[0] + args.nframes*[1 for x in xrange(data_len)])
else:
data_expected = tuple(args.nframes*[1 for x in xrange(data_len)])
data = data_sink.data()
# TODO: args.bits*60=datalen, 60 depends on len of data_carriers
len_data = []
len_data_expected = []
ignored_correlations = []
biterrors_complete = []
biterrorrate_complete = []
biterror = 0
index = 0
for i,j in itertools.izip_longest(data, data_expected):
index += 1
if i != j:
if i is None or j is None:
biterror += 8
else:
#print("Index: {} data: {} expected: {}".format(index, i, j))
biterror += bin(i^j).count("1")
bits = args.bits
len_data = len(data)
len_data_expected = len(data_expected)
ignored_correlations = (len(data_expected)-len(data))/(args.bits*60)
biterrors_complete = float(biterror)/float(8*len(data_expected))
biterrorrate_complete = float(biterror)/float(8*len(data_expected))
print 'len(data): {}'.format(len_data)
print 'len(data_expected): {}'.format(len_data_expected)
print 'Correlation didn\'t work: {} times'.format(ignored_correlations)
print 'biterrors: {} biterrorrate: {}'.format(biterror, biterrorrate_complete)
def bpsk_constellation(m=_def_constellation_points):
if m != _def_constellation_points:
raise ValueError("BPSK can only have 2 constellation points.")
return digital_swig.constellation_bpsk()
"""
print "CDMA PARAMETERS : for adaptive modulation"
prefix = "/home/anastas/gr-cdma" # put the prefix of your gr-cdma trunk
length_tag_name = "packet_len"
num_tag_name = "packet_num"
# header info
bits_per_header = 12 + 12 + 8 + 4
#Zhe Changed 12+16+8 to 12+12+8 because only 12 bits not 16 bits are needed. 4 bits indicating modulation and code mode.
header_mod = digital.constellation_bpsk()
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 = 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 "CDMA PARAMETERS : for adaptive modulation"
prefix = "/home/anastas/gr-cdma" # put the prefix of your gr-cdma trunk
length_tag_name = "packet_len"
num_tag_name = "packet_num"
# header info
bits_per_header = 12 + 12 + 8 + 4
#Zhe Changed 12+16+8 to 12+12+8 because only 12 bits not 16 bits are needed. 4 bits indicating modulation and code mode.
header_mod = digital.constellation_bpsk()
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 = 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
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, "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,
fft_len = 64,
cp_len = 16,
nofdm_symbols = 10,
nofdm_frames = 1,
ofdm_symbol_scale = 1,
constellation = digital.constellation_bpsk(),
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),]),
seq_seed = 42,
debug = False
):
gr.hier_block2.__init__(self,
"ofdm_receive_frames_cb",
gr.io_signature(1, 1, gr.sizeof_gr_complex), # Input signature
gr.io_signature(1, 1, gr.sizeof_char)) # Output signature
# =====================================================================
# Generate class-members
# =====================================================================
self._def_occupied_carriers = occupied_carriers
self._def_pilot_carriers = pilot_carriers
self._def_pilot_symbols = pilot_symbols
self._seq_seed = seq_seed
self.fft_len = fft_len
self.cp_len = cp_len
self.ofdm_symbol_scale = ofdm_symbol_scale
self.constellation = constellation
self.packet_len_tag = "packet_length"
self.frame_len_tag_key = "frame_length"
self.nofdm_symbols = nofdm_symbols
self.nofdm_frames = nofdm_frames
self.debug = debug
# =====================================================================
# Create all blocks
# =====================================================================
# TODO: rename some blocks
sync_detect = digital.ofdm_sync_sc_cfb(
fft_len = fft_len,
cp_len = cp_len
)
delay = blocks.delay(gr.sizeof_gr_complex, self.fft_len+self.cp_len)
oscillator = analog.frequency_modulator_fc(-2.0 / self.fft_len)
mixer = blocks.multiply_cc()
frames = mimoots.ofdm_extract_frame_cvc(
fft_len = self.fft_len,
cp_len = self.cp_len,
nsymbols_per_ofdmframe = self.nofdm_symbols+2 # +2 Sync-Words
)
fft_payload = fft.fft_vcc(
fft_size = self.fft_len,
forward = True,
window = (),
shift = True
)
chanest = digital.ofdm_chanest_vcvc(
sync_symbol1 = utils.ofdm_make_sync_word1(self.fft_len,
self._def_occupied_carriers,
self._def_pilot_carriers),
sync_symbol2 = utils.ofdm_make_sync_word2(self.fft_len,
self._def_occupied_carriers,
self._def_pilot_carriers),
n_data_symbols = self.nofdm_symbols
)
payload_equalizer = digital.ofdm_equalizer_simpledfe(
fft_len = self.fft_len,
constellation = self.constellation.base(),
occupied_carriers = self._def_occupied_carriers,
pilot_carriers = self._def_pilot_carriers,
pilot_symbols = self._def_pilot_symbols,
symbols_skipped = 0,
)
payload_eq = digital.ofdm_frame_equalizer_vcvc(
equalizer = payload_equalizer.base(),
cp_len = cp_len,
len_tag_key = self.frame_len_tag_key,
propagate_channel_state = True,
fixed_frame_len = self.nofdm_symbols
)
# doesn't accept names of parameters
payload_serializer = digital.ofdm_serializer_vcc(
self.fft_len, # fft_len =
self._def_occupied_carriers, # occupied_carriers =
self.frame_len_tag_key, # len_tag_key =
self.packet_len_tag, # packet_len_tag =
0 # symbolsskipped =
)
payload_demod = digital.constellation_decoder_cb(
constellation = self.constellation.base()
)
payload_pack = blocks.repack_bits_bb(
k = self.constellation.bits_per_symbol(),
l = 8,
len_tag_key = self.packet_len_tag,
align_output = True
)
# =====================================================================
# Connect all blocks
# =====================================================================
# TODO: Clean up graph
self.connect(self, sync_detect)
self.connect((sync_detect,0), oscillator, (mixer,0))
self.connect((self,0), delay, (mixer,1))
self.connect((sync_detect,1), (frames,1))
self.connect(mixer, (frames,0), fft_payload, chanest,
payload_eq, payload_serializer, payload_demod,
payload_pack, self)
# =====================================================================
# Debug-Output
# =====================================================================
if self.debug == True:
self.connect(self,
blocks.file_sink(gr.sizeof_gr_complex,
'receive-self.dat'))
self.connect((sync_detect,0),
blocks.file_sink(gr.sizeof_float,
'receive-sync_detect-0.dat'))
self.connect((sync_detect,1),
blocks.file_sink(gr.sizeof_char,
'receive-sync_detect-1.dat'))
self.connect(mixer,
blocks.file_sink(gr.sizeof_gr_complex,
'receive-mixer.dat'))
self.connect(frames,
blocks.file_sink(fft_len*gr.sizeof_gr_complex,
'receive-frames.dat'))
self.connect(fft_payload,
blocks.file_sink(self.fft_len*gr.sizeof_gr_complex,
'receive-fft_payload.dat'))
self.connect(chanest,
blocks.file_sink(self.fft_len*gr.sizeof_gr_complex,
'receive-chanest.dat'))
self.connect(payload_eq,
blocks.file_sink(self.fft_len*gr.sizeof_gr_complex,
'receive-payload_eq.dat'))
self.connect(payload_serializer,
blocks.file_sink(gr.sizeof_char*8,
'receive-payload_serializer.dat'))
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):
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, 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))
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, channel_noise, img_path):
gr.top_block.__init__(self, "Top Block")
##################################################
# Variables
##################################################
self.taps = taps = [1]
self.size_packed_file = size_packed_file = os.path.getsize(
"/home/rcampello/Main/FixedPath/OOT Gnuradio/Material/GnuRadio/Scripts/BPSK/{}"
.format(str(img_path)))
self.samp_rate = samp_rate = 32000 * 20
self.rrc_taps = rrc_taps = firdes.root_raised_cosine(
32, 32, 1.0 / float(2), 0.35, 11 * 2 * 32)
self.noise = noise = channel_noise
self.eq_gain = eq_gain = 0
self.delay2 = delay2 = int(8)
self.delay = delay = int(2)
self.BPSK = BPSK = digital.constellation_bpsk().base()
##################################################
# Blocks
##################################################
self.digital_pfb_clock_sync_xxx_0 = digital.pfb_clock_sync_ccf(
4, 62.8e-3, (rrc_taps), 32, 16, 1.5, 2)
self.digital_costas_loop_cc_0 = digital.costas_loop_cc(
62.8e-3, 2, False)
self.digital_constellation_modulator_0 = digital.generic_mod(
constellation=BPSK,
differential=False,
samples_per_symbol=4,
pre_diff_code=True,
excess_bw=0.35,
verbose=False,
log=False,
)
self.digital_constellation_decoder_cb_0 = digital.constellation_decoder_cb(
BPSK)
self.digital_cma_equalizer_cc_0 = digital.cma_equalizer_cc(
15, 1, eq_gain, 2)
self.channels_channel_model_0 = channels.channel_model(
noise_voltage=noise,
frequency_offset=0.0,
epsilon=1.0,
taps=(taps),
noise_seed=0,
block_tags=False)
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_char * 1, samp_rate,
True)
self.blocks_skiphead_0_0 = blocks.skiphead(gr.sizeof_char * 1,
int(size_packed_file * 0.1))
self.blocks_pack_k_bits_bb_1 = blocks.pack_k_bits_bb(8)
self.blocks_null_sink_0 = blocks.null_sink(gr.sizeof_char * 1)
self.blocks_head_0_0 = blocks.head(gr.sizeof_char * 1,
int(size_packed_file * 2))
self.blocks_file_source_0 = blocks.file_source(
gr.sizeof_char * 1,
'/home/rcampello/Main/FixedPath/OOT Gnuradio/Material/GnuRadio/Scripts/BPSK/{}'
.format(str(img_path)), True)
self.sink = blocks.vector_sink_b()
##################################################
# Connections
##################################################
self.connect((self.blocks_file_source_0, 0),
(self.blocks_throttle_0, 0))
self.connect((self.blocks_head_0_0, 0),
(self.digital_constellation_modulator_0, 0))
self.connect((self.blocks_pack_k_bits_bb_1, 0), (self.sink, 0))
self.connect((self.blocks_pack_k_bits_bb_1, 0),
(self.blocks_null_sink_0, 0))
self.connect((self.blocks_skiphead_0_0, 0), (self.blocks_head_0_0, 0))
self.connect((self.blocks_throttle_0, 0),
(self.blocks_skiphead_0_0, 0))
self.connect((self.channels_channel_model_0, 0),
(self.digital_pfb_clock_sync_xxx_0, 0))
self.connect((self.digital_cma_equalizer_cc_0, 0),
(self.digital_costas_loop_cc_0, 0))
self.connect((self.digital_constellation_decoder_cb_0, 0),
(self.blocks_pack_k_bits_bb_1, 0))
self.connect((self.digital_constellation_modulator_0, 0),
(self.channels_channel_model_0, 0))
self.connect((self.digital_costas_loop_cc_0, 0),
(self.digital_constellation_decoder_cb_0, 0))
self.connect((self.digital_pfb_clock_sync_xxx_0, 0),
(self.digital_cma_equalizer_cc_0, 0))
def __init__(self,
data=[
[],
],
fft_len=64,
cp_len=16,
nofdm_symbols=10,
nofdm_frames=1,
ofdm_symbol_scale=1,
constellation=digital.constellation_bpsk(),
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),
]),
scale=1.0,
seq_seed=42,
debug=False):
gr.hier_block2.__init__(
self,
"ofdm_create_frames_bc",
gr.io_signature(0, 0, 0), # Input signature
gr.io_signature(1, 1, gr.sizeof_gr_complex)) # Output signature
# =====================================================================
# Generate class-members
# =====================================================================
self._def_occupied_carriers = occupied_carriers
self._def_pilot_carriers = pilot_carriers
self._def_pilot_symbols = pilot_symbols
self._seq_seed = seq_seed
self.data = data
self.fft_len = fft_len
self.cp_len = cp_len
self.ofdm_symbol_scale = ofdm_symbol_scale
self.constellation = constellation
self.packet_len_tag = "packet_length"
self.frame_length_tag_key = "frame_length"
self.nofdm_symbols = nofdm_symbols
self.nofdm_frames = nofdm_frames
self.scale = scale
self.debug = debug
# =====================================================================
# Create data to convert into OFDM-Frames
# =====================================================================
random.seed(self._seq_seed)
self.data_len = utils.ofdm_get_data_len(
nofdm_symbols=self.nofdm_symbols,
noccupied_carriers=len(self._def_occupied_carriers[0]),
constellation=constellation)
if debug == True:
print 'Frames: {}'.format(self.nofdm_frames)
print 'Länge: {}'.format(self.data_len)
(data_tosend,
tags) = packet_utils.packets_to_vectors(self.data,
self.packet_len_tag)
# ===================================================================
# Create all blocks
# ===================================================================
self.data_source = blocks.vector_source_b(data=data_tosend,
vlen=1,
tags=tags,
repeat=False)
self.payload_unpack = blocks.repack_bits_bb(
k=8,
l=self.constellation.bits_per_symbol(),
len_tag_key=self.packet_len_tag)
self.payload_mod = digital.chunks_to_symbols_bc(
symbol_table=self.constellation.points())
self.allocator = digital.ofdm_carrier_allocator_cvc(
fft_len=self.fft_len,
occupied_carriers=self._def_occupied_carriers,
pilot_carriers=self._def_pilot_carriers,
pilot_symbols=self._def_pilot_symbols,
sync_words=[
utils.ofdm_make_sync_word1(self.fft_len,
self._def_occupied_carriers,
self._def_pilot_carriers),
utils.ofdm_make_sync_word2(self.fft_len,
self._def_occupied_carriers,
self._def_pilot_carriers),
],
len_tag_key=self.packet_len_tag)
self.ffter = fft.fft_vcc(fft_size=self.fft_len,
forward=False,
window=(),
shift=True)
self.scale = mimoots.ofdm_scale_symbol_vcvc(symbol_len=self.fft_len,
scale=self.scale)
self.cyclic_prefixer = digital.ofdm_cyclic_prefixer(
input_size=self.fft_len,
output_size=self.fft_len + cp_len,
rolloff_len=0,
len_tag_key=self.packet_len_tag)
# ===================================================================
# Connect all blocks
# ===================================================================
self.connect(self.data_source, self.payload_unpack, self.payload_mod,
self.allocator, self.ffter, self.scale,
self.cyclic_prefixer, self)
#self.connect(
# self.data_source, self.payload_mod, self.payload_unpack,
# self.allocator, self.ffter, self.scale, self.cyclic_prefixer,
# self
#)
# ===================================================================
# Debug-Output
# ===================================================================
if self.debug == True:
self.connect(
self.data_source,
blocks.file_sink(gr.sizeof_char, 'create-data_source.dat'))
self.connect(
self.payload_unpack,
blocks.file_sink(gr.sizeof_char, 'create-payload_unpack.dat'))
self.connect(
self.payload_mod,
blocks.file_sink(gr.sizeof_gr_complex,
'create-payload_mod.dat'))
self.connect(
self.allocator,
blocks.file_sink(self.fft_len * gr.sizeof_gr_complex,
'create-allocator.dat'))
self.connect(
self.ffter,
blocks.file_sink(self.fft_len * gr.sizeof_gr_complex,
'create-ffter.dat'))
self.connect(
self.scale,
blocks.file_sink(self.fft_len * gr.sizeof_gr_complex,
'create-scale.dat'))
self.connect(
self.cyclic_prefixer,
blocks.file_sink(gr.sizeof_gr_complex,
'create-cyclic_prefixer.dat'))
def __init__(self):
gr.top_block.__init__(self, "OFDM Transceiver")
self._lock = threading.RLock()
##################################################
# Variables
##################################################
self.sent_pkt = 0
self.rcv_pkt = 0
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 = -80
self.waterfall_max = waterfall_max = -20
self.tx_gain = tx_gain = 0.03
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 = 20e6
self.payload_equalizer = payload_equalizer = digital.ofdm_equalizer_simpledfe(fft_len, payload_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols, 1)
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.osmosdr_source_0 = osmosdr.source( args="numchan=" + str(1) + " " + "bladerf=1" )
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.osmosdr_sink_0 = osmosdr.sink( args="numchan=" + str(1) + " " + "bladerf=0" )
self.osmosdr_sink_0.set_sample_rate(samp_rate)
self.osmosdr_sink_0.set_center_freq(center_freq, 0)
self.osmosdr_sink_0.set_freq_corr(10, 0)
self.osmosdr_sink_0.set_gain(20, 0)
self.osmosdr_sink_0.set_if_gain(20, 0)
self.osmosdr_sink_0.set_bb_gain(0, 0)
self.osmosdr_sink_0.set_antenna("", 0)
self.osmosdr_sink_0.set_bandwidth(samp_rate, 0)
self.digital_ofdm_tx_0 = digital.ofdm_tx(
fft_len=64, 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=1,
rolloff=0,
debug_log=True,
scramble_bits=False
)
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(1024, False)
self.blocks_tagged_stream_to_pdu_0 = blocks.tagged_stream_to_pdu(blocks.byte_t, packet_length_tag_key)
self.blocks_pdu_to_tagged_stream_0 = blocks.pdu_to_tagged_stream(blocks.byte_t, packet_length_tag_key)
self.blocks_multiply_const_vxx_0 = blocks.multiply_const_vcc((tx_gain, ))
self.blocks_message_debug_0 = blocks.message_debug()
self.pdu_block = messaging.pdu_block(self)
##################################################
# Connections
##################################################
self.msg_connect((self.pdu_block, 'out_pdu'), (self.blocks_pdu_to_tagged_stream_0, 'pdus'))
self.msg_connect((self.blocks_tagged_stream_to_pdu_0, 'pdus'), (self.pdu_block, 'in_pdu'))
self.connect((self.blocks_multiply_const_vxx_0, 0), (self.osmosdr_sink_0, 0))
self.connect((self.blocks_pdu_to_tagged_stream_0, 0), (self.digital_ofdm_tx_0, 0))
self.connect((self.dc_blocker_xx_0, 0), (self.digital_ofdm_rx_0, 0))
self.connect((self.digital_ofdm_rx_0, 0), (self.blocks_tagged_stream_to_pdu_0, 0))
self.connect((self.digital_ofdm_tx_0, 0), (self.blocks_multiply_const_vxx_0, 0))
self.connect((self.osmosdr_source_0, 0), (self.dc_blocker_xx_0, 0))
def __init__(self, options):
grc_wxgui.top_block_gui.__init__(self, title="DSSDR")
self.initialized = False
self.stopped = False
self.options = copy.copy(options)
self._constellation = digital.constellation_bpsk()
self._excess_bw = options.excess_bw
self._phase_bw = options.phase_bw
self._freq_bw = options.freq_bw
self._timing_bw = options.timing_bw
self._if_freq = options.if_freq
self._timing_max_dev= 1.5
self._demod_class = digital.bpsk_demod # the demodulator_class we're using
self._chbw_factor = options.chbw_factor # channel filter bandwidth factor
self._samples_per_second = 2e6
self._nav_samples_per_second = 16e6
self._down_decim = 1
self._down_samples_per_second = self._scope_sample_rate = self._samples_per_second/self._down_decim
self._up_samples_per_second = 1e6
self._asm_threshold = 0
self._access_code = None
self._tm_packet_id = 4
self._timestamp_id = 5
self._down_bitrate = options.bitrate
self._up_bitrate = options.up_bitrate
self._up_samples_per_symbol = self._up_samples_per_second/self._up_bitrate
self._samples_per_symbol = self._samples_per_second/self._down_decim/self._down_bitrate
self._down_sub_freq = options.down_sub_freq
self._up_sub_freq = 25e3
self._tm_len = 8920
self._up_tm_len = 8920
self._coding_method = options.coding_method
self._up_coding_method = 'None'
self._up_subcarrier = 'Square'
self._rs_i = 1
self._ccsds_channel = 38
self._uhd_carrier_offset = 10e3
self._turn_div = 749
self._turn_mult = 880
self._modulation_index = 'pi/3'
self._up_modulation_index = 1.047
self._max_carrier_offset = 0.1
self._dssdr_mixer_freq = options.rf_freq
self._up_coding_rate = '1'
self._down_coding_rate = '1'
self._down_conv_en = "False"
self._down_randomizer_en = options.down_randomizer_en
self._down_manchester_en = options.down_manchester_en
self._up_conv_en = "False"
self._up_idle_sequence = "\\x55"
self._down_default_gain = 64
self._up_default_gain = 44
self._up_en = True
if self._access_code is None:
self._access_code = packet_utils.default_access_code
#Construct the lookup table for parameter-setting functions
self.param_setters = {
"DSSDR_CHANNEL": self.setChannel,
"DSSDR_LO_FREQ": self.setLOFreq,
"DSSDR_REF_FREQ": self.setRefFreq,
"DSSDR_TURN_MULT": self.setTurnMult,
"DSSDR_TURN_DIV": self.setTurnDiv,
"DSSDR_UP_GAIN": self.setUpGain,
"DSSDR_DOWN_GAIN": self.setDownGain,
"DSSDR_UP_BITRATE": self.setUpBitrate,
"DSSDR_DOWN_BITRATE": self.setDownBitrate,
"DSSDR_DOWN_SAMPLE_RATE": self.setSampRate,
"DSSDR_UP_SUB_FREQ": self.setUpSubFreq,
"DSSDR_DOWN_SUB_FREQ": self.setDownSubFreq,
"DSSDR_DOPPLER_REPORT": self.dopplerReport,
"DSSDR_PN_RANGE": self.rangePN,
"DSSDR_SEQUENTIAL_RANGE": self.rangeSequential,
"DSSDR_DOWN_CODING_METHOD": self.setDownCodingMethod,
"DSSDR_UP_CODING_METHOD": self.setUpCodingMethod,
"DSSDR_DOWN_TM_LEN": self.setTMLen,
"DSSDR_UP_TM_LEN": self.setUpTMLen,
"DSSDR_DOWN_MOD_IDX": self.setDownModulationIndex,
"DSSDR_UP_MOD_IDX": self.setUpModulationIndex,
"DSSDR_DOWN_CONV_EN": self.setDownConvEn,
"DSSDR_UP_CONV_EN": self.setUpConvEn,
"DSSDR_ASM_TOL": self.setASMThreshold,
"DSSDR_UP_SWEEP": self.freqSweep,
"DSSDR_UP_IDLE": self.setUpIdleSequence,
"DSSDR_UP_EN": self.setUpEn,
"DSSDR_SYNC_TIME": self.syncSDRTime,
"DSSDR_UP_SWEEP": self.freqSweep,
"DSSDR_DOWN_ACQUIRE": self.acquireCarrier,
"DSSDR_INPUT_SELECT": self.setPanelSelect,
"DSSDR_REF_SELECT": self.setRefSelect,
"DSSDR_PPS_SELECT": self.setPPSSelect
}
#TODO:Add status fields for things like DSSDR_REF_LOCK
self._dssdr_channels = {
3: [7149597994, 8400061729],
4: [7150753857, 8401419752],
5: [7151909723, 8402777779],
6: [7153065586, 8404135802],
7: [7154221449, 8405493825],
8: [7155377316, 8406851853],
9: [7156533179, 8408209877],
10: [7157689045, 8409567903],
11: [7158844908, 8410925927],
12: [7160000771, 8412283950],
13: [7161156637, 8413641977],
14: [7162312500, 8415000000],
15: [7163468363, 8416358023],
16: [7164624229, 8417716050],
17: [7165780092, 8419074073],
18: [7166935955, 8420432097],
19: [7168091821, 8421790123],
20: [7169247684, 8423148147],
21: [7170403551, 8424506175],
22: [7171559414, 8425864198],
23: [7172715277, 8427222221],
24: [7173871143, 8428580248],
25: [7175027006, 8429938271],
26: [7176182869, 8431296295],
27: [7177338735, 8432654321],
28: [7178494598, 8434012345],
29: [7179650464, 8435370372],
30: [7180806327, 8436728395],
31: [7181962190, 8438086418],
32: [7183118057, 8439444446],
33: [7184273920, 8440802469],
34: [7185429783, 8442160493],
35: [7186585649, 8443518520],
36: [7187741512, 8444876543],
37: [7188897378, 8446234570],
38: [7190000000, 8450000000],
}
#FLOWGRAPH STUFF
if options.test == True:
self.u = blks2.tcp_source(
itemsize=gr.sizeof_gr_complex*1,
addr="",
port=12905,
server=True
)
elif options.fromfile == True:
self.u2 = blocks.file_meta_source("iq_in.dat")
self.u = blocks.throttle(gr.sizeof_gr_complex*1, self._samples_per_second)
elif options.frombitlog == True:
self.u3 = blocks.file_source(gr.sizeof_char, "bitstream_recording.in", True)
self.u2 = blocks.uchar_to_float()
self.u1 = blocks.throttle(gr.sizeof_float*1, self._down_bitrate)
self.u = blocks.add_const_ff(-0.5)
else:
self.u = uhd.usrp_source(device_addr=options.args, stream_args=uhd.stream_args('fc32'))
self.u.set_clock_source("external")
self.u.set_time_source("external")
self.u.set_samp_rate(self._samples_per_second)
self.u.set_antenna("RX2")
self.u.set_gain(self._down_default_gain)
self.frontend = dfi.dssdrFrontendInterface(self.u)
if options.debug_pps == True:
self.debug_pps = blocks.tag_debug(gr.sizeof_gr_complex, "debug-pps", "rx_time")
if options.tofile == True:
self.u_tx = blocks.file_meta_sink(gr.sizeof_gr_complex, "iq_out.dat", self._up_samples_per_second)
elif options.tonull == True:
self.u_tx = blocks.null_sink(gr.sizeof_gr_complex)
else:
self.u_tx = uhd.usrp_sink(device_addr=options.args, stream_args=uhd.stream_args('fc32'))
self.u_tx.set_clock_source("external")
self.u_tx.set_time_source("external")
self.u_tx.set_samp_rate(self._up_samples_per_second)
self.u_tx.set_antenna("TX/RX")
self.u_tx.set_gain(self._up_default_gain)
#GUI STUFF
if options.graphics == True:
self.nb0 = wx.Notebook(self.GetWin(), style=wx.NB_TOP)
self.nb0.AddPage(grc_wxgui.Panel(self.nb0), "RX")
self.nb0.AddPage(grc_wxgui.Panel(self.nb0), "TX")
self.nb0.AddPage(grc_wxgui.Panel(self.nb0), "Nav")
self.Add(self.nb0)
self.constellation_scope = scopesink2.scope_sink_c(
self.nb0.GetPage(0).GetWin(),
title="Scope Plot",
sample_rate=self._scope_sample_rate,
v_scale=0,
v_offset=0,
t_scale=0,
ac_couple=False,
xy_mode=True,
num_inputs=1,
trig_mode=wxgui.TRIG_MODE_AUTO,
y_axis_label="Counts",
)
self.nb0.GetPage(0).Add(self.constellation_scope.win)
#self.constellation_scope.win.set_marker('plus')
self._scope_is_fft = False
self.time_scope = scopesink2.scope_sink_f(
self.nb0.GetPage(0).GetWin(),
title="Scope Plot",
sample_rate=self._scope_sample_rate,
v_scale=0,
v_offset=0,
t_scale=.005,
ac_couple=False,
xy_mode=False,
num_inputs=1,
trig_mode=wxgui.TRIG_MODE_AUTO,
y_axis_label="Counts",
)
self.nb0.GetPage(0).Add(self.time_scope.win)
self.nb0.GetPage(0).GetWin()._box.Hide(self.time_scope.win)
self.fft_scope = fftsink2.fft_sink_c(
self.nb0.GetPage(0).GetWin(),
baseband_freq=0,
y_per_div=10,
y_divs=10,
ref_level=0,
ref_scale=2.0,
sample_rate=self._scope_sample_rate,
fft_size=1024,
fft_rate=15,
average=False,
avg_alpha=None,
title="FFT Plot",
peak_hold=False,
)
self.nb0.GetPage(0).Add(self.fft_scope.win)
self.nb0.GetPage(0).GetWin()._box.Hide(self.fft_scope.win)
self.row1_sizer = wx.BoxSizer(wx.HORIZONTAL)
self.recording_onoff_chooser = forms.radio_buttons(
parent=self.nb0.GetPage(0).GetWin(),
value='Off',
callback=self.setRecording,
label="IQ Recording",
choices=['Off','On'],
labels=[],
style=wx.RA_HORIZONTAL,
)
self.front_panel_chooser = forms.radio_buttons(
parent=self.nb0.GetPage(0).GetWin(),
value='RF',
callback=self.setPanelSelect,
label="Input Select",
choices=['RF','IF'],
labels=[],
style=wx.RA_HORIZONTAL,
)
self.ref_chooser = forms.radio_buttons(
parent=self.nb0.GetPage(0).GetWin(),
value='Internal',
callback=self.setRefSelect,
label="Ref Select",
choices=['Internal','External'],
labels=[],
style=wx.RA_HORIZONTAL,
)
self.pps_chooser = forms.radio_buttons(
parent=self.nb0.GetPage(0).GetWin(),
value='Internal',
callback=self.setPPSSelect,
label="PPS Select",
choices=['Internal','External'],
labels=[],
style=wx.RA_HORIZONTAL,
)
self.sync_button = forms.button(
parent=self.nb0.GetPage(0).GetWin(),
value='Sync to PPS',
callback=self.syncSDRTime,
choices=['Sync to PPS'],
style=wx.RA_HORIZONTAL,
)
self.ref_locked_text = forms.static_text(
parent=self.nb0.GetPage(0).GetWin(),
value="",
callback=self.setRefLocked,
label="",
converter=forms.str_converter(),
)
self.row1_sizer.Add(self.recording_onoff_chooser, flag=wx.ALIGN_CENTER)
self.row1_sizer.Add(self.front_panel_chooser, flag=wx.ALIGN_CENTER)
self.row1_sizer.Add(self.ref_chooser, flag=wx.ALIGN_CENTER)
self.row1_sizer.Add(self.pps_chooser, flag=wx.ALIGN_CENTER)
self.row1_sizer.Add(self.sync_button, flag=wx.ALIGN_CENTER)
self.row1_sizer.Add(self.ref_locked_text, flag=wx.ALIGN_CENTER)
self.nb0.GetPage(0).Add(self.row1_sizer)
self.complex_scope_chooser = forms.radio_buttons(
parent=self.nb0.GetPage(0).GetWin(),
value='Constellation',
callback=self.setComplexScopeStyle,
label="Complex Scope",
choices=['Constellation','FFT'],
labels=[],
style=wx.RA_HORIZONTAL,
)
self.nb0.GetPage(0).Add(self.complex_scope_chooser)
self.scope_chooser = forms.radio_buttons(
parent=self.nb0.GetPage(0).GetWin(),
value='USRP',
callback=self.setScopePoint,
label="Scope Probe Point",
choices=['USRP','Carrier Tracking','Sub-Carrier Costas','Sub-Carrier Sync','Data Sync'],
labels=[],
style=wx.RA_HORIZONTAL,
)
self.nb0.GetPage(0).Add(self.scope_chooser)
self._bitrate_text_box = forms.text_box(
parent=self.nb0.GetPage(0).GetWin(),
value=self._down_bitrate,
callback=self.setDownBitrate,
label="Symbol Rate",
converter=forms.float_converter(),
)
self.nb0.GetPage(0).Add(self._bitrate_text_box)
self._samprate_text_box = forms.text_box(
parent=self.nb0.GetPage(0).GetWin(),
value=self._samples_per_second,
callback=self.setSampRate,
label="Sampling Rate",
converter=forms.float_converter(),
)
self.nb0.GetPage(0).Add(self._samprate_text_box)
self._subcfreq_text_box = forms.text_box(
parent=self.nb0.GetPage(0).GetWin(),
value=self._down_sub_freq,
callback=self.setDownSubFreq,
label="Downlink Subcarrier Frequency",
converter=forms.float_converter(),
)
self.nb0.GetPage(0).Add(self._subcfreq_text_box)
self._mod_index_chooser = forms.radio_buttons(
parent=self.nb0.GetPage(0).GetWin(),
value=self._modulation_index,
callback=self.setDownModulationIndex,
label="Modulation Index",
choices=['pi/2', 'pi/3'],
labels=[],
style=wx.RA_HORIZONTAL,
)
self.nb0.GetPage(0).Add(self._mod_index_chooser)
self._pktlen_text_box = forms.text_box(
parent=self.nb0.GetPage(0).GetWin(),
value=self._tm_len,
callback=self.setTMLen,
label="Downlink Packet Length (bits)",
converter=forms.float_converter(),
)
self.nb0.GetPage(0).Add(self._pktlen_text_box)
self._coding_chooser = forms.radio_buttons(
parent=self.nb0.GetPage(0).GetWin(),
value=self._coding_method,
callback=self.setDownCodingMethod,
label="Coding",
choices=['None', 'RS', 'Turbo 1/2', 'Turbo 1/3', 'Turbo 1/4', 'Turbo 1/6'],
labels=[],
style=wx.RA_HORIZONTAL,
)
self.nb0.GetPage(0).Add(self._coding_chooser)
self._down_conv_check_box = forms.check_box(
parent=self.nb0.GetPage(0).GetWin(),
value=self._down_conv_en,
callback=self.setDownConvEn,
label="Convolutional Decode",
true="True",
false="False",
)
self.nb0.GetPage(0).Add(self._down_conv_check_box)
self._down_randomizer_check_box = forms.check_box(
parent=self.nb0.GetPage(0).GetWin(),
value=self._down_randomizer_en,
callback=self.setDownRandomizerEn,
label="De-randomizer",
true=True,
false=False,
)
self.nb0.GetPage(0).Add(self._down_randomizer_check_box)
self._down_manchester_check_box = forms.check_box(
parent=self.nb0.GetPage(0).GetWin(),
value=self._down_manchester_en,
callback=self.setDownManchesterEn,
label="Manchester Decode",
true=True,
false=False,
)
self.nb0.GetPage(0).Add(self._down_manchester_check_box)
self._pktlen_text_box = forms.text_box(
parent=self.nb0.GetPage(0).GetWin(),
value=self._asm_threshold,
callback=self.setASMThreshold,
label="ASM Error Tolerance (bits)",
converter=forms.float_converter(),
)
self.nb0.GetPage(0).Add(self._pktlen_text_box)
self._coding_chooser = forms.radio_buttons(
parent=self.nb0.GetPage(0).GetWin(),
value=self._rs_i,
callback=self.setRSI,
label="Reed-Solomon Interleaving Depth",
choices=[1,5],
labels=[],
style=wx.RA_HORIZONTAL,
)
self.nb0.GetPage(0).Add(self._coding_chooser)
self._ccsds_chan_text_box = forms.text_box(
parent=self.nb0.GetPage(0).GetWin(),
value=self._ccsds_channel,
callback=self.setChannel,
label="CCSDS Channel",
converter=forms.int_converter(),
)
self.nb0.GetPage(0).Add(self._ccsds_chan_text_box)
self.setChannel(self._ccsds_channel)
if options.test == True or options.fromfile == True or options.frombitlog == True:
glow = 0.0
ghigh = 1.0
cur_g = 0.5
else:
g = self.u.get_gain_range()
cur_g = self._down_default_gain
# some configurations don't have gain control
if g.stop() <= g.start():
glow = 0.0
ghigh = 1.0
else:
glow = g.start()
ghigh = g.stop()
self._uhd_gain_slider = wx.BoxSizer(wx.HORIZONTAL)
form.slider_field(
parent=self.nb0.GetPage(0).GetWin(),
sizer=self._uhd_gain_slider,
label="USRP RX Gain",
weight=3,
min=int(glow),
max=int(ghigh),
value=cur_g,
callback=self.setDownGain
)
self.nb0.GetPage(0).Add(self._uhd_gain_slider)
#TX chain GUI components
if options.test == True or options.tofile == True or options.tonull == True:
gtxlow = 0.0
gtxhigh = 1.0
cur_gtx = 0.5
else:
gtx = self.u_tx.get_gain_range()
cur_gtx = self._up_default_gain
# some configurations don't have gain control
if gtx.stop() <= gtx.start():
gtxlow = 0.0
gtxhigh = 1.0
else:
gtxlow = gtx.start()
gtxhigh = gtx.stop()
self._up_en_chooser = forms.check_box(
parent=self.nb0.GetPage(1).GetWin(),
value='True',
callback=self.setUpEn,
label="TX Enable",
true='True',
false='False',
)
self.nb0.GetPage(1).Add(self._up_en_chooser)
self._uhd_tx_gain_slider = wx.BoxSizer(wx.HORIZONTAL)
form.slider_field(
parent=self.nb0.GetPage(1).GetWin(),
sizer=self._uhd_tx_gain_slider,
label="USRP TX Gain",
weight=3,
min=int(gtxlow),
max=int(gtxhigh),
value=cur_gtx,
callback=self.setUpGain
)
self.nb0.GetPage(1).Add(self._uhd_tx_gain_slider)
self._subcfreq_up_text_box = forms.text_box(
parent=self.nb0.GetPage(1).GetWin(),
value=self._up_sub_freq,
callback=self.setUpSubFreq,
label="Uplink Subcarrier Frequency",
converter=forms.float_converter(),
)
self.nb0.GetPage(1).Add(self._subcfreq_up_text_box)
self._up_bitrate_text_box = forms.text_box(
parent=self.nb0.GetPage(1).GetWin(),
value=self._up_bitrate,
callback=self.setUpBitrate,
label="Uplink Bitrate",
converter=forms.float_converter(),
)
self.nb0.GetPage(1).Add(self._up_bitrate_text_box)
self._up_data_text_box = forms.text_box(
parent=self.nb0.GetPage(1).GetWin(),
value="1234ABCD",
callback=self.txData,
label="TX Data",
converter=forms.str_converter(),
)
self.nb0.GetPage(1).Add(self._up_data_text_box)
self._up_mod_index_chooser = forms.text_box(
parent=self.nb0.GetPage(1).GetWin(),
value=self._up_modulation_index,
callback=self.setUpModulationIndex,
label="Uplink Modulation Index",
converter=forms.float_converter(),
)
self.nb0.GetPage(1).Add(self._up_mod_index_chooser)
self._up_coding_chooser = forms.radio_buttons(
parent=self.nb0.GetPage(1).GetWin(),
value=self._up_coding_method,
callback=self.setUpCodingMethod,
label="Coding",
choices=['None', 'RS'],
labels=[],
style=wx.RA_HORIZONTAL,
)
self.nb0.GetPage(1).Add(self._up_coding_chooser)
self._subcarrier_chooser = forms.radio_buttons(
parent=self.nb0.GetPage(1).GetWin(),
value=self._up_subcarrier,
callback=self.setUpSubcarrier,
label="Subcarrier Type",
choices=['Square','Sine'],
labels=[],
style=wx.RA_HORIZONTAL,
)
self.nb0.GetPage(1).Add(self._subcarrier_chooser)
self._up_conv_check_box = forms.check_box(
parent=self.nb0.GetPage(1).GetWin(),
value=self._up_conv_en,
callback=self.setUpConvEn,
label="Convolutional Encode",
true="True",
false="False",
)
self.nb0.GetPage(1).Add(self._up_conv_check_box)
self._up_pktlen_text_box = forms.text_box(
parent=self.nb0.GetPage(1).GetWin(),
value=self._up_tm_len,
callback=self.setUpTMLen,
label="Uplink Packet Length (bits)",
converter=forms.float_converter(),
)
self.nb0.GetPage(1).Add(self._up_pktlen_text_box)
self._uhd_offset_text_box = forms.text_box(
parent=self.nb0.GetPage(1).GetWin(),
value=self._uhd_carrier_offset,
callback=self.setUHDCarrierOffset,
label="USRP Offset Frequency (Hz)",
converter=forms.float_converter(),
)
self.nb0.GetPage(1).Add(self._uhd_offset_text_box)
self._sweep_gen_text_box = forms.text_box(
parent=self.nb0.GetPage(1).GetWin(),
value="rf2_1",
callback=self.freqSweep,
label="Frequency Sweep Profile",
converter=forms.str_converter(),
)
self.nb0.GetPage(1).Add(self._sweep_gen_text_box)
self._idle_sequence_text_box = forms.text_box(
parent=self.nb0.GetPage(1).GetWin(),
value=self._up_idle_sequence,
callback=self.setUpIdleSequence,
label="Uplink Idle Sequence",
converter=forms.str_converter(),
)
self.nb0.GetPage(1).Add(self._idle_sequence_text_box)
self._pn_ranging_text_box = forms.text_box(
parent=self.nb0.GetPage(2).GetWin(),
value="",
callback=self.rangePN,
label="Queue PN Ranging",
converter=forms.str_converter(),
)
self.nb0.GetPage(2).Add(self._pn_ranging_text_box)
self._sequential_ranging_text_box = forms.text_box(
parent=self.nb0.GetPage(2).GetWin(),
value="",
callback=self.rangeSequential,
label="Queue Sequential Ranging",
converter=forms.str_converter(),
)
self.nb0.GetPage(2).Add(self._sequential_ranging_text_box)
self.row2_sizer = wx.BoxSizer(wx.HORIZONTAL)
self.freq_acq_button = forms.button(
parent=self.nb0.GetPage(2).GetWin(),
value='Acquire Carrier Offset',
callback=self.acquireCarrier,
choices=['Acquire Carrier Offset'],
style=wx.RA_HORIZONTAL,
)
self.carrier_offset_text = forms.static_text(
parent=self.nb0.GetPage(2).GetWin(),
value="",
label="",
converter=forms.str_converter(),
)
self.row2_sizer.Add(self.freq_acq_button, flag=wx.ALIGN_CENTER)
self.row2_sizer.Add(self.carrier_offset_text, flag=wx.ALIGN_CENTER)
self.nb0.GetPage(2).Add(self.row2_sizer)
self.file_sink = blocks.file_meta_sink(gr.sizeof_gr_complex, "iq_recording.dat", self._samples_per_second)
self.file_sink.close()
self.iq_recording_ctr = 0
# Selection logic to switch between recording and normal flowgraph routes
# NOTE: u_valve logic is implemented backwards in GNURadio....
#self.u_valve = blks2.valve(
# item_size=gr.sizeof_gr_complex,
# open=False
#)
# Temporary code used to verify coherent turnaround
self.turnaround_mixer = blocks.multiply_cc()
self.turnaround_mixer_source = analog.sig_source_c(self._down_samples_per_second, analog.GR_SIN_WAVE, -25e3, 1.0)
self.turnaround_iir = filter.single_pole_iir_filter_cc(0.0001)
self.turnaround_null = blocks.null_sink(gr.sizeof_float)
# PLL and associated carrier for tracking carrier frequency if residual carrier is used
self.carrier_tracking = sdrp.pll_freq_acq_cc(math.pi/2000, math.pi, -math.pi, int(options.acq_samples))
self.imag_to_float = blocks.complex_to_imag()
#Suppressed carrier requires costas after subcarrier mixer
self.subcarrier_costas = digital.costas_loop_cc(0.001, 2)
self.real_to_float = blocks.complex_to_real()
#Square wave subcarrier sync
self.subcarrier_sync = sdrp.square_sub_tracker_ff(0.001, 2*self._down_sub_freq/self._down_samples_per_second*1.0001, 2*self._down_sub_freq/self._down_samples_per_second*0.9999)
#Data sync
self.data_sync = sdrp.square_data_tracker_ff(0.001, self._down_bitrate/self._down_samples_per_second*1.001, self._down_bitrate/self._down_samples_per_second*0.999)
#Data framing
self.soft_correlator = sdrp.correlate_soft_access_tag_ff(conv_packed_binary_string_to_1_0_string('\x1A\xCF\xFC\x1D'), self._asm_threshold, "asm_corr")
self.conv_decoder = sdrp.ccsds_tm_conv_decoder("asm_corr")
self.de_randomizer = sdrp.ccsds_tm_derandomizer("asm_corr")
self.tm_framer = sdrp.ccsds_tm_framer(self._tm_packet_id, self._timestamp_id, "asm_corr", "rx_time", self._down_bitrate)
self.tm_framer.setFrameLength(self._tm_len)
self._current_scope_block = None
self._current_scoped_block = self.u
self._current_scoped_block_port = 0
self._recording = 'Off'
#TX path in flowgraph
self.pkt_gen_msgq = gr.msg_queue(10)
self.pkt_gen = sdrp.ccsds_tm_tx(self._tm_packet_id, self._timestamp_id, 1.0, 16, self.pkt_gen_msgq)
self.conj = blocks.conjugate_cc()
#Sweep generator for transponder lock
self.sweep_gen = sdrp.sweep_generator_cc(self._up_samples_per_second)
# DSSDR subcarrier mixer (either 25 kHz or 0 kHz depending on baud rate)
self.up_subcarrier_mixer = blocks.multiply_ff()
self.subcarrier_mixer_source_tx = analog.sig_source_f(self._up_samples_per_second, analog.GR_SQR_WAVE, 25e3, 2.0, -1.0)
self.phase_mod_tx = analog.phase_modulator_fc(self._up_modulation_index)
self.tx_attenuator = blocks.multiply_const_cc((0.1+0.0j))
#Add in bit recorder if needed
if self.options.bitlog:
self.bit_slicer = digital.binary_slicer_fb()
self.bit_recorder = blocks.file_sink(1, "bitstream_recording.out")
self.setDownCodingMethod(self._coding_method)
self.setUpCodingMethod("None")
self.setUpSubcarrier(self._up_subcarrier)
self.setDownBitrate(self._down_bitrate)
self.setUpBitrate(self._up_bitrate)
self.setDownModulationIndex(self._modulation_index)
self.setUpModulationIndex(self._up_modulation_index)
self.setDownConvEn(self._down_conv_en)
self.setUpConvEn(self._up_conv_en)
self.setUpIdleSequence(self._up_idle_sequence)
self.setDownRandomizerEn(self._down_randomizer_en)
self.setDownManchesterEn(self._down_manchester_en)
#Connection to outside world
self.socket_pdu = blocks.socket_pdu("TCP_SERVER", "127.0.0.1", "12902", 10000)
self.sdrp_interpreter = sdrp.sdrp_packet_interpreter()
self.msg_connect(self.tm_framer, "tm_frame_out", self.sdrp_interpreter, "sdrp_pdu_in")
self.msg_connect(self.sdrp_interpreter, "socket_pdu_out", self.socket_pdu, "pdus")
self.msg_connect(self.socket_pdu, "pdus", self.sdrp_interpreter, "socket_pdu_in")
self.msg_connect(self.sdrp_interpreter,"sdrp_pdu_out", self.pkt_gen, "ccsds_tx_msg_in")
if options.test == False and options.fromfile == False and options.frombitlog == False:
_threading.Thread(target=self.watchRef).start()
self.initialized = True
print "DS-SDR Initialized"
def __init__(self, samp_rate=10000):
gr.hier_block2.__init__(
self, "Sync Radio Hier Grc",
gr.io_signaturev(2, 2, [gr.sizeof_char*1, gr.sizeof_gr_complex*1]),
gr.io_signaturev(2, 2, [gr.sizeof_char*1, gr.sizeof_gr_complex*1]),
)
##################################################
# Parameters
##################################################
self.samp_rate = samp_rate
##################################################
# Variables
##################################################
self.sync_word2 = sync_word2 = [0j,0j,0j,0j,0j,0j,0j,0j,0j,0j,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),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,0j,0j,0j,0j,0j,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,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,1.42,0.0,-1.42,0.0,1.42,0.0,1.42,0.0,1.42,0.0,1.42,0.0,-1.42,0.0,1.42,0.0,1.42,0.0,-1.42,0.0,1.42,0.0,1.42,0.0,1.42,0.0,-1.42,0.0,1.42,0.0,1.42,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0]
self.pilot_symbols = pilot_symbols = ((1, -1,),)
self.pilot_carriers = pilot_carriers = ((-13, 12,),)
self.payload_mod = payload_mod = digital.constellation_qpsk()
self.packet_length_tag_key = packet_length_tag_key = "packet_len"
self.occupied_carriers = occupied_carriers = ([-16, -15, -14, -12, -11, -10, -9, -8, -7, -6, -5, -4, -3, -2, -1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15],)
self.length_tag_key = length_tag_key = "frame_len"
self.header_mod = header_mod = digital.constellation_bpsk()
self.fft_len = fft_len = (len(sync_word1)+len(sync_word2))/2
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.len_ocup_carr = len_ocup_carr = len(occupied_carriers[0])
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=True)
self.header_equalizer = header_equalizer = digital.ofdm_equalizer_simpledfe(fft_len, header_mod.base(), occupied_carriers, pilot_carriers, pilot_symbols)
self.forward_OOB = forward_OOB = [0.005277622700213007, 0.03443705907448985, 0.1214101788557494, 0.29179662246081545, 0.52428014905364, 0.7350677973792328, 0.8210395030022875, 0.7350677973792348, 0.5242801490536404, 0.291796622460816, 0.1214101788557501, 0.03443705907448997, 0.005277622700213012]
self.feedback_OOB = feedback_OOB = [1.0, -1.0455317889337852, 3.9201525346250072, -3.9114761684448958, 6.54266144224035, -5.737287389902878, 5.820328302284336, -4.134700802700442, 2.7949972248757664, -1.4584448495689168, 0.6358650797085171, -0.19847981428665007, 0.04200458351675313]
self.cp_len = cp_len = fft_len/4
self.active_carriers = active_carriers = len(occupied_carriers[0])+4
##################################################
# Blocks
##################################################
self.iir_filter_xxx_1 = filter.iir_filter_ccd((forward_OOB), (feedback_OOB), False)
self.fft_vxx_txpath = fft.fft_vcc(fft_len, False, (()), True, 1)
self.fft_vxx_2_rxpath = fft.fft_vcc(fft_len, True, (), True, 1)
self.fft_vxx_1_rxpath = fft.fft_vcc(fft_len, True, (()), True, 1)
self.digital_packet_headerparser_b_rxpath = digital.packet_headerparser_b(header_formatter.base())
self.digital_packet_headergenerator_bb_txpath = digital.packet_headergenerator_bb(header_formatter.formatter(), "packet_len")
self.digital_ofdm_sync_sc_cfb_rxpath = digital.ofdm_sync_sc_cfb(fft_len, fft_len/4, False)
self.digital_ofdm_serializer_vcc_payload_rxpath = digital.ofdm_serializer_vcc(fft_len, occupied_carriers, length_tag_key, packet_length_tag_key, 1, "", True)
self.digital_ofdm_serializer_vcc_header_rxpath = digital.ofdm_serializer_vcc(fft_len, occupied_carriers, length_tag_key, "", 0, "", True)
self.digital_ofdm_frame_equalizer_vcvc_2_rxpath = digital.ofdm_frame_equalizer_vcvc(payload_equalizer.base(), cp_len, length_tag_key, True, 0)
self.digital_ofdm_frame_equalizer_vcvc_1_rxpath = digital.ofdm_frame_equalizer_vcvc(header_equalizer.base(), cp_len, length_tag_key, True, 1)
self.digital_ofdm_cyclic_prefixer_txpath = digital.ofdm_cyclic_prefixer(fft_len, fft_len+cp_len, rolloff, packet_length_tag_key)
(self.digital_ofdm_cyclic_prefixer_txpath).set_min_output_buffer(24000)
self.digital_ofdm_chanest_vcvc_rxpath = digital.ofdm_chanest_vcvc((sync_word1), (sync_word2), 1, 0, 3, False)
self.digital_ofdm_carrier_allocator_cvc_txpath = digital.ofdm_carrier_allocator_cvc(fft_len, occupied_carriers, pilot_carriers, pilot_symbols, (sync_word1, sync_word2), packet_length_tag_key)
(self.digital_ofdm_carrier_allocator_cvc_txpath).set_min_output_buffer(16000)
self.digital_header_payload_demux_rxpath = digital.header_payload_demux(
3,
fft_len,
cp_len,
length_tag_key,
"",
True,
gr.sizeof_gr_complex,
"rx_time",
samp_rate,
(),
)
#self.digital_crc32_bb_txpath = digital.crc32_bb(False, packet_length_tag_key)
#self.digital_crc32_bb_rxpath = digital.crc32_bb(True, packet_length_tag_key)
self.digital_constellation_decoder_cb_1_rxpath = digital.constellation_decoder_cb(payload_mod.base())
self.digital_constellation_decoder_cb_0 = digital.constellation_decoder_cb(header_mod.base())
self.digital_chunks_to_symbols_x_txpath = digital.chunks_to_symbols_bc((payload_mod.points()), 1)
self.digital_chunks_to_symbols_txpath = digital.chunks_to_symbols_bc((header_mod.points()), 1)
self.blocks_tagged_stream_mux_txpath = blocks.tagged_stream_mux(gr.sizeof_gr_complex*1, packet_length_tag_key, 0)
(self.blocks_tagged_stream_mux_txpath).set_min_output_buffer(16000)
self.blocks_tag_gate_txpath = blocks.tag_gate(gr.sizeof_gr_complex * 1, False)
self.blocks_repack_bits_bb_txpath = blocks.repack_bits_bb(8, payload_mod.bits_per_symbol(), packet_length_tag_key, False)
self.blocks_repack_bits_bb_rxpath = blocks.repack_bits_bb(payload_mod.bits_per_symbol(), 8, packet_length_tag_key, True)
self.blocks_multiply_xx_rxpath = blocks.multiply_vcc(1)
self.blocks_multiply_const_vxx_0 = blocks.multiply_const_vcc((.01, ))
self.blocks_delay_rxpath = blocks.delay(gr.sizeof_gr_complex*1, fft_len+fft_len/4)
self.analog_frequency_modulator_fc_rxpath = analog.frequency_modulator_fc(-2.0/fft_len)
self.analog_agc2_xx_0 = analog.agc2_cc(1e-1, 1e-2, 1.0, 1.0)
self.analog_agc2_xx_0.set_max_gain(65536)
##################################################
# Connections
##################################################
self.connect((self.digital_ofdm_sync_sc_cfb_rxpath, 0), (self.analog_frequency_modulator_fc_rxpath, 0))
self.connect((self.analog_agc2_xx_0, 0), (self.digital_ofdm_sync_sc_cfb_rxpath, 0))
self.connect((self.analog_agc2_xx_0, 0), (self.blocks_delay_rxpath, 0))
self.connect((self.digital_packet_headergenerator_bb_txpath, 0), (self.digital_chunks_to_symbols_txpath, 0))
self.connect((self.blocks_repack_bits_bb_txpath, 0), (self.digital_chunks_to_symbols_x_txpath, 0))
self.connect((self.digital_chunks_to_symbols_txpath, 0), (self.blocks_tagged_stream_mux_txpath, 0))
self.connect((self.digital_chunks_to_symbols_x_txpath, 0), (self.blocks_tagged_stream_mux_txpath, 1))
self.connect((self.digital_ofdm_cyclic_prefixer_txpath, 0), (self.blocks_tag_gate_txpath, 0))
self.connect((self.fft_vxx_txpath, 0), (self.digital_ofdm_cyclic_prefixer_txpath, 0))
self.connect((self.blocks_tagged_stream_mux_txpath, 0), (self.digital_ofdm_carrier_allocator_cvc_txpath, 0))
self.connect((self.digital_ofdm_carrier_allocator_cvc_txpath, 0), (self.fft_vxx_txpath, 0))
self.connect((self.digital_header_payload_demux_rxpath, 0), (self.fft_vxx_1_rxpath, 0))
self.connect((self.fft_vxx_1_rxpath, 0), (self.digital_ofdm_chanest_vcvc_rxpath, 0))
self.connect((self.digital_ofdm_frame_equalizer_vcvc_1_rxpath, 0), (self.digital_ofdm_serializer_vcc_header_rxpath, 0))
self.connect((self.digital_ofdm_chanest_vcvc_rxpath, 0), (self.digital_ofdm_frame_equalizer_vcvc_1_rxpath, 0))
self.connect((self.digital_header_payload_demux_rxpath, 1), (self.fft_vxx_2_rxpath, 0))
self.connect((self.digital_ofdm_frame_equalizer_vcvc_2_rxpath, 0), (self.digital_ofdm_serializer_vcc_payload_rxpath, 0))
self.connect((self.fft_vxx_2_rxpath, 0), (self.digital_ofdm_frame_equalizer_vcvc_2_rxpath, 0))
self.connect((self.digital_ofdm_serializer_vcc_payload_rxpath, 0), (self.digital_constellation_decoder_cb_1_rxpath, 0))
self.connect((self.digital_constellation_decoder_cb_1_rxpath, 0), (self.blocks_repack_bits_bb_rxpath, 0))
self.connect((self.digital_ofdm_serializer_vcc_header_rxpath, 0), (self.digital_constellation_decoder_cb_0, 0))
self.connect((self.analog_frequency_modulator_fc_rxpath, 0), (self.blocks_multiply_xx_rxpath, 0))
self.connect((self.digital_ofdm_sync_sc_cfb_rxpath, 1), (self.digital_header_payload_demux_rxpath, 1))
self.connect((self.blocks_multiply_xx_rxpath, 0), (self.digital_header_payload_demux_rxpath, 0))
self.connect((self.blocks_delay_rxpath, 0), (self.blocks_multiply_xx_rxpath, 1))
self.connect((self.digital_constellation_decoder_cb_0, 0), (self.digital_packet_headerparser_b_rxpath, 0))
self.connect((self, 1), (self.analog_agc2_xx_0, 0))
self.connect((self.blocks_tag_gate_txpath, 0), (self.blocks_multiply_const_vxx_0, 0))
self.connect((self.iir_filter_xxx_1, 0), (self, 1))
'''
self.connect((self, 0), (self.digital_crc32_bb_txpath, 0))
self.connect((self.digital_crc32_bb_txpath, 0), (self.blocks_repack_bits_bb_txpath, 0))
self.connect((self.digital_crc32_bb_txpath, 0), (self.digital_packet_headergenerator_bb_txpath, 0))
'''
self.connect((self, 0), (self.blocks_repack_bits_bb_txpath, 0))
self.connect((self, 0), (self.digital_packet_headergenerator_bb_txpath, 0))
'''
self.connect((self.blocks_repack_bits_bb_rxpath, 0), (self.digital_crc32_bb_rxpath, 0))
self.connect((self.digital_crc32_bb_rxpath, 0), (self, 0))
'''
self.connect((self.blocks_repack_bits_bb_rxpath, 0), (self, 0))
self.connect((self.blocks_multiply_const_vxx_0, 0), (self.iir_filter_xxx_1, 0))
##################################################
# Asynch Message Connections
##################################################
self.msg_connect(self.digital_packet_headerparser_b_rxpath, "header_data", self.digital_header_payload_demux_rxpath, "header_data")
def __init__(self,
antenna=satnogs.not_set_antenna,
baudrate=9600.0,
bb_gain=satnogs.not_set_rx_bb_gain,
decoded_data_file_path='/tmp/.satnogs/data/data',
dev_args=satnogs.not_set_dev_args,
doppler_correction_per_sec=1000,
enable_iq_dump=0,
excess_bw=0.35,
file_path='test.wav',
if_gain=satnogs.not_set_rx_if_gain,
iq_file_path='/tmp/iq.dat',
lo_offset=100e3,
max_cfo=1000.0,
ppm=0,
rf_gain=satnogs.not_set_rx_rf_gain,
rigctl_port=4532,
rx_freq=100e6,
rx_sdr_device='usrpb200',
samp_rate_rx=satnogs.not_set_samp_rate_rx,
udp_IP='127.0.0.1',
udp_port=16887,
waterfall_file_path='/tmp/waterfall.dat'):
gr.top_block.__init__(self, "satnogs_bpsk_ax25")
##################################################
# Parameters
##################################################
self.antenna = antenna
self.baudrate = baudrate
self.bb_gain = bb_gain
self.decoded_data_file_path = decoded_data_file_path
self.dev_args = dev_args
self.doppler_correction_per_sec = doppler_correction_per_sec
self.enable_iq_dump = enable_iq_dump
self.excess_bw = excess_bw
self.file_path = file_path
self.if_gain = if_gain
self.iq_file_path = iq_file_path
self.lo_offset = lo_offset
self.max_cfo = max_cfo
self.ppm = ppm
self.rf_gain = rf_gain
self.rigctl_port = rigctl_port
self.rx_freq = rx_freq
self.rx_sdr_device = rx_sdr_device
self.samp_rate_rx = samp_rate_rx
self.udp_IP = udp_IP
self.udp_port = udp_port
self.waterfall_file_path = waterfall_file_path
##################################################
# Variables
##################################################
self.sps = sps = 4
self.nfilts = nfilts = 32
self.rrc_taps = rrc_taps = firdes.root_raised_cosine(
nfilts, nfilts, 1.0 / float(sps), excess_bw, 11 * sps * nfilts)
self.bpsk_constellation = bpsk_constellation = digital.constellation_bpsk(
).base()
self.audio_samp_rate = audio_samp_rate = 48000
##################################################
# Blocks
##################################################
self.satnogs_waterfall_sink_0 = satnogs.waterfall_sink(
audio_samp_rate, 0.0, 10, 1024, waterfall_file_path, 1)
self.satnogs_udp_msg_sink_0_0 = satnogs.udp_msg_sink(
udp_IP, udp_port, 1500)
self.satnogs_tcp_rigctl_msg_source_0 = satnogs.tcp_rigctl_msg_source(
"127.0.0.1", rigctl_port, False, 1000, 1500)
self.satnogs_ogg_encoder_0 = satnogs.ogg_encoder(
file_path, audio_samp_rate, 1.0)
self.satnogs_iq_sink_0 = satnogs.iq_sink(16768, iq_file_path, False,
enable_iq_dump)
self.satnogs_frame_file_sink_0_1_0 = satnogs.frame_file_sink(
decoded_data_file_path, 0)
self.satnogs_coarse_doppler_correction_cc_0 = satnogs.coarse_doppler_correction_cc(
rx_freq, satnogs.handle_samp_rate_rx(rx_sdr_device, samp_rate_rx))
self.satnogs_ax25_decoder_bm_0_0 = satnogs.ax25_decoder_bm(
'GND', 0, True, False, 1024)
self.satnogs_ax25_decoder_bm_0 = satnogs.ax25_decoder_bm(
'GND', 0, True, True, 1024)
self.pfb_arb_resampler_xxx_0_0 = pfb.arb_resampler_ccf(
(1.0 * sps * baudrate) /
satnogs.handle_samp_rate_rx(rx_sdr_device, samp_rate_rx),
taps=None,
flt_size=32)
self.pfb_arb_resampler_xxx_0_0.declare_sample_delay(0)
self.pfb_arb_resampler_xxx_0 = pfb.arb_resampler_ccf(
audio_samp_rate /
satnogs.handle_samp_rate_rx(rx_sdr_device, samp_rate_rx),
taps=None,
flt_size=32)
self.pfb_arb_resampler_xxx_0.declare_sample_delay(0)
self.osmosdr_source_0 = osmosdr.source(
args="numchan=" + str(1) + " " +
satnogs.handle_rx_dev_args(rx_sdr_device, dev_args))
self.osmosdr_source_0.set_sample_rate(
satnogs.handle_samp_rate_rx(rx_sdr_device, samp_rate_rx))
self.osmosdr_source_0.set_center_freq(rx_freq - lo_offset, 0)
self.osmosdr_source_0.set_freq_corr(ppm, 0)
self.osmosdr_source_0.set_dc_offset_mode(2, 0)
self.osmosdr_source_0.set_iq_balance_mode(0, 0)
self.osmosdr_source_0.set_gain_mode(False, 0)
self.osmosdr_source_0.set_gain(
satnogs.handle_rx_rf_gain(rx_sdr_device, rf_gain), 0)
self.osmosdr_source_0.set_if_gain(
satnogs.handle_rx_if_gain(rx_sdr_device, if_gain), 0)
self.osmosdr_source_0.set_bb_gain(
satnogs.handle_rx_bb_gain(rx_sdr_device, bb_gain), 0)
self.osmosdr_source_0.set_antenna(
satnogs.handle_rx_antenna(rx_sdr_device, antenna), 0)
self.osmosdr_source_0.set_bandwidth(
satnogs.handle_samp_rate_rx(rx_sdr_device, samp_rate_rx), 0)
self.low_pass_filter_0_0 = filter.fir_filter_ccf(
1,
firdes.low_pass(1, audio_samp_rate,
((1.0 + excess_bw) * baudrate / 2.0) +
min(baudrate, abs(max_cfo)), baudrate / 10.0,
firdes.WIN_HAMMING, 6.76))
self.low_pass_filter_0 = filter.fir_filter_ccf(
1,
firdes.low_pass(1, sps * baudrate,
((1.0 + excess_bw) * baudrate / 2.0) +
min(baudrate, abs(max_cfo)), baudrate / 10.0,
firdes.WIN_HAMMING, 6.76))
self.digital_pfb_clock_sync_xxx_0 = digital.pfb_clock_sync_ccf(
sps, 2.0 * math.pi / 100.0, (rrc_taps), nfilts, nfilts / 2, 1.5, 1)
self.digital_costas_loop_cc_0_0 = digital.costas_loop_cc(
2.0 * math.pi / 100.0, 2, True)
self.digital_constellation_receiver_cb_0 = digital.constellation_receiver_cb(
bpsk_constellation, 2.0 * math.pi / 100.0, -0.25, 0.25)
self.blocks_rotator_cc_0_0 = blocks.rotator_cc(
2.0 * math.pi * (1200.0 / audio_samp_rate))
self.blocks_rotator_cc_0 = blocks.rotator_cc(
-2.0 * math.pi *
(lo_offset /
satnogs.handle_samp_rate_rx(rx_sdr_device, samp_rate_rx)))
self.blocks_complex_to_real_0 = blocks.complex_to_real(1)
self.analog_agc2_xx_0_0 = analog.agc2_cc(0.01, 0.001, 0.015, 1.0)
self.analog_agc2_xx_0_0.set_max_gain(65536)
self.analog_agc2_xx_0 = analog.agc2_cc(1e-3, 1e-3, 0.5, 1.0)
self.analog_agc2_xx_0.set_max_gain(65536)
##################################################
# Connections
##################################################
self.msg_connect((self.satnogs_ax25_decoder_bm_0, 'pdu'),
(self.satnogs_frame_file_sink_0_1_0, 'frame'))
self.msg_connect((self.satnogs_ax25_decoder_bm_0, 'pdu'),
(self.satnogs_udp_msg_sink_0_0, 'in'))
self.msg_connect((self.satnogs_ax25_decoder_bm_0_0, 'pdu'),
(self.satnogs_frame_file_sink_0_1_0, 'frame'))
self.msg_connect((self.satnogs_ax25_decoder_bm_0_0, 'pdu'),
(self.satnogs_udp_msg_sink_0_0, 'in'))
self.msg_connect((self.satnogs_tcp_rigctl_msg_source_0, 'freq'),
(self.satnogs_coarse_doppler_correction_cc_0, 'freq'))
self.connect((self.analog_agc2_xx_0, 0),
(self.pfb_arb_resampler_xxx_0_0, 0))
self.connect((self.analog_agc2_xx_0_0, 0),
(self.low_pass_filter_0_0, 0))
self.connect((self.blocks_complex_to_real_0, 0),
(self.satnogs_ogg_encoder_0, 0))
self.connect((self.blocks_rotator_cc_0, 0),
(self.satnogs_coarse_doppler_correction_cc_0, 0))
self.connect((self.blocks_rotator_cc_0_0, 0),
(self.blocks_complex_to_real_0, 0))
self.connect((self.digital_constellation_receiver_cb_0, 0),
(self.satnogs_ax25_decoder_bm_0, 0))
self.connect((self.digital_constellation_receiver_cb_0, 0),
(self.satnogs_ax25_decoder_bm_0_0, 0))
self.connect((self.digital_costas_loop_cc_0_0, 0),
(self.digital_pfb_clock_sync_xxx_0, 0))
self.connect((self.digital_pfb_clock_sync_xxx_0, 0),
(self.digital_constellation_receiver_cb_0, 0))
self.connect((self.low_pass_filter_0, 0),
(self.digital_costas_loop_cc_0_0, 0))
self.connect((self.low_pass_filter_0_0, 0),
(self.blocks_rotator_cc_0_0, 0))
self.connect((self.osmosdr_source_0, 0), (self.blocks_rotator_cc_0, 0))
self.connect((self.pfb_arb_resampler_xxx_0, 0),
(self.analog_agc2_xx_0_0, 0))
self.connect((self.pfb_arb_resampler_xxx_0, 0),
(self.satnogs_iq_sink_0, 0))
self.connect((self.pfb_arb_resampler_xxx_0, 0),
(self.satnogs_waterfall_sink_0, 0))
self.connect((self.pfb_arb_resampler_xxx_0_0, 0),
(self.low_pass_filter_0, 0))
self.connect((self.satnogs_coarse_doppler_correction_cc_0, 0),
(self.analog_agc2_xx_0, 0))
self.connect((self.satnogs_coarse_doppler_correction_cc_0, 0),
(self.pfb_arb_resampler_xxx_0, 0))
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,
baudrate,
samp_rate,
iq,
f_offset=None,
differential=False,
manchester=False,
dump_path=None,
options=None):
gr.hier_block2.__init__(
self, "bpsk_demodulator",
gr.io_signature(1, 1,
gr.sizeof_gr_complex if iq else gr.sizeof_float),
gr.io_signature(1, 1, gr.sizeof_float))
options_block.__init__(self, options)
if dump_path is not None:
dump_path = pathlib.Path(dump_path)
if manchester:
baudrate *= 2
# prevent problems due to baudrate too high
if baudrate >= samp_rate / 4:
print(
f'Sample rate {samp_rate} sps insufficient for {baudrate} baud BPSK demodulation. Demodulator will not work.',
file=sys.stderr)
baudrate = samp_rate / 4
sps = samp_rate / baudrate
max_sps = 10
if sps > max_sps:
decimation = ceil(sps / max_sps)
else:
decimation = 1
sps /= decimation
filter_cutoff = baudrate * 2.0
filter_transition = baudrate * 0.2
if f_offset is None:
f_offset = self.options.f_offset
if f_offset is None:
if iq:
f_offset = 0
elif baudrate <= 2400:
f_offset = 1500
else:
f_offset = 12000
taps = firdes.low_pass(1, samp_rate, filter_cutoff, filter_transition)
f = filter.freq_xlating_fir_filter_ccf if iq else filter.freq_xlating_fir_filter_fcf
self.xlating = f(decimation, taps, f_offset, samp_rate)
agc_constant = 2e-2 / sps # This gives a time constant of 50 symbols
self.agc = rms_agc(agc_constant, 1)
if dump_path is not None:
self.agc_in = blocks.file_sink(gr.sizeof_gr_complex,
str(dump_path / 'agc_in.c64'),
False)
self.agc_out = blocks.file_sink(gr.sizeof_gr_complex,
str(dump_path / 'agc_out.c64'),
False)
self.connect(self.xlating, self.agc_in)
self.connect(self.agc, self.agc_out)
if not self.options.disable_fll:
fll_bw = 2 * pi * decimation / samp_rate * self.options.fll_bw
self.fll = digital.fll_band_edge_cc(sps, self.options.rrc_alpha,
100, fll_bw)
if dump_path is not None:
self.fll_freq = blocks.file_sink(
gr.sizeof_float, str(dump_path / 'fll_freq.f32'), False)
self.fll_phase = blocks.file_sink(
gr.sizeof_float, str(dump_path / 'fll_phase.f32'), False)
self.fll_error = blocks.file_sink(
gr.sizeof_float, str(dump_path / 'fll_error.f32'), False)
self.connect((self.fll, 1), self.fll_freq)
self.connect((self.fll, 2), self.fll_phase)
self.connect((self.fll, 3), self.fll_error)
nfilts = 16
rrc_taps = firdes.root_raised_cosine(nfilts, nfilts, 1.0 / float(sps),
self.options.rrc_alpha,
int(ceil(11 * sps * nfilts)))
ted_gain = 0.5 # "empiric" formula for TED gain of a PFB MF TED for complex BPSK 0.5 sample^{-1}
damping = 1.0
self.clock_recovery = digital.symbol_sync_cc(
digital.TED_SIGNAL_TIMES_SLOPE_ML, sps, self.options.clk_bw,
damping, ted_gain, self.options.clk_limit * sps, 1,
digital.constellation_bpsk().base(), digital.IR_PFB_MF, nfilts,
rrc_taps)
if dump_path is not None:
self.clock_recovery_out = blocks.file_sink(
gr.sizeof_gr_complex,
str(dump_path / 'clock_recovery_out.c64'), False)
self.clock_recovery_err = blocks.file_sink(
gr.sizeof_float, str(dump_path / 'clock_recovery_err.f32'),
False)
self.clock_recovery_T_inst = blocks.file_sink(
gr.sizeof_float, str(dump_path / 'clock_recovery_T_inst.f32'),
False)
self.clock_recovery_T_avg = blocks.file_sink(
gr.sizeof_float, str(dump_path / 'clock_recovery_T_avg.f32'),
False)
self.connect(self.clock_recovery, self.clock_recovery_out)
self.connect((self.clock_recovery, 1), self.clock_recovery_err)
self.connect((self.clock_recovery, 2), self.clock_recovery_T_inst)
self.connect((self.clock_recovery, 3), self.clock_recovery_T_avg)
self.connect(self, self.xlating, self.agc)
if self.options.disable_fll:
self.connect(self.agc, self.clock_recovery)
else:
self.connect(self.agc, self.fll, self.clock_recovery)
self.complex_to_real = blocks.complex_to_real(1)
if manchester:
self.manchester = manchester_sync(self.options.manchester_history)
self.connect(self.clock_recovery, self.manchester)
else:
self.manchester = self.clock_recovery
if differential:
self.delay = blocks.delay(gr.sizeof_gr_complex, 1)
self.multiply_conj = blocks.multiply_conjugate_cc(1)
sign = -1 if manchester else 1
self.multiply_const = blocks.multiply_const_ff(
sign, 1) # take care about inverion in Manchester
self.connect(self.manchester, (self.multiply_conj, 0))
self.connect(self.manchester, self.delay, (self.multiply_conj, 1))
self.connect(self.multiply_conj, self.complex_to_real,
self.multiply_const, self)
else:
costas_bw = 2 * pi / baudrate * self.options.costas_bw
self.costas = digital.costas_loop_cc(costas_bw, 2, False)
if dump_path is not None:
self.costas_out = blocks.file_sink(
gr.sizeof_gr_complex, str(dump_path / 'costas_out.c64'),
False)
self.costas_frequency = blocks.file_sink(
gr.sizeof_float, str(dump_path / 'costas_frequency.f32'),
False)
self.costas_phase = blocks.file_sink(
gr.sizeof_float, str(dump_path / 'costas_phase.f32'),
False)
self.costas_error = blocks.file_sink(
gr.sizeof_float, str(dump_path / 'costas_error.f32'),
False)
self.connect(self.costas, self.costas_out)
self.connect((self.costas, 1), self.costas_frequency)
self.connect((self.costas, 2), self.costas_phase)
self.connect((self.costas, 3), self.costas_error)
self.connect(self.manchester, self.costas, self.complex_to_real,
self)
def main():
args = get_arguments()
constellation = {
1:digital.constellation_bpsk(),
2:digital.constellation_qpsk(),
3:digital.constellation_8psk(),
}
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),])
tb1 = gr.top_block()
tb2 = gr.top_block()
if args.freq == None:
data_source = mimoots.file_source2(
itemsize=(gr.sizeof_gr_complex, gr.sizeof_gr_complex),
filename=(
'1.'.join(args.from_file.rsplit('.',1)),
'2.'.join(args.from_file.rsplit('.',1))
)
)
else:
data_source = mimoots.uhd_source2(freq=args.freq, gain=args.gain)
skip1 = blocks.skiphead(
itemsize=gr.sizeof_gr_complex,
nitems_to_skip=args.skiphead
)
skip2 = blocks.skiphead(
itemsize=gr.sizeof_gr_complex,
nitems_to_skip=args.skiphead
)
ofdm_frames1 = mimoots.ofdm_basebandsignal_to_frames_cvc(
fft_len=fft_len,
cp_len=cp_len,
nofdm_symbols=args.nsymbols
)
ofdm_frames2 = mimoots.ofdm_basebandsignal_to_frames_cvc(
fft_len=fft_len,
cp_len=cp_len,
nofdm_symbols=args.nsymbols
)
buffer1 = blocks.vector_sink_c(vlen=fft_len)
buffer2 = blocks.vector_sink_c(vlen=fft_len)
tb1.connect((data_source,0), skip1, ofdm_frames1, buffer1)
tb1.connect((data_source,1), skip2, ofdm_frames2, buffer2)
tb1.connect(ofdm_frames1, blocks.file_sink(gr.sizeof_gr_complex*fft_len, 'd1.gr'))
tb1.connect(ofdm_frames2, blocks.file_sink(gr.sizeof_gr_complex*fft_len, 'd2.gr'))
tb1.run()
data1 = buffer1.data()
data2 = buffer2.data()
data1 = extract_channeldata(data1, 0, fft_len)
data2 = extract_channeldata(data2, 1, fft_len)
buf_source1 = blocks.vector_source_c(data1)
buf_source2 = blocks.vector_source_c(data2)
data_sink = mimoots.file_sink2(
itemsize=(gr.sizeof_gr_complex, gr.sizeof_gr_complex),
filename=('s1r1.gr', 's2r2.gr')
)
tb2.connect(buf_source1, (data_sink, 0))
tb2.connect(buf_source2, (data_sink, 1))
tb2.run()