線性卷积的FFT算法及其matlab实现.docxVIP

/jpkc2/ipkc/signal/new/course/two/2_3_31.htm线性卷积的FFT算法 线性卷积是求离散系统响应的主要方法之一,许多重要应用都建立在这一理论基础上,如卷积滤波等。 以前曾讨论了用循环卷积计算线性卷积的方法归纳如下: 将长为N2的序列x(n)延长到L,补L-N2个零 将长为N1的序列h(n)延长到L,补L-N1个零 如果L≥N1+N2-1,则循环卷积与线性卷积相等,此时,可有FFT计算线性卷积,方法如下: a.计算X(k)=FFT[x(n)] b.求H(k)=FFT[h(n)] c.求Y(k)=H(k)Y(k) k=0～L-1 d.求y(n)=IFFT[Y(k)] n=0～L-1 可见，只要进行二次FFT,一次IFFT就可完成线性卷积计算。计算表明,L32时,上述计算线性卷积的方法比直接计算线卷积有明显的优越性,因此,也称上述循环卷积方法为快速卷积法。 上述结论适用于x(n),h(n)两序列长度比较接近或相等的情况,如果x(n),h(n)长度相差较多,例如,h(n)为某滤波器的单位脉冲响应,长度有限,用来处理一个很长的输入信号x(n),或者处理一个连续不断的信号，按上述方法,h(n)要补许多零再进行计算,计算量有很大的浪费，或者根本不能实现。为了保持快速卷积法的优越性,可将x(n)分为许多段后处