An efficient method for the numerical integration of measured variable dependent ordinary differential equations

• The proposed numerical integration method can efficiently integrate ODEs which incorporate frequently measured variables.
• The method only evaluates ODEs at time instance at which the measurements are available, such avoiding interpolation.
• For the integration of dynamic neural networks or hybrid models the method can significantly reduce the number of function evaluations as compared to fixed step-size methods.

The Ordinary Differential Equations (ODEs) of dynamic models that are used in process monitoring, control or optimization, are not only functions of time and states, but also of measured variables. So far two possibilities for the numerical integration of such ODEs were given: (i) a fixed step size integration schema could be applied, matching the step size to the time instances of the measurements; or (ii) using an adaptive step size method while interpolating the measurements. While fixed step size methods are computationally expensive, the repetitive interpolation of measurements for the application of adaptive step size methods is prone to errors and time prohibitive, especially for great numbers of measured variables.In this paper, an adaptive step size numerical integration method is proposed and evaluated for dynamic neural network/hybrid semi-parametric models. The method evaluates the ODEs only at time instances at which online measurements are available and adapts the step size according to those time instances. The numerical solution of the ODEs is provided at time instances which are specified by the user, i.e. at time instances of offline measured states. The rationale behind the proposed method is carefully analyzed, and it is demonstrated that its application along with a hybrid model/dynamic neural network model can result into a significant reduction of number of function evaluations, in the studied cases about 50%, while adhering user specified error tolerances for the numerical integration. In addition, the mutual interference between step-size adaption, parameter identification, coping of the neural network and model performance is studied, a fact that other studies have paid little to no attention.