Articulo de referencia

Ajuste de curvas

Ajuste de una curva ruidosa mediante un modelo de pico asimétrico, con un proceso iterativo ( algoritmo de Gauss-Newton con factor de amortiguación variable α). El ajuste de cur...

Ajuste de una curva ruidosa mediante un modelo de pico asimétrico, con un proceso iterativo ( algoritmo de Gauss-Newton con factor de amortiguación variable α).

El ajuste de curvas [ 1 ] [ 2 ] es el proceso de construir una curva , o función matemática , que tenga el mejor ajuste a una serie de puntos de datos , [ 3 ] posiblemente sujeto a restricciones. [ 4 ] [ 5 ] El ajuste de curvas puede implicar interpolación , [ 6 ] [ 7 ] donde se requiere un ajuste exacto a los datos, o suavizado , [ 8 ] [ 9 ] en el que se construye una función "suave" que se ajusta aproximadamente a los datos. Un tema relacionado es el análisis de regresión , [ 10 ] [ 11 ] que se centra más en cuestiones de inferencia estadística como cuánta incertidumbre está presente en una curva que se ajusta a datos observados con errores aleatorios. Las curvas ajustadas se pueden utilizar como ayuda para la visualización de datos, [ 12 ] [ 13 ] para inferir valores de una función cuando no hay datos disponibles, [ 14 ] y para resumir las relaciones entre dos o más variables. [ 15 ] La extrapolación se refiere al uso de una curva ajustada más allá del rango de los datos observados, [ 16 ] y está sujeta a un grado de incertidumbre [ 17 ] ya que puede reflejar el método utilizado para construir la curva tanto como refleja los datos observados.

Para el análisis lineal-algebraico de datos, el "ajuste" generalmente significa intentar encontrar la curva que minimiza el desplazamiento vertical (eje y ) de un punto con respecto a la curva (por ejemplo, mínimos cuadrados ordinarios ). Sin embargo, para aplicaciones gráficas y de imágenes, el ajuste geométrico busca proporcionar el mejor ajuste visual; lo que generalmente significa intentar minimizar la distancia ortogonal a la curva (por ejemplo, mínimos cuadrados totales ), o incluir ambos ejes de desplazamiento de un punto con respecto a la curva. Los ajustes geométricos no son populares porque generalmente requieren cálculos no lineales y/o iterativos, aunque tienen la ventaja de un resultado más estético y geométricamente preciso. [ 18 ] [ 19 ] [ 20 ]

Ajuste algebraico de funciones a puntos de datos

Lo más común es ajustar una función de la forma y = f ( x ) .

Ajuste de líneas y funciones polinómicas a puntos de datos

Curvas polinómicas que se ajustan a una función sinusoidal
Curvas polinómicas que ajustan puntos generados con una función seno. La línea punteada negra representa los datos "reales", la línea roja es un polinomio de primer grado , la línea verde es de segundo grado , la línea naranja es de tercer grado y la línea azul es de cuarto grado .

La ecuación polinómica de primer grado

y=aincógnita+b{\displaystyle y=ax+b\;}

es una línea con pendiente a . Una línea conecta dos puntos cualesquiera, por lo que una ecuación polinómica de primer grado es un ajuste exacto que pasa por dos puntos cualesquiera con coordenadas x distintas.

Si se aumenta el orden de la ecuación a un polinomio de segundo grado, se obtienen los siguientes resultados:

y=aincógnita2+bincógnita+do.{\displaystyle y=ax^{2}+bx+c\;.}

Esto ajustará perfectamente una curva simple a tres puntos.

Si se aumenta el orden de la ecuación a un polinomio de tercer grado, se obtiene lo siguiente:

y=aincógnita3+bincógnita2+doincógnita+d.{\displaystyle y=ax^{3}+bx^{2}+cx+d\;.}

Esto encajará exactamente en cuatro puntos.

