请选择 进入手机版 | 继续访问电脑版

硬汉嵌入式论坛

 找回密码
 立即注册
查看: 9487|回复: 19
收起左侧

[DSP] 【安富莱DSP教程】第32章 实数FFT的实现

  [复制链接]

740

主题

1326

回帖

3546

积分

管理员

春暖花开

Rank: 9Rank: 9Rank: 9

积分
3546
QQ
发表于 2015-4-17 10:17:07 | 显示全部楼层 |阅读模式
特别说明:完整45期数字信号处理教程,原创高性能示波器代码全开源地址:链接
第32章 实数FFT的实现


    本章主要讲解实数的浮点和定点Q31,Q15的实现。关于这部分的知识点和函数的计算结果上,官方的文档有一些小错误,在章节中会跟大家详细讲述,还有一个要注意的问题,调用实数FFT函数一定要使用CMSIS-DSP V1.4.4及其以上版本,以前的版本有bug。
    本章节使用的复数FFT函数来自ARM官方库的TransformFunctions部分
    32.1 复数FFT
    32.2 复数FFT-基2算法
    32.3 复数FFT-基4算法
    32.4 总结

32.1 实数FFT

32.1.1 描述

    CMSIS DSP库里面包含一个专门用于计算实数序列的FFT库,很多情况下,用户只需要计算实数序列即可。计算同样点数FFT的实数序列要比计算同样点数的虚数序列有速度上的优势。
    快速的rfft算法是基于混合基cfft算法实现的。
    一个N点的实数序列FFT正变换采用下面的步骤实现:

    由上面的框图可以看出,实数序列的FFT是先计算N/2个实数的CFFT,然后再重塑数据进行处理从而获得半个FFT频谱即可(利用了FFT变换后频谱的对称性)。
    一个N点的实数序列FFT逆变换采用下面的步骤实现:

    实数FFT支持浮点,Q31和Q15三种数据类型。

32.2 实数FFT

32.2.1 arm_rfft_fast_f32

函数定义如下:
    void arm_rfft_fast_f32(
      arm_rfft_fast_instance_f32 * S,
      float32_t * p, float32_t * pOut,
      uint8_t ifftFlag)
参数定义:
     [in]  *S          points to an arm_rfft_fast_instance_f32 structure.
     [in]  *p          points to the input buffer.
     [in]  *pOut      points to the output buffer.
     [in]  ifftFlag     RFFT if flag is 0, RIFFT if flag is 1
