blochank.m

剑桥大学用于线性和双线性动力学系统辨识的工具箱

function Y = blochank(DATA,M,N)
%BLOCHANK   Block Hankel matrix.
%
%   BLOCHANK(C,M,N) is a block Hankel matrix with block dimension
%      MxN assembled from the block matrix DATA.
%
%   Usage:
%      Y = blochank(C,M,N)
%
%   Inputs: 
%      C    := [C(1);C(2);...;C(N)] Block matrix, i.e. all the C(i) are
%                matrices with the same dimensions.
%      M,N  := Positive integers, with M <= N.
%   
%   Outputs:  
%      Y := Block Hankel matrix, defined as
%
%                [  C(1)    C(2)    C(3)  ... C(N-M+1) ]
%                [  C(2)    C(3)    C(4)  ... C(N-M+2) ]
%                [  C(3)    C(4)    C(5)  ... C(N-M+3) ]
%           Y := [   :       :       :    ...    :     ]
%                [ C(M-2)  C(M-1)   C(M)  ...  C(N-2)  ]
%                [ C(M-1)   C(M)   C(M+1) ...  C(N-1)  ]
%                [  C(M)   C(M+1)  C(M+2) ...   C(N)   ]
%   
%   See also HANKEL, TOEPLITZ, BLOCTOEP.
%
% CUED System Identification Toolbox.
% Cambridge University Engineering Department.
% Copyright (C) 1998-2002. All Rights Reserved.

% Version 1.00, Date: 01/06/2002
% Created by H. Chen and E.C. Kerrigan.

[nrowDATA,ncolDATA] = size(DATA);

if M <= 0 | N <= 0 
  error('M and N are not positive integers.')
end

if M > N
  error('M is larger than N.')
end

if rem(nrowDATA,N) ~= 0
   error('N does not divide into number of rows in DATA.')
end    
 
rows = nrowDATA/N;
  
Y = zeros(M*rows,(N-M+1)*ncolDATA);
 
for j=1:N-M+1
        
  Y(:,(j-1)*ncolDATA+1:j*ncolDATA)=DATA((j-1)*rows+1:(j-1+M)*rows,:);
	 
end
 
% *** last line of blochank.m ***