Una afirmación más general sería decir que se ajustará exactamente a cuatro restricciones . Cada restricción puede ser un punto, un ángulo o una curvatura (que es el recíproco del radio de un círculo osculador ). Las restricciones de ángulo y curvatura se suelen añadir a los extremos de una curva y, en tales casos, se denominan condiciones de extremo . Las condiciones de extremo idénticas se utilizan con frecuencia para garantizar una transición suave entre curvas polinómicas contenidas en una sola spline . También se podrían añadir restricciones de orden superior, como "el cambio en la tasa de curvatura". Esto, por ejemplo, sería útil en el diseño de cruces de autopistas para comprender la tasa de cambio de las fuerzas aplicadas a un automóvil (véase la sacudida ) mientras sigue el cruce, y para establecer límites de velocidad razonables en consecuencia.

La ecuación polinómica de primer grado también podría ajustarse con precisión a un punto y un ángulo, mientras que la ecuación polinómica de tercer grado podría ajustarse con precisión a dos puntos, una restricción angular y una restricción de curvatura. Son posibles muchas otras combinaciones de restricciones para estas ecuaciones y para ecuaciones polinómicas de orden superior.

Si existen más de n  +  1 restricciones ( donde n es el grado del polinomio), la curva polinómica aún puede ajustarse a dichas restricciones. Un ajuste exacto a todas las restricciones no es seguro (aunque podría ocurrir, por ejemplo, en el caso de un polinomio de primer grado que se ajuste con precisión a tres puntos colineales ). En general, sin embargo, se requiere algún método para evaluar cada aproximación. El método de mínimos cuadrados es una forma de comparar las desviaciones.

Se aducen varias razones para obtener un ajuste aproximado cuando es posible simplemente aumentar el grado de la ecuación polinómica y obtener un ajuste exacto.

  • Aunque exista una coincidencia exacta, no significa necesariamente que pueda descubrirse fácilmente. Dependiendo del algoritmo utilizado, puede darse el caso de que no se pueda calcular la coincidencia exacta o que el cálculo requiera demasiado tiempo de procesamiento. En tal caso, podría ser necesaria una solución aproximada.
  • El efecto de promediar los puntos de datos cuestionables en una muestra, en lugar de distorsionar la curva para que se ajuste exactamente a ellos, puede resultar deseable.
  • Fenómeno de Runge : los polinomios de alto grado pueden ser altamente oscilatorios. Si una curva pasa por dos puntos A y B , se esperaría que también pasara cerca del punto medio entre A y B. Esto puede no ocurrir con las curvas de polinomios de alto grado; incluso pueden tener valores muy grandes, tanto positivos como negativos . Con los polinomios de bajo grado, es más probable que la curva pase cerca del punto medio (incluso se garantiza que pase exactamente por el punto medio en un polinomio de primer grado).
  • Los polinomios de bajo grado tienden a ser suaves, mientras que las curvas polinómicas de alto grado tienden a ser irregulares. Para definir esto con mayor precisión, el número máximo de puntos de inflexión posibles en una curva polinómica es n-2 , donde n es el grado de la ecuación polinómica. Un punto de inflexión es un punto en la curva donde cambia de radio positivo a negativo. También podemos decir que es donde pasa de "retener agua" a "desprender agua". Cabe destacar que solo es "posible" que los polinomios de alto grado sean irregulares; también podrían ser suaves, pero no hay garantía de ello, a diferencia de las curvas polinómicas de bajo grado. Un polinomio de grado quince podría tener, como máximo, trece puntos de inflexión, pero también podría tener once, nueve o cualquier número impar hasta uno. (Los polinomios de grado par podrían tener cualquier número par de puntos de inflexión desde n  -  2 hasta cero).

Un grado de la curva polinómica superior al necesario para un ajuste exacto resulta indeseable por todas las razones mencionadas anteriormente para polinomios de alto grado, y además conduce a un número infinito de soluciones. Por ejemplo, un polinomio de primer grado (una línea) restringido por un solo punto, en lugar de los dos habituales, daría lugar a un número infinito de soluciones. Esto plantea el problema de cómo comparar y elegir una sola solución, lo cual puede ser complicado tanto para el software como para los humanos. Por ello, lo más recomendable suele ser elegir el grado más bajo posible para lograr un ajuste exacto en todas las restricciones, e incluso un grado menor si se acepta un ajuste aproximado.

