系统分析:线性规划、非线性规划与排队模型 · 第 11 / 11 篇(展开目录)

从线性规划建模、单纯形法、对偶与灵敏度、运输与指派,到非线性规划、随机过程、马尔可夫链、排队模型和决策分析,每章都配手算例题和自测,最后一章用综合例题把方法串起来。

  1. 第1章 线性规划建模:从文字题到标准型
  2. 第2章 单纯形法:单纯形表、大 M 法与两阶段法
  3. 第3章 对偶与灵敏度分析:最终单纯形表还能回答什么
  4. 第4章 运输问题与指派问题:Vogel 法、位势法与匈牙利法
  5. 第5章 非线性规划(一):一维搜索与多元无约束优化
  6. 第6章 非线性规划(二):拉格朗日乘子与 KKT 条件
  7. 第7章 随机过程基础:几何分布、指数分布与泊松过程
  8. 第8章 马尔可夫链:离散时间与连续时间
  9. 第9章 排队模型:从 M/M/1 到 M/M/m
  10. 第10章 决策分析:不确定性准则、决策树与效用
  11. 第11章 综合例题讲解与解析

约 60 分钟 · 需要前十章的方法

怎么用这一章: 先只看题目,自己写出解题框架(建什么模型、用哪一章的方法),再对照解析。提示框会提醒建模和计算中容易忽略的细节。

例1 线性规划:建模、求解与预算削减

题目。 冰淇淋店做三种口味,每杯的配料成本、人工成本和包装材料如下;每天包装材料最多 10 kg,配料支出不超过 100 元,人工成本不超过 30 元。三种口味每杯利润分别为 2.5、2.3、2.0 元。(a) 建立利润最大的模型;(b) 每天各做多少杯;(c) 老板想削减人工预算,最多能削减多少,而(b)的方案仍可行?

口味 配料(元/杯) 人工(元/杯) 包装(kg/杯)
香草 x1x_1 1.50 0.50 0.20
巧克力 x2x_2 1.70 0.35 0.18
草莓 x3x_3 1.80 0.60 0.16

(a)

