matlab实现FFT正变换和逆变换, matlab的逆变换效果真是完美

eric2013

代码：

1. Fs = 1024;              % 采样率
   
2. N  = 1024; % 采样点数
   
3. n  = 0:N-1;           % 采样序列
   
4. t  = 0:1/Fs:1-1/Fs;     % 时间序列
   
5. f = n * Fs / N;          %真实的频率
   
6. 
   
7. x = 1.5*sin(2*pi*20*t+pi/3) ;  %原始信号
   
8. y = fft(x, N);    %对原始信号做FFT变换
   
9. z = conj(y);     %对转换结果取共轭
   
10. 
    
11. subplot(2,1,2);
    
12. z = fft(z, N);     %再次做FFT
    
13. k = conj(z);      %对转换结果去共轭
    
14. plot(f,  real(k));  %绘制转换后的波形
    
15. title('FFT逆变换后的波形');
    
16. 
    
17. subplot(2,1,1);
    
18. plot(f,  x);       %绘制原始波形
    
19. title('原始波形');

复制代码

eric2013

傅里叶反变化：调用FFT的例子计算IFFT
将FFT的结果取共轭，然后再次调用FFT，将结果取共轭。将所得结果除以N就是所要
void conj(complex *xx)                           //将数据取共轭
{
        int i;
        for(i=0;i<size_x;i++)
        xx[i].img=-(xx[i].img);   
}   
void ifft(complex *xx)
{
        int i;   
        conj(xx);                                          //FFT变换后的X(k)取共轭
        fft(xx);                                          //FFT变换
        conj(xx);                                          //将变换后的结果再次取共轭
        for(i=0;i<size_x;i++)                  //将所得到的结果除以N 就是所要的原始数据
    {
          xx[i].real=xx[i].real/size_x;
          xx[i].img=xx[i].img/size_x;
        }
}