数字信号处理第2章-课件.ppt

1、第2章 离散傅里叶变换 第2章 离散傅里叶变换 2.1 引言引言 2.2 周期序列的离散傅里叶级数（周期序列的离散傅里叶级数（DFS）2.3 离散傅里叶级数（离散傅里叶级数（DFS）的性质）的性质 2.4 有限长序列离散傅里叶变换（有限长序列离散傅里叶变换（DFT）2.5 离散傅里叶变换的性质离散傅里叶变换的性质 2.6 频域采样理论频域采样理论 第2章 离散傅里叶变换 2.1 引引 言言 在第1章中讨论了序列的傅里叶变换和Z变换。由于数字计算机只能计算有限长离散序列，因此有限长序列在数字信号处理中就显得很重要，当然可以用Z变换和傅里叶变换来研究它，但是，这两种变换无法直接利用计算机进行数值计

2、算。针对序列“有限长”这一特点，可以导出一种更有用的变换：离散傅里叶变换（Discrete Fourier Transform，简写为DFT）。它本身也是有限长序列。第2章 离散傅里叶变换 作为有限长序列的一种傅里叶表示法，离散傅里叶变换除了在理论上相当重要之外，而且由于存在有效的快速算法快速离散傅里叶变换，因而在各种数字信号处理的算法中起着核心作用。有限长序列的离散傅里叶变换（DFT）和周期序列的离散傅里叶级数（DFS）本质上是一样的。为了讨论离散傅里叶级数与离散傅里叶变换，我们首先来回顾并讨论傅里叶变换的几种可能形式，见图2-1所示。第2章 离散傅里叶变换 图 2-1 各种形式的傅里叶变换

3、 xa(t)txp(t)ootTpx(nT)oN点xp(n)oN点nTn(a)(b)(c)(d)|Xa(j)|10o0|Xp(jk)|ok|X(ej)|1/T|X(ejk)|sooN点sT第2章 离散傅里叶变换 一个非周期实连续时间信号xa(t)的傅里叶变换，即频谱Xa(j)是一个连续的非周期函数，这一变换对的示意图见图2-1(a)。该变换关系与第1章“连续时间信号的采样”中所涉及到的非周期连续时间信号xa(t)的情况相同。一个周期性连续时间信号xp(t)，其周期为Tp，该信号可展成傅里叶级数，其傅里叶级数的系数为 ，即xp(t)的傅里叶变换或频谱Xp(jk)是由各次谐波分量组成的，并且是非周

