DSP4_快速傅里叶变换综述.pptVIP

* 则： 式中： * * CZT 的实际计算方法： 1. 是 点系列，由 所决定： 2. 是双边无穷长序列，由定义所决定： * 3. 是 点序列，由需要所决定。 希望用DFT实现 截短 延拓 * * 点序列 （取前M个点） * 与本章有关的MATLAB文件主要是fft, ifft和 czt.m。fft实现快速傅立叶变换，ifft实现快速傅立叶反变换，czt.m实现线性调频Z变换 fft的调用格式是：X=fft(x), 或 X=fft(x，N) czt.m 调用格式是： X＝czt(x, M, W, A) 。x是待变换的时域信号，其长度设为N，M是变换的长度，W确定变换的步长，A确定变换的起点。若M=N, A=1, 则CZT变成DFT * 例 设模拟信号 ，以 t = 0.01n (n=0: N-1) 进行取样，试用fft函数对其做频谱分析。N分别为： (1) N=45；(2) N=50；(3) N=55；(2) N=60。 程序清单如下 %计算N=45的FFT并绘出其幅频曲线 N=45;n=0:N-1;t=0.01*n; q=n*2*pi/N; x=2*sin(4*pi*t)+5*cos(8*pi*t); y=fft(x,N); figure(1) subplot(2,2,1) plot(q,abs(y)) title(FFT N=45) * %计算N=50的FFT并绘出其幅频曲线 N=50;n=0:N-1;t=0.01*n; q=n*2*pi/N; x=2*sin(4*pi*t)+5*cos(8*pi*t); y=fft(x,N); figure(1) subplot(2,2,2) plot(q,abs(y)) title(FFT N=50) * %计算N=55的FFT并绘出其幅频曲线 N=55;n=0:N-1;t=0.01*n; q=n*2*pi/N; x=2*sin(4*pi*t)+5*cos(8*pi*t); y=fft(x,N); figure(1) subplot(2,2,3) plot(q,abs(y)) title(FFT N=55) * %计算N=60的FFT并绘出其幅频曲线 N=60;n=0:N-1;t=0.01*n; q=n*2*pi/N; x=2*sin(4*pi*t)+5*cos(8*pi*t); y=fft(x,N); figure(1) subplot(2,2,4) plot(q,abs(y)) title(FFT N=60) * 程序运行结果 从图中可以看出，这几种情况下均有较好的精度。 * 程序运行结果分析 分析：由t=0.01n进行取样可得，采样频率fs=100Hz。而连续信号的最高模拟角频率为Ω＝8 π ，由Ω＝2 πf可得，最高频率为8 π /2 π=4Hz。因此，满足采样定理的要求。 采样序列为 即 为周期序列，周期N=50。 将程序中plot改为stem函数，则可以更清楚地看出频谱。 * 修改程序运行结果 16乘，12加 →1乘，8加 * * * * * DIT与DIF的异同 1. 相同运算量，相同分级数、每级蝶形数， 同址运算特点 2. DIT输入码位倒置，DIF输出码位倒置 ——不是本质区别：输入输出可以重排 3. 蝶形单元不同：DIT先复乘后加减， DIF先减再复乘 —— 本质区别！！ 4. DIF、DIT蝶形单元互为易位（转置） （交换输入输出，所有支路反向） * 作业 画出DIF蝶形单元，写出迭代公式 * 第0级 第1级 第2级 4组 2组 1组 4.4 进一步减少运算量的措施 * FFT中乘法运算主要来自和复指数相乘： （1 组） （2 组） （4 组） （N/2 组） 复数乘法数 （N/4 组） 4.4.1 多类蝶形单元 * 不需要乘法： 平凡（trivial ）旋转因子 旋转因子（twiddle factor) * M 级，前两级都是 ，去除之： 后 M-2 级，含有 个 再去除之： （复乘） （ 如第2级共有 N/8 组， 每组2个平凡因子 ） * 两个复数相乘，需要四次实乘、两次实加。实现和 的相乘，需两次实乘，两次实加。 N点FFT中，有多少个 虚实部相等，trivial ~ 将所有平凡旋转因子去除，或单独考虑 个 第2级起，每组各有1 * 实乘 实加 各种算法比较的基础 以上称为多蝶形单元运算 计