第3章离散傅里叶变换2.pptVIP

算法原理 设序列点数N=2M，M为整数。 将X(k)按k的奇偶分组前，先将输入x(n)按n的 顺序分成前后两半。 按k的奇偶将X(k)分成两部分： 令 则X(2r)和X(2r+1)分别是x1(n)和x2(n)的 N / 2点DFT，记为X1(k)和X2(k) x1(0) x1(1) -1 x1(2) x1(3) -1 x2(0) x2(1) -1 x2(2) x2(3) -1 N/2点 DFT N/2点 DFT x(0) x(7) x(1) x(2) x(3) x(4) x(5) x(6) X1(0)=X(0) X2(0)=X(1) X1(1)=X(2) X1(2)=X(4) X1(3)=X(6) X2(1)=X(3) X2(2)=X(5) X2(3)=X(7) N /2仍为偶数，进一步分解：N /2 N /4 x3(0) x3(1) -1 -1 x4(0) x4(1) N/4点 DFT N/4点 DFT x1(0) x1(1) x1(2) x1(3) X3(0)=X1(0)=X(0) X4(0)=X1(1)=X(2) X3(1)=X1(2)=X(4) X4(1)=X1(3)=X(6) 同理： 其中： 逐级分解，直到2点DFT 当N=8时，即分解到x3(n)，x4(n)，x5(n)，x6(n)，n=0,1 DIF与DIT比较1 相同之处 DIF与DIT两种算法均为原位运算。 DIF与DIT运算量相同。 运算量 不同之处 DIF与DIT两种算法结构倒过来 DIF为输入顺序，输出乱序。运算完毕再运行“二进制倒读”程序。 DIT为输入乱序，输出顺序。先运行“二进制倒读”程序，再进行求DFT。 蝶形结不同 DIT的复数相乘出现在减法之前。 DIF的复数相乘出现在减法之后。 DIF与DIT比较2 快速傅里叶逆变换(IFFT) 以上所讨论的FFT的运算方法同样可用于IDFT的运算，简称为IFFT。即快速付里叶反变换。从IDFT的定义出发，可以导出下列二种利用FFT来计算IFFT的方法。 利用FFT计算IFFT的思路1 把FFT的时间抽取法，用于IDFT运算时，由于输入变量由时间序列x(n)改成频率序列X(k)，原来按x(n)的奇、偶次序分组的时间抽取法FFT，现在就变成了按X(k)的奇偶次序抽取了。 同样，频率抽取的FFT运算用于IDFT运算时，也应改变为时间抽取的IFFT。 利用FFT计算IFFT的思路2 此为DFT可用FFT程序 利用FFT计算IFFT的思路3 共轭 FFT 共轭 乘1/ N 直接调用FFT子程序计算IFFT的方法： 快速傅里叶的应用 快速卷积 快速相关 FFT FFT 序列相乘 IFFT 快速卷积 重叠相加法 FFT FFT 序列相乘 IFFT 共轭 快速相关 fftfilt fft ifft 3.4 与本章内容有关的MATLAB函数 例 已知信号 ，求N点DFT的幅值谱和相位谱。 ％M文件如下： N=64; fs=100; dt=1/fs; n=0:N-1; x=2*sin(5*pi*n*dt)+5*cos(18*pi*n*dt); y=fft(x,N); mag=2*abs(y)/N; pha=angle(y); f=n*fs/N; subplot(121); plot(f,mag); title(Magnitude) subplot(122); plot(f,pha); title(Phase) 小 结 连续时间信号的傅立叶级数和傅立叶变换 序列的傅里叶变换(DTFT) 离散傅里叶变换(DFT) 及其性质 DFT在实际应用中存在的问题 快速算法FFT * 常系数线性差分方程 * 常系数线性差分方程 * 只须将频域成份一个求共轭变换，即首先将X(k)的虚部乘以-1，（先取X(k)的共轭）得X*(k)；然后将X*(k)直接送入FFT程序即可得出Nx*(n)；最后再对运算结果取一次共轭变换，并乘以常数1/N，即可以求出IFFT变换的x(n)的值。 3.3 快速傅里叶变换 FFT 并 不 是 一 种 新 的 变 换 形 式，它 只 是 DFT 的 一 种 快 速 算 法 。并 且 根 据 对 序 列 分 解 与