我使用了系统识别工具箱中的ARX函数和 RESID函数,但得到的残差是:
0
0
0
5
6
8
7
8
零的数量=滞后的数量,我需要一个完整的残差向量
N 阶AR 模型需要前 N 个值来预测下一个值,这就是为什么不预测前 N 个值的原因。您始终可以在开头填充向量(通过复制或零),例如:
load twotankdata
order = 5;
m = arx(y, order);
r = resid([y(1:order);y], m);
r = r(order+1:end);