lagx.m (316B)
1 % LagX.m lagx(x,lag): Returns [x x(-1) ... x(-lag)] 2 % Note: Of course, this matrix is (T-Lag)x(Lag+1) 3 4 function [x,xlag]=lagx(x,lag); 5 if (lag~=0)&(~isempty(x)); 6 n=length(x); 7 xlag=zeros(n,lag); 8 for i=1:lag; 9 xlag(:,i)=[zeros(i,1); x(1:n-i)]; 10 end % i 11 x = x(1+lag:n); 12 xlag = xlag(1+lag:n,:); 13 end % if