diff options
| author | Jacob Gilbert <mrjacobagilbert@gmail.com> | 2020-04-24 12:33:15 -0600 |
|---|---|---|
| committer | mormj <34754695+mormj@users.noreply.github.com> | 2020-04-26 11:17:00 -0400 |
| commit | 44c0882682f2c9110f3a880b86faf2c9252a7c58 (patch) | |
| tree | 734445a155aef0bf5afd9189bee3a2948118c274 /gr-fft/lib/window.cc | |
| parent | cfdf6546024944b8cbbf2bd24b01c5efa5a6c237 (diff) | |
gr-fft: added tukey window generation
Diffstat (limited to 'gr-fft/lib/window.cc')
| -rw-r--r-- | gr-fft/lib/window.cc | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/gr-fft/lib/window.cc b/gr-fft/lib/window.cc index 8638d3738f..8b0a07f430 100644 --- a/gr-fft/lib/window.cc +++ b/gr-fft/lib/window.cc @@ -316,6 +316,29 @@ std::vector<float> window::riemann(int ntaps) return taps; } +std::vector<float> window::tukey(int ntaps, float a) +{ + if ((a < 0) || (a > 1)) + throw std::out_of_range("window::tukey: alpha must be between 0 and 1"); + + float N = static_cast<float>(ntaps - 1); + + float aN = a * N; + float p1 = aN / 2.0; + float mid = midn(ntaps); + std::vector<float> taps(ntaps); + for (int i = 0; i < mid; i++) { + if (abs(i) < p1) { + taps[i] = 0.5 * (1.0 - cos((2 * GR_M_PI * i) / (aN))); + taps[ntaps - 1 - i] = taps[i]; + } else { + taps[i] = 1.0; + taps[ntaps - i - 1] = 1.0; + } + } + return taps; +} + std::vector<float> window::build(win_type type, int ntaps, double beta) { switch (type) { |