[논문 리뷰] Spectral features and asymptotic properties for alpha-circulants and alpha-Toeplitz sequences: theoretical results and examples

Eric Ngondiep, Stefano Serra‐Capizzano외 1명|ArXiv.org|2009. 06. 11.

Advanced Numerical Analysis Techniques참고 문헌 25인용 수 6

한 줄 요약

이 논문은 α-토플리츠 및 α-순환 행렬 수열의 특이값의 渐近 분포를 확립하며, 이들이 적분 가능 함수 f의 푸리에 계수에 의해 지배되는 스펙트럼 분포로 수렴함을 보여준다. α=1인 고전적 결과를 임의의 음이 아닌 정수 α로 일반화하여 다중 격자 방법과 웨이블릿 응용의 이론적 기초를 제공한다.

ABSTRACT

For a given nonnegative integer alpha, a matrix A_{n} of size n is called alpha-Toeplitz if its entries obey the rule A_{n}=[a_{r-alpha*s}]_{r,s=0}^{n-1}. Analogously, a matrix A_{n} again of size n is called alpha-circulant if A_{n}= [a_{(r-alpha*s)mod n}]_{r,s=0}^{n-1}. Such kind of matrices arises in wavelet analysis, subdivision algorithms and more generally when dealing with multigrid/multilevel methods for structured matrices and approximations of boundary value problems. In this paper we study the singular values of alpha-circulants and we provide an asymptotic analysis of the distribution results for the singular values of alpha-Toeplitz sequences in the case where {a_{k}} can be interpreted as the sequence of Fourier coeffcients of an integrable function f over the domain (-pi;pi). Some generalizations to the block, multilevel case, amounting to choose f multivariate and matrix valued, are briefly considered.

연구 동기 및 목표

• 임의의 음이 아닌 정수 α에 대해 α-순환 행렬의 특이값 분포를 특성화하는 것.
• 고전적 α=1의 경우를 초월해 α-토플리츠 수열의 스펙트럼 분포 이론을 확장하는 것.
• 생성 수열 {a_k}가 (−π, π)에서 적분 가능한 함수 f의 푸리에 계수에서 유래할 경우 특이값의 渐近 결과를 확립하는 것.
• 행렬 값 및 다변수 기호를 고려해 다수준 및 블록 구조화된 경우로 분석을 일반화하는 것.
• 구조화된 선형 방정식계에 대한 다중 격자 방법의 수렴 분석과 이론적 결과를 연결하는 것.

제안 방법

• α-토플리츠 행렬을 a_{r−αs}로, α-순환 행렬을 a_{(r−αs) mod n}으로 정의하여 표준 토플리츠 및 순환 구조를 일반화한다.
• 특이값 분해(SVD)를 사용해 행렬의 특이값을 B*B의 고유값과 연결함으로써 스펙트럼 분석을 가능하게 한다.
• 웨일의 스펙트럼 분포 이론을 적용하여 α-토플리츠 수열의 특이값이 기호 f에 의해 결정되는 |f(θ)|의 분포를 渐近적으로 따름을 보인다.
• 특히 DFT 행렬을 통한 대각화를 통해 DFT 기법과 푸리에 분석을 활용해 α-순환 행렬을 분석한다.
• 다변수 기호 f: (−π, π)^d → ℂ^{p×q}를 고려하고 텐서 곱 구조를 사용함으로써 다수준 및 블록 케이스로 결과를 일반화한다.
• 크로네cker 곱 및 텐서 곱 항등식을 활용해 블록 구조의 α-순환 행렬의 곱의 합에 대한 고유값 표현을 유도하며, 특히 J 행렬과 대각 블록을 다룬다.

실험 결과

연구 질문

• RQ1 고정된 α에 대해 n → ∞일 때 α-순환 행렬의 특이값은 어떻게 渐近적으로 행동하는가?
• RQ2 생성 수열 {a_k}가 (−π, π)에서 적분 가능한 함수 f의 푸리에 계수 수열일 경우 α-토플리츠 수열의 특이값은 어떻게 渐近적으로 분포하는가?
• RQ3 다수준 및 블록 구조화된 설정으로 일반화되었을 때 α-토플리츠 및 α-순환 행렬의 스펙트럼 성질은 어떻게 되는가?
• RQ4 α-토플리츠 행렬의 스펙트럼 분포와 다중 격자 방법의 수렴 행동 사이에 어떤 연결 고리가 존재하는가?
• RQ5 DFT 및 고유값 분해 기법을 사용해 α-순환 행렬의 특이값 구조를 완전히 특성화할 수 있는가?

주요 결과

• α-순환 행렬의 특이값은 생성 수열의 푸리에 변환의 절댓값에 따라 渐近적으로 분포하며, DFT와 대각화를 통해 명시적인 고유값 표현이 도출된다.
• α-토플리츠 수열의 경우 특이값은 기호 |f(θ)|를 가진 웨일 스펙트럼 분포를 따르며, f는 (−π, π)에서 적분 가능한 생성 함수이다.
• α-토플리츠 행렬의 특이값 분포는 α에 대해 극한에서 독립적이며 오직 기호 f에만 의존하며, 표준 L1 적분 가능성 조건 하에서 성립한다.
• 블록 구조의 α-순환 행렬에 대해 ∑_{j} C_j^* C_j의 고유값은 n_d,α를 기준으로 군집된 제곱 DFT 계수의 합에 비례하며, n_d,α를 초과하면 질서가 떨어진 구조를 가진다.
• 유도된 스펙트럼 분포 결과는 다수준 및 블록 케이스로 자연스럽게 일반화되며, 기호 f가 다변수 및 행렬 값이 되어도 渐近적 행동을 유지한다.
• 분석 결과는 α > 1에 대해서도 α-토플리츠 행렬의 특이값 분포가 안정함을 드러내며, 수렴 속도는 f의 미세함에 의해 결정되나, 본 논문에서는 이에 대한 정량적 분석은 상세히 제공하지 않는다.