C++: comment transformer un tableau en un autre tableau de taille différente en utilisant fftw3

Question

J'ai besoin de transformer un tableau avec des valeurs entières en domaine fréquentiel (afin de le multiplier par un autre plus tard). Le tableau de sortie doit avoir la taille de 44100 mais le tableau d'entrée variera. Je pense que fftw3 est un bon outil pour cela. Mais comment puis-je créer un «plan» avec différentes tailles de tableaux pour l'entrée et la sortie?

est ici la fonction que j'ai écrit:

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; 
} 

merci beaucoup pour répondre à

Answer 1

La transformée de Fourier discrète sera toujours produire une sortie avec la même taille que son entrée, en tant que principe mathématique de base l'algorithme, et cela est reflété dans FFTW en prenant seulement dans un paramètre de taille. En effet, pour un tableau de N valeurs discrètes dans le domaine temporel, il y a exactement N fréquences uniques possibles en jeu, de 0 à N-1 cycles par unité de temps (cette unité étant la durée de votre entrée).

Vous pouvez augmenter ou réduire la taille de votre entrée à 44100 valeurs, en étirant essentiellement l'entrée et en interpolant ses valeurs pour créer une nouvelle matrice qui correspond à la longueur requise, puis effectuer la transformée de Fourier. Je vous recommande d'effectuer cet étirement sur le domaine de temps d'entrée et non sur le domaine fréquentiel pour éviter les artefacts indésirables.