Eq. (11) — closed form of the cost per period of a stationary policy
ProvedVeinottWagnerSS.RenewalCost.eq11_aCost_closed_formLet the demands be i.i.d. with distribution on , let be the discount factor, the set-up cost and the one-period cost, and let be integers with . Let be the equivalent cost per period of the stationary policy started from stock , where is its total expected discounted cost. Then
with given by (7), the discount renewal function and .
This is the criterion minimized by the paper's algorithm for computing an optimal policy (Section 4).
Formalization Note The model is the paper's reduced model (Eq. (2)): unit purchase cost and one-period cost . is the expected cost of the controlled Markov chain, not a closed form. and are defined by their series ((7) and the first line of (9)). and are passed as natural numbers via toNat. The paper's standing assumptions that is convex and tends to at are not assumed; for the identity holds for every .
import Mathlib import Definitions.Def_VeinottWagnerSS_RenewalCost_Demand import Definitions.Def_VeinottWagnerSS_RenewalCost_RenewalFunctions import Definitions.Def_VeinottWagnerSS_RenewalCost_PolicyCost
namespace VeinottWagnerSS.RenewalCost
/-- Veinott & Wagner (1965), Eq. (11), p. 534: for `0 ≤ α < 1`, `K ≥ 0`, integers `s ≤ S` and
`D = S − s`, the equivalent cost per period `a_α(x) = (1 − α) f(x | s, S)` of the stationary
`(s, S)` policy is
`a_α(x) = (L_α(S, D) + K)/(1 + M_α(D))` for `x < s` and
`a_α(x) = (1 − α) L_α(x, x − s) + [(L_α(S, D) + K)/(1 + M_α(D))] r_α(x − s)` for `x ≥ s`. -/
theorem eq11_aCost_closed_form (D : DemandDist) (α K : ℝ) (hα0 : 0 ≤ α) (hα1 : α < 1)
(hK : 0 ≤ K) (G : ℤ → ℝ) (s S : ℤ) (hsS : s ≤ S) (x : ℤ) :
aCost D α K G s S x =
if x < s then
(LAlpha D.φ α G S (S - s).toNat + K) / (1 + MAlpha D.φ α (S - s).toNat)
else
(1 - α) * LAlpha D.φ α G x (x - s).toNat
+ (LAlpha D.φ α G S (S - s).toNat + K) / (1 + MAlpha D.φ α (S - s).toNat)
* rAlpha D.φ α (x - s).toNat := by sorry
end VeinottWagnerSS.RenewalCost
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.