---
title: Toeplitz-to-Circulant Aliasing Identity
url: https://www.emergentmind.com/topics/toeplitz-to-circulant-aliasing-identity
type: topic
---

# Toeplitz-to-Circulant Aliasing Identity

The Toeplitz-to-Circulant Aliasing Identity is the family of exact and asymptotic relations that convert the nonperiodic, linear-convolution structure of a Toeplitz matrix into periodic or anti-periodic convolution encoded by circulant-type matrices. In its most basic form, Toeplitz coefficients are folded modulo the matrix size, so wrap-around images such as $t_{k+mn}$ or $a_{i-j+\ell N}$ become circulant coefficients; in exact finite-dimensional formulations, the wrap-around contribution is either isolated as a skew-circulant or Hankel-like term, or eliminated by embedding into a larger circulant. These identities are central to FFT diagonalization, circulant and block-circulant preconditioning, CSCS iterations for Sylvester equations, multigrid coarse-grid symbol formation, QFT block-encoding on bounded domains, and several operator-theoretic decompositions used in random matrix theory [2108.03312], [1608.04820], [2605.16749].

## 1. Exact finite-dimensional identities

For a Toeplitz matrix $A \in \mathbb{C}^{n\times n}$ with entries $A_{j,k}=a_{j-k}$, one exact formulation is the circulant-and-skew-circulant splitting
$$
A = A_c + A_s.
$$
In the CSCS framework, the splitting is implemented coefficient-wise through the first column and first row. If $c=(a_0,a_1,\dots,a_{n-1})^T$ is the first column and $r=(a_0,a_{-1},\dots,a_{-(n-1)})$ the first row, then Algorithm 2 constructs
$$
\mathrm{CAc} = \frac{c + [0,\operatorname{reverse}(r(2\!:\!n))]}{2}, \qquad
\mathrm{CAr} = [r(1)/2,\operatorname{reverse}(\mathrm{CAc}(2\!:\!n))],
$$
and
$$
\mathrm{SAc} = \frac{c - [0,\operatorname{reverse}(r(2\!:\!n))]}{2}, \qquad
\mathrm{SAr} = [r(1)/2,-\operatorname{reverse}(\mathrm{SAc}(2\!:\!n))].
$$
Then $A_c=\operatorname{toeplitz}(\mathrm{CAc},\mathrm{CAr})$ and $A_s=\operatorname{toeplitz}(\mathrm{SAc},\mathrm{SAr})$, giving $A=A_c+A_s$ exactly. The “reverse-and-pad-0” operation is the finite-dimensional implementation of the wrap-around indices $j-k\pm n$ [2108.03312].

This exactness is not confined to the CSCS notation. A separate formulation writes a Toeplitz matrix $T_n(t)$ as
$$
T_n(t)=\frac{C_n(c)+S_n(s)}{\sqrt{2}},
$$
where $C_n(c)$ is circulant and $S_n(s)$ is skew-circulant, provided the sequences satisfy
$$
t_k=\frac{c_k+s_k}{\sqrt{2}}, \quad 0\le k\le n-1,
\qquad
t_{-k}=\frac{c_{n-k}-s_{n-k}}{\sqrt{2}}, \quad 1\le k\le n-1.
$$
Entrywise, positive and negative Toeplitz lags are recovered by modulo-$n$ indexing in the circulant term plus the sign-corrected wrap-around in the skew-circulant term [2605.16160].

A third exact form is the larger circulant embedding. For Hermitian Toeplitz $H_N$, if $S=[I_N;0_{N\times N}]$ and $C_{2N}(h^{ext})$ is the $2N\times 2N$ circulant built from
$$
c^{ext}[0]=h[0], \quad c^{ext}[k]=h[-k]\ (1\le k\le N-1), \quad c^{ext}[N]=0, \quad c^{ext}[N+k]=h[k]\ (1\le k\le N-1),
$$
then
$$
H_N x = S^T C_{2N}(h^{ext}) Sx.
$$
Here there is no aliasing correction term because the embedding size is large enough to reproduce linear convolution exactly [1608.04820].

