设为首页收藏本站
开启辅助访问 切换到宽版

 找回密码
 立即注册

QQ登录

只需一步，快速开始

热门： 51单片机 | 24小时必答区 | 单片机教程 | 单片机DIY制作 | STM32 | Cortex M3 | 模数电子 | 电子DIY制作 | 音响/功放 | 拆机乐园 | Arduino | 嵌入式OS | 程序设计

搜索
»论坛 › 嵌入式/单片机论坛 › STM32/8 › STM32学习日志--使用dsp库的FFT函数
返回列表 发新帖
查看: 4996|回复: 2
打印 上一主题 下一主题
收起左侧

STM32学习日志--使用dsp库的FFT函数

[复制链接]
跳转到指定楼层
楼主
ID:224401 发表于 2018-5-21 14:02 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
STM32  dsp  FFT FFT函数单片机源程序如下:
  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                // 幅值系数,最大为32767
  29. #define PERCENT_DC                  0.06                // 直流分量对应的百分比,对应直流分量lBUFMAG[0]=1200
  30. #define PI2  6.28318530717959            // 2*3.14
  31. #define NPT 64                                    /* 采样点数，必须为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，频率步进值1
  60.   *   Freq2: frequency increment for 2nd sweep，频率步进值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，频率步进值为freqinc1
  68.   for (freq=FreqSample/NPT; freq <FreqSample*6/NPT; freq+=freqinc1)  // 40--4000，频率步进值为freqinc1
  69.   {
  70.             MygSin(NPT, FreqSample, freq, 0, ScalingFactor);// 单独的频率，频率步进值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数组中，每个单元数据高字(高16位)中存储采样数据的实部，低字(低16位)存储采样数据的虚部(总是为0)
  127.     lBUFIN[i]= ((short)fZ) << 16 ;  /* sine_cosine  (cos=0x0) */
  128.                                                                         // 高16bit为实部，低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  进行FFT变换，并计算各次谐波幅值
  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黑下载附件…………
复制代码

所有资料51hei提供下载:
STM32学习日志(23)----使用dsp库的FFT函数.zip (446.86 KB, 下载次数: 73)



STM32, dsp, FFT函数
分享到:  QQ好友和群QQ好友和群 QQ空间QQ空间 腾讯微博腾讯微博 腾讯朋友腾讯朋友
收藏收藏2 分享淘帖 顶顶 踩踩

相关帖子
回复

使用道具 举报

沙发
ID:730683 发表于 2020-4-24 18:21 | 只看该作者
您好，我使用STM32F103的官方DSP库，对一个2HZ的标准正弦信号进行FFT变换，出来的幅值信号波动大，有30-50个数，您那边有这种问题吗？？？想问问如何解决的
回复

使用道具 举报

板凳
ID:730683 发表于 2020-4-24 18:23 | 只看该作者
您好，我使用STM32F103的官方DSP库对一个2HZ的标准正弦信号进行FFT变换，但是变换后的幅值波动大，在30-50个数之间，不知道您有没有遇到这种情况，怎么解决的？？？
回复

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

手机版|小黑屋|51黑电子论坛 |51黑电子论坛6群 QQ 管理员QQ:125739409;技术交流QQ群281945664

Powered by 单片机教程网

快速回复 返回顶部 返回列表