STM32F407实现傅里叶变换的三种方法【附源码】

2024-05-02 17:20
文章标签 实现 源码 方法 三种 变换 傅里叶 stm32f407

本文主要是介绍STM32F407实现傅里叶变换的三种方法【附源码】，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

一、浅谈傅里叶变换（Fourier Transformation，FT）

1、傅里叶级数

想要了解傅里叶变换，就要先了解一下什么是傅里叶级数。

 

 如图所示，通过不断合成不同频率的正弦波（频率分量），合成后的波形会逐渐形成类似方波的图形。若叠加的正弦波足够多，那么可以认为最终得到的波形就是方波。

换句话说，方波是由许许多多不同频率的正弦波组成的。

不仅如此，我们能见到的，能听到的，能想到的所有波，不管高矮胖瘦，也不论是否奇形怪状，都可以分解为不同频率的正弦波。

可以从三个维度去观察“波”的特性，他们分别是幅度、频域和时域。

在频域中，频率是从0Hz开始的。0Hz有什么意义呢？

在信号处理中，0Hz的正弦波分量实际上对应了信号的直流（DC）分量。正弦波在0Hz时并不是一个振荡的信号，而是一个恒定的值，即信号的均值或直流偏置。

对于任何实数信号，在频率0处（或称为DC分量）通常是非零的，除非信号的平均值为0。这个DC分量表示了信号中恒定的部分，即不随时间变化的部分。

2、傅里叶变换

不难看出，傅里叶级数处理的是周期的，连续的信号，频域上表现为离散的非周期函数。

而傅里叶变换处理的是非周期的，连续的信号，频域上表现为连续的非周期函数。

可以把傅里叶变换看做是处理周期无限大信号的一种方法。

3、离散傅里叶变换（Discrete Fourier Transform，DFT）

我们都知道，计算机采集的信号都是离散的，不可能是连续的。

这时候我们就需要用到离散傅里叶变换——DFT。

4、快速傅里叶变换（Fast Fourier Transform，FFT）

通俗点来讲，快速傅里叶变换就是高效版的离散傅里叶变换。

二、软件实现

软件有三种方法实现傅里叶变换。

在此之前，要先确定采样点数N和采样频率Fs。

详情请看

【电赛2020E题】从硬件到软件icon-default.png?t=N7T8https://mp.csdn.net/mp_blog/creation/editor/137835864里面有介绍FS和N。

方法一：STM32cubeMX 调用 DSP库

STM32CubeMX关于添加DSP库的使用icon-default.png?t=N7T8https://blog.csdn.net/WandZ123/article/details/125593908

方法二：调用FFT.c和FFT.h文件

方法三：巧用DFT

以电赛2020E题为例，题中要求对1KHz的原信号取五次谐波，那么我们只需关心1KHz到5KHz的的傅里叶变换。即仅需计算1KHz、2KHz...5KHz的点即可（代码中1KHz对应100个点）。

三、原码

这里只展示main.c，需要完整文件请看文章最后。