Relación entre el rendimiento del trigo y la salinidad del suelo [ 21 ]

Ajustar otras funciones a los puntos de datos

En ciertos casos, también se pueden utilizar otros tipos de curvas, como las funciones trigonométricas (por ejemplo, el seno y el coseno).

En espectroscopia, los datos pueden ajustarse mediante funciones gaussianas , lorentzianas , de Voigt y funciones relacionadas.

En biología, ecología, demografía, epidemiología y muchas otras disciplinas, el crecimiento de una población , la propagación de enfermedades infecciosas, etc., se pueden ajustar utilizando la función logística .

En agricultura, la función sigmoide logística invertida (curva en S) se utiliza para describir la relación entre el rendimiento de los cultivos y los factores de crecimiento. La figura azul se obtuvo mediante una regresión sigmoide de datos medidos en tierras agrícolas. Se observa que, inicialmente, es decir, con baja salinidad del suelo, el rendimiento de los cultivos disminuye lentamente al aumentar la salinidad, mientras que posteriormente la disminución se acelera.

Ajuste geométrico de curvas planas a puntos de datos

Si una función de la formay=F(incógnita){\displaystyle y=f(x)}Aunque no se puede postular, aún se puede intentar ajustar una curva plana .

En ciertos casos, también se pueden utilizar otros tipos de curvas, como secciones cónicas (arcos circulares, elípticos, parabólicos e hiperbólicos) o funciones trigonométricas (como el seno y el coseno). Por ejemplo, las trayectorias de los objetos bajo la influencia de la gravedad siguen una trayectoria parabólica cuando se ignora la resistencia del aire. Por lo tanto, sería lógico que los puntos de datos de la trayectoria se ajustaran a una curva parabólica. Las mareas siguen patrones sinusoidales, por lo que los puntos de datos de las mareas deberían ajustarse a una onda sinusoidal, o a la suma de dos ondas sinusoidales de diferentes periodos, si se consideran los efectos de la Luna y el Sol.

Para una curva paramétrica , es efectivo ajustar cada una de sus coordenadas como una función separada de la longitud del arco ; suponiendo que los puntos de datos se pueden ordenar, se puede utilizar la distancia de la cuerda . [ 22 ]

Ajuste de un círculo mediante ajuste geométrico

Ajuste de círculo con el método de Coope, los puntos que describen un arco de círculo, centro (1  ; 1), radio 4
diferentes modelos de ajuste de elipses
Ajuste de elipses minimizando la distancia algebraica (método de Fitzgibbon)

Coope [ 23 ] aborda el problema de encontrar el mejor ajuste visual de un círculo a un conjunto de puntos de datos 2D. El método transforma elegantemente el problema, normalmente no lineal, en un problema lineal que puede resolverse sin utilizar métodos numéricos iterativos, y por lo tanto es mucho más rápido que las técnicas anteriores.

Ajuste de una elipse mediante ajuste aritmético

La ecuación de una sección cónica (incluidas las elipses) tiene la forma [ 24 ].

Aincógnita2+Bincógnitay+doy2+Dincógnita+miy+F=0{\displaystyle Ax^{2}+Bxy+Cy^{2}+Dx+Ey+F=0}

para puntos con coordenadas cartesianas (x,y) en el borde. Para un conjunto de datos de puntos dispersos (x i ,y i ) y una elipse fija, los segundos miembros de esta ecuación ya no son cero. Un enfoque para ajustar los puntos a una elipse es encontrar los seis parámetros A a F tales que la suma de los (cuadrados de) los segundos miembros residuales se minimice. Esto conduce a un problema de ajuste lineal por mínimos cuadrados que se resuelve básicamente encontrando los autovectores de una matriz asociada de 6 x 6 que contiene sumas de potencias mixtas de (x i ,y i ) . [ 25 ] [ 26 ] [ 27 ] [ 28 ]

Ajuste de una elipse mediante ajuste geométrico

