资源简介
CCS 3.3下编译通过,FFT算法的设计过程与验证程序,采用C语言实现。
代码片段和文件信息
#include “math.h“
#define PI 3.1415926
#define SAMPLENUMBER 128
void InitForFFT();
void MakeWave();
void FFT();
int INPUT[SAMPLENUMBER]DATA[SAMPLENUMBER];
float fWaveR[SAMPLENUMBER]fWaveI[SAMPLENUMBER]w[SAMPLENUMBER];
float sin_tab[SAMPLENUMBER]cos_tab[SAMPLENUMBER];
main()
{
int i;
InitForFFT();
MakeWave();
for ( i=0;i {
fWaveR[i]=INPUT[i];
fWaveI[i]=0.0f;
w[i]=0.0f;
}
FFT(fWaveRfWaveI);
for ( i=0;i {
DATA[i]=w[i];
}
while ( 1 ); // break point
}
void FFT(float dataR[SAMPLENUMBER]float dataI[SAMPLENUMBER])
{
int x0x1x2x3x4x5x6xx;
int ijkbpL;
float TRTItemp;
/********** following code invert sequence ************/
for ( i=0;i {
x0=x1=x2=x3=x4=x5=x6=0;
x0=i&0x01; x1=(i/2)&0x01; x2=(i/4)&0x01; x3=(i/8)&0x01;x4=(i/16)&0x01; x5=(i/32)&0x01; x6=(i/64)&0x01;
xx=x0*64+x1*32+x2*16+x3*8+x4*4+x5*2+x6;
dataI[xx]=dataR[i];
}
for ( i=0;i {
dataR[i]=dataI[i]; dataI[i]=0;
}
/************** following code FFT *******************/
for ( L=1;L<=7;L++ )
{ /* for(1) */
b=1; i=L-1;
while ( i>0 )
{
b=b*2; i--;
} /* b= 2^(L-1) */
for ( j=0;j<=b-1;j++ ) /* for (2) */
{
p=1; i=7-L;
while ( i>0 ) /* p=pow(27-L)*j; */
{
p=p*2; i--;
}
p=p*j;
for ( k=j;k<128;k=k+2*b ) /* for (3) */
{
TR=dataR[k]; TI=dataI[k]; temp=dataR[k+b];
dataR[k]=dataR[k]+dataR[k+b]*cos_tab[p]+dataI[k+b]*sin_tab[p];
dataI[k]=dataI[k]-dataR[k+b]*sin_tab[p]+dataI[k+b]*cos_tab[p];
dataR[k+b]=TR-dataR[k+b]*cos_tab[p]-dataI[k+b]*sin_tab[p];
dataI[k+b]=TI+temp*sin_tab[p]-dataI[k+b]*cos_tab[p];
} /* END for (3) */
} /* END for (2) */
} /* END for (1) */
for ( i=0;i {
w[i]=sqrt(dataR[i]*dataR[i]+dataI[i]*dataI[i]);
}
} /* END FFT */
void InitForFFT()
{
int i;
for ( i=0;i {
sin_tab[i]=sin(PI*2*i/SAMPLENUMBER);
cos_tab[i]=cos(PI*2*i/SAMPLENUMBER);
}
}
void MakeWave()
{
int i;
for ( i=0;i {
INPUT[i]=sin(PI*2*i/SAMPLENUMBER*3)*1024;
}
}
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 260 2005-07-18 20:17 Lab0503-FFT\Debug.lkf
文件 260 2005-07-18 20:17 Lab0503-FFT\Debug.lkv
文件 2207 2005-07-18 20:17 Lab0503-FFT\FFT.c
文件 1167 2007-05-09 14:14 Lab0503-FFT\FFT.paf
文件 870 2007-05-09 14:14 Lab0503-FFT\FFT.pjt
文件 13579 2005-07-18 20:48 Lab0503-FFT\FFT.wks
文件 334 2007-05-09 14:14 Lab0503-FFT\c6713.cmd
文件 134 2005-07-18 20:35 Lab0503-FFT\cc_build_Debug.log
文件 5794 2005-07-18 20:17 Lab0503-FFT\Debug\FFT.map
文件 7686 2005-07-18 20:17 Lab0503-FFT\Debug\FFT.obj
文件 17569 2005-07-18 20:17 Lab0503-FFT\Debug\FFT.out
目录 0 2008-04-06 08:57 Lab0503-FFT\Debug
目录 0 2008-04-06 08:57 Lab0503-FFT
----------- --------- ---------- ----- ----
49860 13
- 上一篇:使用C语言编写的病房管理系统
- 下一篇:没有了
相关资源
- FFT电流检测 stm32f4
- 基于STM32的音频信号采样与FFT分析
- 高效FFT的C/C++代码实现包括基2的DIF和
- FFT(快速傅里叶变换)C++源程序
- 数值分析算法程序设计 C++实现
- FFT算法c语言代码
- 信号FFT变换及绘图
- 利用FFT计算频谱图
- FFT的C语言实现
- FFTW3.0快速傅立叶变换库
- 利用c语言实现FFT运算
- 基于VC++的FFT频谱分析
- 数据结构算法设计C++,乐学答案
- 图像FFT卷积
- Frogs_MFCCs.csv
- C语言实现2的指数次点的FFT
- DSP 77个C语言程序 包括FFTFIR IIR等的D
- fftw-3.3.7快速傅里叶变换计算库
- 郑宗汉著算法设计与分析第2版
- picc8.05编译器+注册机完全版+PIC C语言
- 算法设计实验报告-求最大子段和问题
- FFT和DFT完整算法C语言实现
- VB.NET写的FFT算法
- 最大团问题C语言算法设计与分析
- 算法设计实验报告-快速排序和归并排
- 算法设计分析中 图的m色的着色问题
- 自制FFT频谱分析软件
- 快速傅里叶变换的详细源代码C++
- C++ 分治法解决邮局选址问题
- 码头扩建问题
评论
共有 条评论