Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Exercicios Pol Interp2

Download as pdf or txt
Download as pdf or txt
You are on page 1of 4

Exercises

273

M. Abramowitz and I.A. Stegun (Editors), Orthogonal polynomials, Ch. 22 in Handbook of Mathematical Functions with Formulas,
Graphs, and Mathematical Tables, ninth printing, Dover, New York,
pp. 771802, 1972.
Computational aspects of the theory of orthogonal polynomials are discussed in the edited volume
W. Gautschi, G.H. Golub, and G. Opfer (Editors), Applications and Computation of Orthogonal Polynomials, Conference at the
Mathematical Research Institute, Oberwolfach, Germany, March 22
28, 1998, Birkh
auser, Basel, 1999.
A recent survey of the theory and application of orthogonal polynomials
in numerical computations is contained in
W. Gautschi, Orthogonal polynomials: applications and computation, Acta Numerica 5 (A. Iserles, ed.), Cambridge University Press,
Cambridge, pp. 45119, 1996.
Finally, we refer to the books of Powell and Cheney, cited in the
Notes at the end of the previous chapter, concerning the application of
orthogonal polynomials in the eld of best least squares approximation.
Exercises
9.1
9.2

9.3

Construct orthogonal polynomials of degrees 0, 1 and 2 on the


interval (0, 1) with the weight function w(x) = ln x.
Let the polynomials j , j = 0, 1, . . . , form an orthogonal system on the interval (1, 1) with respect to the weight function
w(x) 1. Show that the polynomials j ((2x a b)/(b a)),
j = 0, 1, . . ., represent an orthogonal system for the interval
(a, b) and the same weight function. Hence obtain the polynomials in Example 9.5 from the Legendre polynomials in Example
9.6.
Suppose that the polynomials j , j = 0, 1, . . . , form an orthogonal system on the interval (0, 1) with respect to the weight
function w(x) = x , > 0. Find, in terms of j , a system
of orthogonal polynomials for the interval (0, b) and the same
weight function.

274
9.4

9 Approximation in the 2-norm


Show, by induction or otherwise, that, for 0 k n,
d
dx

(1 x2 )n = (1 x2 )nk qk (x) ,

where qk is a polynomial of degree k. Deduce that all the derivatives of the function (1 x2 )n of order less than n vanish at
x = 1.
Dene j (x) = (d/dx)j (1 x2 )j , and show by repeated integration by parts that
1

k (x)j (x)dx = 0 ,

9.5

0k <j.

Hence verify the expressions in Example 9.6 for the Legendre


polynomials of degrees 0, 1, 2 and 3.
Show, by induction or otherwise, that, for 0 k j,
d
dx

xj ex = xjk qk (x)ex ,

where qk (x) is a polynomial of degree k.


The function j is dened for j 0 by

dj j x
(x e ) .
dxj
Show that, for each j 0, j is a polynomial of degree j,
and that these polynomials form an orthogonal system on the
interval (0, ) with respect to the weight function w(x) = ex .
Write down the polynomials with j = 0, 1, 2 and 3.
Suppose that j , j = 0, 1, . . ., form a system of orthogonal polynomials with weight function w(x) on the interval (a, b). Show
that, for some value of the constant Cj , j+1 (x) Cj xj (x) is
a polynomial of degree j, and hence that
j (x) = ex

9.6

j+1 (x) Cj xj (x) =

jk k (x) ,
k=0

jk R .

Use the orthogonality properties to show that jk = 0 for


k < j 1, and deduce that the polynomials satisfy a recurrence
relation of the form
j+1 (x) (Cj x + Dj )j (x) + Ej j1 (x) = 0 ,

j 1.

Exercises
9.7

275

In the notation of Exercise 6 suppose that the normalisation of


the polynomials is so chosen that for each j the coecient of xj
in j (x) is positive. Show that Cj > 0 for all j. By considering
b
a

w(x)j (x)[j (x) Cj1 xj1 (x)]dx

show that
b

w(x)xj1 (x)j (x)dx > 0 ,

9.8

and deduce that Ej > 0 for all j. Hence show that for all
positive values of j the zeros of j and j1 interlace. (See the
proof of Theorem 5.8.)
Using the weight function w on the interval (a, b) apply a similar
argument to that for Theorem 8.6 to nd the best polynomial
approximation pn of degree n in the 2-norm to the function
xn+1 . Show that
xn+1 pn

9.9

2
2

=
a

2
w(x)2n+1 dx/[cn+1
n+1 ] ,

n+1
in n+1 (x).
where cn+1
n+1 is the coecient of x
Write down the best polynomial approximation of degree 2
to the function x3 in the 2-norm with w(x) 1 on the interval
(1, 1), and evaluate the 2-norm of the error.
Suppose that the weight w is an even function on the interval
(a, a), and that a system of orthogonal polynomials j , j =
0, . . . , n, on the interval (a, a) is constructed by the Gram
Schmidt process. Show that, if j is even, then j is an even
function, and that, if j is odd, then j is an odd function.
Now suppose that the best polynomial approximation of degree n in the 2-norm to the function f on the interval (a, a) is
expressed in the form

pn (x) = 0 0 (x) + + n n (x) .

9.10

Show that if f is an even function, then all the odd coecients


2j1 are zero, and that if f is an odd function, then all the
even coecients 2j are zero.
The function H(x) is dened by H(x) = 1 if x > 0, and
H(x) = H(x). Construct the best polynomial approximations of degrees 0, 1 and 2 in the 2-norm to this function over
the interval (1, 1) with weight function w(x) 1. (It may not

276

9 Approximation in the 2-norm


appear very useful to consider a polynomial approximation to
a discontinuous function, but representations of such functions
by Fourier series will be familiar to most readers. Note that the
function H belongs to L2w (1, 1).)

You might also like