2018年FFT变换相关公式、IFFT变换(FFT逆变换).docVIP

FFT变换相关公式、IFFT变换(FFT逆变换) 　　FFT 变换公式 　　(A )利用 Bulestein 布鲁斯坦所提出的等式: 　　nk =1 　　[n2+k 2? k ? n 2] 　　则: 　　X Z k = x(n)N ? 1 　　n=0A 　　? n 　　ω1[n2+k2 　　? k ? n 2 　　] 　　= 　　x(n)N ? 1 　　n=0 　　A ? n ω12 　　ω? (k? n) 2ω12 = 　　ω12 　　x(n)N ? 1 　　n=0 　　A ? n ω1n 2 　　ω? 　　(k? n) 2 　　(B )令 g n =x(n)A? n ωn 22 　　, n=0,1,…, N-1 　　h n =ω 　　? 　　n 2 , n=0,1,…, M-1 　　可知:Z k 点的 Z 变换值 X Z k 可以通过 g n 和 h n 的离散卷积值并乘 ω1k 2得 　　到,即: 　　X Z k =ωk 2[g n ? h n ] , k=0,1,…, M-1 　　g n 与 h n 的离散卷积,可以用 g n 和 h n 的圆周卷积来实现,而圆周卷积 可用 FFT 的方法来求得。 　　(C ) Z k =A ω? k , k=0,1,…, M-1 　　A =A 0e j θ0 ω=ω0e 　　? j φ0 = Z k =A 0e j θ0? ω0? k e j φ0k =A 0ω0? k e j(θ0+φ0k) =Z 0=A 0e j θ0 , Z 1=A 0ω0? 1e j(φ0+θ0) Z k =A 0ω0? k e j(θ0+φ0k) 　　Z M ? 1=A 0ω0? (M? 1) e j(θ0+φ0(M? 1)) (D )当满足下面特殊条件: 　　(a) M=N 　　(b) A =A 0e j θ0=1 ,即 A 0=1, θ0=0° 　　(c) ω=ω0e 　　? j φ0 　　=e 　　? j 　　2π ,即 ω0=1, φ0= 　　2πN 　　此时 Z k 为均匀分布在单位圆周上,即由 CZT 变换求出该序列的 DFT ,此时: 　　X Z k = x(n)Zk ? n N ? 1 　　n=0= x n [Aω? k ]? n N ? 1 　　n=0 　　　　= 　　x(n)(e? j 　　2π 　　) kn 　　N ? 1 　　n=0 　　= x(n)e 　　? j 　　2πN ? 1 　　n=0 　　　　(E ) h n =ω 　　? 　　n 2可以想象为频率随时间(n )成线性增长的复指数序列: 　　∵ h n = 　　[ω0e ? j φ0]? n 2 　　2 　　=ω0 　　? 　　n 22 　　? e 　　j φ0 　　n 22 　　　　(F )计算 DFT 的具体步骤: 　　①加长补零将 g(n)变成列长为 L 的序列:L =2m g n =x n A ? n ω 　　n 2=x n e 　　? j 2π 　　n 2=x(n)e 　　? j 　　πn 2 0≤ n ≤ N ? 1 　　=g n = x n e 　　? j πn 2 0≤ n ≤ N ? 10 N≤ n ≤ L ? 1 L ≥ 2N ? 1且 L =2m ②求 g(n)的 DFT 　　利用 FFT 方法求 g (n )序列的 DFT 　　G r = g(n)e 　　? j 　　2πrn L ? 1 　　n=0 　　0≤ r ≤ L ? 1 　　③ h(n)补零加长,周期延拓成 L 点序列: 　　h n = ω? 　　n 22 　　0≤ n ≤ N ? 1 [ω=e ? j 　　2πN 　　]? 　　n 22 　　=e 　　j 　　πn 2N 　　0 N≤ n ≤ L ? N ω 　　? 　　(L? N ) 2 　　2 　　L? N +1≤ n ≤ L ? 1 　　　　④求 h(n)的 DFT 　　利用 FFT 方法求 h(n)序列的 DFT 　　H r = h n e ? j 2π 　　rn 0≤ r ≤ L ? 1L ? 1 　　n