summaryrefslogtreecommitdiff
path: root/gr-fft/include/fft/goertzel.h
blob: ff37355e64a1d8289e44bfe66b4469dbf67f4139 (plain)
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
/* -*- c++ -*- */
/*
 * Copyright 2006,2011,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 INCLUDED_FFT_IMPL_GOERTZEL_H
#define INCLUDED_FFT_IMPL_GOERTZEL_H

#include <fft/api.h>
#include <gr_types.h>

namespace gr {
  namespace fft {

    /*!
     * \brief Implements Goertzel single-bin DFT calculation
     * \ingroup misc
     */
    class FFT_API goertzel
    {
    public:
      goertzel(){}
      goertzel(int rate, int len, float freq);
      
      void set_params(int rate, int len, float freq);
      
      // Process a input array
      gr_complex batch(float *in);
      
      // Process sample by sample
      void input(const float &in);
      gr_complex output();
      bool ready() const { return d_processed == d_len; }
      
    private:
      float d_d1;
      float d_d2;
      float d_wr;
      float d_wi;
      int   d_len;
      int   d_processed;
    };
    
  } /* namespace fft */
} /* namespace gr */

#endif /* INCLUDED_FFT_IMPL_GOERTZEL_H */