数字信号处理第5章 信号处理的效率.pptxVIP

第5章 信号处理的效率 数学公式的DFT还有经过计算才能吃完现实。DFT的计算方法优劣直接影响信号处理的速度。5.1 直接计算离散傅里叶变换简写为当n和k从0~N-1变化时，旋转因子WNkn 都在极坐标上绕单位圆旋转。X(k)和x(n)的计算形式相同。 5.1 直接计算5.1.1 直接计算频谱（1）不考虑旋转因子设旋转因子事先算好，并存在计算机存储器中。信号x(n)是N个复数的数组。计算全部频谱需复数乘N2次。计算全部频谱需复数加N(N-1)次。还有，每个k的X(k)都要用到全部x(n)；在算出全部X(k)前，x(n)要用2N个存储单元；X(k)要用2N个存储单元；整个计算过程至少需要4N个存储单元。5.1 直接计算（2）考虑旋转因子计算离散傅里叶变换需要旋转因子计算全部频谱需计算N×N个旋转因子。旋转因子要计算余弦和正弦。 从极坐标看，旋转因子是周期序列。利用周期性，计算离散傅里叶变换时，只需计算旋转因子的N个独立值。5.1 直接计算5.1.2 直接计算卷积设信号x(n)和系统h(n)的长等于N，则系统的输出它的长度2N-1。直接计算卷积需乘N(2N-1)次，加(N-1)(2N-1)次。 5.1 直接计算利用卷积定理也能计算y(n)，条件是循环卷积的长≥2N-1。当循环卷积的长=2N-1时，y(n)的频谱若先算好H(k)，用卷积定理求解y(n)的运算量是：5.1 直接计算复数乘4N(2N-1)次，复数加2(2N-1)(2N-2)次。 相比之下，直接卷积优于卷积定理。例5.1 设语音信号的采样率为6kHz，记录时间为1s，计算机复数乘1次需3μs，复数加1次需1μs。请问：该信号均分为6段，直接计算其频谱要多少时间？若分为3段，要多少时间？5.1 直接计算解 总样本为6000，信号分为6段，直接计算频谱的时间为若把信号分为3段，则直接计算频谱的时间为若信号样本逐秒连续输入，这两种算法不能实时分析。5.2 间接计算直接按定义计算离散傅里叶变换，工作量与N2成正比，还与旋转因子的独立值有关。这是一种启示：缩短DFT长度和减少旋转因子独立值，可以降低DFT的计算量。如果把N点DFT的长度缩短一半，变成两个N/2点DFT的组合，那么，DFT的复乘可从N2变成N2/2，复加可从N(N-1)≈N2变成约N2/2次。常用的分解有：1、按时序的奇偶数将序列分成两段，2、按时序的前后将序列分为两段。5.3 时域抽取的快速算法时域抽取的基本做法是，按时序的奇数和偶数分解序列成两段长度相同的子序列。这种算法要求序列的长度N满足5.3.1 时域抽取的原理基本时域抽取法分两步：第一，将序列分成两段；第二，整理频谱表达式。（1）分解序列 按时序n的偶奇数将N点序列x(n)分解成两个子序列，5.3 时域抽取的快速算法用x0(m)和x1(m)组成x(n)的N点DFT，5.3 时域抽取的快速算法（2）整理频谱为使X0(k)和X1(k)满足N/2点DFT的规定，同时又能反映X(k)的N个值，需对X(k)修改。当k=0~N/2-1时，当k=N/2~N-1时，r=0~N/2-1。5.3 时域抽取的快速算法利用旋转因子的周期和对称性，并将符号r换为k，得到 修改后的DFT基本分解式该式的好处：k值减半，DFT的乘加法都减半，旋转因子也减半。5.3 时域抽取的快速算法5.3.2 原理的推广对N/2点X0(k)和X1(k)继续分解。将X0(k)分解为两个N/4点DFT，同理，X1(k)分解为两个N/4点DFT，5.3 时域抽取的快速算法为了方便分解，将分解公式变为信号流图，称蝶形图。有更简洁的形式，一个碟形运算需复乘1次复加2次。5.3 时域抽取的快速算法例5.2 有一个8点离散傅里叶变换，请用蝶形图表示其时域抽取算法。解 遵循时域抽取法则，8点DFT可分解3次。5.3 时域抽取的快速算法蝶形运算有两个重要特点。（1）反序输入序列的时序等于1点长DFT的下标，下标用二进制表示。它按从左到右递增的二进制顺序，称反序。（2）原位运算蝶形的输入数据在后面的蝶形都不再出现，蝶形运算结果的位置和输入的位置相同。称为原位运算，它能够节省计算机的存储器。5.3 时域抽取的快速算法运用反序和原位运算，蝶形运算图的中间过程符号可略，变为这种时序的奇偶分解DFT方法称时域抽取基2快速傅里叶变换，简称时域抽取快速傅里叶变换5.3 时域抽取的快速算法5.3.3 时域抽取的计算量把每次分解DFT当作一级，时域抽取的蝶形运算共有M级，每级有N/2个蝶，每个蝶复数乘1次复数加2次。所以，时域抽取算法需复数乘MN/2次，复数加MN次。比直接计算法的N2 少许多。5.3 时域抽取的快速算法例5.3 设语音信号的采样率为6kHz，记录时间为1s，