STM32CubeIDE STM32F03C6T6开发板 模拟数据生成 FFT傅里叶运算

已于 2023-09-13 13:09:56 修改
阅读量924 收藏 6
点赞数 8
文章标签： stm32 嵌入式硬件 单片机
于 2023-09-13 11:07:09 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_43476450/article/details/132828887
版权

开发环境：WIN10

开发软件：STM32CubeIDE

使用硬件：STM32F03C6T6开发板

准备TTF软件：是 GCC版本的。

1.新建项目：TTF_Simulate_data_f103c6

       打开STM32CubeIDE-->文件-->新建-->STM32 project

选择器件： 搜索103c6，选择STM32F03A，-->下一步

--- 输入项目名称，点击“完成”，自动进入STM32CubeMX界面

 二，coudMX设置：

 1.项目使用单片机资源：

（1）外部时钟72MH 、

（2）TIM3 

（3）USART1

（4）PC13

---  外部时钟72MH

RCC设置：

 ---  HCLK(MHz) 位置输入 72  ，回车，点击OK， 自动完成设置。

---  TIM3 设置

---  USART1设置：

--- 中断设置：打开TIM3，USART中断

 ---PC13设置：

--- 保存，确认

至此，项目初始化完成。

三，导入FFT软件，选择文件系统导入，浏览到FFT软件所在文件夹

 --- stm32_dsp.h，table_fft.h 导入 Inc 文件夹

 --- cr4_fft_256_stm32.s，cr4_fft_1024_stm32.s 导入Src 文件夹

  

 也可以分别直接复制，粘贴到 Inc 和 Src 这两个文件夹内。

四、添加自己的程序到main.c  找到绿色代码对应位置，复制范围内的代码到自己的main.c。

添加头文件声明

/* USER CODE END Header */

/* Includes ------------------------------------------------------------------*/

#include "main.h"

/* Private includes ----------------------------------------------------------*/

/* USER CODE BEGIN Includes */

#include "stdio.h"

#include "math.h"

#include "stm32_dsp.h"

#include "table_fft.h"

/* USER CODE END Includes */

添加变量声明：

/* Private variables ---------------------------------------------------------*/

TIM_HandleTypeDef htim3;

UART_HandleTypeDef huart1;

/* USER CODE BEGIN PV */

//串口printf重定向

#ifdef __GNUC__

/* With GCC/RAISONANCE, small printf (option LD Linker->Libraries->Small printf

set to 'Yes') calls __io_putchar() */

#define PUTCHAR_PROTOTYPE int __io_putchar(int ch)

#else

#define PUTCHAR_PROTOTYPE int fputc(int ch, FILE *f)

#endif /* __GNUC__ */

/**

* @brief Retargets the C library printf function to the USART.

* @param None

* @retval None

*/

PUTCHAR_PROTOTYPE {

/* Place your implementation of fputc here */

/* e.g. write a character to the EVAL_COM1 and Loop until the end of transmission */

HAL_UART_Transmit(&huart1, (uint8_t*) &ch, 1, 0xFFFF);

return ch;

}

int display_on = 1;

char frequency_domain_value = 0;

char frequency_domain = 1;

char simulate_data = 2;

//串口输出使用

#define FQV frequency_domain_value //频谱

//串口绘图使用

#define FQD frequency_domain //频谱

#define SMD simulate_data //仿真数据

#define SPT 256 //Sample points 采样点

long simulate_data_array[SPT]; //模拟数据数组simulate_data_array

long FFT_out_array[SPT / 2];     //FFT运算结果

long Frequency_amplitude_array[SPT / 2];  //频谱 FFT运算结果转换

float PI2 = 6.28318530717959;

uint32_t SPF = 10240; //采样率

/* USER CODE END PV */

添加函数声明：

/* Private function prototypes -----------------------------------------------*/

void SystemClock_Config(void);

static void MX_GPIO_Init(void);

static void MX_TIM3_Init(void);

static void MX_USART1_UART_Init(void);

/* USER CODE BEGIN PFP */

void simulate_data_generation(void); //Generate analog sampling data 生成模拟采样数据

void FFT_256(void);

void get_Frequency_amplitude(void); Frequency amplitude 频率幅度

//打印绘图

char print_drawing(char str); //频谱值

/* USER CODE END PFP */

 主程序初始化 ：复位后只执行一次

