Lab for Chebyshev Polynomial Approximation

Module for Chebyshev Polynomials

   Check out the new Numerical Analysis Projects page.

Background. Chebyshev Approximation. To construct and evaluate the Chebyshev interpolating polynomial for f(x),
of degree n over the interval [Graphics:ch.txtgr1.gif], where [Graphics:ch.txtgr2.gif] is based on the nodes
[Graphics:ch.txtgr3.gif].

In practice, all we need to know is that the Chebyshev polynomial uses the interpolation nodes
[Graphics:ch.txtgr4.gif]. It will suffice to construct these data points and
use Mathematica's built in "Fit" procedure.

Load in Mathematica's graphics package "Colors".

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr5.gif]
 
 
 
 

Report to be handed in.

Computer Exercises
 

 

Exercise 1. Investigate the Chebyshev polynomial of degree n = 4 for [Graphics:ch.txtgr7.gif].

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr8.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr9.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr10.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr11.gif]

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr12.gif]

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr13.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr14.gif]

Be assured that our method of constructing the Chebyshev polynomial is consistent with the traditional method.
Indeed, Mathematica has a built in procedure for this construction. It is invoked with the following command:

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr15.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr16.gif]

Find the error for the Chebyshev approximation.

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr17.gif]

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr18.gif]

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr19.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr20.gif]

Find the error bound for the Chebyshev approximation.
The error at the right endpoint is:

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr21.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr22.gif]

The error at the local minimum is:

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr23.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr24.gif]

The error bound for the Chebyshev polynomial is:

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr25.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr26.gif]

Construct the Taylor polynomial approximation for f[x] of degree n = 4,
and compare the error bounds.

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr27.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr28.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr29.gif]

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr30.gif]

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr31.gif]

Find the error for the Taylor approximation.

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr32.gif]

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr33.gif]

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr34.gif]

Find the error bound for the Taylor approximation.
The error at the left endpoint is:

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr35.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr36.gif]

The error at the right endpoint is:

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr37.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr38.gif]

The error bound for the Taylor polynomial is:

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr39.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr40.gif]

The error for the Chebyshev polynomial is about 6% of the error for the Taylor polynomial,
as shown by the following computation.

[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr41.gif]
[Graphics:ch.txtgr6.gif][Graphics:ch.txtgr42.gif]

What do you conclude ?

 

 

 

(c) John H. Mathews, 1998