使用STM32提供的DSP库进行FFT(包含板载ADC的例程,AD9220的例程和AD7606的例程)

news/2024/11/7 18:03:09/

使用的芯片为STM32F103,使用的stm32官方提供的FFT库进行FFT处理。
使用的开发板是野火的指南者

一:添加DSP库到自己的工程中
下载得到STM32的DSP库之后,就可以将其添加到自己的工程项目中了。

其中,inc文件夹下的stm32_dsp.h和table_fft.h两个文件是必须添加的。stm32_dsp.h是STM32的DSP库的头文件。

src文件夹下的文件可以有选择的添加(用到那个添加那个即可)。因为我只用到了256点的FFT,所以这里我只添加了cr4_fft_256_stm32.s文件。添加完成后的项目框架如图所示。
  在这里插入图片描述
二.使用DSP库的几点问题
进行256点的FFT,只需要调用STM32 DSP库函数中的cr4_fft_256_stm32()函数即可。该函数的原型为:

void cr4_fft_256_stm32(void *pssOUT, void *pssIN, uint16_t Nbin);

其中,参数pssOUT表示FFT输出数组指针,参数pssIN表示要进行FFT运算的输入数组指针,参数Nbin表示了点数。
  下面是具体的调用实例:

cr4_fft_256_stm32(lBufOutArray, lBufInArray, NPT);

其中,参数lBufOutArray同样是一个long类型的数组,参数lBufInArray就是存放模拟采样数据的采样数组,同样也是long型的数组,NPT为采样点数256。
  
  在输入参数的时候,lBufInArray这个数组的高十六位是信号时域采样值的实数,低十六位是虚部,一般只需要使用实部即可。

调用该函数之后,在lBufOutArray数组中就存放了进行FFT运算之后的结果数据。该数组中每个元素的数据格式为;高16位存储虚部,低16位存储实部。

以上文字来源:https://www.cnblogs.com/menlsh/p/4154070.html

有了上面的准备FFT就很简单了,只要把ADC的数据存到所需要的数组里面即可。

第一个ADC使用的是STM32自带的ADC采集数据,采样速率根据的是
在这里插入图片描述
配置ADC时钟,用DMA存入数组之后再对数组进行统一转换。
输出结果的程序是:
在这里插入图片描述
第i个点对应的频率点可以使用公式:
片上ADC转换一次的时间(根据数据手册)为:
1总采样周期:设置的采样周期(1.5或7.5或13.5或28.5或41.5或55.5或71.5或239.5)+12.5个采样周期
2.转换时间=总采样周期/ADC时钟分频后的频率
3.【256/(1/转换时间)】× i = 所对应的频率(单位为赫兹)

源码下载:https://download.csdn.net/download/asukadesu/12924414

AD9220+FFT源码下载(基于康威科技):https://download.csdn.net/download/asukadesu/12929968

AD7606的使用:

基于康威科技的例程,但是由于康威科技的例程是基于八通道的采样,使用起来不方便进行采样时间的控制。因此进行了单通道的微调。

1.7606的使用:
根据数据手册
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
因此从程序上,先开机让他一直进行转换,然后设定定时器中断读取数据存入数组里即可

void TIM4_IRQHandler(void)
{ad7606_IRQSrc();
}
void ad7606_IRQSrc(void)
{uint8_t i;uint16_t usReadValue;int32_t temp;TIM_ClearFlag(TIM4, TIM_FLAG_Update);AD_CS_LOW();usReadValue = ad7606_ReadBytes();	AD_CS_HIGH();temp = ((int32_t)10000)*((float)((short)usReadValue)/65536);	//单位1mv,除以65536是±5V,除以32768是±10Vprintf( "Range%d = %d\n", i, temp);ad7606_StartConv();
}

对于这条公式

temp = ((int32_t)10000)*((float)((short)usReadValue)/65536);	

在数据手册中,采样值用带符号的16位来存储,因此有效的数位是15位,即32768为最大值,前面的【(int32_t)10 000】实际就是最大量程10V,要调成5V既可以改这里10 000变成5000,也可以把65536改成32768。

此外,在非实验室条件下要产生负压其实也很简单,只要用单片机随便给一个单片机供电,然后把AD7606的地接3.3,输入接地即可。

AD7606测试结果:在这里插入图片描述
在这里插入图片描述

AD7606的源码下载:https://download.csdn.net/download/asukadesu/12931633


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

相关文章

lan 9220

http://www.alldatasheet.com/datasheet-pdf/pdf/312187/SMSC/LAN9220-ABZJ.html gpmc和fpga通信: http://blog.csdn.net/zhangxiao93/article/details/40650829

AT32(六)TMR输出时钟和DMA请求实现两片AD9220同步采样

前言 有个小项目需要用到两片AD9220——10MSPS并口12位ADC——做同步采样,想到AT32F435可以通过定时器实现让DMA按照特定频率从GPIO数据寄存器搬运数据到内存,于是把这个想法实现一下。 软件/硬件环境 VScode ( EIDE Cortex Debug )Open On-Chip Deb…

C++ 第一弹入门基础

目录 目录 1.关键字 2.命名空间 3.标准IO输入输出 4.缺省参数 5.函数重载 6.引用 7.内联函数 1.关键字 c98的关键字一共有63个 在之后用的时候都会再次详细介绍 2.命名空间 2.1什么是命名空间,为什么要有他? 在c中为了避免相同名字的变量&am…

Win 8 RT: 你应该知道的15件事

这次的Win 8直接来了一个Win 8 RT, Win 8 Pro兄弟版。很多用户只知道Win 8 RT为 ARM 硬体而设,比 Win 8 Pro价格便宜,比 Win 8 Pro 薄一些,但并不明白两者的本质区别。而下面我们将用15个问题来让你了解什么是Win 8 RT。 1. Win 8…

关于 Win 8 RT 你应该知道的 15 件事

2019独角兽企业重金招聘Python工程师标准>>> 这次的Win 8直接来了一个Win 8 RT, Win 8 Pro兄弟版。很多用户只知道Win 8 RT为 ARM 硬体而设,比 Win 8 Pro价格便宜,比 Win 8 Pro 薄一些,但并不明白两者的本质区别。而下…

专项练习9

目录 一、选择题 1、在 JavaScript 中,用于阻止默认事件的默认操作的方法是 2、以下代码执行后,result 的值为() 3、不能从字符串 const str qwbewrbbeqqbbbweebbbbqee;中能得到结果 ["b", "bb", "bbb…

网络安装Fuji Xerox DocuCentre-III 2007打印机

新装了Window2012Server,准备装上打印机做打印服务器,结果怎么安装都不能打印,看来是出问题了,折腾了好多次,原来是我选错了,必须是TCP/IP端口。 在这里输入打印机的IP地址,如果不知道IP地址要通…

松下dp-8016p打印机-网络打印设置

第一步:设置打印机上IP地址 按打印机面板的“功能”键,再按“1”一般功能。然后按“设定”,再按数字键9,一般参数,要求输入密码屏幕提示,键操作员模式输入000,之后按设定,进入键操作…