z=rho*z(-1)+eps_z;
z=rho*z(-1)-eps_z;
% Solutions
w_solution = NaN(10,T);
x_solution = NaN(10,T);
w_solution(:,1) = [0.1; 0.1; 0.1; 0.1; 0.1; 0.1; 0.1; 0.1; 0.1; 0.1];
x_solution(:,1) = inv(Pi_t_11) * w_solution(:,1);
for t = 2:T
x_solution(:,t) = inv(Upsilon_11)*Psi_11*x_solution(:,t-1);
w_solution(:,t) = Pi_t_11*x_solution(:,t);
end
Users browsing this forum: No registered users and 4 guests