《算法设计与分析动态规划实例讲解.pptx》由会员分享,可在线阅读,更多相关《算法设计与分析动态规划实例讲解.pptx(82页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、动态规划是用来解决多阶段决策过程最优化的一种数量方法。其特点在于,它可以把一个n维决策问题变换为几个一维最优化问题,从而一个一个地去解决。需指出:动态规划是求解某类问题的一种方法,是考察问题的一种途径,而不是一种算法。必须对具体问题进行具体分析,运用动态规划的原理和方法,建立相应的模型,然后再用动态规划方法去求解。第1页/共82页即在系统发展的不同时刻(或阶段)根据系统所处的状态,不断地做出决策;每个阶段都要进行决策,目的是使整个过程的决策 达到最优效果。动态决策问题的特点:系统所处的状态和时刻是进行决策的重要因素;找到不同时刻的最优决策以及整个过程的最优策略。多阶段决策问题:在多阶段决策过程
2、中,系统的动态过程可以按照时间进程分为状态相互联系而又相互区别的各个阶段;第2页/共82页多阶段决策问题的典型例子:1.生产决策问题:企业在生产过程中,由于需求是随时间变化的,因此企业为了获得全年的最佳生产效益,就要在整个生产过程中逐月或逐季度地根据库存和需求决定生产计划。2.机器负荷分配问题:某种机器可以在高低两种不同的负荷下进行生产。在高负荷下进行生产时,产品的年产量g和投入生产的机器数量u1的关系为g=g(u1)12n状态决策状态决策状态状态决策第3页/共82页 这时,机器的年完好率为a,即如果年初完好机器的数量为u,到年终完好的机器就为au,0a1。在低负荷下生产时,产品的年产量h和投
3、入生产的机器数量u2的关系为 h=h(u2)假定开始生产时完好的机器数量为s1。要求制定一个五年计划,在每年开始时,决定如何重新分配完好的机器在两种不同的负荷下生产的数量,使在五年内产品的总产量达到最高。相应的机器年完好率b,0 b1。第4页/共82页 3.航天飞机飞行控制问题:由于航天飞机的运动的环境是不断变化的,因此就要根据航天飞机飞行在不同环境中的情况,不断地决定航天飞机的飞行方向和速度(状态),使之能最省燃料和实现目的(如软着落问题)。4.不包含时间因素的线性规划、非线性规划等静态决策问题(本质上是一次决策问题)也可以适当地引入阶段的概念,作为多阶段的决策问题用动态规划方法来解决。第5
4、页/共82页 5.最短路问题:给定一个交通网络图如下,其中两点之间的数字表示距离(或花费),试求从A点到G点的最短距离(总费用最小)。123456AB1B2C1C2C3C4D1D2D3E1E2E3F1F2G531368763685338422213335256643第6页/共82页(一)、基本概念1、阶段:把一个问题的过程,恰当地分为若干个相互联系的阶段,以便于按一定的次序去求解。描述阶段的变量称为阶段变量(k)。k=1,2,3,,n阶段的划分,一般是根据时间和空间的自然特征来进行的,但要便于问题转化为多阶段决策。2、状态:表示每个阶段开始所处的自然状况或客观条件。通常一个阶段有若干个状态,描
5、述过程状态的变量称为状态变量sk(表示第k阶段的状态变量)。年、月、路段一个数、一组数、一个向量状态变量的取值有一定的允许集合或范围,此集合称为状态允许集合SK=s1,s2,sk,。一、动态规划的基本思想一、动态规划的基本思想第7页/共82页3、决策:表示当过程处于某一阶段的某个状态时,可以作出不同的决定,从而确定下一阶段的状态,这种决定称为决策。描述决策的变量,称为决策变量。常用uk(sk)表示第k阶段当状态为sk时的决策变量。决策变量是状态变量的函数。可用一个数、一组数或一向量(多维情形)来描述。在实际问题中决策变量的取值往往在某一范围之内,此范围称为允许决策集合。常用Dk(sk)表示第k
6、阶段从状态sk出发的允许决策集合,显然uk(sk)Dk(sk)。第8页/共82页4、多阶段决策过程 可以在各个阶段进行决策,去控制过程发展的多段过程;其发展是通过一系列的状态转移来实现的;系统在某一阶段的状态转移不但与系统的当前的状态和决策有关,而且还与系统过去的历史状态和决策有关。第9页/共82页图示如下:状态转移方程是确定过程由一个状态到另一个状态的演变过程。如果第k阶段状态变量sk的值、该阶段的决策变量一经确定,第k+1阶段状态变量sk+1的值也就确定。其状态转移方程如下(一般形式)12ks1u1s2u2s3skuksk+1 能用动态规划方法求解的多阶段决策过程是一类特殊的多阶段决策过程
7、,即具有无后效性的多阶段决策过程。第10页/共82页 如果状态变量不能满足无后效性的要求,应适当地改变状态的定义或规定方法。动态规划中能处理的状态转移方程的形式。状态具有无后效性的多阶段决策过程的状态转移方程如下无后效性(马尔可夫性)如果某阶段状态给定后,则在这个阶段以后过程的发展不受这个阶段以前各段状态的影响;过程的过去历史只能通过当前的状态去影响它未来的发展;构造动态规划模型时,要充分注意是否满足无后效性的要求;状态变量要满足无后效性的要求;第11页/共82页 5、策略:相互连接的决策序列称为一个策略。从第k阶段开始到第n阶段结束称为一个子策略。Pk,n ,全策略 P1,n .所有策略当中
8、有最优值的策略称为最优策略。6、状态转移方程:是确定过程由一个状态到另一个状态的演变过程,描述了状态转移规律。第12页/共82页 7、指标函数和最优值函数:用来衡量所实现过程优劣的一种数量指标,为指标函数。阶段指标函数:Vk(sk,uk)表示第 k 阶段位于sk 状态、决策为 uk 的指标值。策略指标函数:各决策序列指标值之和。(个别情况为乘积)指标函数的最优值,称为最优值函数。在不同的问题中,指标函数的含义是不同的,它可能是距离、利润、成本、产量或资源消耗等。动态规划模型的指标函数,应具有可分离性,并满足递推关系。第13页/共82页小结:方程:状态转移方程概念:阶段变量k状态变量sk决策变量
9、uk;指标:动态规划本质上是多阶段决策过程;效益指标函数形式:和、积无后效性可递推第14页/共82页解多阶段决策过程问题,求出 最优策略,即最优决策序列f1(s1)最优轨线,即执行最优策略时的状态序列 最优目标函数值从 k 到终点最优策略子策略的最优目标函数值第15页/共82页1、动态规划方法的关键在于正确地写出基本的递推关系式和恰当的边界条件(简称基本方程)。要做到这一点,就必须将问题的过程分成几个相互联系的阶段,恰当的选取状态变量和决策变量及定义最优值函数,从而把一个大问题转化成一组同类型的子问题,然后逐个求解。即从边界条件开始,逐段递推寻优,在每一个子问题的求解中,均利用了它前面的子问题
10、的最优化结果,依次进行,最后一个子问题所得的最优解,就是整个问题的最优解。(二)、动态规划的基本思想(二)、动态规划的基本思想第16页/共82页2、在多阶段决策过程中,动态规划方法是既把当前一段和未来一段分开,又把当前效益和未来效益结合起来考虑的一种最优化方法。因此,每段决策的选取是从全局来考虑的,与该段的最优选择答案一般是不同的.最优化原理:作为整个过程的最优策略具有这样的性质:无论过去的状态和决策如何,相对于前面的决策所形成的状态而言,余下的决策序列必然构成最优子策略。”也就是说,一个最优策略的子策略也是最优的。3、在求整个问题的最优策略时,由于初始状态是已知的,而每段的决策都是该段状态的
11、函数,故最优策略所经过的各段状态便可逐段变换得到,从而确定了最优路线。第17页/共82页(三)、建立动态规划模型的步骤 1、划分阶段k划分阶段是运用动态规划求解多阶段决策问题的第一步,在确定多阶段特性后,按时间或空间先后顺序,将过程划分为若干相互联系的阶段。对于静态问题要人为地赋予“时间”概念,以便划分阶段。2、正确选择状态变量sk选择变量既要能确切描述过程演变又要满足无后效性,而且各阶段状态变量的取值能够确定。一般地,状态变量的选择是从过程演变的特点中寻找。3、确定决策变量uk(sk)及允许决策集合Dk(sk)通常选择所求解问题的关键变量作为决策变量,同时要给出决策变量的取值范围,即确定允许
12、决策集合。第18页/共82页 4、确定状态转移方程根据k 阶段状态变量和决策变量,写出k+1阶段状态变量,状态转移方程应当具有递推关系。sk+1=Tk (sk,uk)Tk 函数关系 5、确定阶段指标函数和最优指标函数,建立动态规划基本方程 阶段指标函数是指第k 阶段的收益,最优指标函数是指从第k 阶段状态出发到第n 阶段末所获得收益的最优值,最后写出动态规划基本方程。fk(sk)=OptVk(sk,uk)+fk+1(sk+1)fn+1(sn+1)=0Opt最优化(max,min)第19页/共82页 以上五步是建立动态规划数学模型的一般步骤。由于动态规划模型与线性规划模型不同,动态规划模型没有统
13、一的模式,建模时必须根据具体问题具体分析,只有通过不断实践总结,才能较好掌握建模方法与技巧。f1(s1)是整个问题的最优策略,最优值。fk(sk)表示从第k阶段(状态sk)到终点的最优指标值。(距离、利润、成本等)第20页/共82页例一、从A地到D地要铺设一条煤气管道,其中需经过两级中间站,两点之间的连线上的数字表示距离,如图所示。问应该选择什么路线,使总距离最短?AB1B2C1C2C3D24333321114二、最短路径问题二、最短路径问题第21页/共82页解:整个计算过程分三个阶段,从最后一个阶段开始。第三阶段(C D):C 有三条路线到终点D。AB1B2C1C2C3D2433332111
14、4DC1C2C3显然有f3(C1)=1;f3(C2)=3;f3(C3)=4第22页/共82页 d(B1,C1)+f3(C1)3+1 f2(B1)=mind(B1,C2)+f3(C2)=min3+3d(B1,C3)+f3(C3)1+44=min6=45第二阶段(B C):B 到C 有六条路线。AB1B2C1C2C3D24333321114DC1C2C3B1B2(最短路线为B1C1D)第23页/共82页 d(B2,C1)+f3(C1)2+1 f2(B2)=mind(B2,C2)+f3(C2)=min3+3d(B2,C3)+f3(C3)1+43=min6=35AB1B2C1C2C3D24333321
15、114DC1C2C3B1B2(最短路线为B2C1D)第24页/共82页第一阶段(A B):A 到B 有二条路线。f1(A)1=d(A,B1)f2(B1)246f1(A)2=d(A,B2)f2(B2)437f1(A)=min=min6,7=6d(A,B1)f2(B1)d(A,B2)f2(B2)(最短路线为AB1C1D)AB1B2C1C2C3D24333321114DC1C2C3B1B2A第25页/共82页AB1B2C1C2C3D24333321114DC1C2C3B1B2A最短路线为AB1C1D路长为6第26页/共82页三、非线性规划问题三、非线性规划问题【例例7-47-4】用动态规划方法解下列
16、非线性规划问题用动态规划方法解下列非线性规划问题第27页/共82页解解:解解决决这这一一类类静静态态规规划划问问题题,需需要要人人为为地地赋赋予予时时间间概概念念,从从而而将将该该问问题题转转化化为为多多阶阶段决策过程。段决策过程。按问题的变量个数划分阶段,把它看作一个三阶段决策问题,按问题的变量个数划分阶段,把它看作一个三阶段决策问题,k=1k=1,2 2,3 3设状态变量为设状态变量为s s1 1,s s2 2,s s3 3,s s4 4并记并记s s1 1cc取问题中的变量取问题中的变量x x1 1,x x2 2,x x3 3为决策变量为决策变量第28页/共82页状态转移方程为:状态转移
17、方程为:s3=x3s3+x2=s2s2+x1=s1c允许决策集合为:允许决策集合为:x3=s30 x2s20 x1s1各阶段指标函数为:各阶段指标函数为:v1(x1)=x1v2(x2)=x22v3(x3)=x3,各指标函数以乘积方式结合,最优指标函数各指标函数以乘积方式结合,最优指标函数fk(sk)表示从第表示从第k k阶段初始状态阶段初始状态sk出发到第出发到第3 3阶段阶段所得到的最大值,则动态规划基本方程为:所得到的最大值,则动态规划基本方程为:第29页/共82页用逆序解法由后向前依次求解:用逆序解法由后向前依次求解:k=3k=3时时,x3*=s3 k=2 k=2时时,第30页/共82页
18、令h2(s2,x2)=x22(s2x2)用用经典解析法经典解析法求极值点:求极值点:解得:解得:x2=0(舍)(舍)所以所以是极大值点。是极大值点。第31页/共82页k=1时,令解得:x1=s1(舍)所以所以 是极大值点。是极大值点。第32页/共82页由于由于s1未知,所以对未知,所以对s1再求极值,再求极值,显然显然s s1 1=c c时,时,f f1 1(s s1 1)取得最大值取得最大值 反向追踪得各阶段最优决策及最优值:反向追踪得各阶段最优决策及最优值:s1=c所以最优解为:所以最优解为:第33页/共82页一一般般地地,如如果果阶阶段段指指标标函函数数vk(sk,uk)是是线线性性函函
19、数数或或凸凸函函数数时时,最最优优指指标标函函数数fk(sk)的的表表达达式式比比较较容容易易得得到到,但但是是当当vk(sk,uk)不不具具备备上上述述特特性性时时,最最优优指指标标函函数数fk(sk)的的表表达达式式不不易易得得到到,就就需需要要采采用用数数值值法法,即即对对连连续续变变量量进进行行离离散散化化处处理理,再再分散求解。分散求解。例如静态规划模型其动态规划基本方程为:其动态规划基本方程为:状态转移方程为状态转移方程为sk+1=skxks1=a第34页/共82页状状态态变变量量sk及及决决策策变变量量xk都都是是连连续续变变量量,对对其其进进行行离离散化处理,具体做法是:散化处
20、理,具体做法是:1.1.对区间对区间00,aa进行分割,分割数进行分割,分割数m=m=,其中,其中是是分割后的小区间的长度,其大小可以根据所求解问分割后的小区间的长度,其大小可以根据所求解问题要求的精度及计算机运算能力而定,分割点为题要求的精度及计算机运算能力而定,分割点为0 0,22,m=am=a。2.2.规定状态变量规定状态变量sk及决策变量及决策变量xk仅在离散点仅在离散点0 0,22,mm处取值,最优指标函数处取值,最优指标函数fk(sk)也定也定义在这些离散点上。动态规划基本方程可以写为:义在这些离散点上。动态规划基本方程可以写为:其中其中s sk k=q q,x xk k=p p。
21、3.3.由后向前逐段递推,直至求出整个过程最优解。由后向前逐段递推,直至求出整个过程最优解。第35页/共82页【例7-5】解解 按按变变量量个个数数将将原原问问题题分分为为三三个个阶阶段段,阶阶段段变变量量k=1,2,3;选择选择xk为决策变量;为决策变量;状态变量状态变量sk表示第表示第k阶段至第阶段至第3阶段决策变量之和;阶段决策变量之和;取小区间长度取小区间长度=1=1,小区间数目,小区间数目m=6/1=6m=6/1=6,状态,状态变量变量sk的取值点为:的取值点为:状态转移方程:状态转移方程:sk+1=skxk;允许决策集合:允许决策集合:Dk(sk)=xk|0 xkskk=1,2,3
22、xk,sk均在分割点上取值;均在分割点上取值;第36页/共82页阶段指标函数分别为:阶段指标函数分别为:g1(x1)=x12g2(x2)=x2g3(x3)=x33,最优指标函数最优指标函数fk(sk)表示从第表示从第k阶段状态阶段状态sk出发到第出发到第3阶段所得到的最大值,动态规阶段所得到的最大值,动态规划的基本方程为:划的基本方程为:k=3时,时,s3及及x3取值点较多,计算结果以表格形式给出,见表取值点较多,计算结果以表格形式给出,见表7-17-1所示。所示。第37页/共82页表71取值第38页/共82页k=2k=2时,时,计算结果见表计算结果见表7-27-2第39页/共82页k=1k=
23、1时,时,其中其中s1=6,计算结果见表计算结果见表7-37-3所示。所示。由由表表7-37-3知知,x1*=2,s1=6,则则s2=s1x1*=62=4,查查 表表 7-27-2得得:x2*=1,则则s3=s2x2*=41=3,查查表表7-17-1得得:x3*=3,所所以以最最优优解解为为:x1*=2,x2*=1,x3*=3,f1(s1)=108。本例也可用经典解析法求得各段的极值,读者本例也可用经典解析法求得各段的极值,读者可自行求解,二者结论完全相同。需要指出的是当连可自行求解,二者结论完全相同。需要指出的是当连续变量离散化处理以后,由于状态变量和决策变量只续变量离散化处理以后,由于状态
24、变量和决策变量只在给定的离散点上取值,故有可能漏掉最优解,因此在给定的离散点上取值,故有可能漏掉最优解,因此需要慎重选择参数需要慎重选择参数m m与与。第40页/共82页资源分配问题就是将一定数量的一种或若干种资源(原材料、资金、设备等)合理分配给若干使用者,使得资源分配后总结果最优。一种资源的分配问题称为一维资源分配问题,两种资源的分配问题称为二维资源分配问题。四、四、资源分配问题资源分配问题第41页/共82页假设有一种资源,数量为a a,将其分配给n n个使用者,分配给第i个使用者数量xi时,相应的收益为gi(xi),问如何分配使得总收入最大?这就是一维资源分配问题,该问题的数学模型为:这
25、是一个静态规划问题,应用动态规划方法求解时人为赋予时间概念,将其看作是一个多阶段决策问题。第42页/共82页按变量个数划分阶段,k=1,2,n;设决策变量uk=xk,表示分配给第k个使用者的资源数量;设状态变量为sk,表示分配给第k个至第n个使用者的总资源数量;状态转移方程:sk+1=skxk,其中s1=a;允许决策集合:Dk(sk)=xk|0 xksk阶段指标函数:vk(sk,uk)=gk(xk)表示分配给第k个使用者数量xk时的收益;最优指标函数fk(sk)表示以数量sk的资源分配给第k个至第n个使用者所得到的最大收益,则动态规划基本方程为:由后向前逐段递推,f1(a)即为所求问题的最大收
26、益。第43页/共82页【例例7-67-6】某某公公司司打打算算在在3 3个个不不同同的的地地区区设设置置4 4个个销销售售点点,根根据据市市场场部部门门估估计计,在在不不同同地地区区设设置置不不同同数数量量的的销销售售点点每每月月可可得得到到的的利利润润如如表表7-47-4所所示示。试试问问在在各各地地区区如如何何设设置置销销售点可使每月总利润最大。售点可使每月总利润最大。表表7-47-4第44页/共82页解 如前所述,建立动态规划数学模型:将问题分为3 3个阶段,k=1,2,3;决策变量xk表示分配给第k个地区的销售点数;状态变量为sk表示分配给第k个至第3个地区的销售点总数;状态转移方程:
27、sk+1=skxk,其中s1=4;允许决策集合:Dk(sk)=xk|0 xksk阶段指标函数:gk(xk)表示xk个销售点分配给第k个地区所获得的利润;最优指标函数fk(sk)表示将数量为sk的销售点分配给第k个至第3个地区所得到的最大利润,动态规划基本方程为:第45页/共82页k=3时,数值计算如下表7-57-5表7-5 7-5 第46页/共82页k=2k=2时,计算结果见下表7-67-6表7-67-6第47页/共82页k=1k=1时,k=1k=1时,只有s s1 1=4=4的情况。计算结果如表7-77-7所示。所以最优解为:x x1 1*=2=2,x x2 2*=1=1,x x3 3*=1
28、=1,f f1 1(4)=47(4)=47,即在第1 1个地区设置2 2个销售点,第2 2个地区设置1 1个销售点,第3 3个地区设置1 1个销售点,每月可获利润4747。表7-77-7第48页/共82页【例例7-7】机器负荷问题机器负荷问题某工厂有100100台机器,拟分四期使用,每一期都可在高、低两种不同负荷下进行生产。若把x x台机器投入高负荷下进行生产,则在本期结束时将有1/3x1/3x台机器损坏报废;余下的机器全部投入低负荷下进行生产,则在期末有1/101/10的机器报废。如果高负荷下生产时每台机器可获利润为1010,低负荷下生产时每台机器可获利润为7 7,问怎样分配机器使四期的总利
29、润最大?解 将问题按周期分为4 4个阶段,k=1,2,3,4;状 态 变 量sk表 示 第k阶 段 初 完 好 的 机 器 数,s1=100,0sk100;决策变量xk表示第k阶段投入高负荷下生产的机器数,则skxk表示第k阶段投入低负荷下生产的机器数;允许决策集合:Dk(sk)=xk|0 xksk状态转移方程:sk+1=Tk(sk,xk),即第k+1阶段初拥有的完好机器数sk+1为:第49页/共82页阶段指标函数:vk(sk,xk)=10 xk+7(skxk)表示第k阶段所获得的利润;最优指标函数f fk k(s sk k)表示从第k k阶段初完好机器数为s sk k至第四阶段的最大利润,动
30、态规划基本方程为:k=4k=4时,x x4 4*=s s4 4 k=3k=3时,x x3 3*=s s3 3第50页/共82页k=2时,x2*=0k=1时,x1*=0第51页/共82页因为s1=100,所以f1(100)=2680,逆向追踪得:s1=100,x1*=0 x2*=0 x3*=s3=81x4*=s4=54即,第1,2期把全部完好机器投入低负荷下生产,第3,4期把全部完好机器投入高负荷下生产所得利润最大。第52页/共82页五、生产计划问题五、生产计划问题 在企业生产经营活动中,经常会遇到如何合理安排生产、库存及销售计划,使总效益最高的问题,这一类问题统称为生产计划问题。第53页/共8
31、2页【例例例例7-87-8】(生产(生产库存问题)库存问题)某工厂要对一种产品制定今后四个时期的生产计划,据估计在今后四个时期内,市场对该产品的需求量分别为2,3,2,4单位,假设每批产品固定成本为3千元,若不生产为0,每单位产品成本为1千元,每个时期最大生产能力不超过6个单位,每期期末未出售产品,每单位需付存贮费0.5千元,假定第1期初和第4期末库存量均为0,问该厂如何安排生产与库存,可在满足市场需求的前提下总成本最小。解解 以以每每个个时时期期作作为为一一个个阶阶段段,该该问问题题分分为为4 4个个阶阶段段,k k=1=1,2 2,3 3,4 4;决策变量决策变量x xk k表示第表示第k
32、 k阶段生产的产品数;阶段生产的产品数;状态变量状态变量s sk k表示第表示第k k阶段初的库存量;阶段初的库存量;第54页/共82页以dk表示第k阶段的需求,则状态转移方程:sk+1=sk+xkdkk=4,3,2,1由于期初及期末库存为0,所以s1=0,s5=0;允许决策集合Dk(sk)的确定:当skdk时,xk可以为0,当skdk时,至少应生产dksk,故xk的下限为max(0,dksk);每期最大生产能力为6,xk最大不超过6,由于期末库存为0,xk还应小于本期至4期需求之和减去本期的库存量,所以xk的上限为min(,6),故有:Dk(sk)=xk|max(0,dksk)xkmin(,
33、6)第55页/共82页阶段指标函数rk(sk,xk)表示第k期的生产费用与存贮费用之和:最优指标函数fk(sk)表示第k期库存为sk到第4期末的生产与存贮最低费用,动态规划基本方程为:第56页/共82页先求出各状态允许状态集合及允许决策集合,如表7-8所示。表表7-8第57页/共82页k=4时,计算结果见表7-9所示。表表7-9第58页/共82页k=3时,计算结果如下表:第59页/共82页k=2时,计算结果如下表第60页/共82页k=1时,计算结果见表7-12所示 逆向追踪可得:x1*=5,s2=3,x2*=0,s3=0,x3*=6,s4=4,x4*=0,即第1时期生产5个单位,第3时期生产6
34、个单位,第2,4时期不生产,可使总费用最小,最小费用为20.5千元。第61页/共82页【例7-9】(库存销售问题)设某公司计划在1至4月份从事某种商品经营。已知仓库最多可存储600件这种商品,已知1月初存货200件,根据预测知1至4月份各月的单位购货成本及销售价格如表7-13所示,每月只能销售本月初的库存,当月进货供以后各月销售,问如何安排进货量和销售量,使该公司四个月获得利润最大(假设四月底库存为零)。表7-13第62页/共82页解按月份划分阶段,k=1,2,3,4;状态变量sk表示第k月初的库存量,s1=200,s5=0;决策变量xk表示第k月售出的货物数量,yk表示第k月购进的货物数量;
35、状态转移方程:sk+1=sk+ykxk;允许决策集合:0 xksk,0yk600(skxk);阶段指标函数为:pkxkckyk表示k月份的利润,其中pk为第k月份的单位销售价格,ck为第k月份的单位购货成本;最优指标函数fk(sk)表示第k月初库存为sk时从第k月至第4月末的最大利润,则动态规划基本方程为:第63页/共82页k=4时,x4*=s4y4*=0k=3时,为求出使44s35x3+4y3最大的x3,y3,须求解线性规划问题:第64页/共82页只有两个变量x3,y3,可用图解法也可用单纯形法求解,图解法求解示意图如图7-5所示:在A点处取得最优解,x3*=0,y3*=600s3,f3(s
36、3)=40s3+2400As3600y3x30600s3图7-5第65页/共82页k=2时,类似地求得:x2*=s2,y2*=600,f2(s2)=42s2+3600k=1时,类似地求得:x1*=s1,y1*=600,f1(s1)=45s1+4800=13800第66页/共82页逆向追踪得各月最优购货量及销售量:x1*=s1=200y1*=600;x2*=s2=s1+y1*x1*=600y2*=600;x3*=0y3*=600s3=600(s2+y2*x2*)=0 x4*=s4=(s3+y3*x3*)=600y4*=0即1月份销售200件,进货600件,2月份销售600件,进货600件,3月份
37、销售量及进货量均为0,4月份销售600件,不进货,可获得最大总利润13800。第67页/共82页六、背包问题六、背包问题有人携带背包上山,其可携带物品的重量限度为a公斤,现有n种物品可供选择,设第i种物品的单件重量为ai公斤,其在上山过程中的价值是携带数量xi的函数ci(xi),问应如何安排携带各种物品的数量,使总价值最大。这就是背包问题,类似的货物装载问题,下料问题都等同于背包问题。背包问题的数学模型为:第68页/共82页下面用动态规划方法求解:按照装入物品的种类划分阶段,k=1,2,n;状态变量sk表示装入第k种至第n种物品的总重量;决策变量xk表示装入第k种物品的件数;状态转移方程为:s
38、k+1=skakxk允许决策集合为:其中表示不超过的最大整数;阶段指标函数ck(xk)表示第k阶段装入第k种商品xk件时的价值;最优指标函数fk(sk)表示第k阶段装入物品总重量为sk时的最大价值,动态规划基本方程为:第69页/共82页【例7-10】某工厂生产三种产品,各产品重量与利润关系如表7-14所示,现将此三种产品运往市场销售,运输能力总重量不超过6吨,问如何安排运输使总利润最大?表7-14种类123单位重量(吨)234单位利润(元)80130180解设xi为装载第i种货物的件数,i=1,2,3,该问题数学模型为:第70页/共82页按前述方法建立动态规划模型;k=3时,计算结果如表7-1
39、5所示。第71页/共82页k=2时,计算结果如表7-16所示。表7-16第72页/共82页k=1时,计算结果如表7-17所示。表7-17反向追踪得最优方案:x1*=0,x2*=2,x3*=0;最优方案:x1*=1,x2*=0,x3*=1;最大总利润为260元。第73页/共82页七、复合系统工作可靠性问题七、复合系统工作可靠性问题某个机器工作系统由n个部件串联而成,其中只要有一个部件失效,则整个系统不能正常工作,因此为了提高系统工作的可靠性,在设计时,每个主要部件上都装有备用元件,一旦某个主要部件失效,备用元件会自动投入系统工作,显然备用元件越多,系统工作可靠性越大,但是备用元件越多,系统的成本
40、、重量、体积相应增大,工作精度降低,因此在上述限制条件下,应选择合理的备用元件数,使整个系统的工作可靠性最大。第74页/共82页设第i(i=1,2,n)个部件上装有ui个备用元件,正常工作的概率为pi(ui),则整个系统正常工作的可靠性为,装第i个部件的费用为ci,重量为wi,要求总费用不超过c,总重量不超过w,则静态规划数学模型为:第75页/共82页按部件个数划分阶段,k=1,2,n;决策变量uk表示部件k上的备用元件数;状态变量xk表示从第k个到第n个部件的总费用,yk表示从第k个到第n个部件的总重量;状态转移方程为:xk+1=xkckukyk+1=ykwkuk允许决策集合为:阶段指标函数
41、为pk(uk),表示第k个部件的正常工作概率;最优指标函数fk(xk,yk)表示由状态xk,yk出发,从部件k到部件n的系统工作最大可靠性,则动态规划基本方程为:f1(c,w)即为整个系统工作的最大可靠性。第76页/共82页【例7-11】某厂设计的一种电子设备由三种元件A、B、C串联而成,已知三种元件的价格及可靠性如表7-18所示,要求设计中使用元件的总费用不超过10万元,问如何设计使设备的可靠性达到最大(不考虑重量限制)。表7-18第77页/共82页解如前所述建立动态规划数学模型;按元件种类划分为3个阶段,k=1,2,3;决策变量xk表示第k个部件配备的元件数;状态变量sk表示从第k阶段到第
42、3阶段配备元件的总费用;状态转移方程为:sk+1=skckxk其中ck表示第k种部件的元件单价;允许决策集合为:以pk表示第k个部件中的1个元件的正常工作概率,假定部件k的xk个元件是并联的,则为xk个元件均不正常工作的概率,fk(sk)表示由状态sk开始从第k个到第3个部件的设备最大可靠性。第78页/共82页k=3时,由于A、B至少要购置1件,用于购置C的最高金额为s3=1023=5万元,计算结果如表7-19所示。表7-19第79页/共82页k=2时,计算结果如表7-20所示。表7-20第80页/共82页k=1时,计算结果如表7-21所示。表7-21逆向追踪得:x1*=2,s2=6,x2*=1,s3=3,x3*=3,即A元件用2个,B元件用1个,C元件用3个,最高可靠性为0.682。第81页/共82页谢谢您的观看!第82页/共82页