Problema inverso

Proceso de cálculo de los factores causales que produjeron un conjunto de observaciones.

Un problema inverso en ciencia es el proceso de calcular a partir de un conjunto de observaciones los factores causales que las produjeron: por ejemplo, calcular una imagen en tomografía computarizada de rayos X , la reconstrucción de la fuente en acústica o calcular la densidad de la Tierra a partir de mediciones de su campo gravitatorio . Se llama problema inverso porque comienza con los efectos y luego calcula las causas. Es el inverso de un problema directo, que comienza con las causas y luego calcula los efectos.

Los problemas inversos son algunos de los problemas matemáticos más importantes en ciencia y matemáticas porque nos informan sobre parámetros que no podemos observar directamente. Tienen una amplia aplicación en identificación de sistemas , óptica , radar , acústica , teoría de la comunicación , procesamiento de señales , imágenes médicas , visión artificial , ^{[1] [2]} geofísica , oceanografía , astronomía , teledetección , procesamiento del lenguaje natural , aprendizaje automático , ^[3] pruebas no destructivas , análisis de estabilidad de taludes ^[4] y muchos otros campos. ^{[ cita requerida ]}

Historia

Comenzar por los efectos para descubrir las causas ha preocupado a los físicos durante siglos. Un ejemplo histórico son los cálculos de Adams y Le Verrier que llevaron al descubrimiento de Neptuno a partir de la trayectoria perturbada de Urano . Sin embargo, el estudio formal de los problemas inversos no se inició hasta el siglo XX.

Uno de los primeros ejemplos de una solución a un problema inverso fue descubierto por Hermann Weyl y publicado en 1911, describiendo el comportamiento asintótico de los valores propios del operador de Laplace-Beltrami . ^[5] Hoy conocida como la ley de Weyl , tal vez se comprenda más fácilmente como una respuesta a la pregunta de si es posible escuchar la forma de un tambor . Weyl conjeturó que las frecuencias propias de un tambor estarían relacionadas con el área y el perímetro del tambor mediante una ecuación particular, un resultado mejorado por matemáticos posteriores.

El campo de los problemas inversos fue abordado posteriormente por el físico soviético - armenio Viktor Ambartsumian . ^{[6] [7]}

Siendo aún estudiante, Ambartsumian estudió a fondo la teoría de la estructura atómica, la formación de los niveles de energía y la ecuación de Schrödinger y sus propiedades, y cuando dominó la teoría de los valores propios de las ecuaciones diferenciales , señaló la aparente analogía entre los niveles de energía discretos y los valores propios de las ecuaciones diferenciales. Entonces preguntó: dada una familia de valores propios, ¿es posible encontrar la forma de las ecuaciones cuyos valores propios son? Esencialmente, Ambartsumian estaba examinando el problema inverso de Sturm-Liouville , que trataba de determinar las ecuaciones de una cuerda vibrante. Este artículo fue publicado en 1929 en la revista de física alemana Zeitschrift für Physik y permaneció en la oscuridad durante bastante tiempo. Al describir esta situación después de muchas décadas, Ambartsumian dijo: "Si un astrónomo publica un artículo con contenido matemático en una revista de física, entonces lo más probable es que caiga en el olvido".

Sin embargo, hacia el final de la Segunda Guerra Mundial, este artículo, escrito por Ambartsumian, de 20 años, fue encontrado por matemáticos suecos y formó el punto de partida para toda un área de investigación sobre problemas inversos, convirtiéndose en la base de toda una disciplina.

Después, Gelfand y Levitan, en la Unión Soviética, dedicaron importantes esfuerzos a la "solución directa" del problema de dispersión inversa . ^[8] Propusieron un método analítico constructivo para determinar la solución. Cuando las computadoras estuvieron disponibles, algunos autores investigaron la posibilidad de aplicar su enfoque a problemas similares, como el problema inverso en la ecuación de onda 1D. Pero rápidamente se demostró que la inversión es un proceso inestable: el ruido y los errores pueden amplificarse enormemente, haciendo que una solución directa sea difícilmente practicable. Luego, alrededor de los años setenta, aparecieron los métodos de mínimos cuadrados y probabilísticos, que resultaron muy útiles para la determinación de parámetros involucrados en varios sistemas físicos. Este enfoque tuvo mucho éxito. Hoy en día, los problemas inversos también se investigan en campos ajenos a la física, como la química, la economía y la informática. Con el tiempo, a medida que los modelos numéricos se vuelvan predominantes en muchas partes de la sociedad, podemos esperar un problema inverso asociado con cada uno de estos modelos numéricos.

Comprensión conceptual

Desde Newton, los científicos han intentado extensamente modelar el mundo. En particular, cuando se dispone de un modelo matemático (por ejemplo, la ley de la gravedad de Newton o la ecuación de Coulomb para la electrostática), podemos prever, dados algunos parámetros que describen un sistema físico (como una distribución de masa o una distribución de cargas eléctricas), el comportamiento del sistema. Este enfoque se conoce como modelado matemático y los parámetros físicos mencionados anteriormente se denominan parámetros del modelo o simplemente el modelo . Para ser precisos, introducimos la noción de estado del sistema físico : es la solución de la ecuación del modelo matemático. En la teoría del control óptimo , estas ecuaciones se denominan ecuaciones de estado . En muchas situaciones, no estamos realmente interesados ​​en conocer el estado físico, sino solo sus efectos sobre algunos objetos (por ejemplo, los efectos que tiene el campo gravitatorio sobre un planeta específico). Por lo tanto, tenemos que introducir otro operador, llamado operador de observación , que convierte el estado del sistema físico (aquí el campo gravitatorio predicho) en lo que queremos observar (aquí los movimientos del planeta considerado). Ahora podemos introducir el llamado problema hacia adelante , que consta de dos pasos:

• determinación del estado del sistema a partir de los parámetros físicos que lo describen
• aplicación del operador de observación al estado estimado del sistema para predecir el comportamiento de lo que queremos observar.

Esto lleva a introducir otro operador ( F significa "adelante") que asigna los parámetros del modelo a los datos que el modelo predice que son el resultado de este procedimiento de dos pasos. El operador se denomina operador de avance o mapa de avance . En este enfoque, básicamente, intentamos predecir los efectos conociendo las causas. ${\displaystyle F}$ ${\displaystyle p}$ ${\displaystyle F(p)}$ ${\displaystyle p}$ ${\displaystyle F}$

La siguiente tabla muestra, considerando la Tierra como sistema físico y para diferentes fenómenos físicos, los parámetros del modelo que describen el sistema, la cantidad física que describe el estado del sistema físico y las observaciones que se hacen comúnmente sobre el estado del sistema.

En el enfoque del problema inverso, en términos generales, tratamos de conocer las causas dados los efectos.

Enunciado general del problema inverso

El problema inverso es el "inverso" del problema directo: en lugar de determinar los datos producidos por parámetros de modelo particulares, queremos determinar los parámetros de modelo que producen los datos que son la observación que hemos registrado (el subíndice obs significa observado). Nuestro objetivo, en otras palabras, es determinar los parámetros de modelo de manera que (al menos aproximadamente) donde sea el mapa directo. Denotamos por el número (posiblemente infinito) de parámetros de modelo y por el número de datos registrados. Presentamos algunos conceptos útiles y las notaciones asociadas que se utilizarán a continuación: ${\displaystyle d_{\text{obs}}}$ ${\displaystyle p}$ $${\displaystyle d_{\text{obs}}=F(p)}$$ ${\displaystyle F}$ ${\displaystyle M}$ ${\displaystyle N}$

