【STM32H7的DSP教程】第31章 STM32H7实数浮点FFT(支持单精度和双精度)

news/2024/10/30 11:31:34/

完整版教程下载地址:http://www.armbbs.cn/forum.php?mod=viewthread&tid=94547

第31章       STM32H7实数浮点FFT(支持单精度和双精度)

本章主要讲解实数浮点FTT,支持单精度和双精度。

目录

31.1 初学者重要提示

31.2 实数浮点FFT说明

31.3 单精度函数arm_rfft_fast_f32的使用(含幅频和相频)

31.3.1 函数说明

31.3.2 使用举例并和Matlab比较

31.4 双精度函数arm_rfft_fast_f64的使用(含幅频和相频)

31.4.1 函数说明

31.4.2 使用举例并和Matlab比较

31.5 实验例程说明(MDK)

31.6 实验例程说明(IAR)

31.7 总结


 

31.1 初学者重要提示

  1.   与上一章节的复数FFT相比,实数FFT仅需用户输入实部即可。输出结果根据FFT的对称性,也仅输出一半的频谱。

31.2 实数浮点FFT说明

CMSIS DSP库里面包含一个专门用于计算实数序列的FFT库,很多情况下,用户只需要计算实数序列即可。计算同样点数FFT的实数序列要比计算同样点数的虚数序列有速度上的优势。

快速的rfft算法是基于混合基cfft算法实现的。

一个N点的实数序列FFT正变换采用下面的步骤实现:

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

一个N点的实数序列FFT逆变换采用下面的步骤实现:

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

31.3 单精度函数arm_rfft_fast_f32的使用(含幅频和相频)

31.3.1 函数说明

函数原型:

void arm_rfft_fast_f32(const arm_rfft_fast_instance_f32 * S,float32_t * p,float32_t * pOut,uint8_t ifftFlag)

函数描述:

这个函数用于单精度浮点实数FFT。

函数参数:

  •   第1个参数是封装好的浮点FFT例化,需要用户先调用函数arm_rfft_fast_init_f32初始化,然后供此函数arm_rfft_fast_f32调用。支持32, 64, 128, 256, 512, 1024, 2048, 4096点FFT。

比如做1024点FFT,代码如下:

arm_rfft_fast_instance_f32 S;

arm_rfft_fast_init_f32(&S, 1024);

arm_rfft_fast_f32(&S, testInput_f32, testOutput_f32, ifftFlag);

  •   第2个参数是实数地址,比如我们要做1024点实数FFT,要保证有1024个缓冲。
  •   第3个参数是FFT转换结果,转换结果不是实数了,而是复数,按照实部,虚拟,实部,虚部,依次排列。比如做1024点FFT,这里的输出也会有1024个数据,即512个复位。
  •   第4个参数用于设置正变换和逆变换,ifftFlag=0表示正变换,ifftFlag=1表示逆变换。

31.3.2 使用举例并和Matlab比较

下面通过在开发板上运行这个函数并计算幅频相应,然后再与Matlab计算的结果做对比。

/*
*********************************************************************************************************
*    函 数 名: arm_rfft_f32_app
*    功能说明: 调用函数arm_rfft_fast_f32计算幅频和相频
*    形    参:无
*    返 回 值: 无
*********************************************************************************************************
*/
static void arm_rfft_f32_app(void)
{uint16_t i;arm_rfft_fast_instance_f32 S;/* 正变换 */ifftFlag = 0; /* 初始化结构体S中的参数 */arm_rfft_fast_init_f32(&S, TEST_LENGTH_SAMPLES);for(i=0; i<1024; i++){/* 波形是由直流分量,50Hz正弦波组成,波形采样率1024,初始相位60° */testInput_f32[i] = 1 + cos(2*3.1415926f*50*i/1024 + 3.1415926f/3);}/* 1024点实序列快速FFT */ arm_rfft_fast_f32(&S, testInput_f32, testOutput_f32, ifftFlag);/* 为了方便跟函数arm_cfft_f32计算的结果做对比,这里求解了1024组模值,实际函数arm_rfft_fast_f32只求解出了512组  */ arm_cmplx_mag_f32(testOutput_f32, testOutputMag_f32, TEST_LENGTH_SAMPLES);printf("=========================================\r\n");    /* 求相频 */PowerPhaseRadians_f32(testOutput_f32, Phase_f32, TEST_LENGTH_SAMPLES, 0.5f);/* 串口打印求解的幅频和相频 */for(i=0; i<TEST_LENGTH_SAMPLES; i++){printf("%f, %f\r\n", testOutputMag_f32[i], Phase_f32[i]);}
}