注意事项:
结构arm_rfft_fast_instance_f32的定义如下(在文件arm_math.h文件):
      typedef struct
      {
          arm_cfft_instance_f32 Sint;      /**< Internal CFFT structure. */
          uint16_t fftLenRFFT;            /**< length of the real sequence */
      float32_t * pTwiddleRFFT;       /**< Twiddle factors real stage  */
      } arm_rfft_fast_instance_f32 ;

    下面通过在开发板上运行函数arm_rfft_fast_f32和arm_cfft_f32计算幅频响应,然后将相应的频率响应结果在Matlab上面绘制出来。
  1. /*
  2. *********************************************************************************************************
  3. *    函 数 名: arm_rfft_fast_f32_app
  4. *    功能说明: 调用函数arm_rfft_fast_f32计算1024点实数序列的幅频响应并跟使用函数arm_cfft_f32计算结果做对比
  5. *    形    参:无
  6. *    返 回 值: 无
  7. *********************************************************************************************************
  8. */
  9. static void arm_rfft_fast_f32_app(void)
  10. {
  11. uint16_t i;
  12. arm_rfft_fast_instance_f32 S;
  13. /* 实数序列FFT长度 */
  14. fftSize = 1024;
  15. /* 正变换 */
  16.     ifftFlag = 0;
  17. /* 初始化结构体S中的参数 */
  18.      arm_rfft_fast_init_f32(&S, fftSize);
  19. /* 按照实部,虚部,实部,虚部..... 的顺序存储数据 */
  20. for(i=0; i<1024; i++)
  21. {
  22. /* 50Hz正弦波,采样率1KHz */
  23. testInput_f32_10khz[i] = 1.2f*arm_sin_f32(2*3.1415926f*50*i/1000)+1;
  24. }
  25. /* 1024点实序列快速FFT */
  26. arm_rfft_fast_f32(&S, testInput_f32_10khz, testOutput_f32_10khz, ifftFlag);
  27. /* 为了方便跟函数arm_cfft_f32计算的结果做对比,这里求解了1024组模值,实际函数arm_rfft_fast_f32
  28.    只求解出了512组  
  29. */
  30.      arm_cmplx_mag_f32(testOutput_f32_10khz, testOutput, fftSize);
  31. /* 串口打印求解的模值 */
  32. for(i=0; i<fftSize; i++)
  33. {
  34. printf("%f\r\n", testOutput[i]);
  35. }
  36. printf("****************************分割线***************************************\r\n");
  37. for(i=0; i<1024; i++)
  38. {
  39. /* 虚部全部置零 */
  40. testInput_f32_10khz[i*2+1] = 0;
  41. /* 50Hz正弦波,采样率1KHz ,作为实部 */
  42. testInput_f32_10khz[i*2] = 1.2f*arm_sin_f32(2*3.1415926f*50*i/1000)+1;
  43. }
  44. arm_cfft_f32(&arm_cfft_sR_f32_len1024, testInput_f32_10khz, ifftFlag, doBitReverse);
  45. /* 求解模值  */
  46.      arm_cmplx_mag_f32(testInput_f32_10khz, testOutput, fftSize);
  47. /* 串口打印求解的模值 */
  48. for(i=0; i<fftSize; i++)
  49. {
  50. printf("%f\r\n", testOutput[i]);
  51. }
  52. }
复制代码
运行如上函数可以通过串口打印出函数arm_rfft_fast_f32和arm_cfft_f32计算的幅频模值,下面通过Matlab绘制波形来对比这两种模值。
    对比前需要先将串口打印出的两组数据加载到Matlab中,arm_rfft_fast_f32的计算结果起名signal,arm_cfft_f32的计算结果起名sampledata,
加载方法在前面的教程中已经讲解过,这里不做赘述了。Matlab中运行的代码如下:
Fs = 1000;              % 采样率
N  = 1024;            % 采样点数

n  = 0:N-1;            % 采样序列
f = n * Fs / N;          %真实的频率

subplot(3,1,1);
plot(f,  signal);        %绘制RFFT结果
title('实数FFT');
xlabel('时间');
ylabel('幅值');

subplot(3,1,2);
plot(f,  sampledata);   %CFFT结果
title('复数FFT');
xlabel('时间');
ylabel('幅值');

Matlab运行结果如下:

    从上面的前512点对比中,我们可以看出两者的计算结果是相符的。这里有一点要特别注意,官方文档中对于函数arm_rfft_fast_f32输出结果的实部,虚部排列顺序说明是错误的。函数arm_rfft_fast_f32的输出结果仍然是实部,虚部,实部,虚部.. 依次排列下去。
    函数arm_rfft_fast_f32在计算直流分量(也就是频率为0的值)的虚部上是有错误的。关于这点大家可以将实际的实部和虚部输出结果打印出来做对比,但差别很小,基本可以忽略。

32.2.2 arm_rfft_q15

函数定义如下:
    void arm_rfft_q15(
          const arm_rfft_instance_q15 * S,
          q15_t * pSrc,
          q15_t * pDst)
参数定义:
    [in]  *S    points to an instance of the Q15 RFFT/RIFFT structure.   
    [in]  *pSrc points to the input buffer.   
    [out] *pDst points to the output buffer.   
    return       none.  