/* USER CODE BEGIN 2 */

HAL_TIM_Base_Start_IT(&htim3);      //启动定时器3

simulate_data_generation(); //仿真数据生成 //Generate analog sampling data 生成模拟采样数据

FFT_256(); //

get_Frequency_amplitude(); Frequency amplitude 频率幅度

//打印绘图

print_drawing(FQV); //频谱值

/* USER CODE END 2 */

主程序循环

/* USER CODE END WHILE */

/* USER CODE BEGIN 3 */

if (display_on == 1) {

HAL_GPIO_TogglePin(GPIOC, GPIO_PIN_13);   //LED灯状态反转

//打印绘图

print_drawing(FQD); //频谱 

//打印绘图

//需要看频谱，将注释掉下一行代码

// print_drawing(SMD); //仿真数据

display_on = 0;

}

//----------------

}

/* USER CODE END 3 */

添加自己的函数：

//仿真数据生成 SPT是采样点:256,SPF是采样频率:10.240KHz

/* USER CODE BEGIN 4 */

//生成FFT需要的仿真数据

void simulate_data_generation(void) {

unsigned short i;

float Signal;

for (i = 0; i < SPT; i++) {

Signal = 100 * sin(PI2 * i * 60.0 / SPF)    //频率60

+ 100 * sin(PI2 * i * 1280.0 / SPF)      //频率1280

+ 100 * sin(PI2 * i * 5000.0 / SPF);     //频率5000

simulate_data_array[i] = ((signed short) Signal) << 16;

}

}

//FFT

void FFT_256() {

cr4_fft_256_stm32(FFT_out_array, simulate_data_array, SPT);

}

//频谱

void get_Frequency_amplitude() {

signed short lX, lY;

float X, Y, Mag;

for (uint16_t i = 0; i < SPT / 2; i++) {

lX = (FFT_out_array[i] << 16) >> 16;

lY = (FFT_out_array[i] >> 16);

X = SPT * ((float) lX) / 32768;

Y = SPT * ((float) lY) / 32768;

Mag = sqrt(X * X + Y * Y) / SPT;

if (i == 0)

Frequency_amplitude_array[i] = (unsigned long) (Mag * 32768);

else

Frequency_amplitude_array[i] = (unsigned long) (Mag * 65536);

}

}

//打印绘图

char print_drawing(char str) {

for (int i = 0; i < SPT / 2; i++) {

switch (str) {

case 0: //FQV串口输出使用 频谱

printf("i = %d", i); //\r\n

printf(" F = %ld", ((SPF / 256) * i)); //175\r\n

printf(" V = %ld\r\n", Frequency_amplitude_array[i]); //\r\n

break;

//串口绘图使用

case 1: //FQD频谱

printf("%ld\r\n", Frequency_amplitude_array[i]); //

break;

case 2: //SMD仿真数据

printf("%ld\r\n", simulate_data_array[i]); //

break;

}

printf("\r\n");

}

return 0;

}

/* USER CODE END 4 */

中断设置：打开 stm32f1xx_it.c 找到对应位置

添加变量：

/* Private variables ---------------------------------------------------------*/

/* USER CODE BEGIN PV */

extern int display_on;

/* USER CODE END PV */

定时器3

void TIM3_IRQHandler(void) {

/* USER CODE BEGIN TIM3_IRQn 0 */

/* USER CODE END TIM3_IRQn 0 */

HAL_TIM_IRQHandler(&htim3);

/* USER CODE BEGIN TIM3_IRQn 1 */

display_on = 1;

/* USER CODE END TIM3_IRQn 1 */

}

 程序部分结束。

构建项目成功：

下载程序：

STM32F03C6T6开发板 运行：

代码：print_drawing(FQV); //频谱值

//FQV串口输出  打印结果：i = 系列 与FFT采样点有关，F = 频率，V = 频率对应的值。

FFT采样点，采样率确定后，频谱的频率就固定下来了，去看SPF = 10240; //采样率，函数 simulate_data_generation(void) 中的生成频率，60Hz，1280Hz，5000Hz。

10240(采样率)  / 256(采样点) = 40  ，很明显生成频率60Hz不在40的整数倍上，看不到60Hz的峰值点，但是邻居频率点幅值是递减的；1280Hz，5000Hz能看到峰值。

i = 0   F = 0   V = 17

i = 1   F = 40   V = 70

i = 2   F = 80   V = 60

