matlab建立ar时间序列模型_时序分析基本理论与方法——时间序列模型及其性质...

本文介绍了时间序列模型的种类,如AR、MA和ARMA模型,并重点讨论了AR模型的Yule-Walker方程及自相关、偏自相关函数的性质。此外,提到了ARMA模型的自相关系数特征。
摘要由CSDN通过智能技术生成

时间序列模型及其性质

为了能够描述时间序列之间的内在联系及其发展规律,需要建立时序模型,它能够有效的用于时间序列的预测和控制分析。目前,常用的时间序列模型有自回归模型(Auto Regressive Model,简称 AR 模型),移动平均模型(Moving Average Model,简称 MA 模型),自回归移动平均模型(Auto Regressive Moving Average Model,简称 ARMA 模型 )以及自回归求和移动平均模型(Integrate Auto Regressive Moving Average Model,简称 ARIMA 模型)[37]。考虑到非平稳序列的 ARIMA 模型可通过对非平稳序列做一定阶次的差分,转化为平稳序列的ARMA 模型,下文将不再对非平稳序列的 ARIMA 模型作介绍。

1. 自回归模型—AR 模型

f48a34c419e459e0d32925353c20fa32.png

a316289ae503d28bec6625ade28e15ac.png

212b9250f3bf52c71a98a022414a0418.png

时序模型的性质可从自相关函数(ACF)、偏相关函数(PACF)角度分析。AR模型的自相关系数可由下列方程组求得:

c3224128b79a8caec20b23baf7fd9b12.png

上述联立的方程组称为 Yule-Walker 方程,解方程得:

2713c04f552c4f1ddfa7e61949bbc0e7.png

通常,AR 模型的自相关系数具有拖尾性,即其自相关系数(ACF)在某一步之后不为零,而是按某一形式缓慢衰减;而 AR 模型的偏自相关系数(PACF)不同于其自相关系数的变化趋势,它将在 p 步以后以某种方式迅速衰减为零,这个变化过程可视其具有截尾性。

2 移动平均模型—MA 模型

60fb4daf57ad72545a9a3136e531bbcc.png

MA(q)模型的自相关系数表示如下:

75284d17e547fd6b826718e9a33beaaf.png

531fbd0d940cbd8b4b5f4fbbef39094b.png

3 自回归移动平均模型—ARMA 模型

e2c0ebb3e7f2b6f107fc9cce609c6168.png

da0a6cb0fbde7556a3aa3aec70b486f1.png

557de7a5025e4ac084f86eb2cf40cc52.png

对 ARMA(p、q)而言,其自相关系数和偏自相关系数都具有拖尾性,可根据这一特性判断时间序列是否具有ARMA(p、q)模型结构。

《来源科技文献,经本人分析整理,以技术会友,广交天下朋友》

MATLABAR模型功率谱估计中AR阶次估计的实现-psd_my.rar (最近看了几个关于功率谱的问题,有关AR模型的谱估计,在此分享一下,希望大家不吝指正) (声明:本文内容摘自我的毕业论文——心率变异信号的预处理及功率谱估计) (按:AR模型功率谱估计是对非平稳随机信号功率谱估计的常用方法,但是其模型阶次的估计,除了HOSA工具箱里的arorder函数外,没有现成的函数可用,arorder函数是基于矩阵SVD分解的阶次估计方法,为了比较各种阶次估计方法的区别,下面的函数使用了'FPE', 'AIC', 'MDL', 'CAT'集中准则一并估计,并采用试验方法确定那一个阶次更好。) ………………………………以上省略…………………………………………………………………… 假设原始数据序列为x,那么n阶参数使用最小二乘估计在MATLAB中实现如下: Y = x; Y(1:n) = []; m = N-n; X = [];% 构造系数矩阵 for i = 1:m     for j = 1:n         X(i,j) = xt(n i-j);     end end beta = inv(X'*X)*X'*Y'; 复制代码 beta即为用最小二乘法估计出的模型参数。 此外,还有估计AR模型参数的Yule-Walker方程法、基于线性预测理论的Burg算法和修正的协方差算法等[26]。相应的参数估计方法MATLAB中都有现成的函数,比如aryule、arburg以及arcov等。 4.3.3 AR模型阶次的选择及实验设计 文献[26]中介绍了五种不同的AR模型定阶准则,分别为矩阵奇异值分解(Singular Value Decomposition, SVD)定阶法、最小预测定误差阶准则(Final Prediction Error Criterion, FPE)、AIC定阶准则(Akaika’s Information theoretic Criterion, AIC)、MDL定阶准则以及CAT定阶准则。文献[28]中还介绍了一种BIC定阶准则。SVD方法是对Yule-Walker方程中的自相关矩阵进行SVD分解来实现的,在MATLAB工具箱中arorder函数就是使用的该算法。其他五种算法的基本思想都是建立目标函数,阶次估计的标准是使目标函数最小化。 以上定阶准则在MATLAB中也可以方便的实现,下面是本文实现FPE、AIC、MDL、CAT定阶准则的程序(部分): for m = 1:N-1    ……       % 判断是否达到所选定阶准则的要求    if strcmp(criterion,'FPE')        objectfun(m 1) = (N (m 1))/(N-(m 1))*E(m 1);    elseif strcmp(criterion,'AIC')        objectfun(m 1) = N*log(E(m 1)) 2*(m 1);    elseif strcmp(criterion,'MDL')        objectfun(m 1) = N*log(E(m 1)) (m 1)*log(N);    elseif strcmp(criterion,'CAT')        for index = 1:m 1            temp = temp (N-index)/(N*E(index));        end        objectfun(m 1) = 1/N*temp-(N-(m 1))/(N*E(m 1));    end        if objectfun(m 1) >= objectfun(m)        orderpredict = m;        break;    end end 复制代码 orderpredict变量即为使用相应准则预测的AR模型阶次。 (注:以上代码为结合MATLAB工具箱函数pburg,arburg两个功率谱估计函数增加而得,修改后的pburg等函数会在附件中示意,名为pburgwithcriterion) 登录/注册后可看大图 程序1.JPG (35.14 KB, 下载次数: 20352) 下载附件  保存到相册 2009-8-28 20:54 上传 登录/注册后可看大图 程序2.JPG (51.78 KB, 下载次数: 15377) 下载附件  保存到相册 2009-8-28 20:54 上传
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值