饮酒驾车模型-数学建模论文及论大学生写作能力

饮酒驾车模型摘要本模型对饮酒驾车问题的研究，根据饮酒方式的不同，分别给出了两个血液中酒精含量的微分模型，模型Ⅰ模型Ⅱ对题目中的参考数据和网上资料，利用剩余法，最小二乘法及应用软件得出几组估计值。选取两个参数估计值，对文中的各个问题都得出较好的结论。关键词：酒精浓度模型最小二乘法剩余法时间1：问题重述据报载，2003年全国道路交通事故死亡人数为10.4372万，其中因饮酒驾车造成的占有相当的比例。针对这种严重的道路交通情况，国家质量监督检验检疫局2004年5月31日发布了新的《车辆驾驶人员血液、呼气酒精含量阈值与检验》国家标准，新标准规定，车辆驾驶人员血液中的酒精含量大于或等于20毫克／百毫升，小于80毫克／百毫升为饮酒驾车（原标准是小于100毫克／百毫升），血液中的酒精含量大于或等于80毫克／百毫升为醉酒驾车（原标准是大于或等于100毫克／百毫升）。大李在中午12点喝了一瓶啤酒，下午6点检查时符合新的驾车标准，紧接着他在吃晚饭时又喝了一瓶啤酒，为了保险起见他呆到凌晨2点才驾车回家，又一次遭遇检查时却被定为饮酒驾车，这让他既懊恼又困惑，为什么喝同样多的酒，两次检查结果会不一样呢？请你参考下面给出的数据（或自己收集资料）建立饮酒后血液中酒精含量的数学模型，并讨论以下问题：1.对大李碰到的情况做出解释；2.在喝了3瓶啤酒或者半斤低度白酒后多长时间内驾车就会违反上述标准，在以下情况下回答：酒是在很短时间内喝的；酒是在较长一段时间（比如2小时）内喝的。3.怎样估计血液中的酒精含量在什么时间最高。4.根据你的模型论证：如果天天喝酒，是否还能开车？ 5.根据你做的模型并结合新的国家标准写一篇短文，给想喝一点酒的司机如何驾车提出忠告。参考数据1.人的体液占人的体重的65%至70%，其中血液只占体重的7%左右；而药物（包括酒精）在血液中的含量与在体液中的含量大体是一样的。2.体重约70kg的某人在短时间内喝下2瓶啤酒后，隔一定时间测量他的血液中酒精含量（毫克／百毫升），得到数据如下：时间(小时)511.522.533.544.55酒精含量306875828277686858515041时间(小时)678910111213141516酒精含量38352825181512107742：模型假设和符号说明人体分为体液室和血液室，体液和血液体积都不变，分别为，；酒精从体液室向血液室转移速率与该室的酒精含量成正比，血液中酒精的排除率与血液室的酒精含量成正比；体液室与血液室酒精浓度含量大致相同；本模型中啤酒的容积，酒精含量为；5、为某参数估计值；6、（）表示为酒精浓度；7、为一瓶啤酒的酒精质量；8、为时刻血液中酒精剩余浓度。3：模型的建立与分析模型分析：根据药物动力学原理，建立单一的房室模型，把人体的血液流动系统看作一个血液室，其它体液系统作为吸收室，模型示意图如下：体液室血液室排除 建立模型 ，根据以上分析，得快速饮酒模型Ⅰ如下： 模型Ⅰ 较长时间饮酒（间隔时间为）模型假设该段时间内饮酒方式是匀速，平均速率为，并且是以同样的速率进入血液室模型Ⅱ4：模型求解快速饮酒型在模型Ⅰ中，短时内将酒精全部吸入血液室，取解得实验表明，通常情况,即酒精的吸收速率一般要大于酒精的消除速率，此时不考虑这种情况.若假定吸收率不等于排除率，模型简化为模型Ⅰ’： 较长时间饮酒模型Ⅱ:模型参数估计快速饮酒模型的参数估计在一般情况下，人体血液对酒精吸收率大于其代谢率，可假定，当时间足够大时，有，此时模型Ⅱ=，从表选取部分数据，即（），利用最小二乘法来拟合，求得，，求得剩余酒精浓度函数，利用剩余法来求相应的，对表的前行数据按不同时间不同的酒精含量分三段，具体表。求得每段吸入速率分别为，，，用加权平均求得。（详见表1）因此，根据参考数据得出拟合函数：由表1最后一列算出的观测值与表中的数据走势一致,但峰值有一定的误差，对此我们再利用更多数据进行拟合，取用同样的方法求得，，。相应的拟合函数：用编写画图程序如下：Plot[{146*(Exp[-0.214*t]-Exp[-0.7*t])},{t,0,16}]的图示从图示可知，前面的数据拟合值偏低，但后面数据拟合值较好。Plot[{121.51*(Exp[-0.196*t]-Exp[-1.57*t])},{t,0,16}]的图示从图示可知，前面的数据拟合值较好，但后面数据拟合值便高。 :模型的检验 6.1问题1的解答在参考数据表中，以简单平均计算，大李喝两瓶酒后，到凌晨点时，大李的血液酒精度应为，显然不符合事实。我们以上得出每一对的估计值为所给参考数表中的拟合函数的参数，用其来计算大李当时的酒精含量显然也不合事实。故拟合的较好的那组数据不符合实际，此时的吸入酒精的速率是前6个数据剩余浓度得出估计量，显然数据越多越反映实际吸收速率而且从网上得知，人体的代谢率一般要小于，我们取值。对此做如下的猜想，人的酒精正常代谢率一般为，因此，实际上的应在与我们拟合出来的之间，用试探法取其平均值。即求出,即得.假设大李的重量为，则血液重量约为:，，血液的体积为：，代入上式得，即大李喝一瓶酒到晚上点时，故符合新的标准。下午点又喝一瓶，到凌晨点，其,酒精浓度内，属于饮酒驾车，6.2问题2的解答6.2.1在快速饮酒的情况下，假定是大李喝下瓶啤酒，，利用模型Ⅱ，将，，，代入，，令，利用计算如下FindRoot[{150.03*(Exp[-0.15*t]-Exp[-0.7*t])==20},{t,0.1}]得,FindRoot[{150.03*(Exp[-0.15*t]-Exp[-0.7*t])==20},{t,4}]得。结论：在较短时间内喝酒，将在喝完后的16.2分钟到13.43小时内是属于饮酒驾车。6.2.2在较长时间饮酒的情况下，（以2小时为例）假定大李是喝下3瓶啤酒，，假定吸入体液室的酒精速率是匀速的，设在小时内他分次将瓶酒喝完的。即每次用分钟喝了的酒精，即得到他吸入体液室内的酒精速率，，，，代入模型Ⅲ得从该函数的图形易知：令，用编程求解如下：FindRoot[{392.93*(1-Exp[-0.15*t])==20},{t,0.2}]解得：，即他还没有和完酒，故舍去。FindRoot[{392.93*(1-Exp[-0.15*2])*Exp[-0.15*(t-2)]==20},{t,4}]解得：结论：在较长时间内喝酒，喝完后12.85-2=10.85小时以内开车属于饮酒驾车。6.3问题3的解答6.3.1快速喝酒，利用函数的极值法，对模型Ⅱ进行求导数得峰值，，只与血液室的吸收率和血液室的排除率有关，将估计值，代入，得出小时。6.3.2慢速喝酒，从问题4.2.2中看出，为增函数，为减函数，在两个小时后达到峰值。假定在小时喝完3杯啤酒，所以在第个小时达到峰值。结论：快速喝酒时血液中的酒精含量在小时后达到最高，慢速喝酒时血液中的酒精含量在小时后达到最高。6.4问题4的解答6.4.1假设大李每天在中午点和下午点各饮酒一瓶，每隔小时测一次他的酒精含量，在天后他的酒精含量为一瓶啤酒的体积为，将估计值，，，代入，解得，令趋向于无穷，得出，所以在这种情况下，已经超出饮酒驾车的标准。6.4.2进一步讨论，大李在同样的假设条件下，如果每天每次只喝（瓶）或更少，这种情况下就不会超出饮酒驾车的标准。7问题5：为你我的幸福,大家共同来关注饮酒驾车有句俗话：无酒不成宴席,千百年来，酒在中国的饮食业中形成一种文化，源远流长。家有婚庆、生子、上大学等乐事，会聚朋友欢歌畅饮；这事本身并无不对，但过量贪杯却使美酒变了味儿。若不加节制的狂酗滥饮，不仅对身体无益，耽误事，出洋相,丢面子,更甚者会给你的一生带来刻骨铭心的祸事。据报载，2003年全国道路交通事故死亡人数为10.4372万，其中因饮酒驾车造成的占有相当的比例。针对这种严重的道路交通情况，国家质量监督检验检疫局2004年5月31日发布了新的《车辆驾驶人员血液、呼气酒精含量阈值与检验》国家标准，新标准规定，车辆驾驶人员血液中的酒精含量大于或等于20毫克／百毫升，小于80毫克／百毫升为饮酒驾车（原标准是小于100毫克／百毫升），血液中的酒精含量大于或等于80毫克／百毫升为醉酒驾车（原标准是大于或等于100毫克／百毫升）。大李在中午12点喝了一瓶啤酒，下午6点检查时符合新的驾车标准，紧接着他在吃晚饭时又喝了一瓶啤酒，为了保险起见他呆到凌晨2点才驾车回家，又一次遭遇检查时却被定为饮酒驾车，这让他既懊恼又困惑，为什么喝同样多的酒，两次检查结果会不一样呢？又是如何控制好酒量以及超标准的时间段?很多司机因为应酬多，经常被这些问题捆扰着，而又无法合理地解决。在这里我们提供数据和模型，对想喝一点酒的司机提供参考。在问题4中，进一步得出结果，一个人每天在中午12点和下午6点各喝小于或等于0.6瓶，可以天天非饮酒开车。我国法律明文规定：醉酒的人犯罪，应当负刑事责任。酒精在为人们生活助兴的同时，也无情地吞噬着人们的心灵。为了您和家人的节日欢乐、幸福，千万别踏这根“高压线”啊！7：模型的优缺点与改进方向本建模思想易理解，结果结合实际，有指导意义。2、本文在解决长时间饮酒方式的时候，假设酒精进入血液为匀速的，速率有待改进。3、参数的拟合值不是很理想，有待于改进。8：附表110.250.0625303.40.8511.5620.50.25684.222.1117.8130.750.5625754.323.2418.66411824.414.4119.4551.52.25824.416.6219.4624774.438.6818.8472.56.25684.2210.5517.81839684.2212.6617.8193.512.25584.0614.2116.4810416513.9315.7215.44114.520.25503.9117.615.2912525413.7118.5513.7613636383.6421.8413.2514749353.5624.912.6715864283.3326.6411.0916981253.2228.9810.3681710100182.8928.98.3521811121152.7129.817.34441912144122.4829.766.152013169102.3029.95.29211419671.9527.33.8221522571.9529.253.8231625641.3122.241.932表1取后8项用剩余法及最小二乘法拟合得排除速率。1138.4108.44.691.17222131.2717.23124.4215.24117.935.93.583.5811.825105.9610.565.88.41785.517.52.867.168.2876.524.71969.0411.061062.0311.06147.8311114.8813611.143.461432.6-2.41526.4-1.61621.3-3.71717.2-0.81813.9-1.11911.2-0.8209.04-0.96217.30.3-1.2-16.86225.89-1.11234.80.8-0.22-3.57参考文献：[1]刘来福，曾文艺，数学模型与数学建模，北京：北京师范大学出版社，。[2]李天然，工程数学，湖南：高等教育出版社，[3]杨珏何旭洪赵昊彤，应用指南，北京：人民邮电出版社。[4]car,net/news/show/shownews.asp?nid=72D题之一（全国一等奖）公务员招聘模型河池学院，莫明丽黎川韦继乐指导老师：邹永福摘要：本文中，我们将公务员的招聘问题看作是一个分派问题，通过对应聘者的笔试成绩等级化和各方面能力的量化，找到了有利于发挥应聘者个人能力和特长并满足各部门需求的目标函数，建立了公务员招聘问题的线性规划模型:s.t.且讨论了因笔试成绩的比重变化而导致招聘结果变化的几种情形，得出了在招聘过程中，我们必须给出各能力评价因素的限定，这样才能更公平合理的实现招聘。关键词：决策变量，目标函数，分派问题，线性规划一、问题的重述与分析1.1问题的重述目前,我国招聘公务员的程序一般分三步进行：公开考试（笔试）、面试考核、择优录取。现有某市直属单位因工作需要，拟向社会公开招聘8名公务员，具体的招聘办法和程序如下：（一）公开考试（笔试）。考试科目有：综合基础知识、专业知识和“行政职业能力测验”三个部分，每科满分为100分。根据考试总分的高低排序按1:2的比例（共16人）选择进入第二阶段的面试考核。（二）面试考核：主要考核应聘人员的知识面、对问题的理解能力、应变能力、表达能力等综合素质。按照一定的标准，面试专家组对每个应聘人员的各个方面都给出一个等级评分，从高到低分成A/B/C/D四个等级，具体结果见赛题附录中的表１所示。（三）由招聘领导小组综合专家组的意见、笔初试成绩以及各用人部门需求确定录用名单,并分配到各用人部门。该单位拟将录用的8名公务员安排到所属的7个部门，并且要求每个部门至少安排一名公务员。这7个部门按工作性质可分为四类：(1)行政管理、(2)技术管理、(3)行政执法、(4)公共事业。见赛题附录中的表2所示。招聘领导小组在确定录用名单的过程中，本着公平、公开的原则，同时考虑录用人员的合理分配和使用，有利于发挥个人的特长和能力。招聘领导小组将7个用人单位的基本情况和四类工作对聘用公务员的具体条件的希望达到的要求都向所有应聘人员公布（见赛题附录中的表2）。每一位参加面试人员都可以申报两个自己的工作类别志愿（见赛题附录中的表1）。请研究下列问题：(1)如果不考虑应聘人员的意愿，择优按需录用，试帮助招聘领导小组设计一种录用分配方案；(2)在考虑应聘人员意愿和用人部门的希望要求的情况下，请你帮助招聘领导小组设计一种分配方案；(3)你的方法对于一般情况，即N个应聘人员M个用人单位时，是否可行？(4)你对上述招聘公务员过程认为还有哪些地方值得改进，给出你的建议。1.2问题的分析公务员招聘通过公开考试，面试考核，择优录取三步进行。即应聘者参加笔试后，根据考试总分的高低按一定比例进入第二阶段的面试考核。面试专家组对应聘者的能力给出一个等级评分后，从高到低分成A/B/C/D四个等级。最后由招聘领导小组按笔试成绩，面试成绩和各部门需求择优确定录用人选，并分配到各个部门。因此公务员的招聘问题可看作是一个分派问题。我们把应聘者的笔试成绩看作是应聘者的一种能力的表现，笔试能力、知识面、对问题的理解能力、应变能力、表达能力这五项可作为衡量应聘者能力强弱的因素，但由于笔试能力给定的是数字成绩，而其它各项能力给定的是符号等级，因此我们需要把笔试成绩的等级化（见附录1）。通常笔试成绩在招聘过程中占有一定比重，因此我们在该模型中考虑笔试成绩占的比重为，而其它四项各占。这样有便于我们对考试成绩的所占比重的考虑。另外，若不考虑笔试成绩，只要即可；若只考虑笔试成绩，只要。二、模型的假设(1)、假设应聘者笔试成绩及专家组对应聘人员的知识面，对问题的理解能力，应变能力，表达能力的评定都是合理正确的(2)、假设应聘者的能力和特长的发挥，主要根据其被评定了的各方面能力等级来确定。即应聘者在知识面、理解能力、应变能力、表达能力各方面的等级比其所去的部门相应方面的要求能力等级越高，则应聘者更能发挥其能力和特长。(3)、假设各部门对应聘者应试能力(即笔试能力)的要求是相同的。(4)、假设应聘者的知识面、理解能力、应变能力、表达能力这四方面在招聘的过程中处于相同的地位，即比重相同。三、符号定义i——第i名应聘者j——第j个部门k——应聘者i能力的第k项或部门j对应聘者要求的能力的第k项——表示应聘者i的第k项能力，其中分别表示应聘者在知识面，对问题的理解能力，应变能力，表达能力和笔试能力，。——表示各部门j对应聘者要求的第k项能力，其中分别表示各部门对应聘者的知识面，对问题的理解能力，应变能力，表达能力和笔试能力的要求能力，。——笔试成绩在招聘过程中占的比重。——应聘者的其它四种能力在招聘过程中各占的比重，且。四、模型的建立由于公务员的招聘问题可看作是一个分派问题，因此可用0～1变量表示分派的决策变量，令应聘者被分配到某部门工作，其个人能力与该部门的要求能力必然存在能力差距。由于题中给出的应聘者各方面能力等级和各部门对应聘者各方面能力的要求等级都是A、B、C、D四个等级，为了更好计算能力的等级差，我们需要对能力等级的量化（见附录2）。我们把应聘者i在某几方面个人能力超出部门j要求的能力值之和，定义为剩余能力，即。同时也把应聘者i在某几方面个人能力低于该部门j要求的能力值的绝对值之和，定义为欠缺能力，即。因此应聘者i被分派去部门j工作产生的过剩能力的总和，产生的欠缺能力总和为。应聘者i被分配到部门j工作，如果其各方面能力等级比部门j相应方面的要求能力等级越高则更能发挥其个人特长和能力，即越大，越大，也就越能发挥i的能力。当然有时可能出现应聘者i去部门j的很大，但是应聘者i可能有某些项能力比部门j相应的要求低很多，这样就会阻碍应聘者去部门j后的个人能力的发挥，因此我们还应要求应聘者i去部门j的尽可能的小，即尽可能的小，应聘者i的各项能力与部门j相应的要求相差不大。于是我们得到了公务员招聘的两个需要满足的目标函数：,。这两个目标函数可化为下面的单目标函数：。（1）由于该市将拟录用8名公务员安排到所属的7个部门，并要求门各部门至少安排一名公务员，所以目标函数（1）应满足下面的条件于是公务员招聘问题等价于下面的线性规划模型：s.t.（2）五、模型的求解5.1回答问题1如果不考虑应聘者的意愿，择优按需录用。成绩在招聘过程中所占的比重不同，招聘结果也有所不同，因此我们考虑成绩所占比重为：1)，2)，3)，4),5)，6)，7)。结果具体如下：表1成绩所占比重为的招聘方案人员12345678部门71123456表2成绩所占比重为的招聘方案人员12345689部门72435612表3成绩所占比重为的招聘方案人员12345689部门27163254表4成绩所占比重为的招聘方案人员124568912部门74251637表5成绩所占比重为的招聘方案人员124578912部门74257631表6成绩所占比重为的招聘方案人员124578912部门73761542表7成绩所占比重为的招聘方案人员12489121516部门22534671问题1的替换原则：由于部门2与部门3对应聘者各方面的能力要求相同，所以其能力差距也相应同，因此上面的分配方案中的部门2可用部门3替换，同理分配方案中的部门3可用部门2替换。类似的部门4与部门5，部门6与部门7也可以互换。上面情形是在不考虑应聘者志愿的前提下，招聘方案随成绩所占比重不同而不同。若想得到其它不同成绩比重下的招聘方案，只要通过附录3中给出的程序，先用Matlab程序段得到相应成绩比重下的系数矩阵，再通过Lingo程序算出相应下的招聘方案，然后得到的方案中若有满足问题1的替换原则中情形，用相应的部门进行替换。从上面的各表数据得出：当成绩所占的比重在区间内录用的人员不变，在内录用的人员不变，在点，，时录用的人员各不相同，发生变化。5.2回答问题2考虑应聘者的意愿，择优按需录用。这时相当于应聘者i不被分配到没有报意愿的部门j工作，即应聘者与该部门的决策变量，我们只要把这些为零的决策变量作为约束条件加入到(2)即可得到考虑应聘者意愿的招聘问题的线性规划方法。同时我们仍考虑成绩所占比重为：1)，2)，3),4)，5)的情形。结果具体如下：表8成绩所占比重为的招聘方案人员12345678部门54162673表9成绩所占比重为的招聘方案人员12345689部门54263671表10成绩所占比重为的招聘方案人员12345689部门54372761表11成绩所占比重为的招聘方案人员124578912部门25436716表12成绩所占比重为的招聘方案人员124578912部门35621747同样的，部门2可与部门3、部门4与部门5、部门6与部门7可以互换。其它成绩所占比重下招聘方案也可以通过附录3中的程序得到，从上面的结果中，我们发现：当成绩所占的比重在区间内录用的人员不变，在内录用的人员不变，在点时录用的人员与前两个区间各不相同，发生变化。5.3回答问题3当有N个应聘人员M个用人单位时，上述方法是可行的。只要我们将人员数，单位数代入（2）即可，于是有N个应聘人员，M个用人单位组成的公务员招聘问题相当于下面的线性规划模型（其中K为要招聘的人员总数）：s.t.（3）六、模型的改进与推广在该模型中，我们对成绩进行了量化。我们也可以考虑对应聘者各方面的能力进行量化来建立模型，如将应聘者在知识面、理解能力、应变能力、表达能力的等级转化成100分制来处理。而在招聘的过程中，招聘领导小组并没有给出笔试成绩及各项能力指标在招聘过程中所占的比重，我们从本模型中的结果发现，成绩的比重不同得到的招聘结果也会不同，因此公务员的招聘应该对各项指标所占的比重应给出必要的衡量标准。本模型对公务员招聘的问题用数学规划的方法进行了处理，该模型还可推广到学生成绩的评定，会议安排，生产计划等等一系列问题。七、参考文献[1]姜启源、谢金星、叶俊编，数学模型（第三版），北京：高等教育出版社，2003年8月。[2]ShoichiroNakamura著，科学计算引论-基于MATLAB的数值分析，北京：电子工业出版社，2002年6月。[3]叶其孝主编，大学生数学建模竞赛辅导教材，湖南：湖南教育出版社，1997年6月。[4]谢兆鸿、范正森、王艮远编著，数学建模技术，北京：中国水利水电出版社，2003年9月。[5]程理民、吴江、张玉林编著，运筹学模型与方法教程，北京：清华大学出版社，2000年1月。[6]周义仓、赫孝良编，数学建模实验，西安：西安交通大学出版社，1999年10月。八、附录8.1附录1笔试成绩的等级化我们把应聘者的笔试成绩看作是应聘的应试能力表现值，把成绩转化为A、B、C、D四个等级的情形，具体方式如下：用应试者中的成绩最高分T与最低分L的差的1/4，即称为分级区间长度。于是介于T与T-d之间的应聘者的成绩转化为A等级；介于T-d与T-2d之间的应聘者的成绩转化为B等级；介于T-2d与T-3d之间的应聘者的成绩转化为C等级；介于T-3d与L之间的应聘者的成绩转化为D等级。参加面试人员的成绩的等级转化结果如下表1所示。表1面试人员成绩等级表人员12345678910111213141516成绩290288288285283283280280280280278277275275274273等级AAABBBCCCCCDDDDD8.2附录2能力和特长等级的量化在招聘过程中，招聘领导小组根据应聘者的笔试成绩，面试成绩和各部门的实际需求来录用和分配人员。由于应聘者的各方面能力和特长以及各部门对应聘者各方面能力的要求分为A、B、C、D四个等级。因此应聘者被分配到某部门工作，必然存在能力等级差距。为了便于计算这些等级差距，我们需要对应聘者的各方面能力等级和各部门对应聘者各方面能力的要求等级进行量化。设，这是由于A与B相差一个等级，因而A-B=4-3=1，A与D相差三个等级，因而A-D=4-1=3。于是应聘者各方面的能力等级量化后的数据如下表2。各部门对应聘者各方面的要求等级量化后的数据如下表3，其中各部门对应聘者的笔试能力的要求在该题中并未给出，我们可以假设为0，即任何的笔试成绩都满足要求。当然这在实际中是不会出现的，因为各用人部门都会给出对笔试成绩不同的要求。表2应聘者各方面的能力值人员12345678910111213141516知识面4434334331143143理解能力4343413433232334应变能力3413342444321423表达能力3223233232444332笔试能力4443332222211111表3各部门对应聘者各方面要求的能力值部门1234567知识面3442222理解能力4332233应变能力2334433表达能力4224444笔试能力00000008.3附录3问题处理的源程序如下：%%%%%%%%%%%%%%%%%%%%下面的程序求应聘者与部门要求的能力的差距矩阵。%%%%%%%%%%%%%%%%%%%%d2004.mclearalllanda1=0.19;landa2=0.81;landa2=landa2/4;%%%%%%%%%%%%%%%%%%%%%%%考试成绩转换程序cj=[290,288,288,285,283,283,280,280,280,280,278,277,275,275,274,273];p=(max(cj)-min(cj))/4;j=0;d=max(cj);fori=1:length(cj)if(cj(i)>=d-p)cj(i)=4-j;elsej=j+1;cj(i)=4-j;d=d-p;endend%%%%%%%%%%%%%%%%%%%%%%%a=[4433;4342;3412;4333;3432;3143;4323;3442;3343;1342;1234;4324;3214;1343;4323;3432];a=[a,cj'];d(1:7)=0;b=[3424;4332;4332;2244;2244;2334;2334];b=[b,d'];%%%%%%%%%%%%%%%%%%%%%%%%%下面的程序求应聘者与部门要求的能力的相对差距矩阵t1(1:7,1:5)=0;fori=1:16forj=1:7t1(j,:)=a(i,:)-b(j,:);t((i-1)*7+j,:)=t1(j,:);endenda=t;%%%%%%%%%tt1(1:112)=0;tt2(1:112)=0;forj=1:112fori=1:5ifi<5ifa(j,i)>=0tt1(j)=tt1(j)+landa2*a(j,i);elsett2(j)=tt2(j)+landa2*abs(a(j,i));endelseifa(j,i)>=0tt1(j)=tt1(j)+landa1*a(j,i);elsett2(j)=tt2(j)+landa1*abs(a(j,i));endendendendd(1:16,1:7)=0;fori=1:16forj=1:7d(i,j)=tt1((i-1)*7+j);endendd1=d;d(1:16,1:7)=0;fori=1:16forj=1:7d(i,j)=tt2((i-1)*7+j);endendd2=d;d=d1-d2%%%%%%%%%%%%%%%%%%%%%%%%%上面的程序求应聘者与部门要求的能力的差距矩阵!!!!!!!!!!!!!!!!!!!!!!!下面是Lingo程序，求解招聘方案的解。!!!!!!!!!!!!!!!!!!!!!!!d2004.lg4MODEL:SETS:w/1..16/:c;v/1..7/:d;links(w,v):a1981,x;ENDSETSDATA:a1981=0.96251.16501.16501.16501.16501.16501.16500.76000.96250.96250.96250.96250.96250.96250.15250.35500.35500.35500.35500.35500.35500.57000.77250.77250.77250.77250.77250.77250.36750.57000.57000.57000.57000.57000.57000.16500.36750.36750.36750.36750.36750.36750.17750.38000.38000.38000.38000.38000.38000.38000.58250.58250.58250.58250.58250.58250.38000.58250.58250.58250.58250.58250.5825-0.2275-0.0250-0.0250-0.0250-0.0250-0.0250-0.0250-0.2275-0.0250-0.0250-0.0250-0.0250-0.0250-0.02500.19000.39250.39250.39250.39250.39250.3925-0.4175-0.2150-0.2150-0.2150-0.2150-0.2150-0.2150-0.2150-0.0125-0.0125-0.0125-0.0125-0.0125-0.0125-0.01250.19000.19000.19000.19000.19000.1900-0.01250.19000.19000.19000.19000.19000.1900;ENDDATA！！！！！！！！！！！！！！！！！！！！！！！！！！！考虑应聘者意愿的！！！！！！！！！！！！！！！！！！！！！！！！！！！也就是设置某些不可能的决策变量为0x(1,1)=0;x(1,6)=0;x(1,7)=0;x(2,2)=0;x(2,3)=0;x(2,6)=0;x(2,7)=0;x(3,4)=0;x(3,5)=0;x(3,6)=0;x(3,7)=0;x(4,1)=0;x(4,2)=0;x(4,3)=0;x(5,1)=0;x(5,6)=0;x(5,7)=0;x(6,1)=0;x(6,2)=0;x(6,3)=0;x(7,2)=0;x(7,3)=0;x(7,4)=0;x(7,5)=0;x(8,1)=0;x(8,4)=0;x(8,5)=0;x(9,2)=0;x(9,3)=0;x(9,6)=0;x(9,7)=0;x(10,2)=0;x(10,3)=0;x(10,6)=0;x(10,7)=0;x(11,2)=0;x(11,3)=0;x(11,4)=0;x(11,5)=0;x(12,1)=0;x(12,2)=0;x(12,3)=0;x(13,4)=0;x(13,5)=0;x(13,6)=0;x(13,7)=0;x(14,2)=0;x(14,3)=0;x(14,6)=0;x(14,7)=0;x(15,2)=0;x(15,3)=0;x(15,4)=0;x(15,5)=0;x(16,2)=0;x(16,3)=0;x(16,4)=0;x(16,5)=0;！！！！！！！！！！！！！！！！！！！！！！！！！！上面的这些决策变量不设置，则相当于不考虑应聘意愿的招聘问题。max=@SUM(links(i,j):a1981(i,j)*x(i,j));！！！！！！！！！！！！！！！！线性规划的目标函数@sum(w(i):@sum(v(j):x(i,j)))=8;!应聘人员总数为8名;@for(w(i):@sum(v(j):x(i,j))<=1);!应聘人员最多去一个部门;@for(v(j):@sum(w(i):x(i,j))<=2);!每个部门最多要两名应聘者;@for(v(j):@sum(w(i):x(i,j))>=1);!每个部门至少要一名应聘者;@for(w(i):@for(v(j):@bin(x(i,j))));END