i = 3   F = 120   V = 24

i = 4   F = 160   V = 13

i = 5   F = 200   V = 11

i = 6   F = 240   V = 10

i = 7   F = 280   V = 8

i = 8   F = 320   V = 8

i = 9   F = 360   V = 8

i = 10   F = 400   V = 10

i = 11   F = 440   V = 10

i = 12   F = 480   V = 8

i = 13   F = 520   V = 7

i = 14   F = 560   V = 10

i = 15   F = 600   V = 10

i = 16   F = 640   V = 7

i = 17   F = 680   V = 4

i = 18   F = 720   V = 2

i = 19   F = 760   V = 6

i = 20   F = 800   V = 2

i = 21   F = 840   V = 6

i = 22   F = 880   V = 6

i = 23   F = 920   V = 6

i = 24   F = 960   V = 6

i = 25   F = 1000   V = 6

i = 26   F = 1040   V = 2

i = 27   F = 1080   V = 6

i = 28   F = 1120   V = 4

i = 29   F = 1160   V = 4

i = 30   F = 1200   V = 6

i = 31   F = 1240   V = 2

i = 32   F = 1280   V = 96

i = 33   F = 1320   V = 4

i = 34   F = 1360   V = 4

i = 35   F = 1400   V = 6

i = 36   F = 1440   V = 4

i = 37   F = 1480   V = 6

i = 38   F = 1520   V = 4

i = 39   F = 1560   V = 4

i = 40   F = 1600   V = 4

i = 41   F = 1640   V = 4

i = 42   F = 1680   V = 6

i = 43   F = 1720   V = 6

i = 44   F = 1760   V = 4

i = 45   F = 1800   V = 6

i = 46   F = 1840   V = 6

i = 47   F = 1880   V = 6

i = 48   F = 1920   V = 2

i = 49   F = 1960   V = 2

i = 50   F = 2000   V = 2

i = 51   F = 2040   V = 2

i = 52   F = 2080   V = 2

i = 53   F = 2120   V = 2

i = 54   F = 2160   V = 2

i = 55   F = 2200   V = 2

i = 56   F = 2240   V = 4

i = 57   F = 2280   V = 4

i = 58   F = 2320   V = 5

i = 59   F = 2360   V = 4

i = 60   F = 2400   V = 2

i = 61   F = 2440   V = 4

i = 62   F = 2480   V = 2

i = 63   F = 2520   V = 2

i = 64   F = 2560   V = 0

i = 65   F = 2600   V = 0

i = 66   F = 2640   V = 2

i = 67   F = 2680   V = 2

i = 68   F = 2720   V = 2

i = 69   F = 2760   V = 2

i = 70   F = 2800   V = 2

i = 71   F = 2840   V = 0

i = 72   F = 2880   V = 2

i = 73   F = 2920   V = 2

i = 74   F = 2960   V = 2

i = 75   F = 3000   V = 4

i = 76   F = 3040   V = 0

i = 77   F = 3080   V = 2

i = 78   F = 3120   V = 4

i = 79   F = 3160   V = 2

i = 80   F = 3200   V = 2

i = 81   F = 3240   V = 4

i = 82   F = 3280   V = 2

i = 83   F = 3320   V = 4

i = 84   F = 3360   V = 0

i = 85   F = 3400   V = 2

i = 86   F = 3440   V = 4

i = 87   F = 3480   V = 2

i = 88   F = 3520   V = 2

i = 89   F = 3560   V = 4

i = 90   F = 3600   V = 4

i = 91   F = 3640   V = 2

i = 92   F = 3680   V = 2

i = 93   F = 3720   V = 4

i = 94   F = 3760   V = 4

i = 95   F = 3800   V = 2

i = 96   F = 3840   V = 2

i = 97   F = 3880   V = 4

i = 98   F = 3920   V = 4

i = 99   F = 3960   V = 4

i = 100   F = 4000   V = 0

i = 101   F = 4040   V = 4

i = 102   F = 4080   V = 2

i = 103   F = 4120   V = 0

i = 104   F = 4160   V = 2

i = 105   F = 4200   V = 4

i = 106   F = 4240   V = 4

i = 107   F = 4280   V = 2

i = 108   F = 4320   V = 0

i = 109   F = 4360   V = 4

i = 110   F = 4400   V = 4

i = 111   F = 4440   V = 4

