WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Gewöhnliche DGL divergiert

 

Ondy
Forum-Newbie

Forum-Newbie


Beiträge: 6
Anmeldedatum: 11.02.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 08.01.2014, 17:21     Titel: Gewöhnliche DGL divergiert
  Antworten mit Zitat      
Hallo Leute,

Habe ein Problem mit der Simulation eines Solarkollektors, bei dem Globalstrahlungskonstante, Windgeschwindigkeit, FLuideintrittstemperatur, Umgebungstemperatur und Luftfeuchte initiert werden und daraus der Wärmenutzenstrom berechnet werden soll. Dabei divergiert die Lösung gegen Ende und ich bekomme unbrauchbare Werte.



[code]
y0=[t_abs;t_abd;t_wd1;t_wd2;t_wd3];
[~,Y]=ode23s(@fdtdtau3,[Tau;Tau2],y0,[],C_abd,C_abs,C_wd,Q_zu,Q_absabd,Q_abdamb,Q_abswd,Q_absfluid,Q_wdamb,Q_wd1,Q_wd2);
t_abs=Y(end,1);
t_abd=Y(end,2);
t_wd1=Y(end,3);
t_wd2=Y(end,4);
t_wd3=Y(end,5);
[/code]

wobei:

[code]
function dy=fdtdtau3(~,y,C_abd,C_abs,C_wd,Q_zu,Q_absabd,Q_abdamb,Q_abswd,Q_absfluid,Q_wdamb,Q_wd1,Q_wd2)
dy=zeros(5,1);
dy(1)=(Q_zu-Q_absabd-Q_absfluid-Q_abswd)/C_abs;
dy(2)=(Q_absabd-Q_abdamb)/C_abd;
dy(3)=(Q_abswd-Q_wd1)/(C_wd/3);
dy(4)=(Q_wd1-Q_wd2)/(C_wd/3);
dy(5)=(Q_wd2-Q_wdamb)/(C_wd/3);
end
[/code]

Habe alle Mögliche Solver ausprobiert (ode23s, ode23t,ode23tb,..), doch alle divergieren ab einer gewissen Stelle im Programm.
Ob Divergenz auftritt, ist eigentlich nur von den gewählten Initialwerten abhängig.

t_abs (=dy(1) in fdtdtau3) scheint hauptverantwortlich dafür zu sein und "schwingt" nicht immer so ein wie es sollte...

Kann Mir bitte jemand helfen? Falls noch weitere Informationen benötigt werden, einfach sagen. Denke es handelt sich nicht um ein zu komplexen Sachverhalt, aber ich weiß einfach nicht weiter...

Grüße
Ondi
Private Nachricht senden Benutzer-Profile anzeigen


Thomas84
Forum-Meister

Forum-Meister


Beiträge: 546
Anmeldedatum: 10.02.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 09.01.2014, 13:48     Titel:
  Antworten mit Zitat      
In der Funktion fdttau3 taucht gar kein y auf. Ich vermute daher das die DGL noch nicht richtig ist. Was sind denn C_abd, etc? Ich vermute diese Werte ändern sich zeitlich?

viele Grüße
Thomas
Private Nachricht senden Benutzer-Profile anzeigen
 
Ondy
Themenstarter

Forum-Newbie

Forum-Newbie


Beiträge: 6
Anmeldedatum: 11.02.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 09.01.2014, 15:00     Titel:
  Antworten mit Zitat      
Ja das mit dem y hat mich auch schon stutzig gemacht, hab ich vom Vorgängermodell einfach so übernommen ehrlich gesagt, weils funktioniert hat. Wenn ich es weglasse, meldet er:

Too many input arguments

Wenn etwas anderes als y dasteht, läuft das Skript wenigstens...
denke ich hab irgendwo nen zentralen denkfehler, kann des sein?

C_abd,C_abs,C_wd sind Wärmekapazitäten und somit zeitlich konstant. Q_zu ist ebenso zeitlich konstant...

Viele Grüße!
Private Nachricht senden Benutzer-Profile anzeigen
 
Ondy
Themenstarter

Forum-Newbie

Forum-Newbie


Beiträge: 6
Anmeldedatum: 11.02.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 10.01.2014, 19:54     Titel:
  Antworten mit Zitat      
Hallo Leute,
Hab den Fehler selbst gefunden. Lag an einer Konstante, die ich nicht entsprechend angepasst habe!
Danke für die Mühen!

Grüße
Ondy
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

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 | goMatlab RSS Button 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.