《【教学课件】第三章离散付里叶变换(DFT)DiscreteFourierTransform.ppt》由会员分享,可在线阅读,更多相关《【教学课件】第三章离散付里叶变换(DFT)DiscreteFourierTransform.ppt(169页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第三章第三章离散付里叶变换离散付里叶变换(DFT)Discrete Fourier Transform第一节第一节 引引 言言一、序列分类一、序列分类对一个序列长度未加以任何限制,则一个序列可分为:无限长序列:n=-或n=0或n=-0有限长序列:0nN-1有限长序列在数字信号处理是很重要的一种序列。由于计算机容量的限制,只能对过程进行逐段分析。二、二、DFT引入引入由于有限长序列,引入DFT(离散付里叶变换)。DFT它是反映了“有限长有限长”这一特点的一种有用工具。DFT变换除了作为有限长序列的一种付里叶表示,在理论上重要之外,而且由于存在着计算计算机机DFT的有效快速算法的有效快速算法-FF
2、T,因而使离散付里叶变换(DFT)得以实现,它使DFT在各种数字信号处理的算法中起着核心的作用。三、本章主要讨论三、本章主要讨论离散付里叶变换的推导离散付里叶变换的有关性质离散付里叶变换逼近连续时间信号的问题序列的抽取与插值第二节付里叶变换的几种形式一、引言傅 里 叶 变 换:建 立 以 时时 间间t 为 自 变 量 的“信 号”与 以 频频 率率 f为 自 变 量 的“频 率 函 数”(频谱)之 间 的 某 种 变 换 关 系.所 以“时 间”或“频 率”取 连 续 还 是 离 散 值,就 形 成 各 种 不 同 形 式 的 傅 里 叶 变 换 对。在 深 入 讨 论 离 散 傅 里 叶 变
3、 换 D F T 之 前,先 概 述 四种 不 同 形式 的 傅 里 叶 变 换 对 二、四种不同付里叶变换对1、傅 里 叶 级 数(FS):连 续 时 间,离 散 频 率 的 傅 里 叶 变 换。2、傅 里 叶 变 换(FT):连 续 时 间,连 续 频 率 的 傅 里 叶 变 换。3、序 列 的 傅 里 叶 变 换(DTFT):离 散 时 间,连 续 频 率 的 傅 里 叶 变 换.4、离 散 傅 里 叶 变 换(DFT):离 散 时 间,离 散 频 率 的 傅 里 叶 变 换假定数字频率为w,模拟频率为。1.傅 里 叶 级 数(FS)周期连续时间信号 非周期离散频谱函数。周期为T0的周期
4、性连续时间函数 x(t)可展成傅里叶级数X(jk0),是离散非周期性频谱,表 示为:FS例子通过以下 变 换 对 可 以 看 出 时 域 的 连 续 函 数 造 成 频 域 是 非 周 期 的 频 谱 函 数,而 频 域 的 离 散 频 谱 就 与 时 域 的 周 期 时 间 函 数 对 应.(频域采样,时域周期延 拓)2.傅 里 叶 变 换(FT)非周期连续时间信号通过连续付里叶变换(FT)得到非周期连续频谱密度函数。例子以下变换对可以看出时域 连 续 函 数 造成频域是非周期的谱,而时域的非周期造成频域是连续的谱.3.序 列 的 傅 里 叶 变 换(DTFT)非周期离散的时间信号(单位圆上
5、的Z变换(DTFT)得到周期性连续的频率函数。例子同样可看出,时域的离散造成频域的周期延拓,而时域的非周期对应于频域的连续.4.离 散 傅 里 叶 变 换(DFT)上面讨论的三种傅里叶变换对,都不适用在计算机上运算,因为至少在一个域(时 域 或 频 域)中,函 数 是 连 续 的.因 为 从 数 字 计 算 角 度,我 们 感 兴 趣 的 是 时 域 及 频 域 都 是 离 散 的 情 况,这 就 是 我 们 这 里 要 谈 到 的 离 散 傅 里 叶 变 换.周期性离散时间信号从上可以推断:周期性时间信号可以产生频谱是离散的 离散时间信号可以产生频谱是周期性的。得出其频谱为周期性离散的。也即
6、我们所希望的。总之,一个域的离散必然造成另一个域的周期延拓。其中正变换:反变换:二、四种付里叶变换形式的归纳周期()和连续离散(T)和非周期周期()和离散离散(T)和周期非周期和连续连续和非周期非周期和离散连续和周期(T)频谱频谱函数时间时间函数作业参看程佩青的光盘第三章的第一节付里叶变换的四种可能形式的测验题。第三节离散付里叶级数(DFS)我 们 先 从 周 期 性 序 列 的 离 散 傅 里 叶 级 数(DFS)开 始 讨 论,然 后 再 讨 论 可 作 为 周 期 函 数 一 个 周 期 的 有 限 长 序 列 的 离 散 傅 里 叶 变 换(DFT).一、引言二、DFS定义反变换:三、
7、DFS离散付里级数的推导意义用数字计算机对信号进行频谱分析时,要求信号必须以离散值作为输入,而且上面讨论可知:只有第四种形式(DFS)对数字信号处理有实用价值。但如果将前三种形式要么在时域上采样,要么在频域上采样,变成离散函数,就可以在计算机上应用。所以我们要先了解如何从以上三种形式推出DFS.1.由非周期连续时间信号推出DFSx(t)经过抽样为x(nT),对离散的时间信号进行DTFT得到周期连续频谱密度函数。再经过抽样,得到周期性离散频谱密度函数即为DFS.x(t)t取样x(t)tDTFTX(ejT)采样X(ejw)w2.周期性连续时间信号函数周期性连续时间信号函数经采样后,得到周期性的离散
8、时间函数(DFS)。x(t)X(ejw)tw采样x(n)nDFS3.非周期离散时间信号非周期离散时间信号经过序列付里时变换(即单位园上的Z变换)DTFT,得到周期连续谱密度函数,再经采样为周期离散频谱密度函数(DFS)。x(t)tX(ejT)wX(ejw)DTFT采样四、推导DFS正变换 以下由第三种付里叶级数形式为例推导出离散付里叶级数变换。非周期信号x(n),其DTFT(单位园上Z变换)为 其为周期连续频谱密度函数,对其频域进行采样,使其成为周期性离散频谱函数。设在一周期内采样N个点,则两采样点间距为即得出DFS的正变换:得到各抽样频点频率为:代入DTFT式子中,这时由于抽样,信号变成周期
9、离散信号 ,得五、DFS的反变换解:已知两边同乘以 ,并对一个周期求和根据正交定理用n替换r,可得:即得:六、回顾DFS设 x(n)为周 期 为 N 的 周 期 序 列,则 其 离 散 傅 里 叶 级 数(DFS)变 换 对 为:正 变 换反变 换第四节离散付里叶级数的性质一、引言可以由抽样Z变换来解析DFS,它的许多性质与Z变换性质类似。它们与Z变换主要区别为:(1)与 两者具有周期性,与Z变换不同。(2)DFS在时域和频域之间具有严格的对偶关系。它们主要性质分为:线性、序列移位(循环移位)、调制性、周期卷积和假设令 和 皆是周期为N的周期序列,它们各自的DFS为(1)线性 其中a,b为任意
10、常数,所得到的频域序列也是周期序列,周期为N。(2)序列移位(循环、移位)时域证明:令 i=n+m,得证毕(3)调制性频域证明:(4)时域卷积 周期卷积和与以前卷积不同,它的卷积过程限在一个周期内称为周期卷积。频域相乘等于时域卷积(指周期卷积)。频域相乘等于时域卷积(指周期卷积)。频域:则有:相乘时域卷积 证明:代入:则:(5)频域卷积时域:相乘周期卷积频域:作业第133页:第1、2题第五节离散付里叶变换DFT一、由DFS引出DFT的定义周期序列实际上只有有限个序列值才有意义,因 而它的离散傅里叶级数表示式也适用于有限长 序列,这就得到有限长序列的傅里叶变换(DFT).具体而言,我们把:(1)
11、时域周期序列看作是有限 长序列x(n)的周期延拓;(2)把频域周期序列看作是有限长序列X(k)的周期延拓.(3)这样我们只要把DFS的定义式两边(时域、频域)各取主值区间,就得到关于有限长序列的时频域的对应变换对.这就是数字信号处理课程里最重要的变换-离 散 傅 里 叶 变 换(DFT).二、DFT定义正变换反变换X(k)、x(n)为有限长序列的离散付里叶变换对,已知其中一个序列就能确定另一个序列。注意在 离 散 傅 里 叶 变 换 关 系 中,有 限 长 序 列 都 作 为 周 期 序 列 的 一 个 周 期 来 表 示,都 隐 含 有 周 期 性 意 义.三、DFT涉及的基本概念1.主 值
12、(主值区间、主值序列)2.移 位(线性移位、圆周移位)3.卷 积(线性卷积、圆周卷积)4.对 称(序列的对称性、序列的对称分量)5.相 关(线性相关、圆周相关)1.主 值(主值区间、主值序列)2.移位线线 性性 移移 位:位:序 列 沿 坐 标 轴 的 平 移.圆周移位:圆周移位:将 有 限 长 序 列 x(n)以 长 度 N 为 周 期,延 拓 为 周 期 序 列,并 加 以 线 性 移 位 后,再 取 它 的 主 值 区 间 上 的 序 列 值,m 点 圆 周 移 位 记 作:其 中(.)N 表 示 N 点 周 期 延 拓.(1)有 限 长 序 列 圆 周 移 位 的 实 现 步 骤(2)
13、例子12131 0.5(1)周期延拓:N=5时nx(n)2131x(n)0.521310.51120.5n(2)周期延拓:N=6 时,补零加长2131x(n)0.521310.51123n32131 0.5nx(n)(3)m=1时,左移(取主值)131x(n)0.52(4)m=-2时,右移(取主值)2131nx(n)0.5n 3.卷 积卷积在此我们主要介绍:(1)线性卷积(2)圆周卷积(3)圆周卷积与线性卷积的性质对比(1)线性卷积线 性 卷 积 定 义:有 限 长 序 列 x1(n),0nN1-1;x2(n),0nN2-1则 线 性 卷 积 为 注意:线 性 卷 积 结 果 长 度 变 为
14、N1+N2-1.(2)圆周卷积令则圆 周 卷 积 结 果 长 度 不 变,为 N.圆 周 卷 积 的 实 现 步 骤例子线性卷积与圆周卷积步骤比较231x(n)54n0N1=5213h(n)n0N2=3得到线性卷积线性卷积结果的示意图14265ny(n)2014830 5 4 3 2 1 1 2 3 15 12 9 6 3 10 8 6 4 2 5 4 3 2 1 5 14 26 20 14 8 3N=7231x(n)54n0N1=5213h(n)n0N2=3(1)圆周卷积圆周卷积:(N=7)补零加长 231x(k)540N=7k213h(k)k0N2=3231h(k)0k(2)圆周卷积需进行
15、周期延拓,而线卷积无需周期延拓:圆卷积的反折(并取主值区间):231231231h(-k)k0(3)平移0231h(1-k)k(4)相乘x(k)h(-k)=51=5x(k)h(1-k)=5*2+4*1=14x(k)h(2-k)=5*3+4*2+3*1=26231x(k)540N=7kx(k)h(3-k)=4*3+3*2+2*1=20 x(k)h(4-k)=3*3+2*2+1*1=14x(k)h(5-k)=2*3+1*2=8x(k)h(6-k)=1*3=3231h(-k)k(5)相加得到圆周卷积的示意图14265ny(n)2014830可见可见,线性卷积与圆周卷积相同线性卷积与圆周卷积相同(当当
16、NN1(5)+N2(3)-1=7时时)课后练习用图表求解圆卷积 x(k)=5,4,3,2,1,h(n)=1,2,3,同上求N=7点的圆卷积。解:(1)将x(n)补零加长为x(k)=5,4,3,2,1,0,0,(2)将h(n)补零加长至N=7,并周期延拓,(3)反折得到:h(-k)=1,0,0,0,0,3,2(4)作图表结果同上。若圆周卷积取长度为N=5,则求圆周卷积231x(k)540N=5k231h(-k)k0求得圆周卷积x(k)h(-k)=5*1+2*3+1*2=13x(k)h(1-k)=5*2+4*1+1*3=17x(k)h(2-k)=5*3+4*2+3*1=26x(k)h(3-k)=4
17、*3+3*2+2*1=20 x(k)h(4-k)=3*3+2*2+1*1=14看出圆卷积与线卷积不同.171326y(n)n02014用图表求解圆卷积 x(k)=5,4,3,2,1,h(n)=1,2,3,同上求N=5点的圆卷积。解:(1)x(n)无需补零加长x(k)=5,4,3,2,1,(2)将h(n)补零加长至N=5,并周期延拓,(3)反折得到:h(-k)=1,0,0,3,2(4)作图表1713262014y(n)n0作业2P133 第3,4,7,8,9,10题参看程佩青的光盘中第三章的离散付里叶图形的测验第1第2题(3)圆 周 卷 积 与 线 性 卷 积 的 性 质 对 比 圆周卷积线性卷
18、积是针对FFT引出的一种表示方法表示方法信号通过线性系统时,信号输出等于输入与系统单位冲输入与系统单位冲激响应的卷积激响应的卷积两序列长度必须相等相等,不等时按要求补足零值点补足零值点。两序列长度可以不等不等。如x1(n)为 N1点,x2(n)为 N2点卷积结果长度与两信号长度相等皆为N卷积结果长度为N=N1+N2-14.对称分为:(1)序列的对称性(2)序列的对称分量(1)序列的对称性(a)奇 对 称(序 列)和 偶 对 称(序 列)(b)圆 周 奇 对 称(序 列)和 圆 周 偶 对 称(序 列)(c)共 轭 对 称(序列)和 共 轭 反 对 称(序 列)(d)圆 周 共 轭 对 称(序列
19、)和 圆 周 共 轭 反 对 称(序 列)(a)奇 对 称(序 列)和 偶 对 称(序 列)4、满 足xe(n)=xe(-n)的 序 列 xe(n)称 为 偶偶 对对 称称 序序 列列1、x(n)与-x(-n)互称为奇对称奇对称。2、满足x0(n)=-x0(-n)的序列x0(n)称为称为奇对称序列奇对称序列。3、x(n)与 x(-n)互称为 偶偶 对对 称称;例子0 xe(n)n0 x(n)n0y(n)=x(-n)nx(n)与y(n)互为偶对称为偶对称序列0 x(n)n0 x(-n)n互为奇对称0 xo(n)n为奇对称序列(b)圆 周 奇 对 称(序 列)和 圆 周 偶 对 称(序 列)1、长
20、 度 为N的 有 限 长 序 列 x(n)与y(n)=-x(-n)NRN(n)互 为 圆圆 周周 奇奇 对对 称称.2、长 度 为 N 的 有 限 长 序 列x(n)若 满 足 x(n)=-x(-n)NRN(n)则x(n)是 圆圆 周周 奇奇 对对 称称 序序 列列.x(n)y(n)=-x(-n)NRN(n)x(n)与y(n)互 为 圆圆 周周 奇奇 对对 称称.圆 周 奇 对 称圆 周 奇 对 称(序 列)x(n)4、长 度 为 N 的 有 限 长 序 列 xe(n),若 满 足 x(n)=x(-n)NRN(n)则 x(n)是 圆圆 周周 偶偶 对对 称称 序序 列列.3、长 度 为 N 的
21、 有 限 长 序 列 x(n)与y(n)=x(-n)NRN(n)互 为 圆圆 周周 偶偶 对对 称称.圆 周 偶 对 称(序 列)周期延拓判断 序列的圆周奇偶对称性的简便方法在n=N处补上补上与n=0处相同的序列值:(1)如果此新的序列对n=N/2是偶对称,则原序列一定为圆周偶对称序列。(2)如果此新的序列对n=N/2是奇对称,则原序列一定为圆周奇对称序列。(c)共 轭 对 称(序列)和 共 轭 反 对 称(序 列)1、序 列 x(n)与y(n)=x*(-n)互 为 共共 轭轭 对对 称称.2、共共 轭轭 对对 称称 序序 列列:一个序列x(n),其满足xe(n)=x*e(-n),即称此序列为
22、共轭对称序列共轭对称序列。对 于 实 序 列 来 说,这 一 条 件 变 成 xe(n)=xe(-n),即 为 偶偶 对对 称称 序序 列列.(c)共 轭 对 称(序列)和 共 轭 反 对 称(序 列)4、共共 轭轭 反反 对对 称称 序序 列列:若一序列x(n),其满足xo(n)=-x*o(-n),称此序列为共共 轭轭 反反 对对 称称 序序 列列 对 于 实 序列 来 说 ,即 为 xo(n)=-xo(-n)奇奇 对对 称称 序序 列列.3、两序列 x(n)与y(n)若满足y(n)=-x*(-n)则互为共共 轭轭 反反 对对 称称.(d)圆 周 共 轭 对 称(序列)和 圆 周 共 轭 反
23、 对 称(序 列)1、N 点 有 限 长 序 列 x(n)与x*(-n)NRN(n)互 为 圆圆 周周 共共 轭轭 对对 称称.2、圆 周 共 轭 对 称 序 列 是 满 足 xep(n)=xep*(-n)NRN(n)即 xep(n)的 模是 圆 周 偶 对 称,辐 角是 圆 周 奇 对 称(或 说 实 部 圆 周 偶 对 称,虚 部 圆 周 奇 对 称).即把xep(n)看成分布在 N等分的圆上,在 n=0 的左半圆与右半 圆上,序列是共轭对称的。圆圆 周周 共共 轭轭 对对 称称(序列序列)的例子的例子虚部虚部实部实部实实 部部 圆圆 周周 偶偶 对对 称称,虚虚 部部 圆圆 周周 奇奇
24、对对 称称圆 周 共 轭 反 对 称(序 列)3、N 点 有 限 长 序 列 x(n)与-x*(-n)NRN(n)互 为 圆圆 周周 共共 轭轭 反反 对对 称称.4、圆、圆 周周 共共 轭轭 反对反对 称称 序序 列:列:序 列 满 足 xop(n)=-xop*(-n)NRN(n)即 xop(n)的 模是 圆 周 奇 对 称,辐 角是 圆 周 偶 对 称(或 说 实 部 圆 周 奇 对 称,虚 部 圆 周 偶 对 称).即 把 xop(n)看 成 分布 在 N 等 分 的 圆 上,在 n=0 的 左 半 圆 与 右半 圆上,序 列 是 共 轭 反 对 称 的。圆 周 共 轭 反 对 称(序
25、列)例子实实 部部 圆圆 周奇周奇 对对 称称,虚虚 部部 圆圆 周周 偶偶 对对 称称实实部部虚虚部部(2)序列的对称分量(a)奇 对 称 分 量 和 偶 对 称 分 量(b)圆 周 奇 对 称 分 量 和 圆 周 偶 对 称 分 量 (c)共 轭 对 称 分 量 和 共 轭 反 对 称 分 量(d)圆 周 共 轭 对 称 分 量 和 圆 周 共 轭 反 对 称 分 量(a)奇 对 称 分 量 和 偶 对 称 分 量说明若 x(n)为 有 限 长 序 列 且0nN-1 ,则 xo(n)与xe(n)点 数长 度 均 为(2N-1).区 别 于 奇 对 称(序列)和 偶 对 称(序列).(b)圆
26、周奇对称分量和圆周偶对称分量 1、x(n)是 长 度 N 的 有有 限限 长长 序 列,可表示成:一个 圆周奇对称序列xop(n)+一个圆周偶对称序列xep(n)即x(n)=xep(n)+xop(n).其中xop(n)称为 x(n)的圆周奇对称分量圆周奇对称分量;xep(n)称 为 x(n)的 圆周偶对称分量圆周偶对称分量.(c)共轭对称分量和共轭反对称分量 1、x(n)=共轭对称序列 xo(n)+共轭反对称序列xe(n)即x(n)=xo(n)+xe(n).其中,xo(n)又称为x(n)的 共轭反共轭反 对对 称分量称分量;xe(n)又称 为 x(n)的 共轭共轭 对对 称称 分分 量量.看出
27、xo(n)和xe(n)分 别 满 足 奇 对 称 和 偶 对 称 的 条 件,且 二 者 之 和 为 x(n)。(d)圆 周 共 轭 对 称 分 量 和 圆 周 共 轭 反 对 称 分 量 1、x(n)是长度为N的有限长序的有限长序 列列,可表示成一圆周共轭反 对称序列xop(n)+一圆周共轭对称序列xep(n).即 x(n)=xep(n)+xop(n)看 出 满 足 圆 周 奇 对 称 和 圆 周 偶 对 称 的 条 件,且 二 者 之 和 为 x(n).其中:xop(n)称为 x(n)的圆周共轭反对称分量;xep(n)称为x(n)的圆周共轭对称分量 2、x(n)是长度为N的有限长序列,可表
28、示:4.相关(1)线性相关(2)圆周相关(1)线性相关(2)圆周相关注:圆周相关结果长度不变为N。相关通信中很重要。第六节离散付里叶变换的性质一、引入在由DFS引出DFT的过程中我们知道,DFT本质上是和周期序列的DFS概念紧密相关的,因而它们在性质上有着极大的相似,并由DFT隐含周期性(对应于DFS的显式周期性)所保证。二、DFT的性质和定理分类(1)线 性(2)时 间 移 位(3)频 率 移 位(4)圆 周 卷 积 定 理(5)圆 周 相 关 定 理(6)对 称 性 质 (7)DFT形式的帕赛瓦尔定理能量计算公式 (8)DFT 的 奇,偶,虚,实 关 系三、假设条件设x1(n),x2(n)
29、都是两个长度为N的有限长序列,它们的离散付里时变换分别为四、性质(1)线性x1(n),x2(n)的线性组合有:其中a,b为任一常数,本性质可由定义直接证明。证:线性说明如果x1(n)和x2(n)长度皆为N,即0nN-1范围有值,则aX1(k)+bX2(k)的长 度也是N;若x1(n)和x2(n)长度不等长度不等,设x1(n)长度为N1,x2(n)长度为N2,则ax1(n)+bx2(n)的长度应为N=maxN1,N2,故DFT必须按长度N计算。若N1频域相乘频域卷积-时域相乘(4)圆 周 卷 积 定 理-2-说明时 域 卷 积 对 应 于 频 域 相 乘,而 时 域 相 乘 对 应 于 频 域
30、卷 积.这 与 我 们 曾 学 过 的 其 他 变 换(FT/L/Z)的 卷 积 定 理 是 相 似 的.但 注 意,由 于 DFT 隐 含 的 周 期 性,卷 积 必 须 是 圆圆 周周 卷卷 积积 才 有 此 性 质.注 意 第 二 个 关 系 中 的 系 数,不 要 忽 略。(4)圆 周 卷 积 定 理-3线卷积和圆卷积步骤比较线卷积:反折、平移、相乘、积分(或相加)圆卷积:周期化周期化、反折、平移、相乘、相加(5)圆 周 相 关 定 理有限长序列的相关运算可分为圆相关(循环相关)与线相关两种形式通常,可借助于圆相关求线相关。(复习)卷积离散线卷积:离散圆卷积:离散线相关:离散圆相关:卷
31、积与相关不同:y是共轭且y中为m-n,卷积与次序无关而相关与次序有关。作业P133页第5,6,12题参看程佩青的光盘中第三章的离散付里叶图形的测验第3题(是本书第134页第11题。)五、对 称 性 质DFT 的 对称性质较为复杂,归为以下三类:1.共轭与圆周共轭对称 在 时 频 域 的 对 应 关 系;2.实(虚)部 与 圆 周 共 轭 对 称(反 对 称)分 量 在 时,频 域 的 对 应 关 系;3.时 域 为 实 序 列 时 对 应 DFT 特 征;另 外,在 以 上 对 称 性 质 的 基 础 上,可 归 纳 总 结 出 x(n)与 X(k)的 奇,偶,虚,实 关 系,利 用 这 些
32、关 系,可 减 少 计 算 DFT 时 的 运 算 量。1.共轭与圆周共轭对称在时频域 的 对 应 关 系设 x(n)为 N 点 有 限 长 序 列,0nN-1则有:如下关系1,关系2 和关系3.(1)关系1时 域 x(n)取 共 轭,对 应 于 频 域 X(k)取 圆 周 共 轭 对 称.若 x(n)本 身 是 实 序 列,对 应 于 频 域 X(k)就 是 圆 周 共 轭 对 称 序 列;反 之 亦 然.原序列序列共轭原序列频域圆周共轭原序列为实序列,其频域为圆周共轭对称序列证明(2)关系2时 域 x(n)取 圆圆 周周 偶偶 对对 称称,对 应 于 频 域 X(k)也 取 圆圆 周周 偶
33、偶 对对 称称.若 x(n)本 身 是 圆 周 偶 对 称 序 列,对 应 频 域 X(k)也 是 圆 周 偶 对 称 序 列;反 之 亦 然.证明解释:设有限长N序列为y(n)=ye(n)+yo(n)已知时 域 x(n)取圆周偶对称取圆周偶对称,则有:对 应 于 频 域 X(k)也 取 圆圆 周周 偶偶 对对 称称.如果y(n)是圆周偶对称序列,即只有ye(n)分量,则Y(k)当然也是圆周偶对称序列。原序列序列时域圆周偶对称序列频域圆周偶对称(3)关系3此 关 系 与 关 系1成 对 偶关系.频 域 X(k)取 共 轭,对 应 于 时 域 x(n)取 圆 周 共 轭 对 称.若 X(k)是
34、实 序 列,则 对 应 时 域 x(n)是 圆 周 共 轭 对 称 序 列;反 之 亦 然.2.实(虚)部与圆周共轭对称(反对称)分量 在 时 频 域 的 对 应关系设x(n)为N点有限长序列0nN-1则有关系1,关系2,关系3:关系1时 域 x(n)取 实 部,对 应 频 域 取 X(k)的 圆 周 共 轭 对 称 分 量.若 x(n)本 身 是 实 序 列,那 么 由 于因 而 对 应 频 域 X(k)是 圆 周 共 轭 对 称 序 列;反 之 亦 然.关系2时 域 x(n)取 虚 部 并 加 权 j,对 应 频 域 取 X(k)的 圆 周 共 轭 反 对 称 分 量.若x(n)本身是纯虚
35、序列,那么X(k)关系3说明:(1)对时域x(n)取圆周共轭对称分量(xep(n),即对频域X(k)取实部;对时域x(n)取圆周共轭反对称分量(xop(n),即对频域X(k)取虚部加权j;若X(k)本身是实序列,则时域x(n)是圆周共轭对称序列;若X(k)本身是纯虚序列,则时域x(n)是圆周共轭反对称序列;反之亦然。3.时域是实序列时对应DFT特征设x(n)为长度为N的有限长实序列,0 nN-1,DFTx(n)=X(k)有以下几个特征:(5个)(1)特征1x(n)为实序列,显然X(k)只有圆周共轭对称分量即X(k)=X*(N-k)NRN(k)说明:(1)x(n)的DFT即X(k)是圆周共轭对称
36、序列。(2)是实(虚)部与圆周共轭对称(反对称)分量在时域、频域的对应关系。(2)特征2ReX(k)=ReX(N-k)NRN(k)说明:X(k)的实部是圆周偶对称序列。(3)特征3ImX(k)=-ImX(N-k)NRN(k)说明:X(k)的虚部是圆周奇对称序列。(4)特征4|X(k)|=|X(N-k)N|RN(k)说明:X(k)的模是圆周偶对称序列。(5)特征5argX(k)=-argX(N-k)NRN(k)说明:X(k)的相角是圆周奇对称序列。4.序列及其DFT的奇偶虚实关系由上对称性质基础上,可归纳总结出x(n)与X(k)的奇、偶;虚、实关系,利用这些关系,可以减少计算DFT的运算量。下面
37、总结归纳出有限长序列及其DFT的奇、偶;虚、实关系。这一关系清晰地展示了时域序列的奇、偶;虚、实特性与频域序列的奇、偶;虚、实特性是如何对应的。(1)奇、偶;虚、实的含义所 谓 奇,偶,虚,实 的 含 义 如 下:奇 -指 序 列 是 圆 周 奇 对 称 序 列 偶 -指 序 列 是 圆 周 偶 对 称 序 列虚 -指 序 列 是 纯 虚 序 列实 -指 序 列 是 实 序 列(2)奇偶虚实关系表六、DFT形式下的帕塞瓦尔定理(Parsevals Theorem)说明:(1)这是DFT形式下的帕塞瓦尔定理(Parsevals,Theorem)(2)只需令y(n)=x(n),再两边取模,便得到明
38、确物理意义的能量计算公式。证明Parseval定理七、DFT性质一览表1七、DFT性质一览表2作业3第133页,第15,16,17,18,19题参看程佩青的光盘中第三章的有关DFT性质的测验题。测验:离散付里叶变换性质(对称性)1.时域序列取共轭共轭,则频域序列有什么关系?2.时域序列取圆周共轭反折圆周共轭反折,则频域序列有什么关系?3.时域序列取实部实部,则频域序列有什么关系?4.时域序列取虚部虚部,则频域序列有什么关系?5.时域序列取圆周共轭对称分量圆周共轭对称分量,则频域序列有什么关系?6.时域序列取圆周共轭反对称分量圆周共轭反对称分量,则频域序列有什么关系?7.时域序列取圆周反折圆周反
39、折,则频域序列有什么关系?第七节抽样z变换频率抽样理论复习:时域抽样定理奈奎斯特抽样定理:要想抽样后能够不失真的还原出原信号,则抽样频率必须大于两倍信号谱的最高频率。或抽样内插公式即由信号的抽样值xa(mT)经此公式而得到连续信号xa(t).主要内容阐明:(1)z变换与DFT的关系(抽样z变换),并引出抽样z变换的概念,讨论频域抽样不失真条件。(2)频域抽样理论(频域抽样不失真条件)(3)频域内插公式一、z变换与DFT关系(1)引入FT引出DFT定义式。DFT看作是DTFT在 频 域 抽 样 后 的 变 换 对.在Z变换与L变换中,可了解到DTFT是单位圆上的Z 变 换.所以对DTFT进行频域
40、抽样时,自 然可以看作是对单位圆上的 Z变换进行抽样.(2)推 导w是 单 位 圆 上 各 点 的 数 字 角 频 率.Z 变 换 的 定 义 式(正 变 换):取z=ejw 代 入,得 到 单 位 圆 上 Z 变 换 为则这正是离散傅里叶变换(DFT)正变换定义式.再 抽 样-N 等 分抽样间隔w=2k/N,即w值为0,2/N,4/N,考虑x(n)是N点有限长序列,n只需0N-1即可。将w=2k/N代入并改变上下限,得(3)结论1 看 出:有 限 长 序 列 x(n)的 DFT的 X(k)序 列 的 各 点 值 等于 x(n)Z 变 换 后 在 单 位 圆 上 N 等 分 抽 样 的 各 点
41、 处 所 得 的 Z 变 换 值,即 这 就 是 Z 变 换 与 DFT 的 关 系.(4)结论2有限长序列补零加长 N增加,求其DFT。发现频 谱包络不变,只是抽样点更密.原因:即N补零加长并不改变有限长序列本身,因而其 Z变换不变,而只是增加了N值。根 据 每个 X(k)仍 等 于X(ejw)这 一 包 络.由于0kN-1,X(k)值的个数增加了,谱线变密.二、频率抽样理论(频域抽样不失真条件)(1)问题引入1、是否任何一序列(或说任何一个频率特性)都能用频域抽样的办法去逼近呢?2、其 限 制 条 件 是 什 么?(2)分析即 频 域 按 每 周 期 N 点 抽 样,时 域 便 按 N 点
42、 周 期 延 拓.(3)结论长度为M的有限长序列,频域抽样不失真的条件:频域抽样点数N要大于或等于序列长度M,即满足NM.此时可得到表明长度为N(或小于N)的有限长序列可用它的z变换在单位圆上的N个均分点上的抽样值精确地表示.(4)抽样后序列能否无失真恢复原时域信号有限长序列非非有限长序列1.长度为M2.当频域抽样不够密,即NM或=M时,可利用其z变换在单位圆上的N个均分点上的抽样值精确地表示。1.时域周期延拓后必然必然造成混叠现象,因而一定会产生误差。2.当n增加时信号衰减得越快,或抽样越密(即抽样点数N越多)则误差越小。(5)注意点 DFT 变 换 对 的 一 一 对 应 关 系 也 是
43、由 此 而 得 到 保 证 的.实 际 上,在 我 们 从 连 续 傅 里 叶 变 换 引 出 DFT 时,也 只 有 按 此 条 件 对 频 域 进 行 抽 样,才 能 在 最 后 正 确 导 出 DFT 变 换 对 定 义式.(6)例子频域抽样:看一个矩形序列,频域抽样是指对时域已是离散,频域仍是连续信号。现在频域上进行抽样处理,使其频域也离散化。解:频域抽样,按N=5点,频域抽样,时域延拓相加,时域延拓的周期个数等于频域的抽样点数N=5,由于N=M,所以时域延拓恰好无混叠现象。按N=4时进行抽样,由于N=4,而序列长度为M=5,NM,时域延拓后产生混叠现象。(原信号为红色,延拓取主值区间
44、后的恢复信号为兰色。)三、频域内插公式从 频 域 抽 样 不 失 真 条 件 可 以 知 道:N 个 频 域 抽 样 X(k)能 不 失 真 的 还 原 出 长 度 为 N 的 有 限 长 序 列 x(n)。那 么 用 N 个 X(k)也 一 定 能 完 整 地 表 示 出 X(z)以 及 频 率 响 应 即 单 位 圆 上 的 X(z).过 程 很 简 单,先 把 N 个 X(k)作 IDFT 得 到 x(n),再 把 x(n)作 Z 变 换 便 得 到 X(z).(1)内插公式(2)内插函数推导频域采样X(k)表示X(z)的内插公式和内插函数(3)频域响应的内插公式(4)频域响应的内插函数推导频域响应的内插函数(5)说明从 公 式 中看 出:在 每 个 抽 样 点 上X(ejw)就 精 确 地 等 于 X(k)(因 为 其 他 的 内 插 函 数 在 这 一 点 上 的 值 为 零,无 影 响),即各 抽 样 点 之 间 的X(ejw)值,则 由 各 抽 样 点 的 加 权 内 插 函 数 在 所 求 点 上 的 值 的 叠 加 而 得到.频 率 响 应 的 内 插 函 数 具 有 线 性 相 位.