精 品00数学建模论文——图书馆的馆藏图书分配

2011 精品
2010 高教社杯全国大学生数学建模竞赛
承
诺
书
我们仔细阅读了中国大学生数学建模竞赛的竞赛规则． 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮 件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问 题． 我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他 公开的资料（包括网上查到的资料） ，必须按照规定的参考文献的表述方式在正 文引用处和参考文献中明确列出． 我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性．如有违反 竞赛规则的行为，我们将受到严肃处理．
我们参赛选择的题号是（从 A/B/C/D 中选择一项填写） ： 我们的参赛报名号为（如果赛区设置报名号的话） ： 所属学校（请填写完整的全名） ： 参赛队员 (打印并签名) ：1． 2． 3． 指导教师或指导教师组负责人 宁 波 工 程 学 院 李瑜苗 杨路捷 吴建明 (打印并签名)： 数 模 组
A
日期：
2010 年 8 月 7 日
赛区评阅编号（由赛区组委会评阅前进行编号） ：
ksdowe
2011 精品
2010 高教社杯全国大学生数学建模竞赛 高教社杯全国大学生数学建模竞赛
编 号 专 用 页
赛区评阅编号（由赛区组委会评阅前进行编号）：
赛区评阅记录（可供赛区评阅时使用）： 评 阅 人 评 分 备 注
全国统一编号（由赛区组委会送交全国前编号）：
全国评阅编号（由全国组委会评阅前进行编号）：
ksdowe
2011 精品
图书馆的馆藏图书分配摘 要
本文根据题意将数据进行分析处理，建立三个模型得到各类图书的综合权 重，以综合权重来确定应有的册数比例，给图书馆提供合理的图书采购方案． 对于问题一， 利用已知的信息， 采用 AHP 层次分析法[1]和模糊综合评价法[2]， 综合考虑各类图书的重点实验室和重点学科建设的需求与否、常用和热门程度、 重要公共课与否、技能课图书与否、上年的出借册次分布五个指标，确定各指标 对各类图书隶属度关系和模糊矩阵，从而得到相对重要程度的权重． 即各类书的相对重要程度的权向量分别为： 0.0484 0.0508 0.0367 0.0418 0.0405 0.0564 0.0368 0.0621 0.0554 0.0418 0.0375 0.0472 0.0357 0.0556 0.0495 0.0362 0.0414 0.0448 0.0673 0.0521 0.0357 对于问题二，采用和问题一类似的方法，对确定了的一年内的相对流通量、 图书平均借用时间、图书利用率三个指标进行处理，从而得到书籍在该校的实际 使用价值的相应权重．但值得注意的是，在问题二的处理过程中，采用了熵值法 对三个指标赋权，减少了主观因素对反映实际使用价值过程的影响． 即各类书的实际使用价值的权向量

分别为： 0.0415 0.0533 0.0402 0.0432 0.0506 0.0542 0.0389 0.0582 0.0497 0.0468 0.0421 0.0581 0.0436 0.0564 0.0549 0.0354 0.0473 0.0447 0.0585 0.0535 0.0288 由上数据可知对应的权值越大相应的价值越大． 对于问题三，引进了一种“席位分配[4]”的数学模型，在满足各类图书的最 低更新率的基础上，使得结果的相对不公平指标最小，从而得到了最优化的购书 分配方法． A、 类书增加 270 册， N 类书增加 520 册,D 类书增加 164 册,E、 即 B C、 R、S、U、V 类书增加 218 册,X 类书增加 1773 册,F 类书增加 1058 册,G 类书增 加 191 册,H3 类书增加 599 册,I 类书增加 898 册,J 类书增加 213 册,K 类书增加 234 册,Q 类书增加 247 册,P 类书增加 1987 册,O 类书增加 966 册,TH 类书增加 1068 册,TM、TS 类书增加 109 册,TN 类书增加 143 册,TQ 类书增加 483 册,TP 类书增 加 693 册,TU、TV 类书增加 315 册,Z 类书增加 965 册． 对于问题四，通过所编写的 MATLAB 的程序，不善于数学建模的人只需要 输入上一年该普通高校图书馆的馆藏图书的分布及流通情况和下一年计划投入 的总资金额这两个数据，就可以得到图书馆下一年的购书资金分配方案．
关键词： 关键词：AHP 层次分析法
熵值法
图书分配
席位分配
1
2011 精品
一、问题的提出现代化图书馆馆藏图书，主要目的不是为了收藏而是为了使用．除了国家图 书馆等特大型的图书馆以外，一般图书馆都有特定的服务群体，办馆宗旨就是要 尽量好地为这些特定群体服务，提高馆藏资源的利用率、读者文献信息需求的满 足率以及对图书馆服务功能的满意率．图书馆每年用于购书的经费是有限的，如 何合理分配使用， 以便使有限的购书经费最大限度地发挥其特定的经济效益是图 书馆工作的重要环节之一． 以某学校图书馆为例，要实现办馆效益，必须做到入藏文献合乎本校教师、 学生（有时也兼顾社会）的需求，使图书馆藏书结构（学科结构、文种结构、文 献类型结构等）能满足本校教学科研的要求，以求藏书体系与本校专业设置相适 应．所购图书要能够真实地反映读者的实际需要，使读者结构和藏书结构尽量吻 合，以便减少读者借不到图书的现象，即降低读者拒借的比率、增加满足率．文 献只有在流通中才能传播信息，产生效益．文献资料得不到利用，购置文献资料 所耗费的资金就体现不出其价值．因此，图书馆在增加藏书规模的同时，要千方 百计地把文献提供给读者，以增加图书的出借次数、出借时间以及在借图书的数 量等，力求使有限的价值投入获得最大的办馆效益． 该校图书馆每学年都要投入大量资金购置图书，图书覆盖全院各学科专业、 具有较完整的

