11.08.2013 Views

Apunte

Apunte

Apunte

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.

108 5. INTERPOLACI ÓN<br />

Teorema 5.19. Dada f ∈ C[a, b] y a = x0 < x1 < x2 . . . < xn = b, existe un única S ∈ C 2 [a, b]<br />

tal que<br />

con S ′′ (a) = S ′′ (b) = 0.<br />

⎧<br />

⎨<br />

⎩<br />

S(xj) = f(xj) 0 ≤ j ≤ n<br />

S | [xj,xj+1]∈ P3<br />

Demostración. Para j = 0, . . . , n − 1 usaremos la notación<br />

Sj = S | [xj,xj+1] y hj = xj+1 − xj.<br />

La función S buscada debe cumplir que S ′′ es una poligonal. Por lo tanto, si S ′′ (xj) = yj, S ′′<br />

j se<br />

escribe como<br />

S ′′ xj+1 − x x − xj<br />

j (x) = yj + yj+1<br />

hj<br />

hj<br />

0 ≤ j ≤ n − 1.<br />

Veremos que es posible encontrar valores yj de tal forma que se cumplan las condiciones requeridas<br />

para S. Integrando dos veces obtenemos, para x ∈ [xj, xj+1], con 0 ≤ j ≤ n − 1<br />

Sj(x) = yj<br />

6hj<br />

(xj+1 − x) 3 + yj+1<br />

(x − xj)<br />

6hj<br />

3 + cj(x − xj) + dj(xj+1 − x) (5.7)<br />

donde cj, dj son constantes a determinar que provienen de la integración.<br />

Observemos que para cualquier elección de yj, cj y dj, S ′′ resulta continua por ser poligonal.<br />

Por lo tanto resta ver que esas constantes pueden elegirse de manera que se verifiquen las otras<br />

condiciones requeridas sobre S y S ′ .<br />

Para que S sea continua e interpole a f tenemos que elegir cj, dj tal que<br />

Sj(xj) = f(xj) y Sj(xj+1) = f(xj+1), 0 ≤ j ≤ n − 1<br />

de lo que, reemplazando en (5.7), obtenemos<br />

cj = f(xj+1)<br />

hj<br />

y por lo tanto, para cada 0 ≤ j ≤ n − 1<br />

Sj(x) = yj<br />

6hj<br />

− yj+1hj<br />

6<br />

(xj+1 − x) 3 + yj+1<br />

(x − xj)<br />

6hj<br />

3 +<br />

y dj = f(xj)<br />

hj<br />

− yjhj<br />

6<br />

<br />

f(xj+1)<br />

+ −<br />

hj<br />

yj+1hj<br />

<br />

f(xj)<br />

(x − xj) + −<br />

6<br />

hj<br />

yjhj<br />

<br />

(xj+1 − x).<br />

6

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

Saved successfully!

Ooh no, something went wrong!