注意事项:
结构arm_rfft_instance_q15的定义如下(在文件arm_math.h文件):
      typedef struct
       {
      uint32_t fftLenReal;   
      uint8_t ifftFlagR;   
           uint8_t bitReverseFlagR;                 
           uint32_t twidCoefRModifier;   
           q15_t *pTwiddleAReal;                    
          q15_t *pTwiddleBReal;                     
          const arm_cfft_instance_q15 *pCfft;      
      } arm_rfft_instance_q15;

    下面通过在开发板上运行函数arm_rfft_q15和arm_cfft_f32计算幅频响应,然后将相应的频率响应结果在Matlab上面绘制出来。
  1. /*
  2. *********************************************************************************************************
  3. *    函 数 名: arm_rfft_q15_app
  4. *    功能说明: 调用函数arm_rfft_q15计算1024点实数序列的幅频响应并跟使用函数arm_cfft_f32计算的结果做对比。
  5. *    形    参:无
  6. *    返 回 值: 无
  7. *********************************************************************************************************
  8. */
  9. static void arm_rfft_q15_app(void)
  10. {
  11. uint16_t i,j;
  12. arm_rfft_instance_q15 S;
  13. /* 实数序列FFT长度 */
  14. fftSize = 1024;
  15. /* 正变换 */
  16.     ifftFlag = 0;
  17. /* 码位倒序 */
  18.     doBitReverse = 1;
  19. /* 初始化结构体S */
  20.      arm_rfft_init_q15(&S, fftSize, ifftFlag, doBitReverse);
  21. /* 按照实部,虚部,实部,虚部..... 的顺序存储数据 */
  22. for(i=0; i<1024; i++)
  23. {
  24. /* 51.2Hz正弦波,采样率1024Hz。
  25.    arm_sin_q15输入参数的范围[0, 32768), 这里每20次为一个完整的正弦波,
  26.    32768 / 20 = 1638.4
  27. */
  28. j = i % 20;
  29.       testInput_q15_50hz[i] = arm_sin_q15(1638*j);
  30. }
  31. /* 1024点实序列快速FFT */
  32. arm_rfft_q15(&S, testInput_q15_50hz, testOutput_q15_50hz);
  33. /* 由于输出结果的格式是Q5,所以这里将定点数转换为浮点数 */
  34. for(i = 0; i < fftSize; i++)
  35. {
  36. testOutput_f32_10khz[i] = (float32_t)testOutput_q15_50hz[i]/32;
  37. }
  38. /* 为了方便对比,这里求解了1024组复数,实际上面的变化只有512组  
  39.    实际函数arm_rfft_q15只求解出了512组  */
  40.      arm_cmplx_mag_f32(testOutput_f32_10khz, testOutput, fftSize);
  41. /* 串口打印求解的模值 */
  42. for(i=0; i<fftSize; i++)
  43. {
  44. printf("%f\r\n", testOutput[i]);
  45. }
  46. printf("****************************分割线***************************************\r\n");
  47. for(i=0; i<1024; i++)
  48. {
  49. /* 51.2Hz正弦波,采样率1024Hz。
  50.    arm_sin_q15输入参数的范围[0, 32768), 这里每20次为一个完整的正弦波,
  51.    32768 / 20 = 1638.4
  52. */
  53. j = i % 20;
  54. testInput_f32_10khz[i*2] = (float32_t) arm_sin_q15(1638*j)/32768;
  55. /* 虚部全部置零 */
  56. testInput_f32_10khz[i*2+1] = 0;
  57. }
  58. arm_cfft_f32(&arm_cfft_sR_f32_len1024, testInput_f32_10khz, ifftFlag, doBitReverse);
  59. /* 求解模值  */
  60.      arm_cmplx_mag_f32(testInput_f32_10khz, testOutput, fftSize);
  61. /* 串口打印求解的模值 */
  62. for(i=0; i<fftSize; i++)
  63. {
  64. printf("%f\r\n", testOutput[i]);
  65. }
  66. }
