Section 13.8

The Least Squares Method and Interpolation of Functions

[13/54%]
Show
LaTeX
Problem 1395

Consider some vector xx and a subspace LL in the Euclidean space Rn\mathbb {R}^{n}. Prove that among the vectors of the subspace LL, the orthogonal projection x∥x^{\| } of the vector xx onto LL ``deviates least'' from the given vector xx, i.e.,

∣x−x∥∣⩽∣x−l∣∀l∈L \left|x-x^{\| }\right| \leqslant \left|x-l\right| \quad \forall l \in L
?
Problem 1396

Prove that a pseudosolution can be found as the solution of a new inhomogeneous system of equations formed using the Gram matrix G(a1,a2,…,an)G\left(a_{1}, a_{2}, \ldots , a_{n}\right):

((a1,a1)(a1,a2)…(a1,an)⋮⋮⋮(an,a1)(an,a2)…(an,an))(x1⋮xn)=((b,a1)⋮(b,an)) \left(\begin{smallmatrix} \left(a_{1}, a_{1}\right) & \left(a_{1}, a_{2}\right) & \ldots & \left(a_{1}, a_{n}\right) \\ \vdots & \vdots & & \vdots \\ \left(a_{n}, a_{1}\right) & \left(a_{n}, a_{2}\right) & \ldots & \left(a_{n}, a_{n}\right) \end{smallmatrix}\right)\left(\begin{smallmatrix} x_{1} \\ \vdots \\ x_{n} \end{smallmatrix}\right) =\left(\begin{smallmatrix} \left(b, a_{1}\right) \\ \vdots \\ \left(b, a_{n}\right) \end{smallmatrix}\right)

What is the geometric meaning of the constants x1,x2,…,xnx_{1}, x_{2}, \ldots , x_{n} found in this way?

?
Problem 1397

Using the method of least squares, find the pseudosolution of the following inconsistent systems of linear equations:

?
(1)

