fft（快速傅里叶变换）学习草稿，逆dft证明

最新推荐文章于 2024-09-06 17:21:28 发布

置顶 YxuanwKeith

最新推荐文章于 2024-09-06 17:21:28 发布

阅读量2.5k

点赞数 9

分类专栏：算法-快速傅里叶变换

本文链接：https://blog.csdn.net/yxuanwkeith/article/details/56510757

版权

算法-快速傅里叶变换专栏收录该内容

0 篇文章 0 订阅

订阅专栏

先上最原始的式子（模 $n$ 意义下）：

c_{i} = \sum_{j = 0}^{i} a_{j} * b_{i - j}

$c_i=\sum_{j=0}^{i}a_j*b_{i-j}$
变形一下：

c i = \sum j \sum k a j * b k * [- i + j + k = = 0] (1)

$c_i=\sum_{j}\sum_{k}a_j*b_k*[-i+j+k==0](1)$
考虑怎么判断一个数是否等于0，观察下面的式子：

[p = = 0] = \sum n - 1 i = 0 w i p n (w 是 单 位 复 数 根)

$[p==0]=\frac{\sum_{i=0}^{n-1}w^{ip}}{n}(w是单位复数根)$
现在要证明它是对的
1.当

p=0 p = 0 $p=0$ 时，

i∗p=0 i ∗ p = 0 $i*p=0$ ，所以

wip=1 w i p = 1 $w^{ip}=1$ ，而总共有

n n $n$ 项，除

n

$n$ 后恰好等于1。
2.当

p≠0 p ≠ 0 $p\not=0$ 时，直接用等比数列求和。

S n = a 1 * 1 - q n 1 - q ， 代 数 进 去 S n = w 0 * 1 - w p n 1 - w p ， 由 于 w n = 1 ， 所 以 S n = 0

$S_n=a1*\frac{1-q^n}{1-q}，代数进去 \\ S_n=w^0*\frac{1-w^{pn}}{1-w^p}，由于w^n=1，所以 \\ S_n=0$
得证。

那么把 $[p==0]=\frac{\sum_{i=0}^{n-1}w^{ip}}{n}$ 代入 $(1)$

c i = \sum j \sum k a j * b k * [- i + j + k = = 0] c i = 1 n \sum j \sum k a j * b k \sum l = 0 n - 1 w - i l * w j l * w k l c i = 1 n \sum l = 0 n - 1 w - i l * (\sum j = 0 n - 1 a j * w j l) * (\sum k = 0 n - 1 b k * w k l)

$c_i=\sum_{j}\sum_{k}a_j*b_k*[-i+j+k==0] \\ c_i=\frac{1}{n}\sum_{j}\sum_{k}a_j*b_k\sum_{l=0}^{n-1}w^{-il}*w^{jl}*w^{kl} \\ c_i=\frac{1}{n}\sum_{l=0}^{n-1}w^{-il}*(\sum_{j=0}^{n-1}a_j*w^{jl})*(\sum_{k=0}^{n-1}b_k*w^{kl})$
对应到dft上，

∑n−1l=0∑n−1j=0aj∗wjl ∑ l = 0 n − 1 ∑ j = 0 n − 1 a j ∗ w j l $\sum_{l=0}^{n-1}\sum_{j=0}^{n-1}a_j*w^{jl}$ 和

∑n−1l=0∑n−1k=0bk∗wkl ∑ l = 0 n − 1 ∑ k = 0 n − 1 b k ∗ w k l $\sum_{l=0}^{n-1}\sum_{k=0}^{n-1}b_k*w^{kl}$ 相当于对

a,b a , b $a,b$ 数组各做一次dft，而为了最后得到

ci c i $c_i$ 还要把

a,b a , b $a,b$ 对应的点值乘起来，设

dl=(∑n−1j=0aj∗wjl)∗(∑n−1k=0bk∗wkl) d l = ( ∑ j = 0 n − 1 a j ∗ w j l ) ∗ ( ∑ k = 0 n − 1 b k ∗ w k l ) $d_l=(\sum_{j=0}^{n-1}a_j*w^{jl})*(\sum_{k=0}^{n-1}b_k*w^{kl})$ ，即