• El espacio de modelos denotado por : el espacio vectorial abarcado por los parámetros del modelo; tiene dimensiones; ${\displaystyle P}$ ${\displaystyle M}$
• El espacio de datos denotado por : si organizamos las muestras medidas en un vector con componentes (si nuestras mediciones consisten en funciones, es un espacio vectorial con dimensiones infinitas); ${\displaystyle D}$ ${\displaystyle D=\mathbb {R} ^{N}}$ ${\displaystyle N}$ ${\displaystyle D}$
• ${\displaystyle F(p)}$:la respuesta del modelo ; consiste en los datos predichos por el modelo ; ${\displaystyle p}$ ${\displaystyle p}$
• ${\displaystyle F(P)}$:la imagen del mapa hacia adelante, es un subconjunto de (pero no un subespacio a menos que sea lineal) hecho de respuestas de todos los modelos; ${\displaystyle P}$ ${\displaystyle D}$ ${\displaystyle F}$
• ${\displaystyle d_{\text{obs}}-F(p)}$: los desajustes de datos (o residuos) asociados con el modelo : se pueden organizar como un vector, un elemento de . ${\displaystyle p}$ ${\displaystyle D}$

El concepto de residuos es muy importante: en el marco de la búsqueda de un modelo que coincida con los datos, su análisis revela si el modelo considerado puede considerarse realista o no . Las discrepancias sistemáticas no realistas entre los datos y las respuestas del modelo también revelan que el mapa de avance es inadecuado y pueden brindar información sobre un mapa de avance mejorado.

Cuando el operador es lineal, el problema inverso es lineal. De lo contrario, que es lo más frecuente, el problema inverso es no lineal. Además, los modelos no siempre se pueden describir mediante un número finito de parámetros. Es el caso cuando buscamos parámetros distribuidos (una distribución de velocidades de onda, por ejemplo): en tales casos, el objetivo del problema inverso es recuperar una o varias funciones. Tales problemas inversos son problemas inversos con dimensión infinita. ${\displaystyle F}$

Problemas de inversa lineal

En el caso de un mapa directo lineal y cuando tratamos con un número finito de parámetros del modelo, el mapa directo se puede escribir como un sistema lineal donde es la matriz que caracteriza el mapa directo. $${\displaystyle d=Fp}$$ ${\displaystyle F}$

Un ejemplo elemental: el campo gravitacional de la Tierra

Solo unos pocos sistemas físicos son realmente lineales con respecto a los parámetros del modelo. Uno de esos sistemas de la geofísica es el del campo gravitatorio de la Tierra . El campo gravitatorio de la Tierra está determinado por la distribución de la densidad de la Tierra en el subsuelo. Debido a que la litología de la Tierra cambia de manera bastante significativa, podemos observar diferencias minúsculas en el campo gravitatorio de la Tierra en la superficie de la Tierra. A partir de nuestra comprensión de la gravedad (Ley de gravitación de Newton), sabemos que la expresión matemática para la gravedad es: aquí es una medida de la aceleración gravitatoria local, es la constante gravitatoria universal , es la masa local (que está relacionada con la densidad) de la roca en el subsuelo y es la distancia desde la masa hasta el punto de observación. $${\displaystyle d={\frac {Gp}{r^{2}}};}$$ ${\displaystyle d}$ ${\displaystyle G}$ ${\displaystyle p}$ ${\displaystyle r}$

Al discretizar la expresión anterior, podemos relacionar las observaciones de datos discretos en la superficie de la Tierra con los parámetros discretos del modelo (densidad) en el subsuelo sobre los que deseamos saber más. Por ejemplo, considere el caso en el que tenemos mediciones realizadas en 5 ubicaciones en la superficie de la Tierra. En este caso, nuestro vector de datos es un vector columna de dimensión (5×1): su componente -ésimo está asociado con la ubicación de observación -ésima. También sabemos que solo tenemos cinco masas desconocidas en el subsuelo (poco realista, pero se usa para demostrar el concepto) con ubicación conocida: denotamos por la distancia entre la ubicación de observación -ésima y la masa -ésima. Por lo tanto, podemos construir el sistema lineal que relaciona las cinco masas desconocidas con los cinco puntos de datos de la siguiente manera: ${\displaystyle d}$ ${\displaystyle i}$ ${\displaystyle i}$ ${\displaystyle p_{j}}$ ${\displaystyle r_{ij}}$ ${\displaystyle i}$ ${\displaystyle j}$ $${\displaystyle d=Fp,}$$ $${\displaystyle d={\begin{bmatrix}d_{1}\\d_{2}\\d_{3}\\d_{4}\\d_{5}\end{bmatrix}},\quad p={\begin{bmatrix}p_{1}\\p_{2}\\p_{3}\\p_{4}\\p_{5}\end{bmatrix}},}$$ $${\displaystyle F={\begin{bmatrix}{\frac {G}{r_{11}^{2}}}&{\frac {G}{r_{12}^{2}}}&{\frac {G}{r_{13}^{2}}}&{\frac {G}{r_{14}^{2}}}&{\frac {G}{r_{15}^{2}}}\\{\frac {G}{r_{21}^{2}}}&{\frac {G}{r_{22}^{2}}}&{\frac {G}{r_{23}^{2}}}&{\frac {G}{r_{24}^{2}}}&{\frac {G}{r_{25}^{2}}}\\{\frac {G}{r_{31}^{2}}}&{\frac {G}{r_{32}^{2}}}&{\frac {G}{r_{33}^{2}}}&{\frac {G}{r_{34}^{2}}}&{\frac {G}{r_{35}^{2}}}\\{\frac {G}{r_{41}^{2}}}&{\frac {G}{r_{42}^{2}}}&{\frac {G}{r_{43}^{2}}}&{\frac {G}{r_{44}^{2}}}&{\frac {G}{r_{45}^{2}}}\\{\frac {G}{r_{51}^{2}}}&{\frac {G}{r_{52}^{2}}}&{\frac {G}{r_{53}^{2}}}&{\frac {G}{r_{54}^{2}}}&{\frac {G}{r_{55}^{2}}}\end{bmatrix}}}$$

Para calcular los parámetros del modelo que se ajustan a nuestros datos, podríamos invertir la matriz para convertir directamente las mediciones en los parámetros de nuestro modelo. Por ejemplo: un sistema con cinco ecuaciones y cinco incógnitas es una situación muy específica: nuestro ejemplo fue diseñado para terminar con esta especificidad. En general, las cantidades de datos y de incógnitas son diferentes, por lo que la matriz no es cuadrada. ${\displaystyle F}$ $${\displaystyle p=F^{-1}d_{\text{obs}}}$$ ${\displaystyle F}$

Sin embargo, incluso una matriz cuadrada no puede tener inversa: la matriz puede ser deficiente en rango (es decir, tener cero valores propios) y la solución del sistema no es única. Entonces, la solución del problema inverso será indeterminada. Esta es una primera dificultad. Los sistemas sobredeterminados (más ecuaciones que incógnitas) tienen otros problemas. Además, el ruido puede corromper nuestras observaciones, haciendo que posiblemente estén fuera del espacio de posibles respuestas a los parámetros del modelo, de modo que la solución del sistema puede no existir. Esta es otra dificultad. ${\displaystyle F}$ ${\displaystyle p=F^{-1}d_{\text{obs}}}$ ${\displaystyle d}$ ${\displaystyle F(P)}$ ${\displaystyle p=F^{-1}d_{\text{obs}}}$

Herramientas para superar la primera dificultad

La primera dificultad refleja un problema crucial: nuestras observaciones no contienen suficiente información y se requieren datos adicionales. Los datos adicionales pueden provenir de información física previa sobre los valores de los parámetros, sobre su distribución espacial o, de manera más general, sobre su dependencia mutua. También pueden provenir de otros experimentos: por ejemplo, podemos pensar en integrar datos registrados por gravímetros y sismógrafos para una mejor estimación de las densidades. La integración de esta información adicional es básicamente un problema de estadística . Esta disciplina es la que puede responder a la pregunta: ¿Cómo mezclar cantidades de diferente naturaleza? Seremos más precisos en la sección "Enfoque bayesiano" a continuación.

En lo que respecta a los parámetros distribuidos, la información previa sobre su distribución espacial a menudo consiste en información sobre algunas derivadas de estos parámetros distribuidos. Además, es una práctica común, aunque un tanto artificial, buscar el modelo "más simple" que coincida razonablemente con los datos. Esto generalmente se logra penalizando la norma del gradiente (o la variación total ) de los parámetros (este enfoque también se conoce como maximización de la entropía). También se puede simplificar el modelo mediante una parametrización que introduzca grados de libertad solo cuando sea necesario. ${\displaystyle L^{1}}$