运行函数arm_rfft_f32_app可以通过串口打印出计算的模值和相角,下面我们就通过Matlab计算的模值和相角跟arm_rfft_fast_f32计算的做对比。

对比前需要先将串口打印出的数据加载到Matlab中,并给这个数组起名sampledata,加载方法在前面的教程的第13章13.6小结已经讲解,这里不做赘述了。Matlab中运行的代码如下:

Fs = 1024;               % 采样率
N  = 1024;               % 采样点数
n  = 0:N-1;              % 采样序列
t  = 0:1/Fs:1-1/Fs;      % 时间序列
f = n * Fs / N;          %真实的频率%波形是由直流分量,50Hz正弦波正弦波组成
x = 1 + cos(2*pi*50*t + pi/3)   ;  
y = fft(x, N);               %对原始信号做FFT变换
Mag = abs(y);subplot(2,2,1);
plot(f, Mag); 
title('Matlab计算幅频响应');
xlabel('频率');
ylabel('赋值');subplot(2,2,2);
realvalue = real(y);
imagvalue = imag(y);
plot(f, atan2(imagvalue, realvalue)*180/pi.*(Mag>=200)); 
title('Matlab计算相频响应');
xlabel('频率');
ylabel('相角');subplot(2,2,3);
plot(f, sampledata1);  %绘制STM32计算的幅频相应
title('STM32计算幅频响应');
xlabel('频率');
ylabel('赋值');subplot(2,2,4);
plot(f, sampledata2);   %绘制STM32计算的相频相应
title('STM32计算相频响应');
xlabel('频率');
ylabel('相角');

运行Matlab后的输出结果如下:

 

从上面的对比结果中可以看出,从上面的前512点对比中,我们可以看出两者的计算结果是相符的Matlab和函数arm_rfft_fast_f32计算的结果基本是一直的。幅频响应求出的幅值和相频响应中的求出的初始相角都是没问题的。

31.4 双精度函数arm_rfft_fast_f64的使用(含幅频和相频)

31.4.1 函数说明

函数原型:

void arm_rfft_fast_f64(arm_rfft_fast_instance_f64 * S,float64_t * p,float64_t * pOut,uint8_t ifftFlag)

函数描述:

这个函数用于双精度浮点实数FFT。

函数参数:

  •   第1个参数是封装好的浮点FFT例化,需要用户先调用函数arm_rfft_fast_init_f64初始化,然后供此函数arm_rfft_fast_f64调用。支持32, 64, 128, 256, 512, 1024, 2048, 4096点FFT。

比如做1024点FFT,代码如下:

arm_rfft_fast_instance_f64 S;

arm_rfft_fast_init_f64(&S, 1024);

arm_rfft_fast_f64(&S, testInput_f64, testOutput_f64, ifftFlag);

  •   第2个参数是实数地址,比如我们要做1024点实数FFT,要保证有1024个缓冲。
  •   第3个参数是FFT转换结果,转换结果不是实数了,而是复数,按照实部,虚拟,实部,虚部,依次排列。比如做1024点FFT,这里的输出也会有1024个数据,即512个复位。
  •   第4个参数用于设置正变换和逆变换,ifftFlag=0表示正变换,ifftFlag=1表示逆变换

31.4.2 使用举例并和Matlab比较

下面通过在开发板上运行这个函数并计算幅频相应,然后再与Matlab计算的结果做对比。

