北京邮电大学《数字信号处理》门爱东-DSP03-离散傅里叶变换及其快速计算方法.pptVIP

门爱东教授menad@ 数字信号处理 Digital Signal Processing 3.1 问题的提出：连续信号的傅里叶变换 非周期连续信号 xa(t)，其傅里叶变换为 xa(t) 为时域连续信号。 Xa(Ω) 为频域连续信号。 3.1 问题的提出：离散信号的变换 离散信号在两种变换域中的表示方法： 1）离散时间傅里叶变换 DTFT -- 提供了绝对可加的离散时间序列在频域（ω）中的表示方法。 2）Z 变换 -- 提供任意序列的 z 域表示。 3.1 问题的提出：可计算性 问题：X(z)，X(ejω) 都是连续的，利用计算机处理有困难，例如使用 Matlab，因此 提出了在频域内取样，使频谱离散化的问题； 必须截断序列，得到有限个点的序列。 目标：我们需要得到一个可进行数值计算的变换 方法： 1）DTFT -? 频域中原始信号频谱的周期拓展 2）对 DTFT 在频域中采样 --? DFS。 3）将 DFS 推广到有限持续时间序列 ? DFT。 （DFT 避免了前面提到的那两个问题，并且它是计算机可实现的变换方式。） DFT 已成为 DSP 算法中的核心变换，原因： 成为有限长序列傅里叶变换的重要方法。 有快速算法。 3.1 问题的提出：傅里叶变换的四种形式 非周期连续时间—傅里叶变换(FT)－连续频率 周期连续时间—傅里叶级数(FS)－离散频率 非周期离散时间—离散时间傅里叶变换(DTFT)－连续频率 周期离散时间—离散傅里叶级数(DFS)－离散频率 3.1 问题的提出：傅里叶变换的四种形式 2. 周期连续时间信号：傅里叶级数 FS 时域连续函数造成频域是非周期的谱。 频域的离散对应时域是周期函数。 3. 非周期离散信号：离散时间傅里叶变换 DTFT 时域的离散化造成频域的周期延拓 时域的非周期对应于频域的连续 4. 周期离散时间信号：离散傅里叶级数 DFS 一个域的离散造成另一个域的周期延拓 离散傅里叶级数的时域和频域都是离散的和周期的 四种傅里叶变换形式的归纳总结： 3.2.1 DFS 定义：预备知识 3.2.1 DFS 定义：正变换 为了推导 的关系，作下列变量代换： 时域： 频域： 则得： 3.2.1 DFS 定义：正变换 周期离散序列的 Z 变换存在（收敛）的问题 因为周期离散序列， 而对于周期信号，严格数学意义上讲，其 Z 变换不收敛，因为： 而对于 找不到衰减因子使它绝对可和（收敛）。为此，定义新函数，其 Z 变换： 3.2.1 DFS 定义：正变换 其频谱： 3.2.1 DFS 定义：正变换 频域取样 X(ejω) 是连续变量 ω 的周期函数，周期为2π。把ω 离散化，即在0～2π区间内等间隔取 N 个点，取样间隔为 2π/N。 另一个角度看， X(ejω) 是 Z 平面单位圆上的 Z 变换。连续变量 ω 的离散化也可以认为是把单位圆分 N 等份，每份为 2π/N 。 其中： 称为频域中的取样间隔， 也称为频率分辨率。 3.2.1 DFS 定义：正变换 3.2.1 DFS 定义：正变换 3.2.1 DFS 定义：反变换 反变换IDFS 正变换两端乘以 ，m=0,1,…,N-1 然后令 k=0，1，…，N-1 求和，得： 3.2.1 DFS 定义：反变换 在什么条件下不产生混迭失真？ —频率取样 频率取样：若时间信号有限长，当满足下列条件时，X(ejω) 的样本值 X(k) 能不失真的恢复成原信号。 为了避免时间上的混迭： （1）必须是时间限制（有限时宽，有 M 个样值） （2）为了防止时域混迭，频域取样造成的时域周期 T 的持续时间应该大于等于信号时宽 M，假设频域 2π周期内取 N 点，则 ，即取样频率间隔 频率分量 如果信号 DFS 可表示为： 因此，时域 n 及频域 k 都是有物理意义的。 更具体地，傅里叶系数的标号 k 和频率 f 的关系为： 所以： 对应关系： 傅里叶系数标号k ：0～N 数字频率ω：0～2π