También se puede integrar información adicional mediante restricciones de desigualdad sobre los parámetros del modelo o algunas funciones de ellos. Estas restricciones son importantes para evitar valores irreales para los parámetros (valores negativos, por ejemplo). En este caso, el espacio abarcado por los parámetros del modelo ya no será un espacio vectorial sino un subconjunto de modelos admisibles denotados por en la secuela. ${\displaystyle P_{\text{adm}}}$

Herramientas para superar la segunda dificultad

Como se mencionó anteriormente, el ruido puede ser tal que nuestras mediciones no sean la imagen de ningún modelo, de modo que no podemos buscar un modelo que produzca los datos sino más bien buscar el mejor (u óptimo) modelo : es decir, el que mejor se ajusta a los datos. Esto nos lleva a minimizar una función objetivo , es decir, una funcional que cuantifica qué tan grandes son los residuos o qué tan lejos están los datos predichos de los datos observados. Por supuesto, cuando tenemos datos perfectos (es decir, sin ruido), entonces el modelo recuperado debería ajustarse perfectamente a los datos observados. Una función objetivo estándar, , es de la forma: donde es la norma euclidiana (será la norma cuando las mediciones sean funciones en lugar de muestras) de los residuos. Este enfoque equivale a hacer uso de mínimos cuadrados ordinarios , un enfoque ampliamente utilizado en estadística. Sin embargo, se sabe que la norma euclidiana es muy sensible a los valores atípicos: para evitar esta dificultad podemos pensar en usar otras distancias, por ejemplo la norma, en reemplazo de la norma. ${\displaystyle \varphi }$ ${\displaystyle \varphi (p)=\|Fp-d_{\text{obs}}\|^{2}}$ ${\displaystyle \|\cdot \|}$ ${\displaystyle L^{2}}$ ${\displaystyle L^{1}}$ ${\displaystyle L^{2}}$

Enfoque bayesiano

Muy similar al enfoque de mínimos cuadrados es el enfoque probabilístico: si conocemos las estadísticas del ruido que contamina los datos, podemos pensar en buscar el modelo más probable m, que es el modelo que cumple el criterio de máxima verosimilitud . Si el ruido es gaussiano , el criterio de máxima verosimilitud aparece como un criterio de mínimos cuadrados, siendo reemplazado el producto escalar euclidiano en el espacio de datos por un producto escalar que involucra la covarianza del ruido. Además, si se dispone de información previa sobre los parámetros del modelo, podríamos pensar en utilizar la inferencia bayesiana para formular la solución del problema inverso. Este enfoque se describe en detalle en el libro de Tarantola. ^[9]

Solución numérica de nuestro ejemplo elemental

Aquí utilizamos la norma euclidiana para cuantificar los desajustes de los datos. Como se trata de un problema inverso lineal, la función objetivo es cuadrática. Para su minimización, es clásico calcular su gradiente utilizando el mismo razonamiento (como lo haríamos para minimizar una función de una sola variable). En el modelo óptimo , este gradiente se desvanece, lo que se puede escribir como: donde F ^T denota la matriz transpuesta de F . Esta ecuación se simplifica a: ${\displaystyle p_{\text{opt}}}$ $${\displaystyle \nabla _{p}\varphi =2(F^{\mathrm {T} }Fp_{\text{opt}}-F^{\mathrm {T} }d_{\text{obs}})=0}$$ $${\displaystyle F^{\mathrm {T} }Fp_{\text{opt}}=F^{\mathrm {T} }d_{\text{obs}}}$$

Esta expresión se conoce como ecuación normal y nos da una posible solución al problema inverso. En nuestro ejemplo, la matriz resulta ser, en general, de rango completo, de modo que la ecuación anterior tiene sentido y determina de forma única los parámetros del modelo: no necesitamos integrar información adicional para llegar a una solución única. ${\displaystyle F^{\mathrm {T} }F}$

Aspectos matemáticos y computacionales

Los problemas inversos suelen estar mal planteados, a diferencia de los problemas bien planteados que suelen cumplirse en el modelado matemático. De las tres condiciones para un problema bien planteado sugeridas por Jacques Hadamard (existencia, unicidad y estabilidad de la solución o soluciones), la condición de estabilidad es la que se viola con mayor frecuencia. En el sentido del análisis funcional , el problema inverso se representa mediante una aplicación entre espacios métricos . Si bien los problemas inversos a menudo se formulan en espacios de dimensión infinita, las limitaciones a un número finito de mediciones y la consideración práctica de recuperar solo un número finito de parámetros desconocidos pueden llevar a que los problemas se reformulen en forma discreta. En este caso, el problema inverso normalmente estará mal condicionado . En estos casos, la regularización se puede utilizar para introducir suposiciones suaves sobre la solución y evitar el sobreajuste . Muchos casos de problemas inversos regularizados se pueden interpretar como casos especiales de inferencia bayesiana . ^[10]

Solución numérica del problema de optimización

Algunos problemas inversos tienen una solución muy simple, por ejemplo, cuando uno tiene un conjunto de funciones unisolventes , es decir , un conjunto de funciones ${\displaystyle n}$ tales que evaluarlas en puntos distintos ${\displaystyle n}$ produce un conjunto de vectores linealmente independientes . Esto significa que dada una combinación lineal de estas funciones, los coeficientes se pueden calcular organizando los vectores como las columnas de una matriz y luego invirtiendo esta matriz. El ejemplo más simple de funciones unisolventes son los polinomios construidos, utilizando el teorema de unisolvencia , de modo que sean unisolventes. Concretamente, esto se hace invirtiendo la matriz de Vandermonde . Pero esta es una situación muy específica.

En general, la solución de un problema inverso requiere algoritmos de optimización sofisticados. Cuando el modelo se describe mediante un gran número de parámetros (el número de incógnitas implicadas en algunas aplicaciones de tomografía por difracción puede alcanzar los mil millones), la resolución del sistema lineal asociado a las ecuaciones normales puede resultar engorrosa. El método numérico que se utilice para resolver el problema de optimización depende en particular del coste necesario para calcular la solución del problema directo. Una vez elegido el algoritmo adecuado para resolver el problema directo (una simple multiplicación matriz-vector puede no ser adecuada cuando la matriz es enorme), el algoritmo adecuado para llevar a cabo la minimización se puede encontrar en libros de texto que tratan de métodos numéricos para la solución de sistemas lineales y para la minimización de funciones cuadráticas (véase, por ejemplo, Ciarlet ^[11] o Nocedal ^[12] ). ${\displaystyle Fp}$ ${\displaystyle F}$

Además, el usuario puede desear añadir restricciones físicas a los modelos: en este caso, deben estar familiarizados con los métodos de optimización con restricciones , un tema en sí mismo. En todos los casos, calcular el gradiente de la función objetivo suele ser un elemento clave para la solución del problema de optimización. Como se mencionó anteriormente, la información sobre la distribución espacial de un parámetro distribuido se puede introducir a través de la parametrización. También se puede pensar en adaptar esta parametrización durante la optimización. ^[13]

Si la función objetivo se basa en una norma distinta a la norma euclidiana, debemos abandonar el área de optimización cuadrática. Como resultado, el problema de optimización se vuelve más difícil. En particular, cuando la norma se utiliza para cuantificar el desajuste de los datos, la función objetivo ya no es diferenciable: su gradiente ya no tiene sentido. En este caso, se utilizan métodos específicos (véase, por ejemplo, Lemaréchal ^[14] ) de optimización no diferenciable. ${\displaystyle L^{1}}$

