数学建模 自习室管理

1、一问题重述：近年来,大学用电浪费比较严重，集中体现在学生上晚自习上，一种情况是去某个教室上自习的人比较少,但是教室内的灯却全部打开,第二种情况是晚上上自习的总人数比较少,但是开放的教室比较多，这要求提供一种最节约、最合理的管理方法。根据题目所给出的数据，有以下问题。数据见表。1 假如学校有8000名同学,每个同学是否上自习相互独立，上自习的可能性为0。7。要使需要上自习的同学满足程度不低于95，开放的教室满座率不低于4/5，同时尽量不超过90。问该安排哪些教室开放，能达到节约用电的目的。2 在第一问基础上，假设这8000名同学分别住在10个宿舍区，现有的45个教室分为9个自习区,按顺序5个教室

2、为1个区，即1,2，3，4，5为第1区，41，42，43,44，45为第9区.这10个宿舍区到9个自习区的距离见表2.学生到各教室上自习的满意程度与到该教室的距离有关系，距离近则满意程度高,距离远则满意程度降低.假设学生从宿舍区到一个自习区的距离与到自习区任何教室的距离相同。请给出合理的满意程度的度量，并重新考虑如何安排教室，既达到节约用电目的，又能提高学生的满意程度。另外尽量安排开放同区的教室.3 假设临近期末，上自习的人数突然增多,每个同学上自习的可能性增大为0。85,要使需要上自习的同学满足程度不低于99%,开放的教室满座率不低于4/5,同时尽量不超过95.这时可能出现教室不能满足需要，

