转自:http://blog.renren.com/share/408963653/15068964503
其实我感觉这个学期算法最难最搞不懂的绝对不是动态规划啊!绝对是快速傅里叶变换啊!最近才弄懂有木有。
有不少人问我,于是干脆就写成日志吧。
首先明确一下基本概念吧,就三点,DFT,FFT,蝴蝶操作。
DFT(离散傅里叶变换):书上写的最清楚的一句话叫做,向量y=(y0,y1,……yn-1)是系数向量a=(a0,a1,,……,an-1)的离散傅里叶变换,也写作y=DFTn(a)。说白了,就是求n个y值,但是n个自变量x的取值很特殊
FFT(快速傅里叶变换):这个名词好理解,就是很快地算出这n个y值。一般我们计算n个n次的多项式值需要O(n^2)的时间。现在用FFT可以减少到O(n*logn)。具体原理一两句话说不完……要考的话也太理论性了…
蝴蝶操作:这个应该是让大家最费解的…他就是一个FFT得实际应用…因为这是实践的东西所以要掌握。书上有一个图,虽然感觉看书不一定看的怎么懂,但是这个图还是要记住。不过我觉得从下往上写看起来更舒服。
转过来以后,向左边是加,向右边是减。谨记啊!
突然发现好难解释,先把宋老师课件里面上次坑爹的没有答案的课后练习拿来看看。
我这里a0,a1,a2……就不调换顺序了,其实是一样的,格式不同罢了。
a=(0,1)
经过FFT,Y=(1,-1),这个就是书上最简单的应用了。
再来一个简单的 a=(1,0).
y=(1,1).
我觉得要考试一般就是考四个的吧,两个太简单了,8个的太复杂了(我后面有写)
例子
a=(1,1,0,1)
y=(3,1,-1,1) 注意算出来是 y0,y2,y1,y3的顺序。
a=(0,1,2,3)
y=(6,-2-2i, -2 , -2+2i).
注意右上角的蝶形运算中的旋转因子变成了(恩!?怎么不能插公式啊!)w(1/4)=i。。。那是因为在a0,a2的蝶形运算中我们多乘了(n/2),这里n=2.(这个我解释不来…不过这个数是固定的,背就行了)。
那我上个八个的,有点复杂,估计可能不会考
a=(0,1,2,3,4,5,6,7).
y的顺序如图
告诉大家一个比较好的办法判断是否算对了,直接人肉使用O(n^2)算法就好了,举四个的例子。就是把x得值带入多项式 y=a3*x^3+a2*x^2+a1*x+a0算出y值什么的。x0=1,x1=i,x2=-1,x3=-i带入。这样的话,其实考察算没算对完全是个 伪命题吗……
再提醒一次,左加右减哦!
最后再说一下这玩意有什么用,什么信号学的废话就不说了。无非是FFT用N*LogN时间算出点值方便用点值法算出多项式乘法的系数结果,这个大家看书上510的图我觉得就能理解了。插值神马的,背公式吧……虽然我也不记得了。
时间有点紧迫,感觉写的巨烂无比,大神们如果发现写错了留言给我立马改,关键不要误人子弟……