[論文レビュー] Inner product quadratures

Yu Chen|arXiv (Cornell University)|May 3, 2012

Mathematical Analysis and Transform Methods参考文献 1被引用数 3

ひとこと要約

本稿では、n個の関数の内積を統合するための新規n項四則積分法、すなわち内積四則積分法を提案する。これは、関数の統合ではなく内積統合を直接標的としているため、古典的ガウス四則積分法を一般化するものである。基底関数の直交列と四則積分ノードの間に双対性を確立し、商行列の固有値分解を用いてグラミアン行列の正確な統合を可能にする。画像処理、積分方程式、および非定値重みを伴う高次元四則積分設計への応用に応用可能である。

ABSTRACT

We introduce a n-term quadrature to integrate inner products of n functions, as opposed to a Gaussian quadrature to integrate 2n functions. We will characterize and provide computational tools to construct the inner product quadrature, and establish its connection to the Gaussian quadrature.

研究の動機と目的

• 古典的ガウス四則積分法が2n個の関数を統合するのではなく、n個の関数の内積を直接統合する四則積分法を開発すること。
• 関数基底行列における直交列構造と正の重みを持つn項四則積分規則の存在との間の双対性を確立すること。
• 非正定値重みを扱えるようにフレームワークを拡張し、積分方程式および画像処理問題のための四則積分を構築すること。
• グラミアンから導出される商行列の固有値分解を用いた、このような四則積分を生成する構成的アルゴリズムを提供すること。
• 非定値重みを伴う高次元への一般化を図り、縮退技術を用いて2変数および3変数多項式空間における四則積分規則を構築すること。

提案手法

• n次元関数空間T(n,x)に含まれるすべての内積を統合するタイプ2四則積分を提案し、直交基底行列Q(n, {x_j})のn列が直交することを保証する。
• 内積条件とモーメント条件をそれぞれB(n,n) = T(n,x) · T(x,n)およびA(n,n) = T(n,x) · μ(x)T(x,n)として定義する。
• 両方のグラミアンが正確に統合される場合、四則積分ノードx_jは行列商AB⁻¹の固有値であることを確立する。
• 双対性条件を適用：Q(n, {x_j})の列がノルム1/√w_jで直交するならば、四則積分規則{x_j, w_j}は単位グラミアンを正確に統合する。
• 既知のノードとその寄与をグラミアンAとBから反復的に除去するための縮退技術を導入し、縮退した商行列の固有値分解により追加のノードを回復可能にする。
• 複数のモーメント行列（例：A_x, A_y）を用いて次元を拡張し、縮退を適用することで2変数および3変数多項式空間におけるノードを回復する。

実験結果

リサーチクエスチョン

• RQ1 与えられた関数空間内にn個のノードを持つ四則積分規則が、n個の関数のすべての内積を正確に統合可能であるか。このような規則が存在する条件は何か。
• RQ2 古典的ガウス四則積分法の枠組みが非正定値重みのため適用不能となる場合、四則積分ノードと重みをどのように計算できるか。
• RQ3 基底関数の列の直交性と有効なn項四則積分規則の存在との間の関係は何か。
• RQ4 関数空間の次元を超える必要ノード数を要する高次元統合問題に、この手法をどのように拡張できるか。
• RQ5 縮退技術を用いて、商行列の固有値分解によって得られるノードに加えて、追加の四則積分ノードを反復的に回復できるか。

主な発見

• n項内積四則積分の存在の必要十分条件は、直交基底行列Q(n, {x_j})の列が互いに直交することである。
• 四則積分ノードx_jは、重みuおよび単純関数μの下での関数空間のグラミアンAとBの商行列AB⁻¹の固有値である。
• 多項式基底の場合、μ(x) = xのとき、この手法は古典的ガウス四則積分ノードをAB⁻¹の固有値として回復する。
• n未満の次数の2変数多項式空間では、固有値分解により最大n(n+1)/2個のノードを生成可能であり、追加のノードは縮退により回復可能である。
• 縮退により、2次元では最大n(n−1)/6個の追加ノードを回復可能であり、この手法は高次元におけるガウス・ラダウおよびガウス・ロバート型規則の構築に有効である。
• Coulomb型ポテンシャルの急速な減衰に起因する理論的収束性が、初期ノード推定値からの誤差伝搬を最小限に抑える。