第三章离散傅里叶变换(数字信号处理).pptVIP

DSP 本章内容 一、引言 傅立叶变换 傅立叶级数 序列的傅立叶变换 离散傅立叶变换 四种形式归纳 二、离散傅立叶级数（DFS） 说明 DFS变换对公式表明，一个周期序列虽然是无穷长序列，但是只要知道它一个周期的内容（一个周期内信号的变化情况），其它的内容也就都知道了，所以这种无穷长序列实际上只有N个序列值的信息是有用的，因此周期序列与有限长序列有着本质的联系。 三、离散傅立叶变换DFT 定 义 定 义 有限长序列的DFT是有限长的 DFT与DFS无本质区别，DFT是DFS的主值 DFT和DTFT，ZT，DFS的关系 设序列x(n)的长度为M 例 周期延拓 DFT的性质 1、DFT的线性 3、共轭对称性 圆周共轭对称 圆周共轭反对称 几种特例 4、DFT形式下的帕塞瓦定理 5、循环卷积 时域循环卷积过程： 1）补零 2）周期延拓 3）翻褶，取主值序列 4）循环移位 5）相乘相加 图解 循环卷积程序 x1 = [1,2,2]; x2 = [1,2,3,4]; y = circonvt(x1,x2,4) Stem(y) 线性卷积 function [y,ny]=conv_improve(x,nx,h,nh) %[x,nx]为第一个信号;%[h,nh]为第二个信号 %conv(x,h)可以实现两个有限长度序列的卷积 ny1=nx(1)+nh(1); ny2=nx(length(x))+nh(length(h)); ny=[ny1:ny2]; y=conv(x,h); 在命令窗口调用卷积函数。x=[3 4 0 -2 2 3 5]; nx=[-3:3]; h=[1 4 5 6 0 1]; nh=[N:N+5]; 例 四、频域采样定理 四、频域采样定理 采样Z变换 采样Z变换 频域采样恢复 内插函数的特性 频响特性 续1 总 结 总 结 五、利用DFT对连续时间信号的逼近（选讲） 图例1 图例2 图例3 续2 利用DFT计算连续时间信号的几个问题 混叠现象 例 频率泄漏 续1 栅栏效应 频率分辨力 分析 当x(n)为线性时不变系统的单位冲激响应h(n)时： 即H(z)可看做是由N个采样带通滤波器并联构成。 对连续信号，能用解析式精确表达其频谱信息。而实际应用中，信号受各种环境噪声或异常事件影响，不能明确知道信号的数学解析表达式，只能作数值分析。 利用计算机用DFT对包罗万象的信号进行分析和合成是当前主要的应用方法。 %Matlab程序 clear;clc; N=32; n=0:N-1; x=1.0*sin(1*n*2*pi/N) +0.1*sin(2*n*2*pi/N) +0.3*sin(3*n*2*pi/N); y=fft(x); subplot(2,1,1); plot(x); hold on; stem(x); subplot(2,1,2); stem(abs(y)) %Matlab程序 clear;clc; N=32; n=0:2*N-1; x=1.0*sin(1*n*2*pi/N) +0.1*sin(2*n*2*pi/N) +0.3*sin(3*n*2*pi/N); y=fft(x); subplot(2,1,1); plot(x); hold on; stem(x); subplot(2,1,2); stem(abs(y)) %Matlab程序 clear; N=32; n=0:N-1; x=power(0.9,n); y=fft(x); subplot(2,1,1); plot(x); hold on; stem(x); subplot(2,1,2); stem(abs(y)) 利用DFT对连续时间信号的逼近 1 2 3 0 0 0 0 1 2 3 0 0 0 0 1 2 3 0 0 0 0 h((m))NRN Y(6)=3 0 0 0 0 3 2 1 0 0 0 0 3 2 1 0 0 0 0 3 2 1 h((6-m))NRN Y(5)=8 0 0 0 3 2 1 0 0 0 0 3 2 1 0 0 0 0 3 2 1 0 h((5-m))NRN Y(4)=14 Y(3)=20 Y(2)=26 Y(1)=14 Y(0)=5 y(n) 0 0 3 2 1 0 0 0 0 3 2 1 0 0 0 0 3 2 1 0 0 h((4-m))NRN 0 3 2 1 0 0 0 0 3 2 1 0 0 0 0 3 2 1 0 0 0 h((3-m))NRN 3