網络与规划模型实验.docVIP

网络与规划模型实验 [学习目标] 能进行机器负荷分配问题中的动态规划运算； 会进行食谱问题中的线性规划运算； 会进行污水处理控制问题中的线性规划运算。 机器负荷分配问题 问题 设某工厂有1000台机器，生产两种产品A、B，若投入y台机器生产A产品，则纯收入为5y，若投入y台机器生产b产品，则纯收入为4y，又知：生产A产品机器的年折损率为20%,生产B产品机器的年折损率为10%，问在5年内如何安排各年度的生产计划，才能使总收入最高？ 问题分析与建立模型 问题可分为5个阶段（k=1，2，3，4，5）。 第k个阶段——第k年初到k+1年初 令 xk——第k年初完好机器数 uk——第k年安排生产A种产品的机器数 则显然有： xk - uk——第k年安排生产B种产品的机器数，且0≤uk≤xk 在以上假设的基础上，则知： 第k + 1年初完好的机器数=（1 - 生产A产品的机器折旧率（20%））×第k年安排生产A产品的机器数+（1 - 生产B产品的机器折旧率（10%））×第k年安排生产B产品的机器数，即 xk +1 =（1 - 0.2）uk+（1 - 0.1）（xk - uk） =0.8 uk + 0.9（xk - uk） =0.9xk - 0.1uk 又令：L（xk，uk）——第k年的纯收入； vk（xk）——第k年初往后各年的最大利润之和。 计算过程 算法 第1步 定义 X=Table{x（i），{i，5}} U=Table{u（i），{i，5}} V=Table{v（i），{i，5}} 第2步 v6（x6）=0 第3步 f（x，y，z）=5y + 4（x - y）+ z， g（x，y）=0.9x - 0.1y 第4步 k = 5 第5步 Price=f（xk，uk，vk +1） 其中vk +1= g（xk，uk） 若：Price关于uk 的导数大于0，则：uk = xk 否则uk =0 vk（xk）= Price 第6步 k = k - 1 第7步 判断k 1？若是则转5，否则，转8 第8步 x1 = 1000 第9步 i = 2 第10步 xi = 0.9xi -1 - 0.1ui -1 第11步 i = i + 1 第12步 若i = 6，则转13，否则，转10 第13步 输出X，U，X - U，最大利润。 源程序 Clear[X，U，V，x，y，z]； X=Table[x[i]，{i，5}]；U=Table[u[i]，{i，5}]；V=Table[v[i]，{i，5}]； v[6]=0； f[x_，y_，z_]：=5y+4（x-y）+z；g[x_，y_]：=0.9x - 0.1y； k=5； While[k=1，xh=f[x[k]，u[k]，v[k+1]]/.x[k+1]→g[x[k]，u[k]]； If[D[xh，u[k]]0， u[k]= x[k]；v[k]=xh/.u[k]→x[k]，u[k]=0；v[k]= xh/.u[k]→0；]；k――；]； x[1]=1000；i=2； While[i=5， x[i]=N[0.9x[i-1]-0.1u[i-1]]；i++；]； Print[″X=″，X] Print[″U=″，U] Print[″X-U=″，X-U] Print[″MAXPRICE=″，V[[1]]] 运行结果为： X={1000，900，810，648，518.4} U={0，0，810，648，518.4} X-U={1000，900，0，0，0} MAXPRICE=17482 结果分析 这是一个很不错的数学实验素材，它用一般人最关心的问题——利润最大，来引发读者探寻如何才能达到最优，在此基础上，使读者充分了解数学在实际工作中的重要作用，激发起学习数学的兴趣。 由结果看出，这种资源分配方式出乎许多人意料之外，它不是按一般人的想法，先投入大量机器生产A产品而盲目追逐“高”利润，而是在前几年先保护机器，然后再投入所有的力量去获得最大利润，这是一个优秀的企业管理者应该具有的素质。