离散傅里叶变换及其快速算法1.pptVIP

根据x(n)与X(k)的对称性，同样可找到X(k)的实部、虚部与x(n)的共轭偶部与共轭奇部的关系。 分别以xe(n)及x0(n)表示序列x(n)的循环（折叠时用到循环折叠）共轭偶对称部分与循环共轭奇部： 同样应从圆周意义上理解 x(N-0)=x(0)，或者x*(N-n)。 可证明: DFT[xe(n)]=Re[X(k)] DFT[x0(n)]=jIm[X(k)] ? 第32页/共53页 例： 设N=8，定义在主值区间的序列x=[ones(1,6),zeros(1,2)],试将其分解为偶对称序列和奇对称序列，求它们的DFT并验证本节所述的特性。 N=8; x=[ones(1,6),zeros(1,2)]; %给定实序列 nx=0:N-1; xr =x(mod(-nx,N)+1); % 求x的循环折叠xr xe=0.5*(x+xr); xo=0.5*(x-xr); %序列分解 第33页/共53页 X=fft(x,N), Xe=fft(xe,N), Xo=fft(xo,N) % 求x,xe,xo的DFT DFT[xe(n)]=Re[X(k)] DFT[x0(n)]=jIm[X(k)] 第34页/共53页 实序列的DFT 第35页/共53页 （6）选频性 对复指数函数 进行采样得复序列 x(n) 0≤n≤N-1 其中q为整数。当ω0=2π/N时，x(n)=ej2πnq/N,其离散傅里叶变换为 写成闭解形式 可见，当输入频率为qω0时，变换X(k)的N个值中只有 X(q)=N，其余皆为零，如果输入信号为若干个不同频率的信号的组合，经离散傅里叶变换后，不同的k上，X(k)将有一一对应的输出，因此，离散傅里叶变换算法实质上对频率具有选择性。 第36页/共53页 （7）DFT与Z变换和DTFT 有限长序列可以进行z变换 比较z变换与DFT变换，可见，当 时， 即 是z平面单位圆上幅角为 的点，即将z平面上的单位圆N等分后的第k点。 1）X(k)也就是z变换在单位圆上等间隔的采样值。 2）X(k)也可看作是对序列付氏变换X(ejω)的采样，采样间隔为：ωN=2π/N。 即 结论： 第37页/共53页 图 DFT与z变换 o o o o o o o o o o o X(ejω) X(k) o Re[z] jIm[z] o 第38页/共53页 变量 周期 分辨率 采样定律告诉我们，一个频带有限的信号，可以对它进行时域采样而不丢失任何信息； DFT变换进一步告诉我们，对于时间有限的信号（有限长序列），也可以对其进行频域采样，而不丢失任何信息，这正反映了傅立叶变换中时域、频域的对称关系。它有十分重要的意义，由于时域上的采样，使我们能够采用数字技术来处理这些时域上的信号（序列），而DFT的理论不仅在时域，而且在频域也离散化，因此使得在频域采用数字技术处理成为可能。 FFT就是频域数字处理中最有成效的一例。 第39页/共53页 （8）DFT形式下的Parseval定理 利用循环卷积和共轭对称特性，可证明DFT形式下的Parseval定律： 第40页/共53页 各种傅里叶变换及其相互关系 时域的周期性对应于频域的离散化，时域的离散化对应于频域的周期性 第41页/共53页 2.2 利用DFT做连续信号的频谱分析 利用DFT计算连续信号的频谱 采样 截短 DFT 第42页/共53页 (1) 混叠 对连续信号x(t)进行数字处理前，要进行采样 采样序列的频谱是连续信号频谱的周期延拓，周期为fs，如采样率过低，不满足采样定理，fs2fh，则导致频谱混迭，使一个周期内的谱对原信号谱产生失真，无法恢复原信号，进一步的数字处理失去依据。 ? （2）? 泄漏