频率域采样.ppt

第 * 页 X 第 * 页 X 3.3 频率域采样 青海民族大学物电学院 引入：DFT的定义 表示在单位圆周上对x(n)的Fourier变换N点等间隔采样。 长度为M的有限长序列 长度为N的有限长序列 则： 的离散Fourier级数： 频域采样定理： 时域离散周期信号 频域离散周期信号 采样 离散化 周期化 采样 离散化 周期化 时域采样定理是保证频域周期延拓时不混叠 频域采样定理是保证时域周期延拓时不混叠 DFT的实质 DFT的定义 用X(K)表示X(Z)和X(ejw)内插公式和内插函数 将上式代入X(z)的表示式中得： 式中WN-kN=1 内插公式 内插函数 当z=ejω时，代入上式整理得： 频域内插公式 频域内插函数 【例3.3.1】( ep331.m) 长度为26的三角形序列x(n)如图3.3.1(a)所示。编写MATLAB程序验证频域采样理论。 解题思想： 序列的长度M=26,要验证采样定理， 即取采样点数N满足采样定理和不 满足采样定理的点数，看是否发生 时域的混叠。从而验证采样定理。 所以N取32和16。 最后分别对X16(k)和X32(k)求IDFT MATLAB求解程序ep331.m如下 M=26; N=32; n=0:M; xa=0:M/2; xb= ceil(M/2)-1:-1:0; xn=［xa, xb］; 　 %产生M长三角波序列x(n) Xk=fft(xn, 512); 　 %512点FFT［x(n)］ X32k=fft(xn, 32); %32点FFT［x(n)］ X32n=ifft(X32k); %32点IFFT［X32(k)］得到x32(n) X16k=X32k(1:2:N); %隔点抽取X32k得到X16(k) X16n=ifft(X16k, N/2); %16点IFFT［X16(k)］得到x16(n) 绘图部分省略 图3.3.1 频域采样定理验证 512点FFT 16点FFT 32点FFT 16点IFFT 32点IFFT 结论：当采样点数N=16M（26）时，由于存在时域混叠失真，x16(n)≠x(n)； 当采样点数N=32M（26）时，无时域混叠失真，x32(n)=IDFT［X32(k)］=x(n)。 第 * 页 X 第 * 页 X