<p>This paper investigates quasi-Monte Carlo (QMC) integration of Lebesgue integrable functions with respect to a density function over <InlineEquation ID="IEq4"> <EquationSource Format="TEX">\(\mathbb {R}^s\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mrow> <mi mathvariant="double-struck">R</mi> </mrow> <mi>s</mi> </msup> </math></EquationSource> </InlineEquation>. We extend the construction-free median QMC rule proposed by Goda and L’Ecuyer (SIAM J Sci Comput, 2022) to the weighted unanchored Sobolev space of functions defined over <InlineEquation ID="IEq5"> <EquationSource Format="TEX">\(\mathbb {R}^s\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mrow> <mi mathvariant="double-struck">R</mi> </mrow> <mi>s</mi> </msup> </math></EquationSource> </InlineEquation> introduced by Nichols and Kuo (J Complexity 2014). By taking the median of <InlineEquation ID="IEq6"> <EquationSource Format="TEX">\(k = \mathcal {O}(\log N)\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>k</mi> <mo>=</mo> <mi mathvariant="script">O</mi> <mo stretchy="false">(</mo> <mo>log</mo> <mi>N</mi> <mo stretchy="false">)</mo> </mrow> </math></EquationSource> </InlineEquation> independent randomized QMC estimators, we prove that for any <InlineEquation ID="IEq7"> <EquationSource Format="TEX">\(\epsilon \in (0,r-\frac{1}{2}]\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>ϵ</mi> <mo>∈</mo> <mo stretchy="false">(</mo> <mn>0</mn> <mo>,</mo> <mi>r</mi> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo stretchy="false">]</mo> </mrow> </math></EquationSource> </InlineEquation>, our method achieves a mean absolute error bound of <InlineEquation ID="IEq8"> <EquationSource Format="TEX">\(\mathcal {O}(N^{-r+\epsilon })\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi mathvariant="script">O</mi> <mo stretchy="false">(</mo> <msup> <mi>N</mi> <mrow> <mo>-</mo> <mi>r</mi> <mo>+</mo> <mi>ϵ</mi> </mrow> </msup> <mo stretchy="false">)</mo> </mrow> </math></EquationSource> </InlineEquation>, where <i>N</i> is the number of points and <InlineEquation ID="IEq9"> <EquationSource Format="TEX">\(r&gt;\frac{1}{2}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>r</mi> <mo>&gt;</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </math></EquationSource> </InlineEquation> is a parameter determined by the function space. This rate matches the rate of randomly shifted lattice rules obtained via a component-by-component (CBC) construction, while our approach requires no specific CBC constructions or prior knowledge of the space’s weight structure. Numerical experiments demonstrate that our method attains an accuracy comparable to the CBC construction based method, and outperforms the Monte Carlo method.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

A median QMC method for unbounded integrands over \(\mathbb {R}^{s}\) in weighted unanchored Sobolev spaces

  • Ziyang Ye,
  • Josef Dick,
  • Xiaoqun Wang

摘要

This paper investigates quasi-Monte Carlo (QMC) integration of Lebesgue integrable functions with respect to a density function over \(\mathbb {R}^s\) R s . We extend the construction-free median QMC rule proposed by Goda and L’Ecuyer (SIAM J Sci Comput, 2022) to the weighted unanchored Sobolev space of functions defined over \(\mathbb {R}^s\) R s introduced by Nichols and Kuo (J Complexity 2014). By taking the median of \(k = \mathcal {O}(\log N)\) k = O ( log N ) independent randomized QMC estimators, we prove that for any \(\epsilon \in (0,r-\frac{1}{2}]\) ϵ ( 0 , r - 1 2 ] , our method achieves a mean absolute error bound of \(\mathcal {O}(N^{-r+\epsilon })\) O ( N - r + ϵ ) , where N is the number of points and \(r>\frac{1}{2}\) r > 1 2 is a parameter determined by the function space. This rate matches the rate of randomly shifted lattice rules obtained via a component-by-component (CBC) construction, while our approach requires no specific CBC constructions or prior knowledge of the space’s weight structure. Numerical experiments demonstrate that our method attains an accuracy comparable to the CBC construction based method, and outperforms the Monte Carlo method.