Verfasst am: 21.08.2015, 18:09
Titel: Gleichungssystem abspeichern und iterieren
Hallo,
ich versuche momentan ein Gleichungssystem über mehrere Perioden zu lösen. Was ich am Ende raushaben will sollte ungefähr so aussehen:
1. Ich setzte die Anfangswerte von 4 Zustandsvariablen: Omeg0, E0, uM0, uI0
2. Die Zustandsvariablen lösen ein Gleichungssystem in einer m-Funktion nach 3 Variablen auf: lambda, thetaM, thetaI
3. Mit den Variablen definiere ich die Werte der Zustandsvariablen in der Folgeperiode, also vereinfacht Omeg1=lambd*Omeg0
4. Ich speichere alle neuen Zustandsvariablen als Startwerte und starte den Prozess von vorne, also z.B. Omeg1 = Omeg0
Leider weiß ich im Moment noch nicht, wie ich den Vektor der Zustandsvariablen mit der m-Funktion verbinde und wie ich die Resultate am besten abspeichere. Anhand von verschiedenen Tutorials habe ich bisher diesen Code geschrieben (3 Dateien, ein Mainfile mit der Iteration, 1 Funktionsdatei und eine Datei mit Parametern die ich für die beiden anderen verwende).
Verfasst am: 22.08.2015, 23:05
Titel: Re: Gleichungssystem abspeichern und iterieren
Hallo ouou,
Das ist ein großes Stück Code und auf Anhieb sehe ich nicht, zu welchem Teil genau diese Aussage passt:
Zitat:
Leider weiß ich im Moment noch nicht, wie ich den Vektor der Zustandsvariablen mit der m-Funktion verbinde und wie ich die Resultate am besten abspeichere.
Ich kann mir unter "Vektor mit M-File verbinden" nichts vorstellen.
Wie Du die Resultate "am besten" abspeicherst, hängt von Deinen Wünschen ab. Was wäre denn da wichtig?
Ich rate Dir, zunächst viel konkretere Fragen zu formulieren.
vielen Dank für deine Antwort. Ich weiß das es sehr viel Code ist, ich wollte vollständigkeitshalber alles posten, das File mit den Parametern und der Teil mit den Hilfsgleichungen kann getrost ignoriert werden.
Mit
Zitat:
Leider weiß ich im Moment noch nicht, wie ich den Vektor der Zustandsvariablen mit der m-Funktion verbinde und wie ich die Resultate am besten abspeichere.
meine ich konkret zwei Sachen:
1. Ich möchte folgende Startwerte für die Funktion equilbiriumfun.m einfügen:
Omeg0 = 0.15;
E0 = 0.007;
uM0 = 0.4;
uI0 = 0.393;
2. Nach der ersten Iteration sollen diese überschrieben werden mit den errechneten Werten Omeg1, E1, uM1 und uI1.
Daraus folgen die 2 Fragen:
1. Wie kann ich die Funktion so umschreiben, dass sie sich die Startwerte aus dem Mainfile holt?
2. Wie kann ich dem Mainfile sagen, dass es die Startwerte nach der ersten Berechnung mit den neuen Werten ersetzt?
1. Genau, ist eigentlich das selbe Problem, allerdings dachte ich, dass ich y0 jetzt als anonyme Funktion festgesetzt habe, so wie du es dort vorgeschlagen hattest. Im Beispiel in der Dokumentation http://de.mathworks.com/help/optim/.....parameters.html#brhkghv-9 sieht es so aus, als könnte die Funktion die Sachen selbstständig aus dem Masterfile extrahieren, aber d.h. ich muss das noch extra festlegen? In der Dokumentation zu http://de.mathworks.com/help/matlab.....cripts-and-functions.html sieht es so aus als könnte ich es mit GLOBAL machen, da kommt bei mir aber auch nur wieder "Undefined function or variable 'uI0'." raus...
Und es ist kein Problem für die späteren Iterationen, dass ich am Anfang des Skripts sowas wie "Omeg0=0.15" stehen habe?
hast du meien Vorschläge umzusetzen versucht? Mit welchem Ergebnis?
Zitat:
y0 jetzt als anonyme Funktion festgesetzt habe
y0 ist ein ganz normaler Vektor.
Zitat:
sieht es so aus, als könnte die Funktion die Sachen selbstständig aus dem Masterfile extrahieren, aber d.h. ich muss das noch extra festlegen?
Ich habe keine Ahnung, was du meinst. Man muss lediglich darauf achten, was die festen und was die variablen Parameter sein sollen. Anhand des Beispiels in der Doku und deines eigenen funktioniernden Beispiels sollte das doch nun wirklich kein Problem sein?
Zitat:
sieht es so aus als könnte ich es mit GLOBAL machen, da kommt bei mir aber auch nur wieder "Undefined function or variable 'uI0'." raus
Man kann es mit global machen, es ist aber unübersichtlich, schwierig zu debuggen, und daher nicht empfehlenswert. Ohne zu wissen, was du nun wo mit global machst, ist es nicht möglich, einen Korrekturvorschlag zu machen.
Zitat:
Und es ist kein Problem für die späteren Iterationen, dass ich am Anfang des Skripts sowas wie "Omeg0=0.15" stehen habe?
Nein, weil du das außerhalb der for-Schleife und damit nur einmalig setzt. Wenn du den Debugger nutzt, kannst du dich auch wunderbar selbst davon überzeugen.
vielen Dank für deine Hilfe, es tut mir Leid, wenn meine Fragen überflüssig erscheinen oder ich den Anschein erwecken, dass ich mich nicht ausreichend mit deinen Antworten befasse, ich versuche bestmöglich deine Tipps umzusetzen. Leider habe ich wirklich nicht viel Basiswissen, ich habe deshalb auch gerade erst verstanden, was du im anderen Thread meintest. Ich dachte die "anonymous function" wäre schon ausreichend mit einem Vektor im Mainfile definiert (hatte damals direkt das mit der Struktur ausprobiert und die anonymous function nicht weiter verfolgt, deswegen habe ich kein funktionierendes Beispiel).
Nach erneuter Lektüre der Dokumentation, habe ich jetzt den Teil mit dem normalen Vektor y0 = [...] ersetzt durch:
Ist dass ungefähr, was du meintest? Im Moment bekomme ich die Fehlermeldung:
Code:
Warning: Struct field assignment overwrites a value with class "double". See MATLAB R14SP2 Release Notes, Assigning
Nonstructure Variables As Structures Displays Warning, for details.
> In createOptionFeedback (line33)
In prepareOptionsForSolver (line31)
In fsolve(line141)
In Equilibrium (line19) Error using @(x)equilibriumfun(x,y0)
Too many input arguments.
Der Startvektor x0 muss transponiert werden, weil du in equilibriumfun von einem Spaltenvektor ausgehst.
Jetzt gibt es noch eine Fehlermeldung wegen mu, da es als Funktionsaufruf interpretiert wird, obwohl es als Variable gedacht ist. Man sollte also mu durch einen anderen Variablennamen (z.B. mue) ersetzen.
Wie angehängt läuft der Code durch. Bitte die Änderungen nachvollziehen und überprüfen.
Schließlich gibt es im Editor noch den Code Analyzer, der per roter / oranger Markierung Verbesserungsvorschläge im Code macht. Diese solltest du dir noch ansehen.
vielen Dank für deine Hilfe und deine Gedulg, ich hatte ursprünglich etwas ganz ähnliches probiert, aber es hat nicht auf Anhieb funktioniert und ich dachte dann, dass es zu leicht wäre. Falls ich eine Dankeswidmung für meine Masterarbeit schreibe, kriegst du auf alle Fälle einen Platz
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