Matlab：索引向量的正确方法;索引时间序列的矩阵维度错误

Question

我试图实现此公式：

其中y（t）是自回归过程AR（2）的输出。以相等的间隔t = 1,2，...，5000个数据采样。 p = 2。该代码引发错误

Error using - 
Matrix dimensions must agree. 

这是因为我不知道如何表示存储在变量L1索引符号的滞后。有人可以帮助缓解这个错误，通过显示正确的方式来索引？这就是我所做的：

var_eps = 1; 
epsilon = sqrt(var_eps)*randn(5000,1); % Gaussian signal exciting the AR model 
Y(1) = 0.0; 
    Y(2) = 0.0; 
    for n= 3:5000 
    Y(n)= 0.1950*Y(n-1) -0.9500*Y(n-2)+ epsilon(n); %AR(1) model 
    end 
p=2; 
t= p+1; 
rho0 = randn(); %c 
phi1 = 0.1950; 
phi2 = -0.9500; 
L1 = (Y(t:end) - rho0 - phi1.*Y(2:end) - phi2.*Y(1:end)).^2; 

L11 = sum(L1); 

Answer 1

这是因为Y(t:end)，Y(2:end)和Y(1:end)都是不同的长度，所以你不能增加或减少它们。可能更易于使用for循环：

L1 = 0; 
for t=p+1:5000 
    L1 = L1 + (Y(t) -rho0 - phi1*Y(t-1) - phi2*Y(t-2))^2; 
end