i x i y i

Transcription

1 Department of Mathematics MTL107: Numerical Methods and Computations Exercise Set 8: Approximation-Linear Least Squares Polynomial approximation, Chebyshev Polynomial approximation. 1. Compute the linear least square polynomial for the data: i x i y i Find the least square polynomials of degrees 1,2 and 3 for the data in the following talbe. Compute the error E in each case. Graph the data and the polynomials. : x i y i Given the data: x i y i a. Construct the least squared polynomial of degree 1, and compute the error. b. Construct the least squared polynomial of degree 2, and compute the error. c. Construct the least squared polynomial of degree 3, and compute the error. d. Construct the least squares approximation of the form be ax, and compute the error. e. Construct the least squares approximation of the form bx a, and compute the error. 4. The following table lists the college grade-point averages of 20 mathematics and computer science majors, together with the scores that these students received on the mathematics portion of the ACT (Americal College Testing Program) test while in high school. Plot these data, and find the equation of the least squares line for this data:

2 : ACT Grade-point ACT Grade-point score average score average Find the linear least squares polynomial approximation to f(x) on the indicated interval if a. f(x) = x 2 + 3x + 2, [0, 1]; b. f(x) = x 3, [0, 2]; c. f(x) = 1 x, [1, 3]; d. f(x) = ex, [0, 2]; e. f(x) = 1 cos x + 1 sin 2x, [0, 1]; f. f(x) = x ln x, [1, 3]; Find the least square polynomial approximation of degrees 2 to the functions and intervals in Exercise Compute the error E for the approximations in Exercise Use the Gram-Schmidt process to construct φ 0 (x), φ 1 (x), φ 2 (x) and φ 3 (x) for the following intervals. a. [0,1] b. [0,2] c. [1,3] 9. Obtain the least square approximation polynomial of degree 3 for the functions in Exercise 5 using the results of Exercise Use the Gram-Schmidt procedure to calculate L 1, L 2, L 3 where {L 0 (x), L 1 (x), L 2 (x), L 3 (x)} is an orthogonal set of polynomials on (0, ) with respect to the weight functions w(x) = e x and L 0 (x) = 1. The polynomials obtained from this procedure are called the Laguerre polynomials. 11. Use the zeros of T 3, to construct an interpolating polynomial of degree 2 for the following functions on the interval [-1,1]: a. f(x) = e x, b. f(x) = sin x, c. f(x) = ln(x + 2), d. f(x) = x 4.

3 12. Find a bound for the maximum error of the approximation in Exercise 1 on the interval [-1,1]. 13. Use the zeros of T 3 and transformations of the given interval to construct an interpolating polynomial of degree 2 for the following functions on the indicated intervals: a. f(x) = 1 x, [1, 3]; b. f(x) = e x, [0, 2]; c. f(x) = 1 cos x + 1 sin 2x, [0, 1]; d. f(x) = x ln x, [1, 3]; Find the sixth Maclaurin polynomial for sin x, and use Chebyshev economization to obtain a lesser degree polynomial approximation while keeping the error less than 0.01 on [-1,1]. 15. Show that for each Chebyshev polynomial T n (x), we have 1 [T n (x)] 2 1 x 2 dx = π Show that for each n, the derivative of the Chebyshev polynomial T n (x), has n distinct zeros in (-1,1). ANSWERS 1. The linear least-squares polynomial is x The least-squares polynomials with their errors are, respectively, x, with E = ; x x 2, with E = ; and x x x 3, with E = a. The linear least-squares polynomial is x , with error 329. b. The least-squares polynomial of degree two is x x , with error c. The least-squares polynomial of degree three is x x x , with error d. The least-squares polynomial of the form be ax is e x, with error 418. e. The least-squares polynomial of the form bx a is x , with error The least squares line for the point average is (ACT score)

4 5. The linear least-sqaures approximations are: a. P 1 (x) = x b. P 1 (x) = x c. P 1 (x) = x d. P 1 (x) = x e. P 1 (x) = x f. P 1 (x) = x. 6. The linear least-sqaures approximations of degree two are: a. P 2 (x) = x x 2 b. P 2 (x) = x x 2 c. P 2 (x) = x x 2 d. P 2 (x) = x x 2 e. P 2 (x) = x x 2 f. P 2 (x) = x x a b c d e f The Gram-Schmidt process produces the following collections of polynomials: a. φ 0 (x) = 1, φ 1 (x) = x 0.5, φ 2 (x) = x 2 x + 1 6, and φ 3(x) = x 3 1.5x x 0.05 b. φ 0 (x) = 1, φ 1 (x) = x 1, φ 2 (x) = x 2 2x + 2 3, and φ 3(x) = x 3 3x x 2 5 c. φ 0 (x) = 1, φ 1 (x) = x 2, φ 2 (x) = x 2 4x , and φ 3(x) = x 3 6x x The least-squares polynomial of degree two are: a. P 2 (x) = φ 0 (x) + 4φ 1 (x) φ 2 (x), b P 2 (x) = 2φ 0 (x) + 3.6φ 1 (x) + 3φ 2 (x), c. P 2 (x) = φ 0 (x) φ 1 (x) φ 2 (x), d. P 2 (x) = φ 0 (x) + 3φ 1 (x) φ 2 (x), e. P 2 (x) = φ 0 (x) φ 1 (x) φ 2 (x), f. P 2 (x) = φ 0 (x) φ 1 (x) φ 2 (x). 10. The Laguerre polynomials are L 1 (x) = x 1, L 1 (x) = x 2 4x + 2, and L 3 (x) = x 3 9x x The interpolating polynomials of degree two are: a. P 2 (x) = (x ) x(x ) b. P 2 (x) = (x ) c. P 2 (x) = (x ) x(x ) d. P 2 (x) = (x ) x(x ) 12. Bounds for the maximum errors of polynomials in Exercise 11 are: a b c d The zeros of T 3 produce the following interpolating polynomials of degree two: a. P 2 (x) = (x ) (x )(x 2)

5 b. P 2 (x) = (x ) (x )(x 1) c. P 2 (x) = (x ) (x )(x 0.5) d. P 2 (x) = (x ) (x )(x 2) 14. The cubic polynomial x 5 32 x3 approximates sin x with error at most The change of variable x = cos θ produces 1 T n 2(x) 1 x 2 dx = 1 [cos(n arccos x)] 2 1 x 2 dx = π 0 (cos(nθ))2 dθ = π It was shown in text that the zeros of T n(x) occur at x k = cos(kπ/n) for k = 1,.., n 1. Because x 0 = cos(0), x n = cos(π/) =, and all values of the cosine lie in the interval [-1,1] it reamins only to show that the zeros are distinct. This follows from the fact that for each k = 1,.., n 1, we have x k in the interval (0, π) and on this interval D x cos(x) = sin x < 0. As a consequence, T n(x) is one-to-one on (0, π), and these n 1 zeros of T n(x) are distinct.