def __init__(self, vlen_in=1, vlen_out=1, n_tailbits=1, denom_mother_code_rate=1, gen_poly=0, bits_per_symbol=1, map_tab=0, pp_0=0, interl_seq_0_2=(0,1), interl_seq_1_2=(0,1), pp_1_tail=0, pp_1=0, pp_0_tail=0, M_total=0, part_len_top=3, part_len_bot=3):
gr.hier_block2.__init__(
self, "MLC 16QAM",
gr.io_signature(1, 1, gr.sizeof_char*1),
gr.io_signature(1, 1, gr.sizeof_gr_complex*1),
)
##################################################
# Parameters
##################################################
self.vlen_in = vlen_in
self.vlen_out = vlen_out
self.n_tailbits = n_tailbits
self.denom_mother_code_rate = denom_mother_code_rate
self.gen_poly = gen_poly
self.bits_per_symbol = bits_per_symbol
self.map_tab = map_tab
self.pp_0 = pp_0
self.interl_seq_0_2 = interl_seq_0_2
self.interl_seq_1_2 = interl_seq_1_2
self.pp_1_tail = pp_1_tail
self.pp_1 = pp_1
self.pp_0_tail = pp_0_tail
self.M_total = M_total
self.part_len_top = part_len_top
self.part_len_bot = part_len_bot
##################################################
# Blocks
##################################################
self.trellis_encoder_xx_top = trellis.encoder_bb(trellis.fsm(1, denom_mother_code_rate, gen_poly), 0, 0) if False else trellis.encoder_bb(trellis.fsm(1, denom_mother_code_rate, gen_poly), 0)
self.trellis_encoder_xx_bot = trellis.encoder_bb(trellis.fsm(1, denom_mother_code_rate, gen_poly), 0, 0) if False else trellis.encoder_bb(trellis.fsm(1, denom_mother_code_rate, gen_poly), 0)
self.drm_qam_map_bc_0 = drm.qam_map_bc(map_tab, bits_per_symbol, vlen_out, 2)
self.drm_punct_bb_top = drm.punct_bb(pp_0, pp_0_tail, (part_len_top + n_tailbits) * denom_mother_code_rate, vlen_out * 2, n_tailbits * denom_mother_code_rate)
self.drm_punct_bb_bot = drm.punct_bb(pp_1, pp_1_tail, (part_len_bot + n_tailbits) * denom_mother_code_rate, vlen_out * 2, n_tailbits * denom_mother_code_rate)
self.drm_partitioning_16_bb_0 = drm.partitioning_bb(vlen_in, M_total)
self.drm_interleaver_bb_top = drm.interleaver_bb((interl_seq_0_2))
self.drm_interleaver_bb_bot = drm.interleaver_bb((interl_seq_1_2))
self.drm_add_tailbits_bb_top = drm.add_tailbits_bb(part_len_top, n_tailbits)
self.drm_add_tailbits_bb_bot = drm.add_tailbits_bb(part_len_bot, n_tailbits)
self.blocks_unpack_k_bits_bb_top = blocks.unpack_k_bits_bb(denom_mother_code_rate)
self.blocks_unpack_k_bits_bb_bot = blocks.unpack_k_bits_bb(denom_mother_code_rate)
##################################################
# Connections
##################################################
self.connect((self.blocks_unpack_k_bits_bb_bot, 0), (self.drm_punct_bb_bot, 0))
self.connect((self.blocks_unpack_k_bits_bb_top, 0), (self.drm_punct_bb_top, 0))
self.connect((self.drm_add_tailbits_bb_bot, 0), (self.trellis_encoder_xx_bot, 0))
self.connect((self.drm_add_tailbits_bb_top, 0), (self.trellis_encoder_xx_top, 0))
self.connect((self.drm_interleaver_bb_bot, 0), (self.drm_qam_map_bc_0, 1))
self.connect((self.drm_interleaver_bb_top, 0), (self.drm_qam_map_bc_0, 0))
self.connect((self.drm_partitioning_16_bb_0, 1), (self.drm_add_tailbits_bb_bot, 0))
self.connect((self.drm_partitioning_16_bb_0, 0), (self.drm_add_tailbits_bb_top, 0))
self.connect((self.drm_punct_bb_bot, 0), (self.drm_interleaver_bb_bot, 0))
self.connect((self.drm_punct_bb_top, 0), (self.drm_interleaver_bb_top, 0))
self.connect((self.drm_qam_map_bc_0, 0), (self, 0))
self.connect((self, 0), (self.drm_partitioning_16_bb_0, 0))
self.connect((self.trellis_encoder_xx_bot, 0), (self.blocks_unpack_k_bits_bb_bot, 0))
self.connect((self.trellis_encoder_xx_top, 0), (self.blocks_unpack_k_bits_bb_top, 0))
def __init__(self): gr.hier_block2.__init__(self, "dvb_convolutional_encoder_bb", gr.io_signature(1, 1, gr.sizeof_char), # Input signature gr.io_signature(1, 1, gr.sizeof_char)) # Output signature self.encoder = trellis.encoder_bb(trellis.fsm(dvb_swig.k, dvb_swig.n, G), dvb_swig.dimensionality) self.unpack = gr.unpack_k_bits_bb(dvb_swig.dimensionality) self.connect(self, self.encoder, self.unpack, self)
def __init__(self):
gr.hier_block2.__init__(
self,
"dvb_convolutional_encoder_bb",
gr.io_signature(1, 1, gr.sizeof_char), # Input signature
gr.io_signature(1, 1, gr.sizeof_char)) # Output signature
self.encoder = trellis.encoder_bb(
trellis.fsm(dvb_swig.k, dvb_swig.n, G), dvb_swig.dimensionality)
self.unpack = gr.unpack_k_bits_bb(dvb_swig.dimensionality)
self.connect(self, self.encoder, self.unpack, self)
def setup_test06(self):
print "... benchmarking QPSK + 1/2 mapper"
self.mode = 2
bitspermode = [0.5, 1, 1.5, 2, 3, 4, 4.5, 5, 6]
self.nobits = (bitspermode[self.mode - 1])
self.data_subcarriers = 200
self.blks = self.N * (int(10 * self.nobits) + 1)
self.tb = gr.top_block()
self.encoder = self._encoder = trellis.encoder_bb(self.fo, 0)
self.unpack = self._unpack = blocks.unpack_k_bits_bb(2)
self.puncturing = puncture_bb(self.data_subcarriers)
bmaptrig_stream_puncturing = [1] + [0] * (10 / 2 - 1)
self.bitmap_trigger_puncturing = blocks.vector_source_b(
bmaptrig_stream_puncturing, True)
if self.interleave:
int_object = trellis.interleaver(2000, 666)
self.interlv = trellis.permutation(int_object.K(),
int_object.INTER(), 1,
gr.sizeof_char)
#self.bitmap = [self.nobits]*self.data_subcarriers
self.bitmap = [self.mode] * self.data_subcarriers + [
self.mode
] * self.data_subcarriers
#self.bmaptrig_stream = [1, 1]+[0]*(11-2)
self.bitdata = [
randint(0, 1) for i in range(self.blks * self.data_subcarriers)
]
#self.data = numpy.array(self.bitdata)*(-2)+1
self.src = blocks.vector_source_b(self.bitdata)
self.bitmap_src = blocks.vector_source_b(self.bitmap, True,
self.data_subcarriers)
#self.bitmap_trigger = blocks.vector_source_b(self.bmaptrig_stream, True)
self.modulator = generic_mapper_bcv(self.data_subcarriers, self.coding,
10)
self.snk = blocks.null_sink(gr.sizeof_gr_complex *
self.data_subcarriers)
#Connect blocks
self.tb.connect(self.src, self.encoder, self.unpack,
self.puncturing) #, self.modulator, self.snk)
if self.interleave:
self.tb.connect(self.puncturing, self.interlv, self.modulator,
self.snk)
else:
self.tb.connect(self.puncturing, self.modulator, self.snk)
self.tb.connect(self.bitmap_src, (self.modulator, 1))
self.tb.connect(self.bitmap_src, (self.puncturing, 1))
self.tb.connect(self.bitmap_trigger_puncturing, (self.puncturing, 2))
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 __init__(self, vlen_in=1, vlen_out=1, n_tailbits=6, denom_mother_code_rate=6, gen_poly=(91, 121, 101, 91, 121, 101), N=1, bits_per_symbol=0, pp=0, pp_tail=0, interl_seq=range(2), map_tab=0): gr.hier_block2.__init__( self, "DRM MLC 4-QAM", gr.io_signature(1, 1, gr.sizeof_char*vlen_in), gr.io_signature(1, 1, gr.sizeof_gr_complex*vlen_out), ) ################################################## # Parameters ################################################## self.vlen_in = vlen_in self.vlen_out = vlen_out self.n_tailbits = n_tailbits self.denom_mother_code_rate = denom_mother_code_rate self.gen_poly = gen_poly self.N = N self.bits_per_symbol = bits_per_symbol self.pp = pp self.pp_tail = pp_tail self.interl_seq = interl_seq self.map_tab = map_tab ################################################## # Blocks ################################################## self.trellis_encoder_xx_0 = trellis.encoder_bb(trellis.fsm(1, denom_mother_code_rate, gen_poly), 0) self.gr_vector_to_stream_0 = gr.vector_to_stream(gr.sizeof_char*1, vlen_in + n_tailbits) self.gr_unpack_k_bits_bb_0 = gr.unpack_k_bits_bb(denom_mother_code_rate) self.gr_stream_to_vector_0 = gr.stream_to_vector(gr.sizeof_char*1, (vlen_in + n_tailbits) * denom_mother_code_rate) self.drm_qam_map_vbvc_0 = drm.qam_map_vbvc(map_tab, bits_per_symbol, vlen_out, 1) self.drm_punct_vbvb_0 = drm.punct_vbvb(pp, pp_tail, (vlen_in + n_tailbits) * denom_mother_code_rate, vlen_out * 2, n_tailbits * denom_mother_code_rate) self.drm_interleaver_vbvb_0 = drm.interleaver_vbvb((interl_seq)) self.add_tailbits_vbvb_0 = drm.add_tailbits_vbvb(vlen_in, n_tailbits) ################################################## # Connections ################################################## self.connect((self, 0), (self.add_tailbits_vbvb_0, 0)) self.connect((self.add_tailbits_vbvb_0, 0), (self.gr_vector_to_stream_0, 0)) self.connect((self.gr_vector_to_stream_0, 0), (self.trellis_encoder_xx_0, 0)) self.connect((self.trellis_encoder_xx_0, 0), (self.gr_unpack_k_bits_bb_0, 0)) self.connect((self.gr_unpack_k_bits_bb_0, 0), (self.gr_stream_to_vector_0, 0)) self.connect((self.gr_stream_to_vector_0, 0), (self.drm_punct_vbvb_0, 0)) self.connect((self.drm_punct_vbvb_0, 0), (self.drm_interleaver_vbvb_0, 0)) self.connect((self.drm_interleaver_vbvb_0, 0), (self.drm_qam_map_vbvc_0, 0)) self.connect((self.drm_qam_map_vbvc_0, 0), (self, 0))
def setup_test06(self):
print "... benchmarking QPSK + 1/2 mapper"
self.mode = 2
bitspermode= [0.5,1,1.5,2,3,4,4.5,5,6]
self.nobits = (bitspermode[self.mode-1])
self.data_subcarriers = 200
self.blks = self.N*(int(10*self.nobits) +1)
self.tb = gr.top_block()
self.encoder = self._encoder = trellis.encoder_bb(self.fo,0)
self.unpack = self._unpack = blocks.unpack_k_bits_bb(2)
self.puncturing = puncture_bb(self.data_subcarriers)
bmaptrig_stream_puncturing = [1]+[0]*(10/2-1)
self.bitmap_trigger_puncturing = blocks.vector_source_b(bmaptrig_stream_puncturing, True)
if self.interleave:
int_object=trellis.interleaver(2000,666)
self.interlv = trellis.permutation(int_object.K(),int_object.INTER(),1,gr.sizeof_char)
#self.bitmap = [self.nobits]*self.data_subcarriers
self.bitmap = [self.mode]*self.data_subcarriers + [self.mode]*self.data_subcarriers
#self.bmaptrig_stream = [1, 1]+[0]*(11-2)
self.bitdata = [randint(0,1) for i in range(self.blks*self.data_subcarriers)]
#self.data = numpy.array(self.bitdata)*(-2)+1
self.src = blocks.vector_source_b(self.bitdata)
self.bitmap_src = blocks.vector_source_b(self.bitmap,True, self.data_subcarriers)
#self.bitmap_trigger = blocks.vector_source_b(self.bmaptrig_stream, True)
self.modulator = generic_mapper_bcv(self.data_subcarriers,self.coding,10)
self.snk = blocks.null_sink(gr.sizeof_gr_complex*self.data_subcarriers)
#Connect blocks
self.tb.connect(self.src,self.encoder,self.unpack,self.puncturing)#, self.modulator, self.snk)
if self.interleave:
self.tb.connect(self.puncturing,self.interlv, self.modulator, self.snk)
else:
self.tb.connect(self.puncturing, self.modulator, self.snk)
self.tb.connect(self.bitmap_src, (self.modulator,1))
self.tb.connect(self.bitmap_src, (self.puncturing,1))
self.tb.connect(self.bitmap_trigger_puncturing, (self.puncturing,2))
def __init__(self, vlen_in=1, vlen_out=1, n_tailbits=1, denom_mother_code_rate=1, map_tab=0, interl_seq=(0,1), bits_per_symbol=1, pp=0, pp_tail=0, gen_poly=0):
gr.hier_block2.__init__(
self, "MLC 4QAM",
gr.io_signature(1, 1, gr.sizeof_char*1),
gr.io_signature(1, 1, gr.sizeof_gr_complex*1),
)
##################################################
# Parameters
##################################################
self.vlen_in = vlen_in
self.vlen_out = vlen_out
self.n_tailbits = n_tailbits
self.denom_mother_code_rate = denom_mother_code_rate
self.map_tab = map_tab
self.interl_seq = interl_seq
self.bits_per_symbol = bits_per_symbol
self.pp = pp
self.pp_tail = pp_tail
self.gen_poly = gen_poly
##################################################
# Blocks
##################################################
self.trellis_encoder_xx_0 = trellis.encoder_bb(trellis.fsm(1, denom_mother_code_rate, gen_poly), 0, 0) if False else trellis.encoder_bb(trellis.fsm(1, denom_mother_code_rate, gen_poly), 0)
self.drm_qam_map_bc_0 = drm.qam_map_bc(map_tab, bits_per_symbol, vlen_out, 1)
self.drm_punct_bb_0 = drm.punct_bb(pp, pp_tail, (vlen_in + n_tailbits) * denom_mother_code_rate, vlen_out * 2, n_tailbits * denom_mother_code_rate)
self.drm_interleaver_bb_0 = drm.interleaver_bb(((interl_seq)))
self.blocks_unpack_k_bits_bb_0 = blocks.unpack_k_bits_bb(denom_mother_code_rate)
self.add_tailbits_bb_0 = drm.add_tailbits_bb(vlen_in, n_tailbits)
##################################################
# Connections
##################################################
self.connect((self.add_tailbits_bb_0, 0), (self.trellis_encoder_xx_0, 0))
self.connect((self.blocks_unpack_k_bits_bb_0, 0), (self.drm_punct_bb_0, 0))
self.connect((self.drm_interleaver_bb_0, 0), (self.drm_qam_map_bc_0, 0))
self.connect((self.drm_punct_bb_0, 0), (self.drm_interleaver_bb_0, 0))
self.connect((self.drm_qam_map_bc_0, 0), (self, 0))
self.connect((self, 0), (self.add_tailbits_bb_0, 0))
self.connect((self.trellis_encoder_xx_0, 0), (self.blocks_unpack_k_bits_bb_0, 0))
def __init__(self):
grc_wxgui.top_block_gui.__init__(self, title="Dvbs Tx")
##################################################
# Variables
##################################################
self.symbol_rate = symbol_rate = 4500000
self.samp_rate = samp_rate = symbol_rate * 2
self.rrc_taps = rrc_taps = 20
##################################################
# Blocks
##################################################
self.wxgui_fftsink2_0 = fftsink2.fft_sink_c(
self.GetWin(),
baseband_freq=435000000,
y_per_div=10,
y_divs=10,
ref_level=0,
ref_scale=2.0,
sample_rate=samp_rate,
fft_size=1024,
fft_rate=15,
average=True,
avg_alpha=None,
title="FFT Plot",
peak_hold=False,
)
self.Add(self.wxgui_fftsink2_0.win)
self.trellis_encoder_xx_0 = trellis.encoder_bb(trellis.fsm(1, 2, (0171, 0133)), 0)
self.osmosdr_sink_0 = osmosdr.sink( args="numchan=" + str(1) + " " + "" )
self.osmosdr_sink_0.set_sample_rate(samp_rate)
self.osmosdr_sink_0.set_center_freq(435e6, 0)
self.osmosdr_sink_0.set_freq_corr(0, 0)
self.osmosdr_sink_0.set_gain(2, 0)
self.osmosdr_sink_0.set_if_gain(0, 0)
self.osmosdr_sink_0.set_bb_gain(-4, 0)
self.osmosdr_sink_0.set_antenna("", 0)
self.osmosdr_sink_0.set_bandwidth(6000000, 0)
self.fft_filter_xxx_0 = filter.fft_filter_ccc(1, (firdes.root_raised_cosine(1.79, samp_rate, samp_rate/2, 0.35, rrc_taps)), 1)
self.fft_filter_xxx_0.declare_sample_delay(0)
self.dvbs_reed_solomon_enc_bb_0 = dvbs.reed_solomon_enc_bb()
self.dvbs_randomizer_bb_0 = dvbs.randomizer_bb()
self.dvbs_puncture_bb_0 = dvbs.puncture_bb(dvbs.C1_2)
self.dvbs_modulator_bc_0 = dvbs.modulator_bc()
self.dvbs_interleaver_bb_0 = dvbs.interleaver_bb()
self.blocks_unpack_k_bits_bb_0 = blocks.unpack_k_bits_bb(2)
self.blocks_packed_to_unpacked_xx_0 = blocks.packed_to_unpacked_bb(1, gr.GR_MSB_FIRST)
self.blocks_pack_k_bits_bb_0 = blocks.pack_k_bits_bb(2)
self.blocks_file_source_0 = blocks.file_source(gr.sizeof_char*1, "/home/re/xfer/adv.ts", False)
##################################################
# Connections
##################################################
self.connect((self.blocks_file_source_0, 0), (self.dvbs_randomizer_bb_0, 0))
self.connect((self.dvbs_randomizer_bb_0, 0), (self.dvbs_reed_solomon_enc_bb_0, 0))
self.connect((self.dvbs_reed_solomon_enc_bb_0, 0), (self.dvbs_interleaver_bb_0, 0))
self.connect((self.dvbs_interleaver_bb_0, 0), (self.blocks_packed_to_unpacked_xx_0, 0))
self.connect((self.blocks_packed_to_unpacked_xx_0, 0), (self.trellis_encoder_xx_0, 0))
self.connect((self.trellis_encoder_xx_0, 0), (self.blocks_unpack_k_bits_bb_0, 0))
self.connect((self.blocks_unpack_k_bits_bb_0, 0), (self.dvbs_puncture_bb_0, 0))
self.connect((self.dvbs_puncture_bb_0, 0), (self.blocks_pack_k_bits_bb_0, 0))
self.connect((self.dvbs_modulator_bc_0, 0), (self.fft_filter_xxx_0, 0))
self.connect((self.blocks_pack_k_bits_bb_0, 0), (self.dvbs_modulator_bc_0, 0))
self.connect((self.fft_filter_xxx_0, 0), (self.osmosdr_sink_0, 0))
self.connect((self.fft_filter_xxx_0, 0), (self.wxgui_fftsink2_0, 0))
def __init__(self, options):
gr.hier_block2.__init__(
self, "fbmc_transmit_path", gr.io_signature(0, 0, 0), gr.io_signature(1, 1, gr.sizeof_gr_complex)
)
common_options.defaults(options)
config = self.config = station_configuration()
config.data_subcarriers = options.subcarriers
config.cp_length = 0
config.frame_data_blocks = options.data_blocks
config._verbose = options.verbose
config.fft_length = options.fft_length
config.dc_null = options.dc_null
config.training_data = default_block_header(config.data_subcarriers, config.fft_length, config.dc_null, options)
config.coding = options.coding
config.fbmc = options.fbmc
config.adaptive_fbmc = options.adaptive_fbmc
config.frame_id_blocks = 1 # FIXME
# digital rms amplitude sent to USRP
rms_amp = options.rms_amplitude
self._options = copy.copy(options)
config.block_length = config.fft_length + config.cp_length
config.frame_data_part = config.frame_data_blocks + config.frame_id_blocks
config.frame_length = config.training_data.fbmc_no_preambles + 2 * config.frame_data_part
config.subcarriers = config.data_subcarriers + config.training_data.pilot_subcarriers
config.virtual_subcarriers = config.fft_length - config.subcarriers - config.dc_null
# default values if parameters not set
if rms_amp is None:
rms_amp = math.sqrt(config.subcarriers)
config.rms_amplitude = rms_amp
# check some bounds
if config.fft_length < config.subcarriers:
raise SystemError, "Subcarrier number must be less than FFT length"
if config.fft_length < config.cp_length:
raise SystemError, "Cyclic prefix length must be less than FFT length"
## shortcuts
blen = config.block_length
flen = config.frame_length
dsubc = config.data_subcarriers
vsubc = config.virtual_subcarriers
# Adaptive Transmitter Concept
used_id_bits = config.used_id_bits = 8 # TODO: no constant in source code
rep_id_bits = config.rep_id_bits = config.data_subcarriers / used_id_bits # BPSK
if config.data_subcarriers % used_id_bits <> 0:
raise SystemError, "Data subcarriers need to be multiple of %d" % (used_id_bits)
## Allocation Control
self.allocation_src = allocation_src(
config.data_subcarriers,
config.frame_data_blocks,
config.coding,
"tcp://*:3333",
"tcp://" + options.rx_hostname + ":3322",
)
if options.static_allocation: # DEBUG
# how many bits per subcarrier
if options.coding:
mode = 1 # Coding mode 1-9
bitspermode = [0.5, 1, 1.5, 2, 3, 4, 4.5, 5, 6] # Information bits per mode
modulbitspermode = [1, 2, 2, 4, 4, 6, 6, 6, 8] # Coding bits per mode
bitcount_vec = [(int)(config.data_subcarriers * config.frame_data_blocks * bitspermode[mode - 1])]
modul_bitcount_vec = [config.data_subcarriers * config.frame_data_blocks * modulbitspermode[mode - 1]]
bitcount_src = blocks.vector_source_i(bitcount_vec, True, 1)
modul_bitcount_src = blocks.vector_source_i(modul_bitcount_vec, True, 1)
bitloading = mode
else:
bitloading = 1
bitcount_vec = [config.data_subcarriers * config.frame_data_blocks * bitloading]
bitcount_src = blocks.vector_source_i(bitcount_vec, True, 1)
modul_bitcount_src = bitcount_src
# id's for frames
id_vec = range(0, 256)
id_src = blocks.vector_source_s(id_vec, True, 1)
# bitloading for ID symbol and then once for data symbols
# bitloading_vec = [1]*dsubc+[0]*(dsubc/2)+[2]*(dsubc/2)
test_allocation = (
[bitloading] * (int)(config.data_subcarriers / 8)
+ [0] * (int)(config.data_subcarriers / 4 * 3)
+ [bitloading] * (int)(config.data_subcarriers / 8)
)
# bitloading_vec = [1]*dsubc+[bitloading]*dsubc
bitloading_vec = [1] * dsubc + test_allocation
bitloading_src = blocks.vector_source_b(bitloading_vec, True, dsubc)
# bitcount for frames
# bitcount_vec = [config.data_subcarriers*config.frame_data_blocks*bitloading]
bitcount_vec = [config.frame_data_blocks * sum(test_allocation)]
bitcount_src = blocks.vector_source_i(bitcount_vec, True, 1)
# power loading, here same for all symbols
power_vec = (
[1] * (int)(config.data_subcarriers / 8)
+ [0] * (int)(config.data_subcarriers / 4 * 3)
+ [1] * (int)(config.data_subcarriers / 8)
)
power_src = blocks.vector_source_f(power_vec, True, dsubc)
# mux control stream to mux id and data bits
mux_vec = [0] * dsubc + [1] * bitcount_vec[0]
mux_ctrl = blocks.vector_source_b(mux_vec, True, 1)
else:
id_src = (self.allocation_src, 0)
bitcount_src = (self.allocation_src, 4)
bitloading_src = (self.allocation_src, 2)
power_src = (self.allocation_src, 1)
if options.coding:
modul_bitcount_src = (self.allocation_src, 5)
else:
modul_bitcount_src = bitcount_src
mux_ctrl = ofdm.tx_mux_ctrl(dsubc)
self.connect(modul_bitcount_src, mux_ctrl)
# Initial allocation
self.allocation_src.set_allocation([4] * config.data_subcarriers, [1] * config.data_subcarriers)
if options.benchmarking:
self.allocation_src.set_allocation([4] * config.data_subcarriers, [1] * config.data_subcarriers)
if options.lab_special_case:
self.allocation_src.set_allocation(
[0] * (config.data_subcarriers / 4)
+ [2] * (config.data_subcarriers / 2)
+ [0] * (config.data_subcarriers / 4),
[1] * config.data_subcarriers,
)
if options.log:
log_to_file(self, id_src, "data/id_src.short")
log_to_file(self, bitcount_src, "data/bitcount_src.int")
log_to_file(self, bitloading_src, "data/bitloading_src.char")
log_to_file(self, power_src, "data/power_src.cmplx")
## GUI probe output
zmq_probe_bitloading = zeromq.pub_sink(gr.sizeof_char, dsubc, "tcp://*:4445")
# also skip ID symbol bitloading with keep_one_in_n (side effect)
# factor 2 for bitloading because we have two vectors per frame, one for id symbol and one for all payload/data symbols
# factor config.frame_data_part for power because there is one vector per ofdm symbol per frame
self.connect(bitloading_src, blocks.keep_one_in_n(gr.sizeof_char * dsubc, 2 * 40), zmq_probe_bitloading)
zmq_probe_power = zeromq.pub_sink(gr.sizeof_float, dsubc, "tcp://*:4444")
# self.connect(power_src, blocks.keep_one_in_n(gr.sizeof_gr_complex*dsubc,40), blocks.complex_to_real(dsubc), zmq_probe_power)
self.connect(power_src, blocks.keep_one_in_n(gr.sizeof_float * dsubc, 40), zmq_probe_power)
## Workaround to avoid periodic structure
seed(1)
whitener_pn = [randint(0, 1) for i in range(used_id_bits * rep_id_bits)]
## ID Encoder
id_enc = self._id_encoder = repetition_encoder_sb(used_id_bits, rep_id_bits, whitener_pn)
self.connect(id_src, id_enc)
if options.log:
id_enc_f = gr.char_to_float()
self.connect(id_enc, id_enc_f)
log_to_file(self, id_enc_f, "data/id_enc_out.float")
## Reference Data Source
ber_ref_src = ber_reference_source(self._options)
self.connect(id_src, (ber_ref_src, 0))
self.connect(bitcount_src, (ber_ref_src, 1))
if options.log:
log_to_file(self, ber_ref_src, "data/ber_rec_src_tx.char")
if options.log:
log_to_file(self, btrig, "data/bitmap_trig.char")
## Frame Trigger
ftrig_stream = [1] + [0] * (config.frame_data_part - 1)
ftrig = self._frame_trigger = blocks.vector_source_b(ftrig_stream, True)
## Data Multiplexer
# Input 0: control stream
# Input 1: encoded ID stream
# Inputs 2..n: data streams
dmux = self._data_multiplexer = stream_controlled_mux_b()
self.connect(mux_ctrl, (dmux, 0))
self.connect(id_enc, (dmux, 1))
if options.coding:
fo = trellis.fsm(1, 2, [91, 121])
encoder = self._encoder = trellis.encoder_bb(fo, 0)
unpack = self._unpack = blocks.unpack_k_bits_bb(2)
self.connect(ber_ref_src, encoder, unpack)
if options.interleave:
int_object = trellis.interleaver(2000, 666)
interlv = trellis.permutation(int_object.K(), int_object.INTER(), 1, gr.sizeof_char)
if not options.nopunct:
bmaptrig_stream_puncturing = [1] + [0] * (config.frame_data_blocks / 2 - 1)
btrig_puncturing = self._bitmap_trigger_puncturing = blocks.vector_source_b(
bmaptrig_stream_puncturing, True
)
puncturing = self._puncturing = puncture_bb(config.data_subcarriers)
self.connect(bitloading_src, (puncturing, 1))
self.connect(self._bitmap_trigger_puncturing, (puncturing, 2))
self.connect(unpack, puncturing)
last_block = puncturing
if options.interleave:
self.connect(last_block, interlv)
last_block = interlv
if options.benchmarking:
self.connect(last_block, blocks.head(gr.sizeof_char, options.N), (dmux, 2))
else:
self.connect(last_block, (dmux, 2))
else:
if options.benchmarking:
self.connect(unpack, blocks.head(gr.sizeof_char, options.N), (dmux, 2))
else:
self.connect(unpack, (dmux, 2))
else:
if options.benchmarking:
self.connect(ber_ref_src, blocks.head(gr.sizeof_char, options.N), (dmux, 2))
else:
self.connect(ber_ref_src, (dmux, 2))
if options.log:
dmux_f = gr.char_to_float()
self.connect(dmux, dmux_f)
log_to_file(self, dmux_f, "data/dmux_out.float")
## Modulator
mod = self._modulator = generic_mapper_bcv(config.data_subcarriers, config.coding, config.frame_data_part)
self.connect(dmux, (mod, 0))
self.connect(bitloading_src, (mod, 1))
if options.log:
log_to_file(self, mod, "data/mod_out.compl")
modi = blocks.complex_to_imag(config.data_subcarriers)
modr = blocks.complex_to_real(config.data_subcarriers)
self.connect(mod, modi)
self.connect(mod, modr)
log_to_file(self, modi, "data/mod_imag_out.float")
log_to_file(self, modr, "data/mod_real_out.float")
## Power allocator
pa = self._power_allocator = multiply_frame_fc(config.frame_data_part, config.data_subcarriers)
self.connect(mod, (pa, 0))
self.connect(power_src, (pa, 1))
if options.log:
log_to_file(self, pa, "data/pa_out.compl")
if options.fbmc:
psubc = pa
else:
psubc = self._pilot_subcarrier_inserter = pilot_subcarrier_inserter()
self.connect(pa, psubc)
if options.log:
log_to_file(self, psubc, "data/psubc_out.compl")
subcarriers = config.subcarriers
# fbmc_pblocks_timing = self._fbmc_timing_pilot_block_inserter = fbmc_timing_pilot_block_inserter(5,False)
oqam_prep = self._oqam_prep = fbmc_oqam_preprocessing_vcvc(config.subcarriers, 0, 0)
self.connect(psubc, oqam_prep)
fbmc_pblocks = self._fbmc_pilot_block_inserter = fbmc_pilot_block_inserter(5, False)
self.connect(oqam_prep, fbmc_pblocks)
# log_to_file(self, fbmc_pblocks, "data/fbmc_pblocks_out.compl")
# fbmc_insert_pream = self._fbmc_insert_pream = fbmc_insert_preamble_vcvc(M, syms_per_frame, preamble)
# log_to_file(self, oqam_prep, "data/oqam_prep.compl")
# log_to_file(self, psubc, "data/psubc_out.compl")
# fbmc_pblocks = fbmc_pblocks_timing
# log_to_file(self, fbmc_pblocks, "data/fbmc_pblocks_out.compl")
beta_mult = self._beta_mult = fbmc_beta_multiplier_vcvc(config.subcarriers, 4, 4 * config.fft_length - 1, 0)
self.connect(fbmc_pblocks, beta_mult)
log_to_file(self, beta_mult, "data/beta_mult.compl")
## Add virtual subcarriers
if config.fft_length > subcarriers:
vsubc = self._virtual_subcarrier_extender = vector_padding_dc_null(
config.subcarriers, config.fft_length, config.dc_null
)
self.connect(beta_mult, vsubc)
else:
vsubc = self._virtual_subcarrier_extender = beta_mult
if options.log:
log_to_file(self, vsubc, "data/vsubc_out.compl")
## IFFT, no window, block shift
ifft = self._ifft = fft_blocks.fft_vcc(config.fft_length, False, [], True)
self.connect(vsubc, ifft)
if options.log:
log_to_file(self, ifft, "data/ifft_out.compl")
# FBMC separate stream + filterbanks
separate_oqam = self._separate_oqam = fbmc_separate_vcvc(config.fft_length, 2)
poly_netw_1 = self._poly_netw_1 = fbmc_polyphase_network_vcvc(
config.fft_length, 4, 4 * config.fft_length - 1, False
)
poly_netw_2 = self._poly_netw_2 = fbmc_polyphase_network_vcvc(
config.fft_length, 4, 4 * config.fft_length - 1, False
)
overlap_p2s = self._overlap_p2s = fbmc_overlapping_parallel_to_serial_vcc(config.fft_length)
self.connect(ifft, (separate_oqam, 0), poly_netw_1)
self.connect((separate_oqam, 1), poly_netw_2)
self.connect(poly_netw_1, (overlap_p2s, 0))
self.connect(poly_netw_2, (overlap_p2s, 1))
## Pilot blocks (preambles)
# pblocks = self._pilot_block_inserter = pilot_block_inserter2(5,False)
# self.connect( overlap_p2s, blocks.stream_to_vector(gr.sizeof_gr_complex,config.fft_length/2), pblocks )
# log_to_file(self, pblocks, "data/fbmc_pilot_block_ins_out.compl")
if options.log:
log_to_file(self, pblocks, "data/pilot_block_ins_out.compl")
## Cyclic Prefix
# cp = self._cyclic_prefixer = cyclic_prefixer(config.fft_length,
# config.block_length)
# cp= blocks.vector_to_stream(gr.sizeof_gr_complex, config.fft_length/2)
# self.connect(pblocks, cp )
# self.connect( overlap_p2s,blocks.stream_to_vector(gr.sizeof_gr_complex,config.fft_length/2), cp )
lastblock = overlap_p2s
if options.log:
log_to_file(self, overlap_p2s, "data/overlap_p2s_out.compl")
# Digital Amplifier for resource allocation
if config.adaptive_fbmc:
rep = blocks.repeat(gr.sizeof_gr_complex, config.frame_length * config.block_length)
amp = blocks.multiply_cc()
self.connect(lastblock, (amp, 0))
self.connect((self.allocation_src, 3), rep, (amp, 1))
lastblock = amp
else:
self.connect((self.allocation_src, 3), blocks.null_sink(gr.sizeof_gr_complex))
## Digital Amplifier
# amp = self._amplifier = gr.multiply_const_cc(1)
amp = self._amplifier = ofdm.multiply_const_ccf(1.0)
self.connect(lastblock, amp)
self.set_rms_amplitude(rms_amp)
# log_to_file(self, amp, "data/amp_tx_out.compl")
if options.log:
log_to_file(self, amp, "data/amp_tx_out.compl")
## Tx parameters
bandwidth = options.bandwidth or 2e6
bits = 8 * config.data_subcarriers * config.frame_data_blocks # max. QAM256
samples_per_frame = config.frame_length * config.block_length
tb = samples_per_frame / bandwidth
# set dummy carrier frequency if none available due to baseband mode
if options.tx_freq is None:
options.tx_freq = 0.0
self.tx_parameters = {
"carrier_frequency": options.tx_freq / 1e9,
"fft_size": config.fft_length,
"cp_size": config.cp_length,
"subcarrier_spacing": options.bandwidth / config.fft_length / 1e3,
"data_subcarriers": config.data_subcarriers,
"bandwidth": options.bandwidth / 1e6,
"frame_length": config.frame_length,
"symbol_time": (config.cp_length + config.fft_length) / options.bandwidth * 1e6,
"max_data_rate": (bits / tb) / 1e6,
}
## Setup Output
self.connect(amp, self)
# Display some information about the setup
if config._verbose:
self._print_verbage()
def __init__(self, vlen_in=1, vlen_out=1, n_tailbits=6, denom_mother_code_rate=6, gen_poly=(91, 121, 101, 91, 121, 101), bits_per_symbol=0, N=1, map_tab=0, pp_0=0, pp_0_tail=0, pp_1=0, pp_1_tail=0, part_len_top=1, part_len_bot=1, M_total=0, interl_seq_0_2=range(2), interl_seq_1_2=range(2)): gr.hier_block2.__init__( self, "DRM MLC 16-QAM", gr.io_signature(1, 1, gr.sizeof_char*vlen_in), gr.io_signature(1, 1, gr.sizeof_gr_complex*vlen_out), ) ################################################## # Parameters ################################################## self.vlen_in = vlen_in self.vlen_out = vlen_out self.n_tailbits = n_tailbits self.denom_mother_code_rate = denom_mother_code_rate self.gen_poly = gen_poly self.bits_per_symbol = bits_per_symbol self.N = N self.map_tab = map_tab self.pp_0 = pp_0 self.pp_0_tail = pp_0_tail self.pp_1 = pp_1 self.pp_1_tail = pp_1_tail self.part_len_top = part_len_top self.part_len_bot = part_len_bot self.M_total = M_total self.interl_seq_0_2 = interl_seq_0_2 self.interl_seq_1_2 = interl_seq_1_2 ################################################## # Blocks ################################################## self.trellis_encoder_xx_0_0 = trellis.encoder_bb(trellis.fsm(1, denom_mother_code_rate, gen_poly), 0) self.trellis_encoder_xx_0 = trellis.encoder_bb(trellis.fsm(1, denom_mother_code_rate, gen_poly), 0) self.gr_vector_to_stream_1_0 = gr.vector_to_stream(gr.sizeof_char*1, part_len_bot+ n_tailbits) self.gr_vector_to_stream_1 = gr.vector_to_stream(gr.sizeof_char*1, part_len_top + n_tailbits) self.gr_unpack_k_bits_bb_0_0 = gr.unpack_k_bits_bb(denom_mother_code_rate) self.gr_unpack_k_bits_bb_0 = gr.unpack_k_bits_bb(denom_mother_code_rate) self.gr_stream_to_vector_0_0 = gr.stream_to_vector(gr.sizeof_char*1, (part_len_bot + n_tailbits) * denom_mother_code_rate) self.gr_stream_to_vector_0 = gr.stream_to_vector(gr.sizeof_char*1, (part_len_top + n_tailbits) * denom_mother_code_rate) self.drm_qam_map_vbvc_0 = drm.qam_map_vbvc(map_tab, bits_per_symbol, vlen_out, 2) self.drm_punct_vbvb_0_0 = drm.punct_vbvb(pp_1, pp_1_tail, (part_len_bot + n_tailbits) * denom_mother_code_rate, vlen_out * 2, n_tailbits * denom_mother_code_rate) self.drm_punct_vbvb_0 = drm.punct_vbvb(pp_0, pp_0_tail, (part_len_top + n_tailbits) * denom_mother_code_rate, vlen_out * 2, n_tailbits * denom_mother_code_rate) self.drm_partitioning_16_vbvb_0 = drm.partitioning_vbvb(vlen_in, M_total) self.drm_interleaver_vbvb_0_0 = drm.interleaver_vbvb((interl_seq_1_2)) self.drm_interleaver_vbvb_0 = drm.interleaver_vbvb((interl_seq_0_2)) self.add_tailbits_vbvb_0_0 = drm.add_tailbits_vbvb(part_len_bot, n_tailbits) self.add_tailbits_vbvb_0 = drm.add_tailbits_vbvb(part_len_top, n_tailbits) ################################################## # Connections ################################################## self.connect((self.drm_interleaver_vbvb_0, 0), (self.drm_qam_map_vbvc_0, 0)) self.connect((self.drm_qam_map_vbvc_0, 0), (self, 0)) self.connect((self.trellis_encoder_xx_0, 0), (self.gr_unpack_k_bits_bb_0, 0)) self.connect((self.gr_stream_to_vector_0, 0), (self.drm_punct_vbvb_0, 0)) self.connect((self.gr_unpack_k_bits_bb_0, 0), (self.gr_stream_to_vector_0, 0)) self.connect((self.drm_punct_vbvb_0, 0), (self.drm_interleaver_vbvb_0, 0)) self.connect((self.trellis_encoder_xx_0_0, 0), (self.gr_unpack_k_bits_bb_0_0, 0)) self.connect((self.gr_stream_to_vector_0_0, 0), (self.drm_punct_vbvb_0_0, 0)) self.connect((self.gr_unpack_k_bits_bb_0_0, 0), (self.gr_stream_to_vector_0_0, 0)) self.connect((self.drm_punct_vbvb_0_0, 0), (self.drm_interleaver_vbvb_0_0, 0)) self.connect((self.drm_interleaver_vbvb_0_0, 0), (self.drm_qam_map_vbvc_0, 1)) self.connect((self.gr_vector_to_stream_1, 0), (self.trellis_encoder_xx_0, 0)) self.connect((self.gr_vector_to_stream_1_0, 0), (self.trellis_encoder_xx_0_0, 0)) self.connect((self, 0), (self.drm_partitioning_16_vbvb_0, 0)) self.connect((self.add_tailbits_vbvb_0, 0), (self.gr_vector_to_stream_1, 0)) self.connect((self.add_tailbits_vbvb_0_0, 0), (self.gr_vector_to_stream_1_0, 0)) self.connect((self.drm_partitioning_16_vbvb_0, 1), (self.add_tailbits_vbvb_0_0, 0)) self.connect((self.drm_partitioning_16_vbvb_0, 0), (self.add_tailbits_vbvb_0, 0))
def __init__(self, options):
gr.hier_block2.__init__(self, "transmit_path",
gr.io_signature(0, 0, 0),
gr.io_signature(1, 1, gr.sizeof_gr_complex))
common_options.defaults(options)
config = self.config = station_configuration()
config.data_subcarriers = options.subcarriers
config.cp_length = options.cp_length
config.frame_data_blocks = options.data_blocks
config._verbose = options.verbose
config.fft_length = options.fft_length
config.dc_null = options.dc_null
config.training_data = default_block_header(config.data_subcarriers,
config.fft_length,
config.dc_null, options)
config.coding = options.coding
config.bandwidth = options.bandwidth
config.gui_frame_rate = options.gui_frame_rate
config.fbmc = options.fbmc
config.frame_id_blocks = 1 # FIXME
# digital rms amplitude sent to USRP
rms_amp = options.rms_amplitude
self._options = copy.copy(options)
config.block_length = config.fft_length + config.cp_length
config.frame_data_part = config.frame_data_blocks + config.frame_id_blocks
config.frame_length = config.frame_data_part + \
config.training_data.no_pilotsyms
config.subcarriers = config.data_subcarriers + \
config.training_data.pilot_subcarriers
config.virtual_subcarriers = config.fft_length - config.subcarriers - config.dc_null
# default values if parameters not set
if rms_amp is None:
rms_amp = math.sqrt(config.subcarriers)
config.rms_amplitude = rms_amp
# check some bounds
if config.fft_length < config.subcarriers:
raise SystemError, "Subcarrier number must be less than FFT length"
if config.fft_length < config.cp_length:
raise SystemError, "Cyclic prefix length must be less than FFT length"
## shortcuts
blen = config.block_length
flen = config.frame_length
dsubc = config.data_subcarriers
vsubc = config.virtual_subcarriers
# Adaptive Transmitter Concept
used_id_bits = config.used_id_bits = 8 #TODO: no constant in source code
rep_id_bits = config.rep_id_bits = config.data_subcarriers / used_id_bits #BPSK
if config.data_subcarriers % used_id_bits <> 0:
raise SystemError, "Data subcarriers need to be multiple of %d" % (
used_id_bits)
# adapt OFDM frame rate and GUI display frame rate
self.keep_frame_n = int(
1.0 / (config.frame_length *
(config.cp_length + config.fft_length) / config.bandwidth) /
config.gui_frame_rate)
## Allocation Control
self.allocation_src = allocation_src(
config.data_subcarriers, config.frame_data_blocks, config.coding,
"tcp://*:3333", "tcp://" + options.rx_hostname + ":3322")
if options.static_allocation: #DEBUG
# how many bits per subcarrier
if options.coding:
mode = 1 # Coding mode 1-9
bitspermode = [0.5, 1, 1.5, 2, 3, 4, 4.5, 5,
6] # Information bits per mode
modulbitspermode = [1, 2, 2, 4, 4, 6, 6, 6,
8] # Coding bits per mode
bitcount_vec = [
(int)(config.data_subcarriers * config.frame_data_blocks *
bitspermode[mode - 1])
]
modul_bitcount_vec = [
config.data_subcarriers * config.frame_data_blocks *
modulbitspermode[mode - 1]
]
bitcount_src = blocks.vector_source_i(bitcount_vec, True, 1)
modul_bitcount_src = blocks.vector_source_i(
modul_bitcount_vec, True, 1)
bitloading = mode
else:
bitloading = 1
bitcount_vec = [
config.data_subcarriers * config.frame_data_blocks *
bitloading
]
bitcount_src = blocks.vector_source_i(bitcount_vec, True, 1)
modul_bitcount_src = bitcount_src
# id's for frames
id_vec = range(0, 256)
id_src = blocks.vector_source_s(id_vec, True, 1)
# bitloading for ID symbol and then once for data symbols
#bitloading_vec = [1]*dsubc+[0]*(dsubc/2)+[2]*(dsubc/2)
#test_allocation = [bitloading]*(int)(config.data_subcarriers/8)+ [0]*(int)(config.data_subcarriers/4*3) + [bitloading]*(int)(config.data_subcarriers/8)
#bitloading_vec = [1]*dsubc+[bitloading]*dsubc
test_allocation = [bitloading] * dsubc
bitloading_vec = [bitloading] * dsubc + test_allocation
bitloading_src = blocks.vector_source_b(bitloading_vec, True,
dsubc)
# bitcount for frames
#bitcount_vec = [config.data_subcarriers*config.frame_data_blocks*bitloading]
bitcount_vec = [config.frame_data_blocks * sum(test_allocation)]
bitcount_src = blocks.vector_source_i(bitcount_vec, True, 1)
# power loading, here same for all symbols
#power_vec = [1]*(int)(config.data_subcarriers/8)+ [0]*(int)(config.data_subcarriers/4*3) + [1]*(int)(config.data_subcarriers/8)
power_vec = [1] * config.data_subcarriers
power_src = blocks.vector_source_f(power_vec, True, dsubc)
# mux control stream to mux id and data bits
mux_vec = [0] * dsubc + [1] * bitcount_vec[0]
mux_ctrl = blocks.vector_source_b(mux_vec, True, 1)
else:
id_src = (self.allocation_src, 0)
bitcount_src = (self.allocation_src, 4)
bitloading_src = (self.allocation_src, 2)
power_src = (self.allocation_src, 1)
if options.coding:
modul_bitcount_src = (self.allocation_src, 5)
else:
modul_bitcount_src = bitcount_src
mux_ctrl = ofdm.tx_mux_ctrl(dsubc)
self.connect(modul_bitcount_src, mux_ctrl)
#Initial allocation
self.allocation_src.set_allocation([2] * config.data_subcarriers,
[1] * config.data_subcarriers)
self.allocation_src.set_allocation_scheme(0)
if options.benchmarking:
self.allocation_src.set_allocation(
[4] * config.data_subcarriers,
[1] * config.data_subcarriers)
if options.lab_special_case:
self.allocation_src.set_allocation(
[0] * (config.data_subcarriers / 4) + [2] *
(config.data_subcarriers / 2) + [0] *
(config.data_subcarriers / 4), [1] * config.data_subcarriers)
if options.log:
log_to_file(self, id_src, "data/id_src.short")
log_to_file(self, bitcount_src, "data/bitcount_src.int")
log_to_file(self, bitloading_src, "data/bitloading_src.char")
log_to_file(self, power_src, "data/power_src.cmplx")
## GUI probe output
zmq_probe_bitloading = zeromq.pub_sink(gr.sizeof_char, dsubc,
"tcp://*:4445")
# also skip ID symbol bitloading with keep_one_in_n (side effect)
# factor 2 for bitloading because we have two vectors per frame, one for id symbol and one for all payload/data symbols
# factor config.frame_data_part for power because there is one vector per ofdm symbol per frame
self.connect(bitloading_src,
blocks.keep_one_in_n(gr.sizeof_char * dsubc, 2 * 40),
zmq_probe_bitloading)
zmq_probe_power = zeromq.pub_sink(gr.sizeof_float, dsubc,
"tcp://*:4444")
#self.connect(power_src, blocks.keep_one_in_n(gr.sizeof_gr_complex*dsubc,40), blocks.complex_to_real(dsubc), zmq_probe_power)
self.connect(power_src,
blocks.keep_one_in_n(gr.sizeof_float * dsubc, 40),
zmq_probe_power)
## Workaround to avoid periodic structure
seed(1)
whitener_pn = [
randint(0, 1) for i in range(used_id_bits * rep_id_bits)
]
## ID Encoder
id_enc = self._id_encoder = repetition_encoder_sb(
used_id_bits, rep_id_bits, whitener_pn)
self.connect(id_src, id_enc)
if options.log:
id_enc_f = gr.char_to_float()
self.connect(id_enc, id_enc_f)
log_to_file(self, id_enc_f, "data/id_enc_out.float")
## Reference Data Source
ber_ref_src = ber_reference_source(self._options)
self.connect(id_src, (ber_ref_src, 0))
self.connect(bitcount_src, (ber_ref_src, 1))
if options.log:
log_to_file(self, ber_ref_src, "data/ber_rec_src_tx.char")
if options.log:
log_to_file(self, btrig, "data/bitmap_trig.char")
## Bitmap Update Trigger for puncturing
if not options.nopunct:
bmaptrig_stream_puncturing = [
1
] + [0] * (config.frame_data_blocks / 2 - 1)
btrig_puncturing = self._bitmap_trigger_puncturing = blocks.vector_source_b(
bmaptrig_stream_puncturing, True)
bmapsrc_stream_puncturing = [1] * dsubc + [2] * dsubc
bsrc_puncturing = self._bitmap_src_puncturing = blocks.vector_source_b(
bmapsrc_stream_puncturing, True, dsubc)
if options.log and options.coding and not options.nopunct:
log_to_file(self, btrig_puncturing,
"data/bitmap_trig_puncturing.char")
## Frame Trigger
ftrig_stream = [1] + [0] * (config.frame_data_part - 1)
ftrig = self._frame_trigger = blocks.vector_source_b(
ftrig_stream, True)
## Data Multiplexer
# Input 0: control stream
# Input 1: encoded ID stream
# Inputs 2..n: data streams
dmux = self._data_multiplexer = stream_controlled_mux_b()
self.connect(mux_ctrl, (dmux, 0))
self.connect(id_enc, (dmux, 1))
if options.coding:
fo = trellis.fsm(1, 2, [91, 121])
encoder = self._encoder = trellis.encoder_bb(fo, 0)
unpack = self._unpack = blocks.unpack_k_bits_bb(2)
self.connect(ber_ref_src, encoder, unpack)
if options.interleave:
int_object = trellis.interleaver(2000, 666)
interlv = trellis.permutation(int_object.K(),
int_object.INTER(), 1,
gr.sizeof_char)
if not options.nopunct:
bmaptrig_stream_puncturing = [
1
] + [0] * (config.frame_data_blocks / 2 - 1)
btrig_puncturing = self._bitmap_trigger_puncturing = blocks.vector_source_b(
bmaptrig_stream_puncturing, True)
puncturing = self._puncturing = puncture_bb(
config.data_subcarriers)
self.connect(bitloading_src, (puncturing, 1))
self.connect(self._bitmap_trigger_puncturing, (puncturing, 2))
self.connect(unpack, puncturing)
last_block = puncturing
if options.interleave:
self.connect(last_block, interlv)
last_block = interlv
if options.benchmarking:
self.connect(last_block,
blocks.head(gr.sizeof_char, options.N),
(dmux, 2))
else:
self.connect(last_block, (dmux, 2))
else:
if options.benchmarking:
self.connect(unpack, blocks.head(gr.sizeof_char,
options.N), (dmux, 2))
else:
self.connect(unpack, (dmux, 2))
else:
if options.benchmarking:
self.connect(ber_ref_src, blocks.head(gr.sizeof_char,
options.N), (dmux, 2))
else:
self.connect(ber_ref_src, (dmux, 2))
if options.log:
dmux_f = gr.char_to_float()
self.connect(dmux, dmux_f)
log_to_file(self, dmux_f, "data/dmux_out.float")
## Modulator
mod = self._modulator = generic_mapper_bcv(config.data_subcarriers,
config.coding,
config.frame_data_part)
self.connect(dmux, (mod, 0))
self.connect(bitloading_src, (mod, 1))
#log_to_file(self, mod, "data/mod_out.compl")
if options.log:
log_to_file(self, mod, "data/mod_out.compl")
modi = blocks.complex_to_imag(config.data_subcarriers)
modr = blocks.complex_to_real(config.data_subcarriers)
self.connect(mod, modi)
self.connect(mod, modr)
log_to_file(self, modi, "data/mod_imag_out.float")
log_to_file(self, modr, "data/mod_real_out.float")
## Power allocator
pa = self._power_allocator = multiply_frame_fc(config.frame_data_part,
config.data_subcarriers)
self.connect(mod, (pa, 0))
self.connect(power_src, (pa, 1))
if options.log:
log_to_file(self, pa, "data/pa_out.compl")
# Standard Transmitter Parts
## Pilot subcarriers
psubc = self._pilot_subcarrier_inserter = pilot_subcarrier_inserter()
self.connect(pa, psubc)
if options.log:
log_to_file(self, psubc, "data/psubc_out.compl")
## Add virtual subcarriers
if config.fft_length > config.subcarriers:
vsubc = self._virtual_subcarrier_extender = \
vector_padding_dc_null(config.subcarriers, config.fft_length,config.dc_null)
self.connect(psubc, vsubc)
else:
vsubc = self._virtual_subcarrier_extender = psubc
if options.log:
log_to_file(self, vsubc, "data/vsubc_out.compl")
## IFFT, no window, block shift
ifft = self._ifft = fft_blocks.fft_vcc(config.fft_length, False, [],
True)
self.connect(vsubc, ifft)
if options.log:
log_to_file(self, ifft, "data/ifft_out.compl")
## Pilot blocks (preambles)
pblocks = self._pilot_block_inserter = pilot_block_inserter(5, False)
self.connect(ifft, pblocks)
if options.log:
log_to_file(self, pblocks, "data/pilot_block_ins_out.compl")
## Cyclic Prefix
cp = self._cyclic_prefixer = cyclic_prefixer(config.fft_length,
config.block_length)
self.connect(pblocks, cp)
lastblock = cp
if options.log:
log_to_file(self, cp, "data/cp_out.compl")
#Digital Amplifier for resource allocation
#if not options.coding:
rep = blocks.repeat(gr.sizeof_gr_complex,
config.frame_length * config.block_length)
amp = blocks.multiply_cc()
self.connect(lastblock, (amp, 0))
self.connect((self.allocation_src, 3), rep, (amp, 1))
lastblock = amp
## Digital Amplifier
#amp = self._amplifier = gr.multiply_const_cc(1)
amp = self._amplifier = ofdm.multiply_const_ccf(1.0)
self.connect(lastblock, amp)
self.set_rms_amplitude(rms_amp)
if options.log:
log_to_file(self, amp, "data/amp_tx_out.compl")
## Tx parameters
bandwidth = options.bandwidth or 2e6
bits = 8 * config.data_subcarriers * config.frame_data_blocks # max. QAM256
samples_per_frame = config.frame_length * config.block_length
tb = samples_per_frame / bandwidth
# set dummy carrier frequency if none available due to baseband mode
if (options.tx_freq is None):
options.tx_freq = 0.0
self.tx_parameters = {'carrier_frequency':options.tx_freq/1e9,'fft_size':config.fft_length, 'cp_size':config.cp_length \
, 'subcarrier_spacing':options.bandwidth/config.fft_length/1e3 \
, 'data_subcarriers':config.data_subcarriers, 'bandwidth':options.bandwidth/1e6 \
, 'frame_length':config.frame_length \
, 'symbol_time':(config.cp_length + config.fft_length)/options.bandwidth*1e6, 'max_data_rate':(bits/tb)/1e6}
## Setup Output
self.connect(amp, self)
# Display some information about the setup
if config._verbose:
self._print_verbage()
def __init__(self,
vlen_in=1,
vlen_out=1,
n_tailbits=6,
denom_mother_code_rate=6,
gen_poly=(91, 121, 101, 91, 121, 101),
bits_per_symbol=0,
N=1,
map_tab=0,
pp_0=0,
pp_0_tail=0,
pp_1=0,
pp_1_tail=0,
part_len_top=1,
part_len_bot=1,
M_total=0,
interl_seq_0_2=range(2),
interl_seq_1_2=range(2)):
gr.hier_block2.__init__(
self,
"DRM MLC 16-QAM",
gr.io_signature(1, 1, gr.sizeof_char * vlen_in),
gr.io_signature(1, 1, gr.sizeof_gr_complex * vlen_out),
)
##################################################
# Parameters
##################################################
self.vlen_in = vlen_in
self.vlen_out = vlen_out
self.n_tailbits = n_tailbits
self.denom_mother_code_rate = denom_mother_code_rate
self.gen_poly = gen_poly
self.bits_per_symbol = bits_per_symbol
self.N = N
self.map_tab = map_tab
self.pp_0 = pp_0
self.pp_0_tail = pp_0_tail
self.pp_1 = pp_1
self.pp_1_tail = pp_1_tail
self.part_len_top = part_len_top
self.part_len_bot = part_len_bot
self.M_total = M_total
self.interl_seq_0_2 = interl_seq_0_2
self.interl_seq_1_2 = interl_seq_1_2
##################################################
# Blocks
##################################################
self.trellis_encoder_xx_0_0 = trellis.encoder_bb(
trellis.fsm(1, denom_mother_code_rate, gen_poly), 0)
self.trellis_encoder_xx_0 = trellis.encoder_bb(
trellis.fsm(1, denom_mother_code_rate, gen_poly), 0)
self.gr_vector_to_stream_1_0 = gr.vector_to_stream(
gr.sizeof_char * 1, part_len_bot + n_tailbits)
self.gr_vector_to_stream_1 = gr.vector_to_stream(
gr.sizeof_char * 1, part_len_top + n_tailbits)
self.gr_unpack_k_bits_bb_0_0 = gr.unpack_k_bits_bb(
denom_mother_code_rate)
self.gr_unpack_k_bits_bb_0 = gr.unpack_k_bits_bb(
denom_mother_code_rate)
self.gr_stream_to_vector_0_0 = gr.stream_to_vector(
gr.sizeof_char * 1,
(part_len_bot + n_tailbits) * denom_mother_code_rate)
self.gr_stream_to_vector_0 = gr.stream_to_vector(
gr.sizeof_char * 1,
(part_len_top + n_tailbits) * denom_mother_code_rate)
self.drm_qam_map_vbvb_0 = drm.qam_map_vbvb(map_tab, bits_per_symbol,
vlen_out, 2)
self.drm_punct_vbvb_0_0 = drm.punct_vbvb(
pp_1, pp_1_tail,
(part_len_bot + n_tailbits) * denom_mother_code_rate, vlen_out * 2,
n_tailbits * denom_mother_code_rate)
self.drm_punct_vbvb_0 = drm.punct_vbvb(
pp_0, pp_0_tail,
(part_len_top + n_tailbits) * denom_mother_code_rate, vlen_out * 2,
n_tailbits * denom_mother_code_rate)
self.drm_partitioning_16_vbvb_0 = drm.partitioning_vbvb(
vlen_in, M_total)
self.drm_interleaver_vbvb_0_0 = drm.interleaver_vbvb((interl_seq_1_2))
self.drm_interleaver_vbvb_0 = drm.interleaver_vbvb((interl_seq_0_2))
self.add_tailbits_vbvb_0_0 = drm.add_tailbits_vbvb(
part_len_bot, n_tailbits)
self.add_tailbits_vbvb_0 = drm.add_tailbits_vbvb(
part_len_top, n_tailbits)
##################################################
# Connections
##################################################
self.connect((self.drm_interleaver_vbvb_0, 0),
(self.drm_qam_map_vbvb_0, 0))
self.connect((self.drm_qam_map_vbvb_0, 0), (self, 0))
self.connect((self.trellis_encoder_xx_0, 0),
(self.gr_unpack_k_bits_bb_0, 0))
self.connect((self.gr_stream_to_vector_0, 0),
(self.drm_punct_vbvb_0, 0))
self.connect((self.gr_unpack_k_bits_bb_0, 0),
(self.gr_stream_to_vector_0, 0))
self.connect((self.drm_punct_vbvb_0, 0),
(self.drm_interleaver_vbvb_0, 0))
self.connect((self.trellis_encoder_xx_0_0, 0),
(self.gr_unpack_k_bits_bb_0_0, 0))
self.connect((self.gr_stream_to_vector_0_0, 0),
(self.drm_punct_vbvb_0_0, 0))
self.connect((self.gr_unpack_k_bits_bb_0_0, 0),
(self.gr_stream_to_vector_0_0, 0))
self.connect((self.drm_punct_vbvb_0_0, 0),
(self.drm_interleaver_vbvb_0_0, 0))
self.connect((self.drm_interleaver_vbvb_0_0, 0),
(self.drm_qam_map_vbvb_0, 1))
self.connect((self.gr_vector_to_stream_1, 0),
(self.trellis_encoder_xx_0, 0))
self.connect((self.gr_vector_to_stream_1_0, 0),
(self.trellis_encoder_xx_0_0, 0))
self.connect((self, 0), (self.drm_partitioning_16_vbvb_0, 0))
self.connect((self.add_tailbits_vbvb_0, 0),
(self.gr_vector_to_stream_1, 0))
self.connect((self.add_tailbits_vbvb_0_0, 0),
(self.gr_vector_to_stream_1_0, 0))
self.connect((self.drm_partitioning_16_vbvb_0, 1),
(self.add_tailbits_vbvb_0_0, 0))
self.connect((self.drm_partitioning_16_vbvb_0, 0),
(self.add_tailbits_vbvb_0, 0))
def __init__(self,
vlen_in=1,
vlen_out=1,
n_tailbits=1,
denom_mother_code_rate=1,
map_tab=0,
interl_seq=(0, 1),
bits_per_symbol=1,
pp=0,
pp_tail=0,
gen_poly=0):
gr.hier_block2.__init__(
self,
"MLC 4QAM",
gr.io_signature(1, 1, gr.sizeof_char * 1),
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
)
##################################################
# Parameters
##################################################
self.vlen_in = vlen_in
self.vlen_out = vlen_out
self.n_tailbits = n_tailbits
self.denom_mother_code_rate = denom_mother_code_rate
self.map_tab = map_tab
self.interl_seq = interl_seq
self.bits_per_symbol = bits_per_symbol
self.pp = pp
self.pp_tail = pp_tail
self.gen_poly = gen_poly
##################################################
# Blocks
##################################################
self.trellis_encoder_xx_0 = trellis.encoder_bb(
trellis.fsm(1, denom_mother_code_rate, gen_poly), 0,
0) if False else trellis.encoder_bb(
trellis.fsm(1, denom_mother_code_rate, gen_poly), 0)
self.drm_qam_map_bc_0 = drm.qam_map_bc(map_tab, bits_per_symbol,
vlen_out, 1)
self.drm_punct_bb_0 = drm.punct_bb(
pp, pp_tail, (vlen_in + n_tailbits) * denom_mother_code_rate,
vlen_out * 2, n_tailbits * denom_mother_code_rate)
self.drm_interleaver_bb_0 = drm.interleaver_bb(((interl_seq)))
self.blocks_unpack_k_bits_bb_0 = blocks.unpack_k_bits_bb(
denom_mother_code_rate)
self.add_tailbits_bb_0 = drm.add_tailbits_bb(vlen_in, n_tailbits)
##################################################
# Connections
##################################################
self.connect((self.add_tailbits_bb_0, 0),
(self.trellis_encoder_xx_0, 0))
self.connect((self.blocks_unpack_k_bits_bb_0, 0),
(self.drm_punct_bb_0, 0))
self.connect((self.drm_interleaver_bb_0, 0),
(self.drm_qam_map_bc_0, 0))
self.connect((self.drm_punct_bb_0, 0), (self.drm_interleaver_bb_0, 0))
self.connect((self.drm_qam_map_bc_0, 0), (self, 0))
self.connect((self, 0), (self.add_tailbits_bb_0, 0))
self.connect((self.trellis_encoder_xx_0, 0),
(self.blocks_unpack_k_bits_bb_0, 0))
def setUp (self):
self.tb = gr.top_block ()
# Read in successfully decoded live data from Matlab
linf=open('/home/demel/exchange/matlab_d.txt')
lintu=range(120)
for i in lintu:
lintu[i]=float(linf.readline())
#print lintu
# source for live data
self.srcl = gr.vector_source_f(lintu,False,120)
# Read in .txt file with example MIB encoded + CRC checksum
inf=open('/home/demel/exchange/crc.txt')
self.intu=range(40)
for i in self.intu:
self.intu[i]=float(inf.readline())
#inf=open('/home/demel/exchange/matlab_d.txt')
#intu=range(120)
#for i in range(120):
# intu[i]=float(inf.readline())
# Source and conversions
self.src = gr.vector_source_f(self.intu,False,40)
self.conv = gr.float_to_char(40,1)
# Resize vector with repetition of last part
# Vector to stream for encoder
my_map1=range(46)
for i in range(40):
my_map1[i+6]=i
for i in range(6):
my_map1[i]=i+40
self.map1 = lte.vector_resize_vbvb(my_map1,40,46)
self.vtos = gr.vector_to_stream(1*gr.sizeof_char,46)
# Encoding of input data
self.fsm = trellis.fsm(1,3,[91,121,117])
self.enc = trellis.encoder_bb(self.fsm,0)
# unpack packed bits from encoder
self.unp = gr.unpack_k_bits_bb(3)
# stream to vector
self.stov = gr.stream_to_vector(1*gr.sizeof_char,138)
# Remove first part which contains tail-biting init stuff
map2 = range(120)
for i in map2:
map2[i]= i+18
self.map2 = lte.vector_resize_vbvb(map2,138,120)
# conversion from char to float to match input of decoder
self.conv2= gr.char_to_float(120,1)
###############################################
# From here on only "receiver side" processing
###############################################
# like QPSK demodulation: NRZ coding.
vec2=range(120)
for i in vec2:
vec2[i]=float(-2.0)
self.mult = gr.multiply_const_vff(vec2)
vec=range(120)
for i in vec:
vec[i]=1
self.add = gr.add_const_vff(vec)
# this is the actual unit under test
self.vit = lte.viterbi_vfvb()
# Sinks
self.snk = gr.vector_sink_b(40)
self.snk2 = gr.vector_sink_f(120)
# connecting blocks
self.tb.connect(self.src,self.conv,self.map1,self.vtos,self.enc,self.unp)
self.tb.connect(self.unp,self.stov,self.map2,self.conv2)
self.tb.connect(self.conv2,self.mult,self.add)
self.tb.connect(self.srcl,self.vit,self.snk)
self.tb.connect(self.add,self.snk2)
def __init__(self):
gr.top_block.__init__(self, "Superposition Coding")
Qt.QWidget.__init__(self)
self.setWindowTitle("Superposition Coding")
qtgui.util.check_set_qss()
try:
self.setWindowIcon(Qt.QIcon.fromTheme('gnuradio-grc'))
except:
pass
self.top_scroll_layout = Qt.QVBoxLayout()
self.setLayout(self.top_scroll_layout)
self.top_scroll = Qt.QScrollArea()
self.top_scroll.setFrameStyle(Qt.QFrame.NoFrame)
self.top_scroll_layout.addWidget(self.top_scroll)
self.top_scroll.setWidgetResizable(True)
self.top_widget = Qt.QWidget()
self.top_scroll.setWidget(self.top_widget)
self.top_layout = Qt.QVBoxLayout(self.top_widget)
self.top_grid_layout = Qt.QGridLayout()
self.top_layout.addLayout(self.top_grid_layout)
self.settings = Qt.QSettings("GNU Radio", "int_cancellation")
self.restoreGeometry(self.settings.value("geometry").toByteArray())
##################################################
# Variables
##################################################
self.snr_db = snr_db = 10
self.prefix = prefix = "/opt/usr/src/gnuradio/gr-trellis/examples/python/fsm_files/"
self.noisevar = noisevar = 10**(-snr_db / 10)
self.fsm2 = fsm2 = "awgn1o2_4.fsm"
self.fsm1 = fsm1 = "awgn1o2_16.fsm"
self.block = block = 1000
self.alpha = alpha = .1
self.R = R = 100e3
##################################################
# Blocks
##################################################
self._alpha_range = Range(0, 1, .01, .1, 200)
self._alpha_win = RangeWidget(self._alpha_range, self.set_alpha,
'P1/P', "counter_slider", float)
self.top_layout.addWidget(self._alpha_win)
self.trellis_viterbi_combined_xx_2 = trellis.viterbi_combined_cb(
trellis.fsm(prefix + fsm2), block, -1, -1, 1,
((1 - alpha)**0.5 * 1, (1 - alpha)**0.5 * 1j,
(1 - alpha)**0.5 * (-1j), (1 - alpha)**0.5 * (-1)),
digital.TRELLIS_EUCLIDEAN)
self.trellis_viterbi_combined_xx_1 = trellis.viterbi_combined_cb(
trellis.fsm(prefix + fsm1), block, -1, -1, 1,
(alpha**0.5 * 1, alpha**0.5 * 1j, alpha**0.5 * (-1j), alpha**0.5 *
(-1)), digital.TRELLIS_EUCLIDEAN)
self.trellis_viterbi_combined_xx_0_0 = trellis.viterbi_combined_cb(
trellis.fsm(prefix + fsm1), block, -1, -1, 1,
(alpha**0.5 * 1, alpha**0.5 * 1j, alpha**0.5 * (-1j), alpha**0.5 *
(-1)), digital.TRELLIS_EUCLIDEAN)
self.trellis_viterbi_combined_xx_0 = trellis.viterbi_combined_cb(
trellis.fsm(prefix + fsm2), block, -1, -1, 1,
((1 - alpha)**0.5 * 1, (1 - alpha)**0.5 * 1j,
(1 - alpha)**0.5 * (-1j), (1 - alpha)**0.5 * (-1)),
digital.TRELLIS_EUCLIDEAN)
self.trellis_encoder_xx_2_0 = trellis.encoder_bb(
trellis.fsm(prefix + fsm2), 0, 0) if False else trellis.encoder_bb(
trellis.fsm(prefix + fsm2), 0)
self.trellis_encoder_xx_2 = trellis.encoder_bb(
trellis.fsm(prefix + fsm1), 0, 0) if False else trellis.encoder_bb(
trellis.fsm(prefix + fsm1), 0)
self.trellis_encoder_xx_1 = trellis.encoder_bb(
trellis.fsm(prefix + fsm2), 0, 0) if False else trellis.encoder_bb(
trellis.fsm(prefix + fsm2), 0)
self.trellis_encoder_xx_0 = trellis.encoder_bb(
trellis.fsm(prefix + fsm1), 0, 0) if False else trellis.encoder_bb(
trellis.fsm(prefix + fsm1), 0)
self._snr_db_range = Range(0, 20, 1, 10, 200)
self._snr_db_win = RangeWidget(self._snr_db_range, self.set_snr_db,
'P/sigma^2 (dB)', "counter_slider",
float)
self.top_layout.addWidget(self._snr_db_win)
self.qtgui_number_sink_0_0_0_0 = qtgui.number_sink(
gr.sizeof_float, 0, qtgui.NUM_GRAPH_HORIZ, 1)
self.qtgui_number_sink_0_0_0_0.set_update_time(0.10)
self.qtgui_number_sink_0_0_0_0.set_title(
'BER 1 (after cancelling user 2)')
labels = ['BER', '', '', '', '', '', '', '', '', '']
units = ['', '', '', '', '', '', '', '', '', '']
colors = [("black", "black"), ("black", "black"), ("black", "black"),
("black", "black"), ("black", "black"), ("black", "black"),
("black", "black"), ("black", "black"), ("black", "black"),
("black", "black")]
factor = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
for i in xrange(1):
self.qtgui_number_sink_0_0_0_0.set_min(i, 0)
self.qtgui_number_sink_0_0_0_0.set_max(i, 1)
self.qtgui_number_sink_0_0_0_0.set_color(i, colors[i][0],
colors[i][1])
if len(labels[i]) == 0:
self.qtgui_number_sink_0_0_0_0.set_label(
i, "Data {0}".format(i))
else:
self.qtgui_number_sink_0_0_0_0.set_label(i, labels[i])
self.qtgui_number_sink_0_0_0_0.set_unit(i, units[i])
self.qtgui_number_sink_0_0_0_0.set_factor(i, factor[i])
self.qtgui_number_sink_0_0_0_0.enable_autoscale(False)
self._qtgui_number_sink_0_0_0_0_win = sip.wrapinstance(
self.qtgui_number_sink_0_0_0_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_number_sink_0_0_0_0_win, 1,
0, 1, 1)
self.qtgui_number_sink_0_0_0 = qtgui.number_sink(
gr.sizeof_float, 0, qtgui.NUM_GRAPH_HORIZ, 1)
self.qtgui_number_sink_0_0_0.set_update_time(0.10)
self.qtgui_number_sink_0_0_0.set_title(
"BER 2 (after cancelling user 1)")
labels = ['BER', '', '', '', '', '', '', '', '', '']
units = ['', '', '', '', '', '', '', '', '', '']
colors = [("black", "black"), ("black", "black"), ("black", "black"),
("black", "black"), ("black", "black"), ("black", "black"),
("black", "black"), ("black", "black"), ("black", "black"),
("black", "black")]
factor = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
for i in xrange(1):
self.qtgui_number_sink_0_0_0.set_min(i, 0)
self.qtgui_number_sink_0_0_0.set_max(i, 1)
self.qtgui_number_sink_0_0_0.set_color(i, colors[i][0],
colors[i][1])
if len(labels[i]) == 0:
self.qtgui_number_sink_0_0_0.set_label(i, "Data {0}".format(i))
else:
self.qtgui_number_sink_0_0_0.set_label(i, labels[i])
self.qtgui_number_sink_0_0_0.set_unit(i, units[i])
self.qtgui_number_sink_0_0_0.set_factor(i, factor[i])
self.qtgui_number_sink_0_0_0.enable_autoscale(False)
self._qtgui_number_sink_0_0_0_win = sip.wrapinstance(
self.qtgui_number_sink_0_0_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_number_sink_0_0_0_win, 1, 1,
1, 1)
self.qtgui_number_sink_0_0 = qtgui.number_sink(gr.sizeof_float, 0,
qtgui.NUM_GRAPH_HORIZ,
1)
self.qtgui_number_sink_0_0.set_update_time(0.10)
self.qtgui_number_sink_0_0.set_title("BER 2 (Raw)")
labels = ['BER', '', '', '', '', '', '', '', '', '']
units = ['', '', '', '', '', '', '', '', '', '']
colors = [("black", "black"), ("black", "black"), ("black", "black"),
("black", "black"), ("black", "black"), ("black", "black"),
("black", "black"), ("black", "black"), ("black", "black"),
("black", "black")]
factor = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
for i in xrange(1):
self.qtgui_number_sink_0_0.set_min(i, 0)
self.qtgui_number_sink_0_0.set_max(i, 1)
self.qtgui_number_sink_0_0.set_color(i, colors[i][0], colors[i][1])
if len(labels[i]) == 0:
self.qtgui_number_sink_0_0.set_label(i, "Data {0}".format(i))
else:
self.qtgui_number_sink_0_0.set_label(i, labels[i])
self.qtgui_number_sink_0_0.set_unit(i, units[i])
self.qtgui_number_sink_0_0.set_factor(i, factor[i])
self.qtgui_number_sink_0_0.enable_autoscale(False)
self._qtgui_number_sink_0_0_win = sip.wrapinstance(
self.qtgui_number_sink_0_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_number_sink_0_0_win, 0, 1,
1, 1)
self.qtgui_number_sink_0 = qtgui.number_sink(gr.sizeof_float, 0,
qtgui.NUM_GRAPH_HORIZ, 1)
self.qtgui_number_sink_0.set_update_time(0.10)
self.qtgui_number_sink_0.set_title("BER1 (Raw)")
labels = ['BER', '', '', '', '', '', '', '', '', '']
units = ['', '', '', '', '', '', '', '', '', '']
colors = [("black", "black"), ("black", "black"), ("black", "black"),
("black", "black"), ("black", "black"), ("black", "black"),
("black", "black"), ("black", "black"), ("black", "black"),
("black", "black")]
factor = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
for i in xrange(1):
self.qtgui_number_sink_0.set_min(i, 0)
self.qtgui_number_sink_0.set_max(i, 1)
self.qtgui_number_sink_0.set_color(i, colors[i][0], colors[i][1])
if len(labels[i]) == 0:
self.qtgui_number_sink_0.set_label(i, "Data {0}".format(i))
else:
self.qtgui_number_sink_0.set_label(i, labels[i])
self.qtgui_number_sink_0.set_unit(i, units[i])
self.qtgui_number_sink_0.set_factor(i, factor[i])
self.qtgui_number_sink_0.enable_autoscale(False)
self._qtgui_number_sink_0_win = sip.wrapinstance(
self.qtgui_number_sink_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_number_sink_0_win, 0, 0, 1,
1)
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.top_layout.addWidget(self._qtgui_const_sink_x_0_win)
self.digital_chunks_to_symbols_xx_0_0_1 = digital.chunks_to_symbols_bc(
(1, 1j, -1j, -1), 1)
self.digital_chunks_to_symbols_xx_0_0_0 = digital.chunks_to_symbols_bc(
(1, 1j, -1j, -1), 1)
self.digital_chunks_to_symbols_xx_0_0 = digital.chunks_to_symbols_bc(
(1, 1j, -1j, -1), 1)
self.digital_chunks_to_symbols_xx_0 = digital.chunks_to_symbols_bc(
(1, 1j, -1j, -1), 1)
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_char * 1, R, True)
self.blocks_sub_xx_2_0 = blocks.sub_cc(1)
self.blocks_sub_xx_2 = blocks.sub_cc(1)
self.blocks_multiply_const_vxx_2_0 = blocks.multiply_const_vcc(
((1 - alpha)**0.5, ))
self.blocks_multiply_const_vxx_2 = blocks.multiply_const_vcc(
(alpha**0.5, ))
self.blocks_multiply_const_vxx_1 = blocks.multiply_const_vcc(
((1 - alpha)**0.5, ))
self.blocks_multiply_const_vxx_0 = blocks.multiply_const_vcc(
(alpha**0.5, ))
self.blocks_add_xx_0 = blocks.add_vcc(1)
self.blks2_error_rate_0_0_0_0 = grc_blks2.error_rate(
type='BER',
win_size=block * 100,
bits_per_symbol=2,
)
self.blks2_error_rate_0_0_0 = grc_blks2.error_rate(
type='BER',
win_size=block * 100,
bits_per_symbol=2,
)
self.blks2_error_rate_0_0 = grc_blks2.error_rate(
type='BER',
win_size=block * 100,
bits_per_symbol=2,
)
self.blks2_error_rate_0 = grc_blks2.error_rate(
type='BER',
win_size=block * 100,
bits_per_symbol=2,
)
self.analog_random_source_x_1 = blocks.vector_source_b(
map(int, numpy.random.randint(0, 2, block)), True)
self.analog_random_source_x_0 = blocks.vector_source_b(
map(int, numpy.random.randint(0, 2, block)), True)
self.analog_noise_source_x_0 = analog.noise_source_c(
analog.GR_GAUSSIAN, noisevar, -42)
##################################################
# Connections
##################################################
self.connect((self.analog_noise_source_x_0, 0),
(self.blocks_add_xx_0, 2))
self.connect((self.analog_random_source_x_0, 0),
(self.blks2_error_rate_0, 0))
self.connect((self.analog_random_source_x_0, 0),
(self.blks2_error_rate_0_0_0_0, 0))
self.connect((self.analog_random_source_x_0, 0),
(self.trellis_encoder_xx_0, 0))
self.connect((self.analog_random_source_x_1, 0),
(self.blks2_error_rate_0_0, 0))
self.connect((self.analog_random_source_x_1, 0),
(self.blks2_error_rate_0_0_0, 0))
self.connect((self.analog_random_source_x_1, 0),
(self.trellis_encoder_xx_1, 0))
self.connect((self.blks2_error_rate_0, 0),
(self.qtgui_number_sink_0, 0))
self.connect((self.blks2_error_rate_0_0, 0),
(self.qtgui_number_sink_0_0, 0))
self.connect((self.blks2_error_rate_0_0_0, 0),
(self.qtgui_number_sink_0_0_0, 0))
self.connect((self.blks2_error_rate_0_0_0_0, 0),
(self.qtgui_number_sink_0_0_0_0, 0))
self.connect((self.blocks_add_xx_0, 0), (self.blocks_sub_xx_2, 0))
self.connect((self.blocks_add_xx_0, 0), (self.blocks_sub_xx_2_0, 0))
self.connect((self.blocks_add_xx_0, 0), (self.qtgui_const_sink_x_0, 0))
self.connect((self.blocks_add_xx_0, 0),
(self.trellis_viterbi_combined_xx_1, 0))
self.connect((self.blocks_add_xx_0, 0),
(self.trellis_viterbi_combined_xx_2, 0))
self.connect((self.blocks_multiply_const_vxx_0, 0),
(self.blocks_add_xx_0, 0))
self.connect((self.blocks_multiply_const_vxx_1, 0),
(self.blocks_add_xx_0, 1))
self.connect((self.blocks_multiply_const_vxx_2, 0),
(self.blocks_sub_xx_2, 1))
self.connect((self.blocks_multiply_const_vxx_2_0, 0),
(self.blocks_sub_xx_2_0, 1))
self.connect((self.blocks_sub_xx_2, 0),
(self.trellis_viterbi_combined_xx_0, 0))
self.connect((self.blocks_sub_xx_2_0, 0),
(self.trellis_viterbi_combined_xx_0_0, 0))
self.connect((self.blocks_throttle_0, 0),
(self.digital_chunks_to_symbols_xx_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0, 0),
(self.blocks_multiply_const_vxx_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0_0, 0),
(self.blocks_multiply_const_vxx_1, 0))
self.connect((self.digital_chunks_to_symbols_xx_0_0_0, 0),
(self.blocks_multiply_const_vxx_2, 0))
self.connect((self.digital_chunks_to_symbols_xx_0_0_1, 0),
(self.blocks_multiply_const_vxx_2_0, 0))
self.connect((self.trellis_encoder_xx_0, 0),
(self.blocks_throttle_0, 0))
self.connect((self.trellis_encoder_xx_1, 0),
(self.digital_chunks_to_symbols_xx_0_0, 0))
self.connect((self.trellis_encoder_xx_2, 0),
(self.digital_chunks_to_symbols_xx_0_0_0, 0))
self.connect((self.trellis_encoder_xx_2_0, 0),
(self.digital_chunks_to_symbols_xx_0_0_1, 0))
self.connect((self.trellis_viterbi_combined_xx_0, 0),
(self.blks2_error_rate_0_0_0, 1))
self.connect((self.trellis_viterbi_combined_xx_0_0, 0),
(self.blks2_error_rate_0_0_0_0, 1))
self.connect((self.trellis_viterbi_combined_xx_1, 0),
(self.blks2_error_rate_0, 1))
self.connect((self.trellis_viterbi_combined_xx_1, 0),
(self.trellis_encoder_xx_2, 0))
self.connect((self.trellis_viterbi_combined_xx_2, 0),
(self.blks2_error_rate_0_0, 1))
self.connect((self.trellis_viterbi_combined_xx_2, 0),
(self.trellis_encoder_xx_2_0, 0))
def __init__(self, pkt_len, nb_pkt, EbN0dB):
gr.top_block.__init__(self, "BER vs Eb/N0 (7,5) AWGN")
##################################################
# Variables
##################################################
self.prefix = prefix = os.getcwd()
self.pkt_len = pkt_len
#self.fsm = fsm = trellis.fsm(prefix + "/fsm/229_159.fsm")
#self.fsm = fsm = trellis.fsm(prefix + "/fsm/rsc_15_13.fsm")
#self.fsm = fsm = trellis.fsm(prefix + "/fsm/171_133.fsm")
self.fsm = fsm = trellis.fsm(prefix + "/fsm/5_7.fsm")
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(
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.pack_src = blocks.pack_k_bits_bb(8)
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.trellis_viterbi = trellis.viterbi_b(trellis.fsm(fsm), pkt_len, 0,
-1)
self.viterbi = lv.viterbi(trellis.fsm(fsm), pkt_len, 0, -1)
self.lazy = lv.lazy_viterbi(trellis.fsm(fsm), pkt_len, 0, -1)
self.viterbi_vb = lv.viterbi_volk_branch(trellis.fsm(fsm), pkt_len, 0,
-1)
self.viterbi_vs = lv.viterbi_volk_state(trellis.fsm(fsm), pkt_len, 0,
-1)
self.pack_trellis_viterbi = blocks.pack_k_bits_bb(8)
self.pack_viterbi = blocks.pack_k_bits_bb(8)
self.pack_lazy = blocks.pack_k_bits_bb(8)
self.pack_viterbi_vb = blocks.pack_k_bits_bb(8)
self.pack_viterbi_vs = blocks.pack_k_bits_bb(8)
self.ber_computer_trellis_viterbi = fec.ber_bf(False)
self.ber_computer_viterbi = fec.ber_bf(False)
self.ber_computer_lazy = fec.ber_bf(False)
self.ber_computer_viterbi_vb = fec.ber_bf(False)
self.ber_computer_viterbi_vs = fec.ber_bf(False)
self.vector_sink_trellis_viterbi = blocks.vector_sink_f()
self.vector_sink_viterbi = blocks.vector_sink_f()
self.vector_sink_lazy = blocks.vector_sink_f()
self.vector_sink_viterbi_vb = blocks.vector_sink_f()
self.vector_sink_viterbi_vs = 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.trellis_viterbi, 0))
self.connect((self.metrics_computer, 0), (self.viterbi, 0))
self.connect((self.metrics_computer, 0), (self.lazy, 0))
self.connect((self.metrics_computer, 0), (self.viterbi_vb, 0))
self.connect((self.metrics_computer, 0), (self.viterbi_vs, 0))
self.connect((self.bits_src, 0), (self.pack_src, 0))
self.connect((self.trellis_viterbi, 0), (self.pack_trellis_viterbi, 0))
self.connect((self.viterbi, 0), (self.pack_viterbi, 0))
self.connect((self.lazy, 0), (self.pack_lazy, 0))
self.connect((self.viterbi_vb, 0), (self.pack_viterbi_vb, 0))
self.connect((self.viterbi_vs, 0), (self.pack_viterbi_vs, 0))
self.connect((self.pack_src, 0),
(self.ber_computer_trellis_viterbi, 0))
self.connect((self.pack_src, 0), (self.ber_computer_viterbi, 0))
self.connect((self.pack_src, 0), (self.ber_computer_lazy, 0))
self.connect((self.pack_src, 0), (self.ber_computer_viterbi_vb, 0))
self.connect((self.pack_src, 0), (self.ber_computer_viterbi_vs, 0))
self.connect((self.pack_trellis_viterbi, 0),
(self.ber_computer_trellis_viterbi, 1))
self.connect((self.pack_viterbi, 0), (self.ber_computer_viterbi, 1))
self.connect((self.pack_lazy, 0), (self.ber_computer_lazy, 1))
self.connect((self.pack_viterbi_vb, 0),
(self.ber_computer_viterbi_vb, 1))
self.connect((self.pack_viterbi_vs, 0),
(self.ber_computer_viterbi_vs, 1))
self.connect((self.ber_computer_trellis_viterbi, 0),
(self.vector_sink_trellis_viterbi, 0))
self.connect((self.ber_computer_viterbi, 0),
(self.vector_sink_viterbi, 0))
self.connect((self.ber_computer_lazy, 0), (self.vector_sink_lazy, 0))
self.connect((self.ber_computer_viterbi_vb, 0),
(self.vector_sink_viterbi_vb, 0))
self.connect((self.ber_computer_viterbi_vs, 0),
(self.vector_sink_viterbi_vs, 0))
def __init__(self,
vlen_in=1,
vlen_out=1,
n_tailbits=1,
denom_mother_code_rate=1,
gen_poly=0,
bits_per_symbol=1,
map_tab=0,
pp_0=0,
interl_seq_0_2=(0, 1),
interl_seq_1_2=(0, 1),
pp_1_tail=0,
pp_1=0,
pp_0_tail=0,
M_total=0,
part_len_top=3,
part_len_bot=3):
gr.hier_block2.__init__(
self,
"MLC 16QAM",
gr.io_signature(1, 1, gr.sizeof_char * 1),
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
)
##################################################
# Parameters
##################################################
self.vlen_in = vlen_in
self.vlen_out = vlen_out
self.n_tailbits = n_tailbits
self.denom_mother_code_rate = denom_mother_code_rate
self.gen_poly = gen_poly
self.bits_per_symbol = bits_per_symbol
self.map_tab = map_tab
self.pp_0 = pp_0
self.interl_seq_0_2 = interl_seq_0_2
self.interl_seq_1_2 = interl_seq_1_2
self.pp_1_tail = pp_1_tail
self.pp_1 = pp_1
self.pp_0_tail = pp_0_tail
self.M_total = M_total
self.part_len_top = part_len_top
self.part_len_bot = part_len_bot
##################################################
# Blocks
##################################################
self.trellis_encoder_xx_top = trellis.encoder_bb(
trellis.fsm(1, denom_mother_code_rate, gen_poly), 0,
0) if False else trellis.encoder_bb(
trellis.fsm(1, denom_mother_code_rate, gen_poly), 0)
self.trellis_encoder_xx_bot = trellis.encoder_bb(
trellis.fsm(1, denom_mother_code_rate, gen_poly), 0,
0) if False else trellis.encoder_bb(
trellis.fsm(1, denom_mother_code_rate, gen_poly), 0)
self.drm_qam_map_bc_0 = drm.qam_map_bc(map_tab, bits_per_symbol,
vlen_out, 2)
self.drm_punct_bb_top = drm.punct_bb(
pp_0, pp_0_tail,
(part_len_top + n_tailbits) * denom_mother_code_rate, vlen_out * 2,
n_tailbits * denom_mother_code_rate)
self.drm_punct_bb_bot = drm.punct_bb(
pp_1, pp_1_tail,
(part_len_bot + n_tailbits) * denom_mother_code_rate, vlen_out * 2,
n_tailbits * denom_mother_code_rate)
self.drm_partitioning_16_bb_0 = drm.partitioning_bb(vlen_in, M_total)
self.drm_interleaver_bb_top = drm.interleaver_bb((interl_seq_0_2))
self.drm_interleaver_bb_bot = drm.interleaver_bb((interl_seq_1_2))
self.drm_add_tailbits_bb_top = drm.add_tailbits_bb(
part_len_top, n_tailbits)
self.drm_add_tailbits_bb_bot = drm.add_tailbits_bb(
part_len_bot, n_tailbits)
self.blocks_unpack_k_bits_bb_top = blocks.unpack_k_bits_bb(
denom_mother_code_rate)
self.blocks_unpack_k_bits_bb_bot = blocks.unpack_k_bits_bb(
denom_mother_code_rate)
##################################################
# Connections
##################################################
self.connect((self.blocks_unpack_k_bits_bb_bot, 0),
(self.drm_punct_bb_bot, 0))
self.connect((self.blocks_unpack_k_bits_bb_top, 0),
(self.drm_punct_bb_top, 0))
self.connect((self.drm_add_tailbits_bb_bot, 0),
(self.trellis_encoder_xx_bot, 0))
self.connect((self.drm_add_tailbits_bb_top, 0),
(self.trellis_encoder_xx_top, 0))
self.connect((self.drm_interleaver_bb_bot, 0),
(self.drm_qam_map_bc_0, 1))
self.connect((self.drm_interleaver_bb_top, 0),
(self.drm_qam_map_bc_0, 0))
self.connect((self.drm_partitioning_16_bb_0, 1),
(self.drm_add_tailbits_bb_bot, 0))
self.connect((self.drm_partitioning_16_bb_0, 0),
(self.drm_add_tailbits_bb_top, 0))
self.connect((self.drm_punct_bb_bot, 0),
(self.drm_interleaver_bb_bot, 0))
self.connect((self.drm_punct_bb_top, 0),
(self.drm_interleaver_bb_top, 0))
self.connect((self.drm_qam_map_bc_0, 0), (self, 0))
self.connect((self, 0), (self.drm_partitioning_16_bb_0, 0))
self.connect((self.trellis_encoder_xx_bot, 0),
(self.blocks_unpack_k_bits_bb_bot, 0))
self.connect((self.trellis_encoder_xx_top, 0),
(self.blocks_unpack_k_bits_bb_top, 0))
def setUp (self):
self.tb = gr.top_block ()
# Read in successfully decoded live data from Matlab
linf=open('/home/demel/exchange/matlab_d.txt')
lintu=range(120)
for i in lintu:
lintu[i]=float(linf.readline())
#print lintu
# source for live data
self.srcl = gr.vector_source_f(lintu,False,120)
# Read in .txt file with example MIB encoded + CRC checksum
inf=open('/home/demel/exchange/crc.txt')
self.intu=range(40)
for i in self.intu:
self.intu[i]=float(inf.readline())
#inf=open('/home/demel/exchange/matlab_d.txt')
#intu=range(120)
#for i in range(120):
# intu[i]=float(inf.readline())
# Source and conversions
self.src = gr.vector_source_f(self.intu,False,40)
self.conv = gr.float_to_char(40,1)
# Resize vector with repetition of last part
# Vector to stream for encoder
my_map1=range(46)
for i in range(40):
my_map1[i+6]=i
for i in range(6):
my_map1[i]=i+40
self.map1 = lte.vector_resize_vbvb(my_map1,40,46)
self.vtos = gr.vector_to_stream(1*gr.sizeof_char,46)
# Encoding of input data
self.fsm = trellis.fsm(1,3,[91,121,117])
self.enc = trellis.encoder_bb(self.fsm,0)
# unpack packed bits from encoder
self.unp = gr.unpack_k_bits_bb(3)
# stream to vector
self.stov = gr.stream_to_vector(1*gr.sizeof_char,138)
# Remove first part which contains tail-biting init stuff
map2 = range(120)
for i in map2:
map2[i]= i+18
self.map2 = lte.vector_resize_vbvb(map2,138,120)
# conversion from char to float to match input of decoder
self.conv2= gr.char_to_float(120,1)
###############################################
# From here on only "receiver side" processing
###############################################
# like QPSK demodulation: NRZ coding.
vec2=range(120)
for i in vec2:
vec2[i]=float(-2.0)
self.mult = gr.multiply_const_vff(vec2)
vec=range(120)
for i in vec:
vec[i]=1
self.add = gr.add_const_vff(vec)
# this is the actual unit under test
self.vit = lte.viterbi_vfvb()
# Sinks
self.snk = gr.vector_sink_b(40)
self.snk2 = gr.vector_sink_f(120)
# connecting blocks
self.tb.connect(self.src,self.conv,self.map1,self.vtos,self.enc,self.unp)
self.tb.connect(self.unp,self.stov,self.map2,self.conv2)
self.tb.connect(self.conv2,self.mult,self.add)
self.tb.connect(self.srcl,self.vit,self.snk)
self.tb.connect(self.add,self.snk2)
def main(args):
nargs = len(args)
if nargs == 1:
infile = args[0]
outfile = None
elif nargs == 2:
infile = args[0]
outfile = args[1]
else:
sys.stderr.write("Usage: dvbs-blade.py input_file [output_file]\n")
sys.exit(1)
symbol_rate = 4500000
samp_rate = symbol_rate * 2
code_rate = dvbs.C1_2
rrc_taps = 20
center_freq = 435000000
txvga1_gain = -4
txvga2_gain = 1
bandwidth = 6000000
tb = gr.top_block()
src = blocks.file_source(gr.sizeof_char, infile, True)
dvbs_randomizer = dvbs.randomizer_bb()
dvbs_reed_solomon_enc = dvbs.reed_solomon_enc_bb()
dvbs_interleaver = dvbs.interleaver_bb()
blocks_packed_to_unpacked = blocks.packed_to_unpacked_bb(1, gr.GR_MSB_FIRST)
trellis_encoder = trellis.encoder_bb(trellis.fsm(1, 2, (0171, 0133)), 0)
blocks_unpack_k_bits = blocks.unpack_k_bits_bb(2)
dvbs_puncture = dvbs.puncture_bb(code_rate)
blocks_pack_k_bits = blocks.pack_k_bits_bb(2)
dvbs_modulator = dvbs.modulator_bc()
fft_filter = filter.fft_filter_ccc(
1, (firdes.root_raised_cosine(1.79, samp_rate, samp_rate / 2, 0.35, rrc_taps)), 1
)
fft_filter.declare_sample_delay(0)
out = osmosdr.sink(args="bladerf=0,buffers=128,buflen=32768")
out.set_sample_rate(samp_rate)
out.set_center_freq(center_freq, 0)
out.set_freq_corr(0, 0)
out.set_gain(txvga2_gain, 0)
out.set_bb_gain(txvga1_gain, 0)
out.set_bandwidth(bandwidth, 0)
tb.connect(src, dvbs_randomizer)
tb.connect(dvbs_randomizer, dvbs_reed_solomon_enc)
tb.connect(dvbs_reed_solomon_enc, dvbs_interleaver)
tb.connect(dvbs_interleaver, blocks_packed_to_unpacked)
tb.connect(blocks_packed_to_unpacked, trellis_encoder)
tb.connect(trellis_encoder, blocks_unpack_k_bits)
tb.connect(blocks_unpack_k_bits, dvbs_puncture)
tb.connect(dvbs_puncture, blocks_pack_k_bits)
tb.connect(blocks_pack_k_bits, dvbs_modulator)
tb.connect(dvbs_modulator, fft_filter)
tb.connect(fft_filter, out)
if outfile:
dst = blocks.file_sink(gr.sizeof_gr_complex, outfile)
tb.connect(fft_filter, dst)
tb.run()
