在拉格朗日插值法与牛顿插值多项式
中有说明当给定n+1个点值序列 ( x 0 , y 0 ) ， ( x 1 , y 1 ) ， ⋅ ⋅ ⋅ ， ( x n , y n ) (x_0,y_0)，(x_1,y_1)，···，(x_n,y_n) (x0,y0)，(x1,y1)，⋅⋅⋅，(xn,yn)，其中任意两个 x x x都互不相同时，有且只有一个n次多项式函数包含这些点值，换句话说，一个n次多项式可以由 n + 1 n+1 n+1个点值来表示，如果我们需要获取一个多项式的点值表示，只需要选取 n + 1 n+1 n+1个互不相同的 x x x代入即可，这样做的复杂度是O(n^2)，有没有可能更快呢？考虑代入单位复数根 e 2 π i / n e^{2\pi i/n} e2πi/n，选择单位复数根代入获取点值表示的这样一种方式，和离散傅里叶变换(DFT)完全相同，或者说这就是离散傅里叶变换，采用快速傅里叶变换(FFT)可以将DFT复杂度降低至O(nlogn)，也就是说可以采用快速傅里叶变换将多项式点值表示的复杂度降至O(nlogn)
快速傅里叶变换
设已有的一个多项式函数为 P ( x ) = a 0 + a 1 x + a 2 x 2 + ⋅ ⋅ ⋅ + a n x n P(x)=a_0+a_1x+a_2x^2+···+a_nx^n P(x)=a0+a1x+a2x2+⋅⋅⋅+anxn， w n i = e 2 π i / n w_n^i=e^{2\pi i/n} wni=e2πi/n，n是已知的要变换的序列的长度，并且一般限制为2的幂次，如果多项式的次数不够，就往后补一些系数为0的幂次。现在我们的目的是获取点值序列： ( P ( w n 0 ) , P ( w n 1 ) , P ( w n 2 ) , ⋅ ⋅ ⋅ , P ( w n n − 1 ) ) (P(w_n^0),P(w_n^1),P(w_n^2),···,P(w_n^{n-1})) (P(wn0),P(wn1),P(wn2),⋅⋅⋅,P(wnn−1))
P ( w n i ) = ∑ j = 0 n a j ( w n i ) j P(w_n^i) = \sum_{j=0}^{n}a_j(w_n^i)^j P(wni)=∑j=0naj(wni)j，这样算完整个序列显然复杂度为 O ( n 2 ) O(n^2) O(n2)，考虑将这个求和式按照奇偶性分组： P ( w n i ) = ∑ j = 0 n 2 − 1 a 2 j ( w n i ) 2 j + ∑ j = 0 n 2 − 1 a 2 j + 1 ( w n i ) 2 j + 1 P(w_n^i) = \sum_{j=0}^{\frac{n}{2}-1}a_{2j}(w_n^i)^{2j}+\sum_{j=0}^{\frac{n}{2}-1}a_{2j+1}(w_n^i)^{2j+1} P(wni)=j=0∑2n−1a2j(wni)2j+j=0∑2n−1a2j+1(wni)2j+1
由于复数根有性质： w n i = w n 2 i 2 w_n^i=w_{\frac{n}{2}}^{\frac{i}{2}} wni=w2n2i， ( w n i ) 2 j + 1 = ( w n i ) ( w n i ) 2 j = w n i w n 2 i j (w_n^i)^{2j+1}=(w_n^i)(w_n^i)^{2j}=w_n^iw_{\frac{n}{2}}^{ij} (wni)2j+1=(wni)(wni)2j=wniw2n
最低0.47元/天 解锁文章
扫一扫
827
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