max⁡ z=2.5x1+2.3x2+2.0x3s.t.{0.2x1+0.18x2+0.16x3≤10(包装)1.5x1+1.7x2+1.8x3≤100(配料)0.5x1+0.35x2+0.6x3≤30(人工)x1,x2,x3≥0\max\ z=2.5x_1+2.3x_2+2.0x_3\quad\text{s.t.}\quad \begin{cases} 0.2x_1+0.18x_2+0.16x_3\le10&(\text{包装})\\ 1.5x_1+1.7x_2+1.8x_3\le100&(\text{配料})\\ 0.5x_1+0.35x_2+0.6x_3\le30&(\text{人工})\\ x_1,x_2,x_3\ge0 \end{cases}

(b) 加松弛变量 s1,s2,s3s_1,s_2,s_3,化为最小化 −z-z。

表 0:检验行 (−2.5, −2.3, −2, 0,0,0 ∣ 0)(-2.5,\ -2.3,\ -2,\ 0,0,0\ |\ 0)。x1x_1 进基,比值 10/0.2=5010/0.2=50、100/1.5=66.7100/1.5=66.7、30/0.5=6030/0.5=60,s1s_1 出基。

表 1:

基 x1x_1 x2x_2 x3x_3 s1s_1 s2s_2 s3s_3 右端
x1x_1 1 9/10 4/5 5 0 0 50
s2s_2 0 7/20 3/5 −15/2 1 0 25
s3s_3 0 −1/10 1/5 −5/2 0 1 5
检验行 0 −1/20 0 25/2 0 0 125

x2x_2 的检验数 −1/20 仍为负,x2x_2 进基;比值 50÷910=55.650\div\frac{9}{10}=55.6、25÷720=71.425\div\frac{7}{20}=71.4,x1x_1 出基。

表 2:

基 x1x_1 x2x_2 x3x_3 s1s_1 s2s_2 s3s_3 右端
x2x_2 10/9 1 8/9 50/9 0 0 500/9
s2s_2 −7/18 0 13/45 −85/9 1 0 50/9
s3s_3 1/9 0 13/45 −35/18 0 1 95/9
检验行 1/18 0 2/45 115/9 0 0 1150/9

检验行全部非负。线性规划最优解:只做巧克力 x2=500/9≈55.6x_2=500/9\approx55.6 杯,利润 1150/9≈127.81150/9\approx127.8 元。只有包装材料是紧约束(影子价格 115/9≈12.8115/9\approx12.8 元/kg)。

上述是允许产量连续变化的线性规划解。实际按整杯生产时,需要另加整数约束:简单取整为 55 杯巧克力,利润为 126.5 元;精确的整数最优方案是 52 杯巧克力 + 4 杯草莓,利润 127.6 元(包装恰好 9.36+0.64=109.36+0.64=10 kg)。直接取整不保证整数最优。

(c) 对连续产量的最优解,人工成本只用了 0.35×500/9=175/9≈19.440.35\times500/9=175/9\approx19.44 元,人工约束的松弛 s3=95/9≈10.56s_3=95/9\approx10.56。人工预算从 30 降到 175/9175/9 之前,当前解都可行,最多削减约 10.56 元。用第3章的语言:右端变化 Δ\Delta 时 s3s_3 变成 95/9+Δ95/9+\Delta,要求 ≥0\ge0,即 Δ≥−95/9\Delta\ge-95/9。若采用整杯最优方案 (0,52,4)(0,52,4),人工成本为 0.35×52+0.6×4=20.60.35\times52+0.6\times4=20.6 元,可削减的预算相应为 9.4 元。

要点

非紧约束能削减的量就是它的松弛量,直接从最终表右端读。“紧”还是“松”要看最终表的基变量:s3s_3 在基中且为正,说明人工预算没用完。

例2 含等式约束的单纯形法与目标系数灵敏度

题目。 (a) 用单纯形法求 min⁡ z=2x1+3x2\min\ z=2x_1+3x_2,s.t. 3x1−2x2≤33x_1-2x_2\le3,8x1+4x2=58x_1+4x_2=5,x1,x2≥0x_1,x_2\ge0。(b) 目标变为 2x1+(3+Δ)x22x_1+(3+\Delta)x_2 时,Δ\Delta 在什么范围内(a)的解仍最优?(c) 若问题是最大化型,第一阶段中人工变量的目标系数应如何取?

(a) 3x1−2x2+s1=33x_1-2x_2+s_1=3,8x1+4x2+a1=58x_1+4x_2+a_1=5,初始基 (s1,a1)(s_1,a_1)。检验行拆成常数行和 M 行:

基 x1x_1 x2x_2 s1s_1 a1a_1 右端
s1s_1 3 −2 1 0 3
a1a_1 8 4 0 1 5
常数行 2 3 0 0 0
M 行 −8 −4 0 0 −5

M 行最负 −8,x1x_1 进基;比值 3/3=13/3=1、5/85/8,a1a_1 出基:

基 x1x_1 x2x_2 s1s_1 右端
s1s_1 0 −7/2 1 9/8
x1x_1 1 1/2 0 5/8
常数行 0 2 0 −5/4

人工变量已出基(删去该列和 M 行),常数行全部非负,最优:**x1=5/8, x2=0x_1=5/8,\ x_2=0,z=5/4z=5/4**。

(b) x2x_2 是非基变量,只影响它自己的检验数:2+Δ≥02+\Delta\ge0,即 **Δ≥−2\Delta\ge-2**(c2≥1c_2\ge1)时解不变;c2c_2 增大不会改变最优解。

(c) 最大化问题中人工变量取 −M-M(目标越大越好,人工变量越大越被惩罚);若先把目标乘 −1-1 化为最小化,则取 +M+M。第一阶段等价于最小化人工变量之和,与原问题是最大还是最小无关。

例3 运输问题与运价灵敏度

题目。 三个工厂向三个地区供货(箱),单位运费如下。(a) 写成平衡运输问题;(b) 用 Vogel 法求初始解并求最优解;(c) 工厂 3 到地区 3 的运费变化 Δ\Delta,求当前解保持最优的范围。

地区 1 地区 2 地区 3 供应
工厂 1 20 24 24 60
工厂 2 18 23 20 70
工厂 3 20 25 26 50
需求 40 65 75

(a) 总供应 180 = 总需求 180,已平衡:min⁡∑cijxij\min\sum c_{ij}x_{ij},∑jxij=si\sum_jx_{ij}=s_i,∑ixij=dj\sum_ix_{ij}=d_j,xij≥0x_{ij}\ge0。

(b) Vogel 法:

步 行罚数 列罚数 选中 分配
1 4, 2, 5 2, 1, 4 行 3 格(3,1)=40
2 0, 3, 1 —, 1, 4 列 3 格(2,3)=70
3 0, —, 1 —, 1, 2 列 3 格(1,3)=5
4 24, —, 25 —, 1, — 行 3 格(3,2)=10
5 24 24 行 1 格(1,2)=55

基格 5 个 =3+3−1=3+3-1。z=24⋅55+24⋅5+20⋅70+20⋅40+25⋅10=3890z=24\cdot55+24\cdot5+20\cdot70+20\cdot40+25\cdot10=3890。

位势(u1=0u_1=0):v2=24, v3=24, u2=−4, u3=1, v1=19v_2=24,\ v_3=24,\ u_2=-4,\ u_3=1,\ v_1=19。非基格检验数:格(1,1) =1=1,(2,1) =3=3,(2,2) =3=3,(3,3) =1=1,全部为正。Vogel 初始解已经最优,最小运费 3890。

(c) 格(3,3) 是非基格,位势不变,只需 1+Δ≥01+\Delta\ge0:**Δ≥−1\Delta\ge-1**,即运费不低于 25 时方案不变(低于 25 改走这条线路更便宜)。

例4 生产—库存问题

题目。 某器材 11 月—2 月的需求为 100、200、180、300 套,各月产能 50、180、280、270 套,生产成本 40 元/套。提前生产要付保管费,每套每月 0.5 元;延期交货要付罚金,每套每月 2 元。求最优生产计划。

产能合计 780 = 需求合计 780,平衡。运价 cijc_{ij} = 40 + 保管或罚金:

生产 \ 交货 11 月 12 月 1 月 2 月 产能
11 月 40 40.5 41 41.5 50
12 月 42 40 40.5 41 180
1 月 44 42 40 40.5 280
2 月 46 44 42 40 270
需求 100 200 180 300

Vogel 法依次分配:格(4,4)=270、(1,1)=50、(2,1)=50、(2,2)=130、(3,2)=70、(3,4)=30、(3,3)=180,共 7 个基格。位势检验所有非基格检验数 ≥0\ge0(格(3,1) 为 0,存在另一个同样最优的方案),最小总成本 31455 元。

计划解读:

月份 生产 用途
11 月 50 全部用于 11 月
12 月 180 50 套补 11 月的缺口(延期一个月),130 套用于 12 月
1 月 280 70 套补 12 月(延期一个月),180 套用于 1 月,30 套存到 2 月
2 月 270 用于 2 月

要点

先判断“提前生产”和“延期交货”哪个便宜:这里保管费 0.5 < 罚金 2,但 11 月产能太小,延期交货不可避免。运价表写对了,后面就是标准的运输算法。

例5 不平衡且有禁止格的指派问题

题目。 4 人 5 项任务,“—”表示此人不能做此任务。求总成本最小的指派。

J1 J2 J3 J4 J5
P1 20 15 16 17 —
P2 — — 16 20 17
P3 12 15 11 — 17
P4 13 — 14 20 12

“—”记为 MM,加虚拟人 P5(成本全 0)。行化简(减 15、16、11、12;P5 不变),列最小值全为 0:

J1 J2 J3 J4 J5
P1 5 0 1 2 M
P2 M M 0 4 1
P3 1 4 0 M 6
P4 1 M 2 8 0
P5 0 0 0 0 0

P2、P3 都只在 J3 有 0,最多 4 个独立的 0。最少 4 条线:P1、P4、P5 三行和 J3 列。未覆盖的最小元素为 1(P2→J5、P3→J1),未覆盖减 1、交叉加 1:

J1 J2 J3 J4 J5
P1 5 0 2 2 M
P2 M M 0 3 0
P3 0 3 0 M 5
P4 1 M 3 8 0
P5 0 0 1 0 0

最优指派:P1→J2,P2→J3,P3→J1,P4→J5,J4 无人承担。总成本 15+16+12+12=5515+16+12+12=55。

例6 拉格朗日乘子:两个等式约束与右端扰动

题目一。 min⁡ f=x12+x22+x32\min\ f=x_1^2+x_2^2+x_3^2,s.t. x1+x2+3x3=2x_1+x_2+3x_3=2,5x1+2x2+x3=55x_1+2x_2+x_3=5。

2x1+λ1+5λ2=0,2x2+λ1+2λ2=0,2x3+3λ1+λ2=0.2x_1+\lambda_1+5\lambda_2=0,\quad 2x_2+\lambda_1+2\lambda_2=0,\quad 2x_3+3\lambda_1+\lambda_2=0 .

解出 xx 用 λ\lambda 表示后代入两条约束:

x=−12(λ1+5λ2, λ1+2λ2, 3λ1+λ2),{11λ1+10λ2=−410λ1+30λ2=−10x=-\tfrac12\big(\lambda_1+5\lambda_2,\ \lambda_1+2\lambda_2,\ 3\lambda_1+\lambda_2\big),\qquad \begin{cases}11\lambda_1+10\lambda_2=-4\\10\lambda_1+30\lambda_2=-10\end{cases}

得 λ1=−223\lambda_1=-\tfrac{2}{23},λ2=−723\lambda_2=-\tfrac{7}{23},于是

x∗=(3746, 823, 1346)≈(0.804, 0.348, 0.283),f∗=3946≈0.848.x^*=\Big(\tfrac{37}{46},\ \tfrac{8}{23},\ \tfrac{13}{46}\Big)\approx(0.804,\ 0.348,\ 0.283),\qquad f^*=\tfrac{39}{46}\approx0.848 .

∇XX2L=2I\nabla^2_{XX}L=2I 正定,是极小点(也是全局最小:凸目标 + 线性约束)。

题目二。 min⁡ f=3x12+2x22+12x32\min\ f=3x_1^2+2x_2^2+\tfrac12x_3^2,s.t. x1+x2=15x_1+x_2=15,x3=1x_3=1。(a) 求最优解;(b) 第一条约束改为 x1+x2=15+Δx_1+x_2=15+\Delta,最优值约变化多少?

(a) 6x1+λ1=06x_1+\lambda_1=0,4x2+λ1=04x_2+\lambda_1=0,x3+λ2=0x_3+\lambda_2=0。由前两式 x2=1.5x1x_2=1.5x_1,代入 x1+x2=15x_1+x_2=15 得 x1=6, x2=9x_1=6,\ x_2=9,x3=1x_3=1;λ1=−36, λ2=−1\lambda_1=-36,\ \lambda_2=-1。f∗=108+162+0.5=270.5f^*=108+162+0.5=270.5。

(b) Δf∗≈−λ1Δ=36Δ\Delta f^*\approx-\lambda_1\Delta=36\Delta。精确值 f∗(Δ)=270.5+36Δ+1.2Δ2f^*(\Delta)=270.5+36\Delta+1.2\Delta^2,小 Δ\Delta 时二次项可忽略。

例7 容器设计

题目。 设计有盖长方体容器装 98 cm³ 化学粉末,长 x1x_1 与高 x2x_2 之比必须为 2.5,宽为 x3x_3,材料 10 元/cm²。(a) 建立成本最小的模型;(b) 用拉格朗日乘子求尺寸和最低成本。

(a)

min⁡ C=20(x1x2+x1x3+x2x3)s.t.x1x2x3=98,x1=2.5x2,x1,x2,x3>0.\min\ C=20(x_1x_2+x_1x_3+x_2x_3)\quad\text{s.t.}\quad x_1x_2x_3=98,\quad x_1=2.5x_2,\quad x_1,x_2,x_3>0 .

(b) 用 f=x1x2+x1x3+x2x3f=x_1x_2+x_1x_3+x_2x_3(成本 =20f=20f)。L=f+λ1(x1x2x3−98)+λ2(x1−2.5x2)L=f+\lambda_1(x_1x_2x_3-98)+\lambda_2(x_1-2.5x_2):

x2+x3+λ1x2x3+λ2=0x1+x3+λ1x1x3−2.5λ2=0x1+x2+λ1x1x2=0\begin{aligned} x_2+x_3+\lambda_1x_2x_3+\lambda_2&=0\\ x_1+x_3+\lambda_1x_1x_3-2.5\lambda_2&=0\\ x_1+x_2+\lambda_1x_1x_2&=0 \end{aligned}

令 x2=tx_2=t,x1=2.5tx_1=2.5t。第三式:3.5t+2.5λ1t2=0⇒λ1=−1.4t3.5t+2.5\lambda_1t^2=0\Rightarrow\lambda_1=-\tfrac{1.4}{t}。第一式乘 2.5 加第二式消去 λ2\lambda_2:

2.5t+2.5x3−3.5x3+2.5t+x3−3.5x3=0 ⇒ 5t=3.5x3 ⇒ x3=107t.2.5t+2.5x_3-3.5x_3+2.5t+x_3-3.5x_3=0\ \Rightarrow\ 5t=3.5x_3\ \Rightarrow\ x_3=\tfrac{10}{7}t .

容积:2.5t⋅t⋅107t=257t3=982.5t\cdot t\cdot\tfrac{10}{7}t=\tfrac{25}{7}t^3=98,t3=27.44t^3=27.44,t=27.443≈3.016t=\sqrt[3]{27.44}\approx3.016。

尺寸:长 x1≈7.54x_1\approx7.54 cm,高 x2≈3.02x_2\approx3.02 cm,宽 x3≈4.31x_3\approx4.31 cm。表面积 2f≈136.462f\approx136.46 cm²,最低成本约 1364.6 元。

要点

先用比例约束把变量减到一个参数 tt,再用不含 λ2\lambda_2 的那个方程求出 λ1\lambda_1,最后用另外两个方程消去 λ2\lambda_2。数字不整齐时保留根式或三位小数即可。

例8 离散时间马尔可夫链:只在售罄时补货

题目。 店里每次进 3 块豪华手表,只有全部卖完才在下月初再进 3 块。每月需求 DkD_k 独立,服从 λ=1\lambda=1 的泊松分布;LkL_k 为第 kk 月底剩余块数,L0=3L_0=3。(a) 写出状态空间和转移概率矩阵(用 ee 表示);(b) 求各状态的平均逗留时间(3 位小数),哪个状态停留最短;(c) 若每次进 mm 块(m>3m>3),写出状态空间和转移概率。

(a) 状态空间 {0,1,2,3}\{0,1,2,3\}。月初库存:月底为 0 → 3;为 1 → 1;为 2 → 2;为 3 → 3。用“从库存 ss 卖到剩 jj”:j≥1j\ge1 时 P(D=s−j)P(D=s-j),j=0j=0 时 P(D≥s)P(D\ge s)。

P=(1−52e−112e−1e−1e−11−e−1e−1001−2e−1e−1e−101−52e−112e−1e−1e−1).P=\begin{pmatrix} 1-\frac52e^{-1}&\frac12e^{-1}&e^{-1}&e^{-1}\\ 1-e^{-1}&e^{-1}&0&0\\ 1-2e^{-1}&e^{-1}&e^{-1}&0\\ 1-\frac52e^{-1}&\frac12e^{-1}&e^{-1}&e^{-1} \end{pmatrix}.

(b) E(Ti)=1/(1−pii)E(T_i)=1/(1-p_{ii}):

状态 0 1 2 3
piip_{ii} 1−52e−1=0.0801-\frac52e^{-1}=0.080 e−1=0.368e^{-1}=0.368 e−1=0.368e^{-1}=0.368 e−1=0.368e^{-1}=0.368
平均逗留(月) 1.087 1.582 1.582 1.582

状态 0(售罄)停留最短:一旦售罄就会补货,下个月几乎总有剩余。

(c) 状态空间 {0,1,…,m}\{0,1,\dots,m\}。i≥1i\ge1 时月初库存为 ii:

pij=e−1(i−j)! (1≤j≤i),pi0=1−∑k=0i−1e−1k!,pij=0 (j>i).p_{ij}=\frac{e^{-1}}{(i-j)!}\ (1\le j\le i),\qquad p_{i0}=1-\sum_{k=0}^{i-1}\frac{e^{-1}}{k!},\qquad p_{ij}=0\ (j>i).

i=0i=0 时月初库存为 mm,第 0 行与第 mm 行相同:p0j=e−1(m−j)! (1≤j≤m)p_{0j}=\dfrac{e^{-1}}{(m-j)!}\ (1\le j\le m),p00=1−∑k=0m−1e−1k!p_{00}=1-\sum_{k=0}^{m-1}\dfrac{e^{-1}}{k!}。

例9 连续时间马尔可夫链:工作流程的速率

题目。 化妆师的工作:去客户家(速率 2/小时)→ 给新娘化妆(速率 3)→ 给新娘母亲化妆(速率 3)→ 去下一位客户。每次化妆中眼线可能花掉(速率 0.5),需要重画(速率 1),重画后继续原来的化妆。各时间独立且服从指数分布。状态:0 路上,1 新娘化妆,2 母亲化妆,3 新娘重画,4 母亲重画。求稳态概率和平均每小时完成的化妆数,并与只给新娘化妆的情况比较。

转移速率:0→10\to1:2;1→21\to2:3,1→31\to3:0.5;3→13\to1:1;2→02\to0:3,2→42\to4:0.5;4→24\to2:1。

速率平衡:

状态 0: 2π0=3π2状态 3: π3=0.5π1状态 4: π4=0.5π2状态 1: 3.5π1=2π0+π3 ⇒ 3π1=2π0\begin{aligned} \text{状态 0:}&\ 2\pi_0=3\pi_2\\ \text{状态 3:}&\ \pi_3=0.5\pi_1\\ \text{状态 4:}&\ \pi_4=0.5\pi_2\\ \text{状态 1:}&\ 3.5\pi_1=2\pi_0+\pi_3\ \Rightarrow\ 3\pi_1=2\pi_0\\ \end{aligned}

于是 π1=π2=23π0\pi_1=\pi_2=\tfrac23\pi_0,π3=π4=13π0\pi_3=\pi_4=\tfrac13\pi_0,π0(1+23+23+13+13)=1\pi_0(1+\tfrac23+\tfrac23+\tfrac13+\tfrac13)=1,π0=13\pi_0=\tfrac13:

π=(13, 29, 29, 19, 19).\pi=\Big(\tfrac13,\ \tfrac29,\ \tfrac29,\ \tfrac19,\ \tfrac19\Big).

平均化妆速率 =3π1+3π2=43=3\pi_1+3\pi_2=\tfrac43 次/小时。只给新娘化妆时(三状态:路上、化妆、重画),稳态为 (12,13,16)(\tfrac12,\tfrac13,\tfrac16),化妆速率 3×13=13\times\tfrac13=1 次/小时。一次出行做两个人,路上的时间被分摊,化妆速率提高到 4/34/3(但每小时服务的客户家庭数从 1 降到 2/32/3)。

例10 排队:三种设计与汽修店

题目一。 食堂卖三类饭菜,到达为泊松过程、服务为指数分布。系统 1:三个窗口各卖一类,各自排队,每队到达率 λ/3\lambda/3、服务率 μ\mu;系统 2:一个窗口卖全部,到达率 λ\lambda、服务率 3μ3\mu;系统 3:三个窗口都能卖全部,共用一条队(取号),到达率 λ\lambda、每窗口服务率 μ\mu。(a) 写出模型和 WiW_i,稳定条件;(b) λ=μ=1\lambda=\mu=1 时排序;(c) 一般情况下排序。

(a) 系统 1:三个 M/M/1,W1=1μ−λ/3=33μ−λW_1=\dfrac{1}{\mu-\lambda/3}=\dfrac{3}{3\mu-\lambda}。系统 2:M/M/1,W2=13μ−λW_2=\dfrac{1}{3\mu-\lambda}。系统 3:M/M/3,代入第9章公式化简得

W3=18μ2+6λμ−λ2(3μ−λ)(λ2+4λμ+6μ2).W_3=\frac{18\mu^2+6\lambda\mu-\lambda^2}{(3\mu-\lambda)(\lambda^2+4\lambda\mu+6\mu^2)} .

三者稳定条件都是 0<λ<3μ0<\lambda<3\mu(μ>0\mu>0)。

(b) λ=μ=1\lambda=\mu=1:W1=32=1.5W_1=\tfrac32=1.5,W2=12=0.5W_2=\tfrac12=0.5,W3=18+6−12(1+4+6)=2322≈1.045W_3=\dfrac{18+6-1}{2(1+4+6)}=\tfrac{23}{22}\approx1.045。排序 W2<W3<W1W_2<W_3<W_1。

(c) 两个比值的分母在稳定范围内都为正。比较分子与分母之差:

W3W2−1=2(3μ−λ)(2μ+λ)λ2+4λμ+6μ2>0,\frac{W_3}{W_2}-1=\frac{2(3\mu-\lambda)(2\mu+\lambda)}{\lambda^2+4\lambda\mu+6\mu^2}>0, W1W3−1=4λ2+6λμ18μ2+6λμ−λ2>0.\frac{W_1}{W_3}-1=\frac{4\lambda^2+6\lambda\mu}{18\mu^2+6\lambda\mu-\lambda^2}>0 .

所以 W2<W3<W1W_2<W_3<W_1:一个快窗口的总停留时间最短,分开排队最长。这里的比较依赖“合力服务速率确实达到 3μ3\mu”这一假设。

题目二。 汽修店每天来 10 辆车(泊松),3 名技师。A:每人独立修一辆,速率 6/天,共用一条队;B:三人合修一辆,速率 18/天。比较 L、Q、W、D 并给出建议。

A 是 M/M/3(ρ=5/9\rho=5/9),B 是 M/M/1(ρ=5/9\rho=5/9):

L Q W(天) D(天)
A:M/M/3 2.04 0.37 0.204 0.037
B:M/M/1 1.25 0.69 0.125 0.069

B 的总停留时间短(合修更快),A 的排队等候短(更少需要排队)。车主最在意“多久能取车”,即 WW,所以建议 B——前提是三人合作确实能把速度提高到 3 倍。

例11 决策树:最终资产及其概率

题目。 钻石贸易公司现有资产 2000(千元),一批原石要卖给客户,运费 50。不切割直接运:53% 被接受,利润 300(扣成本前)。切割(成本 75):切得好的概率 70%,此时接受率 80%;切得不好时接受率只有 12%,但可以再花 25 抛光,抛光后接受率 90%。被拒收时只承担已发生的成本。求使期望最终资产最大的策略、可能的最终资产及概率、最可能的最终资产。

回推:

  • 不切割:0.53×(2000+300−50)+0.47×(2000−50)=0.53×2250+0.47×1950=21090.53\times(2000+300-50)+0.47\times(2000-50)=0.53\times2250+0.47\times1950=2109。
  • 切割、切得好:0.8×(2300−75−50)+0.2×(2000−75−50)=0.8×2175+0.2×1875=21150.8\times(2300-75-50)+0.2\times(2000-75-50)=0.8\times2175+0.2\times1875=2115。
  • 切割、切得不好:直接运 0.12×2175+0.88×1875=19110.12\times2175+0.88\times1875=1911;抛光 0.9×(2175−25)+0.1×(1875−25)=0.9×2150+0.1×1850=21200.9\times(2175-25)+0.1\times(1875-25)=0.9\times2150+0.1\times1850=2120。选抛光。
  • 切割的期望:0.7×2115+0.3×2120=2116.5>21090.7\times2115+0.3\times2120=2116.5>2109。

最优策略:切割;切得不好就抛光。 期望最终资产 2116.5。

最终资产 路径 概率
2175 切得好 → 接受 0.7×0.8=0.560.7\times0.8=0.56
1875 切得好 → 拒收 0.7×0.2=0.140.7\times0.2=0.14
2150 切得不好 → 抛光 → 接受 0.3×0.9=0.270.3\times0.9=0.27
1850 切得不好 → 抛光 → 拒收 0.3×0.1=0.030.3\times0.1=0.03

最可能的最终资产是 2175(概率 0.56)。

例12 效用与 Hurwicz 准则

题目一。 某人有 10000 元,要买机票去开会。廉价票 300 元,有 5% 的概率严重延误,错过会议造成额外损失 2500 元;灵活票 300+p300+p 元,保证不会错过会议。效用函数 u(x)=x/10000u(x)=\sqrt{x/10000},xx 为最终财富。他最多愿付的升级费 pp 是多少?

E[u廉价]=0.950.97+0.050.72=0.95×0.98489+0.05×0.84853=0.97807.E[u_{\text{廉价}}]=0.95\sqrt{0.97}+0.05\sqrt{0.72}=0.95\times0.98489+0.05\times0.84853=0.97807 . 9700−p10000=0.97807 ⇒ 9700−p=9566.2 ⇒ p≈133.8.\sqrt{\frac{9700-p}{10000}}=0.97807\ \Rightarrow\ 9700-p=9566.2\ \Rightarrow\ p\approx133.8 .

最多愿付约 133.8 元,高于期望损失 0.05×2500=1250.05\times2500=125 元,差额 8.8 元是风险溢价。

题目二。 四种生意在三种发展情景下的回报率(%):便利店 (8, 12, 20),理发店 (7, 11, 21),咖啡店 (9, 13, 25),补习中心 (6, 12, 30)。用悲观、乐观、Laplace、后悔和 Hurwicz 准则决策。

准则 计算 选择
悲观 最小值 8, 7, 9, 6 咖啡店
乐观 最大值 20, 21, 25, 30 补习中心
Laplace 平均 13.3, 13, 15.7, 16 补习中心
后悔 列最好值 9, 13, 30;最大后悔 10, 9, 5, 3 补习中心

Hurwicz:H便利店=8+12KH_{\text{便利店}}=8+12K,H理发店=7+14KH_{\text{理发店}}=7+14K,H咖啡店=9+16KH_{\text{咖啡店}}=9+16K,H补习=6+24KH_{\text{补习}}=6+24K。咖啡店在 [0,1][0,1] 上始终高于便利店和理发店,只需比较咖啡店与补习中心:9+16K=6+24K⇒K=3/89+16K=6+24K\Rightarrow K=3/8。

KK 选择
0≤K<3/80\le K<3/8 咖啡店
K=3/8K=3/8 两者并列
3/8<K≤13/8<K\le1 补习中心

总结:拿到题先想三件事

  1. 这是哪一类模型? 线性目标 + 线性约束 → 线性规划;供需表 → 运输;一人一事 → 指派;有平方、乘积 → 非线性规划;“每期”“每天” + 随机需求 → 离散时间马尔可夫链;指数时间、速率 → 连续时间马尔可夫链或排队;收益表 + 情景 → 决策分析。
  2. 状态、变量、约束写全了吗? 隐含的非负、平衡、容量限制最容易漏。
  3. 结果合理吗? 代回约束、检查概率和为 1、用 Little 公式交叉验证、看最优值是否随参数朝合理方向变化。