Skip to content

定义Definition

概率模型检查

Probabilistic model checking · Quantitative model checking

在Markov链上计算事件概率,并以保持边缘的联合转移核和关系不变量证明概率比较与误差界。

形式陈述 ​

概率转移模型 ​

本页先固定有限状态、时间齐次的离散时间 Markov 链,写作

D=(S,P,μ0,L),

其中 S 有限,μ0 是初始分布,P(s,s′)∈[0,1] 且对每个 s 有 ∑s′P(s,s′)=1。它复用Markov 链的路径测度,并以状态标记 L 解释逻辑命题。

路径公式 ψ 诱导可测路径事件,概率算子可写成

P≥p[ψ].

状态 s 满足该公式,当从 s 出发的路径中满足 ψ 的概率至少为 p。

这扩展普通模型检查的布尔答案:工具可以计算精确或有误差界的概率,再与阈值比较。

到达概率方程 ​

目标集合为 G。最终到达 G 的概率 xs 满足

xs=1(s∈G),xs=∑s′P(s,s′)xs′(s∉G).

对不可能到达 G 的状态置零,在剩余状态上解线性方程或用 value iteration。若没有先处理不进入 G 的封闭强连通类,从全一向量迭代甚至可能停在非最小不动点;到达概率是上述方程的最小非负解,不能把任意解当作概率。

线性方程的解是语义对象;浮点求解的舍入和停止误差还需给上下界,尤其当结果接近阈值 p 时不能用未经证明的近似直接判真。

用联合核比较两个有限概率程序 ​

固定有限状态集 S,T、两个随机矩阵 P,Q 和初始分布 μ0,ν0。这对应已明确随机选择概率、没有未解决调度选择的两个离散时间程序。选联合初始分布 λ0,使其边缘为 μ0,ν0;再给出乘积空间上的随机矩阵 K,对每个 (x,y) 满足

(边缘)∑y′∈TK((x,y),(x′,y′))=P(x,x′),∑x′∈SK((x,y),(x′,y′))=Q(y,y′).

还须 K≥0、每行和为一。这就是耦合从一次抽样扩展为逐步联合转移的条件。[4, §3.1] 联合过程按 λn+1=λnK 演化;它是证明中构造的共同实验,不要求两个实际程序共享随机源。

路径边缘定理。 任意固定有限路径 x0,…,xn 的联合过程第一边缘概率恰为

∑y0,…,ynλ0(x0,y0)∏t=0n−1K((xt,yt),(xt+1,yt+1))=μ0(x0)∏t=0n−1P(xt,xt+1).

证明从最内层 yn 求和开始:最后一个 K 按边缘等式变成不依赖 yn−1 的 P(xn−1,xn),可提出求和;再消去 yn−1,依次重复,最后 ∑y0λ0(x0,y0)=μ0(x0)。第二边缘对称。因此不仅每个时刻的状态分布正确,整段有限路径的联合分布也正确。仅检查几个时刻的一维分布不足以替代这个结论。

给定关系 R⊆S×T,若

λ0(R)=1,(x,y)∈R⟹K((x,y),R)=1,

则 R 是联合过程的概率一归纳不变式。初始成立;若第 t 步质量全在 R,每个有质量的行都把全部质量送回 R,所以第 t+1 步仍成立。对有限多个时刻取并,整段路径逐点位于 R 的概率为一。

事件比较定理。 若对每对逐点满足 R 的长度 n+1 路径,都有

A(x0:n)⟹B(y0:n),

则 PrP,μ0(A)≤PrQ,ν0(B)。在联合实验中,指标满足 1A≤1B 几乎处处;取期望并使用路径边缘定理即得结论。关系保持与事件蕴含是两项不同义务:任意不变量都不自动比较任意事件。

直觉

冗余服务轨迹 ​

两个独立组件每轮各以概率 0.01 失效,系统在两个都失效时进入 down。若模型只做一轮,失败概率是 10−4;多轮允许修复或持久故障时,状态必须记录每个组件当前健康性,不能把单轮概率重复相乘而忽略状态。

从 (up,up) 到四种健康组合的概率由乘积给出,下一轮再按当前组合转移。目标 reachability down 的概率通过上述方程累计所有不同时间的首次失败路径。

