傅里叶变换及逆变换

     关于FFT，书上已经给出了实现方法；曾在研2时也使用迭代法实现了自己的FFT，速度上要慢一些，但是理解起来要容易一些；

     最近看书，发现了一些以前没有注意到的问题；比如，FFT产生是到底是什么呢？是频率的信息吗？完整吗？程序表现出来的结果到底正确吗？等等一些问题；以前没有考虑过。

     今天来给出答案，当然是本人的一些个人理解，不一定正确!

一，FFT产生的到底是什么？

书上曾经把FFT后的信息用一幅图像来表示，其实这很勉强；目的只是让大家了解频域到底反应了一些什么东西，比如，边缘上在频域上就可以很直观的反应出来；别说一幅图像了，就是多幅图像也不一定能表示得了，因为频域变换出来的数据不是整型，如果强制为整型，则会造成数据的丢失；除此之外，更重要的一点，一幅图像的频域只反应了一个频谱，也就是相位并没有考虑进去，这样，通过FFT产生的频域来还原时域图像是不可能的。

二，这些频率信息完整吗？

如上所述，这些信息是不完整的，我曾试图使用一幅24位的bmp来表示频率的信息，但是最终因为图像的整数表示而放弃，当然也是可以做的，只是效果要差很多；以后有机会我会试一试的。频域应该是一个复数；在我看来，频域的的完整信息应该有两部分表示，一部分是实部，一部分是虚部；

三，以前图像处理程序中给出的实现正确吗？

首先，我应该说是正确的，因为频率上确实是这样的；但是，以前程序在为了更加细节描述频域的信息，并没有除M*N，取代除了个100；这样更能突出细节；

 

     关于反向FFT的算法其实很简单，使用前向变换来计算就很方便!我试着用频域信息还原时域图像，有很少的地方失真，失真的原因是由于计算机在计算时浮点数的省略；

     我还原时域图像的方法是这样的：使用FFT对时域图像进行变换，得到频域的信息，保存在两个txt文件当中，这两个文件分别是频域的实部和虚部!然后使用前向变换，用txt文件来构建时域图像。

代码如下：

.h文件

 // FFT.h: interface for the CFFT class. // // #if !defined(AFX_FFT_H__C3BE732F_8D82_4870_B422_25580A5BABC9__INCLUDED_) #define AFX_FFT_H__C3BE732F_8D82_4870_B422_25580A5BABC9__INCLUDED_ #if _MSC_VER > 1000 #pragma once #endif // _MSC_VER > 1000 /********************************************** *类名:FFT及反向FFT *目的:时域到频域，频域到时域的变换 *作者:无忧狂澜 *日期:2010-02-19 **********************************************/ class CFFT { public: BOOL ComputeConverseFFT(LONG lWidth, LONG lHeight, double *fFReal, double *fFImag, double *fTReal, double *fTImag, bool bCenter); BOOL ConverseFourier(LPSTR lpDIBBits, LONG lWidth, LONG lHeight, double *fFReal, double *fFImag); BOOL Fourier(LPSTR lpDIBBits, LONG lWidth, LONG lHeight,double *fReal,double *fImag,bool bCenter); CFFT(); virtual ~CFFT(); }; #endif // !defined(AFX_FFT_H__C3BE732F_8D82_4870_B422_25580A5BABC9__INCLUDED_)

 

.cpp文件

 // FFT.cpp: implementation of the CFFT class. // // #include "stdafx.h" #include "mydip.h" #include "FFT.h" #include "math.h" #include <complex> using namespace std; #define PI 3.1415926 #define TWOPI 6.2831852 #ifdef _DEBUG #undef THIS_FILE static char THIS_FILE[]=__FILE__; #define new DEBUG_NEW #endif // // Construction/Destruction // /************************************************************************* * * 函数名称： * FFT() * * 参数: * complex<double> * TD - 指向时域数组的指针 * complex<double> * FD - 指向频域数组的指针 * r - 2的幂数，即迭代次数 * * 返回值: * 无。 * * 说明: * 用迭代法实现的快速付立叶变换。 * 速度比蝶形算法要慢一些，但理解起来容易 ************************************************************************/ /* void init(int r) { f=new complex<double>[1<<r]; int tempIndex,i,j; double angle; W=new complex<double>*[1<<r]; for(tempIndex=1;tempIndex<=r;tempIndex++) { W[1<<tempIndex]=new complex<double>[(1<<(r-1))-1]; }