Symmetric convolution

In order to compute symmetric convolution effectively, one must know which particular frequency domains (which are reachable by transforming real data through DSTs or DCTs) the inputs and outputs to the convolution can be and then tailor the symmetries of the transforms to the required symmetries of the convolution.

The following table documents which combinations of the domains from the main eight commonly used DST I-IV and DCT I-IV satisfy ${\displaystyle f*g=h}$ where ${\displaystyle *}$ represents the symmetric convolution operator. Convolution is a commutative operator, and so ${\displaystyle f}$ and ${\displaystyle g}$ are interchangeable.

Forward transforms of ${\displaystyle f}$, ${\displaystyle g}$ and ${\displaystyle h}$, through the transforms specified should allow the symmetric convolution to be computed as a pointwise multiplication, with any excess undefined frequency amplitudes set to zero. Possibilities for symmetric convolutions involving DSTs and DCTs V-VIII derived from the discrete Fourier transforms (DFTs) of odd logical order can be determined by adding four to each type in the above tables.

There are a number of advantages to computing symmetric convolutions in DSTs and DCTs in comparison with the more common circular convolution with the Fourier transform.

Most notably the implicit symmetry of the transforms involved is such that only data unable to be inferred through symmetry is required. For instance using a DCT-II, a symmetric signal need only have the positive half DCT-II transformed, since the frequency domain will implicitly construct the mirrored data comprising the other half. This enables larger convolution kernels to be used with the same cost as smaller kernels circularly convolved on the DFT. Also the boundary conditions implicit in DSTs and DCTs create edge effects that are often more in keeping with neighbouring data than the periodic effects introduced by using the Fourier transform.

• Martucci, S. A. (1994). "Symmetric convolution and the discrete sine and cosine transforms". IEEE Trans. Signal Process. SP-42 (5): 1038–1051. Bibcode:1994ITSP...42.1038M. doi:10.1109/78.295213.