Robust and Efficient Numerical Methods for Matrix Problems with Singularity
Project/Area Number |
20K14356
|
Research Category |
Grant-in-Aid for Early-Career Scientists
|
Allocation Type | Multi-year Fund |
Review Section |
Basic Section 12040:Applied mathematics and statistics-related
|
Research Institution | University of Tsukuba |
Principal Investigator |
保國 惠一 筑波大学, システム情報系, 助教 (90765934)
|
Project Period (FY) |
2020-04-01 – 2024-03-31
|
Project Status |
Granted (Fiscal Year 2022)
|
Budget Amount *help |
¥4,290,000 (Direct Cost: ¥3,300,000、Indirect Cost: ¥990,000)
Fiscal Year 2023: ¥780,000 (Direct Cost: ¥600,000、Indirect Cost: ¥180,000)
Fiscal Year 2022: ¥1,040,000 (Direct Cost: ¥800,000、Indirect Cost: ¥240,000)
Fiscal Year 2021: ¥1,170,000 (Direct Cost: ¥900,000、Indirect Cost: ¥270,000)
Fiscal Year 2020: ¥1,300,000 (Direct Cost: ¥1,000,000、Indirect Cost: ¥300,000)
|
Keywords | 数値線形代数 / 数値計算 / 非適切問題 / 固有値問題 / 特異性 / 射影法 / 最小二乗問題 / 摂動解析 / 丸め誤差 / 行列解析 / 機械学習 / 最適化 |
Outline of Research at the Start |
本課題の目的は、数値的に解くことが難しい特異や特異に近い行列問題を頑健に求解できるようにすることである。具体的には、(I)既存の特異なモデルや行列問題及び正則なモデルや行列問題を特異に拡張した場合の機序を数値解析的な視点から明らかにすること、(II)従来より広いクラスの特異性をもった問題に対して適用できるように拡張した頑健な手法を設計すること、(III)設計した手法が破綻することなく頑健に求解できるための条件を与えること、(IV)実問題・実データに対して従来法に比べて設計した手法が有効であることを計算機実験で示すこと等である。
|
Outline of Annual Research Achievements |
本課題の目的は、数値的な求解が困難な特異及び特異に近い性質をもつ行列問題に対する頑健な数理アルゴリズムを設計し、それが従来法に比べて優れることを示すことである。 非適切問題に対する非一致な投影・逆投影を用いた一般化最小残差法(GMRES法)を考案した。これまで研究してきた、係数行列の転置とは異なる行列を用いることで行列の対称性を敢えて崩して前処理するGMRES法を用いる。方程式を離散化して得られる行列およびその転置は、必ずしも互いに転置の関係となるとは限らないが、GPUを用いるような特定の計算機環境において行列ベクトル積の効率がよい実装が可能である。一次の摂動解析により、投影・逆投影および観測誤差が解に与える影響を定量的に評価した。特異空間の解析により、解が求めるべき特異空間を復元していることを示した。数値実験により、提案法が従来法と同等の精度で求解できることを示した。 また、一般エルミート行列固有値問題に対する、複素モーメント行列を使ったレイリーリッツの技法に基づく固有値・固有ベクトルの新しい精度保証付き数値計算を考案した。数値積分や線形方程式の求解等に生じるすべての誤差を評価するための区間演算に適した式を得た。近接固有値や重複固有値をもつような悪条件の問題に対しても頑健に精度保証できることを数値実験により示した。 さらに、複数右辺ベクトルをもつ線形方程式に対する非対称行列向けLanczos法が数値的に到達可能な限界精度を改善するため、残差ノルムの平滑化を行う新しい更新式を考案した。従来の平滑化手法は単調な収束を保証できる一方、到達可能精度の改善は不可能であることが理論で示されていた。本研究は、主反復と双対反復を相互に伝搬させることでこの問題を克服できることを示した。 他に、無限次元における微分作用素の新しい固有値計算手法も考案した。
|
Current Status of Research Progress |
Current Status of Research Progress
2: Research has progressed on the whole more than it was originally planned.
Reason
当初の計画したテーマによる成果は、学術論文誌に掲載された。当初の計画から派生した成果もあり、その発表準備にもあたっているが、当初の計画に近い状況にある。
|
Strategy for Future Research Activity |
当初計画していた実問題向けの実装開発を実施する。必要があれば理論的な評価を改善し、実用性の評価を示すため従来法との比較評価を進める。
|
Report
(3 results)
Research Products
(44 results)