For diffusions, a well-developed approach in rare event estimation is to introduce a suitable factorization of the reach probability and then to estimate these factors through simulation of an Interacting Particle System (IPS). This paper studies IPS based reach probability estimation for General Stochastic Hybrid Systems (GSHS). The continuous-time executions of a GSHS evolve in a hybrid state space under influence of combinations of diffusions, spontaneous jumps and forced jumps. In applying IPS to a GSHS, simulation of the GSHS execution plays a central role. From literature, two basic approaches in simulating GSHS execution are known. One approach is direct simulation of a GSHS execution. An alternative is to first transform the spontaneous jumps of a GSHS to forced transitions, and then to simulate executions of this transformed version. This paper will show that the latter transformation yields an extra Markov state component that should be treated as being unobservable for the IPS process. To formally make this state component unobservable for IPS, this paper also develops an enriched GSHS transformation prior to transforming spontaneous jumps to forced jumps. The expected improvements in IPS reach probability estimation are also illustrated through simulation results for a simple GSHS example.

This paper focuses on estimating reach probability of a closed unsafe set by a stochastic process. A well-developed approach is to make use of multi-level MC simulation, which consists of encapsulating the unsafe set by a sequence of increasing closed sets and conducting a sequence of MC simulations to estimate the reach probability of each inner set from the previous set. An essential step is to copy (split) particles that have reached the next level (inner set) prior to conducting a MC simulation to the next level. The aim of this paper is to prove that the variance of the multi-level MC estimated reach probability under fixed assignment splitting is smaller or equal than under random assignment splitting methods. The approaches are illustrated for a geometric Brownian motion example.

This paper studies estimation of reach probability for a generalized stochastic hybrid system (GSHS). For diffusion processes a well-developed approach in reach probability estimation is to introduce a suitable factorization of the reach probability and then to estimate these factors through simulation of an Interacting Particle System (IPS). The theory of this IPS approach has been extended to arbitrary strong Markov processes, which includes GSHS executions. Because Monte Carlo simulation of GSHS particles involves sampling of Brownian motion as well as sampling of random discontinuities, the practical elaboration of the IPS approach for GSHS is not straightforward. The aim of this paper is to elaborate the IPS approach for GSHS by using complementary Monte Carlo sampling techniques. For a simple GSHS example, it is shown that and why the specific technique selected for sampling discontinuities can have a major influence on the effectiveness of IPS in reach probability estimation.