stringtranslate.com

Ajuste de mínimos cuadrados

El ajuste por mínimos cuadrados es un modelo para la solución de un sistema de ecuaciones sobredeterminado basado en el principio de mínimos cuadrados de los residuos de observación . Se utiliza ampliamente en las disciplinas de topografía , geodesia y fotogrametría (el campo de la geomática , en conjunto).

Formulación

Hay tres formas de ajuste de mínimos cuadrados: paramétrico , condicional y combinado :

Claramente, los ajustes paramétricos y condicionales corresponden al caso combinado más general cuando f ( X , Y ) = h ( X ) - Y y f ( X , Y ) = g ( Y ) , respectivamente. Sin embargo, los casos especiales justifican soluciones más simples, como se detalla a continuación. A menudo en la literatura, Y puede denotarse L .

Solución

Las igualdades anteriores solo se cumplen para los parámetros estimados y las observaciones , por lo tanto . Por el contrario, las observaciones medidas y los parámetros aproximados producen un cierre erróneo distinto de cero : Se puede proceder a la expansión de las ecuaciones en serie de Taylor, lo que da como resultado los jacobianos o matrices de diseño : el primero, y el segundo, El modelo linealizado se lee entonces: donde son correcciones de parámetros estimados a los valores a priori , y son residuos de observación posteriores al ajuste .

En el ajuste paramétrico, la segunda matriz de diseño es una identidad, B = - I , y el vector de cierre incorrecto puede interpretarse como los residuos previos al ajuste, , por lo que el sistema se simplifica a: que está en forma de mínimos cuadrados ordinarios . En el ajuste condicional, la primera matriz de diseño es nula, A = 0 . Para los casos más generales, se introducen multiplicadores de Lagrange para relacionar las dos matrices jacobianas y transformar el problema de mínimos cuadrados restringido en uno sin restricciones (aunque más grande). En cualquier caso, su manipulación conduce a los vectores y , así como a los respectivos parámetros y observaciones, matrices de covarianza a posteriori .

Cálculo

Dadas las matrices y vectores anteriores, su solución se encuentra a través de métodos estándar de mínimos cuadrados; por ejemplo, formando la matriz normal y aplicando la descomposición de Cholesky , aplicando la factorización QR directamente a la matriz jacobiana, métodos iterativos para sistemas muy grandes, etc.

Ejemplos resueltos

Aplicaciones

Conceptos relacionados

Extensiones

Si se encuentra una deficiencia de rango , a menudo se puede rectificar mediante la inclusión de ecuaciones adicionales que imponen restricciones a los parámetros y/o las observaciones, lo que conduce a mínimos cuadrados restringidos .

Referencias

  1. ^ Kotz, Samuel; Read, Campbell B.; Balakrishnan, N.; Vidakovic, Brani; Johnson, Norman L. (15 de julio de 2004). "Modelo de Gauss-Helmert". Enciclopedia de ciencias estadísticas . Hoboken, NJ, EE. UU.: John Wiley & Sons, Inc. doi :10.1002/0471667196.ess0854.pub2. ISBN 978-0-471-66719-3.
  2. ^ Förstner, Wolfgang; Wrobel, Bernhard P. (2016). "Estimación". Visión artificial fotogramétrica . Geometría y computación. Vol. 11. Cham: Springer International Publishing. págs. 75–190. doi :10.1007/978-3-319-11550-4_4. ISBN 978-3-319-11549-8. ISSN  1866-6795.
  3. ^ Schaffrin, Burkhard; Snow, Kyle (2010). "Regularización de mínimos cuadrados totales del tipo Tykhonov y una antigua pista de carreras en Corinto". Álgebra lineal y sus aplicaciones . 432 (8). Elsevier BV: 2061–2076. doi : 10.1016/j.laa.2009.09.014 . ISSN  0024-3795.
  4. ^ Neitzel, Frank (17 de septiembre de 2010). "Generalización de mínimos cuadrados totales en el ejemplo de transformación de similitud 2D ponderada y no ponderada". Journal of Geodesy . 84 (12). Springer Science and Business Media LLC: 751–762. Bibcode :2010JGeod..84..751N. doi :10.1007/s00190-010-0408-0. ISSN  0949-7714. S2CID  123207786.

Bibliografía

Apuntes de clase e informes técnicos
Libros y capítulos