fillb.m

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

function [B,k] = fillb(Bind,m,par,ki)
%FILLB   Fills the B matrix with estimates.
%
%   [B,k] = FILLB(Bind,m,par,ki) fills matrix B with estimates beginning in
%   par(ki) and ending in par(k-1). B is a matrix with m columns and as many
%   rows as elements in Bind. Bind(j) contains the index of the first
%   element in the j'th row of B which is non-zero. B is such that the first
%   non-zero element of each row is strictly positive. It is assumed that
%   each row of B is non-zero.
%
%   See also SIGNEXTB.
%
% 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.

if any(Bind > m)
  disp('Indices in Bind are not allowed to be larger than number of columns m in B.')
  error(' ')
end

if size(Bind,2) > 1
  disp('Bind has to be a column vector.')
  error(' ')
end

total = sum(m-Bind+1);

if total > length(par) - ki + 1
  disp('Not enough elements in par to fill B.')
  error(' ')
end

n = length(Bind);
  
B = zeros(n,m);
k = ki;

for j = 1:n % Extract the elements row-wise
  B(j,Bind(j)) = exp(par(k)); % Guarantees that the first non-zero element in
							  % row j is strictly positive. See SIGNEXTB
  k = k + 1;
  for i = Bind(j)+1:m  % Go through the rest of the columns
    B(j,i) = par(k);   % Any sign allowed
	k = k + 1;
  end
end