复制代码
运行如上函数可以通过串口打印出函数arm_rfft_q15和arm_cfft_f32计算的幅频模值,下面通过Matlab绘制波形来对比这两种模值。
    对比前需要先将串口打印出的两组数据加载到Matlab中,arm_rfft_q15的计算结果起名signal,arm_cfft_f32的计算结果起名sampledata,
加载方法在前面的教程中已经讲解过,这里不做赘述了。Matlab中运行的代码如下:
Fs = 1000;             % 采样率
N  = 1024;            % 采样点数

n  = 0:N-1;            % 采样序列
f = n * Fs / N;          %真实的频率

subplot(3,1,1);
plot(f,  signal);        %绘制RFFT结果
title('实数FFT');
xlabel('时间');
ylabel('幅值');

subplot(3,1,2);
plot(f,  sampledata);   %CFFT结果
title('复数FFT');
xlabel('时间');
ylabel('幅值');

Matlab运行结果如下:

    从上面的前512点对比中,我们可以看出两者的计算结果是相符的。这里有一点要特别注意,官方文档中对于函数arm_rfft_q31输出结果的实部,虚部排列顺序说明是错误的。函数arm_rfft_q31的输出结果仍然是实部,虚部,实部,虚部.. 依次排列下去。

32.2.3 arm_rfft_q31

函数定义如下:
    void arm_rfft_q31(
          const arm_rfft_instance_q31 * S,
          q31_t * pSrc,
          q31_t * pDst)
参数定义:
    [in]  *S    points to an instance of the Q31 RFFT/RIFFT structure.   
    [in]  *pSrc points to the input buffer.   
    [out] *pDst points to the output buffer.   
    return       none.  
注意事项:
结构arm_rfft_instance_q31的定义如下(在文件arm_math.h文件):
      typedef struct
      {
          uint32_t fftLenReal;                     
          uint8_t ifftFlagR;                          
          uint8_t bitReverseFlagR;                  
          uint32_t twidCoefRModifier;                 
          q31_t *pTwiddleAReal;                       
          q31_t *pTwiddleBReal;                       
          const arm_cfft_instance_q31 *pCfft;      
      } arm_rfft_instance_q31;

    下面通过在开发板上运行函数arm_rfft_q31和arm_cfft_f32计算幅频响应,然后将相应的频率响应结果在Matlab上面绘制出来。
  1. /*
  2. *********************************************************************************************************
  3. *    函 数 名: arm_rfft_q31_app
  4. *    功能说明: 调用函数arm_rfft_q31计算1024点实数序列的幅频响应并跟使用函数arm_cfft_f32计算的结果做对比。
  5. *    形    参:无
  6. *    返 回 值: 无
  7. *********************************************************************************************************
  8. */
  9. static void arm_rfft_q31_app(void)
  10. {
  11. uint16_t i,j;
  12. arm_rfft_instance_q31 S;
  13. /* 实数序列FFT长度 */
  14. fftSize = 1024;
  15. /* 正变换 */
  16.     ifftFlag = 0;
  17. /* 码位倒序 */
  18.     doBitReverse = 1;
  19. /* 初始化结构体S */
  20.      arm_rfft_init_q31(&S, fftSize, ifftFlag, doBitReverse);
  21. /* 按照实部,虚部,实部,虚部..... 的顺序存储数据 */
  22. for(i=0; i<1024; i++)
  23. {
  24. /* 51.2Hz正弦波,采样率1024Hz。
  25.    arm_sin_q31输入参数的范围0-2^31, 这里每20次为一个完整的正弦波,
  26.    2^31 / 20 = 107374182.4
  27. */
  28. j = i % 20;
  29.       testInput_q31_50hz[i] = arm_sin_q31(107374182*j);
  30. }
  31. /* 1024点实序列快速FFT */
  32. arm_rfft_q31(&S, testInput_q31_50hz, testOutput_q31_50hz);
  33. /* 由于输出结果的格式是Q21,所以这里将定点数转换为浮点数 */
  34. for(i = 0; i < fftSize; i++)
  35. {
  36. /* 输出的数据是11.21格式,2^21 = 4194304*/
  37. testOutput_f32_10khz[i] = (float32_t)testOutput_q31_50hz[i]/2097152;
  38. }
  39. /* 为了方便对比,这里求解了1024组复数,实际上面的变化只有512组  
  40.    实际函数arm_rfft_q31只求解出了512组  */
  41.      arm_cmplx_mag_f32(testOutput_f32_10khz, testOutput, fftSize);
  42. /* 串口打印求解的模值 */
  43. for(i=0; i<fftSize; i++)
  44. {
  45. printf("%f\r\n", testOutput[i]);
  46. }
  47. printf("****************************分割线***************************************\r\n");
  48. for(i=0; i<1024; i++)
  49. {
  50. /* 51.2Hz正弦波,采样率1024Hz。
  51.    arm_sin_q31输入参数的范围0-2^31, 这里每20次为一个完整的正弦波,
  52.    2^31 / 20 = 107374182.4
  53. */
  54. j = i % 20;
  55.       testInput_f32_10khz[i*2] = (float32_t)arm_sin_q31(107374182*j)/2147483648;
  56. /* 虚部全部置零 */
  57. testInput_f32_10khz[i*2+1] = 0;
  58. }
  59. arm_cfft_f32(&arm_cfft_sR_f32_len1024, testInput_f32_10khz, ifftFlag, doBitReverse);
  60. /* 求解模值  */
  61.      arm_cmplx_mag_f32(testInput_f32_10khz, testOutput, fftSize);
  62. /* 串口打印求解的模值 */
  63. for(i=0; i<fftSize; i++)
  64. {
  65. printf("%f\r\n", testOutput[i]);
  66. }
  67. }