中外文文献资源． 如何合理分配资金用于各种图书的购置成为一个 非常有价值的问题．
二、基本假设1、假设该校各专业学生比例大致不会变且总人数相对稳定； 2、假设所借的图书没有不归还或丢失； 3、假设每年借书人数相对稳定； 4、假设题目给的数据真实有效．
三、定义符号说明符号 ω 、ω ' ： 说明 指标权重向量 各个权向量对应的权重 模糊关系矩阵 最大特征根 一致性指标 随机一致性指标 模糊综合评价结果的向量 标准矩阵 各指标同度量化 熵值
ωi ：R：
λmax ：CI： RI： B 、B'： Y： p： ωe ：
2
2011 精品
g： R'： ω '' ：
差异性系数 各指标对各类图书模糊矩阵 三项指标的权向量 每个综合评价值所对应的图书册数 剩余的资金 按每个综合评价值分配得到的图书对应的资金 综合评价值
ci ei M
pi
e
四、模型的分析、模型的建立及求解4．1 问题一的分析及模型建立： ． 问题一的分析及模型建立： 4．1．1 问题一的提出和分析 ． ． 已知： 已知：①该普通高校的重点学科、重点专业的设置情况；②上一年该普通高 校图书馆的馆藏图书的分布及流通情况表． 要求： 要求：同时考虑重点实验室和重点学科建设的需求、常用书籍和流行热门书 籍、重要公共课、技能课图书（如英语、计算机类）的普通要求等几个方面，以 确定各类图书对于该校图书馆的相对重点程度（即相对权重） ． 分析： 我们以各类图书的重点实验室和重点学科建设的需求 分析：根据题目要求， 与否、常用和热门程度、重要公共课与否、技能课图书与否、上年的出借册次分 布这五个指标来衡量各类图书对于该校图书馆的相对重点程度（即相对权重） ． [1] 首先，通过层次分析法 确定个指标之间的权重，利用方根法计算出反映各 指标相对权重的权向量，并进行一致性检验；然后，确定各类图书与各指标之间 的隶属度关系，最后，利用加权平均型合成算子确定模糊综合评价[2]结果，得到 各类图书对于该学校图书馆的相对重要程度（即相对权重） ． 4.1.2 模型一的建立及求解 模型一 首先，用层次分析法确定权重,判断矩阵由 A．L．Saaty 的 1～9 比率标度方 法确定,结果见下关系表如下表 1： 重点学科 常用热门 公共课 技能课 出借册数 重点 1 4/3 2 4 5 常用 3/4 1 2 3 4 公共 1/2 1/2 1 2/1 3/1 技能 1/4 1/3 1/2 1 3/2 出借册数 1/5 1/4 1/3 2/3 1 表一 由表可得
3
2011 精品
? 1 ? ?3/ 4 R = ? 1/ 2 ? ? 1/ 4 ? 1/ 5 ?
4/3 2 4 5 ? ? 1 2 3 4 ? 1/ 2 1 2 3 ? ? 1/ 3 1/ 2 1 3 / 2 ? 1/ 4 1/ 3 2 / 3 1 ? ?
用方根法计算得出权向量： ω = {ω1 , ω2 , ω3 , ω4 , ω5 } = ( 0.3647 0.2941 0.1795 0.0946 0.0672 )CI =

