def __init__(self, size, factor, itemsize=gr.sizeof_gr_complex):
"""
size: (int) vector size (FFT size) of next block
factor: (int) output will have this many more samples than input
If size is not divisible by factor, then the output will necessarily have jitter.
"""
size = int(size)
factor = int(factor)
# assert size % factor == 0
offset = size // factor
gr.hier_block2.__init__(
self, type(self).__name__,
gr.io_signature(1, 1, itemsize),
gr.io_signature(1, 1, itemsize),
)
if factor == 1:
# No duplication needed; simplify flowgraph
# GR refused to connect self to self, so insert a dummy block
self.connect(self, blocks.copy(itemsize), self)
else:
interleave = blocks.interleave(itemsize * size)
self.connect(
interleave,
blocks.vector_to_stream(itemsize, size),
self)
for i in xrange(0, factor):
self.connect(
self,
blocks.delay(itemsize, (factor - 1 - i) * offset),
blocks.stream_to_vector(itemsize, size),
(interleave, i))
def __init__(self, fft_len, cp_len, nofdm_symbols):
gr.hier_block2.__init__(
self, "ofdm_basebandsignal_to_frames_cvc",
gr.io_signature(1, 1, gr.sizeof_gr_complex),
gr.io_signature(1, 1, fft_len * gr.sizeof_gr_complex))
self.fft_len = fft_len
self.cp_len = cp_len
self.nofdm_symbols = nofdm_symbols
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
)
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))
self.connect(frames, self)
def __init__(self,subcarriers, frame_length):
#config = station_configuration()
total_subc = subcarriers
vlen = total_subc
gr.hier_block2.__init__(self,"ofdm_frame_sampler_grc",
gr.io_signature2(2,2,gr.sizeof_gr_complex*vlen,
gr.sizeof_char),
gr.io_signature2(2,2,gr.sizeof_gr_complex*vlen,
gr.sizeof_char))
ft = [0] * frame_length
ft[0] = 1
# The next block ensures that only complete frames find their way into
# the old outer receiver. The dynamic frame start trigger is hence
# replaced with a static one, fixed to the frame length.
frame_sampler = vector_sampler( gr.sizeof_gr_complex * total_subc,
frame_length )
symbol_output = blocks.vector_to_stream( gr.sizeof_gr_complex * total_subc,
frame_length )
delayed_frame_start = blocks.delay( gr.sizeof_char, frame_length - 1 )
damn_static_frame_trigger = blocks.vector_source_b( ft, True )
self.connect( self, frame_sampler, symbol_output, self )
self.connect( (self,1), delayed_frame_start, ( frame_sampler, 1 ) )
self.connect( damn_static_frame_trigger, (self,1) )
def __init__(self):
gr.top_block.__init__(self, "Top Block")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 32000
self.delay = delay = 0
##################################################
# Blocks
##################################################
self.blocks_unpack_k_bits_bb_0 = blocks.unpack_k_bits_bb(8)
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_char*1, samp_rate * 20,True)
self.blocks_pack_k_bits_bb_1 = blocks.pack_k_bits_bb(8)
self.blocks_head_1 = blocks.head(gr.sizeof_char*1, int(7 * 200e3 * 15/15))
self.blocks_file_source_0 = blocks.file_source(gr.sizeof_char*1, '/home/rcampello/Main/FixedPath/OOT Gnuradio/gr-ITpp/examples/BPSK Ladams/27-06-2019 (BER Sem Diff)/Com C\xc3\xb3digo/output2.data', True)
self.blocks_file_sink_0 = blocks.file_sink(gr.sizeof_char*1, '/home/rcampello/Main/FixedPath/OOT Gnuradio/gr-ITpp/examples/BPSK Ladams/27-06-2019 (BER Sem Diff)/Com C\xc3\xb3digo/output3.data', False)
self.blocks_file_sink_0.set_unbuffered(False)
self.blocks_delay_0 = blocks.delay(gr.sizeof_char*1, int(delay))
self.ITpp_BCH_Decoder_0 = ITpp.BCH_Decoder(15, 2)
##################################################
# Connections
##################################################
self.connect((self.ITpp_BCH_Decoder_0, 0), (self.blocks_pack_k_bits_bb_1, 0))
self.connect((self.blocks_delay_0, 0), (self.ITpp_BCH_Decoder_0, 0))
self.connect((self.blocks_file_source_0, 0), (self.blocks_throttle_0, 0))
self.connect((self.blocks_head_1, 0), (self.blocks_file_sink_0, 0))
self.connect((self.blocks_pack_k_bits_bb_1, 0), (self.blocks_head_1, 0))
self.connect((self.blocks_throttle_0, 0), (self.blocks_unpack_k_bits_bb_0, 0))
self.connect((self.blocks_unpack_k_bits_bb_0, 0), (self.blocks_delay_0, 0))
def __init__(self, code_index, code_size=31, DataRate_Kbps=850, Nburst=32, Ncpb=16, MeanPRF_KHz=15600, Nsync=64, deltaL=16, Nsdf=8, bypass_conv_enc=0, isRangingPacket=0, msgq_limit=2):
"""
Hierarchical block for the IEEE 802.15.4a UWB modulation.
Packets to be sent are enqueued by calling send_pkt.
The output is the complex modulated signal at baseband.
@param fg: flow graph
@type fg: flow graph
@param msgq_limit: maximum number of messages in message queue
@type msgq_limit: int
@param pad_for_usrp: If true, packets are padded such that they end up a multiple of 128 samples
@param bypass_conv_enc: bypass the convolutional encoding (1) or not (0)[default]
@type bypass_conv_enc: int
"""
gr.hier_block2.__init__(self, "ieee804154a_uwb_mod_pkt",
gr.io_signature(0, 0, 0), # Input signature
gr.io_signature(1, 1, gr.sizeof_char)) # Output signature
# setting parameters
self.code_index = code_index
self.code_size = code_size
self.DataRate_Kbps = DataRate_Kbps
self.Nburst = Nburst
self.Ncpb = Ncpb
self.MeanPRF_KHz = MeanPRF_KHz
self.Nsync = Nsync
self.deltaL = deltaL
self.Nsdf = Nsdf
self.bypass_conv_enc = bypass_conv_enc
self.isRangingPacket = bypass_conv_enc
############ THE SYNCH PATH
# synchronization header
self.synch_seq = ieee802_15_4a.msg_formatter.set_synch(self.code_size, self.code_index, self.Nsync, self.deltaL, self.Nsdf)
self.synch_msg = gr.message_from_string(self.synch_seq)
self._synch_input = blocks.message_source(gr.sizeof_char, msgq_limit)
############ THE MAIN PATH
# accepts messages from the outside world
self._pkt_input = blocks.message_source(gr.sizeof_char, msgq_limit)
# convolutional encoding
self._convolutional_encoder = ieee802_15_4a.conv_encoder (self.bypass_conv_enc);
# BPSK BPM modulator
self._modulator = ieee802_15_4a.bpsk_bpm_modulator(self.code_index, self.Nburst, self.Ncpb)
# connect the blocks
self.connect(self._pkt_input, self._convolutional_encoder, self._modulator)
# THE CONNECTION BETWEEN THE PATHS
self.delay = blocks.delay (gr.sizeof_char, len(self.synch_seq))
self.sum2 = blocks.or_bb(1);
self.connect (self._synch_input, (self.sum2, 0))
self.connect (self._modulator, self.delay)
self.connect (self.delay, (self.sum2, 1))
self.connect (self.sum2, self)
def create_demux_input_blocks(self):
bw = (float(self.used_carriers) / float(self.fft_len)) / 2.0
tbw = bw * 0.08
filter_taps = firdes.low_pass(1.0, 1, bw + tbw, tbw,
firdes.WIN_HAMMING, 6.76)
self.fft_filter = filter.fft_filter_ccc(1, filter_taps, 1)
self.delay = blocks.delay(gr.sizeof_gr_complex,
self.fft_len + self.cp_len)
self.ofdm_sync = digital.ofdm_sync_sc_cfb(self.fft_len, self.cp_len,
False)
self.freq_modulator = analog.frequency_modulator_fc(-2.0 /
self.fft_len)
self.multiply = blocks.multiply_vcc(1)
header_len = 3 #In OFDM Symbols
self.header_payload_demux = digital.header_payload_demux(
header_len,
self.fft_len,
self.cp_len,
self.rx_length_tag_key,
"",
True,
gr.sizeof_gr_complex,
"rx_time",
self.samp_rate,
(),
0,
)
def __init__(self, delay=0, taps=[]):
gr.hier_block2.__init__(
self,
"Conj FS IQBal",
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
)
##################################################
# Parameters
##################################################
self.delay = delay
self.taps = taps
##################################################
# Blocks
##################################################
self.filter_fir_filter_xxx_0 = filter.fir_filter_ccc(1, (taps))
self.delay_0 = blocks.delay(gr.sizeof_gr_complex * 1, delay)
self.blocks_conjugate_cc_0 = blocks.conjugate_cc()
self.blocks_add_xx_0 = blocks.add_vcc(1)
##################################################
# Connections
##################################################
self.connect((self.blocks_add_xx_0, 0), (self, 0))
self.connect((self, 0), (self.blocks_conjugate_cc_0, 0))
self.connect((self.filter_fir_filter_xxx_0, 0), (self.blocks_add_xx_0, 1))
self.connect((self.blocks_conjugate_cc_0, 0), (self.filter_fir_filter_xxx_0, 0))
self.connect((self, 0), (self.delay_0, 0))
self.connect((self.delay_0, 0), (self.blocks_add_xx_0, 0))
def __init__(self, delay):
gr.top_block.__init__(self, "Top Block")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 32000
self.delay = delay = delay
##################################################
# Blocks
##################################################
self.blocks_unpack_k_bits_bb_0 = blocks.unpack_k_bits_bb(8)
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_char*1, samp_rate * 20,True)
self.blocks_pack_k_bits_bb_1 = blocks.pack_k_bits_bb(8)
self.blocks_head_1 = blocks.head(gr.sizeof_char*1, 393460)
self.blocks_file_source_0 = blocks.file_source(gr.sizeof_char*1, '/home/rcampello/Main/FixedPath/OOT Gnuradio/gr-ITpp/examples/Simula\xc3\xa7\xc3\xa3o BPSK/Transmiss\xc3\xa3o RedDot Usando ConstellationMod/output.data', True)
self.blocks_file_sink_0 = blocks.file_sink(gr.sizeof_char*1, '/home/rcampello/Main/FixedPath/OOT Gnuradio/gr-ITpp/examples/Simula\xc3\xa7\xc3\xa3o BPSK/Transmiss\xc3\xa3o RedDot Usando ConstellationMod/output2.data', False)
self.blocks_file_sink_0.set_unbuffered(False)
self.blocks_delay_0 = blocks.delay(gr.sizeof_char*1, int(delay))
self.ITpp_BCH_Decoder_0 = ITpp.BCH_Decoder(15, 2)
##################################################
# Connections
##################################################
self.connect((self.ITpp_BCH_Decoder_0, 0), (self.blocks_pack_k_bits_bb_1, 0))
self.connect((self.blocks_delay_0, 0), (self.ITpp_BCH_Decoder_0, 0))
self.connect((self.blocks_file_source_0, 0), (self.blocks_throttle_0, 0))
self.connect((self.blocks_head_1, 0), (self.blocks_file_sink_0, 0))
self.connect((self.blocks_pack_k_bits_bb_1, 0), (self.blocks_head_1, 0))
self.connect((self.blocks_throttle_0, 0), (self.blocks_unpack_k_bits_bb_0, 0))
self.connect((self.blocks_unpack_k_bits_bb_0, 0), (self.blocks_delay_0, 0))
def __init__(self):
gr.top_block.__init__(self, "Acpqheadphones")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 48000
##################################################
# Blocks
##################################################
self.blocks_wavfile_sink_0 = blocks.wavfile_sink('pq.wav', 2, samp_rate, 16)
self.blocks_multiply_const_vxx_0 = blocks.multiply_const_vff((-1, ))
self.blocks_delay_0 = blocks.delay(gr.sizeof_float*1, 800)
self.blocks_add_xx_0 = blocks.add_vff(1)
self.audio_source_0 = audio.source(samp_rate, 'hw:1', True)
self.audio_sink_0 = audio.sink(samp_rate, 'hw:1', True)
self.analog_rail_ff_0 = analog.rail_ff(-0.8, 0.8)
self.analog_agc2_xx_0 = analog.agc2_ff(1e-1, 0.5, 0.01, 1.0)
self.analog_agc2_xx_0.set_max_gain(65536)
##################################################
# Connections
##################################################
self.connect((self.analog_agc2_xx_0, 0), (self.analog_rail_ff_0, 0))
self.connect((self.analog_rail_ff_0, 0), (self.audio_sink_0, 0))
self.connect((self.audio_source_0, 0), (self.blocks_add_xx_0, 0))
self.connect((self.audio_source_0, 0), (self.blocks_delay_0, 0))
self.connect((self.audio_source_0, 0), (self.blocks_wavfile_sink_0, 0))
self.connect((self.blocks_add_xx_0, 0), (self.analog_agc2_xx_0, 0))
self.connect((self.blocks_add_xx_0, 0), (self.blocks_wavfile_sink_0, 1))
self.connect((self.blocks_delay_0, 0), (self.blocks_multiply_const_vxx_0, 0))
self.connect((self.blocks_multiply_const_vxx_0, 0), (self.blocks_add_xx_0, 1))
def __init__(self, freq_sample_delay_samps, freq_samps_to_avg,
mag_samps_to_avg, thresh):
gr.hier_block2.__init__(
self, "Sample and Hold Detector",
gr.io_signaturev(2, 2, [gr.sizeof_float * 1, gr.sizeof_float * 1]),
gr.io_signaturev(4, 4, [
gr.sizeof_float * 1, gr.sizeof_float * 1, gr.sizeof_float * 1,
gr.sizeof_float * 1
]))
'''
Constructor
@param freq_sample_delay_samps -
@param freq_samps_to_avg -
@param mag_samps_to_avg -
@param thresh -
'''
##################################################
# Parameters
##################################################
self.freq_sample_delay_samps = freq_sample_delay_samps
self.freq_samps_to_avg = freq_samps_to_avg
self.mag_samps_to_avg = mag_samps_to_avg
self.thresh = thresh
##################################################
# Blocks
##################################################
self.edge_detector = timing_utils.edge_detector_bb(
timing_utils.RISING_EDGE)
self.threshold = blocks.threshold_ff(thresh / 4.0, thresh, 0)
self.samp_hold = blocks.sample_and_hold_ff()
self.mag_avg = blocks.moving_average_ff(int(mag_samps_to_avg),
1.0 / (mag_samps_to_avg), 4000)
self.freq_avg = blocks.moving_average_ff(int(freq_samps_to_avg),
1.0 / (freq_samps_to_avg),
4000)
self.f2c = blocks.float_to_char(1, 1)
self.delay = blocks.delay(
gr.sizeof_float * 1,
int(freq_samps_to_avg - mag_samps_to_avg +
freq_sample_delay_samps))
##################################################
# Connections
##################################################
self.connect((self.delay, 0), (self.mag_avg, 0))
self.connect((self.f2c, 0), (self.edge_detector, 0))
self.connect((self.freq_avg, 0), (self.samp_hold, 0))
self.connect((self.freq_avg, 0), (self, 1))
self.connect((self.mag_avg, 0), (self.threshold, 0))
self.connect((self.mag_avg, 0), (self, 3))
self.connect((self.samp_hold, 0), (self, 0))
self.connect((self.threshold, 0), (self.f2c, 0))
self.connect((self.threshold, 0), (self, 2))
self.connect((self, 0), (self.delay, 0))
self.connect((self, 1), (self.freq_avg, 0))
self.connect((self.edge_detector, 0), (self.samp_hold, 1))
def __init__(self):
gr.top_block.__init__(self, "Not titled yet")
##################################################
# Variables
##################################################
self.trigger_cmd = trigger_cmd = gr.tag_utils.python_to_tag(
(0, pmt.intern("TRIGGER"), pmt.PMT_T, pmt.intern("src")))
self.samp_rate = samp_rate = 10000
self.cmd_msg = cmd_msg = pmt.to_pmt({'lo_offset': 100})
##################################################
# Blocks
##################################################
self.pdu_utils_tag_message_trigger_0 = pdu_utils.tag_message_trigger_b(
pmt.intern("TRIGGER"), pmt.PMT_NIL, cmd_msg, 0, samp_rate, 0.0, 5,
False)
self.blocks_vector_source_x_0_0 = blocks.vector_source_b(
(0, 0), False, 1, [trigger_cmd])
self.blocks_throttle_0_0 = blocks.throttle(gr.sizeof_char * 1,
samp_rate, True)
self.blocks_message_debug_0 = blocks.message_debug()
self.blocks_delay_0 = blocks.delay(gr.sizeof_char * 1, 5 * samp_rate)
##################################################
# Connections
##################################################
self.msg_connect((self.pdu_utils_tag_message_trigger_0, 'msg'),
(self.blocks_message_debug_0, 'print'))
self.connect((self.blocks_delay_0, 0), (self.blocks_throttle_0_0, 0))
self.connect((self.blocks_throttle_0_0, 0),
(self.pdu_utils_tag_message_trigger_0, 0))
self.connect((self.blocks_vector_source_x_0_0, 0),
(self.blocks_delay_0, 0))
def __init__(self, fft_length, pn_weights):
gr.hier_block2.__init__(self, "modified_timing_metric",
gr.io_signature(1,1,gr.sizeof_gr_complex),
gr.io_signature(1,1,gr.sizeof_float))
assert(len(pn_weights) == fft_length)
self.input = gr.kludge_copy(gr.sizeof_gr_complex)
self.connect(self,self.input)
# P(d) = sum(0 to L-1, conj(delayed(r)) * r)
conj = gr.conjugate_cc()
mixer = gr.multiply_cc()
nominator = gr.fir_filter_ccf(1,[pn_weights[fft_length-i-1]*pn_weights[fft_length/2-i-1] for i in range(fft_length/2)])
self.connect(self.input, delay(gr.sizeof_gr_complex,fft_length/2), conj, (mixer,0))
self.connect(self.input, (mixer,1))
self.connect(mixer, nominator)
# moving_avg = P(d)
# R(d)
denominator = schmidl_denominator(fft_length)
# |P(d)| ** 2 / (R(d)) ** 2
p_mag_sqrd = gr.complex_to_mag_squared()
r_sqrd = gr.multiply_ff()
self.timing_metric = gr.divide_ff()
self.connect(nominator, p_mag_sqrd, (self.timing_metric,0))
self.connect(self.input, denominator, (r_sqrd,0))
self.connect(denominator, (r_sqrd,1))
self.connect(r_sqrd, (self.timing_metric,1))
self.connect(self.timing_metric, self)
def __init__(self, sample_rate): super(fsk_center_tracking, self).__init__( "fsk_center_tracking", gr.io_signature(1, 1, gr.sizeof_gr_complex*1), gr.io_signature(1, 1, gr.sizeof_gr_complex*1), ) # Arbitrary averaging values that seem to work OK. window_symbols = 20 symbol_rate = 19200 average_window = int(floor(sample_rate * window_symbols / symbol_rate)) self.delay = blocks.delay(gr.sizeof_gr_complex*1, average_window) self.demod = analog.quadrature_demod_cf(1) self.moving_average = blocks.moving_average_ff(average_window, 1.0 / average_window, 4000) self.vco = blocks.vco_c(sample_rate, -sample_rate, 1) self.multiply = blocks.multiply_vcc(1) self.connect((self, 0), (self.delay, 0)) self.connect((self.delay, 0), (self.multiply, 0)) self.connect((self, 0), (self.demod, 0)) self.connect((self.demod, 0), (self.moving_average, 0)) self.connect((self.moving_average, 0), (self.vco, 0)) self.connect((self.vco, 0), (self.multiply, 1)) self.connect((self.multiply, 0), (self, 0))
def __init__(self):
gr.hier_block2.__init__(
self,
"Amplitude Difference ",
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
gr.io_signature(1, 1, gr.sizeof_float * 1),
)
##################################################
# Blocks
##################################################
self.dpd_signum_func_0 = dpd.signum_func()
self.blocks_sub_xx_0 = blocks.sub_ff(1)
self.blocks_delay_0 = blocks.delay(gr.sizeof_gr_complex * 1, 1)
self.blocks_complex_to_mag_0_0 = blocks.complex_to_mag(1)
self.blocks_complex_to_mag_0 = blocks.complex_to_mag(1)
##################################################
# Connections
##################################################
self.connect((self.blocks_complex_to_mag_0, 0),
(self.blocks_sub_xx_0, 0))
self.connect((self.blocks_complex_to_mag_0_0, 0),
(self.blocks_sub_xx_0, 1))
self.connect((self.blocks_delay_0, 0),
(self.blocks_complex_to_mag_0_0, 0))
self.connect((self.blocks_sub_xx_0, 0), (self.dpd_signum_func_0, 0))
self.connect((self.dpd_signum_func_0, 0), (self, 0))
self.connect((self, 0), (self.blocks_complex_to_mag_0, 0))
self.connect((self, 0), (self.blocks_delay_0, 0))
def __init__(self,
frame_size=223,
differential=False,
dual_basis=False,
use_scrambler=True,
nasa_dsn=False,
syncword_threshold=None,
options=None):
gr.hier_block2.__init__(self, "ccsds_concatenated_deframer",
gr.io_signature(1, 1, gr.sizeof_float),
gr.io_signature(0, 0, 0))
options_block.__init__(self, options)
self.message_port_register_hier_out('out')
self.delay1 = blocks.delay(gr.sizeof_float, 1)
self.viterbi0 = ccsds_viterbi(nasa_dsn=nasa_dsn)
self.viterbi1 = ccsds_viterbi(nasa_dsn=nasa_dsn)
self.char2float0 = blocks.char_to_float(1, 1)
self.char2float1 = blocks.char_to_float(1, 1)
self.addconst0 = blocks.add_const_ff(-0.5)
self.addconst1 = blocks.add_const_ff(-0.5)
self.rs0 = ccsds_rs_deframer(frame_size, differential, dual_basis,
use_scrambler, syncword_threshold,
options)
self.rs1 = ccsds_rs_deframer(frame_size, differential, dual_basis,
use_scrambler, syncword_threshold,
options)
self.connect(self, self.viterbi0, self.char2float0, self.addconst0,
self.rs0)
self.connect(self, self.delay1, self.viterbi1, self.char2float1,
self.addconst1, self.rs1)
self.msg_connect((self.rs0, 'out'), (self, 'out'))
self.msg_connect((self.rs1, 'out'), (self, 'out'))
def __init__(self, sample_rate):
super(fsk_center_tracking, self).__init__(
"fsk_center_tracking",
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
)
# Arbitrary averaging values that seem to work OK.
window_symbols = 20
symbol_rate = 19200
average_window = int(floor(sample_rate * window_symbols / symbol_rate))
self.delay = blocks.delay(gr.sizeof_gr_complex * 1, average_window)
self.demod = analog.quadrature_demod_cf(1)
self.moving_average = blocks.moving_average_ff(average_window,
1.0 / average_window,
4000)
self.vco = blocks.vco_c(sample_rate, -sample_rate, 1)
self.multiply = blocks.multiply_vcc(1)
self.connect((self, 0), (self.delay, 0))
self.connect((self.delay, 0), (self.multiply, 0))
self.connect((self, 0), (self.demod, 0))
self.connect((self.demod, 0), (self.moving_average, 0))
self.connect((self.moving_average, 0), (self.vco, 0))
self.connect((self.vco, 0), (self.multiply, 1))
self.connect((self.multiply, 0), (self, 0))
def __init__(self, M):
gr.hier_block2.__init__(self,
"overlapping_parallel_to_serial_vcc",
gr.io_signature(2, 2, gr.sizeof_gr_complex*M), # Input signature
gr.io_signature(1, 1, gr.sizeof_gr_complex*1)) # Output signature
# Variables
self.M = M
# Assertions
assert (M>0 and int(math.log(M,2))==math.log(M,2)), "M should be of 2^n!"
# Blocks
self.vector_to_stream0 = blocks.vector_to_stream(gr.sizeof_gr_complex*M/2, 2)
self.vector_to_stream1 = blocks.vector_to_stream(gr.sizeof_gr_complex*M/2, 2)
self.delay = blocks.delay(gr.sizeof_gr_complex*M/2, 1)
self.adder = blocks.add_vcc(M/2)
# Connections
self.connect((self,0),self.vector_to_stream0)
self.connect((self,1),self.vector_to_stream1)
self.connect(self.vector_to_stream0, (self.adder,0))
self.connect(self.vector_to_stream1, self.delay)
self.connect(self.delay, (self.adder,1))
self.connect((self.adder,0), blocks.vector_to_stream(gr.sizeof_gr_complex, M/2),(self,0))
class receiver(gr.hier_block2):
def __init__(self):
gr.hier_block2.__init__(self, "ofdm_tx",
gr.io_signature(1, 1, gr.sizeof_gr_complex),
gr.io_signature(1, 1, gr.sizeof_char))
symbol_settings = para.ofdm_symbol()
self.constellationP = helpers.get_constellation(settings.PAYLOAD_BPS)
self.constellationH = helpers.get_constellation(settings.HEADER_BPS)
detector = digital.ofdm_sync_sc_cfb(symbol_settings.get_fft_length(), symbol_settings.get_cp_length())
delayer = blocks.delay(gr.sizeof_gr_complex, symbol_settings.get_time_length_of_symbol)
oscillator = analog.frequency_modulator_fc(-2.0 / symbol_settings.get_fft_length())
splitter = digital.header_payload_demux(
3,
symbol_settings.get_fft_length(), symbol_settings.get_cp_length(),
settings.LENGTH_TAG_KEY,
"",
True
)
mixer = blocks.multiply_cc()
header_fft = fft.fft_vcc(symbol_settings.get_fft_length(),, True, (), True)
chanest = digital.ofdm_chanest_vcvc(symbol_settings._generate_sync_word_one(),symbol_settings._generate_sync_word_two(), 1)
header_equalizer = digital.ofdm_equalizer_simpledfe(
symbol_settings.get_fft_length(),,
self.constellationH.base(),
symbol_settings.get_carrier_tones(),
symbol_settings.get_pilot_tones(),
symbol_settings._generate_pilot_symbols(),
symbols_skipped=0,
)
self.connect(self,detector)
self.connect(self,delay, (mixer,0), (splitter,0))
self.connect((detector,0), oscillator, (mixer,1))
self.connect((detector,1),(splitter,1))
def __init__(self, delay=0, taps=[]):
gr.hier_block2.__init__(
self,
"Conj FS IQBal",
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
)
##################################################
# Parameters
##################################################
self.delay = delay
self.taps = taps
##################################################
# Blocks
##################################################
self.filter_fir_filter_xxx_0 = filter.fir_filter_ccc(1, (taps))
self.delay_0 = blocks.delay(gr.sizeof_gr_complex * 1, delay)
self.blocks_conjugate_cc_0 = blocks.conjugate_cc()
self.blocks_add_xx_0 = blocks.add_vcc(1)
##################################################
# Connections
##################################################
self.connect((self.blocks_add_xx_0, 0), (self, 0))
self.connect((self, 0), (self.blocks_conjugate_cc_0, 0))
self.connect((self.filter_fir_filter_xxx_0, 0),
(self.blocks_add_xx_0, 1))
self.connect((self.blocks_conjugate_cc_0, 0),
(self.filter_fir_filter_xxx_0, 0))
self.connect((self, 0), (self.delay_0, 0))
self.connect((self.delay_0, 0), (self.blocks_add_xx_0, 0))
def __init__(self, fft_len, freq_sample_delay_samps, freq_samps_to_avg, mag_samps_to_avg, thresh):
gr.hier_block2.__init__(self,
"Coarse Dehopper",
gr.io_signature(1, 1, gr.sizeof_gr_complex*1),
gr.io_signature(1, 1, gr.sizeof_gr_complex*1))
'''
Constructor
@param fft_len -
@param freq_sample_delay_samps -
@param freq_samps_to_avg -
@param mag_samps_to_avg -
@param thresh -
'''
##################################################
# Parameters
##################################################
self.fft_len = fft_len
self.freq_sample_delay_samps = freq_sample_delay_samps
self.freq_samps_to_avg = freq_samps_to_avg
self.mag_samps_to_avg = mag_samps_to_avg
self.thresh = thresh
##################################################
# Blocks
##################################################
self.s_and_h_detector = s_and_h_detector(
freq_sample_delay_samps=freq_sample_delay_samps,
freq_samps_to_avg=freq_samps_to_avg,
mag_samps_to_avg=mag_samps_to_avg,
thresh=thresh,
)
self.resamp = pfb.arb_resampler_ccf(1.0 / (fft_len / 4.0), taps=None, flt_size=32)
self.resamp.declare_sample_delay(0)
self.fir = filter.fir_filter_ccc(2, (firdes.low_pass_2(1,1,.30,.05,60)))
self.fir.declare_sample_delay(0)
self.fft_peak = fft_peak(fft_len=fft_len)
self.vco = blocks.vco_c(1, 2.0 * pi / fft_len, 1)
self.mult_conj = blocks.multiply_conjugate_cc(1)
self.delay = blocks.delay(gr.sizeof_gr_complex*1, int(freq_samps_to_avg) + freq_sample_delay_samps)
self.c2mag = blocks.complex_to_mag(1)
##################################################
# Connections
##################################################
self.connect((self.c2mag, 0), (self.s_and_h_detector, 0))
self.connect((self.delay, 0), (self.mult_conj, 0))
self.connect((self.mult_conj, 0), (self.fir, 0))
self.connect((self.vco, 0), (self.mult_conj, 1))
self.connect((self.fft_peak, 0), (self.s_and_h_detector, 1))
self.connect((self.fir, 0), (self.resamp, 0))
self.connect((self, 0), (self.c2mag, 0))
self.connect((self, 0), (self.delay, 0))
self.connect((self, 0), (self.fft_peak, 0))
self.connect((self.resamp, 0), (self, 0))
self.connect((self.s_and_h_detector, 0), (self.vco, 0))
def oqpsktx(self,
carrier=10000,
samp_rate=80000,
bw=4000,
amp=1,
code=codes.mycode,
**kwargs):
code_table, code_len = codes.codes2table(code), len(code)
chunk_len = int(log(code_len, 2))
topblock(self, carrier, samp_rate, bw, amp)
##################################################
# Blocks
##################################################
self.blocks_packed_to_unpacked_xx_0 = blocks.packed_to_unpacked_bb(
chunk_len, gr.GR_LSB_FIRST)
self.digital_chunks_to_symbols_xx_0 = digital.chunks_to_symbols_bc(
(code_table), code_len)
self.blocks_repeat_0 = blocks.repeat(gr.sizeof_gr_complex * 1, 4)
self.blocks_vector_source_x_0 = blocks.vector_source_c(
[0, sin(pi / 4), 1, sin(3 * pi / 4)], True, 1, [])
self.blocks_multiply_xx_0 = blocks.multiply_vcc(1)
self.blocks_complex_to_float_0 = blocks.complex_to_float(1)
self.blocks_delay_0 = blocks.delay(gr.sizeof_float * 1, 2)
self.blocks_float_to_complex_0 = blocks.float_to_complex(1)
self.rational_resampler_xxx_0 = filter.rational_resampler_ccc(
interpolation=int(samp_rate / bw),
decimation=1,
taps=None,
fractional_bw=None,
)
##################################################
# Connections
##################################################
self.connect((self.source, 0), (self.blocks_packed_to_unpacked_xx_0, 0))
self.connect((self.blocks_packed_to_unpacked_xx_0, 0),
(self.digital_chunks_to_symbols_xx_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0, 0),
(self.blocks_repeat_0, 0))
self.connect((self.blocks_vector_source_x_0, 0),
(self.blocks_multiply_xx_0, 0))
self.connect((self.blocks_repeat_0, 0), (self.blocks_multiply_xx_0, 1))
self.connect((self.blocks_multiply_xx_0, 0),
(self.blocks_complex_to_float_0, 0))
self.connect((self.blocks_complex_to_float_0, 0),
(self.blocks_float_to_complex_0, 0))
self.connect((self.blocks_complex_to_float_0, 1), (self.blocks_delay_0, 0))
self.connect((self.blocks_delay_0, 0), (self.blocks_float_to_complex_0, 1))
self.connect((self.blocks_float_to_complex_0, 0),
(self.rational_resampler_xxx_0, 0))
self.connect((self.rational_resampler_xxx_0, 0), (self.sink, 0))
def test_delay():
top = gr.top_block()
src = blocks.null_source(gr.sizeof_gr_complex)
delay = blocks.delay(gr.sizeof_gr_complex, 3000)
probe = blocks.probe_rate(gr.sizeof_gr_complex)
top.connect(src, delay, probe)
return top, probe
def __init__(self, preamble_length):
gr.hier_block2.__init__(self, "fqsweep_corr",
gr.io_signature(1, 1, gr.sizeof_gr_complex),
gr.io_signature(2, 2, gr.sizeof_gr_complex))
despread = self.gen_despread(preamble_length, np.pi * 3.0 / 4.0)
# Blocks
self.equiv_delay = blocks.delay(gr.sizeof_gr_complex,
preamble_length + 2)
self.despread_src = blocks.vector_source_c(despread, True, 1, [])
self.despread_mul = blocks.multiply_conjugate_cc(1)
self.deriv_0_delay = blocks.delay(gr.sizeof_gr_complex, 1)
self.deriv_0_mul = blocks.multiply_conjugate_cc(1)
self.deriv_1_delay = blocks.delay(gr.sizeof_gr_complex, 1)
self.deriv_1_mul = blocks.multiply_conjugate_cc(1)
self.avg = blocks.moving_average_cc(preamble_length,
1.0 / preamble_length, 2**16)
# Connections
# pass output
self.connect((self, 0), (self.equiv_delay, 0))
self.connect((self.equiv_delay, 0), (self, 0))
# despread
self.connect((self, 0), (self.despread_mul, 0))
self.connect((self.despread_src, 0), (self.despread_mul, 1))
# first phase derivate (frequency)
self.connect((self.despread_mul, 0), (self.deriv_0_mul, 0))
self.connect((self.despread_mul, 0), (self.deriv_0_delay, 0))
self.connect((self.deriv_0_delay, 0), (self.deriv_0_mul, 1))
# second phase derivate (change of frequency)
self.connect((self.deriv_0_mul, 0), (self.deriv_1_mul, 0))
self.connect((self.deriv_0_mul, 0), (self.deriv_1_delay, 0))
self.connect((self.deriv_1_delay, 0), (self.deriv_1_mul, 1))
# average
self.connect((self.deriv_1_mul, 0), (self.avg, 0))
self.connect((self.avg, 0), (self, 1))
def test_delay():
top = gr.top_block()
src = blocks.null_source(gr.sizeof_gr_complex)
delay = blocks.delay(gr.sizeof_gr_complex, 3000)
probe = blocks.probe_rate(gr.sizeof_gr_complex)
top.connect(src, delay, probe)
return top, probe
def __init__(self):
gr.hier_block2.__init__(
self, "IEEE802.15.4 OQPSK PHY",
gr.io_signature(1, 1, gr.sizeof_gr_complex*1),
gr.io_signature(1, 1, gr.sizeof_gr_complex*1),
)
self.message_port_register_hier_in("txin")
self.message_port_register_hier_out("rxout")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 4000000
##################################################
# Blocks
##################################################
self.single_pole_iir_filter_xx_0 = filter.single_pole_iir_filter_ff(0.00016, 1)
self.ieee802_15_4_packet_sink_0 = ieee802_15_4.packet_sink(10)
self.ieee802_15_4_access_code_prefixer_0 = ieee802_15_4.access_code_prefixer()
self.foo_burst_tagger_0 = foo.burst_tagger(pmt.intern("pdu_length"), 128)
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(2, 0.000225, 0.5, 0.03, 0.0002)
self.digital_chunks_to_symbols_xx_0 = digital.chunks_to_symbols_bc(([(1+1j), (-1+1j), (1-1j), (-1+1j), (1+1j), (-1-1j), (-1-1j), (1+1j), (-1+1j), (-1+1j), (-1-1j), (1-1j), (-1-1j), (1-1j), (1+1j), (1-1j), (1-1j), (-1-1j), (1+1j), (-1-1j), (1-1j), (-1+1j), (-1+1j), (1-1j), (-1-1j), (-1-1j), (-1+1j), (1+1j), (-1+1j), (1+1j), (1-1j), (1+1j), (-1+1j), (-1+1j), (-1-1j), (1-1j), (-1-1j), (1-1j), (1+1j), (1-1j), (1+1j), (-1+1j), (1-1j), (-1+1j), (1+1j), (-1-1j), (-1-1j), (1+1j), (-1-1j), (-1-1j), (-1+1j), (1+1j), (-1+1j), (1+1j), (1-1j), (1+1j), (1-1j), (-1-1j), (1+1j), (-1-1j), (1-1j), (-1+1j), (-1+1j), (1-1j), (-1-1j), (1-1j), (1+1j), (1-1j), (1+1j), (-1+1j), (1-1j), (-1+1j), (1+1j), (-1-1j), (-1-1j), (1+1j), (-1+1j), (-1+1j), (-1-1j), (1-1j), (-1+1j), (1+1j), (1-1j), (1+1j), (1-1j), (-1-1j), (1+1j), (-1-1j), (1-1j), (-1+1j), (-1+1j), (1-1j), (-1-1j), (-1-1j), (-1+1j), (1+1j), (1+1j), (-1-1j), (-1-1j), (1+1j), (-1+1j), (-1+1j), (-1-1j), (1-1j), (-1-1j), (1-1j), (1+1j), (1-1j), (1+1j), (-1+1j), (1-1j), (-1+1j), (1-1j), (-1+1j), (-1+1j), (1-1j), (-1-1j), (-1-1j), (-1+1j), (1+1j), (-1+1j), (1+1j), (1-1j), (1+1j), (1-1j), (-1-1j), (1+1j), (-1-1j), (1+1j), (1-1j), (1+1j), (-1+1j), (1-1j), (-1+1j), (1+1j), (-1-1j), (-1-1j), (1+1j), (-1+1j), (-1+1j), (-1-1j), (1-1j), (-1-1j), (1-1j), (1-1j), (1+1j), (1-1j), (-1-1j), (1+1j), (-1-1j), (1-1j), (-1+1j), (-1+1j), (1-1j), (-1-1j), (-1-1j), (-1+1j), (1+1j), (-1+1j), (1+1j), (-1-1j), (1+1j), (-1+1j), (-1+1j), (-1-1j), (1-1j), (-1-1j), (1-1j), (1+1j), (1-1j), (1+1j), (-1+1j), (1-1j), (-1+1j), (1+1j), (-1-1j), (-1+1j), (1-1j), (-1-1j), (-1-1j), (-1+1j), (1+1j), (-1+1j), (1+1j), (1-1j), (1+1j), (1-1j), (-1-1j), (1+1j), (-1-1j), (1-1j), (-1+1j), (-1-1j), (1-1j), (-1-1j), (1-1j), (1+1j), (1-1j), (1+1j), (-1+1j), (1-1j), (-1+1j), (1+1j), (-1-1j), (-1-1j), (1+1j), (-1+1j), (-1+1j), (-1+1j), (1+1j), (-1+1j), (1+1j), (1-1j), (1+1j), (1-1j), (-1-1j), (1+1j), (-1-1j), (1-1j), (-1+1j), (-1+1j), (1-1j), (-1-1j), (-1-1j), (1-1j), (-1+1j), (1+1j), (-1-1j), (-1-1j), (1+1j), (-1+1j), (-1+1j), (-1-1j), (1-1j), (-1-1j), (1-1j), (1+1j), (1-1j), (1+1j), (-1+1j), (1+1j), (-1-1j), (1-1j), (-1+1j), (-1+1j), (1-1j), (-1-1j), (-1-1j), (-1+1j), (1+1j), (-1+1j), (1+1j), (1-1j), (1+1j), (1-1j), (-1-1j)]), 16)
self.blocks_vector_source_x_0 = blocks.vector_source_c([0, sin(pi/4), 1, sin(3*pi/4)], True, 1, [])
self.blocks_sub_xx_0 = blocks.sub_ff(1)
self.blocks_repeat_0 = blocks.repeat(gr.sizeof_gr_complex*1, 4)
self.blocks_pdu_to_tagged_stream_0_0_0 = blocks.pdu_to_tagged_stream(blocks.byte_t, 'pdu_length')
self.blocks_packed_to_unpacked_xx_0 = blocks.packed_to_unpacked_bb(4, gr.GR_LSB_FIRST)
self.blocks_multiply_xx_0 = blocks.multiply_vcc(1)
self.blocks_float_to_complex_0 = blocks.float_to_complex(1)
self.blocks_delay_0 = blocks.delay(gr.sizeof_float*1, 2)
self.blocks_complex_to_float_0 = blocks.complex_to_float(1)
self.analog_quadrature_demod_cf_0 = analog.quadrature_demod_cf(1)
##################################################
# Connections
##################################################
self.msg_connect((self.ieee802_15_4_access_code_prefixer_0, 'out'), (self.blocks_pdu_to_tagged_stream_0_0_0, 'pdus'))
self.msg_connect((self.ieee802_15_4_packet_sink_0, 'out'), (self, 'rxout'))
self.msg_connect((self, 'txin'), (self.ieee802_15_4_access_code_prefixer_0, 'in'))
self.connect((self.analog_quadrature_demod_cf_0, 0), (self.blocks_sub_xx_0, 0))
self.connect((self.analog_quadrature_demod_cf_0, 0), (self.single_pole_iir_filter_xx_0, 0))
self.connect((self.blocks_complex_to_float_0, 1), (self.blocks_delay_0, 0))
self.connect((self.blocks_complex_to_float_0, 0), (self.blocks_float_to_complex_0, 0))
self.connect((self.blocks_delay_0, 0), (self.blocks_float_to_complex_0, 1))
self.connect((self.blocks_float_to_complex_0, 0), (self.foo_burst_tagger_0, 0))
self.connect((self.blocks_multiply_xx_0, 0), (self.blocks_complex_to_float_0, 0))
self.connect((self.blocks_packed_to_unpacked_xx_0, 0), (self.digital_chunks_to_symbols_xx_0, 0))
self.connect((self.blocks_pdu_to_tagged_stream_0_0_0, 0), (self.blocks_packed_to_unpacked_xx_0, 0))
self.connect((self.blocks_repeat_0, 0), (self.blocks_multiply_xx_0, 1))
self.connect((self.blocks_sub_xx_0, 0), (self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.blocks_vector_source_x_0, 0), (self.blocks_multiply_xx_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0, 0), (self.blocks_repeat_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.ieee802_15_4_packet_sink_0, 0))
self.connect((self.foo_burst_tagger_0, 0), (self, 0))
self.connect((self, 0), (self.analog_quadrature_demod_cf_0, 0))
self.connect((self.single_pole_iir_filter_xx_0, 0), (self.blocks_sub_xx_0, 1))
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, txstr="Hello World", carrier=10000, samp_rate = 80000, bw=4000, amp=1):
gr.top_block.__init__(self, "Top Block")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate
self.carrier = carrier
self.bw = bw
##################################################
# Blocks
##################################################
self.source = blocks.vector_source_b(tuple(bytearray(txstr)), False, 1, [])
self.rational_resampler_xxx_0 = filter.rational_resampler_ccc(
interpolation=int(samp_rate/bw),
decimation=1,
taps=None,
fractional_bw=None,
)
self.digital_chunks_to_symbols_xx_0 = digital.chunks_to_symbols_bc((CODE_TABLE), CODE_LEN)
self.blocks_vector_source_x_0 = blocks.vector_source_c([0, sin(pi/4), 1, sin(3*pi/4)], True, 1, [])
self.blocks_repeat_0 = blocks.repeat(gr.sizeof_gr_complex*1, 4)
self.blocks_packed_to_unpacked_xx_0 = blocks.packed_to_unpacked_bb(CHUNK_LEN, gr.GR_LSB_FIRST)
self.blocks_multiply_xx_0_0 = blocks.multiply_vcc(1)
self.blocks_multiply_xx_0 = blocks.multiply_vcc(1)
self.blocks_float_to_complex_0 = blocks.float_to_complex(1)
self.blocks_delay_0 = blocks.delay(gr.sizeof_float*1, 2)
self.blocks_complex_to_real_0 = blocks.complex_to_real(1)
self.blocks_complex_to_float_0 = blocks.complex_to_float(1)
self.audio_sink_0 = audio.sink(samp_rate, "")
#XXX Hack: 0.07 should actually be parameter amp, but RPI crashes
self.analog_sig_source_x_0 = analog.sig_source_c(samp_rate, analog.GR_COS_WAVE, carrier, 0.07, 0)
##################################################
# Connections
##################################################
self.connect((self.source, 0), (self.blocks_packed_to_unpacked_xx_0, 0))
self.connect((self.blocks_packed_to_unpacked_xx_0, 0), (self.digital_chunks_to_symbols_xx_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0, 0), (self.blocks_repeat_0, 0))
self.connect((self.blocks_vector_source_x_0, 0), (self.blocks_multiply_xx_0, 0))
self.connect((self.blocks_repeat_0, 0), (self.blocks_multiply_xx_0, 1))
self.connect((self.blocks_multiply_xx_0, 0), (self.blocks_complex_to_float_0, 0))
self.connect((self.blocks_complex_to_float_0, 0), (self.blocks_float_to_complex_0, 0))
self.connect((self.blocks_complex_to_float_0, 1), (self.blocks_delay_0, 0))
self.connect((self.blocks_delay_0, 0), (self.blocks_float_to_complex_0, 1))
self.connect((self.blocks_float_to_complex_0, 0), (self.rational_resampler_xxx_0, 0))
self.connect((self.rational_resampler_xxx_0, 0), (self.blocks_multiply_xx_0_0, 0))
self.connect((self.analog_sig_source_x_0, 0), (self.blocks_multiply_xx_0_0, 1))
self.connect((self.blocks_multiply_xx_0_0, 0), (self.blocks_complex_to_real_0, 0))
self.connect((self.blocks_complex_to_real_0, 0), (self.audio_sink_0, 0))
def __init__(self, syncword_threshold=None, options=None):
gr.hier_block2.__init__(self, "lilacsat_1_deframer",
gr.io_signature(1, 1, gr.sizeof_float),
gr.io_signature(0, 0, 0))
options_block.__init__(self, options)
self.message_port_register_hier_out('out')
self.message_port_register_hier_out('codec2')
if syncword_threshold is None:
syncword_threshold = self.options.syncword_threshold
self.delay1 = blocks.delay(gr.sizeof_float, 1)
self.viterbi0 = ccsds_viterbi()
self.viterbi1 = ccsds_viterbi()
self.differential0 = digital.diff_decoder_bb(2)
self.differential1 = digital.diff_decoder_bb(2)
self.tag0 = digital.correlate_access_code_tag_bb(
_syncword, syncword_threshold, 'syncword')
self.tag1 = digital.correlate_access_code_tag_bb(
_syncword, syncword_threshold, 'syncword')
self.scrambler0 = digital.additive_scrambler_bb(
0xA9, 0xFF, 7, count=0, bits_per_byte=1, reset_tag_key='syncword')
self.scrambler1 = digital.additive_scrambler_bb(
0xA9, 0xFF, 7, count=0, bits_per_byte=1, reset_tag_key='syncword')
self.demux0 = lilacsat1_demux('syncword')
self.demux1 = lilacsat1_demux('syncword')
self.pdu2tag = blocks.pdu_to_tagged_stream(blocks.byte_t, 'packet_len')
self.pack = blocks.unpacked_to_packed_bb(1, gr.GR_MSB_FIRST)
self.multiply_length = blocks.tagged_stream_multiply_length(
gr.sizeof_char * 1, 'packet_len', 1 / 8.0)
self.tag2pdu = blocks.tagged_stream_to_pdu(blocks.byte_t, 'packet_len')
self.pdu2tag_codec2 = blocks.pdu_to_tagged_stream(
blocks.byte_t, 'packet_len')
self.pack_codec2 = blocks.unpacked_to_packed_bb(1, gr.GR_MSB_FIRST)
self.multiply_length_codec2 = blocks.tagged_stream_multiply_length(
gr.sizeof_char * 1, 'packet_len', 1 / 8.0)
self.tag2pdu_codec2 = blocks.tagged_stream_to_pdu(
blocks.byte_t, 'packet_len')
self.connect(self, self.viterbi0, self.differential0, self.tag0,
self.scrambler0, self.demux0)
self.connect(self, self.delay1, self.viterbi1, self.differential1,
self.tag1, self.scrambler1, self.demux1)
self.connect(self.pdu2tag, self.pack, self.multiply_length,
self.tag2pdu)
self.connect(self.pdu2tag_codec2, self.pack_codec2,
self.multiply_length_codec2, self.tag2pdu_codec2)
self.msg_connect((self.demux0, 'kiss'), (self.pdu2tag, 'pdus'))
self.msg_connect((self.demux1, 'kiss'), (self.pdu2tag, 'pdus'))
self.msg_connect((self.demux0, 'codec2'),
(self.pdu2tag_codec2, 'pdus'))
self.msg_connect((self.demux1, 'codec2'),
(self.pdu2tag_codec2, 'pdus'))
self.msg_connect((self.tag2pdu, 'pdus'), (self, 'out'))
self.msg_connect((self.tag2pdu_codec2, 'pdus'), (self, 'codec2'))
def __init__(self, context, mode, angle=0.0):
gr.hier_block2.__init__(
self,
type(self).__name__,
gr.io_signature(1, 1, gr.sizeof_float * 1),
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
)
self.__angle = 0.0 # dummy statically visible value will be overwritten
# TODO: My signal level parameters are probably wrong because this signal doesn't look like a real VOR signal
vor_30 = analog.sig_source_f(self.__audio_rate, analog.GR_COS_WAVE,
self.__vor_sig_freq, 1, 0)
vor_add = blocks.add_cc(1)
vor_audio = blocks.add_ff(1)
# Audio/AM signal
self.connect(
vor_30,
blocks.multiply_const_ff(0.3), # M_n
(vor_audio, 0))
self.connect(
self,
blocks.multiply_const_ff(audio_modulation_index), # M_i
(vor_audio, 1))
# Carrier component
self.connect(analog.sig_source_c(0, analog.GR_CONST_WAVE, 0, 0, 1),
(vor_add, 0))
# AM component
self.__delay = blocks.delay(gr.sizeof_gr_complex,
0) # configured by set_angle
self.connect(
vor_audio,
make_resampler(self.__audio_rate, self.__rf_rate
), # TODO make a complex version and do this last
blocks.float_to_complex(1),
self.__delay,
(vor_add, 1))
# FM component
vor_fm_mult = blocks.multiply_cc(1)
self.connect( # carrier generation
analog.sig_source_f(self.__rf_rate,
analog.GR_COS_WAVE, fm_subcarrier, 1, 0),
blocks.float_to_complex(1), (vor_fm_mult, 1))
self.connect( # modulation
vor_30,
make_resampler(self.__audio_rate, self.__rf_rate),
analog.frequency_modulator_fc(2 * math.pi * fm_deviation /
self.__rf_rate),
blocks.multiply_const_cc(0.3), # M_d
vor_fm_mult,
(vor_add, 2))
self.connect(vor_add, self)
# calculate and initialize delay
self.set_angle(angle)
def __init__(self):
gr.top_block.__init__(self, "Scanner Grc")
##################################################
# Variables
##################################################
self.f_symb = f_symb = 1625000.0/6.0
self.f_900_b = f_900_b = 921.2e6
self.samp_rate = samp_rate = f_symb*4
self.fs = fs = f_900_b
self.f_900_e = f_900_e = 959.8e6
self.f_1800_e = f_1800_e = 1879.8e6
self.f_1800_b = f_1800_b = 1805.2e6
self.OSR = OSR = 4
##################################################
# Blocks
##################################################
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(fs, 0)
self.osmosdr_source_0.set_freq_corr(0, 0)
self.osmosdr_source_0.set_dc_offset_mode(0, 0)
self.osmosdr_source_0.set_iq_balance_mode(2, 0)
self.osmosdr_source_0.set_gain_mode(True, 0)
self.osmosdr_source_0.set_gain(30, 0)
self.osmosdr_source_0.set_if_gain(30, 0)
self.osmosdr_source_0.set_bb_gain(30, 0)
self.osmosdr_source_0.set_antenna("", 0)
self.osmosdr_source_0.set_bandwidth(200000, 0)
self.low_pass_filter_0 = filter.fir_filter_ccf(1, firdes.low_pass(
1, samp_rate, 200e3, 10e3, firdes.WIN_HAMMING, 6.76))
self.threshold_result = blocks.threshold_ff(0, 0.2, 0)
self.blocks_threshold_ff_0_0 = blocks.threshold_ff(0, 0, 0)
self.blocks_threshold_ff_0 = blocks.threshold_ff(int((138)*samp_rate/f_symb), int((138)*samp_rate/f_symb), 0)
self.blocks_null_sink_0 = blocks.null_sink(gr.sizeof_float*1)
self.blocks_multiply_conjugate_cc_0 = blocks.multiply_conjugate_cc(1)
self.blocks_moving_average_xx_0 = blocks.moving_average_ff(int((142)*samp_rate/f_symb), 1, int(1e6))
self.blocks_delay_0 = blocks.delay(gr.sizeof_gr_complex*1, int(OSR))
self.blocks_complex_to_arg_0 = blocks.complex_to_arg(1)
##################################################
# Connections
##################################################
self.connect((self.blocks_complex_to_arg_0, 0), (self.blocks_threshold_ff_0_0, 0))
self.connect((self.blocks_delay_0, 0), (self.blocks_multiply_conjugate_cc_0, 1))
self.connect((self.blocks_moving_average_xx_0, 0), (self.blocks_threshold_ff_0, 0))
self.connect((self.blocks_multiply_conjugate_cc_0, 0), (self.blocks_complex_to_arg_0, 0))
self.connect((self.blocks_threshold_ff_0, 0), (self.threshold_result, 0))
self.connect((self.blocks_threshold_ff_0_0, 0), (self.blocks_moving_average_xx_0, 0))
self.connect((self.threshold_result, 0), (self.blocks_null_sink_0, 0))
self.connect((self.low_pass_filter_0, 0), (self.blocks_delay_0, 0))
self.connect((self.low_pass_filter_0, 0), (self.blocks_multiply_conjugate_cc_0, 0))
self.connect((self.osmosdr_source_0, 0), (self.low_pass_filter_0, 0))
def __init__ ( self, fft_length ):
gr.hier_block2.__init__(self, "recursive_timing_metric",
gr.io_signature(1,1,gr.sizeof_gr_complex),
gr.io_signature(1,1,gr.sizeof_float))
self.input = gr.kludge_copy(gr.sizeof_gr_complex)
self.connect(self, self.input)
# P(d) = sum(0 to L-1, conj(delayed(r)) * r)
conj = gr.conjugate_cc()
mixer = gr.multiply_cc()
mix_delay = delay(gr.sizeof_gr_complex,fft_length/2+1)
mix_diff = gr.sub_cc()
nominator = accumulator_cc()
inpdelay = delay(gr.sizeof_gr_complex,fft_length/2)
self.connect(self.input, inpdelay,
conj, (mixer,0))
self.connect(self.input, (mixer,1))
self.connect(mixer,(mix_diff,0))
self.connect(mixer, mix_delay, (mix_diff,1))
self.connect(mix_diff,nominator)
rmagsqrd = gr.complex_to_mag_squared()
rm_delay = delay(gr.sizeof_float,fft_length+1)
rm_diff = gr.sub_ff()
denom = accumulator_ff()
self.connect(self.input,rmagsqrd,rm_diff,gr.multiply_const_ff(0.5),denom)
self.connect(rmagsqrd,rm_delay,(rm_diff,1))
ps = gr.complex_to_mag_squared()
rs = gr.multiply_ff()
self.connect(nominator,ps)
self.connect(denom,rs)
self.connect(denom,(rs,1))
div = gr.divide_ff()
self.connect(ps,div)
self.connect(rs,(div,1))
self.connect(div,self)
def __init__(self, context, mode, angle=0.0):
gr.hier_block2.__init__(
self, 'SimulatedDevice VOR modulator',
gr.io_signature(1, 1, gr.sizeof_float * 1),
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
)
self.__angle = 0.0 # dummy statically visible value will be overwritten
# TODO: My signal level parameters are probably wrong because this signal doesn't look like a real VOR signal
vor_30 = analog.sig_source_f(self.__audio_rate, analog.GR_COS_WAVE, self.__vor_sig_freq, 1, 0)
vor_add = blocks.add_cc(1)
vor_audio = blocks.add_ff(1)
# Audio/AM signal
self.connect(
vor_30,
blocks.multiply_const_ff(0.3), # M_n
(vor_audio, 0))
self.connect(
self,
blocks.multiply_const_ff(audio_modulation_index), # M_i
(vor_audio, 1))
# Carrier component
self.connect(
analog.sig_source_c(0, analog.GR_CONST_WAVE, 0, 0, 1),
(vor_add, 0))
# AM component
self.__delay = blocks.delay(gr.sizeof_gr_complex, 0) # configured by set_angle
self.connect(
vor_audio,
make_resampler(self.__audio_rate, self.__rf_rate), # TODO make a complex version and do this last
blocks.float_to_complex(1),
self.__delay,
(vor_add, 1))
# FM component
vor_fm_mult = blocks.multiply_cc(1)
self.connect( # carrier generation
analog.sig_source_f(self.__rf_rate, analog.GR_COS_WAVE, fm_subcarrier, 1, 0),
blocks.float_to_complex(1),
(vor_fm_mult, 1))
self.connect( # modulation
vor_30,
make_resampler(self.__audio_rate, self.__rf_rate),
analog.frequency_modulator_fc(2 * math.pi * fm_deviation / self.__rf_rate),
blocks.multiply_const_cc(0.3), # M_d
vor_fm_mult,
(vor_add, 2))
self.connect(
vor_add,
self)
# calculate and initialize delay
self.set_angle(angle)
def __init__(self, fft_length):
gr.hier_block2.__init__(self, "recursive_timing_metric",
gr.io_signature(1, 1, gr.sizeof_gr_complex),
gr.io_signature(1, 1, gr.sizeof_float))
self.input = gr.kludge_copy(gr.sizeof_gr_complex)
self.connect(self, self.input)
# P(d) = sum(0 to L-1, conj(delayed(r)) * r)
conj = gr.conjugate_cc()
mixer = gr.multiply_cc()
mix_delay = delay(gr.sizeof_gr_complex, fft_length / 2 + 1)
mix_diff = gr.sub_cc()
nominator = accumulator_cc()
inpdelay = delay(gr.sizeof_gr_complex, fft_length / 2)
self.connect(self.input, inpdelay, conj, (mixer, 0))
self.connect(self.input, (mixer, 1))
self.connect(mixer, (mix_diff, 0))
self.connect(mixer, mix_delay, (mix_diff, 1))
self.connect(mix_diff, nominator)
rmagsqrd = gr.complex_to_mag_squared()
rm_delay = delay(gr.sizeof_float, fft_length + 1)
rm_diff = gr.sub_ff()
denom = accumulator_ff()
self.connect(self.input, rmagsqrd, rm_diff, gr.multiply_const_ff(0.5),
denom)
self.connect(rmagsqrd, rm_delay, (rm_diff, 1))
ps = gr.complex_to_mag_squared()
rs = gr.multiply_ff()
self.connect(nominator, ps)
self.connect(denom, rs)
self.connect(denom, (rs, 1))
div = gr.divide_ff()
self.connect(ps, div)
self.connect(rs, (div, 1))
self.connect(div, self)
def __init__(self, bias, freq_sample_delay_samps, freq_samps_to_avg, mag_samps_to_avg, resamp_rate, thresh):
gr.hier_block2.__init__(self,
"Fine Dehopper",
gr.io_signature(1, 1, gr.sizeof_gr_complex*1),
gr.io_signature(1, 1, gr.sizeof_gr_complex*1))
##################################################
# Parameters
##################################################
self.bias = bias
self.freq_sample_delay_samps = freq_sample_delay_samps
self.freq_samps_to_avg = freq_samps_to_avg
self.mag_samps_to_avg = mag_samps_to_avg
self.resamp_rate = resamp_rate
self.thresh = thresh
##################################################
# Blocks
##################################################
self.s_and_h_detector = s_and_h_detector(
freq_sample_delay_samps=freq_sample_delay_samps,
freq_samps_to_avg=freq_samps_to_avg,
mag_samps_to_avg=mag_samps_to_avg,
thresh=thresh,
)
self.resamp = pfb.arb_resampler_ccf(resamp_rate * 2.0, taps=None, flt_size=32)
self.resamp.declare_sample_delay(0)
self.fir = filter.fir_filter_ccc(2, (firdes.low_pass_2(1,1,.25,.05,60)))
self.fir.declare_sample_delay(0)
self.vco = blocks.vco_c(1, 1, 1)
self.mult_conj = blocks.multiply_conjugate_cc(1)
self.delay = blocks.delay(gr.sizeof_gr_complex*1, int(freq_samps_to_avg) + freq_sample_delay_samps)
self.c2mag = blocks.complex_to_mag(1)
self.add_const = blocks.add_const_vff((-1.0 * bias * (resamp_rate), ))
self.demod = analog.quadrature_demod_cf(1)
##################################################
# Connections
##################################################
self.connect((self.demod, 0), (self.s_and_h_detector, 1))
self.connect((self.add_const, 0), (self.vco, 0))
self.connect((self.c2mag, 0), (self.s_and_h_detector, 0))
self.connect((self.delay, 0), (self.mult_conj, 0))
self.connect((self.mult_conj, 0), (self.fir, 0))
self.connect((self.vco, 0), (self.mult_conj, 1))
self.connect((self.fir, 0), (self.resamp, 0))
self.connect((self, 0), (self.demod, 0))
self.connect((self, 0), (self.c2mag, 0))
self.connect((self, 0), (self.delay, 0))
self.connect((self.resamp, 0), (self, 0))
self.connect((self.s_and_h_detector, 0), (self.add_const, 0))
def __init__(self, payload_len=96, sync_seq=[complex(1,0),complex(-1,0)], threshold=.8):
gr.hier_block2.__init__(self,
"generic_rx_path",
gr.io_signature(1, 1, gr.sizeof_gr_complex*1),
gr.io_signaturev(3, 3, [gr.sizeof_gr_complex*1, gr.sizeof_gr_complex*1, gr.sizeof_float*1]),
)
##################################################
# Parameters
##################################################
self.threshold = threshold
self.sync_seq = sync_seq
self.sync_header_len = len(sync_seq)
self.items_per_header_symbol = 1
self.payload_len = payload_len
##################################################
# Variables
##################################################
self.tag_len_name = tag_len_name = "packet_len"
##################################################
# Blocks
##################################################
self.frame_sync_0 = dctk.frame_sync (self.sync_seq, self.threshold)
self.digital_header_payload_demux_0 = digital.header_payload_demux(self.sync_header_len, self.items_per_header_symbol, 0, tag_len_name, "", False, gr.sizeof_gr_complex)
self.dctk_header_tag_generator_0 = dctk.header_tag_generator(self.sync_header_len, payload_len, tag_len_name)
self.blocks_float_to_char_0 = blocks.float_to_char(1, 1)
self.blocks_delay_0 = blocks.delay(gr.sizeof_gr_complex*1, self.sync_header_len)
self.blocks_complex_to_float_0 = blocks.complex_to_float(1)
##################################################
# Connections
##################################################
self.connect((self, 0), (self.blocks_delay_0, 0))
self.connect((self, 0), (self.frame_sync_0, 0))
self.connect((self.blocks_delay_0, 0), (self.digital_header_payload_demux_0, 0))
self.connect((self.frame_sync_0, 1), (self.digital_header_payload_demux_0, 1))
self.connect((self.digital_header_payload_demux_0, 0), (self.blocks_complex_to_float_0, 0))
self.connect((self.blocks_complex_to_float_0, 0), (self.blocks_float_to_char_0, 0))
self.connect((self.blocks_float_to_char_0, 0), (self.dctk_header_tag_generator_0, 0))
self.connect((self.frame_sync_0, 0), (self, 2))
self.connect((self.digital_header_payload_demux_0, 1), (self, 1))
self.connect((self.digital_header_payload_demux_0, 0), (self, 0))
##################################################
# Asynch Message Connections
##################################################
self.msg_connect(self.dctk_header_tag_generator_0, "header_data", self.digital_header_payload_demux_0, "header_data")
def __init__(self, OSR=4):
grgsm.hier_block.__init__(
self,
"FCCH bursts detector",
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
)
##################################################
# Parameters
##################################################
self.OSR = OSR
##################################################
# Variables
##################################################
self.f_symb = f_symb = 1625000.0 / 6.0
self.samp_rate = samp_rate = f_symb * OSR
##################################################
# Blocks
##################################################
self.gsm_fcch_burst_tagger_0 = grgsm.fcch_burst_tagger(OSR)
self.blocks_threshold_ff_0_0 = blocks.threshold_ff(0, 0, 0)
self.blocks_threshold_ff_0 = blocks.threshold_ff(
int((138) * samp_rate / f_symb), int((138) * samp_rate / f_symb),
0)
self.blocks_multiply_conjugate_cc_0 = blocks.multiply_conjugate_cc(1)
self.blocks_moving_average_xx_0 = blocks.moving_average_ff(
int((142) * samp_rate / f_symb), 1, int(1e6))
self.blocks_delay_0 = blocks.delay(gr.sizeof_gr_complex * 1, int(OSR))
self.blocks_complex_to_arg_0 = blocks.complex_to_arg(1)
##################################################
# Connections
##################################################
self.connect((self, 0), (self.blocks_multiply_conjugate_cc_0, 0))
self.connect((self.blocks_delay_0, 0),
(self.blocks_multiply_conjugate_cc_0, 1))
self.connect((self.blocks_complex_to_arg_0, 0),
(self.blocks_threshold_ff_0_0, 0))
self.connect((self, 0), (self.blocks_delay_0, 0))
self.connect((self.blocks_multiply_conjugate_cc_0, 0),
(self.blocks_complex_to_arg_0, 0))
self.connect((self.blocks_moving_average_xx_0, 0),
(self.blocks_threshold_ff_0, 0))
self.connect((self.blocks_threshold_ff_0_0, 0),
(self.blocks_moving_average_xx_0, 0))
self.connect((self.gsm_fcch_burst_tagger_0, 0), (self, 0))
self.connect((self, 0), (self.gsm_fcch_burst_tagger_0, 0))
self.connect((self.blocks_threshold_ff_0, 0),
(self.gsm_fcch_burst_tagger_0, 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, preamble_length):
gr.hier_block2.__init__(self, "burstsilence_corr",
gr.io_signature(1, 1, gr.sizeof_gr_complex),
gr.io_signature(2, 2, gr.sizeof_gr_complex))
half_length = int(preamble_length) / 2
# Blocks
self.equiv_delay = blocks.delay(gr.sizeof_gr_complex, preamble_length)
self.to_power = blocks.complex_to_mag_squared(1)
self.pw_delay = blocks.delay(gr.sizeof_float, half_length)
self.sub_sigref = blocks.sub_ff(1)
self.avg = blocks.moving_average_ff(half_length, 1, 4096)
self.null_src = blocks.null_source(gr.sizeof_float)
self.clamp = blocks.max_ff(1, 1)
self.to_complex = blocks.float_to_complex(1)
# Connections
self.connect((self, 0), (self.to_power, 0))
self.connect((self, 0), (self.equiv_delay, 0))
self.connect((self.to_power, 0), (self.pw_delay, 0))
self.connect((self.to_power, 0), (self.sub_sigref, 1))
self.connect((self.pw_delay, 0), (self.sub_sigref, 0))
self.connect((self.sub_sigref, 0), (self.avg, 0))
self.connect((self.avg, 0), (self.clamp, 0))
self.connect((self.null_src, 0), (self.clamp, 1))
self.connect((self.clamp, 0), (self.to_complex, 0))
self.connect((self.null_src, 0), (self.to_complex, 1))
self.connect((self.equiv_delay, 0), (self, 0))
self.connect((self.to_complex, 0), (self, 1))
def test_010(self):
delta_t = 10
tb = self.tb
src_data = [float(x) for x in range(0, 100)]
expected_result = tuple(delta_t * [0.0] + src_data)
src = blocks.vector_source_f(src_data)
op = blocks.delay(gr.sizeof_float, delta_t)
dst = blocks.vector_sink_f()
tb.connect(src, op, dst)
tb.run()
dst_data = dst.data()
self.assertEqual(expected_result, dst_data)
def test_010(self):
delta_t = 10
tb = self.tb
src_data = [float(x) for x in range(0, 100)]
expected_result = tuple(delta_t*[0.0] + src_data)
src = blocks.vector_source_f(src_data)
op = blocks.delay(gr.sizeof_float, delta_t)
dst = blocks.vector_sink_f()
tb.connect(src, op, dst)
tb.run()
dst_data = dst.data()
self.assertEqual(expected_result, dst_data)
def __init__(self, carrier_freq_hi=1750, carrier_freq_lo=1080, input_path_hi='/tmp/remote.bin', input_path_lo='/tmp/local.bin', lo_delay=180):
gr.top_block.__init__(self, "Dual Fsk Modulation")
##################################################
# Parameters
##################################################
self.carrier_freq_hi = carrier_freq_hi
self.carrier_freq_lo = carrier_freq_lo
self.input_path_hi = input_path_hi
self.input_path_lo = input_path_lo
self.lo_delay = lo_delay
##################################################
# Variables
##################################################
self.oversample = oversample = 1
self.samp_rate = samp_rate = 44100
self.baud = baud = 300/oversample
self.fsk_deviation_hz = fsk_deviation_hz = 100
self.SPS = SPS = samp_rate/baud
##################################################
# Blocks
##################################################
self.digital_chunks_to_symbols_xx_0_0_0 = digital.chunks_to_symbols_bf(((2*3.14*carrier_freq_lo-2*3.14*fsk_deviation_hz,2*3.14*carrier_freq_lo+2*3.14*fsk_deviation_hz)), 1)
self.digital_chunks_to_symbols_xx_0_0 = digital.chunks_to_symbols_bf(((2*3.14*carrier_freq_hi-2*3.14*fsk_deviation_hz,2*3.14*carrier_freq_hi+2*3.14*fsk_deviation_hz)), 1)
self.blocks_wavfile_sink_0 = blocks.wavfile_sink('/tmp/data.wav', 1, samp_rate, 16)
self.blocks_vco_f_0_0 = blocks.vco_f(samp_rate*oversample, 1, 0.25)
self.blocks_vco_f_0 = blocks.vco_f(samp_rate*oversample, 1, 0.25)
self.blocks_repeat_0_0 = blocks.repeat(gr.sizeof_float*1, SPS)
self.blocks_repeat_0 = blocks.repeat(gr.sizeof_float*1, SPS)
self.blocks_file_source_0_0 = blocks.file_source(gr.sizeof_char*1, input_path_lo, False)
self.blocks_file_source_0 = blocks.file_source(gr.sizeof_char*1, input_path_hi, False)
self.blocks_delay_0 = blocks.delay(gr.sizeof_float*1, lo_delay * SPS)
self.blocks_add_xx_0 = blocks.add_vff(1)
##################################################
# Connections
##################################################
self.connect((self.blocks_add_xx_0, 0), (self.blocks_wavfile_sink_0, 0))
self.connect((self.blocks_delay_0, 0), (self.blocks_add_xx_0, 1))
self.connect((self.blocks_file_source_0, 0), (self.digital_chunks_to_symbols_xx_0_0, 0))
self.connect((self.blocks_file_source_0_0, 0), (self.digital_chunks_to_symbols_xx_0_0_0, 0))
self.connect((self.blocks_repeat_0, 0), (self.blocks_vco_f_0, 0))
self.connect((self.blocks_repeat_0_0, 0), (self.blocks_vco_f_0_0, 0))
self.connect((self.blocks_vco_f_0, 0), (self.blocks_add_xx_0, 0))
self.connect((self.blocks_vco_f_0_0, 0), (self.blocks_delay_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0_0, 0), (self.blocks_repeat_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0_0_0, 0), (self.blocks_repeat_0_0, 0))
def publish_tm_window(self,unique_id):
"""
corbaname: ofdm_ti.unique_id
"""
raise SystemError,"Bad guy! Obey the gnuradio hierarchy ..."
config = self.config
msgq = gr.msg_queue(10)
msg_sink = gr.message_sink(gr.sizeof_float*config.fft_length,msgq,True)
sampler = vector_sampler(gr.sizeof_float,config.fft_length)
self.connect(self.receiver.timing_metric,(sampler,0))
self.connect(self.receiver.time_sync,delay(gr.sizeof_char,config.fft_length/2-1),(sampler,1))
self.connect(sampler,msg_sink)
self.servants.append(corba_data_buffer_servant(str(unique_id),config.fft_length,msgq))
def add_vor(freq, angle): compensation = math.pi / 180 * -6.5 # empirical, calibrated against VOR receiver (and therefore probably wrong) angle = angle + compensation angle = angle % (2 * math.pi) vor_sig_freq = 30 phase_shift = int(rf_rate / vor_sig_freq * (angle / (2 * math.pi))) vor_dev = 480 vor_channel = make_channel(freq) vor_30 = analog.sig_source_f(audio_rate, analog.GR_COS_WAVE, vor_sig_freq, 1, 0) vor_add = blocks.add_cc(1) vor_audio = blocks.add_ff(1) # Audio/AM signal self.connect( vor_30, blocks.multiply_const_ff(0.3), # M_n (vor_audio, 0)) self.connect(audio_signal, blocks.multiply_const_ff(0.07), # M_i (vor_audio, 1)) # Carrier component self.connect( analog.sig_source_c(0, analog.GR_CONST_WAVE, 0, 0, 1), (vor_add, 0)) # AM component self.connect( vor_audio, blocks.float_to_complex(1), make_interpolator(), blocks.delay(gr.sizeof_gr_complex, phase_shift), (vor_add, 1)) # FM component vor_fm_mult = blocks.multiply_cc(1) self.connect( # carrier generation analog.sig_source_f(rf_rate, analog.GR_COS_WAVE, 9960, 1, 0), blocks.float_to_complex(1), (vor_fm_mult, 1)) self.connect( # modulation vor_30, filter.interp_fir_filter_fff(interp, interp_taps), # float not complex analog.frequency_modulator_fc(2 * math.pi * vor_dev / rf_rate), blocks.multiply_const_cc(0.3), # M_d vor_fm_mult, (vor_add, 2)) self.connect( vor_add, vor_channel) signals.append(vor_channel)
def __init__(self, OSR=4):
gr.hier_block2.__init__(
self,
"FCCH bursts detector",
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
)
##################################################
# Parameters
##################################################
self.OSR = OSR
##################################################
# Variables
##################################################
self.f_symb = f_symb = 1625000.0 / 6.0
self.samp_rate = samp_rate = f_symb * OSR
##################################################
# Blocks
##################################################
self.gsm_fcch_burst_tagger_0 = grgsm.fcch_burst_tagger(OSR)
self.blocks_threshold_ff_0_0 = blocks.threshold_ff(0, 0, 0)
self.blocks_threshold_ff_0 = blocks.threshold_ff(
int((138) * samp_rate / f_symb), int((138) * samp_rate / f_symb), 0
)
self.blocks_multiply_conjugate_cc_0 = blocks.multiply_conjugate_cc(1)
self.blocks_moving_average_xx_0 = blocks.moving_average_ff(int((142) * samp_rate / f_symb), 1, int(1e6))
self.blocks_delay_0 = blocks.delay(gr.sizeof_gr_complex * 1, int(OSR))
self.blocks_complex_to_arg_0 = blocks.complex_to_arg(1)
##################################################
# Connections
##################################################
self.connect((self, 0), (self.blocks_multiply_conjugate_cc_0, 0))
self.connect((self.blocks_delay_0, 0), (self.blocks_multiply_conjugate_cc_0, 1))
self.connect((self.blocks_complex_to_arg_0, 0), (self.blocks_threshold_ff_0_0, 0))
self.connect((self, 0), (self.blocks_delay_0, 0))
self.connect((self.blocks_multiply_conjugate_cc_0, 0), (self.blocks_complex_to_arg_0, 0))
self.connect((self.blocks_moving_average_xx_0, 0), (self.blocks_threshold_ff_0, 0))
self.connect((self.blocks_threshold_ff_0_0, 0), (self.blocks_moving_average_xx_0, 0))
self.connect((self.gsm_fcch_burst_tagger_0, 0), (self, 0))
self.connect((self, 0), (self.gsm_fcch_burst_tagger_0, 0))
self.connect((self.blocks_threshold_ff_0, 0), (self.gsm_fcch_burst_tagger_0, 1))
def oqpsktx(self, carrier=10000, samp_rate = 80000, bw=4000, amp=1, code=codes.mycode, **kwargs):
code_table, code_len = codes.codes2table(code), len(code)
chunk_len = int(log(code_len,2))
topblock(self, carrier, samp_rate, bw, amp)
##################################################
# Blocks
##################################################
self.blocks_packed_to_unpacked_xx_0 = blocks.packed_to_unpacked_bb(chunk_len, gr.GR_LSB_FIRST)
self.digital_chunks_to_symbols_xx_0 = digital.chunks_to_symbols_bc((code_table), code_len)
self.blocks_repeat_0 = blocks.repeat(gr.sizeof_gr_complex*1, 4)
self.blocks_vector_source_x_0 = blocks.vector_source_c([0, sin(pi/4), 1, sin(3*pi/4)], True, 1, [])
self.blocks_multiply_xx_0 = blocks.multiply_vcc(1)
self.blocks_complex_to_float_0 = blocks.complex_to_float(1)
self.blocks_delay_0 = blocks.delay(gr.sizeof_float*1, 2)
self.blocks_float_to_complex_0 = blocks.float_to_complex(1)
self.rational_resampler_xxx_0 = filter.rational_resampler_ccc(
interpolation=int(samp_rate/bw),
decimation=1,
taps=None,
fractional_bw=None,
)
##################################################
# Connections
##################################################
self.connect((self.source, 0), (self.blocks_packed_to_unpacked_xx_0, 0))
self.connect((self.blocks_packed_to_unpacked_xx_0, 0), (self.digital_chunks_to_symbols_xx_0, 0))
self.connect((self.digital_chunks_to_symbols_xx_0, 0), (self.blocks_repeat_0, 0))
self.connect((self.blocks_vector_source_x_0, 0), (self.blocks_multiply_xx_0, 0))
self.connect((self.blocks_repeat_0, 0), (self.blocks_multiply_xx_0, 1))
self.connect((self.blocks_multiply_xx_0, 0), (self.blocks_complex_to_float_0, 0))
self.connect((self.blocks_complex_to_float_0, 0), (self.blocks_float_to_complex_0, 0))
self.connect((self.blocks_complex_to_float_0, 1), (self.blocks_delay_0, 0))
self.connect((self.blocks_delay_0, 0), (self.blocks_float_to_complex_0, 1))
self.connect((self.blocks_float_to_complex_0, 0), (self.rational_resampler_xxx_0, 0))
self.connect((self.rational_resampler_xxx_0, 0), (self.sink, 0))
def __init__(self, src_file):
gr.top_block.__init__(self)
sample_rate = 11025
ampl = 0.1
print src_file
# Audio source (.wav file)
# TODO : Make the filename a VARIABLE
# src_file = input("Enter .wav File PSK31 : ")
src = blocks.wavfile_source(src_file, False)
# Raw float data output file.
# TODO : To make the raw file also a variable, for psk31decoder2.py to run
dst = blocks.file_sink(1, "./output.raw")
# Delay line. This delays the signal by 32ms
dl = blocks.delay(gr.sizeof_float, int(round(sample_rate/31.25)))
# Multiplier
# Multiplying the source and the delayed version will give us
# a negative output if there was a phase reversal and a positive output
# if there was no phase reversal
mul = blocks.multiply_ff(1)
# Low Pass Filter. This leaves us with the envelope of the signal
lpf_taps = filter.firdes.low_pass(
5.0,
sample_rate,
15,
600,
filter.firdes.WIN_HAMMING)
lpf = filter.fir_filter_fff(1, lpf_taps)
# Binary Slicer (comparator)
slc = digital.binary_slicer_fb()
# Connect the blocks.
self.connect(src, dl)
self.connect(src, (mul, 0))
self.connect(dl, (mul, 1))
self.connect(mul, lpf)
self.connect(lpf, slc)
self.connect(slc, dst)
def test_020(self):
tb = self.tb
vector_sink = blocks.vector_sink_f(1)
ref_sink = blocks.vector_sink_f(1)
tags_strobe = blocks.tags_strobe(gr.sizeof_float*1, pmt.intern("TEST"), 100, pmt.intern("strobe"))
head = blocks.head(gr.sizeof_float*1, 10**5)
delay = blocks.delay(gr.sizeof_float*1, 100)
tb.connect((delay, 0), (head, 0))
tb.connect((head, 0), (vector_sink, 0))
tb.connect((tags_strobe, 0), (delay, 0))
tb.connect((tags_strobe, 0), (ref_sink, 0))
tb.run()
tags = vector_sink.tags()
self.assertNotEqual(len(tags), 0)
lastoffset = tags[0].offset - 100
for tag in tags:
newoffset = tag.offset
self.assertEqual(newoffset, lastoffset + 100)
lastoffset = newoffset
def __init__(self,options):
config = station_configuration()
total_subc = config.subcarriers
vlen = total_subc
gr.hier_block2.__init__(self,"fbmc_frame_sampler",
gr.io_signature2(2,2,gr.sizeof_gr_complex*vlen,
gr.sizeof_char),
gr.io_signature2(2,2,gr.sizeof_gr_complex*vlen,
gr.sizeof_char))
frame_size = config.frame_data_part + config.training_data.fbmc_no_preambles/2
print "frame_size: ", frame_size
ft = [0] * frame_size
ft[0] = 1
# The next block ensures that only complete frames find their way into
# the old outer receiver. The dynamic frame start trigger is hence
# replaced with a static one, fixed to the frame length.
frame_sampler = ofdm.vector_sampler( gr.sizeof_gr_complex * total_subc,
frame_size )
symbol_output = blocks.vector_to_stream( gr.sizeof_gr_complex * total_subc,
frame_size )
#delayed_frame_start = blocks.delay( gr.sizeof_char, config.frame_length-1 - config.training_data.fbmc_no_preambles/2 )
delayed_frame_start = blocks.delay( gr.sizeof_char, config.frame_length/2-1)
damn_static_frame_trigger = blocks.vector_source_b( ft, True )
#oqam_postpro = ofdm.fbmc_oqam_postprocessing_vcvc(total_subc,0,0)
self.connect( self, frame_sampler, symbol_output ,self)
#self.connect( (self,1), blocks.keep_m_in_n(gr.sizeof_char,config.frame_data_part + config.training_data.fbmc_no_preambles/2,2*config.frame_data_part + config.training_data.fbmc_no_preambles,0),delayed_frame_start, ( frame_sampler, 1 ) )
self.connect( (self,1), delayed_frame_start, ( frame_sampler, 1 ) )
#self.connect( self, blocks.multiply_const_vcc(([1.0]*total_subc)) ,self)
#terminate_stream(self,frame_sampler)
self.connect( damn_static_frame_trigger, (self,1) )
def __init__(self, fft_length):
gr.hier_block2.__init__(self, "schmidl_nominator",
gr.io_signature(1,1,gr.sizeof_gr_complex),
gr.io_signature(1,1,gr.sizeof_gr_complex))
self.input=gr.kludge_copy(gr.sizeof_gr_complex)
# P(d) = sum(0 to L-1, conj(delayed(r)) * r)
conj = gr.conjugate_cc()
mixer = gr.multiply_cc()
moving_avg = gr.fir_filter_ccf(1,[1.0 for i in range(fft_length/2)])
self.connect(self, self.input, delay(gr.sizeof_gr_complex,fft_length/2), conj, (mixer,0))
self.connect(self.input, (mixer,1))
self.connect(mixer, moving_avg, self)
# moving_avg = P(d)
try:
gr.hier_block.update_var_names(self, "schmidl_nom", vars())
gr.hier_block.update_var_names(self, "schmidl_nom", vars(self))
except:
pass
def __init__(self,options):
config = station_configuration()
total_subc = config.subcarriers
vlen = total_subc
gr.hier_block2.__init__(self,"ofdm_frame_sampler",
gr.io_signature2(2,2,gr.sizeof_gr_complex*vlen,
gr.sizeof_char),
gr.io_signature2(2,2,gr.sizeof_gr_complex*vlen,
gr.sizeof_char))
ft = [0] * config.frame_length
ft[0] = 1
# The next block ensures that only complete frames find their way into
# the old outer receiver. The dynamic frame start trigger is hence
# replaced with a static one, fixed to the frame length.
frame_sampler = ofdm.vector_sampler( gr.sizeof_gr_complex * total_subc,
config.frame_length )
symbol_output = blocks.vector_to_stream( gr.sizeof_gr_complex * total_subc,
config.frame_length )
delayed_frame_start = blocks.delay( gr.sizeof_char, config.frame_length - 1 )
damn_static_frame_trigger = blocks.vector_source_b( ft, True )
#oqam_postpro = ofdm.fbmc_oqam_postprocessing_vcvc(total_subc,0,0)
if options.enable_erasure_decision:
self.frame_gate = vector_sampler(
gr.sizeof_gr_complex * total_subc * config.frame_length, 1 )
self.connect( self, frame_sampler, self.frame_gate,
symbol_output )
else:
self.connect( self, frame_sampler, symbol_output, self )
self.connect( (self,1), delayed_frame_start, ( frame_sampler, 1 ) )
self.connect( damn_static_frame_trigger, (self,1) )
def __init__(self, fft_len, cp_len, nofdm_symbols):
gr.hier_block2.__init__(
self,
"ofdm_basebandsignal_to_frames_cvc",
gr.io_signature(1, 1, gr.sizeof_gr_complex),
gr.io_signature(1, 1, fft_len*gr.sizeof_gr_complex)
)
self.fft_len = fft_len
self.cp_len = cp_len
self.nofdm_symbols = nofdm_symbols
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
)
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))
self.connect(frames, self)
def __init__(self):
gr.top_block.__init__(self, "Wifi Rx")
Qt.QWidget.__init__(self)
self.setWindowTitle("Wifi Rx")
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", "wifi_rx")
self.restoreGeometry(self.settings.value("geometry").toByteArray())
##################################################
# Variables
##################################################
self.window_size = window_size = 48
self.sync_length = sync_length = 320
self.samp_rate = samp_rate = 5e6
self.lo_offset = lo_offset = 0
self.gain = gain = 20
self.freq = freq = 943000000.0
self.chan_est = chan_est = 1
##################################################
# Blocks
##################################################
self._samp_rate_options = [1e6, 5e6, 10e6, 20e6]
self._samp_rate_labels = ["1 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"+": "))
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.top_layout.addWidget(self._samp_rate_tool_bar)
self.qtgui_time_sink_x_1_0_0_0 = qtgui.time_sink_f(
1024*16*4, #size
samp_rate, #samp_rate
"", #name
1 #number of inputs
)
self.qtgui_time_sink_x_1_0_0_0.set_update_time(0.10)
self.qtgui_time_sink_x_1_0_0_0.set_y_axis(-0.1, 1.4)
self.qtgui_time_sink_x_1_0_0_0.set_y_label("Amplitude ----", "")
self.qtgui_time_sink_x_1_0_0_0.enable_tags(-1, True)
self.qtgui_time_sink_x_1_0_0_0.set_trigger_mode(qtgui.TRIG_MODE_NORM, qtgui.TRIG_SLOPE_POS, 0.63, 0, 0, "")
self.qtgui_time_sink_x_1_0_0_0.enable_autoscale(False)
self.qtgui_time_sink_x_1_0_0_0.enable_grid(False)
self.qtgui_time_sink_x_1_0_0_0.enable_control_panel(False)
if not True:
self.qtgui_time_sink_x_1_0_0_0.disable_legend()
labels = ["RATIO", "", "", "", "",
"", "", "", "", ""]
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 = [2.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_0.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_time_sink_x_1_0_0_0.set_line_label(i, labels[i])
self.qtgui_time_sink_x_1_0_0_0.set_line_width(i, widths[i])
self.qtgui_time_sink_x_1_0_0_0.set_line_color(i, colors[i])
self.qtgui_time_sink_x_1_0_0_0.set_line_style(i, styles[i])
self.qtgui_time_sink_x_1_0_0_0.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_1_0_0_0.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_1_0_0_0_win = sip.wrapinstance(self.qtgui_time_sink_x_1_0_0_0.pyqwidget(), Qt.QWidget)
self.top_layout.addWidget(self._qtgui_time_sink_x_1_0_0_0_win)
self.qtgui_time_sink_x_1_0_0 = qtgui.time_sink_f(
1024*16*4, #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(-0.001, 0.2)
self.qtgui_time_sink_x_1_0_0.set_y_label("Amplitude", "")
self.qtgui_time_sink_x_1_0_0.enable_tags(-1, True)
self.qtgui_time_sink_x_1_0_0.set_trigger_mode(qtgui.TRIG_MODE_NORM, qtgui.TRIG_SLOPE_POS, 0.05, 0, 0, "")
self.qtgui_time_sink_x_1_0_0.enable_autoscale(True)
self.qtgui_time_sink_x_1_0_0.enable_grid(False)
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 = ["Power av.", "", "", "", "",
"", "", "", "", ""]
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.top_layout.addWidget(self._qtgui_time_sink_x_1_0_0_win)
self.qtgui_time_sink_x_1_0 = qtgui.time_sink_f(
1024*16*4, #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(-0.001, 0.2)
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_NORM, qtgui.TRIG_SLOPE_POS, 0.05, 0, 0, "")
self.qtgui_time_sink_x_1_0.enable_autoscale(True)
self.qtgui_time_sink_x_1_0.enable_grid(False)
self.qtgui_time_sink_x_1_0.enable_control_panel(False)
if not True:
self.qtgui_time_sink_x_1_0.disable_legend()
labels = ["Correlation av.", "", "", "", "",
"", "", "", "", ""]
widths = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
colors = ["blue", "red", "green", "black", "cyan",
"magenta", "yellow", "dark red", "dark green", "blue"]
styles = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1,
-1, -1, -1, -1, -1]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_time_sink_x_1_0.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_time_sink_x_1_0.set_line_label(i, labels[i])
self.qtgui_time_sink_x_1_0.set_line_width(i, widths[i])
self.qtgui_time_sink_x_1_0.set_line_color(i, colors[i])
self.qtgui_time_sink_x_1_0.set_line_style(i, styles[i])
self.qtgui_time_sink_x_1_0.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_1_0.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_1_0_win = sip.wrapinstance(self.qtgui_time_sink_x_1_0.pyqwidget(), Qt.QWidget)
self.top_layout.addWidget(self._qtgui_time_sink_x_1_0_win)
self.nutaq_rtdex_source_0 = nutaq.rtdex_source("nutaq_carrier_perseus_0",gr.sizeof_short,1,3)
self.nutaq_rtdex_source_0.set_type(0)
self.nutaq_rtdex_source_0.set_packet_size(8192)
self.nutaq_rtdex_source_0.set_channels("1")
self.nutaq_radio420_tx_0_0 = nutaq.radio420_tx("nutaq_carrier_perseus_0", 1, 0)
self.nutaq_radio420_tx_0_0.set_default_enable(1)
self.nutaq_radio420_tx_0_0.set_default_tx_freq(943e6)
self.nutaq_radio420_tx_0_0.set_default_reference(0)
self.nutaq_radio420_tx_0_0.set_default_datarate(samp_rate*2)
self.nutaq_radio420_tx_0_0.set_default_calibrate(1)
self.nutaq_radio420_tx_0_0.set_default_band(0)
self.nutaq_radio420_tx_0_0.set_default_update_rate(1)
self.nutaq_radio420_tx_0_0.set_default_tx_vga1_gain(-10)
self.nutaq_radio420_tx_0_0.set_default_tx_vga2_gain(15)
self.nutaq_radio420_tx_0_0.set_default_tx_gain3(3)
self.nutaq_radio420_tx_0_0.set_default_tx_lpf_bandwidth(6)
self.nutaq_radio420_tx_0_0.set_default_ref_clk_ctrl(0)
self.nutaq_radio420_tx_0_0.set_default_rf_ctrl(0)
self.nutaq_radio420_tx_0_0.set_default_tx_gain_ctrl(0)
self.nutaq_radio420_tx_0_0.set_default_pll_cpld_ctrl(0)
self.nutaq_radio420_rx_0 = nutaq.radio420_rx("nutaq_carrier_perseus_0", 1, 1)
self.nutaq_radio420_rx_0.set_default_enable(1)
self.nutaq_radio420_rx_0.set_default_rx_freq(943e6)
self.nutaq_radio420_rx_0.set_default_reference(0)
self.nutaq_radio420_rx_0.set_default_datarate(samp_rate*2)
self.nutaq_radio420_rx_0.set_default_calibrate(1)
self.nutaq_radio420_rx_0.set_default_band(0)
self.nutaq_radio420_rx_0.set_default_update_rate(1)
self.nutaq_radio420_rx_0.set_default_rx_lna_gain(2)
self.nutaq_radio420_rx_0.set_default_rx_vga1_gain(1)
self.nutaq_radio420_rx_0.set_default_rx_gain2(16)
self.nutaq_radio420_rx_0.set_default_rx_gain3(5)
self.nutaq_radio420_rx_0.set_default_rx_rf_filter(2)
self.nutaq_radio420_rx_0.set_default_rx_lpf_bandwidth(2)
self.nutaq_radio420_rx_0.set_default_ref_clk_ctrl(0)
self.nutaq_radio420_rx_0.set_default_rf_ctrl(0)
self.nutaq_radio420_rx_0.set_default_rx_gain_ctrl(0)
self.nutaq_radio420_rx_0.set_default_pll_cpld_ctrl(0)
self.nutaq_custom_register_0_0 = nutaq.custom_register("nutaq_carrier_perseus_0",5)
self.nutaq_custom_register_0_0.set_index(4)
self.nutaq_custom_register_0_0.set_update_rate(1)
self.nutaq_custom_register_0 = nutaq.custom_register("nutaq_carrier_perseus_0",4)
self.nutaq_custom_register_0.set_index(1)
self.nutaq_custom_register_0.set_default_value(6)
self.nutaq_custom_register_0.set_update_rate(1)
self.nutaq_carrier_perseus_0 = nutaq.carrier(0,"nutaq_carrier_perseus_0", "192.168.0.103")
self.ieee802_11_ofdm_sync_short_0 = ieee802_11.ofdm_sync_short(0.85 * 0.75, 2, False, False)
self.ieee802_11_ofdm_sync_long_0 = ieee802_11.ofdm_sync_long(sync_length, False, False)
self.ieee802_11_ofdm_equalize_symbols_0 = ieee802_11.ofdm_equalize_symbols(ieee802_11.LMS, False)
self.ieee802_11_ofdm_decode_signal_0 = ieee802_11.ofdm_decode_signal(False, False)
self.ieee802_11_ofdm_decode_mac_0 = ieee802_11.ofdm_decode_mac(True, True)
self.ieee802_11_moving_average_xx_1 = ieee802_11.moving_average_ff(window_size + 16)
self.ieee802_11_moving_average_xx_0 = ieee802_11.moving_average_cc(window_size)
self._gain_range = Range(0, 100, 1, 20, 200)
self._gain_win = RangeWidget(self._gain_range, self.set_gain, "gain", "slider")
self.top_layout.addWidget(self._gain_win)
self._freq_options = [943000000.0, 2412000000.0, 2417000000.0, 2422000000.0, 2427000000.0, 2432000000.0, 2437000000.0, 2442000000.0, 2447000000.0, 2452000000.0, 2457000000.0, 2462000000.0, 2467000000.0, 2472000000.0, 2484000000.0, 5170000000.0, 5180000000.0, 5190000000.0, 5200000000.0, 5210000000.0, 5220000000.0, 5230000000.0, 5240000000.0, 5260000000.0, 5280000000.0, 5300000000.0, 5320000000.0, 5500000000.0, 5520000000.0, 5540000000.0, 5560000000.0, 5580000000.0, 5600000000.0, 5620000000.0, 5640000000.0, 5660000000.0, 5680000000.0, 5700000000.0, 5745000000.0, 5765000000.0, 5785000000.0, 5805000000.0, 5825000000.0, 5860000000.0, 5870000000.0, 5880000000.0, 5890000000.0, 5900000000.0, 5910000000.0, 5920000000.0]
self._freq_labels = [' 0 | 943.0 | ??', ' 1 | 2412.0 | 11g', ' 2 | 2417.0 | 11g', ' 3 | 2422.0 | 11g', ' 4 | 2427.0 | 11g', ' 5 | 2432.0 | 11g', ' 6 | 2437.0 | 11g', ' 7 | 2442.0 | 11g', ' 8 | 2447.0 | 11g', ' 9 | 2452.0 | 11g', ' 10 | 2457.0 | 11g', ' 11 | 2462.0 | 11g', ' 12 | 2467.0 | 11g', ' 13 | 2472.0 | 11g', ' 14 | 2484.0 | 11g', ' 34 | 5170.0 | 11a', ' 36 | 5180.0 | 11a', ' 38 | 5190.0 | 11a', ' 40 | 5200.0 | 11a', ' 42 | 5210.0 | 11a', ' 44 | 5220.0 | 11a', ' 46 | 5230.0 | 11a', ' 48 | 5240.0 | 11a', ' 52 | 5260.0 | 11a', ' 56 | 5280.0 | 11a', ' 58 | 5300.0 | 11a', ' 60 | 5320.0 | 11a', '100 | 5500.0 | 11a', '104 | 5520.0 | 11a', '108 | 5540.0 | 11a', '112 | 5560.0 | 11a', '116 | 5580.0 | 11a', '120 | 5600.0 | 11a', '124 | 5620.0 | 11a', '128 | 5640.0 | 11a', '132 | 5660.0 | 11a', '136 | 5680.0 | 11a', '140 | 5700.0 | 11a', '149 | 5745.0 | 11a', '153 | 5765.0 | 11a', '157 | 5785.0 | 11a', '161 | 5805.0 | 11a', '165 | 5825.0 | 11a', '172 | 5860.0 | 11p', '174 | 5870.0 | 11p', '176 | 5880.0 | 11p', '178 | 5890.0 | 11p', '180 | 5900.0 | 11p', '182 | 5910.0 | 11p', '184 | 5920.0 | 11p']
self._freq_tool_bar = Qt.QToolBar(self)
self._freq_tool_bar.addWidget(Qt.QLabel("Channel"+": "))
self._freq_combo_box = Qt.QComboBox()
self._freq_tool_bar.addWidget(self._freq_combo_box)
for label in self._freq_labels: self._freq_combo_box.addItem(label)
self._freq_callback = lambda i: Qt.QMetaObject.invokeMethod(self._freq_combo_box, "setCurrentIndex", Qt.Q_ARG("int", self._freq_options.index(i)))
self._freq_callback(self.freq)
self._freq_combo_box.currentIndexChanged.connect(
lambda i: self.set_freq(self._freq_options[i]))
self.top_layout.addWidget(self._freq_tool_bar)
self.foo_wireshark_connector_0 = foo.wireshark_connector(127, False)
self.fft_vxx_0 = fft.fft_vcc(64, True, (window.rectangular(64)), True, 1)
self.blocks_stream_to_vector_0 = blocks.stream_to_vector(gr.sizeof_gr_complex*1, 64)
self.blocks_short_to_float_0_0 = blocks.short_to_float(1, 2**11-1)
self.blocks_short_to_float_0 = blocks.short_to_float(1, 2**11-1)
self.blocks_null_sink_0 = blocks.null_sink(gr.sizeof_gr_complex*48)
self.blocks_multiply_xx_0 = blocks.multiply_vcc(1)
self.blocks_float_to_complex_0 = blocks.float_to_complex(1)
self.blocks_file_sink_0 = blocks.file_sink(gr.sizeof_char*1, "/tmp/ofdm.pcap", True)
self.blocks_file_sink_0.set_unbuffered(True)
self.blocks_divide_xx_0 = blocks.divide_ff(1)
self.blocks_delay_0_0 = blocks.delay(gr.sizeof_gr_complex*1, 16)
self.blocks_delay_0 = blocks.delay(gr.sizeof_gr_complex*1, sync_length)
self.blocks_deinterleave_0 = blocks.deinterleave(gr.sizeof_short*1, 1)
self.blocks_conjugate_cc_0 = blocks.conjugate_cc()
self.blocks_complex_to_mag_squared_0 = blocks.complex_to_mag_squared(1)
self.blocks_complex_to_mag_0 = blocks.complex_to_mag(1)
##################################################
# Connections
##################################################
self.msg_connect((self.ieee802_11_ofdm_decode_mac_0, 'out'), (self.foo_wireshark_connector_0, 'in'))
self.connect((self.blocks_complex_to_mag_0, 0), (self.blocks_divide_xx_0, 0))
self.connect((self.blocks_complex_to_mag_0, 0), (self.qtgui_time_sink_x_1_0, 0))
self.connect((self.blocks_complex_to_mag_squared_0, 0), (self.ieee802_11_moving_average_xx_1, 0))
self.connect((self.blocks_conjugate_cc_0, 0), (self.blocks_multiply_xx_0, 1))
self.connect((self.blocks_deinterleave_0, 0), (self.blocks_short_to_float_0, 0))
self.connect((self.blocks_deinterleave_0, 1), (self.blocks_short_to_float_0_0, 0))
self.connect((self.blocks_delay_0, 0), (self.ieee802_11_ofdm_sync_long_0, 1))
self.connect((self.blocks_delay_0_0, 0), (self.blocks_conjugate_cc_0, 0))
self.connect((self.blocks_delay_0_0, 0), (self.ieee802_11_ofdm_sync_short_0, 0))
self.connect((self.blocks_divide_xx_0, 0), (self.ieee802_11_ofdm_sync_short_0, 2))
self.connect((self.blocks_divide_xx_0, 0), (self.qtgui_time_sink_x_1_0_0_0, 0))
self.connect((self.blocks_float_to_complex_0, 0), (self.blocks_complex_to_mag_squared_0, 0))
self.connect((self.blocks_float_to_complex_0, 0), (self.blocks_delay_0_0, 0))
self.connect((self.blocks_float_to_complex_0, 0), (self.blocks_multiply_xx_0, 0))
self.connect((self.blocks_multiply_xx_0, 0), (self.ieee802_11_moving_average_xx_0, 0))
self.connect((self.blocks_short_to_float_0, 0), (self.blocks_float_to_complex_0, 0))
self.connect((self.blocks_short_to_float_0_0, 0), (self.blocks_float_to_complex_0, 1))
self.connect((self.blocks_stream_to_vector_0, 0), (self.fft_vxx_0, 0))
self.connect((self.fft_vxx_0, 0), (self.ieee802_11_ofdm_equalize_symbols_0, 0))
self.connect((self.foo_wireshark_connector_0, 0), (self.blocks_file_sink_0, 0))
self.connect((self.ieee802_11_moving_average_xx_0, 0), (self.blocks_complex_to_mag_0, 0))
self.connect((self.ieee802_11_moving_average_xx_0, 0), (self.ieee802_11_ofdm_sync_short_0, 1))
self.connect((self.ieee802_11_moving_average_xx_1, 0), (self.blocks_divide_xx_0, 1))
self.connect((self.ieee802_11_moving_average_xx_1, 0), (self.qtgui_time_sink_x_1_0_0, 0))
self.connect((self.ieee802_11_ofdm_decode_signal_0, 0), (self.blocks_null_sink_0, 0))
self.connect((self.ieee802_11_ofdm_decode_signal_0, 0), (self.ieee802_11_ofdm_decode_mac_0, 0))
self.connect((self.ieee802_11_ofdm_equalize_symbols_0, 0), (self.ieee802_11_ofdm_decode_signal_0, 0))
self.connect((self.ieee802_11_ofdm_sync_long_0, 0), (self.blocks_stream_to_vector_0, 0))
self.connect((self.ieee802_11_ofdm_sync_short_0, 0), (self.blocks_delay_0, 0))
self.connect((self.ieee802_11_ofdm_sync_short_0, 0), (self.ieee802_11_ofdm_sync_long_0, 0))
self.connect((self.nutaq_rtdex_source_0, 0), (self.blocks_deinterleave_0, 0))
def __init__(self, fft_length, cp_length, occupied_tones, snr, ks, carrier_map_bin, nc_filter, logging=False):
"""
Hierarchical block for receiving OFDM symbols.
The input is the complex modulated signal at baseband.
Synchronized packets are sent back to the demodulator.
@param fft_length: total number of subcarriers
@type fft_length: int
@param cp_length: length of cyclic prefix as specified in subcarriers (<= fft_length)
@type cp_length: int
@param occupied_tones: number of subcarriers used for data
@type occupied_tones: int
@param snr: estimated signal to noise ratio used to guide cyclic prefix synchronizer
@type snr: float
@param ks: known symbols used as preambles to each packet
@type ks: list of lists
@param logging: turn file logging on or off
@type logging: bool
"""
gr.hier_block2.__init__(self, "ofdm_receiver",
gr.io_signature(1, 1, gr.sizeof_gr_complex), # Input signature
gr.io_signature2(2, 2, gr.sizeof_gr_complex*occupied_tones, gr.sizeof_char)) # Output signature
bw = (float(occupied_tones) / float(fft_length)) / 2.0
tb = bw*0.04
print "ofdm_receiver:__init__:occupied_tones %s fft_length %d " % (occupied_tones, fft_length)
chan_coeffs = filter.firdes.low_pass (1.0, # gain
1.0, # sampling rate
bw+tb, # midpoint of trans. band
tb, # width of trans. band
filter.firdes.WIN_HAMMING) # filter type
self.chan_filt = filter.fft_filter_ccc(1, chan_coeffs)
# linklab, get ofdm parameters
self._fft_length = fft_length
self._occupied_tones = occupied_tones
self._cp_length = cp_length
self._nc_filter = nc_filter
self._carrier_map_bin = carrier_map_bin
win = [1 for i in range(fft_length)]
# linklab, initialization function
self.initialize(ks, self._carrier_map_bin)
zeros_on_left = int(math.ceil((fft_length - occupied_tones)/2.0))
ks0 = fft_length*[0,]
ks0[zeros_on_left : zeros_on_left + occupied_tones] = ks[0]
ks0 = np_fft.ifftshift(ks0)
ks0time = np_fft.ifft(ks0)
# ADD SCALING FACTOR
ks0time = ks0time.tolist()
SYNC = "pn"
if SYNC == "ml":
nco_sensitivity = -1.0/fft_length # correct for fine frequency
self.ofdm_sync = ofdm_sync_ml(fft_length,
cp_length,
snr,
ks0time,
logging)
elif SYNC == "pn":
nco_sensitivity = -2.0/fft_length # correct for fine frequency
self.ofdm_sync = ofdm_sync_pn(fft_length,
cp_length,
logging)
elif SYNC == "pnac":
nco_sensitivity = -2.0/fft_length # correct for fine frequency
self.ofdm_sync = ofdm_sync_pnac(fft_length,
cp_length,
ks0time,
logging)
# for testing only; do not user over the air
# remove filter and filter delay for this
elif SYNC == "fixed":
self.chan_filt = gr.multiply_const_cc(1.0)
nsymbols = 18 # enter the number of symbols per packet
freq_offset = 0.0 # if you use a frequency offset, enter it here
nco_sensitivity = -2.0/fft_length # correct for fine frequency
self.ofdm_sync = ofdm_sync_fixed(fft_length,
cp_length,
nsymbols,
freq_offset,
logging)
# Set up blocks
# Create a delay line, linklab
self.delay = blocks.delay(gr.sizeof_gr_complex, fft_length)
self.nco = analog.frequency_modulator_fc(nco_sensitivity) # generate a signal proportional to frequency error of sync block
self.sigmix = blocks.multiply_cc()
self.sampler = gr_papyrus.ofdm_sampler(fft_length, fft_length+cp_length)
self.fft_demod = gr_fft.fft_vcc(fft_length, True, win, True)
self.ofdm_frame_acq = gr_papyrus.ofdm_frame_acquisition(occupied_tones,
fft_length,
cp_length, ks[0])
# linklab, check current mode: non-contiguous OFDM or not
if self._nc_filter:
print '\nMulti-band Filter Turned ON!'
# linklab, non-contiguous filter
self.ncofdm_filt = ncofdm_filt(self._fft_length, self._occupied_tones, self._carrier_map_bin)
self.connect(self, self.chan_filt, self.ncofdm_filt)
self.connect(self.ncofdm_filt, self.ofdm_sync) # into the synchronization alg.
self.connect((self.ofdm_sync,0), self.nco, (self.sigmix,1)) # use sync freq. offset output to derotate input signal
self.connect(self.ncofdm_filt, self.delay, (self.sigmix,0)) # signal to be derotated
else :
print '\nMulti-band Filter Turned OFF!'
self.connect(self, self.chan_filt)
self.connect(self.chan_filt, self.ofdm_sync) # into the synchronization alg.
self.connect((self.ofdm_sync,0), self.nco, (self.sigmix,1)) # use sync freq. offset output to derotate input signal
self.connect(self.chan_filt, self.delay, (self.sigmix,0)) # signal to be derotated
self.connect(self.sigmix, (self.sampler,0)) # sample off timing signal detected in sync alg
self.connect((self.ofdm_sync,1), (self.sampler,1)) # timing signal to sample at
self.connect((self.sampler,0), self.fft_demod) # send derotated sampled signal to FFT
self.connect(self.fft_demod, (self.ofdm_frame_acq,0)) # find frame start and equalize signal
self.connect((self.sampler,1), (self.ofdm_frame_acq,1)) # send timing signal to signal frame start
self.connect((self.ofdm_frame_acq,0), (self,0)) # finished with fine/coarse freq correction,
self.connect((self.ofdm_frame_acq,1), (self,1)) # frame and symbol timing, and equalization
if logging:
self.connect(self.chan_filt, gr.file_sink(gr.sizeof_gr_complex, "ofdm_receiver-chan_filt_c.dat"))
self.connect(self.fft_demod, gr.file_sink(gr.sizeof_gr_complex*fft_length, "ofdm_receiver-fft_out_c.dat"))
self.connect(self.ofdm_frame_acq,
gr.file_sink(gr.sizeof_gr_complex*occupied_tones, "ofdm_receiver-frame_acq_c.dat"))
self.connect((self.ofdm_frame_acq,1), gr.file_sink(1, "ofdm_receiver-found_corr_b.dat"))
self.connect(self.sampler, gr.file_sink(gr.sizeof_gr_complex*fft_length, "ofdm_receiver-sampler_c.dat"))
self.connect(self.sigmix, gr.file_sink(gr.sizeof_gr_complex, "ofdm_receiver-sigmix_c.dat"))
self.connect(self.nco, gr.file_sink(gr.sizeof_gr_complex, "ofdm_receiver-nco_c.dat"))
def __init__( self, options, log = False ):
## Read configuration
config = station_configuration()
fft_length = config.fft_length
cp_length = config.cp_length
block_header = config.training_data
data_subc = config.data_subcarriers
virtual_subc = config.virtual_subcarriers
total_subc = config.subcarriers
block_length = config.block_length
frame_length = config.frame_length
dc_null = config.dc_null
L = block_header.mm_periodic_parts
## Set Input/Output signature
gr.hier_block2.__init__( self,
"ofdm_inner_receiver",
gr.io_signature(
1, 1,
gr.sizeof_gr_complex ),
gr.io_signaturev(
4, 4,
[gr.sizeof_gr_complex * total_subc, # OFDM blocks
gr.sizeof_char, # Frame start
gr.sizeof_float * total_subc,
gr.sizeof_float] ) ) # Normalized |CTF|^2
## Input and output ports
self.input = rx_input = self
out_ofdm_blocks = ( self, 0 )
out_frame_start = ( self, 1 )
out_disp_ctf = ( self, 2 )
out_disp_cfo = ( self, 3 )
## pre-FFT processing
if options.ideal is False and options.ideal2 is False:
if options.old_receiver is False:
## Compute autocorrelations for S&C preamble
## and cyclic prefix
self._sc_metric = sc_metric = autocorrelator( fft_length/2, fft_length/2 )
self._gi_metric = gi_metric = autocorrelator( fft_length, cp_length )
self.connect( rx_input, sc_metric )
self.connect( rx_input, gi_metric )
## Sync. Output contains OFDM blocks
sync = ofdm.time_sync( fft_length, cp_length )
self.connect( rx_input, ( sync, 0 ) )
self.connect( sc_metric, ( sync, 1 ) )
self.connect( gi_metric, ( sync, 2 ) )
ofdm_blocks = ( sync, 0 )
frame_start = ( sync, 1 )
#log_to_file( self, ( sync, 1 ), "data/peak_detector.char" )
else:
#Testing old/new metric
self.tm = schmidl.recursive_timing_metric(fft_length)
self.connect( self.input, self.tm)
#log_to_file( self, self.tm, "data/rec_sc_metric_ofdm.float" )
timingmetric_shift = -2#int(-cp_length/4)# 0#-2 #int(-cp_length * 0.8)
tmfilter = filter.fft_filter_fff(1, [1./cp_length]*cp_length)
self.connect( self.tm, tmfilter )
self.tm = tmfilter
self._pd_thres = 0.3
self._pd_lookahead = fft_length / 2 # empirically chosen
peak_detector = ofdm.peak_detector_02_fb(self._pd_lookahead, self._pd_thres)
self.connect(self.tm, peak_detector)
#log_to_file( self, peak_detector, "data/rec_peak_detector.char" )
frame_start = [0]*frame_length
frame_start[0] = 1
frame_start = self.frame_trigger_old = blocks.vector_source_b(frame_start,True)
delayed_timesync = blocks.delay(gr.sizeof_char,
(frame_length-1)*block_length + timingmetric_shift)
self.connect( peak_detector, delayed_timesync )
self.block_sampler = ofdm.vector_sampler(gr.sizeof_gr_complex,block_length*frame_length)
self.discard_cp = ofdm.vector_mask(block_length,cp_length,fft_length,[])
self.connect(self.input,self.block_sampler)
self.connect(delayed_timesync,(self.block_sampler,1))
# TODO: dynamic solution
vt2s = blocks.vector_to_stream(gr.sizeof_gr_complex*block_length,
frame_length)
self.connect(self.block_sampler,vt2s,self.discard_cp)
#terminate_stream(self,ofdm_blocks)
ofdm_blocks = self.discard_cp
# else:
# serial_to_parallel = blocks.stream_to_vector(gr.sizeof_gr_complex,block_length)
# discard_cp = ofdm.vector_mask(block_length,cp_length,fft_length,[])
# ofdm_blocks = discard_cp
# self.connect( rx_input, serial_to_parallel, discard_cp )
# frame_start = [0]*frame_length
# frame_start[0] = 1
# frame_start = blocks.vector_source_b(frame_start,True)
#
# print "Disabled time synchronization stage"
## Compute autocorrelations for S&C preamble
## and cyclic prefix
#log_to_file( self, sc_metric, "data/sc_metric_ofdm.float" )
#log_to_file(self, frame_start, "data/frame_start.compl")
# log_to_file(self,ofdm_blocks,"data/ofdm_blocks_original.compl")
if options.disable_time_sync or options.ideal or options.ideal2:
if options.ideal is False and options.ideal2 is False:
terminate_stream(self, ofdm_blocks)
terminate_stream(self, frame_start)
serial_to_parallel = blocks.stream_to_vector(gr.sizeof_gr_complex,block_length)
discard_cp = ofdm.vector_mask_dc_null(block_length,cp_length,fft_length,dc_null, [])
ofdm_blocks = discard_cp
self.connect( rx_input, serial_to_parallel, discard_cp )
frame_start = [0]*frame_length
frame_start[0] = 1
frame_start = blocks.vector_source_b(frame_start,True)
print "Disabled time synchronization stage"
print"\t\t\t\t\tframe_length = ",frame_length
if options.ideal is False and options.ideal2 is False:
## Extract preamble, feed to Morelli & Mengali frequency offset estimator
assert( block_header.mm_preamble_pos == 0 )
morelli_foe = ofdm.mm_frequency_estimator( fft_length, L,1,0 )
sampler_preamble = ofdm.vector_sampler( gr.sizeof_gr_complex * fft_length,
1 )
self.connect( ofdm_blocks, ( sampler_preamble, 0 ) )
self.connect( frame_start, ( sampler_preamble, 1 ) )
self.connect( sampler_preamble, morelli_foe )
freq_offset = morelli_foe
## Adaptive LMS FIR filtering of frequency offset
lms_fir = ofdm.lms_fir_ff( 20, 1e-3 ) # TODO: verify parameter choice
self.connect( freq_offset, lms_fir )
freq_offset = lms_fir
#self.zmq_probe_freqoff = zeromq.pub_sink(gr.sizeof_float, 1, "tcp://*:5557")
self.connect(lms_fir, blocks.keep_one_in_n(gr.sizeof_float,20) ,out_disp_cfo)
else:
self.connect(blocks.vector_source_f ([1]) ,out_disp_cfo)
#log_to_file(self, lms_fir, "data/lms_fir.float")
if options.disable_freq_sync or options.ideal or options.ideal2:
if options.ideal is False and options.ideal2 is False:
terminate_stream(self, freq_offset)
freq_offset = blocks.vector_source_f([0.0],True)
print "Disabled frequency synchronization stage"
if options.ideal is False and options.ideal2 is False:
## Correct frequency shift, feed-forward structure
frequency_shift = ofdm.frequency_shift_vcc( fft_length, -1.0/fft_length,
cp_length )
self.connect( ofdm_blocks, ( frequency_shift, 0 ) )
self.connect( freq_offset, ( frequency_shift, 1 ) )
self.connect( frame_start, ( frequency_shift, 2 ) )
ofdm_blocks = frequency_shift
## FFT
fft = fft_blocks.fft_vcc( fft_length, True, [], True )
self.connect( ofdm_blocks, fft )
ofdm_blocks = fft
#log_to_file( self, fft, "data/compen.float" )
## Remove virtual subcarriers
if fft_length > data_subc:
subcarrier_mask = ofdm.vector_mask_dc_null( fft_length, virtual_subc/2,
total_subc, dc_null, [] )
self.connect( ofdm_blocks, subcarrier_mask )
ofdm_blocks = subcarrier_mask
#log_to_file(self, ofdm_blocks, "data/vec_mask.compl")
## Least Squares estimator for channel transfer function (CTF)
if options.logcir:
log_to_file( self, ofdm_blocks, "data/OFDM_Blocks.compl" )
inv_preamble_fd = numpy.array( block_header.pilotsym_fd[
block_header.channel_estimation_pilot[0] ] )
inv_preamble_fd = numpy.concatenate([inv_preamble_fd[:total_subc/2],inv_preamble_fd[total_subc/2+dc_null:]])
#print "Channel estimation pilot: ", inv_preamble_fd
inv_preamble_fd = 1. / inv_preamble_fd
LS_channel_estimator0 = ofdm.multiply_const_vcc( list( inv_preamble_fd ) )
self.connect( ofdm_blocks, LS_channel_estimator0, gr.null_sink(gr.sizeof_gr_complex*total_subc))
log_to_file( self, LS_channel_estimator0, "data/OFDM_Blocks_eq.compl" )
## post-FFT processing
## extract channel estimation preamble from frame
if options.ideal is False and options.ideal2 is False:
chest_pre_trigger = blocks.delay( gr.sizeof_char, 1)
sampled_chest_preamble = ofdm.vector_sampler( gr.sizeof_gr_complex * total_subc, 1)
self.connect( frame_start, chest_pre_trigger )
self.connect( chest_pre_trigger, ( sampled_chest_preamble, 1 ) )
self.connect( ofdm_blocks, ( sampled_chest_preamble, 0 ) )
## Least Squares estimator for channel transfer function (CTF)
inv_preamble_fd = numpy.array( block_header.pilotsym_fd[
block_header.channel_estimation_pilot[0] ] )
inv_preamble_fd = numpy.concatenate([inv_preamble_fd[:total_subc/2],inv_preamble_fd[total_subc/2+dc_null:]])
#print "Channel estimation pilot: ", inv_preamble_fd
inv_preamble_fd = 1. / inv_preamble_fd
LS_channel_estimator = ofdm.multiply_const_vcc( list( inv_preamble_fd ) )
self.connect( sampled_chest_preamble, LS_channel_estimator )
estimated_CTF = LS_channel_estimator
if options.logcir:
log_to_file( self, sampled_chest_preamble, "data/PREAM.compl" )
if not options.disable_ctf_enhancer:
if options.logcir:
ifft1 = fft_blocks.fft_vcc(total_subc,False,[],True)
self.connect( estimated_CTF, ifft1,gr.null_sink(gr.sizeof_gr_complex*total_subc))
summ1 = ofdm.vector_sum_vcc(total_subc)
c2m =gr.complex_to_mag(total_subc)
self.connect( estimated_CTF,summ1 ,gr.null_sink(gr.sizeof_gr_complex))
self.connect( estimated_CTF, c2m,gr.null_sink(gr.sizeof_float*total_subc))
log_to_file( self, ifft1, "data/CIR1.compl" )
log_to_file( self, summ1, "data/CTFsumm1.compl" )
log_to_file( self, estimated_CTF, "data/CTF1.compl" )
log_to_file( self, c2m, "data/CTFmag1.float" )
## MSE enhancer
ctf_mse_enhancer = ofdm.CTF_MSE_enhancer( total_subc, cp_length + cp_length)
self.connect( estimated_CTF, ctf_mse_enhancer )
# log_to_file( self, ctf_mse_enhancer, "data/ctf_mse_enhancer_original.compl")
#ifft3 = fft_blocks.fft_vcc(total_subc,False,[],True)
#null_noise = ofdm.noise_nulling(total_subc, cp_length + cp_length)
#ctf_mse_enhancer = fft_blocks.fft_vcc(total_subc,True,[],True)
#ctf_mse_enhancer = ofdm.vector_mask( fft_length, virtual_subc/2,
# total_subc, [] )
#self.connect( estimated_CTF, ifft3,null_noise,ctf_mse_enhancer )
estimated_CTF = ctf_mse_enhancer
print "Disabled CTF MSE enhancer"
if options.logcir:
ifft2 = fft_blocks.fft_vcc(total_subc,False,[],True)
self.connect( estimated_CTF, ifft2,gr.null_sink(gr.sizeof_gr_complex*total_subc))
summ2 = ofdm.vector_sum_vcc(total_subc)
c2m2 =gr.complex_to_mag(total_subc)
self.connect( estimated_CTF,summ2 ,gr.null_sink(gr.sizeof_gr_complex))
self.connect( estimated_CTF, c2m2,gr.null_sink(gr.sizeof_float*total_subc))
log_to_file( self, ifft2, "data/CIR2.compl" )
log_to_file( self, summ2, "data/CTFsumm2.compl" )
log_to_file( self, estimated_CTF, "data/CTF2.compl" )
log_to_file( self, c2m2, "data/CTFmag2.float" )
## Postprocess the CTF estimate
## CTF -> inverse CTF (for equalizer)
## CTF -> norm |.|^2 (for CTF display)
ctf_postprocess = ofdm.postprocess_CTF_estimate( total_subc )
self.connect( estimated_CTF, ctf_postprocess )
inv_estimated_CTF = ( ctf_postprocess, 0 )
disp_CTF = ( ctf_postprocess, 1 )
# if options.disable_equalization or options.ideal:
# terminate_stream(self, inv_estimated_CTF)
# inv_estimated_CTF_vec = blocks.vector_source_c([1.0/fft_length*math.sqrt(total_subc)]*total_subc,True,total_subc)
# inv_estimated_CTF_str = blocks.vector_to_stream(gr.sizeof_gr_complex, total_subc)
# self.inv_estimated_CTF_mul = ofdm.multiply_const_ccf( 1.0/config.rms_amplitude )
# #inv_estimated_CTF_mul.set_k(1.0/config.rms_amplitude)
# inv_estimated_CTF = blocks.stream_to_vector(gr.sizeof_gr_complex, total_subc)
# self.connect( inv_estimated_CTF_vec, inv_estimated_CTF_str, self.inv_estimated_CTF_mul, inv_estimated_CTF)
# print "Disabled equalization stage"
'''
## LMS Phase tracking
## Track residual frequency offset and sampling clock frequency offset
nondata_blocks = []
for i in range(config.frame_length):
if i in config.training_data.pilotsym_pos:
nondata_blocks.append(i)
print"\t\t\t\t\tnondata_blocks=",nondata_blocks
pilot_subc = block_header.pilot_tones
pilot_subcarriers = block_header.pilot_subc_sym
print "PILOT SUBCARRIERS: ", pilot_subcarriers
phase_tracking = ofdm.lms_phase_tracking_03( total_subc, pilot_subc,
nondata_blocks, pilot_subcarriers,0 )
self.connect( ofdm_blocks, ( phase_tracking, 0 ) )
self.connect( inv_estimated_CTF, ( phase_tracking, 1 ) )
self.connect( frame_start, ( phase_tracking, 2 ) ) ##
if options.scatter_plot_before_phase_tracking:
self.before_phase_tracking = equalizer
if options.disable_phase_tracking or options.ideal:
terminate_stream(self, phase_tracking)
print "Disabled phase tracking stage"
else:
ofdm_blocks = phase_tracking
'''
## Channel Equalizer
if options.disable_equalization or options.ideal or options.ideal2:
print "Disabled equalization stage"
if options.ideal is False and options.ideal2 is False:
terminate_stream(self, inv_estimated_CTF)
else:
equalizer = ofdm.channel_equalizer( total_subc )
self.connect( ofdm_blocks, ( equalizer, 0 ) )
self.connect( inv_estimated_CTF, ( equalizer, 1 ) )
self.connect( frame_start, ( equalizer, 2 ) )
ofdm_blocks = equalizer
#log_to_file(self, equalizer,"data/equalizer_siso.compl")
#log_to_file(self, ofdm_blocks, "data/equalizer.compl")
## LMS Phase tracking
## Track residual frequency offset and sampling clock frequency offset
if options.ideal is False and options.ideal2 is False:
nondata_blocks = []
for i in range(config.frame_length):
if i in config.training_data.pilotsym_pos:
nondata_blocks.append(i)
print"\t\t\t\t\tnondata_blocks=",nondata_blocks
pilot_subc = block_header.pilot_tones
pilot_subcarriers = block_header.pilot_subc_sym
print "PILOT SUBCARRIERS: ", pilot_subcarriers
phase_tracking2 = ofdm.lms_phase_tracking_dc_null( total_subc, pilot_subc,
nondata_blocks, pilot_subcarriers, dc_null )
self.connect( ofdm_blocks, ( phase_tracking2, 0 ) )
self.connect( frame_start, ( phase_tracking2, 1 ) ) ##
if options.disable_phase_tracking or options.ideal or options.ideal2:
if options.ideal is False and options.ideal2 is False:
terminate_stream(self, phase_tracking2)
print "Disabled phase tracking stage"
else:
ofdm_blocks = phase_tracking2
if options.scatter_plot_before_phase_tracking:
self.before_phase_tracking = equalizer
## Output connections
self.connect( ofdm_blocks, out_ofdm_blocks )
self.connect( frame_start, out_frame_start )
if options.ideal is False and options.ideal2 is False:
self.connect( disp_CTF, out_disp_ctf )
else:
self.connect( blocks.vector_source_f([1.0]*total_subc),blocks.stream_to_vector(gr.sizeof_float,total_subc), out_disp_ctf )
if log:
log_to_file( self, sc_metric, "data/sc_metric.float" )
log_to_file( self, gi_metric, "data/gi_metric.float" )
log_to_file( self, morelli_foe, "data/morelli_foe.float" )
log_to_file( self, lms_fir, "data/lms_fir.float" )
log_to_file( self, sampler_preamble, "data/preamble.compl" )
log_to_file( self, sync, "data/sync.compl" )
log_to_file( self, frequency_shift, "data/frequency_shift.compl" )
log_to_file( self, fft, "data/fft.compl")
log_to_file( self, fft, "data/fft.float", mag=True )
if vars().has_key( 'subcarrier_mask' ):
log_to_file( self, subcarrier_mask, "data/subcarrier_mask.compl" )
log_to_file( self, ofdm_blocks, "data/ofdm_blocks_out.compl" )
log_to_file( self, frame_start, "data/frame_start.float",
char_to_float=True )
log_to_file( self, sampled_chest_preamble,
"data/sampled_chest_preamble.compl" )
log_to_file( self, LS_channel_estimator,
"data/ls_channel_estimator.compl" )
log_to_file( self, LS_channel_estimator,
"data/ls_channel_estimator.float", mag=True )
if "ctf_mse_enhancer" in locals():
log_to_file( self, ctf_mse_enhancer, "data/ctf_mse_enhancer.compl" )
log_to_file( self, ctf_mse_enhancer, "data/ctf_mse_enhancer.float",
mag=True )
log_to_file( self, (ctf_postprocess,0), "data/inc_estimated_ctf.compl" )
log_to_file( self, (ctf_postprocess,1), "data/disp_ctf.float" )
log_to_file( self, equalizer, "data/equalizer.compl" )
log_to_file( self, equalizer, "data/equalizer.float", mag=True )
log_to_file( self, phase_tracking, "data/phase_tracking.compl" )
def __init__(self):
gr.top_block.__init__(self, "Simulation for SDR TCC")
Qt.QWidget.__init__(self)
self.setWindowTitle("Simulation for SDR TCC")
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", "sdrSim2")
self.restoreGeometry(self.settings.value("geometry").toByteArray())
##################################################
# Variables
##################################################
self.sps = sps = 4
self.nfilts = nfilts = 32
self.eb = eb = 0.35
self.timing_loop_bw = timing_loop_bw = 0.005
self.time_offset = time_offset = 1.00
self.taps = taps = [1.0, 0.25-0.25j, 0.50 + 0.10j, -0.3 + 0.2j]
self.samp_rate = samp_rate = 32000
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.phase_bw = phase_bw = 6.28/100.0
self.payload_size = payload_size = 992
self.noise_volt = noise_volt = 0.0001
self.matched_filter = matched_filter = firdes.root_raised_cosine(nfilts, nfilts, 1, eb, int(11*sps*nfilts))
self.freq_offset = freq_offset = 0
self.freq_bw = freq_bw = 6.28/100.0
self.excess_bw = excess_bw = 0.35
self.eq_gain = eq_gain = 0.01
self.delay = delay = 0
self.arity = arity = 4
##################################################
# Blocks
##################################################
self.controls = Qt.QTabWidget()
self.controls_widget_0 = Qt.QWidget()
self.controls_layout_0 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.controls_widget_0)
self.controls_grid_layout_0 = Qt.QGridLayout()
self.controls_layout_0.addLayout(self.controls_grid_layout_0)
self.controls.addTab(self.controls_widget_0, "Channel")
self.controls_widget_1 = Qt.QWidget()
self.controls_layout_1 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.controls_widget_1)
self.controls_grid_layout_1 = Qt.QGridLayout()
self.controls_layout_1.addLayout(self.controls_grid_layout_1)
self.controls.addTab(self.controls_widget_1, "Receiver")
self.top_grid_layout.addWidget(self.controls, 0,0,1,2)
self._timing_loop_bw_layout = Qt.QVBoxLayout()
self._timing_loop_bw_label = Qt.QLabel("Time: BW")
self._timing_loop_bw_slider = Qwt.QwtSlider(None, Qt.Qt.Horizontal, Qwt.QwtSlider.BottomScale, Qwt.QwtSlider.BgSlot)
self._timing_loop_bw_slider.setRange(0.0, 0.2, 0.005)
self._timing_loop_bw_slider.setValue(self.timing_loop_bw)
self._timing_loop_bw_slider.setMinimumWidth(200)
self._timing_loop_bw_slider.valueChanged.connect(self.set_timing_loop_bw)
self._timing_loop_bw_label.setAlignment(Qt.Qt.AlignBottom | Qt.Qt.AlignHCenter)
self._timing_loop_bw_layout.addWidget(self._timing_loop_bw_label)
self._timing_loop_bw_layout.addWidget(self._timing_loop_bw_slider)
self.controls_grid_layout_1.addLayout(self._timing_loop_bw_layout, 0,0,1,1)
self._time_offset_layout = Qt.QVBoxLayout()
self._time_offset_tool_bar = Qt.QToolBar(self)
self._time_offset_layout.addWidget(self._time_offset_tool_bar)
self._time_offset_tool_bar.addWidget(Qt.QLabel("Timing Offset"+": "))
class qwt_counter_pyslot(Qwt.QwtCounter):
def __init__(self, parent=None):
Qwt.QwtCounter.__init__(self, parent)
@pyqtSlot('double')
def setValue(self, value):
super(Qwt.QwtCounter, self).setValue(value)
self._time_offset_counter = qwt_counter_pyslot()
self._time_offset_counter.setRange(0.999, 1.001, 0.0001)
self._time_offset_counter.setNumButtons(2)
self._time_offset_counter.setValue(self.time_offset)
self._time_offset_tool_bar.addWidget(self._time_offset_counter)
self._time_offset_counter.valueChanged.connect(self.set_time_offset)
self._time_offset_slider = Qwt.QwtSlider(None, Qt.Qt.Horizontal, Qwt.QwtSlider.BottomScale, Qwt.QwtSlider.BgSlot)
self._time_offset_slider.setRange(0.999, 1.001, 0.0001)
self._time_offset_slider.setValue(self.time_offset)
self._time_offset_slider.setMinimumWidth(200)
self._time_offset_slider.valueChanged.connect(self.set_time_offset)
self._time_offset_layout.addWidget(self._time_offset_slider)
self.controls_grid_layout_0.addLayout(self._time_offset_layout, 0,2,1,1)
self.received = Qt.QTabWidget()
self.received_widget_0 = Qt.QWidget()
self.received_layout_0 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.received_widget_0)
self.received_grid_layout_0 = Qt.QGridLayout()
self.received_layout_0.addLayout(self.received_grid_layout_0)
self.received.addTab(self.received_widget_0, "Constellation")
self.received_widget_1 = Qt.QWidget()
self.received_layout_1 = Qt.QBoxLayout(Qt.QBoxLayout.TopToBottom, self.received_widget_1)
self.received_grid_layout_1 = Qt.QGridLayout()
self.received_layout_1.addLayout(self.received_grid_layout_1)
self.received.addTab(self.received_widget_1, "Symbols")
self.top_grid_layout.addWidget(self.received, 2,0,1,1)
self._phase_bw_layout = Qt.QVBoxLayout()
self._phase_bw_label = Qt.QLabel("Phase: Bandwidth")
self._phase_bw_slider = Qwt.QwtSlider(None, Qt.Qt.Horizontal, Qwt.QwtSlider.BottomScale, Qwt.QwtSlider.BgSlot)
self._phase_bw_slider.setRange(0.0, 1.0, 0.01)
self._phase_bw_slider.setValue(self.phase_bw)
self._phase_bw_slider.setMinimumWidth(200)
self._phase_bw_slider.valueChanged.connect(self.set_phase_bw)
self._phase_bw_label.setAlignment(Qt.Qt.AlignBottom | Qt.Qt.AlignHCenter)
self._phase_bw_layout.addWidget(self._phase_bw_label)
self._phase_bw_layout.addWidget(self._phase_bw_slider)
self.controls_grid_layout_1.addLayout(self._phase_bw_layout, 0,2,1,1)
self._noise_volt_layout = Qt.QVBoxLayout()
self._noise_volt_tool_bar = Qt.QToolBar(self)
self._noise_volt_layout.addWidget(self._noise_volt_tool_bar)
self._noise_volt_tool_bar.addWidget(Qt.QLabel("Noise Voltage"+": "))
class qwt_counter_pyslot(Qwt.QwtCounter):
def __init__(self, parent=None):
Qwt.QwtCounter.__init__(self, parent)
@pyqtSlot('double')
def setValue(self, value):
super(Qwt.QwtCounter, self).setValue(value)
self._noise_volt_counter = qwt_counter_pyslot()
self._noise_volt_counter.setRange(0, 1, 0.01)
self._noise_volt_counter.setNumButtons(2)
self._noise_volt_counter.setValue(self.noise_volt)
self._noise_volt_tool_bar.addWidget(self._noise_volt_counter)
self._noise_volt_counter.valueChanged.connect(self.set_noise_volt)
self._noise_volt_slider = Qwt.QwtSlider(None, Qt.Qt.Horizontal, Qwt.QwtSlider.BottomScale, Qwt.QwtSlider.BgSlot)
self._noise_volt_slider.setRange(0, 1, 0.01)
self._noise_volt_slider.setValue(self.noise_volt)
self._noise_volt_slider.setMinimumWidth(200)
self._noise_volt_slider.valueChanged.connect(self.set_noise_volt)
self._noise_volt_layout.addWidget(self._noise_volt_slider)
self.controls_grid_layout_0.addLayout(self._noise_volt_layout, 0,0,1,1)
self._freq_offset_layout = Qt.QVBoxLayout()
self._freq_offset_tool_bar = Qt.QToolBar(self)
self._freq_offset_layout.addWidget(self._freq_offset_tool_bar)
self._freq_offset_tool_bar.addWidget(Qt.QLabel("Frequency Offset"+": "))
class qwt_counter_pyslot(Qwt.QwtCounter):
def __init__(self, parent=None):
Qwt.QwtCounter.__init__(self, parent)
@pyqtSlot('double')
def setValue(self, value):
super(Qwt.QwtCounter, self).setValue(value)
self._freq_offset_counter = qwt_counter_pyslot()
self._freq_offset_counter.setRange(-0.1, 0.1, 0.001)
self._freq_offset_counter.setNumButtons(2)
self._freq_offset_counter.setValue(self.freq_offset)
self._freq_offset_tool_bar.addWidget(self._freq_offset_counter)
self._freq_offset_counter.valueChanged.connect(self.set_freq_offset)
self._freq_offset_slider = Qwt.QwtSlider(None, Qt.Qt.Horizontal, Qwt.QwtSlider.BottomScale, Qwt.QwtSlider.BgSlot)
self._freq_offset_slider.setRange(-0.1, 0.1, 0.001)
self._freq_offset_slider.setValue(self.freq_offset)
self._freq_offset_slider.setMinimumWidth(200)
self._freq_offset_slider.valueChanged.connect(self.set_freq_offset)
self._freq_offset_layout.addWidget(self._freq_offset_slider)
self.controls_grid_layout_0.addLayout(self._freq_offset_layout, 0,1,1,1)
self._eq_gain_layout = Qt.QVBoxLayout()
self._eq_gain_label = Qt.QLabel("Equalizer: rate")
self._eq_gain_slider = Qwt.QwtSlider(None, Qt.Qt.Horizontal, Qwt.QwtSlider.BottomScale, Qwt.QwtSlider.BgSlot)
self._eq_gain_slider.setRange(0.0, 0.1, 0.001)
self._eq_gain_slider.setValue(self.eq_gain)
self._eq_gain_slider.setMinimumWidth(200)
self._eq_gain_slider.valueChanged.connect(self.set_eq_gain)
self._eq_gain_label.setAlignment(Qt.Qt.AlignBottom | Qt.Qt.AlignHCenter)
self._eq_gain_layout.addWidget(self._eq_gain_label)
self._eq_gain_layout.addWidget(self._eq_gain_slider)
self.controls_grid_layout_1.addLayout(self._eq_gain_layout, 0,1,1,1)
self._delay_layout = Qt.QVBoxLayout()
self._delay_tool_bar = Qt.QToolBar(self)
self._delay_layout.addWidget(self._delay_tool_bar)
self._delay_tool_bar.addWidget(Qt.QLabel("Delay"+": "))
class qwt_counter_pyslot(Qwt.QwtCounter):
def __init__(self, parent=None):
Qwt.QwtCounter.__init__(self, parent)
@pyqtSlot('double')
def setValue(self, value):
super(Qwt.QwtCounter, self).setValue(value)
self._delay_counter = qwt_counter_pyslot()
self._delay_counter.setRange(0, 200, 1)
self._delay_counter.setNumButtons(2)
self._delay_counter.setValue(self.delay)
self._delay_tool_bar.addWidget(self._delay_counter)
self._delay_counter.valueChanged.connect(self.set_delay)
self._delay_slider = Qwt.QwtSlider(None, Qt.Qt.Horizontal, Qwt.QwtSlider.BottomScale, Qwt.QwtSlider.BgSlot)
self._delay_slider.setRange(0, 200, 1)
self._delay_slider.setValue(self.delay)
self._delay_slider.setMinimumWidth(200)
self._delay_slider.valueChanged.connect(self.set_delay)
self._delay_layout.addWidget(self._delay_slider)
self.top_grid_layout.addLayout(self._delay_layout, 1,0,1,1)
self.tutorial_my_qpsk_demod_cb_1 = tutorial.my_qpsk_demod_cb(True)
self.qtgui_time_sink_x_0_0 = qtgui.time_sink_f(
500, #size
samp_rate, #samp_rate
"", #name
2 #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, 2)
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(False)
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):
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.top_grid_layout.addWidget(self._qtgui_time_sink_x_0_0_win, 2,1,1,1)
self.qtgui_time_sink_x_0 = qtgui.time_sink_f(
500, #size
samp_rate, #samp_rate
"", #name
1 #number of inputs
)
self.qtgui_time_sink_x_0.set_update_time(0.10)
self.qtgui_time_sink_x_0.set_y_axis(-1, 4)
self.qtgui_time_sink_x_0.enable_tags(-1, True)
self.qtgui_time_sink_x_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, 0, "")
self.qtgui_time_sink_x_0.enable_autoscale(False)
labels = ["", "", "", "", "",
"", "", "", "", ""]
widths = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
colors = ["blue", "red", "green", "black", "cyan",
"magenta", "yellow", "dark red", "dark green", "blue"]
styles = [1, 1, 1, 1, 1,
1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1,
-1, -1, -1, -1, -1]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_time_sink_x_0.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_time_sink_x_0.set_line_label(i, labels[i])
self.qtgui_time_sink_x_0.set_line_width(i, widths[i])
self.qtgui_time_sink_x_0.set_line_color(i, colors[i])
self.qtgui_time_sink_x_0.set_line_style(i, styles[i])
self.qtgui_time_sink_x_0.set_line_marker(i, markers[i])
self.qtgui_time_sink_x_0.set_line_alpha(i, alphas[i])
self._qtgui_time_sink_x_0_win = sip.wrapinstance(self.qtgui_time_sink_x_0.pyqwidget(), Qt.QWidget)
self.received_grid_layout_1.addWidget(self._qtgui_time_sink_x_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.enable_autoscale(False)
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.received_grid_layout_0.addWidget(self._qtgui_const_sink_x_0_win, 0,0,1,1)
self._freq_bw_layout = Qt.QVBoxLayout()
self._freq_bw_label = Qt.QLabel("Frequency Bandwidth")
self._freq_bw_slider = Qwt.QwtSlider(None, Qt.Qt.Horizontal, Qwt.QwtSlider.BottomScale, Qwt.QwtSlider.BgSlot)
self._freq_bw_slider.setRange(0.0, 1.0, 0.01)
self._freq_bw_slider.setValue(self.freq_bw)
self._freq_bw_slider.setMinimumWidth(200)
self._freq_bw_slider.valueChanged.connect(self.set_freq_bw)
self._freq_bw_label.setAlignment(Qt.Qt.AlignBottom | Qt.Qt.AlignHCenter)
self._freq_bw_layout.addWidget(self._freq_bw_label)
self._freq_bw_layout.addWidget(self._freq_bw_slider)
self.controls_grid_layout_1.addLayout(self._freq_bw_layout, 0,3,1,1)
self.digital_pfb_clock_sync_xxx_0 = digital.pfb_clock_sync_ccf(sps, timing_loop_bw, (rrc_taps), nfilts, nfilts/2, 1.5, 2)
self.digital_map_bb_0 = digital.map_bb(([0,1,3,2]))
self.digital_diff_decoder_bb_0 = digital.diff_decoder_bb(4)
self.digital_costas_loop_cc_0 = digital.costas_loop_cc(phase_bw, arity)
self.digital_correlate_and_sync_cc_0 = digital.correlate_and_sync_cc((preamble), (matched_filter), sps)
self.digital_constellation_modulator_0 = digital.generic_mod(
constellation=qpsk,
differential=True,
samples_per_symbol=sps,
pre_diff_code=True,
excess_bw=excess_bw,
verbose=False,
log=False,
)
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_volt,
frequency_offset=freq_offset,
epsilon=time_offset,
taps=(taps),
noise_seed=0,
block_tags=False
)
self.blocks_vector_source_x_0_0 = blocks.vector_source_b(map(lambda x: (-x+1)/2, preamble), True, 1, [])
self.blocks_unpack_k_bits_bb_0_0 = blocks.unpack_k_bits_bb(8)
self.blocks_unpack_k_bits_bb_0 = blocks.unpack_k_bits_bb(2)
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_gr_complex*1, samp_rate,True)
self.blocks_stream_mux_0_0_0 = blocks.stream_mux(gr.sizeof_char*1, (len(preamble)/8,payload_size))
self.blocks_pack_k_bits_bb_1 = blocks.pack_k_bits_bb(8)
self.blocks_pack_k_bits_bb_0 = blocks.pack_k_bits_bb(8)
self.blocks_null_sink_0 = blocks.null_sink(gr.sizeof_gr_complex*1)
self.blocks_file_source_0 = blocks.file_source(gr.sizeof_char*1, "/home/franchz/tcc/gnuradio/tx_teste.txt", True)
self.blocks_file_sink_0 = blocks.file_sink(gr.sizeof_char*1, "/home/franchz/tcc/gnuradio/rx_teste.txt", False)
self.blocks_file_sink_0.set_unbuffered(False)
self.blocks_delay_0 = blocks.delay(gr.sizeof_float*1, int(delay))
self.blocks_char_to_float_0_0_0 = blocks.char_to_float(1, 1)
self.blocks_char_to_float_0_0 = blocks.char_to_float(1, 1)
self.blocks_char_to_float_0 = blocks.char_to_float(1, 1)
##################################################
# Connections
##################################################
self.connect((self.blocks_throttle_0, 0), (self.channels_channel_model_0, 0))
self.connect((self.digital_costas_loop_cc_0, 0), (self.qtgui_const_sink_x_0, 0))
self.connect((self.blocks_char_to_float_0, 0), (self.qtgui_time_sink_x_0, 0))
self.connect((self.digital_cma_equalizer_cc_0, 0), (self.digital_costas_loop_cc_0, 0))
self.connect((self.digital_pfb_clock_sync_xxx_0, 0), (self.digital_cma_equalizer_cc_0, 0))
self.connect((self.blocks_unpack_k_bits_bb_0_0, 0), (self.blocks_char_to_float_0_0_0, 0))
self.connect((self.blocks_unpack_k_bits_bb_0, 0), (self.blocks_char_to_float_0_0, 0))
self.connect((self.blocks_vector_source_x_0_0, 0), (self.blocks_pack_k_bits_bb_0, 0))
self.connect((self.blocks_pack_k_bits_bb_0, 0), (self.blocks_stream_mux_0_0_0, 0))
self.connect((self.channels_channel_model_0, 0), (self.digital_correlate_and_sync_cc_0, 0))
self.connect((self.digital_correlate_and_sync_cc_0, 0), (self.digital_pfb_clock_sync_xxx_0, 0))
self.connect((self.digital_correlate_and_sync_cc_0, 1), (self.blocks_null_sink_0, 0))
self.connect((self.blocks_file_source_0, 0), (self.blocks_stream_mux_0_0_0, 1))
self.connect((self.blocks_file_source_0, 0), (self.blocks_unpack_k_bits_bb_0_0, 0))
self.connect((self.blocks_char_to_float_0_0, 0), (self.qtgui_time_sink_x_0_0, 1))
self.connect((self.blocks_char_to_float_0_0_0, 0), (self.blocks_delay_0, 0))
self.connect((self.blocks_delay_0, 0), (self.qtgui_time_sink_x_0_0, 0))
self.connect((self.blocks_pack_k_bits_bb_1, 0), (self.blocks_file_sink_0, 0))
self.connect((self.blocks_unpack_k_bits_bb_0, 0), (self.blocks_pack_k_bits_bb_1, 0))
self.connect((self.blocks_stream_mux_0_0_0, 0), (self.digital_constellation_modulator_0, 0))
self.connect((self.digital_constellation_modulator_0, 0), (self.blocks_throttle_0, 0))
self.connect((self.tutorial_my_qpsk_demod_cb_1, 0), (self.digital_map_bb_0, 0))
self.connect((self.digital_map_bb_0, 0), (self.digital_diff_decoder_bb_0, 0))
self.connect((self.digital_costas_loop_cc_0, 0), (self.tutorial_my_qpsk_demod_cb_1, 0))
self.connect((self.digital_diff_decoder_bb_0, 0), (self.blocks_char_to_float_0, 0))
self.connect((self.digital_diff_decoder_bb_0, 0), (self.blocks_unpack_k_bits_bb_0, 0))
def __init__(self, fft_length, cp_length, snr, kstime, logging):
''' Maximum Likelihood OFDM synchronizer:
J. van de Beek, M. Sandell, and P. O. Borjesson, "ML Estimation
of Time and Frequency Offset in OFDM Systems," IEEE Trans.
Signal Processing, vol. 45, no. 7, pp. 1800-1805, 1997.
'''
gr.hier_block2.__init__(self, "ofdm_sync_ml",
gr.io_signature(1, 1, gr.sizeof_gr_complex), # Input signature
gr.io_signature2(2, 2, gr.sizeof_float, gr.sizeof_char)) # Output signature
self.input = blocks.add_const_cc(0)
SNR = 10.0**(snr / 10.0)
rho = SNR / (SNR + 1.0)
symbol_length = fft_length + cp_length
# ML Sync
# Energy Detection from ML Sync
self.connect(self, self.input)
# Create a delay line
self.delay = blocks.delay(gr.sizeof_gr_complex, fft_length)
self.connect(self.input, self.delay)
# magnitude squared blocks
self.magsqrd1 = blocks.complex_to_mag_squared()
self.magsqrd2 = blocks.complex_to_mag_squared()
self.adder = blocks.add_ff()
moving_sum_taps = [rho / 2 for i in range(cp_length)]
self.moving_sum_filter = filter.fir_filter_fff(1,moving_sum_taps)
self.connect(self.input,self.magsqrd1)
self.connect(self.delay,self.magsqrd2)
self.connect(self.magsqrd1,(self.adder,0))
self.connect(self.magsqrd2,(self.adder,1))
self.connect(self.adder,self.moving_sum_filter)
# Correlation from ML Sync
self.conjg = blocks.conjugate_cc();
self.mixer = blocks.multiply_cc();
movingsum2_taps = [1.0 for i in range(cp_length)]
self.movingsum2 = filter.fir_filter_ccf(1,movingsum2_taps)
# Correlator data handler
self.c2mag = blocks.complex_to_mag()
self.angle = blocks.complex_to_arg()
self.connect(self.input,(self.mixer,1))
self.connect(self.delay,self.conjg,(self.mixer,0))
self.connect(self.mixer,self.movingsum2,self.c2mag)
self.connect(self.movingsum2,self.angle)
# ML Sync output arg, need to find maximum point of this
self.diff = blocks.sub_ff()
self.connect(self.c2mag,(self.diff,0))
self.connect(self.moving_sum_filter,(self.diff,1))
#ML measurements input to sampler block and detect
self.f2c = blocks.float_to_complex()
self.pk_detect = blocks.peak_detector_fb(0.2, 0.25, 30, 0.0005)
self.sample_and_hold = blocks.sample_and_hold_ff()
# use the sync loop values to set the sampler and the NCO
# self.diff = theta
# self.angle = epsilon
self.connect(self.diff, self.pk_detect)
# The DPLL corrects for timing differences between CP correlations
use_dpll = 0
if use_dpll:
self.dpll = gr.dpll_bb(float(symbol_length),0.01)
self.connect(self.pk_detect, self.dpll)
self.connect(self.dpll, (self.sample_and_hold,1))
else:
self.connect(self.pk_detect, (self.sample_and_hold,1))
self.connect(self.angle, (self.sample_and_hold,0))
################################
# correlate against known symbol
# This gives us the same timing signal as the PN sync block only on the preamble
# we don't use the signal generated from the CP correlation because we don't want
# to readjust the timing in the middle of the packet or we ruin the equalizer settings.
kstime = [k.conjugate() for k in kstime]
kstime.reverse()
self.kscorr = filter.fir_filter_ccc(1, kstime)
self.corrmag = blocks.complex_to_mag_squared()
self.div = blocks.divide_ff()
# The output signature of the correlation has a few spikes because the rest of the
# system uses the repeated preamble symbol. It needs to work that generically if
# anyone wants to use this against a WiMAX-like signal since it, too, repeats.
# The output theta of the correlator above is multiplied with this correlation to
# identify the proper peak and remove other products in this cross-correlation
self.threshold_factor = 0.1
self.slice = blocks.threshold_ff(self.threshold_factor, self.threshold_factor, 0)
self.f2b = blocks.float_to_char()
self.b2f = blocks.char_to_float()
self.mul = blocks.multiply_ff()
# Normalize the power of the corr output by the energy. This is not really needed
# and could be removed for performance, but it makes for a cleaner signal.
# if this is removed, the threshold value needs adjustment.
self.connect(self.input, self.kscorr, self.corrmag, (self.div,0))
self.connect(self.moving_sum_filter, (self.div,1))
self.connect(self.div, (self.mul,0))
self.connect(self.pk_detect, self.b2f, (self.mul,1))
self.connect(self.mul, self.slice)
# Set output signals
# Output 0: fine frequency correction value
# Output 1: timing signal
self.connect(self.sample_and_hold, (self,0))
self.connect(self.slice, self.f2b, (self,1))
if logging:
self.connect(self.moving_sum_filter, blocks.file_sink(gr.sizeof_float, "ofdm_sync_ml-energy_f.dat"))
self.connect(self.diff, blocks.file_sink(gr.sizeof_float, "ofdm_sync_ml-theta_f.dat"))
self.connect(self.angle, blocks.file_sink(gr.sizeof_float, "ofdm_sync_ml-epsilon_f.dat"))
self.connect(self.corrmag, blocks.file_sink(gr.sizeof_float, "ofdm_sync_ml-corrmag_f.dat"))
self.connect(self.kscorr, blocks.file_sink(gr.sizeof_gr_complex, "ofdm_sync_ml-kscorr_c.dat"))
self.connect(self.div, blocks.file_sink(gr.sizeof_float, "ofdm_sync_ml-div_f.dat"))
self.connect(self.mul, blocks.file_sink(gr.sizeof_float, "ofdm_sync_ml-mul_f.dat"))
self.connect(self.slice, blocks.file_sink(gr.sizeof_float, "ofdm_sync_ml-slice_f.dat"))
self.connect(self.pk_detect, blocks.file_sink(gr.sizeof_char, "ofdm_sync_ml-peaks_b.dat"))
if use_dpll:
self.connect(self.dpll, blocks.file_sink(gr.sizeof_char, "ofdm_sync_ml-dpll_b.dat"))
self.connect(self.sample_and_hold, blocks.file_sink(gr.sizeof_float, "ofdm_sync_ml-sample_and_hold_f.dat"))
self.connect(self.input, blocks.file_sink(gr.sizeof_gr_complex, "ofdm_sync_ml-input_c.dat"))