Una vez calculado el modelo óptimo, debemos responder a la pregunta: "¿Podemos confiar en este modelo?" La pregunta puede formularse de la siguiente manera: ¿Qué tan grande es el conjunto de modelos que coinciden con los datos "casi tan bien" como este modelo? En el caso de las funciones objetivo cuadráticas, este conjunto está contenido en un hiperelipsoide, un subconjunto de ( es el número de incógnitas), cuyo tamaño depende de lo que entendamos por "casi tan bien", es decir, del nivel de ruido. La dirección del eje más grande de este elipsoide ( vector propio asociado con el valor propio más pequeño de la matriz ) es la dirección de los componentes mal determinados: si seguimos esta dirección, podemos introducir una fuerte perturbación en el modelo sin cambiar significativamente el valor de la función objetivo y, por lo tanto, terminar con un modelo cuasi-óptimo significativamente diferente. Vemos claramente que la respuesta a la pregunta "¿podemos confiar en este modelo?" está determinada por el nivel de ruido y por los valores propios de la hessiana de la función objetivo o, equivalentemente, en el caso en que no se haya integrado ninguna regularización, por los valores singulares de la matriz . Por supuesto, el uso de la regularización (u otros tipos de información previa) reduce el tamaño del conjunto de soluciones casi óptimas y, a su vez, aumenta la confianza que podemos depositar en la solución calculada. ${\displaystyle R^{M}}$ ${\displaystyle M}$ ${\displaystyle F^{T}F}$ ${\displaystyle F}$

Estabilidad, regularización y discretización de modelos en dimensión infinita

Aquí nos centramos en la recuperación de un parámetro distribuido. Cuando buscamos parámetros distribuidos tenemos que discretizar estas funciones desconocidas. Al hacerlo, reducimos la dimensión del problema a algo finito. Pero ahora, la pregunta es: ¿hay algún vínculo entre la solución que calculamos y la del problema inicial? Luego, otra pregunta: ¿qué queremos decir con la solución del problema inicial? Dado que un número finito de datos no permite la determinación de una infinidad de incógnitas, la función de desajuste de los datos originales tiene que ser regularizada para asegurar la unicidad de la solución. Muchas veces, reducir las incógnitas a un espacio de dimensión finita proporcionará una regularización adecuada: la solución calculada se verá como una versión discreta de la solución que estábamos buscando. Por ejemplo, una discretización ingenua a menudo funcionará para resolver el problema de deconvolución : funcionará siempre que no permitamos que aparezcan frecuencias faltantes en la solución numérica. Pero muchas veces, la regularización tiene que ser integrada explícitamente en la función objetivo.

Para entender lo que puede suceder, debemos tener en cuenta que resolver un problema inverso lineal de este tipo equivale a resolver una ecuación integral de Fredholm de primer tipo: $${\displaystyle d(x)=\int _{\Omega }K(x,y)p(y)dy}$$

donde es el núcleo, y son vectores de , y es un dominio en . Esto es válido para una aplicación 2D. Para una aplicación 3D, consideramos . Nótese que aquí los parámetros del modelo consisten en una función y que la respuesta de un modelo también consiste en una función denotada por . Esta ecuación es una extensión a dimensión infinita de la ecuación matricial dada en el caso de problemas discretos. ${\displaystyle K}$ ${\displaystyle x}$ ${\displaystyle y}$ ${\displaystyle R^{2}}$ ${\displaystyle \Omega }$ ${\displaystyle R^{2}}$ ${\displaystyle x,y\in R^{3}}$ ${\displaystyle p}$ ${\displaystyle d(x)}$ ${\displaystyle d=Fp}$

Para que sea suficientemente suave, el operador definido anteriormente es compacto en espacios de Banach razonables como el . La teoría de F. Riesz establece que el conjunto de valores singulares de dicho operador contiene cero (de ahí la existencia de un espacio nulo), es finito o, como mucho, numerable y, en el último caso, constituyen una secuencia que tiende a cero. En el caso de un núcleo simétrico, tenemos una infinidad de valores propios y los vectores propios asociados constituyen una base hilbertiana de . Por lo tanto, cualquier solución de esta ecuación está determinada hasta una función aditiva en el espacio nulo y, en el caso de infinidad de valores singulares, la solución (que implica el recíproco de valores propios pequeños arbitrarios) es inestable: ¡dos ingredientes que hacen que la solución de esta ecuación integral sea un típico problema mal planteado! Sin embargo, podemos definir una solución a través de la pseudoinversa de la función directa (de nuevo hasta una función aditiva arbitraria). Cuando la función forward es compacta, la regularización clásica de Tikhonov funcionará si la usamos para integrar información previa que indique que la norma de la solución debe ser lo más pequeña posible: esto hará que el problema inverso esté bien planteado. Sin embargo, como en el caso de dimensión finita, tenemos que cuestionar la confianza que podemos depositar en la solución calculada. Nuevamente, básicamente, la información se encuentra en los valores propios del operador hessiano. Si se exploran subespacios que contienen vectores propios asociados con valores propios pequeños para calcular la solución, entonces la solución difícilmente será confiable: algunos de sus componentes estarán mal determinados. El valor propio más pequeño es igual al peso introducido en la regularización de Tikhonov. ${\displaystyle K}$ ${\displaystyle L^{2}}$ ${\displaystyle L^{2}}$ ${\displaystyle L^{2}}$

Los núcleos irregulares pueden producir una función hacia delante que no sea compacta e incluso ilimitada si equipamos ingenuamente el espacio de modelos con la norma. En tales casos, el hessiano no es un operador acotado y la noción de valor propio ya no tiene sentido. Se requiere un análisis matemático para convertirlo en un operador acotado y diseñar un problema bien planteado: se puede encontrar una ilustración en ^[15] . Nuevamente, tenemos que cuestionar la confianza que podemos depositar en la solución calculada y tenemos que generalizar la noción de valor propio para obtener la respuesta. ^[16] ${\displaystyle L^{2}}$

El análisis del espectro del operador hessiano es, por tanto, un elemento clave para determinar la fiabilidad de la solución calculada. Sin embargo, un análisis de este tipo suele ser una tarea muy pesada. Esto ha llevado a varios autores a investigar enfoques alternativos en el caso de que no estemos interesados ​​en todos los componentes de la función desconocida, sino solo en las subincógnitas que son las imágenes de la función desconocida por un operador lineal. Estos enfoques se conocen como el "método de Backus y Gilbert ^[17] ", el enfoque de los centinelas de Lions ^[18] y el método SOLA ^[19] : estos enfoques resultaron estar fuertemente relacionados entre sí, como se explica en Chavent ^[20]. Finalmente, el concepto de resolución limitada , a menudo invocado por los físicos, no es más que una visión específica del hecho de que algunos componentes mal determinados pueden corromper la solución. Pero, en términos generales, estos componentes mal determinados del modelo no están necesariamente asociados a altas frecuencias.

Algunos problemas inversos lineales clásicos para la recuperación de parámetros distribuidos

Los problemas mencionados a continuación corresponden a diferentes versiones de la integral de Fredholm: cada una de ellas está asociada a un núcleo específico . ${\displaystyle K}$

Desconvolución

El objetivo de la deconvolución es reconstruir la imagen o señal original que aparece ruidosa y borrosa en los datos . ^[21] Desde un punto de vista matemático, el núcleo aquí solo depende de la diferencia entre y . ${\displaystyle p(x)}$ ${\displaystyle d(x)}$ ${\displaystyle K(x,y)}$ ${\displaystyle x}$ ${\displaystyle y}$

Métodos tomográficos

En estos métodos intentamos recuperar un parámetro distribuido, la observación consiste en la medición de las integrales de este parámetro realizada a lo largo de una familia de líneas. Denotamos por la línea en esta familia asociada con el punto de medición . La observación en puede escribirse así como: donde es la longitud de arco a lo largo de y una función de ponderación conocida. Comparando esta ecuación con la integral de Fredholm anterior, notamos que el núcleo es una especie de función delta que alcanza su pico en la línea . Con un núcleo de este tipo, la función directa no es compacta. ${\displaystyle \Gamma _{x}}$ ${\displaystyle x}$ ${\displaystyle x}$ $${\displaystyle d(x)=\int _{\Gamma _{x}}w(x,y)p(y)\,dy}$$ ${\displaystyle s}$ ${\displaystyle {\Gamma _{x}}}$ ${\displaystyle w(x,y)}$ ${\displaystyle K(x,y)}$ ${\displaystyle {\Gamma _{x}}}$

Tomografía computarizada

