ffmpeg做fft

最新推荐文章于 2023-02-22 13:37:47 发布

qq229873466

最新推荐文章于 2023-02-22 13:37:47 发布

阅读量857

点赞数

分类专栏： ffmpeg

本文链接：https://blog.csdn.net/qq229873466/article/details/78093806

版权

ffmpeg 专栏收录该内容

2 篇文章 1 订阅

订阅专栏

#define inline _inline

#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#include "libavcodec/fft.h"

int main()
{
    float a[8] = { 0, 1, 2, 3, 4, 5, 6, 7 };

    FFTComplex *tab;
    FFTContext s;

    int i;
    int do_inverse = 0;
    int fft_nbits = 3, fft_size;

    fft_size = 1 << fft_nbits;
    tab = av_malloc_array(fft_size, sizeof(FFTComplex));

    if (!(tab))
        goto cleanup;

    ff_fft_init(&s, fft_nbits, do_inverse);

    for (i = 0; i < fft_size; i++) {
        tab[i].re = a[i];
        tab[i].im = 0;
    }

    s.fft_permute(&s, tab);
    s.fft_calc(&s, tab);
    ff_fft_end(&s);

    for (i = 0; i < fft_size; i++)
    {
        printf("%f %fi\n", tab[i].re, tab[i].im);
    }
    
cleanup:
    av_free(tab);

    return 0;
}