La técnica anterior se extiende a elipses generales [ 29 ] mediante la adición de un paso no lineal, lo que da como resultado un método rápido que, además, encuentra elipses visualmente agradables de orientación y desplazamiento arbitrarios. El ajuste geométrico minimiza la suma de las distancias ortogonales al cuadrado de los puntos de datos al borde de la elipse. [ 30 ] [ 31 ] Un método similar que aproxima estas distancias ortogonales se basa en la métrica de Sampson. [ 32 ] [ 33 ]

Superficies de ajuste

Cabe señalar que, si bien esta discusión se centró en curvas bidimensionales, gran parte de esta lógica también se aplica a superficies tridimensionales, cada una de las cuales se define mediante una red de curvas en dos direcciones paramétricas, generalmente denominadas u y v . Una superficie puede estar compuesta por uno o más parches de superficie en cada dirección.

Software

Muchos paquetes estadísticos como R y software numérico como gnuplot , GNU Scientific Library , Igor Pro , MLAB , Maple , MATLAB , TK Solver 6.0, Scilab , Mathematica , GNU Octave y SciPy incluyen comandos para realizar ajustes de curvas en diversos escenarios. También existen programas específicamente diseñados para el ajuste de curvas; estos se pueden encontrar en las listas de programas de análisis estadístico y numérico, así como en la categoría: Software de regresión y ajuste de curvas .

Véase también

