系统分析:线性规划、非线性规划与排队模型 · 第 7 / 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章 综合例题讲解与解析

约 25 分钟 · 需要条件概率

读完你能: 用几何分布和指数分布描述寿命、加工时间、修理时间;用无记忆性计算条件概率;求两个独立指数随机变量谁先发生、最小值的分布;区分四类随机过程;计算泊松过程在一段时间内到达 kk 个的概率,并处理泊松流的合并与分流。

1. 几何分布:第几次才成功

每次试验独立,成功概率都是 pp。XX = 第一次成功发生在第几次:

P(X=n)=(1−p)n−1p,n=1,2,… ;P(X>n)=(1−p)n;E(X)=1p.P(X=n)=(1-p)^{n-1}p,\qquad n=1,2,\dots;\qquad P(X>n)=(1-p)^n;\qquad E(X)=\frac1p .

P(X>n)P(X>n) 的含义是“前 nn 次都失败”,所以是 (1−p)n(1-p)^n。在系统分析里,“成功”常常是某个事件发生:一个时间单位里工件加工完、刀具损坏、故障修好……

例1: 刀具每次加工有 2% 的概率损坏,XX 是它坚持到第几次加工时损坏。

  • 平均寿命 E(X)=1/0.02=50E(X)=1/0.02=50 次加工。
  • 至少还能用 20 次以上:P(X>20)=0.9820=0.6676P(X>20)=0.98^{20}=0.6676。

无记忆性

P(X>m+n∣X>m)=P(X>n),P(X=m+n∣X>m)=P(X=n).P(X>m+n\mid X>m)=P(X>n),\qquad P(X=m+n\mid X>m)=P(X=n).

证明只需一行:P(X>m+n∣X>m)=(1−p)m+n(1−p)m=(1−p)nP(X>m+n\mid X>m)=\dfrac{(1-p)^{m+n}}{(1-p)^m}=(1-p)^n。

例1 续: 已经用了 10 次还没坏,再用 30 次以上不坏的概率是多少?P(X>40∣X>10)=P(X>30)=0.9830=0.5455P(X>40\mid X>10)=P(X>30)=0.98^{30}=0.5455。已经用过的 10 次不影响将来——“旧刀”和“新刀”一样。

几何分布是唯一具有无记忆性的离散分布。

2. 指数分布:连续时间里的“等多久”

f(x)=λe−λx,P(X>x)=e−λx,P(X≤x)=1−e−λx,E(X)=1λ(x≥0).f(x)=\lambda e^{-\lambda x},\qquad P(X>x)=e^{-\lambda x},\qquad P(X\le x)=1-e^{-\lambda x},\qquad E(X)=\frac1\lambda\qquad(x\ge0).

λ\lambda 叫速率(每单位时间发生多少次),1/λ1/\lambda 是平均时间。它常用来描述到达间隔、加工时间、准备时间、故障前时间、修理时间等。

无记忆性:

P(X>x+y∣X>x)=P(X>y),P(X≤x+y∣X>x)=P(X≤y).P(X>x+y\mid X>x)=P(X>y),\qquad P(X\le x+y\mid X>x)=P(X\le y).

指数分布是唯一具有无记忆性的连续分布。这一点非常关键:正是因为它,连续时间马尔可夫链在某个状态里“已经待了多久”不影响将来,才能只用当前状态描述系统。

例2: 修理时间服从指数分布,平均 20 分钟(λ=3\lambda=3 次/小时)。

  • 超过 30 分钟的概率:P(X>0.5)=e−3×0.5=e−1.5=0.2231P(X>0.5)=e^{-3\times0.5}=e^{-1.5}=0.2231。
  • 已经修了 20 分钟,再过 10 分钟内修好的概率:P(X≤30∣X>20)=P(X≤10 分钟)=1−e−0.5=0.3935P(X\le30\mid X>20)=P(X\le10\text{ 分钟})=1-e^{-0.5}=0.3935。

3. 两个指数时间“赛跑”

设 X1∼EXP(λ1)X_1\sim\mathrm{EXP}(\lambda_1)、X2∼EXP(λ2)X_2\sim\mathrm{EXP}(\lambda_2) 相互独立。

