Approximation Exercises

Ex1-3: Find the polynomial of least degree that interpolates the following data 𝒟={(xi,yi,yi',yi''),i=0,1,,n}.

Exercise 1. Find the polynomial of least degree that interpolates the data 𝒟={(xi,yi),i=0,1,,n}={(3,10),(7,146),(1,2),(2,1)}.

Exercise 2. Find the polynomial of least degree that interpolates the data 𝒟={(xi,yi,yi'),i=0,1,,n}=𝒟={(3,10,14),(1,2,-6)}.

Exercise 3. Find the polynomial of least degree that interpolates the data 𝒟={(xi,yi,yi',yi''),i=0,1,,n}=𝒟={(1,2,-6,10)}.

Exercise 4. Consider ={f|f:[a,b]}, 𝒫n={p|p.(t)=a0+a1t++antn}, 𝒙n+1, xi[a,b], ijxixj. Define the operators S:(n+1)×2, I:(n+1)×2𝒫n,

f,𝒚=f(𝒙)n+1,S(f)=[ 𝒙 𝒚 ],

I([ 𝒙 𝒚 ])=pn(t),pn(𝒙)=𝒚=f(𝒙).

For each of the mappings, S, I, P=IS, determine if the mapping is linear, and if so write its matrix representation. Also determine if the mappings are surjective, bijective, one-to-one.

Exercise 5. Recall the eigenvalue problem for finite-dimensional linear operators, 𝑨:mm, 𝑨𝑿=𝑿𝚲, 𝑿m×m, 𝚲m×m, 𝚲=diag(λ1,,λm). The generalization of the finite-dimensional problem to linear operators on vector space 𝒰, L:𝒰𝒰 is immediate, Lu=λu, u𝒰, λ. What is the appropriate generalization of C(𝑿) for the operator P:?

Exercise 6. Provide a relative error bound for the polynomial interpolation of f:[-1,1], f(t)=cosh(t) by a polynomial p22(t) of degree n=22.

Exercise 7. Prove that the above bound ε5×10-16, irrespective of the choice of distinct interpolation nodes.

Exercise 8. In the limit x1x0 the divided difference

f[x0,x1]=[y0,y1]=y1-y0x1-x0,yi=f(xi),i=0,1,

has limit f[x0,x1]f'(x0). Write and establish the validity of the finite difference form of the product rule (fg)'=f'g+fg'.

Exercise 9. Repeat the above for second order finite differences and (fg)''=f''g+2f'g'+fg''.

Exercise 10. A natural cubic spline has zero curvature at the end points. Prove that of all cubic spline interpolations of data 𝒟={(xi,yi=f(xi)),i=0,1,,n}, the natural spline S(t) curvature two-norm is bounded by the function curvature two-norm

x0xn[S''(t)]2dtx0xn[f''(t)]2dt.

Exercise 11. Recall the behavior of global polynomial interpolation for the Runge function, and interpret the significance of the above inequality for spline interpolation. In particular study the bound for the curvature of the global polynomial interpolant

x0xn[S''(t)]2dt.

Exercise 12. Prove that best inf-norm approximant of f:[-1,1], f(t)=cosh(t) by a quadratic polynomial has form p2(t)=a+bt2, with b=cosh1-1. Compute a.

Exercise 13. Compare the error bound for p2(t) from Ex 12 with that for p22(t) from Ex. 6.

Exercise 14. Prove that

-11f(t)dt190[7f(0)+32f(14)+12f(12)+32f(34)+7f(1)]

is exact for f𝒫4, (polynomials of degree up to 4).

Exercise 15. Determine the quadrature formula

01f(t)dtc0f(0)+c1f(1)

that is exact for all functions of form f(t)=aet+bcos(πt/2).

Exercise 16. Determine the quadrature formula

02πf(t)dtc0f(0)+c1f(1)

that is exact for all functions of form f(t)=a+bcos(t).

Exercise 17. Verify that the above formula is exact for

f(t)=k=0n[akcos(2k+1)t+bksinkt]

Exercise 18. Can the quadrature

01f(t)dtα[f(x0)+f(x1)]

be exact for all quadratic polynomials?

Exercise 19. Let the quadrature

-11f(t)dti=0ncif(xi)

be exact for all polynomials of degree n=2k, k. Show that if the quadrature nodes xi are symmetric around the origin the formula is also exact for polynomials of degree n+1.

Exercise 20. Determine the minimum number of subintervals needed to approximate

12(t+e-t2)dt

to absolute accuracy e<0.5×10-7 using the composite trapezoid rule.

Exercise 21. Repeat above for the composite Simpson rule.

Exercise 22. Consider the equation f(x)=0, f:, fC(2)(). Generate an approximating sequence {xn}n, xnr as n, f(r)=0, based upon constructing the linear interpolant of data 𝒟={(xn-2,fn-2),(xn-1,fn-1)} to obtain xn (fk=f(xk)).

Exercise 23. Consider the equation f(x)=0, f:, fC(2)(). Generate an approximating sequence {xn}n, xnr as n, f(r)=0, based upon constructing the linear interpolant of data 𝒟={(xn-1,fn-1,fn-1')} to obtain xn (fk=f(xk),fk'=f(xk)).

Exercise 24. Consider the equation f(x)=0, f:, fC(2)(). Generate an approximating sequence {xn}n, xnr as n, f(r)=0, based upon constructing the linear interpolant of data 𝒟={(xn-3,fn-3),(xn-2,fn-2),(xn-1,fn-1)} to obtain xn (fk=f(xk)).

Exercise 25. Consider the equation f(x)=0, f:, fC(2)(). Generate an approximating sequence {xn}n, xnr as n, f(r)=0, based upon constructing the linear interpolant of data 𝒟={(xn-1,fn-1,fn-1',fn-1'')} to obtain xn (fk=f(xk),fk'=f(xk), fk''=f''(xk)).