En la tomografía computarizada de rayos X, las líneas en las que se integra el parámetro son líneas rectas: la reconstrucción tomográfica de la distribución de parámetros se basa en la inversión de la transformada de Radon . Aunque desde un punto de vista teórico se entienden bien muchos problemas de inversión lineal, los problemas que involucran la transformada de Radon y sus generalizaciones aún presentan muchos desafíos teóricos con cuestiones de suficiencia de datos aún sin resolver. Tales problemas incluyen datos incompletos para la transformada de rayos X en tres dimensiones y problemas que involucran la generalización de la transformada de rayos X a campos tensoriales. Las soluciones exploradas incluyen la Técnica de Reconstrucción Algebraica , la retroproyección filtrada y, a medida que el poder de cómputo ha aumentado, los métodos de reconstrucción iterativos como la Varianza Mínima Asintótica Dispersa iterativa . ^[22]

Tomografía por difracción

La tomografía de difracción es un problema inverso lineal clásico en sismología de exploración: la amplitud registrada en un momento dado para un par fuente-receptor es la suma de las contribuciones que surgen de puntos tales que la suma de las distancias, medidas en tiempos de viaje, desde la fuente y el receptor, respectivamente, es igual al tiempo de registro correspondiente. En 3D, el parámetro no se integra a lo largo de líneas sino sobre superficies. Si la velocidad de propagación es constante, dichos puntos se distribuyen en un elipsoide. El problema inverso consiste en recuperar la distribución de puntos de difracción de los sismogramas registrados a lo largo del estudio, siendo conocida la distribución de velocidad. Una solución directa ha sido propuesta originalmente por Beylkin y Lambaré et al. ^[23]: estos trabajos fueron los puntos de partida de los enfoques conocidos como migración de amplitud preservada (ver Beylkin ^{[24] [25]} y Bleistein ^[26] ). Si se utilizan técnicas de óptica geométrica (es decir, rayos) para resolver la ecuación de onda, estos métodos resultan estar estrechamente relacionados con los llamados métodos de migración de mínimos cuadrados ^[27] derivados del enfoque de mínimos cuadrados (ver Lailly, ^[28] Tarantola ^[29] ).

Tomografía Doppler (astrofísica)

Si consideramos un objeto estelar en rotación, las líneas espectrales que podemos observar en un perfil espectral se desplazarán debido al efecto Doppler. La tomografía Doppler tiene como objetivo convertir la información contenida en el seguimiento espectral del objeto en una imagen 2D de la emisión (en función de la velocidad radial y de la fase en el movimiento periódico de rotación) de la atmósfera estelar. Como explica Tom Marsh ^[30], este problema inverso lineal es similar a la tomografía: tenemos que recuperar un parámetro distribuido que se ha integrado a lo largo de las líneas para producir sus efectos en los registros.

Conducción de calor inversa

Las primeras publicaciones sobre conducción térmica inversa surgieron a partir de la determinación del flujo de calor superficial durante el reingreso atmosférico a partir de sensores de temperatura enterrados. ^{[31] [32]} Otras aplicaciones en las que se necesita el flujo de calor superficial pero los sensores de superficie no son prácticos incluyen: dentro de motores alternativos, dentro de motores de cohetes; y, pruebas de componentes de reactores nucleares. ^[33] Se han desarrollado diversas técnicas numéricas para abordar la mala postura y la sensibilidad al error de medición causado por la amortiguación y el retraso en la señal de temperatura. ^{[34] [35] [36]}

Problemas inversos no lineales

Los problemas inversos no lineales constituyen una familia inherentemente más difícil de problemas inversos. En este caso, la función directa es un operador no lineal. La modelización de fenómenos físicos a menudo se basa en la solución de una ecuación diferencial parcial (véase la tabla anterior, excepto para la ley de la gravedad): aunque estas ecuaciones diferenciales parciales suelen ser lineales, los parámetros físicos que aparecen en estas ecuaciones dependen de forma no lineal del estado del sistema y, por tanto, de las observaciones que hagamos sobre él. ${\displaystyle F}$

Algunos problemas inversos no lineales clásicos

Problemas de dispersión inversa

Mientras que los problemas inversos lineales se resolvieron completamente desde el punto de vista teórico a finales del siglo XIX ^{[ cita requerida ]} , solo una clase de problemas inversos no lineales lo fue antes de 1970, los problemas espectrales inversos y de dispersión inversa (de una dimensión espacial) , después del trabajo seminal de la escuela matemática rusa ( Krein , Gelfand , Levitan, Marchenko ). Chadan y Sabatier han realizado una amplia revisión de los resultados en su libro "Inverse Problems of Quantum Scattering Theory" (dos ediciones en inglés, una en ruso).

En este tipo de problemas, los datos son propiedades del espectro de un operador lineal que describen la dispersión. El espectro está formado por valores propios y funciones propias , que forman juntos el "espectro discreto", y generalizaciones, llamadas espectro continuo. El punto físico más notable es que los experimentos de dispersión dan información sólo sobre el espectro continuo, y que conocer su espectro completo es necesario y suficiente para recuperar el operador de dispersión. Por lo tanto, tenemos parámetros invisibles, mucho más interesantes que el espacio nulo que tiene una propiedad similar en problemas lineales inversos. Además, hay movimientos físicos en los que el espectro de dicho operador se conserva como consecuencia de dicho movimiento. Este fenómeno está gobernado por ecuaciones de evolución diferencial parcial no lineales especiales, por ejemplo la ecuación de Korteweg-de Vries . Si el espectro del operador se reduce a un único valor propio, su movimiento correspondiente es el de una única protuberancia que se propaga a velocidad constante y sin deformación, una onda solitaria llamada " solitón ".

La señal perfecta y sus generalizaciones para la ecuación de Korteweg-de Vries u otras ecuaciones diferenciales parciales no lineales integrables son de gran interés y tienen muchas aplicaciones posibles. Esta área se ha estudiado como una rama de la física matemática desde la década de 1970. Los problemas inversos no lineales también se estudian actualmente en muchos campos de la ciencia aplicada (acústica, mecánica, mecánica cuántica, dispersión electromagnética, en particular sondeos de radar, sondeos sísmicos y casi todas las modalidades de imagen).

Un último ejemplo relacionado con la hipótesis de Riemann fue dado por Wu y Sprung: la idea es que en la antigua teoría cuántica semiclásica , la inversa del potencial dentro del hamiltoniano es proporcional a la semiderivada de la función de conteo de valores propios (energías)  n ( x ).

Adaptación de la permeabilidad en yacimientos de petróleo y gas

El objetivo es recuperar el coeficiente de difusión en la ecuación diferencial parcial parabólica que modela flujos de fluidos monofásicos en medios porosos. Este problema ha sido objeto de numerosos estudios desde un trabajo pionero realizado a principios de los años setenta. ^[37] En el caso de los flujos bifásicos, un problema importante es la estimación de las permeabilidades relativas y las presiones capilares. ^[38]

Problemas inversos en las ecuaciones de onda

El objetivo es recuperar las velocidades de las ondas (ondas P y S) y las distribuciones de densidad a partir de sismogramas . Estos problemas inversos son de gran interés en sismología y geofísica de exploración . Básicamente, podemos considerar dos modelos matemáticos:

• La ecuación de onda acústica (en la que se ignoran las ondas S cuando las dimensiones del espacio son 2 o 3)
• La ecuación de elastodinámica en la que las velocidades de las ondas P y S pueden derivarse de los parámetros de Lamé y la densidad.

Estas ecuaciones hiperbólicas básicas se pueden mejorar incorporando atenuación , anisotropía , ...

La solución del problema inverso en la ecuación de onda 1D ha sido objeto de muchos estudios. Es uno de los pocos problemas inversos no lineales para los que podemos probar la unicidad de la solución. ^[8] El análisis de la estabilidad de la solución fue otro desafío. ^[39] Se desarrollaron aplicaciones prácticas, utilizando el enfoque de mínimos cuadrados. ^{[39] [40]} La extensión a problemas 2D o 3D y a las ecuaciones de elastodinámica se intentó desde los años 80, ¡pero resultó ser muy difícil! Este problema, a menudo denominado inversión de forma de onda completa (FWI), aún no está completamente resuelto: entre las principales dificultades se encuentran la existencia de ruido no gaussiano en los sismogramas, problemas de salto de ciclo (también conocido como ambigüedad de fase) y el comportamiento caótico de la función de desajuste de datos. ^[41] Algunos autores han investigado la posibilidad de reformular el problema inverso para hacer que la función objetivo sea menos caótica que la función de desajuste de datos. ^{[42] [43]}