## 2. Spectral diagonalization and convolutional meaning

The operational force of the identity comes from the fact that circulant and skew-circulant components are diagonalizable by Fourier transforms. For the unitary DFT matrix $F$ of order $n$,
$$
A_c = F^* \Lambda F, \qquad
\Lambda = \operatorname{diag}(\lambda_0,\dots,\lambda_{n-1}), \qquad
\lambda_k=\sum_{m=0}^{n-1} c_m \omega_n^{mk},
$$
where $(c_0,\dots,c_{n-1})$ is the first column of $A_c$. For the skew-circulant part, if
$$
D=\operatorname{diag}(1,e^{\pi i/n},e^{2\pi i/n},\dots,e^{(n-1)\pi i/n}), \qquad \Phi=FD,
$$
then
$$
A_s=\Phi^* E \Phi, \qquad
E=\operatorname{diag}(\theta_0,\dots,\theta_{n-1}),
$$
with eigenvalues computable by FFTs. The paper emphasizes that both $\Lambda$ and $E$ are computable in $O(n\log n)$ [2108.03312].

At the operator level, this is a decomposition of linear convolution into periodic and anti-periodic convolution. Toeplitz multiplication $y=Ax$ is linear convolution truncated by boundaries. Periodizing the domain introduces wrap-around terms, and the CSCS split partitions them into a symmetric periodic part, carried by $A_c x$, and an antisymmetric anti-periodic part, carried by $A_s x$. In the coefficient formulas these wrap terms appear explicitly as $a_{j-k\pm n}$; the identity therefore isolates finite-$n$ boundary effects into two FFT-diagonalizable pieces rather than discarding them.

The same viewpoint clarifies why Fourier methods work so broadly for Toeplitz problems and also why they can fail if the wrong geometry is imposed. A DFT diagonalizes circulants because its basis is periodic; the skew-circulant phase shift implements the corresponding anti-periodic twist. Toeplitz structure is therefore not “approximately Fourier diagonal” in a universal sense: the exact statement is that Toeplitz structure can be recoded into circulant-type pieces whose wrap-around images are explicit and algebraically controlled.

## 3. Weighted folding, optimal circulants, and asymptotic equivalence

A second major line of work treats aliasing not as an exact decomposition into two pieces, but as a one-piece circulant surrogate. In its simplest form, the first column of the $n\times n$ circulant induced by a Toeplitz sequence $\{t_k\}_{k\in\mathbb{Z}}$ is
$$
c_k = \sum_{m\in\mathbb{Z}} t_{k+mn}, \qquad k=0,1,\dots,n-1.
$$
This is the canonical folding of linear-convolution coefficients modulo $n$. In the frequency domain, the DFT of $\{c_k\}$ produces circulant eigenvalues that approximate symbol samples, with aliasing error controlled by the tails $\sum_{|k|\ge n}|t_k|$. The same mechanism extends to block Toeplitz with Toeplitz blocks, where the BCCB embedding uses the two-dimensional wrap-around
$$
c_{k,\ell}=\sum_{p,q\in\mathbb{Z}} t_{k+pn_1,\ell+qn_2}
$$
[1609.00573].

Several standard circulant constructions are weighted variants of this folding. Strang’s preconditioner keeps the nearest Toeplitz diagonals and performs a single wrap. Chan’s Frobenius-optimal circulant uses
$$
c_k^{(C)}=\frac{n-k}{n}t_k+\frac{k}{n}t_{k-n}, \qquad k=0,1,\dots,n-1.
$$
For real symmetric Toeplitz matrices, the nearest circulant in the Frobenius norm is likewise symmetric, with
$$
c_r=\frac{(n-r)t_r + r\,t_{n-r}}{n}, \qquad r=1,\dots,n-1, \qquad c_0=t_0.
$$
Its eigenvalues satisfy
$$
\lambda_m(C_n)=\sum_{d=-(n-1)}^{n-1}\left(1-\frac{|d|}{n}\right)t_d e^{-2\pi i m d/n},
$$
so the spectrum is the Fejér-kernel smoothing of the Toeplitz symbol sampled on the $n$-point grid. The same Fejér-averaged aliasing formula appears in the optimal circulant projection used for preconditioning functions of Toeplitz matrices [2208.05771], [1807.10929].

