A geometric approach to approximating the limit set of eigenvalues for banded Toeplitz matrices
Fuente:
arXiv
Guardado en:
| Autores principales: | , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2023
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866912016590962688 |
|---|---|
| author | Bucht, Teodor Christiansen, Jacob S. |
| author_facet | Bucht, Teodor Christiansen, Jacob S. |
| contents | This article is about finding the limit set for banded Toeplitz matrices. Our main result is a new approach to approximate the limit set $Λ(b)$ where $b$ is the symbol of the banded Toeplitz matrix. The new approach is geometrical and based on the formula $Λ(b) = \cap_{ρ\in (0, \infty)} \text{sp } T(b_ρ)$, where $ρ$ is a scaling factor, i.e. $b_ρ(t) := b(ρt)$, and $\text{sp }(\cdot)$ denotes the spectrum. We show that the full intersection can be approximated by the intersection for a finite number of $ρ$'s, and that the intersection of polygon approximations for $\text{sp } T(b_ρ)$ yields an approximating polygon for $Λ(b)$ that converges to $Λ(b)$ in the Hausdorff metric. Further, we show that one can slightly expand the polygon approximations for $\text{sp } T(b_ρ)$ to ensure that they contain $\text{sp } T(b_ρ)$. Then, taking the intersection yields an approximating superset of $Λ(b)$ which converges to $Λ(b)$ in the Hausdorff metric, and is guaranteed to contain $Λ(b)$. Combining the established algebraic (root-finding) method with our approximating superset, we are able to give an explicit bound on the Hausdorff distance to the true limit set. We implement the algorithm in Python and test it. It performs on par to and better in some cases than existing algorithms. We argue, but do not prove, that the average time complexity of the algorithm is $O(n^2 + mn\log m)$, where $n$ is the number of $ρ$'s and $m$ is the number of vertices for the polygons approximating $\text{sp } T(b_ρ)$. Further, we argue that the distance from $Λ(b)$ to both the approximating polygon and the approximating superset decreases as $O(1/\sqrt{k})$ for most of $Λ(b)$, where $k$ is the number of elementary operations required by the algorithm. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2308_00829 |
| institution | arXiv |
| publishDate | 2023 |
| record_format | arxiv |
| spellingShingle | A geometric approach to approximating the limit set of eigenvalues for banded Toeplitz matrices Bucht, Teodor Christiansen, Jacob S. Numerical Analysis Spectral Theory 15B05, 47B06, 47B35, 65F15 This article is about finding the limit set for banded Toeplitz matrices. Our main result is a new approach to approximate the limit set $Λ(b)$ where $b$ is the symbol of the banded Toeplitz matrix. The new approach is geometrical and based on the formula $Λ(b) = \cap_{ρ\in (0, \infty)} \text{sp } T(b_ρ)$, where $ρ$ is a scaling factor, i.e. $b_ρ(t) := b(ρt)$, and $\text{sp }(\cdot)$ denotes the spectrum. We show that the full intersection can be approximated by the intersection for a finite number of $ρ$'s, and that the intersection of polygon approximations for $\text{sp } T(b_ρ)$ yields an approximating polygon for $Λ(b)$ that converges to $Λ(b)$ in the Hausdorff metric. Further, we show that one can slightly expand the polygon approximations for $\text{sp } T(b_ρ)$ to ensure that they contain $\text{sp } T(b_ρ)$. Then, taking the intersection yields an approximating superset of $Λ(b)$ which converges to $Λ(b)$ in the Hausdorff metric, and is guaranteed to contain $Λ(b)$. Combining the established algebraic (root-finding) method with our approximating superset, we are able to give an explicit bound on the Hausdorff distance to the true limit set. We implement the algorithm in Python and test it. It performs on par to and better in some cases than existing algorithms. We argue, but do not prove, that the average time complexity of the algorithm is $O(n^2 + mn\log m)$, where $n$ is the number of $ρ$'s and $m$ is the number of vertices for the polygons approximating $\text{sp } T(b_ρ)$. Further, we argue that the distance from $Λ(b)$ to both the approximating polygon and the approximating superset decreases as $O(1/\sqrt{k})$ for most of $Λ(b)$, where $k$ is the number of elementary operations required by the algorithm. |
| title | A geometric approach to approximating the limit set of eigenvalues for banded Toeplitz matrices |
| topic | Numerical Analysis Spectral Theory 15B05, 47B06, 47B35, 65F15 |
| url | https://arxiv.org/abs/2308.00829 |