TI DSP的 FFT与IFFT库函数的使用测试

ToneChip

2022-03-25 14:39:47

关注

VIPTI DSP的 FFT与IFFT库函数的使用测试

/** * Copyright (C) 2013 Guangzhou Tronlong Electronic Technology Co., Ltd. - www.tronlong.com * * @file dsplib_fft.c * * @brief Example application main file. * This application will run FFT performance test. * * @author Tronlong <support@tron.

开通可读全文 ·2.6k 人阅读

最低0.47元/天解锁文章

/**
 * Copyright (C) 2013 Guangzhou Tronlong Electronic Technology Co., Ltd. - www.tronlong.com
 *
 * @file dsplib_fft.c
 *
 * @brief Example application main file.
 * This application will run FFT performance test.
 *
 * @author Tronlong <support@tronlong.com>
 *
 * @version V1.0
 *
 * @date 2019-09-12
 *
 **/

/* C Standard Header files */
#include <stdio.h>
#include <stdlib.h>

/* Compiler Header files */
#include <c6x.h>

/* CSL Header files */
#include <ti/csl/csl_cacheAux.h>

/* DSPLIB Header files */
#include <inc/dsplib.h>

/* MATHLIB Header files */
#include <inc/mathlib.h>

/* Local Header files */
#include "system/platform.h"

/* MAX length of FFT in complex samples */
#define MAXN    4096

/* Sampling Rate */
#define Fs      1000.0

#define PI      3.141592654

/* bit reverse table containing 64 entries  */
unsigned char brev[64] =
{
    0x00, 0x20, 0x10, 0x30, 0x08, 0x28, 0x18, 0x38,
    0x04, 0x24, 0x14, 0x34, 0x0c, 0x2c, 0x1c, 0x3c,
    0x02, 0x22, 0x12, 0x32, 0x0a, 0x2a, 0x1a, 0x3a,
    0x06, 0x26, 0x16, 0x36, 0x0e, 0x2e, 0x1e, 0x3e,
    0x01, 0x21, 0x11, 0x31, 0x09, 0x29, 0x19, 0x39,
    0x05, 0x25, 0x15, 0x35, 0x0d, 0x2d, 0x1d, 0x3d,
    0x03, 0x23, 0x13, 0x33, 0x0b, 0x2b, 0x1b, 0x3b,
    0x07, 0x27, 0x17, 0x37, 0x0f, 0x2f, 0x1f, 0x3f
};

/**
 * @details generate specialized sequence of twiddle factors for FFT
 *
 * @param w pointer to complex twiddle factor
 *
 * @param n length of FFT in complex samples
 *
 * @return NULL
 */
void tw_fft_gen(float *w, int n)
{
    int i, j, k;

    for(j = 1, k = 0; j <= n >> 2; j = j << 2) {
        for(i = 0; i < n >> 2; i += j) {
            w[k]     = (float)sinsp(2 * PI * i / n);