Asymptotic equivalence results show how far these surrogates can be pushed. For uniformly absolutely bounded Hermitian Toeplitz sequences $\{H_N\}$ and derived circulants $\{C_N\}$, square summability of the Toeplitz coefficients implies
$$
\frac{1}{N}\|H_N-C_N\|_F^2 \to 0,
$$
hence $\|H_N-C_N\|_F=o(\sqrt{N})$. This yields equal distribution of eigenvalues in the collective sense. Individual eigenvalue convergence is stronger: it holds for the three standard derived circulants when $h\in\ell^1(\mathbb{Z})$, and for the Cesàro circulant under weaker assumptions when the symbol is bounded, Riemann integrable, and has connected essential range. For band Toeplitz matrices the individual error is $O(1/N)$, and the largest and smallest eigenvalues converge to the essential supremum and infimum of the symbol [1608.04820].

A frequent misunderstanding is to treat these statements as interchangeable. They are not. Exact finite-$n$ identities, Frobenius-optimal projections, and asymptotic equivalence belong to different levels of the theory; the existence of a small Frobenius error does not by itself imply individual eigenvalue convergence.

## 4. Numerical linear algebra: Sylvester solvers, preconditioners, and multigrid

The CSCS iterative method for the Sylvester equation
$$
AX + XB = C
$$
assumes $A$ and $B$ are Toeplitz and splits them as
$$
A=C_A+S_A, \qquad B=C_B+S_B.
$$
With positive shift parameters $\alpha,\beta>0$, the method performs two substeps per iteration:
$$
(\alpha I + C_A)Z + Z(\beta I + C_B)=R,
$$
followed by
$$
(\alpha I + S_A)Z + Z(\beta I + S_B)=R'.
$$
Because $C_A,C_B,S_A,S_B$ are diagonalizable by FFTs or phase-shifted FFTs, each Sylvester subproblem reduces to diagonal linear systems. The iteration matrix for the Kronecker-sum formulation is
$$
M_y = (yI+\widehat S)^{-1}(yI-\widehat C)(yI+\widehat C)^{-1}(yI-\widehat S), \qquad y=\alpha+\beta,
$$
with spectral radius bounded by
$$
\rho(M_y)\le \theta_y
= \max_{\lambda\in X(\widehat C)}\left|\frac{y-\lambda}{y+\lambda}\right|
  \max_{\mu\in X(\widehat S)}\left|\frac{y-\mu}{y+\mu}\right|.
$$
If one of $C_A,C_B,S_A,S_B$ is positive definite and the remaining three are positive semidefinite, the CSCS iteration converges. Theorem 4.3 also gives a closed-form parameter choice $y^\*$ minimizing the bound, with the natural practical choice $\alpha=\beta=y^\*/2$ when $A$ and $B$ have comparable norms. The implementation uses spectral residual computation to avoid explicit dense multiplications, giving $O(n\log n)$ complexity per FFT per column [2108.03312].