Tomografía de tiempo de viaje

Al darse cuenta de la dificultad del problema inverso en la ecuación de onda, los sismólogos investigaron un enfoque simplificado que hace uso de la óptica geométrica. En particular, se propusieron invertir la distribución de la velocidad de propagación, conociendo los tiempos de llegada de los frentes de onda observados en los sismogramas. Estos frentes de onda pueden estar asociados a llegadas directas o a reflexiones asociadas a reflectores cuya geometría se debe determinar, junto con la distribución de la velocidad.

La distribución del tiempo de llegada ( es un punto en el espacio físico) de un frente de onda emitido desde una fuente puntual, satisface la ecuación de Eikonal : donde denota la distribución de lentitud (recíproca de la velocidad). La presencia de hace que esta ecuación no sea lineal. Se resuelve clásicamente disparando rayos (trayectorias sobre las cuales el tiempo de llegada es estacionario) desde la fuente puntual. ${\displaystyle {\tau }(x)}$ ${\displaystyle x}$ $${\displaystyle \|\nabla \tau (x)\|=s(x),}$$ ${\displaystyle s(x)}$ ${\displaystyle \|\cdot \|}$

Este problema es similar a la tomografía: los tiempos de llegada medidos son la integral a lo largo de la trayectoria del rayo de la lentitud. Pero este problema similar a la tomografía es no lineal, principalmente porque la geometría desconocida de la trayectoria del rayo depende de la distribución de la velocidad (o lentitud). A pesar de su carácter no lineal, la tomografía de tiempo de viaje resultó ser muy eficaz para determinar la velocidad de propagación en la Tierra o en el subsuelo, siendo este último aspecto un elemento clave para la obtención de imágenes sísmicas, en particular utilizando los métodos mencionados en la Sección "Tomografía por difracción".

Aspectos matemáticos: las preguntas de Hadamard

Las preguntas se refieren a la precisión del planteamiento: ¿el problema de mínimos cuadrados tiene una solución única que depende continuamente de los datos (problema de estabilidad)? Es la primera pregunta, pero también es difícil debido a la no linealidad de . Para ver de dónde surgen las dificultades, Chavent ^[44] propuso dividir conceptualmente la minimización de la función de desajuste de datos en dos pasos consecutivos ( es el subconjunto de modelos admisibles): ${\displaystyle F}$ ${\displaystyle P_{\text{adm}}}$

• Paso de proyección: dada, encuentre una proyección sobre (punto más cercano sobre de acuerdo con la distancia involucrada en la definición de la función objetivo) ${\displaystyle d_{\text{obs}}}$ ${\displaystyle F(P_{\text{adm}})}$ ${\displaystyle F(P_{\text{adm}})}$
• Dada esta proyección, encuentre una preimagen que sea un modelo cuya imagen por operador sea esta proyección. ${\displaystyle F}$

Pueden surgir dificultades (y normalmente surgirán) en ambos pasos:

1. No es probable que el operador sea uno a uno, por lo tanto, puede haber más de una preimagen, ${\displaystyle F}$
2. Incluso cuando es uno a uno, su inversa puede no ser continua en , ${\displaystyle F}$ ${\displaystyle F(P)}$
3. La proyección sobre puede no existir, si este conjunto no es cerrado, ${\displaystyle F(P_{\text{adm}})}$
4. La proyección sobre puede no ser única ni continua, ya que puede ser no convexa debido a la no linealidad de . ${\displaystyle F(P_{\text{adm}})}$ ${\displaystyle F}$

Nos remitimos a Chavent ^[44] para un análisis matemático de estos puntos.

Aspectos computacionales

Una función de desajuste de datos no convexos

Como el mapa de avance no es lineal, es probable que la función de desajuste de datos no sea convexa, lo que hace que las técnicas de minimización local sean ineficientes. Se han investigado varios enfoques para superar esta dificultad:

• uso de técnicas de optimización global como el muestreo de la función de densidad posterior y el algoritmo Metropolis en el marco probabilístico del problema inverso, ^[45] algoritmos genéticos (solos o en combinación con el algoritmo Metropolis: ver ^[46] para una aplicación a la determinación de permeabilidades que coinciden con los datos de permeabilidad existentes), redes neuronales, técnicas de regularización incluyendo análisis de múltiples escalas;
• reformulación de la función objetivo de mínimos cuadrados para hacerla más suave (véase ^{[42] [43]} para el problema inverso en las ecuaciones de onda).

Cálculo del gradiente de la función objetivo

Los problemas inversos, especialmente en dimensión infinita, pueden ser de gran tamaño, por lo que requieren un tiempo de cálculo importante. Cuando la función directa no es lineal, las dificultades computacionales aumentan y minimizar la función objetivo puede resultar difícil. Al contrario de lo que ocurre en la situación lineal, un uso explícito de la matriz de Hesse para resolver las ecuaciones normales no tiene sentido en este caso: la matriz de Hesse varía con los modelos. La evaluación del gradiente de la función objetivo para algunos modelos es mucho más eficaz. Se puede ahorrar un importante esfuerzo computacional cuando podemos evitar el cálculo muy pesado del jacobiano (a menudo llamado " derivadas de Fréchet "): el método de estados adjuntos, propuesto por Chavent y Lions, ^[47] tiene como objetivo evitar este cálculo muy pesado. En la actualidad, se utiliza ampliamente. ^[48]

Un algoritmo de inversión (publicado bajo una licencia Creative Commons, CC BY-NC-ND por Elsevier) ^[49]

Aplicaciones

La teoría del problema inverso se utiliza ampliamente en predicciones meteorológicas, oceanografía, hidrología e ingeniería petrolera. ^{[50] [51] [52]} Otra aplicación es la inversión de ondas elásticas para la caracterización no destructiva de estructuras de ingeniería. ^[49]

También se encuentran problemas inversos en el campo de la transferencia de calor, donde se estima un flujo de calor superficial ^[53] saliente a partir de datos de temperatura medidos dentro de un cuerpo rígido; y en la comprensión de los controles sobre la descomposición de la materia vegetal. ^[54] El problema inverso lineal también es fundamental para la estimación espectral y la estimación de la dirección de llegada (DOA) en el procesamiento de señales .

La litografía inversa se utiliza en el diseño de fotomáscaras para la fabricación de dispositivos semiconductores .

Véase también

• Problema inverso atmosférico  : medición de la distribución vertical de las propiedades físicas de la columna atmosféricaPages displaying wikidata descriptions as a fallback
• Método Backus-Gilbert
• Tomografía computarizada  : procedimiento de obtención de imágenes médicas que utiliza rayos X para producir imágenes transversales.Pages displaying short descriptions of redirect targets
  • Técnica de reconstrucción algebraica  – Técnica en tomografía computarizada
  • Retroproyección filtrada  – Transformación integralPages displaying short descriptions of redirect targets
  • Reconstrucción iterativa
• Asimilación de datos  – Técnica para actualizar el modelo numérico con datos observados
• Optimización de ingeniería  – Técnicas de optimización
• Modelo de caja gris  : modelo de producción de datos matemáticos con estructura limitada
• Geofísica matemática  – Aplicación de las matemáticas en el campo geográficoPages displaying short descriptions of redirect targets
• Estimación óptima  : estadística aplicada, un método inverso de matriz regularizada basado en el teorema de BayesPages displaying wikidata descriptions as a fallback
• Inversión sísmica
• Regularización de Tikhonov  : técnica de regularización para problemas mal planteadosPages displaying short descriptions of redirect targets
• Detección comprimida  : técnica de procesamiento de señales
• Problema de inducción  – Cuestión de si el razonamiento inductivo conduce al conocimiento definitivo

Revistas académicas

Cuatro revistas académicas principales cubren problemas inversos en general:

• Problemas inversos
• Revista de problemas inversos y mal planteados ^[55]
• Problemas inversos en ciencia e ingeniería ^[56]
• Problemas inversos y formación de imágenes ^[57]

Muchas revistas sobre imágenes médicas, geofísica, pruebas no destructivas, etc. están dominadas por problemas inversos en esas áreas.

Referencias

1. Mohamad-Djafari, Ali (29 de enero de 2013). Problemas inversos en la visión y la tomografía 3D. John Wiley & Sons. ISBN 978-1-118-60046-7.
2. Pizlo, Zygmunt. "La percepción vista como un problema inverso". Vision research 41.24 (2001): 3145-3161.
3. Vito, Ernesto De, et al. "Aprender a partir de ejemplos como un problema inverso". Journal of Machine Learning Research 6 de mayo (2005): 883-904.
4. Cardenas, IC (2019). "Sobre el uso de redes bayesianas como un enfoque de metamodelado para analizar incertidumbres en el análisis de estabilidad de taludes". Georisk: evaluación y gestión del riesgo para sistemas de ingeniería y geopeligros . 13 (1): 53–65. Bibcode :2019GAMRE..13...53C. doi :10.1080/17499518.2018.1498524. S2CID  216590427.
5. Weyl, Hermann (1911). "Über die asymptotische Verteilung der Eigenwerte". Nachrichten der Königlichen Gesellschaft der Wissenschaften zu Göttingen : 110-117. Archivado desde el original el 1 de agosto de 2013 . Consultado el 14 de mayo de 2018 .
6. » Epílogo — El artículo de Ambartsumian Viktor Ambartsumian
7. Ambartsumian, Rouben V. (1998). "Una vida en astrofísica. Artículos seleccionados de Viktor A. Ambartsumian". Astrofísica . 41 (4): 328–330. doi :10.1007/BF02894658. S2CID  118952753.
8. Burridge, Robert (1980). "Las ecuaciones integrales de Gelfand-Levitan, Marchenko y Gopinath-Sondhi de la teoría de dispersión inversa, consideradas en el contexto de problemas de impulso-respuesta inversa". Wave Motion . 2 (4): 305–323. Bibcode :1980WaMot...2..305B. doi :10.1016/0165-2125(80)90011-6.
9. Tarantola, Albert (1987). Teoría del problema inverso (1.ª ed.). Elsevier. ISBN 9780444599674.
10. Tarantola, Albert (2005). "Front Matter" (PDF) . Teoría del problema inverso y métodos para la estimación de parámetros de modelos . SIAM. pp. i–xii. doi :10.1137/1.9780898717921.fm. ISBN. 978-0-89871-572-9.
11. Ciarlet, Philippe (1994). Introducción al análisis numérico matricielle y a la optimización . París: Masson. ISBN 9782225688935.
12. Nocedal, Jorge (2006). Optimización numérica . Springer.
13. Ben Ameur, Hend; Chavent, Guy; Jaffré, Jérôme (2002). "Indicadores de refinamiento y engrosamiento para la parametrización adaptativa: aplicación a la estimación de transmisividades hidráulicas" (PDF) . Problemas inversos . 18 (3): 775–794. Bibcode :2002InvPr..18..775B. doi :10.1088/0266-5611/18/3/317. S2CID  250892174.
14. Lemaréchal, Claude (1989). Optimización, Manuales de investigación de operaciones y ciencia de la gestión . Elsevier. págs. 529–572.
15. Delprat-Jannaud, Florence; Lailly, Patrick (1993). "Formulaciones mal planteadas y bien planteadas del problema de la tomografía del tiempo de recorrido de la reflexión". Revista de investigación geofísica . 98 (B4): 6589–6605. Código Bibliográfico :1993JGR....98.6589D. doi :10.1029/92JB02441.
16. Delprat-Jannaud, Florence; Lailly, Patrick (1992). "¿Qué información sobre el modelo terrestre proporcionan los tiempos de viaje de reflexión?". Journal of Geophysical Research . 98 (B13): 827–844. Bibcode :1992JGR....9719827D. doi :10.1029/92JB01739.
17. Backus, George; Gilbert, Freeman (1968). "El poder de resolución de los datos brutos de la Tierra". Revista geofísica de la Royal Astronomical Society . 16 (10): 169–205. Código Bibliográfico :1968GeoJ...16..169B. doi : 10.1111/j.1365-246X.1968.tb00216.x .
18. Leones, Jacques Louis (1988). "Sur les sentinelles des systèmes distribués". CR Acad. Ciencia. París . I Matemáticas: 819–823.
19. Pijpers, Frank; Thompson, Michael (1993). "El método SOLA para la inversión heliosísmica". Astronomía y Astrofísica . 281 (12): 231–240. Bibcode :1994A&A...281..231P.
20. Chavent, chico (1998). Mínimos cuadrados, centinelas y promedio sustractivo óptimamente localizado en ecuaciones aux derivadas parciales y aplicaciones. París: Gauthier Villars. págs. 345–356.
21. Kaipio, J., y Somersalo, E. (2010). Problemas inversos estadísticos y computacionales. Nueva York, NY: Springer.
22. Abeida, Habti; Zhang, Qilin; Li, Jian; Merabtine, Nadjim (2013). "Enfoques iterativos dispersos asintóticos basados ​​en varianza mínima para procesamiento de matrices" (PDF) . IEEE Transactions on Signal Processing . 61 (4): 933–944. arXiv : 1802.03070 . Bibcode :2013ITSP...61..933A. doi :10.1109/tsp.2012.2231676. ISSN  1053-587X. S2CID  16276001.
23. Lambaré, Gilles; Virieux, Jean; Madariaga, Raul; Jin, Side (1992). "Inversión asintótica iterativa en la aproximación acústica". Geofísica . 57 (9): 1138–1154. Bibcode :1992Geop...57.1138L. doi :10.1190/1.1443328. S2CID  55836067.
24. Beylkin, Gregory (1984). "El problema de la inversión y las aplicaciones de la transformada generalizada de Radon" (PDF) . Comunicaciones sobre Matemática Pura y Aplicada . XXXVII (5): 579–599. doi :10.1002/cpa.3160370503.
25. Beylkin, Gregory (1985). "Imágenes de discontinuidades en el problema de dispersión inversa mediante la inversión de una transformada causal generalizada de Radon". J. Math. Phys . 26 (1): 99–108. Bibcode :1985JMP....26...99B. doi :10.1063/1.526755.
26. Bleistein, Norman (1987). "Sobre la formación de imágenes de reflectores en la tierra". Geofísica . 52 (7): 931–942. Bibcode :1987Geop...52..931B. doi :10.1190/1.1442363. S2CID  5095133.
27. Nemeth, Tamas; Wu, Chengjun; Schuster, Gerard (1999). "Migración de datos de reflexión incompletos por mínimos cuadrados" (PDF) . Geofísica . 64 (1): 208–221. Bibcode :1999Geop...64..208N. doi :10.1190/1.1444517.
28. Lailly, Patrick (1983). El problema sísmico inverso como una secuencia de migraciones previas a la pila . Filadelfia: SIAM. pp. 206–220. ISBN 0-89871-190-8.
29. Tarantola, Albert (1984). "Inversión de datos de reflexión sísmica en la aproximación acústica". Geofísica . 49 (8): 1259–1266. Bibcode :1984Geop...49.1259T. doi :10.1190/1.1441754. S2CID  7596552.
30. Marsh, Tom (2005). "Tomografía Doppler". Astrofísica y ciencia espacial . 296 (1–4): 403–415. arXiv : astro-ph/0011020 . Código Bibliográfico :2005Ap&SS.296..403M. doi :10.1007/s10509-005-4859-3. S2CID  15334110.
31. Shumakov, NV (1957). "Un método para el estudio experimental del proceso de calentamiento de un cuerpo sólido". Física soviética – Física técnica (traducido por el Instituto Americano de Física) . 2 : 771.
32. Stolz, G. Jr. (1960). "Soluciones numéricas a un problema inverso de conducción de calor para formas simples". Journal of Heat Transfer . 82 : 20–26. doi :10.1115/1.3679871.
33. Beck, JV; Blackwell, B.; St. Clair, CR Jr. (1985). Conducción inversa del calor. Problemas mal planteados . Nueva York: J. Wiley & Sons. ISBN 0471083194.
34. Beck, JV; Blackwell, B.; Haji-Sheikh, B. (1996). "Comparación de algunos métodos de conducción de calor inversa utilizando datos experimentales". Revista internacional de transferencia de calor y masa . 39 (17): 3649–3657. Código Bibliográfico :1996IJHMT..39.3649B. doi :10.1016/0017-9310(96)00034-8.
35. Ozisik, MN; Orlande, HRB (2021). Transferencia de calor inversa, fundamentos y aplicaciones (2.ª ed.). CRC Press. ISBN 9780367820671.
36. Manual de ingeniería inversa, editado por KA Woodbury . CRC Press. 2002. ISBN 9780849308611.
37. Chavent, Guy; Lemonnier, Patrick; Dupuy, Michel (1975). "Coincidencia histórica mediante el uso de la teoría de control óptimo". Revista de la Sociedad de Ingenieros Petroleros . 15 (2): 74–86. doi :10.2118/4627-PA.
38. Chavent, Guy; Cohen, Gary; Espy, M. (1980). "Determinación de permeabilidades relativas y presiones capilares mediante un método de ajuste automático". Sociedad de Ingenieros Petroleros (enero). doi :10.2118/9237-MS.
39. Bamberger, Alain; Chavent, Guy; Lailly, Patrick (1979). "Acerca de la estabilidad del problema inverso en la ecuación de onda 1D, aplicación a la interpretación de los perfiles sísmicos". Revista de Matemática Aplicada y Optimización . 5 : 1–47. doi :10.1007/bf01442542. S2CID  122428594.
40. Macé, Danièle; Lailly, Patrick (1986). "Solución del problema inverso unidimensional VSP". Prospección geofísica . 34 (7): 1002–1021. Código Bibliográfico :1986GeopP..34.1002M. doi :10.1111/j.1365-2478.1986.tb00510.x. OSTI  6901651.
41. Virieux, Jean; Operto, Stéphane (2009). "Una visión general de la inversión de forma de onda completa en la geofísica de exploración". Geofísica . 74 (6): WCC1–WCC26. doi :10.1190/1.3238367.
42. Clément, François; Chavent, Guy; Gomez, Suzana (2001). "Inversión de la forma de onda del tiempo de viaje basada en la migración de estructuras simples en 2-D: un ejemplo sintético". Geofísica . 66 (3): 845–860. Bibcode :2001Geop...66..845C. doi :10.1190/1.1444974.
43. Symes, William; Carrazone, Jim (1991). "Inversión de velocidad mediante optimización de semblanza diferencial". Geofísica . 56 (5): 654–663. Bibcode :1991Geop...56..654S. doi :10.1190/1.1443082.
44. Chavent, Guy (2010). Mínimos cuadrados no lineales para problemas inversos . Springer. ISBN 978-90-481-2785-6.
45. Koren, Zvi; Mosegaard, Klaus; Landa, Evgeny; Thore, Pierre; Tarantola, Albert (1991). "Estimación de Monte Carlo y análisis de resolución de velocidades sísmicas de fondo". Revista de investigación geofísica . 96 (B12): 20289–20299. Código Bibliográfico :1991JGR....9620289K. doi :10.1029/91JB02278.
46. Tahmasebi, Pejman; Javadpour, Farzam; Sahimi, Muhammad (agosto de 2016). "Correspondencia estocástica de la permeabilidad de esquisto: caracterización y modelado tridimensional". Revista internacional de geología del carbón . 165 : 231–242. Código Bibliográfico :2016IJCG..165..231T. doi :10.1016/j.coal.2016.08.024.
47. Chavent, chico (1971). Identificación de coeficientes repartidos en ecuaciones aux derivadas partielles . Université Paris 6: Thèse d'Etat.{{cite book}}: CS1 maint: location (link)
48. Plessix, René (2006). "Una revisión del método de estados adjuntos para calcular el gradiente de un funcional con aplicaciones geofísicas". Revista Geofísica Internacional . 167 (2): 495–503. Código Bibliográfico :2006GeoJI.167..495P. doi : 10.1111/j.1365-246X.2006.02978.x .
49. Taslimian, Rohollah; Jaganathan, Arun P. (2024). "Inversión de ondas elásticas circunferenciales para la caracterización de tuberías de hormigón". Wave Motion . 127 : 103272. Código Bibliográfico :2024WaMot.12703272T. doi : 10.1016/j.wavemoti.2024.103272 .
50. Carl Wunsch (13 de junio de 1996). El problema inverso de la circulación oceánica. Cambridge University Press. pp. 9–. ISBN 978-0-521-48090-1.
51. Tahmasebi, Pejman; Javadpour, Farzam; Sahimi, Muhammad (agosto de 2016). "Correspondencia estocástica de la permeabilidad de esquisto: caracterización y modelado tridimensional". Revista internacional de geología del carbón . 165 : 231–242. Código Bibliográfico :2016IJCG..165..231T. doi :10.1016/j.coal.2016.08.024.
52. Knighton, James; Singh, Kanishka; Evaristo, Jaivime (2020). "Comprensión de las estrategias de captación de agua de las raíces forestales a escala de cuenca en los Estados Unidos continentales mediante el modelado ecohidrológico inverso". Geophysical Research Letters . 47 (1): e2019GL085937. Código Bibliográfico :2020GeoRL..4785937K. doi : 10.1029/2019GL085937 . ISSN  1944-8007. S2CID  213914582.
53. Patric Figueiredo (diciembre de 2014). Desarrollo de un método iterativo para resolver problemas multidimensionales de conducción de calor inversa. Lehrstuhl für Wärme- und Stoffübertragung RWTH Aachen.^{[ enlace muerto permanente ]}
54. Forney, David C.; Rothman, Daniel H. (7 de septiembre de 2012). "Estructura común en la heterogeneidad de la descomposición de la materia vegetal". Journal of the Royal Society Interface . 9 (74): 2255–2267. doi :10.1098/rsif.2012.0122. PMC 3405759 . PMID  22535699. 
55. "Revista de problemas inversos y mal planteados". Archivado desde el original el 1 de febrero de 2013.
56. "Problemas inversos en ciencia e ingeniería: vol. 25, n.º 4".
57. "IPI". Archivado desde el original el 11 de octubre de 2006.

