快速傅里叶变换_快速傅里叶变换的计算机实现基于实数

最新推荐文章于 2023-04-13 19:16:28 发布

weixin_39807541

最新推荐文章于 2023-04-13 19:16:28 发布

阅读量313

点赞数

文章标签：快速傅里叶变换快速傅里叶变换fft

本文链接：https://blog.csdn.net/weixin_39807541/article/details/111299139

版权

本文提供了一个基于实数的快速傅里叶变换(FFT)算法代码，输入和输出均为实数。代码简洁高效，适用于8192个点的正变换和逆变换，计算时间小于1毫秒。作者分享了邮箱和微信号，鼓励读者转发分享，以支持其继续贡献技术力量。

摘要由CSDN通过智能技术生成

上回文章发布的快速傅里叶变换是基于查表法计算的，本次发布一个自带蝶形变换的算法代码。

本次发布的代码，函数的输入和输出均是实数。

由于增加了碟形变换单元准备，因此代码所需的单词量从85个增加到94个。

对的，总共94个单词

//基2的fft代码
//res为输出数组，data为输入数组，size为数据长度,sign=1为FFT,-1为IFFT

void fft(float *res,float *data,int size,int sign)
{
int n=size;
float *x,*y;
int i,j,k,l,m,n1,n2;
float c,c1,e,s,s1,t,tr,ti;

m=log((float)n)/log(2.0f)+0.2;

x=new float[n];
y=new float[n];
memset(x,0,n*4);
memset(y,0,n*4);
memcpy(x,data,size*4);

n1=n-1;
for(j=0,i=0;i {
  if(i  {
   tr=x[j];
   ti=y[j];
   x[j]=x[i];
   y[j]=y[i];
   x[i]=tr;
   y[i]=ti;
  }
  k=n/2;
  while(k  {
   j=j-k;
   k=k/2;
  }
  j=j+k;
}

n1=1;
for(l=1;l<=m;l++)
{
  n1=2*n1;
  n2=n1/2;
  e=3.14159265359/n2;
  c=1.0;
  s=0.0;
  c1=cos(e);
  s1=-sign*sin(e);

  for(j=0;j  {
   for(i=j;i   {
    k=i+n2;
    tr=c*x[k]-s*y[k];
    ti=c*y[k]+s*x[k];
    x[k]=x[i]-tr;
    y[k]=y[i]-ti;
    x[i]=x[i]+tr;
    y[i]=y[i]+ti;
   }
   t=c;
   c=c*c1-s*s1;
   s=t*s1+s*c1;
  }
}
if(sign==-1)
{
  for(i=0;i  {
   res[i]=x[i]/2;
  }
}
else
{
  int len=n;
  float adjust=2.0f/size;
  for(i=0;i   res[i]=sqrt(x[i]*x[i]+y[i]*y[i])*adjust;
}
delete x;
delete y;
x=NULL;
y=NULL;
return;
}