i = 112   F = 4480   V = 2

i = 113   F = 4520   V = 2

i = 114   F = 4560   V = 2

i = 115   F = 4600   V = 2

i = 116   F = 4640   V = 2

i = 117   F = 4680   V = 2

i = 118   F = 4720   V = 0

i = 119   F = 4760   V = 2

i = 120   F = 4800   V = 2

i = 121   F = 4840   V = 2

i = 122   F = 4880   V = 4

i = 123   F = 4920   V = 2

i = 124   F = 4960   V = 0

i = 125   F = 5000   V = 96

i = 126   F = 5040   V = 2

i = 127   F = 5080   V = 0

函数 get_Frequency_amplitude() 转换的数据，打印绘图结果：

代码：print_drawing(FQD); //频谱

函数 simulate_data_generation(void) 生成的数据，打印绘图结果：

代码：print_drawing(SMD); //仿真数据

不同采样率，不同生成频率：

参考链接：

【玩转单片机系列002】 如何使用STM32提供的DSP库进行FFT - 知乎 (zhihu.com)

串口 printf 重定向参考：

STM32CubeIDE 利用自带HAL库 串口收发_stm32cube实现串口接收数据_zateper的博客-CSDN博客

qq_666076
关注
STM32的FFT程序
08-30
基于STM32的FFT算法程序，可以做FFT程序设计的参考资料
STM32F103 FFT算法的实现
12-15
1024点的FFT算法实现
STM32 CubeMX配置ADC+DMA进行FFT（1）
ANNOLRA的博客
07-29 9060
STM32 cubeMX配置ADC+DMA，计算采样率和ADC电压 最近，因为项目需要处理音频信号，对AD采集的音频信号进行FFT运算，记录一下学习过程。 主要内容： 1.ADC和DMA各配置有什么实际意义 2.配置好的ADC如何计算采样率和电压值 cubeMX配置ADC+DMA 这里我用的是STM32F407处理器，与其他的处理器思路大体一致 ADC配置 Mode是ADC工作模式 我们设置为independent mode（独立模式）因为只需要单独ADC使用，如果是多个ADC联合使用才需要更改。 Clo
用STM32实现FFT
qq_44164303的博客
11-11 3万+
前言 电子设计大赛已经过去很久了，一直想写一篇关于FFT的文章也没有来得及，现在写一下来记录分享一下。 本篇文章不讲复杂的FFT原理，只讲如何在stm32里面怎么实现FFT 一、FFT是什么，能干啥？ FFT(fast Fourier transform)快速傅里叶变换，把时域信号变换到频域，至于换到频域后能做啥，那就多了 二、先用matlab来认识一下FFT 1.上代码 代码如下： fs = 1024; %采样频率 NPT = 1024; %样点个数 t=0:1/f
【STM32】使用STM32提供的DSP库进行FFT（附详细代码）
最新发布
Symean的博客
09-03 433
最近，因为项目需要在STM32F103系列处理器上，对采集的音频信号进行FFT运算，然而STM32F103毕竟不是STM32F4系列的处理器，对于一般的FFT运算程序还是比较缓慢的。幸亏官方提供了针对FFT的官方库，但是去官网找了半天居然找不到那个库的下载，好像官方早就把那个库下架了，估计是为了给带DSP指令集的F4系列让路。然后就只好从别人的项目中把这个官方库给扒出来了……提取码: dn2d。
关于STM32CubeIDE 使用汇编FFT库cr4_fft_1024_stm32.s编译报错
cjx_tom的博客
11-24 2146
近期开始使用STM32CubeIDE，想要将keil上写的标准库程序移植到Cube里用HAL库实现。 最近移植DSP库FFT
STM32F030-Demo.rar_DEMO_STM32F030_STM32F030工程_STM32F03x_FWLib_st
07-15
stm32f030工程源码，原子风格代码
STM32F103C8T6资料
02-21
6. **开发板与硬件资源**：描述中的“STM32F103C8T6核心板资料”可能包括了开发板的原理图和PCB布局，这些资料对初学者和开发者来说非常有价值，他们可以了解如何连接外部元件，以及如何在实际电路中使用STM32F103C8...
STM32F03系列 STM32F030K6T6核心板资料
03-17
STM32F03系列 STM32F030K6T6核心板资料， 1、核心板硬件资料 2、核心板特性与使用说明 3、芯片资料 4、编程外设库 5、其他元器件资料 6、软件安装包 7、工程模板(库函数、寄存器、HAL库)
YD-STM32F030（051）C8T6最小系统资料.rar_DEMO_STM32F030_YD-LED实验_stm32F03
07-14
STM32F030是意法半导体（STMicroelectronics）推出的一款基于ARM Cortex-M0内核的微控制器，属于STM32系列的超低功耗产品。在本压缩包中，我们主要关注的是STM32F030C8T6型号，这是一款具有32KB闪存和2KB RAM的芯片...
FFT在STM32处理器上的实现完整代码
10-28
FFT在STM32处理器上的实现完整代码,可以直接运行在STM32系列芯片上
STM32的傅里叶算法
09-03
用STM32系列芯片做傅里叶变换的算法和程序，希望对大家有帮助
STM32进行FFT傅里叶变换 CUBEMX
02-28
STM32进行FFT傅里叶变换 CUBEMX
stm32f103上实现的一个fft算法
04-17
以定时器配合DMA ADC模式，形成20KHz的adc采样率，再进行1024个点的fft产生频谱数组，usart打印出来，可供学习或项目中参考。
STM32 +差分ad纯C实现2048点FFT，采集信号频谱图
01-04
自己写，用于学习交流，采样用ad在时域以一定频率等间隔采样，完成一组数据采集存储后，MCU进行FFT输出频谱数组
STM32F030 - BH1750_stm32f030iic_stm32f030模拟iic_传感器模拟IIC_stm32f03
09-11
在本文中，我们将深入探讨如何使用STM32F030单片机通过模拟IIC（Inter-Integrated Circuit）通信协议来驱动BH1750光照传感器。STM32F030是一款基于ARM Cortex-M0内核的微控制器，具有高效能、低功耗的特点，适用于...
STM32CubeMx移植DSP库 傅立叶变化(FFT)测试
岸上游泳
10-25 4426
基于STM32移植DSP库实现FFT变换示例及程序编写
stm32f1单片机上用FFT测量信号频率（高精度、过程详细）
热门推荐
weixin_43368814的博客
01-13 9万+
在STM32F1系列单片机上面实现FFT 最近需要做一台基于stm32的示波器，测量信号参数，用单片机上面一些传统的方法局限性还是比较大，故开始琢磨直接上fft。 ...
手把手STM32G431CubeMx+keil添加DSP库并实现FFT简单实验
m0_74923693的博客
03-19 3174
三天前的积分赛需要用到FFT来测量信号的频谱，于是乎在网上找找教程。发现之前已经有学长分享过经验，但是自己跟着添加的过程中还是遇到了一些keil里面的报错，在网上又找了找相关文章，最后终于是解决了问题，并且在校初选中运用FFT测量了信号的频谱，于是想写下这篇文章当作备忘。贴下学长的帖子的链接：https://blog.csdn.net/qq_34022877/article/details/123190943。
stm32cube fft
12-03
STM32Cube是STMicroelectronics公司推出的一套基于STM32微控制器的软件生态系统。FFT（Fast Fourier Transform，快速傅里叶变换）是一种用于信号处理和频谱分析的算法，能够将时域内的信号转换为频域内的频谱图。在STM32Cube软件生态系统中，STMicroelectronics提供了基于STM32微控制器的FFT库，这个库可以方便地在STM32微控制器上实现FFT算法，用于对信号进行频谱分析和处理。 利用STM32Cube FFT库，开发者可以通过简单的API调用实现对信号的FFT变换，并获取其频谱信息。这对于需要实时信号频谱分析的应用场景非常有用，比如音频处理、振动检测、频谱分析等领域。开发者可以通过STM32CubeMX工具配置对应的STM32微控制器，然后利用CubeIDE集成开发环境进行FFT库的集成和开发，从而实现对信号的频域分析。 通过STM32Cube FFT库，开发者可以充分利用STM32微控制器的性能优势，实现高效的信号处理和频谱分析，提供更加丰富的功能和应用体验。同时，STM32Cube还提供了丰富的示例代码和文档，对开发者进行指导和帮助，使得FFT算法的应用更加便捷和高效。 总之，STM32Cube FFT库为基于STM32微控制器的应用提供了方便、高效的信号处理和频谱分析解决方案，为开发者在嵌入式系统中实现信号处理提供了强大的支持和工具。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值