FFTW3 вычисляет взаимную корреляцию в том же сигнале

В настоящее время я создаю код C, который принимает в качестве входных данных wav файл (в частности, только один канал оригинала wav файл), и он выполняет кратковременное преобразование Фурье.
Основная часть кода это:

stft_data = (fftw_complex*)(fftw_malloc(sizeof(fftw_complex)*windowSize));

fft_result= (fftw_complex*)(fftw_malloc(sizeof(fftw_complex)*windowSize));

storage = (fftw_complex*)(fftw_malloc(sizeof(fftw_complex)*storage_capacity));

//define the fftw plane
fftw_plan plan_forward;
plan_forward = fftw_plan_dft_1d(windowSize, stft_data, fft_result, FFTW_FORWARD, FFTW_ESTIMATE);

//integer indexes
int i,counter ;
counter = 0 ;
//create a Hamming window
double hamming_result[windowSize];
hamming(windowSize, hamming_result);

//implement the stft position indexes
int chunkPosition = 0; //actual chunk position
int readIndex ; //read the index of the wav file

while (chunkPosition < wav_length ){
//read the window
for(i=0; i<windowSize; i++){

readIndex = chunkPosition + i;

if (readIndex < wav_length){
stft_data[i] = wav_data[readIndex]*hamming_result[i]*_Complex_I  + 0.0*I;
}
else{
//if we are beyond the wav_length
stft_data[i] = 0.0*_Complex_I + 0.0*I;//padding
break;
}
}
//compute the fft
fftw_execute(plan_forward);
//store the stft in a data structure
for (i=0; i<windowSize;i++)
{
//printf("RE: %.2f  IM: %.2f\n", creal(fft_result[i]),cimag(fft_result[i]));
storage[counter] = creal(fft_result[i]) + cimag(fft_result[i]);
counter+=1;
}

//update indexes
chunkPosition += hop_size;
printf("Chunk Position %d\n", chunkPosition);
printf("Counter position %d\n", counter);
printf("Fourier transform done\n");

}

Как только БПФ вычислено для выбранного окна, я сохраняю реальную и мнимую части БПФ в storage переменная.

После этого я хотел бы вычислить взаимную корреляцию между точками данных в каждом из N окон, которые у меня есть в конце.
В качестве примера я хотел бы вычислить корреляцию между первой точкой данных первого окна ( storage[0] ) с первым элементом второго окна (storage[windowSize+1]).
Однако я сталкиваюсь с некоторыми проблемами, и у меня нет разумных ценностей. Согласно тому, что я изучал, корреляция в пространстве Фурье это просто сложное умножение между двумя членами Фурье. Таким образом,
что я делаю, это что-то вроде:

correlation = storage[0]*conj(storage[windowSize+1]);

Однако я получил очень большие значения, что заставляет меня задуматься, действительно ли я вычисляю корреляцию.

Где я не прав?
Как мне масштабировать результаты корреляции?
Как мне вычислить корреляцию со значениями Фурье?
и затем, как я должен построить значения Фурье, которые я получаю из расчетов FFTW3? мне сдвинуть все значения или они уже сдвинуты?

Спасибо большое