Ich habe folgendes Optimierungsproblem:
Gesucht wird der Winkel alpha, welcher von theta (gegeben) und p und q abhängt.
alpha = p*theta + q
Wie eine Optimierung in Matlab mit einem Parameter funktioniert, weiss ich. Doch nun habe ich ja zwei Parameter. Hier bin ich irgendwie überfordert. Ich benutze in Matlab lsqnonlin für die Optimierung und habe auch schon einige m-files geschrieben, doch irgendwo ist der Wurm drin. Kann sonst die m-files auch mal hier posten, falls das weiter helfen könnte.
lsqnonlin kann genausogut Funktionen von mehreren Variablen (bei Dir: Parametern) verarbeiten. Zum Beispiel ist dann das zweite Input-Argument x0 ein Vektor. Du müsstest also p und q zu einem Vektor zusammenfassen.
Wenn Dir das noch nicht hilft, poste am besten wirklich den (lauffähigen) Code, den Du für eine Variable verwendet hast.
Danke für die Antwort. Das Ganze ist leider nur ein Teil einer grösseren Optimierungsaufgabe. Ich habe mal alles zusammen geschrieben. Doch beim Ausführen tritt der folgende Fehler auf:
??? Error using ==> lsqncommon at 101
LSQNONLIN cannot continue because user supplied objective function failed with the following error:
Error using ==> mrdivide
Matrix dimensions must agree.
Error in ==> lsqnonlin at 182
[x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...
Error in ==> optimise_myfun at 24
[x, resnorm] = lsqnonlin(@myfun,x0,xl,xu,options)
Woran kann das liegen? Welche Matrize stimmt nicht?
Ich poste hier mal meine Codes und die verwendeten Daten sind im Excel-Dokument.
d1_zf ist eine Zahl, also hier unproblematisch. alpha hat die Größe 1045x1 und l1_zf die Größe 1x256. Deshalb gibt es Probleme beim /-Operator, der auch mrdivide heißt (siehe Fehlermeldung).
Wenn die Vektoren alpha und l1_zf gleich lang wären und Du eine elemntweise Division wolltest, ginge das so:
also mit einem Punkt vor /. Da die Vektoren aber gar nicht gleich lang sind, ist mir völlig unklar, was da genau geteilt werden soll.
Ich hab Dein myfun zum Testen übrigens folgendermaßen aufgedröselt. So kann man erstens sehen, was eigentlich gerechnet wird, und andererseits werden die häufig verwendeten Terme auch nur einmal berechnet, was etwas Zeit spart.
Danke für die Antwort. Ich habe später dann auch festgestellt, dass die Grössen der einzelnen Vektoren nicht übereinstimmen. Ich habe es jetzt geändert und es funktioniert. Und deine Aufteilung der einzelnen Terme ist gut. Werde ich noch so verwenden. Danke für die Hilfe
Gruss Alexa
Einstellungen und Berechtigungen
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
RSS