/* USER CODE BEGIN Header */
/********************************************************************************* @file           : main.c* @brief          : Main program body******************************************************************************* @attention** Copyright (c) 2024 STMicroelectronics.* All rights reserved.** This software is licensed under terms that can be found in the LICENSE file* in the root directory of this software component.* If no LICENSE file comes with this software, it is provided AS-IS.********************************************************************************/
/* USER CODE END Header */
/* Includes ------------------------------------------------------------------*/
#include "main.h"
#include "adc.h"
#include "tim.h"
#include "usart.h"
#include "gpio.h"/* Private includes ----------------------------------------------------------*/
/* USER CODE BEGIN Includes */
#include "stdio.h"
#include "arm_math.h"
#include "fft.h"
/* USER CODE END Includes *//* Private typedef -----------------------------------------------------------*/
/* USER CODE BEGIN PTD */
#define N 1024
/* USER CODE END PTD *//* Private define ------------------------------------------------------------*/
/* USER CODE BEGIN PD */float FFT_INPUT[N];
float FFT_OUTPUT[N];
float FFT_1024POINT[2*N];
complex FFT_POINT[N];
_Bool Flag=0;
arm_cfft_radix4_instance_f32 scfft;
/* USER CODE END PD *//* Private macro -------------------------------------------------------------*/
/* USER CODE BEGIN PM *//* USER CODE END PM *//* Private variables ---------------------------------------------------------*//* USER CODE BEGIN PV *//* USER CODE END PV *//* Private function prototypes -----------------------------------------------*/
void SystemClock_Config(void);
/* USER CODE BEGIN PFP *//* USER CODE END PFP *//* Private user code ---------------------------------------------------------*/
/* USER CODE BEGIN 0 */
float abscomplex(complex *in)
{float out;out=in->real*in->real+in->imag*in->imag;return sqrt(out);
}
void DSP_FFTTC(float* IN)//IN是指主函数的FFT_INPUT[i]
{u32 i;for(i=0;i<N;i++){FFT_1024POINT[2*i]=*(IN+i);//存入实部，存入该数组偶数位置FFT_1024POINT[2*i+1]=0;//虚部置零，存入该数组奇数位置}arm_cfft_radix4_init_f32(&scfft,N,0,1);//初始化FFT结构体scfftarm_cfft_radix4_f32(&scfft,FFT_1024POINT);//进行FFT运算arm_cmplx_mag_f32(FFT_1024POINT,FFT_OUTPUT,N);//计算FFT的幅度printf("用DSP库的FFT\n");for(i=0;i<N/2;i++){if(i==0)//在FFT的结果中，索引为0的元素对应的是直流分量，即频率为0的分量{FFT_OUTPUT[i]=FFT_OUTPUT[i]/N;//对于索引为0的元素，对应直流分量，需要除以N进行修正}else{FFT_OUTPUT[i]=FFT_OUTPUT[i]*2/N;//对于其他元素，乘以2并除以N进行修正}printf("第%d个FFT:%.2f\n",i,FFT_OUTPUT[i]);//测试}
}
void L_FFT(float* IN)
{u32 i;for(i=0;i<N;i++){FFT_POINT[i].real=*(IN+i);FFT_POINT[i].imag=*(IN+i);/*fft.h中定义了结构体
typedef struct complex //复数类型
{float real;		//实部float imag;		//虚部
}complex;*/}fft(N,FFT_POINT);//傅里叶变换函数c_abs(FFT_POINT,FFT_OUTPUT,N);//求所有复数的模printf("普通的FFT\n");for(i=0;i<N/2;i++){/*详情注释见DSP库FFT*/if(i==0){FFT_OUTPUT[i]=FFT_OUTPUT[i]/N;}else{FFT_OUTPUT[i]=FFT_OUTPUT[i]*2/N;}printf("第%d个FFT:%.2f\n",i,FFT_OUTPUT[i]);//测试}
}void DFT(float* IN)
{u32 i,k;complex y={.real=0,.imag=0};printf("DFT\n");for(k=0;k<6;k++){for(i=0;i<N;i++){y.real+=*(IN+i)*cos(-2*100*PI*i*k/N);y.imag+=*(IN+i)*sin(-2*100*PI*i*k/N);}y.real=2*y.real/N;y.imag=2*y.imag/N;printf("第%d个点的值为%.2f\n",k,abscomplex(&y));}
}
/* USER CODE END 0 *//*** @brief  The application entry point.* @retval int*/
int main(void)
{/* USER CODE BEGIN 1 */u16 i;/* USER CODE END 1 *//* MCU Configuration--------------------------------------------------------*//* Reset of all peripherals, Initializes the Flash interface and the Systick. */HAL_Init();/* USER CODE BEGIN Init *//* USER CODE END Init *//* Configure the system clock */SystemClock_Config();/* USER CODE BEGIN SysInit *//* USER CODE END SysInit *//* Initialize all configured peripherals */MX_GPIO_Init();MX_ADC1_Init();MX_USART1_UART_Init();MX_TIM2_Init();MX_TIM3_Init();/* USER CODE BEGIN 2 */printf("开始\n");HAL_TIM_Base_Start_IT(&htim2);HAL_TIM_Base_Start(&htim3);HAL_TIM_PWM_Start(&htim3,TIM_CHANNEL_1);//PA6输出PWM/* USER CODE END 2 *//* Infinite loop *//* USER CODE BEGIN WHILE */while (1){if(Flag){HAL_TIM_Base_Stop(&htim2);for(i=0;i<1024;i++){FFT_INPUT[i]=FFT_INPUT[i]*3.3f/4096;//除以4096是因为单片机采到的值经过转换后数据寄存器最高位为2的12次方，正好是4096，乘3.3v是为了转化为电压printf("ADC_VALUE:%.2f\n",FFT_INPUT[i]);}DSP_FFTTC(FFT_INPUT);L_FFT(FFT_INPUT);DFT(FFT_INPUT);HAL_Delay(1000);Flag=0;HAL_TIM_Base_Start(&htim2);}/* USER CODE END WHILE *//* USER CODE BEGIN 3 */}/* USER CODE END 3 */
}/*** @brief System Clock Configuration* @retval None*/
void SystemClock_Config(void)
{RCC_OscInitTypeDef RCC_OscInitStruct = {0};RCC_ClkInitTypeDef RCC_ClkInitStruct = {0};/** Configure the main internal regulator output voltage*/__HAL_RCC_PWR_CLK_ENABLE();__HAL_PWR_VOLTAGESCALING_CONFIG(PWR_REGULATOR_VOLTAGE_SCALE1);/** Initializes the RCC Oscillators according to the specified parameters* in the RCC_OscInitTypeDef structure.*/RCC_OscInitStruct.OscillatorType = RCC_OSCILLATORTYPE_HSE;RCC_OscInitStruct.HSEState = RCC_HSE_ON;RCC_OscInitStruct.PLL.PLLState = RCC_PLL_ON;RCC_OscInitStruct.PLL.PLLSource = RCC_PLLSOURCE_HSE;RCC_OscInitStruct.PLL.PLLM = 4;RCC_OscInitStruct.PLL.PLLN = 168;RCC_OscInitStruct.PLL.PLLP = RCC_PLLP_DIV2;RCC_OscInitStruct.PLL.PLLQ = 4;if (HAL_RCC_OscConfig(&RCC_OscInitStruct) != HAL_OK){Error_Handler();}/** Initializes the CPU, AHB and APB buses clocks*/RCC_ClkInitStruct.ClockType = RCC_CLOCKTYPE_HCLK|RCC_CLOCKTYPE_SYSCLK|RCC_CLOCKTYPE_PCLK1|RCC_CLOCKTYPE_PCLK2;RCC_ClkInitStruct.SYSCLKSource = RCC_SYSCLKSOURCE_PLLCLK;RCC_ClkInitStruct.AHBCLKDivider = RCC_SYSCLK_DIV1;RCC_ClkInitStruct.APB1CLKDivider = RCC_HCLK_DIV4;RCC_ClkInitStruct.APB2CLKDivider = RCC_HCLK_DIV2;if (HAL_RCC_ClockConfig(&RCC_ClkInitStruct, FLASH_LATENCY_5) != HAL_OK){Error_Handler();}
}/* USER CODE BEGIN 4 */void HAL_TIM_PeriodElapsedCallback(TIM_HandleTypeDef *htim)
{static u16 counter=0;if(htim==&htim2){if(!Flag){HAL_GPIO_TogglePin(GPIOB,GPIO_PIN_15);//翻转PB15 用示波器测量以判断定时器频率FFT_INPUT[counter++]=get_adcvalue();if(counter==1024){counter=0;Flag=1;}}}
}/* USER CODE END 4 *//*** @brief  This function is executed in case of error occurrence.* @retval None*/
void Error_Handler(void)
{/* USER CODE BEGIN Error_Handler_Debug *//* User can add his own implementation to report the HAL error return state */__disable_irq();while (1){}/* USER CODE END Error_Handler_Debug */
}#ifdef  USE_FULL_ASSERT
/*** @brief  Reports the name of the source file and the source line number*         where the assert_param error has occurred.* @param  file: pointer to the source file name* @param  line: assert_param error line source number* @retval None*/
void assert_failed(uint8_t *file, uint32_t line)
{/* USER CODE BEGIN 6 *//* User can add his own implementation to report the file name and line number,ex: printf("Wrong parameters value: file %s on line %d\r\n", file, line) *//* USER CODE END 6 */
}
#endif /* USE_FULL_ASSERT */

链接：https://pan.baidu.com/s/1DyjH0DaRx3qYEF7I_QdbpA?pwd=hk0x 
提取码：hk0x

参考链接

通俗易懂的理解傅里叶变换（一）[收藏]

这篇关于STM32F407实现傅里叶变换的三种方法【附源码】的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！


http://www.chinasem.cn/article/954776。

相关文章

hdu1043（八数码问题，广搜 + hash（实现状态压缩） ）

hdu1043（八数码问题，广搜 + hash（实现状态压缩） ）

利用康拓展开将一个排列映射成一个自然数，然后就变成了普通的广搜题。 #include<iostream>#include<algorithm>#include<string>#include<stack>#include<queue>#include<map>#include<stdio.h>#include<stdlib.h>#include<ctype.h>#inclu
阅读更多...
JAVA智听未来一站式有声阅读平台听书系统小程序源码

JAVA智听未来一站式有声阅读平台听书系统小程序源码

智听未来，一站式有声阅读平台听书系统 🌟&nbsp;开篇：遇见未来，从“智听”开始 在这个快节奏的时代，你是否渴望在忙碌的间隙，找到一片属于自己的宁静角落？是否梦想着能随时随地，沉浸在知识的海洋，或是故事的奇幻世界里？今天，就让我带你一起探索“智听未来”——这一站式有声阅读平台听书系统，它正悄悄改变着我们的阅读方式，让未来触手可及！ 📚&nbsp;第一站：海量资源，应有尽有 走进“智听
阅读更多...
【C++】_list常用方法解析及模拟实现

【C++】_list常用方法解析及模拟实现

相信自己的力量，只要对自己始终保持信心，尽自己最大努力去完成任何事，就算事情最终结果是失败了，努力了也不留遗憾。💓💓💓 目录   ✨说在前面 🍋知识点一：什么是list？ •🌰1.list的定义 •🌰2.list的基本特性 •🌰3.常用接口介绍 🍋知识点二：list常用接口 •🌰1.默认成员函数 🔥构造函数（⭐） 🔥析构函数 •🌰2.list对象
阅读更多...
【Prometheus】PromQL向量匹配实现不同标签的向量数据进行运算

【Prometheus】PromQL向量匹配实现不同标签的向量数据进行运算

✨✨ 欢迎大家来到景天科技苑✨✨ 🎈🎈 养成好习惯，先赞后看哦~🎈🎈 🏆 作者简介：景天科技苑 🏆《头衔》：大厂架构师，华为云开发者社区专家博主，阿里云开发者社区专家博主，CSDN全栈领域优质创作者，掘金优秀博主，51CTO博客专家等。 🏆《博客》：Python全栈，前后端开发，小程序开发，人工智能，js逆向，App逆向，网络系统安全，数据分析，Django，fastapi
阅读更多...
让树莓派智能语音助手实现定时提醒功能

让树莓派智能语音助手实现定时提醒功能

最初的时候是想直接在rasa 的chatbot上实现，因为rasa本身是带有remindschedule模块的。不过经过一番折腾后，忽然发现，chatbot上实现的定时，语音助手不一定会有响应。因为，我目前语音助手的代码设置了长时间无应答会结束对话，这样一来，chatbot定时提醒的触发就不会被语音助手获悉。那怎么让语音助手也具有定时提醒功能呢？ 我最后选择的方法是用threading.Time
阅读更多...
Android实现任意版本设置默认的锁屏壁纸和桌面壁纸（两张壁纸可不一致）

Android实现任意版本设置默认的锁屏壁纸和桌面壁纸（两张壁纸可不一致）

客户有些需求需要设置默认壁纸和锁屏壁纸  在默认情况下 这两个壁纸是相同的  如果需要默认的锁屏壁纸和桌面壁纸不一样 需要额外修改 Android13实现 替换默认桌面壁纸： 将图片文件替换frameworks/base/core/res/res/drawable-nodpi/default_wallpaper.*  （注意不能是bmp格式） 替换默认锁屏壁纸： 将图片资源放入vendo
阅读更多...
C#实战|大乐透选号器[6]：实现实时显示已选择的红蓝球数量

C#实战|大乐透选号器[6]：实现实时显示已选择的红蓝球数量

哈喽，你好啊，我是雷工。 关于大乐透选号器在前面已经记录了5篇笔记，这是第6篇； 接下来实现实时显示当前选中红球数量，蓝球数量； 以下为练习笔记。 01 效果演示 当选择和取消选择红球或蓝球时，在对应的位置显示实时已选择的红球、蓝球的数量； 02 标签名称 分别设置Label标签名称为：lblRedCount、lblBlueCount
阅读更多...
浅谈主机加固，六种有效的主机加固方法

浅谈主机加固，六种有效的主机加固方法

在数字化时代，数据的价值不言而喻，但随之而来的安全威胁也日益严峻。从勒索病毒到内部泄露，企业的数据安全面临着前所未有的挑战。为了应对这些挑战，一种全新的主机加固解决方案应运而生。 MCK主机加固解决方案，采用先进的安全容器中间件技术，构建起一套内核级的纵深立体防护体系。这一体系突破了传统安全防护的局限，即使在管理员权限被恶意利用的情况下，也能确保服务器的安全稳定运行。 普适主机加固措施：
阅读更多...
webm怎么转换成mp4？这几种方法超多人在用！

webm怎么转换成mp4？这几种方法超多人在用！

webm怎么转换成mp4？WebM作为一种新兴的视频编码格式，近年来逐渐进入大众视野，其背后承载着诸多优势，但同时也伴随着不容忽视的局限性，首要挑战在于其兼容性边界，尽管WebM已广泛适应于众多网站与软件平台，但在特定应用环境或老旧设备上，其兼容难题依旧凸显，为用户体验带来不便，再者，WebM格式的非普适性也体现在编辑流程上，由于它并非行业内的通用标准，编辑过程中可能会遭遇格式不兼容的障碍，导致操
阅读更多...
透彻！驯服大型语言模型（LLMs）的五种方法，及具体方法选择思路

透彻！驯服大型语言模型（LLMs）的五种方法，及具体方法选择思路

引言 随着时间的发展，大型语言模型不再停留在演示阶段而是逐步面向生产系统的应用，随着人们期望的不断增加，目标也发生了巨大的变化。在短短的几个月的时间里，人们对大模型的认识已经从对其zero-shot能力感到惊讶，转变为考虑改进模型质量、提高模型可用性。 「大语言模型(LLMs)其实就是利用高容量的模型架构（例如Transformer）对海量的、多种多样的数据分布进行建模得到，它包含了大量的先验
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2