当前位置：文档库 › 基于STM32芯片的128点FFT

基于STM32芯片的128点FFT

Main函数

* Author : MCD Application Team

* Version : V2.0.1

* Date : 06/13/2008

* Description : Main program body

******************************************************************************* *

* THE PRESENT FIRMWARE WHICH IS FOR GUIDANCE ONL Y AIMS A T PROVIDING CUSTOMERS

* WITH CODING INFORMA TION REGARDING THEIR PRODUCTS IN ORDER FOR THEM TO SA VE TIME.

* AS A RESULT, STMICROELECTRONICS SHALL NOT BE HELD LIABLE FOR ANY DIRECT,

* INDIRECT OR CONSEQUENTIAL DAMAGES WITH RESPECT TO ANY CLAIMS ARISING FROM THE

* CONTENT OF SUCH FIRMW ARE AND/OR THE USE MADE BY CUSTOMERS OF THE CODING

* INFORMATION CONTAINED HEREIN IN CONNECTION WITH THEIR PRODUCTS.

******************************************************************************* /

/* Includes ------------------------------------------------------------------*/

#include "stm32f10x_lib.h"

#include "stdio.h"

#include "math.h"

#include "stm32_dsp.h"

/* Private typedef -----------------------------------------------------------*/

/* Private define ------------------------------------------------------------*/

#define ADC1_DR_Address ((u32)0x4001244C)

#define USART1_DR_Base ((u32)0x40013804)

#define LEN 128

/* Private macro -------------------------------------------------------------*/

/* Private variables ---------------------------------------------------------*/

ErrorStatus HSEStartUpStatus;

/* Private function prototypes -----------------------------------------------*/

void RCC_Configuration(void);

void GPIO_Configuration(void);

void NVIC_Configuration(void);

void TIM_BaseConfiguration(void);

void DMA_Configuration(void);

void USART_Configuration(void);

void ADC_Configuration(void);

void ADC_ALL_Init(void);

void FFT_Configuration(double *pr,double *pi,int n,int k,double*fr,double*fi,int l,int il);

void FFT_solve(void);

void Delay(u32 counter);

/* Private functions ---------------------------------------------------------*/

u16 data_buff[LEN];

u16 px_buff[LEN];

u16 pz_buff[LEN];

double x[128];

double pr[128];

double pi[128];

double fr[128];

double fi[128];

double mo[128];

//int *px = (int*)0x1f00;

//int *pz = (int*)0x1f80;

int xm,zm,i,t = 0;

u8 Flag = 0;

/****************************************************************************** *

* Function Name : main

* Description : Main program

* Input : None

* Output : None

* Return : None

******************************************************************************* /

int main(void)

{

#ifdef DEBUG

debug();

#endif

/* System clocks configuration ---------------------------------------------*/

RCC_Configuration();

/* NVIC configuration ------------------------------------------------------*/

NVIC_Configuration();

/* GPIO configuration ------------------------------------------------------*/

GPIO_Configuration();

/* ADC1 configuration ------------------------------------------------------*/

ADC_Configuration();

/*USART Configuration--------------------------------------------------------*/

USART_Configuration();

ADC_ALL_Init();

while (1)

{

FFT_solve();

}

/****************************************************************************** *

* Function Name : RCC_Configuration

* Description : Configures the different system clocks.

* Input : None

* Output : None

* Return : None

******************************************************************************* /

void RCC_Configuration(void)

{

/* RCC system reset(for debug purpose) */

RCC_DeInit();

/* Enable HSE */

RCC_HSEConfig(RCC_HSE_ON);

/* Wait till HSE is ready */

HSEStartUpStatus = RCC_WaitForHSEStartUp();

if(HSEStartUpStatus == SUCCESS)

