Difference between revisions of "Rectangle rule"
(TeX) |
m (fix TeX) |
||
Line 2: | Line 2: | ||
A formula for calculating an integral over a finite interval $[a,b]$: | A formula for calculating an integral over a finite interval $[a,b]$: | ||
− | + | \begin{equation}\label{eq:1} | |
+ | \int\limits_a^bf(x)dx\cong h\sum_{k=1}^Nf(\alpha+(k-1)h), | ||
+ | \end{equation} | ||
where $h=(b-a)/N$ and $\alpha\in[a,a+h]$. Its algebraic degree of accuracy is 1 if $\alpha=a+h/2$ and 0 otherwise. | where $h=(b-a)/N$ and $\alpha\in[a,a+h]$. Its algebraic degree of accuracy is 1 if $\alpha=a+h/2$ and 0 otherwise. | ||
− | The [[ | + | The [[quadrature formula]] \eqref{eq:1} is exact for the trigonometric functions |
$$\cos\frac{2\pi}{b-a}kx,\quad\sin\frac{2\pi}{b-a}kx,\quad k=0,\dots,N-1.$$ | $$\cos\frac{2\pi}{b-a}kx,\quad\sin\frac{2\pi}{b-a}kx,\quad k=0,\dots,N-1.$$ | ||
− | If $b-a=2\pi$, then \ | + | If $b-a=2\pi$, then \eqref{eq:1} is exact for all trigonometric polynomials of order at most $N-1$; moreover, its trigonometric degree of accuracy is $N-1$. No other quadrature formula with $N$ real nodes can have trigonometric degree of accuracy larger than $N-1$, so that the rectangle rule with $b-a=2\pi$ has the highest trigonometric degree of accuracy. |
− | Let $R(f,\alpha)$ be the error of the rectangle rule, i.e. the difference between the left- and right-hand sides of \ | + | Let $R(f,\alpha)$ be the error of the rectangle rule, i.e. the difference between the left- and right-hand sides of \eqref{eq:1}. If the integrand $f$ is twice continuously differentiable on $[a,b]$, then for $\alpha=a+h/2$ one has |
$$R\left(f,a+\frac h2\right)=\frac{b-a}{24}h^2f''(\xi),$$ | $$R\left(f,a+\frac h2\right)=\frac{b-a}{24}h^2f''(\xi),$$ | ||
Line 20: | Line 22: | ||
$$R(f,\alpha)=-(b-a)B_{2k}\frac{h^{2k}}{(2k)!}f^{(2k)}(\eta),$$ | $$R(f,\alpha)=-(b-a)B_{2k}\frac{h^{2k}}{(2k)!}f^{(2k)}(\eta),$$ | ||
− | for some $\eta\in[a,b]$, where $B_{2k}$ is the Bernoulli number (cf. [[ | + | for some $\eta\in[a,b]$, where $B_{2k}$ is the Bernoulli number (cf. [[Bernoulli numbers]]). |
Latest revision as of 17:35, 24 March 2018
A formula for calculating an integral over a finite interval $[a,b]$:
\begin{equation}\label{eq:1} \int\limits_a^bf(x)dx\cong h\sum_{k=1}^Nf(\alpha+(k-1)h), \end{equation}
where $h=(b-a)/N$ and $\alpha\in[a,a+h]$. Its algebraic degree of accuracy is 1 if $\alpha=a+h/2$ and 0 otherwise.
The quadrature formula \eqref{eq:1} is exact for the trigonometric functions
$$\cos\frac{2\pi}{b-a}kx,\quad\sin\frac{2\pi}{b-a}kx,\quad k=0,\dots,N-1.$$
If $b-a=2\pi$, then \eqref{eq:1} is exact for all trigonometric polynomials of order at most $N-1$; moreover, its trigonometric degree of accuracy is $N-1$. No other quadrature formula with $N$ real nodes can have trigonometric degree of accuracy larger than $N-1$, so that the rectangle rule with $b-a=2\pi$ has the highest trigonometric degree of accuracy.
Let $R(f,\alpha)$ be the error of the rectangle rule, i.e. the difference between the left- and right-hand sides of \eqref{eq:1}. If the integrand $f$ is twice continuously differentiable on $[a,b]$, then for $\alpha=a+h/2$ one has
$$R\left(f,a+\frac h2\right)=\frac{b-a}{24}h^2f''(\xi),$$
for some $\xi\in[a,b]$. If $f$ is a periodic function with period $b-a$ and has a continuous derivative of order $2k$ (where $k$ is a natural number) on the entire real axis, then for any $\alpha\in[a,a+h]$,
$$R(f,\alpha)=-(b-a)B_{2k}\frac{h^{2k}}{(2k)!}f^{(2k)}(\eta),$$
for some $\eta\in[a,b]$, where $B_{2k}$ is the Bernoulli number (cf. Bernoulli numbers).
Comments
References
[a1] | D.M. Young, R.T. Gregory, "A survey of numerical mathematics" , Dover, reprint (1988) pp. 362ff |
Rectangle rule. Encyclopedia of Mathematics. URL: http://encyclopediaofmath.org/index.php?title=Rectangle_rule&oldid=33490