Theorem 2.2, pp. 4–5 — chains to the unit rows give ρ(F) ≤ (1 − (1 − α₁)α₂^{|C*|})^{1/(|C*|+1)} < 1
ProvedSubSuperStoch.SubSpec.theorem_2_2Let be an sub-stochastic matrix: and for all . Let
be non-empty. Suppose that for each there is a non-zero element chain
whose entries are all non-zero. Put
where is the number of elements of . Then the spectral radius of satisfies
The theorem gives a checkable sufficient condition for a sub-stochastic matrix to be strictly stable, , with an explicit bound: it suffices that every row of sum one can be reached by a chain of non-zero entries from a row of sum less than one. It is the matrix tool behind the convergence results for asynchronous bipartite tracking in the same paper.
Formalization Note Indices are 0-based. The chain for is a function c with c 0 = i ∈ S₁, c (len k) = k, 1 ≤ len k, entries F (c (t+1)) (c t) ≠ 0 and c (t+1) ≠ c t; the hypothesis provides one chain per , of length len k, and is the largest of these chosen lengths (IsGreatest {m | ∃ k, rowSum F k = 1 ∧ m = len k} C). and are passed as values with IsGreatest/IsLeast hypotheses over exactly the printed sets. is stated for every complex eigenvalue (specC F), and the root is the real power ; its radicand lies in .
import Mathlib import Definitions.Def_SubSuperStoch_SubSpec_Setting
namespace SubSuperStoch.SubSpec
theorem theorem_2_2 {n : ℕ} (F : Matrix (Fin n) (Fin n) ℝ) (hF : IsSubStochastic F)
(hS₁ : ∃ s, rowSum F s < 1) (hS₂ : ∃ s, rowSum F s = 1)
(len : Fin n → ℕ)
(hchain : ∀ k, rowSum F k = 1 → ∃ i, rowSum F i < 1 ∧
∃ c : ℕ → Fin n, c 0 = i ∧ c (len k) = k ∧ 1 ≤ len k ∧ IsChain F c (len k))
(α₁ : ℝ) (hα₁ : IsGreatest {x | ∃ s, rowSum F s < 1 ∧ x = rowSum F s} α₁)
(α₂ : ℝ) (hα₂ : IsLeast {x | ∃ i j, 0 < F i j ∧ x = F i j} α₂)
(C : ℕ) (hC : IsGreatest {m | ∃ k, rowSum F k = 1 ∧ m = len k} C) :
(∀ μ ∈ specC F, ‖μ‖ ≤ (1 - (1 - α₁) * α₂ ^ C) ^ (1 / ((C : ℝ) + 1))) ∧
(1 - (1 - α₁) * α₂ ^ C) ^ (1 / ((C : ℝ) + 1)) < 1 := by sorry
end SubSuperStoch.SubSpec
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.