小弟刚开始学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 d