离散傅里叶变换数字信号处理陈后金课件.pptxVIP

离散傅里叶变换数字信 号处理陈后金 05 十月 2023 • 了解四种信号的傅里叶变换的数学概念及特点。 • 深刻理解有限长序列DFT的定义及概念。 • 掌握序列DFT与序列DTFT和z变换的相互关系 。 • 掌握利用DFT分析任意信号频谱的原理和方法 。 • 理解DFT分析信号频谱中出现的现象以及改善这些 现象的方法 。 • 掌握利用DFT实现序列线性卷积的原理和方法。 学 习 要 求 重 点 和 难 点 信号的频域分析在信息技术领域广泛应用 为什么进行信号频谱的数值化分析? 1.许多实际信号不存在数学解析式 2.利用计算机数值计算,简单快捷 问题的提出 四种信号傅里叶表示 有限长序列离散傅里叶变换 DFT矩阵表示 利用nVTLVB计算DFT 有限长序列的傅里叶分析 频谱特点： 离散非周期谱 1. 周期为T0的连续时间周期信号 四种信号傅里叶表示 频谱特点： 连续非周期谱 四种信号傅里叶表示 2. 连续时间非周期信号 频谱特点： 周期为2冗的连续谱 四种信号傅里叶表示 3. 离散非周期信号 频谱特点： 周期为N的离散谱 4. 周期为N的离散周期信号 四种信号傅里叶表示 有限长序列离散傅里叶变换 符号表示 IDFT DFT 有限长序列x[k]离散傅里叶变换X[m]是其离散时 间傅里叶变换X(ejW)在一个周期[0,2p)的等间隔抽样 有限长序列DFT与DTFT关系 DFT可以看成是截取DFS的主值区间构成的变换对 DFT与DFS关系 例3: 求有限长4点序列 的DFT。 例4 有限长4点序列DFT矩阵表示。 例题： 例4: 求有限长4点序列 的DFT。 解: 如果序列后补零，其DFT有何变化? X[m]={2 ，2 ，-2 ，2}, m=0,1,2,3 有限长4点序列DFT矩阵表示: DFT矩阵表示 DFT矩阵表示 fft(x) 计算M点的DFT 。M是序列x 的长度。 fft(x,N) 计算N点的DFT。 MN ，将原序列裁为N点计算N点的DFT； MN ，将原序列补零至N点，然后计算N点DFT。 fft(x) fft(x,N) ifft(x) ifft(x,N) 利用MATLAB计算DFT 利用MATLAB计算DFT---N=4、8、16 stem(0:15,abs(xm2)); axis([0 16 -1 3]); stem(0:7,abs(xm1)); axis([0 8 -1 3]); stem(0:3,abs(xm)); axis([0 4 -1 3]); xm2=fft(x,16); subplot(313); xm=fft(x,4); subplot(311); xm1=fft(x,8); subplot(312); x=[1 1 -1 1]; 利用MATLAB计算DFT---N=16、64、256 x=[0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0]; xm=fft(x,16); subplot(311); stem(0:15,abs(xm)); subplot(312); xm1=fft(x,64); stem(0:63,abs(xm1)); axis([0 64 -1 7]); subplot(313); xm2=fft(x,256); stem(0:255,abs(xm2)); axis([0 256 -1 7]); axis([0 16 -1 7]); x=[0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0]; x1=[0 0 0 0 0 1];N1=6; xm1=fft(x1);subplot(211);stem(0:N1-1,abs(xm1)); xm16=fft(x1,16);subplot(212);stem(0:15,abs(xm16)); x2=[0 0 0 0 0 1 1];N2=7; xm2=fft(x2);subplot(211);stem(0:N2-1,abs(xm2));grid xm16=fft(x2,16);subplot(212);stem(0:15,abs(xm16));grid x3=[0 0 0 0 0 1 1 1];N3=8; xm3=fft(x3);subplot(211);stem(0:N3-1,abs(xm3)); xm16=fft(x3,16);subplot(212);stem(0:15,abs(xm16));