Kleinste absolute Abweichungen, die mit dem Barrodale-Roberts-Algorithmus gelöst werden: Vorzeitige Beendigung?

9

Bitte entschuldigen Sie die längere Frage, es bedarf nur einer Erklärung, um zum eigentlichen Problem zu gelangen. Diejenigen, die mit den genannten Algorithmen vertraut sind, könnten wahrscheinlich direkt zum ersten Simplex-Tablau springen.


Um Probleme mit der geringsten absoluten Abweichung (auch bekannt als -Optimierung) zu lösen , ist der Barrodale-Roberts-Algorithmus eine spezielle Simplex-Methode, die weitaus weniger Speicher- und Rechenaufwand benötigt, um ein geeignetes Minimum zu finden.L1

Meine Implementierung des Algorithmus endet an einem einfachen Beispiel, bevor ein angemessenes Minimum erreicht ist. Lassen Sie mich das Problem jedoch wahrscheinlich zuerst ausführlicher darlegen:

Bei gegebenen Daten versucht die L 1 -Optimierung, c m zu finden , das n i = 1 | minimiert y i - f ( x i ) |(xi,yi)L1cm wobei A x ist eine n × m - Matrixdie in gewisser Weise davon abhängt , x . Dieses Problem kann als lineares Programm angegeben und daher unter anderem mit simplexartigen Methoden gelöst werden.

i=1n|yif(xi)|withf(x):=Axϕ
Axn×mx

L1rank(A)

Lei und Anderson schlugen 2002 eine kleine Modifikation vor , die die numerische Stabilität erhöhen und daher bekannte Probleme mit dem Simplex-Algorithmus überwinden soll.

Grundsätzlich geht dieser Algorithmus davon aus, dass Sie mit einem bestimmten Satz von Punkten beginnen, die interpoliert werden müssen, die angegebenen Prozeduren verwenden, um ein Simplex-Tableau zu erstellen, und dann anhand der Regeln von Barrodale und Roberts entscheiden, welche Basisvariablen geändert und daher modifiziert werden sollen Satz von Datenpunkten, der angenähert wird.

{(1,1),(2,1),(3,2),(4,3),(5,2)}a1+a2x

BasisRu1u3b11/23/21/2v21/21/21/2b21/21/21/2u41/21/23/2v5112Marginal cost210

2

Da alle nichtbasischen Vektoren nicht positive Grenzkosten haben [...]

Die Iteration ist beendet und das Optimum ist erreicht.

{2,5}

BasisRu2u5u11/34/31/3b11/35/32/3u32/32/31/3u44/31/32/3b21/31/31/3Marginal cost7/310/35/3

u2u1

Zusätzliche Informationen: Wenn ich mit dem von Tablrodale und Roberts angegebenen anfänglichen Tableau beginne, kann ich das obige Tableau auch mit gewöhnlichen Simplex-Schritten reproduzieren. Daher bin ich ziemlich sicher, dass die tatsächlichen numerischen Werte korrekt sind und meine Interpretation der Pivot-Auswahlregel ist fehlerhaft.

Irgendwelche Gedanken dazu?

Mir ist klar, dass die Frage selbst ziemlich kompliziert ist und wahrscheinlich Kenntnisse über mindestens den Barrodale- und Roberts-Algorithmus erfordert, um ausreichend beantwortet zu werden. Der Algorithmus als Ganzes ist zu lang, um ihn hier im Detail zu wiederholen. Wenn Sie jedoch weitere Fragen zu den von mir unternommenen Schritten oder zu fehlenden Informationen haben, können Sie diese gerne stellen, und ich werde die Frage gerne ergänzen.

Thilo
quelle
Wenn jemand mit ausreichendem Ruf ein Tag nach dem Vorbild von "Least-Absolute-Deviations" oder "L1-Approximation" erstellen könnte, wäre ich dankbar.
Thilo
Die Optimalitätsbedingung ist, dass die Basislösung machbar sein muss (in Bezug auf ihre Nicht-Negativitätsbeschränkungen) und dass die reduzierten Kosten kleiner oder gleich 0 sein müssen. Wenn Ihre Basislösung nicht realisierbar ist, sind alle Wetten ungültig.
Brian Borchers
Die Grundlösung ist konstruktionsbedingt realisierbar. Somit sollte es kein Problem geben. Ich habe jedoch eine erste Idee, wo das Problem liegen könnte. Ich werde eine entsprechende Antwort hinzufügen, wenn ich recht habe.
Thilo

Antworten:

4

Ich habe es gelöst. Eigentlich haben Barrodale und Roberts es gelöst und ich habe es einfach nicht sorgfältig gelesen.

uiiui=0vi

bjcjuivi

Daher muss mein Simplex-Tableau oben wie folgt aussehen:

BasisRu2u5v2v5u11/34/31/34/31/3b11/35/32/35/32/3u32/32/31/32/31/3u44/31/32/31/32/3b21/31/31/31/31/3Marginal cost7/310/35/34/31/3

v2

Vielen Dank, dass Sie gelesen haben und mir einen Platz zum Aufschreiben meines Problems gegeben haben, was normalerweise dazu beiträgt, die Lösung erheblich einzugrenzen. Hoffentlich kann diese Antwort manchmal für andere hilfreich sein, die versuchen, Barrodale & Roberts zu implementieren.

Thilo
quelle