[其他ST产品] STM32F4 Timer+ADC+DMA+FFT的理解与应用

[复制链接]
4448|77
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:00 | 显示全部楼层 |阅读模式
ST, DMA, DM, ADC, AD
STM32F407 利用Timer+ADC+DMA在FFT上的理解与运用
最近遇到有关信号的总谐波失真THD值的测量,根据THD的计算公式可知我们需要知道该信号各次谐波分量的幅值,所以需要将ADC采样到的信号进行FFT,将时域上的信号转换到频域上,获得其幅值。

这里我们假设要采集的信号为1Khz。
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:00 | 显示全部楼层
方案
利用STM32F407上的定时器Timer来触发ADC采样,并利用DMA搬运采样到的AD值,最后用dsp库里的有关FFT运算的函数进行各次谐波幅值的获取。 6036764eb1eb56a652.png
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:00 | 显示全部楼层
原理
设定数据
  1. 被采信号频率:f
  2. 被采信号周期:T = 1/f
  3. 采样频率:fs
  4. 采样周期:Ts = 1/fs

  5. 采样总点数:NPT
  6. 采样总时间:t = NPT * Ts

  7. 频谱图频率分辨率:f0 = fs/NPT

  8. 采到的被采信号周期数:NT = t/T
  9. 一个被采信号周期的采样点数:fs/f
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:00 | 显示全部楼层
带入数据
我们需要采的信号频率为1KHz,此次准备以32KHz的采样频率进行采集,并且一共只采集256个点进行FFT运算。

故:
  1. 被采信号频率:f = 1000Hz
  2. 被采信号周期:T = 1/f  = 0.001s
  3. 采样频率:fs = 32 000Hz
  4. 采样周期:Ts = 1/fs =0.000 031 25s

  5. 采样总点数:NPT = 256
  6. 采样总时间:t = NPT * Ts = 0.008s

  7. 频谱图频率分辨率:f0 = fs/NPT = 32 000Hz/256 = 125Hz

  8. 采到的被采信号周期数:NT = t/T = 0.008s/0.001s = 8
  9. 一个被采信号周期的被采到的点数:fs/f = 32 000Hz/1000Hz = 32
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:00 | 显示全部楼层
分析数据
FFT后的幅度谱的横坐标是频率,并且是离散的。是以频率f0的n倍展开的。故横坐标为