/*
*********************************************************************************************************
*    函 数 名: arm_rfft_f64_app
*    功能说明: 调用函数arm_rfft_fast_f64计算幅频和相频
*    形    参:无
*    返 回 值: 无
*********************************************************************************************************
*/
static void arm_rfft_f64_app(void)
{uint16_t i;float64_t lX,lY;arm_rfft_fast_instance_f64 S;/* 正变换 */ifftFlag = 0; /* 初始化结构体S中的参数 */arm_rfft_fast_init_f64(&S, TEST_LENGTH_SAMPLES);for(i=0; i<1024; i++){/* 波形是由直流分量,50Hz正弦波组成,波形采样率1024,初始相位60° */testInput_f64[i] = 1 + cos(2*3.1415926*50*i/1024 + 3.1415926/3);}/* 1024点实序列快速FFT */ arm_rfft_fast_f64(&S, testInput_f64, testOutput_f64, ifftFlag);/* 求解模值  */ for (i =0; i < TEST_LENGTH_SAMPLES; i++){lX = testOutput_f64[2*i];                    /* 实部*/lY = testOutput_f64[2*i+1];                   /* 虚部 */  testOutputMag_f64[i] = sqrt(lX*lX+ lY*lY);   /* 求模 */}printf("=========================================\r\n");    /* 求相频 */PowerPhaseRadians_f64(testOutput_f64, Phase_f64, TEST_LENGTH_SAMPLES, 0.5);/* 串口打印幅值和相频 */for(i=0; i<TEST_LENGTH_SAMPLES; i++){printf("%.11f, %.11f\r\n", testOutputMag_f64[i], Phase_f64[i]);}    }

运行函数arm_rfft_f64_app可以通过串口打印出计算的模值和相角,下面我们就通过Matlab计算的模值和相角跟arm_rfft_fast_f32计算的做对比。

对比前需要先将串口打印出的数据加载到Matlab中,并给这个数组起名sampledata,加载方法在前面的教程的第13章13.6小结已经讲解,这里不做赘述了。Matlab中运行的代码如下:

Fs = 1024;               % 采样率
N  = 1024;               % 采样点数
n  = 0:N-1;              % 采样序列
t  = 0:1/Fs:1-1/Fs;      % 时间序列
f = n * Fs / N;          %真实的频率%波形是由直流分量,50Hz正弦波正弦波组成
x = 1 + cos(2*pi*50*t + pi/3)   ;  
y = fft(x, N);               %对原始信号做FFT变换
Mag = abs(y);subplot(2,2,1);
plot(f, Mag); 
title('Matlab计算幅频响应');
xlabel('频率');
ylabel('赋值');subplot(2,2,2);
realvalue = real(y);
imagvalue = imag(y);
plot(f, atan2(imagvalue, realvalue)*180/pi.*(Mag>=200)); 
title('Matlab计算相频响应');
xlabel('频率');
ylabel('相角');subplot(2,2,3);
plot(f, sampledata1);  %绘制STM32计算的幅频相应
title('STM32计算幅频响应');
xlabel('频率');
ylabel('赋值');subplot(2,2,4);
plot(f, sampledata2);   %绘制STM32计算的相频相应
title('STM32计算相频响应');
xlabel('频率');
ylabel('相角');

运行Matlab后的输出结果如下:

从上面的对比结果中可以看出,从上面的前512点对比中,我们可以看出两者的计算结果是相符的Matlab和函数arm_rfft_fast_f64计算的结果基本是一直的。幅频响应求出的幅值和相频响应中的求出的初始相角都是没问题的。

31.5 实验例程说明(MDK)

配套例子:

V7-221_实数浮点FTT(支持单精度和双精度)

实验目的:

  1. 学习实数浮点FFT,支持单精度浮点和双精度浮点

实验内容:

  1. 启动一个自动重装软件定时器,每100ms翻转一次LED2。
  2. 按下按键K1,串口打印1024点实数单精度FFT的幅频响应和相频响应。
  3. 按下按键K2,串口打印1024点实数双精度FFT的幅频响应和相频响应。

使用AC6注意事项

特别注意附件章节C的问题

上电后串口打印的信息:

波特率 115200,数据位 8,奇偶校验位无,停止位 1。

RTT方式打印信息:

程序设计:

  系统栈大小分配:

  RAM空间用的DTCM:

  硬件外设初始化

硬件外设的初始化是在 bsp.c 文件实现:

