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.
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 ) | 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.
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.
Da alle nichtbasischen Vektoren nicht positive Grenzkosten haben [...]
Die Iteration ist beendet und das Optimum ist erreicht.
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.
Antworten:
Ich habe es gelöst. Eigentlich haben Barrodale und Roberts es gelöst und ich habe es einfach nicht sorgfältig gelesen.
Daher muss mein Simplex-Tableau oben wie folgt aussehen:
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.
quelle