… , -nf0, …-2f0, -f0, 0, f0, 2f0, … ,nf0, … 而纵坐标为对应频点的幅值信息。
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:01 | 显示全部楼层
本次是设定的f0 = 125Hz, 理想情况下FFT后的幅度谱看起来应该会是下面这个样子: 1523764eb1ee75b535.png
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:01 | 显示全部楼层
其中横坐标0刻度的地方即频率为0的点,为直流分量 。如果被采信号没有偏置的话这一点幅度应该为0。而横坐标其它点皆为125Hz的倍数。
由于被采信号的频率为1000Hz,故横坐标的正半轴和负半轴的1000Hz处都会有“擎天柱”。但在运用中我们只会取正半轴部分,毕竟正半轴知道了也就知道了负半轴,所以负半轴可以说是没什么用的。
值得注意的是:对于直流分量来说,纵坐标对应的值就是直流分量的幅度,而其他频率对应的纵坐标值只是其实际幅值的1/2
关于这一点有一个简单的理解方式,就是它的幅度平均分到了正负半轴,导致只有一半。
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:01 | 显示全部楼层
如果从公式的角度,可以这样理解: 6338864eb1f046af76.png 可以看到,正负两边均只有信号峰值A的一半,对于幅度而言也是一样的也是一半。这是本人个人的理解方式,可能会不严谨。
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:01 | 显示全部楼层
代码
ADC+DMA部分:
ADC注意配置成外部时钟触发,不连续转换,单通道不扫描。
DMA需要外设不自增,内存自增的方式来存储采到的连续的256个点,非循环模式。
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:02 | 显示全部楼层
  1. #include "ADCDMA.h"

  2. //PF3 ADC3:IN9

  3. uint16_t AD3_Value[AD3_Value_Length];

  4. void AD3_Init(){
  5.         //结构体
  6.         GPIO_InitTypeDef  GPIO_InitStructure;                                 //GPIO结构体       
  7.         DMA_InitTypeDef DMA_InitStructure;                                        //DMA结构体       
  8.         ADC_CommonInitTypeDef ADC_CommonInitStructure;                 //ADCcommon结构体
  9.         ADC_InitTypeDef ADC_InitStructure;                                        //ADC结构体
  10.        
  11.         //时钟开启
  12.         RCC_AHB1PeriphClockCmd(RCC_AHB1Periph_GPIOF, ENABLE);        //使能GPIOF时钟
  13.         RCC_AHB1PeriphClockCmd(RCC_AHB1Periph_DMA2,ENABLE);                //使能DMA2时钟
  14.         RCC_APB2PeriphClockCmd(RCC_APB2Periph_ADC3, ENABLE);        //使能ADC3时钟
  15.        
  16.         //GPIO配置
  17.         GPIO_InitStructure.GPIO_Pin = GPIO_Pin_3;               
  18.         GPIO_InitStructure.GPIO_Mode = GPIO_Mode_AN;                        //模拟模式
  19. //        GPIO_InitStructure.GPIO_OType = GPIO_OType_PP;                        //推挽输出
  20.         GPIO_InitStructure.GPIO_Speed = GPIO_Speed_100MHz;                //100MHz
  21.         GPIO_InitStructure.GPIO_PuPd = GPIO_PuPd_NOPULL;                //浮空
  22.         GPIO_Init(GPIOF, &GPIO_InitStructure);                                        //GPIO初始化
  23.        
  24.         //ADCcommon配置
  25.         ADC_CommonInitStructure.ADC_Mode = ADC_Mode_Independent;                                        //独立模式
  26.         ADC_CommonInitStructure.ADC_DMAAccessMode = ADC_DMAAccessMode_Disabled;                //非多重模式,多重模式下才开启此配置的DMA
  27.         ADC_CommonInitStructure.ADC_Prescaler = ADC_Prescaler_Div4;                                        //采样频率4分频                84MHz/4 = 21MHz
  28.         ADC_CommonInitStructure.ADC_TwoSamplingDelay = ADC_TwoSamplingDelay_5Cycles;
  29.         ADC_CommonInit(&ADC_CommonInitStructure);                                                                        //ADCcommon结构体初始化
  30.        
  31.         //ADC配置
  32.         ADC_InitStructure.ADC_ContinuousConvMode = DISABLE;                                                        //是否连续转换
  33.         ADC_InitStructure.ADC_DataAlign = ADC_DataAlign_Right;                                                //数据对齐:右对齐
  34.         ADC_InitStructure.ADC_ExternalTrigConv = ADC_ExternalTrigConv_T2_TRGO;//
  35.         ADC_InitStructure.ADC_ExternalTrigConvEdge = ADC_ExternalTrigConvEdge_RisingFalling   ;        //
  36.         ADC_InitStructure.ADC_NbrOfConversion = AD3_Value_Length;                                        //转换数量:一波采集采集的AD值个数,多通道时一般为通道数量
  37.         ADC_InitStructure.ADC_Resolution = ADC_Resolution_12b;                                                //12bit
  38.         ADC_InitStructure.ADC_ScanConvMode = DISABLE;                                                                  //是否扫描:多通道需扫描
  39.         ADC_Init(ADC3, &ADC_InitStructure);                                                                                        //ADC3初始化


  40.         //ADC序列 转换顺序
  41.         ADC_RegularChannelConfig(ADC3,ADC_Channel_9, 1,ADC_SampleTime_84Cycles);


  42.         //DMA配置
  43.         DMA_InitStructure.DMA_Channel = DMA_Channel_2;                                                        //DMA_CH2       
  44.         DMA_InitStructure.DMA_Mode = DMA_Mode_Normal;                                                         //非循环
  45.         DMA_InitStructure.DMA_DIR = DMA_DIR_PeripheralToMemory;                                        //外设到内存
  46.         DMA_InitStructure.DMA_BufferSize = AD3_Value_Length;                                        //传输次数,数组长度
  47.         DMA_InitStructure.DMA_Priority = DMA_Priority_Medium;                                        //优先级
  48.                 //DMA_FIFO
  49.         DMA_InitStructure.DMA_FIFOMode = DMA_FIFOMode_Disable;
  50.         DMA_InitStructure.DMA_FIFOThreshold = DMA_FIFOThreshold_HalfFull;
  51.        
  52.         DMA_InitStructure.DMA_Memory0BaseAddr =        (uint32_t)AD3_Value;                                 //内存地址:收集AD值得数组               
  53.         DMA_InitStructure.DMA_MemoryBurst = DMA_MemoryBurst_Single;
  54.         DMA_InitStructure.DMA_MemoryDataSize = DMA_MemoryDataSize_HalfWord;                        //16 bit
  55.         DMA_InitStructure.DMA_MemoryInc = DMA_MemoryInc_Enable;                                                //内存自增

  56.         DMA_InitStructure.DMA_PeripheralBaseAddr = (uint32_t)&(ADC3->DR);                        //外设地址,ADC3地址,多通道但仅有一个寄存器
  57.         DMA_InitStructure.DMA_PeripheralBurst = DMA_PeripheralBurst_Single;
  58.         DMA_InitStructure.DMA_PeripheralDataSize = DMA_PeripheralDataSize_HalfWord;        //16 bit
  59.         DMA_InitStructure.DMA_PeripheralInc = DMA_PeripheralInc_Disable;                         //外设不自增,始终为ADC3地址
  60.         DMA_Init(DMA2_Stream0,&DMA_InitStructure);                                                                        //DMA2_Stream0 初始化

  61.         DMA_Cmd(DMA2_Stream0,ENABLE);                                                //DMA2_Stream0 使能       
  62.         ADC_DMARequestAfterLastTransferCmd(ADC3, ENABLE);        //源数据变化时开启DMA传输
  63.         ADC_DMACmd(ADC3, ENABLE);                                                        //ADC3_DMA使能
  64.         ADC_Cmd(ADC3, ENABLE);                                                                //ADC3 使能
  65.        
  66. //        ADC_SoftwareStartConv(ADC3);                                                //软件触发ADC转换
  67. }

  68. //******清空ADC与DMA的中断标志位,起到再次触发ADC与DMA的作用,否则ADC只会采集一波数据然后DMA搬运,若需多次采集必须使用次函数,可在数据处理完成后再次使用***//
  69. void ADC_DMA_Trigger(){

  70.         DMA_Cmd(DMA2_Stream0,DISABLE);//若用循环模式就可不用disable再enable,关掉再重启主要起重装NDTR和保护数据的作用
  71. //        DMA_SetCurrDataCounter(DMA2_Stream0,AD3_Value_Length);
  72. //        DMA_ClearITPendingBit( DMA2_Stream0 ,DMA_IT_TCIF0|DMA_IT_DMEIF0|DMA_IT_TEIF0|DMA_IT_HTIF0|DMA_IT_TCIF0 );
  73.         DMA_ClearITPendingBit( DMA2_Stream0 ,DMA_IT_TCIF0);       
  74.         ADC_ClearITPendingBit(ADC3,ADC_IT_OVR);//ADC3->SR = 0;
  75.         DMA_Cmd(DMA2_Stream0,ENABLE);
  76. }
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:02 | 显示全部楼层
  1. #ifndef __ADCDMA_H__
  2. #define __ADCDMA_H__

  3. #include "stm32f4xx.h"  

  4. #define AD3_Value_Length 256 //因为采样点数是256,故需要长度为256的数组

  5. extern uint16_t AD3_Value[AD3_Value_Length];

  6. void AD3_Init();
  7. void ADC_DMA_Trigger();

  8. #endif
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:02 | 显示全部楼层
Timer部分:
由于是定时器触发ADC所以采样率由定时器控制,即:

