06.09.2021 Views

First Semester in Numerical Analysis with Julia, 2020a

First Semester in Numerical Analysis with Julia, 2020a

First Semester in Numerical Analysis with Julia, 2020a

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

CHAPTER 5. APPROXIMATION THEORY 203<br />

Section 4.3 3 . They can be obta<strong>in</strong>ed from the follow<strong>in</strong>g recursion<br />

n =1, 2,..., and they satisfy<br />

L n+1 (x) =<br />

2n +1<br />

n +1 xL n(x) −<br />

n<br />

n +1 L n−1(x),<br />

〈L n ,L n 〉 = 2<br />

2n +1 .<br />

Exercise 5.3-1:<br />

L 2 (x) are orthogonal.<br />

Show, by direct <strong>in</strong>tegration, that the Legendre polynomials L 1 (x) and<br />

Example 91 (Chebyshev polynomials). If we take w(x) =(1− x 2 ) −1/2 and [a, b] =[−1, 1],<br />

and aga<strong>in</strong> drop the orthonormal requirement <strong>in</strong> Gram-Schmidt, we obta<strong>in</strong> the follow<strong>in</strong>g<br />

orthogonal polynomials:<br />

T 0 (x) =1,T 1 (x) =x, T 2 (x) =2x 2 − 1,T 3 (x) =4x 3 − 3x, ...<br />

These polynomials are called Chebyshev polynomials and satisfy a curious identity:<br />

T n (x) =cos(n cos −1 x),n≥ 0.<br />

Chebyshev polynomials also satisfy the follow<strong>in</strong>g recursion:<br />

T n+1 (x) =2xT n (x) − T n−1 (x)<br />

for n =1, 2,...,and<br />

⎧<br />

0 if j ≠ k<br />

⎪⎨<br />

〈T j ,T k 〉 = π if j = k =0<br />

⎪⎩ π/2 if j = k>0.<br />

If we take the first n +1Legendre or Chebyshev polynomials, call them φ 0 , ..., φ n , then<br />

these polynomials form a basis for the vector space P n . In other words, they form a l<strong>in</strong>early<br />

<strong>in</strong>dependent set of functions, and any polynomial from P n can be written as a unique l<strong>in</strong>ear<br />

comb<strong>in</strong>ation of them. These statements follow from the follow<strong>in</strong>g theorem, which we will<br />

leave unproved.<br />

3 The Legendre polynomials <strong>in</strong> Section 4.3 differ from these by a constant factor. For example, <strong>in</strong> Section<br />

4.3 the third polynomial was L 2 (x) =x 2 − 1 3 , but here it is L 2(x) = 3 2 (x2 − 1 3<br />

). Observe that multiply<strong>in</strong>g these<br />

polynomials by a constant does not change their roots (what we were <strong>in</strong>terested <strong>in</strong> Gaussian quadrature),<br />

or their orthogonality.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!