DSP28335做简单的频谱分析fft算法

DSP28335AD采样后做fft变换，得出频谱 根据在RAM中调试的需要，这个项目配置成"boot to SARAM".2833x引导模式 表如下显示. 常用的还有"boot to Flash"模式，当程序在RAM调试完善后就 可以将代码烧进Flash中并使用"boot to Flash"引导模式.

#include “DSP2833x_Device.h“     // DSP2833x Headerfile Include File
#include “DSP2833x_Examples.h“   // DSP2833x Examples Include File
#include  
// Determine when the shift to right justify the data takes place
// Only one of these should be defined as 1.
// The other two should be defined as 0.
#define POST_SHIFT   0  // Shift results after the entire sample table is full
#define INLINE_SHIFT 1  // Shift results as the data is taken from the results regsiter
#define NO_SHIFT     0  // Do not shift the results

#define pi 3.141593 // float小数点后6位
#define NL 256  // NL为合成信号点数，与N的值必须是一致的，即NL=N
    int N=NL;   //FFT点数

// ADC start parameters
#if （CPU_FRQ_150MHZ）     // Default - 150 MHz SYSCLKOUT
  #define ADC_MODCLK 0x3 // HSPCLK = SYSCLKOUT/2*ADC_MODCLK2 = 150/（2*3）   = 25.0 MHz
#endif
#if （CPU_FRQ_100MHZ）
  #define ADC_MODCLK 0x2 // HSPCLK = SYSCLKOUT/2*ADC_MODCLK2 = 100/（2*2）   = 25.0 MHz
#endif
#define ADC_CKPS   0x0   // ADC module clock = HSPCLK/1      = 25.5MHz/（1）   = 25.0 MHz
#define ADC_SHCLK  0x1   // S/H width in ADC module periods                  = 2 ADC cycle
#define AVG        1000  // Average sample limit
#define ZOFFSET    0x00  // Average Zero offset
#define BUF_SIZE   512  // Sample buffer size
/**************************************函数声明************************************************/
void scia_echoback_init（void）;  //声明SCI-A工作方式和参数配置
void scia_xmit（int a）;          //声明发送字节的函数
void scia_msg（char *msg）;       //声明发送数组的函数
/**********************************************************************************************/
// Global variable for this example
Uint16 SampleTable[BUF_SIZE];
struct Complex		// 定义复数结构体
{
   float realimag;
};
struct Complex Wn;//定义旋转因子
struct Complex Vn;//每一级第一个旋转因子虚部为0，实部为1
struct Complex T;//存放旋转因子与X（k+B）的乘积

//float Realin[NL]={0};// 采样输入的实数
float output[NL]={0};// 输出的FFT幅值（复数的模）
struct Complex Sample[NL];// 采样输入的实数转化为复数

struct Complex MUL（struct Complex astruct Complex b）//定义复乘
{
   struct Complex c;
   c.real=a.real*b.real-a.imag*b.imag;
   c.imag=a.real*b.imag+a.imag*b.real;
   return（c）;
}

void MYFFT（struct Complex *xinint N）//输入为复数指针*xin，做N点FFT
{
   int L=0; // 级间运算层
   int J=0; // 级内运算层
   int K=0KB=0; // 蝶形运算层
   int M=1Nn=0;// N=2^M
   float B=0; // 蝶形运算两输入数据间隔
   /* 以下是为倒序新建的局部变量*/
   int LH=0J2=0N1=0IK2=0;
   struct Complex T;
   /*以下是倒序*/
   LH=N/2; // LH=N/2
   J2=LH;
   N1=N-2;
   for（I=1;I<=N1;I++）
    {
     if（I	 {
       T=xin[I];
       xin[I]=xin[J2];
       xin[J2]=T;
	 }
	 K2=LH;
	 while（J2>=K2）
	  {
        J2-=K2;
        K2=K2/2;// K2=K2/2
	  }
	  J2+=K2;
    }
   /* 以下为计算出M */
   Nn=N;
   while（Nn!=2）// 计算出N的以2为底数的幂M
   {
     M++;
	 Nn=Nn/2;
   }

   /* 蝶形运算 */
   for（L=1;L<=M;L++）  // 级间
	{
	  B=pow（2（L-1））;
      Vn.real=1;
	  Vn.imag=0;
      Wn.real=cos（pi/B）;
      Wn.imag=-sin（pi/B）;
	  for（J=0;J	   {
		 for（K=J;K


 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-10-19 15:31  lab25-AD\
     文件         329  2018-10-19 15:31  lab25-AD\.ccsproject
     文件       17159  2015-04-20 14:12  lab25-AD\.cdtbuild
     文件       17159  2018-10-19 15:31  lab25-AD\.cdtbuild_initial
     文件         553  2015-04-20 14:12  lab25-AD\.cdtproject
     文件       31099  2018-10-19 15:31  lab25-AD\.cproject
     文件         838  2018-10-19 15:31  lab25-AD\.project
     文件         588  2018-10-19 15:31  lab25-AD\.project_initial
     目录           0  2018-10-19 15:31  lab25-AD\.settings\
     文件          62  2018-10-19 15:31  lab25-AD\.settings\org.eclipse.cdt.codan.core.prefs
     文件         123  2018-10-19 15:31  lab25-AD\.settings\org.eclipse.cdt.debug.core.prefs
     文件         606  2018-10-19 15:31  lab25-AD\.settings\org.eclipse.cdt.managedbuilder.core.prefs
     文件         301  2018-10-19 15:31  lab25-AD\.settings\org.eclipse.core.resources.prefs
     目录           0  2018-10-19 15:31  lab25-AD\AD.CS_\
     文件        3072  2015-04-20 14:12  lab25-AD\AD.CS_\FILE.CDX
     文件        1492  2015-04-20 14:12  lab25-AD\AD.CS_\FILE.DBF
     文件        2256  2015-04-20 14:12  lab25-AD\AD.CS_\FILE.FPT
     文件      323072  2015-04-20 14:12  lab25-AD\AD.CS_\SYMBOL.CDX
     文件      238414  2015-04-20 14:12  lab25-AD\AD.CS_\SYMBOL.DBF
     文件      412331  2015-04-20 14:12  lab25-AD\AD.CS_\SYMBOL.FPT
     文件        9265  2015-04-20 14:12  lab25-AD\AD.paf2
     文件        6688  2015-04-20 14:12  lab25-AD\AD.sbl
     目录           0  2018-10-19 15:31  lab25-AD\CMD\
     文件        3941  2015-04-20 14:12  lab25-AD\CMD\28335_RAM_lnk.cmd
     文件        8438  2015-04-20 14:12  lab25-AD\CMD\DSP2833x_Headers_nonBIOS.cmd
     目录           0  2018-11-23 15:06  lab25-AD\Debug\
     文件       30284  2018-11-23 15:06  lab25-AD\Debug\AD.map
     文件       27525  2018-11-23 15:06  lab25-AD\Debug\AD.obj
     文件      144727  2018-11-23 15:06  lab25-AD\Debug\AD.out
     文件      155372  2018-11-23 15:06  lab25-AD\Debug\AD_linkInfo.xml
     目录           0  2018-10-19 15:31  lab25-AD\Debug\CMD\
............此处省略84个文件信息