Алгоритм Герцеля (Goertzel algorithm)

#include < stdio.h > 
#include < stdlib.h > 
#define _USE_MATH_DEFINES 
#include < math.h >


int main(){
    // исходные данные    
    int N = 8; //точек ДПФ
    int k = 1; //номер спектрального отсчта

    //исходный сигнал
    double s[8] = {3, 2, 1, -1, 1, -2, -3, -2};

    //массив промежуточных результатов
    double v[8] = {0, 0, 0, 0, 0, 0, 0, 0};

    //параметр альфа
    double a = 2*cos(2*M_PI*k/N);

    //поворотный к-т W_N^-k
    double wr =    cos(2*M_PI*k/N); //реальная часть
    double wi =    sin(2*M_PI*k/N); //мнимая часть

    //учет v[-1] = v[-2] = 0
    v[0] = s[0];    
    v[1] = s[1]+a*v[0];

    // итерационный расчет массива v согласно (14)
    for(int i = 2; i < N; i++)
        v[i] = s[i] + a*v[i-1] - v[i-2];

    //реальная и мнимая части спектрального отсчета S(1) согласно (13)
    double SR = v[N-1]*wr-v[N-2];
    double SI = v[N-1]*wi;

    //печать результата
    printf("S(%d) = %.4f%+.4fj\n",k,SR,SI);

    system("Pause");
    return 0;
}