P(X1<X2)=λ1λ1+λ2,min⁡{X1,X2}∼EXP(λ1+λ2).P(X_1<X_2)=\frac{\lambda_1}{\lambda_1+\lambda_2},\qquad \min\{X_1,X_2\}\sim\mathrm{EXP}(\lambda_1+\lambda_2).

第二个结论:P(min⁡>x)=P(X1>x)P(X2>x)=e−(λ1+λ2)xP(\min>x)=P(X_1>x)P(X_2>x)=e^{-(\lambda_1+\lambda_2)x}。第一个结论:

P(X1<X2)=∫0∞λ1e−λ1x e−λ2x dx=λ1λ1+λ2.P(X_1<X_2)=\int_0^\infty\lambda_1e^{-\lambda_1x}\,e^{-\lambda_2x}\,dx=\frac{\lambda_1}{\lambda_1+\lambda_2}.

推广到多个:nn 个独立指数时间的最小值速率是 ∑λi\sum\lambda_i,第 ii 个最先发生的概率是 λi/∑λj\lambda_i/\sum\lambda_j。

例3: 机器加工一个零件的时间 EXP(0.4)\mathrm{EXP}(0.4)(每小时),故障前时间 EXP(0.1)\mathrm{EXP}(0.1)。

  • 加工过程中“先完成还是先故障”:零件先完成的概率 0.4/(0.4+0.1)=0.80.4/(0.4+0.1)=0.8。
  • 机器在“加工”状态停留的时间:min⁡∼EXP(0.5)\min\sim\mathrm{EXP}(0.5),平均 2 小时——不论以哪种方式离开。

这正是第8章连续时间马尔可夫链里“一个状态有几个出口”的计算方法。

几何分布的对应结论

离散情形也有类似结论:独立几何随机变量 G1,G2G_1,G_2(参数 p1,p2p_1,p_2)中 G1G_1 严格先发生的概率是 p1(1−p2)p1+p2−p1p2\dfrac{p_1(1-p_2)}{p_1+p_2-p_1p_2}。离散时间里两者可能在同一步发生,所以公式比连续情形复杂。

4. 随机过程的四种类型

随机过程 {X(t):t∈T}\{X(t):t\in T\} 是一族随机变量,tt 通常是时间。所有可能取值组成状态空间 SS。按时间和状态是离散还是连续分成四类:

状态离散(称为“链”) 状态连续
时间离散 每小时检查一次机器是好是坏 每天中午测量的温度
时间连续 任意时刻机器的状态;[0,t][0,t] 内到达的零件数 机械臂在 tt 时刻走过的距离

接下来讨论两类“链”:离散时间马尔可夫链(DTMC)和连续时间马尔可夫链(CTMC),见第8章;以及一种特殊的连续时间链——泊松过程。

5. 泊松过程

X(t)X(t) 表示 [0,t][0,t] 内到达的个数。齐次泊松过程从 X(0)=0X(0)=0 开始,不相交时间区间内的增量相互独立,长度为 tt 的区间内到达数服从同一个均值为 λt\lambda t 的泊松分布:

P(X(t)=k)=e−λt(λt)kk!,k=0,1,2,…,P(X(t)=k)=\frac{e^{-\lambda t}(\lambda t)^k}{k!},\qquad k=0,1,2,\dots,

此时 {X(t):t≥0}\{X(t):t\ge0\} 是速率为 λ\lambda 的泊松过程,E[X(t)]=λtE[X(t)]=\lambda t。它描述零件到达车间、任务到达计算中心、电话到达交换机、顾客到达柜台。只有单个时刻的泊松分布还不够,独立增量与平稳增量同样是定义的一部分。

令 k=0k=0:P([0,t] 内没有到达)=e−λtP([0,t]\ \text{内没有到达})=e^{-\lambda t}。这正是“下一次到达的等待时间 TT 超过 tt”的概率,所以

P(T>t)=e−λt⇒T∼EXP(λ).P(T>t)=e^{-\lambda t}\quad\Rightarrow\quad T\sim\mathrm{EXP}(\lambda).