/*
*********************************************************************************************************
*    函 数 名: bsp_Init
*    功能说明: 初始化所有的硬件设备。该函数配置CPU寄存器和外设的寄存器并初始化一些全局变量。只需要调用一次
*    形    参:无
*    返 回 值: 无
*********************************************************************************************************
*/
void bsp_Init(void)
{/* 配置MPU */MPU_Config();/* 使能L1 Cache */CPU_CACHE_Enable();/* STM32H7xx HAL 库初始化,此时系统用的还是H7自带的64MHz,HSI时钟:- 调用函数HAL_InitTick,初始化滴答时钟中断1ms。- 设置NVIC优先级分组为4。*/HAL_Init();/* 配置系统时钟到400MHz- 切换使用HSE。- 此函数会更新全局变量SystemCoreClock,并重新配置HAL_InitTick。*/SystemClock_Config();/* Event Recorder:- 可用于代码执行时间测量,MDK5.25及其以上版本才支持,IAR不支持。- 默认不开启,如果要使能此选项,务必看V7开发板用户手册第8章*/    
#if Enable_EventRecorder == 1  /* 初始化EventRecorder并开启 */EventRecorderInitialize(EventRecordAll, 1U);EventRecorderStart();
#endifbsp_InitKey();        /* 按键初始化,要放在滴答定时器之前,因为按钮检测是通过滴答定时器扫描 */bsp_InitTimer();      /* 初始化滴答定时器 */bsp_InitUart();    /* 初始化串口 */bsp_InitExtIO();    /* 初始化FMC总线74HC574扩展IO. 必须在 bsp_InitLed()前执行 */    bsp_InitLed();        /* 初始化LED */    
}

  MPU配置和Cache配置:

数据Cache和指令Cache都开启。配置了AXI SRAM区(本例子未用到AXI SRAM),FMC的扩展IO区。

/*
*********************************************************************************************************
*    函 数 名: MPU_Config
*    功能说明: 配置MPU
*    形    参: 无
*    返 回 值: 无
*********************************************************************************************************
*/
static void MPU_Config( void )
{MPU_Region_InitTypeDef MPU_InitStruct;/* 禁止 MPU */HAL_MPU_Disable();/* 配置AXI SRAM的MPU属性为Write back, Read allocate,Write allocate */MPU_InitStruct.Enable           = MPU_REGION_ENABLE;MPU_InitStruct.BaseAddress      = 0x24000000;MPU_InitStruct.Size             = MPU_REGION_SIZE_512KB;MPU_InitStruct.AccessPermission = MPU_REGION_FULL_ACCESS;MPU_InitStruct.IsBufferable     = MPU_ACCESS_BUFFERABLE;MPU_InitStruct.IsCacheable      = MPU_ACCESS_CACHEABLE;MPU_InitStruct.IsShareable      = MPU_ACCESS_NOT_SHAREABLE;MPU_InitStruct.Number           = MPU_REGION_NUMBER0;MPU_InitStruct.TypeExtField     = MPU_TEX_LEVEL1;MPU_InitStruct.SubRegionDisable = 0x00;MPU_InitStruct.DisableExec      = MPU_INSTRUCTION_ACCESS_ENABLE;HAL_MPU_ConfigRegion(&MPU_InitStruct);/* 配置FMC扩展IO的MPU属性为Device或者Strongly Ordered */MPU_InitStruct.Enable           = MPU_REGION_ENABLE;MPU_InitStruct.BaseAddress      = 0x60000000;MPU_InitStruct.Size             = ARM_MPU_REGION_SIZE_64KB;    MPU_InitStruct.AccessPermission = MPU_REGION_FULL_ACCESS;MPU_InitStruct.IsBufferable     = MPU_ACCESS_BUFFERABLE;MPU_InitStruct.IsCacheable      = MPU_ACCESS_NOT_CACHEABLE;    MPU_InitStruct.IsShareable      = MPU_ACCESS_NOT_SHAREABLE;MPU_InitStruct.Number           = MPU_REGION_NUMBER1;MPU_InitStruct.TypeExtField     = MPU_TEX_LEVEL0;MPU_InitStruct.SubRegionDisable = 0x00;MPU_InitStruct.DisableExec      = MPU_INSTRUCTION_ACCESS_ENABLE;HAL_MPU_ConfigRegion(&MPU_InitStruct);/*使能 MPU */HAL_MPU_Enable(MPU_PRIVILEGED_DEFAULT);
}/*
*********************************************************************************************************
*    函 数 名: CPU_CACHE_Enable
*    功能说明: 使能L1 Cache
*    形    参: 无
*    返 回 值: 无
*********************************************************************************************************
*/
static void CPU_CACHE_Enable(void)
{/* 使能 I-Cache */SCB_EnableICache();/* 使能 D-Cache */SCB_EnableDCache();
}

  主功能:

主程序实现如下操作:

  •   启动一个自动重装软件定时器,每100ms翻转一次LED2。
  •   按下按键K1,串口打印1024点实数单精度FFT的幅频响应和相频响应。
  •   按下按键K2,串口打印1024点实数双精度FFT的幅频响应和相频响应。
/*
*********************************************************************************************************
*    函 数 名: main
*    功能说明: c程序入口
*    形    参: 无
*    返 回 值: 错误代码(无需处理)
*********************************************************************************************************
*/
int main(void)
{uint8_t ucKeyCode;        /* 按键代码 */bsp_Init();        /* 硬件初始化 */PrintfLogo();    /* 打印例程信息到串口1 */PrintfHelp();    /* 打印操作提示信息 */bsp_StartAutoTimer(0, 100);    /* 启动1个100ms的自动重装的定时器 *//* 进入主程序循环体 */while (1){bsp_Idle();        /* 这个函数在bsp.c文件。用户可以修改这个函数实现CPU休眠和喂狗 */if (bsp_CheckTimer(0))    /* 判断定时器超时时间 */{/* 每隔100ms 进来一次 */bsp_LedToggle(4);    /* 翻转LED2的状态 */   }ucKeyCode = bsp_GetKey();    /* 读取键值, 无键按下时返回 KEY_NONE = 0 */if (ucKeyCode != KEY_NONE){switch (ucKeyCode){case KEY_DOWN_K1:            /* K1键按下 */arm_rfft_f32_app();break;case KEY_DOWN_K2:            /* K2键按下 */arm_rfft_f64_app();break;default:/* 其它的键值不处理 */break;}}}
}

31.6 实验例程说明(IAR)

配套例子:

V7-221_实数浮点FTT(支持单精度和双精度)

实验目的:

  1. 学习实数浮点FFT,支持单精度浮点和双精度浮点

实验内容:

  1. 启动一个自动重装软件定时器,每100ms翻转一次LED2。
  2. 按下按键K1,串口打印1024点实数单精度FFT的幅频响应和相频响应。
  3. 按下按键K2,串口打印1024点实数双精度FFT的幅频响应和相频响应。

上电后串口打印的信息:

波特率 115200,数据位 8,奇偶校验位无,停止位 1。

RTT方式打印信息:

程序设计:

  系统栈大小分配:

  RAM空间用的DTCM:

  硬件外设初始化

硬件外设的初始化是在 bsp.c 文件实现:

/*
*********************************************************************************************************
*    函 数 名: bsp_Init
*    功能说明: 初始化所有的硬件设备。该函数配置CPU寄存器和外设的寄存器并初始化一些全局变量。只需要调用一次
*    形    参:无
*    返 回 值: 无
*********************************************************************************************************
*/
void bsp_Init(void)
{/* 配置MPU */MPU_Config();/* 使能L1 Cache */CPU_CACHE_Enable();/* STM32H7xx HAL 库初始化,此时系统用的还是H7自带的64MHz,HSI时钟:- 调用函数HAL_InitTick,初始化滴答时钟中断1ms。- 设置NVIC优先级分组为4。*/HAL_Init();/* 配置系统时钟到400MHz- 切换使用HSE。- 此函数会更新全局变量SystemCoreClock,并重新配置HAL_InitTick。*/SystemClock_Config();/* Event Recorder:- 可用于代码执行时间测量,MDK5.25及其以上版本才支持,IAR不支持。- 默认不开启,如果要使能此选项,务必看V7开发板用户手册第8章*/    
#if Enable_EventRecorder == 1  /* 初始化EventRecorder并开启 */EventRecorderInitialize(EventRecordAll, 1U);EventRecorderStart();
#endifbsp_InitKey();        /* 按键初始化,要放在滴答定时器之前,因为按钮检测是通过滴答定时器扫描 */bsp_InitTimer();      /* 初始化滴答定时器 */bsp_InitUart();    /* 初始化串口 */bsp_InitExtIO();    /* 初始化FMC总线74HC574扩展IO. 必须在 bsp_InitLed()前执行 */    bsp_InitLed();        /* 初始化LED */    
}

 MPU配置和Cache配置:

数据Cache和指令Cache都开启。配置了AXI SRAM区(本例子未用到AXI SRAM),FMC的扩展IO区。

/*
*********************************************************************************************************
*    函 数 名: MPU_Config
*    功能说明: 配置MPU
*    形    参: 无
*    返 回 值: 无
*********************************************************************************************************
*/
static void MPU_Config( void )
{MPU_Region_InitTypeDef MPU_InitStruct;/* 禁止 MPU */HAL_MPU_Disable();/* 配置AXI SRAM的MPU属性为Write back, Read allocate,Write allocate */MPU_InitStruct.Enable           = MPU_REGION_ENABLE;MPU_InitStruct.BaseAddress      = 0x24000000;MPU_InitStruct.Size             = MPU_REGION_SIZE_512KB;MPU_InitStruct.AccessPermission = MPU_REGION_FULL_ACCESS;MPU_InitStruct.IsBufferable     = MPU_ACCESS_BUFFERABLE;MPU_InitStruct.IsCacheable      = MPU_ACCESS_CACHEABLE;MPU_InitStruct.IsShareable      = MPU_ACCESS_NOT_SHAREABLE;MPU_InitStruct.Number           = MPU_REGION_NUMBER0;MPU_InitStruct.TypeExtField     = MPU_TEX_LEVEL1;MPU_InitStruct.SubRegionDisable = 0x00;MPU_InitStruct.DisableExec      = MPU_INSTRUCTION_ACCESS_ENABLE;HAL_MPU_ConfigRegion(&MPU_InitStruct);/* 配置FMC扩展IO的MPU属性为Device或者Strongly Ordered */MPU_InitStruct.Enable           = MPU_REGION_ENABLE;MPU_InitStruct.BaseAddress      = 0x60000000;MPU_InitStruct.Size             = ARM_MPU_REGION_SIZE_64KB;    MPU_InitStruct.AccessPermission = MPU_REGION_FULL_ACCESS;MPU_InitStruct.IsBufferable     = MPU_ACCESS_BUFFERABLE;MPU_InitStruct.IsCacheable      = MPU_ACCESS_NOT_CACHEABLE;    MPU_InitStruct.IsShareable      = MPU_ACCESS_NOT_SHAREABLE;MPU_InitStruct.Number           = MPU_REGION_NUMBER1;MPU_InitStruct.TypeExtField     = MPU_TEX_LEVEL0;MPU_InitStruct.SubRegionDisable = 0x00;MPU_InitStruct.DisableExec      = MPU_INSTRUCTION_ACCESS_ENABLE;HAL_MPU_ConfigRegion(&MPU_InitStruct);/*使能 MPU */HAL_MPU_Enable(MPU_PRIVILEGED_DEFAULT);
}/*
*********************************************************************************************************
*    函 数 名: CPU_CACHE_Enable
*    功能说明: 使能L1 Cache
*    形    参: 无
*    返 回 值: 无
*********************************************************************************************************
*/
static void CPU_CACHE_Enable(void)
{/* 使能 I-Cache */SCB_EnableICache();/* 使能 D-Cache */SCB_EnableDCache();
}

  主功能:

主程序实现如下操作:

  •   启动一个自动重装软件定时器,每100ms翻转一次LED2。
  •   按下按键K1,串口打印1024点实数单精度FFT的幅频响应和相频响应。
  •   按下按键K2,串口打印1024点实数双精度FFT的幅频响应和相频响应。
