spectral_irregularity.hpp

Go to the documentation of this file.

#ifndef EDSP_SPECTRAL_IRREGULARITY_HPP
#define EDSP_SPECTRAL_IRREGULARITY_HPP

#include <iterator>
namespace edsp { namespace feature { inline namespace spectral {

    template <typename ForwardIt>
    constexpr auto spectral_irregularity(ForwardIt first, ForwardIt last) {
        using value_type = typename std::iterator_traits<ForwardIt>::value_type;
        value_type square_diff = 0, square_ampl = 0;
        for (auto until = (last - 1); first != until; ++first) {
            const value_type diff = (*first + 1) - (*first);
            square_diff           += diff * diff;
            square_ampl           += (*first) * (*first);
        }
        return square_diff / square_ampl;
    }

}}}


#endif //EDSP_SPECTRAL_IRREGULARITY_HPP