[論文レビュー] Spectral features and asymptotic properties for alpha-circulants and alpha-Toeplitz sequences: theoretical results and examples

Eric Ngondiep, Stefano Serra‐Capizzano+1 more|ArXiv.org|Jun 11, 2009

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} を考慮し、テンソル積構造を用いることで、多段階およびブロックケースに結果を一般化する。
• 特に 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 の滑らかさに依存するが、本稿では詳細な定量的評価は行われていない。