全国大学生数学建模竞赛A题一等奖论文.docx

1、全国大学生数学建模竞赛A题一等奖论文出版社的资源配置模型摘 要本文讨论出版社的资源优化问题。根据出版社的工作流程，我们将问题分为两个阶段。第一阶段是总出版社如何将总数一定的书号分配给各分出版社；第二阶段是分出版社如何将分得的书号数分配到具体的课程上，以实现利润的最大化。在建立模型确定第一阶段的书号分配方案时，本文侧重于体现长远发展战略和增加强势出版社支持力度的原则，为此我们引入强势的概念，并以此作为目标函数。强势是反映各分出版社的市场占有率、满意度、市场排位等的一个综合指标。我们首先对附件2所给数据提取市场占有率、满意度、市场排位等影响书号数分配的因素，统计出各因素历年的数据，并采用熵权法得到

2、相应的指标权重，然后通过TOPSIS方法得到各分社在总社中的排名强势系数。最后我们将所得到的强势系数带入目标函数，利用Lingo软件计算出各分社应分配的书号数。为了取得更好更贴近实际的结果，我们对模型进行优化，通过引入稳定性的概念来约束分配方案中的奇异现象，最后得到更好的分配方案（表4.6）。 在第二阶段的书号分配过程中，我们以各分社利润最大化为目标又建立了一个优化模型。这里需要解决的难点是预测当年各课程的单位书号的销售量。通过对附件3，4的分析处理，得到各课程往年的单位书号的销售量，并以此为基础运用灰色预测的方法预测出2006年单位书号的销售量。最后用Lingo软件包求解得到结果(表4.8与

3、附录3)。最后我们根据得出的结果，对出版社提出了相应的建议，给出了出版社在分配书号的过程中兼顾短期效益和长远利益时应该考虑的影响因素。关键词：资源优化，熵权法，TOPSIS方法，灰色预测，强势值。1 问题的重述出版社资源配置的好坏直接决定着出版社的经济效益和长远的发展战略，所以如何合理的分配出版社的资源，以达到出版社每年获得的利润最大，而且有利于出版社的长远发展，这就是本题所要解决的问题。出版社最重要的资源就是书号，书号就包括了一个出版社的人力资源、生产资源、资金和管理资源等信息，所以对出版社资源的合理分配就是对出版社的书号进行合理的分配。书号的分配在每个出版社都有一定的程序，以A 出版社为例

4、，假设A出版社主要出版教材类书本，出版社在机构上分为总出版社和分出版社，其中分出版社的划分是根据学科来划分，例如出版计算机类的书为一个分社，出版英语类书本的为另外一个分社，依此类推将A出版社分为9个分社，其关系如图1.1，分社又按课程的不同进行了细分，总社在整个的过程当中只起一个领导规划的作用，对分社的具体资源分配不参与策划。书号的具体分配分为两个步骤，首先就是总社根据各分社提出的书号数申请、人力资源状况和历年的市场信息，在综合考虑当年效益和长远规划的前提下将定量的书号数分给其隶属的9个分社，其中的分配还要遵循以下原则，就是总社要加强对9个分社当中的强势产品的支持力度，优化书号的配置。总社的书

5、号分配完毕之后，各分社再根据各自所分得的书号数按课程进行具体的定量分配，也就是将书号分给每一个课程，其中分配的原则就是要使自身分社在当年获得的利润要最大，分配好之后再安排具体的出版计划进行书本的出版，在分社的具体分配书号的过程当中，总社不参与策划，而且各分社之间的书号分配也是独立的，相互书号的分配没有影响。从出版行业的实际情况出发，通常市场的信息是不完整的，而且各出版社对资料信息的积累和集也是不完善的，也就是说不管从市场角度来看，还是从出版社自身的角度来考虑，信息量都是不全面的，所以这对书号的分配带来了问题，这在实际当中也是一个比较普遍的问题。现在要解决的问题就是在给定一定的市场信息和出版社自