若两组件共享电源,独立性假设错误,乘积模型会严重低估风险。概率模型检查忠实计算模型分布,不验证分布参数和独立假设来自现实。

共用随机数是构造方法,边缘等式才是证书 ​

设一个组件尚未失效时下一步失效概率为 p,另一个为 q≥p。可共同抽 U∼Unif[0,1]:第一个在 U≤p 时失效,第二个在 U≤q 时失效。较可靠组件失效时,另一组件也失效,因而能逐条路径比较“某时刻前已失效”的事件。

但“使用同一随机数”本身没有验证力;要明确映射和概率。若强行让两条链始终取相同的新状态,却只按 p 抛硬币,第二条链的失效概率会被改成 p。这种伪证书看似保持相等,实际不满足第二边缘。

例子与边界

两条吸收失效链的完整证书 ​

状态 0 表示正常,1 表示已失效且不再修复。令 0≤p≤q≤1,

P=(1−pp01),Q=(1−qq01),λ0=δ(0,0).

在全部四个配对状态上给出联合核;表外没有其他转移:

当前配对 下一配对及概率
00 00:1−q,01:q−p,11:p
01 01:1−p,11:p
10 10:1−q,11:q
11 11:1

例如 00 行中,第一边缘转到 0 的概率为 (1−q)+(q−p)=1−p,转到 1 为 p;第二边缘转到 0 为 1−q,转到 1 为 (q−p)+p=q。01 行第一边缘为 (1−p,p)、第二边缘恒为 1;10 行相反地让第一边缘恒为 1、第二边缘为 (1−q,q);11 行两个边缘都吸收。非负性和每行归一化也逐项成立。

关系 R={00,01,11} 即 x≤y。初态在 R,上述三行都不进入 10,故关系保持。取 A=“第一条链截至第 n 步(含)已失效”、B=“第二条链截至第 n 步(含)已失效”,逐点次序直接给出事件蕴含,所以前者概率不超过后者。

令 p=1/4,q=1/2,完整可达分布为:

时刻 λn(00) λn(01) λn(11) Pr(Xn=1) Pr(Yn=1)
0 1 0 0 0 0
1 1/2 1/4 1/4 1/4 1/2
2 1/4 5/16 7/16 7/16 3/4

第二步的 01 质量为 (1/2)(1/4)+(1/4)(3/4)=5/16,11 质量为 (1/2)(1/4)+(1/4)(1/4)+1/4=7/16。独立计算单链也得到 1−(3/4)2=7/16 与 1−(1/2)2=3/4,核验两个边缘。

保持边缘的联合核与两种不等概率

终态不等与路径曾经不等 ​

两状态空间相同时,任意终态事件 A⊆S 满足

|Pr(Xn∈A)−Pr(Yn∈A)|≤Pr(Xn≠Yn),

因为在相等时两个指标一致。对所有事件取上确界就是总变差上界。[5, §2.2] 本例只有 01 不相等,由 Pr(Xn=0)=(1−p)n 和 Pr(Yn=0)=(1−q)n 得

Pr(Xn≠Yn)=(1−p)n−(1−q)n.

在 n=2,失效事件的差本身就是 5/16,故终态上界取到等号。

若 A 依赖整个有限前缀,直接的耦合论证必须改成

|Pr(X0:n∈A)−Pr(Y0:n∈A)|≤Pr(∃t≤n:Xt≠Yt).

在本例中,首次分歧必须在仍停于 00 时走 00→01,这些首次时刻互斥,故右边为

(q−p)∑t=1n(1−q)t−1.

n=2 时是 1/4+(1/2)(1/4)=3/8。其中路径 00→01→11 的概率为 1/16:它曾经不等,终态却相等,正好解释 3/8−5/16。这是本联合构造给出的路径上界,不宣称它是路径分布间最小可能距离。

两链从 00 出发已经相等,之后仍可分离,因为 P≠Q。因此不能把“首次相遇时间”代入通常相遇后锁合的耦合界。另若使用独立转移耦合,边缘依然正确,但在 00 会以 p(1−q) 的概率进入 10,在本例中它等于 1/8>0,不能用它证明次序不变量。

MDP 与调度器 ​

Markov decision process 同时含非确定动作与概率后继。固定 scheduler 后才得到路径概率测度,因此性质通常询问

