图书馆借阅量预测数学建模.docVIP

一、问题重述： 某单位资料室近11周图书借阅量如下表所示： 周 周 次 借 阅 量 星 期 1 2 3 4 5 6 7 8 9 10 11 一 942 1142 821 1166 900 731 1098 807 810 862 二 978 854 716 823 735 648 904 760 931 715 三 1039 710 683 697 719 525 821 770 265 788 四 1077 640 419 602 550 651 564 574 560 204 482 五 1291 891 523 1081 992 303 1094 1073 887 411 1063 1.请建立马氏链模型来对12-15周图书借阅量进行预测； 2.请自己再选择另外一类模型进行预测； 3.对两个模型的预测结果进行比较和评价。 二、问题分析： 针对问题一，要建立马氏链模型，必须对所给数据进行状态划分，我们可以对1-11周的数据进行处理，求其每周的图书借阅总量，因所得数据波动围较大，对其划分状态可能会导致模型不准确，因此，求其每周的天平均图书借阅量来进行状态划分并建立状态转移矩阵进行预测。 对于问题二，我们可以对数据进行稳定处理，然后建立一元线性回归模型，分别以周次和经稳定处理后的图书借阅量为自变量和因变量建立一元线性方程，用Matlab软件求解进行预测。 问题三，对以上两种模型预测的结果进行比较和评价，分析其优缺点并将两种模型进行推广。 三、模型假设： 1.新增图书的种类及数量对图书借阅量无影响。 2.不考虑图书馆的制度调整。 3.不考虑借阅图书岗位职工人数变动。 4.特殊天气因素对图书借阅量无影响或影响很小。 5.不考虑单位的特殊安排。 四、定义与符号说明： 1.表示从状态i到状态j的转移概率。 2.是系统从状态i到转移到其他状态的次数。 3.是由状态i到状态j的转移次数。 4.是残差。 5.是相对误差。 7.R是回归方程的决定系数(R是相关系数)。 8.F是统计量。 9．P是与统计量对应的概率值。 10.S是剩余方差。 五、模型的建立与求解： 问题一：基于马氏链模型预测。 1.模型建立： 按照系统的发展，时间离散化为n=0,1,2，...，对每个n，系统的状态用随机变量表示，设可以取k个离散值=1,2，...，k，从=i到=j的概率及=P(=j| =i),即转移概率。如果的取值只取决于的取值及转移概率，而与,...的取值无关，那么这种离散状态按照离散时间的随机转移过程称为马氏链。其具体方法步骤如下: （1）建立指标值的分级标准。 （2）按步骤①所建立的分级标准确定资料序列中各时段指标所对应的状态。 （3）对步骤②所得的结果进行统计，可得不同滞时（步长）马尔科夫链的转移概率矩阵，它决定了指标值状态转移过程的概率法则。 （4）以（m,m+k）表示马氏链“在时刻出现=的条件下，在时刻出现=条件的概率，即： （m,m+k）=P（==） （i,j=1,2,…，N；m,k都是正整数。） 并称为转移概率。 当转移概率（m,m+k）的k=1时，有： =（m,m+k）=P（==）称之为一步转移概率，表示马氏链由状态经过一次转移到达状态转移概率。所有的转移概率可以构成一个转移概率矩阵: （5）分别以前面若干时段的指标值为初始状态，结合其相应的各阶转移概率矩阵即可预测出该时段指标值的状态概率（i∈E）k为滞时（步长），k=1,2,…，m。 （6）将同一状态的各预测概率加权和作为指标值处于该状态的预测概率，即 =（i∈E）max{，i∈E }所对应的i即为该时段指标值的预测状态。从而确定预测值取值模糊区间[,],预测值该区间的中间值，即=（+）/2。 （7）计算预测值与实际值的相对误差，并以此对新模型进行预测精度评估，其计算公式为： =（-）/ 2.模型求解： 记随机变量（n=1,2,...，11）为第n周图书借阅量，表示从状态i到状态j的转移概率，因的取值只取决于的取值及转移概率，故需要对所给数据进行状态划分。 将1～11周的图书借阅量进行统计，求出每周的图书借阅总量为 5327,4237,3162,4369,3896,2858,4481,1647,3784,2621,3910。但第8周的周一至周三无数据，可能是清明节放假，图书馆关闭，并且每周图书借阅总量的数 波动围较大，故可再求出每周平均每天的图书借阅量，数据如下：1065,847,632,874,779,572,896,824,757,524,782。这些数据是位于524～1065之间，可以将这些数据进行如下分组：500～600,600～700,700～800,800～900,900～1000,1000～110