数学建模算法与应用1
1 绪论。1.1 时间序列的概念和分类。
时间序列是按时间顺序排列的、随时间变化且相互关联的数据序列。分析时间序列的方法构成数据分析的一个重要领域,即时间序列分析。
时间序列根据所研究的依据不同,可有不同的分类。
1.按所研究的对象的多少分,有一元时间序列和多元时间序列。
2.按时间的连续性可将时间序列分为离散时间序列和连续时间序列两种。
3.按序列的统计特性分,有平稳时间序列和非平稳时间序列。如果一个时间序列的概率分布与时间t无关,则称该序列为严格的(狭义的)平稳时间序列。如果序列的。
一、二阶矩存在,而且对任意时刻t满足:
1)均值为常数。
2)协方差为时间间隔的函数。
则称该序列为宽平稳时间序列,也叫广义平稳时间序列。我们以后所研究的时间序列主要是宽平稳时间序列。
4.按时间序列的分布规律来分,有高斯型时间序列和非高斯型时间序列。
1.2 时间序列分析方法概述。
时间序列**技术就是通过对**目标自身时间序列的处理,来研究其变化趋势的。一个时间序列往往是以下几类变化形式的叠加或耦合。
1)长期趋势变动。它是指时间序列朝着一定的方向持续上升或下降,或停留在某一水平上的倾向,它反映了客观事物的主要变化趋势。
2)季节变动。
3)循环变动。通常是指周期为一年以上,由非季节因素引起的涨落起伏波形相似的波动。
4)不规则变动。通常它分为突然变动和随机变动。
通常用表示长期趋势项,表示季节变动趋势项,表示循环变动趋势项,表示随机干扰项。常见的时间序列模型有以下几种类型:
1)加法模型。
2)乘法模型。
3)混合模型。
其中:是观测目标的观测记录,均值,方差。如果在**时间范围以内,无突然变动且随机变动的方差较小,并且有理由认为过去和现在的演变趋势将继续发展到未来时,可用一些经验方法进行**。
下面介绍移动平均法。
2 移动平均法。
移动平均法是根据时间序列资料逐渐推移,依次计算包含一定项数的时序平均数,以反映长期趋势的方法。当时间序列的数值由于受周期变动和不规则变动的影响,起伏较大,不易显示出发展趋势时,可用移动平均法,消除这些因素的影响,分析、**序列的长期趋势。 移动平均法有简单移动平均法,加权移动平均法,趋势移动平均法等。
2.1 简单移动平均法。
设观测序列为,取移动平均的项数n2.1.1)
当**目标的基本趋势是在某一水平上下波动时,可用一次简单移动平均方法建立**模型:
其**标准误差为:
最近n期序列值的平均值作为未来各期的测试结果。一般n取值范围:
当历史序列的基本趋势变化不大且序列中随机变动成分较多时,n的取值应该大一些。否则n的取值应小一些。在有确定的季节变动周期的资料中,移动平均的项数应取周期长度。
选择最佳n值的一个有效方法是,比较若干模型的**误差。**标准误差最小者为好。
简单移动平均法只适合做近期**,而且是**目标的发展趋势变化不大的情况。如果目标的发展趋势存在其它的变化,采用简单移动平均法就会产生较大的**偏差和滞后。
2.2 加权移动平均法。
在简单移动平均公式中,每期数据在求平均时的作用是等同的。但是,每期数据所包含的信息量不一样,近期数据包含着更多关于未来情况的信息。因此,把各期数据等同看待是不尽合理的,应考虑各期数据的重要性,对近期数据给予较大的权重,这就是加权移动平均法的基本思想。
设时间序列为;加权移动平均公式为
式中为t期加权移动平均数;为的权数,它体现了相应的在加权平均数中的重要性。
利用加权移动平均数来做**,其**公式为。
即以第t期加权移动平均数作为第t+1期的**值。
在加权移动平均法中,的选择,同样具有一定的经验性。一般的原则是:近期数据的权数大,远期数据的权数小。至于大到什么程度和小到什么程度,则需要按照**者对序列的了解和分析来确定。
2.3 趋势移动平均法。
简单移动平均法和加权移动平均法,在时间序列没有明显的趋势变动时,能够准确反映实际情况。但当时间序列出现直线增加或减少的变动趋势时,用简单移动平均法和加权移动平均法来**就会出现滞后偏差。因此,需要进行修正,修正的方法是作二次移动平均,利用移动平均滞后偏差的规律来建立直线趋势的**模型。
这就是趋势移动平均法。
一次移动的平均数为。
在一次移动平均的基础上再进行一次移动平均就是二次移动平均,其计算公式为。
下面讨论如何利用移动平均的滞后偏差建立直线趋势**模型。
设时间序列从某时期开始具有直线趋势,且认为未来时期也按此直线趋势变化,则可设此直线趋势**模型为。
其中为当前时期数;为由至**期的时期数;为截距;为斜率。两者又称为平滑系数。
现在,我们根据移动平均值来确定平滑系数。由上式可知。
所以。因此。
可得平滑系数的计算公式为。
趋势移动平均法对于同时存在直线趋势与周期波动的序列,是一种既能反映趋势变化,又可以有效地分离出来周期变动的方法。
3案例分析。
我国1965~2023年的发电总量如下表所示,试**2023年和2023年的发电总量。
解:将上表中的数据导入matlab中,得到发电总量趋势散点图如下,从图中可看出发电总量基本呈直线上升趋势,可用趋势移动平均法来**。
取n=6,分别计算一次和二次移动平均值并列于上表中。
由公式(2.3.7),可得t=21时,有。
于是,得t=21时直线趋势**模型为。
**2023年和2023年的发电总量为。
计算的matlab程序如下:
clc,clear
load(''把原始数据保存在纯文本文件中。
m1=length(y);
n=6; %n为移动平均的项数。
for i=1:m1-n+1
yhat1(i)=sum(y(i:i+n-1))/n;
endyhat1
m2=length(yhat1);
for i=1:m2-n+1
yhat2(i)=sum(yhat1(i:i+n-1))/n;
endyhat2
plot(1:21,y,'*g')
xlabel('t');ylabel('发电总量')
title('发电总量趋势散点图')
a21=2*yhat1(end)-yhat2(end)
b21=2*(yhat1(end)-yhat2(end))/n-1)
y1986=a21+b21
y1987=a21+2*b21
程序运行结果如下:
4实验总结
通过本次实验我进一步熟悉了时间序列的相关知识。移动平均法是时间序列分析法中重要的方法之一,但移动平均法只适合做近期**,而且只能**目标的发展趋势变化不大的情况。如果目标的发展趋势存在其他的变化,采用简单移动平均法就会产生较大的**偏差和滞后。
此次实验我还在一定程度上锻炼了我的编程能力及对matlab的使用。
5参考文献。
1]张树京,齐立心。时间序列分析简明教程[m].北京:清华大学出版社,2003
2]吴怀宇。时间序列分析与综合[m].武汉大学出版社。 2004:10-98.
3]王燕。应用时间序列分析[m].北京:中国人民大学出版社,2005:1-239.
5]陈桂明,戚红雨,潘伟,matlab数理统计(北京:科学出版社,2002.
6]韩中庚。数学建模方法及应用。北京:高等教育出版社,2005.
某集团下设两个子公司:子公司a、子公司b。各子公司财务分别独立核算。
每个子公司都实施了对雇员的医疗保障计划,由各子公司自行承担雇员的全部医疗费用。过去的统计数据表明,每个子公司的雇员人数以及每一年龄段的雇员比例,在各年度都保持相对稳定。各子公司各年度的医疗费用支出见下表(附录1)。
数学建模课程设计
高速公路修建费用问题。摘要。本文研究了高速公路修建费用问题。在建设高速公路时,要求建造费用最小,这是个关于最优化求解的问题。根据不同地貌上建造成本的不同,找出a b两地高速路的最低建造成本路线。以确定最便宜的路线为目标建立了非线性规划模型来确定在各个地貌交界处的汇合点。建造总费用为各个地貌中建造公路...
数学建模课程设计
题目 服务机构劳务安排的优化设计。服务机构劳务安排的优化设计。摘要。生活中,我们会接触到各种各样的服务机构,就服务机构自身而言,不同的时间段内需要的服务量和需求量有显著的不同,而且不同时段对劳务的支付工资也有不同。本文主要站在管理者的角度,考虑在既要满足需要,又要尽量节约劳务开支的基础上如何设计才能...
数学建模课程设计
0840503220 苏阳。0840503224 张明。0840503226 郑景旻。影院座位设计。问题回顾 影院座位的满意程度主要取决于视角和仰角,视角是观众眼睛到屏幕上下边缘的视线的夹角,越大越好 仰角是观众眼睛到屏幕上边缘视线与水平线的夹角,太大使人的头部过分上仰,引起不适,一般要求仰角不超过...