4、期离散频率函数，xp(t)和Xp(jk)的示意图见图2-1(b)。其中，离散频谱相邻两谱线之间的角频率间隔为=2F=2/Tp，k为谱谐波序号。)(kXp第2章 离散傅里叶变换 在第1章里讨论了一个非周期连续时间信号xa(t)经过等间隔采样的信号（x(nT)），即离散时间信号序列x(n)，其傅里叶变换X(ej)是以2为周期的连续函数，振幅特性如图2-1(c)所示。这里的是数字频率，它和模拟角频率的关系为=T。若振幅特性的频率轴用表示，则周期为s=2/T。比较图2-1(a)、（b）和(c)可发现有以下规律：如果信号频域是离散的，表现为周期性的时间函数。相反，在时域上是离散的，则该信号在频域必然表现

5、为周期性的频率函数。不难设想，一个离散周期序列，它一定具有既是周期又是离散的频谱，其振幅特性如图2-1(d)所示。第2章 离散傅里叶变换 表表2-1 四种傅里叶变换形式的归纳四种傅里叶变换形式的归纳 时间函数 频率函数 连续和非周期 非周期和连续 连续和周期 非周期和离散 离散和非周期 周期和连续 散和周期 周期和离散 第2章 离散傅里叶变换 可以得出一般的规律：一个域的离散对应另一个域的周期延拓，一个域的连续必定对应另一个域的非周期。表2-1对这四种傅里叶变换形式的特点作了简要归纳。下面我们先从周期性序列的离散傅里叶级数开始讨论，然后讨论可作为周期函数一个周期的有限长序列的离散傅里叶变换。第

6、2章 离散傅里叶变换 2.2 周期序列的离散傅里叶级数周期序列的离散傅里叶级数(DFS)设 是一个周期为N的周期序列，即)()(rNnxnx)(nx r为任意整数 周期序列不是绝对可和的，所以不能用Z变换表示，因为在任何z值下，其Z变换都不收敛，也就是 nnznx|)(|第2章 离散傅里叶变换 但是，正如连续时间周期信号可以用傅里叶级数表示一样，周期序列也可以用离散傅里叶级数来表示，该级数相当于成谐波关系的复指数序列（正弦型序列）之和。也就是说，复指数序列的频率是周期序列 的基频（2/N）的整数倍。这些复指数序列ek(n)的形式为)(nx)()(2neenerNkknNjk（2-1）式中,k,

7、r为整数。第2章 离散傅里叶变换 由式（2-1）可见，复指数序列ek(n)对k呈现周期性，周期也为N。也就是说,离散傅里叶级数的谐波成分只有N个独立量，这是和连续傅里叶级数的不同之处（后者有无穷多个谐波成分），因而对离散傅里叶级数，只能取k=0 到N-1的N个独立谐波分量，不然就会产生二义性。因而 可展成如下的离散傅里叶级数，即)(nx102)(1)(NkknNjekXNnx（2-2）式中，求和号前所乘的系数1/N是习惯上已经采用的常数，是k次谐波的系数。)(kX第2章 离散傅里叶变换 下面我们来求解系数 ，这要利用复正弦序列的正交特性，即)(kX01111122102rNjrNNjNnrnN

8、jeeNeNr=mN,m为整数 其他r（2-3）将式（2-2）两端同乘以 ，然后从n=0 到N-1的一个周期内求和，则得到 rnNje2第2章 离散傅里叶变换)(1)()(1)(10)(210)(21010102rXeNkXekXNenxNnnrkNjNknrkNjNkNnNnrnNj把r换成k可得 102)()(NnknNjenxkX（2-4）这就是求k=0 到N-1的N个谐波系数的公式。同时看出 也是一个以N为周期的周期序列，即)(kX)(kX第2章 离散傅里叶变换)()()()(21010)(2kXenxenxmNkXknNjNnNnnmNkNj 这和离散傅里叶级数只有N个不同的系数 的

9、说法是一致的。可以看出，时域周期序列 的离散傅里叶级数在频域（即其系数 也是一个周期序列。因而 与 是频域与时域的一个周期序列对，式（2-2）与式（2-4）一起可看作是一对相互表达周期序列的离散傅里叶级数（DFS）对。为了表示方便，常常利用复数量WN来写这两个式子。WN定义为)(kX)(nx)(kX)(kX)(nxNjNeW2（2-5）第2章 离散傅里叶变换 使用WN，式（2-4）及式（2-2）可表示为:nkNNknkNjNknkNNnnkNjNnWkXNekXNkXIDFSnxWnxenxnxDFSkX)(1)(1)()()()()()(1021010210（2-6）（2-7）式中，DFS表

10、示离散傅里叶级数正变换，IDFS表示离散傅里叶级数反变换。从上面看出，只要知道周期序列一个周期的内容，其他的内容也都知道了。所以，这种无限长序列实际上只有一个周期中的N个序列值有信息。因而周期序列和有限长序列有着本质的联系。第2章 离散傅里叶变换 例例2-1 设 为周期脉冲串)(nx)()(rNnnxr（2-8）因为对于0nN-1，,所以利用式（2-6）求出 的DFS系数为)()(nnx)(nx1)()()(1010nkNNnnkNNnWnWnxkX（2-9）在这种情况下，对于所有的k值 均相同。于是，将式（2-9）代入式（2-7）可以得出表示式)(kX1021011)()(NknkNjnkN

11、NkreNWNrNnnx（2-10）第2章 离散傅里叶变换 例例2-2 已知周期序列 如图2-2所示，其周期N=10,试求解它的傅里叶级数系数 。)(kX)(kX图2-2 例2-2的周期序列 （周期N=10）)(nx100 1 2 3 4 5 6 7 8 910n)(nx第2章 离散傅里叶变换 由式（2-6）40102101100)()(nnkjnkneWnxkX（2-11）这一有限求和有闭合形式 40102101100)()(nnkjnkneWnxkX（2-12）图 2-3 图2-2所示序列的傅里叶级数系数 的幅值)(kX10 1 2 3 4 5 6 7 8 9101520k|)(|kx5第

12、2章 离散傅里叶变换 式（2-6）中的周期序列 可看成是对 的第一个周期x(n)作Z变换，然后将Z变换在Z平面单位圆上按等间隔角2/N采样而得到的。令 )(kX)(nx0)()()()(nxnRnxnxN0nN-1 其他n 通常称x(n)为 的主值区序列，则x(n)的Z变换为)(nx10)()()(NnnnnznxznxzX（2-13）第2章 离散傅里叶变换 把式（2-13）与式（2-6）比较可知 kKGNjkNeWzzXkX4*2)()(（2-14）可以看出，当0kN-1 时，是对X(z)在Z平面单位圆上的N点等间隔采样，在此区间之外随着k的变化，的值呈周期变化。图2-4画出了这些特点。)(

13、kX)(kX第2章 离散傅里叶变换 jImz234567(N1)k02/NRezo图2-4|z|11第2章 离散傅里叶变换 由于单位圆上的Z变换即为序列的傅里叶变换，所以周期序列 也可以解释为的一个周期x(n)的傅里叶变换的等间隔采样。因为)(kX)(nx1010)()()(NnnjNnnjjenxenxeX（2-15）比较式（2-15）和式（2-6），可以看出 NkjeXkX/2)()(这相当于以2/N的频率间隔对傅里叶变换进行采样。（2-16）第2章 离散傅里叶变换 例例2-3 为了举例说明傅里叶级数系数 和周期信号 的一个周期的傅里叶变换之间的关系，我们再次研究图2-2所示的序列 。在序

14、列 的一个周期中：)(kX)(nx)(nx)(nx01)(nx0n4 其他（2-17）则 的一个周期的傅里叶变换是)(nx)2/5sin()2/5sin(11)(2405jnjjnjjeeeeeX)10/sin()10/5sin()()(10410/2kkeeXkXkjkj（2-18）可以证明，若将=2k/10 代入式（2-18），即 第2章 离散傅里叶变换 图 2-5 对图2-2所示序列的一个周期作傅里叶变换的幅值|X(ej)|5o234第2章 离散傅里叶变换 图 2-6 图2-3和图2-5的重叠图，它表明一个周期序列的DFS系数等于主值区序列的傅里叶变换的采样|X(ej)|,|X(k)|o

15、2341020k50第2章 离散傅里叶变换 2.3 离散傅里叶级数（离散傅里叶级数（DFS）的性质）的性质 由于可以用采样变换来解释DFS，因此它的许多性质与变换性质非常相似。但是，由于 和 两者都具有周期性，这就使它与Z变换性质还有一些重要差别。此外，DFS在时域和频域之间具有严格的对偶关系，这是序列的Z变换表示所不具有的。设 和皆是周期为N的周期序列，们各自的DFS分别为：)(nx)(kX)(1nx)(2nx)()()()(2211nxDFSkXnxDFSkX第2章 离散傅里叶变换 2.3.1 线性线性)()()()(2121kXbkXanxbnxaDFS（2-19）式中，a和b为任意常数

16、，所得到的频域序列也是周期序列，周期为N。这一性质可由DFS定义直接证明，留给读者自己去做。第2章 离散傅里叶变换 2.3.2 序列的移位序列的移位)()()()()(2lkXnxWDFSkXekXWmnxDFSnlNmkNjmkN（2-20）（2-21a）或)()()(2nxenxWlkXIDFSnlNjnlN证（2-21b）mkNkiNmNminkNNnWWixWmnxmnxDFS110)()()(i=n+m 第2章 离散傅里叶变换 由于 都是以N为周期的周期函数，故 kiNWix及)()()()(10kXWWixWmnxDFSmkNkiNNimkN 由于 与 的对称特点，可以用相似的方法

17、证明式（2-21a）：)(nx)(kX)()()()()(1010lkXWnxWnxWnxWDFSnklNNnknNNinlNnlN第2章 离散傅里叶变换 2.3.3 周期卷积周期卷积 如果)()()(21kXkXkY则）（mnxmxkYIDFSnyNm2101)()()(或）（mnxmxnyNm1102)()(证证 knNNkWkXkXNkXkXIDFSny）（210121)(1)()()(代入 mkNNmWmxnX1011)()(2-22)第2章 离散傅里叶变换 得)()(1)()(1)(2101)(210101)(210101mnxmxWkXNmxWkXmxNnyNmkmnNNkNmkm

18、nNNkNm）（）（将变量进行简单换元，即可得等价的表示式）（mnxmxnyNm1102)()(第2章 离散傅里叶变换 式（2-22）是一个卷积公式，但是它与非周期序列的线性卷积不同。首先，和（或 和 都是变量m的周期序列，周期为N，故乘积也是周期为N的周期序列;其次，求和只在一个周期上进行，即m=0到N-1，所以称为周期卷积。)(1mx)(2mnx)(2mx)(1mnx第2章 离散傅里叶变换 周期卷积的过程可以用图2-7来说明，这是一个N=7 的周期卷积。每一个周期里 有一个宽度为4的矩形脉冲，有一个宽度为3的矩形脉冲，图中画出了对应于n=0,1,2 时的 。周期卷积过程中一个周期的某一序列

19、值移出计算区间时，相邻的同一位置的序列值就移入计算区间。运算在m=0到N-1区间内进行，即在一个周期内将与 逐点相乘后求和，先计算出n=0,1,N-1的结果，然后将所得结果周期延拓，就得到所求的整个周期序列 。)(1nx)(2nx)(2mnx)(2mnx)(1mx)(ny第2章 离散傅里叶变换 图 2-7 两个周期序列（N=7）的周期卷积 n0n(a)(c)(1mxm)(1nx NN)(2nx10 NN0N1(d)第2章 离散傅里叶变换 图 2-7 两个周期序列（N=7）的周期卷积(d)(e)mn 1N10)1(2mx)0(2mxN10mn 0第2章 离散傅里叶变换 图 2-7 两个周期序列（

20、N=7）的周期卷积(f)(g)N0Nn)(ny1123 320)2(2mxmn 2N1第2章 离散傅里叶变换 由于DFS和IDFS变换的对称性，可以证明（请读者自己证明）时域周期序列的乘积对应着频域周期序列的周期卷积。即，如果)()()(21nxnxny则)()(1)()(1)()()(1102210110lkXlXNlkXlXNWnynyDFSkYNlNlNnnkN(2-23)第2章 离散傅里叶变换 2.4 有限长序列离散傅里叶变换（有限长序列离散傅里叶变换（DFT）2.4.1 DFT的定义的定义 上一节我们讨论的周期序列实际上只有有限个序列值有意义，因而它和有限长序列有着本质的联系。本节将

21、根据周期序列和有限长序列之间的关系，由周期序列的离散傅里叶级数表示式推导得到有限长序列的离散频域表示即离散傅里叶变换（DFT）。设x(n)为有限长序列，长度为N，即x(n)只在n=0到N-1点上有值，其他n时，x(n)=0。即 nNnnxnx其他010)()(第2章 离散傅里叶变换 为了引用周期序列的概念，我们把它看成周期为N的周期序列 的一个周期，而把 看成x(n)的以N为周期的周期延拓，即表示成:)(nx)(nxnNnnxnx其他010)()(rrNnxnx)()(这个关系可以用图2-8来表明。通常把 的第一个周期n=0 到n=N-1 定义为“主值区间”，故x(n)是 的“主值序列”，即主

22、值区间上的序列。而称为x(n)的周期延拓。对不同r值x(n+rN)之间彼此并不重叠，故上式可写成)(nx)(nx)(nxNnxNnxnx)()mod()(2-26)第2章 离散傅里叶变换 图2-80N1n)(nx N0N1n主值区间x(n)第2章 离散傅里叶变换 用(n)N表示（n mod N），其数学上就是表示“n对N取余数”，或称“n对N取模值”。令 mNnn10n1N-1,m为整数 则n1为n对N的余数。例如，是周期为N=9的序列，则有：)(nx)8()1()1()4()22()22()4()13()13()8()8()8(9999xxxxxxxxxxxx第2章 离散傅里叶变换 利用前面

23、的矩形序列RN(n)，式（2-24）可写成)()()(nRnxnxN(2-27)同理，频域的周期序列 也可看成是对有限长序列X(k)的周期延拓，而有限长序列X(k)可看成是周期序列 的主值序列，即:)(kX)(kX)()()()()(kRkXkXkXkXNN(2-28)(2-29)我们再看表达DFS与IDFS的式(2-6)和式(2-7)：1010)(1)()()()()(NknkNNnnkNWkXNkXIDFSnxWnxnxDFSkX第2章 离散傅里叶变换 这两个公式的求和都只限定在n=0到N-1和k=0 到N-1 的主值区间进行，它们完全适用于主值序列x(n)与X(k)，因而我们可以得到有限

24、长序列的离散傅里叶变换的定义:1010)(1)()()()()(NknkNNnnkNWkXNkXIDFTnxWnxnxDFTkX0kN-1 0nN-1（2-30）（2-31）第2章 离散傅里叶变换 x(n)和X(k)是一个有限长序列的离散傅里叶变换对。我们称式(2-30)为x(n)的N点离散傅里叶变换(DFT),称式(2-31)为X(k)的N点离散傅里叶反变换(IDFT)。已知其中的一个序列，就能惟一地确定另一个序列。这是因为x(n)与X(k)都是点数为N的序列，都有N个独立值（可以是复数），所以信息当然等量。此外，值得强调得是，在使用离散傅里叶变换时，必须注意所处理的有限长序列都是作为周期序

25、列的一个周期来表示的。换句话说，离散傅里叶变换隐含着周期性。第2章 离散傅里叶变换 例例2-4 已知序列x(n)=(n)，求它的N点DFT。解解 单位脉冲序列的DFT很容易由DFT的定义式（2-30）得到：1001)()(NnNnkNWWnkX k=0,1,N-1 (n)的X(k)如图2-9。这是一个很特殊的例子，它表明对序列(n)来说，不论对它进行多少点的DFT，所得结果都是一个离散矩形序列。第2章 离散傅里叶变换 图2-9 序列(n)及其离散傅里叶变换 10n(n)X(k)10 12N1k第2章 离散傅里叶变换 例例 2-5 已知x(n)=cos(n/6)是一个长度N=12的有限长序列，求

26、它的N点DFT。解解 由DFT的定义式（2-30）110)1(122110)1(122110122661211021216cos)(nknjnknjnnkjnjnjnkneeeeeWnkX 利用复正弦序列的正交特性（2-3）式，再考虑到k的取值区间，可得 11,0,011,16)(kkkkX其他第2章 离散傅里叶变换 图 2-10 有限长序列及其DFT0 1 211x(n)n01X(k)11n第2章 离散傅里叶变换 例例 2-6 已知如下X(k):13)(kXk=0 1k9 求其10点IDFT。解解 X(k)可以表示为 X(k)=1+2(k)0k9 写成这种形式后，就可以很容易确定离散傅里叶反

27、变换。由于一个单位脉冲序列的DFT为常数：1)()()()(111nxDFTkXnnx第2章 离散傅里叶变换 2.4.2 DFT与序列傅里叶变换、与序列傅里叶变换、Z变换的关系变换的关系 若x(n)是一个有限长序列，长度为N，对x(n)进行Z变换 10)()(NnnznxzX比较Z变换与DFT，我们看到，当z=W-kN时)()()(10nxDFTWnxzXNnnkNWzkN即 kNWzzXkX)()(（2-32）第2章 离散傅里叶变换 表明 是Z平面单位圆上幅角为 的点，也即将Z平面单位圆N等分后的第k点，所以X(k)也就是对X(z)在Z平面单位圆上N点等间隔采样值，如图2-11所示。此外，由

28、于序列的傅里叶变换X(ej)即是单位圆上的Z变换，根据式（2-32），DFT与序列傅里叶变换的关系为 kNjkNeWz2kNWkN2NeXeXkXNjkkNjN2)()()(2（2-33）（2-34）第2章 离散傅里叶变换 式（2-33）说明X(k)也可以看作序列x(n)的傅里叶变换X(ej)在区间0,2上的N点等间隔采样，其采样间隔为N=2/N，这就是DFT的物理意义。显而易见，DFT的变换区间长度N不同，表示对X(ej)在区间0,2上的采样间隔和采样点数不同，所以DFT的变换结果也不同。图 2-11 DFT与序列傅里叶变换、Z变换的关系 jIm(z)o2NW1NW0NWk0)2(NNW)3

