A review of numerical methods for digital simulation

This review is based on an extended study of the litera ture and experience in the detailed simulation of guided missiles and similar systems. After a brief recapitulation of the main ideas behind the classical methods for numer ically solving ordinary differential equations, the recent literature on this subject is surveyed and discussed. Pre dictor, predictor-corrector, and Runge-Kutta methods of various orders are then compared experimentally and theoretically. The theoretical comparison is based on the performance in computing the steady-state frequency response of linearized systems and includes consideration of accuracy, numerical stability, recovery of continuous outputs from sampled outputs, and computing time per step. The factors determining the optimum order are brought out. Other procedures such as different step lengths in different parts of the problem, special-purpose difference equations for "stiff" linear sections of the prob lem, and procedures for inserting wideband noise are dis cussed. The methods currently in use with a number of simulation languages are then summarized, and it is found that there is little consensus on "best" methods. It is concluded that further development is needed of meth ods for handling (a) stiff equations, including the non linear case, (b) mixed step lengths, (c) noise insertion, (d) discontinuities, and (e) recovery of continuous out puts. A list of 75 references is given.

[1]  Russell W. Stineman.  Digital Time-Domain Analysis of Systems with Widely Separated Poles , 1965, JACM.

[2]  Maury E. Fowler,et al.  A New Numerical Method for Simulation , 1965 .

[3]  M. Rubinoff,et al.  Numerical solution of differential equations , 1954, AIEE-IRE '54 (Eastern).

[4]  T. E. Hull,et al.  Corrector Formulas for Multi-Step Integration Methods , 1962 .

[5]  G. Dahlquist.  Convergence and stability in the numerical integration of ordinary differential equations , 1956 .

[6]  Morris Rubinoff,et al.  Digital Computers for Real-Time Simulation , 1955, JACM.

[7]  James D. Riley,et al.  Stability Properties of Adams-Moulton Type Methods , 1965 .

[8]  Fred T. Krogh.  Predictor-Corrector Methods of High Order With Improved Stability Characteristics , 1966, JACM.

[9]  Elmer G. Gilbert.  Some critical remarks on a new numerical method for simulation of dynamical systems , 1966 .

[10]  Richard Wesley Hamming,et al.  Stable Predictor-Corrector Methods for Ordinary Differential Equations , 1959, JACM.

[11]  Hans J. Stetter,et al.  Generalized Multistep Predictor-Corrector Methods , 1964, JACM.

[12]  R.W.H. Sargent,et al.  Dynamic behaviour of multi-component multi-stage systems. Numerical methods for the solution , 1962 .

[13]  A. Nordsieck,et al.  On Numerical Integration of Ordinary Differential Equations , 1953 .

[14]  William H. Anderson,et al.  A Numerical Method for Solving Control Differential Equations on Digital Computers , 1960, JACM.

[15]  Peter R. Benyon,et al.  State variable difference methods for digital simulation , 1967 .

[16]  Harry J. Gray.  Propagation of Truncation Errors in the Numerical Solution of Ordinary Differential Equations by Repeated Closures , 1955, JACM.

[17]  Roger A. Gaskill.  Fact and fallacy in digital simulation , 1965 .

[18]  Elmer G. Gilbert.  Dynamic-error analysis of digital and combined analog-digital computer systems , 1966 .

[19]  C. E. Treanor,et al.  A Method for the Numerical Integration of Coupled First-Order Differential Equations with Greatly Different Time Constants , 1966 .

[20]  W. E. Milne,et al.  Numerical Integration of Differential Equations. , 1957 .

[21]  C. V. D. Forrington.  Extensions of the Predictor-corrector Method for the Solution of Systems of ordinary Differential Equations , 1961, Comput. J..

[22]  Harry C. Shaw.  Discrete Analogs for Continuous Filters , 1966, JACM.

[23]  Emil Grosswald.  Transformations Useful in Numerical Integration Methods , 1959 .

[24]  D. W. Martin.  Runge-Kutta Methods for Integrating Differential Equations on High Speed Digital Computers , 1958, Comput. J..

[25]  T. E. Hull,et al.  Efficiency of Predictor-Corrector Procedures , 1963, JACM.

[26]  Lucio Tavernini.  On UNIP and the construction of digital simulation programs , 1966 .

[27]  Herbert M. Gurk.  The use of stability charts in the synthesis of numerical quadrature formulae , 1955 .

[28]  Roger L. Crane,et al.  Stability of a Generalized Corrector Formula , 1962, JACM.

[29]  Hans J. Stetter,et al.  Stabilizing predictors for weakly unstable correctors , 1965 .

[30]  P. E. Chase.  Stability Properties of Predictor-Corrector Methods for Ordinary Differential Equations , 1962, JACM.

[31]  C. W. Gear.  The numerical integration of ordinary differential equations , 1967 .

[32]  H. J. Gray,et al.  Numerical methods in digital real-time simulation , 1954 .

[33]  R. R. Reynolds,et al.  Fifth-Order Methods for the Numerical Solution of Ordinary Differential Equations , 1962, JACM.

[34]  John R. Rice,et al.  Split Runge-Kutta method for simultaneous equations , 1960 .

[35]  A. C. R. Newbery.  Multistep integration formulas , 1963 .

[36]  R. L. Crane,et al.  A Predictor-Corrector Algorithm with an Increased Range of Absolute Stability , 1965, JACM.

[37]  R. R. Reynolds,et al.  Stability of a Numerical Solution of Differential Equations , 1959, JACM.

[38]  S. Gill,et al.  A process for the step-by-step integration of differential equations in an automatic digital computing machine , 1951, Mathematical Proceedings of the Cambridge Philosophical Society.

[39]  John C. Butcher,et al.  On the Convergence of Numerical Solutions to Ordinary Differential Equations , 1966 .

[40]  David A. Pope.  An exponential method of numerical integration of ordinary differential equations , 1963, CACM.