DAC-ADC-DSP库FFT获取正弦波频率

已于 2023-07-26 17:51:38 修改
阅读量1.4k 收藏 27
点赞数 8
文章标签： stm32
于 2023-07-24 20:54:10 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_63636870/article/details/131904634
版权

DSP与FFT

DSP库的添加

cubemx配置
下载DSP库

在这里插入图片描述
在这里插入图片描述

添加DSP库到工程中

在这里插入图片描述在这里插入图片描述

工程配置
添加宏

在这里插入图片描述

ARM_MATH_CM4,__CC_ARM,ARM_MATH_MATRIX_CHECK,ARM_MATH_ROUNDING

第一个宏：

M3:ARM_MATH_CM3

M4:ARM_MATH_CM4

H7:ARM_MATH_CM7

include

在这里插入图片描述

#include "arm_math.h"
#include "math.h"
#include "arm_const_structs.h"

#include "arm_const_structs.h" 如果cubemx生成的库没有，可以找板子的例程，将DSP的lib文件添加到工程中
在这里插入图片描述

测试
/* USER CODE BEGIN PV */
uint8_t lcd_id[12];				//存放LCD ID字符串
float data = 0;//定义一个float型变量
float fft_inputbuf[FFT_LENGTH*2];	//FFT输入数组
float fft_outputbuf[FFT_LENGTH];	//FFT输出数组
/* USER CODE END PV */

/* USER CODE BEGIN 1 */
	arm_cfft_radix4_instance_f32 scfft;
	uint8_t key,t=0;
	float time; 
	uint8_t buf[50]; 
	uint16_t i; 
	char send[15];
/* USER CODE END 1 */

/* USER CODE BEGIN 2 */
	LCD_Init();
	POINT_COLOR=RED;      //画笔颜色：红色
	sprintf((char*)lcd_id,"LCD ID:%04X",lcddev.id);//将LCD ID打印到lcd_id数组。		
	LCD_ShowString(30,40,210,24,24,"mcudev STM32F4");	
	LCD_ShowString(30,70,200,16,16,"TFTLCD TEST");
	LCD_ShowString(30,90,200,16,16,"mcudev.taobao.com");
	arm_cfft_radix4_init_f32(&scfft,FFT_LENGTH,0,1);//初始化scfft结构体，设定FFT相关参数
  /* USER CODE END 2 */