进 行 一 致 性 检 验 ： 最 大 特 征 根 λ max = 5.0166 ， 一 致 性 指 标= 0.0042 ， 平 均 随 机 一 致 性 指 标 RI = 1.12 ， 随 机 一 致 性 比 率 n ?1 CI CR= = 0.00375 < 0.1 ， 因此该判断矩阵具有满意的一致性． 故通过一致性检验． RI 然后，分别给 21 类书在 5 个指标上进行打分 重点学 图书类 重要公 技能课 出借册 常用 科和重 别 共课 图书 数 点专业 A、B 3 3 6 3 7433 C 、N 3 4 3 3 20209 D 3 2 3 3 3773 E、 、 、 R S U、V 3 3 3 3 5900 X 3 3 3 3 1725 F 3 4 3 3 38194 G 3 2 3 3 4007 H3 3 4 6 6 23893 I 3 4 3 3 34937 J 3 3 3 3 5681 K 3 2 3 3 6164 Q 3 4 3 3 8670 P 3 2 3 3 462 O 3 4 6 3 2987 TH 5 4 3 3 2997 TM、TS 3 2 3 3 2236 TN 3 3 3 3 4413 TQ 5 3 3 3 2274 TP 5 4 6 6 27314 TU、TV 5 4 3 3 11435 Z 3 2 3 3 558 表二
λ max ? n
利用 MATLAB 软件对其进行归一化处理：4
2011 精品
图书类 别 A、B C、N D E、 S、 R、 U、V X F G H3 I J K Q P O TH TM、TS TN TQ TP TU、TV Z
重点 0.044944 0.044944 0.044944 0.044944 0.044944 0.044944 0.044944 0.044944 0.044944 0.044944 0.044944 0.044944 0.044944 0.05618 0.05618 0.044944 0.044944 0.05618 0.05618 0.05618 0.044944
常用 0.045455 0.060606 0.030303 0.045455 0.045455 0.060606 0.030303 0.060606 0.060606 0.045455 0.030303 0.060606 0.030303 0.060606 0.060606 0.030303 0.045455 0.045455 0.060606 0.060606 0.030303
公共课 0.068182 0.034091 0.034091 0.034091 0.034091 0.034091 0.034091 0.068182 0.034091 0.034091 0.034091 0.034091 0.034091 0.068182 0.034091 0.034091 0.034091 0.034091 0.068182 0.034091 0.034091 表三
技能书 0.043478 0.043478 0.043478 0.043478 0.043478 0.043478 0.043478 0.086957 0.043478 0.043478 0.043478 0.043478 0.043478 0.043478 0.043478 0.043478 0.043478 0.043478 0.086957 0.043478 0.043478
出借册 数 0.03453 0.093881 0.017527 0.027408 0.008013 0.17743 0.018615 0.110995 0.1623 0.026391 0.028635 0.040277 0.002146 0.013876 0.013923 0.010387 0.020501 0.010564 0.126887 0.053121 0.002592
确定各指标对各类图书的隶属度，得到模糊矩阵 R ? 0.042253521 0.042253521 L 0.070422535 0.042253521? ? ? 0.060606 L 0.060606 0.030303 ? ? 0.045455 R = ? 0.068182 0.034091 L 0.034091 0.034091 ? ? ? ? 0.046512 0.046512 L 0.046512 0.046512 ? ? 0.096422 L 0.038864 0.0206 ? ? ? 0.031669 利用加权平均型合成算子，将权向量 ω 与模糊关系矩阵 R 合成，得到模糊 综合评价结果向量 B ： B = ω ? R = (0.0484 0.0508 L 0.0521 0.0357)
5
2011 精品
各类书的重要性权重 0.07 0.06 0.05 0.04 0.03 0.02 0.01 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 书类 权重 系列1
即各类书的相对重要程度的权向量分别为： 0.0484 0.0508 0.0367 0.0418 0.0405 0.0564 0.0368 0.0621 0.0554 0.0418 0.0375 0.0472 0.0357 0.0556 0.0495 0.0362 0.0414 0.0448 0.0673 0.0521 0.0357 4．2 问题二的分析及模型建立： ． 问题二的分析及模型建立： 4．2．1 问题二的提出及分

