1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
|
#!/usr/bin/env python
#
# Copyright 2007 Free Software Foundation, Inc.
#
# This file is part of GNU Radio
#
# GNU Radio is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2, or (at your option)
# any later version.
#
# GNU Radio is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with GNU Radio; see the file COPYING. If not, write to
# the Free Software Foundation, Inc., 51 Franklin Street,
# Boston, MA 02110-1301, USA.
#
from gnuradio import gr, eng_notation
from gnuradio.eng_option import eng_option
from optparse import OptionParser
from usrp_source import usrp_source_c
from sounder_rx import sounder_rx
n2s = eng_notation.num_to_str
class usrp_sounder_rx(gr.hier_block2):
def __init__(self, options):
gr.hier_block2.__init__(self, "usrp_sounder_rx",
gr.io_signature(0,0,0),
gr.io_signature(0,0,0))
self._options = options
self._u = usrp_source_c(0,
self._options.rx_subdev_spec,
self._options.gain,
self._options.chip_rate,
self._options.freq,
self._options.cal,
self._options.verbose)
self._options.chip_rate = self._u._if_rate
self._length = 2**self._options.degree-1
self._receiver = sounder_rx(self._options.chip_rate,
self._options.degree,
self._options.verbose)
samples = 100 * self._length**2
self.define_component("usrp", self._u)
self.define_component("head", gr.head(gr.sizeof_gr_complex, samples))
self.define_component("rx", self._receiver)
self.define_component("c2m", gr.complex_to_mag())
self.define_component("s2v", gr.stream_to_vector(gr.sizeof_float, self._length))
self.define_component("lpf", gr.single_pole_iir_filter_ff(self._options.alpha, self._length))
self.define_component("v2s", gr.vector_to_stream(gr.sizeof_float, self._length))
self.define_component("sink", gr.file_sink(gr.sizeof_float, "impulse.dat"))
self.connect("usrp", 0, "head", 0)
self.connect("head", 0, "rx", 0)
self.connect("rx", 0, "c2m", 0)
self.connect("c2m", 0, "s2v", 0)
self.connect("s2v", 0, "lpf", 0)
self.connect("lpf", 0, "v2s", 0)
self.connect("v2s", 0, "sink", 0)
if self._options.verbose:
print "Chip rate is", n2s(self._options.chip_rate), "chips/sec"
print "Resolution is", n2s(1.0/self._options.chip_rate), "sec"
print "Using PN code of degree", self._options.degree
def main():
parser = OptionParser(option_class=eng_option)
parser.add_option("-R", "--rx-subdev-spec", type="subdev", default=None,
help="select USRP Rx side A or B (default=first found)")
parser.add_option("-f", "--freq", type="eng_float", default=0.0,
help="set center frequency (default=%default)")
parser.add_option("-c", "--cal", type="eng_float", default=0.0,
help="set frequency calibration offset (default=%default)")
parser.add_option("-v", "--verbose", action="store_true", default=False,
help="print extra debugging info")
parser.add_option("-d", "--degree", type="int", default=10,
help="set PN code degree (length=2**degree-1, default=%default)")
parser.add_option("-r", "--chip-rate", type="eng_float", default=8e6,
help="set sounder chip rate (default=%default)")
parser.add_option("-g", "--gain", type="eng_float", default=None,
help="set receiver gain (default=%default)")
parser.add_option("", "--alpha", type="eng_float", default=1.0,
help="set smoothing constant (default=%default)")
(options, args) = parser.parse_args()
if len(args) != 0:
parser.print_help()
sys.exit(1)
top_block = usrp_sounder_rx(options)
runtime = gr.runtime(top_block)
try:
runtime.run()
except KeyboardInterrupt:
pass
if __name__ == '__main__':
main ()
|