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

v4 set_shocks.m: added M_.maximum_lag to the periods index 

git-svn-id: https://www.dynare.org/svn/dynare/dynare_v4@725 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
michel 2006-04-29 14:40:35 +00:00
parent eca3df3338
commit 7d0b4fa90e
1 changed files with 16 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
matlab/set_shocks.m
View File

@ -1,18 +1,25 @@
% Copyright (C) 2003 Michel Juillard
%
% set
function set_shocks(flag,k,ivar,values)
global oo_
global oo_ M_
n = size(oo_.exo_simul,1);
if k(end) > n
oo_.exo_simul = [oo_.exo_simul; ones(k(end)-n,1)*oo_.exo_steady_state'];
k = k + M_.maximum_lag;
n1 = size(oo_.exo_simul,1);
n2 = size(oo_.exo_det_simul,1);
if k(end) > n1 & flag <= 1
oo_.exo_simul = [oo_.exo_simul; repmat(oo_.exo_steady_state',k(end)-n1,1)];
elseif k(end) > n2 & flag > 1
oo_.exo_det_simul = [oo_.exo_det_simul; repmat(oo_.exo_steady_state',k(end)-n2,1)];
end
if flag == 0
oo_.exo_simul(k,ivar) = ones(length(k),1).*values;
else
switch flag
case 0
oo_.exo_simul(k,ivar) = repmat(values,length(k),1);
case 1
oo_.exo_simul(k,ivar) = oo_.exo_simul(k,ivar).*values;
case 2
oo_.exo_det_simul(k,ivar) = repmat(values,length(k),1);
case 3
oo_.exo_det_simul(k,ivar) = oo_.exo_det_simul(k,ivar).*values;
end
% 05/29/03 MJ