区块折积

区块折积，又称区块卷积、分段折积、分段卷积、分块卷积等，是一种计算线性离散卷积的方法，在信号处理以及线性非时变系统领域的应用层面上有相当高的价值。
区块卷积主要常使用于计算一固定离散信号与另一非固定离散讯号间的线性卷积，且非固定讯号的长度通常较固定方长上许多，一个很具体的例子就是计算一讯号经有限长度滤波器后的输出。
其主要具有两个优点，第一，其计算复杂度较低，在输入讯号长度为
formula_1时，区块卷积的计算复杂度为formula_2，而直接对整段输入讯号进行卷积计算的复杂度为formula_3。
第二，当我们要以硬体进行实作时，区块卷积仅需要单一固定的硬体架构，即可处理不同长度（甚至是无穷长）的输入，可以像工厂的流水生产线一般连续的接受输入并同时连续的吐出结果，
而如果要直接对整段输入讯号进行卷积，则需要对不同的输入长度设计不同的硬体运算架构，在应用上是不切实际并且没有必要的。
目前对于区块卷积，主流使用的英文专有词汇为sectioned convolution，但在一些中文圈作者所撰写的文章中，也可能会采用取「块」字翻英后的block convolution。
计算.
区块卷积的计算核心概念，便是将较长的讯号进行分段，将原先的求卷积问题分解为规模较小的问题，再将结果进行整合，以得到原先问题的结果。
但根据将讯号分段以及将结果整合的方式，可以再将区块卷积细分为两个类别，可见下文说明。
卷积的性质观察.
首先，在对区块卷积的计算法进行细分之前，我们需要先了解一些关于卷积的性质，以及关于切分后子问题的结果与母问题结果的关联。
假设有一长度为formula_1的离散讯号（通常很长）：
formula_5
与一长度为formula_6的离散讯号（通常很短）：
formula_7
将两者进行卷积后的结果formula_8，我们可以知道其长度应该为formula_9：
formula_10
根据卷积的定义，我们可以知道其中formula_11的第formula_12点输出formula_13，
会受到相对应位置的formula_14，以及其向前formula_15点的影响：
formula_16
这时如果我们将formula_17切出长度为formula_18的小段，并假设起始位置为formula_19：
formula_20
我们可以发现，其与formula_21进行卷积后的结果formula_22，和原先母问题在对应区段的结果formula_11存在差异，
具体来说，formula_24中formula_25部分的影响，例如在formula_26
就仅考量一项formula_17的影响，和formula_11相差许多：
formula_29
同理，在formula_30的后段部分也会有类似的问题，少考量了formula_17在段落之后的影响，而只有计算到部分残留影响的结果。
只有在formula_30的中段部分，会与母问题结果formula_11完全相同。
而要如何处理这个差异，并透过多段的formula_30还原出formula_11，便是以下两种方法的最主要区别。
重叠-储存之折积法.
重叠-储存算法，便是舍弃子结果formula_30的前后段，只取中段与母结果formula_11相同的部分并进行拼接的演算法。
具体来说，重叠-储存算法切分问题的主要考量为母结果formula_11，
在实作这个算法的时候必须先决定将母结果formula_11切分的长度。
假设决定将母结果切分为长度formula_40的区间，
根据前面的讨论，我们可以得知：
1.formula_41的长度需要为formula_42才能够确保formula_40长度的有效区段
2.在每次长度为formula_44的计算结果formula_30中，需要舍去前后各formula_46长度的部分
3.在第一个区段，因为formula_17在formula_48的区域是没有资讯的，所以formula_49的长度只需要formula_40即可
综合以上三点，我们可以写出所需要的切分区段描述：
formula_51
formula_52
以及结果的描述：
formula_53
而重叠-储存算法便是得名于其在输入区段的切分上会有重叠的部分，须将前一个子问题的输入后半先行储存再作为下一个子问题的输入前半。
这个算法的好处在于不需要额外的加法，并且如果在子问题求解时是采取直接计算卷积的方式，可以不用在讯号后面补零。
但缺点为同样的分段数量下，每个子问题会需要处理较长的讯号，子问题计算量可能会较大一些。
重叠-相加之折积法.
重叠-相加算法，则是将子结果formula_30的后段，
与下一段子结果formula_30的前段，
像拼图一样根据对应的编号进行相加，
重建出完整母结果formula_11的演算法。
具体来说，重叠-储存算法切分问题的主要考量为输入formula_17，
在实作这个算法的时候必须先决定将输入formula_17切分的长度。
假设决定将输入切分为长度formula_18的区间，
我们可以直接写出所需要的切分区段描述：
formula_60
以及结果的描述：
formula_61
这个算法特别需要注意的便是，重建出的母结果formula_11牵涉到子结果对应项目的相加，要正确地对上才会是正确的结果，例如formula_63就会是两个子结果的相加：
formula_64
而重叠-相加算法便是得名于其在输出的子结果上会有项目编号重叠的部分，须将其进行相加才能得到母结果。
这个算法的好处在于切分的方式相当直觉，而且在同样的分段数量下，每个子问题需要处理的讯号长度较短，子问题计算量较小。
但缺点为可能需要一些额外的加法，而且就算采取直接计算卷积的方式来计算子问题，也仍需要补零在讯号后方。
复杂度与比较.
虽然区块卷积具有两种差异不小的算法，但总体来说，对于每个固定架构的子问题，所需要的计算量都是常数，
而子问题的数量则正比于输入讯号长度formula_1，所以计算复杂度对两个算法来说都是formula_2，
相较于使用快速傅里叶变换技巧进行整个问题的计算复杂度formula_3要来的低。
所以一般来说在formula_1很大时，显然的使用区块卷积会较处理整个问题来的有效率。
不过，有时候输入讯号formula_17并不具有足够的长度可以保证区块卷积较有效率，
又或者在一些极端情况下，直接按定义进行计算卷积可能会更有效率，
所以接著便是要稍微谈论一下不同主流计算卷积方法的计算量，
以及如何根据输入讯号长度formula_1与固定讯号长度formula_6间的关系选择算法。
直接计算卷积.
如果是根据卷积的定义，直接以乘法及加法计算卷积的话，
我们可以知道，每一个formula_72都需要轮流和每一个formula_73进行相乘，
所以总共需要的乘法数量为
formula_74
又因为讯号可能是复数，所以上述的乘法为复数乘法，将其统一以时数乘法进行实作后可以得到需要的实数乘法数为
formula_75
可以发现，其实直接计算卷积对于输入长度formula_1来说，其复杂度也是formula_2，
但其常数为formula_78，当formula_6增大时，运算量会增加很快，运算效率会不及区块卷积，
但相对的当formula_6落在很小的值时，直接计算卷积可能会是较高效率的计算方法。
使用快速傅立叶变换技巧计算卷积.
除了直接计算卷积之外，另外一种主流的卷积计算方法，便是透过卷积定理，
将卷积的计算化为两次傅立叶变换相乘后再进行反傅立叶变换的问题，
实作如下：
formula_81
其中formula_82为formula_83点快速傅立叶变换，
formula_84为逆变换（形式与快速傅立叶变换相同），formula_83根据圆周卷积定理必须大于formula_86，
而formula_72及formula_73需透过补零来补齐长度。
如果采用这个方法，并且假设formula_21为固定不变，可以将formula_90计算一次后储存备用，
那可以得到所需乘法数为：
formula_91
其中formula_92为formula_83点快速傅立叶变换所需的乘法数，
估计为formula_94，
乘两倍是因为需计算formula_95与最后的逆变换，
formula_96则是因为formula_97含有formula_83个复数的乘法，统一转换为实数乘法则是formula_96。
将formula_100与formula_92的估计值带入上面的所需乘法数进行估计，可以得到所需乘法数为：
formula_102
可以发现，使用快速傅立叶变换技巧计算卷积对于输入长度formula_1来说，其复杂度是formula_3，
但当formula_6大约与formula_1落在接近的数量级时，运算量会较其他方法低上许多。
使用区块卷积搭配快速傅立叶算法计算卷积.
最后，另外一种主流的卷积计算法便是先利用区块卷积将问题拆分，再将拆分后的问题透过快速傅立叶算法来处理，
而之所以通常不是搭配直接计算的方式，
是因为我们通常会假设拆分过后的问题中，要进行卷积的两讯号长度会差不多长，
而根据前面的讨论，我们已经知道使用快速傅立叶转换来进行计算会有效率的多。
而使用区块卷积搭配快速傅立叶算法计算卷积的复杂度，
和前面两种方法的计算复杂度只与formula_107有关不同，
分段所使用的段落长度formula_108也会很大程度的影响这个方法的计算量。
以重叠-相加算法作为代表，
假设采用的分段长度为formula_108，
可以得到分段的段数formula_110为：
formula_111
而每个子问题的实际计算量根据上面可得：
formula_112
乘上需要解的子问题数可得实际总计算量：
formula_113
同样根据上面可再进一步估计得：
formula_114
可以发现这个计算法的复杂度为formula_2，
并且常数约为formula_116，
在formula_117的数值略大时会较采取直接计算会较有效率，
但在formula_117与formula_119的数量级接近时，
因为formula_117的大小被迫推升，
导致比起直接对整段讯号进行快速傅立叶算法的表现还要差。
所以根据上述可以估计这个方法可以展现优势的区段，
大概是界在前面两种方法的formula_107比例之间。
而在实际应用时，因为快速傅立叶转换的计算量在局部会有不小的波动，并不如理论值估计那般平滑，所以为了得到较好的计算效率，
除了须将上面的计算量视作formula_108的函数，透过数值绘图的方法或微分，求得在理论估计上最适合的formula_108值之外，
还需在求得的最适合子问题大小附近进行找寻，
挑选数个可能的合适子问题傅立叶转换大小，
再将反求得的分段长度带回计算运算量，
以求得真正的最佳分段长度formula_108。
例如：在
formula_125
的情况下，求得的理论最适合分段长度为
formula_126
理论最适合子问题傅立叶转换大小为
formula_127点
但在formula_128点傅立叶转换的计算量有一个明显的低谷，将由此逆推得到的可能分段长度
formula_129
与其他可能值都带回去计算实际计算量后，可以确认formula_130才是最佳的问题分段长度。
最后，同样是基于快速傅立叶转换的计算量在局部的波动，在极为特殊的情况下，如果formula_131，
我们便有可能使子问题仅需使用formula_132点傅立叶转换，而在这样的情况下虽会使用大量的加法，
但可不必在傅立叶转换上使用任何乘法，可能会较直接计算卷积有效率。
参考资料.
Jian-Jiun Ding, advanced digital signal processing lecture note , the Department of Electrical Engineering, National Taiwan University （NTU）, Taipei, Taiwan, 2022.