def __init__(self,
if_rate=None,
filter_type=None,
excess_bw=_def_excess_bw,
symbol_rate=_def_symbol_rate):
"""
Hierarchical block for P25 demodulation base class
@param if_rate: sample rate of complex input channel
@type if_rate: int
"""
self.if_rate = if_rate
self.symbol_rate = symbol_rate
self.bb_sink = None
self.null_sink = blocks.null_sink(gr.sizeof_float)
self.baseband_amp = blocks.multiply_const_ff(_def_bb_gain)
coeffs = op25_c4fm_mod.c4fm_taps(
sample_rate=self.if_rate,
span=9,
generator=op25_c4fm_mod.transfer_function_rx).generate()
sps = self.if_rate / 4800
if filter_type == 'rrc':
ntaps = 7 * sps
if ntaps & 1 == 0:
ntaps += 1
coeffs = filter.firdes.root_raised_cosine(1.0, if_rate,
symbol_rate, excess_bw,
ntaps)
if filter_type == 'gmsk':
# lifted from gmsk.py
_omega = sps
_gain_mu = _def_gmsk_mu
_mu = _def_mu
if not _gain_mu:
_gain_mu = 0.175
_gain_omega = .25 * _gain_mu * _gain_mu # critically damped
self.symbol_filter = blocks.multiply_const_ff(1.0)
self.fsk4_demod = digital.clock_recovery_mm_ff(
_omega, _gain_omega, _mu, _gain_mu, _def_omega_relative_limit)
self.slicer = digital.binary_slicer_fb()
elif filter_type == 'fsk4mm':
self.symbol_filter = filter.fir_filter_fff(1, coeffs)
_omega = sps
_gain_mu = _def_gmsk_mu
_mu = _def_mu
if not _gain_mu:
_gain_mu = 0.0175
_gain_omega = .25 * _gain_mu * _gain_mu # critically damped
self.fsk4_demod = digital.clock_recovery_mm_ff(
_omega, _gain_omega, _mu, _gain_mu, _def_omega_relative_limit)
levels = [-2.0, 0.0, 2.0, 4.0]
self.slicer = op25_repeater.fsk4_slicer_fb(levels)
else:
self.symbol_filter = filter.fir_filter_fff(1, coeffs)
autotuneq = gr.msg_queue(2)
self.fsk4_demod = op25.fsk4_demod_ff(autotuneq, self.if_rate,
self.symbol_rate)
levels = [-2.0, 0.0, 2.0, 4.0]
self.slicer = op25_repeater.fsk4_slicer_fb(levels)
def __init__(self, baudrate=1200, host_address='127.0.0.1', host_port=7355, output_file="/tmp/rx_data.bin", samp_rate=48000):
gr.top_block.__init__(self, "Udp Decoder")
##################################################
# Parameters
##################################################
self.baudrate = baudrate
self.host_address = host_address
self.host_port = host_port
self.output_file = output_file
self.samp_rate = samp_rate
##################################################
# Blocks
##################################################
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(samp_rate/baudrate, 0.001, 0, 0.25, 0.001)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.blocks_udp_source_0 = blocks.udp_source(gr.sizeof_short*1, host_address, host_port, 1472, True)
self.blocks_short_to_float_0 = blocks.short_to_float(1, 32767)
self.blocks_file_sink_0 = blocks.file_sink(gr.sizeof_char*1, output_file, False)
self.blocks_file_sink_0.set_unbuffered(False)
##################################################
# Connections
##################################################
self.connect((self.blocks_short_to_float_0, 0), (self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.blocks_udp_source_0, 0), (self.blocks_short_to_float_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.blocks_file_sink_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0))
def __init__(self, samplerate, symbolrate = SYMRATE, channel_str = None, sendmsg = True, debug = False, samplepersymbol = SPS, fmdeviation = FM_DEVIATION ): gr.hier_block2.__init__(self, "pocsag", gr.io_signature(1, 1, gr.sizeof_gr_complex), gr.io_signature(1, 1, 1)) self.samplerate = samplerate self.symbolrate = symbolrate self.sendmsg = sendmsg self.debug = debug self.samplepersymbol = samplepersymbol self.fmdeviation = fmdeviation self.fractional_interpolator = gr.fractional_interpolator_cc(0, 1.0 * samplerate / (symbolrate * samplepersymbol)) self.quadrature_demod = gr.quadrature_demod_cf((symbolrate * samplepersymbol) / (fmdeviation * 4.0)) self.low_pass_filter = gr.fir_filter_fff(1, gr.firdes.low_pass(1, symbolrate * samplepersymbol, symbolrate * 2, symbolrate / 2.0, gr.firdes.WIN_HAMMING, 6.76)) self.digital_clock_recovery_mm = digital.clock_recovery_mm_ff(samplepersymbol, 0.03 * 0.03 * 0.3, 0.4, 0.03, 1e-4) self.digital_binary_slicer_fb = digital.binary_slicer_fb() self.pktdecoder = pocsag_pktdecoder(channel_str = channel_str, sendmsg = sendmsg, debug = debug) self.connect(self, self.fractional_interpolator, self.quadrature_demod, self.low_pass_filter, self.digital_clock_recovery_mm, self.digital_binary_slicer_fb, self.pktdecoder, self)
def __init__(self):
gr.top_block.__init__(self)
self.samp_rate = samp_rate = 1000000
self.dump_freq = dump_freq = 2400490000
self._u = uhd.usrp_source(
device_addr="%default",
io_type=uhd.io_type.COMPLEX_FLOAT32,
num_channels=1)
self._u.set_gain(0, 0)
self._u.set_samp_rate(self.samp_rate)
treq = uhd.tune_request(self.dump_freq, 0)
tr = self._u.set_center_freq(treq)
if tr == None:
sys.stderr.write('Failed to set center frequency\n')
raise SystemExit, 1
self.filter1 = gr.interp_fir_filter_ccf(1, firdes.low_pass(1, samp_rate, 500000, 10000, firdes.WIN_HAMMING, 6.76))
self.squelch = gr.pwr_squelch_cc(-100, 0.001, 0, True)
self.demod = gr.quadrature_demod_cf(1)
self.sync = digital.clock_recovery_mm_ff(2, 0.0076562, 0.5, 0.175, 0.005)
self.slicer = digital.binary_slicer_fb()
self.detect_seq = digital.correlate_access_code_bb("01010101010101010101010101010101", 1)
self.dump = flysky.dumpsync()
self.connect(self._u,self.filter1,self.squelch,self.demod,self.sync,self.slicer,self.detect_seq, self.dump)
def __init__(self, dc_block_len=1000, samp_per_sym=10, rx_bit_inv=1): gr.hier_block2.__init__( self, "GMSK Demodulator", gr.io_signature(1, 1, gr.sizeof_float*1), gr.io_signature(1, 1, gr.sizeof_char*1), ) ################################################## # Parameters ################################################## self.dc_block_len = dc_block_len self.samp_per_sym = samp_per_sym self.rx_bit_inv = rx_bit_inv ################################################## # Variables ################################################## self.samp_rate = samp_rate = 48000 ################################################## # Blocks ################################################## self.digital_clock_recovery_mm_xx_0_0 = digital.clock_recovery_mm_ff(samp_per_sym*(1+0.0), 0.25*0.175*0.175, 0.5, 0.175, 0.005) self.digital_binary_slicer_fb_0_0 = digital.binary_slicer_fb() self.dc_blocker_xx_0 = filter.dc_blocker_ff(dc_block_len, False) self.blocks_multiply_const_vxx_0_0 = blocks.multiply_const_vff((rx_bit_inv, )) ################################################## # Connections ################################################## self.connect((self.dc_blocker_xx_0, 0), (self.blocks_multiply_const_vxx_0_0, 0)) self.connect((self.digital_clock_recovery_mm_xx_0_0, 0), (self.digital_binary_slicer_fb_0_0, 0)) self.connect((self.blocks_multiply_const_vxx_0_0, 0), (self.digital_clock_recovery_mm_xx_0_0, 0)) self.connect((self.digital_binary_slicer_fb_0_0, 0), (self, 0)) self.connect((self, 0), (self.dc_blocker_xx_0, 0))
def test02(self):
# Test float/float version
omega = 2
gain_omega = 0.01
mu = 0.5
gain_mu = 0.01
omega_rel_lim = 0.001
self.test = digital.clock_recovery_mm_ff(omega, gain_omega,
mu, gain_mu,
omega_rel_lim)
data = 100*[1,]
self.src = blocks.vector_source_f(data, False)
self.snk = blocks.vector_sink_f()
self.tb.connect(self.src, self.test, self.snk)
self.tb.run()
expected_result = 100*[0.99972, ] # doesn't quite get to 1.0
dst_data = self.snk.data()
# Only compare last Ncmp samples
Ncmp = 30
len_e = len(expected_result)
len_d = len(dst_data)
expected_result = expected_result[len_e - Ncmp:]
dst_data = dst_data[len_d - Ncmp:]
#print expected_result
#print dst_data
self.assertFloatTuplesAlmostEqual(expected_result, dst_data, 5)
def _update_sliced(self, filtered_symbols):
self.slicer_view.data = filtered_symbols
if filtered_symbols is None:
self.sliced_view.data = None
return
omega = float(filtered_symbols.sampling_rate) / self.burst.symbol_rate
mu = 0.5
data_source = filtered_symbols.samples
numpy_source = NumpySource(data_source)
clock_recovery = digital.clock_recovery_mm_ff(omega, self._gain_omega, mu, self._gain_mu, self._omega_relative_limit)
numpy_sink = NumpySink(numpy.float32)
top = gr.top_block()
top.connect(numpy_source, clock_recovery)
top.connect(clock_recovery, numpy_sink)
top.run()
symbol_data = numpy_sink.data
# TODO: Adjust sampling rate
bits = []
for i in range(len(symbol_data)):
if symbol_data[i] >= 0:
symbol_data[i] = 1
bits.append('1')
else:
symbol_data[i] = -1
bits.append('0')
bits = ''.join(bits)
#print(bits)
self.sliced_view.data = TimeData(symbol_data, self.burst.symbol_rate)
def test04(self):
# Test float/float version
omega = 2
gain_omega = 0.01
mu = 0.25
gain_mu = 0.1
omega_rel_lim = 0.001
self.test = digital.clock_recovery_mm_ff(omega, gain_omega, mu,
gain_mu, omega_rel_lim)
data = 1000 * [1, 1, -1, -1]
self.src = blocks.vector_source_f(data, False)
self.snk = blocks.vector_sink_f()
self.tb.connect(self.src, self.test, self.snk)
self.tb.run()
expected_result = 1000 * [-1.2, 1.2]
dst_data = self.snk.data()
# Only compare last Ncmp samples
Ncmp = 100
len_e = len(expected_result)
len_d = len(dst_data)
expected_result = expected_result[len_e - Ncmp:]
dst_data = dst_data[len_d - Ncmp:]
#print expected_result
#print dst_data
self.assertFloatTuplesAlmostEqual(expected_result, dst_data, 1)
def _update_sliced(self, filtered_symbols):
self.slicer_view.data = filtered_symbols
if filtered_symbols is None:
self.sliced_view.data = None
return
omega = float(filtered_symbols.sampling_rate) / self.burst.symbol_rate
mu = 0.5
data_source = filtered_symbols.samples
numpy_source = NumpySource(data_source)
clock_recovery = digital.clock_recovery_mm_ff(omega, self._gain_omega, mu, self._gain_mu, self._omega_relative_limit)
numpy_sink = NumpySink(numpy.float32)
top = gr.top_block()
top.connect(numpy_source, clock_recovery)
top.connect(clock_recovery, numpy_sink)
top.run()
symbol_data = numpy_sink.data
numpy.max(symbol_data)
# TODO: Adjust sampling rate
bits = []
for i in xrange(len(symbol_data)):
if symbol_data[i] >= 0:
symbol_data[i] = 1
bits.append('1')
else:
symbol_data[i] = -1
bits.append('0')
bits = ''.join(bits)
#print(bits)
self.sliced_view.data = TimeData(symbol_data, self.burst.symbol_rate)
def __init__(self, sample_rate, freq_offset, queue, logger = None, group_description_csv = None):
gr.hier_block2.__init__(
self,
"SmartZone Control Channel",
gr.io_signature(1, 1, gr.sizeof_gr_complex), # input signature
gr.io_signature(0, 0, 0) # output signature
)
self._CC_DEVIATION = 4e3 # observed
self._symbol_rate = 3600.0 # control channel rate is 3.6kb/s
self._oversample = 4 # XXX reduce
# get close to the desired sample rate with decimation
channel_decimation = int(sample_rate / (self._oversample * self._symbol_rate)) & ~1
channel_rate = sample_rate / channel_decimation
samples_per_symbol = channel_rate / self._symbol_rate
# channel_bw = self._CC_DEVIATION + self._symbol_rate # from pager source
channel_bw = 3 * self._symbol_rate
# taps = firdes.low_pass(1, sample_rate, int(3.0 * self._symbol_rate), int(3.0 * self._symbol_rate / 10.0), firdes.WIN_HAMMING)
channel_taps = firdes.low_pass(1, sample_rate, channel_bw, channel_bw / 10.0, firdes.WIN_HAMMING)
channel_filter = filter.freq_xlating_fir_filter_ccf(channel_decimation, channel_taps, freq_offset, sample_rate)
#quad_demod = analog.quadrature_demod_cf(1.0)
quad_demod = analog.quadrature_demod_cf(channel_rate / (2 * math.pi * self._CC_DEVIATION))
clock = digital.clock_recovery_mm_ff(omega = samples_per_symbol, gain_omega = 0.001, mu = 0, gain_mu = 0.001, omega_relative_limit = 0.005)
slicer = digital.binary_slicer_fb()
digital_correlate = digital.correlate_access_code_bb("10101100", 0)
cc_sink = control_channel_sink(logger, queue, group_description_csv)
self.connect(self, channel_filter, quad_demod, clock, slicer, digital_correlate, cc_sink)
def set_decoder(self, chan, num):
print("Add decoder for channel %d" % num)
decoder = btbr.btdecode(num)
# Both methods are quite equivalent
if 1:
quad_demod = analog.quadrature_demod_cf(self.samp_per_sym /
(2 * math.pi))
clock_recov = digital.clock_recovery_mm_ff(
omega=self.samp_per_sym,
gain_omega=0.25 * self.gain_mu**2,
mu=0.32,
gain_mu=self.gain_mu,
omega_relative_limit=0.005)
binary_slicer = digital.binary_slicer_fb()
self.connect(chan, quad_demod)
self.connect(quad_demod, clock_recov)
self.connect(clock_recov, binary_slicer)
self.connect(binary_slicer, decoder)
self.msg_connect((decoder, 'out'), (self.bthandler, 'in'))
else:
gfsk_demod = digital.gfsk_demod(
samples_per_symbol=self.samp_per_sym,
sensitivity=2 * math.pi / self.samp_per_sym,
gain_mu=self.gain_mu,
mu=0.32,
omega_relative_limit=0.005,
freq_error=0.0,
verbose=False,
log=False,
)
self.connect(chan, gfsk_demod)
self.connect(gfsk_demod, decoder)
self.msg_connect((decoder, 'out'), (self.bthandler, 'in'))
def __init__(self, options, filename):
gr.top_block.__init__(self)
inf_str = None
symbol_rate = 152.34e3
sample_rate = 1e6
#if len(options) != 0:
# inf_str = args[0]
squelch = analog.pwr_squelch_cc(float(options.squelch), 0.1, 0, True)
demod = analog.quadrature_demod_cf(1.0)
cr = digital.clock_recovery_mm_ff(sample_rate/symbol_rate, 0.00765625, 0, 0.175, 0.005)
slicer = digital.binary_slicer_fb()
corr = digital.correlate_access_code_bb(AC, 3)
sink = sniffer()
if False:
print "Reading from: " + inf_str
src = blocks.file_source(gr.sizeof_gr_complex, inf_str, False)
else:
freqs = {
'AA':917.0e6, 'AB':913.0e6, 'AC':914.0e6, 'AD':915.0e6,
'BA':916.0e6, 'BB':919.0e6, 'BC':920.0e6, 'BD':921.0e6,
'CA':922.0e6, 'CB':923.0e6, 'CC':907.0e6, 'CD':908.0e6,
'DA':905.5e6, 'DB':909.0e6, 'DC':911.0e6, 'DD':910.0e6}
frequency = freqs[options.channel]
print "Channel: " + options.channel + " (" + str(frequency/1e6) + "MHz)"
# Create a UHD device source
src = uhd.usrp_source(device_addr=options.args, stream_args=uhd.stream_args('fc32', "sc16", args=""))
# Set the subdevice spec
if(options.spec):
src.set_subdev_spec(options.spec, 0)
# Set the antenna
if(options.antenna):
src.set_antenna(options.antenna, 0)
# Set receiver sample rate
src.set_samp_rate(options.samp_rate)
# Set receive daughterboard gain
if options.gain is None:
g = src.get_gain_range()
options.gain = float(g.start()+g.stop())/2
print "Using mid-point gain of", options.gain, "(", g.start(), "-", g.stop(), ")"
src.set_gain(options.gain)
# Set frequency (tune request takes lo_offset)
treq = uhd.tune_request(frequency)
tr = src.set_center_freq(treq)
if tr == None:
sys.stderr.write('Failed to set center frequency\n')
raise SystemExit, 1
self.connect(src, squelch, demod, cr, slicer, corr, sink)
def __init__(self, *args, **kwargs):
"""
Hierarchical block for O-QPSK demodulation.
The input is the complex modulated signal at baseband
and the output is a stream of bytes.
@param sps: samples per symbol
@type sps: integer
"""
try:
self.sps = kwargs.pop('sps')
self.log = kwargs.pop('log')
except KeyError:
pass
gr.hier_block2.__init__(
self,
"ieee802_15_4_demod",
gr.io_signature(1, 1, gr.sizeof_gr_complex), # Input
gr.io_signature(1, 1, gr.sizeof_float)) # Output
# Demodulate FM
sensitivity = (pi / 2) / self.sps
#self.fmdemod = gr.quadrature_demod_cf(1.0 / sensitivity)
self.fmdemod = gr.quadrature_demod_cf(1)
# Low pass the output of fmdemod to allow us to remove
# the DC offset resulting from frequency offset
alpha = 0.0008 / self.sps
self.freq_offset = gr.single_pole_iir_filter_ff(alpha)
self.sub = gr.sub_ff()
self.connect(self, self.fmdemod)
self.connect(self.fmdemod, (self.sub, 0))
self.connect(self.fmdemod, self.freq_offset, (self.sub, 1))
# recover the clock
omega = self.sps
gain_mu = 0.03
mu = 0.5
omega_relative_limit = 0.0002
freq_error = 0.0
gain_omega = .25 * gain_mu * gain_mu # critically damped
self.clock_recovery = digital.clock_recovery_mm_ff(
omega, gain_omega, mu, gain_mu, omega_relative_limit)
# Connect
self.connect(self.sub, self.clock_recovery, self)
if self.log:
self.connect(self.fmdemod,
gr.file_sink(gr.sizeof_float, 'rx-fmdemod.dat'))
self.connect(self.freq_offset,
gr.file_sink(gr.sizeof_float, 'rx-fo.dat'))
self.connect(self.sub, gr.file_sink(gr.sizeof_float, 'rx-sub.dat'))
self.connect(self.clock_recovery,
gr.file_sink(gr.sizeof_float, 'rx-recovery.dat'))
def __init__(self):
gr.top_block.__init__(self, "Top Block")
##################################################
# Variables
##################################################
self.variable_slider_0 = variable_slider_0 = .846
self.test = test = .005
self.shift = shift = .906
self.samp_rate_0 = samp_rate_0 = 1.2e6
self.samp_rate = samp_rate = 1.2e6/4
self.pows = pows = 1.3
self.lpf = lpf = .724
self.go = go = 0.564
self.gm = gm = 1.61
self.centre_freq = centre_freq = 439.95e6
##################################################
# Blocks
##################################################
self.rtlsdr_source_0 = osmosdr.source( args="numchan=" + str(1) + " " + "" )
self.rtlsdr_source_0.set_sample_rate(samp_rate_0)
self.rtlsdr_source_0.set_center_freq(439.9e6, 0)
self.rtlsdr_source_0.set_freq_corr(0, 0)
self.rtlsdr_source_0.set_dc_offset_mode(0, 0)
self.rtlsdr_source_0.set_iq_balance_mode(2, 0)
self.rtlsdr_source_0.set_gain_mode(False, 0)
self.rtlsdr_source_0.set_gain(20, 0)
self.rtlsdr_source_0.set_if_gain(10, 0)
self.rtlsdr_source_0.set_bb_gain(10, 0)
self.rtlsdr_source_0.set_antenna("", 0)
self.rtlsdr_source_0.set_bandwidth(0, 0)
self.low_pass_filter_1 = filter.fir_filter_fff(10, firdes.low_pass(
1, samp_rate, 2.56e3*lpf, (2.56e3/2)*lpf, firdes.WIN_HAMMING, 6.76))
self.freq_xlating_fir_filter_xxx_0 = filter.freq_xlating_fir_filter_ccc(4, (firdes.low_pass_2(1,samp_rate_0,100e3,50e3,40)), 0, samp_rate_0)
self.digital_clock_recovery_mm_xx_1 = digital.clock_recovery_mm_ff(11.6439*(1+test), 0.25*0.175*0.175*go, 0.5, 0.175*gm, 0.005*variable_slider_0)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.blocks_add_const_vxx_0 = blocks.add_const_vff((-12*shift, ))
self.blks2_tcp_sink_0 = grc_blks2.tcp_sink(
itemsize=gr.sizeof_char*1,
addr="127.0.0.1",
port=9000,
server=False,
)
self.analog_quadrature_demod_cf_0 = analog.quadrature_demod_cf(10)
self.analog_pwr_squelch_xx_0 = analog.pwr_squelch_cc(-40*pows, .001, 0, False)
##################################################
# Connections
##################################################
self.connect((self.analog_pwr_squelch_xx_0, 0), (self.analog_quadrature_demod_cf_0, 0))
self.connect((self.analog_quadrature_demod_cf_0, 0), (self.low_pass_filter_1, 0))
self.connect((self.blocks_add_const_vxx_0, 0), (self.digital_clock_recovery_mm_xx_1, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.blks2_tcp_sink_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_1, 0), (self.digital_binary_slicer_fb_0, 0))
self.connect((self.freq_xlating_fir_filter_xxx_0, 0), (self.analog_pwr_squelch_xx_0, 0))
self.connect((self.low_pass_filter_1, 0), (self.blocks_add_const_vxx_0, 0))
self.connect((self.rtlsdr_source_0, 0), (self.freq_xlating_fir_filter_xxx_0, 0))
def __init__(self, channel=26):
gr.top_block.__init__(self, "Top Block")
##################################################
# Parameters
##################################################
self.channel = channel
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 4000000
##################################################
# Blocks
##################################################
self.xmlrpc_server_0 = SimpleXMLRPCServer.SimpleXMLRPCServer(('localhost', 8080), allow_none=True)
self.xmlrpc_server_0.register_instance(self)
self.xmlrpc_server_0_thread = threading.Thread(target=self.xmlrpc_server_0.serve_forever)
self.xmlrpc_server_0_thread.daemon = True
self.xmlrpc_server_0_thread.start()
self.single_pole_iir_filter_xx_0 = filter.single_pole_iir_filter_ff(0.00016, 1)
self.rftap_rftap_encap_0 = rftap.rftap_encap(2, 195, '')
self.osmosdr_source_0 = osmosdr.source( args="numchan=" + str(1) + " " + '' )
self.osmosdr_source_0.set_sample_rate(samp_rate)
self.osmosdr_source_0.set_center_freq(1000000 * (2400 + 5 * (channel - 10)), 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(0, 0)
self.osmosdr_source_0.set_gain_mode(False, 0)
self.osmosdr_source_0.set_gain(10, 0)
self.osmosdr_source_0.set_if_gain(20, 0)
self.osmosdr_source_0.set_bb_gain(20, 0)
self.osmosdr_source_0.set_antenna('', 0)
self.osmosdr_source_0.set_bandwidth(0, 0)
self.ieee802_15_4_packet_sink_0 = ieee802_15_4.packet_sink(10)
self.epy_block_0 = epy_block_0.blk()
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(2, 0.000225, 0.5, 0.03, 0.0002)
self.blocks_sub_xx_0 = blocks.sub_ff(1)
self.blocks_socket_pdu_0_0_0 = blocks.socket_pdu("UDP_CLIENT", '127.0.0.1', '52002', 10000, False)
self.blocks_message_debug_0 = blocks.message_debug()
self.analog_quadrature_demod_cf_0 = analog.quadrature_demod_cf(1)
##################################################
# Connections
##################################################
self.msg_connect((self.epy_block_0, 'out'), (self.rftap_rftap_encap_0, 'in'))
self.msg_connect((self.ieee802_15_4_packet_sink_0, 'out'), (self.epy_block_0, 'in'))
self.msg_connect((self.rftap_rftap_encap_0, 'out'), (self.blocks_message_debug_0, 'print_pdu'))
self.msg_connect((self.rftap_rftap_encap_0, 'out'), (self.blocks_socket_pdu_0_0_0, 'pdus'))
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_sub_xx_0, 0), (self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.ieee802_15_4_packet_sink_0, 0))
self.connect((self.osmosdr_source_0, 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 __init__(self):
gr.top_block.__init__(self, "Quetzal-1 Receiver")
##################################################
# Variables
##################################################
self.symb_rate = symb_rate = 4800
self.samp_rate = samp_rate = 48000
self.samp_per_sym = samp_per_sym = int(samp_rate/symb_rate)
self.log = log = 0
self.homedir = homedir = ""
self.hhmmss = hhmmss = ""
self.gain_mu = gain_mu = 0.175*3
self.filename_raw_beacon = filename_raw_beacon = ""
self.filename_parsed_beacon = filename_parsed_beacon = ""
self.filename_image_metadata = filename_image_metadata = ""
self.filename_image = filename_image = ""
self.do_once = do_once = 0
self.data_rate = data_rate = 4800
self.count = count = 0
self.center_freq = center_freq = 437.2e6
##################################################
# Blocks
##################################################
self.zeromq_pub_sink_0 = zeromq.pub_sink(gr.sizeof_char, 1, 'tcp://127.0.0.1:1502', 100, False, -1)
self.satellites_strip_ax25_header_0 = satellites.strip_ax25_header()
self.satellites_nrzi_decode_0 = satellites.nrzi_decode()
self.satellites_hdlc_deframer_0 = satellites.hdlc_deframer(check_fcs=True, max_length=10000)
self.quetzal1_parse = quetzal1_parse.quetzal1_parse(filename_parsed_beacon='', filename_raw_beacon='')
self.low_pass_filter_0_0 = filter.fir_filter_fff(1, firdes.low_pass(
1, 48000, 2400, 2000, firdes.WIN_HAMMING, 6.76))
self.digital_descrambler_bb_0 = digital.descrambler_bb(0x21, 0, 16)
self.digital_clock_recovery_mm_xx_0_0 = digital.clock_recovery_mm_ff(10, 0.25*gain_mu*gain_mu, 0.5, gain_mu, 0.005)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.blocks_wavfile_source_0 = blocks.wavfile_source('/home/dan/Documents/repos/gr-quetzal1/recordings/example_beacon_quetzal1.wav', False)
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_float*1, samp_rate,True)
self.blocks_pdu_to_tagged_stream_0 = blocks.pdu_to_tagged_stream(blocks.byte_t, 'packet_len')
self.blocks_multiply_const_vxx_0 = blocks.multiply_const_vff((5, ))
self.blocks_add_const_vxx_0 = blocks.add_const_vff((0, ))
##################################################
# Connections
##################################################
self.msg_connect((self.satellites_hdlc_deframer_0, 'out'), (self.satellites_strip_ax25_header_0, 'in'))
self.msg_connect((self.satellites_strip_ax25_header_0, 'out'), (self.blocks_pdu_to_tagged_stream_0, 'pdus'))
self.msg_connect((self.satellites_strip_ax25_header_0, 'out'), (self.quetzal1_parse, 'in'))
self.connect((self.blocks_add_const_vxx_0, 0), (self.blocks_multiply_const_vxx_0, 0))
self.connect((self.blocks_multiply_const_vxx_0, 0), (self.low_pass_filter_0_0, 0))
self.connect((self.blocks_pdu_to_tagged_stream_0, 0), (self.zeromq_pub_sink_0, 0))
self.connect((self.blocks_throttle_0, 0), (self.blocks_add_const_vxx_0, 0))
self.connect((self.blocks_wavfile_source_0, 0), (self.blocks_throttle_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.satellites_nrzi_decode_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0_0, 0), (self.digital_binary_slicer_fb_0, 0))
self.connect((self.digital_descrambler_bb_0, 0), (self.satellites_hdlc_deframer_0, 0))
self.connect((self.low_pass_filter_0_0, 0), (self.digital_clock_recovery_mm_xx_0_0, 0))
self.connect((self.satellites_nrzi_decode_0, 0), (self.digital_descrambler_bb_0, 0))
def __init__(self):
grc_wxgui.top_block_gui.__init__(self, title="Top Block")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 1000000
self.dec_rate = dec_rate = 2
##################################################
# Blocks
##################################################
self.uhd_usrp_source_0 = uhd.usrp_source(
device_addr="",
stream_args=uhd.stream_args(
cpu_format="fc32",
channels=range(1),
),
)
self.uhd_usrp_source_0.set_samp_rate(1000000)
self.uhd_usrp_source_0.set_center_freq(2400490000, 0)
self.uhd_usrp_source_0.set_gain(0, 0)
self.low_pass_filter_0_0 = gr.interp_fir_filter_ccf(
1,
firdes.low_pass(1, samp_rate, 500000, 10000, firdes.WIN_HAMMING,
6.76))
self.low_pass_filter_0 = gr.fir_filter_fff(
1,
firdes.low_pass(1, samp_rate, 500000, 10000, firdes.WIN_HAMMING,
6.76))
self.gr_quadrature_demod_cf_0 = gr.quadrature_demod_cf(1)
self.gr_pwr_squelch_xx_0 = gr.pwr_squelch_cc(-100, 0.001, 0, True)
self.flysky_dumpsync_0 = flysky.dumpsync()
self.digital_correlate_access_code_bb_0_1 = digital.correlate_access_code_bb(
"010101010101010101010101010101010101010001110101", 1)
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(
2, 0.0076562, 0.5, 0.175, 0.005)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
##################################################
# Connections
##################################################
self.connect((self.digital_clock_recovery_mm_xx_0, 0),
(self.digital_binary_slicer_fb_0, 0))
self.connect((self.gr_quadrature_demod_cf_0, 0),
(self.low_pass_filter_0, 0))
self.connect((self.low_pass_filter_0, 0),
(self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.gr_pwr_squelch_xx_0, 0),
(self.gr_quadrature_demod_cf_0, 0))
self.connect((self.uhd_usrp_source_0, 0),
(self.low_pass_filter_0_0, 0))
self.connect((self.low_pass_filter_0_0, 0),
(self.gr_pwr_squelch_xx_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0),
(self.digital_correlate_access_code_bb_0_1, 0))
self.connect((self.digital_correlate_access_code_bb_0_1, 0),
(self.flysky_dumpsync_0, 0))
def __init__(self, source_data, sampling_rate, carrier_hz, symbol_rate, deviation, access_code):
super(FSKDemodulator, self).__init__()
self._decoded = {}
self._carrier_hz = carrier_hz
self._deviation = deviation
self._access_code = access_code
samp_rate = sampling_rate
#symbol_rate = 9920
self.samples_per_symbol = float(samp_rate) / symbol_rate
omega = self.samples_per_symbol * 1.0
mu = 0.0
gain_mu = 0.2
gain_omega = 0.25 * gain_mu * gain_mu
omega_relative_limit = 0.001
tap_count = int(math.floor(self.samples_per_symbol))
hz_n = (carrier_hz - deviation)
taps_n = numpy.exp(numpy.arange(tap_count, dtype=numpy.float32) * 2.0j * numpy.pi * hz_n / samp_rate)
hz_p = (carrier_hz + deviation)
taps_p = numpy.exp(numpy.arange(tap_count, dtype=numpy.float32) * 2.0j * numpy.pi * hz_p / samp_rate)
#source = blocks.file_source(gr.sizeof_gr_complex*1, filepath_in, False)
# Concatenate data to compensate for correlate_access_code_bb latency
source_data_padding_count = int(math.ceil(self.samples_per_symbol * 64))
source_data = numpy.concatenate((source_data, numpy.zeros((source_data_padding_count,), dtype=numpy.complex64)))
source = NumpySource(source_data)
filter_n = filter.fir_filter_ccc(1, taps_n.tolist())
self.connect(source, filter_n)
filter_p = filter.fir_filter_ccc(1, taps_p.tolist())
self.connect(source, filter_p)
mag_n = blocks.complex_to_mag(1)
self.connect(filter_n, mag_n)
mag_p = blocks.complex_to_mag(1)
self.connect(filter_p, mag_p)
sub_pn = blocks.sub_ff()
self.connect(mag_p, (sub_pn, 0))
self.connect(mag_n, (sub_pn, 1))
clock_recovery = digital.clock_recovery_mm_ff(omega, gain_omega, mu, gain_mu, omega_relative_limit)
self.connect(sub_pn, clock_recovery)
slicer = digital.binary_slicer_fb()
self.connect(clock_recovery, slicer)
access_code_correlator = digital.correlate_access_code_bb(access_code, 0)
self.connect(slicer, access_code_correlator)
self.packetizer = Packetizer()
self.connect(access_code_correlator, self.packetizer)
def __init__(self):
gr.top_block.__init__(self, "Afsk Csu")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 44100
self.bps = bps = 1200
self.window_len = window_len = (samp_rate / bps) * 2
self.dec = dec = 4
self.window_sync = window_sync = signal.windows.hann(116)
self.window = window = signal.windows.cosine(window_len)
self.dec_samp_rate = dec_samp_rate = samp_rate / dec
self.bpf_width = bpf_width = bps + 200
self.bpf_trans = bpf_trans = 200
##################################################
# Blocks
##################################################
self.freq_xlating_fir_filter_xxx_0_0 = filter.freq_xlating_fir_filter_fcf(
dec, (window), (2200 + 1200) / 2, samp_rate)
self.digital_diff_decoder_bb_0 = digital.diff_decoder_bb(2)
self.digital_clock_recovery_mm_xx_0_0 = digital.clock_recovery_mm_ff(
9.1875, 0.28, 0, 0.024, 0.01)
self.digital_binary_slicer_fb_0_0 = digital.binary_slicer_fb()
self.blocks_wavfile_source_0 = blocks.wavfile_source(
'/tmp/gs-sample.wav', False)
self.blocks_multiply_const_vxx_0 = blocks.multiply_const_vff((1, ))
self.blocks_file_sink_1_0_1_0 = blocks.file_sink(
gr.sizeof_char * 1, '/tmp/gs-sample.binary', False)
self.blocks_file_sink_1_0_1_0.set_unbuffered(True)
self.band_pass_filter_0 = filter.fir_filter_fff(
1,
firdes.band_pass(1, samp_rate, 1700 - bpf_width, 1700 + bpf_width,
bpf_trans, firdes.WIN_HAMMING, 6.76))
self.analog_quadrature_demod_cf_0 = analog.quadrature_demod_cf(
(dec_samp_rate / bps) / (math.pi * 0.5))
##################################################
# Connections
##################################################
self.connect((self.analog_quadrature_demod_cf_0, 0),
(self.blocks_multiply_const_vxx_0, 0))
self.connect((self.band_pass_filter_0, 0),
(self.freq_xlating_fir_filter_xxx_0_0, 0))
self.connect((self.blocks_multiply_const_vxx_0, 0),
(self.digital_clock_recovery_mm_xx_0_0, 0))
self.connect((self.blocks_wavfile_source_0, 0),
(self.band_pass_filter_0, 0))
self.connect((self.digital_binary_slicer_fb_0_0, 0),
(self.digital_diff_decoder_bb_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0_0, 0),
(self.digital_binary_slicer_fb_0_0, 0))
self.connect((self.digital_diff_decoder_bb_0, 0),
(self.blocks_file_sink_1_0_1_0, 0))
self.connect((self.freq_xlating_fir_filter_xxx_0_0, 0),
(self.analog_quadrature_demod_cf_0, 0))
def __init__(self, source_data, sampling_rate, carrier_hz, symbol_rate, deviation, access_code):
super(FSKDemodulator, self).__init__()
self._decoded = {}
self._carrier_hz = carrier_hz
self._deviation = deviation
self._access_code = access_code
samp_rate = sampling_rate
#symbol_rate = 9920
self.samples_per_symbol = float(samp_rate) / symbol_rate
omega = self.samples_per_symbol * 1.0
mu = 0.0
gain_mu = 0.2
gain_omega = 0.25 * gain_mu * gain_mu
omega_relative_limit = 0.001
tap_count = int(math.floor(self.samples_per_symbol))
hz_n = (carrier_hz - deviation)
taps_n = numpy.exp(numpy.arange(tap_count, dtype=numpy.float32) * 2.0j * numpy.pi * hz_n / samp_rate)
hz_p = (carrier_hz + deviation)
taps_p = numpy.exp(numpy.arange(tap_count, dtype=numpy.float32) * 2.0j * numpy.pi * hz_p / samp_rate)
#source = blocks.file_source(gr.sizeof_gr_complex*1, filepath_in, False)
# Concatenate data to compensate for correlate_access_code_bb latency
source_data_padding_count = int(math.ceil(self.samples_per_symbol * 64))
source_data = numpy.concatenate((source_data, numpy.zeros((source_data_padding_count,), dtype=numpy.complex64)))
source = NumpySource(source_data)
filter_n = filter.fir_filter_ccc(1, taps_n.tolist())
self.connect(source, filter_n)
filter_p = filter.fir_filter_ccc(1, taps_p.tolist())
self.connect(source, filter_p)
mag_n = blocks.complex_to_mag(1)
self.connect(filter_n, mag_n)
mag_p = blocks.complex_to_mag(1)
self.connect(filter_p, mag_p)
sub_pn = blocks.sub_ff()
self.connect(mag_p, (sub_pn, 0))
self.connect(mag_n, (sub_pn, 1))
clock_recovery = digital.clock_recovery_mm_ff(omega, gain_omega, mu, gain_mu, omega_relative_limit)
self.connect(sub_pn, clock_recovery)
slicer = digital.binary_slicer_fb()
self.connect(clock_recovery, slicer)
access_code_correlator = digital.correlate_access_code_bb(access_code, 0)
self.connect(slicer, access_code_correlator)
self.packetizer = Packetizer()
self.connect(access_code_correlator, self.packetizer)
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 __init__(self,
samp_rate=1600000,
samp_per_sym=16,
freq_error=-0.0025000):
gr.hier_block2.__init__(
self,
"Wireless M-Bus Demod",
gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
gr.io_signature(1, 1, gr.sizeof_char * 1),
)
##################################################
# Parameters
##################################################
self.samp_rate = samp_rate
self.samp_per_sym = samp_per_sym
self.freq_error = freq_error
##################################################
# Variables
##################################################
self.cutoff = cutoff = 120e3
self.chip_rate = chip_rate = samp_rate / samp_per_sym
##################################################
# Blocks
##################################################
self.low_pass_filter_0 = gnuradio.filter.fir_filter_ccf(
1,
gnuradio.filter.firdes.low_pass(1, samp_rate, cutoff, cutoff / 2,
gnuradio.filter.firdes.WIN_HAMMING,
6.76))
self.gr_sub_xx_0 = blocks.sub_ff(1)
self.gr_single_pole_iir_filter_xx_0 = gnuradio.filter.single_pole_iir_filter_ff(
0.0512 / samp_per_sym, 1)
self.gr_quadrature_demod_cf_0 = analog.quadrature_demod_cf(1)
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(
samp_per_sym * (1 + freq_error), .25 * 0.06 * 0.06 * 4, 0.5,
0.06 * 2, 0.002 * 2)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
##################################################
# Connections
##################################################
self.connect((self.gr_quadrature_demod_cf_0, 0),
(self.gr_single_pole_iir_filter_xx_0, 0))
self.connect((self.gr_single_pole_iir_filter_xx_0, 0),
(self.gr_sub_xx_0, 1))
self.connect((self.gr_quadrature_demod_cf_0, 0), (self.gr_sub_xx_0, 0))
self.connect((self.gr_sub_xx_0, 0),
(self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0),
(self.digital_binary_slicer_fb_0, 0))
self.connect((self.low_pass_filter_0, 0),
(self.gr_quadrature_demod_cf_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self, 0))
self.connect((self, 0), (self.low_pass_filter_0, 0))
def __init__(self, callsign='', invert=1, ip='::', latitude=0, longitude=0, port=7355, recstart=''):
gr.top_block.__init__(self, "TW-1B decoder")
##################################################
# Parameters
##################################################
self.callsign = callsign
self.invert = invert
self.ip = ip
self.latitude = latitude
self.longitude = longitude
self.port = port
self.recstart = recstart
##################################################
# Variables
##################################################
self.threshold = threshold = 6
self.gain_mu = gain_mu = 0.175*3
self.access_code = access_code = "10010011000010110101000111011110"
##################################################
# Blocks
##################################################
self.sync_to_pdu_packed_0 = sync_to_pdu_packed(
packlen=256,
sync=access_code,
threshold=threshold,
)
self.sids_submit_0 = sids.submit('http://tlm.pe0sat.nl/tlmdb/frame_db.php', 40927, callsign, longitude, latitude, recstart)
self.sids_print_timestamp_0 = sids.print_timestamp('%Y-%m-%d %H:%M:%S')
self.low_pass_filter_0 = filter.fir_filter_fff(1, firdes.low_pass(
1, 48000, 2400, 2000, firdes.WIN_HAMMING, 6.76))
self.digital_descrambler_bb_0 = digital.descrambler_bb(0x21, 0x00, 16)
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(10, 0.25*gain_mu*gain_mu, 0.5, gain_mu, 0.005)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.blocks_udp_source_0 = blocks.udp_source(gr.sizeof_short*1, ip, port, 1472, False)
self.blocks_short_to_float_0 = blocks.short_to_float(1, 32767.0)
self.blocks_multiply_const_vxx_0 = blocks.multiply_const_vff((invert*10, ))
self.blocks_message_debug_0 = blocks.message_debug()
self.ax100_gomx3_rs_decode_0 = ax100.gomx3_rs_decode(False)
##################################################
# Connections
##################################################
self.msg_connect((self.ax100_gomx3_rs_decode_0, 'out'), (self.sids_print_timestamp_0, 'in'))
self.msg_connect((self.ax100_gomx3_rs_decode_0, 'out'), (self.sids_submit_0, 'in'))
self.msg_connect((self.sids_print_timestamp_0, 'out'), (self.blocks_message_debug_0, 'print_pdu'))
self.msg_connect((self.sync_to_pdu_packed_0, 'out'), (self.ax100_gomx3_rs_decode_0, 'in'))
self.connect((self.blocks_multiply_const_vxx_0, 0), (self.low_pass_filter_0, 0))
self.connect((self.blocks_short_to_float_0, 0), (self.blocks_multiply_const_vxx_0, 0))
self.connect((self.blocks_udp_source_0, 0), (self.blocks_short_to_float_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.digital_descrambler_bb_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0))
self.connect((self.digital_descrambler_bb_0, 0), (self.sync_to_pdu_packed_0, 0))
self.connect((self.low_pass_filter_0, 0), (self.digital_clock_recovery_mm_xx_0, 0))
def __init__(self, freq, ppm, osmosdr_args):
gr.top_block.__init__(self, "gr_omnicon")
self.msgq_out = blocks_message_sink_0_msgq_out = gr.msg_queue(0)
##################################################
# Variables
##################################################
self.xlate_bandwidth = xlate_bandwidth = 15.0e3
self.samp_rate = samp_rate = 1.2e6
self.xlate_decimation = xlate_decimation = int(samp_rate/(xlate_bandwidth*3.2))
self.baud_rate = baud_rate = 2400
self.lowpass_decimation = lowpass_decimation = int((samp_rate/xlate_decimation)/(baud_rate*4))
self.freq_offset = freq_offset = 250000
self.sps = sps = (samp_rate/xlate_decimation/lowpass_decimation)/baud_rate
self.omega_rel_limit = omega_rel_limit = ((2450.0-2400.0)/2400.0)
self.gain_omega = gain_omega = 0
self.gain_mu = gain_mu = 0.1
self.freq_tune = freq_tune = freq-freq_offset
##################################################
# Blocks
##################################################
self.rtlsdr_source_0 = osmosdr.source(args=osmosdr_args)
self.rtlsdr_source_0.set_sample_rate(samp_rate)
self.rtlsdr_source_0.set_center_freq(freq_tune, 0)
self.rtlsdr_source_0.set_freq_corr(ppm, 0)
self.rtlsdr_source_0.set_dc_offset_mode(0, 0)
self.rtlsdr_source_0.set_iq_balance_mode(0, 0)
self.rtlsdr_source_0.set_gain_mode(True, 0)
self.rtlsdr_source_0.set_gain(42, 0)
self.rtlsdr_source_0.set_if_gain(10, 0)
self.rtlsdr_source_0.set_bb_gain(10, 0)
self.rtlsdr_source_0.set_antenna("", 0)
self.rtlsdr_source_0.set_bandwidth(0, 0)
self.low_pass_filter_0 = filter.fir_filter_fff(lowpass_decimation, firdes.low_pass(
1, samp_rate/xlate_decimation, baud_rate, (baud_rate)/10, firdes.WIN_HAMMING, 6.76))
self.freq_xlating_fir_filter_xxx_0 = filter.freq_xlating_fir_filter_ccc(xlate_decimation, (firdes.low_pass(1, samp_rate, xlate_bandwidth, xlate_bandwidth/20 )), freq_offset, samp_rate)
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(sps, gain_omega, 0.5, gain_mu, omega_rel_limit)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.dc_blocker_xx_0 = filter.dc_blocker_ff(500, True)
self.blocks_message_sink_0 = blocks.message_sink(gr.sizeof_char*1, blocks_message_sink_0_msgq_out, False)
self.analog_quadrature_demod_cf_0_0_0 = analog.quadrature_demod_cf(1)
##################################################
# Connections
##################################################
self.connect((self.analog_quadrature_demod_cf_0_0_0, 0), (self.dc_blocker_xx_0, 0))
#self.connect((self.blocks_message_sink_0, 'msg'), (self, 0))
self.connect((self.dc_blocker_xx_0, 0), (self.low_pass_filter_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.blocks_message_sink_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0))
self.connect((self.freq_xlating_fir_filter_xxx_0, 0), (self.analog_quadrature_demod_cf_0_0_0, 0))
self.connect((self.low_pass_filter_0, 0), (self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.rtlsdr_source_0, 0), (self.freq_xlating_fir_filter_xxx_0, 0))
def __init__(self, *args, **kwargs):
"""
Hierarchical block for O-QPSK demodulation.
The input is the complex modulated signal at baseband
and the output is a stream of bytes.
@param sps: samples per symbol
@type sps: integer
"""
try:
self.sps = kwargs.pop('sps')
self.log = kwargs.pop('log')
except KeyError:
pass
gr.hier_block2.__init__(self, "ieee802_15_4_demod",
gr.io_signature(1, 1, gr.sizeof_gr_complex), # Input
gr.io_signature(1, 1, gr.sizeof_float)) # Output
# Demodulate FM
sensitivity = (pi / 2) / self.sps
#self.fmdemod = gr.quadrature_demod_cf(1.0 / sensitivity)
self.fmdemod = gr.quadrature_demod_cf(1)
# Low pass the output of fmdemod to allow us to remove
# the DC offset resulting from frequency offset
alpha = 0.0008/self.sps
self.freq_offset = gr.single_pole_iir_filter_ff(alpha)
self.sub = gr.sub_ff()
self.connect(self, self.fmdemod)
self.connect(self.fmdemod, (self.sub, 0))
self.connect(self.fmdemod, self.freq_offset, (self.sub, 1))
# recover the clock
omega = self.sps
gain_mu=0.03
mu=0.5
omega_relative_limit=0.0002
freq_error=0.0
gain_omega = .25*gain_mu*gain_mu # critically damped
self.clock_recovery = digital.clock_recovery_mm_ff(omega, gain_omega, mu, gain_mu,
omega_relative_limit)
# Connect
self.connect(self.sub, self.clock_recovery, self)
if self.log:
self.connect(self.fmdemod, gr.file_sink(gr.sizeof_float, 'rx-fmdemod.dat'))
self.connect(self.freq_offset, gr.file_sink(gr.sizeof_float, 'rx-fo.dat'))
self.connect(self.sub, gr.file_sink(gr.sizeof_float, 'rx-sub.dat'))
self.connect(self.clock_recovery, gr.file_sink(gr.sizeof_float, 'rx-recovery.dat'))
def __init__(self):
gr.top_block.__init__(self, "FSK Demod Demo")
# Variables
self.symbol_rate = symbol_rate = 125e3
self.samp_rate = samp_rate = symbol_rate
self.f_center = f_center = 868e6
self.sps = sps = 2
self.sensitivity = sensitivity = (pi / 2) / sps
self.alpha = alpha = 0.0512 / sps
self.bandwidth = bandwidth = 100e3
# Blocks
self.uhd_usrp_source_0 = uhd.usrp_source(
device_addr="",
stream_args=uhd.stream_args(
cpu_format="fc32",
channels=range(1),
),
)
self.uhd_usrp_source_0.set_samp_rate(samp_rate)
self.uhd_usrp_source_0.set_center_freq(f_center, 0)
self.uhd_usrp_source_0.set_gain(0, 0)
self.uhd_usrp_source_0.set_bandwidth(bandwidth, 0)
self.fm_demod = gr.quadrature_demod_cf(1 / sensitivity)
self.freq_offset = gr.single_pole_iir_filter_ff(alpha)
self.sub = gr.sub_ff()
self.add = gr.add_ff()
self.multiply = gr.multiply_ff()
self.invert = gr.multiply_const_vff((-1, ))
# recover the clock
omega = sps
gain_mu = 0.03
mu = 0.5
omega_relative_limit = 0.0002
freq_error = 0.0
gain_omega = .25 * gain_mu * gain_mu # critically damped
self.clock_recovery = digital.clock_recovery_mm_ff(
omega, gain_omega, mu, gain_mu, omega_relative_limit)
self.slice = digital.binary_slicer_fb()
self.sink = gr.vector_sink_b(1)
self.file_sink = gr.file_sink(gr.sizeof_char, 'fsk_dump.log')
# Connections
self.connect(self.fm_demod, (self.add, 0))
self.connect(self.fm_demod, self.freq_offset, (self.add, 1))
self.connect(self.uhd_usrp_source_0, self.fm_demod)
self.connect(self.add, self.clock_recovery, self.invert, self.slice,
self.file_sink)
self.connect(self.slice, self.sink)
def __init__(self):
grc_wxgui.top_block_gui.__init__(self, title="Top Block")
_icon_path = "/usr/share/icons/hicolor/32x32/apps/gnuradio-grc.png"
self.SetIcon(wx.Icon(_icon_path, wx.BITMAP_TYPE_ANY))
##################################################
# Variables
##################################################
self.lowpass2 = lowpass2 = 10e3
self.lowpass = lowpass = 10e3
self.base = base = "/home/ali/Documents/UTAT/GNURadio/2018-12-28/es_pipe_U@?U@?U@?"
self.samp_rate = samp_rate = 4003200
self.filename_2_0 = filename_2_0 = base+"_mmclock"
self.filename_2 = filename_2 = base+"_bits"
self.filename_1 = filename_1 = base+"_quaddemod_lpf_"+str((lowpass/1000.0))+"lpf2_"+str((lowpass2/1000.0))
self.filename_0 = filename_0 = base+"_quaddemod_lpf_"+str((lowpass/1000.0))
self.filename = filename = base+"_lpf_"+str((lowpass/1000.0))
##################################################
# Blocks
##################################################
self.low_pass_filter_1 = filter.fir_filter_fff(1, firdes.low_pass(
1, samp_rate, lowpass2, lowpass2, firdes.WIN_HAMMING, 6.76))
self.low_pass_filter_0 = filter.fir_filter_ccf(1, firdes.low_pass(
1, samp_rate, lowpass, lowpass/3, firdes.WIN_HAMMING, 6.76))
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(417*(1+0.0), 0.25*0.175*0.175, 0.5, 0.175, 0.005)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_gr_complex*1, samp_rate,True)
self.blocks_pack_k_bits_bb_1 = blocks.pack_k_bits_bb(8)
self.blocks_multiply_xx_0 = blocks.multiply_vcc(1)
self.blocks_file_source_0 = blocks.file_source(gr.sizeof_gr_complex*1, base, False)
self.blocks_file_sink_0_2 = blocks.file_sink(gr.sizeof_float*1, filename_2_0, False)
self.blocks_file_sink_0_2.set_unbuffered(False)
self.blocks_file_sink_0_1_0 = blocks.file_sink(gr.sizeof_char*1, filename_2, False)
self.blocks_file_sink_0_1_0.set_unbuffered(False)
self.blocks_add_const_vxx_0 = blocks.add_const_vff((-15, ))
self.analog_sig_source_x_0 = analog.sig_source_c(samp_rate, analog.GR_COS_WAVE, -1e6, 1, 0)
self.analog_quadrature_demod_cf_0 = analog.quadrature_demod_cf(2123.8262)
##################################################
# Connections
##################################################
self.connect((self.analog_quadrature_demod_cf_0, 0), (self.low_pass_filter_1, 0))
self.connect((self.analog_sig_source_x_0, 0), (self.blocks_multiply_xx_0, 1))
self.connect((self.blocks_add_const_vxx_0, 0), (self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.blocks_file_source_0, 0), (self.blocks_throttle_0, 0))
self.connect((self.blocks_multiply_xx_0, 0), (self.low_pass_filter_0, 0))
self.connect((self.blocks_pack_k_bits_bb_1, 0), (self.blocks_file_sink_0_1_0, 0))
self.connect((self.blocks_throttle_0, 0), (self.blocks_multiply_xx_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.blocks_pack_k_bits_bb_1, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.blocks_file_sink_0_2, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0))
self.connect((self.low_pass_filter_0, 0), (self.analog_quadrature_demod_cf_0, 0))
self.connect((self.low_pass_filter_1, 0), (self.blocks_add_const_vxx_0, 0))
def __init__(self):
gr.top_block.__init__(self, "grc_cc1111_test_local_loop")
##################################################
# Variables
##################################################
self.symbole_rate = symbole_rate = 80000
self.samp_rate = samp_rate = 2e6
self.samp_per_sym = samp_per_sym = int(samp_rate / symbole_rate)
self.preamble = preamble = '0101010101010101'
self.myqueue_out = myqueue_out = gr.msg_queue(2)
self.myqueue_in = myqueue_in = gr.msg_queue(2)
self.bit_per_sym = bit_per_sym = 1
self.access_code = access_code = '11010011100100011101001110010001'
##################################################
# Blocks
##################################################
self.digital_gmsk_mod_0 = digital.gmsk_mod(
samples_per_symbol=int(samp_per_sym),
bt=0.5,
verbose=False,
log=False,
)
self.digital_correlate_access_code_bb_0_0 = digital.correlate_access_code_bb(access_code, 1)
self.digital_clock_recovery_mm_xx_0_0 = digital.clock_recovery_mm_ff(samp_per_sym*(1+0.0), 0.25*0.175*0.175, 0.5, 0.175, 0.005)
self.digital_binary_slicer_fb_0_0_0 = digital.binary_slicer_fb()
self.cc1111_cc1111_packet_encoder_0 = cc1111.cc1111_packet_mod_base(cc1111.cc1111_packet_encoder(
samples_per_symbol=samp_per_sym,
bits_per_symbol=bit_per_sym,
preamble=preamble,
access_code=access_code,
pad_for_usrp=True,
do_whitening=True,
add_crc=True
),
source_queue=myqueue_in
)
self.cc1111_cc1111_packet_decoder_0 = cc1111.cc1111_packet_decoder(myqueue_out,True, True, False, True)
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_char*1, samp_rate,True)
self.blocks_null_sink_0 = blocks.null_sink(gr.sizeof_char*1)
self.analog_quadrature_demod_cf_0_0 = analog.quadrature_demod_cf(1)
##################################################
# Connections
##################################################
self.connect((self.digital_clock_recovery_mm_xx_0_0, 0), (self.digital_binary_slicer_fb_0_0_0, 0))
self.connect((self.analog_quadrature_demod_cf_0_0, 0), (self.digital_clock_recovery_mm_xx_0_0, 0))
self.connect((self.digital_binary_slicer_fb_0_0_0, 0), (self.digital_correlate_access_code_bb_0_0, 0))
self.connect((self.digital_gmsk_mod_0, 0), (self.analog_quadrature_demod_cf_0_0, 0))
self.connect((self.digital_correlate_access_code_bb_0_0, 0), (self.cc1111_cc1111_packet_decoder_0, 0))
self.connect((self.cc1111_cc1111_packet_decoder_0, 0), (self.blocks_null_sink_0, 0))
self.connect((self.blocks_throttle_0, 0), (self.digital_gmsk_mod_0, 0))
self.connect((self.cc1111_cc1111_packet_encoder_0, 0), (self.blocks_throttle_0, 0))
def __init__(self, options):
gr.hier_block2.__init__(self, "fsk_demod",
gr.io_signature(1, 1, gr.sizeof_gr_complex), # Input signature
gr.io_signature(1, 1, gr.sizeof_char)) # Output signature
self._syms_per_sec = options.syms_per_sec # ditto
self._samples_per_second = options.samples_per_second
self._gain_mu = options.gain_mu # for the clock recovery block
self._mu = options.mu
self._omega_relative_limit = options.omega_relative_limit
self._freqoffset = options.offset
#first bring that input stream down to a manageable level, let's say 3 samples per bit.
self._clockrec_oversample = 3
self._downsampletaps = gr.firdes.low_pass(1, self._samples_per_second, 10000, 1000, firdes.WIN_HANN)
self._decim = int(self._samples_per_second / (self._syms_per_sec * self._clockrec_oversample))
print "Demodulator decimation: %i" % (self._decim,)
self._downsample = gr.freq_xlating_fir_filter_ccf(self._decim, #decimation
self._downsampletaps, #taps
self._freqoffset, #freq offset
self._samples_per_second) #sampling rate
#using a pll to demod gets you a nice IIR LPF response for free
self._demod = gr.pll_freqdet_cf(2.0 / self._clockrec_oversample, #gain alpha, rad/samp
2*pi/self._clockrec_oversample, #max freq, rad/samp
-2*pi/self._clockrec_oversample) #min freq, rad/samp
self._sps = float(self._samples_per_second)/self._decim/self._syms_per_sec
#band edge filter FLL with a low bandwidth is very good
#at synchronizing to continuous FSK signals
self._carriertrack = digital.fll_band_edge_cc(self._sps,
0.6, #rolloff factor
64, #taps
1.0) #loop bandwidth
print "Samples per symbol: %f" % (self._sps,)
self._softbits = digital.clock_recovery_mm_ff(self._sps,
0.25*self._gain_mu*self._gain_mu, #gain omega, = mu/2 * mu_gain^2
self._mu, #mu (decision threshold)
self._gain_mu, #mu gain
self._omega_relative_limit) #omega relative limit
self._subtract = gr.sub_ff()
self._slicer = digital.binary_slicer_fb()
self.connect(self, self._downsample, self._carriertrack, self._demod, self._softbits, self._slicer, self)
def __init__(self): gr.top_block.__init__(self, "FSK Demod Demo") # Variables self.symbol_rate = symbol_rate = 125e3 self.samp_rate = samp_rate = symbol_rate self.f_center = f_center = 868e6 self.sps = sps = 2 self.sensitivity = sensitivity = (pi / 2) / sps self.alpha = alpha = 0.0512/sps self.bandwidth = bandwidth = 100e3 # Blocks self.uhd_usrp_source_0 = uhd.usrp_source( device_addr="", stream_args=uhd.stream_args( cpu_format="fc32", channels=range(1), ), ) self.uhd_usrp_source_0.set_samp_rate(samp_rate) self.uhd_usrp_source_0.set_center_freq(f_center, 0) self.uhd_usrp_source_0.set_gain(0, 0) self.uhd_usrp_source_0.set_bandwidth(bandwidth, 0) self.fm_demod = gr.quadrature_demod_cf(1 / sensitivity) self.freq_offset = gr.single_pole_iir_filter_ff(alpha) self.sub = gr.sub_ff() self.add = gr.add_ff() self.multiply = gr.multiply_ff() self.invert = gr.multiply_const_vff((-1, )) # recover the clock omega = sps gain_mu = 0.03 mu = 0.5 omega_relative_limit = 0.0002 freq_error = 0.0 gain_omega = .25 * gain_mu * gain_mu # critically damped self.clock_recovery = digital.clock_recovery_mm_ff(omega, gain_omega, mu, gain_mu, omega_relative_limit) self.slice = digital.binary_slicer_fb() self.sink = gr.vector_sink_b(1) self.file_sink = gr.file_sink(gr.sizeof_char, 'fsk_dump.log') # Connections self.connect(self.fm_demod, (self.add, 0)) self.connect(self.fm_demod, self.freq_offset, (self.add, 1)) self.connect(self.uhd_usrp_source_0, self.fm_demod) self.connect(self.add, self.clock_recovery, self.invert, self.slice, self.file_sink) self.connect(self.slice, self.sink)
def __init__(self, sps, gain_mu):
gr.hier_block2.__init__(self, "fsk_demod",
gr.io_signature(1, 1, gr.sizeof_gr_complex), # Input signature
gr.io_signature(1, 1, gr.sizeof_char)) # Output signature
self._sps = float(sps)
self._gain_mu = gain_mu # for the clock recovery block
self._mu = 0.5
self._omega_relative_limit = 0.35
#first bring that input stream down to a manageable level
self._clockrec_oversample = 3.0
self._decim = self._sps / self._clockrec_oversample
print "Demodulator decimation: %f" % self._decim
self._downsampletaps = filter.firdes.low_pass(1.0/self._decim,
1.0, 0.4,
0.05, filter.firdes.WIN_HANN)
# self._downsample = filter.fft_filter_ccc(self._decim,
# self._downsampletaps)
#sure this works but it's a little heavy on the CPU at high rates
self._downsample = pfb.arb_resampler_ccf(1/self._decim)
self._clockrec_sps = self._sps / self._decim
#using a pll to demod gets you a nice IIR LPF response for free
self._demod = analog.pll_freqdet_cf(2.0 / self._clockrec_sps, #gain alpha, rad/samp
2*pi/self._clockrec_sps, #max freq, rad/samp
-2*pi/self._clockrec_sps) #min freq, rad/samp
#band edge filter FLL with a low bandwidth is very good
#at synchronizing to continuous FSK signals
self._carriertrack = digital.fll_band_edge_cc(self._clockrec_sps,
0.6, #rolloff factor
64, #taps
1.0) #loop bandwidth
print "Samples per symbol: %f" % (self._clockrec_sps,)
self._softbits = digital.clock_recovery_mm_ff(self._clockrec_sps,
0.25*self._gain_mu*self._gain_mu, #gain omega, = mu/2 * mu_gain^2
self._mu, #mu (decision threshold)
self._gain_mu, #mu gain
self._omega_relative_limit) #omega relative limit
self._slicer = digital.binary_slicer_fb()
if self._decim > 1:
self.connect(self, self._downsample, self._carriertrack, self._demod, self._softbits, self._slicer, self)
else:
self.connect(self, self._carriertrack, self._demod, self._softbits, self._slicer, self)
def __init__(self, bt = 0.3, samples_per_symbol = 2, ti_adj=False):
gr.hier_block2.__init__(self, "msk_demod",
gr.io_signature(1, 1, gr.sizeof_char),
gr.io_signature(1, 1, gr.sizeof_gr_complex))
self.sps = 2
self.bt = 0.35
self.mu = 0.5
self.gain_mu = 0.175
self.freq_error = 0.0
self.omega_relative_limit = 0.005
self.omega = self.sps * (1 + self.freq_error)
self.gain_omega = .25 * self.gain_mu * self.gain_mu # critically damped
ntaps = 4 * samples_per_symbol # up to 3 bits in filter at once
sensitivity = (pi / 2) / samples_per_symbol # phase change per bit = pi / 2
# Turn it into NRZ data.
self.unpack = gr.packed_to_unpacked_bb(1, gr.GR_MSB_FIRST)
self.nrz = digital.chunks_to_symbols_bf([-1, 1], 1) # note could also invert bits here
# Form Gaussian filter
# Generate Gaussian response (Needs to be convolved with window below).
self.gaussian_taps = gr.firdes.gaussian(1, samples_per_symbol, bt, ntaps)
self.sqwave = (1,) * samples_per_symbol # rectangular window
self.taps = numpy.convolve(numpy.array(self.gaussian_taps),numpy.array(self.sqwave))
self.gaussian_filter = filter.interp_fir_filter_fff(samples_per_symbol, self.taps)
# the clock recovery block tracks the symbol clock and resamples as needed.
# the output of the block is a stream of soft symbols (float)
self.clock_recovery = digital.clock_recovery_mm_ff(self.omega, self.gain_omega,
self.mu, self.gain_mu,
self.omega_relative_limit)
# FM modulation
self.fmmod = gr.frequency_modulator_fc(sensitivity)
# TODO: this is hardcoded, how to figure out this value?
self.offset = gr.add_const_vff((-.5, ))
# CC430 RF core is inverted with respect to USRP for some reason
self.invert = gr.multiply_const_vff((-1, ))
# Connect & Initialize base class
if ti_adj:
self.connect(self, self.unpack, self.nrz, self.invert, self.offset, self.gaussian_filter, self.fmmod, self)
else:
self.connect(self, self.unpack, self.nrz, self.gaussian_filter, self.fmmod, self)
def __init__(self, input_filename="usrp.dat", output_filename="frames.dat", decim=32, pll_alpha=0.01, clock_alpha=0.01): gr.top_block.__init__(self, "Hrpt Demod") ################################################## # Parameters ################################################## self.input_filename = input_filename self.output_filename = output_filename self.decim = decim self.pll_alpha = pll_alpha self.clock_alpha = clock_alpha ################################################## # Variables ################################################## self.sym_rate = sym_rate = 600*1109 self.sample_rate = sample_rate = 4e6 self.sps = sps = sample_rate/sym_rate self.max_clock_offset = max_clock_offset = 100e-6 self.max_carrier_offset = max_carrier_offset = 2*math.pi*100e3/sample_rate self.hs = hs = int(sps/2.0) ################################################## # Blocks ################################################## self.pll = noaa.hrpt_pll_cf(pll_alpha, pll_alpha**2/4.0, max_carrier_offset) self.noaa_hrpt_deframer_0 = noaa.hrpt_deframer() self.gr_moving_average_xx_0 = gr.moving_average_ff(hs, 1.0/hs, 4000) self.gr_file_sink_0 = gr.file_sink(gr.sizeof_short*1, output_filename) self.gr_file_sink_0.set_unbuffered(False) self.file_source = gr.file_source(gr.sizeof_short*1, input_filename, False) self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(sps/2.0, clock_alpha**2/4.0, 0.5, clock_alpha, max_clock_offset) self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb() self.decoder = noaa.hrpt_decoder(True,False) self.cs2cf = gr.interleaved_short_to_complex() self.agc = gr.agc_cc(1e-5, 1.0, 1.0/32768.0, 1.0) ################################################## # Connections ################################################## self.connect((self.file_source, 0), (self.cs2cf, 0)) self.connect((self.pll, 0), (self.gr_moving_average_xx_0, 0)) self.connect((self.cs2cf, 0), (self.agc, 0)) self.connect((self.agc, 0), (self.pll, 0)) self.connect((self.noaa_hrpt_deframer_0, 0), (self.gr_file_sink_0, 0)) self.connect((self.noaa_hrpt_deframer_0, 0), (self.decoder, 0)) self.connect((self.gr_moving_average_xx_0, 0), (self.digital_clock_recovery_mm_xx_0, 0)) self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0)) self.connect((self.digital_binary_slicer_fb_0, 0), (self.noaa_hrpt_deframer_0, 0))
def __init__(self, input_rate, symbol_rate): super(clock_recovery, self).__init__( "clock_recovery", gr.io_signature(1, 1, gr.sizeof_float*1), gr.io_signature(1, 1, gr.sizeof_char*1), ) samples_per_symbol = float(input_rate) / symbol_rate omega_relative_limit = 0.02 gain_mu = 0.4 / samples_per_symbol self.clock_recovery = digital.clock_recovery_mm_ff(samples_per_symbol*(1+0.00), 0.25*gain_mu*gain_mu, 0.5, gain_mu, omega_relative_limit) self.slicer = digital.binary_slicer_fb() self.connect((self, 0), (self.clock_recovery, 0)) self.connect((self.clock_recovery, 0), (self.slicer, 0)) self.connect((self.slicer, 0), (self, 0))
def __init__(self):
gr.top_block.__init__(self, "Elster Rx Nogui")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 2400000
self.rx_gain = rx_gain = 45
self.corr = corr = 0
self.channel_rate = channel_rate = 400000
self.channel_decimation = channel_decimation = 4
self.center_freq = center_freq = 904600000
##################################################
# Blocks
##################################################
self.osmosdr_source_0 = osmosdr.source( args="numchan=" + str(1) + " " + "" )
self.osmosdr_source_0.set_sample_rate(samp_rate)
self.osmosdr_source_0.set_center_freq(center_freq, 0)
self.osmosdr_source_0.set_freq_corr(corr, 0)
self.osmosdr_source_0.set_dc_offset_mode(0, 0)
self.osmosdr_source_0.set_iq_balance_mode(0, 0)
self.osmosdr_source_0.set_gain_mode(0, 0)
self.osmosdr_source_0.set_gain(rx_gain, 0)
self.osmosdr_source_0.set_if_gain(20, 0)
self.osmosdr_source_0.set_bb_gain(20, 0)
self.osmosdr_source_0.set_antenna("", 0)
self.osmosdr_source_0.set_bandwidth(0, 0)
self.low_pass_filter_1 = filter.fir_filter_fff(channel_decimation, firdes.low_pass(
1, channel_rate, 20000, 5000, firdes.WIN_HAMMING, 6.76))
self.elster_packetize_0 = elster.packetize(1)
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(channel_rate * 56.48E-6 / 2 / channel_decimation, 0.25*(0.05*0.05), 0.5, 0.05, 0.005)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.blocks_keep_one_in_n_0 = blocks.keep_one_in_n(gr.sizeof_gr_complex*1, samp_rate/channel_rate)
self.analog_quadrature_demod_cf_0 = analog.quadrature_demod_cf(-channel_rate/(115000*2*3.1416))
##################################################
# Connections
##################################################
self.connect((self.osmosdr_source_0, 0), (self.blocks_keep_one_in_n_0, 0))
self.connect((self.blocks_keep_one_in_n_0, 0), (self.analog_quadrature_demod_cf_0, 0))
self.connect((self.analog_quadrature_demod_cf_0, 0), (self.low_pass_filter_1, 0))
self.connect((self.low_pass_filter_1, 0), (self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.elster_packetize_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0))
def __init__(self):
grc_wxgui.top_block_gui.__init__(self, title="Top Block")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 1000000
self.dec_rate = dec_rate = 2
##################################################
# Blocks
##################################################
self.uhd_usrp_source_0 = uhd.usrp_source(
device_addr="", stream_args=uhd.stream_args(cpu_format="fc32", channels=range(1))
)
self.uhd_usrp_source_0.set_samp_rate(1000000)
self.uhd_usrp_source_0.set_center_freq(2400490000, 0)
self.uhd_usrp_source_0.set_gain(0, 0)
self.low_pass_filter_0_0 = gr.interp_fir_filter_ccf(
1, firdes.low_pass(1, samp_rate, 500000, 10000, firdes.WIN_HAMMING, 6.76)
)
self.low_pass_filter_0 = gr.fir_filter_fff(
1, firdes.low_pass(1, samp_rate, 500000, 10000, firdes.WIN_HAMMING, 6.76)
)
self.gr_quadrature_demod_cf_0 = gr.quadrature_demod_cf(1)
self.gr_pwr_squelch_xx_0 = gr.pwr_squelch_cc(-100, 0.001, 0, True)
self.flysky_dumpsync_0 = flysky.dumpsync()
self.digital_correlate_access_code_bb_0_1 = digital.correlate_access_code_bb(
"010101010101010101010101010101010101010001110101", 1
)
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(2, 0.0076562, 0.5, 0.175, 0.005)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
##################################################
# Connections
##################################################
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0))
self.connect((self.gr_quadrature_demod_cf_0, 0), (self.low_pass_filter_0, 0))
self.connect((self.low_pass_filter_0, 0), (self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.gr_pwr_squelch_xx_0, 0), (self.gr_quadrature_demod_cf_0, 0))
self.connect((self.uhd_usrp_source_0, 0), (self.low_pass_filter_0_0, 0))
self.connect((self.low_pass_filter_0_0, 0), (self.gr_pwr_squelch_xx_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.digital_correlate_access_code_bb_0_1, 0))
self.connect((self.digital_correlate_access_code_bb_0_1, 0), (self.flysky_dumpsync_0, 0))
def __init__(self,
baudrate=1200,
client_address='127.0.0.1',
client_port=4336,
host_address='127.0.0.1',
host_port=4326,
samp_rate=48000):
gr.top_block.__init__(self, "Udp Decoder")
##################################################
# Parameters
##################################################
self.baudrate = baudrate
self.client_address = client_address
self.client_port = client_port
self.host_address = host_address
self.host_port = host_port
self.samp_rate = samp_rate
##################################################
# Blocks
##################################################
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(
samp_rate / baudrate, 0.001, 0, 0.25, 0.001)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.blocks_udp_source_0 = blocks.udp_source(gr.sizeof_short * 1,
host_address, host_port,
1472, True)
self.blocks_udp_sink_0 = blocks.udp_sink(gr.sizeof_char * 1,
client_address, client_port,
1472, True)
self.blocks_short_to_float_0 = blocks.short_to_float(1, 32767)
##################################################
# Connections
##################################################
self.connect((self.blocks_short_to_float_0, 0),
(self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.blocks_udp_source_0, 0),
(self.blocks_short_to_float_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0),
(self.blocks_udp_sink_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0),
(self.digital_binary_slicer_fb_0, 0))
def __init__(self, input_rate, symbol_rate):
super(clock_recovery, self).__init__(
"clock_recovery",
gr.io_signature(1, 1, gr.sizeof_float * 1),
gr.io_signature(1, 1, gr.sizeof_char * 1),
)
samples_per_symbol = float(input_rate) / symbol_rate
omega_relative_limit = 0.02
gain_mu = 0.4 / samples_per_symbol
self.clock_recovery = digital.clock_recovery_mm_ff(
samples_per_symbol * (1 + 0.00), 0.25 * gain_mu * gain_mu, 0.5,
gain_mu, omega_relative_limit)
self.slicer = digital.binary_slicer_fb()
self.connect((self, 0), (self.clock_recovery, 0))
self.connect((self.clock_recovery, 0), (self.slicer, 0))
self.connect((self.slicer, 0), (self, 0))
def _update_filtered(self, translated_burst):
if self._rrc:
translated_data_1 = translate_burst(translated_burst, self._deviation)
filtered_data_1 = TimeData(scipy.signal.lfilter(self._taps, 1, translated_data_1.samples), translated_data_1.sampling_rate)
translated_data_2 = translate_burst(translated_burst, -self._deviation)
filtered_data_2 = TimeData(scipy.signal.lfilter(self._taps, 1, translated_data_2.samples), translated_data_2.sampling_rate)
self.eye_view.data = (filtered_data_1.abs, filtered_data_2.abs)
else:
filtered_data_1 = TimeData(numpy.complex64(scipy.signal.lfilter(self._taps_n, 1, translated_burst.samples)), translated_burst.sampling_rate)
filtered_data_2 = TimeData(numpy.complex64(scipy.signal.lfilter(self._taps_p, 1, translated_burst.samples)), translated_burst.sampling_rate)
self.eye_view.data = (filtered_data_1.abs, filtered_data_2.abs)
filtered_diff = TimeData(filtered_data_2.abs.samples - filtered_data_1.abs.samples, filtered_data_1.sampling_rate)
self.slicer_view.data = filtered_diff
#print('sliced abs sum: %s' % sum(abs(filtered_diff.samples)))
omega = self._samples_per_symbol
mu = 0.5
data_source = filtered_diff.samples
numpy_source = NumpySource(data_source)
clock_recovery = digital.clock_recovery_mm_ff(omega, self._gain_omega, mu, self._gain_mu, self._omega_relative_limit)
#clock_recovery = digital.pfb_clock_sync_fff(self._samples_per_symbol, 1.0, self._taps)
numpy_sink = NumpySink(numpy.float32)
top = gr.top_block()
top.connect(numpy_source, clock_recovery)
top.connect(clock_recovery, numpy_sink)
top.run()
data_sink = numpy_sink.data
# TODO: Adjust sampling rate
bits = []
for i in range(len(data_sink)):
if data_sink[i] >= 0:
data_sink[i] = 1
bits.append('1')
else:
data_sink[i] = -1
bits.append('0')
bits = ''.join(bits)
#print(bits)
self.sliced_view.data = TimeData(data_sink, filtered_diff.sampling_rate)
def __init__(self, samp_rate=1600000, samp_per_sym=16, freq_error=-0.0025000): gr.hier_block2.__init__( self, "Wireless M-Bus Demod", gr.io_signature(1, 1, gr.sizeof_gr_complex*1), gr.io_signature(1, 1, gr.sizeof_char*1), ) ################################################## # Parameters ################################################## self.samp_rate = samp_rate self.samp_per_sym = samp_per_sym self.freq_error = freq_error ################################################## # Variables ################################################## self.cutoff = cutoff = 120e3 self.chip_rate = chip_rate = samp_rate/samp_per_sym ################################################## # Blocks ################################################## self.low_pass_filter_0 = gnuradio.filter.fir_filter_ccf(1, gnuradio.filter.firdes.low_pass( 1, samp_rate, cutoff, cutoff/2, gnuradio.filter.firdes.WIN_HAMMING, 6.76)) self.gr_sub_xx_0 = blocks.sub_ff(1) self.gr_single_pole_iir_filter_xx_0 = gnuradio.filter.single_pole_iir_filter_ff(0.0512/samp_per_sym, 1) self.gr_quadrature_demod_cf_0 = analog.quadrature_demod_cf(1) self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(samp_per_sym*(1+freq_error), .25 *0.06*0.06*4, 0.5, 0.06*2, 0.002*2) self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb() ################################################## # Connections ################################################## self.connect((self.gr_quadrature_demod_cf_0, 0), (self.gr_single_pole_iir_filter_xx_0, 0)) self.connect((self.gr_single_pole_iir_filter_xx_0, 0), (self.gr_sub_xx_0, 1)) self.connect((self.gr_quadrature_demod_cf_0, 0), (self.gr_sub_xx_0, 0)) self.connect((self.gr_sub_xx_0, 0), (self.digital_clock_recovery_mm_xx_0, 0)) self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0)) self.connect((self.low_pass_filter_0, 0), (self.gr_quadrature_demod_cf_0, 0)) self.connect((self.digital_binary_slicer_fb_0, 0), (self, 0)) self.connect((self, 0), (self.low_pass_filter_0, 0))
def __init__(self):
"""
Hierarchical block for FSK demodulation.
The input is the complex modulated signal at baseband
and the output is a stream of floats.
"""
# Initialize base class
gr.hier_block2.__init__(self, "fsk_demod",
gr.io_signature(1, 1, gr.sizeof_gr_complex),
gr.io_signature(1, 1, gr.sizeof_float))
# Variables
self.sps = sps = 2
self.sensitivity = sensitivity = (pi / 2) / sps
self.alpha = alpha = 0.0512 / sps
self.fm_demod = gr.quadrature_demod_cf(1 / sensitivity)
self.freq_offset = gr.single_pole_iir_filter_ff(alpha)
self.sub = gr.sub_ff()
self.add = gr.add_ff()
self.multiply = gr.multiply_ff()
self.invert = gr.multiply_const_vff((-1, ))
# recover the clock
omega = sps
gain_mu = 0.03
mu = 0.5
omega_relative_limit = 0.0002
freq_error = 0.0
gain_omega = .25 * gain_mu * gain_mu # critically damped
self.clock_recovery = digital.clock_recovery_mm_ff(
omega, gain_omega, mu, gain_mu, omega_relative_limit)
self.slice = digital.binary_slicer_fb()
# Connections
self.connect(self.fm_demod, (self.add, 0))
self.connect(self.fm_demod, self.freq_offset, (self.add, 1))
self.connect(self, self.fm_demod)
self.connect(self.add, self.clock_recovery, self.invert, self)
def __init__(self):
"""
Hierarchical block for FSK demodulation.
The input is the complex modulated signal at baseband
and the output is a stream of floats.
"""
# Initialize base class
gr.hier_block2.__init__(
self, "fsk_demod", gr.io_signature(1, 1, gr.sizeof_gr_complex), gr.io_signature(1, 1, gr.sizeof_float)
)
# Variables
self.sps = sps = 2
self.sensitivity = sensitivity = (pi / 2) / sps
self.alpha = alpha = 0.0512 / sps
self.fm_demod = gr.quadrature_demod_cf(1 / sensitivity)
self.freq_offset = gr.single_pole_iir_filter_ff(alpha)
self.sub = gr.sub_ff()
self.add = gr.add_ff()
self.multiply = gr.multiply_ff()
self.invert = gr.multiply_const_vff((-1,))
# recover the clock
omega = sps
gain_mu = 0.03
mu = 0.5
omega_relative_limit = 0.0002
freq_error = 0.0
gain_omega = 0.25 * gain_mu * gain_mu # critically damped
self.clock_recovery = digital.clock_recovery_mm_ff(omega, gain_omega, mu, gain_mu, omega_relative_limit)
self.slice = digital.binary_slicer_fb()
# Connections
self.connect(self.fm_demod, (self.add, 0))
self.connect(self.fm_demod, self.freq_offset, (self.add, 1))
self.connect(self, self.fm_demod)
self.connect(self.add, self.clock_recovery, self.invert, self)
def __init__(self, dc_block_len=1000, samp_per_sym=10, rx_bit_inv=1):
gr.hier_block2.__init__(
self,
"GMSK Demodulator",
gr.io_signature(1, 1, gr.sizeof_float * 1),
gr.io_signature(1, 1, gr.sizeof_char * 1),
)
##################################################
# Parameters
##################################################
self.dc_block_len = dc_block_len
self.samp_per_sym = samp_per_sym
self.rx_bit_inv = rx_bit_inv
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 48000
##################################################
# Blocks
##################################################
self.digital_clock_recovery_mm_xx_0_0 = digital.clock_recovery_mm_ff(
samp_per_sym * (1 + 0.0), 0.25 * 0.175 * 0.175, 0.5, 0.175, 0.005)
self.digital_binary_slicer_fb_0_0 = digital.binary_slicer_fb()
self.dc_blocker_xx_0 = filter.dc_blocker_ff(dc_block_len, False)
self.blocks_multiply_const_vxx_0_0 = blocks.multiply_const_vff(
(rx_bit_inv, ))
##################################################
# Connections
##################################################
self.connect((self.dc_blocker_xx_0, 0),
(self.blocks_multiply_const_vxx_0_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0_0, 0),
(self.digital_binary_slicer_fb_0_0, 0))
self.connect((self.blocks_multiply_const_vxx_0_0, 0),
(self.digital_clock_recovery_mm_xx_0_0, 0))
self.connect((self.digital_binary_slicer_fb_0_0, 0), (self, 0))
self.connect((self, 0), (self.dc_blocker_xx_0, 0))
def __init__(self):
gr.top_block.__init__(self, "TW-1 test decoder")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 48000
self.samp_per_sym = samp_per_sym = 10
self.gain_mu = gain_mu = 0.175*3
##################################################
# Blocks
##################################################
self.synctags_fixedlen_tagger_0 = synctags.fixedlen_tagger("syncword", "packet_len", 256*8, numpy.byte)
self.low_pass_filter_0 = filter.fir_filter_fff(1, firdes.low_pass(
1, samp_rate, 2400, 2000, firdes.WIN_HAMMING, 6.76))
self.digital_descrambler_bb_0 = digital.descrambler_bb(0x21, 0x00, 16)
self.digital_correlate_access_code_tag_bb_0 = digital.correlate_access_code_tag_bb("10010011000010110101000111011110", 4, "syncword")
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(samp_per_sym*(1+0.0), 0.25*gain_mu*gain_mu, 0.5, gain_mu, 0.005)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.blocks_wavfile_source_0 = blocks.wavfile_source("/tmp/tw-1c.wav", False)
self.blocks_unpacked_to_packed_xx_0 = blocks.unpacked_to_packed_bb(1, gr.GR_MSB_FIRST)
self.blocks_tagged_stream_to_pdu_0 = blocks.tagged_stream_to_pdu(blocks.byte_t, "packet_len")
self.blocks_tagged_stream_multiply_length_0 = blocks.tagged_stream_multiply_length(gr.sizeof_char*1, "packet_len", 1/8.0)
self.blocks_message_debug_1 = blocks.message_debug()
self.ax100_gomx3_rs_decode_0 = ax100.gomx3_rs_decode(True)
##################################################
# Connections
##################################################
self.msg_connect((self.ax100_gomx3_rs_decode_0, 'out'), (self.blocks_message_debug_1, 'print_pdu'))
self.msg_connect((self.blocks_tagged_stream_to_pdu_0, 'pdus'), (self.ax100_gomx3_rs_decode_0, 'in'))
self.connect((self.blocks_tagged_stream_multiply_length_0, 0), (self.blocks_tagged_stream_to_pdu_0, 0))
self.connect((self.blocks_unpacked_to_packed_xx_0, 0), (self.blocks_tagged_stream_multiply_length_0, 0))
self.connect((self.blocks_wavfile_source_0, 0), (self.low_pass_filter_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.digital_descrambler_bb_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0))
self.connect((self.digital_correlate_access_code_tag_bb_0, 0), (self.synctags_fixedlen_tagger_0, 0))
self.connect((self.digital_descrambler_bb_0, 0), (self.digital_correlate_access_code_tag_bb_0, 0))
self.connect((self.low_pass_filter_0, 0), (self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.synctags_fixedlen_tagger_0, 0), (self.blocks_unpacked_to_packed_xx_0, 0))
def __init__(self):
"""
Hierarchical block for MSK demodulation.
The input is the complex modulated signal at baseband
and the output is a stream of floats.
"""
# Initialize base class
gr.hier_block2.__init__(self, "msk_demod",
gr.io_signature(1, 1, gr.sizeof_gr_complex),
gr.io_signature(1, 1, gr.sizeof_float))
self.sps = 2
self.bt = 0.35
self.mu = 0.5
self.gain_mu = 0.175
self.freq_error = 0.0
self.omega_relative_limit = 0.005
self.omega = self.sps * (1 + self.freq_error)
self.gain_omega = .25 * self.gain_mu * self.gain_mu # critically damped
# Demodulate FM
sensitivity = (pi / 2) / self.sps
self.fmdemod = gr.quadrature_demod_cf(1.0 / sensitivity)
self.invert = gr.multiply_const_vff((-1, ))
# TODO: this is hardcoded, how to figure out this value?
self.offset = gr.add_const_vff((-1.4, ))
# the clock recovery block tracks the symbol clock and resamples as needed.
# the output of the block is a stream of soft symbols (float)
self.clock_recovery = digital.clock_recovery_mm_ff(
self.omega, self.gain_omega, self.mu, self.gain_mu,
self.omega_relative_limit)
self.slicer = digital.binary_slicer_fb()
self.connect(self, self.fmdemod, self.invert, self.clock_recovery,
self.offset, self)
def __init__(self):
"""
Hierarchical block for MSK demodulation.
The input is the complex modulated signal at baseband
and the output is a stream of floats.
"""
# Initialize base class
gr.hier_block2.__init__(
self, "msk_demod", gr.io_signature(1, 1, gr.sizeof_gr_complex), gr.io_signature(1, 1, gr.sizeof_float)
)
self.sps = 2
self.bt = 0.35
self.mu = 0.5
self.gain_mu = 0.175
self.freq_error = 0.0
self.omega_relative_limit = 0.005
self.omega = self.sps * (1 + self.freq_error)
self.gain_omega = 0.25 * self.gain_mu * self.gain_mu # critically damped
# Demodulate FM
sensitivity = (pi / 2) / self.sps
self.fmdemod = gr.quadrature_demod_cf(1.0 / sensitivity)
self.invert = gr.multiply_const_vff((-1,))
# TODO: this is hardcoded, how to figure out this value?
self.offset = gr.add_const_vff((-1.2,))
# the clock recovery block tracks the symbol clock and resamples as needed.
# the output of the block is a stream of soft symbols (float)
self.clock_recovery = digital.clock_recovery_mm_ff(
self.omega, self.gain_omega, self.mu, self.gain_mu, self.omega_relative_limit
)
self.slicer = digital.binary_slicer_fb()
self.connect(self, self.fmdemod, self.invert, self.clock_recovery, self.offset, self)
def __init__(self,
samplerate,
symbolrate=SYMRATE,
channel_str=None,
sendmsg=True,
debug=False,
samplepersymbol=SPS,
fmdeviation=FM_DEVIATION):
gr.hier_block2.__init__(self, "pocsag",
gr.io_signature(1, 1, gr.sizeof_gr_complex),
gr.io_signature(1, 1, 1))
self.samplerate = samplerate
self.symbolrate = symbolrate
self.sendmsg = sendmsg
self.debug = debug
self.samplepersymbol = samplepersymbol
self.fmdeviation = fmdeviation
self.fractional_interpolator = gr.fractional_interpolator_cc(
0, 1.0 * samplerate / (symbolrate * samplepersymbol))
self.quadrature_demod = gr.quadrature_demod_cf(
(symbolrate * samplepersymbol) / (fmdeviation * 4.0))
self.low_pass_filter = gr.fir_filter_fff(
1,
gr.firdes.low_pass(1, symbolrate * samplepersymbol, symbolrate * 2,
symbolrate / 2.0, gr.firdes.WIN_HAMMING, 6.76))
self.digital_clock_recovery_mm = digital.clock_recovery_mm_ff(
samplepersymbol, 0.03 * 0.03 * 0.3, 0.4, 0.03, 1e-4)
self.digital_binary_slicer_fb = digital.binary_slicer_fb()
self.pktdecoder = pocsag_pktdecoder(channel_str=channel_str,
sendmsg=sendmsg,
debug=debug)
self.connect(self, self.fractional_interpolator, self.quadrature_demod,
self.low_pass_filter, self.digital_clock_recovery_mm,
self.digital_binary_slicer_fb, self.pktdecoder, self)
def __init__(self):
gr.top_block.__init__(self, "Top Block")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 1024000
self.decimation = decimation = 10
self.symbol_rate = symbol_rate = 38400
self.input_rate = input_rate = samp_rate / decimation
self.window_symbols = window_symbols = 1
self.symbol_taps_length = symbol_taps_length = int(
(float(input_rate) / symbol_rate))
self.samples_per_symbol = samples_per_symbol = float(
input_rate) / symbol_rate
self.grab_freq = grab_freq = 868200000
self.symbol_taps = symbol_taps = (1, ) * symbol_taps_length
self.offset_sign2 = offset_sign2 = (-34780 * grab_freq) / 868200000
self.offset_sign1 = offset_sign1 = (11600 * grab_freq) / 868200000
self.gain_mu = gain_mu = 0.4 / samples_per_symbol
self.fsk_deviation_hz = fsk_deviation_hz = 5200
self.average_window = average_window = int(
(input_rate * window_symbols / symbol_rate))
##################################################
# Blocks
##################################################
self.urmetEasyRead_urmetEasyRead_0 = urmetEasyRead.urmetEasyRead()
self.digital_correlate_access_code_bb_0 = digital.correlate_access_code_bb(
'11010011100100011101001110010001', 1)
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(
samples_per_symbol, 0.25 * gain_mu * gain_mu, 0.5, gain_mu, 0.02)
self.digital_binary_slicer_fb_0_0 = digital.binary_slicer_fb()
self.cc1111_cc1111_packet_decoder_0 = cc1111.cc1111_packet_decoder(
gr.msg_queue(1), False, True, True, True)
self.blocks_uchar_to_float_0 = blocks.uchar_to_float()
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_gr_complex * 1,
samp_rate, True)
self.blocks_null_sink_1 = blocks.null_sink(gr.sizeof_char * 1)
self.blocks_float_to_complex_0 = blocks.float_to_complex(1)
self.blocks_file_source_0 = blocks.file_source(gr.sizeof_char * 1,
sys.argv[1], False)
self.blocks_deinterleave_0 = blocks.deinterleave(
gr.sizeof_float * 1, 1)
self.blocks_add_const_vxx_0 = blocks.add_const_vcc((-128 - 128j, ))
self.analog_simple_squelch_cc_0 = analog.simple_squelch_cc(40, 1)
self.analog_quadrature_demod_cf_0 = analog.quadrature_demod_cf(1)
self.FXFIR1 = filter.freq_xlating_fir_filter_ccc(
decimation, (1, ), -14e3, samp_rate)
##################################################
# Connections
##################################################
self.connect((self.FXFIR1, 0), (self.analog_simple_squelch_cc_0, 0))
self.connect((self.analog_quadrature_demod_cf_0, 0),
(self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.analog_simple_squelch_cc_0, 0),
(self.analog_quadrature_demod_cf_0, 0))
self.connect((self.blocks_add_const_vxx_0, 0),
(self.blocks_throttle_0, 0))
self.connect((self.blocks_deinterleave_0, 0),
(self.blocks_float_to_complex_0, 0))
self.connect((self.blocks_deinterleave_0, 1),
(self.blocks_float_to_complex_0, 1))
self.connect((self.blocks_file_source_0, 0),
(self.blocks_uchar_to_float_0, 0))
self.connect((self.blocks_float_to_complex_0, 0),
(self.blocks_add_const_vxx_0, 0))
self.connect((self.blocks_throttle_0, 0), (self.FXFIR1, 0))
self.connect((self.blocks_uchar_to_float_0, 0),
(self.blocks_deinterleave_0, 0))
self.connect((self.cc1111_cc1111_packet_decoder_0, 0),
(self.urmetEasyRead_urmetEasyRead_0, 0))
self.connect((self.digital_binary_slicer_fb_0_0, 0),
(self.digital_correlate_access_code_bb_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0),
(self.digital_binary_slicer_fb_0_0, 0))
self.connect((self.digital_correlate_access_code_bb_0, 0),
(self.cc1111_cc1111_packet_decoder_0, 0))
self.connect((self.urmetEasyRead_urmetEasyRead_0, 0),
(self.blocks_null_sink_1, 0))
def __init__(self):
gr.top_block.__init__(self, "Uhd Ais 3")
Qt.QWidget.__init__(self)
self.setWindowTitle("Uhd Ais 3")
qtgui.util.check_set_qss()
try:
self.setWindowIcon(Qt.QIcon.fromTheme('gnuradio-grc'))
except:
pass
self.top_scroll_layout = Qt.QVBoxLayout()
self.setLayout(self.top_scroll_layout)
self.top_scroll = Qt.QScrollArea()
self.top_scroll.setFrameStyle(Qt.QFrame.NoFrame)
self.top_scroll_layout.addWidget(self.top_scroll)
self.top_scroll.setWidgetResizable(True)
self.top_widget = Qt.QWidget()
self.top_scroll.setWidget(self.top_widget)
self.top_layout = Qt.QVBoxLayout(self.top_widget)
self.top_grid_layout = Qt.QGridLayout()
self.top_layout.addLayout(self.top_grid_layout)
self.settings = Qt.QSettings("GNU Radio", "uhd_ais_3")
self.restoreGeometry(self.settings.value("geometry").toByteArray())
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 250e3
self.decim = decim = 5
self.baud = baud = 9600
self.samp_per_sym = samp_per_sym = (samp_rate / decim / 50 * 48) / baud
self.rx_gain = rx_gain = 45
self.fsk_deviation = fsk_deviation = 10e3
self.freq = freq = 162e6
##################################################
# Blocks
##################################################
self._rx_gain_tool_bar = Qt.QToolBar(self)
self._rx_gain_tool_bar.addWidget(Qt.QLabel("rx_gain" + ": "))
self._rx_gain_line_edit = Qt.QLineEdit(str(self.rx_gain))
self._rx_gain_tool_bar.addWidget(self._rx_gain_line_edit)
self._rx_gain_line_edit.returnPressed.connect(lambda: self.set_rx_gain(
eng_notation.str_to_num(
str(self._rx_gain_line_edit.text().toAscii()))))
self.top_grid_layout.addWidget(self._rx_gain_tool_bar, 8, 0, 1, 2)
self.rational_resampler_xxx_0_0 = filter.rational_resampler_fff(
interpolation=48,
decimation=50,
taps=None,
fractional_bw=None,
)
self.rational_resampler_xxx_0 = filter.rational_resampler_fff(
interpolation=48,
decimation=50,
taps=None,
fractional_bw=None,
)
self.qtgui_waterfall_sink_x_0_0 = qtgui.waterfall_sink_c(
1024, #size
firdes.WIN_BLACKMAN_hARRIS, #wintype
0, #fc
samp_rate, #bw
"AIS-B", #name
1 #number of inputs
)
self.qtgui_waterfall_sink_x_0_0.set_update_time(0.010)
self.qtgui_waterfall_sink_x_0_0.enable_grid(True)
self.qtgui_waterfall_sink_x_0_0.enable_axis_labels(True)
if not True:
self.qtgui_waterfall_sink_x_0_0.disable_legend()
if "complex" == "float" or "complex" == "msg_float":
self.qtgui_waterfall_sink_x_0_0.set_plot_pos_half(not True)
labels = ['', '', '', '', '', '', '', '', '', '']
colors = [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_waterfall_sink_x_0_0.set_line_label(
i, "Data {0}".format(i))
else:
self.qtgui_waterfall_sink_x_0_0.set_line_label(i, labels[i])
self.qtgui_waterfall_sink_x_0_0.set_color_map(i, colors[i])
self.qtgui_waterfall_sink_x_0_0.set_line_alpha(i, alphas[i])
self.qtgui_waterfall_sink_x_0_0.set_intensity_range(-60, 10)
self._qtgui_waterfall_sink_x_0_0_win = sip.wrapinstance(
self.qtgui_waterfall_sink_x_0_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_waterfall_sink_x_0_0_win, 2,
4, 2, 4)
self.qtgui_waterfall_sink_x_0 = qtgui.waterfall_sink_c(
1024, #size
firdes.WIN_BLACKMAN_hARRIS, #wintype
0, #fc
samp_rate, #bw
"AIS-A", #name
1 #number of inputs
)
self.qtgui_waterfall_sink_x_0.set_update_time(0.010)
self.qtgui_waterfall_sink_x_0.enable_grid(True)
self.qtgui_waterfall_sink_x_0.enable_axis_labels(True)
if not True:
self.qtgui_waterfall_sink_x_0.disable_legend()
if "complex" == "float" or "complex" == "msg_float":
self.qtgui_waterfall_sink_x_0.set_plot_pos_half(not True)
labels = ['', '', '', '', '', '', '', '', '', '']
colors = [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_waterfall_sink_x_0.set_line_label(
i, "Data {0}".format(i))
else:
self.qtgui_waterfall_sink_x_0.set_line_label(i, labels[i])
self.qtgui_waterfall_sink_x_0.set_color_map(i, colors[i])
self.qtgui_waterfall_sink_x_0.set_line_alpha(i, alphas[i])
self.qtgui_waterfall_sink_x_0.set_intensity_range(-60, 10)
self._qtgui_waterfall_sink_x_0_win = sip.wrapinstance(
self.qtgui_waterfall_sink_x_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_waterfall_sink_x_0_win, 0,
4, 2, 4)
self.qtgui_freq_sink_x_0 = qtgui.freq_sink_c(
2048, #size
firdes.WIN_BLACKMAN_hARRIS, #wintype
0, #fc
samp_rate / decim, #bw
"AIS", #name
2 #number of inputs
)
self.qtgui_freq_sink_x_0.set_update_time(0.01)
self.qtgui_freq_sink_x_0.set_y_axis(-60, 10)
self.qtgui_freq_sink_x_0.set_y_label('Relative Gain', 'dB')
self.qtgui_freq_sink_x_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, 0.0, 0,
"")
self.qtgui_freq_sink_x_0.enable_autoscale(False)
self.qtgui_freq_sink_x_0.enable_grid(True)
self.qtgui_freq_sink_x_0.set_fft_average(1.0)
self.qtgui_freq_sink_x_0.enable_axis_labels(True)
self.qtgui_freq_sink_x_0.enable_control_panel(False)
if not True:
self.qtgui_freq_sink_x_0.disable_legend()
if "complex" == "float" or "complex" == "msg_float":
self.qtgui_freq_sink_x_0.set_plot_pos_half(not True)
labels = ['', '', '', '', '', '', '', '', '', '']
widths = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
colors = [
"green", "red", "green", "black", "cyan", "magenta", "yellow",
"dark red", "dark green", "dark blue"
]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(2):
if len(labels[i]) == 0:
self.qtgui_freq_sink_x_0.set_line_label(
i, "Data {0}".format(i))
else:
self.qtgui_freq_sink_x_0.set_line_label(i, labels[i])
self.qtgui_freq_sink_x_0.set_line_width(i, widths[i])
self.qtgui_freq_sink_x_0.set_line_color(i, colors[i])
self.qtgui_freq_sink_x_0.set_line_alpha(i, alphas[i])
self._qtgui_freq_sink_x_0_win = sip.wrapinstance(
self.qtgui_freq_sink_x_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_freq_sink_x_0_win, 0, 0, 4,
4)
self.pyqt_text_output_0_0 = pyqt.text_output()
self._pyqt_text_output_0_0_win = self.pyqt_text_output_0_0
self.top_grid_layout.addWidget(self._pyqt_text_output_0_0_win, 4, 0, 4,
4)
self.pyqt_text_output_0 = pyqt.text_output()
self._pyqt_text_output_0_win = self.pyqt_text_output_0
self.top_grid_layout.addWidget(self._pyqt_text_output_0_win, 4, 4, 4,
4)
self.low_pass_filter_0_0 = filter.fir_filter_ccf(
decim,
firdes.low_pass(1, samp_rate, 7e3, 1e3, firdes.WIN_HAMMING, 6.76))
self.low_pass_filter_0 = filter.fir_filter_ccf(
decim,
firdes.low_pass(1, samp_rate, 7e3, 1e3, firdes.WIN_HAMMING, 6.76))
self.digital_hdlc_deframer_bp_0_0 = digital.hdlc_deframer_bp(11, 1000)
self.digital_hdlc_deframer_bp_0 = digital.hdlc_deframer_bp(11, 1000)
self.digital_diff_decoder_bb_0_0 = digital.diff_decoder_bb(2)
self.digital_diff_decoder_bb_0 = digital.diff_decoder_bb(2)
self.digital_clock_recovery_mm_xx_0_0 = digital.clock_recovery_mm_ff(
samp_per_sym * (1 + 0.0), 0.25 * 0.175 * 0.175, 0.5, 0.175, 0.005)
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(
samp_per_sym * (1 + 0.0), 0.25 * 0.175 * 0.175, 0.5, 0.175, 0.005)
self.digital_binary_slicer_fb_0_0 = digital.binary_slicer_fb()
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_gr_complex * 1,
samp_rate * 4, True)
self.blocks_socket_pdu_0 = blocks.socket_pdu("TCP_SERVER", '', '52001',
10000, False)
self.blocks_multiply_xx_1 = blocks.multiply_vcc(1)
self.blocks_multiply_xx_0 = blocks.multiply_vcc(1)
self.blocks_message_debug_0_1 = blocks.message_debug()
self.blocks_file_source_0 = blocks.file_source(
gr.sizeof_gr_complex * 1,
'/home/zleffke/workspace/captures/ais/ais_20161218_250k_2.32fc',
True)
self.analog_sig_source_x_1 = analog.sig_source_c(
samp_rate, analog.GR_COS_WAVE, 25e3 + 400, 1, 0)
self.analog_sig_source_x_0 = analog.sig_source_c(
samp_rate, analog.GR_COS_WAVE, -25e3 + 400, 1, 0)
self.analog_quadrature_demod_cf_0_0 = analog.quadrature_demod_cf(
(samp_rate / decim) / (2 * math.pi * fsk_deviation / 8.0))
self.analog_quadrature_demod_cf_0 = analog.quadrature_demod_cf(
(samp_rate / decim) / (2 * math.pi * fsk_deviation / 8.0))
self.analog_agc2_xx_0_0 = analog.agc2_cc(1e-3, 1e-1, 1.0, 1.0)
self.analog_agc2_xx_0_0.set_max_gain(65536)
self.analog_agc2_xx_0 = analog.agc2_cc(1e-3, 1e-1, 1.0, 1.0)
self.analog_agc2_xx_0.set_max_gain(65536)
self.ais_pdu_to_nmea_0_0 = ais.pdu_to_nmea('B')
self.ais_pdu_to_nmea_0 = ais.pdu_to_nmea('A')
self.ais_invert_0_0 = ais.invert()
self.ais_invert_0 = ais.invert()
##################################################
# Connections
##################################################
self.msg_connect((self.ais_pdu_to_nmea_0, 'out'),
(self.blocks_message_debug_0_1, 'print'))
self.msg_connect((self.ais_pdu_to_nmea_0, 'out'),
(self.blocks_message_debug_0_1, 'print_pdu'))
self.msg_connect((self.ais_pdu_to_nmea_0, 'out'),
(self.blocks_socket_pdu_0, 'pdus'))
self.msg_connect((self.ais_pdu_to_nmea_0, 'out'),
(self.pyqt_text_output_0_0, 'pdus'))
self.msg_connect((self.ais_pdu_to_nmea_0_0, 'out'),
(self.blocks_message_debug_0_1, 'print'))
self.msg_connect((self.ais_pdu_to_nmea_0_0, 'out'),
(self.blocks_message_debug_0_1, 'print_pdu'))
self.msg_connect((self.ais_pdu_to_nmea_0_0, 'out'),
(self.blocks_socket_pdu_0, 'pdus'))
self.msg_connect((self.ais_pdu_to_nmea_0_0, 'out'),
(self.pyqt_text_output_0, 'pdus'))
self.msg_connect((self.digital_hdlc_deframer_bp_0, 'out'),
(self.ais_pdu_to_nmea_0, 'to_nmea'))
self.msg_connect((self.digital_hdlc_deframer_bp_0_0, 'out'),
(self.ais_pdu_to_nmea_0_0, 'to_nmea'))
self.connect((self.ais_invert_0, 0),
(self.digital_hdlc_deframer_bp_0, 0))
self.connect((self.ais_invert_0_0, 0),
(self.digital_hdlc_deframer_bp_0_0, 0))
self.connect((self.analog_agc2_xx_0, 0),
(self.analog_quadrature_demod_cf_0_0, 0))
self.connect((self.analog_agc2_xx_0, 0), (self.qtgui_freq_sink_x_0, 0))
self.connect((self.analog_agc2_xx_0, 0),
(self.qtgui_waterfall_sink_x_0, 0))
self.connect((self.analog_agc2_xx_0_0, 0),
(self.analog_quadrature_demod_cf_0, 0))
self.connect((self.analog_agc2_xx_0_0, 0),
(self.qtgui_freq_sink_x_0, 1))
self.connect((self.analog_agc2_xx_0_0, 0),
(self.qtgui_waterfall_sink_x_0_0, 0))
self.connect((self.analog_quadrature_demod_cf_0, 0),
(self.rational_resampler_xxx_0, 0))
self.connect((self.analog_quadrature_demod_cf_0_0, 0),
(self.rational_resampler_xxx_0_0, 0))
self.connect((self.analog_sig_source_x_0, 0),
(self.blocks_multiply_xx_0, 1))
self.connect((self.analog_sig_source_x_1, 0),
(self.blocks_multiply_xx_1, 0))
self.connect((self.blocks_file_source_0, 0),
(self.blocks_throttle_0, 0))
self.connect((self.blocks_multiply_xx_0, 0),
(self.low_pass_filter_0_0, 0))
self.connect((self.blocks_multiply_xx_1, 0),
(self.low_pass_filter_0, 0))
self.connect((self.blocks_throttle_0, 0),
(self.blocks_multiply_xx_0, 0))
self.connect((self.blocks_throttle_0, 0),
(self.blocks_multiply_xx_1, 1))
self.connect((self.digital_binary_slicer_fb_0, 0),
(self.digital_diff_decoder_bb_0, 0))
self.connect((self.digital_binary_slicer_fb_0_0, 0),
(self.digital_diff_decoder_bb_0_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0),
(self.digital_binary_slicer_fb_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0_0, 0),
(self.digital_binary_slicer_fb_0_0, 0))
self.connect((self.digital_diff_decoder_bb_0, 0),
(self.ais_invert_0, 0))
self.connect((self.digital_diff_decoder_bb_0_0, 0),
(self.ais_invert_0_0, 0))
self.connect((self.low_pass_filter_0, 0), (self.analog_agc2_xx_0, 0))
self.connect((self.low_pass_filter_0_0, 0),
(self.analog_agc2_xx_0_0, 0))
self.connect((self.rational_resampler_xxx_0, 0),
(self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.rational_resampler_xxx_0_0, 0),
(self.digital_clock_recovery_mm_xx_0_0, 0))
def __init__(self):
gr.top_block.__init__(self, "OPS-SAT UHF demodulator/decoder")
Qt.QWidget.__init__(self)
self.setWindowTitle("OPS-SAT UHF demodulator/decoder")
qtgui.util.check_set_qss()
try:
self.setWindowIcon(Qt.QIcon.fromTheme('gnuradio-grc'))
except:
pass
self.top_scroll_layout = Qt.QVBoxLayout()
self.setLayout(self.top_scroll_layout)
self.top_scroll = Qt.QScrollArea()
self.top_scroll.setFrameStyle(Qt.QFrame.NoFrame)
self.top_scroll_layout.addWidget(self.top_scroll)
self.top_scroll.setWidgetResizable(True)
self.top_widget = Qt.QWidget()
self.top_scroll.setWidget(self.top_widget)
self.top_layout = Qt.QVBoxLayout(self.top_widget)
self.top_grid_layout = Qt.QGridLayout()
self.top_layout.addLayout(self.top_grid_layout)
self.settings = Qt.QSettings("GNU Radio", "os_demod_decode")
self.restoreGeometry(self.settings.value("geometry").toByteArray())
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 57600
self.baud_rate = baud_rate = 9600
self.gaussian_taps = gaussian_taps = firdes.gaussian(
1.5, 2 * (samp_rate / baud_rate), 0.5, 12)
self.gain_mu = gain_mu = 0.175
##################################################
# Blocks
##################################################
self.zeromq_sub_source_0 = zeromq.sub_source(gr.sizeof_gr_complex, 1,
'tcp://127.0.0.1:5555',
100, False, -1)
self.zeromq_pub_sink_0 = zeromq.pub_sink(gr.sizeof_char, 1,
'tcp://127.0.0.1:38211', 100,
False, -1)
self.satellites_strip_ax25_header_0 = satellites.strip_ax25_header()
self.satellites_nrzi_decode_0 = satellites.nrzi_decode()
self.satellites_hdlc_deframer_0_0 = satellites.hdlc_deframer(
check_fcs=True, max_length=1000)
self.satellites_decode_rs_0 = satellites.decode_rs(True, 0)
self.satellites_check_address_0 = satellites.check_address(
'DP0OPS', "from")
self.qtgui_waterfall_sink_x_0 = qtgui.waterfall_sink_c(
512, #size
firdes.WIN_HAMMING, #wintype
0, #fc
samp_rate, #bw
"OPS-SAT UHF BEACON", #name
1 #number of inputs
)
self.qtgui_waterfall_sink_x_0.set_update_time(0.03)
self.qtgui_waterfall_sink_x_0.enable_grid(False)
self.qtgui_waterfall_sink_x_0.enable_axis_labels(True)
if not True:
self.qtgui_waterfall_sink_x_0.disable_legend()
if "complex" == "float" or "complex" == "msg_float":
self.qtgui_waterfall_sink_x_0.set_plot_pos_half(not True)
labels = ['', '', '', '', '', '', '', '', '', '']
colors = [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_waterfall_sink_x_0.set_line_label(
i, "Data {0}".format(i))
else:
self.qtgui_waterfall_sink_x_0.set_line_label(i, labels[i])
self.qtgui_waterfall_sink_x_0.set_color_map(i, colors[i])
self.qtgui_waterfall_sink_x_0.set_line_alpha(i, alphas[i])
self.qtgui_waterfall_sink_x_0.set_intensity_range(-140, 10)
self._qtgui_waterfall_sink_x_0_win = sip.wrapinstance(
self.qtgui_waterfall_sink_x_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_waterfall_sink_x_0_win)
self.qtgui_freq_sink_x_0 = qtgui.freq_sink_c(
512, #size
firdes.WIN_BLACKMAN_hARRIS, #wintype
0, #fc
samp_rate, #bw
"", #name
1 #number of inputs
)
self.qtgui_freq_sink_x_0.set_update_time(0.10)
self.qtgui_freq_sink_x_0.set_y_axis(-140, 10)
self.qtgui_freq_sink_x_0.set_y_label('Relative Gain', 'dB')
self.qtgui_freq_sink_x_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, 0.0, 0,
"")
self.qtgui_freq_sink_x_0.enable_autoscale(False)
self.qtgui_freq_sink_x_0.enable_grid(False)
self.qtgui_freq_sink_x_0.set_fft_average(0.1)
self.qtgui_freq_sink_x_0.enable_axis_labels(True)
self.qtgui_freq_sink_x_0.enable_control_panel(False)
if not True:
self.qtgui_freq_sink_x_0.disable_legend()
if "complex" == "float" or "complex" == "msg_float":
self.qtgui_freq_sink_x_0.set_plot_pos_half(not True)
labels = ['', '', '', '', '', '', '', '', '', '']
widths = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
colors = [
"blue", "red", "green", "black", "cyan", "magenta", "yellow",
"dark red", "dark green", "dark blue"
]
alphas = [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0]
for i in xrange(1):
if len(labels[i]) == 0:
self.qtgui_freq_sink_x_0.set_line_label(
i, "Data {0}".format(i))
else:
self.qtgui_freq_sink_x_0.set_line_label(i, labels[i])
self.qtgui_freq_sink_x_0.set_line_width(i, widths[i])
self.qtgui_freq_sink_x_0.set_line_color(i, colors[i])
self.qtgui_freq_sink_x_0.set_line_alpha(i, alphas[i])
self._qtgui_freq_sink_x_0_win = sip.wrapinstance(
self.qtgui_freq_sink_x_0.pyqwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_freq_sink_x_0_win)
self.fir_filter_xxx_0 = filter.fir_filter_fff(1, (gaussian_taps))
self.fir_filter_xxx_0.declare_sample_delay(0)
self.digital_descrambler_bb_0_0 = digital.descrambler_bb(0x21, 0, 16)
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(
(samp_rate / baud_rate) * (1 + 0.0), 0.25 * gain_mu * gain_mu, 0.5,
gain_mu, 0.005)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.digital_additive_scrambler_bb_0_0 = digital.additive_scrambler_bb(
0xA9,
0xFF,
7,
count=0,
bits_per_byte=1,
reset_tag_key="packet_len")
self.blocks_unpacked_to_packed_xx_0_0_0_0 = blocks.unpacked_to_packed_bb(
1, gr.GR_MSB_FIRST)
self.blocks_tagged_stream_to_pdu_0_0_0_0_0 = blocks.tagged_stream_to_pdu(
blocks.byte_t, 'packet_len')
self.blocks_pdu_to_tagged_stream_1 = blocks.pdu_to_tagged_stream(
blocks.byte_t, 'packet_len')
self.blocks_pdu_to_tagged_stream_0 = blocks.pdu_to_tagged_stream(
blocks.byte_t, 'packet_len')
self.blocks_packed_to_unpacked_xx_0 = blocks.packed_to_unpacked_bb(
1, gr.GR_MSB_FIRST)
self.blocks_message_debug_0 = blocks.message_debug()
self.analog_quadrature_demod_cf_0 = analog.quadrature_demod_cf(
2 * (samp_rate / baud_rate) / (math.pi))
##################################################
# Connections
##################################################
self.msg_connect((self.blocks_tagged_stream_to_pdu_0_0_0_0_0, 'pdus'),
(self.satellites_decode_rs_0, 'in'))
self.msg_connect((self.satellites_check_address_0, 'ok'),
(self.satellites_strip_ax25_header_0, 'in'))
self.msg_connect((self.satellites_decode_rs_0, 'out'),
(self.blocks_message_debug_0, 'print_pdu'))
self.msg_connect((self.satellites_decode_rs_0, 'out'),
(self.blocks_pdu_to_tagged_stream_1, 'pdus'))
self.msg_connect((self.satellites_hdlc_deframer_0_0, 'out'),
(self.satellites_check_address_0, 'in'))
self.msg_connect((self.satellites_strip_ax25_header_0, 'out'),
(self.blocks_pdu_to_tagged_stream_0, 'pdus'))
self.connect((self.analog_quadrature_demod_cf_0, 0),
(self.fir_filter_xxx_0, 0))
self.connect((self.blocks_packed_to_unpacked_xx_0, 0),
(self.digital_additive_scrambler_bb_0_0, 0))
self.connect((self.blocks_pdu_to_tagged_stream_0, 0),
(self.blocks_packed_to_unpacked_xx_0, 0))
self.connect((self.blocks_pdu_to_tagged_stream_1, 0),
(self.zeromq_pub_sink_0, 0))
self.connect((self.blocks_unpacked_to_packed_xx_0_0_0_0, 0),
(self.blocks_tagged_stream_to_pdu_0_0_0_0_0, 0))
self.connect((self.digital_additive_scrambler_bb_0_0, 0),
(self.blocks_unpacked_to_packed_xx_0_0_0_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0),
(self.satellites_nrzi_decode_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0),
(self.digital_binary_slicer_fb_0, 0))
self.connect((self.digital_descrambler_bb_0_0, 0),
(self.satellites_hdlc_deframer_0_0, 0))
self.connect((self.fir_filter_xxx_0, 0),
(self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.satellites_nrzi_decode_0, 0),
(self.digital_descrambler_bb_0_0, 0))
self.connect((self.zeromq_sub_source_0, 0),
(self.analog_quadrature_demod_cf_0, 0))
self.connect((self.zeromq_sub_source_0, 0),
(self.qtgui_freq_sink_x_0, 0))
self.connect((self.zeromq_sub_source_0, 0),
(self.qtgui_waterfall_sink_x_0, 0))
def __init__(self):
grc_wxgui.top_block_gui.__init__(self, title="Keyfob Decode Wav")
##################################################
# Variables
##################################################
self.sps = sps = 5
self.samp_rate = samp_rate = 48000
##################################################
# Blocks
##################################################
self.wxgui_scopesink2_1 = scopesink2.scope_sink_f(
self.GetWin(),
title="Scope Plot",
sample_rate=3400,
v_scale=0,
v_offset=0,
t_scale=0,
ac_couple=False,
xy_mode=False,
num_inputs=1,
trig_mode=wxgui.TRIG_MODE_AUTO,
y_axis_label="Counts",
)
self.Add(self.wxgui_scopesink2_1.win)
self.wxgui_scopesink2_0 = scopesink2.scope_sink_f(
self.GetWin(),
title="Scope Plot",
sample_rate=3400 * 5,
v_scale=0,
v_offset=0,
t_scale=0,
ac_couple=False,
xy_mode=False,
num_inputs=1,
trig_mode=wxgui.TRIG_MODE_AUTO,
y_axis_label="Counts",
)
self.Add(self.wxgui_scopesink2_0.win)
self.pfb_arb_resampler_xxx_0 = pfb.arb_resampler_fff(
3400.0*sps/48000,
taps=None,
flt_size=32)
self.pfb_arb_resampler_xxx_0.declare_sample_delay(0)
self.keyfob_parse_packet_0 = keyfob.parse_packet()
self.keyfob_manchester_decode_0 = keyfob.manchester_decode()
self.digital_correlate_access_code_tag_bb_0 = digital.correlate_access_code_tag_bb("10101000", 0, "packet_start")
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(5*(1+0.0), 0.25*0.175*0.175, 0.5, 0.175, 0.05)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.blocks_wavfile_source_0 = blocks.wavfile_source("../gqrx_20150306_154200_434400000.wav", False)
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_float*1, samp_rate*3,True)
self.blocks_moving_average_xx_1 = blocks.moving_average_ff(sps, 1, 4000)
##################################################
# Connections
##################################################
self.connect((self.blocks_moving_average_xx_1, 0), (self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.blocks_moving_average_xx_1, 0), (self.wxgui_scopesink2_0, 0))
self.connect((self.blocks_throttle_0, 0), (self.pfb_arb_resampler_xxx_0, 0))
self.connect((self.blocks_wavfile_source_0, 0), (self.blocks_throttle_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.keyfob_manchester_decode_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.wxgui_scopesink2_1, 0))
self.connect((self.digital_correlate_access_code_tag_bb_0, 0), (self.keyfob_parse_packet_0, 0))
self.connect((self.keyfob_manchester_decode_0, 0), (self.digital_correlate_access_code_tag_bb_0, 0))
self.connect((self.pfb_arb_resampler_xxx_0, 0), (self.blocks_moving_average_xx_1, 0))
def __init__(self):
grc_wxgui.top_block_gui.__init__(self, title="Top Block")
_icon_path = "/usr/share/icons/hicolor/32x32/apps/gnuradio-grc.png"
self.SetIcon(wx.Icon(_icon_path, wx.BITMAP_TYPE_ANY))
##################################################
# Variables
##################################################
self.initpathprefix = initpathprefix = "/home/user/alarm-fingerprint/AlarmGnuRadioFiles/"
self.pathprefix = pathprefix = "/home/user/alarm-fingerprint/AlarmGnuRadioFiles/Captured/"
self.finput = finput = initpathprefix + "Capture_init.cap"
self.foutput = foutput = pathprefix + finput.rsplit("/", 1)[1]
self.symb_rate = symb_rate = 4000
self.samp_rate = samp_rate = 10e6
self.decimation = decimation = 100
self.channel_spacing = channel_spacing = 2000000 + 1000000
self.addconst = addconst = 0
self.symb_rate_slider = symb_rate_slider = 4000
self.samp_per_sym = samp_per_sym = int((samp_rate / decimation) / symb_rate)
self.recfile4 = recfile4 = initpathprefix + "/init/_AddConst" + str(addconst) + "_DSC.dat"
self.freq_offset = freq_offset = (channel_spacing / 2) + (channel_spacing * 0.1)
self.freq = freq = 433.92e6
self.channel_trans = channel_trans = 1.2e6
##################################################
# Blocks
##################################################
_channel_trans_sizer = wx.BoxSizer(wx.VERTICAL)
self._channel_trans_text_box = forms.text_box(
parent=self.GetWin(),
sizer=_channel_trans_sizer,
value=self.channel_trans,
callback=self.set_channel_trans,
label="channel_trans",
converter=forms.float_converter(),
proportion=0,
)
self._channel_trans_slider = forms.slider(
parent=self.GetWin(),
sizer=_channel_trans_sizer,
value=self.channel_trans,
callback=self.set_channel_trans,
minimum=0,
maximum=1.8e6,
num_steps=10,
style=wx.SL_HORIZONTAL,
cast=float,
proportion=1,
)
self.Add(_channel_trans_sizer)
_symb_rate_slider_sizer = wx.BoxSizer(wx.VERTICAL)
self._symb_rate_slider_text_box = forms.text_box(
parent=self.GetWin(),
sizer=_symb_rate_slider_sizer,
value=self.symb_rate_slider,
callback=self.set_symb_rate_slider,
label="symb_rate_slider",
converter=forms.float_converter(),
proportion=0,
)
self._symb_rate_slider_slider = forms.slider(
parent=self.GetWin(),
sizer=_symb_rate_slider_sizer,
value=self.symb_rate_slider,
callback=self.set_symb_rate_slider,
minimum=0,
maximum=10e3,
num_steps=100,
style=wx.SL_HORIZONTAL,
cast=float,
proportion=1,
)
self.Add(_symb_rate_slider_sizer)
self.low_pass_filter_0 = filter.fir_filter_fff(
decimation, firdes.low_pass(1, samp_rate, 8e3, 1.8e6, firdes.WIN_BLACKMAN, 6.76)
)
self.freq_xlating_fir_filter_xxx_0 = filter.freq_xlating_fir_filter_ccc(
1,
(firdes.low_pass(1, samp_rate, channel_spacing, channel_trans, firdes.WIN_BLACKMAN, 6.76)),
-freq_offset,
samp_rate,
)
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(
samp_per_sym * (1 + 0.0), 0.25 * 0.175 * 0.175, 0.5, 0.175, 0.005
)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_gr_complex * 1, samp_rate, True)
self.blocks_file_source_0 = blocks.file_source(gr.sizeof_gr_complex * 1, finput, False)
self.blocks_file_sink_1 = blocks.file_sink(gr.sizeof_char * 1, recfile4, False)
self.blocks_file_sink_1.set_unbuffered(False)
self.blocks_complex_to_mag_squared_0 = blocks.complex_to_mag_squared(1)
self.blocks_add_const_vxx_0 = blocks.add_const_vff((addconst,))
##################################################
# Connections
##################################################
self.connect((self.low_pass_filter_0, 0), (self.digital_clock_recovery_mm_xx_0, 0))
self.connect((self.freq_xlating_fir_filter_xxx_0, 0), (self.blocks_complex_to_mag_squared_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0))
self.connect((self.blocks_complex_to_mag_squared_0, 0), (self.blocks_add_const_vxx_0, 0))
self.connect((self.blocks_add_const_vxx_0, 0), (self.low_pass_filter_0, 0))
self.connect((self.blocks_file_source_0, 0), (self.blocks_throttle_0, 0))
self.connect((self.blocks_throttle_0, 0), (self.freq_xlating_fir_filter_xxx_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.blocks_file_sink_1, 0))
def __init__(self):
gr.top_block.__init__(self, "grc_cc1111_hackrf_receiver")
##################################################
# Variables
##################################################
self.symbole_rate = symbole_rate = 40000
self.samp_rate = samp_rate = 2e06
self.rat_interop = rat_interop = 8
self.rat_decim = rat_decim = 5
self.firdes_transition_width = firdes_transition_width = 15000
self.firdes_decim = firdes_decim = 4
self.firdes_cuttoff = firdes_cuttoff = 21e3
self.samp_per_sym = samp_per_sym = ((samp_rate/2/firdes_decim)*rat_interop/rat_decim) / symbole_rate
self.myqueue_out = myqueue_out = gr.msg_queue(2)
self.frequency_shift = frequency_shift = 520000
self.frequency_center = frequency_center = 869.02e06
self.freq_adjust1 = freq_adjust1 = 0
self.firdes_filter = firdes_filter = firdes.low_pass(1,samp_rate/2, firdes_cuttoff, firdes_transition_width)
self.access_code = access_code = '11010011100100011101001110010001'
##################################################
# Blocks
##################################################
self.rational_resampler_xxx_0_0 = filter.rational_resampler_ccc(
interpolation=rat_interop,
decimation=rat_decim,
taps=None,
fractional_bw=None,
)
self.osmosdr_source_0 = osmosdr.source( args="numchan=" + str(1) + " " + "hackrf" )
self.osmosdr_source_0.set_sample_rate(samp_rate)
self.osmosdr_source_0.set_center_freq(frequency_center-frequency_shift, 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(0, 0)
self.osmosdr_source_0.set_gain_mode(False, 0)
self.osmosdr_source_0.set_gain(10, 0)
self.osmosdr_source_0.set_if_gain(20, 0)
self.osmosdr_source_0.set_bb_gain(20, 0)
self.osmosdr_source_0.set_antenna("", 0)
self.osmosdr_source_0.set_bandwidth(0, 0)
self.freq_xlating_fir_filter_xxx_1 = filter.freq_xlating_fir_filter_ccc(2, (1, ), frequency_shift, samp_rate)
self.freq_xlating_fir_filter_xxx_0_0 = filter.freq_xlating_fir_filter_ccc(firdes_decim, (firdes_filter), freq_adjust1, samp_rate/2)
self.digital_correlate_access_code_bb_0_0 = digital.correlate_access_code_bb(access_code, 1)
self.digital_clock_recovery_mm_xx_0_0 = digital.clock_recovery_mm_ff(samp_per_sym*(1+0.0), 0.25*0.175*0.175, 0.5, 0.175, 0.005)
self.digital_binary_slicer_fb_0_0_0 = digital.binary_slicer_fb()
self.cc1111_cc1111_packet_decoder_0 = cc1111.cc1111_packet_decoder(myqueue_out,True, True, False, True)
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_gr_complex*1, samp_rate/2,True)
self.blocks_null_sink_0_0 = blocks.null_sink(gr.sizeof_char*1)
self.analog_quadrature_demod_cf_0_0 = analog.quadrature_demod_cf(2)
##################################################
# Connections
##################################################
self.connect((self.blocks_throttle_0, 0), (self.freq_xlating_fir_filter_xxx_0_0, 0))
self.connect((self.digital_binary_slicer_fb_0_0_0, 0), (self.digital_correlate_access_code_bb_0_0, 0))
self.connect((self.analog_quadrature_demod_cf_0_0, 0), (self.digital_clock_recovery_mm_xx_0_0, 0))
self.connect((self.rational_resampler_xxx_0_0, 0), (self.analog_quadrature_demod_cf_0_0, 0))
self.connect((self.digital_clock_recovery_mm_xx_0_0, 0), (self.digital_binary_slicer_fb_0_0_0, 0))
self.connect((self.freq_xlating_fir_filter_xxx_0_0, 0), (self.rational_resampler_xxx_0_0, 0))
self.connect((self.osmosdr_source_0, 0), (self.freq_xlating_fir_filter_xxx_1, 0))
self.connect((self.freq_xlating_fir_filter_xxx_1, 0), (self.blocks_throttle_0, 0))
self.connect((self.digital_correlate_access_code_bb_0_0, 0), (self.cc1111_cc1111_packet_decoder_0, 0))
self.connect((self.cc1111_cc1111_packet_decoder_0, 0), (self.blocks_null_sink_0_0, 0))
def __init__(self, filepath_in):
gr.top_block.__init__(self)
#grc_wxgui.top_block_gui.__init__(self, title="Top Block")
##################################################
# Variables
##################################################
self.samp_rate = samp_rate = 200e3
self.bb_interpolation = bb_interpolation = 100
self.bb_decimation = bb_decimation = 612
self.samples_per_symbol = samples_per_symbol = 4
self.gain_mu = gain_mu = 0.03
self.bb_rate = bb_rate = float(samp_rate) * bb_interpolation / bb_decimation
self.bb_filter_freq = bb_filter_freq = 10e3
self.omega = omega = samples_per_symbol
self.mu = mu = 0.5
self.gain_omega = gain_omega = 0.25 * gain_mu * gain_mu
self.bb_taps = bb_taps = gr.firdes.low_pass(1.0, samp_rate, bb_filter_freq, bb_filter_freq * 0.1)
self.baud_rate = baud_rate = bb_rate / samples_per_symbol
#self.average = average = 64
##################################################
# Blocks
##################################################
# self.wxgui_scopesink2_1_0_0 = scopesink2.scope_sink_f(
# self.GetWin(),
# title="Scope Plot",
# sample_rate=baud_rate,
# v_scale=0,
# v_offset=0,
# t_scale=0,
# ac_couple=False,
# xy_mode=False,
# num_inputs=1,
# trig_mode=gr.gr_TRIG_MODE_NORM,
# y_axis_label="Counts",
# )
# self.Add(self.wxgui_scopesink2_1_0_0.win)
#self.freq_xlating_fir_filter_xxx_0 = filter.freq_xlating_fir_filter_ccc(1, (bb_taps), 6e3, samp_rate)
self.digital_correlate_access_code_bb_0 = digital.correlate_access_code_bb("10101010101010101010101010101", 1)
self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(omega, gain_omega, mu, gain_mu, 0.0002)
self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb()
self.dc_blocker_xx_0 = filter.dc_blocker_ff(64, True)
#self.blocks_uchar_to_float_0_0 = blocks.uchar_to_float()
#self.blocks_throttle_0 = blocks.throttle(gr.sizeof_gr_complex*1, samp_rate)
#self.blocks_file_source_0_0 = blocks.file_source(gr.sizeof_gr_complex*1, "/mnt/hgfs/tmp/rf_captures/315.000m_200.000k_20130623_133451_extract_am_2.cfile", True)
self.blocks_file_source_0_0 = blocks.file_source(gr.sizeof_gr_complex*1, filepath_in, False)
self.blocks_complex_to_mag_0 = blocks.complex_to_mag(1)
self.blks2_rational_resampler_xxx_0 = blks2.rational_resampler_fff(
interpolation=bb_interpolation,
decimation=bb_decimation,
taps=None,
fractional_bw=None,
)
# _bb_filter_freq_sizer = wx.BoxSizer(wx.VERTICAL)
# self._bb_filter_freq_text_box = forms.text_box(
# parent=self.GetWin(),
# sizer=_bb_filter_freq_sizer,
# value=self.bb_filter_freq,
# callback=self.set_bb_filter_freq,
# label="BB Freq",
# converter=forms.int_converter(),
# proportion=0,
# )
# self._bb_filter_freq_slider = forms.slider(
# parent=self.GetWin(),
# sizer=_bb_filter_freq_sizer,
# value=self.bb_filter_freq,
# callback=self.set_bb_filter_freq,
# minimum=5e3,
# maximum=30e3,
# num_steps=250,
# style=wx.SL_HORIZONTAL,
# cast=int,
# proportion=1,
# )
# self.Add(_bb_filter_freq_sizer)
# _average_sizer = wx.BoxSizer(wx.VERTICAL)
# self._average_text_box = forms.text_box(
# parent=self.GetWin(),
# sizer=_average_sizer,
# value=self.average,
# callback=self.set_average,
# label="Average Length",
# converter=forms.int_converter(),
# proportion=0,
# )
# self._average_slider = forms.slider(
# parent=self.GetWin(),
# sizer=_average_sizer,
# value=self.average,
# callback=self.set_average,
# minimum=0,
# maximum=256,
# num_steps=256,
# style=wx.SL_HORIZONTAL,
# cast=int,
# proportion=1,
# )
# self.Add(_average_sizer)
##################################################
# Connections
##################################################
self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0))
self.connect((self.digital_binary_slicer_fb_0, 0), (self.digital_correlate_access_code_bb_0, 0))
#self.connect((self.digital_correlate_access_code_bb_0, 0), (self.blocks_uchar_to_float_0_0, 0))
#self.connect((self.blocks_throttle_0, 0), (self.freq_xlating_fir_filter_xxx_0, 0))
#self.connect((self.freq_xlating_fir_filter_xxx_0, 0), (self.blocks_complex_to_mag_0, 0))
self.connect((self.blocks_complex_to_mag_0, 0), (self.blks2_rational_resampler_xxx_0, 0))
self.connect((self.blks2_rational_resampler_xxx_0, 0), (self.dc_blocker_xx_0, 0))
self.connect((self.dc_blocker_xx_0, 0), (self.digital_clock_recovery_mm_xx_0, 0))
#self.connect((self.blocks_uchar_to_float_0_0, 0), (self.wxgui_scopesink2_1_0_0, 0))
#self.connect((self.blocks_file_source_0_0, 0), (self.blocks_throttle_0, 0))
#self.connect((self.blocks_file_source_0_0, 0), (self.freq_xlating_fir_filter_xxx_0, 0))
self.connect((self.blocks_file_source_0_0, 0), (self.blocks_complex_to_mag_0, 0))
self.packetizer = Packetizer(82)
self.connect((self.digital_correlate_access_code_bb_0, 0), (self.packetizer, 0))
def __init__(self): grc_wxgui.top_block_gui.__init__(self, title="APRS Receiver") ################################################## # Variables ################################################## self.space = space = 1200 self.mark = mark = 2200 self.xlate_decim = xlate_decim = 8 self.xlate_bandwidth = xlate_bandwidth = 1200*6 self.sym_dev = sym_dev = (mark-space)/2 self.samp_rate = samp_rate = 1e6 self.quad_rate = quad_rate = 96000 self.gain = gain = 10 self.freq_offset = freq_offset = 390e3 self.freq = freq = 144e6 self.baud = baud = 1200 self.audio_rate = audio_rate = 48000 self.audio_mul = audio_mul = 1 self.aprs_rate = aprs_rate = 12000 self.ant = ant = 'TX/RX' ################################################## # Message Queues ################################################## ax25_hdlc_framer_b_0_msgq_out = ax25_print_frame_0_msgq_in = gr.msg_queue(2) ################################################## # Blocks ################################################## self.nb = self.nb = wx.Notebook(self.GetWin(), style=wx.NB_TOP) self.nb.AddPage(grc_wxgui.Panel(self.nb), "Baseband") self.nb.AddPage(grc_wxgui.Panel(self.nb), "Waterfall") self.nb.AddPage(grc_wxgui.Panel(self.nb), "Signal") self.nb.AddPage(grc_wxgui.Panel(self.nb), "Slicer") self.nb.AddPage(grc_wxgui.Panel(self.nb), "Eye") self.Add(self.nb) _gain_sizer = wx.BoxSizer(wx.VERTICAL) self._gain_text_box = forms.text_box( parent=self.GetWin(), sizer=_gain_sizer, value=self.gain, callback=self.set_gain, label="RF Gain", converter=forms.float_converter(), proportion=0, ) self._gain_slider = forms.slider( parent=self.GetWin(), sizer=_gain_sizer, value=self.gain, callback=self.set_gain, minimum=0, maximum=50, num_steps=50, style=wx.SL_HORIZONTAL, cast=float, proportion=1, ) self.Add(_gain_sizer) _freq_offset_sizer = wx.BoxSizer(wx.VERTICAL) self._freq_offset_text_box = forms.text_box( parent=self.GetWin(), sizer=_freq_offset_sizer, value=self.freq_offset, callback=self.set_freq_offset, label="Freq Offset", converter=forms.float_converter(), proportion=0, ) self._freq_offset_slider = forms.slider( parent=self.GetWin(), sizer=_freq_offset_sizer, value=self.freq_offset, callback=self.set_freq_offset, minimum=-500e3, maximum=500e3, num_steps=1000, style=wx.SL_HORIZONTAL, cast=float, proportion=1, ) self.Add(_freq_offset_sizer) self._freq_text_box = forms.text_box( parent=self.GetWin(), value=self.freq, callback=self.set_freq, label="Freq", converter=forms.float_converter(), ) self.Add(self._freq_text_box) _audio_mul_sizer = wx.BoxSizer(wx.VERTICAL) self._audio_mul_text_box = forms.text_box( parent=self.GetWin(), sizer=_audio_mul_sizer, value=self.audio_mul, callback=self.set_audio_mul, label="Audio", converter=forms.float_converter(), proportion=0, ) self._audio_mul_slider = forms.slider( parent=self.GetWin(), sizer=_audio_mul_sizer, value=self.audio_mul, callback=self.set_audio_mul, minimum=0, maximum=10, num_steps=1000, style=wx.SL_HORIZONTAL, cast=float, proportion=1, ) self.Add(_audio_mul_sizer) self._ant_chooser = forms.drop_down( parent=self.GetWin(), value=self.ant, callback=self.set_ant, label="Antenna", choices=['TX/RX', 'RX2'], labels=[], ) self.Add(self._ant_chooser) self.wxgui_waterfallsink2_0 = waterfallsink2.waterfall_sink_c( self.nb.GetPage(1).GetWin(), baseband_freq=0, dynamic_range=50, ref_level=-65, ref_scale=2.0, sample_rate=aprs_rate, fft_size=512, fft_rate=15, average=False, avg_alpha=None, title="Waterfall Plot", ) self.nb.GetPage(1).Add(self.wxgui_waterfallsink2_0.win) self.wxgui_scopesink2_0_0_0 = scopesink2.scope_sink_f( self.nb.GetPage(4).GetWin(), title="Scope Plot", sample_rate=aprs_rate/10, v_scale=0.5, v_offset=0, t_scale=0.002, ac_couple=False, xy_mode=False, num_inputs=1, trig_mode=gr.gr_TRIG_MODE_AUTO, y_axis_label="Counts", ) self.nb.GetPage(4).Add(self.wxgui_scopesink2_0_0_0.win) self.wxgui_scopesink2_0_0 = scopesink2.scope_sink_f( self.nb.GetPage(3).GetWin(), title="Scope Plot", sample_rate=aprs_rate, v_scale=0.5, v_offset=0, t_scale=0.002, ac_couple=False, xy_mode=False, num_inputs=1, trig_mode=gr.gr_TRIG_MODE_AUTO, y_axis_label="Counts", ) self.nb.GetPage(3).Add(self.wxgui_scopesink2_0_0.win) self.wxgui_scopesink2_0 = scopesink2.scope_sink_f( self.nb.GetPage(2).GetWin(), title="Scope Plot", sample_rate=aprs_rate, v_scale=0.05, v_offset=0, t_scale=0.002, ac_couple=False, xy_mode=False, num_inputs=1, trig_mode=gr.gr_TRIG_MODE_AUTO, y_axis_label="Counts", ) self.nb.GetPage(2).Add(self.wxgui_scopesink2_0.win) self.wxgui_fftsink2_0 = fftsink2.fft_sink_c( self.nb.GetPage(0).GetWin(), baseband_freq=0, y_per_div=10, y_divs=10, ref_level=-20, ref_scale=2.0, sample_rate=samp_rate, fft_size=1024, fft_rate=15, average=True, avg_alpha=0.5, title="FFT Plot", peak_hold=False, ) self.nb.GetPage(0).Add(self.wxgui_fftsink2_0.win) def wxgui_fftsink2_0_callback(x, y): self.set_freq_offset(x) self.wxgui_fftsink2_0.set_callback(wxgui_fftsink2_0_callback) self.uhd_usrp_source_0 = uhd.usrp_source( device_addr="", stream_args=uhd.stream_args( cpu_format="fc32", channels=range(1), ), ) self.uhd_usrp_source_0.set_samp_rate(samp_rate) self.uhd_usrp_source_0.set_center_freq(freq, 0) self.uhd_usrp_source_0.set_gain(gain, 0) self.uhd_usrp_source_0.set_antenna(ant, 0) self.low_pass_filter_0 = gr.fir_filter_ccf(1, firdes.low_pass( 1, aprs_rate, 2e3, 600, firdes.WIN_HAMMING, 6.76)) self.gr_single_pole_iir_filter_xx_0 = gr.single_pole_iir_filter_ff(0.0001, 1) self.gr_null_sink_0 = gr.null_sink(gr.sizeof_float*1) self.gr_multiply_xx_0 = gr.multiply_vcc(1) self.gr_multiply_const_vxx_0 = gr.multiply_const_vff((audio_mul, )) self.gr_agc_xx_1 = gr.agc_ff(1e-3, 0.8, 0.1, 10.0) self.freq_xlating_fir_filter_xxx_0 = filter.freq_xlating_fir_filter_ccc(xlate_decim, (firdes.low_pass(1, samp_rate, xlate_bandwidth/2, 1000)), freq_offset, samp_rate) self.digital_clock_recovery_mm_xx_0 = digital.clock_recovery_mm_ff(10, .25 * (0.05)**2, 0.5, 0.005, 0.005) self.digital_binary_slicer_fb_0 = digital.binary_slicer_fb() self.blocks_sub_xx_0 = blocks.sub_ff(1) self.blocks_float_to_complex_0 = blocks.float_to_complex(1) self.blks2_rational_resampler_xxx_0_0 = blks2.rational_resampler_ccc( interpolation=quad_rate, decimation=int(samp_rate/xlate_decim), taps=None, fractional_bw=None, ) self.blks2_rational_resampler_xxx_0 = blks2.rational_resampler_ccc( interpolation=aprs_rate, decimation=quad_rate, taps=None, fractional_bw=None, ) self.blks2_nbfm_rx_0_0 = blks2.nbfm_rx( audio_rate=audio_rate, quad_rate=quad_rate, tau=75e-6, max_dev=25000, ) self.blks2_nbfm_rx_0 = blks2.nbfm_rx( audio_rate=aprs_rate, quad_rate=quad_rate, tau=75e-6, max_dev=3e3, ) self.ax25_print_frame_0 = packetradio.queue_watcher_thread(ax25_print_frame_0_msgq_in) self.ax25_hdlc_framer_b_0 = packetradio.hdlc_framer(ax25_hdlc_framer_b_0_msgq_out, False) self.analog_sig_source_x_0 = analog.sig_source_c(aprs_rate, analog.GR_SIN_WAVE, -(min(mark,space)+sym_dev), 1, 0) self.analog_quadrature_demod_cf_0 = analog.quadrature_demod_cf(aprs_rate/(2*math.pi*sym_dev)) self.analog_pwr_squelch_xx_0_0_0 = analog.pwr_squelch_cc(-70, 1e-1, 0, False) self.analog_pwr_squelch_xx_0_0 = analog.pwr_squelch_cc(-70, 1e-1, 0, False) ################################################## # Connections ################################################## self.connect((self.uhd_usrp_source_0, 0), (self.wxgui_fftsink2_0, 0)) self.connect((self.blks2_rational_resampler_xxx_0_0, 0), (self.blks2_rational_resampler_xxx_0, 0)) self.connect((self.freq_xlating_fir_filter_xxx_0, 0), (self.blks2_rational_resampler_xxx_0_0, 0)) self.connect((self.blks2_rational_resampler_xxx_0, 0), (self.wxgui_waterfallsink2_0, 0)) self.connect((self.uhd_usrp_source_0, 0), (self.freq_xlating_fir_filter_xxx_0, 0)) self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.digital_binary_slicer_fb_0, 0)) self.connect((self.blocks_sub_xx_0, 0), (self.digital_clock_recovery_mm_xx_0, 0)) self.connect((self.gr_single_pole_iir_filter_xx_0, 0), (self.blocks_sub_xx_0, 1)) self.connect((self.analog_quadrature_demod_cf_0, 0), (self.gr_single_pole_iir_filter_xx_0, 0)) self.connect((self.analog_quadrature_demod_cf_0, 0), (self.blocks_sub_xx_0, 0)) self.connect((self.low_pass_filter_0, 0), (self.analog_quadrature_demod_cf_0, 0)) self.connect((self.gr_multiply_xx_0, 0), (self.low_pass_filter_0, 0)) self.connect((self.analog_sig_source_x_0, 0), (self.gr_multiply_xx_0, 1)) self.connect((self.blocks_float_to_complex_0, 0), (self.gr_multiply_xx_0, 0)) self.connect((self.blks2_nbfm_rx_0, 0), (self.blocks_float_to_complex_0, 0)) self.connect((self.blks2_nbfm_rx_0, 0), (self.blocks_float_to_complex_0, 1)) self.connect((self.digital_clock_recovery_mm_xx_0, 0), (self.wxgui_scopesink2_0_0_0, 0)) self.connect((self.digital_binary_slicer_fb_0, 0), (self.ax25_hdlc_framer_b_0, 0)) self.connect((self.blks2_nbfm_rx_0, 0), (self.wxgui_scopesink2_0, 0)) self.connect((self.blks2_rational_resampler_xxx_0_0, 0), (self.analog_pwr_squelch_xx_0_0, 0)) self.connect((self.analog_pwr_squelch_xx_0_0_0, 0), (self.blks2_nbfm_rx_0, 0)) self.connect((self.blks2_rational_resampler_xxx_0_0, 0), (self.analog_pwr_squelch_xx_0_0_0, 0)) self.connect((self.analog_quadrature_demod_cf_0, 0), (self.wxgui_scopesink2_0_0, 0)) self.connect((self.blks2_nbfm_rx_0_0, 0), (self.gr_agc_xx_1, 0)) self.connect((self.analog_pwr_squelch_xx_0_0, 0), (self.blks2_nbfm_rx_0_0, 0)) self.connect((self.gr_agc_xx_1, 0), (self.gr_multiply_const_vxx_0, 0)) self.connect((self.gr_multiply_const_vxx_0, 0), (self.gr_null_sink_0, 0))
