论文信息 - Solving 0 = F(t, y(t), y′(t)) in Matlab

Solving 0 = F(t, y(t), y′(t)) in Matlab

Abstract Important algorithms and design decisions of a new code, ode15i, for solving 0 = F(t, y(t), y′(t)) are presented. They were developed to exploit Matlab, a popular problem solving environment (PSE). Codes widely-used in general scientific computation (GSC) compute consistent initial conditions only for restricted forms of the differential equations. ode15i does this for the general problem, a task that is qualitatively different. Unlike popular codes in GSC, ode15i saves partial derivatives because this is more efficient in the PSE. A new representation of fixed leading coefficient BDFs is based on the Lagrangian form of polynomial interpolation. Basic computations are both clear and efficient in the PSE when using this form. Some numerical experiments show that ode15i is an effective way to solve fully implicit ODEs and DAEs of index 1 in Matlab.

Lawrence F. Shampine | L. Shampine

[1] Kevin Burrage,et al. An implementation of singly-implicit Runge-Kutta methods , 1980 .

[2] Linda R. Petzold,et al. Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

[3] H. Keller,et al. Analysis of Numerical Methods , 1967 .

[4] C. Kelley. Iterative Methods for Linear and Nonlinear Equations , 1987 .

[5] Linda R. Petzold,et al. Consistent Initial Condition Calculation for Differential-Algebraic Systems , 1998, SIAM J. Sci. Comput..

[6] L. Shampine,et al. Numerical Solution of Ordinary Differential Equations. , 1995 .

[7] E. Hairer,et al. Solving Ordinary Differential Equations I , 1987 .

[8] Ralph E. White,et al. An initialization subroutine for DAEs solvers: DAEIS , 2001 .

[9] J. M. Watt. Numerical Initial Value Problems in Ordinary Differential Equations , 1972 .

[10] Ron Sacks-Davis,et al. An Alternative Implementation of Variable Step-Size Multistep Formulas for Stiff ODEs , 1980, TOMS.

[11] Stig Skelboe,et al. The control of order and steplength for backward differentiation methods , 1977 .

[12] R. Brayton,et al. A new efficient algorithm for solving differential-algebraic systems using implicit backward differentiation formulas , 1972 .

[13] M. Powell,et al. On the Estimation of Sparse Jacobian Matrices , 1974 .

[14] D. Arnold,et al. Computer Solution of Ordinary Differential Equations. , 1981 .