系统分析:线性规划、非线性规划与排队模型 · 第 11 / 11 篇(展开目录)
从线性规划建模、单纯形法、对偶与灵敏度、运输与指派,到非线性规划、随机过程、马尔可夫链、排队模型和决策分析,每章都配手算例题和自测,最后一章用综合例题把方法串起来。
- 第1章 线性规划建模:从文字题到标准型
- 第2章 单纯形法:单纯形表、大 M 法与两阶段法
- 第3章 对偶与灵敏度分析:最终单纯形表还能回答什么
- 第4章 运输问题与指派问题:Vogel 法、位势法与匈牙利法
- 第5章 非线性规划(一):一维搜索与多元无约束优化
- 第6章 非线性规划(二):拉格朗日乘子与 KKT 条件
- 第7章 随机过程基础:几何分布、指数分布与泊松过程
- 第8章 马尔可夫链:离散时间与连续时间
- 第9章 排队模型:从 M/M/1 到 M/M/m
- 第10章 决策分析:不确定性准则、决策树与效用
- 第11章 综合例题讲解与解析
约 60 分钟 · 需要前十章的方法
怎么用这一章: 先只看题目,自己写出解题框架(建什么模型、用哪一章的方法),再对照解析。提示框会提醒建模和计算中容易忽略的细节。
例1 线性规划:建模、求解与预算削减
题目。 冰淇淋店做三种口味,每杯的配料成本、人工成本和包装材料如下;每天包装材料最多 10 kg,配料支出不超过 100 元,人工成本不超过 30 元。三种口味每杯利润分别为 2.5、2.3、2.0 元。(a) 建立利润最大的模型;(b) 每天各做多少杯;(c) 老板想削减人工预算,最多能削减多少,而(b)的方案仍可行?
| 口味 |
配料(元/杯) |
人工(元/杯) |
包装(kg/杯) |
| 香草 |
1.50 |
0.50 |
0.20 |
| 巧克力 |
1.70 |
0.35 |
0.18 |
| 草莓 |
1.80 |
0.60 |
0.16 |
(a)
(b) 加松弛变量 ,化为最小化 。
表 0:检验行 。 进基,比值 、、, 出基。
表 1:
| 基 |
|
|
|
|
|
|
右端 |
|
1 |
9/10 |
4/5 |
5 |
0 |
0 |
50 |
|
0 |
7/20 |
3/5 |
−15/2 |
1 |
0 |
25 |
|
0 |
−1/10 |
1/5 |
−5/2 |
0 |
1 |
5 |
| 检验行 |
0 |
−1/20 |
0 |
25/2 |
0 |
0 |
125 |
的检验数 −1/20 仍为负, 进基;比值 、, 出基。
表 2:
| 基 |
|
|
|
|
|
|
右端 |
|
10/9 |
1 |
8/9 |
50/9 |
0 |
0 |
500/9 |
|
−7/18 |
0 |
13/45 |
−85/9 |
1 |
0 |
50/9 |
|
1/9 |
0 |
13/45 |
−35/18 |
0 |
1 |
95/9 |
| 检验行 |
1/18 |
0 |
2/45 |
115/9 |
0 |
0 |
1150/9 |
检验行全部非负。线性规划最优解:只做巧克力 杯,利润 元。只有包装材料是紧约束(影子价格 元/kg)。
上述是允许产量连续变化的线性规划解。实际按整杯生产时,需要另加整数约束:简单取整为 55 杯巧克力,利润为 126.5 元;精确的整数最优方案是 52 杯巧克力 + 4 杯草莓,利润 127.6 元(包装恰好 kg)。直接取整不保证整数最优。
(c) 对连续产量的最优解,人工成本只用了 元,人工约束的松弛 。人工预算从 30 降到 之前,当前解都可行,最多削减约 10.56 元。用第3章的语言:右端变化 时 变成 ,要求 ,即 。若采用整杯最优方案 ,人工成本为 元,可削减的预算相应为 9.4 元。
要点
非紧约束能削减的量就是它的松弛量,直接从最终表右端读。“紧”还是“松”要看最终表的基变量: 在基中且为正,说明人工预算没用完。
例2 含等式约束的单纯形法与目标系数灵敏度
题目。 (a) 用单纯形法求 ,s.t. ,,。(b) 目标变为 时, 在什么范围内(a)的解仍最优?(c) 若问题是最大化型,第一阶段中人工变量的目标系数应如何取?
(a) ,,初始基 。检验行拆成常数行和 M 行:
| 基 |
|
|
|
|
右端 |
|
3 |
−2 |
1 |
0 |
3 |
|
8 |
4 |
0 |
1 |
5 |
| 常数行 |
2 |
3 |
0 |
0 |
0 |
| M 行 |
−8 |
−4 |
0 |
0 |
−5 |
M 行最负 −8, 进基;比值 、, 出基:
| 基 |
|
|
|
右端 |
|
0 |
−7/2 |
1 |
9/8 |
|
1 |
1/2 |
0 |
5/8 |
| 常数行 |
0 |
2 |
0 |
−5/4 |
人工变量已出基(删去该列和 M 行),常数行全部非负,最优:**,**。
(b) 是非基变量,只影响它自己的检验数:,即 ****()时解不变; 增大不会改变最优解。
(c) 最大化问题中人工变量取 (目标越大越好,人工变量越大越被惩罚);若先把目标乘 化为最小化,则取 。第一阶段等价于最小化人工变量之和,与原问题是最大还是最小无关。
例3 运输问题与运价灵敏度
题目。 三个工厂向三个地区供货(箱),单位运费如下。(a) 写成平衡运输问题;(b) 用 Vogel 法求初始解并求最优解;(c) 工厂 3 到地区 3 的运费变化 ,求当前解保持最优的范围。
|
地区 1 |
地区 2 |
地区 3 |
供应 |
| 工厂 1 |
20 |
24 |
24 |
60 |
| 工厂 2 |
18 |
23 |
20 |
70 |
| 工厂 3 |
20 |
25 |
26 |
50 |
| 需求 |
40 |
65 |
75 |
|
(a) 总供应 180 = 总需求 180,已平衡:,,,。
(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 个 。。
位势():。非基格检验数:格(1,1) ,(2,1) ,(2,2) ,(3,3) ,全部为正。Vogel 初始解已经最优,最小运费 3890。
(c) 格(3,3) 是非基格,位势不变,只需 :****,即运费不低于 25 时方案不变(低于 25 改走这条线路更便宜)。
例4 生产—库存问题
题目。 某器材 11 月—2 月的需求为 100、200、180、300 套,各月产能 50、180、280、270 套,生产成本 40 元/套。提前生产要付保管费,每套每月 0.5 元;延期交货要付罚金,每套每月 2 元。求最优生产计划。
产能合计 780 = 需求合计 780,平衡。运价 = 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 个基格。位势检验所有非基格检验数 (格(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 |
“—”记为 ,加虚拟人 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 无人承担。总成本 。
例6 拉格朗日乘子:两个等式约束与右端扰动
题目一。 ,s.t. ,。
解出 用 表示后代入两条约束:
得 ,,于是
正定,是极小点(也是全局最小:凸目标 + 线性约束)。
题目二。 ,s.t. ,。(a) 求最优解;(b) 第一条约束改为 ,最优值约变化多少?
(a) ,,。由前两式 ,代入 得 ,;。。
(b) 。精确值 ,小 时二次项可忽略。
例7 容器设计
题目。 设计有盖长方体容器装 98 cm³ 化学粉末,长 与高 之比必须为 2.5,宽为 ,材料 10 元/cm²。(a) 建立成本最小的模型;(b) 用拉格朗日乘子求尺寸和最低成本。
(a)
(b) 用 (成本 )。:
令 ,。第三式:。第一式乘 2.5 加第二式消去 :
容积:,,。
尺寸:长 cm,高 cm,宽 cm。表面积 cm²,最低成本约 1364.6 元。
要点
先用比例约束把变量减到一个参数 ,再用不含 的那个方程求出 ,最后用另外两个方程消去 。数字不整齐时保留根式或三位小数即可。
例8 离散时间马尔可夫链:只在售罄时补货
题目。 店里每次进 3 块豪华手表,只有全部卖完才在下月初再进 3 块。每月需求 独立,服从 的泊松分布; 为第 月底剩余块数,。(a) 写出状态空间和转移概率矩阵(用 表示);(b) 求各状态的平均逗留时间(3 位小数),哪个状态停留最短;(c) 若每次进 块(),写出状态空间和转移概率。
(a) 状态空间 。月初库存:月底为 0 → 3;为 1 → 1;为 2 → 2;为 3 → 3。用“从库存 卖到剩 ”: 时 , 时 。
(b) :
| 状态 |
0 |
1 |
2 |
3 |
|
|
|
|
|
| 平均逗留(月) |
1.087 |
1.582 |
1.582 |
1.582 |
状态 0(售罄)停留最短:一旦售罄就会补货,下个月几乎总有剩余。
(c) 状态空间 。 时月初库存为 :
时月初库存为 ,第 0 行与第 行相同:,。
例9 连续时间马尔可夫链:工作流程的速率
题目。 化妆师的工作:去客户家(速率 2/小时)→ 给新娘化妆(速率 3)→ 给新娘母亲化妆(速率 3)→ 去下一位客户。每次化妆中眼线可能花掉(速率 0.5),需要重画(速率 1),重画后继续原来的化妆。各时间独立且服从指数分布。状态:0 路上,1 新娘化妆,2 母亲化妆,3 新娘重画,4 母亲重画。求稳态概率和平均每小时完成的化妆数,并与只给新娘化妆的情况比较。
转移速率::2;:3,:0.5;:1;:3,:0.5;:1。
速率平衡:
于是 ,,,:
平均化妆速率 次/小时。只给新娘化妆时(三状态:路上、化妆、重画),稳态为 ,化妆速率 次/小时。一次出行做两个人,路上的时间被分摊,化妆速率提高到 (但每小时服务的客户家庭数从 1 降到 )。
例10 排队:三种设计与汽修店
题目一。 食堂卖三类饭菜,到达为泊松过程、服务为指数分布。系统 1:三个窗口各卖一类,各自排队,每队到达率 、服务率 ;系统 2:一个窗口卖全部,到达率 、服务率 ;系统 3:三个窗口都能卖全部,共用一条队(取号),到达率 、每窗口服务率 。(a) 写出模型和 ,稳定条件;(b) 时排序;(c) 一般情况下排序。
(a) 系统 1:三个 M/M/1,。系统 2:M/M/1,。系统 3:M/M/3,代入第9章公式化简得
三者稳定条件都是 ()。
(b) :,,。排序 。
(c) 两个比值的分母在稳定范围内都为正。比较分子与分母之差:
所以 :一个快窗口的总停留时间最短,分开排队最长。这里的比较依赖“合力服务速率确实达到 ”这一假设。
题目二。 汽修店每天来 10 辆车(泊松),3 名技师。A:每人独立修一辆,速率 6/天,共用一条队;B:三人合修一辆,速率 18/天。比较 L、Q、W、D 并给出建议。
A 是 M/M/3(),B 是 M/M/1():
|
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 的排队等候短(更少需要排队)。车主最在意“多久能取车”,即 ,所以建议 B——前提是三人合作确实能把速度提高到 3 倍。
例11 决策树:最终资产及其概率
题目。 钻石贸易公司现有资产 2000(千元),一批原石要卖给客户,运费 50。不切割直接运:53% 被接受,利润 300(扣成本前)。切割(成本 75):切得好的概率 70%,此时接受率 80%;切得不好时接受率只有 12%,但可以再花 25 抛光,抛光后接受率 90%。被拒收时只承担已发生的成本。求使期望最终资产最大的策略、可能的最终资产及概率、最可能的最终资产。
回推:
- 不切割:。
- 切割、切得好:。
- 切割、切得不好:直接运 ;抛光 。选抛光。
- 切割的期望:。
最优策略:切割;切得不好就抛光。 期望最终资产 2116.5。
| 最终资产 |
路径 |
概率 |
| 2175 |
切得好 → 接受 |
|
| 1875 |
切得好 → 拒收 |
|
| 2150 |
切得不好 → 抛光 → 接受 |
|
| 1850 |
切得不好 → 抛光 → 拒收 |
|
最可能的最终资产是 2175(概率 0.56)。
例12 效用与 Hurwicz 准则
题目一。 某人有 10000 元,要买机票去开会。廉价票 300 元,有 5% 的概率严重延误,错过会议造成额外损失 2500 元;灵活票 元,保证不会错过会议。效用函数 , 为最终财富。他最多愿付的升级费 是多少?
最多愿付约 133.8 元,高于期望损失 元,差额 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:,,,。咖啡店在 上始终高于便利店和理发店,只需比较咖啡店与补习中心:。
|
选择 |
|
咖啡店 |
|
两者并列 |
|
补习中心 |
总结:拿到题先想三件事
- 这是哪一类模型? 线性目标 + 线性约束 → 线性规划;供需表 → 运输;一人一事 → 指派;有平方、乘积 → 非线性规划;“每期”“每天” + 随机需求 → 离散时间马尔可夫链;指数时间、速率 → 连续时间马尔可夫链或排队;收益表 + 情景 → 决策分析。
- 状态、变量、约束写全了吗? 隐含的非负、平衡、容量限制最容易漏。
- 结果合理吗? 代回约束、检查概率和为 1、用 Little 公式交叉验证、看最优值是否随参数朝合理方向变化。