Métodos de Spline Cúbico en Cálculo
Métodos de Spline Cúbico en Cálculo
TRABAJO DE INVESTIGACIÓN
SPLINE CÚBICO
Profesor: CARBAJAL ORDINOLA DAVID
Asignatura: Cálculo I
Resumen
El término "spline" hace referencia a una amplia clase de funciones que son utilizadas en
aplicaciones que requieren la interpolación de datos, o un suavizado de curvas. Los splines son
utilizados para trabajar tanto en una como en varias dimensiones. Las funciones para la
interpolación por splines normalmente se determinan como minimizadores de la aspereza
sometidas a una serie de restricciones. En este artículo nos referiremos con el término "spline" a su
versión restringida en una dimensión y polinomial, que es la más comúnmente utilizada.
Es nuestro día a día, la forma de las acciones que uno realiza no sigue una ley o regla general
puesto que no hay o existen superficies rectas ni mucho menos líneas rectas. Tanto los caminos
como las trayectorias y formas toman como ejemplo a las superficies curvas y líneas que por su
naturaleza no proceden de una ley o regla en específico.
Nuestro trabajo de investigación tiene el propósito de dar a conocer un método cuyo fin es
obtener funciones que se parecen o asemejan a la correspondencia que tendrían que tener un
determinado grupo de puntos, ello parte de la técnica denominada interpolación cuya
conceptualización se dará más adelante conforme se va desarrollando nuestro trabajo.
Tenemos la esperanza como grupo que este trabajo tenga la posibilidad de ayudar a nuestra
sociedad, sobre todo contribuir en el buen provecho de los métodos tratados y que estos sean de
mayor alcance para toda nuestra comunidad universitaria y allegados al estudio académico-
profesional de diversas ramas de nuestra actualida
1 OBJETIVOS
1.1 GENERAL
● Explicar el método de splines cúbicos y como estas sirven para obtener funciones que se
parecen o se asemejan a la correspondencia que tendrían que tener un determinado grupo
de puntos.
1.2 ESPECÍFICOS
● Calcular los coeficientes de las funciones planteadas, así como su diseño de cada una de
ellas.
1
2 LA INTERPOLACIÓN
2.1 Generalidades:
Teorema:
El problema de interpolación de Lagrange tiene una solución única.
2
Demostración:
Escribiendo pN (x) = a0 +a1x+· · ·+aN x N entonces los coeficientes son la solución del
sistema lineal a0 + a1x0 + · · · + aN x N 0 = f(x0), a0 + a1x1 + · · · + aN x N 1 = f(x1), . . . . .
. . . . a0 + a1
coincide con el valor f(x0), luego p1(x) de grado ≤ 1 que interpola en x0 y x1, por lo que p1(x)
= p0 + q1(x), tal que q1(x0) = 0, por lo que q1(x) = c1(x − x0) con c1 calculado en la igualdad
f(x0) + c1(x1 − x0) = f(x1).
Después p2(x) = p1(x) + q2(x), tal que q2(x0) = q2(x1) = 0, por lo que q2(x) = c2(x − x0)(x −
x1)...
Por consiguiente pN (x) = c0 + c1(x − x0) + c2(x − x0)(x − x1) . . . + cN (x − x0)(x − x1)· · ·
(x − xN−1), donde los coeficientes ci son las diferencias divididas de la función con la
notación
pN (x) = f[x0]
+ f[x0, x1] (x − x0)
+f[x0, x1, x2] (x − x0)(x − x1) . . .
3
Es decir, que la diferencia dividida de orden n es el cociente de la potencia n − sima del
polinomio de interpolación. Está claro entonces que la diferencia dividida de orden cero es
el valor de la función, es decir, f[xi ] = f(xi) y las diferencias de orden N ≥ 1 pueden calcularse
de las de orden N − 1 con el siguiente teorema.
Teorema
Si x0, . . . , xN son N + 1 puntos distintos dos a dos, entonces
f[x0, x1, · · · xN ] = f[x1, · · · xN ] − f[x0, · · · xN−1]
xN − x0
.
Demostración:
Utilizando este teorema se puede construir la siguiente tabla. Las dos primeras
columnas son las abscisas y los valores de la función, después cada elemento se
Se obtiene a partir de los que están a su oeste y noroeste.
x0 f[x0]
x1 f[x1] f[x0, x1]
x2 f[x2] f[x1, x2] f[x0, x1, x2]
.
.
.
.
.
.
.
.
xN f[xN ] f[xN−1, xN ] f[xN−2, xN−1, xN ] · · · f[x0, · · · , xN ] Finalmente, los elementos
diagonales son los coeficiente del polinomio interpolador en la forma de Newton.
4
3 LOS SPLINES CÚBICOS
En el caso del spline cúbico cada polinomio P(x) que construimos en [m,n] tiene grado 3.
Esto quiere decir, que va a tener la forma P(x) = ax³ + bx² + cx + d. En este caso vamos a
tener cuatro incógnitas por cada intervalo (a,b,c,d)}, y una nueva condición para cada punto
común a dos intervalos, respecto a la derivada segunda:
● Que las partes de la función a trozos P(x) pasen por ese punto. Es decir, que las dos
Pn(x) que rodean al f(x) que queremos aproximar, sean igual a f(x) en cada uno de estos
puntos.
● Que la derivada en un punto siempre coincida para ambos "lados" de la función definida
a trozos que pasa por tal punto común.
● Que la derivada segunda en un punto siempre coincida para ambos "lados" de la función
definida a trozos que pasa por tal punto común.
Como puede deducirse al compararlo con el caso de splines cuadráticos, ahora no nos va a
faltar una sino dos ecuaciones (condiciones) para el número de incógnitas que tenemos.
La forma de solucionar esto, determina el carácter de los splines cúbicos. Así, podemos
usar:
5
3.1 Método de los spline:
1. La restricción Si de S en cada intervalo [xi, xi+1], para i = 0,1, ... , n − 1, es un polinomio de grado
no superior a tres.
2. S(x) ∈ C1 ([a, b]), es decir, S = S(x) es una función diferenciable, con derivada continua, en el
intervalo [a, b].
1. La restricción Si de S en cada intervalo [xi, xi+], para i = 0,1, ... , n − 1, es un polinomio de grado
no superior a tres.
2. S(x) ∈ C2([a, b]), es decir, S = S(x) es una función dos veces diferenciable, con las dos derivadas
continuas, en el intervalo [a, b].
6
3.1.4 Spline cardinal
Este tipo de spline se le puede definir como la secuencia de las curvas individuales unidas por
una mayor a esta. El spline se especifica mediante una matriz de puntos y un parámetro de tensión.
Un spline cardinal pasa suavemente por cada punto del array; no hay esquinas afiladas ni cambios
bruscos en la tensión de la curva. Además tengamos en cuenta que los spline cardinales son solo
un subconjunto de las curvas de Hermite donde no requieren los puntos tangentes porque se
calculan de los puntos de control. La formula de las tangentes para los splines cardinales:
Ti = a * ( Pi+1 - Pi-1 )
El parámetro τ se conoce como "tensión" y afecta a la forma en que la curva presenta una
curvatura en los puntos de control interpolados de la función polinómica como podemos apreciar en
la figura 2. Nos facilita el trabajo de encontrar las tangentes en la función, entonces si damos como
función 4 puntos de control como pi-2, pi-1, pi y pi+1 y al ser polinómica presentaría este
comportamiento al derivar.
7
3.1.6 B-Spline
Es una función spline que tiene el mínimo soporte con respecto a un determinado grado,
suavidad y partición del dominio. Un teorema fundamental establece que cada función spline de un
determinado grado, suavidad y partición del dominio, se puede representar como una combinación
lineal de B-splines del mismo grado y suavidad, y sobre la misma partición.
con
8
Propiedades de B-spline
● Toda B-spline va a pasar por el primer y último punto
Nota:
Cuando el número de puntos de control de Boor es el mismo que el grado, la B-Spline degenera en
una curva de Bézier. La forma de las funciones base es determinada por la posición de los nodos.
Escalar o trasladar el vector de nodo no altera las funciones de base.
El spline está contenido en el casco convexo de sus puntos de control. Una B-spline básica de grado
n
9
3.1.6.2.2 B-spline lineal
La B-spline lineal se define en dos tramos de nodo consecutivos y es continua sobre los nodos,
pero no diferenciable.
donde i es el número de puntos de control y t es un parámetro global dados los valores de los nodos.
Esta formulación expresa una curva B-spline como una combinación lineal de funciones B-spline
básicas, de ahí el nombre.
Hay dos tipos de B-spline - uniforme y no uniforme. Una B-spline no uniforme es una curva donde
los intervalos entre los puntos sucesivos de control no son, o no necesariamente son, iguales (el
vector de nodos de espacios de nodo interiores no son iguales). Una forma común es donde los
intervalos se reducen sucesivamente a cero, interpolando los puntos de control.
10
3.1.6.2.5 B-spline cúbica uniforme
La B-spline cúbica con vector-nodo uniforme es la forma más usual de B-spline. La función base
puede ser fácilmente calculada, y es igual para cada segmento, en este caso. Puesto en forma de
matriz, esto es:
4 TRABAJO DE INVESTIGACIÓN:
Dados los siguientes puntos:
Determinaremos el Spline cúbico natural de tipo C2 interpolante de los pares de puntos anteriores.
11
4.1 Ecuaciones para resolver el spline cúbico:
(Lo señalado con una flecha son las ecuaciones planteadas para hallar todas las variables)
12
4.2 Primeras ecuaciones:
13
𝒂𝟒 (𝟏. 𝟓𝟎)𝟑 + 𝒃𝟒 (𝟏. 𝟓𝟎)𝟐 + 𝒄𝟒 (𝟏. 𝟓𝟎) + 𝒅𝟒 = 𝟎. 𝟓𝟒𝟕𝟒, 𝒙 = 𝟏. 𝟓𝟎
→ 𝟑. 𝟑𝟕𝟓𝒂𝟒 + 𝟐. 𝟐𝟓𝒃𝟒 + 𝟏. 𝟓𝒄𝟒 + 𝒅𝟒 = 𝟎. 𝟓𝟒𝟕𝟒
𝒂𝟓 (𝟏. 𝟓𝟎)𝟑 + 𝒃𝟓 (𝟏. 𝟓𝟎)𝟐 + 𝒄𝟓 (𝟏. 𝟓𝟎) + 𝒅𝟓 = 𝟎. 𝟓𝟒𝟕𝟒, 𝒙 = 𝟏. 𝟓𝟎
→ 𝟑. 𝟑𝟕𝟓𝒂𝟓 + 𝟐. 𝟐𝟓𝒃𝟓 + 𝟏. 𝟓𝒄𝟓 + 𝒅𝟓 = 𝟎. 𝟓𝟒𝟕𝟒
𝒂𝟏𝟎 (𝟐. 𝟏𝟎)𝟑 + 𝒃𝟏𝟎 (𝟐. 𝟏𝟎)𝟐 + 𝒄𝟏𝟎 (𝟐. 𝟏𝟎) + 𝒅𝟏𝟎 = 𝟎. 𝟓𝟏𝟕𝟎, 𝒙 = 𝟐. 𝟏𝟎
14
→ 𝟗. 𝟐𝟔𝟏𝒂𝟏𝟎 + 𝟒. 𝟒𝟏𝒃𝟏𝟎 + 𝟐. 𝟏𝒄𝟏𝟎 + 𝒅𝟏𝟎= 𝟎. 𝟓𝟏𝟕𝟎
𝟑 𝟐
𝒂𝟏𝟏 (𝟐. 𝟏𝟎) + 𝒃𝟏𝟏 (𝟐. 𝟏𝟎) + 𝒄𝟏𝟏 (𝟐. 𝟏𝟎) + 𝒅𝟏𝟏 = 𝟎. 𝟓𝟏𝟕𝟎, 𝒙 = 𝟐. 𝟏𝟎
→ 𝟗. 𝟐𝟔𝟏𝒂𝟏𝟏 + 𝟒. 𝟒𝟏𝒃𝟏𝟏 + 𝟐. 𝟏𝒄𝟏𝟏 + 𝒅𝟏𝟏 = 𝟎. 𝟓𝟏𝟕𝟎
𝒂𝟏𝟏 (𝟐. 𝟐𝟎)𝟑 + 𝒃𝟏𝟏 (𝟐. 𝟐𝟎)𝟐 + 𝒄𝟏𝟏 (𝟐. 𝟐𝟎) + 𝒅𝟏𝟏 = 𝟎. 𝟓𝟏𝟒𝟕, 𝒙 = 𝟐. 𝟐𝟎
→ 𝟏𝟎. 𝟔𝟒𝟖𝒂𝟏𝟏 + 𝟒. 𝟖𝟒𝒃𝟏𝟏 + 𝟐. 𝟐𝒄𝟏𝟏 + 𝒅𝟏𝟏 = 𝟎. 𝟓𝟏𝟒𝟕
𝒂𝟏𝟐 (𝟐. 𝟐𝟎)𝟑 + 𝒃𝟏𝟐 (𝟐. 𝟐𝟎)𝟐 + 𝒄𝟏𝟐 (𝟐. 𝟐𝟎) + 𝒅𝟏𝟐 = 𝟎. 𝟓𝟏𝟒𝟕, 𝒙 = 𝟐. 𝟐𝟎
→ 𝟏𝟎. 𝟔𝟒𝟖𝒂𝟏𝟐 + 𝟒. 𝟖𝟒𝒃𝟏𝟐 + 𝟐. 𝟐𝒄𝟏𝟐 + 𝒅𝟏𝟐 = 𝟎. 𝟓𝟏𝟒𝟕
𝒂𝟏𝟐 (𝟐. 𝟑𝟎)𝟑 + 𝒃𝟏𝟐 (𝟐. 𝟑𝟎)𝟐 + 𝒄𝟏𝟐 (𝟐. 𝟑𝟎) + 𝒅𝟏𝟐 = 𝟎. 𝟓𝟏𝟐𝟕, 𝒙 = 𝟐. 𝟑𝟎
→ 𝟏𝟐. 𝟏𝟔𝟕𝒂𝟏𝟐 + 𝟓. 𝟐𝟗𝒃𝟏𝟐 + 𝟐. 𝟑𝒄𝟏𝟐 + 𝒅𝟏𝟐 = 𝟎. 𝟓𝟏𝟐𝟕
𝒂𝟏𝟑 (𝟐. 𝟑𝟎) + 𝒃𝟏𝟑 (𝟐. 𝟑𝟎)𝟐 + 𝒄𝟏𝟑 (𝟐. 𝟑𝟎) + 𝒅𝟏𝟑 = 𝟎. 𝟓𝟏𝟐𝟕,
𝟑
𝒙 = 𝟐. 𝟑𝟎
→ 𝟏𝟐. 𝟏𝟔𝟕𝒂𝟏𝟑 + 𝟓. 𝟐𝟗𝒃𝟏𝟑 + 𝟐. 𝟑𝒄𝟏𝟑 + 𝒅𝟏𝟑 = 𝟎. 𝟓𝟏𝟐𝟕
𝒂𝟏𝟑 (𝟐. 𝟒𝟎)𝟑 + 𝒃𝟏𝟑 (𝟐. 𝟒𝟎)𝟐 + 𝒄𝟏𝟑 (𝟐. 𝟒𝟎) + 𝒅𝟏𝟑 = 𝟎. 𝟓𝟏𝟎𝟗, 𝒙 = 𝟐. 𝟒𝟎
→ 𝟏𝟑. 𝟖𝟐𝟒𝒂𝟏𝟑 + 𝟓. 𝟕𝟔𝒃𝟏𝟑 + 𝟐. 𝟒𝒄𝟏𝟑 + 𝒅𝟏𝟑 = 𝟎. 𝟓𝟏𝟎𝟗
𝒂𝟏𝟒 (𝟐. 𝟒𝟎) + 𝒃𝟏𝟒 (𝟐. 𝟒𝟎)𝟐 + 𝒄𝟏𝟒 (𝟐. 𝟒𝟎) + 𝒅𝟏𝟒 = 𝟎. 𝟓𝟏𝟎𝟗,
𝟑
𝒙 = 𝟐. 𝟒𝟎
→ 𝟏𝟑. 𝟖𝟐𝟒𝒂𝟏𝟒 + 𝟓. 𝟕𝟔𝒃𝟏𝟒 + 𝟐. 𝟒𝒄𝟏𝟒 + 𝒅𝟏𝟒 = 𝟎. 𝟓𝟏𝟎𝟗
𝒂𝟏𝟒 (𝟐. 𝟓𝟎)𝟑 + 𝒃𝟏𝟒 (𝟐. 𝟓𝟎)𝟐 + 𝒄𝟏𝟒 (𝟐. 𝟓𝟎) + 𝒅𝟏𝟒 = 𝟎. 𝟓𝟎𝟗𝟓, 𝒙 = 𝟐. 𝟓𝟎
→ 𝟏𝟓. 𝟔𝟐𝟓𝒂𝟏𝟒 + 𝟔. 𝟐𝟓𝒃𝟏𝟒 + 𝟐. 𝟓𝒄𝟏𝟒 + 𝒅𝟏𝟒 = 𝟎. 𝟓𝟎𝟗𝟓
𝒂𝟏𝟓 (𝟐. 𝟓𝟎) + 𝒃𝟏𝟓 (𝟐. 𝟓𝟎)𝟐 + 𝒄𝟏𝟓 (𝟐. 𝟓𝟎) + 𝒅𝟏𝟓 = 𝟎. 𝟓𝟎𝟗𝟓,
𝟑
𝒙 = 𝟐. 𝟓𝟎
→ 𝟏𝟓. 𝟔𝟐𝟓𝒂𝟏𝟓 + 𝟔. 𝟐𝟓𝒃𝟏𝟓 + 𝟐. 𝟓𝒄𝟏𝟓 + 𝒅𝟏𝟓 = 𝟎. 𝟓𝟎𝟗𝟓
𝒂𝟏𝟓 (𝟐. 𝟔𝟎)𝟑 + 𝒃𝟏𝟓 (𝟐. 𝟔𝟎)𝟐 + 𝒄𝟏𝟓 (𝟐. 𝟔𝟎) + 𝒅𝟏𝟓 =, 𝟎. 𝟓𝟎𝟖𝟒 𝒙 = 𝟐. 𝟔𝟎
→ 𝟏𝟕. 𝟓𝟕𝟔𝒂𝟏𝟓 + 𝟔. 𝟕𝟔𝒃𝟏𝟓 + 𝟐. 𝟔𝒄𝟏𝟓 + 𝒅𝟏𝟓 = 𝟎. 𝟓𝟎𝟖𝟒
𝒂𝟏𝟔 (𝟐. 𝟔𝟎) + 𝒃𝟏𝟔 (𝟐. 𝟔𝟎)𝟐 + 𝒄𝟏𝟔 (𝟐. 𝟔𝟎) + 𝒅𝟏𝟔 =, 𝟎. 𝟓𝟎𝟖𝟒
𝟑
𝒙 = 𝟐. 𝟔𝟎
→ 𝟏𝟕. 𝟓𝟕𝟔𝒂𝟏𝟔 + 𝟔. 𝟕𝟔𝒃𝟏𝟔 + 𝟐. 𝟔𝒄𝟏𝟔 + 𝒅𝟏𝟔 = 𝟎. 𝟓𝟎𝟖𝟒
𝒂𝟏𝟔 (𝟐. 𝟕𝟎)𝟑 + 𝒃𝟏𝟔 (𝟐. 𝟕𝟎)𝟐 + 𝒄𝟏𝟔 (𝟐. 𝟕𝟎) + 𝒅𝟏𝟔 = 𝟎. 𝟓𝟎𝟕𝟐, 𝒙 = 𝟐. 𝟕𝟎
→ 𝟏𝟗. 𝟔𝟖𝟑𝒂𝟏𝟔 + 𝟕. 𝟐𝟗𝒃𝟏𝟔 + 𝟐. 𝟕𝒄𝟏𝟔 + 𝒅𝟏𝟔 = 𝟎. 𝟓𝟎𝟕𝟐
15
𝒂𝟏𝟕 (𝟐. 𝟕𝟎)𝟑 + 𝒃𝟏𝟕 (𝟐. 𝟕𝟎)𝟐 + 𝒄𝟏𝟕 (𝟐. 𝟕𝟎) + 𝒅𝟏𝟕 = 𝟎. 𝟓𝟎𝟕𝟐, 𝒙 = 𝟐. 𝟕𝟎
→ 𝟏𝟗. 𝟔𝟖𝟑𝒂𝟏𝟕 + 𝟕. 𝟐𝟗𝒃𝟏𝟕 + 𝟐. 𝟕𝒄𝟏𝟕 + 𝒅𝟏𝟕 = 𝟎. 𝟓𝟎𝟕𝟐
𝒂𝟏𝟕 (𝟐. 𝟖𝟎)𝟑 + 𝒃𝟏𝟕 (𝟐. 𝟖𝟎)𝟐 + 𝒄𝟏𝟕 (𝟐. 𝟖𝟎) + 𝒅𝟏𝟕 = 𝟎. 𝟓𝟎𝟔𝟐, 𝒙 = 𝟐. 𝟖𝟎
→ 𝟐𝟏. 𝟗𝟓𝟐𝒂𝟏𝟕 + 𝟕. 𝟖𝟒𝒃𝟏𝟕 + 𝟐. 𝟖𝒄𝟏𝟕 + 𝒅𝟏𝟕 = 𝟎. 𝟓𝟎𝟔𝟐
𝒂𝟏𝟖 (𝟐. 𝟖𝟎)𝟑 + 𝒃𝟏𝟖 (𝟐. 𝟖𝟎)𝟐 + 𝒄𝟏𝟖 (𝟐. 𝟖𝟎) + 𝒅𝟏𝟖 = 𝟎. 𝟓𝟎𝟔𝟐, 𝒙 = 𝟐. 𝟖𝟎
→ 𝟐𝟏. 𝟗𝟓𝟐𝒂𝟏𝟖 + 𝟕. 𝟖𝟒𝒃𝟏𝟖 + 𝟐. 𝟖𝒄𝟏𝟖 + 𝒅𝟏𝟖 = 𝟎. 𝟓𝟎𝟔𝟐
𝒂𝟏𝟖 (𝟐. 𝟗𝟎)𝟑 + 𝒃𝟏𝟖 (𝟐. 𝟗𝟎)𝟐 + 𝒄𝟏𝟖 (𝟐. 𝟗𝟎) + 𝒅𝟏𝟖 = 𝟎. 𝟓𝟎𝟓𝟕, 𝒙 = 𝟐. 𝟗𝟎
→ 𝟐𝟒. 𝟑𝟖𝟗𝒂𝟏𝟖 + 𝟖. 𝟒𝟏𝒃𝟏𝟖 + 𝟐. 𝟗𝒄𝟏𝟖 + 𝒅𝟏𝟖 = 𝟎. 𝟓𝟎𝟓𝟕
𝒂𝟏𝟗 (𝟐. 𝟗𝟎) + 𝒃𝟏𝟗 (𝟐. 𝟗𝟎)𝟐 + 𝒄𝟏𝟗 (𝟐. 𝟗𝟎) + 𝒅𝟏𝟗 = 𝟎. 𝟓𝟎𝟓𝟕,
𝟑
𝒙 = 𝟐. 𝟗𝟎
→ 𝟐𝟒. 𝟑𝟖𝟗𝒂𝟏𝟗 + 𝟖. 𝟒𝟏𝒃𝟏𝟗 + 𝟐. 𝟗𝒄𝟏𝟗 + 𝒅𝟏𝟗 = 𝟎. 𝟓𝟎𝟓𝟕
𝒂𝟏𝟗 (𝟑. 𝟎𝟎)𝟑 + 𝒃𝟏𝟗 (𝟑. 𝟎𝟎)𝟐 + 𝒄𝟏𝟗 (𝟑. 𝟎𝟎) + 𝒅𝟏𝟗 = 𝟎. 𝟓𝟎𝟓𝟐, 𝒙 = 𝟑. 𝟎𝟎
→ 𝟐𝟕𝒂𝟏𝟗 + 𝟗𝒃𝟏𝟗 + 𝟑𝒄𝟏𝟗 + 𝒅𝟏𝟗 = 𝟎. 𝟓𝟎𝟓𝟐
𝒂𝟐𝟎 (𝟑. 𝟎𝟎) + 𝒃𝟐𝟎 (𝟑. 𝟎𝟎)𝟐 + 𝒄𝟐𝟎 (𝟑. 𝟎𝟎) + 𝒅𝟐𝟎 = 𝟎. 𝟓𝟎𝟓𝟐,
𝟑
𝒙 = 𝟑. 𝟎𝟎
→ 𝟐𝟕𝒂𝟐𝟎 + 𝟗𝒃𝟐𝟎 + 𝟑𝒄𝟐𝟎 + 𝒅𝟐𝟎 = 𝟎. 𝟓𝟎𝟓𝟐
𝒂𝟐𝟎 (𝟑. 𝟏𝟎)𝟑 + 𝒃𝟐𝟎 (𝟑. 𝟏𝟎)𝟐 + 𝒄𝟐𝟎 (𝟑. 𝟏𝟎) + 𝒅𝟐𝟎 = 𝟎. 𝟓𝟎𝟒𝟓, 𝒙 = 𝟑. 𝟏𝟎
→ 𝟐𝟗. 𝟕𝟗𝟏𝒂𝟐𝟎 + 𝟗. 𝟔𝟏𝒃𝟐𝟎 + 𝟑. 𝟏𝒄𝟐𝟎 + 𝒅𝟐𝟎 = 𝟎. 𝟓𝟎𝟒𝟓
𝒂𝟐𝟏 (𝟑. 𝟏𝟎)𝟑 + 𝒃𝟐𝟏 (𝟑. 𝟏𝟎)𝟐 + 𝒄𝟐𝟏 (𝟑. 𝟏𝟎) + 𝒅𝟐𝟏 = 𝟎. 𝟓𝟎𝟒𝟓, 𝒙 = 𝟑. 𝟏𝟎
→ 𝟐𝟗. 𝟕𝟗𝟏𝒂𝟐𝟏 + 𝟗. 𝟔𝟏𝒃𝟐𝟏 + 𝟑. 𝟏𝒄𝟐𝟏 + 𝒅𝟐𝟏 = 𝟎. 𝟓𝟎𝟒𝟓
𝒂𝟐𝟏 (𝟑. 𝟐𝟎)𝟑 + 𝒃𝟐𝟏 (𝟑. 𝟐𝟎)𝟐 + 𝒄𝟐𝟏 (𝟑. 𝟐𝟎) + 𝒅𝟐𝟏 = 𝟎. 𝟓𝟎𝟑𝟕, 𝒙 = 𝟑. 𝟐𝟎
→ 𝟑𝟐. 𝟕𝟔𝟖𝒂𝟐𝟏 + 𝟏𝟎. 𝟐𝟒𝒃𝟐𝟏 + 𝟑. 𝟐𝒄𝟐𝟏 + 𝒅𝟐𝟏 = 𝟎. 𝟓𝟎𝟑𝟕
𝒂𝟐𝟐 (𝟑. 𝟐𝟎)𝟑 + 𝒃𝟐𝟐 (𝟑. 𝟐𝟎)𝟐 + 𝒄𝟐𝟐 (𝟑. 𝟐𝟎) + 𝒅𝟐𝟐 = 𝟎. 𝟓𝟎𝟑𝟕, 𝒙 = 𝟑. 𝟐𝟎
→ 𝟑𝟐. 𝟕𝟔𝟖𝒂𝟐𝟐 + 𝟏𝟎. 𝟐𝟒𝒃𝟐𝟐 + 𝟑. 𝟐𝒄𝟐𝟐 + 𝒅𝟐𝟐 = 𝟎. 𝟓𝟎𝟑𝟕
𝒂𝟐𝟐 (𝟑. 𝟑𝟎) + 𝒃𝟐𝟐 (𝟑. 𝟑𝟎)𝟐 + 𝒄𝟐𝟐 (𝟑. 𝟑𝟎) + 𝒅𝟐𝟐 = 𝟎. 𝟓𝟎𝟑𝟕,
𝟑
𝒙 = 𝟑. 𝟑𝟎
→ 𝟑𝟓. 𝟗𝟑𝟕𝒂𝟐𝟐 + 𝟏𝟎. 𝟖𝟗𝒃𝟐𝟐 + 𝟑. 𝟑𝒄𝟐𝟐 + 𝒅𝟐𝟐 = 𝟎. 𝟓𝟎𝟑𝟕
16
𝒂𝟐𝟑 (𝟑. 𝟑𝟎)𝟑 + 𝒃𝟐𝟑 (𝟑. 𝟑𝟎)𝟐 + 𝒄𝟐𝟑 (𝟑. 𝟑𝟎) + 𝒅𝟐𝟑 = 𝟎. 𝟓𝟎𝟑𝟕, 𝒙 = 𝟑. 𝟑𝟎
→ 𝟑𝟓. 𝟗𝟑𝟕𝒂𝟐𝟑 + 𝟏𝟎. 𝟖𝟗𝒃𝟐𝟑 + 𝟑. 𝟑𝒄𝟐𝟑 + 𝒅𝟐𝟑 = 𝟎. 𝟓𝟎𝟑𝟕
𝒂𝟐𝟒 (𝟑. 𝟒𝟎)𝟑 + 𝒃𝟐𝟒 (𝟑. 𝟒𝟎)𝟐 + 𝒄𝟐𝟒 (𝟑. 𝟒𝟎) + 𝒅𝟐𝟒 = 𝟎. 𝟓𝟎𝟑𝟒, 𝒙 = 𝟑. 𝟒𝟎
→ 𝟑𝟗. 𝟑𝟎𝟒𝒂𝟐𝟑 + 𝟏𝟏. 𝟓𝟔𝒃𝟐𝟑 + 𝟑. 𝟒𝒄𝟐𝟑 + 𝒅𝟐𝟑 = 𝟎. 𝟓𝟎𝟑𝟒
𝒂𝟐𝟒 (𝟑. 𝟒𝟎)𝟑 + 𝒃𝟐𝟒 (𝟑. 𝟒𝟎)𝟐 + 𝒄𝟐𝟒 (𝟑. 𝟒𝟎) + 𝒅𝟐𝟒 = 𝟎. 𝟓𝟎𝟑𝟒, 𝒙 = 𝟑. 𝟒𝟎
→ 𝟑𝟗. 𝟑𝟎𝟒𝒂𝟐𝟒 + 𝟏𝟏. 𝟓𝟔𝒃𝟐𝟒 + 𝟑. 𝟒𝒄𝟐𝟒 + 𝒅𝟐𝟒 = 𝟎. 𝟓𝟎𝟑𝟒
𝒂𝟐𝟒 (𝟑. 𝟓𝟎)𝟑 + 𝒃𝟐𝟒 (𝟑. 𝟓𝟎)𝟐 + 𝒄𝟐𝟒 (𝟑. 𝟓𝟎) + 𝒅𝟐𝟒 = 𝟎. 𝟓𝟎𝟐𝟗, 𝒙 = 𝟑. 𝟓𝟎
→ 𝟒𝟐. 𝟖𝟕𝟓𝒂𝟐𝟒 + 𝟏𝟐. 𝟐𝟓𝒃𝟐𝟒 + 𝟑. 𝟓𝒄𝟐𝟒 + 𝒅𝟐𝟒 = 𝟎. 𝟓𝟎𝟐𝟗
𝒂𝟐𝟓 (𝟑. 𝟓𝟎) + 𝒃𝟐𝟓 (𝟑. 𝟓𝟎)𝟐 + 𝒄𝟐𝟓 (𝟑. 𝟓𝟎) + 𝒅𝟐𝟓 = 𝟎. 𝟓𝟎𝟐𝟗,
𝟑
𝒙 = 𝟑. 𝟓𝟎
→ 𝟒𝟐. 𝟖𝟕𝟓𝒂𝟐𝟓 + 𝟏𝟐. 𝟐𝟓𝒃𝟐𝟓 + 𝟑. 𝟓𝒄𝟐𝟓 + 𝒅𝟐𝟓 = 𝟎. 𝟓𝟎𝟐𝟗
𝒂𝟐𝟓 (𝟑. 𝟔𝟎)𝟑 + 𝒃𝟐𝟓 (𝟑. 𝟔𝟎)𝟐 + 𝒄𝟐𝟓 (𝟑. 𝟔𝟎) + 𝒅𝟐𝟓 = 𝟎. 𝟓𝟎𝟐𝟒, 𝒙 = 𝟑. 𝟔𝟎
→ 𝟒𝟔. 𝟔𝟓𝟔𝒂𝟐𝟓 + 𝟏𝟐. 𝟗𝟔𝒃𝟐𝟓 + 𝟑. 𝟔𝒄𝟐𝟓 + 𝒅𝟐𝟓 = 𝟎. 𝟓𝟎𝟐𝟒
𝒂𝟐𝟔 (𝟑. 𝟔𝟎) + 𝒃𝟐𝟔 (𝟑. 𝟔𝟎)𝟐 + 𝒄𝟐𝟔 (𝟑. 𝟔𝟎) + 𝒅𝟐𝟔 = 𝟎. 𝟓𝟎𝟐𝟒,
𝟑
𝒙 = 𝟑. 𝟔𝟎
→ 𝟒𝟔. 𝟔𝟓𝟔𝒂𝟐𝟔 + 𝟏𝟐. 𝟗𝟔𝒃𝟐𝟔 + 𝟑. 𝟔𝒄𝟐𝟔 + 𝒅𝟐𝟔 = 𝟎. 𝟓𝟎𝟐𝟒
𝒂𝟐𝟔 (𝟑. 𝟕𝟎)𝟑 + 𝒃𝟐𝟔 (𝟑. 𝟕𝟎)𝟐 + 𝒄𝟐𝟔 (𝟑. 𝟕𝟎) + 𝒅𝟐𝟔 = 𝟎. 𝟓𝟎𝟐𝟐, 𝒙 = 𝟑. 𝟕𝟎
→ 𝟓𝟎. 𝟔𝟓𝟑𝒂𝟐𝟔 + 𝟏𝟑. 𝟔𝟗𝒃𝟐𝟔 + 𝟑. 𝟕𝒄𝟐𝟔 + 𝒅𝟐𝟔 = 𝟎. 𝟓𝟎𝟐𝟐
𝒂𝟐𝟕 (𝟑. 𝟕𝟎) + 𝒃𝟐𝟕 (𝟑. 𝟕𝟎)𝟐 + 𝒄𝟐𝟕 (𝟑. 𝟕𝟎) + 𝒅𝟐𝟕 = 𝟎. 𝟓𝟎𝟐𝟐,
𝟑
𝒙 = 𝟑. 𝟕𝟎
→ 𝟓𝟎. 𝟔𝟓𝟑𝒂𝟐𝟕 + 𝟏𝟑. 𝟔𝟗𝒃𝟐𝟕 + 𝟑. 𝟕𝒄𝟐𝟕 + 𝒅𝟐𝟕 = 𝟎. 𝟓𝟎𝟐𝟐
𝒂𝟐𝟕 (𝟑. 𝟖𝟎)𝟑 + 𝒃𝟐𝟕 (𝟑. 𝟖𝟎)𝟐 + 𝒄𝟐𝟕 (𝟑. 𝟖𝟎) + 𝒅𝟐𝟕 = 𝟎. 𝟓𝟎𝟐𝟏, 𝒙 = 𝟑. 𝟖𝟎
→ 𝟓𝟒. 𝟖𝟕𝟐𝒂𝟐𝟕 + 𝟏𝟒. 𝟒𝟒𝒃𝟐𝟕 + 𝟑. 𝟖𝒄𝟐𝟕 + 𝒅𝟐𝟕 = 𝟎. 𝟓𝟎𝟐𝟏
𝒂𝟐𝟖 (𝟑. 𝟖𝟎)𝟑 + 𝒃𝟐𝟖 (𝟑. 𝟖𝟎)𝟐 + 𝒄𝟐𝟖 (𝟑. 𝟖𝟎) + 𝒅𝟐𝟖 = 𝟎. 𝟓𝟎𝟐𝟏, 𝒙 = 𝟑. 𝟖𝟎
→ 𝟓𝟒. 𝟖𝟕𝟐𝒂𝟐𝟖 + 𝟏𝟒. 𝟒𝟒𝒃𝟐𝟖 + 𝟑. 𝟖𝒄𝟐𝟖 + 𝒅𝟐𝟖 = 𝟎. 𝟓𝟎𝟐𝟏
𝒂𝟐𝟖 (𝟑. 𝟗𝟎)𝟑 + 𝒃𝟐𝟖 (𝟑. 𝟗𝟎)𝟐 + 𝒄𝟐𝟖 (𝟑. 𝟗𝟎) + 𝒅𝟐𝟖 = 𝟎. 𝟓𝟎𝟐𝟏, 𝒙 = 𝟑. 𝟗𝟎
→ 𝟓𝟗. 𝟑𝟏𝟗𝒂𝟐𝟖 + 𝟏𝟓. 𝟐𝟏𝒃𝟐𝟖 + 𝟑. 𝟗𝒄𝟐𝟖 + 𝒅𝟐𝟖 = 𝟎. 𝟓𝟎𝟐𝟏
17
𝒂𝟐𝟗 (𝟑. 𝟗𝟎)𝟑 + 𝒃𝟐𝟗 (𝟑. 𝟗𝟎)𝟐 + 𝒄𝟐𝟗 (𝟑. 𝟗𝟎) + 𝒅𝟐𝟗 = 𝟎. 𝟓𝟎𝟐𝟏, 𝒙 = 𝟑. 𝟗𝟎
→ 𝟓𝟗. 𝟑𝟏𝟗𝒂𝟐𝟗 + 𝟏𝟓. 𝟐𝟏𝒃𝟐𝟗 + 𝟑. 𝟗𝒄𝟐𝟗 + 𝒅𝟐𝟗 = 𝟎. 𝟓𝟎𝟐𝟏
𝒂𝟐𝟗 (𝟒. 𝟎𝟎)𝟑 + 𝒃𝟐𝟗 (𝟒. 𝟎𝟎)𝟐 + 𝒄𝟐𝟗 (𝟒. 𝟎𝟎) + 𝒅𝟐𝟗 = 𝟎. 𝟓𝟎𝟏𝟕, 𝒙 = 𝟒. 𝟎𝟎
→ 𝟔𝟒𝒂𝟐𝟗 + 𝟏𝟔𝒃𝟐𝟗 + 𝟒𝒄𝟐𝟗 + 𝒅𝟐𝟗 = 𝟎. 𝟓𝟎𝟏𝟕
18
𝟑𝒂𝟔 (𝟏. 𝟕𝟎)𝟐 + 𝟐𝒃𝟔 (𝟏. 𝟕𝟎) + 𝒄𝟔
= 𝟑𝒂𝟕 (𝟏. 𝟕𝟎)𝟐 + 𝟐𝒃𝟕 (𝟏. 𝟕𝟎) + 𝒄𝟕 [𝟏. 𝟕𝟎]
→ 𝟖. 𝟔𝟕𝒂𝟔 + 𝟑. 𝟒𝒃𝟔 + 𝒄𝟔 = 𝟖. 𝟔𝟕𝒂𝟕 + 𝟑. 𝟒𝒃𝟕 + 𝒄𝟕
19
𝟑𝒂𝟏𝟓 (𝟐. 𝟔𝟎)𝟐 + 𝟐𝒃𝟏𝟓 (𝟐. 𝟔𝟎) + 𝒄𝟏𝟓
= 𝟑𝒂𝟏𝟔 (𝟐. 𝟔𝟎)𝟐 + 𝟐𝒃𝟏𝟔 (𝟐. 𝟔𝟎) + 𝒄𝟏𝟔 [𝟐. 𝟔𝟎]
→ 𝟐𝟎. 𝟐𝟖𝒂𝟏𝟓 + 𝟓. 𝟐𝒃𝟏𝟓 + 𝒄𝟏𝟓 = 𝟐𝟎. 𝟐𝟖𝒂𝟏𝟔 + 𝟓. 𝟐𝒃𝟏𝟔 + 𝒄𝟏𝟔
20
𝟑𝒂𝟐𝟒 (𝟑. 𝟓𝟎)𝟐 + 𝟐𝒃𝟐𝟒 (𝟑. 𝟓𝟎) + 𝒄𝟐𝟒
= 𝟑𝒂𝟐𝟓 (𝟑. 𝟓𝟎)𝟐 + 𝟐𝒃𝟐𝟓 (𝟑. 𝟓𝟎) + 𝒄𝟐𝟓 [𝟑. 𝟓𝟎]
→ 𝟑𝟔. 𝟕𝟓𝒂𝟐𝟒 + 𝟕𝒃𝟐𝟒 + 𝒄𝟐𝟒 = 𝟑𝟔. 𝟕𝟓𝒂𝟐𝟓 + 𝟕𝒃𝟐𝟓 + 𝒄𝟐𝟓
𝑺´´(𝒙) = 𝟔𝒂𝒙 + 𝟐𝒃
21
𝟔𝒂𝟑 (𝟏. 𝟒𝟎) + 𝟐𝒃𝟑 = 𝟔𝒂𝟒 (𝟏. 𝟒𝟎) + 𝟐𝒃𝟒
→ 𝟖. 𝟒𝒂𝟑 + 𝟐𝒃𝟑 = 𝟖. 𝟒𝒂𝟒 + 𝟐𝒃𝟒
22
→ 𝟏𝟓𝒂𝟏𝟒 + 𝟐𝒃𝟏𝟒 = 𝟏𝟓𝒂𝟏𝟓 + 𝟐𝒃𝟏𝟓
23
𝟔𝒂𝟐𝟔 (𝟑. 𝟕𝟎) + 𝟐𝒃𝟐𝟔 = 𝟔𝒂𝟐𝟕 (𝟑. 𝟕𝟎) + 𝟐𝒃𝟐𝟕
→ 𝟐𝟐. 𝟐𝒂𝟐𝟔 + 𝟐𝒃𝟐𝟔 = 𝟐𝟐. 𝟐𝒂𝟐𝟕 + 𝟐𝒃𝟐𝟕
𝑺´´(𝒙𝟎 ) = 𝟎
𝟔𝒂𝟎 (𝟏. 𝟎𝟎) + 𝟐𝒃𝟎 = 𝟎
→ 𝟔𝒂𝟎 + 𝟐𝒃𝟎 = 𝟎
𝑺´´(𝒙𝟐𝟗 ) = 𝟎
𝟔𝒂𝟐𝟗 (𝟒. 𝟎𝟎) + 𝟐𝒃𝟐𝟗 = 𝟎
→ 𝟐𝟒𝒂𝟐𝟗 + 𝟐𝒃𝟐𝟗 = 𝟎
24
Al ponerlo en una matriz 120x120 y otra de 1x120, usar la eliminación de Gauss
en Matlab para calcular las incógnitas:
a0 0.867506642149556
b0 -2.602519926448730
c0 2.366844860027300
d0 -0.006831575728055
a1 0.062466789252610
b1 0.054111588111260
c1 -0.555449805988756
d1 1.064676468477790
a2 -0.494026155822210
b2 2.057486190380620
c2 -2.975732805145340
d2 2.045776449287090
a3 -0.209709809302328
b3 0.948652438953057
c3 -1.534248928289530
d3 1.421133435982890
a4 0.050055638102474
b4 -0.142362440147111
c4 -0.006828097549295
d4 0.708337048304156
a5 0.052434130533805
b5 -0.153065656088092
c5 0.009226726362178
d5 0.701915118739567
a6 0.102743357268419
b6 -0.394549944414027
c6 0.395601587683718
25
d6 0.495848526034744
a7 -0.586080195825631
b7 3.118450176365630
c7 -5.576498617641110
d7 3.880038642385790
a8 0.441577426034422
b8 -2.430900981678290
c8 4.412333466839360
d8 -2.113260608302360
a9 -0.180229508310703
b9 1.113398544088290
c9 -2.321835632118470
d9 2.151713154370550
a10 -0.020659392793055
b10 0.155977850982042
c10 -0.406994245907527
d10 0.875152230229677
a11 -0.037132920516382
b11 0.259761075638484
c11 -0.624939017687054
d11 1.027713570475550
a12 -0.030808925142413
b12 0.218022706170192
c12 -0.533114604858095
d12 0.960375667733051
a13 0.060368621086136
26
b13 -0.411102362806060
c13 0.913873053788946
d13 -0.148981537227981
a14 -0.010665559201712
b14 0.100343735266155
c14 -0.313597581583826
d14 0.832994971069888
a15 -0.117706384279221
b15 0.903149923347922
c15 -2.320613051785450
d15 2.505507862907740
a16 0.081491096318107
b16 -0.650590425311003
c16 1.719111854728430
d16 -0.995587056071390
a17 0.091741999005889
b17 -0.733622737082450
c17 1.943299096505840
d17 -1.197355573677230
a18 -0.148459092340618
b18 1.284066430229020
c18 -3.706230571961660
d18 4.075538783563990
a19 0.002094370357483
b19 -0.025748695244591
c19 0.092233291912327
27
d19 0.403690381817136
a20 -0.059918389089864
b20 0.532366139781370
c20 -1.582111213167470
d20 2.078034886893430
a21 0.337579186003643
b21 -3.164361308590620
c21 9.877743876782910
d21 -9.763815372717620
a22 -0.390398354935213
b22 3.824223084424570
c22 -12.485726180862100
d22 14.090552688768400
a23 0.124014233752400
b23 -1.268461543581590
c23 4.320133091562740
d23 -4.395892510893910
a24 -0.005658580079891
b24 0.054201157506717
c24 -0.176920092141704
d24 0.700767763970362
a25 0.098620086567359
b25 -1.040724842288780
c25 3.655320907143910
d25 -3.770180068531730
a26 -0.088821766188648
b26 0.983647167475564
c26 -3.632418328009360
28
d26 4.975107013649780
a27 0.056666978184561
b27 -0.631277895066812
c27 2.342804403393840
d27 -2.394334355084420
a28 -0.137846146550419
b28 1.586171726911120
c28 -6.083504160119400
d28 8.278989825371750
a29 -0.005282391978239
b29 0.035175798416903
c29 -0.034620038992117
d29 0.415440467905423
𝑆(𝑥) =
𝟎. 𝟖𝟔𝟕𝟓𝟎𝟔𝟔𝟒𝟐𝟏𝟒𝟗𝟓𝟓𝟔𝒙𝟑 − 𝟐. 𝟔𝟎𝟐𝟓𝟏𝟗𝟗𝟐𝟔𝟒𝟒𝟖𝟕𝟑𝟎𝒙𝟐 + 𝟐. 𝟑𝟔𝟔𝟖𝟒𝟒𝟖𝟔𝟎𝟎𝟐𝟕𝟑𝟎𝟎𝒙
−𝟎. 𝟎𝟎𝟔𝟖𝟑𝟏𝟓𝟕𝟓𝟕𝟐𝟖𝟎𝟓𝟓, 𝒙 = [𝟏. 𝟎𝟎, 𝟏. 𝟏𝟎]
𝟎. 𝟎𝟔𝟐𝟒𝟔𝟔𝟕𝟖𝟗𝟐𝟓𝟐𝟔𝟏𝟎𝒙𝟑 + 𝟎. 𝟎𝟓𝟒𝟏𝟏𝟏𝟓𝟖𝟖𝟏𝟏𝟏𝟐𝟔𝟎𝒙𝟐
− 𝟎. 𝟓𝟓𝟓𝟒𝟒𝟗𝟖𝟎𝟓𝟗𝟖𝟖𝟕𝟓𝟔𝒙
+𝟏. 𝟎𝟔𝟒𝟔𝟕𝟔𝟒𝟔𝟖𝟒𝟕𝟕𝟕𝟗𝟎, 𝒙 = [𝟏. 𝟏𝟎, 𝟏. 𝟐𝟎]
29
𝟎. 𝟎𝟓𝟎𝟎𝟓𝟓𝟔𝟑𝟖𝟏𝟎𝟐𝟒𝟕𝟒𝒙𝟑 − 𝟎. 𝟏𝟒𝟐𝟑𝟔𝟐𝟒𝟒𝟎𝟏𝟒𝟕𝟏𝟏𝟏𝒙𝟐
− 𝟎. 𝟎𝟎𝟔𝟖𝟐𝟖𝟎𝟗𝟕𝟓𝟒𝟗𝟐𝟗𝟓𝒙
+𝟎. 𝟕𝟎𝟖𝟑𝟑𝟕𝟎𝟒𝟖𝟑𝟎𝟒𝟏𝟓𝟔, 𝒙 = [𝟏. 𝟒𝟎, 𝟏. 𝟓𝟎]
𝟎. 𝟎𝟓𝟐𝟒𝟑𝟒𝟏𝟑𝟎𝟓𝟑𝟑𝟖𝟎𝟓𝒙𝟑 − 𝟎. 𝟏𝟓𝟑𝟎𝟔𝟓𝟔𝟓𝟔𝟎𝟖𝟖𝟎𝟗𝟐𝒙𝟐
+ 𝟎. 𝟎𝟎𝟗𝟐𝟐𝟔𝟕𝟐𝟔𝟑𝟔𝟐𝟏𝟕𝟖𝒙
𝟎. 𝟕𝟎𝟏𝟗𝟏𝟓𝟏𝟏𝟖𝟕𝟑𝟗𝟓𝟔𝟕, 𝒙 = [𝟏. 𝟓𝟎, 𝟏. 𝟔𝟎]
𝟎. 𝟏𝟎𝟐𝟕𝟒𝟑𝟑𝟓𝟕𝟐𝟔𝟖𝟒𝟏𝟗𝒙𝟑 − 𝟎. 𝟑𝟗𝟒𝟓𝟒𝟗𝟗𝟒𝟒𝟒𝟏𝟒𝟎𝟐𝟕𝒙𝟐
+ 𝟎. 𝟑𝟗𝟓𝟔𝟎𝟏𝟓𝟖𝟕𝟔𝟖𝟑𝟕𝟏𝟖𝒙
+𝟎. 𝟒𝟗𝟓𝟖𝟒𝟖𝟓𝟐𝟔𝟎𝟑𝟒𝟕𝟒𝟒, 𝒙 = [𝟏. 𝟔𝟎, 𝟏. 𝟕𝟎]
𝟎. 𝟒𝟒𝟏𝟓𝟕𝟕𝟒𝟐𝟔𝟎𝟑𝟒𝟒𝟐𝟐𝒙𝟑 − 𝟐. 𝟒𝟑𝟎𝟗𝟎𝟎𝟗𝟖𝟏𝟔𝟕𝟖𝟐𝟗𝟎𝒙𝟐
+ 𝟒. 𝟒𝟏𝟐𝟑𝟑𝟑𝟒𝟔𝟔𝟖𝟑𝟗𝟑𝟔𝟎𝒙
−𝟐. 𝟏𝟏𝟑𝟐𝟔𝟎𝟔𝟎𝟖𝟑𝟎𝟐𝟑𝟔𝟎, 𝒙 = [𝟏. 𝟖𝟎, 𝟏. 𝟗𝟎]
30
𝟎. 𝟎𝟔𝟎𝟑𝟔𝟖𝟔𝟐𝟏𝟎𝟖𝟔𝟏𝟑𝟔𝒙𝟑 − 𝟎. 𝟒𝟏𝟏𝟏𝟎𝟐𝟑𝟔𝟐𝟖𝟎𝟔𝟎𝟔𝟎𝒙𝟐
+ 𝟎. 𝟗𝟏𝟑𝟖𝟕𝟑𝟎𝟓𝟑𝟕𝟖𝟖𝟗𝟒𝟔𝒙
−𝟎. 𝟏𝟒𝟖𝟗𝟖𝟏𝟓𝟑𝟕𝟐𝟐𝟕𝟗𝟖𝟏, 𝒙 = [𝟐. 𝟑𝟎, 𝟐. 𝟒𝟎]
𝟎. 𝟎𝟖𝟏𝟒𝟗𝟏𝟎𝟗𝟔𝟑𝟏𝟖𝟏𝟎𝟕𝒙𝟑 − 𝟎. 𝟔𝟓𝟎𝟓𝟗𝟎𝟒𝟐𝟓𝟑𝟏𝟏𝟎𝟎𝟑𝒙𝟐
+ 𝟏. 𝟕𝟏𝟗𝟏𝟏𝟏𝟖𝟓𝟒𝟕𝟐𝟖𝟒𝟑𝟎𝒙
−𝟎. 𝟗𝟗𝟓𝟓𝟖𝟕𝟎𝟓𝟔𝟎𝟕𝟏𝟑𝟗𝟎, 𝒙 = [𝟐. 𝟔𝟎, 𝟐. 𝟕𝟎]
𝟎. 𝟎𝟗𝟏𝟕𝟒𝟏𝟗𝟗𝟗𝟎𝟎𝟓𝟖𝟖𝟗𝒙𝟑 − 𝟎. 𝟕𝟑𝟑𝟔𝟐𝟐𝟕𝟑𝟕𝟎𝟖𝟐𝟒𝟓𝟎𝒙𝟐
+ 𝟏. 𝟗𝟒𝟑𝟐𝟗𝟗𝟎𝟗𝟔𝟓𝟎𝟓𝟖𝟒𝟎𝒙
−𝟏. 𝟏𝟗𝟕𝟑𝟓𝟓𝟓𝟕𝟑𝟔𝟕𝟕𝟐𝟑𝟎, 𝒙 = [𝟐. 𝟕𝟎, 𝟐. 𝟖𝟎]
𝟎. 𝟎𝟎𝟐𝟎𝟗𝟒𝟑𝟕𝟎𝟑𝟓𝟕𝟒𝟖𝟑𝒙𝟑 − 𝟎. 𝟎𝟐𝟓𝟕𝟒𝟖𝟔𝟗𝟓𝟐𝟒𝟒𝟓𝟗𝟏𝒙𝟐
+ 𝟎. 𝟎𝟗𝟐𝟐𝟑𝟑𝟐𝟗𝟏𝟗𝟏𝟐𝟑𝟐𝟕𝒙
+𝟎. 𝟒𝟎𝟑𝟔𝟗𝟎𝟑𝟖𝟏𝟖𝟏𝟕𝟏𝟑𝟔, 𝒙 = [𝟐. 𝟗𝟎, 𝟑. 𝟎𝟎]
𝟎. 𝟑𝟑𝟕𝟓𝟕𝟗𝟏𝟖𝟔𝟎𝟎𝟑𝟔𝟒𝟑𝒙𝟑 − 𝟑. 𝟏𝟔𝟒𝟑𝟔𝟏𝟑𝟎𝟖𝟓𝟗𝟎𝟔𝟐𝟎𝒙𝟐
+ 𝟗. 𝟖𝟕𝟕𝟕𝟒𝟑𝟖𝟕𝟔𝟕𝟖𝟐𝟗𝟏𝟎𝒙
−𝟗. 𝟕𝟔𝟑𝟖𝟏𝟓𝟑𝟕𝟐𝟕𝟏𝟕𝟔𝟐𝟎, 𝒙 = [𝟑. 𝟏𝟎, 𝟑. 𝟐𝟎]
31
−𝟎. 𝟑𝟗𝟎𝟑𝟗𝟖𝟑𝟓𝟒𝟗𝟑𝟓𝟐𝟏𝟑𝒙𝟑 + 𝟑. 𝟖𝟐𝟒𝟐𝟐𝟑𝟎𝟖𝟒𝟒𝟐𝟒𝟓𝟕𝟎𝒙𝟐
− 𝟏𝟐. 𝟒𝟖𝟓𝟕𝟐𝟔𝟏𝟖𝟎𝟖𝟔𝟐𝟏𝟎𝟎𝒙
𝟏𝟒. 𝟎𝟗𝟎𝟓𝟓𝟐𝟔𝟖𝟖𝟕𝟔𝟖𝟒𝟎𝟎, 𝒙 = [𝟑. 𝟐𝟎, 𝟑. 𝟑𝟎]
𝟎. 𝟏𝟐𝟒𝟎𝟏𝟒𝟐𝟑𝟑𝟕𝟓𝟐𝟒𝟎𝟎𝒙𝟑 − 𝟏. 𝟐𝟔𝟖𝟒𝟔𝟏𝟓𝟒𝟑𝟓𝟖𝟏𝟓𝟗𝟎𝒙𝟐
+ 𝟒. 𝟑𝟐𝟎𝟏𝟑𝟑𝟎𝟗𝟏𝟓𝟔𝟐𝟕𝟒𝟎𝒙
−𝟒. 𝟑𝟗𝟓𝟖𝟗𝟐𝟓𝟏𝟎𝟖𝟗𝟑𝟗𝟏𝟎, 𝒙 = [𝟑. 𝟑𝟎, 𝟑. 𝟒𝟎]
𝟎. 𝟎𝟗𝟖𝟔𝟐𝟎𝟎𝟖𝟔𝟓𝟔𝟕𝟑𝟓𝟗𝒙𝟑 − 𝟏. 𝟎𝟒𝟎𝟕𝟐𝟒𝟖𝟒𝟐𝟐𝟖𝟖𝟕𝟖𝟎𝒙𝟐
+ 𝟑. 𝟔𝟓𝟓𝟑𝟐𝟎𝟗𝟎𝟕𝟏𝟒𝟑𝟗𝟏𝟎𝒙
−𝟑. 𝟕𝟕𝟎𝟏𝟖𝟎𝟎𝟔𝟖𝟓𝟑𝟏𝟕𝟑𝟎, 𝒙 = [𝟑. 𝟓𝟎, 𝟑. 𝟔𝟎]
𝟎. 𝟎𝟓𝟔𝟔𝟔𝟔𝟗𝟕𝟖𝟏𝟖𝟒𝟓𝟔𝟏𝒙𝟑 − 𝟎. 𝟔𝟑𝟏𝟐𝟕𝟕𝟖𝟗𝟓𝟎𝟔𝟔𝟖𝟏𝟐𝒙𝟐
+ 𝟐. 𝟑𝟒𝟐𝟖𝟎𝟒𝟒𝟎𝟑𝟑𝟗𝟑𝟖𝟒𝟎𝒙
−𝟐. 𝟑𝟗𝟒𝟑𝟑𝟒𝟑𝟓𝟓𝟎𝟖𝟒𝟒𝟐𝟎, 𝒙 = [𝟑. 𝟕𝟎, 𝟑. 𝟖𝟎]
32
5 Conclusiones:
● También se concluye que si bien es cierto que el hallar los coeficientes de cada una de las
ecuaciones puede llegar a ser complicado pero con la utilización de los software (Matlab)
que nos brinda el avance de la tecnología nos facilita enormemente ya que si no tendríamos
que usar un sistema de matrices y aún más el procedimiento se vuelve más tedioso cuando
existen varios intervalos.
● Por otro lado, que el método de spline cúbico nos facilita las cosas a la hora de graficar las
ecuaciones lo que al principio se podría pensar que las ecuaciones nos darían gráficas
complicadas pero por medio de este método hace que se aproximen gráficas sencillas fáciles
de construir como lo podemos ver en los resultados de nuestra gráfica, está presenta en la
mayor parte de los intervalos funciones decrecientes y cóncava hacia arriba en su mayoría.
● Podemos decir también que si bien es cierto que el método de spline cúbicos es una área
muy interesante y aún más en áreas donde se recolectan datos para hacer simulaciones
como es el caso de ingeniería civil donde se tiene que diseñar columnas puentes , el usar
esta herramientas facilita las cosas en cuanto al diseño de dichas construcciones de manera
eficiente.
33
6 Anexos
La cinemática de algunas piezas robóticas se basa en una trayectoria curvilínea, pero esta
curvatura no debe de ser muy brusca, es decir, curvas muy inclinadas o muy onduladas, sino más
bien curvas con pendiente suave. Justamente por estos últimos detalles de la curvatura adecuada
que debe de realizarse para poder realizar un óptimo movimiento de alguna pieza de un robot es el
método de interpolación por spline cúbicos. Una vez obtenido dicha gráfica de una función
interpolada por este método se le va a llevar a un programa que se encargue de llevar dicha
información de la gráfica a la pieza del robot al que se le va a emplear dicho movimiento. Pero antes,
dicha gráfica tiene que pasar por uno de los dos métodos de operación de variables e interpretación
de gráficas para que se pueda realizar dicho movimiento. Estos son:
*Método práctico (2D): Este método es fácil, pero puede resultar tedioso o impreciso. Se usa más
en figuras simétricas, en curvas de pendiente suave y que el trabajo sea realizado en superficies
2D y liso.
Figura 1.Curvas unidas obtenidas de algunas funciones interpoladas a través del procesamiento de imágenes (3D).Recuperado de
“Aplicación de la interpolación por spline cúbicos en definir trayectorias para movimientos de robots”.Documento académico
Se puede realizar las ecuaciones de la función interpolada a través de Matlab, C++ u otros que nos
ayuden con esta parte del cálculo.
En síntesis, el spline cúbico para el movimiento de los robots ayuda para poder trazar letras en
placas de automóviles o en algunos procesos industriales, etc.
34
Aplicación de los spline cúbicos naturales en las grúas de los buques
Para poder realizar el traslado de la grúa, los buques requieren del uso de medidas realizadas
de obtención de puntos de control in-situ para el desarrollo de la curva que la grúa va a recorrer,
esta curva debe ser continua y suave. En la obtención de los puntos de control que se desean
buscar tendrán un ruido asociado o formar parte de curvas muy pronunciadas, esto se puede
solucionar usando el método de interpolación por partes: el spline cúbico natural.
Figura 2. Presentación de una curva producida por una función interpolada mediante algunos puntos de control para el recorrido del
carguero.(https://repositorio.upct.es/xmlui/bitstream/handle/10317/8267/pfc-ege-spl.pdf?sequence=1&isAllowed=y
35
6.2 Anexo 2: Muestras de los resultados.
36
0.406994245907527*x+0.875152230229677)+((x>=2.1)&(x<2.2)).*(-
0.037132920516382*x.^3+0.259761075638484*x.^2-
0.624939017687054*x+1.027713570475550)+((x>=2.2)&(x<2.3)).*(-
0.030808925142413*x.^3+0.218022706170192*x.^2-
0.533114604858095*x+0.960375667733051)+((x>=2.3)&(x<2.4)).*(0.060368621086136*x.^3-
0.411102362806060*x.^2+0.913873053788946*x-0.148981537227981)+((x>=2.4)&(x<2.5)).*(-
0.010665559201712*x.^3+0.100343735266155*x.^2-
0.313597581583826*x+0.832994971069888)+((x>=2.5)&(x<2.6)).*(-
0.117706384279221*x.^3+0.903149923347922*x.^2-
2.320613051785450*x+2.505507862907740)+((x>=2.6)&(x<2.7)).*(0.081491096318107*x.^3-
0.650590425311003*x.^2+1.719111854728430*x-
0.995587056071390)+((x>=2.7)&(x<2.8)).*(0.091741999005889*x.^3-
0.733622737082450*x.^2+1.943299096505840*x-1.197355573677230)+((x>=2.8)&(x<2.9)).*(-
0.148459092340618*x.^3+1.284066430229020*x.^2-
3.706230571961660*x+4.075538783563990)+((x>=2.9)&(x<3)).*(0.002094370357483*x.^3-
0.025748695244591*x.^2+0.092233291912327*x+0.403690381817136)+((x>=3)&(x<3.1)).*(-
0.059918389089864*x.^3+0.532366139781370*x.^2-
1.582111213167470*x+2.078034886893430)+((x>=3.1)&(x<3.2)).*(0.337579186003643*x.^3-
3.164361308590620*x.^2+9.877743876782910*x-9.763815372717620)+((x>=3.2)&(x<3.3)).*(-
0.390398354935213*x.^3+3.824223084424570*x.^2-
12.485726180862100*x+14.090552688768400)+((x>=3.3)&(x<3.4)).*(0.124014233752400*x.^3-
1.268461543581590*x.^2+4.320133091562740*x-4.395892510893910)+((x>=3.4)&(x<3.5)).*(-
0.005658580079891*x.^3+0.054201157506717*x.^2-
0.176920092141704*x+0.700767763970362)+((x>=3.5)&(x<3.6)).*(0.098620086567359*x.^3-
1.040724842288780*x.^2+3.655320907143910*x-3.770180068531730)+((x>=3.6)&(x<3.7)).*(-
0.088821766188648*x.^3+0.983647167475564*x.^2+-
3.632418328009360*x+4.975107013649780)+((x>=3.7)&(x<3.8)).*(0.056666978184561*x.^3-
0.631277895066812*x.^2+2.342804403393840*x-2.394334355084420)+((x>=3.8)&(x<3.9)).*(-
0.137846146550419*x.^3+1.586171726911120*x.^2-
6.083504160119400*x+8.278989825371750)+((x>=3.9)&(x<=4)).*(-
0.005282391978239*x.^3+0.035175798416903*x.^2-
0.034620038992117*x+0.415440467905423);
plot(x,f,'r'),grid
end
37
7 Referencias:
● Molano, L. A. (2012). Ciencia y desarrollo. Interpolación polinomial, algunas técnicas y su
programación. Recuperado de: [Link].
● Vadillo, F. (2021). Una introducción a la interpolación polinomial. Recuperado de: [Link]
● González, C. Spline: Curvas y superficies. [Link]. ¿Cuál es la diferencia entre interpolación
y regresión? (s/f). Eluniversaledomex.mx. Recuperado el 17 de agosto de 2022, de [Link]
● Diferencia Fundamental Entre Regresión e Interpolación. (s/f). Scribd. Recuperado el 17 de
agosto de 2022, de [Link]
● Interpolación. (s/f). Ehu.es. Recuperado el 17 de agosto de 2022, de [Link]
● Del Rosario, E. [Link]
● Twigg, C. (2003). Catmull-Rom splines. Recuperado de: [Link]
● Gonzalo, E. (2017). Funciones spline. [Link].
● Pipenbrinck, N. (1998). Hermite Curve Interpolation. Recuperado de: [Link].
● Spline de Hermite Cúbico. (s.f.). Wikipedia. t. Recuperado de [Link]
● Interpolación Temario. (s.f.). Aragón, Universidad Nacional Autónoma de México. t.
Recuperado de [Link].
● Interpolación. (s.f.). Wikidot. Recuperado de [Link]
● Interpolación. (s.f.). Ecured.t. Recuperado de [Link]