Daniel@0: function net = rbfunpak(net, w)
Daniel@0: %RBFUNPAK Separates a vector of RBF weights into its components.
Daniel@0: %
Daniel@0: %	Description
Daniel@0: %	NET = RBFUNPAK(NET, W) takes an RBF network data structure NET and  a
Daniel@0: %	weight vector W, and returns a network data structure identical to
Daniel@0: %	the input network, except that the centres C, the widths WI, the
Daniel@0: %	second-layer weight matrix W2 and the second-layer bias vector B2
Daniel@0: %	have all been set to the corresponding elements of W.
Daniel@0: %
Daniel@0: %	See also
Daniel@0: %	RBFPAK, RBF
Daniel@0: %
Daniel@0: 
Daniel@0: %	Copyright (c) Ian T Nabney (1996-2001)
Daniel@0: 
Daniel@0: % Check arguments for consistency
Daniel@0: errstring = consist(net, 'rbf');
Daniel@0: if ~errstring
Daniel@0:   error(errstring);
Daniel@0: end
Daniel@0: 
Daniel@0: if net.nwts ~= length(w)
Daniel@0:   error('Invalid length of weight vector')
Daniel@0: end
Daniel@0: 
Daniel@0: nin 	= net.nin;
Daniel@0: nhidden = net.nhidden;
Daniel@0: nout 	= net.nout;
Daniel@0: 
Daniel@0: mark1 = nin*nhidden;
Daniel@0: net.c = reshape(w(1:mark1), nhidden, nin);
Daniel@0: if strcmp(net.actfn, 'gaussian')
Daniel@0:   mark2 = mark1 + nhidden;
Daniel@0:   net.wi = reshape(w(mark1+1:mark2), 1, nhidden);
Daniel@0: else
Daniel@0:   mark2 = mark1;
Daniel@0:   net.wi = [];
Daniel@0: end
Daniel@0: mark3 = mark2 + nhidden*nout;
Daniel@0: net.w2 = reshape(w(mark2+1:mark3), nhidden, nout);
Daniel@0: mark4 = mark3 + nout;
Daniel@0: net.b2 = reshape(w(mark3+1:mark4), 1, nout);