第3章 离散傅利叶变换(DFT).ppt

第三章 离散傅利叶变换（DFT） 由DFS导出DFT的过程 离散傅里叶变换(DFT)： 有限长序列的傅里叶变换称为离散傅里叶变换，简写为DFT。 由此得DFT定义： 正变换 反变换 注意 在 离 散 傅 里 叶 变 换 关 系 中 , 有 限 长 序 列 都 作 为 周 期 序 列 的 一 个 周 期 来 表 示 , 都 隐 含 有 周 期 性 意 义 . P78 频域循环卷积定理 序列及其DFT的实、虚、偶、奇关系 （3）结论 长度为M的有限长序列，频域抽样不失真的条件：频域抽样点数N要大于或等于序列长度M，即满足N≥M。此时可得到： （4）例子 P 89 频域抽样：看一个矩形序列，频域抽样是指对时域已是离散，频域仍是连续信号。现在频域上进行抽样处理，使其频域也离散化。 解：频域抽样，按N=5点，频域抽样，时域延拓相加……,时域延拓的周期个数等于频域的抽样点数N=5，由于N=M，所以时域延拓恰好无混叠现象。 按N=4时进行抽样，由于N=4，而序列长度为M=5，NM,时域延拓后产生混叠现象。（原信号为红色，延拓取主值区间后的恢复信号为兰色。） 3.3.3 频域内插公式 从频域抽样不失真条件可以知道： N个频域抽样X(k)能不失真的还原出长度为N的有限长序列x(n)。 那么用N个X(k)也一定能完整地表示出x(n)整个X(z)以及频率响应X(ejω)。 过程很简单,先把N个X(k)作IDFT得到x(n),再把x(n)作Z变换便得到X(z)。 频域响应（傅利叶变换）的内插公式（P 90） 各抽样点之间的X(ejw)值,则由各抽样点的加权内插函数在所求点上的值的叠加而得到。 信号的频谱分析：计算信号的傅里叶变换 例子： 有一频谱分析仪用的FFT处理器，其抽样点数必须是2的整数幂。假定没有采用任何特殊的数据处理措施，已给条件为： (1)频率分辨率≤10Hz (2)信号的最高频率≤2.5kHz 试确定以下参量： （1）最小记录长度Tp； （2）抽样点的最大时间间隔T； （3）在一个记录中的最少点数N。 解：（1）由分辨力的要求确定最小记录长度Tp Tp=1/F=1/10=0.1(s) 故最小记录长度为0.1秒。 （2）从信号的最高频率确定最大的抽样时间间隔T fs≥2fc， T=1/fs≤1/2fc=0.2*10-3 (s) （3）最小记录点数N，它应满足： N≥2fc /F = 500 该处理器所需最少采样点数为N=29=512点。 频谱泄漏 改善方法： 栅栏效应 改善方法： 增加频域抽样点数N（时域补零），使谱线更密 1. 用DFT对连续非周期信号进行谱分析 P 97 在工程应用中经常碰到的是连续非周期信号的傅里叶变换，即连续非周期性的频谱函数，数字计算机难于处理的,因而我们采用DFT对其进行逼近。 推导： 设连续信号xa(t)持续时间为Tp，最高频率为fc， xa(t)的傅里叶变换为： X( jf ) 是f 的连续周期函数 对X(jf )在[0，fs]上等间隔采样N点，采样间隔为F 设共采样N点 同理，由 可以推出： 连续信号的频谱特性可以通过对连续信号采样并进行DFT再乘以T的近似方法得到； 时域采样信号可通过（3.4.8）得出。 (3.4.7) (3.4.8) （1）栅栏效应 利用DFT逼近连续时间信号的傅里叶变换,其频谱将不再是连续函数而是离散点的频谱，谱线仅出现在基频F的整数倍处。 用DFT计算频谱,就如通过一个栅栏观看一个景色,只能在离散点的地方看到真实的景象，从而产生栅栏效应。 如果在两离散的谱线间,频谱有很大变化,不作特殊处理,则无法将其检测出来。 减小栅栏效应方法 使频域抽样更密。即增加频域抽样点数N,在不改变时域数据的情况下，必然是在时域数据末端增加一些零值点，使一个周期内的点数增加，但并不改变原有的记录数据。频域抽样为 ，N增加，必然使样点间距更近，谱线更密，谱线变密后原来看不到的谱分量就可能看到了。 0 k |X(k)| 0 x(n) n 0 x(n) n 谱线细化 （2）频率响应的混叠失真及参数选择 信号最高频率为fc，按抽样定理，抽样频率应满足（抽样前用前置滤波器将高于fc频率的信号分量滤除，防止混叠） 采样间隔： 频率采样间隔: 频率分辨率用频率采样间隔F描述，表示谱分析中能够分辨的两个频谱分量最小间隔。（P98） 信号的最高频率分量fc与频率分辨率F间有着矛盾 关系，要想fc增加，为