-
Notifications
You must be signed in to change notification settings - Fork 17
/
mlpderiv.m
47 lines (41 loc) · 1.13 KB
/
mlpderiv.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
function g = mlpderiv(net, x)
%MLPDERIV Evaluate derivatives of network outputs with respect to weights.
%
% Description
% G = MLPDERIV(NET, X) takes a network data structure NET and a matrix
% of input vectors X and returns a three-index matrix G whose I, J, K
% element contains the derivative of network output K with respect to
% weight or bias parameter J for input pattern I. The ordering of the
% weight and bias parameters is defined by MLPUNPAK.
%
% See also
% MLP, MLPPAK, MLPGRAD, MLPBKP
%
% Copyright (c) Ian T Nabney (1996-2001)
% Check arguments for consistency
errstring = consist(net, 'mlp', x);
if ~isempty(errstring);
error(errstring);
end
[y, z] = mlpfwd(net, x);
ndata = size(x, 1);
if isfield(net, 'mask')
nwts = size(find(net.mask), 1);
temp = zeros(1, net.nwts);
else
nwts = net.nwts;
end
g = zeros(ndata, nwts, net.nout);
for k = 1 : net.nout
delta = zeros(1, net.nout);
delta(1, k) = 1;
for n = 1 : ndata
if isfield(net, 'mask')
temp = mlpbkp(net, x(n, :), z(n, :), delta);
g(n, :, k) = temp(logical(net.mask));
else
g(n, :, k) = mlpbkp(net, x(n, :), z(n, :),...
delta);
end
end
end