|
|
|
System Identification Toolbox: Parameterverlauf bei pem |
|
| EliteTUM |
Forum-Fortgeschrittener
|
 |
Beiträge: 70
|
|
|
|
Anmeldedatum: 21.04.11
|
|
|
|
Wohnort: München
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 09.02.2012, 12:02
Titel: System Identification Toolbox: Parameterverlauf bei pem
|
 |
| |
 |
 |
Hallo zusammen,
ich verwende pem zur Parameterschätzung verschiedener Dynamikmodelle. Die Ergebnisse sind recht zufriedenstellend.
Zur Gegenüberstellung mit Kalman-Filtern, welche ich ebenfalls implementiert habe, würde ich gerne den Parameterverlauf darstellen. D.h. wie die Schätzung der Parameter von Schätzschritt zu Schätzschritt sich entwickelt. Ich weiß, dass die entsprechende Toolbox in Simulink dies machen kann, jedoch steht mir diese hier nicht zur Verfügung.
Nun würde ich dies gerne mit dem pem-Befehl machen.
Gibt es irgendwie eine Möglichkeit, die Parameterverläufe zu erhalten? Auf die einzelnen Parameterwerte kann man am Ende ja mit
zugreifen. Meine einzige, aber zugegeben sehr unschöne Idee wäre:
die Anzahl maximaler Iterationen "MaxIter" in einer while-Schleife hochzählen lassen und jedesmal das Endergebnis für die Parameter abspeichern. Als Abbruchkriterium kann man ja entweder eine vorher definierte maximale Iterationszahl wählen und/oder einen minimalen Gradienten aller Parameter. Aber das ist denkbar höchst ineffizient. Andere Ideen?
_________________
- EliteTUM
_____________________________________
|
|
|
|
|
|
| User_Max |
Gast
|
|
Beiträge: ---
|
|
|
|
Anmeldedatum: ---
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 31.07.2012, 16:59
Titel: Parameteroptimierung mit PEM
|
|
Hey,
Ich beschäftige mich gerade mit dem PEM und dem GREYEST Befehlen in Matlab. Ein sehr einfaches dynamisches Modell (Zweimassenschwinger) soll angepasst werden. Das Modell kommt mit der System Identifikation Toolbox (wenn ich die Option State Space Modell [6x6] verwende) auf eine 100% Übereinstimmung, was zu erwarten war.
Jedoch komm ich mit dem PEM Befehl nur auf 78.31%.
Also irgendwas scheine ich falsch zu machen, kannst du mir da vll. weiterhelfen ?
Hier der Code:
data = frd(freqresp(sys_init(1,1),w),w,0);
m_idgrey = idgrey('KfE21_return_ss',par,'c',aux,T);
m_est = pem(data,m_idgrey);
compare(data,m_est)
|
|
|
|
| EliteTUM |
Themenstarter
Forum-Fortgeschrittener
|
|
Beiträge: 70
|
|
|
|
Anmeldedatum: 21.04.11
|
|
|
|
Wohnort: München
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 05.08.2012, 19:59
Titel:
|
|
ich würde mal darauf tippen, dass unterschiedliche Standard-Werte für die Schätzalgorithmen verwendet werden.
Lass dir ausgeben, weshalb der Befehl PEM beendet wurde. Das lässt sich im Struct des Typs IDGREY, also bei dir der Variable "m_est" herausfinden (siehe HELP PEM oder HELP IDGREY). Evtl. solltest du einen anderen Algorithmus (z.B. Levenberg-Marquardt) oder niedrigere Toleranzen für Fehler verwenden. Evtl. auch mehr Iterationsschritte?!?
_________________
- EliteTUM
_____________________________________
|
|
|
|
|
|
|
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
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
|  RSS
Hosted by:
Copyright © 2007 - 2025
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
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.
|
|
