Resolver dos problemas inversos con la misma solución


8

Tengo dos problemas inversos,

A1 x=b1A2 x=b2

Hasta ahora los he estado resolviendo de forma independiente usando la regularización de Tikhonov y obteniendo dos estimaciones para . Sin embargo, en mi caso, representa la misma solución en ambas ecuaciones. ¿Es posible hacer una solución 'simultánea'? Lo ideal sería encontrar la respuesta paraxx

min(A1xb12+A2xb22+Γx2)

Donde Γ=α I e I es la matriz de identidad como en la Regularización de Tikhonov (también conocida como regresión de cresta). Supongo que podría tomar el promedio de ambas soluciones, preguntándome si hay una forma estadísticamente más poderosa de abordar esto.


1
¿Cuál es la precisión relativa de las mediciones en y ? Es posible que necesite escalar para ajustar esto. ¿Son independientes todas las medidas? La correlación puede complicar las cosas. b1b2
Brian Borchers, el

En este momento todo está modelado, así que sé y perfectamente, pero en la práctica sabré con tal vez 10 veces más precisión. Sin embargo, en este paso quiero asumir que los conozco a ambos por igual y que son independientes. b1b2b1
abnowack

Entonces, ¿cuál es tu pregunta aquí? Es fácil resolver el problema de los mínimos cuadrados de tres términos que ha dado en su pregunta.
Brian Borchers, el

¿Lo es? Si explica en una respuesta, lo marcaré como correcto. Solo uso rutinas básicas como el solucionador menos cuadrado de numpy . No soy de un fondo CS, por lo que podría estar perdiendo algo obvio.
abnowack

Respuestas:


9

Puedes escribir tu problema como

minFmg22

dónde

F=[A1A2αI]

y

g=[b1b20].

Puede usar cualquier solucionador lineal de mínimos cuadrados que desee para resolver este problema.

Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.