FFT快速傅立叶变换-Matlab和c语言实现

毕业设计做的是基于SOPC的数字频谱分析仪，所以用了一下Matlab。其实数字信号处理学的很差。

FFT来历很复杂，但是他的结果很简单。但凡xxx变换，其实就是一堆数经过运算的到另外一堆数。

对于FFT，输入的数的个数和输出的数的个数相等，这就更妙了。

我们把采样到的数据放入一个虚数数组，数组下标是采样的时间点，计算完成后结果还是存储在这个数组中，

但是下标的含义变成了频率。

Matlab自带的fft函数功能很全很强大，但是这上一个自制的m程序：

function y=myfft(xr,n)

%信号长度为2的整数次方，基于时间抽选的FFT

p=0:n-1; %开始倒位序
nu=log2(n);
p1=p;
b=zeros(1,n);
for t=1:nu;
    p2=floor(p1/2);
    b=b*2+(p1-2*p2);
    p1=p2;
end;
yr(p+1)=xr(b+1);
xr=yr;   %倒位序结束
t=0:n/2-1; %计算因子w开始（只计算w0到w n/2-1）
for v=0:n/2-1;
    w=exp(-2*i*pi*t/n);
end; %计算因子w结束
for m=1:nu;% 计算x（k）开始
    h=2^(m-1);
    k=1;
    while(k<n+1)
        for t=1:h;
            y=bitshift(k-1,nu-m,nu)+1;%求w的幂次数
            xch(k)=xr(k)+w(y)*xr(k+h);
            k=k+1;
        end;
        for t=1:h;
            y=bitshift(k-1-h,nu-m,nu)+1;%求w的幂次数