matlab建立ar时间序列模型_请教matlab下建立AR模型的问题(附代码)

小弟刚开始学ar模型,打算用此模型根据已有的101时间序列个数据,预测接下来的数据,在网上查了许多资料,但是都没有相对比较完整的代码,我自己连拼带凑写了几行,不知对不对,附带把不清楚的地方用红色字体标注出来了,麻烦各位大侠指点一下,十分感谢!

x=load('d:/traffic.txt');

y=aryule(x,4);%建立4阶模型,求出  y[n] = a1y[n-1] + a2y[n-2] + .....+ap[yn-p]   的预测系数a(p)

%这一步在matlab命令窗口执行后显示如下信息:

%y =1.0000   -1.0202    0.0799    0.1174   -0.1629

z=idpoly([y],[]);%把多项式转换成模型

%这一步在matlab命令窗口执行后显示如下信息:

%Discrete-time IDPOLY model: A(q)y(t) = e(t)

%A(q) = 1 - 1.02 q^-1 + 0.07993 q^-2 + 0.1174 q^-3 - 0.1629 q^-4

%

%This model was not estimated from data. (Question1:这是什么意思?难道我的数据没有用到吗?)

%Sampling interval: 1

m=iddata([x],[]);%Question2:这一行是什么用处?

%这一步在matlab命令窗口执行后显示如下信息:

%Time domain data set with 101 samples.

%Sampling interval: 1

%

%Outputs      Unit (if specified)

%  y1

n=ar(x,4,'yw');%Question3:这样就把模型建立好了吗?

%这一步在matlab命令窗口执行后显示如下信息:

%Discrete-time IDPOLY model: A(q)y(t) = e(t)

%A(q) = 1 - 1.02 q^-1 + 0.07993 q^-2 + 0.1174 q^-3 - 0.1629 q^-4

%

%Estimated using AR ('yw'/'ppw') from data set x

%Loss function 3.08869e+008 and FPE 3.32402e+008

%Sampling interval: 1

compare(n,m,1);%向前预测一个时间单位的数据。

运行结束后显示图片:

untitled.jpg (73.98 KB, 下载次数: 130)

2009-10-6 16:43 上传

对于图片的显示,Question4:只显示了101个时间单位,也就是没有显示出来预测的第102个数据,怎么才能让它显示出来?

还有就是上边两条曲线相差一个时间单位,Question5:怎么把结果的原数据和预测数据在相同的时间单位上显示出来?体现出二者的对比?

下边附件是实验用的数据文件

traffic.txt

(797 Bytes, 下载次数: 8265)

2009-10-6 16:50 上传

点击文件名下载附件

这是实验用的数据文件

[本帖最后由 rockleo 于 2009-10-6 16:51 编辑]

  • 0
    点赞
  • 0
    评论
  • 0
    收藏
  • 一键三连
    一键三连
  • 扫一扫,分享海报

表情包
插入表情
评论将由博主筛选后显示,对所有人可见 | 还能输入1000个字符
相关推荐
©️2020 CSDN 皮肤主题: 数字20 设计师:CSDN官方博客 返回首页
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值