STM32F7 adc+tim+dma 进行FFT计算的问题

啾啾啾

大家好！想请教一下大家。我在使用DSP库进行FFT处理风速传感器的波形信号。代码如下：
ADC+DMA:

1. ADC_ChannelConfTypeDef hadc2_channel;
   
2. void MY_ADC_Init(void)
   
3. {
   
4.     __HAL_RCC_DMA2_CLK_ENABLE();
   
5. 
   
6.     hdma_adc2.Instance = DMA2_Stream2;
   
7.     hdma_adc2.Init.Channel = DMA_CHANNEL_1;
   
8.     hdma_adc2.Init.Direction = DMA_PERIPH_TO_MEMORY;
   
9.     hdma_adc2.Init.PeriphInc = DMA_PINC_DISABLE;               //???????
   
10.     hdma_adc2.Init.MemInc = DMA_MINC_ENABLE;                   //???????
    
11.     hdma_adc2.Init.PeriphDataAlignment = DMA_PDATAALIGN_HALFWORD;  //??????:16?
    
12.     hdma_adc2.Init.MemDataAlignment = DMA_MDATAALIGN_HALFWORD;     //???????:16?
    
13.     hdma_adc2.Init.Mode = DMA_CIRCULAR;                          //???????
    
14.     hdma_adc2.Init.Priority = DMA_PRIORITY_LOW;             //?????
    
15.     hdma_adc2.Init.FIFOMode = DMA_FIFOMODE_DISABLE;            /* ??FIFO*/
    
16. 
    
17.     HAL_DMA_Init(&hdma_adc2);
    
18. 
    
19.     __HAL_LINKDMA(&hadc2, DMA_Handle, hdma_adc2);                //?DMA?ADC????
    
20. 
    
21. 
    
22.     hadc2.Instance = ADC2;
    
23.     hadc2.Init.ClockPrescaler = ADC_CLOCK_SYNC_PCLK_DIV4; //4??,ADCCLK=PCLK2/4=108/4=27MHZ
    
24.     hadc2.Init.Resolution = ADC_RESOLUTION_12B;           //12???
    
25.     hadc2.Init.ScanConvMode = DISABLE;                    //?????
    
26.     hadc2.Init.ContinuousConvMode = DISABLE;              //??????
    
27.     hadc2.Init.DiscontinuousConvMode = DISABLE;           //?????????
    
28.     hadc2.Init.ExternalTrigConvEdge = ADC_EXTERNALTRIGCONVEDGE_RISING; //??????
    
29.     hadc2.Init.ExternalTrigConv = ADC_EXTERNALTRIGCONV_T2_CC2;     //????
    
30.     hadc2.Init.DataAlign = ADC_DATAALIGN_RIGHT;           //???
    
31.     hadc2.Init.NbrOfConversion = 1;                       //1????????? ??????????1
    
32.     hadc2.Init.DMAContinuousRequests = ENABLE;           //??DMA??
    
33.     hadc2.Init.EOCSelection = ADC_EOC_SINGLE_CONV;
    
34.     HAL_ADC_Init(&hadc2);
    
35. 
    
36.     hadc2_channel.Channel = ADC_CHANNEL_3;                                 //??
    
37.     hadc2_channel.Rank = 1;                                     //??1
    
38.     hadc2_channel.SamplingTime = ADC_SAMPLETIME_3CYCLES;      //????
    
39.     hadc2_channel.Offset = 0;
    
40.     HAL_ADC_ConfigChannel(&hadc2, &hadc2_channel);       //????
    
41.                
    
42. }

复制代码

TIM:

1. TIM_HandleTypeDef htim2;
   
2. TIM_OC_InitTypeDef TIM2_CH2Handler;
   
3. void MY_TIM2_Init(void)
   
