Koksma-Hlawka型不等式を礎とする準モンテカルロ法の研究
Project/Area Number |
20K14326
|
Research Category |
Grant-in-Aid for Early-Career Scientists
|
Allocation Type | Multi-year Fund |
Review Section |
Basic Section 12010:Basic analysis-related
|
Research Institution | Hiroshima University |
Principal Investigator |
鈴木 航介 広島大学, 先進理工系科学研究科(理), 助教 (20868674)
|
Project Period (FY) |
2020-04-01 – 2023-03-31
|
Project Status |
Discontinued (Fiscal Year 2022)
|
Budget Amount *help |
¥4,160,000 (Direct Cost: ¥3,200,000、Indirect Cost: ¥960,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,170,000 (Direct Cost: ¥900,000、Indirect Cost: ¥270,000)
|
Keywords | 準モンテカルロ法 / 乱択化準モンテカルロ法 / ソボル列 / スクランブル / デジタルネット / 優良格子点法 / 数値積分 / Koksma-Hlawka型不等式 |
Outline of Research at the Start |
高次元数値積分のアルゴリズムとして、サンプル点を一様ランダムに選ぶモンテカルロ法がよく知られている。しかしモンテカルロ法の収束オーダーは、サンプル点の個数をNとしたときO(1/sqrt(N))と決して速くない。ここで超一様点集合という均等に散らばる点集合をサンプルに用いる準モンテカルロ法を用いると、性質のよい関数の積分誤差が1/N、もしくはそれ以上に高速なオーダーで収束する。本研究では、準モンテカルロ法の理論的根拠となるKoksma-Hlawka型不等式に基づき、不等式の離散化、積分領域の一般化、幅広い関数に対応できる積分則の構築などを研究する。
|
Outline of Annual Research Achievements |
高次元数値積分のアルゴリズムとして、サンプル点を一様ランダムに選ぶモンテカルロ法がよく知られているが、収束オーダーはサンプル点の個数の平方根の逆数となり決して速くない。準モンテカルロ法では、超一様点集合という注意深く設計された点集合をサンプルに用いることでモンテカルロ法よりも高速な積分誤差の収束を目指す。応用上では、両アルゴリズムの利点を得るために点集合の一様性を壊さないようにランダマイズする乱択化準モンテカルロ法(RQMC)という手法がよく用いられる。ここ2年ほど、RQMCにおいて平均値ではなく中央値を推定値とする手法(ここではmedian-RQMCと呼ぶ)が提案されてきた。この手法では、関数の性質を知ることなく、滑らかさなどの関数の良い性質を自動的に反映した高速な誤差収束が高い確率で得られる。本研究の最終年度では、scrambled sobol' (スクランブルソボル)列による median-RQMCについて研究し、以下の結果を得た。 一つ目の研究では、スクランブルソボル列 の gain coefficient と呼ばれる量を調べ、既存のものよりも良い評価式を得た。本研究では、ソボル列のデジタルネットとしての構造に注目して双対性を使い、議論を双対デジタルネット上での数え上げに帰着した。 二つ目の研究では、スクランブルソボル列による median-RQMC に対して、滑らかさなどの関数の良い性質を自動的に反映した高速な誤差収束を証明した。既存の結果を超え、本研究では関数の周期性を要求せず、非常に広いクラスの関数に対してほぼ最適な収束を示した。そのため多くの例に対し応用可能であり、意義のある結果だと考えている。
|
Report
(3 results)
Research Products
(5 results)