/* USER CODE BEGIN WHILE */
  while (1)
  {
	  data = arm_sin_f32(PI/6);					//计算sin30°
	  sprintf((char*)lcd_id,"LCD ID:%.5f",data);//0.49999
	  LCD_ShowString(30,110,200,16,16,lcd_id);
		for(i=0;i<FFT_LENGTH;i++)//生成信号序列
		{
			 fft_inputbuf[2*i]=100+
							   10*arm_sin_f32(2*PI*i/FFT_LENGTH)+
							   30*arm_sin_f32(2*PI*i*4/FFT_LENGTH)+
							   50*arm_cos_f32(2*PI*i*8/FFT_LENGTH);	//生成输入信号实部
			//printf("%.2f\r\n",fft_inputbuf[2*i]);
			sprintf(send,"%.2f\r\n",fft_inputbuf[i*2]);
			HAL_UART_Transmit(&huart1,(uint8_t*)send,8,0xff);
			 fft_inputbuf[2*i+1]=0;//虚部全部为0
		}
		printf("\r\n");
		arm_cfft_radix4_f32(&scfft,fft_inputbuf);	//FFT计算（基4）
		printf("after:");
		for(i=0;i<FFT_LENGTH;i++)//生成信号序列
		{
			sprintf(send,"%.2f\r\n",fft_inputbuf[i*2]);
			HAL_UART_Transmit(&huart1,(uint8_t*)send,8,0xff);
		}
		arm_cmplx_mag_f32(fft_inputbuf,fft_outputbuf,FFT_LENGTH);	//把运算结果复数求模得幅值
		printf("\r\n%d point FFT runtime:%0.3fms\r\n",FFT_LENGTH,time/1000);
		printf("FFT Result:\r\n");
		for(i=0;i<FFT_LENGTH;i++)
		{
			printf("fft_outputbuf[%d]:%f\r\n",i,fft_outputbuf[i]);
		}		
    /* USER CODE END WHILE */

如果没有LCD

/* USER CODE BEGIN PV */
float data = 0;//定义一个float型变量
float fft_inputbuf[FFT_LENGTH*2];	//FFT输入数组
float fft_outputbuf[FFT_LENGTH];	//FFT输出数组
/* USER CODE END PV */

/* USER CODE BEGIN 1 */
	arm_cfft_radix4_instance_f32 scfft;
	uint8_t key,t=0;
	float time; 
	uint8_t buf[50]; 
	uint16_t i; 
	char send[15];
/* USER CODE END 1 */

/* USER CODE BEGIN 2 */
	arm_cfft_radix4_init_f32(&scfft,FFT_LENGTH,0,1);//初始化scfft结构体，设定FFT相关参数
  /* USER CODE END 2 */

/* USER CODE BEGIN WHILE */
  while (1)
  {
	  data = arm_sin_f32(PI/6);					//计算sin30°
      //开调试看data值或串口输出
		for(i=0;i<FFT_LENGTH;i++)//生成信号序列
		{
			 fft_inputbuf[2*i]=100+
							   10*arm_sin_f32(2*PI*i/FFT_LENGTH)+
							   30*arm_sin_f32(2*PI*i*4/FFT_LENGTH)+
							   50*arm_cos_f32(2*PI*i*8/FFT_LENGTH);	//生成输入信号实部
			//printf("%.2f\r\n",fft_inputbuf[2*i]);
			sprintf(send,"%.2f\r\n",fft_inputbuf[i*2]);
			HAL_UART_Transmit(&huart1,(uint8_t*)send,8,0xff);
			 fft_inputbuf[2*i+1]=0;//虚部全部为0
		}
		printf("\r\n");
		arm_cfft_radix4_f32(&scfft,fft_inputbuf);	//FFT计算（基4）
		printf("after:");
		for(i=0;i<FFT_LENGTH;i++)//生成信号序列
		{
			sprintf(send,"%.2f\r\n",fft_inputbuf[i*2]);
			HAL_UART_Transmit(&huart1,(uint8_t*)send,8,0xff);
		}
		arm_cmplx_mag_f32(fft_inputbuf,fft_outputbuf,FFT_LENGTH);	//把运算结果复数求模得幅值
		printf("\r\n%d point FFT runtime:%0.3fms\r\n",FFT_LENGTH,time/1000);
		printf("FFT Result:\r\n");
		for(i=0;i<FFT_LENGTH;i++)
		{
			printf("fft_outputbuf[%d]:%f\r\n",i,fft_outputbuf[i]);
		}		
    /* USER CODE END WHILE */

DAC

cubemx配置
DAC

在这里插入图片描述

  • Trigger:Timer2 Trigger Out event

在这里插入图片描述

DMA点ADD添加，Circular模式

在这里插入图片描述

确认配置无误

TIM2

在这里插入图片描述

  • 时钟频率为168MHz
  • PSC取168-1
  • ARR取5-1

输出频率为：168000000/168/5 = 200000Hz

工程配置
BEGIN0 产生正弦波函数

在这里插入图片描述

/* USER CODE BEGIN 0 */
void SineWave_Data( uint16_t num,uint16_t *D,float U)
{
    uint16_t i;
    for( i=0;i<num;i++)
    {
        D[i]=(uint16_t)((U*sin(( 1.0*i/(num-1))*2*3.14159265358979)+U)*4095/3.3);
    }
}

/* USER CODE END 0 */
BEGIN1 正弦波数组

在这里插入图片描述

 /* USER CODE BEGIN 1 */
	uint16_t Sine12bit[32] = {
		2048	, 2460	, 2856	, 3218	, 3532	, 3786	, 3969	, 4072	,
		4093	, 4031	, 3887	, 3668	, 3382	, 3042	,2661	, 2255	, 
		1841	, 1435	, 1054	, 714	, 428	, 209	, 65	, 3		,
		24		, 127	, 310	, 564	, 878	, 1240	, 1636	, 2048
	};
	
  /* USER CODE END 1 */
BEGIN2 初始化

在这里插入图片描述

//DAC
	TIM2->ARR = 9;						//修改ARR 可不加
	HAL_TIM_Base_Start(&htim2);			//打开定时器中断
	//HAL_DAC_Start(&hdac,DAC_CHANNEL_1);
	SineWave_Data(n,DualSine12bit,1.6);	//产生正弦波函数
	HAL_DAC_Start_DMA(&hdac, DAC_CHANNEL_1, (uint32_t *)Sine12bit, 32, DAC_ALIGN_12B_R);
输出信号

在这里插入图片描述

ADC

cubemx配置
ADC1

在这里插入图片描述
在这里插入图片描述

在这里插入图片描述

  • 定时器3触发
  • 添加DMA
  • 打开中断
TIM3

在这里插入图片描述

在这里插入图片描述

  • 时钟频率168MHz
  • PSC:168-1
  • ARR:2-1
  • Update Event
  • 打开中断
工程配置
BEGIN PD
#define fft_adc_n 1024    // 采1024个点
BEGIN PV
float adc_data[fft_adc_n*2]={0};        // 存ADC值
uint8_t ad_flag = 0;
BEGIN2

在这里插入图片描述

HAL_TIM_Base_Start_IT(&htim3);				//打开定时器3中断
HAL_ADC_Start_DMA(&hadc1,(uint32_t*)&adc_buffer,fft_adc_n);	//ADC采集一组数据
BEGIN4
void HAL_ADC_ConvCpltCallback(ADC_HandleTypeDef* hadc)
{
    ad_flag=1;
}

ADC中断，采集完成后标志位置1

while(1)
while (1)
  {
	  if(ad_flag == 1)
	  {
          //数据处理
          //......
          
          ad_flag = 0;	//标志位置0
      }
}

FFT

float fft_in_adc_data[fft_adc_n*2]={0}; // FFT输入，实部是ADC值，虚部补0
float fft_out_adc_data[fft_adc_n]={0};  // FFT输出
uint16_t flag[fft_adc_n] = {0};			// 数据处理时暂时存放数据
int k = 1;  							// 找第k大值
float32_t maxValue = 0;					// ADC最大值
uint32_t index = 0;						// 最大值下标
int i,j;
float freq = 0;							// 频率

arm_cfft_radix4_instance_f32 scfft;        //FFT对应结构体变量
arm_cfft_radix4_init_f32(&scfft,fft_adc_n,0,1); //初始化scfft结构体，设置FFT相关参数

if(ad_flag == 1)
{
    printf("finish\r\n");
    ad_flag = 0;
    for(i=0;i<fft_adc_n;i++)
    {
        fft_in_adc_data[i*2] = adc_buffer[i]*3.3/4095;		//实部
        fft_in_adc_data[i*2+1] = 0;							//虚部
    }
    printf("\r\n");
    arm_cfft_f32(&arm_cfft_sR_f32_len1024,fft_in_adc_data,0,1);
    arm_cmplx_mag_f32(fft_in_adc_data,fft_out_adc_data,fft_adc_n);
		  
    fft_out_adc_data[0]/=1024;
    for(i=1;i<fft_adc_n;i++)
        fft_out_adc_data[i]/=512;
    for(i=0;i<fft_adc_n;i++)
    {
        printf("%.5f\t",fft_out_adc_data[i]);
    }
    printf("\r\n");
    //arm_max_f32(&fft_out_adc_data[1],fft_adc_n/2-1,&maxValue,&index);

    /****查找最大值及其下标*********/
	memset(flag, 0, sizeof(flag));
  	for (i = 1; i <= k; i++)
  	{
    	maxValue = 0.0;
    	for (j = 1; j < fft_adc_n / 2; j++)
    	{
            if (flag[j] == 1)
                break;
            if (fft_out_adc_data[j] > maxValue)
            {
                maxValue = fft_out_adc_data[j];
                index = j;
            }
		}
    	flag[index] = 1;
  	}
    printf("index:%d\r\n",index);
    printf("Vmax:%.2f\r\n",maxValue);
    freq = (float)500000/fft_adc_n * (index) /2;		//500000为采样率
    
    printf("freq: %.3f",freq);
    ad_flag = 0;
}
  • arm_cfft_sR_f32_len1024在arm_const_structs.h里面，如果报错，检查头文件
  • 如果更改采样点数，注意修改arm_cfft_sR_f32_len1024，如果采样512个点，则改为arm_cfft_sR_f32_len512，以此类推
  • fft_out_adc_data[0]为直流分量电压值
  • maxValue为基波电压幅值
  • fft_out_adc_data[0]+maxValue=正弦波最大值
  • fft_out_adc_data[0]-maxValue=正弦波最小值
参考文章

https://blog.csdn.net/Nothing_To_Say_/article/details/123606260

https://blog.csdn.net/qq_34022877/article/details/123190943

https://blog.csdn.net/qq_24426625/article/details/129738537

等

m0_63636870
关注
DAC、ADC、FFT使用总结
及时总结
07-28 3289
计算公式， 波形生成， DAC波形频率， ADC采样时间， 离散傅里叶变换DFT， FFT
烧脑长文！近乎完美的DDS正弦波信号音生成器！
qq_43416206的博客
03-16 996
在标准评估板上实施的初级和开箱即用测试显示，用于传统正弦波CW生成的基于处理器的DDS技术要实现高性能指日可待。通过精心设计重构滤波器和模拟输出缓冲级，可以实现–120dBc谐波失真系数。基于DSP的NCO/DDS不只受到单信号音正弦波生成限制。通过使用具备合适的截止频率，且无其他硬件变更的优化AAF（贝塞尔或巴特沃兹），同样的DSP和DAC组合可用作高性能AWG来生成任何类型的波形，例如，完全合成可设置参数的多信号音正弦波（可以完全控制每个分量的相位和幅度）来实施IMD测试。
ADC采样——fft运算的前提
qq_37203178的博客
07-28 3955
1、基本介绍 2、ADC引脚输入一般控制在0~3.3V，超过3.3V可能回烧掉。 VSSA= Vref- =GND 2.4V<= Vref+ <=VDDA=3.3V 3、注入通道 和 规则通道 如果只有规则通道就按照通道排序进行运行，但是若是在规则通道中添加了注入通道，此时运行到插入注入通道的位置时，必须先执行注入通道，注入通道执行结束后再回到之前执行后面的规则通道。如下所示， 4、注意1：STM32F1的ADC是支持单次，连续，扫描或者间断模...
dsp输出正弦波 c语言,快速计算正弦波 - DSP 与 ARM 处理器 - 电子工程世界-论坛...
weixin_42512699的博客
05-22 2318
快速计算正弦波[复制链接]在DSP运用中，经常需要产生正弦波。如果直接用c的数学函数sin，当然可以产生正弦波，但是由于sin函数本身的效率很低，产生正弦波所需要的MIPS就会占去DSP处理能力的相当大的一部分。本章介用递推数列算正弦波的方法，先介绍原理，推导出递推公式，然后用浮点小数实现计算，再用定点小数进一步优化算法，最后进行误差分析，并提出更精确的定点小数算法。先来看看如何推导出递推数列的公...
总结——TI_音频信号分析仪
最新发布
m0_74349248的博客
07-23 1433
使用MSPM0G3507开发音频信号分析仪，介绍了基本的操作
利用adc采样正弦波，计算频率的c语言程序
weixin_42601608的博客
02-17 1339
我不了解你问的问题，但是我可以提供一些有关编程的关键思想。ADC采样正弦波的关键就是要在一定的时间内记录每一次采样的结果，然后再根据这些采样结果计算出频率。你可以使用C语言中的循环结构来实现这一过程，比如使用for循环来完成采样结果的记录，然后再使用函数计算出频率。 ...
如何 FFT(快速傅里叶变换) 求幅度、频率（超详细 含推导过程）
热门推荐
Xav Pun的博客
11-30 10万+
保姆级教学➕全部代码➕原理解释➕推导过程
用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
stm32adc采集正弦波并用stm32dsp库的fft处理数据计算失真度
06-25
本文旨在介绍如何使用STM32的ADC采集正弦波并使用STM32的DSP库的FFT处理数据计算失真度。 第一步，需要确定要采集的正弦波的频率和振幅范围。可以使用STM32的定时器和DAC模块生成一个可调节的正弦波信号。将它的...
基于MATLAB和DSP的数字信号处理课程实验设计-综合文档
05-19
1. **信号生成**：MATLAB提供了丰富的函数库来创建各种类型的数字信号，如正弦波、方波、白噪声等。实验者可以利用这些函数生成模拟信号，以便进行后续的处理和分析。 2. **滤波器设计**：在数字信号处理中，滤波器...
STM32F4 ADC采样FFT运算测试代码
07-29
模拟信号经过ADC采样后变成数字信号，数字信号可以进行FFT运算，在频域中更容易分析信号的特征。此代码用STM32F407的ADC-DMA模式采集4096个点的数据，利用DSP库里的FFT算法进行快速傅里叶变换，经实测可以使用。
基于DSP芯片设计的一种波形发生器-综合文档
05-22
1. **信号生成模块**：DSP芯片内部的高速计算能力使得它可以执行复杂的数学运算，例如傅立叶变换，用于生成正弦、方波、三角波等基本波形，甚至复杂的脉冲序列或任意波形。这通常涉及到查表法（Look-up Table, LUT）...
STM32实现FFT，求取幅度频谱
qq_34022877的博客
02-28 4万+
STM32实现FFT，求取幅度频谱 FFT不太对劲的理解 FFT的原理比较复杂，因为32使用FFT不用去管算法是如何运作的，我在这里就进行简单的介绍了。 因为是简单介绍，就只介绍下幅度频谱图，不考虑相位频谱图。 ​ FFT可以将一个信号从时域变换到频域，比如一个1VPP的1k的正弦信号，它的时域和频域的示意图如下： ​ 频域为我们观察信号提供了一个新的视角。比如下面是1k和2k信号的叠加。 ​ 从时域上看，1k+2k的波形不容易进行处理，也不好猜出来这个波形到底有什么特性（当然这个例子其实还是
STM32F407FFT识别正弦波信号的频率
GYX0518的博客
09-23 3127
正弦信号进行FFT变换后会有较为明显的特征，在所对应的频率处会有一个较大的值，这个值所对应的频率即为正弦信号的频率，这里是以10K的正弦信号为例。运用单片机的DAC口或DDS对识别的信号进行复现，使用DAC口输出信号需要用硬件电路进行处理才能得到稳定的信号。在STM官网下载407系列单片机的DSP固件库，在Keil中添加相关的库。主函数中的max1、max2的大小乘以1K即为输入正弦信号的频率大小。3.对频谱分析得的结果进一步分析，得到输入信号正弦波的频率大小。这里用的TFT屏进行显示。
STM32F1 标准库 ADC正弦波采样+DMA+FFT+输出正弦波频率（超详细版）
weixin_75086069的博客
10-22 1万+
本文代码可以实现用单片机对一个未知频率的正弦波进行采样，经过运算后得到该频率。实现方法为：通过设置采样频率，对采样到的数据进行传输，使用对采样到的波形数据进行傅里叶变换，得到傅里叶序列，对其进行计算得到，当幅值最大时对应的序列就是该未知波形的频率（序列需要乘分辨率）。
开源小项目-基于STM32F103的频率测量实现 V1.1（ADC+TIM+DMA+FFT）
JJ_KING123的博客
05-09 1万+
该项目是基于STM32F103的频率测量，包括频率、振幅、相位、总谐波失真等信息，采用FFT快速傅里叶变换，配合ADC、DMA、TIM、FFT实现。 快速傅里叶变换（Fast Fourier Transform，简称FFT）是一种高效的算法，用于计算离散傅里叶变换（Discrete Fourier Transform，简称DFT）及其逆变换。傅里叶变换是一种将信号从时域（或空间域）转换到频域的数学方法。这意味着它可以将复杂的波形分解成一系列不同频率的正弦波和余弦波的组合，每个波都有其对应的振幅和相位。
【嵌入式】利用arm-DSP库进行FFT计算，获得信号的频谱、幅值及相位(上)
spiremoon的博客
08-24 2万+
电力系统中往往掺杂谐波，而FFT可以将谐波检测出来，具有较大的实用价值。今天主要讲一下在STM32中如何利用dsp库进行快速傅里叶计算，从而得出信号的频谱幅值以及相位。 一、Matlab简单搭建 1.谐波检测搭建 通过50Hz正弦波叠加100Hz、200Hz的正弦波得到叠加后畸变的波形，这个波形主要用来进行FFT运算进行谐波检测。模型很简单如图1所示，三个正弦波相位相同、频率不同，幅值分别为5、1...
添加DSP库启用arm_math.h
ittuann的博客
06-23 4599
开启STM32的FPU硬件单精度浮点数计算单元，并启用CMSIS封装的STM32的DSP静态库。像是卡尔曼滤波的矩阵运算，FIR和IIR滤波还有很多高级数学功能有了DSP库都会很方便和快速。FPU：Float Point Unit。全称是浮点计算单元，用于浮点计算的协处理器。一般单片机带有硬件FPU并且支持浮点指令集时，工程会默认开启了FPU功能。STM32F4就带有32位单精度硬件FPU，支持浮点指令集。在头文件中 中定义了。其值为时可用。在头文件中 中定义了。其值为时可用。在中初始化函数会根据上面两个
STM32CubeMX关于添加DSP库的使用
W_Oilpicture的博客
07-04 8929
数字信号处理（，简称 ） 是一门涉及许多学科而又广泛应用于许多领 域的新兴学科，通过利用计算机或专用处理设备，以数字形式对信号进行采集、变换、滤波、估值、增强、 压缩、识别等处理，以得到符合人们需要的信号形式。 库中包含常用数学运算，复数，矩阵，三角函数，插值、求最大值求最小值、还有重要的滤波器和 ，非常实用，其可以调用的函数十分之多。 完成上述步骤即完成 库的添加。之后就可以生成代码了。 5.代码中添加头文件 以上步骤已完成 库的添加。下面进行调用一个 库的函数进行测试代码： 7.测试
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值