Referencias

• Chadan, Khosrow y Sabatier, Pierre Célestin (1977). Problemas inversos en la teoría de la dispersión cuántica . Springer-Verlag. ISBN 0-387-08092-9 
• Aster, Richard; Borchers, Brian y Thurber, Clifford (2018). Estimación de parámetros y problemas inversos , tercera edición, Elsevier. ISBN 9780128134238 , ISBN 9780128134238  
• Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007). "Sección 19.4. Problemas inversos y el uso de información a priori". Recetas numéricas: el arte de la computación científica (3.ª ed.). Nueva York: Cambridge University Press. ISBN 978-0-521-88068-8Archivado desde el original el 11-08-2011 . Consultado el 17-08-2011 .

Lectura adicional

• CW Groetsch (1999). Problemas inversos: actividades para estudiantes universitarios . Cambridge University Press. ISBN 978-0-88385-716-8.

Enlaces externos

• Asociación Internacional de Problemas Inversos Archivado el 15 de junio de 2017 en Wayback Machine.
• Asociación Euroasiática de Problemas Inversos
• Sociedad finlandesa de problemas inversos
• Red de problemas inversos
• El sitio web de Albert Tarantola incluye una versión PDF gratuita de su libro Teoría de problemas inversos y algunos artículos en línea sobre problemas inversos.
• Página de problemas inversos de la Universidad de Alabama Archivado el 5 de abril de 2014 en Wayback Machine.
• Proyecto de problemas inversos y geoestadística Archivado el 2 de noviembre de 2017 en Wayback Machine , Instituto Niels Bohr, Universidad de Copenhague
• Página de recursos de la teoría inversa geofísica de Andy Ganse
• Centro finlandés de excelencia en investigación de problemas inversos