[論文レビュー] Numerical computation of convolutions in free probability theory

Sheehan Olver, Raj Rao Nadakuditi|arXiv (Cornell University)|Mar 8, 2012

Random Matrices and Applications参考文献 26被引用数 15

ひとこと要約

本稿では、自由確率論における自由畳み込み（加法的、乗法的、圧縮的）を計算する数値的手法を提示する。アドミッサブルな測度の逆コーシー変換を活用することで、可逆な測度（例：半円分布、ガウス分布、平方根減衰型）に対してスペクトル精度を達成する。これは、良好に条件付けられたバナッハ型バナッハ行列系を用いることで実現され、入力測度に広い条件下で収束が保証される。

ABSTRACT

We develop a numerical approach for computing the additive, multiplicative and compressive convolution operations from free probability theory. We utilize the regularity properties of free convolution to identify (pairs of) `admissible' measures whose convolution results in a so-called `invertible measure' which is either a smoothly-decaying measure supported on the entire real line (such as the Gaussian) or square-root decaying measure supported on a compact interval (such as the semi-circle). This class of measures is important because these measures along with their Cauchy transforms can be accurately represented via a Fourier or Chebyshev series expansion, respectively. Thus, knowledge of the functional inverse of their Cauchy transform suffices for numerically recovering the invertible measure via a non-standard yet well-behaved Vandermonde system of equations. We describe explicit algorithms for computing the inverse Cauchy transform alluded to and recovering the associated measure with spectral accuracy. Convergence is guaranteed under broad assumptions on the input measures.

研究の動機と目的

• 自由確率論における自由畳み込み（加法的、乗法的、圧縮的）を計算するための数値フレームワークの構築を目的とする。
• 代数的でない測度に対して、記号的計算を超えた効率的な数値ツールの欠如に対処することを目的とする。
• 半円分布やガウス分布のような一般的な分布を含む自由畳み込みの高精度計算を可能とすることを目的とする。
• 逆コーシー変換と測度回復がスペクトル精度で計算可能となる条件を確立することを目的とする。
• 研究者や実務家が自由確率論の畳み込みを数値的に探索できる、実装可能な実用的アルゴリズムの提供を目的とする。

提案手法

• 本手法は、自由畳み込み演算を入力測度の逆コーシー変換の点ごとの計算に還元する。
• 逆コーシー変換を高精度に計算可能な「アドミッサブルな」測度（滑らかに減衰する、ヤコビ型、半平方根型、点質量、またはそれらの有限組み合わせ）を同定する。
• 「可逆な」測度（シュワーツ型、平方根減衰型、または半平方根／シュワーツ型）に対しては、非標準的だが良好に条件付けられたバナッハ型バナッハ行列系を用いて測度を回復する。
• 逆コーシー変換は数値的台形則と解析接続を用いて計算され、コーシー変換の正則性を活用する。
• 逆コーシー変換をフーリエまたはチェビシェフ級数展開を用いて表現することで、スペクトル精度を達成する。
• バナッハ型行列系の条件数の上限に基づき、台形則の点が一様に分布し、関数が解析的であることに依存して、収束が保証される。

実験結果

リサーチクエスチョン

• RQ1 非代数的確率測度に対して、自由畳み込み演算をスペクトル精度で数値的に計算可能か？
• RQ2 コンパクトまたは無限大に台を持つ測度に対して、逆コーシー変換を高精度に計算するにはどうすればよいか？
• RQ3 どのような測度のクラスにおいて、逆コーシー変換からの出力測度の安定的かつ高精度な回復が可能か？
• RQ4 測度回復に用いるバナッハ型行列系が良好に条件付けられるのはどのような条件下か？
• RQ5 本手法は、たとえばランダム行列の和や積の極限スペクトル分布を計算するといった、ランダム行列理論の実用的問題に応用可能か？

主な発見

• 本手法により、広範な測度クラスに対して自由加法的・乗法的・圧縮的畳み込みのスペクトル精度が達成される。
• アドミッサブルな測度（半円分布やガウス分布を含む）に対して、逆コーシー変換は点ごとの高精度計算が可能である。
• バナッハ型行列系による測度回復は、弱い仮定のもとで安定的かつ収束的であり、大まかな m に対して条件数が √n + δ で有界である。
• 入力関数が解析的で f(−1) = 0 を満たす場合、最小二乗近似の収束速度は超代数的になる（変形されたバナッハ型行列系において）。
• 数値例により、半円分布とガウス分布の自由加法的畳み込み、およびガウス分布の自由圧縮が、高精度で正しく計算されていることが示された。
• 解析関数のトレUNCTION誤差の減衰とモア・ペンローズ一般逆行列の上限に基づき、収束に関する理論的保証が確立された。