{

/* Enable Prefetch Buffer */

FLASH_PrefetchBufferCmd(FLASH_PrefetchBuffer_Enable);

/* Flash 2 wait state */

FLASH_SetLatency(FLASH_Latency_2);

/* HCLK = SYSCLK */

RCC_HCLKConfig(RCC_SYSCLK_Div1);

/* PCLK2 = HCLK */

RCC_PCLK2Config(RCC_HCLK_Div1);

/* PCLK1 = HCLK/2 */

RCC_PCLK1Config(RCC_HCLK_Div2);

/* ADCCLK = PCLK2/4 */

RCC_ADCCLKConfig(RCC_PCLK2_Div4);

/* PLLCLK = 8MHz * 7 = 56 MHz */

RCC_PLLConfig(RCC_PLLSource_HSE_Div1, RCC_PLLMul_7);

/* Enable PLL */

RCC_PLLCmd(ENABLE);

/* Wait till PLL is ready */

while(RCC_GetFlagStatus(RCC_FLAG_PLLRDY) == RESET)

{

}

/* Select PLL as system clock source */

RCC_SYSCLKConfig(RCC_SYSCLKSource_PLLCLK);

/* Wait till PLL is used as system clock source */

while(RCC_GetSYSCLKSource() != 0x08)

{

}

/* Enable peripheral clocks --------------------------------------------------*/

/* Enable DMA1 clock */

RCC_AHBPeriphClockCmd(RCC_AHBPeriph_DMA1, ENABLE);

/* Enable USART2 clock */

RCC_APB1PeriphClockCmd(RCC_APB2Periph_USART1, ENABLE);

/*Enable TIM2 clock*/

RCC_APB1PeriphClockCmd(RCC_APB1Periph_TIM2, ENABLE);

/* Enable ADC1 and GPIOC clock */

RCC_APB2PeriphClockCmd(RCC_APB2Periph_ADC1 | RCC_APB2Periph_GPIOC | RCC_APB2Periph_GPIOA, ENABLE);

}

/****************************************************************************** *

* Function Name : GPIO_Configuration

* Description : Configures the different GPIO ports.

* Input : None

* Output : None

* Return : None

******************************************************************************* /

void GPIO_Configuration(void)

{

GPIO_InitTypeDef GPIO_InitStructure;

/* Configure PC.04 (ADC Channel14) as analog input -------------------------*/

GPIO_InitStructure.GPIO_Pin = GPIO_Pin_4;

GPIO_InitStructure.GPIO_Mode = GPIO_Mode_AIN;

GPIO_Init(GPIOC, &GPIO_InitStructure);

/* Configure PA.03 (ADC Channel14) as analog input -------------------------*/

GPIO_InitStructure.GPIO_Pin = GPIO_Pin_3;

GPIO_InitStructure.GPIO_Mode = GPIO_Mode_AIN;

GPIO_Init(GPIOA, &GPIO_InitStructure);

/* Configure USART1 Rx (PA.10) as input floating */

GPIO_InitStructure.GPIO_Pin = GPIO_Pin_10;

GPIO_InitStructure.GPIO_Mode = GPIO_Mode_IN_FLOATING;

GPIO_Init(GPIOA, &GPIO_InitStructure);

/* Configure USART1 Tx (PA.09) as alternate function push-pull */

GPIO_InitStructure.GPIO_Pin = GPIO_Pin_9;

GPIO_InitStructure.GPIO_Speed = GPIO_Speed_50MHz;

GPIO_InitStructure.GPIO_Mode = GPIO_Mode_AF_PP;

GPIO_Init(GPIOA, &GPIO_InitStructure);

}

/****************************************************************************** *

* Function Name : NVIC_Configuration

* Description : Configures Vector Table base location.

* Input : None

* Output : None

* Return : None

*******************************************************************************

void NVIC_Configuration(void)

{

NVIC_InitTypeDef NVIC_InitStructure;

#ifdef VECT_TAB_RAM

/* Set the Vector Table base location at 0x20000000 */

NVIC_SetVectorTable(NVIC_VectTab_RAM, 0x0);

#else /* VECT_TAB_FLASH */

/* Set the Vector Table base location at 0x08000000 */

NVIC_SetVectorTable(NVIC_VectTab_FLASH, 0x0);

#endif

/* Enable the TIM2 gloabal Interrupt */

NVIC_InitStructure.NVIC_IRQChannel = TIM2_IRQChannel;

NVIC_InitStructure.NVIC_IRQChannelPreemptionPriority = 0;

NVIC_InitStructure.NVIC_IRQChannelSubPriority = 2;

NVIC_InitStructure.NVIC_IRQChannelCmd = ENABLE;

NVIC_Init(&NVIC_InitStructure);

/* Enable the ADC1_EOC Interrupt */

NVIC_InitStructure.NVIC_IRQChannel = ADC1_2_IRQChannel;

NVIC_InitStructure.NVIC_IRQChannelPreemptionPriority = 0;

NVIC_InitStructure.NVIC_IRQChannelSubPriority = 1;

NVIC_InitStructure.NVIC_IRQChannelCmd = ENABLE;

NVIC_Init(&NVIC_InitStructure);

}

/****************************************************************************** *

* Function Name : TIM_BaseConfiguration

* Description : Configures Vector Table base location.

* Input : None

* Output : None

* Return : None

******************************************************************************* /

void TIM_BaseConfiguration(void)

{

TIM_TimeBaseInitTypeDef TIM_TimeBaseStructure;

/* Time base configuration */

TIM_TimeBaseStructure.TIM_Period = 65535;

TIM_TimeBaseStructure.TIM_Prescaler = 0;