两句话记住泊松过程:

  • 到达间隔服从速率 λ\lambda 的指数分布;
  • 长度为 tt 的时间内到达个数服从均值 λt\lambda t 的泊松分布。

例4: 顾客按每小时 6 人的泊松过程到达。

  • 10 分钟内恰好来 2 人:λt=6×16=1\lambda t=6\times\frac16=1,P=e−1⋅12/2!=0.1839P=e^{-1}\cdot1^2/2!=0.1839。
  • 15 分钟内没人来:e−6×0.25=e−1.5=0.2231e^{-6\times0.25}=e^{-1.5}=0.2231。
  • 30 分钟内至少来 2 人:λt=3\lambda t=3,1−e−3(1+3)=0.80091-e^{-3}(1+3)=0.8009。
  • 相邻两位顾客平均间隔 1/61/6 小时 = 10 分钟。

6. 泊松流的合并与分流

合并(叠加): nn 个独立的泊松流,速率 λ1,…,λn\lambda_1,\dots,\lambda_n,合在一起仍是泊松流,速率 ∑λi\sum\lambda_i。例如两条生产线分别以每小时 4 件、2 件送来零件,检验台收到的是每小时 6 件的泊松流。

分流(分解): 速率 λ\lambda 的泊松流中,每个到达独立地以概率 pip_i 分到第 ii 类(∑pi=1\sum p_i=1),则第 ii 类是速率 piλp_i\lambda 的泊松流,且各类之间相互独立。例如每小时 10 位顾客,30% 买 A 套餐,则 A 套餐的点单是每小时 3 单的泊松流。

排队比较题的关键一步

“三个窗口各卖一种饭,顾客按需求分流”的系统,就是把速率 λ\lambda 的总流分成三条速率 λ/3\lambda/3 的泊松流,每个窗口各自是一个单服务台排队。第9章比较几种排队设计时会反复用到这一点。

7. 小结

几何分布 指数分布 泊松分布
类型 离散,取 1,2,…1,2,\dots 连续,取 x≥0x\ge0 离散,取 0,1,2,…0,1,2,\dots
含义 第几次才发生 要等多久 一段时间内发生几次
尾概率/概率 P(X>n)=(1−p)nP(X>n)=(1-p)^n P(X>x)=e−λxP(X>x)=e^{-\lambda x} P(k)=e−λt(λt)k/k!P(k)=e^{-\lambda t}(\lambda t)^k/k!
均值 1/p1/p 1/λ1/\lambda λt\lambda t
无记忆性 有(离散中唯一) 有(连续中唯一) —
  • 独立指数时间的最小值:速率相加;谁先发生:速率占比。
  • 泊松过程 = 到达间隔为指数分布的计数过程;合并速率相加,分流速率按比例。

自测

1. 机器故障前时间服从指数分布,平均 100 小时。已知它已经运行了 40 小时,接下来 10 小时内发生故障的概率是多少?

λ=0.01\lambda=0.01/小时。由无记忆性,P(X≤50∣X>40)=P(X≤10)=1−e−0.1=0.0952P(X\le50\mid X>40)=P(X\le10)=1-e^{-0.1}=0.0952。与已经运行了多久无关。

2. 设备在“工作”状态时,可能因为加工完成(速率 3/小时)、发生故障(速率 0.5/小时)或被调去做别的任务(速率 0.5/小时)而离开。它在“工作”状态平均停留多久?离开时是因为故障的概率是多少?

停留时间是三个指数时间的最小值,速率 3+0.5+0.5=43+0.5+0.5=4/小时,平均 15 分钟。因故障离开的概率 0.5/4=0.1250.5/4=0.125。

3. 网站每分钟平均收到 2 个请求(泊松)。其中 25% 是图片请求。1 分钟内恰好收到 1 个图片请求的概率是多少?

图片请求是速率 2×0.25=0.52\times0.25=0.5/分钟的泊松流。P=e−0.5×0.51/1!=0.3033P=e^{-0.5}\times0.5^1/1!=0.3033。