29、(NNWRezoX(ej)X(k)第2章 离散傅里叶变换 例例 2-7 有限长序列x(n)为 01)(nx0n4 其余n 求其N=5 点离散傅里叶变换X(k)。解解 序列x(n)如图2-12（a）所示。在确定DFT时，我们可以将x(n)看作是一个长度N5的任意有限长序列。首先我们以N=5 为周期将x(n)延拓成周期序列 ，如图2-12（b），的DFS与x(n)的DFT相对应。因为在图2-12（b）中的序列在区间0nN-1 上为常数值，所以可以得出)(nx011)()/2210)/2(NeeekXNkjkjNnnNkj（k=0,N,2N,其他 第2章 离散傅里叶变换 也就是说，只有在k=0 和k

30、=N 的整数倍处才有非零的DFS系数 值。这些DFS系数如图2-12（c）所示。为了说明傅里叶级数 与x(n)的频谱X(ej)间的关系，在图2-12（c）中也画出了傅里叶变换的幅值|X（ej)|。显然，就是X(ej)在频率k=2k/N 处的样本序列。按照式（2-29），x(n)的DFT对应于取 的一个周期而得到的有限长序列X(k)。这样，x(n)的5点DFT如图2-12（d）所示。)(kx)(kX)(kX)(kX0511)()(52215052kjkjnnkjeeenxkX k=0,1,2,3,4 k=0 k=0,1,2,3,4 第2章 离散傅里叶变换 x(n)(nx(a)(b)(c)40nn

31、(d)(kX50101234567891011k|X(ej)|5X(k)k01234O24图 2-12 DFT的举例说明（a）有限长序列x(n);（b）由x(n)形成的周期N=5的周期序列;（c）对应于 的傅里叶级数 和x(n)的傅里叶变换的幅度特性|X(ej)|;（d）x(n)的DFT X(k)(nx)(kX第2章 离散傅里叶变换 图 2-13 DFT的举例说明（a）有限长序列x(n);（b）由x(n)形成的周期N=10的周期序列 ;（c）DFT的幅值 )(nxx(n)104nn1)(nx0410100101053.241.2411.243.24k|X(k)|(a)(b)(c)第2章 离散傅

32、里叶变换 同样，一个常数的DFT是一个单位脉冲序列：1)(2nxX2(k)=DFTx2(n)=N(k)所以)(51)(nnx第2章 离散傅里叶变换 通过式（2-26）和式（2-27）联系起来的有限长序列x(n)和周期序列 之间的差别似乎很小，因为利用这两个关系式可以直接从一个构造出另一个。然而在研究DFT的性质以及改变x(n)对X(k)的影响时，这种差别是很重要的。信号时域采样理论实现了信号时域的离散化，使我们能用数字技术在时域对信号进行处理。而离散傅里叶变换理论实现了频域离散化，因而开辟了用数字技术在频域处理信号的新途径，从而推进了信号的频谱分析技术向更深更广的领域发展。)(nx第2章 离散

33、傅里叶变换 2.5 离散傅里叶变换的性质离散傅里叶变换的性质 本节讨论DFT的一些性质，它们本质上和周期序列的DFS概念有关，而且是由有限长序列及其DFT表示式隐含的周期性得出的。以下讨论的序列都是N点有限长序列，用DFT表示N点DFT，且设:DFTx1(n)=X1(k)DFTx2(n)=X2(k)第2章 离散傅里叶变换 2.5.1 线性线性)()()()(2121kbXkaXnbxnaxDFT式中，a,b为任意常数。该式可根据DFT定义证明。第2章 离散傅里叶变换 2.5.2 圆周移位圆周移位 1.定义定义 一个长度为N的有限长序列x(n)的圆周移位定义为 y(n)=x(n+m)NRN(n)

34、（2-37）我们可以这样来理解上式所表达的圆周移位的含义。首先，将x(n)以N为周期进行周期延拓得到周期序列 ;再将 加以移位：Nnxnx)()()(nx)()(mnxmnxN（2-38）第2章 离散傅里叶变换 然后，再对移位的周期序列 取主值区间（n=0 到N-1）上的序列值，即x(n+m)NRN(n)。所以，一个有限长序列x(n)的圆周移位序列y(n)仍然是一个长度为N的有限长序列，这一过程可用图2-14（a）、（b）、（c）、（d）来表达。从图上可以看出，由于是周期序列的移位，当我们只观察 0nN-1 这一主值区间时，某一采样从该区间的一端移出时，与其相同值的采样又从该区间的另一端循环移

35、进。因而，可以想象x(n)是排列在一个N等分的圆周上，序列x(n)的圆周移位，就相当于x(n)在此圆周上旋转，如图2-14（e）、（f）、(g)所示，因而称为圆周移位。若将x(n)向左圆周移位时，此圆是顺时针旋转;将x(n)向右圆周移位时，此圆是逆时针旋转。此外，如果围绕圆周观察几圈，那么看到的就是周期序列 。)(mnx)(nx第2章 离散傅里叶变换 图 2-14 圆周移位过程示意图(e)x(n)21n 0N1N2on 0N1N221n 0N2N1(f)(g)210 x(n)n0n)(nxNnxnx)2()2(0n)()2(nRnxNN0N1n(a)(b)(c)(d)N1N1N1第2章 离散傅

36、里叶变换 2.时域圆周移位定理时域圆周移位定理设x(n)是长度为N的有限长序列，y(n)为x(n)圆周移位，即)()()(nRmnxnyNN则圆周移位后的DFT为)()()()()(kXWnRmnxDFTnyDFTkYmkNNN证证 利用周期序列的移位性质加以证明。)()()(kXWmnxDFSmnxDFSmkNN第2章 离散傅里叶变换 再利用DFS和DFT关系)()()()()()()(kXWkRkXWnRmnxDFTnRmnxDFTmkNNmkNNNN这表明，有限长序列的圆周移位在离散频域中引入一个和频率成正比的线性相移 ，而对频谱的幅度没有影响。mkNjknNeW2第2章 离散傅里叶变换

37、 3.频域圆周移位定理频域圆周移位定理 对于频域有限长序列X(k)，也可看成是分布在一个N等分的圆周上，所以对于X(k)的圆周移位，利用频域与时域的对偶关系，可以证明以下性质：若)()(nxDFTkX则)()()()(2nxenxWkRlkXIDFTnlNjnlNNN这就是调制特性。它说明，时域序列的调制等效于频域的圆周移位。第2章 离散傅里叶变换 2.5.3 圆周卷积圆周卷积设x1(n)和x2(n)都是点数为N的有限长序列（0nN-1），且有：)()()()(2211kXnxDFTkXnxDFT若)()()(21kXkXkY则 10121021)()()()()()()()(NmNNNmNN

38、nRmnxmxnRmnxmxkYIDFTny第2章 离散傅里叶变换 一般称式(2-41)所表示的运算为x1(n)和x2(n)的N点圆周卷积。下面先证明式(2-41)，再说明其计算方法。证证 这个卷积相当于周期序列 和 作周期卷积后再取其主值序列。先将Y(k)周期延拓，即)(1nx)(2nx)()()(21kXkXkY根据DFS的周期卷积公式 NNmNNmmnxmxmnxmxny)()()()()(10212101第2章 离散傅里叶变换 由于0mN-1 为主值区间，因此)()(11mxmxN1021)()()()()()(NmNNNnRmnxmxnRnyny 将 式经过简单换元，也可证明)(ny

39、1012)()()()(NmNNnRmnxmxny第2章 离散傅里叶变换 卷积过程可以用图2-15来表示。圆周卷积过程中，求和变量为m,n为参变量。先将x2(m)周期化，形成x2(m)N，再反转形成x2(-m)N，取主值序列则得到x2(-m)NRN(m)，通常称之为x2(m)的圆周反转。对x2(m)的圆周反转序列圆周右移n，形成x2(n-m)NRN(m)，当n=0,1,2,N-1时，分别将x1(m)与x2(n-m)NRN(m)相乘，并在m=0 到N-1 区间内求和，便得到圆周卷积y(n)。可以看出,它和周期卷积过程是一样的，只不过这里要取主值序列。特别要注意，两个长度小于等于N的序列的N点圆周

40、卷积长度仍为N，这与一般的线性卷积不同。圆周卷积用符号来表示。圆周内的N表示所作的是N点圆周卷积。N第2章 离散傅里叶变换 图图 2-15 圆周卷积过程示意图圆周卷积过程示意图 x1(n)1N1nx2(n)1N1nx2(0 m)NRN(m)1N1mooo第2章 离散傅里叶变换 图图 2-15 圆周卷积过程示意图圆周卷积过程示意图 x2(1 m)NRN(m)1N1mx2(2 m)NRN(m)1N1my(n)x1(n)x2(n)233211N1noooN第2章 离散傅里叶变换)()()()()()(210121nRmnxmxnxnxnyNNNmN或)()()()()()(110212nRmnxmx

41、nxnxnyNNNmN第2章 离散傅里叶变换 N 利用时域与频域的对称性，可以证明频域圆周卷积定理（请读者自己证明）：若)()()(21nxnxnyx1(n)，x2(n)皆为N点有限长序列，则)()(1)()()(1)()()(1)()(2111022101kXkXNkRlkXlXNkRlkXlXNnyDFTkYNNNlNNNl即时域序列相乘，乘积的DFT等于各个DFT的圆周卷积再乘以1/N。第2章 离散傅里叶变换 2.5.4 有限长序列的线性卷积与圆周卷积有限长序列的线性卷积与圆周卷积 时域圆周卷积在频域上相当于两序列的DFT的乘积，而计算DFT可以采用它的快速算法快速傅里叶变换（FFT）（

42、见第3章），因此圆周卷积与线性卷积相比，计算速度可以大大加快。但是实际问题大多总是要求解线性卷积。例如，信号通过线性时不变系统，其输出就是输入信号与系统的单位脉冲响应的线性卷积，如果信号以及系统的单位脉冲响应都是有限长序列，那么是否能用圆周卷积运算来代替线性卷积运算而不失真呢？下面就来讨论这个问题。设x1(n)是N1点的有限长序列（0nN1-1），x2(n)是N2点的有限长序列（0nN2-1）。第2章 离散傅里叶变换（1）先看它们的线性卷积 1021212111)()()()()()()(Nmmmnxmxmnxmxnxnxnyx1(m)的非零区间为 0mN1-1 x2(n-m)的非零区间为 0

43、n-mN2-1(2-43)第2章 离散傅里叶变换 将两个不等式相加，得到 0nN1+N2-2 在上述区间外，不是x1(m)=0就是x2(n-m)=0，因而y1(n)=0。所以y1(n)是N1+N2-1 点有限长序列，即线性卷积的长度等于参与卷积的两序列的长度之和减1。例如，图2-16 中，x1(n)为N1=4 的矩形序列（图2-16（a），x2(n)为N2=5 的矩形序列（图2-16（b），则它们的线性卷积y1(n)为N=N1+N2-1=8 点的有限长序列（图 2-16（c）。第2章 离散傅里叶变换 （2）我们再来看x1(n)与x2(n)的圆周卷积。先假设进行L点的圆周卷积，再讨论L取何值时，

44、圆周卷积才能代表线性卷积。设y(n)=x1(n)x2(n)是两序列的L点圆周卷积，LmaxN1,N2，这就要将x1(n)与x2(n)都看成是L点的序列。在这L个序列值中，x1(n)只有前N1个是非零值，后L-N1个均为补充的零值。同样，x2(n)只有前N2个是非零值，后L-N2个均为补充的零值。则 L102121)()()()()()(LmLLnRmnxmxnxnxnyL（2-44）第2章 离散傅里叶变换 为了分析其圆周卷积，我们先将序列x1(n)与x2(n)以L为周期进行周期延拓)()()()()()(222111rLnxnxnxkLnxnxnxrLkL它们的周期卷积序列为)()()()()

45、()()()(1210121012101rLnymrLnxmxmrLnxmxmnxmxnyrLmrrLmLLm（2-45）第2章 离散傅里叶变换 前面已经分析了y1(n)具有N1+N2-1个非零值。因此可以看到，如果周期卷积的周期LN1+N2-1，那么y1(n)的周期延拓就必然有一部分非零序列值要交叠起来，从而出现混叠现象。只有在LN1+N2-1 时，才没有交叠现象。这时,在y1(n)的周期延拓 中，每一个周期L内，前N1+N2-1个序列值正好是y1(n)的全部非零序列值，而剩下的L-(N1+N2-1)个点上的序列值则是补充的零值。圆周卷积正是周期卷积取主值序列)(1ny)()()(1nRrL

46、nynyLr)()()()()(21nRnynxnxnyLL因此 第2章 离散傅里叶变换 所以要使圆周卷积等于线性卷积而不产生混叠的必要条件为 121NNL（2-47）满足此条件后就有)()(1nyny即 x1(n)x2(n)=x1(n)*x2(n)L 图2-16（d）、（e）、（f）正反映了（2-46）式的圆周卷积与线性卷积的关系。在图2-16（d）中，L=6小于N1+N2-1=8，这时产生混叠现象，其圆周卷积不等于线性卷积；而在图2-16（e）、（f）中,L=8和L=10，这时圆周卷积结果与线性卷积相同，所得y(n)的前8点序列值正好代表线性卷积结果。所以只要LN1+N2-1，圆周卷积结果

47、就能完全代表线性卷积。第2章 离散傅里叶变换 x1(n)n1N141230 x2(n)n112340N25y1(n)N1 N2 18n123405 6789 1011234(a)(b)(c)图2-16 线性卷积与圆周卷积第2章 离散傅里叶变换 图2-16 线性卷积与圆周卷积x1(n)x2(n)L6n12340 x1(n)x2(n)L8n1234x1(n)x2(n)L10n1234(d)(e)(f)12345012345670123456789LLL第2章 离散傅里叶变换 例例 2-8 一个有限长序列为)5(2)()(nnnx（1）计算序列x(n)的10点离散傅里叶变换。（2）若序列y(n)的D

48、FT为)()(1022kXekYkj式中，X(k)是x(n)的10点离散傅里叶变换，求序列y(n)。第2章 离散傅里叶变换（3）若10点序列y(n)的10点离散傅里叶变换是)()()(kWkXkY式中,X(k)是序列x(n)的10点DFT，W(k)是序列w(n)的10点DFT 01)(nw0n6 其他 求序列y(n)。第2章 离散傅里叶变换 解解（1）由式（2-30）可求得x(n)的10点DFT kkjknkNnnnkNeWWnnWnxkX)1(212121)5(2)()()(510251010101100 （2）X(k)乘以一个WNkm形式的复指数相当于是x(n)圆周移位m点。本题中m=-2

49、,x(n)向左圆周移位了2点，就有 y(n)=x(n+2)10R10(n)=2(n-3)+(n-8)第2章 离散傅里叶变换 （3）X(k)乘以W(k)相当于x(n)与w(n)的圆周卷积。为了进行圆周卷积，可以先计算线性卷积再将结果周期延拓并取主值序列。x(n)与w(n)的线性卷积为z(n)=x(n)*w(n)=1,1,1,1,1,3,3,2,2,2,2,2圆周卷积为)()10()(10nRrnznyr 在 0n9 求和中，仅有序列z(n)和z(n+10)有非零值，用表列出z(n)和z(n+10)的值，对n=0,1,2,9求和，得到：第2章 离散傅里叶变换 n0 1 2 3 4 5 6 7 8

50、910 11Z(n)z(n+10)1 1 1 1 1 3 3 2 2 22 2 0 0 0 0 0 0 0 02 20 0y(n)3 3 1 1 1 3 3 2 2 2_ _所以10点圆周卷积为 y(n)=3,3,1,1，1,3,3,2,2,2 第2章 离散傅里叶变换 2.5.5 共轭对称性共轭对称性 设x*(n)为x(n)的共轭复序列，则 DFTx*(n)=X*(-k)NRN(k)=X*(N-k)NRN(k)=X*(N-k)0kN-1 且 X(N)=X(0)（2-48）证证)(*)()(*)()()()(*)()()()(*)(*10)(10*10kNXkRkNXkRWnxkRkXkRWnx