北航计算机研究生课程 算法设计与分析Assignment1.docx

1、北航计算机研究生课程 算法设计与分析 Assignment1一、解：设第k月的需求量为Nk(k=1,2,3,4)状态变量Xk：第k月初的库存量，X1=X5=0，0XkNk+N4决策变量Uk：第k月的生产量，max0，Nk-XkUkmin6，Nk+N4 - Xk状态转移方程：Xk+1 = Uk + Xk Nk第k月的成本Vk = 0.5*(Xk - Nk) Uk=0 3 + Uk + 0.5*(Uk + Xk - Nk) Uk0设Fk(Xk)是由第k月初的库存量Xk开始到第4月份结束这段时间的最优成本则Fk(Xk) = minVk + Fk+1(X k+1) 1k4= min 3 + Uk +

2、0.5*(Uk + Xk - Nk) + Fk+1(Uk + Xk - Nk) Uk0 min 0.5*(Xk - Nk) + Fk+1(Xk - Nk) Uk=0F5(X5)=0四个月内的最优成本为F1(X1)=F1(0)详细计算步骤如下：（1）k=4时0X44，max0，4 - X4U4min6，4-X4X4U4X5V4F5(X5)V4 + F5(X5)040707=F4(0)130606=F4(1)220505=F4(2)310404=F4(3)400000=F4(4)即对于状态X4的每个取值，都有唯一确定的决策变量U4使得F4(X4)最优（2）k=3时0X36，max0，2 - X3U

3、3min6，6-X3X3U3X4V3F4(X4)V3 + F4(X4)0205712316.5612.5428513539.5413.56411011=F3(0)1104711215.5611.5327512438.5412.55410010=F3(1)200077=F3(2)114.5610.5226511337.5411.5449093010.566.5=F3(3)125510236.5410.534808402156=F3(4)135.549.5247075031.545.5=F3(5)14606604302=F3(6)（3）k=2时0X29，max0，3 - X2U2min6，9-X2

4、X2U2X3V2F3(X3)V2 + F3(X3)03061117417.51017.5529716=F2(0)6310.56.51712051116316.51016.5428715=F2(1)539.56.516641161721041115215.51015.5327714=F2(2)438.56.51554106166511.55.51730001111=F2(3)114.51014.5226713337.56.5144496155510.55.51666122144010.51010.5=F2(4)125713236.56.513348614459.55.51556112135021

5、78=F2(5)135.56.512247613358.55.51446102126031.56.58=F2(6)146612257.55.513369211704268=F2(7)156.55.5122682108052.55.58=F2(8)16729906325=F2(9)（4）k=1时X1=0，max0，2U1min6，11X1U1X2V1F2(X2)V1 + F2(X2)02051621316.51521.54281422539.51120.5=F1(0)641110.521.5由以上计算可得，4个月的总最优成本为F1(0) = 20.5(千元)从k=1回溯，可得最优结果中各阶段的状

6、态变量Xk和决策变量Uk如下表：月份k产量Uk月初库存量Xk需求量Nk每月成本Vk15029.52033036021140440二、解：1、变量设定阶段k：已遍历过k个结点，k=1,26,7。 K=1表示刚从V1出发，k=7表示已回到起点V1状态变量Xk=(i，Sk)：已遍历k个结点，当前位于i结点，还未遍历的结点集合为Sk。则X1=(1，2,3,4,5,6)，X6=(i，)，X7=(1，)决策变量Uk=(i，j)：已遍历k个结点，当前位于i结点，下一个结点选择j。状态转移方程：Xk+1 = T(Xk，Uk) = (j，Sk-j)第k阶段的指标函数Vk = Di,j。最优指标函数Fk(Xk)

