|
|
| papteo |
Forum-Newbie
|
 |
Beiträge: 2
|
|
|
|
Anmeldedatum: 29.10.12
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 29.10.2012, 18:24
Titel: Hilfe mit ODE
|
 |
Hallo,
ich will zwei kombinierte differentielle Gleichung mit ode45 zeichen. Ich will das nur für die Temperatur 1300 ° C eine Gleichung gilt und für die restliche Temperatur gilt es die andere Gleichung. ich habe diese folgende Code geschrieben, aber die IF function funktioniert nicht. Es sollte bei T=1300 die Temperature kurz (0,2 s) konstant bleiben. Haben Sie Vorchlagen?
[code]
function dTdt=f(t,T)
a=-0.96*10^(-10);
b=-1.78;
c=522.27;
d=22876.7;
if (T==1300)
dTdt=a*T^4+b*T+c+d*(exp(-2643.68/T)-exp(-4165.83/T));
else
dTdt=a*T^4+b*T+c;
end
>> clear all
>> time=[0.01 1];
>> iT=1700;
>> [t,T]=ode45(@f,time,iT)
>> plot(t,T)
Vielen Dank im Voraus!! 
|
|
|
|
|
|
| Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
|
Verfasst am: 29.10.2012, 19:52
Titel:
|
|
Hallo,
es ist nicht sinnvoll, die Funktion nur zu einem exakten Wert anders zu definieren. Es ist schließlich nicht davon auszugehe, dass bei einer Simulation mit diskreten Werten (wie von ode45 durchgeführt) überhaupt genau dieser Wert getroffen wird.
Sinnvoller wäre es, wenn eine Bedingung der Form
oder vergleichbares gewählt wird.
Grüße,
Harald
|
|
|
|
| papteo |
Themenstarter
Forum-Newbie
|
|
Beiträge: 2
|
|
|
|
Anmeldedatum: 29.10.12
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 29.10.2012, 20:06
Titel:
|
|
Hallo Harald,
Vielen Dank für deine Antwort. Das habe ich auch schon probiert aber wieder ohne Erfolg. Da ich nicht so viele Erfahrung mit Matlab habe, wollte ich auch fragen, ob meine Code richtig geschrieben ist. Wahrscheinlich liegt dort das Problem.
Vielen Dank und schönen Abend noch!
|
|
|
|
| Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
|
Verfasst am: 29.10.2012, 20:27
Titel:
|
|
Hallo,
der Code an sich sieht gut aus. Die Bedingung greift nur eben nicht, wie oben erklärt.
Grüße,
Harald
|
|
|
|
|
|
|
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 - 2024
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.
|
|
