You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
cuando n → ∞.<br />
6. INTERPOLACIÓN POR POLINOMIOS A TROZOS 107<br />
f − qn∞ ≤ f ′′ ∞<br />
2<br />
h2 4 = f ′′ ∞ (b − a)<br />
8<br />
2<br />
n2 → 0<br />
1<br />
Ejemplo 5.18. Sea f : [−1, 1] → IR, la función de Runge f(x) = . Puede verse que<br />
1 + 25x2 f ′′ ∞ = 50. En consecuencia, aproximando por poligonales obtenemos<br />
f − qn∞ ≤ f ′′ ∞<br />
8<br />
2 2<br />
n<br />
= 25<br />
.<br />
n2 Luego, en este caso, interpolando por poligonales obtenemos una aproximación mucho mejor<br />
que la que se obtiene interpolando con un polinomio, si se utilizan los mismos puntos.<br />
Splines cúbicos.<br />
En muchos problemas interesa aproximar por funciones derivables. Esto no puede lograrse aproximando<br />
por poligonales y por lo tanto es necesario aumentar el grado de los aproximantes. Un<br />
método clásico es el correspondiente a grado tres, splines cúbicos. Vamos a ver que, de esta<br />
manera, puede obtenerse un aproximante C 2 .<br />
Dada f ∈ C[a, b] y a = x0 < x1 < x2 . . . < xn = b buscamos S tal que S, S ′ y S ′′ sean continuas<br />
en [a, b] y además se verifique<br />
S(xj) = f(xj) para 0 ≤ j ≤ n y S | [xj,xj+1]∈ P3.<br />
Por ejemplo si n = 3, como Sj ∈ P3 tenemos 4 × 3 = 12 coeficientes a determinar. Veamos<br />
cuántas condiciones se tienen que satisfacer. Tenemos que verificar,<br />
S0(x0) = f(x0), S1(x1) = f(x1), S2(x2) = f(x2),<br />
S0(x1) = S1(x1), S1(x2) = S1(x2), S2(x3) = f(x3).<br />
es decir, seis condiciones. Si además queremos S ′ y S ′′ continuas en x1, x2 tenemos cuatro<br />
condiciones más. O sea, en total diez condiciones para doce incógnitas. Una cuenta análoga<br />
puede hacerse en el caso general para ver que la cantidad de coeficientes a determinar supera en<br />
dos al número de condiciones.<br />
Luego, si hay solución habrá infinitas pues tenemos dos coeficientes para fijar arbitrariamente.<br />
Lo natural entonces es fijar S ′ (x0) y S ′ (xn) o bien S ′′ (x0) y S ′′ (xn). Elegiremos esta última<br />
opción por ser más simple.