一个关于128点的快速傅立叶的C语言程序
来源：学生作业帮助网 编辑：六六作业网 时间：2024/12/18 17:01:48
一个关于128点的快速傅立叶的C语言程序
一个关于128点的快速傅立叶的C语言程序
一个关于128点的快速傅立叶的C语言程序
这是我写的1024点的快速傅里叶变换程序,下面有验证,你把数组
double A[2049]={0};
double B[1100]={0};
double powerA[1025]={0};
改成 A[256]={0};
B[130]={0};
power[129]={0};就行了,
void FFT(double data[], int nn, int isign)
的程序可以针对任何点数,只要是2的n次方
具体程序如下：
#include
#include "math.h"
#include
#include
#include
#include
#include
void FFT(double data[], int nn, int isign)
{
//复数的快速傅里叶变换
int n,j,i,m,mmax,istep;
double tempr,tempi,theta,wpr,wpi,wr,wi,wtemp;
n = 2 * nn;
j = 1;
for (i = 1; i i)
{
tempr = data[j];
tempi = data[j + 1];
data[j] = data[i];
data[j + 1] = data[i + 1];
data[i] = tempr;
data[i + 1] = tempi;
}
m = n / 2;
while (m >= 2 && j > m)
{
j = j - m;
m = m / 2;
}
j = j + m;
}
mmax = 2;
while( n > mmax )
{
istep = 2 * mmax; //这里表示一次的数字的变化.也体现了级数,若第一级时,也就是书是的第0级,其为两个虚数,所以对应数组应该增加4,这样就可以进入下一组运算
theta = -6.28318530717959 / (isign * mmax);
wpr = -2.0 * sin(0.5 * theta)*sin(0.5 * theta);
wpi = sin(theta);
wr = 1.0;
wi = 0.0;
for( m = 1; m