复制代码
运行如上函数可以通过串口打印出函数arm_rfft_q15和arm_cfft_f32计算的幅频模值,下面通过Matlab绘制波形来对比这两种模值。
    对比前需要先将串口打印出的两组数据加载到Matlab中,arm_rfft_q15的计算结果起名signal,arm_cfft_f32的计算结果起名sampledata,
加载方法在前面的教程中已经讲解过,这里不做赘述了。Matlab中运行的代码如下:
Fs = 1000;             % 采样率
N  = 1024;            % 采样点数

n  = 0:N-1;            % 采样序列
f = n * Fs / N;          %真实的频率

subplot(3,1,1);
plot(f,  signal);        %绘制RFFT结果
title('实数FFT');
xlabel('时间');
ylabel('幅值');

subplot(3,1,2);
plot(f,  sampledata);   %CFFT结果
title('复数FFT');
xlabel('时间');
ylabel('幅值');

Matlab运行结果如下:

    从上面的前512点对比中,我们可以看出两者的计算结果是相符的。这里有一点要特别注意,官方文档中对于函数arm_rfft_q31输出结果的实部,虚部排列顺序说明是错误的。函数arm_rfft_q31的输出结果仍然是实部,虚部,实部,虚部.. 依次排列下去。

32.3 总结
    使用实数FFT计算的时候要特别的注意本章节提到的几个错误点。
32.1.png
32.2.png
32.3.png
32.4.png
32.5.png
努力打造安富莱高质量微信公众号:点击扫描图片关注
回复

使用道具 举报

3

主题

24

回帖

33

积分

新手上路

积分
33
发表于 2018-7-16 16:15:28 | 显示全部楼层
在实数fft例程中,输入序列按照例程输入来看不是按照实部,虚部,实部,虚部..... 存储的吧?
/* 按照实部,虚部,实部,虚部..... 的顺序存储数据 */
for(i=0; i<1024; i++)
{
/* 50Hz正弦波,采样率1KHz */
testInput_f32_10khz[i] = 1.2f*arm_sin_f32(2*3.1415926f*50*i/1000)+1;
}
/
回复

