stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Check that stationary and/or non stationary components exist.

time-shift
Stéphane Adjemian(Charybdis) 2018-04-26 16:35:32 +02:00
parent f5ef781b14
commit e8829eab8e
1 changed files with 10 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
matlab/pac-tools/hVectors.m
View File

@ -24,10 +24,18 @@ n = length(H);
tmp = eye(n*m)-kron(G, transpose(H));
h0 = A_1*A_b*((kron(iota(m, m), H))'*(tmp\kron(iota(m, m), iota(n, ids))));
if isempty(ids)
h0 = [];
else
h0 = A_1*A_b*((kron(iota(m, m), H))'*(tmp\kron(iota(m, m), iota(n, ids))));
end
if nargout>1
h1 = A_1*A_b*(kron(iota(m, m)'*inv(eye(m)-G), H')*(tmp\kron(iota(m, m), iota(n, idns))));
if isempty(idns)
h1 = [];
else
h1 = A_1*A_b*(kron(iota(m, m)'*inv(eye(m)-G), H')*(tmp\kron(iota(m, m), iota(n, idns))));
end
end
if nargout>2