6、身的信息，了解出版社的运做情况下，建立数学模型，将书号进行合理的分配，制定出一个明确的分配方案，使出版社的当年利润最大，对长远的发展有利。2 问题的背景分析随着知识经济社会的到来，知识产业与传媒产业作为新兴产业，正处于飞速发展的进程中。在中国加入WTO之后，出版行业也逐步开始向国内外开放。传统的出版单位也正在进行大规模的转型与改革。巨大的经济效益、良好的投资前景、正在开放的投资机会，使出版发行行业受到投资者的青睐，包括出版发行业在内的传媒产业已经成为最具投资前景的朝阳产业之一。全国图书销售额基本呈现稳步增长趋势。1999-2005年，中国图书纯销售金额从355亿元增长到522亿元左右，年均增长

7、率6.7%。中国图书零售总额从273亿元增长到390亿元左右，年均增长率6.1%。这两个增长率的差异造成库存数量的增加。中国的图书总库存金额从241亿元增长到455亿元，年均增长率高达11%。也就是说，平均每年有30-40亿元的出版新书卖不出去造成库存增加。另一方面，由于库存的增加，最近几年的全球出版业增长趋缓，利润率下降，行业景气程度走低，竞争也是愈演愈烈，所以怎样在这个充满机遇又充满竞争的市场中站稳自己的脚跟，在竞争中处于领先地位是每一个出版社现如今要解决的问题。随着出版行业的对外开放和市场经济的冲击，在市场经济中，市场对资源的配置起基础性作用。根据平均利润率的规律，人们的职业抉择和社会资

8、金会流向获利丰厚、利润率高于社会平均利润率的行业。出版行业虽然走势下降，但是利润率仍然高于社会平均利润。怎样在社会平均利润的基础上，调整好自己内部的各种资源分配关系，以达到最大的利润，并对未来长远的发展有益是出版社的最终要求。3 模型的假设与符号说明3.1模型的假设1.假设问卷调查表的数据客观真实。2.假设一定时期内A出版社的经营正常，不会出现经营不善的情况，而且教材出版行业的市场运行正常。3.假设每年各分社按照总社所分配的实际书号数计算的销售量（附件3）为当年的实际销售量。4.假定同一课程不同书目价格差别不大，同时销售量相近。3.2符号的说明表3.1 符号的说明符 号说 明各个分社的强势总值

9、A出版社书本的利润率 稳定系数每个分社所分得的书号数2006年预测的第i个分社单位书号的销售量第i个分社书本的平均价格第i个分社申请的书号数第i个分社第j个部门的人数 第i个分社第j个部门的工作能力强势系数第j项指标下第i个方案的指标比重指标j的熵值第j项指标的差异系数第j项指标下第i个方案的熵权值第i个分社的利润4 模型的建立与求解4.1 问题的总体分析首先对题目进行分析，可以确定该题是一个优化预测的问题。其中本题给出前5年的相关数据的量相当大，所以数据的统计分析就是一个难点，也是解决本题的一个关键。如何从大量的数据当中获得需要的信息，是本题要做的一个重要工作。本题的总目标是2006年的预计

10、销售利润最大和考虑长远的发展战略。该问题的解决主要分为两个步骤：1）总出版社根据各分出版社提供的信息将总数一定的书号数分配到各个分社，也就是将总社的书号数分为9部分，关键是确立分配的标准与原则；2）分社再根据自己各社的状况对各课程的书号数进行具体的定量分配。 两个步骤的目的就是希望达到利润最大。4.2 步骤1模型的建立与求解总体思路：总社将书号数分配到各个分社，首先就应该以9个分社的总利润为目标，各个分社所分得的书号数为决策变量。利用附件的数据对目标函数进行约束，要求目标所得的利润达到最大，最后用Lingo进行求解。总体的思路比较明了，也由于数据量的庞大，在模型的求解当中会比较麻烦，还需要建立

11、另外的模型对目标函数中的某些量进行预测计算。在分配书号使利润达到最大的情况下，必须要体现的一点就是加强对强势产品的支持力度来优化资源配置。也就要确立强势产品的定义，它与哪些因素有关，定出一个标准来对强势产品进行度量。 4.2.1模型的建立 由上面的分析我们对步骤1首先建立一个单目标的优化模型。在这里我们引入了一个强势的概念。所谓强势是指反映各分出版社的市场占有率、满意度、市场排位等的一个综合指标。目标函数为9个分社的强势值Z，要使得Z的值达到最大，在目标函数当中我们要体现对强势出版社的支持力度来优化资源的配置，在这当中我们引入一个强势系数，以的大小来体现出强势产品，其中，运用熵权法进行确定熵权

