嵌入式开发:傅里叶变换(5):STM32和Matlab联调验证FFT

news/2025/2/26 8:30:09

目录

1. MATLAB获取 STM32 的原始数据

2. 将数据上传到电脑

3. MATLAB 接收数据并验证

STM32进行傅里叶代码

结果分析


STM32 和 MATLAB 联调是嵌入式开发中常见的工作流程,通常目的是将 STM32 采集的数据或控制信号传输到 MATLAB 中进行实时处理、分析和可视化,同时也可以将在开发板上运行算法,然后再与Matlab计算的结果做对比,对我们的算法进行验证

以在项目中对磁通门传感器进行数据处理为例,进行在 STM32 和 MATLAB 联调,以下是三个关键步骤:

1. MATLAB获取 STM32 的原始数据

在 STM32 上,首先需要从传感器或其他外设获取原始数据。这个过程通常涉及:

  • 配置 ADC、传感器接口(如 I2C、SPI 等)来读取传感器数据。
  • 将采集到的数据存储到数组或缓冲区中,准备传输。

示例代码:

// 这是采集到的原始数据:Flugate_frame.dataBuf的缓存数组,包含磁通门传感器X、Y、Z三个轴的数据
// 数据填充过程,读取磁通门传感器X轴数据到数组中,先对X轴的数据做FFT变换
								                                                                                                                         
int groupCount = Flugate_frame.dataLen / 12;  /* 每组12字节(x,y,z) */								                                                                                                                     
for (int i = 0; i < groupCount; i++)								                                                                                                                     
{ 										                                                                                                    
    if(Flugate_Date.Data_idx < 512)										                                                                                                                     
    {	
        //从原始数据缓存数组中提取X轴的数据到数组Flugate_Date.axis_X_Data_FFT_Input中,作为FFT的输入											                        
        memcpy(&Flugate_Date.axis_X_Data_FFT_Input[Flugate_Date.Data_idx], 
        &Flugate_frame.dataBuf[i*12 + 0], 4);										 												                                                                                
        Flugate_Date.Data_idx++;//数据组数										                                                                                 
    }								                                                                                                                         
}

2. 将数据上传到电脑

在 STM32 上,使用串口(USART)、USB、I2C、SPI 或网络通信(如 TCP/IP)将原始数据发送到计算机。

  • 通过串口发送数据:使用 HAL_UART_Transmit() 或直接操作串口接口,发送数据到计算机。
  • 通过 USB:将 STM32 配置为 USB 虚拟串口或 USB 存储设备,将数据通过 USB 上传到计算机。
  • 通过网络(TCP/IP):如果 STM32 具有网络接口(如以太网或 Wi-Fi),则可以通过 TCP/IP 协议将数据传输到计算机。

示例代码(串口发送数据):

comSendBuf(COM4, (uint8_t*)&Flugate_Date.axis_X_Data_FFT_Input[Flugate_Date.Data_idx], 4);

3. MATLAB 接收数据并验证

在 MATLAB 中,接收 STM32 上传的数据,进行验证、处理和可视化:

  • 打开串口:在 MATLAB 中使用 serialport 打开与 STM32 的连接,我这里是通过串口发送数据,所以用到 serialport 来接收数据。
  • 接收数据:使用  read 等函数接收数据。
  • 数据验证与处理:接收的数据可以进行校验、滤波、分析或可视化,确保其正确性。
  • STM2发送的数据是来自传感器,每个数据均为4个字节,共发送了512个数据,对这512个数据进行FFT变换。

示例代码(MATLAB 接收数据):

clear;

serialObj = serialport('COM11',  115200);  % 替换为正确的串口号

% 设置超时时间为 200 秒
serialObj.Timeout = 200;  % 设置超时为 20 秒

% 读取 512 * 4个 字节的数据
numFloats = 512;
% 按字节读取数据
data_UINT8 = read(serialObj, numFloats * 4, "uint8");  

% 将字节数据重新排列为 4 字节一组
% 重塑为 512 行,每行 4 字节
data_UINT8 = reshape(data_UINT8, 4, 512).' ;  

% 更改数据大小端
data_UINT8_change = fliplr(data_UINT8);

% 创建一个 512x1 的全零数组
float_value = zeros(512, 0);  

