6第六篇_快速傅里叶变换.ppt

1234 借助于FFT算法的重叠相加法需要以下步骤： （1）计算 ； （2）计算 ； （3）计算 ； （4）计算 ； （5）将各子序列 重叠相加。 DFT与FFT的计算量 FFT的时间抽取蝶形图的特点和画法 倒位序的获得 IFFT的计算 FFT的编程 为何采用即位运算 快速卷积的获得 如何利用FFT进行频谱分析 请学习第七章 6.4 快速傅里叶反变换IFFT 另： 6.5 快速卷积 线性卷积 循环卷积 5.3 离散傅立叶变换的性质 利用循环卷积计算线性卷积 N2 N1 N1+N2-1 序列补零为N1+N2-1 得到新序列 计算新序列的循环卷积 上述结果即为线性卷积 5.3 离散傅立叶变换的性质 求：N点、2N点循环卷积及线性卷积 N点循环卷积 线性卷积 5.3 离散傅立叶变换的性质 2N点循环卷积 n=0 n=1 5.3 离散傅立叶变换的性质 利用2N点DFT MATLAB: convolution_5_3.m 6.5 快速卷积 利用FFT计算循环卷积, 然后得到线性卷积 序列补零为2M=N1+N2-1 得到新序列 计算新序列的FFT IFFT 将上述结果相乘 取前N1+N2-1位 6.5 快速卷积 线性卷积快速算法的计算量 1. 2. 2次FFT计算和1次IFFT 、一次 点复数乘法 补零后长度均为 直接计算线性卷积需要进行 次实数乘法运算 离散时间系统的 已知或已经事先求出，并已置于存储器中。实际中只需要2次 点FFT的运算量 1 2 N 需要的全部复数乘法运算次数为 6.5 快速卷积 重叠相加法 M N M+N-1 6.5 快速卷积 6.5 快速卷积 在实际中也可以边计算 边进行重叠相加。这样，可以保证信号处理具有最小的延时。重叠相加法已被广泛地应用于处理无限长信号的实时处理，例如，语音信号等。 6.6 FFT的应用 测量线性系统的系统函数 线性系统的零状态响应 电力系统谐波分析 A/D变换时，应该在整个周期T0内采样 kN/2 6.6 FFT的应用 电气工程，采用正弦系统 例1-3：求周期矩形脉冲信号的指数谱，并画出它的幅频特性和相频特性图形 在指数谱中，直流分量满足幅值包络线函数 1.2 周期信号的频谱：二、指数频谱 滤波 电力系统暂态分析中的应用 在电力系统测量中的应用 6.6 FFT的应用 由采样定理和FFT的基2原则，确定采样数 采样时间间隔T=0.25s f=4Hz 例: 试用FFT求其频谱 周期T0=1s, 频率f0=1Hz 6.6 FFT的应用 电子工程，采用指数系统 电气工程，采用正弦系统 离散信号频谱 6.6 FFT的应用 MATLAB: fft_sin_abs.m MATLAB: fft_sin_randn_Power_spectral_density.m MATLAB: fft_sin_stem_spectrum_overlap.m FFT在MATLAB中的信号分析应用 MATLAB: fft_sin_zero_padded.m 作业：6-1、6-2、6-3、6-4、6-5、6-6、6-7、 6-8 * 信号分析与处理 第六章 快速傅里叶变换 电气工程学院电力工程系 安 勃 电力系电信教研室 教一楼335室 752-3031 bo_an_cn@ 第六章 快速傅里叶变换 6.1 提高DFT运算速度的主要方法 6.2 时间抽取算法（N=2M） 10.2 FFT时间抽取算法的程序 6.3 FFT 的频率抽取算法 （N=2M） 6.4 快速傅里叶反变换IFFT 6.5 快速卷积 6.6 FFT的应用 次复数乘法和 次复数加法运算 6.1 提高DFT运算速度的主要方法 ＝ 6.1 提高DFT运算速度的主要方法 周期性 对称性 可以将N点DFT变为N/2点DFT，直至变为2点DFT 提高计算速度的基本原理 时间抽取算法 频率抽取算法 6.2 时间抽取算法（N=2M） 奇偶分解 偶数序列 奇数序列 时域奇偶分解 6.2 时间抽取算法（N=2M） 时域奇偶分解 k’=0,1,2,…,N/2-1 长度为N/2 k’=0,1,2,…,N/2-1 对应的频域顺序排列（前后分解） 6.2 时间抽取算法（N=2M） N