diff options
| -rw-r--r-- | config/grc_gnuradio_examples.m4 | 1 | ||||
| -rw-r--r-- | gnuradio-core/src/lib/general/Makefile.am | 13 | ||||
| -rw-r--r-- | gnuradio-core/src/lib/general/general.i | 4 | ||||
| -rw-r--r-- | gnuradio-core/src/lib/general/gr_probe_density_b.cc | 68 | ||||
| -rw-r--r-- | gnuradio-core/src/lib/general/gr_probe_density_b.h | 72 | ||||
| -rw-r--r-- | gnuradio-core/src/lib/general/gr_probe_density_b.i | 36 | ||||
| -rw-r--r-- | gnuradio-core/src/lib/general/gr_probe_mpsk_snr_c.cc | 85 | ||||
| -rw-r--r-- | gnuradio-core/src/lib/general/gr_probe_mpsk_snr_c.h | 82 | ||||
| -rw-r--r-- | gnuradio-core/src/lib/general/gr_probe_mpsk_snr_c.i | 39 | ||||
| -rw-r--r-- | gnuradio-examples/python/Makefile.am | 1 | ||||
| -rw-r--r-- | gnuradio-examples/python/digital-bert/Makefile.am | 33 | ||||
| -rw-r--r-- | gnuradio-examples/python/digital-bert/README | 62 | ||||
| -rwxr-xr-x | gnuradio-examples/python/digital-bert/benchmark_rx.py | 164 | ||||
| -rwxr-xr-x | gnuradio-examples/python/digital-bert/benchmark_tx.py | 111 | ||||
| -rw-r--r-- | gnuradio-examples/python/digital-bert/receive_path.py | 116 | ||||
| -rw-r--r-- | gnuradio-examples/python/digital-bert/transmit_path.py | 62 |
16 files changed, 945 insertions, 4 deletions
diff --git a/config/grc_gnuradio_examples.m4 b/config/grc_gnuradio_examples.m4 index 451316071..59874ca3c 100644 --- a/config/grc_gnuradio_examples.m4 +++ b/config/grc_gnuradio_examples.m4 @@ -34,6 +34,7 @@ AC_DEFUN([GRC_GNURADIO_EXAMPLES],[ gnuradio-examples/python/dect/Makefile \ gnuradio-examples/python/digital/Makefile \ gnuradio-examples/python/digital_voice/Makefile \ + gnuradio-examples/python/digital-bert/Makefile \ gnuradio-examples/python/mp-sched/Makefile \ gnuradio-examples/python/multi-antenna/Makefile \ gnuradio-examples/python/multi_usrp/Makefile \ diff --git a/gnuradio-core/src/lib/general/Makefile.am b/gnuradio-core/src/lib/general/Makefile.am index 024363253..b46590070 100644 --- a/gnuradio-core/src/lib/general/Makefile.am +++ b/gnuradio-core/src/lib/general/Makefile.am @@ -165,7 +165,9 @@ libgeneral_la_SOURCES = \ malloc16.c \ gr_unpack_k_bits_bb.cc \ gr_descrambler_bb.cc \ - gr_scrambler_bb.cc + gr_scrambler_bb.cc \ + gr_probe_mpsk_snr_c.cc \ + gr_probe_density_b.cc libgeneral_qa_la_SOURCES = \ qa_general.cc \ @@ -324,8 +326,9 @@ grinclude_HEADERS = \ random.h \ gr_unpack_k_bits_bb.h \ gr_descrambler_bb.h \ - gr_scrambler_bb.h - + gr_scrambler_bb.h \ + gr_probe_mpsk_snr_c.h \ + gr_probe_density_b.h noinst_HEADERS = \ qa_general.h \ @@ -454,7 +457,9 @@ swiginclude_HEADERS = \ gri_agc2_cc.i \ gri_agc2_ff.i \ gr_descrambler_bb.i \ - gr_scrambler_bb.i + gr_scrambler_bb.i \ + gr_probe_mpsk_snr_c.i \ + gr_probe_density_b.i CLEANFILES = $(BUILT_SOURCES) *.pyc diff --git a/gnuradio-core/src/lib/general/general.i b/gnuradio-core/src/lib/general/general.i index 8326b27b0..542bf6064 100644 --- a/gnuradio-core/src/lib/general/general.i +++ b/gnuradio-core/src/lib/general/general.i @@ -131,6 +131,8 @@ #include <gr_decode_ccsds_27_fb.h> #include <gr_descrambler_bb.h> #include <gr_scrambler_bb.h> +#include <gr_probe_mpsk_snr_c.h> +#include <gr_probe_density_b.h> %} %include "gr_nop.i" @@ -242,3 +244,5 @@ %include "gr_decode_ccsds_27_fb.i" %include "gr_descrambler_bb.i" %include "gr_scrambler_bb.i" +%include "gr_probe_mpsk_snr_c.i" +%include "gr_probe_density_b.i" diff --git a/gnuradio-core/src/lib/general/gr_probe_density_b.cc b/gnuradio-core/src/lib/general/gr_probe_density_b.cc new file mode 100644 index 000000000..dae4eec74 --- /dev/null +++ b/gnuradio-core/src/lib/general/gr_probe_density_b.cc @@ -0,0 +1,68 @@ +/* -*- c++ -*- */ +/* + * Copyright 2008 Free Software Foundation, Inc. + * + * 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 3, 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. + */ + +#ifdef HAVE_CONFIG_H +#include <config.h> +#endif + +#include <gr_probe_density_b.h> +#include <gr_io_signature.h> +#include <stdexcept> +#include <iostream> + +gr_probe_density_b_sptr +gr_make_probe_density_b(double alpha) +{ + return gr_probe_density_b_sptr(new gr_probe_density_b(alpha)); +} + +gr_probe_density_b::gr_probe_density_b(double alpha) + : gr_sync_block("density_b", + gr_make_io_signature(1, 1, sizeof(char)), + gr_make_io_signature(0, 0, 0)) +{ + set_alpha(alpha); + d_density = 1.0; +} + +gr_probe_density_b::~gr_probe_density_b() +{ +} + +int +gr_probe_density_b::work(int noutput_items, + gr_vector_const_void_star &input_items, + gr_vector_void_star &output_items) +{ + const char *in = (const char *)input_items[0]; + + for (int i = 0; i < noutput_items; i++) + d_density = d_alpha*(double)in[i] + d_beta*d_density; + + return noutput_items; +} + +void +gr_probe_density_b::set_alpha(double alpha) +{ + d_alpha = alpha; + d_beta = 1.0-d_alpha; +} + diff --git a/gnuradio-core/src/lib/general/gr_probe_density_b.h b/gnuradio-core/src/lib/general/gr_probe_density_b.h new file mode 100644 index 000000000..92b98d193 --- /dev/null +++ b/gnuradio-core/src/lib/general/gr_probe_density_b.h @@ -0,0 +1,72 @@ +/* -*- c++ -*- */ +/* + * Copyright 2008 Free Software Foundation, Inc. + * + * 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 3, 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. + */ +#ifndef INCLUDED_GR_PROBE_DENSITY_B_H +#define INCLUDED_GR_PROBE_DENSITY_B_H + +#include <gr_sync_block.h> + +class gr_probe_density_b; + +typedef boost::shared_ptr<gr_probe_density_b> gr_probe_density_b_sptr; + +gr_probe_density_b_sptr gr_make_probe_density_b(double alpha); + +/*! + * This block maintains a running average of the input stream and + * makes it available as an accessor function. The input stream + * is type unsigned char. + * + * If you send this block a stream of unpacked bytes, it will tell + * you what the bit density is. + * + * \param alpha Average filter constant + * + */ + +class gr_probe_density_b : public gr_sync_block +{ +private: + friend gr_probe_density_b_sptr gr_make_probe_density_b(double alpha); + + double d_alpha; + double d_beta; + double d_density; + + gr_probe_density_b(double alpha); + +public: + ~gr_probe_density_b(); + + /*! + * \brief Returns the current density value + */ + double density() const { return d_density; } + + /*! + * \brief Set the average filter constant + */ + void set_alpha(double alpha); + + int work(int noutput_items, + gr_vector_const_void_star &input_items, + gr_vector_void_star &output_items); +}; + +#endif /* INCLUDED_GR_PROBE_DENSITY_B_H */ diff --git a/gnuradio-core/src/lib/general/gr_probe_density_b.i b/gnuradio-core/src/lib/general/gr_probe_density_b.i new file mode 100644 index 000000000..54ed07cf5 --- /dev/null +++ b/gnuradio-core/src/lib/general/gr_probe_density_b.i @@ -0,0 +1,36 @@ +/* -*- c++ -*- */ +/* + * Copyright 2008 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 3, 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. + */ + +GR_SWIG_BLOCK_MAGIC(gr,probe_density_b); + +gr_probe_density_b_sptr gr_make_probe_density_b(double alpha); + +class gr_probe_density_b : public gr_sync_block +{ +public: + double density() const; + + void set_alpha(double alpha); + +private: + gr_probe_density_b(); +}; diff --git a/gnuradio-core/src/lib/general/gr_probe_mpsk_snr_c.cc b/gnuradio-core/src/lib/general/gr_probe_mpsk_snr_c.cc new file mode 100644 index 000000000..a0bb99ce1 --- /dev/null +++ b/gnuradio-core/src/lib/general/gr_probe_mpsk_snr_c.cc @@ -0,0 +1,85 @@ +/* -*- c++ -*- */ +/* + * Copyright 2008 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 3, 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. + */ + +#ifdef HAVE_CONFIG_H +#include <config.h> +#endif + +#include <gr_probe_mpsk_snr_c.h> +#include <gr_io_signature.h> + +gr_probe_mpsk_snr_c_sptr +gr_make_probe_mpsk_snr_c(double alpha) +{ + return gr_probe_mpsk_snr_c_sptr(new gr_probe_mpsk_snr_c(alpha)); +} + +gr_probe_mpsk_snr_c::gr_probe_mpsk_snr_c(double alpha) + : gr_sync_block ("probe_mpsk_snr_c", + gr_make_io_signature(1, 1, sizeof(gr_complex)), + gr_make_io_signature(0, 0, 0)), + d_signal_mean(0.0), + d_noise_variance(0.0) +{ + set_alpha(alpha); +} + +gr_probe_mpsk_snr_c::~gr_probe_mpsk_snr_c() +{ +} + +int +gr_probe_mpsk_snr_c::work(int noutput_items, + gr_vector_const_void_star &input_items, + gr_vector_void_star &output_items) +{ + const gr_complex *in = (const gr_complex *) input_items[0]; + + for (int i = 0; i < noutput_items; i++){ + // Update of signal mean estimate + double mag = abs(in[i]); + d_signal_mean = d_alpha*abs(in[i]) + d_beta*d_signal_mean; + + // Update noise variance estimate + double noise = mag-d_signal_mean; + double var = noise*noise; + d_noise_variance = d_alpha*var + d_beta*d_noise_variance; + } + + return noutput_items; +} + +double +gr_probe_mpsk_snr_c::snr() const +{ + if (d_noise_variance == 0.0) + return 0.0; + else + return 10*log10(d_signal_mean*d_signal_mean/d_noise_variance); +} + +void +gr_probe_mpsk_snr_c::set_alpha(double alpha) +{ + d_alpha = alpha; + d_beta = 1.0-alpha; +} diff --git a/gnuradio-core/src/lib/general/gr_probe_mpsk_snr_c.h b/gnuradio-core/src/lib/general/gr_probe_mpsk_snr_c.h new file mode 100644 index 000000000..67492f686 --- /dev/null +++ b/gnuradio-core/src/lib/general/gr_probe_mpsk_snr_c.h @@ -0,0 +1,82 @@ +/* -*- c++ -*- */ +/* + * Copyright 2008 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 3, 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. + */ +#ifndef INCLUDED_GR_PROBE_MPSK_SNR_C_H +#define INCLUDED_GR_PROBE_MPSK_SNR_C_H + +#include <gr_sync_block.h> + +class gr_probe_mpsk_snr_c; +typedef boost::shared_ptr<gr_probe_mpsk_snr_c> gr_probe_mpsk_snr_c_sptr; + +gr_probe_mpsk_snr_c_sptr +gr_make_probe_mpsk_snr_c(double alpha = 0.0001); + +/*! + * Compute the estimate SNR of an MPSK signal using the Squared Signal + * to Noise Variance (SNV) technique. + * + * This technique assumes an AWGN channel. + * + * \param alpha Mean and variance smoothing filter constant + * \ingroup sink + * + * Compute the running average of the signal mean and noise variance. + * The estimated signal mean, noise variance, and SNR are available + * via accessors. + * + * This SNR estimator is inaccurate below about 7dB SNR. + * + */ +class gr_probe_mpsk_snr_c : public gr_sync_block +{ + double d_alpha; + double d_beta; + double d_signal_mean; + double d_noise_variance; + + // Factory function returning shared pointer of this class + friend gr_probe_mpsk_snr_c_sptr + gr_make_probe_mpsk_snr_c(double alpha); + + // Private constructor + gr_probe_mpsk_snr_c(double alpha); + +public: + ~gr_probe_mpsk_snr_c(); + + int work (int noutput_items, + gr_vector_const_void_star &input_items, + gr_vector_void_star &output_items); + + // Return the estimated signal mean + double signal_mean() const { return d_signal_mean; } + + // Return the estimated noise variance + double noise_variance() const { return d_noise_variance; } + + // Return the estimated signal-to-noise ratio in decibels + double snr() const; + + void set_alpha(double alpha); +}; + +#endif /* INCLUDED_GR_PROBE_MPSK_SNR_C_H */ diff --git a/gnuradio-core/src/lib/general/gr_probe_mpsk_snr_c.i b/gnuradio-core/src/lib/general/gr_probe_mpsk_snr_c.i new file mode 100644 index 000000000..37a86b23d --- /dev/null +++ b/gnuradio-core/src/lib/general/gr_probe_mpsk_snr_c.i @@ -0,0 +1,39 @@ +/* -*- c++ -*- */ +/* + * Copyright 2008 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 3, 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. + */ + +GR_SWIG_BLOCK_MAGIC(gr,probe_mpsk_snr_c); + +gr_probe_mpsk_snr_c_sptr +gr_make_probe_mpsk_snr_c(double alpha = 0.0001); + +class gr_probe_mpsk_snr_c : public gr_sync_block +{ +private: + void gr_probe_mpsk_snr_c(double alpha); + +public: + double signal_mean(); + double noise_variance(); + double snr(); + + void set_alpha (double alpha); +}; diff --git a/gnuradio-examples/python/Makefile.am b/gnuradio-examples/python/Makefile.am index b1780bc6f..2de55d8d8 100644 --- a/gnuradio-examples/python/Makefile.am +++ b/gnuradio-examples/python/Makefile.am @@ -26,6 +26,7 @@ SUBDIRS = \ audio \ dect \ digital \ + digital-bert \ digital_voice \ mp-sched \ multi-antenna \ diff --git a/gnuradio-examples/python/digital-bert/Makefile.am b/gnuradio-examples/python/digital-bert/Makefile.am new file mode 100644 index 000000000..886033792 --- /dev/null +++ b/gnuradio-examples/python/digital-bert/Makefile.am @@ -0,0 +1,33 @@ +# +# Copyright 2008 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 3, 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. +# + +include $(top_srcdir)/Makefile.common + +EXTRA_DIST = \ + benchmark_rx.py \ + benchmark_tx.py \ + receive_path.py \ + transmit_path.py + +ourdatadir = $(exampledir)/digital-bert +ourdata_DATA = $(EXTRA_DIST) + +MOSTLYCLEANFILES = *.pyc *.pyo *~ diff --git a/gnuradio-examples/python/digital-bert/README b/gnuradio-examples/python/digital-bert/README new file mode 100644 index 000000000..74d85b703 --- /dev/null +++ b/gnuradio-examples/python/digital-bert/README @@ -0,0 +1,62 @@ +BERT testing example scripts + + +benchmark_tx.py + +This sets up a BPSK transmitter that is modulated with a pseudorandom +sequence of bits. The PN code is generated by sending an all 1s +sequence through a 7-bit scrambler. The transmitter performs the BPSK +modulation, then passes the complex baseband waveform through a +root-raised-cosine filter and onto the USRP. + +The --sps parameter controls how many baseband samples per symbol +are created and passed through the RRC filter, prior to going to the +USRP over the USB for interpolation to the final DAC rate. + +The baseband bit rate is controlled by -r or --rate. This value, when +multiplied by the --sps parameter, must result in valid interpolation +rate for the USRP. For example, if the baseband rate is 250k bits/sec, +and the samples per symbol is 4, then the final rate is 1M samples/sec, +which results in an interpolation rate of 128. The valid interpolation +rates for the USRP are multiples of 4 between 16 and 512. + +Finally, the RRC excess bandwidth may be specified by --excess-bw. +(See ./benchmark_tx.py -h for additional parameters.) + + +benchmark_rx.py + +This sets up a BPSK receiver to demodulate the received waveform. It +accepts a similar set of parameters as the transmitter, except that one +specifies the USRP decimation rate desired. The resulting sample stream +rate must be an integral number of baseband symbols. For example, the +parameters corresponding to the above transmitter would be to use a +decimation rate of 8 (32 sps), 16 (16 sps), 32 (8 sps), 64, (4 sps), or +128 (2 sps). The lower the USRP decimation, the more CPU is required to +demodulate the signal, so not all valid decimation rates will work. + +The baseband signal from the USRP is first passed through an AGC to +establish an average power of 1.0. It is then passed through a matched +filter (another RRC), a Costas phase-locked loop, and an M&M bit timing +recovery loop. The resulting constellation has an SNR estimation probe +attached, and is then sliced into a bit stream. + +The recovered bits are then passed through a 7-bit descrambler. If +there are no channel errors, the all 1s sequence is recovered. In the +event of a channel error, there will be a 0 in the bit stream for each +feedback tap in the descrambler. In this case, the CCSDS descrambler is +using 3 feedback taps. + +Finally, the signal is passed into a bit density measurement probe. The +channel BER is measured by dividing the 0s density by three. This +measurement is inaccurate at high BER rates (>10%) as the error 0s +begin to overlap. + +The benchmark script will, once per second, output the mean symbol power, +the average noise variance, the estimated SNR, and the average BER. + +NOTE: The particular SNR estimator used is inaccurate below about 7dB, +and will report erroneously high values even for random noise. + +There are a variety of Costas and M&M loop parameters one can adjust. +See ./benchmark_rx.py -h for the full set. diff --git a/gnuradio-examples/python/digital-bert/benchmark_rx.py b/gnuradio-examples/python/digital-bert/benchmark_rx.py new file mode 100755 index 000000000..96bf7fa12 --- /dev/null +++ b/gnuradio-examples/python/digital-bert/benchmark_rx.py @@ -0,0 +1,164 @@ +#!/usr/bin/env python +# +# Copyright 2008 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 3, 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, usrp +from optparse import OptionParser +from gnuradio.eng_option import eng_option +from receive_path import receive_path +import gnuradio.gr.gr_threading as _threading +import sys, time + +n2s = eng_notation.num_to_str + +class status_thread(_threading.Thread): + def __init__(self, tb): + _threading.Thread.__init__(self) + self.setDaemon(1) + self.tb = tb + self.done = False + self.start() + + def run(self): + while not self.done: + print "Magnitude: %5.3f Variance: %5.3f Estimated SNR: %4.1f dB BER: %g" % ( + tb.mag(), tb.var(), tb.snr(), tb.ber()) + try: + time.sleep(1.0) + except KeyboardInterrupt: + self.done = True + +class rx_bpsk_block(gr.top_block): + def __init__(self, options): + + gr.top_block.__init__(self, "rx_mpsk") + + print "USRP decimation rate", options.decim_rate + + # Create a USRP source at desired board, sample rate, frequency, and gain + self._setup_usrp(options.which, + options.decim_rate, + options.rx_subdev_spec, + options.freq, + options.gain) + + # Create the BERT receiver + if_rate = self._usrp.adc_rate()/options.decim_rate + self._receiver = receive_path(if_rate, + options.rate, + options.excess_bw, + options.costas_alpha, + options.costas_beta, + options.costas_max, + options.mm_gain_mu, + options.mm_gain_omega, + options.mm_omega_limit) + + self.connect(self._usrp, self._receiver) + + + def _setup_usrp(self, which, decim, subdev_spec, freq, gain): + self._usrp = usrp.source_c(which=which, decim_rate=decim) + if subdev_spec is None: + subdev_spec = usrp.pick_rx_subdevice(self._usrp) + self._subdev = usrp.selected_subdev(self._usrp, subdev_spec) + mux = usrp.determine_rx_mux_value(self._usrp, subdev_spec) + self._usrp.set_mux(mux) + tr = self._usrp.tune(0, self._subdev, freq) + if not (tr): + print "Failed to tune to center frequency!" + else: + print "Center frequency:", n2s(freq) + if gain is None: + g = self._subdev.gain_range(); + gain = float(g[0]+g[1])/2.0 + self._subdev.set_gain(gain) + print "RX d'board:", self._subdev.side_and_name() + + def snr(self): + return self._receiver.snr() + + def mag(self): + return self._receiver.signal_mean() + + def var(self): + return self._receiver.noise_variance() + + def ber(self): + return self._receiver.ber() + +def get_options(): + parser = OptionParser(option_class=eng_option) + parser.add_option("-w", "--which", type="int", default=0, + help="select which USRP (0, 1, ...) (default is %default)", + metavar="NUM") + parser.add_option("-R", "--rx-subdev-spec", type="subdev", default=None, + help="select USRP Rx side A or B (default=first one with a daughterboard)") + parser.add_option("-f", "--freq", type="eng_float", default=None, + help="set frequency to FREQ", metavar="FREQ") + parser.add_option("-g", "--gain", type="eng_float", default=None, + help="set Rx gain (default is mid-point)") + parser.add_option("-r", "--rate", type="eng_float", default=250e3, + help="Select modulation symbol rate (default=%default)") + parser.add_option("-d", "--decim-rate", type="int", default=8, + help="Select USRP decimation rate (default=%default)") + parser.add_option("", "--excess-bw", type="eng_float", default=0.35, + help="Select RRC excess bandwidth (default=%default)") + parser.add_option("", "--costas-alpha", type="eng_float", default=0.2, + help="set Costas loop 1st order gain, (default=%default)") + parser.add_option("", "--costas-beta", type="eng_float", default=0.01, + help="set Costas loop 2nd order gain, (default=%default)") + parser.add_option("", "--costas-max", type="eng_float", default=0.005, + help="set Costas loop max freq (rad/sample) (default=%default)") + parser.add_option("", "--mm-gain-mu", type="eng_float", default=0.2, + help="set M&M loop 1st order gain, (default=%default)") + parser.add_option("", "--mm-gain-omega", type="eng_float", default=0.01, + help="set M&M loop 2nd order gain, (default=%default)") + parser.add_option("", "--mm-omega-limit", type="eng_float", default=0.005, + help="set M&M max timing error, (default=%default)") + + + (options, args) = parser.parse_args() + if len(args) != 0: + parser.print_help() + sys.exit(1) + + if options.freq == None: + print "You must supply a frequency with -f or --freq" + sys.exit(1) + + return (options, args) + + +if __name__ == "__main__": + (options, args) = get_options() + + tb = rx_bpsk_block(options) + + print "\n*** SNR estimator is inaccurate below about 7dB" + print "*** BER estimator is inaccurate above about 10%\n" + updater = status_thread(tb) + + try: + tb.run() + except KeyboardInterrupt: + updater.done = True + updater = None diff --git a/gnuradio-examples/python/digital-bert/benchmark_tx.py b/gnuradio-examples/python/digital-bert/benchmark_tx.py new file mode 100755 index 000000000..1778a74b5 --- /dev/null +++ b/gnuradio-examples/python/digital-bert/benchmark_tx.py @@ -0,0 +1,111 @@ +#!/usr/bin/env python +# +# Copyright 2008 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 3, 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, usrp +from gnuradio.eng_option import eng_option +from optparse import OptionParser +from transmit_path import transmit_path +import sys + +_dac_rate = 128e6 + +n2s = eng_notation.num_to_str + +class tx_bpsk_block(gr.top_block): + def __init__(self, options): + gr.top_block.__init__(self, "tx_mpsk") + + self._transmitter = transmit_path(options.sps, + options.excess_bw, + options.amplitude) + + if_rate = options.rate*options.sps + interp = int(_dac_rate/if_rate) + + print "Modulation:", n2s(options.rate), "bits/sec" + print "TX IF rate:", n2s(if_rate), "samples/sec" + print "USRP interpolation:", interp + print "DAC amplitude:", options.amplitude + + self._setup_usrp(options.which, + interp, + options.tx_subdev_spec, + options.freq) + + self.connect(self._transmitter, self._usrp) + + + def _setup_usrp(self, which, interp, subdev_spec, freq): + self._usrp = usrp.sink_c(which=which, interp_rate=interp) + if subdev_spec is None: + subdev_spec = usrp.pick_tx_subdevice(self._usrp) + self._usrp.set_mux(usrp.determine_tx_mux_value(self._usrp, subdev_spec)) + self._subdev = usrp.selected_subdev(self._usrp, subdev_spec) + tr = usrp.tune(self._usrp, self._subdev._which, self._subdev, freq) + if not (tr): + print "Failed to tune to center frequency!" + else: + print "Center frequency:", n2s(freq) + gain = float(self._subdev.gain_range()[1]) # Max TX gain + self._subdev.set_gain(gain) + self._subdev.set_enable(True) + print "TX d'board:", self._subdev.side_and_name() + + +def get_options(): + parser = OptionParser(option_class=eng_option) + parser.add_option("-w", "--which", type="int", default=0, + help="select which USRP (0, 1, ...) default is %default", + metavar="NUM") + parser.add_option("-T", "--tx-subdev-spec", type="subdev", default=None, + help="select USRP Tx side A or B (default=first one with a daughterboard)") + parser.add_option("-f", "--freq", type="eng_float", default=None, + help="set frequency to FREQ", metavar="FREQ") + parser.add_option("-a", "--amplitude", type="eng_float", default=2000, + help="set Tx amplitude (0-32767) (default=%default)") + parser.add_option("-r", "--rate", type="eng_float", default=250e3, + help="Select modulation symbol rate (default=%default)") + parser.add_option("", "--sps", type="int", default=2, + help="Select samples per symbol (default=%default)") + parser.add_option("", "--excess-bw", type="eng_float", default=0.35, + help="Select RRC excess bandwidth (default=%default)") + + (options, args) = parser.parse_args() + if len(args) != 0: + parser.print_help() + sys.exit(1) + + if options.freq == None: + print "Must supply frequency as -f or --freq" + sys.exit(1) + + return (options, args) + +if __name__ == "__main__": + (options, args) = get_options() + + tb = tx_bpsk_block(options) + + try: + tb.run() + except KeyboardInterrupt: + pass diff --git a/gnuradio-examples/python/digital-bert/receive_path.py b/gnuradio-examples/python/digital-bert/receive_path.py new file mode 100644 index 000000000..84ff32916 --- /dev/null +++ b/gnuradio-examples/python/digital-bert/receive_path.py @@ -0,0 +1,116 @@ +# +# Copyright 2008 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 3, 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 +import math + +n2s = eng_notation.num_to_str + +class receive_path(gr.hier_block2): + def __init__(self, + if_rate, # Incoming sample rate + symbol_rate, # Original symbol rate + excess_bw, # RRC excess bandwidth, typically 0.35-0.5 + costas_alpha, # Costas loop 1st order gain, typically 0.01-0.2 + costas_beta, # Costas loop 2nd order gain, typically alpha^2/4.0 + costas_max, # Costas loop max frequency offset in radians/sample + mm_gain_mu, # M&M loop 1st order gain, typically 0.001-0.2 + mm_gain_omega, # M&M loop 2nd order gain, typically alpha^2/4.0 + mm_omega_limit, # M&M loop max timing error + ): + + gr.hier_block2.__init__(self, "receive_path", + gr.io_signature(1, 1, gr.sizeof_gr_complex), # Input signature + gr.io_signature(0, 0, 0)) # Output signature + + sps = int(if_rate/symbol_rate) + print "IF sample rate:", n2s(if_rate) + print "Symbol rate:", n2s(symbol_rate) + print "Samples/symbol:", sps + print "RRC bandwidth:", excess_bw + + # Create AGC to scale input to unity + self._agc = gr.agc_cc(1e-5, 1.0, 1.0, 1.0) + + # Create RRC with specified excess bandwidth + taps = gr.firdes.root_raised_cosine(1.0, # Gain + sps, # Sampling rate + 1.0, # Symbol rate + excess_bw, # Roll-off factor + 11*sps) # Number of taps + + self._rrc = gr.fir_filter_ccf(1, taps) + + # Create a Costas loop frequency/phase recovery block + + print "Costas alpha:", costas_alpha + print "Costas beta:", costas_beta + print "Costas max:", costas_max + + self._costas = gr.costas_loop_cc(costas_alpha, # PLL first order gain + costas_beta, # PLL second order gain + costas_max, # Max frequency offset rad/sample + -costas_max, # Min frequency offset rad/sample + 2) # BPSK + + # Create a M&M bit synchronization retiming block + mm_mu = 0.5 + mm_omega = sps + + print "MM gain mu:", mm_gain_mu + print "MM gain omega:", mm_gain_omega + print "MM omega limit:", mm_omega_limit + + self._mm = gr.clock_recovery_mm_cc(mm_omega, # Initial samples/symbol + mm_gain_omega, # Second order gain + mm_mu, # Initial symbol phase + mm_gain_mu, # First order gain + mm_omega_limit) # Maximum timing offset + + # Add an SNR probe on the demodulated constellation + self._snr_probe = gr.probe_mpsk_snr_c(10.0/symbol_rate) + self.connect(self._mm, self._snr_probe) + + # Slice the resulting constellation into bits. + # Get inphase channel and make decision about 0 + self._c2r = gr.complex_to_real() + self._slicer = gr.binary_slicer_fb() + + # Descramble BERT sequence. A channel error will create 3 incorrect bits + self._descrambler = gr.descrambler_bb(0x8A, 0x7F, 7) # CCSDS 7-bit descrambler + + # Measure BER by the density of 0s in the stream + self._ber = gr.probe_density_b(1.0/symbol_rate) + + self.connect(self, self._agc, self._rrc, self._costas, self._mm, + self._c2r, self._slicer, self._descrambler, self._ber) + + def snr(self): + return self._snr_probe.snr() + + def signal_mean(self): + return self._snr_probe.signal_mean() + + def noise_variance(self): + return self._snr_probe.noise_variance() + + def ber(self): + return (1.0-self._ber.density())/3.0 diff --git a/gnuradio-examples/python/digital-bert/transmit_path.py b/gnuradio-examples/python/digital-bert/transmit_path.py new file mode 100644 index 000000000..47b6ab229 --- /dev/null +++ b/gnuradio-examples/python/digital-bert/transmit_path.py @@ -0,0 +1,62 @@ +#!/usr/bin/env python +# +# Copyright 2008 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 3, 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 +from math import pi, log10 +import cmath + +class transmit_path(gr.hier_block2): + """ + This transmits a BERT sequence of bits using filtered BPSK and + outputs the complex baseband waveform. + """ + def __init__(self, + sps, # Samples per symbol + excess_bw, # RRC filter excess bandwidth (typically 0.35-0.5) + amplitude # DAC output level, 0-32767, typically 2000-8000 + ): + + gr.hier_block2.__init__(self, "transmit_path", + gr.io_signature(0, 0, 0), # Input signature + gr.io_signature(1, 1, gr.sizeof_gr_complex)) # Output signature + + # Create BERT data bit stream + self._bits = gr.vector_source_b([1,], True) # Infinite stream of ones + self._scrambler = gr.scrambler_bb(0x8A, 0x7F, 7) # CCSDS 7-bit scrambler + + # Map to constellation + self._constellation = [-1+0j, 1+0j] + self._mapper = gr.chunks_to_symbols_bc(self._constellation) + + # Create RRC with specified excess bandwidth + taps = gr.firdes.root_raised_cosine(sps*amplitude, # Gain + sps, # Sampling rate + 1.0, # Symbol rate + excess_bw, # Roll-off factor + 11*sps) # Number of taps + + self._rrc = gr.interp_fir_filter_ccf(sps, # Interpolation rate + taps) # FIR taps + + # Wire block inputs and outputs + self.connect(self._bits, self._scrambler, self._mapper, self._rrc, self) + |