析 ． ． 已知： 已知：①图书最终的实现价值应取决于图书的被利用率；②上一年该普通高 校图书馆的馆藏图书的分布及流通情况表． 要求： 要求： 评价一本书的真正价值必须考虑到它的流通量大小和借用时间的长短 等，并根据该校上一年各类图书的出借情况，提出一种评价图书在该校实际使用 价值的方法． 分析： 我们同样先给出用来评价图书在该校的实际使用价值 分析：与问题一类似， 的三个指标：一年内的相对流通量、图书平均借用时间、图书利用率．然后，我 们采用熵值法 [3] 对这三个指标进行赋权处理，从而得到指标的权向量；利用 MATLAB 软件对各类图书的三个评价指标数据表中的数据进行归一化，从而得 到各指标对各类图书的隶属度模糊矩阵．最后，利用加权平均型合成算子确定模 糊综合评价结果，得到各类图书在该学校图书馆的实际使用价值． 4．2．2 模型二的建立及求解 ． ． 图书最终的实现价值应取决于图书的被利用率． 因而评价一本书的真正价值 必须考虑到它的流通量大小和借用时间的长短等多方面的指标． 在模型二的评价中，设定了三个指标来评价图书在该校的实际使用价值，三 个指标如下：一年内的相对流通量（出借册数/册数） 、图书平均借用时间（出借 总时间/出借册数） 、图书利用率（出借种类数/内容种类数） （数据见下表） ． 一年内 图书平 图书利 的相对 均借用 用率 （出 流通量 时间 （出 借种类/ （出借 借总时 内容种 册数/册 间/出借 类） 数） 册数） 0.8267 33.998 0.6199536
图书类 别
A、B
2011 精品
C、N D E、 S、 R、 U、V X F G H3 I J K Q P O TH TM、TS TN TQ TP TU、TV Z
1.1667 0.6884 0.812 0.9965 1.0833 0.6285 1.1957 1.1667 0.8 0.7909 1.0547 0.621 1.0481 1.1039 0.6134 0.9257 0.8133 1.1818 1.0887 0.4143
45 0.700035 38.9973 0.590038 42 50.9948 46 41.0015 51.0004 35 47.9982 40.0023 54.997 51.0368 54.0027 51.006 32.0009 40.0029 44.0018 50.0003 45.9989 29.9821 0.580179 0.619741 0.780047 0.549846 0.790039 0.700027 0.649945 0.580017 0.827949 0.597744 0.780022 0.720442 0.549929 0.69961 0.610515 0.828165 0.75 0.44898
采用熵值法对三个指标进行赋权处理： 1．对原始数据进行标准化处理，得到标准矩阵 Y = ( y ij ) 21×3 计算公式为 yij = 值． 得到 Y 矩阵（见附录 1） 2．将各指标同度量化，计算第 j 项指标下第 i 类书指标值的比重 y ij pij = 21 (i ≤ i ≤ 21,1 ≤ j ≤ 3) ∑ yiji =1
xij Mj
(i = 1,2,L ,21; j = 1,2,3) 其中 M j 为第 j 个指标的最大
得到 p 矩阵（见附录 2） 3．计算第 j 项指标的熵值
ω e j = ?k ∑ pij ln piji =1
21
(1 ≤ j ≤ n)
其中 k =
1 ln 21
?1 21 ∑ pij ln pij (1 ≤ j ≤ 3) ． ln 21 i =1 得到 ω e 矩阵（见附录 3） 4．计

算第 j 项指标的差异性系数 g j = 1 ? ω e j (1 ≤ j ≤ 3) ，其值越大，指标就
则 ωe j =
越重要；7
2011 精品
5．确定指标权重，第 j 项指标的权数 ω j ' =
gj
∑gj =1
21
(1 ≤ j ≤ 3) ．j
得到 ω ' 的一个指标权重向量 ω ' = (0.3297 0.3322 0.3382) ． 利用 MATLAB 软件对模型二中各类图书的三个评价指标数据表中的数据进 行归一化处理得到： 图书类 别图书利用 一年内的 借用时间 率 （出借种 相对流通 （出借总 类/内容 量 （出借册 时间/出 数/册数） 种类） 借册数） 0.0434641 0.036754 0.044367 0.0613397 0.048647 0.050098 0.0361929 0.042158 0.042226 图书平均
A、B C、N D
0.045404 0.041521 X 0.0523914 0.055128 0.044352 F 0.0569549 0.049729 0.055824 G 0.0330436 0.044325 0.03935 H3 0.0628644 0.055134 0.056539 I 0.0613397 0.037837 0.050098 J 0.0420603 0.051889 0.046514 K 0.0415819 0.043245 0.041509 Q 0.0554513 0.059455 0.059253 P 0.0326493 0.055174 0.042778 O 0.0551043 0.05838 0.055823 TH 0.058038 0.05514 0.051559 TM、TS 0.0322498 0.034595 0.039356 TN 0.0486691 0.043245 0.050068 TQ 0.0427596 0.047568 0.043692 TP 0.0621336 0.054053 0.059268 TU、TV 0.0572388 0.049727 0.053674 Z 0.021782 0.032412 0.032131 确定各指标对各类图书的隶属度，得到模糊矩阵 R ' ? 0.0434641 0.0613397 L 0.0572388 0.021782 ? ? ? R ' = ? 0.036754 0.048647 L 0.049727 0.032412 ? ? 0.044367 0.050098 L 0.053674 0.032131 ? ? ? 利用加权平均型合成算子，将权向量 ω ' 与模糊关系矩阵 R ' 合成，得到模糊 综合评价结果向量 B ' ： B ' = ω ' R ' = (0.0415 0.0533 L 0.0535 0.0288)
E、R、S、 U、V 0.0426912
8
2011 精品
各类书实用价值权重 0.07 0.06 0.05 权重 0.04 0.03 0.02 0.01 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 书类 系列1
即各类书的实际使用价值的权向量分别为： 0.0415 0.0533 0.0402 0.0432 0.0506 0.0542 0.0389 0.0582 0.0497 0.0468 0.0421 0.0581 0.0436 0.0564 0.0549 0.0354 0.0473 0.0447 0.0585 0.0535 0.0288 由上数据可知对应的权值越大相应的价值越大． 4．3 问题三的分析及模型建立 ． 4．3．1 问题三的提出及分析 ． ． 已知： 已知：①通过前两问研究，我们得到了各类图书的相对重要程度和其在该校 的实际使用价值所对应的权重； ②上一年该普通高校图书馆的馆藏图书的分布及 流通情况表；③图书馆计划投入 100 万元用于购置各种图书． 要求： 在所确定的购书资金分配方案应尽可能符合学校学科发展的需要和教 要求： 学科研需要，又应当尽可能提高读者的满意率，与此同时，图书馆自然还应当注 意到各类馆藏图书的更新率， 以及用于购书的总经费是有限制的这四个条件的约 束下，尽可能满足目标函数（所购图书实际效益最大） ． 分析： 分析： 在前两问的研究基础之上，我们引进两