fs = 84MHz/(arr*psc)
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:02 | 显示全部楼层
  1. #include "Timer.h"

  2. //主频84M

  3. //TIM2 32bit
  4. void TIM2_Init(uint16_t arr, uint16_t psc){
  5.        
  6.         RCC_APB1PeriphClockCmd(RCC_APB1Periph_TIM2,ENABLE);                       //TIM2 时钟使能
  7.        
  8.         TIM_InternalClockConfig(TIM2);       
  9.         TIM_TimeBaseInitTypeDef TIM_TimeBaseInitStructure;                                //时基结构体
  10.         TIM_TimeBaseInitStructure.TIM_Period = arr;                                         //设置自动重装载值
  11.         TIM_TimeBaseInitStructure.TIM_Prescaler =psc;                                         //设置预分频值
  12.         TIM_TimeBaseInitStructure.TIM_ClockDivision = TIM_CKD_DIV1;         //设置时钟分割
  13.         TIM_TimeBaseInitStructure.TIM_CounterMode = TIM_CounterMode_Up; //向上计数模式
  14.        
  15.         TIM_SelectOutputTrigger(TIM2,TIM_TRGOSource_Update);                        //更新溢出向外触发
  16.         TIM_TimeBaseInit(TIM2, &TIM_TimeBaseInitStructure);                         //时基初始化
  17.         TIM_Cmd(TIM2, ENABLE);                                                                                        //定时器使能
  18. }
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:02 | 显示全部楼层
  1. #ifndef __TIMER_H__
  2. #define __TIMER_H__

  3. #include "stm32f4xx.h"                  // Device header

  4. void TIM2_Init(uint16_t arr, uint16_t psc);

  5. #endif
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:02 | 显示全部楼层
main部分:
  1. #include "stm32f4xx.h"                  // Device header
  2. #include "delay.h"
  3. #include "usart.h"
  4. #include "fft_calculate.h"
  5. #include "math.h"
  6. #include "Timer.h"
  7. #include "ADCDMA.h"

  8. /********************************************工程介绍***********************************************************
  9. 此FFT工程,利用Timer定时触发ADC3_IN9触发AD采集并利用DMA2_Stream0搬运至AD3_Value[256] 即采样256个点

  10. /***********************************FFT相关理论计算介绍********************************************************
  11. 设:
  12. 被采目标信号频率:f
  13. 被采目标信号周期:T = 1/f
  14. 采样频率:fs
  15. 采样周期:Ts = 1/fs

  16. 采样点数:NPT
  17. 采样总时间:t = NPT*Ts

  18. 频谱图频率分辨率:f0 = fs/NPT

  19. 采到的被采信号周期数:NT = t/T
  20. 一个被采信号周期的采样点数:fs/f

  21. *************************************************************************************************************/
  22. u16 i;

  23. int main(){
  24.         delay_init(168);
  25.         uart_init(9600);
  26.        
  27.         printf("Start !\r\n");
  28.        
  29.         TIM2_Init(5-1, 525-1);//fs=32KHz        fs=84MHz/(arr*psc)
  30.         AD3_Init();       
  31.         delay_ms(10);


  32.         while(1){
  33.                
  34.                 ADC_DMA_Trigger();                                                                                        //每次都要重新触发,否则只采样一波数据(NPT个AD值)便停止了
  35.                
  36.                 for(i=0; i<NPT; i++){
  37.                        
  38.                         InBufArray[i] = ((signed short)(AD3_Value[i])) << 16;        //将AD值移至实部
  39.                         printf("%d\r\n",AD3_Value[i]);                                                        //打印AD值(片内12位AD:0~4095)
  40.                        
  41.                 }
  42.                
  43.                 cr4_fft_256_stm32(OutBufArray, InBufArray, NPT);                        //FFT运算
  44.                 GetPowerMag();                                                                              //获取信号各次谐波分量的幅值       
  45.                
  46.                 for(i=0; i<NPT/2; i++){                                                                       
  47.                         printf("%d:%d\r\n",i,MagBufArray[i]);                                        //打印幅值
  48.                 }
  49.                 delay_ms(5000);
  50.                
  51. }

  52. }



 楼主| 我爱台妹mmd 发表于 2023-8-27 18:03 | 显示全部楼层
