DVD的在线租赁大学生数学建模论文.doc

DVD在线租赁模型摘 要本文针对DVD的在线租赁服务问题，建立了DVD的需求预测、购买和分配的多种模型。 针对第一个问题建立了一个概率模型。我们从一张DVD被租赁的期望次数去考虑问题，从而避开了从会员租碟去考虑而要分多种情况复杂境地。首先从理想情况出发，建立初步模型I，然后考虑到普遍性，给出了问题更符合实际的模型II及求解（见表5、6），得出应该至少准备6250张DVD1，才能保证希望看到DVD1的会员中至少50%在一个月内能够看到该碟；要保证在三个月内至少95%的会员能够看到DVD1，则需要至少准备3958张。对于第二问，我们首先对数据进行了分析和预处理，推断出对于无论哪种分配方案，网站提供的DVD数是不能满足全部订购会员三张DVD的要求的。接着，本文定义了个体满意程度如下：在个体满意程度定义下，给出了群体满意程度的定义。基于群体满意程度，引入了数据结构中稀疏矩阵的行优先存储方法先建立了一个0-1规划模型III，然后采用贪心算法的思想建立起一个贪心算法模型，求出并列出前30位会员得到的DVD的种类和对应会员的偏爱程度（见表7）。最后，对求解的结果进行分析，求出群体满意度为。在第三个问题的解决过程中，我们避开前两问另辟蹊径，认为一旦会员租赁到他的在线订单中的三种DVD，会员就得到了他想看的DVD。为了使会员的满意度达到最大，我们让对一种DVD满意度为1或2或3的会员都得到满足，这样会员得到的三张DVD就是自己偏爱度为1、2、3的三张。基于此再计算出每种DVD的购买数量（见表8）。在问题四中，考虑到不同的年龄阶段对不同类型的DVD的需求差异一般比较大，所以提出了会员年龄分组问题，经过分析，给出了一种抽样准则，基于这个准则和简单随机抽样性质，建立了一个新的预测模型。此外，还提出了A类和B类会员租赁DVD天数的问题，分析了会员租赁天数对群体满意度的影响，。 在模型改进里，我们把会员租赁DVD的周期设为10天或20天，让会员下订单和网站经营人员分配DVD都有较多的支配时间。而且考虑那些在分配时分到一张或两张DVD的没有分发出去的不成功的DVD分配，再重新进行一次根据偏爱程度进行分配。这样可以有效的提高会员的满意度。 该模型的特点是简单有效，思路清晰，基本能够解决题中的几个问题，但不足之处是模型假设有些理想化，而且缺少对结果验证，这些都有待改进。1问题重述随着信息时代的到来，网络成为人们生活中越来越不可或缺的元素之一。许多网站利用其强大的资源和知名度，面向其会员群提供日益专业化和便捷化的服务。DVD的在线租赁就是其中一种可行的服务。顾客缴纳一定数量的月费成为会员，订购DVD租赁服务。会员对哪些DVD有兴趣，只要在线提交订单，网站就会通过快递的方式尽可能满足要求。会员提交的订单包括多张DVD，这些DVD是基于其偏爱程度排序的。网站会根据手头现有的DVD数量和会员的订单进行分发。每个会员每个月租赁次数不得超过2次，每次获得3张DVD。会员看完3张DVD之后，只需要将DVD放进网站提供的信封里寄回（邮费由网站承担），就可以继续下次租赁。我们需要考虑如下问题：1）网站准备购买一些新的DVD，通过问卷调查1000个会员，得到了愿意观看这些DVD的人数。此外，历史数据显示，60%的会员每月租赁DVD两次，而另外的40%只租一次。网站现有10万个会员，对表1中的每种DVD来说，应该至少准备多少张，才能保证希望看到该DVD的会员中至少50%在一个月内能够看到该DVD？如果要求保证在三个月内至少95%的会员能够看到该DVD应该至少准备多少张？2）表2中列出了网站手上100种DVD的现有张数和当前需要处理的1000位会员的在线订单，如何对这些DVD进行分配，才能使会员获得最大的满意度？并具体列出前30位会员（即C0001C0030）分别获得哪些DVD。3）继续考虑表2，并假设表2中DVD的现有数量全部为0。从网站经营管理人员角度出发，应该如何决定每种DVD的购买量，以及如何对这些DVD进行分配，才能使一个月内95%的会员得到他想看的DVD，并且满意度最大？4）从网站经营管理人员角度出发，考虑并分析在DVD的需求预测、购买和分配中还有哪些重要问题值得研究？并提出问题，尝试建立相应的数学模型。2模型假设1 对1000个会员的调查结果是准确的且具有代表性，能够准确反映10万个会员的租赁愿望；2 网站管理人员完全根据当前的会员在线订单决定每种DVD的购买量，不考虑未来租赁市场的可能变化以及其他因素对管理者决策的影响；3 单个会员对网站的满意度仅与会员对网站租赁给会员的DVD的偏爱程度有关。4 一个会员不会在一次订单中租相同的DVD；5 每月月初会员发定单订购当月的DVD，由网站统一邮寄给客户，在网站第二次统一邮寄前已归还租赁DVD的会员可第二次订购当月的DVD，网站第二次统一邮寄给客户。6 会员在当月租赁的DVD都能在下个月网站发货前成功归还；3符号说明：第i种新购买的DVD数目，i=15；：网站提供的第j种DVD的总量；j=1100；：第i种DVD观众愿意观看的人数，i=15；：一张新购买DVD在一个月内被租赁出去的期望次数；：一张新购买DVD在三个月内被租赁出去的期望次数；：单个个体的满意度；i=11000；：总体的平均满意度；：会员用数字量化的对DVD偏爱程度；i=1，2，3；：表示允许的最大误差，即偏爱程度最低（8、9、10）的三张DVD与分配到偏爱程度最大（1、2、3）的三张DVD时的最大偏爱程度差值；：第i个会员第k次租赁的满意度，i=11000，k=1、2；：的权重，k=1、2；：第i 个年龄区间；：第i个年龄区间占全体会员的比例；：某网站的会员总数。4模型的建立和求解41 问题一411 问题一分析通过问卷调查1000个会员，得到了愿意观看这5种DVD的人数见表1：表1 对1000个会员调查的部分结果：DVD名称DVD1DVD2DVD3DVD4DVD5愿意观看的人数200100502510根据随机抽样的统计性质，随机样本反映了总体的特性，我们可以估计得到10万个会员的结果如表2：表2 10万个会员的估计结果：DVD名称DVD1DVD2DVD3DVD4DVD5愿意观看的人数2000010000500025001000对于客户和网站的订购与经营，基于上面的假设，我们可以做如下推断：在每月月初客户发定单订购当月的DVD，网站统一邮寄给客户，在第二次订购结束前归还已订DVD成功的客户可以参与第二次订购，网站再次统一把DVD邮寄给客户，也就是说，如果第一次订的DVD没在第二次定货结束前归还成功，网站就不予办理本月的第二次订货，这就出现了40%的会员一个月只租赁DVD一次，除此之外，网站有有效的措施使每位客户在当月末归还所有租赁的DVD。根据题意，每个会员每次租赁可获得3张DVD，这里的三张可以是上面5种新购DVD中的1种、2种、3种或是没有，这样就必须对这4种情况分别讨论，从而把问题引进了一个很复杂的境地。但是，我们转换一个角度思考问题，不是从每个会员的租赁DVD的分情况考虑，而是从上面5种新购买DVD的每一张被租赁次数去考虑，即考虑每张新DVD在每个月被租出去的次数，就可以使问题大大简化。412 模型建立和求解历史数据显示，60%的会员每月租赁DVD两次，记为A类会员；而另外的40%每月只租一次，记为B类会员。愿意观看各种新购DVD的会员也理当由A类会员和B类会员组成，且它们的比例为3：2（60%：40%），当满足的50%希望看到该新购DVD的会员在一个月内能够看到该DVD全是A类会员时，即新DVD能在当月第二次订货前归还，那么这张新DVD可以在本月再次被网站租赁出去，这样每张新购DVD在一个月内被租赁出去的次数为2，这样就可以保证用最少的DVD数量来达到50%希望看到该DVD的会员在一个月内能够看到该DVD，此时最少DVD数M1i=，具体的解如下表：表3 5种DVD购买数DVD名称DVD1DVD2DVD3DVD4DVD5购买数量500025001250625250同理可得出要保证在三个月内至少95%的会员能够看到该DVD最理想的情况是租赁该DVD的95%会员是由A 类会员60%和B类会员35%组成的，DVD租给A类会员的在三个月内可以被租出6次，而租给B类会员三个月内只能被租出3次，所以该种DVD的数量M2i与愿意观看该种DVD的人数Ni满足下面关系：具体的解如下表：表4 5种DVD购买数DVD名称DVD1DVD2DVD3DVD4DVD5购买数量37281864932466188 但是，上述情况只是最理想的状况，发生的概率非常小。每张新购买DVD如果在本月第二次订货前归还，则可被租赁两次，而在本月第二次订货前归还的客户都是A类客户，A类客户所占的比率为60%，所以该种新DVD碟的被租赁两次的概率是60%，那么DVD被租赁一次的概率就是40%。这样，一张DVD在一个月内各被租赁出去的期望次数为：a) 保证希望看到该DVD的会员中至少50%在一个月内能够看到该DVD建立目标函数： 约束条件：1）新购买的该种DVD碟的被租出的期望次数要大于愿意观看该碟的人数的50% 约束12）保证DVD的数目非负 约束23）DVD的数目为整数为整数约束3建立模型如下：（I） St 求解模型（I）得出5种DVD的具体购买数量如下表5：表5 5种DVD购买数DVD名称DVD1DVD2DVD3DVD4DVD5购买数量625031251563782391b) 要求保证在三个月内至少95%的会员能够看到该DVD 此时，一张DVD在三个月内被租赁出去的期望次数为：建立目标函数： 约束条件：1）三个月中碟被租出的期望次数要大于愿意观看该碟的人数的95% 约束12）保证DVD的数目非负 约束23）DVD的数目为整数为整数约束3建立模型如下：（II） St 模型（II）求解得出5种DVD的具体购买数量如下表6：表6 5种DVD购买数DVD名称DVD1DVD2DVD3DVD4DVD5购买数量3958197998549324742 问题二421 数据分析和预处理分析题中给出的数据，表格中每一行为一个会员的一个在线订单，每一列为DVD的种类，110的数字表示会员对某种DVD的偏爱程度，数字越小，表示会员对某种DVD的偏爱程度越高，数字0表示该会员对应的DVD不感兴趣，该DVD不在订单中。进一步分析，我们利用Excel的统计、排序和查找功能，发现所有会员的订单中非零数组合由三部分组成：（1）209个会员选择了八种DVD并根据偏爱程度进行了排序；（2）209个会员选择了九种DVD并根据偏爱程度进行了排序；（3）582个会员选择了十种DVD并根据偏爱程度进行了排序。我们又对每一种DVD的数量和对该种DVD偏爱程度为1的会员人数总和进行了统计比较发现，并非所有种类的现有DVD数量都能满足对该种DVD偏爱程度为1的会员人数总和的要求，共有14种DVD不能满足此要求。这就要求我们要对14种DVD偏爱程度为1而没有分配到该DVD的会员在进行偏爱程度为2的DVD分配时优先考虑。在计算过程中我们还发现，第37列的DVD数量为106张，而对D037有偏爱的顾客有91人，因此，D037要剩余15张，而1000个人每人分到三张DVD，总量为3000，网站现有DVD3007张，去除15张后，就不足1000人的需求，因此，不管是那种分配方案都一定有人在此次的分配中得不到自己订单中的三张DVD。此外，在对会员的C0001C1000和DVD的D000D100的编号中，我们发现这些编号有错误，他们并不是按升序排列的，也不是互不相同的，况且会员不会在订货同一时间段里发两次订单。 422 满意度的定义因为会员对某种DVD的偏爱程度都用1-10的数字量化了，数字越小，表示会员对该种DVD的偏爱程度越高，所以我们定义了如下的个体满意度函数3：定义没有分配到3张DVD的会员的。这里为会员所得到的自己订单中的三张DVD的量化的偏爱程度；为会员最理想的三张DVD的偏爱程度， 为偏爱程度最低的三张DVD，表示允许的最大误差，即偏爱程度最低的三张DVD与分配到偏爱程度最高的三张DVD的最大偏爱程度差值。我们量化这个差值，即三张DVD偏爱程度之和的差值，所以有，根据此满意度定义，为（1，2，3）时，该会员的满意度为最大，为（8，9，10）时，该会员的满意度为最小 根据个体满意度函数定义，我们又定义了如下群体满意度函数：423 模型建立和求解4231 0-1规划模型1000个会员的在线订单可以看成是一个1000行100列的稀疏矩阵，里面有大量的零，表示会员没有订购该种DVD，但我们只对里面的非零项（共9373项）感兴趣，我们是要把群体满意度最大作为目标函数，而每一个非零项都是我们的决策变量，我们要决策每一个非零项所对应的DVD是买还是不买，所以它是一个0，1变量，问题也转变成一个0-1规划模型。在这里引入数据结构里稀疏矩阵的按行优先储存方式，只记录下非零项的信息，把它储存在一个9373行3列的矩阵B里，B的第一列是按行的顺序依次存储非零项在订单矩阵的所在的行，第二列是存储该非零项在订单矩阵的所在的列，第三列是存储它在订单矩阵中的值，即它的偏爱程度。分析订单，还可以得到订单矩阵中每一行的非零项在矩阵B的位置，在这里记为temp(i)。基于上，就能建立一个0-1规划模型目标函数：约束条件：1 设对应B矩阵一行的0-1变量为，i=1,2.9372，每个会员的租碟总数不超过3，第i个会员在矩阵B里的顺序是从temp(i)行到temp(i+1)-1行： 2 设为网站提供的第j种DVD的总量，所有会员租赁同一种DVD的总量不超过网站提供的该种DVD的总量； 3 9373个决策标量的0-1约束。建立了一个0-1规划模型：（III）st但是发现目标函数很难用线形的表达式表达出来，而且过多变量都极大的影响的模型的求解可行性。针对该问题，采用了一个简易的贪心算法建立了一个DVD分配模型，来寻找分配的局部最优解。4232 贪心算法模型对于当前1000位会员的在线订单中，有582个会员选择了十种DVD并根据偏爱程度进行了排序，为了让在DVD进行分配时会员的满意度最大，我们可以采取分十步走的策略。即：在每位会员满足分到三张DVD的前提下，所有种类的DVD首先满足对其偏爱程度最高的会员，即：首先满足对其偏爱程度为1的会员，因为当前状态下，这种选择是局部最好的选择，如果该种DVD在满足偏爱程度为1的会员后有剩余则再满足偏爱程度为2的会员，如此依次类推，直至该种DVD被全部分配出去。这种分配只依赖于前面做出的选择，但决不依赖于将来所做出的选择 2。如果该种DVD数量或被分配后剩余的该种DVD在分配中不能满足同一偏爱程度的人数的要求，则根据会员的在线定单中的编号按顺序分配。但是考虑到对于相同偏爱程度下没有分配到DVD的会员，需要对贪心算法做一些改进，在对其紧接下来的较高偏爱程度的DVD分配时给予优先考虑这部分会员（具体程序见附录一）。基于上面思想，我们得出了1000个会员中订单被满足的人数为792人次，还得出了订单被满足的会员所分配到的偏爱的DVD的编号以及对该种DVD的偏爱度。下表是前30个会员分配到的DVD种类及对该种DVD的偏爱程度：表7前30个会员的分配方案会员编号DVD种类及相应偏爱度DVD种类及相应偏爱度DVD种类及相应偏爱度C000181822983C0002614420C000380150243C000471182413C0005661682113C0006191532163C000781182263C000871100C00095311002783C0010601552853C0011591632193C00123112273C0013961782213C0014521232896C0015131853669C0016841972559C0017671472513C0018411602783C0019841862663C0020451892613C002153145224C0022571552383C0023951292813C0024761412374C002591692813C0026221682953C0027581223504C0028813420C0029551302443C0030621372985得出了1000个会员中订单被满足的人数和订单被满足的会员所分配到的偏爱的DVD的编号以及对该种DVD的偏爱度后，利用我们定义的个体满意度函数求出每个（程序见附录），进而根据我们的分配方案求出所有会员的群体满意度= 0.7521 ，这个满意度还是比较满意的，他能够满足792人的3张碟的要求，如表7中所示，表格中DVD种类为0的就说明不能满足三张碟，上面我们已经分析过了，这种情况一定会发生，这样的会员网站都不予以发货。43 问题三431 模型分析与建立 对于问题三，我们还需要做如下说明：一个月内要租赁两次的会员得到DVD后可在随后的115天的任何一天归还；一个月内要租赁一次的会员得到DVD后可在随后的130天的任何一天归还。网站对于第二次收到的订单，再从已经返还的DVD中根据客户的偏爱度进行分配，其中我们约定了A类会员和B类会员租赁DVD的天数为一个区间的均匀分布。我们认为一旦会员租赁到他的在线订单中的三种DVD，会员就得到了他想看的DVD。如果会员想要租赁两次DVD，而网站只能满足他一次的需求，他的满意度将降低。假设会员的第一次租赁都要在每个月的第一天满足，为了满足95%的会员对DVD的需求，网站至少应购买张DVD。考虑使单个会员的满意度最大，网站应租赁给他偏爱程度为1、2和3的DVD。如果网站购买DVD，使得在第一次租赁中95的会员每个会员的满意度都达到最大，那么在第一次租赁中95的会员对网站的满意度就达到了100。这样对于每一种DVD，如果网站购买的量仅仅用来而且刚好满足前950个会员（也可以是其他的任意950个会员）中对该种DVD偏爱程度为1、2或3的会员，那么网站就使一个月内95的会员得到了他想看的DVD，且满意度最大。综合以上两点，网站购买2850张DVD，就既能使购买量最小，又能使第一次租赁中95%会员的满意度最大。网站管理人员决定每种DVD购买量的方法是：表示对该种DVD的偏爱程度为i的会员总数在得到DVD购买量的前提下，我们考虑使需要租赁第二次的会员满意度最大。前面我们定义的单个会员的满意度函数只代表了一次租赁的满意度，租赁二次的满意度定义为其中由于会员希望尽快看到DVD，因此第一次租赁满意度的重要性应大于第二次，应大于0.5。这里我们取为0.75。432模型求解具体的求解流程见下图：nyynyyynnn开始产生1000个随机数赋给hyh(i)<0.6,time(i)取115，sign(i)=2;h(i)>=0.6,time(i)取130，sign(i)=1k=1i=1time(i)=time(i)-1time(i)=0|y(i)=-1?第i个会员归还DVD，DVD库存增加（time(i)=0）sign(i)=2&(y(i)=0|y(i)=-1)?记录当前能租到DVD的种类record和总数sumsum>=3?借到偏爱程度高的三张DVD，y(i)=1y(i)=-1i=i+1k=k+1i<=950?k<=30?计算满意度计算各种DVD的购买量y(i)=0我们再次利用贪心算法的原理（见附录二），让每步的选择是当前最优的选择的情况下，得出100种DVD的购买量如下表：表8 100种DVD的购买量D001D002D003D004D005D006D007D008D009D01020342636192729333423D011D012D013D014D015D016D017D018D019D02027292729253728242637D021D022D023D024D025D026D027D028D029D03032293521282822182540D031D032D033D034D035D036D037D038D039D04029333028343321272627D041D042D043D044D045D046D047D048D049D05049342634342531223032D051D052D053D054D055D056D057D058D059D06038263126313029283133D061D062D063D064D065D066D067D068D069D07025293033302528333229D071D072D073D074D075D076D077D078D079D08033322328242220283028D081D082D083D084D085D086D087D088D089D09028172119332034242326D091D092D093D094D095D096D097D098D099D10039272221382433311834 在这里求解整体满意度时，我们把第二次租赁DVD的会员对DVD的偏爱程度用他在第一次订单中排序但没有分配到的DVD的种类来计算，通过1000次的计算求得群体满意度均值为88.127%，而在第二次的分配中，A类会员达到的整体满意度均值为87.94%。 4 4问题四441 会员年龄分组问题首先在DVD需求预测方面，因为我们必须在每月月初第一次客户订碟之前就购买好本月的碟，所以各种DVD碟的预测就显得格外重要，预测多于实际需求过多，购买的DVD的利用率就低，在会员费固定不变的情况下，网站的收益就变少，反之，预测过少，就会导致供不应求，这就容易导致会员群体的满意度降低，这样的潜在亏损就可能更大。在本题中的实际需求预测是通过问卷调查获得的，题中没有考虑网站会员的年龄组成，而且在确定会员总体的想看某种DVD的数量时只是由样本调查数据线形扩大而成的。其实，会员的年龄组成是一个很值得研究的问题，不同的年龄阶段对不同类型的DVD的需求差异一般比较大，我们可以我们在选择问卷调查时应该把它考虑进去，建立了一个需求预测模型如下：第一步：划分年龄区间按照会员年龄的不同将全体会员划分成K个互不相交的年龄区间：，i=1K求出不同区间各占的比例为；第二步：确定抽样准则因为网站会员的总数较大，设为，想通过问卷来调查网站所有的会员并不现实，所以我们只好通过确定抽样样本的问卷调查来估算整体特性，先确定抽样准则，假设抽样样本数为，则抽样的准则为从K个年龄区间分别随机抽取样本，样本数分别为 ，i=1K；第三步：推出会员总体的愿意观看每种DVD人数处理收回的问卷调查，可以得到第种DVD的各个年龄组的愿意观看人数，再根据简单随机抽样的统计性质，会员总体的第i个年龄组愿意观看某种DVD的人数是，所以对于会员总体，第j种碟的会员观看人数为442 会员租赁DVD天数问题在求解问题三DVD分配时时，我们约定了A类会员和B类会员租赁DVD的天数为一个区间的均匀分布，这有可能与实际情况不符合。因为每个会员租赁DVD的天数对DVD的利用率有影响，会员租赁天数越长，DVD的利用率就越低，会员的群体满意度也会越小，所以，会员租赁DVD的天数也是一个值得研究的重要问题。我们可以通过随机抽样进行问卷调查，分别统计出A类会员总体和B类会员总体的DVD租赁天数的分布规律，从而作为在第三问中的租赁天数随机模拟的分布，这样就更符合实际，从而更正确指导DVD的分配。5、模型改进方向1） 在问题二的模型建立中，我们还可以考虑把客户下订单的时间取消约束，也就是说，在客户是会员的时间内可以归还租赁的DVD的情况下，客户在一个月的任何时间内下订单。而且一个月内要租赁两次的会员得到DVD后可在随后当月的115天的任何一天归还；一个月内要租赁一次的会员得到DVD后可在随后的130天的任何一天归还。网站对于第二次收到的订单再从已经返还的DVD中根据客户的偏爱度进行分配。此外，我们还可以考虑把每月租赁DVD一次或两次的周期设置为10天或20天，而不是理想的半个月和一个月，这样，客户发订单和网站统一邮寄DVD就有很多的机动时间可以安排。2） 我们用贪心算法建立的DVD分配模型中，有一部分的订单没能满足，但他们也分配到了一个或两个的DVD，如果把这部分的DVD也用于能够满足其订单的会员的DVD的分配中，就有可能使客户的整体满意度有所提高。我们在这里考虑的满意度为总体满意度，我们还可以先从每个会员出发来考虑问题，先搜索出对每个会员的分配策略，再从中选择出个体满意度最大的分配方案。6、模型评价模型中我们用到了贪心算法，贪心算法总是做出在当前看来是最好的选择，也就是说贪心算法并不是从整体上考虑最优，它所做出的选择只是在某种意义上的局部最优解。它只会以来前面的决策做出当前的决策，而从不以来将来的决策选择。但它的最终结果却是最优解的很好近似。基于贪心算法的不足，我们又引入了优先考虑的算法原则，在一定程度上弥补了贪心算法的缺陷。在模型三的建立和求解中，我们变换了模型二中会员订购DVD的方式，在模型三中求解B类会员第二次的满意度时，我们直接用其在第一次订购时没有得到的偏爱度为4、5、6的DVD为该会员的最高满意度。这里的假设虽然与模型二中的有些矛盾，但当最佳偏爱度达到满足时，位于其次的偏爱度占据最佳位置，基于此点考虑，模型三的假设还是有很大的可取性的，而且它使得模型三的计算简化了很多。参考文献：1 徐玖平 胡知能 运筹学 （第二版） 北京 科学出版社 2004年5月2 王晓东 算法设计与分析 107页 北京 清华大学出版社 2005年2月3 刘明广 杜刚 基于满意度最大的群决策方法，,2005-9-17附录附录一% 贪心算法程序实现load('data1');data=data1;e=zeros(1000,1);c=zeros(1000,1);d=zeros(1000,3);for m=1:10 for j=1:100 for i=2:1001 if (data(1,j)>0) &(data(i,j)=m) &(c(i-1)<3) c(i-1)=c(i-1)+1; d(i-1,c(i-1)=j; data(1,j)=data(1,j)-1; end; end; endend;附录二sumecnt=0;sumecnt2=0;for i=1:1000 ecnt,ecnt2=dz(); sumecnt=sumecnt+ecnt; sumecnt2=sumecnt2+ecnt2;end;e=sumecnt/1000;e2=sumecnt2/1000;function ecnt,ecnt2=dz()load('data1.mat');goumai=zeros(100,1);left=zeros(100,1);rent2=ones(1000,3);y=zeros(1000,1);r2=zeros(1000,6);content=ones(1000,1);s=0;sumcnt=0;s2=0;sumcnt2=0;for j=1:100 for i=2:951 if data1(i,j)=1 | data1(i,j)=2 | data1(i,j)=3 goumai(j)=goumai(j)+1; end; end;end;%goumai（j）为第j种DVD的购买量hy=rand(1000,1);for i=1:1000 if hy(i)<0.6 sign(i)=2; time(i)=int8(rand()*14)+1; else sign(i)=1; time(i)=int8(rand()*29)+1; end;end;%time(i)表示第i个会员的归还DVD的时间for k=1:30 for i=1:950 time(i)=time(i)-1; if time(i)=0 | y(i)=-1 for j=1:100 if y(i)=0 & (data1(i+1,j)=1 | data1(i+1,j)=2 | data1(i+1,j)=3) left(j)=left(j)+1; end; end; for j=1:3 if y(i)=1 left(rent2(i,j)=left(rent2(i,j)+1; end; end; %第i个会员归还DVD后，DVD库存总量的增加 %y(i)=1表示归还DVD后，又借到三张；y(i)=-1表示归还DVD后，还没有再借 if sign(i)=2 & (y(i)=0 | y(i)=-1) sum=0;record=zeros(1,7); for j=1:100 if left(j)=0 | data1(i+1,j)=0 | data1(i+1,j)=1 | data1(i+1,j)=2 | data1(i+1,j)=3 ; else switch data1(i+1,j) case 4 record(1)=j; case 5 record(2)=j; case 6 record(3)=j; case 7 record(4)=j; case 8 record(5)=j; case 9 record(6)=j; otherwise record(7)=j; end;