快速傅里叶变换fft的c语言算法彻底研究 pdf,快速傅里叶变换fft变换

快速fft变换

快速傅里 叶变换FFT的C语言算法彻底研究LED音乐频谱显示的核心算法就是快速傅里叶变换，FFT的理解和编程还是比较难的，特地撰写此文分享一下研究成果。

一、彻底理解傅里叶变换

快速傅里叶变换(Fast Fourier Transform)是离散傅里叶变换的一种快速算法，简称FFT，通过FFT可以将一个信号从时域变换到频域。模拟信号经过A/D转换变为数字信号的过程称为采样。为保证采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成分的2倍，这称之为采样定理。假设采样频率为fs，采样点数为N，那么FFT结果就是一个N点的复数，每一个点就对应着一个频率点，某一点n(n从1开始)表示的频率为：fn=(n-1)*fs/N。举例说明：用1kHz的采样频率采样128点，则FFT结果的128个数据即对应的频率点分别是0，1k/128，2k/128，3k/128， ，127k/128 Hz。这个频率点的幅值为：该点复数的模值除以N/2(n=1时是直流分量，其幅值是该点的模值除以N)。

二、傅里叶变换的C语言编程

1、对于快速傅里叶变换FFT，第一个要解决的问题就是码位倒序。假设一个N点的输入序列，那么它的序号二进制数位数就是t=log2N.码位倒序要解决两个问题：①将t位二进制数倒序；②将倒序后的两个存储单元进行交换。如果输入序列的自然顺序号i用二进制数表示，例如若最大序号为15，即用4位就可表示n3n2n1n0，则其倒序后j对应的二进制数就是n0n1n2n3，那么怎样才能实现倒序呢？利用C语言的移位功能！

程序如下，我不多说，看不懂者智商一定在180以下！复数类型定义及其运算 #define N 64 //64点

#define log2N 6 //log2N=6

/*复数类型*/

typedef struct

{

float real;

float img;

}complex;

complex xdata x[N]; //输入序列

/*复数加法*/

complex add(complex a,complex b)

{

complex c;

c.real=a.real+b.real;

c.img=a.img+b.img;

return c;

}

/*复数减法*/

complex sub(complex a,complex b)

{

complex c;

c.real=a.real-b.real;