AJUSTE DE CURVAS
INTRODUCCIÓN
Supongamos que tenemos un conjunto de puntos que mostramos en la
siguiente gráfica
De los puntos mostrados nos podemos dar cuenta que parece
tener la forma de un polinomio de segundo grado de la forma:
|
|
(1) |
Esta ecuación (1) puede usarse para representar el conjunto
de valores obtenidos experimentalmente para la cual debemos
determinar los valores
de a
1,
a
2,
a
3, etc.
Para determinar estos valores utilizamos el siguiente procedimiento:
- Establecer el criterio para determinar la ecuación que
represente a los valores (obtenidos experimentalmente).
- Escribir la ecuación que expresa el error
o desviación entre el valor observado y los valores
dados por la ecuación.
- Habiendo obtenido la ecuación del error, minimizar dicho error.
EVALUACIÓN DEL ERROR
Si consideramos las parejas de datos, como se muestra en
la gráfica
donde:
d =
distancia = Yobservada - Y
obtenida por
la ecuación
Yobservada = Valor obtenido experimentalmente.
Y obtenida por la ecuación =
valor de la función evaluada en
cualquier valor X
Observando la gráfica, parece que esta
distancia se puede usar para representar el error,
pero habrá distancias positivas y
negativas, (como se puede observar la
distancia d1 es positiva
y la distancia d2 es negativa)
de modo que el error promedio para los puntos como
los mostrados será pequeño aunque los
errores individuales sean grandes.
Esta dificultad podría ser resuelta usando el
valor absoluto de las distancias, sin embargo al derivar
la función del valor absoluto se generan ciertos problemas.
La solución podría ser definir el error como
el cuadrado de la distancia, esto elimina la dificultad del
signo. Por esta razón el método se
llama: Método de Mínimos Cuadrados.
|
|
(2) |
en donde S es la suma de los cuadrados de las
diferencias entre el valor calculado y el valor observado y por lo
tanto es el valor que se debe minimizar
|
|
(3) |
Siendo el caso de que la curva supuesta es una ecuación de
segundo grado, se tiene la ecuación:
|
|
(4) |
Para minimizar la función anterior, derivando parcialmente
con respecto a
a
1,
a
2 y
a
3 e igualando a cero:
(Obsérvese que las variables son
a
1,
a
2 y
a
3, mientras
que Yi, X
i son constantes)
Las ecuaciones se pueden expresar de acuerdo como sigue:
Lo anterior lo podemos expresar en forma matricial:
|
|
(7) |
La fórmula general para un polinomio de
grado n en donde hay m parejas
de datos es:
|
|
(8) |
Como se puede observar el problema consiste en lo siguiente:
- Obtener la matriz de coeficientes.
- Resolver el sistema de ecuaciones resultantes.
Recordando que:
- Si n es el grado del polinomio,
hay n+1 valores de la matriz de
coeficientes y n+1 ecuaciones.
- El máximo exponenete de X en los términos de
la sumatoria de 2n puede ser que los datos no
representen un polinomio de 2o grado sino que representen uno
de 3o y 4o grados.
El ajuste de curvas es un procedimiento de tanteo y error, si
una curva no representa los datos, entonces se intenta con un
polinomio de grado superior.
|
|
EJEMPLO:
X | Y |
|
X | Y |
0.00 | 0.0000 |
|
0.60 | 0.6367 |
0.10 | 0.1002 |
|
0.70 | 0.7586 |
0.20 | 0.2013 |
|
0.80 | 0.8881 |
0.30 | 0.3045 |
|
0.90 | 1.0265 |
0.40 | 0.4108 |
|
1.00 | 1.1752 |
0.50 | 0.5211 |
|
|
De la tabla de datos, usando
Mínimos Cuadrados, determine
los polinomios de 2o, 3er y 4o grado; graficar para determinar la
curva más aproximada.
SOLUCIÓN
POLINOMIO DE SEGUNDO GRADO
Primero determinamos los coeficientes de la matriz y los elementos
constantes. Los elementos de la matriz son:
M =11
Los términos constantes son:
De acuerdo con esto, el sistema de ecuaciones a resolver es el
siguiente:
Resolviendo por Gauss se obtienen los siguientes resultados.
a1 = 0.006727
a2 = 0.895462
a3 = 0.265963
y el polinomio de segundo grado es:
POLINOMIO DE TERCER GRADO
Para el caso del polinomio de 3er grado se requiere:
y el sistema de ecuaciones a resolver es:
cuya solución es:
a1 = - 0.000112
|
a2 = 1.004150
|
a3 = - 0.019075
|
a4 = 0.190032
|
y el polinomio queda:
POLINOMIO DE CUARTO GRADO
Repitiendo el procedimiento anterior, se obtienen los siguientes
resultados
a1 = - 0.000112
|
a2 = 0.994595
|
a3 = 0.028713
|
a4 = 0.113563
|
a5 = 0.038237
|
quedando el polinomio como se muestra:
No hay comentarios:
Publicar un comentario