DSPFFT试验报告

1、DSP技术课程设计(FFT)作者姓名：钟欣专业名称：电子科学与技术指导教师：张双目录一、实验结果3正弦信号输入的结果3余弦信号输入的结果4心得5参考文献6附件1正弦信号程序7正弦信号文件7文件11附件2余弦信号程序12余弦信号文件12实验结果正弦信号输入的结果余弦信号输入的结果心得通过本次的DSP设计课程，我学会了在CCS软件平台 上运用C语言编程來实现傅里叶运算。在课程中老师讲解了 的基本原理，以及如何通过运用C语言编程來实现。实验过程中，通过对傅里叶变换的深入了解，以及老师 对程序的讲解，我成功的实现了正弦和余弦的傅里叶变换运 算，完成了实现要求。这次课程让我对傅里叶变换有了更加 深入的了

2、解。在使用CCS软件时，要注意加上文件，避免岀错；还有就是文件要放在工程文件夹里，不然要出错，还有就是 iiit和float之间的区分。我经过此次的学习后，无论是傅里叶变换的原理，还是 DSP编程方面，都有了很大的提升。这对以后的课程学习提 供了坚实的基础。我也从中获益良多。参考文献1 杨明广.C程序设计教材.四川科学技术出版社,2012.2 邹彦DSP原理及应用（修订版）.电子工业岀版社,2002.3 吴京.信号分析与处理（修订版）.电子工业出版社,2014.附件1正弦信号程序正弦信号文件* * * * * * * * * * * * * * * * * * */include o #mcl

3、ude o # mclude,HI#defiiie PI#define SAMPLENUMBER 128 void FFT();void IiiitFoiFFTO；void MakeWave();float INPUTSAMPLENUMBER,DATASAMPLENUMBER;float fVaveRSAMPLENUMBER,fWaveISAMPLENUMBER,wSAMPLENUMBER;float sin_tabSANIPLENUlVIBER,cos_tabSAMPLENUMBER;mam()int 1；IiutForFFTO；MakeWave();foi(i=0;iSAMPLENUMBE

4、R;i+ )fVaveRi=INPUTi;fVaveIi=;wi=;FFT(fWaveR,fVaveI);for(i=0;iSAMPLENUMBER;i+)DATAi=wi;wlule(l);void FFT(float dataRSAMPLENUMBER,float dataISAMPLENUMBER)int xO,x 1 ,x2,x3,x4,x5,x6,xx;mt i,j,k,b,p,L；float TR,TI,temp;for( i=O;iSAMPLENUMBER;i-H-)x0=xl=x2=x3=x4=x5=x6=0;x0=i&0x01; xl=(i/2)&0x01 ; x2=(i/4

5、)&0x01 ; x3=(i/8)&0x01; x4=(i/16)&0x01; x5=(i/32)&0x01;x6=(i/64)&0x01;xx=x0*64+xl *32+x2* 16+x3*8+x4*4+x5*2+x6; dataIxx=dataRi;for( i=O;i0 )b=b*2; i-;for(j=0;j0)p=p*2;i-;ip=p*i;for(k=j ;kv 12 8 ;k=k+2 *b)fITR=dataRk ;TI=dataIk ;temp=dataRk+b; dataRk=dataRk+dataRk+b*cos_tabp+dataIk+b*sm_tabp; datalk=

6、datalk -dataRk+b *sin_tabp+datalk+b *cos_tabp; dataRk+b=TR-dataRk+b *cos_tabp -datalk+b *sin_tab p;dataIk+b=TI+temp*sin_tabp-dataIk+b *cos_tabp;iI JfoT(i=0 ；1VS AMPLENUMBER/2 ;i+)wi=sqrt(dataRi *dataRi+dataIi *datali);void IiMtFoiFFT()lilt 1;foi(i=0;iSAMPLENUMBER;i+) sm_tabi=sin(PI*2*i/SAMPLENUMBER)

7、;cos_tabi=cos(PI*2*i/SAMPLENUMBER); void MakeWave()int 1；foi(i=0;iPROGlPAGE0.cinit: PROG1 PAGE 0.switch: PROGlPAGE0.vectois:VECTPAGE0.const: DARAM1PAGE 1.bss :DARAM1 PAGE 1.stack : DARAM2PAGE 1.system: DARAM2PAGE 1.data :DARAM2PAGE 1/* * * * * * * * * * *End of File* * * * * * * * */附件2余弦信号程序余弦信号文件/

8、* * * * * * * * * */#mclude o #mclude o # mcludeHU #define PI #ckflnc SAMPLENUMBER 128void FFT();void IiiitFoiFFTO；void MakeWave();float INPUTSAMPLENUMBER,DATASAMPLENUMBER; float fWavcRS AMPLENUMBER,fWaveISAMPLENUMBER,wSAMPLENUMBER; float sin_tabSAlVIPLENUlVIBER,cos_tabSAMPLENUMBER; inam()int 1；Iiut

9、FoiFFTO；MakeWave();for(i=0;iSAMPLENUMBER;i+ )fWaveRi=INPUTi;fVaveIi=;wi=;FFT(fWaveR,fVaveI); for(i=0;iSAMPLENUMBER;i+) DATAi=wi;wlule(l);void FFT(float dataRSAMPLENUMBER,float dataISAMPLENUMBER) mt x(xl.x2、x3,x4,x5、x6,xx;uit ij,k,bp，L;float TR.TI.temp;for( i=O;iSAMPLENUMBER;i+ )x0=xl=x2=x3=x4=x5=x6=

10、0;x0=i&0x01; xl=(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+xl *32+x2* 16+x3*8+x4*4+x5*2+x6; dataIxx=dataRi;for( i=O;i0 )b=b*2; i-;for(j=0;j0)p=p*2;i-;p=p*j；for(k=j ;kv 12 8 ;k=k+2 *b)TR=dataRk ;TI=dataIk ;temp=dataRk+b;dataRk=dataRk+dataRk+b *c

11、os_tabp+datalk+b *sin_tabp;datalk=datalk -dataRk+b *sin_tabp+datalk+b *cos_tabp;dataRk+b=TR-dataRk+b*cos_tabp-dataIk+b*sin_tabp;dataIk+b=TI+tcmp*sin_tabp -datalk+b *cos_tabp; for(i=0;iSAMPLENUMBER/2;i-l-+)wi=sqrt(dataRi *dataRi+dataIi *datali);void IiiitFoiFFT()lilt 1;foi(i=0;KSAMPLENUMBER;i+) sm_tabi=sin