21.06.2014 Views

Numerical Methods Contents - SAM

Numerical Methods Contents - SAM

Numerical Methods Contents - SAM

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Remark 7.0.2 (Interpolation and approximation: enabling technologies).<br />

Approximation and interpolation are useful for several numerical tasks, like integration, differentiation<br />

and computation of the solutions of differential equations, as well as for computer graphics: smooth<br />

curves and surfaces.<br />

8 Polynomial Interpolation<br />

Remark 7.0.3 (Function representation).<br />

this is a “foundations” part of the course<br />

△<br />

8.1 Polynomials<br />

Notation: Vector space of the polynomials of degree ≤ k, k ∈ N:<br />

P k := {t ↦→ α k t k + α k−1 t k−1 + · · · + α 1 t + α 0 , α j ∈ K} . (8.1.1)<br />

!<br />

General function f : D ⊂ R ↦→ K, D interval, contains an “infinite amount of information”.<br />

Terminology:<br />

the functions t ↦→ t k , k ∈ N 0 , are called monomials<br />

?<br />

How to represent f on a computer?<br />

t ↦→ α k t k + α k−1 t k−1 + · · · + α 0 = monomial representation of a polynomial.<br />

➙ Idea: parametrization, a finite number of parameters α 1,...,α n , n ∈ N, characterizes f.<br />

Ôº¾ º¼<br />

Obvious: P k is a vector space. What is its dimension?<br />

º½<br />

Special case:<br />

Representation with finite linear combination of basis functions<br />

b j : D ⊂ R ↦→ K, j = 1, ...,n:<br />

f = ∑ n<br />

j=1 α jb j , α j ∈ K .<br />

➙ f ∈ finite dimensional function space V n := Span {b 1 , ...,b n }.<br />

△<br />

✬<br />

Theorem 8.1.1 (Dimension of space of polynomials).<br />

dim P k = k + 1 and P k ⊂ C ∞ (R).<br />

✫<br />

Proof. Dimension formula by linear independence of monomials.<br />

✩<br />

✪<br />

Ôº¿½<br />

Why are polynomials important in computational mathematics ?<br />

➙<br />

➙<br />

➙<br />

Easy to compute, integrate and differentiate<br />

Vector space & algebra<br />

Analysis: Taylor polynomials & power series<br />

Remark 8.1.1 (Polynomials in Matlab).<br />

Ôº¿¼ º¼<br />

MATLAB: α k t k + α k−1 t k−1 + · · · + α 0 ➙ Vector (α k , α k−1 , ...,α 0 ) (ordered!).<br />

△<br />

Remark 8.1.2 (Horner scheme).<br />

Ôº¿¾ º½<br />

p(t) = (t · · · t(t(α n t + α n−1 ) + α n−2 ) + · · · + α 1 ) + α 0 . (8.1.2)<br />

Evaluation of a polynomial in monomial representation:<br />

Horner scheme

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

Saved successfully!

Ooh no, something went wrong!