7、= Fk(i，Sk)：已遍历k个结点，当前从i结点出发，访问Sk中的结点一次且仅一次，最后返回起点V1的最短距离。则Fk(i，Sk) = min Di,j + Fk+1(j，Sk-j) 1k6 F7(X7) = F7(1，) = 02、分析：（1）k=6时，F6(i，) = minDi,1 + F7(X7) = Di,1 i=2,3,4,5,6X6=(i，)U6=(i，j)X7=(1，)V6=Di,jF7(1，)V6 + F7(X7)(2，)(2，1)(1，)12012=F6(2,)(3，)(3，1)(1，)23023=F6(3,)(4，)(4，1)(1，)34034=F6(4,)(5，)(5

8、，1)(1，)45045=F6(5,)(6，)(6，1)(1，)56056=F6(6,)即k=6时，对于每一种状态X6，都有唯一的决策U6。（2）k=5时，F5(i，S5) = minDi,j + F6(j，) i=2,3,4,5,6X5=(i,S5)U5=(i,j)X6=(j, )V5=Di,jF6(j,)V5 + F6(X6)(2,6(2,6)(6，)215677=F5(2,6)(2,5(2,5)(5，)254570=F5(2,5)(2,4(2,4)(4，)303464=F5(2,4)(2,3(2,3)(3，)182341=F5(2,3)(3,6)(3,6)(6，)155671=F5(3,

9、6)(3,5)(3,5)(5，)104555=F5(3,5)(3,4)(3,4)(4，)53439=F5(3,4)(3,2)(3,2)(2，)191231=F5(3,2)(4,6)(4,6)(6，)165672=F5(4,6)(4,5)(4,5)(5，)84553=F5(4,5)(4,3)(4,3)(3，)42327=F5(4,3)(4,2)(4,2)(2，)321244=F5(4,2)(5,6)(5,6)(6，)185674=F5(5,6)(5,4)(5,4)(4，)103444=F5(5,4)(5,3)(5,3)(3，)112334=F5(5,3)(5,2)(5,2)(2，)271239=

10、F5(5,2)(6,5)(6,5)(5，)124557=F5(6,5)(6,4)(6,4)(4，)203454=F5(6,4)(6,3)(6,3)(3，)162339=F5(6,3)(6,2)(6,2)(2，)221234=F5(6,2)即k=时，对于每一种状态X5，都有唯一决策U5。（3）k=4时，F4(i,S4) = min(Di,j + F5(j,S5) ) i=2,3,4,5,6X4=(i,S4)U4=(i,j)X5=(j,S5)V4=Di,jF5(j,S5)V4 + F5(j,S5)(2,3,4)(2,3)(3,4)183957=F4(2,3,4)(2,4)(4,3)302757=F

11、4(2,3,4)(2,4,5)(2,4)(4,5)305383(2,5)(5,4)254469=F4(2,4,5)(2,5,6)(2,5)(5,6)257499(2,6)(6,5)215778=F4(2,5,6)(2,3,5)(2,3)(3,5)185573(2,5)(5,3)253459=F4(2,3,5)(2,3,6)(2,3)(3,6)187189(2,6)(6,3)213960=F4(2,3,6)(2,4,6)(2,4)(4,6)3072102(2,6)(6,4)215475=F4(2,4,6)(3,2,4)(3,2)(2,4)196483(3,4)(4,2)54449=F4(3,2,

12、4)(3,2,5)(3,2)(2,5)197089(3,5)(5,2)103949=F4(3,2,5)(3,2,6)(3,2)(2,6)197796(3,6)(6,2)153449=F4(3,2,6)(3,4,5)(3,4)(4,5)55358(3,5)(5,4)104454=F4(3,4,5)(3,4,6)(3,4)(4,6)57277(3,6)(6,4)155469=F4(3,4,6)(3,5,6)(3,5)(5,6)107484(3,6)(6,5)155772=F4(3,5,6)(4,2,3)(4,2)(2,3)324173(4,3)(3,2)43134=F4(4,2,3)(4,2,5)

13、(4,2)(2,5)3270102(4,5)(5,2)83947=F4(4,2,5)(4,2,6)(4,2)(2,6)3277109(4,6)(6,2)163450=F4(4,2,6)(4,3,5)(4,3)(3,5)45559(4,5)(5,3)83442=F4(4,3,5)(4,3,6)(4,3)(3,6)47175(4,6)(6,3)163955=F4(4,3,6)(4,5,6)(4,5)(5,6)87482(4,6)(6,5)165773=F4(4,5,6)(5,2,3)(5,2)(2,3)274168(5,3)(3,2)113142=F4(5,2,3)(5,2,4)(5,2)(2,4

14、)276491(5,4)(4,2)104454=F4(5,2,4)(5,2,6)(5,2)(2,6)2777104(5,6)(6,2)183452=F4(5,2,6)(5,3,4)(5,3)(3,4)113950(5,4)(4,3)102737=F4(5,3,4)(5,3,6)(5,3)(3,6)117182(5,6)(6,3)183957=F4(5,3,6)(5,4,6)(5,4)(4,6)107282(5,6)(6,4)185472=F4(5,4,6)(6,2,3)(6,2)(2,3)224163(6,3)(3,2)163147=F4(6,2,3)(6,2,4)(6,2)(2,4)2264

15、86(6,4)(4,2)204464=F4(6,2,4)(6,2,5)(6,2)(2,5)227092(6,5)(5,2)123951=F4(6,2,5)(6,3,4)(6,3)(3,4)163955(6,4)(4,3)202747=F4(6,3,4)(6,3,5)(6,3)(3,5)165571(6,5)(5,3)123446=F4(6,3,5)(6,4,5)(6,4)(4,5)205373(6,5)(5,4)124456=F4(6,4,5)（4）k=3时，F3(i,S3) = minDi,j + F4(j,S4) i=2,3,4,5,6X3=(i,S3)U3=(i,j)X4=(j,S4)V

16、3=Di,jF4(j,S4)V3 + F4(j,S4)(2,3,4,5)(2,3)(3,4,5)185472(2,4)(4,3,5)304272(2,5)(5,3,4)253762=F3(2,3,4,5)(2,3,4,6)(2,3)(3,4,6)186987(2,4)(4,3,6)305585(2,6)(6,3,4)214768=F3(2,3,4,6)(2,3,5,6)(2,3)(3,5,6)187290(2,5)(5,3,6)255782(2,6)(6,3,5)214667=F3(2,3,5,6)(2,4,5,6)(2,4)(4,5,6)3073103(2,5)(5,4,6)257297(2

17、,6)(6,4,5)215677=F3(2,4,5,6)(3,2,4,5)(3,2)(2,4,5)196988(3,4)(4,2,5)54752=F3(3,2,4,5)(3,5)(5,2,4)105464(3,2,4,6)(3,2)(2,4,6)197594(3,4)(4,2,6)55055=F3(3,2,4,6)(3,6)(6,2,4)156479(3,2,5,6)(3,2)(2,5,6)197897(3,5)(5,2,6)105262=F3(3,2,5,6)(3,6)(6,2,5)155166(3,4,5,6)(3,4)(4,5,6)57378(3,5)(5,4,6)107282(3,6)

18、(6,4,5)155671=F3(3,4,5,6)(4,2,3,5)(4,2)(2,3,5)325991(4,3)(3,2,5)44953(4,5)(5,2,3)84250=F3(4,2,3,5)(4,2,3,6)(4,2)(2,3,6)326092(4,3)(3,2,6)44953=F3(4,2,3,6)(4,6)(6,2,3)164763(4,2,5,6)(4,2)(2,5,6)3278110(4,5)(5,2,6)85260=F3(4,2,5,6)(4,6)(6,2,5)165167(4,3,5,6)(4,3)(3,5,6)47276(4,5)(5,3,6)85765(4,6)(6,3,

19、5)164662=F3(4,3,5,6)(5,2,3,4)(5,2)(2,3,4)275784(5,3)(3,2,4)114960(5,4)(4,2,3)103444=F3(5,2,3,4)(5,2,3,6)(5,2)(2,3,6)276087(5,3)(3,2,6)114960=F3(5,2,3,6)(5,6)(6,2,3)184765(5,2,4,6)(5,2)(2,4,6)2775102(5,4)(4,2,6)105060=F3(5,2,4,6)(5,6)(6,2,4)186482(5,3,4,6)(5,3)(3,4,6)116980(5,4)(4,3,6)105565=F3(5,3,4

20、,6)(5,6)(6,3,4)184765=F3(5,3,4,6)(6,2,3,4)(6,2)(2,3,4)225779(6,3)(3,2,4)164965(6,4)(4,2,3)203454=F3(6,2,3,4)(6,2,3,5)(6,2)(2,3,5)225981(6,3)(3,2,5)164965(6,5)(5,2,3)124254=F3(6,2,3,5)(6,2,4,5)(6,2)(2,4,5)226991(6,4)(4,2,5)204767(6,5)(5,2,4)125466=F3(6,2,4,5)(6,3,4,5)(6,3)(3,4,5)165470(6,4)(4,3,5)204

21、262(6,5)(5,3,4)123749=F3(6,3,4,5)（5）k=2时，F2(i,S2) = minDi,j + F3(j,S3) i=2,3,4,5,6X2=(i,S2)U2=(i,j)X3=(j,S3)V2=Di,jF3(j,S3)V2 + F3(j,S3)(2,3,4,5,6)(2,3)(3,4,5,6)187189(2,4)(4,3,5,6)306292(2,5)(5,3,4,6)256590(2,6)(6,3,4,5)214970=F2(2,3,4,5,6)(3,2,4,5,6)(3,2)(2,4,5,6)197796(3,4)(4,2,5,6)56065=F2(3,2,4,5,6)(3,5)(5,2,4,6)106070(3,6)(6,2,