TIM_TimeBaseStructure.TIM_ClockDivision = 0;

TIM_TimeBaseStructure.TIM_CounterMode = TIM_CounterMode_Up;

TIM_TimeBaseInit(TIM2, &TIM_TimeBaseStructure);

/*Enable TIM2*/

TIM_Cmd(TIM2,ENABLE);

/* TIM IT enable */

TIM_ITConfig(TIM2, TIM_IT_Update , ENABLE);

}

/****************************************************************************** *

* Function Name : ADC_Configuration

* Description : Configures Vector Table base location.

* Input : None

* Output : None

* Return : None

******************************************************************************* /

void ADC_Configuration(void)

{

ADC_InitTypeDef ADC_InitStructure;

ADC_InitStructure.ADC_Mode = ADC_Mode_Independent; /*工作在独立模式*/

ADC_InitStructure.ADC_ScanConvMode = ENABLE; /*规定了模数转换工作在扫描模式（多通道）还是单次（单通道）模式*/

ADC_InitStructure.ADC_ContinuousConvMode = ENABLE; /*规定了模数转换工作在连续还是单次模式*/

ADC_InitStructure.ADC_ExternalTrigConv = ADC_ExternalTrigConv_None; /*转换由软件而不是外部触发启动*/

ADC_InitStructure.ADC_DataAlign = ADC_DataAlign_Right; /*规定了ADC 数据向左边对齐还是向右边对齐*/

ADC_InitStructure.ADC_NbrOfChannel = 1; /*规定了顺序进行规则转换的ADC通道的数目*/

ADC_Init(ADC1, &ADC_InitStructure);

}

void ADC_ALL_Init(void)

{

/* ADC1 regular channel14 configuration */

ADC_RegularChannelConfig(ADC1, ADC_Channel_14, 1, ADC_SampleTime_7Cycles5);

/*设置指定ADC的规则组通道，设置它们的转化顺序和采样时间*/

/* Enable ADC1 EOC interrupts */

ADC_ITConfig(ADC1, ADC_IT_EOC , ENABLE);

/* Enable ADC1 */

ADC_Cmd(ADC1, ENABLE);

/* Enable ADC1 reset calibaration register */

ADC_ResetCalibration(ADC1); /*重置指定的ADC的校准寄存器*/

/* Check the end of ADC1 reset calibration register */

while(ADC_GetResetCalibrationStatus(ADC1));/*获取ADC重置校准寄存器的状态*/

/* Start ADC1 calibaration */

ADC_StartCalibration(ADC1);

/* Check the end of ADC1 calibration */

while(ADC_GetCalibrationStatus(ADC1));

/* Start ADC1 Software Conversion */

ADC_SoftwareStartConvCmd(ADC1, ENABLE);/*使能指定的ADC的软件转换启动功能*/

}

/****************************************************************************** *

* Function Name : USART_Configuration

* Description : Configures Vector Table base location.

* Input : None

* Output : None

* Return : None

******************************************************************************* /

void USART_Configuration(void)

{

USART_InitTypeDef USART_InitStructure;

USART_https://www.wendangku.net/doc/6615347765.html,ART_BaudRate = 9600;

USART_https://www.wendangku.net/doc/6615347765.html,ART_WordLength = USART_WordLength_8b;

USART_https://www.wendangku.net/doc/6615347765.html,ART_StopBits = USART_StopBits_1;

USART_https://www.wendangku.net/doc/6615347765.html,ART_Parity = USART_Parity_No;

USART_https://www.wendangku.net/doc/6615347765.html,ART_HardwareFlowControl = USART_HardwareFlowControl_None;

USART_https://www.wendangku.net/doc/6615347765.html,ART_Mode = USART_Mode_Rx | USART_Mode_Tx;

/* Configure USART1 */

USART_Init(USART1, &USART_InitStructure);

/*Enable USART1 */

USART_Cmd(USART1, ENABLE);

}

/******************************************************************************

* Function Name : FFT_Configuration

* Description :

* Input : None

* Output : None

* Return : None

******************************************************************************* /

void FFT_Configuration(double *pr,double *pi,int n,int k,double*fr,double*fi,int l,int il)

//int n,k,l,il; double pr[],pi[],fr[],fi[];

/*pr(实部)，pi(虚部)，n(点数)，k(阶数)*/

