2019-02-04 21:34:12 +01:00
|
|
|
/*
|
|
|
|
Strawberry Music Player
|
|
|
|
This file was part of Clementine.
|
|
|
|
Copyright 2004, Melchior FRANZ <mfranz@kde.org>
|
|
|
|
Copyright 2010, 2014, John Maguire <john.maguire@gmail.com>
|
|
|
|
Copyright 2014, Krzysztof Sobiecki <sobkas@gmail.com>
|
|
|
|
Copyright 2017, Santiago Gil
|
|
|
|
|
|
|
|
Strawberry 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 of the License, or
|
|
|
|
(at your option) any later version.
|
|
|
|
|
|
|
|
Strawberry 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 Strawberry. If not, see <http://www.gnu.org/licenses/>.
|
|
|
|
*/
|
2018-02-27 18:06:05 +01:00
|
|
|
|
|
|
|
#ifndef FHT_H
|
|
|
|
#define FHT_H
|
|
|
|
|
2019-02-04 21:34:12 +01:00
|
|
|
#include <QVector>
|
|
|
|
|
2018-02-27 18:06:05 +01:00
|
|
|
/**
|
|
|
|
* Implementation of the Hartley Transform after Bracewell's discrete
|
|
|
|
* algorithm. The algorithm is subject to US patent No. 4,646,256 (1987)
|
|
|
|
* but was put into public domain by the Board of Trustees of Stanford
|
|
|
|
* University in 1994 and is now freely available[1].
|
|
|
|
*
|
|
|
|
* [1] Computer in Physics, Vol. 9, No. 4, Jul/Aug 1995 pp 373-379
|
|
|
|
*/
|
|
|
|
class FHT {
|
2019-02-04 21:34:12 +01:00
|
|
|
const int num_;
|
|
|
|
const int exp2_;
|
|
|
|
|
|
|
|
QVector<float> buf_vector_;
|
|
|
|
QVector<float> tab_vector_;
|
|
|
|
QVector<int> log_vector_;
|
|
|
|
|
2021-06-12 20:53:23 +02:00
|
|
|
float *buf_();
|
|
|
|
float *tab_();
|
|
|
|
int *log_();
|
2018-02-27 18:06:05 +01:00
|
|
|
|
|
|
|
/**
|
|
|
|
* Create a table of "cas" (cosine and sine) values.
|
|
|
|
* Has only to be done in the constructor and saves from
|
|
|
|
* calculating the same values over and over while transforming.
|
|
|
|
*/
|
|
|
|
void makeCasTable();
|
|
|
|
|
|
|
|
/**
|
|
|
|
* Recursive in-place Hartley transform. For internal use only!
|
|
|
|
*/
|
|
|
|
void _transform(float*, int, int);
|
|
|
|
|
|
|
|
public:
|
|
|
|
/**
|
|
|
|
* Prepare transform for data sets with @f$2^n@f$ numbers, whereby @f$n@f$
|
|
|
|
* should be at least 3. Values of more than 3 need a trigonometry table.
|
|
|
|
* @see makeCasTable()
|
|
|
|
*/
|
2021-09-12 21:24:22 +02:00
|
|
|
explicit FHT(uint);
|
2018-02-27 18:06:05 +01:00
|
|
|
|
|
|
|
~FHT();
|
2019-02-04 21:34:12 +01:00
|
|
|
int sizeExp() const;
|
|
|
|
int size() const;
|
2021-06-22 13:45:29 +02:00
|
|
|
void scale(float*, float) const;
|
2018-02-27 18:06:05 +01:00
|
|
|
|
|
|
|
/**
|
|
|
|
* Exponentially Weighted Moving Average (EWMA) filter.
|
|
|
|
* @param d is the filtered data.
|
|
|
|
* @param s is fresh input.
|
|
|
|
* @param w is the weighting factor.
|
|
|
|
*/
|
2021-06-22 13:45:29 +02:00
|
|
|
void ewma(float *d, float *s, float w) const;
|
2018-02-27 18:06:05 +01:00
|
|
|
|
|
|
|
/**
|
|
|
|
* Logarithmic audio spectrum. Maps semi-logarithmic spectrum
|
|
|
|
* to logarithmic frequency scale, interpolates missing values.
|
|
|
|
* A logarithmic index map is calculated at the first run only.
|
|
|
|
* @param p is the input array.
|
|
|
|
* @param out is the spectrum.
|
|
|
|
*/
|
2021-06-12 20:53:23 +02:00
|
|
|
void logSpectrum(float *out, float *p);
|
2018-02-27 18:06:05 +01:00
|
|
|
|
|
|
|
/**
|
|
|
|
* Semi-logarithmic audio spectrum.
|
|
|
|
*/
|
|
|
|
void semiLogSpectrum(float*);
|
|
|
|
|
|
|
|
/**
|
|
|
|
* Fourier spectrum.
|
|
|
|
*/
|
|
|
|
void spectrum(float*);
|
|
|
|
|
|
|
|
/**
|
|
|
|
* Calculates a mathematically correct FFT power spectrum.
|
|
|
|
* If further scaling is applied later, use power2 instead
|
|
|
|
* and factor the 0.5 in the final scaling factor.
|
|
|
|
* @see FHT::power2()
|
|
|
|
*/
|
|
|
|
void power(float*);
|
|
|
|
|
|
|
|
/**
|
|
|
|
* Calculates an FFT power spectrum with doubled values as a
|
|
|
|
* result. The values need to be multiplied by 0.5 to be exact.
|
|
|
|
* Note that you only get @f$2^{n-1}@f$ power values for a data set
|
|
|
|
* of @f$2^n@f$ input values. This is the fastest transform.
|
|
|
|
* @see FHT::power()
|
|
|
|
*/
|
|
|
|
void power2(float*);
|
|
|
|
|
|
|
|
/**
|
|
|
|
* Discrete Hartley transform of data sets with 8 values.
|
|
|
|
*/
|
2021-06-22 13:41:38 +02:00
|
|
|
static void transform8(float*);
|
2018-02-27 18:06:05 +01:00
|
|
|
|
|
|
|
void transform(float*);
|
|
|
|
};
|
|
|
|
|
2019-02-04 21:34:12 +01:00
|
|
|
#endif // FHT_H
|