{x+2y=12x+4y=1x+y=0\begin{cases} x+2 y = 1 \\ 2 x+4 y = 1 \\ x+y = 0 \end{cases}

(2)

{x=22x−y=1,x−2y=−1,x+y=2\begin{cases} x = 2 \\ 2 x-y = 1, \\ x-2 y = -1, \\ x+y = 2 \end{cases}

(3)

{2x1−x3=1,x2+x3=−1,x1−x2+x3=0,x1−x3=−1;\begin{cases} 2 x^{1}-x^{3} = 1, \\ x^{2}+x^{3} = -1, \\ x^{1}-x^{2}+x^{3} = 0, \\ x^{1}-x^{3} = -1 ; \end{cases}

(4)

{x1−x3=1,x2+x3−x4=−1,x1−x2=0,−x1+x3+2x4=1,x2+x3−x4=1.\begin{cases} x^{1}-x^{3} = 1, \\ x^{2}+x^{3}-x^{4} = -1, \\ x^{1}-x^{2} = 0, \\ -x^{1}+x^{3}+2 x^{4} = 1, \\ x^{2}+x^{3}-x^{4} = 1. \end{cases}

Problem 1398

Prove that the system of polynomials φk(x),k=0,…,n\varphi_{k}(x), k = 0, \ldots , n, from (13.1) forms an orthonormal basis of the space of polynomials of degree at most nn with respect to the inner product

(f,g)=f(x0)g(x0)+f(x1)g(x1)+…+f(xn)g(xn)−∞<x0<x1<…<xn<+∞ \begin{aligned} (f, g) = f\left(x_{0}\right) g\left(x_{0}\right)+f\left(x_{1}\right) g\left(x_{1}\right) & +\ldots +f\left(x_{n}\right) g\left(x_{n}\right) \\ & -\infty < x_{0} < x_{1} < \ldots < x_{n} < +\infty \end{aligned}
?
Problem 1399

Given 2n+12 n+1 pairwise distinct points xi,i=0,…,2nx_{i}, i = 0, \ldots , 2 n, of the interval [0,2π][0,2 \pi ]. In the space V2n+1V_{2 n+1} of trigonometric polynomials

V2n+1={c0+c1cos⁡x+c2sin⁡x+…+c2n−1cos⁡nx+c2nsin⁡nx} V_{2 n+1} = \left\{ c_{0}+c_{1} \cos x+c_{2} \sin x+\ldots +c_{2 n-1} \cos n x+c_{2 n} \sin n x\right\}

find trigonometric polynomials tkt_{k} such that

tk(xj)={1,j=k0,j≠k t_{k}\left(x_{j}\right) = \begin{cases} 1, & j = k \\ 0, & j \neq k \end{cases}

Show that tkt_{k} form an orthonormal basis of the space V2n+1V_{2 n+1} with respect to the inner product

(f,g)=f(x0)g(x0)+f(x1)g(x1)+…+f(x2n)g(x2n) (f, g) = f\left(x_{0}\right) g\left(x_{0}\right)+f\left(x_{1}\right) g\left(x_{1}\right)+\ldots +f\left(x_{2 n}\right) g\left(x_{2 n}\right)
?
Problem 1400

Prove that in the space V2n+1V_{2 n+1} of trigonometric polynomials there exists a unique solution T2n+1\mathscr {T}_{2 n+1} of the system of equations

T2n+1(xi)=f(xi),i=0,…,2n \mathscr {T}_{2 n+1}\left(x_{i}\right) = f\left(x_{i}\right), \quad i = 0, \ldots , 2 n

where x0,…,x2n−2n+1x_{0}, \ldots , x_{2 n}-2 n+1 pairwise distinct points of the interval [0,2π][0,2 \pi ], and f:[a,b]→Rf:[a, b] \rightarrow \mathbb {R} is some function.

?
Problem 1401

Find a polynomial hh of minimal degree such that, at the given pairwise distinct points x0,x1,…,xnx_{0}, x_{1}, \ldots , x_{n}, its values and derivatives coincide with the values and derivatives of a given function:

h(xi)=f(xi),h′(xi)=f′(xi),0⩽i⩽n h\left(x_{i}\right) = f\left(x_{i}\right), \quad h^{\prime }\left(x_{i}\right) = f^{\prime }\left(x_{i}\right), \quad 0 \leqslant i \leqslant n

Using the Lagrange interpolation polynomial LnfL_{n}^{f}, one can interpolate an arbitrary function ff on the real line from the table of its values f(x0),f(x1),…,f(xn)f\left(x_{0}\right), f\left(x_{1}\right), \ldots , f\left(x_{n}\right), where −∞<x0<……<xn<+∞-\infty < x_{0} < \ldots \ldots < x_{n} < +\infty.

?
Problem 1402

Prove that if the set {xi,i=0,…,n}\left\{ x_{i}, i = 0, \ldots , n\right\} of interpolation nodes is symmetric with respect to zero, then from the evenness (oddness) of the function f(x)f(x) it follows that the function Lnf(x)L_{n}^{f}(x) is even (odd).

?
Problem 1403

Introduce the space Vn+1V_{n+1} of discrete functions f:Mn+1→Rf: M_{n+1} \rightarrow \mathbb {R} on the finite set Mn+1={x0,…,xn}M_{n+1} = \left\{ x_{0}, \ldots , x_{n}\right\}. Show that Vn+1V_{n+1} is a Euclidean space with respect to the inner product (f,g)=f(x0)g(x0)+f(x1)g(x1)+…+f(xn)g(xn)(f, g) = f\left(x_{0}\right) g\left(x_{0}\right)+f\left(x_{1}\right) g\left(x_{1}\right)+\ldots +f\left(x_{n}\right) g\left(x_{n}\right). Exhibit some orthonormal basis of the space Vn+1V_{n+1}.

?
Problem 1404

Prove that the solution Lmf(x)L_{m}^{f}(x) of the least-squares interpolation problem in the space Vn+1V_{n+1} exists and is unique, and that the coefficients ai,i=0,…,ma_{i}, i = 0, \ldots , m, are found from the system of inhomogeneous linear equations

G(φ0,φ1,…,φm)(a0⋮am)=((f,φ0)⋮(f,φm)) G\left(\varphi _{0}, \varphi _{1}, \ldots , \varphi _{m}\right)\left(\begin{smallmatrix} a_{0} \\ \vdots \\ a_{m} \end{smallmatrix}\right) =\left(\begin{smallmatrix} \left(f, \varphi _{0}\right) \\ \vdots \\ \left(f, \varphi _{m}\right) \end{smallmatrix}\right)

where G(φ0,φ1,…,φm)G\left(\varphi_{0}, \varphi_{1}, \ldots , \varphi_{m}\right) is the Gram matrix of the system of functions φ0(x),……,φm(x)\varphi_{0}(x), \ldots \ldots , \varphi_{m}(x) in the Euclidean space Vn+1V_{n+1}.

?
Problem 1405

Using the method of least squares, find the interpolation (best mean-square approximation) of the function ff given by the values f(0)=1,f(1)=2,f(2)=3,f(3)=5f(0) = 1, f(1) = 2, f(2) = 3, f(3) = 5:

?
(1)

by a linear polynomial L1f(x)=b1x+b0L_{1}^{f}(x) = b_{1} x+b_{0} (i.e., interpolate using the system of functions φ0(x)=1,φ1(x)=x\varphi_{0}(x) = 1, \varphi_{1}(x) = x);

(2)

by a quadratic polynomial L2f(x)=b2x2+b1x+b0L_{2}^{f}(x) = b_{2} x^{2}+b_{1} x+b_{0};

(3)

by a polynomial of degree 3, L3f(x)=b3x3+b2x2+b1x+b0L_{3}^{f}(x) = b_{3} x^{3}+b_{2} x^{2}+b_{1} x+b_{0}.

Find the quadratic deviation ∥Lif−f∥4\left\| L_{i}^{f}-f\right\|_{4} in each of these cases.

Problem 1406

Prove that in the case of an orthonormal system of functions φ0(x),…,φn(x)\varphi_{0}(x), \ldots , \varphi_{n}(x), the coefficients ai,i=0,…,na_{i}, i = 0, \ldots , n, equal the Fourier coefficients of the function ff, i.e., the linear combination with the Fourier coefficients of the function ff is the best interpolation of the function ff.

?
Problem 1407

Using the method of least squares, find, in the Euclidean space of continuously differentiable functions with the inner product

(f,g)=∫02πf(x)g(x)dx (f, g) = \int _{0}^{2 \pi } f(x) g(x) d x

the interpolation (best mean-square approximation) of a continuously differentiable function ff using the system of functions φ0(x)=1,φ1(x)=cos⁡x,φ2(x)=sin⁡x\varphi_{0}(x) = 1, \varphi_{1}(x) = \cos x, \varphi_{2}(x) = \sin x.

?