第3章离散变换dft136.pptxVIP

第3章 离散傅里叶变换(DFT);本章重点：;§3.1 离散傅里叶变换的定义 ;IDFT［X(k)］的唯一性证明：;function [Xk] = dft(xn,N) % Computes Discrete Fourier Transform % % [Xk] = dft(xn,N) % Xk = DFT coeff. array over 0 = k = N-1 % xn = N-point finite-duration sequence % N = Length of DFT %----------------------------------- n = [0:1:N-1]; % row vector for n k = [0:1:N-1]; % row vecor for k WN = exp(-j*2*pi/N); % Wn factor nk = n*k; % creates a N by N matrix of nk values WNnk = WN .^ nk; % DFT matrix Xk = xn * WNnk; % row vector for DFT coefficients;function xn = idft(Xk,N) % Computes Inverse Discrete Fourier Transform % xn = idft(Xk,N) % Xk = DFT coeff. array over 0 = k = N-1 % xn = N-point finite-duration sequence % N = Length of DFT %----------------------------------- n = [0:1:N-1]; % row vector for n k = [0:1:N-1]; % row vecor for k WN = exp(-j*2*pi/N); % Wn factor nk = n*k; % creates a N by N matrix of nk values WNnk = WN.^(-nk); % DFT matrix xn =(Xk*WNnk)./N; % xn is N-point sequence suplot(211);stem(k,abs(Xk); suplot(212);stem(n,real(xn);;求序列x(n)={1, 2, 3, 4}4点DFT变换;例3.1.1;DFT结果与变换区间长度N有关;设序列x(n)的长度为N， 其DFT 、Z变换和DTFT分别为：;;例 3.1.1; X(k)与x(n)均为有限长序列，但由于WNkn的周期性， X(k)隐含周期性，且周期均为N。;主值区间和主值序列 任何周期为N的周期序列 可以看作长度为N的有限长序列x(n)的周期延拓序列，x(n)是 的一个周期。; x((n))N表示x(n)以N为周期的周期延拓序列，((n))N表示n对N求余。 如果? n=MN+n1，0≤n1≤N-1，M为整数， 则 ((n))N=n1;图 3.1.2 有限长序列及其周期延拓 ;周期序列的离散傅立叶级数 如果x(n)的长度为N， 且 =x((n))N， 则 的 离散傅里叶级数表示为 上式中的 。 结论：有限长序列的离散傅立叶变换X(k) 是x(n)的周期延拓序 列x((n))N的离散傅立叶级数系数 的主值序列。;关系：;各种傅立叶变换参数的比较;§3.2 离散傅里叶变换的基本性质; 设x(n)为有限长序列，长度为N，则x(n)的循环移 位定义为 y(n)=x((n+m))NRN(n) (3.2.2);图 3.2.1 循环移位过程示意图 ;2. 圆周位移的含义 由于我们取主值序列，即只观察n=0到N-1这一主 值区间，当某一抽样从此区间一端移出时，与它相同 值的抽样又从此区间的另一端进来。如果把x(n) 排列 一个N等分的圆周上，序列的移位就相当于x(n) 在圆 上旋转，故称作圆周移位。当围着圆周观察几圈时， 看到就是周期序列： 。;用MATLAB实现有限长序列的循环移位; 设x(n)是长度为N的有限长序列，y(n)为x(n)的循环移位， 即