|
|
| goldloeckchenMD |
Gast
|
 |
Beiträge: ---
|
|
|
|
Anmeldedatum: ---
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 27.11.2009, 11:42
Titel: Runge-Kutta-Verfahren
|
 |
Hallo!
Ich hänge momentan ziemlich fest und hoffe, dass ihr mir helfen könnt. Wir sollen das Runge-Kutta-Verfahren Implementieren und haben eine DGL vorgegeben.
Ich mein bisher geschriebenes einfach 'reinkopiert. In der Schleife stehe ich auf dem Schlauch. f steht ja für die Funktion, also dxdt, nur orgendwie ist mir der Syntax unklar!
p.mue = 0.6396;
p.k_d = 0.0080;
tspan = [0:5];
laenge_int = 5;
n = 20;
delta_t = laenge_int/n; %Schrittweite
%solver-Aufruf
[t,x] = ode45(@wachtum,tspan,x0,[],p);
figure(2);
plot(t,x)
grid on;
xlabel('t in [1/h]')
ylabel('x in [kg*m^-3]')
%Runge-Kutta-Verfahren
x(1) = 0.05;
for i=2:length(t)
k1=delta_t*f(t(i-1),x(i-1));
k2=delta_t*f(t(i-1)+delta_t/2,x(i-1)+k1/2);
k3=delta_t*f(t(i-1)+deltat/2,x(i-1)+k2/2);
k4=delta_t*f(t(i-1)+delta_t,x(i-1)+k3);
x(i)=x(i-1) + delta_t*(k1+2*k2+2*k3+k4)/6
end
Ich hoffe, ihr könnt mir helfen!
function dxdt=wachtum(t,x,p)
dxdt_neu = p.mue*x-p.k_d*x;
dxdt = [dxdt_neu];
|
|
|
|
|
|
| BSYT'ler_MD |
Gast
|
|
Beiträge: ---
|
|
|
|
Anmeldedatum: ---
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 04.12.2009, 15:41
Titel:
|
|
Hi, Bioprozesstechnik bei Reichl, wa?^^
Das Runge-Kutta-Verfahren ist im ode45 implementiert und muss NICHT extra programmiert werden:
//Aufruf
disp('Runge-Kutta-Verfahren (RKV) mittels ode45:')
Start_b = [0.5]; %Startwert für RKV
Para_b = [mue_b,k_d_b]; %Enthält µ aus Regression und k_d laut Aufgabe
options = odeset;
[T_RKV,X_RKV] = ode45(@RKVb,[0 5],Start_b,options,Para_b);
//Funktion
function Rueckgabe = RKVb(t_span,X,Parameter)
Rueckgabe = Parameter(1)*X-Parameter(2)*X;
end
|
|
|
|
|
|
|
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.
|
|
