Paper 3, Section II, C

Numerical Analysis | Part IB, 2019

(a) Let w(x)w(x) be a positive weight function on the interval [a,b][a, b]. Show that

⟨f,g⟩=∫abf(x)g(x)w(x)dx\langle f, g\rangle=\int_{a}^{b} f(x) g(x) w(x) d x

defines an inner product on C[a,b]C[a, b].

(b) Consider the sequence of polynomials pn(x)p_{n}(x) defined by the three-term recurrence relation

pn+1(x)=(x−αn)pn(x)−βnpn−1(x),n=1,2,…p_{n+1}(x)=\left(x-\alpha_{n}\right) p_{n}(x)-\beta_{n} p_{n-1}(x), \quad n=1,2, \ldots

where

p0(x)=1,p1(x)=x−α0,p_{0}(x)=1, \quad p_{1}(x)=x-\alpha_{0},

and the coefficients αn\alpha_{n} (for n⩾0)\left.n \geqslant 0\right) and βn\beta_{n} (for n⩾1)\left.n \geqslant 1\right) are given by

αn=⟨pn,xpn⟩⟨pn,pn⟩,βn=⟨pn,pn⟩⟨pn−1,pn−1⟩\alpha_{n}=\frac{\left\langle p_{n}, x p_{n}\right\rangle}{\left\langle p_{n}, p_{n}\right\rangle}, \quad \beta_{n}=\frac{\left\langle p_{n}, p_{n}\right\rangle}{\left\langle p_{n-1}, p_{n-1}\right\rangle}

Prove that this defines a sequence of monic orthogonal polynomials on [a,b][a, b].

(c) The Hermite polynomials Hen(x)H e_{n}(x) are orthogonal on the interval (−∞,∞)(-\infty, \infty) with weight function e−x2/2e^{-x^{2} / 2}. Given that

Hen(x)=(−1)nex2/2dndxn(e−x2/2)H e_{n}(x)=(-1)^{n} e^{x^{2} / 2} \frac{d^{n}}{d x^{n}}\left(e^{-x^{2} / 2}\right)

deduce that the Hermite polynomials satisfy a relation of the form (∗)(*) with αn=0\alpha_{n}=0 and βn=n\beta_{n}=n. Show that ⟨Hen,Hen⟩=n!2π\left\langle H e_{n}, H e_{n}\right\rangle=n ! \sqrt{2 \pi}.

(d) State, without proof, how the properties of the Hermite polynomial He⁡N(x)\operatorname{He}_{N}(x), for some positive integer NN, can be used to estimate the integral

∫−∞∞f(x)e−x2/2dx\int_{-\infty}^{\infty} f(x) e^{-x^{2} / 2} d x

where f(x)f(x) is a given function, by the method of Gaussian quadrature. For which polynomials is the quadrature formula exact?

Typos? Please submit corrections to this page on GitHub.