4. {
   
5.     __HAL_RCC_TIM2_CLK_ENABLE();
   
6.        
   
7. 
   
8. 
   
9.     htim2.Instance = TIM2;
   
10.     htim2.Init.Prescaler = 108-1;
    
11.     htim2.Init.CounterMode = TIM_COUNTERMODE_UP;
    
12.     htim2.Init.Period = 100;
    
13.     htim2.Init.ClockDivision = TIM_CLOCKDIVISION_DIV1;
    
14. 
    
15.     HAL_TIM_PWM_Init(&htim2);
    
16. 
    
17. 
    
18.     TIM2_CH2Handler.OCMode = TIM_OCMODE_PWM1;
    
19.     TIM2_CH2Handler.Pulse = 50;
    
20.     TIM2_CH2Handler.OCPolarity = TIM_OCPOLARITY_LOW;
    
21.     TIM2_CH2Handler.OCFastMode = TIM_OCFAST_DISABLE;
    
22.     HAL_TIM_PWM_ConfigChannel(&htim2, &TIM2_CH2Handler, TIM_CHANNEL_2) ;
    
23. 
    
24.     TIM_CCxChannelCmd(htim2.Instance, TIM_CHANNEL_2, TIM_CCx_ENABLE);
    
25.     __HAL_TIM_ENABLE(&htim2);
    
26.                 //__HAL_TIM_DISABLE(&htim2);
    
27. }

复制代码

FFT:

1. uint16_t FFT_ADC(uint16_t* buffer)
   
2. {
   
3.         uint16_t index;
   
4.         for(int j=0;j<COUNT;j++)
   
5.         {
   
6.                 FFT_IN_Buffer[2*j] = buffer[j];            //实部//
   
7.                 FFT_IN_Buffer[2*j+1] = 0;                      //虚部为0//
   
8.         }
   
9.         arm_cfft_f32(&arm_cfft_sR_f32_len1024,FFT_IN_Buffer,0,1);           //1024点FFT//
   
10.         arm_cmplx_mag_f32(FFT_IN_Buffer,FFT_OUT_Buffer,COUNT);              //复数模平方，算幅值//
    
11. //        printf("FFT result is:\r\n");
    
12. //        for (int i = 0; i < COUNT; i++)
    
13. //  {
    
14. //        if (i % 8 == 0) printf("\r\n");
    
15. //        printf("%.1f  ", FFT_OUT_Buffer[i]);
    
16. //        }
    
17.         index = Max_Index(FFT_OUT_Buffer,COUNT);
    
18. //        printf("\r\n****************************************************************\r\n");
    
19.         return index;
    
20. }
    
21. 
    
22. //找出频率最大值对应的下标//
    
23. uint16_t Max_Index(float* buffer,int N)
    
24. {
    
25.         float Max = buffer[1];
    
26.         uint16_t index = 1;
    
27.         for(int i=2;i<N;i++)
    
28.         {
    
29.                 if(Max<buffer[i])
    
30.                 {
    
31.                         Max = buffer[i];
    
32.                         index = i;
    
33.                 }
    
34.         }
    
35.         return index;
    
36. }

复制代码

主程序中

1. int main()
   
2. {
   
3.            MY_ADC_Init();
   
4.         MY_TIM2_Init();
   
5.         HAL_ADC_Start_DMA(&hadc2, (uint32_t *)ADC2_Buffer, sizeof(ADC2_Buffer));
   
6.                 uint16_t index;
   
7.         int Frequency =0;
   
8.          while(1)
   
9.          {
   
10.                                               index = FFT_ADC(ADC2_Buffer);
    
11.                                 Frequency = (Fs/COUNT)*index;
    
12.                                 printf("%d\r\n", Frequency);
    
13.           }
    
14. }

复制代码

我通过调试 观察发现FFT_IN的数组，奇数下标（对应虚部）并没有被赋值为0（如图1）
FFT计算出来的频率也有些不对劲（图2）。请问一下是哪个环节出现了纰漏。谢谢！

bin133

帮顶  

啾啾啾

bin133 发表于 2021-10-21 09:44
帮顶

谢谢大佬

weiwei123321

DMA装满后，停掉DMA和定时器，将数据进行FFT后打开DMAHE 定时器

liaoliaoliao

请问 最后搞定了没 如果要测相位 怎样测呢