Referencias

  1. Sandra Lach Arlinghaus, PHB Manual práctico de ajuste de curvas. CRC Press, 1994.
  2. William M. Kolb. Ajuste de curvas para calculadoras programables . Syntec, Incorporated, 1984.
  3. SS Halli, KV Rao. 1992. Técnicas avanzadas de análisis de población. ISBN 0306439972Página 165 ( cf. ... las funciones se cumplen si tenemos un ajuste de bueno a moderado para los datos observados).
  4. La señal y el ruido: por qué fallan tantas predicciones, pero algunas no. Por Nate Silver
  5. Preparación de datos para la minería de datos : Texto. Por Dorian Pyle.
  6. Métodos numéricos en ingeniería con MATLAB®. Por Jaan Kiusalaas. Página 24.
  7. Métodos numéricos en ingeniería con Python 3. Por Jaan Kiusalaas. Página 21.
  8. Métodos numéricos de ajuste de curvas . Por PG Guest, Philip George Guest. Página 349.
  9. Véase también: Suavizador
  10. Ajuste de modelos a datos biológicos mediante regresión lineal y no lineal . Por Harvey Motulsky y Arthur Christopoulos.
  11. Análisis de regresión por Rudolf J. Freund, William J. Wilson, Ping Sa. Página 269.
  12. Informática Visual. Editado por Halimah Badioze Zaman, Peter Robinson, Maria Petrou, Patrick Olivier, Heiko Schröder. Página 689.
  13. Métodos numéricos para modelos de ingeniería no lineales . Por John R. Hauser. Página 227.
  14. Métodos de física experimental: espectroscopia, volumen 13, parte 1. Por Claire Marton. Página 150.
  15. Enciclopedia del diseño de investigación, volumen 1. Editado por Neil J. Salkind. Página 266.
  16. Técnicas de análisis y planificación comunitaria . Por Richard E. Klosterman. Página 1.
  17. Introducción al riesgo y la incertidumbre en la evaluación de inversiones ambientales. DIANE Publishing. Pág. 69
  18. Ahn, Sung-Joon (diciembre de 2008), "Ajuste geométrico de curvas y superficies paramétricas" (PDF) , Journal of Information Processing Systems , 4 (4): 153–158 , doi : 10.3745/JIPS.2008.4.4.153 , archivado del original (PDF) el 13 de marzo de 2014.
  19. Chernov, N.; Ma, H. (2011), "Ajuste por mínimos cuadrados de curvas y superficies cuadráticas", en Yoshida, Sota R. (ed.), Visión por computadora , Nova Science Publishers, pp. 285–302 , ISBN  9781612093994
  20. Liu, Yang; Wang, Wenping (2008), "Una revisión del ajuste de distancia ortogonal por mínimos cuadrados de curvas y superficies paramétricas", en Chen, F.; Juttler, B. (eds.), Avances en modelado y procesamiento geométrico , Lecture Notes in Computer Science, vol. 4975, pp. 384–397 , CiteSeerX 10.1.1.306.6085 , doi : 10.1007/978-3-540-79246-8_29 , ISBN    978-3-540-79245-1
  21. Calculadora para regresión sigmoide
  22. pág. 51 en Ahlberg y Nilson (1967) La teoría de las splines y sus aplicaciones , Academic Press, 1967
  23. Coope, ID (1993). "Ajuste de círculos mediante mínimos cuadrados lineales y no lineales". Journal of Optimization Theory and Applications . 76 (2): 381– 388. doi : 10.1007/BF00939613 . hdl : 10092/11104 . S2CID 59583785 . 
  24. Rosin, Paul L. (1993). "Una nota sobre el ajuste de elipses por mínimos cuadrados". Pattern Recogn. Lett . 14 (10): 799– 808. Bibcode : 1993PaReL..14..799R . doi : 10.1016/0167-8655(93)90062-I .
  25. Fitzgibbon, Andrew; Pilu, Maurizio; Fisher, Robert B (1999). "Ajuste directo de elipses por mínimos cuadrados". IEEE Transactions on Pattern Analysis and Machine Intelligence . 21 (5): 476– 480. Bibcode : 1999ITPAM..21..476F . doi : 10.1109/34.765658 .
  26. Fitzgibbon, Andrew; Pilu, Maurizio; Fisher, Robert B (1996). "Ajuste directo de elipses por mínimos cuadrados". Actas de la 13.ª Conferencia Internacional sobre Reconocimiento de Patrones . Vol. 1. págs. 253–257 . doi : 10.1109/ICPR.1996.546029 . ISBN   0-8186-7282-X.
  27. Halir, Radim; Flusser, Jan (1998). "Ajuste de elipses mediante mínimos cuadrados directos numéricamente estables" (PDF) . Actas de la 6.ª Conferencia Internacional de Visión Gráfica Computacional de la UE Central .
  28. Al-Sharadqah, A.; Chernov, N. (2012). "Un ajuste elíptico doblemente óptimo". Comput. Stat. Data Anal . 56 (9): 2771– 2781. doi : 10.1016/j.csda.2012.02.028 .
  29. Paul Sheer, Un asistente de software para fotometrología estéreo manual , tesis de maestría, 1997
  30. Gander, Walter; Golub, Gene H.; Strebel, Rolf (1996). "Ajuste por mínimos cuadrados de círculos y elipses" . Bull. Belg. Math. Soc .: 63–84 .
  31. Ahn, Sung Joon; Rau, Wolfgang; Warnecke, Hans-Jurgen (2001). "Ajuste de distancias ortogonales por mínimos cuadrados de círculo, esfera, elipse, hipérbola y parábola". Patt. Recog . 34 (12): 2283– 2303. Bibcode : 2001PatRe..34.2283A . doi : 10.1016/S0031-3203(00)00152-7 .
  32. ^ Szpak, Zygmunt L.; Chojnacki, Wojciech; van den Hengel, Antón (2012). "Ajuste de Elipse garantizado con la distancia Sampson". Lectura. No. comp. Ciencia . Apuntes de conferencias sobre informática. vol. 7576. págs. 87– 100. doi : 10.1007/978-3-642-33715-4_7 . ISBN   978-3-642-33714-7.
  33. Szpak, Zygmunt L.; Chojnacki, Wojciech; van den Hengel, Anton (2015). "Ajuste garantizado de elipse con una región de confianza y una medida de incertidumbre para el centro, los ejes y la orientación". J. Math. Imag. Vis . 52 (2): 173– 199. Bibcode : 2015JMIV...52..173S . doi : 10.1007/s10851-014-0536-x .

Lecturas adicionales

  • N. Chernov (2010), Regresión circular y lineal: Ajuste de círculos y líneas por mínimos cuadrados , Chapman & Hall/CRC, Monografías sobre estadística y probabilidad aplicada, Volumen 117 (256 págs.).