使用道具 举报

1万

主题

6万

回帖

10万

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
105917
QQ
发表于 2018-7-16 16:17:46 | 显示全部楼层
confidentlin 发表于 2018-7-16 16:15
在实数fft例程中,输入序列按照例程输入来看不是按照实部,虚部,实部,虚部..... 存储的吧?
/* 按照实部 ...

是的,注释有误,忘了修改了。
回复

使用道具 举报

3

主题

24

回帖

33

积分

新手上路

积分
33
发表于 2018-7-16 16:20:43 | 显示全部楼层
那输入是按实部,实部,实部,实部.....存储的,输出是实部,虚部,实部,虚部.....存储的是吧?
回复

使用道具 举报

1万

主题

6万

回帖

10万

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
105917
QQ
发表于 2018-7-16 16:22:44 | 显示全部楼层
confidentlin 发表于 2018-7-16 16:20
那输入是按实部,实部,实部,实部.....存储的,输出是实部,虚部,实部,虚部.....存储的是吧?

对。
回复

使用道具 举报

4

主题

11

回帖

23

积分

新手上路

积分
23
发表于 2021-1-12 09:18:08 | 显示全部楼层

请问Q15的实数IFFT算法在使用的过程中,在进行完FFT运算后,直接将FFT的输出结果 testOutput_q15_50hz作为IFFT输入,可以吗?
static void arm_rifft_q15_app(void)
{
                fftSize = 1024;
                /* ??? */
                ifftFlag = 1;
                /* ???? */
                doBitReverse = 1;
       
          uint8_t i;
       
                arm_rfft_instance_q15 S;
       
                arm_rfft_init_q15(&S, fftSize, ifftFlag, doBitReverse);
               
                arm_rfft_q15(&S, testOutput_q15_50hz, testOutput1_q15_50hz);
       
          printf("printf the ifft data.");
                for(i=0; i<fftSize; i++)
        {
                printf("%d\r\n", testOutput1_q15_50hz);
        }

}
回复

使用道具 举报

1万

主题

6万

回帖

10万

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
105917
QQ
发表于 2021-1-12 10:27:31 | 显示全部楼层
liwei0225 发表于 2021-1-12 09:18
请问Q15的实数IFFT算法在使用的过程中,在进行完FFT运算后,直接将FFT的输出结果 testOutput_q15_50hz作 ...

你这个是用arm_rfft_q15做的FFT,然后将输出数据再用arm_rfft_q15做FFT逆变换?

早期我在示波器上玩过,相位上有些错位,正弦波应该好点,你用我的方法玩下
http://www.armbbs.cn/forum.php?m ... 2408&extra=page%3D1

上面是正变化,下面是逆变换

回复

使用道具 举报

0

主题

3

回帖

3

积分

新手上路

积分
3
发表于 2021-3-11 11:10:38 | 显示全部楼层
eric2013 发表于 2021-1-12 10:27
你这个是用arm_rfft_q15做的FFT,然后将输出数据再用arm_rfft_q15做FFT逆变换?

早期我在示波器上玩过 ...

关于 arm_rfft_q15 这个api,感觉很玄乎,自己测试发现有问题,根据官网的图片如下:
arm_rfft_q15正变换:
RFFTQ15.gif
arm_rfft_q15 反变换:
RIFFTQ15.gif

另外参考网上唯一的一个arm_rfft_q15 正变换又反变换的例子如下:
https://stackoverflow.com/questions/46147919/cmsis-real-fft-on-8192-samples-in-q15