论大学生写作能力写作能力是对自己所积累的信息进行选择、提取、加工、改造并将之形成为书面文字的能力。积累是写作的基础，积累越厚实，写作就越有基础，文章就能根深叶茂开奇葩。没有积累，胸无点墨，怎么也不会写出作文来的。写作能力是每个大学生必须具备的能力。从目前高校整体情况上看，大学生的写作能力较为欠缺。一、大学生应用文写作能力的定义那么，大学生的写作能力究竟是指什么呢？叶圣陶先生曾经说过，“大学毕业生不一定能写小说诗歌，但是一定要写工作和生活中实用的文章，而且非写得既通顺又扎实不可。”对于大学生的写作能力应包含什么，可能有多种理解，但从叶圣陶先生的谈话中，我认为：大学生写作能力应包括应用写作能力和文学写作能力，而前者是必须的，后者是“不一定”要具备，能具备则更好。众所周知，对于大学生来说，是要写毕业论文的，我认为写作论文的能力可以包含在应用写作能力之中。大学生写作能力的体现，也往往是在撰写毕业论文中集中体现出来的。本科毕业论文无论是对于学生个人还是对于院系和学校来说，都是十分重要的。如何提高本科毕业论文的质量和水平，就成为教育行政部门和高校都很重视的一个重要课题。如何提高大学生的写作能力的问题必须得到社会的广泛关注，并且提出对策去实施解决。二、造成大学生应用文写作困境的原因：（一）大学写作课开设结构不合理。就目前中国多数高校的学科设置来看，除了中文专业会系统开设写作的系列课程外，其他专业的学生都只开设了普及性的《大学语文》课。学生写作能力的提高是一项艰巨复杂的任务，而我们的课程设置仅把这一任务交给了大学语文教师，可大学语文教师既要在有限课时时间内普及相关经典名著知识，又要适度提高学生的鉴赏能力，且要教会学生写作规律并提高写作能力，任务之重实难完成。（二）对实用写作的普遍性不重视。“大学语文”