[論文レビュー] Recursive bi-orthogonalisation approach and orthogonal projectors

Laura Rebollo‐Neira|ArXiv.org|Sep 13, 2002

Mathematical Analysis and Transform Methods参考文献 15被引用数 5

ひとこと要約

この論文では、線形独立な関数が張る部分空間への直交射影子を表現する双直交関数を構築する再帰的双直交化アルゴリズムを提案する。行列の逆行列計算を回避し、新しい関数が追加された際にも効率的な更新が可能である。この手法は正確な双直交性を保証し、逆演算や明示的な行列の逆行列計算を必要とせず、真の直交射影子を生成する。再帰的更新による収束性と安定性が証明されている。

ABSTRACT

An approach is proposed which, given a family of linearly independent functions, constructs the appropriate biorthogonal set so as to represent the orthogonal projector operator onto the corresponding subspace. The procedure evolves iteratively and it is endowed with the following properties: i) it yields the desired biorthogonal functions avoiding the need of inverse operations ii) it allows to quickly update a whole family of biorthogonal functions each time that a new member is introduced in the given set. The approach is of particular relevance to the approximation problem arising when a function is to be represented as a finite linear superposition of non orthogonal waveforms.

研究の動機と目的

• 線形独立な関数が張る部分空間への直交射影子を生成する数値的に安定した反復的手法を構築すること。
• 非直交ウェーブフォームの双直交基底を構築する際の行列逆行列計算の計算コストを回避すること。
• 部分空間に新しい関数が追加された際に、全体の双直交関数集合を効率的に更新できること。
• 反復処理中に線形従属な関数を検出し、除外するための基準を提供すること。
• 得られた双直交系が部分空間上への有効な直交射影子作用素を生成すること。

提案手法

• アルゴリズムは、新しい関数が現在の部分空間に射影されるのを基に、以前の双直交関数を修正するグラム・シュミットに類似した更新を再帰的に用いて双直交関数 $\tilde{\alpha}_n^{k+1}$ を構築する。
• 各ステップ $k+1$ において、新しい双直交関数は $|\tilde{\alpha}_n^{k+1}\rangle = |\tilde{\alpha}_n^k\rangle - \frac{\langle\alpha_m|\psi_{k+1}\rangle}{|||\psi_{k+1}\rangle||^2}\langle\alpha_{k+1}|\tilde{\alpha}_n^k\rangle$ により更新される。ここで $|\psi_{k+1}\rangle$ は直交化後の残差である。
• 残差 $|\psi_{k+1}\rangle$ は $|\psi_{k+1}\rangle = |\alpha_{k+1}\rangle - \sum_{m=1}^k \frac{\langle\alpha_m|\alpha_{k+1}\rangle}{|||\alpha_m^k\rangle||^2}|\alpha_m^k\rangle$ として定義され、これにより以前の部分空間に対して直交することが保証される。
• 双直交関数は再帰的射影と正規化ステップを直接計算することで、行列の逆行列計算を回避する。
• 射影子は $\hat{P}_{k+1} = \sum_{n=1}^{k+1} |\tilde{\alpha}_n^{k+1}\rangle\langle\alpha_n|$ として表現され、双直交性と射影不変性により、これが直交射影子であることが証明されている。
• ケース解析と数学的帰納法により、すべての $m,n \leq k+1$ に対して双直交性 $\langle\alpha_m|\tilde{\alpha}_n^{k+1}\rangle = \delta_{m,n}$ が保存されることを証明した。

実験結果

リサーチクエスチョン

• RQ1 非直交的で線形独立な関数が張る部分空間への直交射影子を生成する双直交基底を、行列の逆行列計算を伴わずどのように構築できるか？
• RQ2 部分空間に新しい関数が追加された際に、双直交関数を効率的に更新できるか？
• RQ3 どの再帰的メカニズムが、アルゴリズムの各ステップで双直交性と射影子の有効性を保証するか？
• RQ4 構築過程で線形従属な関数をどのように検出し、除外するか？
• RQ5 どのような条件下で、得られた作用素が一般の双直交表現ではなく、真の直交射影子であると保証できるか？

主な発見

• 提案された再帰的アルゴリズムは、与えられた関数が張る部分空間への有効な直交射影子を生成する双直交関数を構築するが、行列の逆行列計算を必要としない。
• 再帰的射影と正規化により、すべての $m,n \leq k+1$ に対して正確な双直交性 $\langle\alpha_m|\tilde{\alpha}_n^{k+1}\rangle = \delta_{m,n}$ が保証される。
• 射影子作用素 $\hat{P}_{k+1} = \sum_{n=1}^{k+1} |\tilde{\alpha}_n^{k+1}\rangle\langle\alpha_n|$ は自己随伴かつ冪等的であることが証明されており、これが直交射影子であることが確認された。
• 新しい関数が追加された際に、全体の双直交関数集合を効率的に更新でき、再計算から完全に回避することができる。
• 線形従属な関数は自動的に検出され、除外される。これは、新しい関数が以前の関数の線形結合に含まれる場合、残差 $|\psi_{k+1}\rangle$ がゼロになるためである。
• この手法は数値的に安定しており、行列の逆行列計算に伴う悪条件化を回避するため、大規模または適応的近似問題に適している。