Determining Approximate Solutions of Nonlinear Ordinary Differential Equations Using Orthogonal Colliding Bodies Optimization

摘要

The solution of nonlinear Ordinary Differential Equations (ODE) finds potential applications in physics, economics, computing and engineering. Conventional approaches used for solving ODE are effective in case of 1st order or 2nd order problems. With increase in order the complexity associated with the problem increases. Thus instead of going for an exact solution, determining an approximate solution is also helpful. In this paper, solving ODE is handled as an optimization problem. Popular Fourier Series expansion is used as an approximation function. A hybrid algorithm Orthogonal Colliding Bodies Optimization (OCBO) is formulated by assembling good features of orthogonal array (exploration of solution in search space) and Colliding Bodies Optimization (bodies after collision quickly move to a position with minimal energy level on the search space). The coefficients of the Fourier series are computed with OCBO. Simulation studies are carried out to determine solution of popular practically used ODEs: Bernoulli Equation for flowing fluids, Integro-Differential equation, Brachistochrone equation for gravity, current response of an oscillatory Tank circuit, voltage and current decay with time in an electrical circuit. Simulations are also carried out on three benchmark ODEs used for modelling the biological processes. Comparative analysis demonstrated the superior approximation of the proposed approach over Orthogonal PSO, water cycle algorithm and Harmonic search.