数学 运筹学

运筹学:第 4 次作业

三阶段试制问题的随机动态规划解答,给出状态、转移概率、Bellman 递推与最优批量;原始题面尚待补录。

本页目录6 节

随机动态规划解答 #

内容状态:解答存档。 当前资料仅包含解答,原始题面尚待补录。以下保留能够确认的参数与计算过程。

已知参数 #

设单件样品合格的概率为

p=0.4,1p=0.6.p=0.4,\qquad 1-p=0.6.

每次试制一批 nn 件样品时的费用为

C(n)=200+100n,n=1,2,,5.C(n)=200+100n,\qquad n=1,2,\dots,5.

其中 200200 元为每批装配费,100n100n 元为本批 nn 件样品的制造费。三个月内未能交出合格样品时,需支付违约金 20002000 元。

1. 状态与决策 #

将三次试制视作三个阶段 k=1,2,3k=1,2,3。在第 kk 次试制前,状态变量记为

sk={0,到目前为止尚未制出合格样品,1,到目前为止已经制出至少一件合格样品.s_k= \begin{cases} 0, & \text{到目前为止尚未制出合格样品},\\[2mm] 1, & \text{到目前为止已经制出至少一件合格样品}. \end{cases}

kk 阶段的决策变量为

xk=n{1,2,3,4,5},x_k=n\in\{1,2,3,4,5\},

表示在本批中生产 nn 件样品。

若在状态 sk=0s_k=0 时选择 xk=nx_k=n,本批中至少出现一件合格样品的概率为

1(1p)n=10.6n,1-(1-p)^n=1-0.6^n,

全部不合格的概率为 0.6n0.6^n。因此

P(sk+1=1sk=0,xk=n)=10.6n,P(sk+1=0sk=0,xk=n)=0.6n.\begin{aligned} P(s_{k+1}=1\mid s_k=0,x_k=n)&=1-0.6^n,\\ P(s_{k+1}=0\mid s_k=0,x_k=n)&=0.6^n. \end{aligned}

若某阶段开始前已经有合格品,即 sk=1s_k=1,则以后保持该状态且不再生产:

P(sk+1=1sk=1)=1,P(sk+1=0sk=1)=0.P(s_{k+1}=1\mid s_k=1)=1,\qquad P(s_{k+1}=0\mid s_k=1)=0.

2. Bellman 递推 #

把第三次试制结束后的时刻记为阶段 k=4k=4。令 fk(s)f_k(s) 表示第 kk 阶段开始时系统处于状态 ss,从当前阶段到合同结束的最小期望总费用。

终端条件为

f4(1)=0,f4(0)=2000.f_4(1)=0,\qquad f_4(0)=2000.

k=1,2,3k=1,2,3,已有合格品时不再发生试制费用,因此

fk(1)=0.f_k(1)=0.

sk=0s_k=0 时选择 xk=nx_k=n,Bellman 方程为

fk(0)=min1n5{C(n)+(10.6n)fk+1(1)+0.6nfk+1(0)}=min1n5{200+100n+0.6nfk+1(0)}.\begin{aligned} f_k(0) &=\min_{1\le n\le 5}\Big\{ C(n)+(1-0.6^n)f_{k+1}(1)+0.6^n f_{k+1}(0) \Big\}\\ &=\min_{1\le n\le 5}\big\{200+100n+0.6^n f_{k+1}(0)\big\}. \end{aligned}

3. 向后递推 #

阶段 3

f4(0)=2000f_4(0)=2000,得到

f3(0)=min1n5{200+100n+0.6n2000}=min{1500, 1120, 932, 859.2, 855.52}.\begin{aligned} f_3(0) &=\min_{1\le n\le 5}\{200+100n+0.6^n\cdot2000\}\\ &=\min\{1500,\ 1120,\ 932,\ 859.2,\ 855.52\}. \end{aligned}

因此 f3(0)=855.52f_3(0)=855.52,最优决策为 n3=5n_3^*=5

阶段 2

f2(0)=min1n5{200+100n+0.6nf3(0)}=min{813.312, 707.9872, 684.79232, 710.875392, 766.5252352}.\begin{aligned} f_2(0) &=\min_{1\le n\le 5}\{200+100n+0.6^n f_3(0)\}\\ &=\min\{813.312,\ 707.9872,\ 684.79232,\ 710.875392,\ 766.5252352\}. \end{aligned}

因此 f2(0)=684.79232f_2(0)=684.79232,最优决策为 n2=3n_2^*=3

阶段 1

f1(0)=min1n5{200+100n+0.6nf2(0)}=min{710.875392, 646.5252352, 647.91514112, 688.749084672, 753.2494508032}.\begin{aligned} f_1(0) &=\min_{1\le n\le 5}\{200+100n+0.6^n f_2(0)\}\\ &=\min\{710.875392,\ 646.5252352,\ 647.91514112,\ 688.749084672,\ 753.2494508032\}. \end{aligned}

因此 f1(0)=646.5252352f_1(0)=646.5252352,最优决策为 n1=2n_1^*=2

结论 #

最优策略为:

{第 1 批生产 2 件样品;若尚无合格品,第 2 批生产 3 件样品;若仍无合格品,第 3 批生产 5 件样品.\begin{cases} \text{第 1 批生产 }2\text{ 件样品};\\ \text{若尚无合格品,第 2 批生产 }3\text{ 件样品};\\ \text{若仍无合格品,第 3 批生产 }5\text{ 件样品}. \end{cases}

从初始状态 s1=0s_1=0 出发,最小期望总费用为

f1(0)646.53 元.f_1(0)\approx 646.53\ \text{元}.

讨论

评论

正在加载评论…

输入关键词开始搜索。