节 离散傅立叶变换.pptx
有限长序列的傅里叶分析有限长序列的傅里叶分析一、四种信号傅里叶表示一、四种信号傅里叶表示1.周期为T0的连续时间周期信号频谱特点：离散非周期谱第1页/共63页2.连续时间非周期信号频谱特点：连续非周期谱第2页/共63页3.离散非周期信号频谱特点：周期为2的连续谱第3页/共63页4.周期为周期为N N 的离散周期信号的离散周期信号频谱特点：周期为N的离散谱第4页/共63页为了便于更好地理解DFT的概念，先讨论周期序列及其离散傅里叶级数（DFS）表示。一个周期为N的周期序列，即，k为任意整数，N为周期周期序列不能进行Z变换，因为其在n=-到+都周而复始永不衰减，即 z 平面上没有收敛域。但是，正象连续时间周期信号可用傅氏级数表达，周期序列也可用离散的傅氏级数来表示，也即用周期为N的正弦序列来表示。离散傅里叶级数（DFS）第5页/共63页周期为N的正弦序列其基频成分为：K次谐波序列为：但离散级数所有谐波成分中只有N个是独立的，这是与连续傅氏级数的不同之处，即因此第6页/共63页将周期序列展成离散傅里叶级数时，只需取k=0到(N-1)这N个独立的谐波分量，所以一个周期序列的离散傅里叶级数只需包含这N个复指数，利用正弦序列的周期性可求解系数。将上式两边乘以，并对一个周期求和第7页/共63页第8页/共63页上式中部分显然只有当k=r时才有值为1,其他任意k值时均为零,所以有或写为1)可求N次谐波的系数2）也是一个由N个独立谐波分量组成的傅立叶级数3）为周期序列，周期为N。第9页/共63页时域上周期序列的离散傅里叶级数在频域上仍是一个周期序列。第10页/共63页是一个周期序列的离散傅里叶级数(DFS)变换对，这种对称关系可表为：习惯上：记，第11页/共63页DFS变换对公式表明，一个周期序列虽然是无穷长序列，但是只要知道它一个周期的内容（一个周期内信号的变化情况），其它的内容也就都知道了，所以这种无穷长序列实际上只有N个序列值的信息是有用的，因此周期序列与有限长序列有着本质的联系。则DFS变换对可写为DFS离散傅里叶级数变换IDFS离散傅里叶级数反变换。第12页/共63页DDFS的几个主要特性：假设都是周期为N的两个周期序列，各自的离散傅里叶级数为：1）线性a,b为任意常数第13页/共63页 2）序列移位证因为及都是以N为周期的函数，所以有第14页/共63页由于与对称的特点，同样可证明第15页/共63页 3）共轭对称性对于复序列其共轭序列满足证:同理:第16页/共63页进一步可得共轭偶对称分量共轭奇对称分量第17页/共63页4）周期卷积若则或第18页/共63页周期卷积第19页/共63页证：这是一个卷积公式，但与前面讨论的线性卷积的差别在于，这里的卷积过程只限于一个周期内（即m=0N-1），称为周期卷积。例：、，周期为N=7,宽度分别为4和3，求周期卷积。结果仍为周期序列，周期为N。第20页/共63页由于DFS与IDFS的对称性，对周期序列乘积，存在着频域的周期卷积公式，若则第21页/共63页我们知道周期序列实际上只有有限个序列值有意义，因此它的许多特性可推广到有限长序列上。一个有限长序列x(n)，长为N，为了引用周期序列的概念，假定一个周期序列，它由长度为N的有限长序列x(n)延拓而成，它们的关系：离散傅里叶变换（DFT）第22页/共63页周期序列的主值区间与主值序列：对于周期序列，定义其第一个周期n=0N-1，为的“主值区间”，主值区间上的序列为主值序列x(n)。x(n)与的关系可描述为：数学表示：RN（n）为矩形序列。符号（n）N是余数运算表达式，表示n对N求余数。第23页/共63页第24页/共63页例：是周期为N=8的序列，求n=11和n=-2对N的余数。因此第25页/共63页频域上的主值区间与主值序列：周期序列的离散付氏级数也是一个周期序列，也可给它定义一个主值区间，以及主值序列X(k)。数学表示：第26页/共63页再看周期序列的离散傅里叶级数变换（DFS）公式：这两个公式的求和都只限于主值区间（0N-1），它们完全适用于主值序列x(n)与X(k)，因而我们可得到一个新的定义有限长序列离散傅里叶变换定义。第27页/共63页长度为N的有限长序列x(n)，其离散傅里叶变换X(k)仍是一个长度为N的有限长序列，它们的关系为：x(n)与X(k)是一个有限长序列离散傅里叶变换对，已知x(n)就能唯一地确定X(k)，同样已知X(k)也就唯一地确定x(n)，实际上x(n)与X(k)都是长度为N的序列（复序列）都有N个独立值，因而具有等量的信息。有限长序列隐含着周期性。第28页/共63页1.线性需将较短序列补零后，再按长序列的点数做DFT2.循环位移(Circular shift of a sequence)循环位移定义为离散傅里叶变换的性质第29页/共63页第30页/共63页DFT频域循环位移特性DFT时域循环位移特性第31页/共63页3.对称性(symmetry)周期共轭对称(Periodicconjugatesymmetry)定义为周期共轭反对称(Periodicconjugateantisymmetry)定义为当序列xk为实序列时，周期偶对称序列满足当序列xk为实序列时，周期奇对称序列满足第32页/共63页对称特性当xk是实序列时第33页/共63页4.循环卷积h(-n)Nh(1-n)Nh(2-n)Nh(3-n)N第34页/共63页卷积定理第35页/共63页序列DFT与z变换的关系xk的 Xm等于其z变换X(z)在单位圆上等间隔取样第36页/共63页设序列xk的长度为NmXIDFTkx变换Z)(zX(内插公式)第37页/共63页问题提出:实际需要：LTI系统响应yk=x khk可否利用DFT计算线性卷积？例：x1k=1,1,1,x2k=1,1,0,1,N=4一、两个有限长序列的线性卷积利用DFT计算线性卷积第38页/共63页第39页/共63页线性卷积的矩阵表示线性卷积的矩阵表示第40页/共63页第41页/共63页循环卷积的矩阵表示循环卷积的矩阵表示第42页/共63页第43页/共63页循环卷积的矩阵表示第44页/共63页若若xk的长度为的长度为N，hk的长度为的长度为M，则，则L=N+M-1点循环卷积等于点循环卷积等于xk 与与hk的线性的线性卷积。卷积。第45页/共63页直接计算与由DFT间接计算结果比较第46页/共63页若x1k为 M 点序列,x2k为L 点序列，L Mx1k L x2k中哪些点不是线性卷积的点?问题讨论第47页/共63页0 k M-2不是线性卷积不是线性卷积的结果，即前的结果，即前(M-1)个点与个点与线性卷积不一样。线性卷积不一样。第48页/共63页线性卷积的矩阵表示线性卷积的矩阵表示第49页/共63页循环卷积的矩阵表示循环卷积的矩阵表示第50页/共63页 x1k L x2kk=0 M-2,前M-1个点不是线性卷积的点k=M-1 L-1,L-M+1个点与线性卷积的点对应线性卷积 L L+M-2 后M-1点没有计算则L点循环卷积结论若x1k为 M 点序列,x2k为L 点序列，L M第51页/共63页长序列和短序列的线性卷积直接利用DFT计算的缺点：(1)信号要全部输入后才能进行计算，延迟太多(2)内存要求大(3)算法效率不高解决问题方法：采用分段卷积分段卷积可采用重叠相加法 和 重叠保留法第52页/共63页1.重叠相加（overlap add）将长序列xk分为若干段长度为L的序列其中第53页/共63页y0k的非零范围y1k-L的非零范围序列y0k,y1k的重叠部分重叠的点数L+M-2-L+1=M-1依次将相邻两段的依次将相邻两段的M-1个重叠点相加，即得到最终的个重叠点相加，即得到最终的线性卷积结果。线性卷积结果。第54页/共63页重叠相加法分段卷积举例第55页/共63页第56页/共63页方法：(1)将xk长序列分段，每段长度为L；(2)各段序列xnk与M点短序列hk循环卷积；(3)从各段循环卷积中提取线性卷积结果。2.重叠保留法(overlap save)前M-1个点不是线性卷积的点因ynk=xn kL hk故分段时，每段与其前一段有M-1个点重叠。第57页/共63页-x k(M 1)M-1-L(M 1)L-1x 0kx 1k2L-Mk第一段前需补M-1个零第58页/共63页记ynk=xn kL hk第59页/共63页y0k中的中的M-1,L-1点对应于线性卷积点对应于线性卷积 xk*hk中的中的0,L-M点点y1k中的M-1,L-1点对应于线性卷积xkhk中的 L-(M-1),2L-M-(M-1)点第60页/共63页例例 已知序列已知序列xk=k+2,0 k 12,hk=1,2,1试分别利用重叠相加和保留法计算试分别利用重叠相加和保留法计算线性卷积线性卷积,取取L=5。yk=2,7,12,16,20,24,28,32,36,40,44,48,52,41,14解:重叠相加法x1k=2,3,4,5,6x2k=7,8,9,10,11x3k=12,13,14y1k=2,7,12,16,20,17,6y2k=7,22,32,36,40,32,11y3k=12,37,52,41,14第61页/共63页解解:重叠保留法重叠保留法yk=2,7,12,16,20,24,28,32,36,40,44,48,52,41,14x1k=0,0,2,3,4x2k=3,4,5,6,7x3k=6,7,8,9,10y1k=x1khk=11,4,2,7,12x4k=9,10,11,12,13y2k=x2khk=23,17,16,20,24y3k=x3khk=35,29,28,32,36y4k=x4khk=47,41,40,44,48x5k=12,13,14,0,0y5k=x5khk=12,37,52,41,14第62页/共63页感谢您的观看。第63页/共63页
