重叠相加法（MATLAB实现）

---

一、重叠相加法原理

对于一个线性时不变系统来说，当给定其离散的输入序列时，要想计算该序列通过该线性时不变系统所引起的零状态响应时，可以使用该序列与该系统的冲激响应序列相卷积进行求解，如下式所示。
$$y(n)=x(n)\ast h(n)=\sum _{k=-\infty }^{\infty }x(k)h(n-k)$$
但当输入的数据序列无限长，而线性时不变系统的冲激响应序列比较短时，该方法会有很大的计算量。尤其是进行实时滤波，难以把所需要使用的输入序列都收集完毕后再进行运算，所以希望能够一边输入序列，一边进行计算。这种情况使用重叠相加法计算线性卷积，首先将输入序列x(n)拆分成许多的小区间， 然后用每一个小区间和冲激序列 h(n)进行卷积计算（卷积计算可以用FFT进行），最后将所得的每一个序列重叠部分进行叠加即可。算法可行性如下式所示：
$$y(n)=x(n)\ast h(n)=h(n)\ast \sum _{k=0}^{\infty }{x}_k(n)=\sum _{k=0}^{\infty }{y}_k(n)$$
示意图如下图所示：
​

二、算法流程

给出$x(n)$及$h(n)$，其中$h(n)$的长度为M。

①将$x(n)$进行分段得到$x_i(n)(i=0,1,2,3,\dots )$，每段长度为N，接着将$x_i(n)$补零延长至L=M+N-1，计算的L点FFT，得到$X_i(K)$。

②将$h(n)$补零延长至L=M+N-1得到$h_e(n)$，计算$h_e(n)$的L点FFT，得到$H_e(K)$。

③计算$Y_i(K)=X_i(K)\ast H_e(K)$。

④求出$Y_i(K)$的L点的IFFT，得到$y_i(n)$。

⑤将各段$y_i(n)$进行相加，但要考虑重叠部分：
$$y_i(M+n)=\{\begin{array}{ll}{y}_{i-1}(M+n)+y_i(n)& 0\le n\le N-1\\ y_i(n)& N-1\le n\le M-1\end{array}$$

三、Matlab仿真

1.h(n)和x(n)序列

$x(n)$为20000点序列，h(n)=[1,2,5,8]，长度为4，两个序列如下图所示。

2.使用重叠相加法计算x(n)与h(n)卷积的结果

3.使用Matlab自带的线性卷积函数的结果

为了验证重叠相加法的正确性，采用Matlab自带的conv函数来计算x(n)与h(n)的线性卷积，结果如下图所示。通过图形比较可以看出，两种方法得到的结果一样。

4.使用两种方法得到的结果的差值

为了更直观，将两种方法得到的结果进行相减，从下图可以看出，两者之差为0，因此重叠相加法的设计正确。

四、线性卷积与重叠相加法的比较

当h(n)的长度较短，而当x(n)序列的长度十分长接近无限长时，使用重叠相加法的运算量比使用线性卷积的运算量少，在实际运算中可考虑使用重叠相加法，但当x(n)序列的长度较短时，使用线性卷积的运算量比使用重叠相加法的运算量少，在实际运算中可考虑使用线性卷积，当然重叠相加法也可以稍作考虑。
（可在Matlab中使用tic和toc函数来比较两个程序的运行时间，从而比较两者运算量）

五、重叠相加法代码

本代码适用于将x(n)刚好分为整数段的情况，因此需设计好每一段长度。
代码链接：https://download.csdn.net/download/m0_46285064/79455268?spm=1001.2014.3001.5503

---

参考文章

重叠相加法计算卷积_鼻始拌饭加个蛋的博客-CSDN博客_重叠相加法计算线性卷积

数字信号处理——重叠相加法计算线性卷积_进修的码农的博客-CSDN博客_重叠相加法

重叠相加法&重叠保留法_373955482的博客-CSDN博客_重叠相加法

fft重叠加法_重叠相加法与重叠保留法_六堡茶之家的博客-CSDN博客