我把主要源码拷贝如下:

  1. if(arm_rfft_init_q15(&S, 8192, 0, 1) != ARM_MATH_SUCCESS)
  2.     return state_error;

  3.   // perform FFT
  4.   arm_rfft_q15(&S, (q15_t*)0xC0000000, (q15_t*)0xC0400000);

  5.   // Post-shift by 12, in place (see doc)
  6.   arm_shift_q15((q15_t*)0xC0400000, 12, (q15_t*)0xC0400000, 16384);

  7.   // Init inverse FFT
  8.   if(arm_rfft_init_q15(&S, 8192, 1, 1) != ARM_MATH_SUCCESS)
  9.     return state_error;

  10.   // Perform iFFT
  11.   arm_rfft_q15(&S, (q15_t*)0xC0400000, (q15_t*)0xC0800000);

  12.   // Post shift
  13.   arm_shift_q15((q15_t*)0xC0800000, 12, (q15_t*)0xC0800000, 8192);
复制代码



自己尝试发现基本上经过正变换又反变换,没法复原原始信号,查看正变换之后经过12移位16bits的数据基本很多都32767了,都饱和了。

不知道FFT 和IFFT是不是上面代码那样使用,如果是为什么信号不能复原,另外IFFT之后是不是还需要shift的操作,能否指教下,谢谢!
回复

使用道具 举报

1万

主题

6万

回帖

10万

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
105917
QQ
发表于 2021-3-11 17:39:51 | 显示全部楼层
joinee0208 发表于 2021-3-11 11:10
关于 arm_rfft_q15 这个api,感觉很玄乎,自己测试发现有问题,根据官网的图片如下:
arm_rfft_q15正变 ...

你这个帖子的楼上回复,我用的是ST的那个汇编FFT做的逆变换,效果还正常。

ARM的这个定点FFT,我还没有测试过,后面我测试下,我们交流下。你可以先用浮点FFT测试下,那个没有定标问题,测试还方便些。
回复

使用道具 举报

0

主题

3

回帖

3

积分

新手上路

积分
3
发表于 2021-3-12 09:52:18 | 显示全部楼层
eric2013 发表于 2021-3-11 17:39
你这个帖子的楼上回复,我用的是ST的那个汇编FFT做的逆变换,效果还正常。

ARM的这个定点FFT,我还没 ...

恩,之所以用q15,是因为数据都是16bits的音频,用q15刚好不用转化,如果用float的就需要每次fft之前从float转化成q15,内存扩大一倍,还消耗计算资源,索性直接用q15,我查看了内网和外网的所有例子(除了上面我提到的stackoverflow的demo),几乎没有q15正变换再反变换回来的成功的例子(基本上正变换之后就不管了),所以觉得很奇怪,也有可能是我理解不对,后续你有测试进展,还望指教下,谢谢了!
回复

使用道具 举报

0

主题

3

回帖

3

积分

新手上路

积分
3
发表于 2021-5-12 08:54:32 | 显示全部楼层
eric2013 发表于 2021-3-11 17:39
你这个帖子的楼上回复,我用的是ST的那个汇编FFT做的逆变换,效果还正常。

ARM的这个定点FFT,我还没 ...

斑竹有没有空啊,能否帮忙测试下啊
回复

使用道具 举报

1万

主题

6万

回帖

10万

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
105917
QQ
发表于 2021-5-12 09:48:46 | 显示全部楼层
joinee0208 发表于 2021-5-12 08:54
斑竹有没有空啊,能否帮忙测试下啊

有时间了试试。
回复

使用道具 举报

1

主题

34

回帖

37

积分

新手上路

积分
37
发表于 2021-8-16 11:18:24 | 显示全部楼层
请问Eric,按您的例程,做1024点FFT,采样率100HZ, 调用函数 arm_rfft_fast_init_f32(&S_FFT, 1024);
arm_rfft_fast_f32(&S_FFT, m_stFFTBuffer.FFTBufIn, m_stFFTBuffer.FFTBufOut, 0);
m_stFFTBuffer.FFTBufIn 的数组长度是1024,实际每次采样256个点的数据存入,后面的全都是补零,每采样完256个数据后执行1024点的FFT, 得到的结果 m_stFFTBuffer.FFTBufOut 使用KEIL watch窗口观察 逐渐变大,25秒钟以后就溢出了,请问 是 arm_rfft_fast_f32 必须要采样满1024个点才能处理吗 ?同样的方式使用cr4_fft_1024_stm32 函数处理就不会有问题
回复

