研究生课程 · 不确定规划

不确定规划

当你掷硬币 10000 次,频率给了概率;当你只有一句专家估计"工期大约 8 到 12 天",频率无从谈起——这是信度的世界。不确定理论(刘宝碇)为后者建立公理体系:不确定测度、不确定分布、期望值模型、机会约束规划,最后用仿真 + 智能算法求解。

1 为什么需要不确定理论

概念两种"不确定"

随机性(概率):事件有明确的频率解释——大量重复试验中比例稳定。适用前提:未来与过去统计同质、样本充足。
不确定性(信度):无法重复、无样本——全新项目的工期、政策的出台时间、某块油田的储量。此时"估计"来自专家经验,硬套概率会得出严重偏离的结论(Liu 的著名反例:桥的寿命)。
判别标准(Liu 建议):有足够样本 → 概率;无样本、仅信度 → 不确定理论;两者混合 → 不确定随机规划。
为什么不用概率硬套:贝叶斯学派也用"主观概率",但其运算基于可加性;而"信度"的实验(Liu 十字模型)表明人对不互斥事件的信心不满足可加性,只满足对偶性 + 次可加性——这正是不确定公理的出发点。

2 不确定测度与公理体系

公理不确定测度 M

集函数 M: 2^Γ → [0,1] 满足:
公理 1(正规性) M{Γ} = 1;
公理 2(对偶性) M{Λ} + M{Λᶜ} = 1 对任意事件 Λ;
公理 3(次可加性) M{∪Λᵢ} ≤ Σ M{Λᵢ}(对可数事件族)。
与概率的根本区别:概率要求可加(等号 + 互斥),不确定测度只要求次可加——"对两个都可能发生的事,整体信心 ≤ 信心之和"。由此可推出:M{∅}=0、单调性、M{Λ₁∩Λ₂} ≥ M{Λ₁}+M{Λ₂}−1(对偶次可加)。

操作对偶性与次可加性演示

调节你对事件 A 的信度 M{A},观察补事件的信度被对偶性"锁死";再看两事件并集的信度范围:

随堂检测公理体系

3 不确定变量与不确定分布

定义Φ(x) = M{ξ ≤ x}

不确定分布是不确定测度作用在 ξ 上得到的单调不减函数,值域 (0,1),不要求左右连续(区别于概率分布函数)。
三大标准分布:
线性 L(a,b):Φ 在 [a,b] 上从 0 线性升到 1;
之字 Z(a,b,c):两段折线,b 处取 ½;
正态 N(e,σ):Φ(x) = [1 + exp(π(e−x)/(√3σ))]⁻¹(与概率正态形似而神异——它来自 Sigmoid)。
经验不确定分布 / 99 表:专家给出 0.01…0.99 分位点的"不确定量表",是工程中采集信度数据的标准格式。

操作不确定分布实验室 ⭐(三大分布 + 99 表)

拖动 α 滑块,虚线显示 α-乐观值与 α-悲观值(Φ⁻¹(1−α) 与 Φ⁻¹(α))——机会约束规划的计算核心。

4 数字特征与关键值

定义期望、方差、乐观值

期望:E[ξ] = ∫₀^{+∞} M{ξ≥r}dr − ∫_{−∞}^0 M{ξ≤r}dr(对一般变量),标准分布有闭式:E[L(a,b)] = (a+b)/2;E[Z(a,b,c)] = (a+2b+c)/4;E[N(e,σ)] = e。
方差:E[(ξ−e)²],Var[L] = (b−a)²/12,Var[N(e,σ)] = σ²。
α-乐观值:ξ_sup(α) = sup{x | M{ξ≥x} ≥ α} = Φ⁻¹(1−α);α-悲观值 = Φ⁻¹(α)。工程上"90% 把握能完成的最短工期"就是乐观值。

随堂检测数字特征

5 不确定规划模型与仿真

模型EVM、CCP 与相关机会

期望值模型(EVM):max E[f(x,ξ)] s.t. E[gᵢ(x,ξ)] ≤ 0——把不确定目标/约束折成期望。
机会约束规划(CCP):max f̄ s.t. M{gᵢ(x,ξ) ≤ 0} ≥ αᵢ——约束只须以置信水平 αᵢ 成立("九成把握按时完工")。
相关机会规划(RCP):max M{事件 h(x,ξ) ≤ 0 发生}——最大化"事情办成"的机会。
求解:难以解析时用不确定仿真(99 表/逆分布采样生成"伪样本")评估目标与约束,外层接遗传算法/粒子群等智能搜索——"仿真 + 智能算法"是不确定规划的标准装备。

操作机会约束求解器(把握多少才算达标)

约束 x + ξ ≥ T,其中 ξ ~ L(1, 3)(工期类变量)。左图:满足测度 M 随决策 x 的变化;右下:99 表与逆分布采样直方图:

操作99 表与不确定仿真(信度 → 样本)

注意:直方图是仿真样本的频率,蓝线是信度曲线 Φ。二者定义不同(频率 vs 信度),仿真只是求解工具,不代表 ξ"真的"服从该频率分布。

随堂检测规划模型