標(biāo)題: STM32學(xué)習(xí)日志--使用dsp庫的FFT函數(shù) [打印本頁]
作者: a353114532    時(shí)間: 2018-5-21 14:02
標(biāo)題: STM32學(xué)習(xí)日志--使用dsp庫的FFT函數(shù)
STM32  dsp  FFT FFT函數(shù)單片機(jī)源程序如下:
  1. /**
  2.   ******************************************************************************
  3.   * @file FFT_Demo/src/main.c
  4.   * @author  MCD Application Team
  5.   * @version  V2.0.0
  6.   * @brief  Main program body
  7.   ******************************************************************************
  8.   * @copy
  9.   *
  10.   * THE PRESENT FIRMWARE WHICH IS FOR GUIDANCE ONLY AIMS AT PROVIDING CUSTOMERS
  11.   * WITH CODING INFORMATION REGARDING THEIR PRODUCTS IN ORDER FOR THEM TO SAVE
  12.   * TIME. AS A RESULT, STMICROELECTRONICS SHALL NOT BE HELD LIABLE FOR ANY
  13.   * DIRECT, INDIRECT OR CONSEQUENTIAL DAMAGES WITH RESPECT TO ANY CLAIMS ARISING
  14.   * FROM THE CONTENT OF SUCH FIRMWARE AND/OR THE USE MADE BY CUSTOMERS OF THE
  15.   * CODING INFORMATION CONTAINED HEREIN IN CONNECTION WITH THEIR PRODUCTS.
  16.   *
  17.   * <h2><center>© COPYRIGHT 2009 STMicroelectronics</center></h2>
  18.   */
  19. /* Includes ------------------------------------------------------------------*/
  20. #include "stm32f10x.h"
  21. #include <math.h>
  22. #include "stm3210b_lcd.h"
  23. #include "stm32_dsp.h"
  24. #include "table_fft.h"   
  25. //------------------------------------------------------------------------------

  27. #define FreqSample           3200                // 50Hz*64=3200
  28. #define ScalingFactor   20000                // 幅值系數(shù),最大為32767
  29. #define PERCENT_DC                  0.06                // 直流分量對應(yīng)的百分比,對應(yīng)直流分量lBUFMAG[0]=1200
  30. #define PI2  6.28318530717959            // 2*3.14
  31. #define NPT 64                                    /* 采樣點(diǎn)數(shù)，必須為64NPT = No of FFT point*/
  32. //------------------------------------------------------------------------------
  33. extern volatile uint32_t TimingDelay ;
  34. long lBUFIN[NPT];         /* Complex input vector */
  35. long lBUFOUT[NPT];        /* Complex output vector */
  36. long lBUFMAG[NPT + NPT/2];/* Magnitude vector */

  38. /* Private function prototypes -----------------------------------------------*/
  39. void MyDualSweep(uint32_t freqinc1,uint32_t freqinc2);
  40. void MygSin(long nfill, long Fs, long Freq1, long Freq2, long Ampli);
  41. void powerMag(long nfill, char* strPara);
  42. void DSPDemoInit(void);
  43. //------------------------------------------------------------------------------
  44. //------------------------------------------------------------------------------
  45. //------------------------------------------------------------------------------
  46. int main(void)
  47. {
  48.   DSPDemoInit();

  50.   while (1)
  51.   {
  52.         MyDualSweep(FreqSample/NPT,6*FreqSample/NPT);     // 測試fft
  53.   }
  54. }
  55. //-----------------------------------------------------------------------------

  57. /**
  58.   * @brief  Produces a combination of two sinewaves as input signal
  59.   * @param freq1: frequency increment for 1st sweep，頻率步進(jìn)值1
  60.   *   Freq2: frequency increment for 2nd sweep，頻率步進(jìn)值2
  61.   * @retval : None
  62.   */
  63. void MyDualSweep(uint32_t freqinc1,uint32_t freqinc2)
  64. {
  65.   uint32_t freq;
  66.   
  67.   //for (freq=FreqSample/NPT; freq <FreqSample/2; freq+=freqinc1)  // 40--4000，頻率步進(jìn)值為freqinc1
  68.   for (freq=FreqSample/NPT; freq <FreqSample*6/NPT; freq+=freqinc1)  // 40--4000，頻率步進(jìn)值為freqinc1
  69.   {
  70.             MygSin(NPT, FreqSample, freq, 0, ScalingFactor);// 單獨(dú)的頻率，頻率步進(jìn)值2=0；
  71.                                                                                            
  72. #if (NPT==64)
  73.             cr4_fft_64_stm32(lBUFOUT, lBUFIN, NPT);
  74. #elif (NPT==256)
  75.             cr4_fft_256_stm32(lBUFOUT, lBUFIN, NPT);
  76. #else
  77.             cr4_fft_1024_stm32(lBUFOUT, lBUFIN, NPT);
  78. #endif

  80.             powerMag(NPT,"2SIDED");
  81.   }

  83. //-----------------------------------------------------------------------------

  85.         for (freq=FreqSample/NPT; freq <FreqSample/2; freq+=freqinc2)
  86.           {
  87.             MygSin(NPT, FreqSample, freq, FreqSample*5/NPT, ScalingFactor);

  89. #if (NPT==64)
  90.             cr4_fft_64_stm32(lBUFOUT, lBUFIN, NPT);
  91. #elif (NPT==256)
  92.             cr4_fft_256_stm32(lBUFOUT, lBUFIN, NPT);
  93. #else
  94.             cr4_fft_1024_stm32(lBUFOUT, lBUFIN, NPT);
  95. #endif
  96.    
  97.             powerMag(NPT,"2SIDED");
  98.   }
  99. }
  100. //------------------------------------------------------------------------------
  101. /**
  102.   * @brief  Produces a combination of two sinewaves as input signal
  103.   * @param ill: length of the array holding input signal
  104.   *   Fs: sampling frequency
  105.   *   Freq1: frequency of the 1st sinewave
  106.   *   Freq2: frequency of the 2nd sinewave
  107.   *   Ampli: scaling factor
  108.   * @retval : None
  109.   */
  110. void MygSin(long nfill, long Fs, long Freq1, long Freq2, long Ampli)
  111. {
  112.   uint32_t i;
  113.   float fFs, fFreq1, fFreq2, fAmpli;
  114.   float fZ,fY;

  116.   fFs = (float) Fs;
  117.   fFreq1 = (float) Freq1;
  118.   fFreq2 = (float) Freq2;
  119.   fAmpli = (float) Ampli;

  121.   for (i=0; i < nfill; i++)
  122.   {
  123.     fY =PERCENT_DC/2+sin(PI2 * i * (fFreq1/fFs)) + 0.2*sin(PI2 * i * (fFreq2/fFs));
  124.     fZ = fAmpli * fY;
  125.        
  126. // lBUFIN數(shù)組中，每個(gè)單元數(shù)據(jù)高字(高16位)中存儲(chǔ)采樣數(shù)據(jù)的實(shí)部，低字(低16位)存儲(chǔ)采樣數(shù)據(jù)的虛部(總是為0)
  127.     lBUFIN[i]= ((short)fZ) << 16 ;  /* sine_cosine  (cos=0x0) */
  128.                                                                         // 高16bit為實(shí)部，低16bit為虛部，沒有的話，就用0代替
  129.   }
  130. }
  131. //------------------------------------------------------------------------------
  132. /**
  133.   * @brief  Removes the aliased part of the spectrum (not tested)
  134.   * @param ill: length of the array holding power mag
  135.   * @retval : None
  136.   */
  137. void onesided(long nfill)
  138. {
  139.   uint32_t i;
  140.   
  141.   lBUFMAG[0] = lBUFMAG[0];
  142.   lBUFMAG[nfill/2] = lBUFMAG[nfill/2];
  143.   for (i=1; i < nfill/2; i++)
  144.   {
  145.     lBUFMAG[i] = lBUFMAG[i] + lBUFMAG[nfill-i];
  146.     lBUFMAG[nfill-i] = 0x0;
  147.   }
  148. }
  149. //------------------------------------------------------------------------------
  150. /**
  151.   * @brief  Compute power magnitude of the FFT transform  進(jìn)行FFT變換，并計(jì)算各次諧波幅值
  152.   * @param ill: length of the array holding power mag
  153.   *   : strPara: if set to "1SIDED", removes aliases part of spectrum (not tested)
  154.   * @retval : None
  155.   */
  156. void powerMag(long nfill, char* strPara)
  157. {
  158.   static int32_t lX,lY;
  159.   uint32_t i;

  161.   for (i=0; i < nfill/2; i++)
  162.   {
  163.     lX= (lBUFOUT[i]<<16)>>16; /* 取低16bit，sine_cosine --> cos */
  164.     lY= (lBUFOUT[i] >> 16);   /* 取高16bit，sine_cosine --> sin */   
  165.     {
  166.       float X=  NPT*((float)lX)/32768;
  167.       float Y = NPT*((float)lY)/32768;
  168.       float Mag = sqrt(X*X+ Y*Y)/nfill;
  169.       lBUFMAG[i] = (uint32_t)(Mag*65536);
  170.     }   
  171.   }
  172. }
  173. //------------------------------------------------------------------------------

  175. /**
  176.   * @brief  Inserts a delay time.
  177.   * @param ount: specifies the delay time length (time base 10 ms).
  178.   * @retval : None
  179.   */
  180. void Delay(uint32_t nCount)
  181. {
  182.   /* Configure the systick */
  183.   __enable_irq();
  184.          
  185.      /* Setup SysTick Timer for 10 msec interrupts */
  186.   if (SysTick_Config(SystemFrequency /100)) /* SystemFrequency is defined in
  187.    搒ystem_stm32f10x.h?and equal to HCLK frequency */
  188.    {
  189.     /* Capture error */
  190.      while (1);
  191.     }     
  192.   /* Enable the SysTick Counter */
  193.   TimingDelay = nCount;
  194.    while (TimingDelay )  {}            
  195.    TimingDelay=0;
  196.    __disable_irq();

  198. }

  200. /**
  201.   * @brief  Initializes the DSP lib demo (clock, peripherals and LCD).
  202.   * @param  None
  203.   * @retval : None
  204.   */
  205. void DSPDemoInit(void)
  206. {
  207.   GPIO_InitTypeDef GPIO_InitStructure;
  208.   ErrorStatus HSEStartUpStatus = ERROR;

  210.   /* SYSCLK, HCLK, PCLK2 and PCLK1 configuration -----------------------------*/
  211.   /* RCC system reset(for debug purpose) */
  212.   RCC_DeInit();

  214.   /* Enable HSE */
  215.   RCC_HSEConfig(RCC_HSE_ON);

  217.   /* Wait till HSE is ready */
  218.   HSEStartUpStatus = RCC_WaitForHSEStartUp();

  220.   if(HSEStartUpStatus == SUCCESS)
  221.   {
  222.     /* Enable Prefetch Buffer */
  223.     FLASH_PrefetchBufferCmd(FLASH_PrefetchBuffer_Enable);

  225.     /* Flash 2 wait state */
  226.     FLASH_SetLatency(FLASH_Latency_2);

  228.     /* HCLK = SYSCLK */
  229.     RCC_HCLKConfig(RCC_SYSCLK_Div1);

  231.     /* PCLK2 = HCLK */
  232.     RCC_PCLK2Config(RCC_HCLK_Div1);

  234.     /* PCLK1 = HCLK/2 */
  235.     RCC_PCLK1Config(RCC_HCLK_Div2);

  237.     /* PLLCLK = 8MHz * 9 = 72 MHz */
  238.     RCC_PLLConfig(RCC_PLLSource_HSE_Div1, RCC_PLLMul_9);

  240.     /* Enable PLL */
  241.     RCC_PLLCmd(ENABLE);

  243.     /* Wait till PLL is ready */
  244.     while(RCC_GetFlagStatus(RCC_FLAG_PLLRDY) == RESET)
  245.     {
  246.     }

  248.     /* Select PLL as system clock source */
  249.     RCC_SYSCLKConfig(RCC_SYSCLKSource_PLLCLK);

  251.     /* Wait till PLL is used as system clock source */
  252.     while(RCC_GetSYSCLKSource() != 0x08)
  253.     {
  254.     }
  255.   }

  257.   /* Enable GPIOA, GPIOB, GPIOC, GPIOD, GPIOE and AFIO clocks */
  258.   RCC_APB2PeriphClockCmd(RCC_APB2Periph_GPIOA | RCC_APB2Periph_GPIOB |RCC_APB2Periph_GPIOC
  259.          | RCC_APB2Periph_GPIOD | RCC_APB2Periph_GPIOE | RCC_APB2Periph_AFIO, ENABLE);

  261.   /* TIM1 Periph clock enable */
  262.   RCC_APB2PeriphClockCmd(RCC_APB2Periph_TIM1, ENABLE);
  263.   RCC_APB2PeriphClockCmd(RCC_APB2Periph_GPIOC, ENABLE);


  266.   /* SPI2 Periph clock enable */
  267.   RCC_APB1PeriphClockCmd(RCC_APB1Periph_SPI2, ENABLE);

  269.   /* TIM2  and TIM4 clocks enable */
  270.   RCC_APB1PeriphClockCmd(RCC_APB1Periph_TIM2 | RCC_APB1Periph_TIM4, ENABLE);

  272.   /*------------------- Resources Initialization -----------------------------*/
  273.   /* GPIO Configuration */
  274.   /* Configure PC.06 and PC.07 as Output push-pull */
  275.   GPIO_InitStructure.GPIO_Pin = GPIO_Pin_6 | GPIO_Pin_7;
  276.   GPIO_InitStructure.GPIO_Speed = GPIO_Speed_50MHz;
  277.   GPIO_InitStructure.GPIO_Mode = GPIO_Mode_Out_PP;
  278.   GPIO_Init(GPIOC, &GPIO_InitStructure);



  282.   /*------------------- Drivers Initialization -------------------------------*/
  283.   /* Initialize the LCD */
  284.   STM3210B_LCD_Init();

  286.   /* Clear the LCD */
  287.   LCD_Clear(White);

  289. ……………………

  291. …………限于本文篇幅 余下代碼請從51黑下載附件…………
復(fù)制代碼

所有資料51hei提供下載:
STM32學(xué)習(xí)日志(23)----使用dsp庫的FFT函數(shù).zip (446.86 KB, 下載次數(shù): 73)



作者: DONGpig    時(shí)間: 2020-4-24 18:21
您好，我使用STM32F103的官方DSP庫，對一個(gè)2HZ的標(biāo)準(zhǔn)正弦信號(hào)進(jìn)行FFT變換，出來的幅值信號(hào)波動(dòng)大，有30-50個(gè)數(shù)，您那邊有這種問題嗎？？？想問問如何解決的
作者: DONGpig    時(shí)間: 2020-4-24 18:23
您好，我使用STM32F103的官方DSP庫對一個(gè)2HZ的標(biāo)準(zhǔn)正弦信號(hào)進(jìn)行FFT變換，但是變換后的幅值波動(dòng)大，在30-50個(gè)數(shù)之間，不知道您有沒有遇到這種情況，怎么解決的？？？



歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1