/*
*********************************************************************************************************
*    函 数 名: main
*    功能说明: c程序入口
*    形    参: 无
*    返 回 值: 错误代码(无需处理)
*********************************************************************************************************
*/
int main(void)
{uint8_t ucKeyCode;        /* 按键代码 */bsp_Init();        /* 硬件初始化 */PrintfLogo();    /* 打印例程信息到串口1 */PrintfHelp();    /* 打印操作提示信息 */bsp_StartAutoTimer(0, 100);    /* 启动1个100ms的自动重装的定时器 *//* 进入主程序循环体 */while (1){bsp_Idle();        /* 这个函数在bsp.c文件。用户可以修改这个函数实现CPU休眠和喂狗 */if (bsp_CheckTimer(0))    /* 判断定时器超时时间 */{/* 每隔100ms 进来一次 */bsp_LedToggle(4);    /* 翻转LED2的状态 */   }ucKeyCode = bsp_GetKey();    /* 读取键值, 无键按下时返回 KEY_NONE = 0 */if (ucKeyCode != KEY_NONE){switch (ucKeyCode){case KEY_DOWN_K1:            /* K1键按下 */arm_rfft_f32_app();break;case KEY_DOWN_K2:            /* K2键按下 */arm_rfft_f64_app();break;default:/* 其它的键值不处理 */break;}}}
}

31.7 总结

本章节设计到实数FFT实现,有兴趣的可以深入了解源码的实现。

 


http://www.ppmy.cn/news/520365.html

相关文章

STM32 实数FFT 极速配置

想要进行FFT操作就得配置DSP环境&#xff0c;操作如下 STM32 keil极速配置DSP环境 关于FFT&#xff0c;先上操作再讲含义 操作 添加定义 #define fftsize 256 //对256个数据进行FFT #define ifftFlag 0float32_t testInput_f32[256]{0}; float32_t testOutput_f32[256…

STM32 FFT算法实现

DSP 库运行环境搭建 在 MDK 里面搭建 STM32F4 的 DSP 运行环境(使用.lib 方式)是很简单的&#xff0c;分为 3 个步骤&#xff1a; 1&#xff0c; 添加文件。 首先&#xff0c;我们在例程工程目录下新建&#xff1a;DSP_LIB 文件夹&#xff0c;存放我们将要添加的文件&#xff…

STM32+CubeMX 通过RMS和FFT进行波形识别

波形识别 本文所展示的程序可以用于分辨正弦波、三角波、方波三种波形。 文章目录 波形识别思路可以判断波形的两个特点时域方面频域方面 外设配置 & DSP库配置代码部分串口重定向时域部分变量定义ADC采集求最大值&#xff0c;最小值&#xff0c;幅值取出波形的一个周期求…

STM32 FIR滤波器使用要点,STM32 DSP篇

1.综述 STM32搭配的CMSIS的DSP库已经写好了FIR&#xff0c;我们做的就是用好用稳定&#xff01; 2.使用流程 2.1使用MATLAB的Filter Designer生成滤波器的参数&#xff08;直接在MATLAB命令中输入filterDesigner即可调出此工具&#xff0c;配置好参数后&#xff0c;选择Targ…

【STM32F407的DSP教程】第30章 STM32F407复数浮点FFT(支持单精度和双精度)

完整版教程下载地址&#xff1a;http://www.armbbs.cn/forum.php?modviewthread&tid94547 第30章 STM32F407复数浮点FFT&#xff08;支持单精度和双精度&#xff09; 本章主要讲解复数浮点FTT&#xff0c;支持单精度和双精度。 目录 30.1 初学者重要提示 30.2 …

C嵌入汇编之vld1.f32和vst1.f32指令理解

想完成类似与memcpy&#xff0c;使用arm的neon指令完成&#xff0c;第一次代码 关键点 第一点&#xff1a; add %0,%0,#32&#xff0c;表示每次移动32/4&#xff1d;&#xff18;个float 第二点&#xff1a; subs r0,r0,#1才能改变状态标志位&#xff0c;sub不行 第三点&a…

%.2f

% .2f 是一个格式化字符串&#xff0c;用于在 Python 程序中格式化浮点数。 其中&#xff0c;% 表示这是一个格式化字符串&#xff0c;.2 表示保留两位小数&#xff0c;f 表示这是一个浮点数。 例如&#xff0c;在 Python 中可以使用如下代码将一个浮点数格式化为字符串&#x…

NEON 常用函数讲解

目录 一、基本的加载存储操作 1. vld1q_f32 2.vst1q_f32 3.vld2q_f32 4.vst2q_f32 5.vld3q_f32 6.vst3q_f32 7.vld4q_f32, vst4q_f32 二、特殊操作 1.vdupq_n_f32 2.vzipq_f32 3.vuzpq_f32 4.vcombine_f32 5. vget_low_f32 6. vget_high_f32 7. vtrnq_f32 8. v…