FFT快速傅里叶变换 C程序 fft快速傅里叶变换

#include "myapp.h"

#include "ICETEK-VC5509-EDU.h"

#include "scancode.h"

#include <math.h>

#define PI 3.1415926

#define SAMPLENUMBER 128

void InitForFFT();

void MakeWave();

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<SAMPLENUMBER;i++ )

{

fWaveR[i]=INPUT[i];

fWaveI[i]=0.0f;

w[i]=0.0f;

}

FFT(fWaveR,fWaveI);

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

{

DATA[i]=w[i];

}

while ( 1 ); // break point

}
FFT(快速傅里叶变换)C程序 fft快速傅里叶变换

void FFT(float dataR[SAMPLENUMBER],float dataI[SAMPLENUMBER])

{

int x0,x1,x2,x3,x4,x5,x6,xx;

int i,j,k,b,p,L;

float TR,TI,temp;

/********** following code invert sequence ************/

for ( i=0;i<SAMPLENUMBER;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<SAMPLENUMBER;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(2,7-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<SAMPLENUMBER/2;i++ )

{

w[i]=sqrt(dataR[i]*dataR[i]+dataI[i]*dataI[i]);

}

} /* END FFT */

void InitForFFT()

{

int i;

for ( i=0;i<SAMPLENUMBER;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<SAMPLENUMBER;i++ )

{

INPUT[i]=sin(PI*2*i/SAMPLENUMBER*3)*1024;

}

}

 

 

这里面有些头文件可能找不到,这是我们老师提供的参考程序,这里只是给个参考,下次有时间我再注释下。

#include"myapp.h"

void CLK_init()

{

ioport unsigned int *clkmd;

clkmd=(unsigned int *)0x1c00;

*clkmd =0x2033; // 0x2033;//0x2413;// 200MHz=0x2513

}

void SetDSPPLL(unsigned int uPLL)

{

ioport unsigned int *clkmd;

clkmd=(unsigned int *)0x1c00;

*clkmd =uPLL;

}

void TMCR_reset( void )

{

ioport unsigned int *TMCR_MGS3=(unsigned int *)0x07FE;

ioport unsigned int *TMCR_MM =(unsigned int *)0x07FF;

*TMCR_MGS3 =0x510;

*TMCR_MM =0x000;

}

  

爱华网本文地址 » http://www.413yy.cn/a/25101013/175264.html

更多阅读

vc6.0使用教程 vc6.0中文版使用教程

vc6.0使用教程——简介  Visual C++ 6.0 是微软在 98 年推出的程序设计套件,由于其设计界面友好,启动快,对机器要求低、所以一直被推荐作为 C++ 程序设计课程的首选工具,本文以简单的 WIN32 终端程序设计过程为例讲解 VC6 的使用。vc

C#程序实现Canny边缘检测算法 边缘检测算法实现

转载自:http://blog.csdn.net/yjz_uestc/article/details/6664937Canny边缘检测是被公认的检测效果最好的边缘检测方法,是由JohnF.Canny于1986年提出,算法目标是找出一个最优的边缘检测的方法,所谓最优即:1.好的检测:算法能够尽可能的标

module_init和module_exit的作用 module init tools

就像你写C程序需要包含C库的头文件那样,Linux内核编程也需要包含Kernel头文件,大多的Linux驱动程序需要包含下面三个头文件:#include <linux/init.h>#include <linux/module.h>#include <linux/kernel.h>其中,init.h 定义了驱动的初始化

C程序设计---练习题四

第9章编译预处理1.以下程序的输出结果是()A、9B、6C、36D、18#define f(x) x*xmain( ){ int a=6,b=2,c;c=f(a) / f(b);

C语言视频教程30课曾怡 c语言教学视频曾怡

最好的c语言视频教程----曾怡副教授讲解30课,csf格式C程序设计视频教程(曾怡):本套视频教程由曾怡副教授讲解,使用教材为:《C语言程序设计》 谭浩强清华大学出版社出版。是难得的C语言学习视频教程。全程共30讲,每讲45分钟左右,CSF视频格

声明:《FFT快速傅里叶变换 C程序 fft快速傅里叶变换》为网友桐哥变砖哥分享!如侵犯到您的合法权益请联系我们删除