Matrix equation representation of the convolution equation and its unique solvability
Yuki Satake, Tomohiro Sogabe, Tomoya Kemmochi, Shao-Liang Zhang
TL;DR
This work shows that the convolution equation $F*X=B$ with a $3\times3$ kernel $F$ and $B\in\mathbb{R}^{m\times n}$ can be reformulated as a generalized Sylvester equation via shift matrices under zero, periodic, or reflexive boundary conditions. For common image-processing filters, the authors derive explicit necessary and sufficient conditions on $m$ and $n$ ensuring a unique solution for every right-hand side $B$, and they further reduce many cases to simpler Sylvester/Lyapunov forms. A comprehensive set of results is provided for box blur, Gaussian blur, edge detection (A, B, C), sharpen, and emboss, with several cases yielding unconditional solvability and others requiring parity constraints on $m$ and $n$. The findings pave the way for efficient solvers that exploit matrix structure, and they outline future work on numerical algorithms and the unresolved reflexive-emboss case, with potential impact on image restoration and related applications.
Abstract
We consider the convolution equation $F*X=B$, where $F\in\mathbb{R}^{3\times 3}$ and $B\in\mathbb{R}^{m\times n}$ are given, and $X\in\mathbb{R}^{m\times n}$ is to be determined. The convolution equation can be regarded as a linear system with a coefficient matrix of special structure. This fact has led to many studies including efficient numerical algorithms for solving the convolution equation. In this study, we show that the convolution equation can be represented as a generalized Sylvester equation. Furthermore, for some realistic examples arising from image processing, we show that the generalized Sylvester equation can be reduced to a simpler form, and analyze the unique solvability of the convolution equation.
