如何在Matlab中解释ARX进程的模拟结果？

Question

我模拟了以下命令：

process_arx= idpoly([1 -0.5],[0 0 0.6 -0.2],1,1,1,'NoiseVariance',0.05,'Ts',1);
uk=idinput(2555,'prbs',[0 0.2], [-1 1]);
yk=sim(process_arx,uk,'Noise');

uk和yk的前10个是：

[uk(1:10) yk(1:10)]

ans =

   -1.0000    0.1202
   -1.0000    0.4702
   -1.0000   -0.8700
   -1.0000   -0.6422
   -1.0000   -0.6498
   -1.0000   -1.0173
   -1.0000   -1.0056
   -1.0000   -0.8262
   -1.0000   -0.0129
   -1.0000    0.2128

据我所知，我正在做的是模拟以下方程：

y(t)= 0.5*y(t-1) + 0.6*u(t-2) -0.2*u(t-3)

例如，如果我想要y(4)，我必须计算：y(4)=0.5*y(3) + 0.6*u(2) -0.2*u(1)

从值: y(3)= -0.87；u(2)= -1和u(1) = -1

y(4)应该是-0.835，而模拟结果是-0.6422！

如何正确解释模拟结果？

Answer 1

我认为你没有考虑到模型创新e(t)是随机的，即

y(t)= 0.5*y(t-1) + 0.6*u(t-2) -0.2*u(t-3) + e(t)