filter: added mmse_fir_interpolator with QA code.

author: Tom Rondeau <trondeau@vt.edu> 2012-06-17 16:44:04 -0400
committer: Tom Rondeau <trondeau@vt.edu> 2012-06-17 16:44:04 -0400
commit: ce3d887a47b47cb2d5351e14066ceb56b2b3c892 (patch)
tree: e27110a4dd06f2c918cef7cc9a3e164010434b51 /gr-filter/include
parent: d647fbbaa07554ad314c9eb2b5c1735b6265d13a (diff)
4 files changed, 158 insertions, 1 deletions
diff --git a/gr-filter/include/filter/CMakeLists.txt b/gr-filter/include/filter/CMakeLists.txt
index 5da171ffa8..3a1b5f150b 100644
--- a/gr-filter/include/filter/CMakeLists.txt
+++ b/gr-filter/include/filter/CMakeLists.txt
@@ -82,6 +82,8 @@ install(FILES
     fft_filter.h
     iir_filter.h
     interpolator_taps.h
+    mmse_fir_interpolator_cc.h
+    mmse_fir_interpolator_ff.h
     pm_remez.h
     polyphase_filterbank.h
     single_pole_iir.h
diff --git a/gr-filter/include/filter/fir_filter.h b/gr-filter/include/filter/fir_filter.h
index 42b4ac2e03..e861112681 100644
--- a/gr-filter/include/filter/fir_filter.h
+++ b/gr-filter/include/filter/fir_filter.h
@@ -42,7 +42,7 @@ namespace gr {
 	std::vector<float> taps() const;
 	unsigned int ntaps() const;
 	
-	inline float filter(const float input[]);
+	float filter(const float input[]);
 	void  filterN(float output[],
 		      const float input[],
 		      unsigned long n);
diff --git a/gr-filter/include/filter/mmse_fir_interpolator_cc.h b/gr-filter/include/filter/mmse_fir_interpolator_cc.h
new file mode 100644
index 0000000000..0436b4a0ff
--- /dev/null
+++ b/gr-filter/include/filter/mmse_fir_interpolator_cc.h
@@ -0,0 +1,79 @@
+/* -*- c++ -*- */
+/*
+ * Copyright 2002,2007,2012 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 _GRI_MMSE_FIR_INTERPOLATOR_CC_H_
+#define _GRI_MMSE_FIR_INTERPOLATOR_CC_H_
+
+#include <filter/api.h>
+#include <filter/fir_filter.h>
+#include <gr_complex.h>
+#include <vector>
+
+namespace gr {
+  namespace filter {
+    
+    /*!
+     * \brief Compute intermediate samples between signal samples x(k*Ts)
+     * \ingroup filter_primitive
+     *
+     * This implements a Mininum Mean Squared Error interpolator with
+     * 8 taps. It is suitable for signals where the bandwidth of
+     * interest B = 1/(4*Ts) Where Ts is the time between samples.
+     *
+     * Although mu, the fractional delay, is specified as a float, it
+     * is actually quantized. 0.0 <= mu <= 1.0. That is, mu is
+     * quantized in the interpolate method to 32nd's of a sample.
+     *
+     * For more information, in the GNU Radio source code, see:
+     * \li gnuradio-core/src/gen_interpolator_taps/README
+     * \li gnuradio-core/src/gen_interpolator_taps/praxis.txt
+     */
+
+    class FILTER_API mmse_fir_interpolator_cc
+    {
+    public:
+      mmse_fir_interpolator_cc();
+      ~mmse_fir_interpolator_cc();
+
+      unsigned ntaps() const;
+      unsigned nsteps() const;
+
+      /*!
+       * \brief compute a single interpolated output value.
+       *
+       * \p input must have ntaps() valid entries and be 8-byte aligned.
+       * input[0] .. input[ntaps() - 1] are referenced to compute the output value.
+       * \throws std::invalid_argument if input is not 8-byte aligned.
+       *
+       * \p mu must be in the range [0, 1] and specifies the fractional delay.
+       *
+       * \returns the interpolated input value.
+       */
+      gr_complex interpolate(const gr_complex input[], float mu) const;
+
+    protected:
+      std::vector<kernel::fir_filter_ccf *> filters;
+    };
+
+  }  /* namespace filter */
+}  /* namespace gr */
+
+#endif /* _MMSE_FIR_INTERPOLATOR_CC_H_ */
diff --git a/gr-filter/include/filter/mmse_fir_interpolator_ff.h b/gr-filter/include/filter/mmse_fir_interpolator_ff.h
new file mode 100644
index 0000000000..4353aa94f1
--- /dev/null
+++ b/gr-filter/include/filter/mmse_fir_interpolator_ff.h
@@ -0,0 +1,76 @@
+/* -*- c++ -*- */
+/*
+ * Copyright 2002,2012 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 _MMSE_FIR_INTERPOLATOR_FF_H_
+#define _MMSE_FIR_INTERPOLATOR_FF_H_
+
+#include <filter/api.h>
+#include <filter/fir_filter.h>
+#include <vector>
+
+namespace gr {
+  namespace filter {
+    
+    /*!
+     * \brief Compute intermediate samples between signal samples x(k*Ts)
+     * \ingroup filter_primitive
+     *
+     * This implements a Mininum Mean Squared Error interpolator with
+     * 8 taps. It is suitable for signals where the bandwidth of
+     * interest B = 1/(4*Ts) Where Ts is the time between samples.
+     *
+     * Although mu, the fractional delay, is specified as a float, it
+     * is actually quantized. 0.0 <= mu <= 1.0. That is, mu is
+     * quantized in the interpolate method to 32nd's of a sample.
+     *
+     * For more information, in the GNU Radio source code, see:
+     * \li gnuradio-core/src/gen_interpolator_taps/README
+     * \li gnuradio-core/src/gen_interpolator_taps/praxis.txt
+     */
+    class FILTER_API mmse_fir_interpolator_ff
+    {
+    public:
+      mmse_fir_interpolator_ff();
+      ~mmse_fir_interpolator_ff();
+
+      unsigned ntaps() const;
+      unsigned nsteps() const;
+
+      /*!
+       * \brief compute a single interpolated output value.
+       * \p input must have ntaps() valid entries.
+       * input[0] .. input[ntaps() - 1] are referenced to compute the output value.
+       *
+       * \p mu must be in the range [0, 1] and specifies the fractional delay.
+       *
+       * \returns the interpolated input value.
+       */
+      float interpolate(const float input[], float mu) const;
+
+    protected:
+      std::vector<kernel::fir_filter_fff*> filters;
+    };
+
+  }  /* namespace filter */
+}  /* namespace gr */
+
+#endif /* _MMSE_FIR_INTERPOLATOR_FF_H_ */
author	Tom Rondeau <trondeau@vt.edu>	2012-06-17 16:44:04 -0400
committer	Tom Rondeau <trondeau@vt.edu>	2012-06-17 16:44:04 -0400
commit	ce3d887a47b47cb2d5351e14066ceb56b2b3c892 (patch)
tree	e27110a4dd06f2c918cef7cc9a3e164010434b51 /gr-filter/include
parent	d647fbbaa07554ad314c9eb2b5c1735b6265d13a (diff)