ladybird/Userland/Libraries/LibDSP/Window.h
kleines Filmröllchen 00e13b5b27 LibDSP: Rename library namespace to DSP
That's the standard naming convention, but I didn't follow it when
originally creating LibDSP and nobody corrected me, so here I am one
year later :^)
2022-07-19 11:17:45 +01:00

71 lines
2.1 KiB
C++

/*
* Copyright (c) 2021, Arne Elster (arne@elster.li)
*
* SPDX-License-Identifier: BSD-2-Clause
*/
#pragma once
#include <AK/FixedArray.h>
#include <AK/Math.h>
namespace DSP {
template<typename T>
class Window final {
public:
template<size_t size>
constexpr static Array<T, size> hamming() { return make_window<size>(calculate_hamming); }
constexpr static FixedArray<T> hamming(size_t size) { return make_window(size, calculate_hamming); }
template<size_t size>
constexpr static Array<T, size> hann() { return make_window<size>(calculate_hann); }
constexpr static FixedArray<T> hann(size_t size) { return make_window(size, calculate_hann); }
template<size_t size>
constexpr static Array<T, size> blackman_harris() { return make_window<size>(calculate_blackman_harris); }
constexpr static FixedArray<T> blackman_harris(size_t size) { return make_window(size, calculate_blackman_harris); }
private:
constexpr static float calculate_hann(size_t index, size_t size)
{
return 0.5f * (1 - AK::cos<float>((2 * AK::Pi<T> * index) / (size - 1)));
}
constexpr static float calculate_hamming(size_t index, size_t size)
{
return 0.54f - 0.46f * AK::cos<float>((2 * AK::Pi<T> * index) / (size - 1));
}
constexpr static float calculate_blackman_harris(size_t index, size_t size)
{
T const a0 = 0.35875;
T const a1 = 0.48829;
T const a2 = 0.14128;
T const a3 = 0.01168;
return a0 - a1 * AK::cos(2 * AK::Pi<T> * index / size) + a2 * AK::cos(4 * AK::Pi<T> * index / size) - a3 * AK::cos(6 * AK::Pi<T> * index / size);
}
template<size_t size>
constexpr static Array<T, size> make_window(auto window_function)
{
Array<T, size> result;
for (size_t i = 0; i < size; i++) {
result[i] = window_function(i, size);
}
return result;
}
constexpr static FixedArray<T> make_window(size_t size, auto window_function)
{
FixedArray<T> result;
result.resize(size);
for (size_t i = 0; i < size; i++) {
result[i] = window_function(i, size);
}
return result;
}
};
}