Download All

Code covered by the BSD License  

Highlights from
Manuscript of 'Solving Index-1 DAEs in MATLAB and Simulink'

  • ex1 This is a singular perturbation problem involving Michaelis-Menten
  • ex1 This is a singular perturbation problem involving Michaelis-Menten
  • ex10 This is the DAS 2 test problem of I. Cameron, Solution of differential-
  • ex10 This is the DAS 2 test problem of I. Cameron, Solution of differential-
  • ex10ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
  • ex11 This is a mechanics problem studied in section 4 of Chapter XII of
  • ex11 This is a mechanics problem studied in section 4 of Chapter XII of
  • ex11ode(t,u,flag) Copyright 1999, The MathWorks, Inc.
  • ex12
  • ex12
  • ex12ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
  • ex13 This is the discharge pressure control problem from E. Hairer, Ch. Lubich,
  • ex13 This is the discharge pressure control problem from E. Hairer, Ch. Lubich,
  • ex13ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
  • ex14 DAE examples of index greater than 1
  • ex14 DAE examples of index greater than 1
  • ex14ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
  • ex1ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
  • ex2 This is an artificial problem with two sets of consistent initial
  • ex2 This is an artificial problem with two sets of consistent initial
  • ex2ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
  • ex3 This is the pendulum problem posed as an index 1 problem. Formulated as
  • ex3 This is the pendulum problem posed as an index 1 problem. Formulated as
  • ex3ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
  • ex4 The Robertson problem is a classic test problem for codes that solve
  • ex4 The Robertson problem is a classic test problem for codes that solve
  • ex4ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
  • ex5 This is the two phase plug flow problem from E. Hairer, Ch. Lubich,
  • ex5 This is the two phase plug flow problem from E. Hairer, Ch. Lubich,
  • ex5ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
  • ex6 This is the two phase plug flow problem from E. Hairer, Ch. Lubich,
  • ex6 This is the two phase plug flow problem from E. Hairer, Ch. Lubich,
  • ex6ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
  • ex7 This is the one transistor amplifier problem displayed on p. 377 of
  • ex7 This is the one transistor amplifier problem displayed on p. 377 of
  • ex7ode(t,u,flag) Copyright 1999, The MathWorks, Inc.
  • ex8 This is the two transistor amplifier problem of E. Hairer, C. Lubich,
  • ex8 This is the two transistor amplifier problem of E. Hairer, C. Lubich,
  • ex8ode(t,u,flag) Copyright 1999, The MathWorks, Inc.
  • ex9 This is the DAS 1 test problem of I. Cameron, Solution of differential-
  • ex9 This is the DAS 1 test problem of I. Cameron, Solution of differential-
  • ex9ode(t,y,flag) Copyright 1999, The MathWorks, Inc.
  • ex1.m
  • ex10.m
  • ex11.m
  • ex12.m
  • ex13.m
  • ex14.m
  • ex2.m
  • ex3.m
  • ex4.m
  • ex5.m
  • ex6.m
  • ex7.m
  • ex8.m
  • ex9.m
  • View all files

3.0

3.0 | 2 ratings Rate this file 25 Downloads (last 30 days) File Size: 308.17 KB File ID: #7481

Manuscript of 'Solving Index-1 DAEs in MATLAB and Simulink'

by Jacek Kierzenka

 

19 Apr 2005 (Updated 12 Oct 2010)

Manuscript of the SIAM Review paper + M-files for all examples.

| Watch this File

File Information
Description

This is a manuscript of the paper

Lawrence F. Shampine, Mark W. Reichelt, and Jacek A. Kierzenka,

Solving Index-1 DAEs in MATLAB and Simulink,
SIAM Review, Vol. 18, No. 3, pp. 538-552, 1999.

The paper describes mathematical and software developments needed for effective solution of differential algebraic equations of index 1. Modifications introduced to ODE15S and ODE23T in MATLAB 5.3 allow for solving linearly-implicit DAEs of index 1.

The M-files for all the original examples described in the paper are available in the directory DAE_examples.

The directory DAE_examples_65 contains the files for the examples updated to take advantage of features available in MATLAB 6.5 (R13). Correspondingly, DAE_examples_70 contains files modified for use with MATLAB 7.0 (R14).

NOTE:
The paper describes solving linearly-implicit DAEs with time-dependent, singular mass matrix
    M(t)*y' = f(t,y)

In MATLAB 6.0, the solvers were extended to handle DAEs with state-dependent mass matrices
    M(t,y)*y'=f(t,y)
MATLAB release MATLAB 5.3 (R11)
Tags for This File  
Everyone's Tags
Tags I've Applied
Add New Tags Please login to tag files.
Comments and Ratings (3)
08 Aug 2005 George Chamoun

the examples only described systems of the form M(t)*y' = f(t,y).

there were no state-dependent mass matrix examples of the form M(t,y)*y'=f(t,y)
21 Nov 2005 valialah zareei

Dear sir
I am simulating packed tower steady-state
by jacobian but don converge and need to
dynamic simulation ,Please help me
  Best regard
27 Oct 2007 Luis Enrique Avendaño  
Please login to add a comment or rating.
Updates
12 Oct 2010

Copyright update
Tag Activity for this File
Tag Applied By Date/Time
differential equations Jacek Kierzenka 22 Oct 2008 07:46:21
dae Jacek Kierzenka 22 Oct 2008 07:46:21
mathematics Jacek Kierzenka 22 Oct 2008 07:46:21
differentialalgebraic equations Jacek Kierzenka 22 Oct 2008 07:46:21
differentialalgebraic equations Saurabh 11 Dec 2008 19:19:45

Contact us at files@mathworks.com