Beispiel - SAM - ETH Zürich

MATLAB-CODE : Lösung einer Index-1-DAE
f = @(t,y) [ ... ];
J = @(t,y) [ ... ];
M = @(t,y) [ ... ];
opts = odeset(’Mass’,M,’Jacobian’,J);
[t,y] = ode15s(f,tspan,y0,opts);
MATLAB-Code zur Lösung allgemeiner
DAEs
M(t,y)ẏ = f(t,y) .
Funktionf
Jacobi-MatrixD y f
Matrix(funktion)M(t,y)
Numerische
Mathemtik
Alternativer Integrator: ode23t (gleicher Aufruf)
(Beachte: Alle MATLAB DAE-Integratoren benutzen adaptive Schrittweitensteuerung, vgl. Sect. 2.6)
△
R. Hiptmair
rev 35327,
24. Juni
2011
Beispiel 3.8.15 (Lösung der Schaltkreis-DAEs mit MATLAB).
DAE (3.8.2) mitC = 1,R 2 = 1,R 1 = 1000,K = 10,I 0 = 10 −4
MATLAB-Integratoren ode23t, ode15s, default Toleranzen
3.8
p. 418