English

Optimal Quadrature Formulas for the Sobolev Space $H^1$

Numerical Analysis 2018-07-17 v1

Abstract

We study optimal quadrature formulas for arbitrary weighted integrals and integrands from the Sobolev space H1([0,1])H^1([0,1]). We obtain general formulas for the worst case error depending on the nodes xjx_j. A particular case is the computation of Fourier coefficients, where the oscillatory weight is given by ρk(x)=exp(2πikx)\rho_k(x) = \exp(- 2 \pi i k x). Here we study the question whether equidistant nodes are optimal or not. We prove that this depends on nn and kk: equidistant nodes are optimal if n2.7k+1n \ge 2.7 |k| +1 but might be suboptimal for small nn. In particular, the equidistant nodes xj=j/kx_j = j/ |k| for j=0,1,,k=n+1j=0, 1, \dots , |k| = n+1 are the worst possible nodes and do not give any useful information. To characterize the worst case function we use certain results from the theory of weak solutions of boundary value problems and related quadratic extremal problems.

Keywords

Cite

@article{arxiv.1609.01146,
  title  = {Optimal Quadrature Formulas for the Sobolev Space $H^1$},
  author = {Erich Novak and Shun Zhang},
  journal= {arXiv preprint arXiv:1609.01146},
  year   = {2018}
}

Comments

24 pages, 1 figure, 1 table