Let (Ω,F,P) be a probability space, let T>0 and Λ≥0 be real numbers, let E be a nonempty finite set (the state space), and let A be a nonempty finite set (the clock labels) with ∣A∣ elements. Write B[0,T] for the trace Borel σ-algebra on [0,T], integrals over compact intervals for the Lebesgue integral over the compact interval (an integral over a degenerate interval [t,t] being read as 0), 1S for the function equal to 1 on a set S and 0 off it, and 1{⋯} for the indicator of the condition inside the braces. The data are:
(D1) (Clocks.) For each a∈A a stochastic process Ya=(Yua)u≥0 on (Ω,F,P), every path of which is a counting path. (No law is prescribed here; the probabilistic input is hypothesis (H2) below. It holds, for instance, when the Ya are independent Poisson clocks with a horizon R>ΛT and, for each r, the levels Tra defined below are Fr-measurable and satisfy the clock-reading bound of Fresh-Start Property for Independent Poisson Clocks Read at Levels Satisfying a Clock-Reading Bound with past Fr: that lemma applied at r with cˉ=Λr gives R−=R−Λr>Λ(T−r).)
(D2) (Rates and transitions.) For each a∈A a rate function ga:[0,T]×E→[0,Λ] such that u↦ga(u,x) is measurable on [0,T] with respect to B[0,T] for every x∈E, and a transition map ϕa:E→E.
(D3) (State process.) A filtration (Ft)t∈[0,T] of sub-σ-algebras of F (time index restricted to [0,T]), an event Ω0∈F with P(Ω0)=1, and a family X=(Xt)t∈[0,T] of maps Xt:Ω→E such that {Xt=x}∈Ft for all t∈[0,T] and x∈E, and such that for every x∈E the map (u,ω)↦1Ω0(ω)1{Xu(ω)=x} is measurable with respect to the product σ-algebra B[0,T]⊗F.
(D4) (Consumed clock times and counters.) For each a∈A and t∈[0,T] a random variable Tta:Ω→[0,Λt] (the consumed clock time) such that for every ω∈Ω0 and every t∈[0,T]
Tta(ω)=∫[0,t]ga(u,Xu(ω))du,
where the integrand u↦1Ω0(ω)ga(u,Xu(ω))=∑x∈Ega(u,x)1Ω0(ω)1{Xu(ω)=x} is, by (D2) and (D3), a [0,Λ]-valued section of a B[0,T]⊗F-measurable map, so that the integral exists for every ω∈Ω0 and is 0 at t=0. The counters are Nta=YTtaa. (For instance Tta=∫[0,t]1Ω0ga(u,Xu)du qualifies: it is F-measurable by the Tonelli theorem applied on [0,T]×Ω to the integrand multiplied by 1[0,t](u), the integral over [0,t] and the integral over [0,T] of the product being the integrals over R of the same zero extension, claim 2 of Restricted Lebesgue Measure and Integral Toolkit on a Compact Interval.) The hypotheses are:
(H1) (Pathwise structure on Ω0.) For every ω∈Ω0: the path t↦Xt(ω) is piecewise constant and right-continuous, that is, there are a count K, either 0 or a natural number, and times 0<t1<⋯<tK≤T such that the path is constant on [0,t1), on [tk,tk+1) for k∈{1,…,K−1}, and on [tK,T] (constant on [0,T] when K=0); for t∈(0,T] write Xt−(ω) for the common value of Xs(ω) over s∈(t−ε,t) for some ε>0, which exists and does not depend on ε by this piecewise-constant structure. Each path t↦Nta(ω) is nondecreasing (the integrand of Ta is nonnegative and every path of Ya is a counting path), so its left limit Nt−a=sups<tNsa exists for t∈(0,T], and likewise for the grand total Nt=∑a∈ANta, whose left limit is Nt−=∑aNt−a (the supremum of a finite sum of nondecreasing functions being the sum of the suprema). It is required that for every t∈(0,T]: if Nt=Nt− then Xt=Xt−, while if Nt=Nt−+1 then Xt=ϕa(Xt−) for the unique label a with Nta=Nt−a (the increments Ntb−Nt−b being nonnegative integers with sum 1). Nothing is required at times t with Nt≥Nt−+2.
(H2) (Fresh start.) For every r∈[0,T], writing ϱr=Λ(T−r) for the residual horizon, the residual clocks Y^ua,r=YTra+ua−YTraa (u≥0, a∈A) satisfy: each Y^ua,r is a random variable; for every a and all real 0=u0<u1<⋯<up≤ϱr the increments Y^u1a,r−Y^u0a,r,…,Y^upa,r−Y^up−1a,r are independent, the q-th having the Poisson distribution with parameter uq−uq−1; and the family of σ-algebras consisting of Fr together with the σ-algebras generated σ(Y^ua,r:0≤u≤ϱr), one for each a∈A, is independent.
Then, writing E for the expectation:
(a) (Forward equation in conditioned form.) For every function F:E→R, all 0≤r≤t≤T, and every event D∈Fr,
E[(F(Xt)−F(Xr))1D]=E[1D∫[r,t]1Ω0a∈A∑ga(u,Xu)(F(ϕa(Xu))−F(Xu))du],
where the inner integral (read as 0 when t=r) is defined for every ω and is F-measurable, its integrand being a bounded section of a B[0,T]⊗F-measurable map as above.
(b) (Solution of the forward equation.) Let r∈[0,T) and D∈Fr, and put μuD(x)=P(D∩{Xu=x}) for u∈[r,T] and x∈E. For u∈[0,T] and x=y in E put qu(x,y)=∑a∈A: ϕa(x)=yga(u,x) (an empty sum being 0). Then u↦qu(x,y) is measurable on [0,T], ∑y=xqu(x,y)≤∣A∣Λ for all u and x, and (μuD)u∈[r,T] is a solution of the forward equation on [r,T] for the rates q in the sense of that lemma, applied with the rate bound there taken to be ∣A∣Λ and with its pairing μ(F)=∑x∈Eμ(x)F(x).