我模拟了以下命令:
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(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!
如何正确解释模拟结果?
答案 0 :(得分:1)
我认为你没有考虑到模型创新e(t)是随机的,即
y(t)= 0.5*y(t-1) + 0.6*u(t-2) -0.2*u(t-3) + e(t)