for index = 1:512
% 将四个字节合并成一个32位无符号整数
int_value = bitshift(data_UINT8_change(index,1), 24) + bitshift(data_UINT8_change(index,2), 16) + bitshift(data_UINT8_change(index,3), 8) + data_UINT8_change(index,4);
% 将该无符号整数转换为float
% int_value 需要先转为int32类型
float_value(index) = typecast(int32(int_value), 'single');  

end
% 关闭串口连接
clear serialObj;

%进行ftft
%采样频率
Fs = 512;
%采样点数
N = 512;
float_value_FFT = fft(float_value,N);

Mag = abs(float_value_FFT);

STM32进行傅里叶代码

具体代码附件可以进行私聊。同时STM32实现FFT可见博客:嵌入式开发傅里叶变换(5):基于STM32-DSP库实现-CSDN博客

结果分析

STM32进行FFT的滤波效果:

Matlab滤波效果:

比较联调可得STM32 FFT变换后的结果和Matlab中基本一致,求出的幅值响应没有问题。


http://www.niftyadmin.cn/n/5868391.html

相关文章

[250224] Yaak 2.0:Git集成、WebSocket支持、OAuth认证等 | Zstandard v1.5.7 发布

目录 Yaak 2.0 发布&#xff1a;Git 集成、WebSocket 支持、OAuth 认证等众多功能&#xff01;Zstandard v1.5.7 发布&#xff1a;性能提升&#xff0c;稳定性增强 Yaak 2.0 发布&#xff1a;Git 集成、WebSocket 支持、OAuth 认证等众多功能&#xff01; Yaak&#xff0c;一款…

vue自定义指令千分位

问题 开发的时候经常会遇到需要在输入框中输入数字转为千分位&#xff0c;点击填写时又转为数字的情况 解决 因此直接在vue中注入自定义指令&#xff0c;通过使用自定义指令达到效果&#xff1b;限制input输入框只能输入数字和一位小鼠带你 自定义指令-千分位 // 自定义指令-千…

双足肌肉骨骼机器人 VS 传统钢铁结构机器人:科技新趋势与跨界创新

引言 近年来&#xff0c;机器人技术突飞猛进&#xff0c;人工智能、自动化、机器人硬件的发展不断刷新我们的认知&#xff0c;尤其是在机器人外骨骼和双足机器人领域。随着波兰某科技公司研发出一款具有肌肉骨骼结构的双足机器人&#xff0c;传统钢铁结构的机器人&#xff0c;尤…

计算机网络:从底层原理到前沿应用,解锁数字世界的连接密码

计算机网络&#xff1a;从底层原理到前沿应用&#xff0c;解锁数字世界的连接密码 在信息如洪流般奔涌的时代&#xff0c;计算机网络宛如无形的脉络&#xff0c;贯穿于我们生活的每一个角落。它不仅是数据传输的通道&#xff0c;更是连接全球、驱动创新的核心力量。从日常的网络…

HWUI 和 Skia

&#x1f4cc; HWUI 和 Skia 的关系 Skia 是 Android 的底层 2D 图形库&#xff0c;提供 CPU 和 GPU 渲染能力&#xff0c;支持 OpenGL、Vulkan、Metal 等后端。HWUI 是 Android UI 组件的 GPU 渲染引擎&#xff0c;主要用于 加速 View、动画、阴影等 UI 元素的绘制。HWUI 依赖…

动态自定义标签属性页面(Tomcat 9)

java文件 &#xff0c;包名org.rain.tag package org.rain.tag; import java.io.IOException; import java.util.HashMap; import java.util.Map; import javax.servlet.jsp.JspException; import javax.servlet.jsp.tagext.DynamicAttributes; import javax.servlet.jsp.…

sage-huga改进SITAN

Sage-Husa自适应滤波算法 Sage-Husa自适应滤波算法是一种在递推滤波过程中实时估计和修正系统噪声和观测噪声统计特性的算法&#xff0c;从而降低系统模型误差&#xff0c;提高滤波精度。该算法基于卡尔曼滤波&#xff0c;并通过自适应调整噪声协方差矩阵来优化滤波效果。 算…

面向对象程序设计-实验十一

&#xff08;给出题目描述&#xff09; 6-1 时钟模拟 代码清单&#xff1a; #include<iostream> using namespace std; class MyTime { private: int m_h; int m_m; int m_s; public: MyTime(int h, int m, int s) { m_h h; m_m m; m_s s; } MyTime() { this-&…