In ill-posed Toeplitz and BTTB systems, aliasing is used more selectively. Dykes, Noschese, and Reichel build the Frobenius-optimal Chan circulant and then replace the $n-p$ smallest-in-magnitude eigenvalues by $1$,
$$
C_p = U\,\operatorname{diag}(\lambda_1,\dots,\lambda_p,1,\dots,1)\,U^*,
$$
so that the invariant subspace associated with the smallest eigenvalues is almost preserved rather than aggressively preconditioned. Their perturbation analysis leads to the criterion
$$
\min_{1\le q<n}\ \frac{1}{|\lambda_q|}\left(\frac{|\lambda_{q+1}|}{|\lambda_1|}+\eta\right),
\qquad \eta=\epsilon/\|b\|,
$$
and then sets $p=\lfloor(3/4)q\rfloor$. The same paper uses the filtered initial approximation
$$
\widetilde C_p^\dagger = U\,\operatorname{diag}(\lambda_1^{-1},\dots,\lambda_p^{-1},0,\dots,0)\,U^*,
\qquad x_0=\widetilde C_p^\dagger b,
$$
and stops the Krylov process by the discrepancy principle. On the reported image restoration and 1D “gravity” tests, the circulant or BCCB preconditioner roughly halves the iteration count while keeping restoration quality comparable to unpreconditioned runs [1609.00573].

In multigrid, aliasing governs coarse-grid symbol formation. For a circulant $A_n=C_n(f)$, a projector $p_{n,g}^k=C_n(p)Z_{n,g}^k$ yields a coarse operator whose symbol is
$$
\widehat f(x)=\frac{1}{g}\sum_{y\in\Omega_g(x/g)} f(y)|p(y)|^2.
$$
Equivalently, if $f|p|^2(\theta)=\sum_{\ell\in\mathbb{Z}} a_\ell e^{i\ell\theta}$, then $\widehat f(\theta)=\sum_{j\in\mathbb{Z}} a_{gj}e^{ij\theta}$. The $g=2$ case suffers from “mirror point” pathology when zeros of $f$ occur at points offset by $\pi$; the paper shows that choosing $g\ge 3$ and designing $p$ to vanish on the mirror set removes this obstruction. With $1<\theta<g$, the resulting multigrid cycle is optimal in the paper’s cost model [1010.5730]. This suggests that Toeplitz-to-circulant aliasing is not merely a device for preconditioning; it is also the algebraic mechanism by which coarse-grid correction is defined.

## 5. Boundary-aware operator compilation and matrix functions

A recent quantum formulation makes the boundary issue explicit. For the semi-discrete fractional Laplacian on a bounded interval with open, zero-extension boundary conditions, the physical operator $A_{\alpha,h}^{(N)}$ is a Toeplitz truncation of an infinite-lattice convolution kernel $a_m$. An $N$-point QFT, however, diagonalizes the circulant operator
$$
\widetilde A_{\alpha,h}^{(N)} = QFT_N^\dagger\, \operatorname{diag}(|\bar\xi_k|^\alpha)\, QFT_N,
$$
not the Toeplitz truncation. The exact aliasing proposition is
$$
(\widetilde A_{\alpha,h}^{(N)})_{ij} = \sum_{\ell\in\mathbb{Z}} a_{i-j+\ell N},
$$
so
$$
\widetilde A_{\alpha,h}^{(N)} = A_{\alpha,h}^{(N)} + E^{(N)},
\qquad
E^{(N)}_{ij}=\sum_{\ell\in\mathbb{Z}\setminus\{0\}} a_{i-j+\ell N}.
$$
Zero-padding to size $M\ge 2N$ and compressing back gives
$$
P_{N\to M}^\dagger \widetilde A_{\alpha,h}^{(M)} P_{N\to M}
= A_{\alpha,h}^{(N)} + E^{(M)},
$$
with
$$
\|E^{(M)}\|_2 \le \sum_{|r|\ge M-N+1}|a_r|
$$
by the Schur test. For the classical second-difference Laplacian, the kernel is finitely supported and the aliasing reduces to corner wraps only; padding to any $M\ge N+2$ eliminates them exactly on the compressed subspace. One call to the resulting block-encoding has cost
$$
G_{BE}=O(m^2)+C_\alpha(m,b)+C_{rot}(b),
\qquad m=\log_2 M
$$
[2605.16749].

