dft.cpp

#include 
#include 
#include 

using namespace std;

const double PI = 3.141592;

vector< complex > DFT（vector< complex >& theData）
{
	// Define the Size of the read in vector
	const int S = theData.size（）;

	// Initalise new vector with size of S
	vector< complex > out（S 0）;
	for（unsigned i=0; （i < S）; i++）
	{
		out[i] = complex（0.0 0.0）;
		for（unsigned j=0; （j < S）; j++）
		{
			out[i] += theData[j] * polar（1.0 - 2 * PI * i * j / S）;
		}
	}

	retur