第6章离散信号与系统的频域分析.ppt

显然，这样的方程式用矩阵表示更为简明方便。一般情况下， 我们可将DFT的定义式写成矩阵的形式: 它们都是列矩阵。WE称为DFT矩阵 . 它们都是列矩阵，WE称为DFT矩阵，当N=8时 在分析一个系统时， 往往是确定的，因而WE矩阵的关键是矩阵中 的指数排列次序。为此，我们将WE矩阵中各元素 的指数单独组成矩阵E，使问题更加简单明了。当N=8时 同理，IDFT也可写成矩阵形式。 设N=22=4，则有 n=0, 1, 2, 3 式中， 令n的排序为0， 2， 1， 3， 则可写出其矩阵表示式为 可见，在n按0, 2, 1, 3排序情况下，DFT矩阵为 每个子矩阵可以写成下面的形式： 矩阵WE可以写为 若令 则 稀疏矩阵 6.7.2 基2FFT概述 对于N=22=4，据DFT定义式(6.7 - 1)有 利用DFT矩阵因子化的结论式(6.7 - 10)，可将上式写成 列矩阵F1(k)： 利用中间矢量F1(k)就可进一步完成式(6.7-12)的计算： 元素F2(0)可用一次复数乘法和一次复数加法确定： 直接算法的计算时间对FFT算法的计算时间之比有下列近似的关系： 图 6.7-1 N=22的FFT信号流程图 例如，对于结点F1(2)有 表 6.2 自然顺序与二进制逆序(N=8) 6.8 离散系统的频域分析 6.8.1 基本信号ejωk激励下的零状态响应 对一任意的周期离散信号f(k)，利用离散傅里叶级数可以将其表示为指数信号  的线性组合，即 式中： 同样，利用离散时间傅里叶变换可以将任一非周期离散信号f(k)表示为指数信号ejωk的线性组合，即 式中： 因此，与连续信号的情况一样，我们将指数信号ejωk称为基本信号。指数信号 实质上与基本信号ejωk一样，它只不过是当 时的特例。 设稳定离散LTI系统的单位响应为h(k)，则据上一章的讨论可知， 系统对基本序列e jωk的零状态响应为 式中的求和项正好是h(k)的离散时间傅里叶变换，记为H(ejω)，即 称H(ejω)为传输函数或频率响应。 一个稳定的离散LTI系统，对ejωk这一基本信号的零状态响应是基本信号ejωk本身乘上一个与时间序数k无关的常系数H(ejω)， 而H(ejω)为系统单位响应h(k)的离散时间傅里叶变换。H(ejω)一般是ω的连续函数， 而且是复函数， 即 称为系统的幅频响应或幅度响应，φ(ω)称为系统的相频响应或相位响应。 设输入周期正弦序列为 应用欧拉公式，可将f(k)写成 f(k)通过一个频响函数为H(ejω)的离散LTI系统的稳态响应可表示为 式中： (6.8 - 10) 由于|H(ejω)|为ω的偶函数，而φ(ω)为ω的奇函数，因而式(6.8 - 10)可写为 例 6.8-1 已知描述一离散LTI系统的差分方程为 若输入正弦序列 求该系统的稳态响应ys(k)。 解 则有 由输入正弦序列的表达式可知，其ω=π/2，所以有 则该离散系统的稳态响应为 若输入正弦序列改 , 即ω改为π，则据 从而其输出的稳态响应为 6.8.2 一般信号f(k)激励下的零状态响应 该系统的零状态响应yf(k)为 应用离散时间傅里叶变换的有关性质，上式在频域中可以表示为 对一个离散LTI系统而言，其输出y(k)与输入f(k)间的N阶线性常系数差分方程一般具有如下形式： 式中，ai和bi都是常数。若系统是稳定的，对上式两端求离散时间傅里叶变换，并应用离散时间傅里叶变换的时移性质， 可以得到 从而得到该系统的频率响应 例 6.8-2 描述一稳定离散LTI系统的差分方程为 解 若给该系统施加信号f(k)=(0.5)kε(k)，则该系统的零状态响应yf(k)据式(6.8-13)即可求得。 由于 则有 同理，将Y(ejω)展开成部分分式有 利用式(6.2 - 13)的变换对， 可得零状态响应 若f(k)是一个有限长序列，系统的h(k)也是一个有限长序列，则在满足由循环卷积求线卷积的条件下，它们的DFT满足 此时，由DFT利用FFT算法可以求得系统的响应y(k)。 周期序列fp(k)的离散时间傅里叶级数表示式为 如果将NFn表示成Fp