形式陈述
概率转移模型
本页先固定有限状态、时间齐次的离散时间 Markov 链,写作
其中 有限, 是初始分布, 且对每个 有 。它复用Markov 链公理库Markov 链Markov chain未来条件分布在给定当前状态后与更早历史无关的随机过程。的路径测度,并以状态标记 解释逻辑命题。
路径公式 诱导可测路径事件,概率算子可写成
状态 满足该公式,当从 出发的路径中满足 的概率公理库概率分布Probability distribution · Law可测空间上总质量为一的测度;随机变量的律是由样本概率推出的一类分布。至少为 。
这扩展普通模型检查公理库模型检查问题Model checking problem · Model checking给定系统模型与形式规格,判定所有指定初始行为是否满足公式,并交付证明结果或诊断见证。的布尔答案:工具可以计算精确或有误差界的概率,再与阈值比较。
到达概率方程
目标集合为 。最终到达 的概率 满足
对不可能到达 的状态置零,在剩余状态上解线性方程或用 value iteration。若没有先处理不进入 的封闭强连通类,从全一向量迭代甚至可能停在非最小不动点;到达概率是上述方程的最小非负解,不能把任意解当作概率。
线性方程的解是语义对象;浮点求解的舍入和停止误差还需给上下界,尤其当结果接近阈值 时不能用未经证明的近似直接判真。
用联合核比较两个有限概率程序
固定有限状态集 、两个随机矩阵 和初始分布 。这对应已明确随机选择概率、没有未解决调度选择的两个离散时间程序。选联合初始分布 ,使其边缘为 ;再给出乘积空间上的随机矩阵 ,对每个 满足
还须 、每行和为一。这就是耦合公理库耦合法Coupling method · Probability coupling在共同概率空间中构造具有指定边缘的随机变量,并用它们相遇的概率比较分布。从一次抽样扩展为逐步联合转移的条件。[4, §3.1] 联合过程按 演化;它是证明中构造的共同实验,不要求两个实际程序共享随机源。
路径边缘定理。 任意固定有限路径 的联合过程第一边缘概率恰为
证明从最内层 求和开始:最后一个 按边缘等式变成不依赖 的 ,可提出求和;再消去 ,依次重复,最后 。第二边缘对称。因此不仅每个时刻的状态分布正确,整段有限路径的联合分布也正确。仅检查几个时刻的一维分布不足以替代这个结论。
给定关系 ,若
则 是联合过程的概率一归纳不变式公理库不变式与归纳不变式Invariant · Inductive invariant · Strengthened invariant区分所有可达状态上成立的性质与由初始性和一步闭包直接证明的归纳不变式。。初始成立;若第 步质量全在 ,每个有质量的行都把全部质量送回 ,所以第 步仍成立。对有限多个时刻取并,整段路径逐点位于 的概率为一。
事件比较定理。 若对每对逐点满足 的长度 路径,都有
则 。在联合实验中,指标满足 几乎处处;取期望并使用路径边缘定理即得结论。关系保持与事件蕴含是两项不同义务:任意不变量都不自动比较任意事件。
直觉
冗余服务轨迹
两个独立组件每轮各以概率 失效,系统在两个都失效时进入 down。若模型只做一轮,失败概率是 ;多轮允许修复或持久故障时,状态必须记录每个组件当前健康性,不能把单轮概率重复相乘而忽略状态。
从 (up,up) 到四种健康组合的概率由乘积给出,下一轮再按当前组合转移。目标 reachability down 的概率通过上述方程累计所有不同时间的首次失败路径。
若两组件共享电源,独立性假设错误,乘积模型会严重低估风险。概率模型检查忠实计算模型分布,不验证分布参数和独立假设来自现实。
共用随机数是构造方法,边缘等式才是证书
设一个组件尚未失效时下一步失效概率为 ,另一个为 。可共同抽 :第一个在 时失效,第二个在 时失效。较可靠组件失效时,另一组件也失效,因而能逐条路径比较“某时刻前已失效”的事件。
但“使用同一随机数”本身没有验证力;要明确映射和概率。若强行让两条链始终取相同的新状态,却只按 抛硬币,第二条链的失效概率会被改成 。这种伪证书看似保持相等,实际不满足第二边缘。
例子与边界
两条吸收失效链的完整证书
状态 表示正常, 表示已失效且不再修复。令 ,
在全部四个配对状态上给出联合核;表外没有其他转移:
| 当前配对 |
下一配对及概率 |
|
|
|
|
|
|
|
|
例如 行中,第一边缘转到 的概率为 ,转到 为 ;第二边缘转到 为 ,转到 为 。 行第一边缘为 、第二边缘恒为 ; 行相反地让第一边缘恒为 、第二边缘为 ; 行两个边缘都吸收。非负性和每行归一化也逐项成立。
关系 即 。初态在 ,上述三行都不进入 ,故关系保持。取 “第一条链截至第 步(含)已失效”、“第二条链截至第 步(含)已失效”,逐点次序直接给出事件蕴含,所以前者概率不超过后者。
令 ,完整可达分布为:
| 时刻 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
第二步的 质量为 , 质量为 。独立计算单链也得到 与 ,核验两个边缘。
保持边缘的联合核与两种不等概率 终态不等与路径曾经不等
两状态空间相同时,任意终态事件 满足
因为在相等时两个指标一致。对所有事件取上确界就是总变差上界。[5, §2.2] 本例只有 不相等,由 和 得
在 ,失效事件的差本身就是 ,故终态上界取到等号。
若 依赖整个有限前缀,直接的耦合论证必须改成
在本例中,首次分歧必须在仍停于 时走 ,这些首次时刻互斥,故右边为
时是 。其中路径 的概率为 :它曾经不等,终态却相等,正好解释 。这是本联合构造给出的路径上界,不宣称它是路径分布间最小可能距离。
两链从 出发已经相等,之后仍可分离,因为 。因此不能把“首次相遇时间”代入通常相遇后锁合的耦合界。另若使用独立转移耦合,边缘依然正确,但在 会以 的概率进入 ,在本例中它等于 ,不能用它证明次序不变量。
MDP 与调度器
Markov decision process 同时含非确定动作与概率后继。固定 scheduler 后才得到路径概率测度,因此性质通常询问
对有限 MDP,先约定每个状态的动作集非空。目标 的最大到达概率在 上为一,在 上满足 Bellman 方程
最大到达概率是带上述目标边界条件的最小非负解;不能仅因某个向量满足方程就认定它是所求概率。
把非确定选择平均成均匀概率擅自加入策略,会得到另一个模型。最坏、最好和某个控制策略下的概率必须分别报告。
对有限 MDP reachability,存在 deterministic memoryless scheduler 达到最优值,但这项结论依赖目标和模型类型。带多目标、部分可观察或历史约束的策略可能需要随机化或记忆,不能把 reachability 特例推广到所有概率规格。
若 scheduler 能观察现实控制器看不到的隐藏状态,计算出的最优概率不可实现。部分可观察模型必须限制策略信息,问题复杂度也会改变。
推论与应用
有限核证书的检查边界
给定有理数表,检查者逐行验证非负性、归一化、两项边缘等式,再验证联合初始分布的非负性、归一化、两项初始边缘等式及 ,并对每个 检查 。事件蕴含仍须由事件定义证明;有限时域可以枚举路径或采用积状态记录事件。接受这些义务后,前述有限和证明给出可靠的概率比较,无须采样猜测。
若允许每步至多 的概率离开关系,必须另外累积坏事件质量;不能继续使用严格不变量的概率一结论。若模型包含调度非确定性,联合核还要与所声称的两侧调度器量词一致。本页证书固定两个随机矩阵,不直接证明任意 MDP 的最优概率比较。
奖励与长期性质
给状态或转移附 reward,可检查到目标前期望成本、有限时域累计奖励或长期平均。期望有限需要到达性和无正成本逃逸循环等条件;存在永不到目标的正概率时,expected time 可能为无穷。
稳态概率只对适当遍历类有唯一极限。多 bottom SCC 的链会让长期分布依赖初态和吸收概率,不能无条件输出一个全局稳态向量。
连续时间 Markov 链用 rate matrix ;离开状态的等待时间呈指数分布,总 exit rate 决定驻留时间。把 rate 直接归一化成离散概率只保留下一跳,丢失时间界性质。CSL 的 bounded until 需结合 uniformization 或瞬态分析。
reward 若在状态驻留期间按时间累计,CTMC 的期望奖励还要乘停留时间;把它当每跳固定奖励会得到另一个量。
PCTL、CSL 等逻辑分别服务离散/连续时间模型;时间界、next 和 until 的语义随模型改变,不能把公式名字相同当成完全可互换。
误差界与阈值判断
value iteration 得到近似向量 时,迭代差小不自动给出到真解的误差界,尤其在近吸收慢混合链中。可靠工具可维护单调下、上界 ,只有当两界都落在阈值同侧才返回布尔结论。
若 落在 内,应继续计算或返回 unknown;用打印小数四舍五入后比较会把数值误差伪装成概率证明。
参数概率模型中的转移含未知参数,结果可能是有理函数或分段区域。代入一个标称参数得到的答案不能代表整个允许区间。
统计估计的不同保证
statistical model checking 通过采样估计满足概率,给置信误差而非穷举式精确结果。它适合状态巨大模拟器,却可能漏极稀有反例;样本独立、停止规则和显著性都要写入保证。不能把 Monte Carlo “未观察到失败”称为传统概率模型检查已证明零概率。
稀有事件 importance sampling 会改变采样分布并用 likelihood ratio 校正。若校正权重错误,估计可有极低表面方差却系统偏置;模拟加速仍需统计证明。
从有限概率模型转向可不终止的程序时,次概率程序语义公理库次概率程序语义Subprobabilistic program semantics将概率程序解释为输入到终止输出的次概率核,用缺失质量记录不终止,并逐步计算顺序组合。用缺失输出质量保存非终止概率。做状态约简时,概率互模拟公理库概率互模拟Probabilistic bisimulation按等价类总转移概率定义有限 Markov 链的强互模拟,算出可合并状态并展示概率差异如何被观察。则要求到每个等价类的总质量一致;只保留相同出边支撑会漏掉概率差异,不能保证约简后的定量性质正确。
参考资料
[1] Christel Baier and Joost-Pieter Katoen, Principles of Model Checking, MIT Press, 2008, Chs. 10–11。
[2] Hans Hansson and Bengt Jonsson, “A Logic for Reasoning about Time and Reliability,” Formal Aspects of Computing 6, 1994, pp. 512–535。
[3] Marta Kwiatkowska, Gethin Norman, and David Parker, “PRISM 4.0: Verification of Probabilistic Real-Time Systems,” CAV, 2011, pp. 585–591。
[4] Lasse Leskelä, “Stochastic relations of random variables and processes”, arXiv:0806.3562,§3.1 核耦合与 Theorem 3.3,§4.3 Theorem 4.8。
[5] Gilles Barthe et al., “Relational reasoning via probabilistic coupling”, 2015,§§2–2.2,关系提升与总变差的耦合界。本文另行展开有限路径边缘、事件蕴含和完整有理数算例。