出版社资源优化配置的线性规划模型（06年）
【摘要】：本文讨论出版社的资源优化配置问题，即出版社在总书号数一定的条件下的书号分配问题，使出版社的利润最大。根据近五年来的数据进行筛选与分析，得出以分社为单位的书号分配方案，和以课程为单位的进一步细化书号分配方案。
由于出版社在定价时保持对所有教材利润率同一，在此原则上制定教材单价，这就意味着利润是由销售总额唯一确定的，销售总额越大，利润也越大。问题转化为求销售总额的最大值的规划问题。
由附件二（读者调查数据表），我们建立了一套市场竞争力评估体系：我们知道市场竞争力的关键因素有市场占有率、出版社在读者心目中的地位（以下简称地位）、读者满意度等，于是把市场竞争力定义为以上三个关键因素的线性组合，这既符合实际又于数据相吻合。通过对历年市场占有率、地位和读者满意度的数据统计分析，我们发现市场占有率是关于年份的线性函数，而地位和读者满意度几年来没有明显变化，故将2006年的地位和读者满意度取为前几年的均值。由此我们得出市场竞争力也应是关于年份的线性函数。
在销量关于书号数和市场竞争力的关系的确定问题上，通过对销量和书号数的统计分析，计算出五年来各个分社单位书号销量，用Matlab作图观察后，我们发现单位书号销量是随年份呈指数增长形势的，因此我们大胆推测单位书号的销量与市场竞争力是成指数型相关的。
销售总额是各分社销售额的和，分社销售额是分社销量与分社均价的乘积，而分社均价我们定义成分社各课程均价的加权平均数，其权就是各课程的销量。2006年的分社均价就取为前五年的分社均价的平均值。这样方便了销售总额的计算和模型的线性化求解，当然也可能产生一定的误差。
为了对这种情况加以改进，我们又建立了更精确的数学模型，站在A出版社的高度直接分配每门课程的书号数，进一步细化分析。
最后我们对模型的可行性，合理性，科学性进行了阐述，得到对模型的整体评价以及需改进之处。
【关键词】：线性规划 海量数据 数据筛选 数据统计分析
数值模拟 市场竞争力评估 单位书号销量 分社均价
1
一. 问题重述
本题研究A出版社,重点关注九个分社72门学科教材的出版销售状况，在以五年市场调研数据为背景的情况下，预测2006年的市场信息,将总量一定的书号数(500)分配给各分社，使得出版社的利润最大化。
其中书号数的分配主要以产生最大经济效益即利润的原则进行，而教材定价是建立在利润率不变的情况下的，因此可以认为利润的大小只决定于销售额的大小，而销售额等于课程销售量乘以课程均价，书本的实际销售量由总社分配给各个分社的书号数决定，但计算方法没有明确给出。经分析，实际销售量是关于书号数和市场竞争力的函数，因此准确的找出这个函数以及市场竞争力的发展规律并正确预测出2006年的市场竞争力是决定分配书号数的关键问题。这需要对前五年的数据进行分析得出函数，以及市场竞争力的发展规律，求出出版社的书号数分配方案，使得2006年的出版社实际销售总额。
二．问题分析
本题的目的是给出版社提出一个最优的方案，即确定分配给各个分社的书号数，使得出版社的利润最高。由于出版社在定价时保持对所有教材利润率同一，在此原则上制定教材单价，这就意味着利润是由销售总额唯一确定的，销售总额越大，利润也越大。问题转化为求销售总额的最大值的规划问题。
以各个分社为单位，销售额是以销售量与价格的乘积来定的。价格在本题中是一个已知量，而销售量是关于书号数与市场竞争力的未知函数，在假设这个未知函数不根据时间的改变而改变的情况下，可根据五年的数据确定这个未知函数。
市场竞争力是影响实际销售量的一个重要因素，它的准确与否与计算的准确度有莫大的关系。这需要基于问卷调查的数据评估出来，而市场占有率，A社在读者心中的位置和其它因素都是市场竞争力影响因素。
市场占有率是在各个分社的基础上，求出每年每个分社的出版物市场占有率，以计算每年各个分社的销售册数与相应类的出版社总销售册数之比，取值范围应在[0,1]之间，分析其市场占有率的发展规律，预计2006年市场占有率。
A社在读者心中的地位（以下简称地位），因地位是数字上的越低相当于其分数值越高，故基于计算和实际意义的考虑，取每年平均地位数据的倒数来计算，取值范围也在为[0,1] 之间，再以此估计第六年的情况。
读者满意度是具体到各个分社的各个学科的，因市场占有率和地位的取值范围都在[0.1]之间，故读者满意度的计算：应先计算每位读者对教材内容、教材的作者、印刷及排版质量和教材价格的平均满意度，范围是[0,5]，再除以5，使得其落在[0,1]之间，再计算所有读者的平均满意度，范围仍在[0,1]之间。以读者满意度的发展规律来预计2006年的读者满意度。
其它因素诸如学校地域性，读者的专业，课程的性质，教材获得方式，新旧书等，在实际计算出版社的市场竞争力问题上影响极其微弱，可忽略，并将大大降低模型的复
2
杂程度。于是得出市场竞争力是关于市场占有率，地位和读者满意度的线性组合。最终预测2006年的市场竞争力。
三．符号说明：
1，Y?(Y1,Y2,Y3,Y4,Y5,Y6)：A出版社每年的销售总额。
?X11?X212，X??Xji???????X51X12?X22?X52X19???X29?： A出版社分配给各个分社的书号数矩阵 。
?????X59?3，X6?(X61,X62,?,X69)：2006年A出版社分配给各个分社的书号数。
'''4，X6'?(X61,X62,?,X69)：2006年各分社申请的书号数。
?Z11?Z215，Z =?Zji?=?????Z51?P11?P216，P??Pji???????P59?R11?R21?7，R??Rji??????R61Z12?Z19??Z22?Z29?：第j 年第i个分社的市场竞争力。 ?????Z52?Z59?P?P1219??P22?P29?：第j 年第i个分社的市场占有率。
?????P59?P59?R12?R19??R22?R29?：第j年第i个分社出版物的平均价格。
?????R62?R69?8，M：A出版社的总书号数，经计算M=500。 9，K??K1,K2,?,K9?：各出版社计划的准确度。 10，
Q：A出版社在读者心目中的地位，用排名平均数的倒数。
3
11，
?B11?B21?B??Bji??????B61?y11?yy??yji???21????y51B12?B19??B22?B29?：读者对各个分社的平均满意度。
?????B62?B69?y12?y22?y52y19???y29?：第j 年第i个分社的实际销售量。 ?????y59?12，
13，
y6??y61,y62,?,y69?：2006年每个分社的销售量。
?y11??yyjk??21???y?51y12?y22???y51?y1,72??y2,72?：第j年第k门课程的实际销售量。 ???y5,72??14，
15， N?(N1,N2,?,N9)：2006年第i个分社的人力资源工作能力。
四．模型假设：
1，读者调查数据基本可信。
2，学校地域性，读者的专业，课程的性质，新旧书等因素对评估出版社的影响微弱，可以不做考虑。
3，2006年的总书号数与前五年是一致的。
4， 分社销量yji与Xji及Zji的函数关系式yji?f(Xji,Zji)每年是不变的。
五．模型的建立：
一， 销售额函数的建立：
由于2001~2005年的各个分社的销售量y已经给出，各个分社的出版物平均价格R可以得到，其计算方法为以分社为单位，相应课程的销售量与课程的均价乘积求平均值，即第j年第i个分社出版物的平均价格
1Rji?yji
i分社的课程k4
?yjkRk，
2006年的平均价格R6在对每年平均价格数据的分析上得出尽可能准确的估计。 故第j年的销售额Yj为第j年各分社销售量yj与第j年分社的课程平均价格Rj的乘积的和，即销售额
?Rj1???Rtj2?Yj=?yj??Rj?=?yj1,yj2,?,yj9??
?????R???j9?从而得出Y?(Y1,Y2,Y3,Y4,Y5)。 本题中2006年的销售总额
?R61???RtY6=?y6??R6?=?y61,y62,?,y69??62?
?????R?69?二， 销量函数的建立：
从销售额的计算中，我们知道分社销量函数的确定是至关重要的。而分社销量
确定分社销量yji是与分社分配得到的书号数Xji以及分社的市场竞争力Zji相关的，
yji与Xji及Zji的函数关系式yji?f(Xji,Zji)也必须基于2001~2005年的分析上。由于分社销量yji与分社分配得到的书号数Xji是已知的，故问题转化为市场竞争力Zji与分社销量yji与分社分配得到的书号数Xji的函数关系式Zji?v(yji,Xji)，由假设4可知Zji?v(yji,Xji)也是不随年份的变化而变化的。 三， 市场竞争力函数的建立：
市场竞争力Zji是根据读者调查数据的分析上，是三个关键因素市场占有率Pji，
读者满意度B， A出版社在读者心目中的地位Q的线性组合，即
Zji?r1iPji?r2iBj?r3iQji,(r1i，r2i，r3i?0)，再由Zji?v(yj,iX)可ji以确定出权r1ji,r2ji,r3ji，得出市场竞争力函数，由于函数不随年份而改变的，所以2006的各个分社的市场竞争力可以做出预测。
四， 模型建立[1]：
5
Max Y6=y6R??X6ieZ6iRit6i?19?9??X6i?500? i?1??s.t ?X6i?Ni,(i?1,2,?,9) ?1''?X6i?X6i?X6i,(i?1,2,?,9)?2??X6i是正整数,(i?1,2,?,9)
六．模型的求解
一：数据的准备与分析：
从附件二读者调查表中得到市场占有率，A社在读者心目中的地位，读者满意度等市场竞争力的关键因素，数据如下：
1， 市场占有率P：
因A社关心的只有1～72门，故市场占有率的计算方法[2]如下：先在“Q2a课程名称（Q2a）”中筛选出1～72门课程的数据，然后在“Q2c学科类别(Q2c)”中分别筛选出九个学科列别，统计“Q2d使用册数(Q2d)”的总册数，这就是所有出版社的此学科的总销售量，再从“Q2g出版社代码(Q2g)”筛选“P115”统计出A出版社的销售量，则
第j年学科类（A出版社第i分社）的销售量， Pji?第j年全部出版社相应学科类的总销售量就得到了第j年第i分社的市场占有率。数据如下：
2001 2002 2003 2004 2005 分社一 分社二 分社三 分社四 分社五 分社六 分社七 分社八 分社九 0.1399 0.1472 0.1644 0.1663 0.1724 0.2619 0.2735 0.2943 0.3059 0.3198 0.5586 0.5735 0.6036 0.6190 0.6366 0.0594 0.0622 0.0616 0.0722 0.0826 0.3350 0.3442 0.3495 0.3684 0.4027 0.5785 0.6099 0.6637 0.6861 0.7758 0.4362 0.4894 0.4894 0.5319 0.5851 0.8598 0.8785 0.8879 0.8785 0.9065 0.7340 0.7660 0.7660 0.8617 0.9574 按列观察发现每个分社的市场占有率基本都是每年递增的。 求其一阶微分如下：
分社一 分社二 分社三 分社四 分社五 分社六 分社七 分社八 0.0074 0.0116 0.0149 0.0028 0.0092 0.0314 0.0532 0.0187 0.0172 0.0209 0.0301 -0.0006 0.0053 0.0538 0 0.0093 0.0018 0.0062
分社九 0.0319 0 0.0957 0.0957 0.0116 0.0139 0.0154 0.0177 0.0106 0.0104 0.0189 0.0343 0.0224 0.0897 0.0426 0.0532 -0.0093 0.0280 6
其二阶微分如下
分社一 分社二 分社三 分社四 分社五 分社六 分社七 分社八 分社九 0.0098 0.0093 0.0152 -0.0035 -0.0039 0.0224 -0.0532 -0.0093 -0.0319 -0.0153 -0.0093 -0.0147 0.0112 0.0137 -0.0314 0.0426 -0.0187 0.0957 0.0043 0.0023 0.0023 -0.0002 0.0153 0.0673 0.0106 0.0374 0.0000 对每个分社分别观察其微分，发现二阶微分都有正负值并且数据接近于0，一阶微分也非常接近，故猜测每个分社的市场占有率是关于年份的线性函数，采用最小二乘法计算求出各个分社关于年份的函数如下：
Pj1?0.008409j?0.1328，Pj2?0.01483j?0.2466，Pj3?0.02015j?0.5378 Pj4?0.005653j?0.05065，Pj5?0.01596j?0.3121，Pj6?0.04709j?0.5215 Pj7?0.03404j?0.4043，Pj8?0.009346j?0.8542，Pj9?0.05426j?0.6543 由此预测2006年的各个分社的市场占有率如下：
P61 P62 P63 P64 P65 P66 P67 P68 P69 0.183254 0.33558 0.6587 0.084568 0.40786 0.80404 0.60854 0.910276 0.97986
2. A社在读者心目中的地位B：
统计附件二问卷调查数据中的“Q1 A出版社在您心中的位置（Q1）”得出每年A社在读者心目中的平均地位数据如下，
2001 2002 2003 2004 2005 2.5522 2.5533 2.5520 2.5518 3．0102 取其倒数得到B值如下：
2001 0.3322
从表中可以明显看出2002~2005年的数据基本持平不变，但2001年的数据有所差别，可能是2001年A出版社刚出道，日后渐进稳定，故预测2006年的B6应取2002~2005年的均值
2002 0.3918 2003 0.3916 2004 0.3919 2005 0.3919 B6?0.39179
3.每年每个分社的读者满意度Q：
统计附件二读者问卷调查表中的“Q2l:对该书的满意度评价”的平均得分再除以总
7
分5，得到Q如下， 分社一 分社二 分社三 分社四 分社五 分社六 分社七 分社八 分社九 2001 0.3957 0.4048 0.4002 0.4500 0.5994 0.6002 0.6016 0.5991 0.5750 2002 0.5752 0.5649 0.5708 0.5667 0.5693 0.5605 0.5628 0.5570 0.6037 2003 0.5654 0.5646 0.5647 0.5670 0.5666 0.5614 0.5824 0.5692 0.5702 2004 0.5646 0.5744 0.5707 0.5674 0.5676 0.5877 0.5447 0.6042 0.5580 2005 0.5674 0.5667 0.5690 0.5680 0.5684 0.5585 0.5452 0.5682 0.5787
对其直接观察可以看出分社一到分社五，2001年的数据与2002~2005年的数据相比有明显异常，而2002~2005年的数据却相差不大，故Q61,Q62,Q63,Q64,Q65取2002~2005年的数据均值。分社六到分社九的数据按年份波动，取Q66,Q67,Q68,Q69为2001~ 2005年的数据均值。2006年各个分社的读者满意度Q6如下： Q61 0.5681 Q62 0.5677 Q63 0.5688 Q64 0.5673 Q65 0.5680 Q66 0.5737 Q67 0.5673 Q68 0.5795 Q69 0.5771 4，分社均价：
由于销售额是各个分社销售量与相应分社的教材均价的乘积，而分社的分社均价应等于课程的销量与课程均价的乘积的和与分社销量的商，根据附件三的课程销量，附件四的课程均价得到每年每个分社的分社均价数据表如下：
分社一 分社二 分社三 分社四 分社五 分社六 分社七 分社八 分社九 2001 25.8942 29.3284 23.8947 30.0625 17.4437 24.4493 24.1634 24.4757 26.2469 2002 25.8989 29.2741 23.9006 29.7631 17.3680 24.1085 24.1491 24.6916 26.5110 2003 25.8410 28.9136 23.8093 29.7738 17.8541 23.7971 24.2340 24.5114 26.6398 2004 25.9321 28.9883 23.8201 28.9501 17.3874 24.5559 24.2928 24.5930 26.9457 2005 25.8339 29.0554 23.8582 29.2945 17.5923 24.3632 24.1021 24.6539 26.7571
取每个分社5年来的课程均价的均值为2006年各个分社的分社均价 分社一 分社二 分社三 分社四 分社五 分社六 分社七 分社八 分社九 25.88 29.119 23.8566 29.5688 17.5291 24.2548 24.1883 24.5851 26.6201
二：销量函数的求解：
8
销量函数
yji?f(Xji,Zji)得出Zji?v(yji,Xji)，因yji,Xji是已知量，从附件3
得到每年每个出版分社的实际销售量表如下：
分社一 分社二 分社三 分社四 分社五 分社六 分社七 分社八 2001 9843 18752 194518 14490 99522 10714 3655 7054 2002 10292 19605 202077 16082 115423 12548 3985 8055 2003 13363 26449 238052 19488 137233 15754 5073 9823 2004 16749 33899 309841 25961 180784 19020 6564 11715 2005 20636 39716 390024 33515 216924 25977 8409 15434 13934 5541 6322 7313 9825 分社九
从附件4得到每年每个出版分社的实际获得的书号数表如下： 分社一 分社二 分社三 分社四 分社五 分社六 分社七 分社八 分社九
2001 64 37 147 95 44 40 24 22 27 2002 69 39 148 87 48 42 25 21 21 2003 77 43 153 77 42 40 20 23 25 2004 69 39 150 83 50 38 17 27 27 2005 73 41 146 85 44 39 19 27 26 我们假设实际销售量和分得的书号数成正比例关系，则
yjiXji反映了单位书号数的销
售量，通过对五年数据的分析以及画图(以下给出其中两个的插值图)
9
6502600t1 vs. tfit 1600240055022005002000450180040016003501400t1 vs. tfit 1300120025011.522.533.544.5511.522.533.544.55
可见，我们发现
yjiXji关于年份成指数增长关系，根据假设
yji?f(Xji,Zji)，这很可
能是因为市场竞争力Zji逐年增长的原因造成的。因为市场竞争力
Zji?r1iPji?r2iBj?r3iQji，由于市场占有率是关于年份的线性函数，故市场竞争力也是关于年份的线性函数，因此Zji?rriQ1iPji?r2iB?3i
?yji?XjieZji??(1)?Zji?r1iPji?r2iB?r3iQi ???Pji?kij?biyjiXji关于Zji成指数关系，据此设
下面确定各个参数：
五年的市场占有率Pji已计算得到，。常量B和Qi取五年数据的平均值得到。由(1)式得
lnyjiXji?r1i(kij?bi)?r2iB?r3iQi
lnyjiXji可由五年数据计算得到，则r1i、r2i和r3i也可由最小二乘法线性拟合得到，经计算
得到第j年每个分社的市场竞争力评估公式如下：
Zj1?20.27Pj1?1.832B?2.369Q1, Zj2?12.43Pj2?2.476B?3.373Q2
10
Zj3?9.027Pj3?4.303B?0.6496Q3, Zj4?42.91Pj4?4.729B?1.317Q4 Zj5?12.32Pj5?3.071B?4.263Q5, Zj6?4.965Pj6?2.47B?3.03Q6 Zj7?8.867Pj7?0.8644B?1.247Q7, Zj8?6.847Pj8?2.563?10?8B?3.052?10?10Q8
Zj9?3.888Pj9?2.753B?2.639Q9;
由这些评估公式来评估2006年各个分社的市场竞争力如下：
分社一 分社二 分社三 分社四 分社五 分社六 分社七 分社八 分社九 5.7782 7.0561 8.0015 6.2287 8.6494 6.698 6.4421 6.2327 6.4113
再由销量函数
?yji?XjieZji???Zji?r1iPji?r2iB?r3iQi?P?kj?bjiii??
可以得到2006年各个分社的销量。
三，2006年的A社的销售总额的确定：
A社的销售总额应等于各个分社的销售额的和，而各个分社的销售额应等于分社销售量与相应的分社均价的乘积，得总销售额为
Y6??Y6i
i?19四，问题的解答：
11
Max Y6=?Y6i??y6iRi??X6ieZ6iRii?1i?1i?19?X6i?500??i?1??s.t ?X6i?Ni,(i?1,2,?,9) ,?1'?X6i?X6i,(i?1,2,?,9)2??999
书号数初值取为2006年申请书号数的乘以各个分社的正确度。
由于市场竞争力Z6i?i=1,2,?,9?由上已确定，是确定的数，故此规划为线性整数规划，通过用Matlab编程[3]（程序见附录）求解得2006年A出版社分配给各个分社的书
号数如下：
分社一 分社二 分社三 分社四 分社五 分社六 分社七 分社八 分社九 55 66 120 59 由此得到的销售额为2.1538?107
72 68 20 20 20 七．对模型的一步细化讨论
由于以上模型是以分社为单位进行建模分析，市场占有率，满意度是以分社为单位统计计算而不是以课程为单位统计计算，课程均价换算为分社均价等这些方法都会产生一定的误差量，实际上最好的情况是A出版社站在各门课程的高度，直接对书号进行配置，出于这样的考虑，我们对模型进行了进一步的细化分析，建立如下模型：
Max Y6=y6R??X6neZ6nRnt6n?172?72?X6n?500??n?1 ?s.t ?X6k?Ni,(i?1,2,?,9) ,??i分社的课程k?1''?X6n?X6n?X6n,(n?1,2,?,72)?2X6n是A出版社分配给第n门课程的书号数， Rn是第n门课程的均价，
12
'X6n是第n门课程申请的书号数，
Z6n是第n门课程的市场竞争力，其计算方法包括市场占有率，读者满意度与前面的计算方法基本一致，区别在于这里是以课程为单位，前面是以分社为单位分析的，这里就不赘述。
导入数据（数据见附件中），在用线性规划求解模型，得到72门课程的分配得到的书号数如下：
课程号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 书号数 9 9 2 3 3 8 6 3 8 4 8 4 2 4 3 3 6 8 6 4 6 19 35 4 课程号 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 书号数 12 17 6 3 12 6 40 2 2 11 4 8 3 6 5 3 4 10 8 10 6 8 12 14 课程号 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 书号数 9 2 5 16 4 8 4 4 2 2 4 5 4 8 4 8 2 4 8 5 4 3 2 4
由此作出的预计销售额2.4749?107。
计算得各个分社得到的书号数如下：
13
分社一 分社二 分社三 分社四 分社五 分社六 分社七 分社八 分社九 55 48
120 84 72 44 21 30 26 八．建议
通过对题目给出的数据的统计分析，我们建模分析，再细化分析，给出如下认为可行的资源优化配置： 分社一 分社二 分社三 分社四 分社五 分社六 分社七 分社八 分社九 55
48 120 84 72 44 21 30 26 如此A出版社将得到的销售总额为：2.4749?107。
九．模型的评价
1．我们在模型中假设市场占有率和市场竞争力是线性变化的，这大大简化了模型的复杂度，使得模型的可操作性大大加强，也方便了模型的求解。但也产生了些误差，有少部分数据拟合效果不太理想。在需要考虑更高精度的预测时可以把占有率和市场竞争力设为非线性变化的函数。
2．误差来源：在对附件二读者调查问卷的数据分析中会产生一定的误差，这是由于调
查资料涵盖内容的局限性，如果能获得更多相关信息，准确度将会有进一步的提高。
十．参考书目
[1] 姜启源等，《数学模型》（第三版），北京：高等教育出版社，2003 [2] 李贤平等，《概率论与数理统计》，上海：复旦大学出版社，2004 [3] 谭永基等，《数学模型》，上海：复旦大学出版社，1997
14
附录
由于程序很多，仅附上几个代表性程序。
1. 原规划模型程序 clear
z=[5.7782 7.0561 8.0015 6.2287 8.6494 6.698 6.4421 6.2327 6.4113]; z=exp(z) load junjia2.mat; z=junjia2.*z'; load nenli1;
x_=[110 66 222 118 72 76 40 40 40]'; load c1; a=x_.*c1(:,1); b=x_.*c1(:,2); clear c1 load chuzi.mat
[x,fval]=linprog(-z,[eye(9);-eye(9);ones(1,9);eye(9)],[nenli1;-0.5.*x_;500;x_],[],[],[],[],chuzi); fval=-fval;
2.改进的规划模型程序 clear load junjia; load xiao; load A111; z=-junjia.*xiao;
b=[10,20,30,40,48,54,60,66,72]; load nenli1;
a=ones(1,72);a1=zeros(1,72);a2=zeros(1,72);
a3=zeros(1,72);a9=zeros(1,72);a4=zeros(1,72);a5=zeros(1,72);a6=zeros(1,72);a7=zeros(1,72);a8=zeros(1,72);
a1(1:b(1))=a(1:b(1));b1=nenli1(1);
a2(b(1)+1:b(2))=a(b(1)+1:b(2));b2=nenli1(2); a3(b(2)+1:b(3))=a(b(2)+1:b(3));b3=nenli1(3); a4(b(3)+1:b(4))=a(b(3)+1:b(4));b4=nenli1(4); a5(b(4)+1:b(5))=a(b(4)+1:b(5));b5=nenli1(5); a6(b(5)+1:b(6))=a(b(5)+1:b(6));b6=nenli1(6); a7(b(6)+1:b(7))=a(b(6)+1:b(7));b7=nenli1(7); a8(b(7)+1:b(8))=a(b(7)+1:b(8));b8=nenli1(8); a9(b(8)+1:b(9))=a(b(8)+1:b(9));b9=nenli1(9);
15
a10=eye(72);b10=0.5*A111; a11=ones(1,72);b11=500; a12=eye(72);b12=A111;
[a,p,c]=linprog(z,[a1;a2;a3;a4;a5;a6;a7;a8;a9;-a10;a11;a12],[b1;b2;b3;b4;b5;b6;b7;b8;b9;-b10;b11;b12]) t=1;c=cell(1,9); for i=1:9
c{i}=a(t:b(i)); t=b(i)+1; end for i=1:9
d(i)=sum(c{i}); end
3．求读者满意度 %按年度分书号求平均 load('ss1.mat'); for i=1:72
n=find(ss1(:,1)==i); temp=ss1(n,2:5); s1(i)=mean(mean(temp)); end
load('ss2.mat'); for i=1:72
n=find(ss2(:,1)==i); temp=ss2(n,2:5); s2(i)=mean(mean(temp)); end
load('ss3.mat'); for i=1:72
n=find(ss3(:,1)==i); temp=ss3(n,2:5); s3(i)=mean(mean(temp)); end
load('ss4.mat'); for i=1:72
n=find(ss4(:,1)==i); temp=ss4(n,2:5); s4(i)=mean(mean(temp)); end
load('ss5.mat');
16
for i=1:72
n=find(ss5(:,1)==i); temp=ss5(n,2:5); s5(i)=mean(mean(temp)); end
s=[s1;s2;s3;s4;s5];s=s/5; %按年度分类求平均
a=[1,11,21,31,41,49,55,61,67,73]; for i=1:9
n=find(ss1(:,1)>=a(i) & ss1(:,1) n=find(ss2(:,1)>=a(i) & ss2(:,1) n=find(ss3(:,1)>=a(i) & ss3(:,1) n=find(ss4(:,1)>=a(i) & ss4(:,1) n=find(ss5(:,1)>=a(i) & ss5(:,1) p=[p1;p2;p3;p4;p5];p=p/5; 4．求可靠度 load('kkd.mat'); b=[11,22,33,44,53,60,67,74,81]; temp=kkd; 17 temp([2*b-1,2*b],:)=[]; kkd1=[];kkd2=[]; for i=1:72 kkd1=[kkd1;temp(2*i,:)./temp(2*i-1,:)]; end temp=kkd(sort([2*b-1,2*b]),:); for i=1:9 kkd2=[kkd2;temp(2*i,:)./temp(2*i-1,:)]; end 5．求各分社均价 load junjia; load shiji; a=[1,11,21,31,41,49,55,61,67,73]; temp=[junjia,shiji1];junjia1=[]; temp=[[1:72]',temp]; for i=1:9 n=find(temp(:,1)>=a(i) & temp(:,1) jia(:,j-1)=Temp(:,1).*Temp(:,j); end junjia1=[junjia1;sum(jia)./sum(Temp(:,2:6))]; end junjia2=mean(junjia1,2) 18