compth.m

新的神经网络算法源程序

function th=compth(t,d,w,c,b, nn)
%Compth: Compose parameter THETA of WaveNet
%
%        THETA=compth(t,d,w,c,b, nn)
%
%        t: translation parameters
%        d: dilation parameters
%        w: linear weights
%        c: direct linear connections
%        b: bias
%        nn = [na nb nk]   (See System Identification Toolbox)
%
%See also decompth.m

%By Qinghua Zhang. March, 1994.

if nargin<5, help compth; error(''); end

[tl, tc] = size(t);
[dl, dc] = size(d);
[wl, wc] = size(w);
[cl, cc] = size(c);
[bl, bc] = size(b);

% TEST consistancy of the dimensions  
if (tl~=dl)|(wl>1)|(cl~=1)|(bl~=1)|(tc~=dc)|(tc~=wc)|(cc~=tl&tl~=0)|(bc~=1)
     error('Inconsistent dimensions of parameters.');
end


nbwavelon = length(w);

nbvar = length(c);
thr = 2*nbvar + 5;

if nargin<6, nn=[0 ones(1,nbvar) zeros(1,nbvar)]; end

nnlen = size(nn,2);
thc = max([3+nnlen nbwavelon  nbvar+1 6]);

th = zeros(thr, thc);

th(1,1:(3+nnlen)) = [nbwavelon nbvar nnlen nn];
if nbwavelon>0
  th(2:nbvar+1,1:nbwavelon) = t;
  th(nbvar+2:2*nbvar+1,1:nbwavelon) = d;
  th(2*nbvar+2,1:nbwavelon) = w;
end
th(2*nbvar+3, 1:nbvar) = c;
th(2*nbvar+3, nbvar+1) = b;
th(2*nbvar+4, 1:6) = clock;
th(2*nbvar+5,1:6) = abs('COMPTH');