Come si imposta un buffer quando si esegue un FFT utilizzando il framework Accelerate?

Question

Utilizzo il framework Accelerate per eseguire una trasformazione Fast Fourier (FFT) e sto cercando di trovare un modo per creare un buffer da utilizzare con una lunghezza di 1024. Ho accesso al picco e al picco medio di un segnale su cui voglio fare la FFT.

Qualcuno può aiutarmi o darmi qualche suggerimento per farlo?

Answer 1

Apple ha alcuni esempi su come impostare FFT nel loro vDSP Programming Guide. Dovresti anche dare un'occhiata all'applicazione di esempio vDSP Examples. Mentre per il Mac, questo codice dovrebbe essere tradotto direttamente su iOS.

Recentemente ho bisogno di fare una semplice FFT di una forma d'onda di ingresso 64 intero, per cui ho usato il codice seguente:

static FFTSetupD fft_weights; 
static DSPDoubleSplitComplex input; 
static double *magnitudes; 

+ (void)initialize 
{ 
    /* Setup weights (twiddle factors) */ 
    fft_weights = vDSP_create_fftsetupD(6, kFFTRadix2); 

    /* Allocate memory to store split-complex input and output data */ 
    input.realp = (double *)malloc(64 * sizeof(double)); 
    input.imagp = (double *)malloc(64 * sizeof(double)); 
    magnitudes = (double *)malloc(64 * sizeof(double)); 
} 

- (CGFloat)performAcceleratedFastFourierTransformAndReturnMaximumAmplitudeForArray:(NSUInteger *)waveformArray; 
{ 
    for (NSUInteger currentInputSampleIndex = 0; currentInputSampleIndex < 64; currentInputSampleIndex++) 
    { 
     input.realp[currentInputSampleIndex] = (double)waveformArray[currentInputSampleIndex]; 
     input.imagp[currentInputSampleIndex] = 0.0f; 
    } 

    /* 1D in-place complex FFT */ 
    vDSP_fft_zipD(fft_weights, &input, 1, 6, FFT_FORWARD); 

    input.realp[0] = 0.0; 
    input.imagp[0] = 0.0; 

    // Get magnitudes 
    vDSP_zvmagsD(&input, 1, magnitudes, 1, 64); 

    // Extract the maximum value and its index 
    double fftMax = 0.0; 
    vDSP_maxmgvD(magnitudes, 1, &fftMax, 64); 

    return sqrt(fftMax); 
} 

Come si può vedere, ho usato solo i valori reali in questo FFT per impostare su i buffer di input, eseguita la FFT, e quindi leggere le grandezze.