Eqs. (29)-(30): gain rate and relative values of a two-state Markov-renewal process
ProvedJewellMRP.GainRate.two_state_gain_rateaverage-rewardmarkov-renewalp2o-batch-p200bp2o-gran-per-chapterp2o-plan-paperp2o-v1
Consider a Markov-renewal program with two states and a stationary policy whose chain is
with mean sojourn times and expected rewards . Let
Then solves the value-determination equations (13) of , and equals the gain rate (B 7) of computed from any stationary probability vector of .
These closed forms let the optimal policy of a two-state program be found by direct evaluation, and they give an explicit, nonvacuous instance of the objects of the algorithm of Fig. 2.
Formalization Note States are indices of Fin 2. The hypotheses make the chain ergodic, as the paper notes below (21); no ergodicity of other policies is assumed.
Preamble
import Mathlib import Definitions.Def_JewellMRP_GainRate_MRP import Definitions.Def_JewellMRP_GainRate_PolicyIteration
Formal statement
namespace JewellMRP.GainRate
theorem two_state_gain_rate {α : Type*} (M : MRP 2 α) (z : Fin 2 → α)
(h12 : 0 < M.p 0 (z 0) 1) (h21 : 0 < M.p 1 (z 1) 0) :
let p12 := M.p 0 (z 0) 1
let p21 := M.p 1 (z 1) 0
let ν1 := M.ν 0 (z 0)
let ν2 := M.ν 1 (z 1)
let ρ1 := M.ρ 0 (z 0)
let ρ2 := M.ρ 1 (z 1)
let g := (p21 * ρ1 + p12 * ρ2) / (ν1 * p21 + ν2 * p12)
M.SolvesValueDetermination z g ![(ν2 * ρ1 - ν1 * ρ2) / (ν1 * p21 + ν2 * p12), 0] ∧
∀ π : Fin 2 → ℝ, IsStationaryDist (M.policyMatrix z) π → M.gainRate z π = g := by sorry
end JewellMRP.GainRate
Source
Jewell, Markov-Renewal Programming. II: Infinite Return Models, Example, Oper. Res. 11 (1963), p. 961, Eqs. (29)-(30), with (21) (p. 960)
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.