Computing bounds for imprecise continuous-time Markov chains using normal cones

The theory of imprecise Markov chains has achieved significant progress in recent years. Its applicability, however, is still very much limited, due in large part to the lack of efficient computational methods for calculating higher-dimensional models. The high computational complexity shows itself especially in the calculation of the imprecise version of the Kolmogorov backward equation. The equation is represented at every point of an interval in the form of a minimization problem, solvable merely with linear programming techniques. Consequently, finding an exact solution on an entire interval is infeasible, whence approximation approaches have been developed. To achieve sufficient accuracy, in general, the linear programming optimization methods need to be used in a large number of time points. The principal goal of this paper is to provide a new, more efficient approach for solving the imprecise Kolmogorov backward equation. It is based on the Lipschitz continuity of the solutions of the equation with respect to time, causing the linear programming problems appearing in proximate points of the time interval to have similar optimal solutions. This property is exploited by utilizing the theory of normal cones of convex sets. The present article is primarily devoted to providing the theoretical basis for the novel technique, yet, the initial testing shows that in most cases it decisively outperforms the existing methods.