个新概念———读者满意度 （（册数/内容种类数）/ 出借册数/出借种类数））和实际效益的综合评价值．但 （ 由于读者满意度越大，图书所需更新的比率应该相对越小，所以用读者满意度的 倒数来作为衡量综合评价值的第三指标．于是，就有了衡量综合评价值的三项指 标，即图书相对重要权重、图书实际使用价值权重、读者满意度的倒数．然后， 利用熵值法来可以确定综合评价实际效益中三个指标之间的相对权重， 继而得到 各类图书综合评价值． 在制定购书资金的分配方案时，我们又引进一个称为“席位分配”的数学模 型．这个模型从每个席位对应的人数出发，定义了一个相对不公平指标，将有限 的代表席位逐个分配到各个小组，结果使得相对不公平指标最小．而不公平程度 c 以各类图书的 i （每个综合评价值所对应的图书册数）的方差值来衡量，方差 ei 值越小，不公平程度也越小．而在图书更新过程中，还应该注意到各类书的更新 率．通过资料查找，我们发现图书馆的各类图书的更新率至少应达到 3%．所以， 采用在满足最低更新率的基础上，进行“席位分配”的方法，利用 MATLAB 软 件所编程序，最终得到最优的购书资金分配方案． 4．3．2 模型三的建立及求解 ． ．9
2011 精品
应用模型二提到的熵值法，确定该三项指标在综合评价中的权重．（附录） 得到三项指标的权向量： ω '' = ( 0.3314 0.3323 0.3363) 再运用模糊综合评价的方法，求得 21 类书的综合评价值 e 为： 0.0461 0.0548 0.0397 0.0452 0.0497 0.0537 0.0391 0.0585 0.0552 0.0444 0.0430 0.0505 0.0390 0.0537 0.0533 0.0374 0.0456 0.0460 0.0593 0.0528 0.0327各类书的综合权重 0.07 0.06 0.05 权重 0.04 0.03 0.02 0.01 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 书类 系列1
在本模型中，以该类书的综合权重反映书的应有册数比重． 要考虑 21 类书的购书比例， 假设有 A、 两类书的综合评价值为 e1 、 2 ， B e c c c c 该类书原有的册数为 c1 、 c 2 ，模型希望达到的是 1 = 2 ，若 1 < 2 ，则认 e1 e2 e1 e2 为 A 书是亏欠的，应补．推广到 21 类书，为了保证每类书都有一定更新率， 在满足各类图书最低分配率的基础上，将剩余的经费分配给其余亏欠的图 书，直到达到资金上限（100 万） ．具体步骤如下：1．先满足每类书 3%的更新率，判断剩余的资金 M 是否大于零，若 M > 0 ，则转步骤 2；若 M <= 0 ，则转步骤 4． c 2．比较所有的 i (i = 1,2,3L 21) ，值最小的那类书，所对应的 ci + 1 ， ei 即购一本书，然后 ci = ci + 1 ， M = M ? p i ，对更新的比值计算方差．3．再判断 M ，若 M > 0 ，则转步骤 2；若 M <= 0 ，则转步骤 4． 4．各类书的购书过