使用道具 举报

1万

主题

6万

回帖

10万

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
105917
QQ
发表于 2021-8-16 14:09:05 | 显示全部楼层
xml2028 发表于 2021-8-16 11:18
请问Eric,按您的例程,做1024点FFT,采样率100HZ, 调用函数 arm_rfft_fast_init_f32(&S_FFT, 1024);
arm_ ...

采集256个点后面点补0没问题。

输出数据逐渐变大直到溢出,不正常。
回复

使用道具 举报

1

主题

34

回帖

37

积分

新手上路

积分
37
发表于 2021-8-16 15:40:13 | 显示全部楼层
eric2013 发表于 2021-8-16 14:09
采集256个点后面点补0没问题。

输出数据逐渐变大直到溢出,不正常。

我用的雅特力的M4内核的AT32F403A单片机,
硬汉哥的DSP教程说:相应文件添加后还有最重要一条,要把stm32_dsp.h文件中的STM32F4的头文件, 我将stm32f4xx.h"替换成了 #include "at32f4xx.h" 这个会不会影响函数arm_rfft_fast_f32的调用 ,我编译也没有报错
回复

使用道具 举报

1

主题

34

回帖

37

积分

新手上路

积分
37
发表于 2021-8-16 15:57:35 | 显示全部楼层
现在发现调用了函数 arm_rfft_fast_f32(&S_FFT, m_stFFTBuffer.FFTBufIn, m_stFFTBuffer.FFTBufOut, 0); 后,      
数组m_stFFTBuffer.FFTBufIn 内的数据发生了很大变化
回复

使用道具 举报

1万

主题

6万

回帖

10万

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
105917
QQ
发表于 2021-8-16 16:13:05 | 显示全部楼层
xml2028 发表于 2021-8-16 15:57
现在发现调用了函数 arm_rfft_fast_f32(&S_FFT, m_stFFTBuffer.FFTBufIn, m_stFFTBuffer.FFTBufOut, 0); 后 ...

这是正常的,此变量空间也要参与运算。

使用DSP库的实数FFT函数arm_rfft_fast_f32,注意输入变量也会参与运算,存入临时数值
http://www.armbbs.cn/forum.php?m ... 7963&fromuid=58
(出处: 硬汉嵌入式论坛)
回复

使用道具 举报

5

主题

31

回帖

46

积分

新手上路

积分
46
发表于 2022-2-24 21:04:57 | 显示全部楼层
请问对于不定长的采样样本FFT的实现,有哪些参考案例或者库吗?
非2的N次方长度的。
感谢。
回复

使用道具 举报

1万

主题

6万

回帖

10万

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
105917
QQ
发表于 2022-2-25 00:24:13 | 显示全部楼层
darkness 发表于 2022-2-24 21:04
请问对于不定长的采样样本FFT的实现,有哪些参考案例或者库吗?
非2的N次方长度的。
感谢。

补零就行。
回复

使用道具 举报

1

主题

8

回帖

11

积分

新手上路

积分
11
发表于 2022-5-12 17:17:28 | 显示全部楼层
joinee0208 发表于 2021-3-11 11:10
关于 arm_rfft_q15 这个api,感觉很玄乎,自己测试发现有问题,根据官网的图片如下:
arm_rfft_q15正变 ...

哥们 你搞定了么  我现在也遇到逆变换的问题了  
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|小黑屋|Archiver|手机版|硬汉嵌入式论坛

GMT+8, 2024-3-29 05:08 , Processed in 0.322962 second(s), 28 queries .

Powered by Discuz! X3.4 Licensed

Copyright © 2001-2023, Tencent Cloud.

快速回复 返回顶部 返回列表