Lösen von zwei inversen Problemen mit derselben Lösung

8

Ich habe zwei umgekehrte Probleme,

A1 x=b1A2 x=b2

Bisher habe ich sie unabhängig voneinander mithilfe der Tikhonov-Regularisierung gelöst und zwei Schätzungen für x . In meinem Fall repräsentiert x in beiden Gleichungen die gleiche Lösung. Ist es möglich, eine "simultane" Lösung durchzuführen? Im Idealfall würde ich die Antwort finden

min(A1xb12+A2xb22+Γx2)

Wobei und I die Identitätsmatrix wie bei der Tikhonov-Regularisierung (auch bekannt als Ridge-Regression) ist. Ich nehme an, ich könnte einfach den Durchschnitt beider Lösungen nehmen und mich fragen, ob es jedoch einen statistisch aussagekräftigeren Weg gibt, dies zu erreichen.Γ=α II

abnowack
quelle
1
Wie ist die relative Genauigkeit der Messungen in und ? Möglicherweise müssen Sie skalieren, um dies anzupassen. Sind alle Messungen unabhängig? Korrelation kann Dinge komplizieren. b1b2
Brian Borchers
Im Moment ist alles modelliert, so dass ich und perfekt kenne , aber in der Praxis werde ich mit vielleicht 10x mehr Genauigkeit kennen. In diesem Schritt möchte ich jedoch davon ausgehen, dass ich beide gleichermaßen kenne und dass sie unabhängig sind. b1b2b1
abnowack
Was ist Ihre Frage hier? Es ist einfach, das Problem der kleinsten Quadrate mit drei Begriffen zu lösen, das Sie in Ihrer Frage angegeben haben.
Brian Borchers
Ist es? Wenn Sie in einer Antwort erklären, werde ich es richtig markieren. Ich benutze nur grundlegende Routinen wie den kleinsten quadratischen Löser von numpy . Ich habe keinen CS-Hintergrund, daher könnte mir etwas Offensichtliches fehlen.
abnowack

Antworten:

9

Sie können Ihr Problem als schreiben

minFmg22

wo

F=[A1A2αI]

und

g=[b1b20].

Sie können einen beliebigen linearen Löser für kleinste Quadrate verwenden, um dieses Problem zu lösen.

Brian Borchers
quelle