12、值。每个分社的强势值为该分社对应的强势系数乘以所分得的书号数。目标函数表达式如下： （1）其中决策变量为每个分社2006年所分得的书号数目。对目标函数的约束条件首先为总的书号数是一定的，为500。 （2）题目要求每个分出版社所分得的书号数不得小于各分社申请总数的一半，即 （3）人力资源对目标的约束，即每个分社每年能完成的书号数有一个上限，上限的确定为各分社三个部门每年能完成书号数的最小值，即 （4）4.2.2模型的求解对于该模型的求解，我们首先要计算出目标函数中的，这样才能运用Lingo软件求出目标函数的最优解。a）附件2的数据处理由于附件2的数据庞大，我们必须对其数据进行压缩、精简，剔除一些

13、和题意无关的数据，提取出我们需要的信息。首先我们要对附件2进行处理，方法如下：1.将附件2中的出版社代号为P115的全部数据提取出来，也就是将A 出版社的全部数据提取出来，得到2001-2005这5年A出版社的全部统计数据。2. Q2f栏为教材的获得方式，由于本题考虑的是从A出版社出版的书本，所以在书的获得方式当中，对于代号3（高年级学生免费赠送或向高年级学生购买二手教材）、代号4（老师复印的教材或讲义）、代号6（其他）这3类对应的数据我们予以剔除。3. Q2k栏为获得的书本是否为新书，由于旧书不能作为当年出版社的一个销售反应量，所以Q2k栏代号为2（旧书）的数据也要剔除。4. 在本题当中，附

14、件3和附件4共统计的课程为72个，也就是A（P115）出版社现在统计的是72类课程书目，在题目给出的统计表中每一个课程名称对应一个编号，但是A出版社所出版的书对应的课程数目却大于72，所以在调查问卷当中，有些课程不在题目所统计的范围之内，也就在调查问卷课程名称一栏当中没有对应的编号，是用中文汉字写出，对于这部分的数据我们也予以剔除。经过上述的处理，我们就将附件2的数据进行了大量的精简，得到了一个数据量相对较小的统计表格，然后再对统计后的表格的每一栏标题作为影响因素，考虑这些因素是否对所建立的目标函数有影响，能确定对书号分配没影响的先进行剔除，其余的再根据计算的数据进行确定。b) 强势系数的确定

15、要确定，首先就要从所给的附件当中提取出与权重相关的指标，经过筛选提取出4类指标，分别为：市场占有率、附件2中的Q1栏（A出版社在你心目中的地位）和Q2l栏（对该书的满意度）。其中对书的满意度有4项，所以一共影响确定的指标数为6个。对于每个指标子我们运用熵权法（参考文献4）进行确定，熵权法是一种决定指标的方法，我们知道，综合指标取决于单个指标数的确定，一般情况下的权重是根据经验来确定的，但是这种确定权重的方法缺少科学根据，也不能保证确立的综合指标能反映原始指标的大部分信息，且权重的确立因人而异，所以其应用受到了限制，而熵权法就能够避免这些问题，使权重的确立具有科学的根据，具有说服力。熵权法的步骤

16、确立如下：1计算第j项指标下第i个方案的指标比重2计算指标j的熵值 （）3计算第j项指标的差异系数 4定义权重则 就为熵权法确定的权重。以下根据上述理论计算每个指标的熵权：1市场占有率的预测从附件2中分别统计出各类学科书本的总量和P115出版社各类学科的书本数，市场占有率方面我们考虑的是A出版社每一学科的书占其对应学科总数的量例如：A出版社计算机类的书市场占有率为A出版社拥有计算机类的书的总数除以市场中计算机类书的总量，各类书的市场占有率计算结果如表4.1表4.1 分社的市场占有率 年份分社20012002200320042005计算机类0.1325150.1411040.1558280.15

