Using integrals of the state transition matrix for efficient transient-response computations

An iterative technique for the computation of lower-order state systems with poles at the origin (such as aircraft lateral dynamics) or with time-polynomial inputs is developed and demonstrated. The integrated input or output signals are computed using the integrals of the state transition matrix, thus reducing the system order and computation time. In an application to the lateral response of a DC-8 aircraft (Grunwald, 1981), this method required only 7824 multiplication or addition operations (at fourth order) to achieve the accuracy obtained in 15,012 operations by a conventional sixth-order computation.