程结束，输出每类书所购的册数，及每循环一次所 得的方差．选择方差最小的所对应的购书比值（即最符合各类书应有的册数 比） ，剩余的资金再按照册数符合综合权重的原则分配． 5．输出最终各类书的所购的册数以及所花费的资金． 利用 MATLAB 软件的程序运行得到最优的分配方案为： 类别 A 、B C 、N D E、 、 、 X R S F G U 、V 520 164 218 1773 1058 191 新增册 270 数 类别 H3 I J K Q P O10
2011 精品
新增册 数 类别 新增册 数
599 TH 1068
898
213
234 TQ 483
247 TP 693
1987
966
TM、TS TN 109 143
TU、TV Z 315 965
共花费总资金 100 万元． 4．4 问题四的分析及求解 ． 问题四是前三问的分析综合得到的权重以确定图书的采购方案， 根据题意只 要简要的阐述馆方输入哪些数字，怎样操作获得合理的购书方案即可． 我们决定编一个小系统，馆方只需两个数据的 txt 文件（其中第一个 txt 数 据为： 重点学科和重点专业 常用 重要公共课 技能课图书 出借册书； 第二个 txt 数据为：出借册数 册数 出借时总时间 出借种类 内容种类）和一个计划投入总 资金，利用 MATLAB 软件读取上面的数据就可以得到比较合适的购书方 案．（matlab 程序见附录问题四）
五、结果分析由综合评价值 e 矩阵可知“自动化技术、计算机技术综合权重”最大，“常 用外国语”“文学”“ 经济”“数理科学和化学” “机械仪表工业” “建筑 科学、水利工程”也较大，将模型出来的结果与题目给出的数据进行比较，可知 权重大的重点类建设对象或是重点学科．而“综合性图书”权重最小，与题目给 的数据也较符合． 根据综合权重的柱状图与原册数的实际比例作对比各类书的综合权重 0.07 0.06 0.05
图书馆现有书册数 40000 35000 30000 25000 20000 15000 10000 5000 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 书类
册数
权重
0.04 0.03 0.02 0.01 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 书类 系列1
系列1
综合权重大的原有册数也多，但有存在例外．并且模型得出的权重大小波动 较小，而原有册数差距较大． 分析原因：1：图书馆现有图书分布并不合理； 2：模型中设置的权重并不很合理．
六、模型推广本题不仅可以用于优化图书的采购方案问题， 也可以用于各类相似的评价问 题中．模型三用到的席位分配模型，可以用以解决生活、工作中可能产生的资源 分配公平与否的问题． 与此同时，本文涉及的模型可以用以评价多种属性（指标）的对象．
11
2011 精品
七、模型的评价与改进优点： 1、利用多种方法确定权重，一定程度上减少了主观因素对结果的影响． 2、基于层次分析、模糊综合评价模型并进行改进

