C语言实现基2FFT算法-频域.doc_第1页
C语言实现基2FFT算法-频域.doc_第2页
C语言实现基2FFT算法-频域.doc_第3页
C语言实现基2FFT算法-频域.doc_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

/*基二频域算法*/#include math.h#include stdio.hstruct compx double real; double imag; compx ;struct compx EE(struct compx b1,struct compx b2)struct compx b3;b3.real=b1.real*b2.real-b1.imag*b2.imag;b3.imag=b1.real*b2.imag+b1.imag*b2.real;return(b3);void FFT(struct compx *xin,int N)int f,m,LH,nm,i,k,j,L;double p , ps ;int le,B,ip;float pi;struct compx v,w,t;LH=N/2; f=N;for(m=1;(f=f/2)!=1;m+); /*2m=N*/for(L=m;L=1;L-) /*这里和时域的也有差别*/ le=pow(2,L);B=le/2; /*每一级碟形运算间隔的点数*/pi=3.14159;for(j=0;j=B-1;j+) p=pow(2,m-L)*j; ps=2*pi/N*p; w.real=cos(ps); w.imag=-sin(ps); for(i=j;i=N-1;i=i+le) ip=i+B; t=xini; xini.real=xini.real+xinip.real; xini.imag=xini.imag+xinip.imag; xinip.real=xinip.real-t.real; xinip.imag=xinip.imag-t.imag; xinip=EE(xinip,w); /*变址运算*/nm=N-2; j=N/2;for(i=1;i=nm;i+)if(i=k)j=j-k;k=k/2;j=j+k; /*main programe*/#include #include #include float result257;struct compx s257;int Num=16;const float pp=3.14159;main()int i;for(i=0;i16;i+)si.real=sin(pp*i/32);si.imag=0;FFT(s,Num);for(i=0;i16;i+)printf(%.4f,si.real);printf(+%.4fj
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论