forked from js05212/CDL
-
Notifications
You must be signed in to change notification settings - Fork 0
/
sdae_get_hidden.m
53 lines (45 loc) · 1.58 KB
/
sdae_get_hidden.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
48
49
50
51
% sdae_get_hidden
% Copyright (C) 2011 KyungHyun Cho, Tapani Raiko, Alexander Ilin
%
%This program is free software; you can redistribute it and/or
%modify it under the terms of the GNU General Public License
%as published by the Free Software Foundation; either version 2
%of the License, or (at your option) any later version.
%
%This program is distributed in the hope that it will be useful,
%but WITHOUT ANY WARRANTY; without even the implied warranty of
%MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
%GNU General Public License for more details.
%
%You should have received a copy of the GNU General Public License
%along with this program; if not, write to the Free Software
%Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
%
function [h_mf] = sdae_get_hidden(my, mask_output, x0, S, target_sparsity)
if nargin < 5
target_sparsity = 0;
end
layers = S.structure.layers;
n_layers = length(layers);
h_mf = x0;
for l = 2:n_layers
h_mf = bsxfun(@plus, h_mf * S.W{l-1}, S.biases{l}');
if my.dropout~=0 && l~=n_layers
% recover from dropout
h_mf = h_mf.*(1-my.dropout);
end
if my.dropout~=0 && l==n_layers && mask_output
% recover from dropout
h_mf = h_mf.*(1-my.dropout);
end
if l < n_layers || S.bottleneck.binary
h_mf = sigmoid(h_mf, S.hidden.use_tanh);
end
end
if S.bottleneck.binary
if target_sparsity > 0
avg_acts = mean(h_mf, 1);
diff_acts = max(avg_acts - (1 - target_sparsity), 0);
h_mf = min(max(bsxfun(@minus, h_mf, diff_acts), 0), 1);
end
end