数字信号处理-第3章.ppt

数字信号处理－离散傅立叶变换 * (2)对周期序列的频谱分析 设序列 (n)=x(n+rN)是周期为N的周期序列，则其 傅立叶变换为： 周期序列的频谱结构可以用离散傅里叶级数系数 表示 取 的主值序列 进行N点DFT，得到 周期序列的频谱结构也可以用其主值序列的离散傅里叶变换X(k)来表示(分析) x ~ ~ x (n) 数字信号处理－离散傅立叶变换 * 令：n=n’+rN，r=0,1,…m-1，n’=0,1…N-1,则 x ~ xM(n)= (n)?RM(n) 即：M=mN，m为整数 截取序列的长度M为 (n)的整数个周期 x ~ 设：n=n’+rN 数字信号处理－离散傅立叶变换 * 3.4 DFT的应用举例 周期序列的频谱结构也可以用xM(k)表示 分析: (1)只有在k=rm时，XM(rm)=m ，表示 (n)的r次谐波谱线，幅度扩大了m倍，在其它k值， XM(k)＝0。 (2)X(r)与XM(rm)对应点的频率相等。 (3)只要截取 (n)整数个周期进行DFT，就可得到它的频谱结构，达到谱分析的目的。 k/m=整数 k/m≠整数 X(r) ~ x ~ x ~ 数字信号处理－离散傅立叶变换 * 3.4 DFT的应用举例 x(n)=cos n?/4 16点相当于取周期序列的两个周期进行DFT 数字信号处理－离散傅立叶变换 * 3.4 DFT的应用举例 若事先不知道 的周期，怎样进行频谱分析： 先截取x (n) M点，则求xM (n)的DFT： xM(n) = x (n)?RM(n)， XM(k) =DFT[xM(n)]，0?k ?M-1; 再截取x (n) 2M点，则求x2M(n)的DFT： x2M(n) = x (n)?R2M(n)， X2M(k)= DFT[x2M(n)]，0?k ?2M-1; 将2次截取序列的频谱进行分析，是否满足误差要求，若不满足，应加大截取窗长度（增加M值），再将结果进行分析。 ~ ~ ~ ~ x(n) ~ 数字信号处理－离散傅立叶变换 * 3.4 DFT的应用举例 3．用DFT进行谱分析的误差问题 DFT(实际中用FFT计算)可用来对连续信号和数字信号进行谱分析。在实际分析过程中，要对连续信号采样和截断，由此可能产生误差分析。 (1)频谱混迭现象： 原因: 不满足时域采样定理 避免措施：采样频率fs?2fc，以避免信号在w=?处附近的混迭。 具体方法是：采样时满足采样定理，采样前对信号进行预滤波，滤去信号中频率高于fs/2的频率分量。 数字信号处理－离散傅立叶变换 * 3.4 DFT的应用举例 (2)栅栏效应： 现象：N点DFT是在区间[0, 2?]上的N点等间隔采样，采样点之间的频谱函数值是不知道的，就好像从N+1个栅栏缝隙中观看信号的频谱特性，得到的是N个缝隙中看到的频谱函数值，这种现象称为栅栏效应。 原因：对信号的频谱进行有限点采样。 后果：栅栏效应可能漏掉(挡住)大的频谱分量 减少栅栏效应的措施：对原序列补0，增大N，以增加采样点； 数字信号处理－离散傅立叶变换 * 3.4 DFT的应用 数字信号处理－离散傅立叶变换 * 3.4 DFT的应用举例 (3)截断效应： 原因: 对序列x(n)截断所引起的。 无限长序列x(n)截短成有限长序列y(n)，即 y(n)=x(n)?RN(n)，则 Y(ejw)=FT[y(n)]=1/(2?)?X(ejw)*RN(ejw) =1/(2?)????? X(ej?)*RN(ej(w-?))?d?， 其中X(ejw)=FT[x(n)] RN(ejw)=FT[RN(n)]=e-jw (N-1)/2?sin(wN/2) ? sin(w/2)=RN(w)ej?(w) RN(w) w 0 2??N N 2??N 矩形窗函数幅度谱 主瓣 旁瓣 数字信号处理－离散傅立叶变换 * 例：x(n)=cos(w0n)，w0=?/4，用DFT分析其频谱特性。 解：序列的幅度谱X(ejw)= ?? [?(w-?/4-2?l)+ ?(w+?/4-2?l)] 加矩形窗截断后,Y(ejw)=1/2?X(ejw)*RN(ejw)，定性图如