11.08.2013 Views

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

46 Polynomial Approximation of Differential Equations<br />

In particular, for the Laguerre case, another set of Lagrange polynomials in Pn−1 will<br />

be useful in what follows. We first define η (n)<br />

0 := 0, and then we construct the Lagrange<br />

basis { ˜ l (n)<br />

j }0≤j≤n−1 w<strong>it</strong>h respect to η (n)<br />

i , 0 ≤ i ≤ n − 1. W<strong>it</strong>h a proof similar to that<br />

of theorem 3.2.1 (this time using (1.6.1) and (3.1.19)), one gets, for n ≥ 1,<br />

(3.2.11)<br />

˜ l (n)<br />

j (x) =<br />

where un = L (α)<br />

n and x ∈ [0,+∞[.<br />

3.3 The interpolation operators<br />

⎧<br />

(n − 1)! Γ(α + 2)<br />

⎪⎨<br />

− u<br />

Γ(n + α + 1)<br />

⎪⎩<br />

′ n(x) if j = 0,<br />

−<br />

x u ′ n(x)<br />

if 1 ≤ j ≤ n − 1,<br />

n un(η (n)<br />

j ) (x − η (n)<br />

j )<br />

Following the guideline of section 2.4, we introduce some new operators defined in the<br />

space of continuous functions. They are named interpolation operators. We first begin<br />

by considering the set of polynomials {uk}k∈N , orthogonal w<strong>it</strong>h respect to the weight<br />

w. Next, we choose n ≥ 1 and evaluate the zeroes ξ (n)<br />

k , 1 ≤ k ≤ n, of un. At this point,<br />

we define Iw,n : C 0 (I) → Pn−1 to be the operator mapping a continuous function f<br />

to the unique polynomial pn ∈ Pn−1 satisfying pn(ξ (n)<br />

j ) = f(ξ (n)<br />

j ), 1 ≤ j ≤ n. This<br />

polynomial is the interpolant of f w<strong>it</strong>h respect to the zeroes of un. The operator Iw,n<br />

is linear. Moreover, one can check that<br />

(3.3.1) Iw,np = p, ∀p ∈ Pn−1.<br />

A question arose in sections 2.3 and 2.4 on how to characterize a certain expression<br />

g(x,p(x)) in terms of a given polynomial p. A very simple answer is obtained when p<br />

is determined by the values attained in a set of given points. For example, we have<br />

(3.3.2) [Iw,ng(x,p(x))] x=ξ (n)<br />

j<br />

for any polynomial p ∈ Pn−1.<br />

= g(ξ (n)<br />

j ,p(ξ (n)<br />

j )), 1 ≤ j ≤ n,

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

Saved successfully!

Ooh no, something went wrong!