第3章 离散傅里叶变换幻灯片.pptVIP

有限长序列的N点离散傅里叶变换X(k)正好是x(n)的周期延拓序列x((n))N的离散傅里叶级数系数的主值序列， ， X(k)实质上是x(n)的周期延拓序列 的频谱特性。 【例2.3.1】设x(n)=R4(n)，将x(n)以N=8为周期进行周期延拓，得到如图2.3.1(a)所示的周期序列，周期为8。求 ，并画出它的幅度谱。 Xk = fft(xn,N) xn = ifft(Xk,N) 频域采样定理： 假设序列x(n)的长度为M，只有当频域采样点数N≥M时，才有xN(n)=IDFT[X(k)]=x(n)，即可由频域采样X(k)恢复原序列x(n)，否则产生时域混叠现象。 用DFT近似计算模拟信号频谱的计算步骤总结如下： (1) 首先确定用DFT对模拟信号频谱进行近似计算的三个参数，即频率分辨率F、 采样频率Fs、 记录时间Tp。 (2) 用已确定的Fs对模拟信号采样。采样后得到时域离散信号为 x(n)=xa(t)|t=nT=xa(nT) (3) 在计算机上调DFT(FFT)函数对信号x(n)进行频谱计算。 六、离散帕塞瓦尔定理 3.3 频率域采样 由于 所以 3.4 DFT的应用举例 一、用DFT计算线性卷积 二、用DFT对信号进行谱分析 一、用DFT计算线性卷积 当L≥N+M-1时，yc(n)=yl(n) 线性卷积与循环卷积的关系： 序列yl(n)的长度为N+M-1 序列yc(n)的长度为L 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分析连续信号频谱的原理示意图 【例3.4.2】对实时信号进行谱分析，要求谱线间距(分辨率)F≤10Hz，信号最高频率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ω)=FT[x(n)]和X(k)=DFT[x(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(2π·50t)，它的幅度曲线是线状谱，如图(a)所示。如果将它截取0.09s的一段，相当于将它乘以长度为0.09 s矩形窗函数，即xa(t)RTp(t)，Tp=0.09s, 对该信号进行谱分析肯定和原信号的谱分析不同，它的幅度谱如图(b)所示。如果Tp=0.18s，幅度谱如图(c)所示。 (a) 50 Hz正弦波的幅频曲线； (b) 50 Hz正弦波加窗后的幅频曲线(Tp=0.09 s); (c) 50 Hz正弦波加窗后的幅频曲线(Tp=0.18 s) 第3章 离散傅里叶变换 (Discrete Fourier Transform—DFT) 3.1 离散傅里叶变换的定义及物理意义 3.2 离散傅里叶变换的性质 3.3 频率域采样 3.4 DFT的应用举例 傅里叶变换的离散性和周期性 1. 连续时间、连续频率——连续时间非周期信号的傅里叶变换 结论：时域非周期-?频域连续；时域连续-?频域非周期 2. 连续时间、离散频率——连续时间周期信号的傅里叶级数 结论：时域周期-?频域离散；时域连续-?频域非周期 3. 离散时间、连续频率——序列的傅里叶变换 结论：时域非周期-?频域连续；时域离散-?频域周期 4. 离散时间、离散频率——离散傅里叶变换 结论：时域周期-?频域离散；时域离散-?频域周期 3.1 离散傅里叶变换的定义及物理意义 一、离散傅里叶变换(DFT)的定义 二、DFT与FT及ZT的关系