Aliasing also underlies preconditioners for functions of Toeplitz matrices. For a Hermitian Toeplitz matrix $A_n$ generated by a positive Wiener-class symbol $f$, the optimal circulant $c(A_n)$ has eigenvalues given by Fejér-aliased samples of $f$, while the superoptimal circulant $T_n$ has eigenvalues
$$
\lambda_i(T_n)=\frac{(U_n A_n^2 U_n^*)_{ii}}{(U_n A_n U_n^*)_{ii}},
$$
that is, gridwise ratios of aliased samples of $f^2$ and $f$. For analytic $h$, the paper proves low-rank-plus-small-norm decompositions of $h(T_n)-h(A_n)$ and shows that the absolute value superoptimal preconditioner $|h(T_n)|$ yields spectra of $|h(T_n)|^{-1}h(A_n)$ clustered around $\pm 1$. The same construction extends to BTTB matrices with optimal BCCB preconditioners, although the multilevel case has $O(n)+O(m)$ outliers rather than a fixed number [1807.10929].

## 6. Parity splittings, random matrices, and conceptual scope

A particularly transparent structured example is Driessel’s near-Toeplitz tridiagonal matrix
$$
R_n = (\pi_n-\pi_n^T)E_+ + (\eta_n-\eta_n^T)E_-,
$$
where $\pi_n$ is the basic circulant, $\eta_n$ the basic skew-circulant, and $E_\pm=\tfrac12(I\pm E_n)$ are the projections onto even and odd subspaces with respect to the exchange matrix $E_n$. On $\operatorname{Ran}(E_+)$, $R_n$ acts exactly as the periodic first-difference operator $\pi_n-\pi_n^T$; on $\operatorname{Ran}(E_-)$, it acts as the anti-periodic first-difference operator $\eta_n-\eta_n^T$. The discrepancy between Toeplitz ends and wrap-around corners vanishes on the corresponding parity subspace [1102.1953].

The same circulant/skew-circulant algebra has recently been used in random matrix theory. For non-symmetric Toeplitz matrices, the identity
$$
T_n = \frac{C_n + S_n}{\sqrt{2}}
$$
expresses Toeplitz structure as a sum of a circulant and a skew-circulant. For Hankel matrices,
$$
H_n = \frac{R_n + L_n}{\sqrt{2}},
$$
with $R_n$ reverse circulant and $L_n$ left skew-circulant. Since the circulant-type components have tractable $*$-limits, the paper derives $*$-convergence results for random Toeplitz and Hankel matrices from these deterministic identities. In particular, symmetric Toeplitz matrices converge in $*$-distribution to the sum of two self-adjoint Gaussian variables, while random Hankel matrices converge to the sum of two self-adjoint symmetrized Rayleigh variables [2605.16160].

Three distinctions delimit the scope of the subject. First, exact finite-$n$ decompositions such as $A=A_c+A_s$ or $T_n=(C_n+S_n)/\sqrt2$ are not the same as optimal one-piece circulant approximations or asymptotic equivalence theorems [1608.04820]. Second, a finite DFT or QFT diagonalizes a circulant surrogate by construction; it does not diagonalize an open-boundary Toeplitz truncation unless an additional boundary adapter, such as zero-padding and compression, is introduced [2605.16749]. Third, in inverse problems the smallest-eigenvalue subspace is often the noise-sensitive part of the model, so an effective circulant preconditioner may need to preserve that subspace rather than flatten it [1609.00573].

Taken together, these results define the Toeplitz-to-Circulant Aliasing Identity as a unifying operator principle: linear-convolution structure can be folded into periodic or anti-periodic convolution, either exactly or asymptotically, provided the wrap-around images are represented explicitly. The resulting diagonalizability by FFT, phase-shifted FFT, or QFT is not an incidental computational convenience; it is the direct algebraic consequence of how Toeplitz boundary data are aliased into circulant geometry.

Source: https://www.emergentmind.com/topics/toeplitz-to-circulant-aliasing-identity