Построение спектров sin и cos сигналов

#include "stdafx.h" 
#include "iostream"
#include "cmath"
 
using namespace std;
 
const float PI = 3.1415926535897932384626433832795;
const int N = 32;
 
int _tmain(int argc, _TCHAR* argv[])
{
    setlocale(LC_ALL, "Russian");
    int s[N] = { 3, 5, 26, -23, -19, -18, 7, -8, 21, -27, 12, 13, 26, 4, 11, -11, -6, -21, 17, -2, 5, 12, 29, 4, -14, -13, -23, 7, 10, -5, 14, -2 };
    float mCos[N];
    float mSin[N];
 
    float sumCos = 0;
    float sumSin = 0;
 
    float ODPF = 0;
 
    float fCos, fSin, fCos2, fSin2;
 
    float amplitude;
    float fase;
 
    cout << endl;
    cout << "ДПФ: " << endl;
    for (int k = 0; k < N; k++)
    {
        for (int n = 0; n < N; n++)
        {
            fCos = (s[n] * (cos((2 * PI * k * n) / N)));
            fSin = (s[n] * (-sin((2 * PI * k * n) / N)));
 
            sumCos += fCos;
            sumSin += fSin;
 
            mCos[k] = sumCos;
            mSin[k] = sumSin;
        }
        cout << "     " << k << " -> " << "(" << sumCos << " | " << sumSin << "i" << ")" << endl;
 
        sumCos = 0;
        sumSin = 0;
    }
 
    cout << endl;
    cout << "ОДПФ: " << endl;
    for (int k = 0; k < N; k++)
    {
        for (int n = 0; n < N; n++)
        {
            fCos2 = (mCos[n] / N) * cos((2 * PI * n * k) / N) - (mSin[n] / N) * sin((2 * PI * n * k) / N);
 
            ODPF += fCos2;  
        }
        cout << "     " << k << " -> " << "(" << ODPF << ")" << endl;
 
        ODPF = 0;
    }
 
    cout << endl;
    cout << "Амплитуда:" << endl;
    for (int i = 0; i < N; i++)
    {
        amplitude = sqrt(pow(abs(mCos[i]), 2) + pow(abs(mSin[i]), 2));
        cout << "     " << i << " -> " << amplitude << endl;
    }
 
    cout << endl;
    cout << "Фаза:" << endl;
    for (int i = 0; i < N; i++)
    {
        fase = atan(mSin[i] / mCos[i]);
        cout << "     " << i << " -> " << fase << endl;
    }
 
    cout << endl;
 
    return 0;
}