FFT应用——傅立叶变换实验报告

《FFT应用——傅立叶变换实验报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、FFT应用——傅立叶变换实验报告1.实验名称：FFT应用——傅立叶变换2.实验目的：1．加深对DFT算法原理和基本性质的理解。2．熟悉FFT算法原理和FFT子程序的应用。3．利用FFT算法对序列信号进行变换及逆变换。3.实验原理：从纯粹的数学意义上看，傅立叶变换是将一个函数转换为一系列周期函数来处理的。从物理效果看，傅立叶变换是将图像从空间域转换到频率域，其逆变换是将图像从频率域转换到空间域。4、实验内容及步骤1.读入长度为N的序列信号。2．调用信号产生子程序产生实验信号。3.利用函数FFT1D,对其进行快速傅立叶变换,F1=fft1d(f).4.显示变换

2、后的实验数据。5.对变换后的信号，利用函数IFFT1D,对其进行傅立叶逆变换,F2=fft1d(f).6.显示变换了中心后的数据，比较和原来的输入信号是否相同。5.程序清单：傅立叶变换函数fft1dvoidfft1d(intflag,intn,doublefr[],doublefi[],doubletblSin[],doubletblCos[]){inti,m,iw,j=0,l,lp,lp2,n2,k;doublec,s,wr,wi,xa,ya;for(i=0;i

3、;ya=fi[i];fi[i]=fi[j];fi[j]=ya;}n2=n/2;while(j>=n2){j=j-n2;n2=n2/2;}j+=n2;}m=0;n2=n;while(n2!=1){m+=1;n2=n2/2;}for(l=1;l<=m;l++){lp=(int)pow(2.0,l);lp2=lp/2.0;k=0;for(j=0;j

4、]*c;fr[iw]=fr[i]-wr;fi[iw]=fi[i]-wi;fr[i]=fr[i]+wr;fi[i]=fi[i]+wi;}}}if(flag==1){for(i=0;i

5、s[i]=cos(arg);}}6.程序运行结果：变换前的信号傅立叶变换后的数据傅立叶逆变换后的数据7.实验分析：信号经过傅立叶变换后，输出的数据实部和虚部分别对称，符合傅立叶变换的性质。再经过逆变换后，得到的数据符合原来的输入数据，说明算法正确。