第3章_离散傅里叶变换.ppt

第3章 离散傅里叶变换 (discrete fourier transformdft),3.1 离散傅里叶变换的定义及物理意义 3.2 离散傅里叶变换的性质 3.3 频率域采样 3.4 dft的应用举例,傅里叶变换的离散性和周期性,1. 连续时间、连续频率连续时间非周期信号的傅里叶变换,结论：时域非周期-频域连续；时域连续-频域非周期,2. 连续时间、离散频率连续时间周期信号的傅里叶级数,结论：时域周期-频域离散；时域连续-频域非周期,3. 离散时间、连续频率序列的傅里叶变换,结论：时域非周期-频域连续；时域离散-频域周期,4. 离散时间、离散频率离散傅里叶变换,结论：时域周期-频域离散；时域离散-频域周期,3.1 离散傅里叶变换的定义及物理意义,一、离散傅里叶变换(dft)的定义 二、dft与ft及zt的关系 三、dft的隐含周期性 四、用matlab计算序列的dft,“借用” 的主值序列x(k)定义为“离散傅里叶变换(dft)”。目的是使傅里叶分析可以利用数字计算机。,一、离散傅里叶变换(discrete fourier transformdft),(1),(4),(2),(3),【例3.1.1】x(n)=r4(n)，求x(n)的4点和8点dft。,dftx(n)结果与n的取值有关,二、dft与傅里叶变换和z变换的关系,dft的物理意义1：序列x(n)的n点dft是x(n)的z变换在单位圆上的n点等间隔采样,dft的物理意义2：x(k)为x(n)的傅里叶变换在区间0, 2上的n点等间隔采样,二、dft与傅里叶变换和z变换的关系,三、dft的隐含周期性,dft物理意义3：有限长序列的n点离散傅里叶变换x(k)正好是x(n)的周期延拓序列x(n)n的离散傅里叶级数系数的主值序列，x(k)实质上是x(n)的周期延拓序列 的频谱特性。,xk = fft(xn,n) xn = ifft(xk,n),四、用matlab计算序列的dft,【例3.1.2】,3.2 离散傅里叶变换的性质,一、线性性质,二、循环移位性质,1. 序列的循环移位,2. 时域循环移位定理,3. 频域循环移位定理,二、循环移位性质,1. 两个有限长序列的循环卷积,三、循环卷积定理,h(n)与x(n)的l点循环卷积,循环倒相序列,x(n)的l点循环卷积矩阵,【例3.2.1】,1. 两个有限长序列的循环卷积,2. 时域循环卷积定理,三、循环卷积定理,3. 频域循环卷积定理,四、复共轭序列的dft,五、dft的共轭对称性,1. 有限长共轭对称序列和共轭反对称序列,任何有限长序列都可以表示成其共轭对称分量和共轭反对称分量之和,2. dft的共轭对称性,(1) 将序列x(n)表示为,有限长序列x(n)的实部的dft为x(k)的共轭对称分量，虚部和j一起的dft为x(k)的共轭反对称分量。,2. dft的共轭对称性,(2)将序列x(n)表示为,有限长序列x(n)的共轭对称分量和共轭反对称分量的dft分别为x(k)的实部和虚部乘以j。,2. dft的共轭对称性,(3) 有限长实序列dft的共轭对称性,【例3.2.2】利用dft的共轭对称性，设计一种高效算法，通过计算一个n点dft，就可以计算出两个实序列x1(n)和x2(n)的n点dft。,六、离散帕塞瓦尔定理,3.3 频率域采样,由于,所以,频域采样定理： 假设序列x(n)的长度为m，只有当频域采样点数nm时，才有xn(n)=idftx(k)=x(n)，即可由频域采样x(k)恢复原序列x(n)，否则产生时域混叠现象。,3.4 dft的应用举例,一、用dft计算线性卷积 二、用dft对信号进行谱分析,一、用dft计算线性卷积,当ln+m-1时，yc(n)=yl(n),线性卷积与循环卷积的关系：,1. 用dft对连续时间非周期信号进行谱分析,(1) 对x(t)以采样间隔t采样得序列 x(n)= x(nt),(2) 将序列x(n)=x(nt)截断成从t=0开始长度为tp的有限长序列，tp=nt,(3) 为了数值计算，在频域也要离散化，一个周期fs等间隔采样n点，每个样点间隔为 f，n=fs/f，f=fs/n=1/tn=1/tp,二、用dft对信号进行谱分析,图3.4.6 用dft分析连续信号频谱的原理示意图,用dft近似计算模拟信号频谱的计算步骤总结如下： (1) 首先确定用dft对模拟信号频谱进行近似计算的三个参数，即频率分辨率f、 采样频率fs、 记录时间tp。 (2) 用已确定的fs对模拟信号采样。采样后得到时域离散信号为 x(n)=xa(t)|t=nt=xa(nt) (3) 在计算机上调dft(fft)函数对信号x(n)进行频谱计算。,【例3.4.2】对实时信号进行谱分析，要求谱线间距(分辨率)f10hz，信号最高频率fc=2.5khz，试确定最小记录时间tpmin、最大采样间隔tmax、最少的采样点数nmin。如果fc不变，要求谱分辨率f增加1倍，nmin和tpmin是多少？,【例3.4.3】用dft对模拟信号进行谱分析，设模拟信号xa(t)的最高频率为200hz，以采样频率fs=400hz采样得到时域离散序列x(n)=xa(nt)，要求频率分辨率为10hz。模拟信号频谱xa(j)如下图所示，试画出x(ej)=ftx(n)和x(k)=dftx(n) 谱线图，并标出每个k=0,20,40对应的数字频率k和模拟频率fk的值。,2. 对连续时间周期信号的傅里叶级数的dfs逼近,(2) 将频域离散序列加以截断，截断长度为一个周期n,(1) 时域抽样,(1),(2),3. 用dft对序列进行谱分析,(2) 对周期序列进行谱分析,(1) 对有限长序列进行谱分析,4. 用dft进行谱分析的误差问题,(1) 频谱混叠,(2) 栅栏效应,(3) 截断效应,【例3.4.5】50hz正弦波xa(t)=sin(250t)，它的幅度曲线是线状谱，如图(a)所示。如果将它截取0.09s的一段，相当于将它乘以长度为0.09 s矩形窗函数，即xa(t)rtp(t)，tp=0.09s, 对该信号进行谱分析