17、88670.165025经管类0.2526070.2398610.2838930.2954810.308227数学类0.555120.5691210.5992160.6145660.631761英语类0.0565970.0592310.0591720.0697920.080208两课类0.387260.3389190.3272160.3545720.377207机械、能源类0.5650220.6008970.6576580.6681610.757848化学、化工类0.372340.436170.4315790.4787230.531915地理、地质类0.9065420.8785050.8878

18、50.8785050.906542环境类0.7021280.7340430.7340430.8297870.925532上述表格当中每个分社的市场占有率有5年的数据，在计算的时候我们需要2006年的市场占有率，由于每一年只有5个数据，属于小样本事件，所以我们用灰色模型GM（1，1）对2006年市场占有率进行预测。运用DPS数据处理软件进行预测得出结果如表4.2表4.2 预测每分社市场占有率分社计算机类经管类数学类英语类两课类机械、能源类化学、化工类地理、地质类环境类占有率0.183840.339150.649480.094360.395850.789360.53910.906770.99133

19、对应的条形图如图4.1图4.1 市场占有率条形图2Q1栏（A出版社在你心目中的地位）的预测首先从附件2中统计出Q1栏中5年历史数据每年每一学科的平均地位，然后建立灰色GM（1，1）模型对Q1栏数据进行预测，平均值和预测值见下表4.3表4.3 Q1栏平均与预测值年份分社200120022003200420052006计算机类2.9443.43.4023.383.3283.31939经管类2.98173.60393.27763.31763.29322.9024数学类3.04151.50281.51331.52291.54151.5799英语类2.98772.32752.35882.3932.419

20、92.3564两课类2.11011.95662.13462.11882.11012.19机械、能源类2.85712.2912.28082.24832.23082.2966化学、化工类3.37142.07322.1222.08892.062.0255地理、地质类2.81442.81912.82112.81912.81442.7871环境类2.93943.71013.71013.73083.71263.80553满意度（4个评价）的预测 从附件2中统计出Q2l栏（对该书的满意度）的数据，求每年每一学科的满意度，运用同样方法建立灰色GM（1，1）模型对满意度进行预测，满意度4个方面的评价平均及预测值

21、见附录。根据上面所得的6组预测值，用Matlab按照熵权法的算法步骤进行编程求解熵权值，可得6个指标的熵权值见表4.4表4.4 各指标的熵权值指标出版社地位评价1评价2评价3评价4市场占有率熵权值0.154950.0052210.004330.001430.0019310.83213将预测的2006年的6个指标值构成一个69的矩阵，如下所示 矩阵1其中为第i个指标第j个学科的2006的预测值，然后将矩阵1矩阵转化成矩阵2，转化方法为，得到矩阵2 ，如下矩阵2选出每一列中的最大元素与最小元素，然后根据马氏距离求出第j列中每一个元素到和的距离，公式如下 （5） （6）这里对前面所定义的作如下表达式

22、： （7）公式7中如果越大，那么相应的就会小，的值就会越大，说明第j个分社的强势越明显。最后根据上面计算的结果运用Lingo编写程序（程序见附录）解得书号的最优分配情况如下表4.5表4.5 各分社分配的书号数分社分配的书号数计算机类55经管类47数学类120英语类60两课类36机械、能源类72化学、化工类30地理、地质类40环境类40 由上述模型解得的结果发现大部分分社所分得的书号数要不是申请数的一半就是申请数，也就是大部分的取值要么为上限，要么为下限，从实际的出版情况来说，由于出版行业没有太大的波动，所以各分社所分得的书号数与往年的书号数相比波动也不大，于是在此基础上，我们要考虑对模型的改进

23、。4.3步骤1模型的改进由实际的出版行业可以得知，在出版行业运行正常，没有什么大的改革措施的情况下，各出版社的书号分配也不会出现太大的波动，而且在模型的假设当中，我们也假设了在一段时间内，出版行业的运行是正常的；另一方面，我们可以从历史数据看出每年的书号分配数目也没有太大的波动，总是保持在一个小区间段内，而且在出版社对书号的分配波动性较大的情况下，其市场对应的风险也相应的越大，所以在模型的改进当中我们引入一个系数，定义为稳定系数， （8）其中为第i个分出版社书号分配的平均值，通过引入系数，步骤1模型的目标函数变为 （9）将改进的模型用Lingo编程解得书号的分配情况如表4.6表4.6 各分社分