，结果符合实际． 3、模型三的算法逻辑清晰、易懂，运用软件，减少大量计算量． 缺点： 1、模型一中，运用层次分析法确定五个指标的权重，带有主观因素，一定程度 上影响结果． 2、多次运用模糊综合评价，算法较单一．
参考文献：[1]吴祈宗，运筹学与最优化方法，北京：机械工业出版社，2003 年； [2]张秀兰，基于模糊综合评判法的研究及应用，科技信息 2008 年 14 期：91— 92，2008 年； [3]沈红丽，因子分析法和熵值法在高校科技创新评价中的应用，河北工业大学 学报 第 38 卷 第 1 期，2009 年 2 月； [4]靖培栋 刘忠厚，图书馆外文核心期刊购买模型探讨，中国图书馆学报（双月 刊）1999 年第 4 期第 25 卷 44-48，1999 年； [5]臧秀平 李萍 张建 丁声铎，改进的模糊综合评判法在评标中的应用，江苏科 技大学学报（自然科学版）第 21 卷 第 6 期，2007 年 12 月．
12
2011 精品
附件代码： 问题一 MATLAB 代码： %Maxlmta.m %和法求最大特征根 clc clear all disp('please choose the filename you want to plot'); %查找数据的文件夹 [filename,pathname]= uigetfile(' *.txt', 'choose the file you want to plot'); if pathname==0 %pathname 返回 0 说明文件打开失败，可能是取消了，或是文件不存在等 等原因 return %return 用于退出整个程序 end name =[pathname filename]; 路径和名字 fid=fopen(name,'r+'); %读取文件 x=fscanf(fid,'%c'); %count 得到数据个数，A 是列向量，用于存放所有数据 A=str2num(x); %下面的 A 是一个测试的程序矩阵 %A=[1 1/2 4 3 3; % 2 1 7 5 5; % 1/4 1/7 1 1/2 1/3; % 1/3 1/5 2 1 1; % 1/3 1/5 3 1 1]; %RI--随机一致性指标 %n--A 的列长度 %w--权向量 %lmta--最大特征根 %CI--一致性指标 %CR--一致性比率13
%文件的
2011 精品
%RIn--A 的一致性指标 %flag--标志变量 %Wij Wi W--临时变量 RI=[0,0,0.58,0.90,1.12,1.24,1.32,1.41,1.45,1.49,1.51]; %将 A 的每一列向量归一化 Asum=sum(A); n=length(A); for j=1:n for i=1:n Wij(i,j)=A(i,j)./Asum(1,j); end end %将 Wij 按行求和 Wij=Wij'; Wi=sum(Wij); Wi=Wi'; %将 Wi 归一化 W=sum(Wi); disp('权向量'); w=Wi./W %计算 lmta disp('最大特征根\n'); lmta=sum(1/n*(A*w)./w) %计算 CR disp('A 的一致性指标'); RIn=RI(1,n); disp('一致性指标'); CI=(lmta-n)/(n-1) CR=CI/RIn; %判断一致性检验 if CR<0.1|n<=2 disp('通过一致性检验') flag=1;14
2011 精品
else disp('不能通过一致性检验') flag=0; end 代码： 问题二 MATLAB 代码： %shuangzhifa.m %层次分析法中的熵值法 function t clear all clc %输入数据 [filename pathname]=uigetfile('*.txt','please choose the file'); name=[pathname filename]; if filename==0 return end fid=fopen (name,'r+'); x=fscanf (fid,'%c'); x=str2num(x); %对 x 标准化得到 y xsize=size(x); temp=max(x); for i=1:xsize(1) y(i,:)=x(i,:)./temp; end clear temp; %计算 p(i,

j) temp=y'; temp=sum(temp); for i=1:xsize(1) for j=1:xsize(2) p(i,j)=y(i,j)/temp(j); end end clear temp; %计算 e(j)15
2011 精品
k=1/log(xsize(1)); for j=1:xsize(2) tempsum=0; for i=1:xsize(1) temp=p(i,j)*log(p(i,j)); tempsum=tempsum+temp; end e(j)=-k*tempsum; end %求差异性系数 g(j) g=1-e; clear temp; temp=sum(g); for j=1:xsize(2) w(j)=g(j)./temp; end w return 问题三 MATLAB 代码： 代码： %shijixiaoyiguihua.m function temp clear all clc c=[8991 17322 5481 7266 1731 35256 6375 19983 29946 7101 7794 8220 744 2850 2715 3645 4767 2796 23112 10503 1347]; e=[0.0461 0.0548 0.0397 0.0452 0.0497 0.0537 0.0391 0.0585 0.0552 0.0444 0.043 0.0505 0.039 0.0537 0.0533 0.0374 0.0456 0.046 0.0593 0.0528 0.0327]; price=[225674.1 486748.2 121130.1 189916 65951.1 1131718 128137.5 561522.3 7339765 199538.1 304745.4 206322 20162.4 100035 98011.5 102424.5 185913 64587.6 834343.2 305637.3 60749.7]; tmp=c./e; ptmp=tmp; pc=c; eachprice=price./c; x=zeros(1,21);16
2011 精品
ppx=round(c*0.03); x=ppx+x; pprice=ppx.*eachprice; px=x; STD=std(ptmp); while sum(eachprice.*px)<=(1000000-sum(pprice)) [pnumber pindex]=min(ptmp); px(pindex)=px(pindex)+1; pc(pindex)=pc(pindex)+1; ptmp=pc./e; ptemp=std(ptmp); STD=[STD ptemp]; end n=size(STD,2); STD=STD(1,[1:n-1]); [number index]=min(STD); for i=1:index [pnumber pindex]=min(tmp); x(pindex)=x(pindex)+1; c(pindex)=c(pindex)+1; tmp=c./e; end x sumprice=sum(x.*eachprice)+sum(pprice) return 代码： 问题四 MATLAB 代码： %wenti4.m function temp clear all clc disp('读入一个 txt 的矩阵数据文件\n');
17
2011 精品
disp('数据的顺序为： 重点学科和重点专业 常用 重要公共课 技能课图书 出借 册书\n'); [filename pathname]=uigetfile('*.txt','请选择数据'); name=[pathname filename]; dif=fopen(name,'r'); x1=fscanf (dif,'%f'); n=size(x,1); w=[0.3647 0.2941 0.1795 0.0946 0.0672];%人为定义的一个奴隶度关系得到的 权向量 xsum=sum(x1); for i=1:n x1(i,:)=x1(i,:)./xsum; end x1=x1'; w1=w*x1; %得到第一个 w1 disp('请输入第二个数的 txt 文件的矩阵数据\n'); [filename pathname]=uigetfile('*.txt','输入的数据格式为：出借册数 册数 出借时总时间 出借种类 内容种类\n'); name=[pathname filename]; if filename==0 return end fid=fopen (name,'r+'); data=fscanf (fid,'%c'); data=str2num(data); x2=[data(1)/data(2) data(3)/data(1) data(4)/data(5)]; %对 x2 标准化得到 y xsize=size(x2); temp=max(x2); for i=1:xsize(1) y(i,:)=x2(i,:)./temp; end clear temp; %计算 p(i,j) temp=y'; temp=sum(temp); for i=1:xsize(1)18
2011 精品
for j=1:xsize(2) p(i,j)=y(i,j)/temp(j); end end clear temp; %计算 e(j) k=1/log(xsize(1)); for j=1:xsize(2) tempsum=0; for i=1:xsize(1) temp=p(i,j)*log(p(i,j)); tempsum=tempsum+temp; end e(j)=-k*tempsum; end %求差异性系数 g(j) g=1-e; clear temp; temp=sum(g); for j=1:xsize(2) w(j)=g(j)./temp; end xxsum=sum(x2); for i=1:xsize(1) x2=x2./xxsum; end x2=x2'; w2=w*x2; %得到第二 w2 数据 w3=(data(1)/data(4))/

