利用DFT/FFT对序列位移—非整数时的陷阱

一、引言

已知，信号$x[n]$及其离散傅里叶变换$X[k]$，则有
$$\begin{array}{rl}X[k]& =\sum _{n=0}^{N-1}x[n]e^{-j2\pi nk/N}\\ x[n]& =\frac{1}{N}\sum _{k=0}^{N-1}X[k]e^{j2\pi nk/N}\end{array}$$

同时，离散傅里叶变换（DFT）及其快速实现（FFT），具有“位移性质”，即，
$$\begin{array}{rl}DFT[x[n-n_0]]& =X[k]e^{-j2\pi n_{0}k/N}\\ IDFT[X[k-k_0]]& =x[n]e^{j2\pi n{k}_0/N}\end{array}$$

本文主要论点：
1、当$n_0,k_0\in Z$时，上式成立。
2、当$n_0,k_0\notin Z$时，上式不一定成立。
3、当$n_0,k_0\notin Z$时，若使上式成立，需要保证信号$x[n]$的频谱$X[k]$不存在“频率泄露（频率分量没有位于采样网格上时产生）”和“频谱折叠（采样率不满足奈奎斯特采样定理时产生）”，即$X[k]$非常理想时，上式才能成立。

二、理论推导

(1) 时域解释
假设位移长度$n_0\notin Z$，设$x[n]$位移后的信号$y[n]$，并设其离散傅里叶变换分别为$Y[k]$，则利用离散傅里叶变换位移性质，可得，
Y [ k ] = X [ k ] e − j 2 π n 0 k / N = ∑ n = 0 N − 1 x [ n ] e − j 2 π ( n + n 0 ) k / N y [ n ] = 1 N ∑ k = 0 N − 1 Y [