24、配的书号数分社分配的书号数计算机类70经管类40数学类120英语类95两课类65机械、能源类40化学、化工类21地理、地质类24环境类25 我们将两次求得的结果与前5年的均值进行比较，如图4.2 图4.2 两次预测结果与均值比较（ 代表原始模型分配书号数；*代表改进后模型的分配书号数；折线代表各年各分社分得书号数的平均值）由图中我们可以很直观的看出改进后的模型的数据与平均值的比较靠近，初始模型与均值的相差较大，也就是波动性较强，与实际情况不是很吻合。所以我们认为改进的模型比较成功。4.4 稳定性分析 对总书号数进行变动在400到560中以20为间隔进行稳定性分析，得到表格和图形如下：表 4.7

25、 书号数的变动量与强势值变动量之间的关系序号书号数变动量强势值变动量1400-100173.347-30.8252420-80193.329-10.8433440-60197.684-6.4884460-40203.200-0.9725480-20204.4340.26265000204.1720752020203.622-0.550854040203.250-0.922956060200.980-3.192图 4.3 总书号数的变化对强势值的影响 从表格和图形中可以看出在强势值随书号数的变动以400分界，在400以下总书号数对强势值的影响较大，在400以上对强势值影响较小。总书号在460时强

26、势值达到最大。在420到540之间是比较稳定的，出版社把每年的书号总数定为500的原因可能在这里。4.5步骤2模型的建立与求解总体思路：步骤2是建立在步骤1的结果之下的。在每个分社的书号分配当中，以利润最大作为每个分社的目标函数建立优化模型，在目标函数的建立当中我们要引入单位书号销售量的概念，对进行灰色预测，与步骤1的改进模型类似，在该模型当中我们同样要加入稳定系数来控制书号数的大幅度波动。对于模型的约束条件我们对书号数确立一个下限，下限为历史书号数据的最低值，上限为所申请的书号数，最终运用Lingo对9个目标函数依次进行求解。4.5.1模型的建立每个分社以利润最大为目标函数建立单目标优化模型

27、，表达式如下 （10）为第i个出版社第j个课程分配的书号数，为2006年第i个分社第j个课程的单位书号的销售量，为第i个分社第j 个课程的均价。约束条件的确立：决策变量的上限为第i个出版社第j个课程的申请书号数。下限为第i个出版社所有课程书号数的最低值。这里我们仅以计算机类为例解答（其它类相似）。由前五年的分配的书号数的可以看出最小值为2，所以有如下的约束条件： (j=1,2,310) 4.5.2模型的求解在求解目标函数时，我们首先要确立的值，其确立方法如下：首先将前5年的每个课程的单位书号销售量进行统计，这里以计算机类的数据为例进行说明，单位书号销售量为每年每门课程的实际销售量比上该课程所分

28、配的书号数，这是因为附件2中假设了同一课程不同书目价格差别不大，同时销售量相近，所以我们这样进行平均的处理。在预测2006年的单位书号销售量的时候我们首先是用回归进行拟合，在处理的过程中我们发现拟合的可信度很低，很难找到一个相关量。这个预测是一个小样本的问题，所以我们选择用灰色模型GM（1，1）进行预测，在这个过程中预测的效果都很好（C值基本上都小于0.5，P值都大于0.7），其统计数据和预测数据如表4.7（其它课程的数据见附录2）。表4.7 计算机类单位书号销售量及预测值年份课程2001200220032004200520061124113154.1667240.0909224.333330

29、3.49332180.9155.0909223.4167274.8333327.25403.5948333.3333392.557.5382.5101.3333238.242446069.7562.66667136.6667202.75289.0524548.6666760.3333349.7597.25139.6667188.97476111.498.18182118.75168.8243.5455327.30547280.25296.125316.7778354458.8889526.1727862.6666761.3333395.7539.666677968.983149223.1205235.3636303.8355431.40710148.25162161.1667243361.6479.5463