{

int it,m,is,i,j,nv,l0;

double p,q,s,vr,vi,poddr,poddi;

/*雷德(Rader)算法，实现数据的倒位序排列-----------*/

for (it=0; it<=n-1; it++)

{

m=it;

is=0;

for (i=0; i<=k-1; i++)

{

j=m/2;

is=2*is+(m-2*j);

m=j;

}

fr[it]=pr[is];

fi[it]=pi[is];

}

pr[0]=1.0;

pi[0]=0.0;

p=6.283185306/(1.0*n);

pr[1]=cos(p);

pi[1]=-sin(p);

if (l!=0)

pi[1]=-pi[1];

for (i=2; i<=n-1; i++)

{

p=pr[i-1]*pr[1];

q=pi[i-1]*pi[1];

s=(pr[i-1]+pi[i-1])*(pr[1]+pi[1]);

pr[i]=p-q;

pi[i]=s-p-q;

}

for (it=0; it<=n-2; it=it+2)

{

vr=fr[it];

vi=fi[it];

fr[it]=vr+fr[it+1];

fi[it]=vi+fi[it+1];

fr[it+1]=vr-fr[it+1];

fi[it+1]=vi-fi[it+1];

}

m=n/2;

nv=2;

for (l0=k-2; l0>=0; l0--)

{

m=m/2;

nv=2*nv;

for (it=0; it<=(m-1)*nv; it=it+nv)

for (j=0; j<=(nv/2)-1; j++)

{

p=pr[m*j]*fr[it+j+nv/2];

q=pi[m*j]*fi[it+j+nv/2];

s=pr[m*j]+pi[m*j];

s=s*(fr[it+j+nv/2]+fi[it+j+nv/2]);

poddr=p-q;

poddi=s-p-q;

fr[it+j+nv/2]=fr[it+j]-poddr;

fi[it+j+nv/2]=fi[it+j]-poddi;

fr[it+j]=fr[it+j]+poddr;

fi[it+j]=fi[it+j]+poddi;

}

if (l!=0)

for (i=0; i<=n-1; i++)

{ fr[i]=fr[i]/(1.0*n);

fi[i]=fi[i]/(1.0*n);

}

if (il!=0)

for (i=0; i<=n-1; i++)

{

pr[i]=sqrt(fr[i]*fr[i]+fi[i]*fi[i]);

if (fabs(fr[i])<0.000001*fabs(fi[i]))

{

if ((fi[i]*fr[i])>0) pi[i]=90.0;

else

pi[i]=-90.0;

}

else

pi[i]=atan(fi[i]/fr[i])*360.0/6.283185306;

}

/****************************************************************************** *

* Function Name : FFT_solve

* Description :

* Input : None

* Output : None

* Return : None

******************************************************************************* /

void FFT_solve(void)

{

if(Flag)

{

/*寄存器接收满后Flag = 1,关闭ADC中断,进行FFT运算--------*/

ADC_Cmd(ADC1, DISABLE);

//px = (int*)0x1f00;

/*转换存储位置-------------*/

for(i=0;i<128;i++)

{

//*px=data_buff[i];

//px++;

px_buff[i] = data_buff[i];

}

//px = (int*)0x1f00;

/*将ADC采来的数据转换为实际的幅度值--------------*/

for (i=0; i<=128; i++)

{

//xm=*px;

xm = px_buff[i];

x[i]=xm/32768.0;

pr[i]=x[i];/*pr数组存放函数点的实部-------------*/

pi[i]=0; /*pi数组存放函数点的虚部-------------*/

//px++;

}

/*进行FFT运算----------------------*/

FFT_Configuration(pr,pi,128,7,fr,fi,0,1);

//pz = (int*)0x1f80;

/*将最终计算出的值放入pz 中----------------------*/

for (i=0;i<=128;i++)

{

mo[i] = sqrt(fr[i]*fr[i]+fi[i]*fi[i]);

zm = (int)(mo[i]*1000.0);

//*pz = zm;

//pz++;

pz_buff[i] = zm;

}

Flag = 0;

}

/****************************************************************************** *

* Function Name : Delay

* Description :

* Input : None

* Output : None

* Return : None

******************************************************************************* /

void Delay(u32 counter)

{

for(;counter > 0;counter--);

}

#ifdef DEBUG

/****************************************************************************** *

* Function Name : assert_failed

* Description : Reports the name of the source file and the source line number

* where the assert_param error has occurred.

* Input : - file: pointer to the source file name

* - line: assert_param error line source number

* Output : None

* Return : None

******************************************************************************* /

void assert_failed(u8* file, u32 line)

{

/* User can add his own implementation to report the file name and line number, ex: printf("Wrong parameters value: file %s on line %d\r\n", file, line) */

/* Infinite loop */

while (1)

{

}

#endif

stm32f10x_it.c

* Author : MCD Application Team

* Version : V2.0.1

* Date : 06/13/2008

* Description : Main Interrupt Service Routines.

* This file provides template for all exceptions handler

* and peripherals interrupt service routine.