C++ : 어떻게 fftw3을 사용하여 다른 크기의 다른 배열로 푸리에 변환 할 수 있습니까?

Question

정수 값을 가진 배열을 주파수 도메인으로 푸리에 변환해야합니다 (나중에 다른 것과 곱하기 위해). 출력 배열의 크기는 44100이어야하지만 입력 배열은 다양합니다.

나는 fftw3이 그 좋은 도구라고 생각한다. 하지만 입출력을위한 다양한 배열 크기의 '계획'을 어떻게 만듭니 까? 여기

은 내가 쓴 기능입니다 :

fftw_complex* fourier(int* samples, int numberOfSamples){ 

    fftw_complex* input; 
    fftw_complex* output; 
    fftw_plan plan; 
    input = (fftw_complex*)fftw_malloc(sizeof(fftw_complex)*numberOfSamples); 
    output = (fftw_complex*)fftw_malloc(sizeof(fftw_complex)*44100);     //cd-quality 

    plan = fftw_plan_dft_1d(44100, input, output, FFTW_FORWARD, FFTW_ESTIMATE); 
    // here is the problem because ^this 
    // array has a different size than 
    //          ^this array 

    for(size_t index = 0; index < numberOfSamples; index++){ 
     input[index][0] = (double)(samples[index]); 
     input[index][1] = 0; 
    } 
    fftw_execute(plan); 
    fftw_destroy_plan(plan); 
    fftw_free(input); 
    return output; 
} 

감사의 기본적인 수학 전제로,

Answer 1

이산 푸리에 항상 입력과 같은 크기의 출력을 생성 변환 응답에 대한 많은 이 알고리즘은 하나의 크기 매개 변수만을 사용하여 FFTW에 반영됩니다. 이것은 시간 도메인에서 N 개의 이산 값 배열에 대해 단위 시간당 0에서 N-1 사이클 (이 단위는 입력 시간 길이 임)부터 정확히 N 개의 가능한 고유 주파수가 재생되기 때문입니다.

기본적으로 입력을 늘리고 값을 보간하여 원하는 길이에 맞는 새 배열을 만든 다음 푸리에 변환을 수행하여 입력을 44100 값까지 위/아래로 확장 할 수 있습니다. 원하지 않는 아티팩트를 피하기 위해 주파수 영역이 아닌 입력 시간 영역에서이 스트레칭을 수행하는 것이 좋습니다.