数字信号处理办法西电版第三章快速傅里叶变换复习.pptVIP

例 用频谱分析下列时域信号的组成( fs=32 Hz ) 对序列作32点DFT，|X(k)|如右图示。频率分辨率 F=fs/N=1Hz 3.8 FFT的其它应用 3.8.1 线性卷积的FFT算法——快速卷积 1. FFT的快速卷积法（利用圆周卷积的FFT计算线性卷积） FFT计算y(n)的步骤 求H(k)=DFT［h(n)］ 求X(k)=DFT［x(n)］ 计算Y(k)=X(k)H(k); 求y(n)=IDFT［Y(k)］ 运算量的比较（直接线性卷积和FFT法计算线性卷积） x(n)与h(n)点数差不多时, 有如下结果(KM为两种卷积运算量比值) M=L 8 16 32 64 128 256 512 1024 2048 4096 Km 0.572 0.941 1.6 2.78 5.92 8.82 16 29.24 53.9 99.9 当M=L且M超过32以后，M越长， FFT法的优势越明显 x(n)与h(n)点数相差很大时 采用FFT运算时，要求长序列全部输入，短序列补充很多零点，这样既不经济，运算时间也长。解决的方法是将长序列分段计算，即重叠相加法与重叠保留法 设h(n)的点数为M，长序列x(n)可被分成为若干长度为L点的段 iL≤n≤(i+1)L-1 其他n i=0, 1, … 输入序列表示 2. 重叠相加法 x(n)和h(n)的线性卷积等于 计算N点FFT， H(k)=DFT［h(n)］ （ N=2m≥L+M-1 ） 计算N点FFT，Xi(k)=DFT［xi(n)］ （ N=2m≥L+M-1 ） 相乘，Yi(k)=Xi(k)H(k) 计算N点IFFT，yi(n)=IDFT［Yi(k)］ 计算步骤 将各段yi(n)（包括重叠部分）相加 xi(n)与yi(n)序列长度不同, 导致相邻两段输出有若干点重叠，应该将重叠部分相加后再和不重叠的部分共同组成输出y(n) 重叠相加法卷积图示 将x(n)分段，每段L=N-M+1个点 (h(n) 长度为M，做N点FFT) 保留第i个分段尾部的(M-1)点作为第(i+1)分段前面的(M-1)个输入值，组成L+M-1点序列xi＋1(n) 将卷积的结果yi+1(n)前面重叠部分(M-1)个点舍去后，剩余的按顺序排列就构成所需要的结果 为对分段卷积结果采取统一的舍去算法，第1段的前面可以补上(M-1)个零 3. 重叠保留法 计算步骤 重叠保留法图示 x(n)每段L=N-M+1个点 保留第i个分段尾部(M-1)点为xi＋1(n) 前(M-1)个输入值 重叠保留法图示 将yi+1(n)前(M-1)个点舍去，剩余的按顺序排列 FFT消噪是对含噪信号的频谱进行处理，即将噪声所在频段的X(k)值全部置零。这种方法要求噪声与信号的频谱不在同一频段 3.8.2 信号消噪 例 语音消噪 信号淹没在啸叫噪声中 信号与噪声的功率谱 去噪后的功率谱 重构的原语音信号 Dual-Tone Multi Frequency是按键式电话信令中的一般名称，应用于需要交互控制的系统中，如语音控制、ATM交换机等 3.8.3 FFT在双音多频（DTMF）信号中的应用 DTMF信令系统的工作原理 DTMF信令系统中共有八个频率，分为四个高频音和低频音，用高频音和低频音的组合来代表数字、字符*和# 按键拨号设计中八个频率的分配方案 按键拨号设计中双音检测方案 FFT算法实现双音检测 首先通过FFT计算DTMF信号的频谱，然后检测八个对应频率点的幅值来确定输入的信号 THE END 两种FFT运算方法等价 DIT-FFT DIF-FFT 3.4.3 IDFT的高效算法 FFT算法同样可以用于 IDFT，称为快速傅里叶反变换(IFFT) 比较DFT和IDFT的运算公式: 左图为在DIF-FFT流图上改动后，得到的DIT-IFFT运算流图 3.4.3 IDFT的高效算法 上述IFFT算法需要改动FFT的程序和参数才能实现。也可以利用DFT的性质，直接用FFT程序来完成IDFT的运算 已知 两边取共轭 两边再取共轭 3.5 N为复合数的FFT算法 3.6 线性调频Z变换（Chirp-Z变换）算法 FFT算法用于计算有限长序列的z变换X(z)在z平面单位圆上N个等间隔抽样点zk上的采样值 实际应用中在很多情况下并不一定需要计算全部频谱值,而仅需对某一频带内的信号频谱作较密集的分析 另外，采样也不一定局限于单位圆上，而需要计算出某一螺旋线上的等角度间隔的采样值 例如语音信号分析时，往往在靠近语音信号序列z变换的极点的螺旋线上进行采样，可以使语音信号的共振峰变得更尖锐，便于精确确定共振峰频率 3.6 线性调频Z变换