3 * Copyright 2002 Free Software Foundation, Inc.
5 * This file is part of GNU Radio
7 * GNU Radio is free software; you can redistribute it and/or modify
8 * it under the terms of the GNU General Public License as published by
9 * the Free Software Foundation; either version 3, or (at your option)
12 * GNU Radio is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 * GNU General Public License for more details.
17 * You should have received a copy of the GNU General Public License
18 * along with GNU Radio; see the file COPYING. If not, write to
19 * the Free Software Foundation, Inc., 51 Franklin Street,
20 * Boston, MA 02110-1301, USA.
25 #include <gri_goertzel.h>
27 gri_goertzel::gri_goertzel(int rate, int len, float freq)
32 float w = 2.0*M_PI*freq/rate;
33 d_wr = 2.0*std::cos(w);
40 gr_complex gri_goertzel::batch(float *in)
45 for(int i = 0; i < d_len; i++)
51 void gri_goertzel::input(const float &input)
53 float y = input + d_wr*d_d1 - d_d2;
59 gr_complex gri_goertzel::output()
61 gr_complex out((0.5*d_wr*d_d1-d_d2)/d_len, (d_wi*d_d1)/d_len);