(data(2)/data(5)); pw=[w1 w2 w3]; ppw=pw; ppwsum=sum(ppw); n=size(ppw,1); for i=1:n ppw(i,:)=pw(i,:)./ppwsum;19
2011 精品
end R=ppw'; %对 pw 标准化得到 y pwsize=size(pw); temp=max(pw); for i=1:pwsize(1) y(i,:)=pw(i,:)./temp; end clear temp; %计算 p(i,j) temp=y'; temp=sum(temp); for i=1:pwsize(1) for j=1:pwsize(2) p(i,j)=y(i,j)/temp(j); end end clear temp; %计算 e(j) k=1/log(pwsize(1)); for j=1:pwsize(2) tempsum=0; for i=1:pwsize(1) temp=p(i,j)*log(p(i,j)); tempsum=tempsum+temp; end e(j)=-k*tempsum; end %求差异性系数 g(j) g=1-e; clear temp; temp=sum(g); for j=1:pwsize(2) ppw(j)=g(j)./temp;20
2011 精品
end pxxsum=sum(ppw); for i=1:xsize(1) ppw=ppw./pxxsum; end ppw=ppw'; w=ppw*pw; %得到 w 数据 e=w*R; c=data(1); disp('请输入一个 txt 的图书总价数据的矩阵'); [filename pathname]=uigetfile('*.txt','请择择文件'); name=[pathname filename]; dif=fopen(name,'r'); price=fscanf (dif,'%f'); %c=[8991 17322 5481 7266 1731 35256 6375 19983 29946 7101 7794 8220 744 2850 2715 3645 4767 2796 23112 10503 1347]; %e=[0.0461 0.0548 0.0397 0.0452 0.0497 0.0537 0.0391 0.0585 0.0552 0.0444 0.043 0.0505 0.039 0.0537 0.0533 0.0374 0.0456 0.046 0.0593 0.0528 0.0327]; %price=[225674.1 486748.2 121130.1 189916 65951.1 1131718 128137.5 561522.3 7339765 199538.1 304745.4 206322 20162.4 100035 98011.5 102424.5 185913 64587.6 834343.2 305637.3 60749.7]; tmp=c./e; ptmp=tmp; pc=c; eachprice=price./c; x=zeros(1,21); ppx=round(c*0.03); x=ppx+x; pprice=ppx.*eachprice; sum(pprice) px=x; STD=std(ptmp); zongzijin=input('请输入总资金数：（单位\元）'); while sum(eachprice.*px)<=(zongzijin-sum(pprice)) [pnumber pindex]=min(ptmp); px(pindex)=px(pindex)+1; pc(pindex)=pc(pindex)+1;21
2011 精品
ptmp=pc./e; ptemp=std(ptmp); STD=[STD ptemp]; end n=size(STD,2); STD=STD(1,[1:n-1]); [number index]=min(STD); for i=1:index [pnumber pindex]=min(tmp); x(pindex)=x(pindex)+1; c(pindex)=c(pindex)+1; tmp=c./e; end x sumprice=sum(x.*eachprice)+sum(pprice) return
22


本文来源：https://www.bwwdw.com/article/gtbq.html