线性卷积matlab代码,线性卷积的FFT算法及其matlab实现

线性卷积的FFT算法

线性卷积是求离散系统响应的主要方法之一,许多重要应用都建立在这一理论基础上,如卷积滤波等。

以前曾讨论了用循环卷积计算线性卷积的方法归纳如下:

将长为N

2的序列x(n)延长到L,补L-N

2

个零

将长为N

1的序列h(n)延长到L,补L-N

1

个零

如果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就可完成线性卷积计算。计算表明,L>32时,上述计算线性卷积的方法比直接计算线卷积有明显的优越性,因此,也称上述循环卷积方法为快速卷积法。

上述结论适用于x(n),h(n)两序列长度比较接近或相等的情况,如果

x(n),h(n)长度相差较多,例如,h(n)为某滤波器的单位脉冲响应,长度有限,用来处理一个很长的输入信号x(n),或者处理一个连续不断的信号，按上述方法,h(n)要补许多零再进行计算,计算量有很大的浪费，或者根本不能实现。为了保持快速卷积法的优越性,可将x(n)分为许多段后处理,每小段的长与h(n)接近,其处理方法有两种：

(1)重叠相加法——由分段卷积的各段相加构成总的卷积输出