在用函数cr4_fft_256_stm32(OutBufArray, InBufArray, NPT)时,注意将AD值送入InBufArray[]的实部,它的高16位表示复数的实部,低16位表示虚部。

而函数GetPowerMag()是对FFT后输出的复数OutBufArray[]求模值。
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:03 | 显示全部楼层
关于函数ADC_DMA_Trigger()的作用,这里本人才疏学浅也很迷惑,但没有这一句有关中断标志位清除,确实只会运作一次,即只会采一次256个点运算,然后就没有然后了(还望高手指点迷津[抱拳])。我通过debug模式观察寄存器,发现ADC并没有在采样,DMA也没有在搬运数据。但这反而更方便了,让每一次采集都变得可控。需要采集时,使用此函数触发即可。
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:03 | 显示全部楼层
实践
做了这么多工作现在来实操验证一下。此次实验输入的被采的信号均为正弦波。
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:03 | 显示全部楼层
输入1KHz 0V到3V(即幅度:3V,偏移:1.5V)
4972264eb1f7bd66d3.png
 楼主| 我爱台妹mmd 发表于 2023-8-27 18:03 | 显示全部楼层
用示波器先浅浅观察一下: 90464eb1f87df0e4.png
您需要登录后才可以回帖 登录 | 注册

本版积分规则

72

主题

648

帖子

0

粉丝
快速回复 在线客服 返回列表 返回顶部