datamat.m

来自「递归贝叶斯估计的工具包」· M 代码 · 共 45 行

% DATAMAT  Packs a vector of data (length N) into a datamatrix of dimension M-by-(N-M+1)
%
%   D = datamat(x,M)
%
%   INPUT
%           x       vector of data
%           M       data matrix window (frame) size
%   OUTPUT
%           D       M-by-(N-M+1) datamatrix
%
%
%   Example    D = datamat([1 2 3 4 5 6 7 8 9],3)
%
%     will generate the following datamatrix,
%
%     D = | 3  4  5  6  7  8  9 |
%         | 2  3  4  5  6  7  8 |
%         | 1  2  3  4  5  6  7 |
%
%

%   Copyright  (c) Rudolph van der Merwe (2002)
%
%   This file is part of the ReBEL Toolkit. The ReBEL Toolkit is available free for
%   academic use only (see included license file) and can be obtained by contacting
%   [e-mail].  Businesses wishing to obtain a copy of the software should
%   contact [e-mail] for commercial licensing information.
%
%   See LICENSE (which should be part of the main toolkit distribution) for more
%   detail.

%===============================================================================================

function dm=datamat(x,M);

N=length(x);

dm=zeros(N,M);

i=0:N-M;
ii=M:-1:1;
im=repmat(ii',1,N-M+1)+repmat(i,M,1);

dm=x(im);