infσPrsσ(ψ),supσPrsσ(ψ).

对有限 MDP,先约定每个状态的动作集非空。目标 G 的最大到达概率在 s∈G 上为一,在 s∉G 上满足 Bellman 方程

xs=maxa∈Act(s)∑s′P(s,a,s′)xs′.

最大到达概率是带上述目标边界条件的最小非负解;不能仅因某个向量满足方程就认定它是所求概率。

把非确定选择平均成均匀概率擅自加入策略,会得到另一个模型。最坏、最好和某个控制策略下的概率必须分别报告。

对有限 MDP reachability,存在 deterministic memoryless scheduler 达到最优值,但这项结论依赖目标和模型类型。带多目标、部分可观察或历史约束的策略可能需要随机化或记忆,不能把 reachability 特例推广到所有概率规格。

若 scheduler 能观察现实控制器看不到的隐藏状态,计算出的最优概率不可实现。部分可观察模型必须限制策略信息,问题复杂度也会改变。

推论与应用

有限核证书的检查边界 ​

给定有理数表,检查者逐行验证非负性、归一化、两项边缘等式,再验证联合初始分布的非负性、归一化、两项初始边缘等式及 λ0(R)=1,并对每个 (x,y)∈R 检查 K((x,y),Rc)=0。事件蕴含仍须由事件定义证明;有限时域可以枚举路径或采用积状态记录事件。接受这些义务后,前述有限和证明给出可靠的概率比较,无须采样猜测。

若允许每步至多 εt 的概率离开关系,必须另外累积坏事件质量;不能继续使用严格不变量的概率一结论。若模型包含调度非确定性,联合核还要与所声称的两侧调度器量词一致。本页证书固定两个随机矩阵,不直接证明任意 MDP 的最优概率比较。

奖励与长期性质 ​

给状态或转移附 reward,可检查到目标前期望成本、有限时域累计奖励或长期平均。期望有限需要到达性和无正成本逃逸循环等条件;存在永不到目标的正概率时,expected time 可能为无穷。

稳态概率只对适当遍历类有唯一极限。多 bottom SCC 的链会让长期分布依赖初态和吸收概率,不能无条件输出一个全局稳态向量。

连续时间 Markov 链用 rate matrix R(s,s′);离开状态的等待时间呈指数分布,总 exit rate 决定驻留时间。把 rate 直接归一化成离散概率只保留下一跳,丢失时间界性质。CSL 的 bounded until 需结合 uniformization 或瞬态分析。

reward 若在状态驻留期间按时间累计,CTMC 的期望奖励还要乘停留时间;把它当每跳固定奖励会得到另一个量。

PCTL、CSL 等逻辑分别服务离散/连续时间模型;时间界、next 和 until 的语义随模型改变,不能把公式名字相同当成完全可互换。

误差界与阈值判断 ​

value iteration 得到近似向量 x(n) 时,迭代差小不自动给出到真解的误差界,尤其在近吸收慢混合链中。可靠工具可维护单调下、上界 ln≤x≤un,只有当两界都落在阈值同侧才返回布尔结论。

若 p 落在 [ln,un] 内,应继续计算或返回 unknown;用打印小数四舍五入后比较会把数值误差伪装成概率证明。

参数概率模型中的转移含未知参数,结果可能是有理函数或分段区域。代入一个标称参数得到的答案不能代表整个允许区间。

统计估计的不同保证 ​

statistical model checking 通过采样估计满足概率,给置信误差而非穷举式精确结果。它适合状态巨大模拟器,却可能漏极稀有反例;样本独立、停止规则和显著性都要写入保证。不能把 Monte Carlo “未观察到失败”称为传统概率模型检查已证明零概率。

稀有事件 importance sampling 会改变采样分布并用 likelihood ratio 校正。若校正权重错误,估计可有极低表面方差却系统偏置;模拟加速仍需统计证明。

从有限概率模型转向可不终止的程序时,次概率程序语义用缺失输出质量保存非终止概率。做状态约简时,概率互模拟则要求到每个等价类的总质量一致;只保留相同出边支撑会漏掉概率差异,不能保证约简后的定量性质正确。

参考资料

[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,关系提升与总变差的耦合界。本文另行展开有限路径边缘、事件蕴含和完整有理数算例。

关系图谱19 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系