数字信号处理(丁玉美版)教案第三章-.pptVIP

第三章 离散傅立叶变换 Discrete Fourier Transform 内容提要 频率域采样 DFT的应用 3.4 频率域采样 1、频率采样定理 对任意 ，有 且 的收敛域包含单位圆，则 表示在区间[0，2?]上对 的 傅立叶变换 的N点等间隔采样。用 有限长序列 表示： 频域采样定理 如果x（n）的长度为M，则只有当频域采样点数N?M时，才有 例 已知 求其FT，再进行4点和8点采样，通过IDFT 求其时域表示； 解： 3.5 DFT的应用 DFT的快速算法FFT的出现， 使DFT在 数字通信、 语言信号处理、 图像处理、 功率谱估计、 仿真、 系统分析、 雷达 理论、 光学、 医学、 地震以及数值分 析等各个领域都得到广泛应用。 3.5.1 用DFT计算卷积 （1）用DFT计算循环卷积 在实际应用中，为了分析时域离散线性移不变系统或者对序列进行滤波处理等，需要计算两个序列的线性卷积。 请问可以使用直接使用DFT计算两个序列线性卷积吗？ （2）循环卷积与线性卷积的关系 循环卷积等于线性卷积以L为周期的 周期延拓序列的主值序列。 只有当循环卷积长度L?N+M-1时,以 L为周期进行周期延拓才无混叠现象。 因此，循环卷积等于线性卷积的条件是：L?N+M-1 For examples h(n)=[1, 2, –1, 1] M=4 x(n)=[1, 1, 2, 1, 2, 2, 1, 1] N=8 1.Compute the linear convolution. 2.For the values L=9, 11, 13, compute the circular convolution of the two signals of x(n) and h(n). Circular convolution(循环卷积）: L=9时 yc(n)=[1,4,3,5,3,7,4,3,3] 比较 yl(n)= [1,3,3,5,3,7,4,3,3,0,1] L=11时 yc(n) =[1,3,3,5,3,7,4,3,3,0,1] L=13时 yc(n) =[1,3,3,5,3,7,4,3,3,0,1,0,0] %Linear convolution h=[1, 2, -1, 1] ; x=[1, 1, 2, 1, 2, 2, 1, 1] ; yl=conv(h,x); subplot(4,1,1) stem(yl); title(‘线性卷积输出结果’); yc9=circonv(h,x,9); subplot(4,1,2) stem(yc9); title(‘9点循环卷积输出结果’); only L ? N+M-1= 8+4-1=11 linear convolution（线性卷积） =circular convolution（循环卷积） 如果 L ? N+M-1，可用DFT（FFT）求线性卷积。 Matlab 应用举例 (演示文件见dft4conv.m) %利用DFT计算线性卷积 x=[1 2 0 1 ];h=[2 2 1 1]; %计算线性卷积的长度 L=length(x)+length(h)-1; X=fft(x,L);H=fft(h,L); %计算频域乘积的IDFT y1=ifft(X.*H); %画出卷积结果及误差 k=0:L-1;subplot(2,1,1); stem(k,real(y1));axis([0 6 07]); title(‘result of linear convoltution’); xlabel(‘Time index k’);ylabel(‘Amplitude’); y2=conv(x,h);error=y1-y2; subplot(2,1,2);stem(k,abs(error)); xlabel(‘Time index k’);ylabel(‘Amplitude’); title(‘幅度误差’);%计算误差 3.5.2 Overlap-Add Methods （重叠相加法） 当输入信号的长度很长或为无限长时， 输出信号也将很长或无限长，那么 L ? N+M-1 的条件就很难能满足。 于是将长序列截短，分段求输出 这就是重叠相加法的基本思路。 设 长度为N， 为无限长序列。将 分为长