3、需要临时搭建几个教室。假设现有的45个教室仍按问题2中要求分为9个区.搭建的教室紧靠在某区，每个区只能搭建一个教室，搭建的教室与该区某教室的规格相同(所有参数相同），学生到该教室的距离与到该区任何教室的距离假设相同.问至少要搭建几个教室,并搭建在什么位置，既达到节约用电目的，又能提高学生的满意程度。表格见附录1.需要研究的问题：1 统计出上自习的人数和所需要的座位数2 把节约用电作为问题一的约束条件求解3 根据宿舍区到自习区的距离（附录1表2）构造学生上自习满意程度的函数4 在解决问题一的基础上，同时考虑节约用电和满意程度配置开放自习教室，进行多目标规划。5 改变约束条件，重新计算上自习人数和

4、所需要的座位数6 考虑搭建若干个教室提供足够座位给期末时上自习人数,同时兼顾提高满意度和节约用电的要求二模型设计和求解：（一）模型假设：1 每个同学上自习相互独立，且概率相同2 每个同学随机选择自习教室，不受距离、楼层等因素的干扰3 计算过程中，座位数和教室数满足整数的要求4 满意度只与学生区到自习区的距离有关5 情况1：学生人数共8000人， 学生区不对总人数进行平均分配即不考虑10个学生区人数的居住分配情况情况2：10个学生区，每个区域平均配置即居住有学生800名 6 问题3中在未搭建临时教室之前10个学生区中没有座位的人数相同7 若某教室开放，则此教室所有灯管全部打开8 不考虑搭建临时教

5、室的成本问题(二）符号说明：符号含义n样本容量上自习事件的概率分布为0或1，分别表示第i号自习教室关闭或开放k上自习的可能人数为第号教室单盏灯管消耗功率（瓦/每只）为第号教室所含灯管数为第号教室消耗总功率为第号教室拥有座位数r所需座位数，满意度函数所用系数从i学生区到j自习区的满意度从i学生区到j自习区的人数从i学生区到j自习区的距离第号自习区提供的座位数可在第号自习区的上自习人数上限目标函数的加权系数Ai设宿舍区依次为A1，A2,A10Bi自习区为B1，B2，B9与第号教室规格相同的搭建教室的选择变量q所有教室总座位数所有教室总的最小功率 总的最大满意度(三）解题思路及过程：问题1基于题目情

6、况,根据题目所给的表格，运用概率统计的相关知识，分析和计算学生上自习的人数以及所需要的座位数目。然后根据节约用电的原则，把耗电最小作为教室选择的约束条件，得到结果。具体步骤如下：(1）计算所需座位数此问题符合概率统计中的二项分布。由于样本值较大，则可以用正态分布对二项分布进行近似计算.应用“棣莫弗一拉普拉斯（ DeMoivreLaplace）定理”进行样本计算。将满足程度不低于95%理解为上自习得同学有95都有自习座位坐.即每个上自习人能够正常上自习的概率为0。95。由此可以计算出上自习所需座位数.再由开放教室的满座率求得座位数的上限和下限。计算过程：样本容量n=8000, 所需座位数为r有

7、（i=1，2,8000） 表示上自习的人数，.要使得由棣莫弗一拉普拉斯( DeMoivreLaplace）中心极限定理,有查正态分布表得, 解得 由满座率介于80至90之间，求得座位上限 座位下限由所有教室总座位数, 所以座位上限是虚约束.（2)优化选择教室方法一：由于教室的选择只有两种方式:选择与不选择。顾此部分采用0，1整数规划方法.为抉择变量,有设目标函数：，即开放教室用电总功率的最小值.构造约束条件：由Lingo软件实现(见附录2)。得到结果：开放2,3,4,5，6,7，8,9，10，12，13，14，17,18，19，20，21，22,23，24，25，26，27,28,29,30,

8、31,32，34,35，36，37,38，39,40，41，42,43号自习教室； 关闭1，11,15，16，33，44,45号自习教室。由题目所提供表1确定:共提供座位6301个，消耗总功率为80577瓦。方法二:方法二采用穷举法.顾名思义，穷举法就是把所有的可能情况一一列出来，进行验算。穷举法用时间上的牺牲换来了解的全面性保证,尤其是随着计算机运算速度的飞速发展，穷举法的形象已经不再是最低等和原始的无奈之举。此题,可以通过穷举法进行计算。根据题意，要使自习教室提供的座位在6298到7085之间，则设变量设单行矩阵Z=单列矩阵C=则耗电总功率让按二进制递增到，取的最小值就是所求的最小耗电总功

9、率.现按照单位座位耗电功率进行排序，见下表：次序教室号座位数教室单位座位耗电功率12416010。12523216010。12534020011。5241421011。90552821011。90563621011。90573921011.90581319012.12692719012。126102919012.126113519012。126123719012。126133819012。126143020512.195151819512.30816319312.43517419312.435181719212.519512812。656201912812。656214318012。82291

10、1013。091232311013。091243111013。091252625613.184263425613528812013。5291012013.5302012013。5312212013。5321224713.66433712014。4342112014。4354215015。3636257015。42937337015。429384115016。667临界点：628339116416。87540447017。8574145120184228819。09143168519。765441570244516426。25其中为已知，Z矩阵具有可能情况。穷举次数过多

11、，计算机运算时间过长，因此需要减少穷举次数.观察此表,当排到如表所示的临界点时,已经包括座位数6283个,要达到6298个的座位下限，12号教室单位座位的消耗功率为13.664W,7号教室单位座位的消耗功率为14。4W,差别较大。因此,只对上表所列的后13个教室进行穷举，即次序为3345号的教室开放情况进行穷举(13个教室已经远大于临界点后的7个教室，完全可以实现对教室耗电总功率的最小配置，即这里的局部最优解就是全局最优解).此时,次序号132号的教室全部座位综合为5603个，消耗功率68557W。即次序号132号教室的全部为1。则单行矩阵Z=单列矩阵C=则耗电总功率(68557+)W，取最小

12、值。由Matlab软件编程得:ticmin=100000;for i=0：213-1 xi=dec2bin(i）； xi=xi148; l=length（xi）; xishu=zeros（1，13）； for j=（13-l+1）:13 if(j=13） xishu(j)=xi（13-j+1）； end end fenliang=120；120；150；70；70；150；64；70;120;88;85；70；64； dan=14。4；14.4;15。36；15。429；15。429；16。667；16.875；17。857；18；19.091;19.765；24；26。25; sum=695

13、57+xishu（fenliang。dan）； zuo=5603+xishufenliang； if（zuo6297&summin） min=sum； z=xishu； azuo=zuo; endendtoc解得 ,则此教室关闭。解得需要关闭的教室号为33,11，44，45，16,15,1.其余教室开放。此时共提供座位6301个，消耗总功率为80577瓦.与方法一同解。问题1两种方法的比较：方法一是在考虑选择问题时的规范解法,考虑抉择变量和约束条件即可,具有通用性和普遍性。方法二在很多领域也都可以用到，但是此题的穷举次数是相当庞大的，因此需要选取好的局部最优解来实现全局最优解.所以第二

14、种方法的重点是放在如何选取好的局部最优解来减少穷举次数。但穷举法具有的通俗性以及其通过计算机的易实现的特点是其优势。问题2方法一：使用分层序列法实现双目标规划由题目可知,问题2可以在问题1的基础上构建，再考虑分区问题和满意度问题。在此，采用多目标规划方法中的分层序列法。所谓分层序列法，就是把多目标规划问题中的p个目标按其重要程度排出一个次序，假设最重要，次之，再次之，最后一个目标为.先求出以第一个目标为目标函数。这里把节约用电作为首要目标，即以问题1作为第一目标函数.问题1解得的所开放教室为最节约用电时的选择。由于该目标函数的最优解已在问题1中求得，在此只考虑满意程度的问题（第二目标函数)，即

15、分配区域让上自习学生获得最大满意度。（1） 由表2构造满意度矩阵由题目可知,自习的满意程度只和从学生区到自习区的距离有关系，则可以利用距离关系构造满意度矩阵，为目标规划确定满意系数.此处采用线性构造方法,能够直观体现满意程度和距离的关系。由表2可以得到,学生区距离自习区最近的路程为305米,最远的路程为696米。设距离305米时，满意度为1 距离696米时,满意度为0构造线性方程组解二元一次方程得 得满意度计算公式： 由此计算公式得到满意度矩阵B1B2B3B4B5B6B7B8B9A10。8721210.0971870。808180.708440。335040。721230.531970。946

16、29A20。0025580。416880。580560。485930.670080。570330。782610。419440。23529A30。470590.358060.797950.624040。212280。317140。54220。432230.19949A40。951410。396420。961640。588240。700770。117650。997440。227620。02046A500.20460.565220.503840。792840。35550.685420.0306910。26854A60。590790.250640。739130.562660。0588240。31458

17、0.79540.153450。36829A70。874680.800510。39130。368290.634270.424550。549870。966750。98465A80.6930910。618930。314580.918160.976980。386190。39130。99744A90。994880.818410。411760.953960。636830。365730。278770。304350.92583A100。547310.56010.652170。856780.322250.296680。268540。52430.44501（2） 采用分层序列法,在问题1结果基础上考虑满意度实现多

18、目标规划设目标函数 即所有上自习学生总的满意程度条件1:情况1：不考虑10个学生区的人数平均分配情况,即可能出现一个学生区上自习的人数为0的情况.情况2：在问题1中,解得上自习人数为5668人,若分为十个学生区，则简化模型，设每个学生区上自习的人数为567人。条件2：由问题1结论可得具体的开放教室，按照题目要求分组,则每个自习区的座位数固定，见下表：自习区提供座位数提供座位数×0。9B1602542B2590531B3647582B4635572B5580522B61051946B7716644B81000900B9480432总和63015671由此表得约束条件: 即从不同的10个

19、学生区到j号自习区的人数不得高于j号自习区的人数上限条件3:由问题1的结论确定 座位上限 座位下限由Lingo软件编程得到如下结果：情况1B1B2B3B4B5B6B7B8B9A1A2A3A4582644A5A6A7900A8531522946432A9542572A10满意度：97.38结论：情况1为一极限情况:开放教室与问题1结论相同，即开放全部9个自习区,以达到最大程度的节约用电的目的。上表所示，在不考虑学生区分配居住人数的情况下，第4学生区的同学有582人到第三自习区上自习,有644人到第七自习区上自习;第7学生区的同学有900人去第八自习区上自习；第8学生区各有531，522，946，

20、432人分别去第二，第五，第六,第九自习区上自习。第9学生区各有542,572人分别去第一和第四自习区上自习。 这种情况下，学生的满意程度为最大,达到97。38。但是，这种情况与实际情况不相符。所以设计了第二种情况： 情况2B1B2B3B4B5B6B7B8B9A153136A2379188A3423144A437615932A552245A6567A7567A8567A91665396A10567542531582572522946644899432满意度:85。04结论：上表所示为假设每个学生区居住有800名学生，根据问题1可知，每个学生区有567名同学上自习.得到以下结果。第1学生区各有5

21、31人，36人分别去第二和第九自习区自习;第2学生区各有379人，188人分别去第六和第八自习区自习；第3学生区各有423人,144人分别去第三和第八自习区自习；第4学生区各有376人，159人，32人分别去第一，第三和第七自习区自习;第5学生区各有522人，45人分别去第五和第七自习区自习；第6学生区有567人去第七自习区自习；第7学生区有567人去第八自习区自习;第8学生区有567人去第六自习区自习；第9学生区各有166人,5人,396人分别去第一，第四和第九自习区自习；第10学生区有567人去第六自习区自习；此时，上自习学生的满意程度为85。04%，此种情况也是在一个理想情况下，首先考虑

22、节约用电，然后考虑学生满意程度的优先顺序下实现的优化分配.由于题目要求尽量把所关闭的教室放到同一个区内,但是在最节电情况下，开放的教室分散于9个自习区中，并且任何一个区关闭的总的座位数与所需的座位数之和都超出了自习室所能提供最大的座位数，因此不可能实现关闭一个区的情况。方法二:采用加权系数法根据题目要求，满意度只与从学生区到自习区的距离有关。构造关于最大满意度和最小消耗功率的目标函数，并从这两个函数中寻找约束关系，应用Lingo软件求解.如果每个学去自习区的路程最短，则所有学生去自习区的总路程也是最短的，同时满意度函数达到最大值.从而对总距离和总消耗的功率进行线性加权求和，通过对加权系数的调整

23、改变功耗和路程在目标函数中所占的比重,考虑约束条件限制，可以求得不同权重下的最优解。建立目标函数其中为号教室的耗电功率，为选择变量，表示从号学生区到号自习区的距离，表示从号学生区到号自习区的人数,为加权系数，其和为1，若改变权重，则意为考虑优先程度的讨论。若则考虑节约用电的目标优先于满意度的要求.越大，则节约用电的考虑程度所占比重越大，相应地,满意度的考虑程度越小：反之越大，则满意度的考虑程度所占比重越大，节约用电的考虑程度越小。如方法一中表格所示，行为学生区，列为自习区，条件约束如下：（1） 行约束：此种方法假设使用方法一中情况2的条件,由于总自习的人数为5668人,则设每个学生区有567名

24、学生上自习,由此得到约束条件:（2） 列约束:设为号自习区的座位数，为号教室的座位数，为选择变量，此条件仍然满足问题1的结论,即开放的教室满座率不低于4/5，同时尽量不超过90%。由号自习区座位数区内被选中教室座位数之和构造约束条件：通式： 即从各个学生区到第号自习区的自习人数位于自习区座位数的80到90区间内.使之符合满座率的要求.（3） 座位总数约束：提供的所有教室座位总和必须大于问题1所求的座位下限，由此构造约束条件：其中为号教室的座位数。由Lingo软件编程实现如下结果:问题2结论表线性系数关闭的的教室号功耗(w）路程（m）满意度：97。91情况一0。10。91637384483649

25、2018763宿舍到对应自习区的人数B1B2B3B4B5B6B7B8B9A153136A2373194A3567A453136A552245A699468A75589A8567A9684495A10567情况二0.30。71151637816832036487满意度：84。27%45宿舍到对应自习区的人数B1B2B3B4B5B6B7B8B9A153136A237837152A3567A44377258A552245A6567A7567A8567A91044459A10567情况三0。70。3121516807172056993满意度:69。18%4144宿舍到对应自习区的人数B1B2B3B4B

26、5B6B7B8B9A153136A2378189A3435132A426820495A552245A6567A7567A8567A91944369A10567现在通过软件编程提供的数据构造满意函数:最远总路程 =2145988最近总路程 2106045满意度 则可计算三种情况的满意度： x=2016045时，y97.91% x=2036487时,y84。27x=2056993时,y69。18% 如表中所示：情况1：时，关闭教室16，37，38，44，其余教室开放，此时满意度为97。91，总消耗功率83649瓦。情况2：时，关闭教室1,15,16,37，45，其余教室开放，此时满意度为84。27

27、,总消耗功率81683瓦.情况3：时，关闭教室1,2，15，16，41，44,其余教室开放，此时满意度为69.18，总消耗功率80717瓦。问题2两种方法的比较：方法一的分层规划实际上就是极限情况下的加权规划构造的满意函数虽然不同,但原理相似,都可以体现满意程度。加权系数法具有更大的灵活性，可以通过改变权重实现目标函数的调整。问题3（1） 重新计算上自习人数以及所需要的座位数由于题设条件改变，上自习人数增多,由问题1所得结论已经不能使用。现重新计算上自习的学生总数以及所需座位数的上下限.解题思路与问题1的步骤1相同，仍然利用概率统计的相关知识,在此不重复赘述，仅列出计算过程:样本容量n=800

28、0， 所需座位数为r有 (i=1，2，8000) 表示上自习的人数,.要使得由棣莫弗一拉普拉斯( DeMoivreLaplace）中心极限定理,有查正态分布表得, 解得 由满座率介于80%至90之间,求得座位上限 座位下限（2）已经解得，上自习的人数为6875人,所需要的座位数下限为7237个。由于要考虑搭建教室的成本问题，所以现有教室必须全部开放才可能达到最小成本. 由已知，开放所有教室提供的座位总数为6844个 所需要的座位数下限为7237个 则搭建教室所需提供的座位数下限为72376844393个在原先自习室自习的人数为6844×0.96502人 则需要在搭建自习教室自习的人数

29、为68756502373人现构造目标函数：其中：为与第号教室规格相同的搭建教室的选择变量则有为第号教室的消耗功率为从号学生区到第号教室的距离为从号学生区到第号教室上自习的学生人数现进行条件约束：（1） 座位限制:搭建教室所提供的座位下限为393个即其中为号教室的座位数.（2） 基于学生区的人数限制：由于需要在搭建自习教室自习的人数为68756502373人，所以设每个学生区都有38个学生无法去已有的自习教室自习。即（3） 基于教室人数及其满座率的限制：由于每个教室的满座率在8095%之间，搭建教室必须也满足次要有，有： 即分别从每个学生区到第个搭建教室的人数和必须符合此教室满座率上下限的要求。

30、（4） 一个区只能搭建一个教室,且搭建的教室与该区某个已有教室所有参数相同的条件限制： 现用Lingo软件编程实现如下结果：问题3结论表方案1教室数功耗（w)路程（m)57440132918学生区到对应自习区的人数（括号里数字代表与搭建教室规格相同的教室号)B1B2B3（11）B4（16）B5(25）B6B7(31)B8B9（45）A138A238A338A41424A538A638A738A81820A938A1038方案2教室数功耗（w）路程（m）46880134022学生区到对应自习区的人数(括号里数字代表与搭建教室规格相同的教室号)B1B2B3（11）B4（16)B5B6B7(32)B

31、8B9(41）A138A238A338A41424A538A638A738A838A938A1038方案3教室数功耗(w）路程（m)35400136344学生区到对应自习区的人数（括号里数字代表与搭建教室规格相同的教室号）B1B2B3B4(19）B5B6B7（32）B8B9（45）A138A238A338A438A538A638A738A838A938A1038方案4教室数功耗（w)路程(m）34680148858学生区到对应自习区的人数（括号里数字代表与搭建教室规格相同的教室号)B1B2（9)B3B4B5（24）B6B7（32）B8B9A138A2434A338A438A538A638A73

32、8A838A92810A1038表格中有四种搭建自习教室的方案,均由调整不同的线性拟合系数获得。在编程解释此问题的过程中，把满意程度直接用距离衡量。如表所示:方案1：分别在B3，B4，B5，B7,B9五个自习区构建与11，16，25，31,45号教室规格相同的教室。新搭建教室消耗功率为7440瓦.方案2：分别在B3，B4，B7，B9五个自习区构建与11，16，32，45号教室规格相同的教室。新搭建教室消耗功率为6880瓦.方案3：分别在B4，B7，B9五个自习区构建与19，32,45号教室规格相同的教室。新搭建教室消耗功率为5400瓦。方案4：分别在B2，B5，B7五个自习区构建与9，24，3

33、2号教室规格相同的教室。新搭建教室消耗功率为4680瓦。三模型评价优点：1、 本模型由于综合考虑了教室的开放与学生满意度和节约用电的关系，从而影响达到二者兼顾下的最优情况,由此可以推广至多重目标的统筹兼顾。2、 本题采用多种解题方案，给出的都是简单易行的最优原则，从不同的角度不同的侧重点来分析问题，对各种实际情况有更强的适应性。3、 充分利用MATLAB和Lingo软件来求解矩阵问题和线性规划问题，达到了实现方法多样化的目的。4、 可拓展性强。该模型的适用范围比较广，在工厂和仓库等方面都有应用。5、 联系所给出的数据,可以观察出于实际符合的很好。缺点:模型中只对极限情况进行约束，没有考虑各个教

34、室的实际满座率。四参考文献：1. 朱德通，最优化模型与实验，上海,同济大学出版社，20032. 陈宝林，最优化理论与算法（第2版），北京，清华大学出版社，20053. 贾秋玲等,基于Matlab 7。x系统仿真、分析及设计，西安，西北工业大学出版社,20064. 徐伟等，概率论与数理统计，西安，西北工业大学出版社,20025. 肖华勇,随机数学基础，西安,高等教育出版社，20056. 维普期刊五附录1. 题目所给表格表1 教室相关数据教室座位数灯管数开关数一个开关控制的灯管数灯管的功率/每只1644231440w2884231440w31934841250w419350510 48w51283

35、6218 45w612036218 45w7120364948w81203631245w91103631240w10120364945w1164273940w122477551545w131904831648w14210505105016854231440w171924841250w1819550510 48w1912836218 45w2012036218 45w21120364948w221203631245w231103631240w24160364945w2570273940w262567551545w271904831648w282105051050w291

36、904831648w302055051050w311103631240w32160364945w3370273940w342567551545w351904831648w362105051050w371904831648w381904831648w392105051050w402004831648w411505051050w421504831648w431804831648w4470255550w451204531548w表2 学生区(标号为A）到自习区(标号为B）的距离（单位：米）B1B2B3B4B5B6B7B8B9A1355305658380419565414488326A26955334

37、69506434473390532604A3512556384452613572484527618A4324541320466422650306607688A5696616475499386557428684591A6465598407476673573385636552A7354383543552448530481318311A8425305454573337314545543306A9307376535323447553587577334A104824774413615705805914915222. 问题1方法一Lingo软件编程源代码：model:sets:shu/1.45/:seat

38、，p,x;endsetsdata:seat=64,88,193,193,128，120，120，120，110，120,64，247,190，210，70，85,192,195，128,120,120，120，110，160，70,256,190，210，190,205,110，160，70，256，190，210，190，190，210，200,150，150，180，70，120；p=1680,1680,2400，2400,1620，1620,1728,1620，1440,1620,1080,3375，2304，2500，1680,1680,2400，2400,1620，1620，1728，1620,1440,1620，1080，3375，2304,2500,2304，2500,1440，1620，1080，3375，2304，2500,2304,2304,2500，2304，2500，2304，2304,1250，2160；enddatamin=sum（shu（i）：p（i)*x（i）；sum（shu（i）：x（i）*seat（i)=6298；for(shu(i）：bin(x（i)）;end问题1方法一由Lingo实现结果：Globa