苏州模板建站定制,网页版传奇怎么开,阿里云备案域名,网站开发的标准流程神经网络向量化 在本节#xff0c;我们将引入神经网络的向量化版本。在前面关于神经网络介绍的章节中#xff0c;我们已经给出了一个部分向量化的实现#xff0c;它在一次输入一个训练样本时是非常有效率的。下边我们看看如何实现同时处理多个训练样本的算法。具体来讲… 神经网络向量化 在本节我们将引入神经网络的向量化版本。在前面关于神经网络介绍的章节中我们已经给出了一个部分向量化的实现它在一次输入一个训练样本时是非常有效率的。下边我们看看如何实现同时处理多个训练样本的算法。具体来讲我们将把正向传播、反向传播这两个步骤以及稀疏特征集学习扩展为多训练样本版本。 Contents [hide] 1 正向传播2 反向传播3 稀疏自编码网络4 中英文对照5 中文译者 正向传播 考虑一个三层网络(一个输入层、一个隐含层、以及一个输出层)并且假定x是包含一个单一训练样本 的列向量。则向量化的正向传播步骤如下 这对于单一训练样本而言是非常有效的一种实现但是当我们需要处理m个训练样本时则需要把如上步骤放入一个for循环中。 更具体点来说参照逻辑回归向量化的例子我们用Matlab/Octave风格变量x表示包含输入训练样本的矩阵x(:,i)代表第个训练样本。则x正向传播步骤可如下实现 % 非向量化实现
for i1:m, z2 W1 * x(:,i) b1;a2 f(z2);z3 W2 * a2 b2;h(:,i) f(z3);
end; 这个for循环能否去掉呢对于很多算法而言我们使用向量来表示计算过程中的中间结果。例如在前面的非向量化实现中z2,a2,z3都是列向量分别用来计算隐层和输出层的激励结果。为了充分利用并行化和高效矩阵运算的优势我们希望算法能同时处理多个训练样本。让我们先暂时忽略前面公式中的b1和b2(把它们设置为0)那么可以实现如下: % 向量化实现 (忽略 b1, b2)
z2 W1 * x;
a2 f(z2);
z3 W2 * a2;
h f(z3) 在这个实现中z2,a2,z3都是矩阵每个训练样本对应矩阵的一列。在对多个训练样本实现向量化时常用的设计模式是虽然前面每个样本对应一个列向量比如z2但我们可把这些列向量堆叠成一个矩阵以充分享受矩阵运算带来的好处。这样在这个例子中a2就成了一个s2 X m的矩阵(s2是网络第二层中的神经元数m是训练样本个数)。矩阵a2的物理含义是当第i个训练样本x(:i)输入到网络中时它的第i列就表示这个输入信号对隐神经元 (网络第二层)的激励结果。 在上面的实现中我们假定激活函数f(z)接受矩阵形式的输入z并对输入矩阵按列分别施以激活函数。需要注意的是你在实现f(z)的时候要尽量多用Matlab/Octave的矩阵操作并尽量避免使用for循环。假定激活函数采用Sigmoid函数则实现代码如下所示: % 低效的、非向量化的激活函数实现
function output unvectorized_f(z)
output zeros(size(z))
for i1:size(z,1), for j1:size(z,2),output(i,j) 1/(1exp(-z(i,j)));end;
end;
end% 高效的、向量化激活函数实现
function output vectorized_f(z)
output 1./(1exp(-z)); % ./ 在Matlab或Octave中表示对矩阵的每个元素分别进行除法操作
end 最后我们上面的正向传播向量化实现中忽略了b1和b2现在要把他们包含进来为此我们需要用到Matlab/Octave的内建函数repmat % 正向传播的向量化实现
z2 W1 * x repmat(b1,1,m);
a2 f(z2);
z3 W2 * a2 repmat(b2,1,m);
h f(z3) repmat(b1,1,m)的运算效果是它把列向量b1拷贝m份然后堆叠成如下矩阵 这就构成一个s2 X m的矩阵。它和W1 * x相加就等于是把W1 * x矩阵译者注这里x是训练矩阵而非向量, 所以W1 * x代表两个矩阵相乘结果还是一个矩阵的每一列加上b1。如果不熟悉的话可以参考Matlab/Octave的帮助文档获取更多信息(输入“help repmat”)。rampat作为Matlab/Octave的内建函数运行起来是相当高效的远远快过我们自己用for循环实现的效果。 反向传播 现在我们来描述反向传播向量化的思路。在阅读这一节之前强烈建议各位仔细阅读前面介绍的正向传播的例子代码确保你已经完全理解。下边我们只会给出反向传播向量化实现的大致纲要而由你来完成具体细节的推导见向量化练习。 对于监督学习我们有一个包含m个带类别标号样本的训练集。 (对于自编码网络我们只需令y(i) x(i)即可, 但这里考虑的是更一般的情况。) 假定网络的输出有s3维因而每个样本的类别标号向量就记为。在我们的Matlab/Octave数据结构实现中把这些输出按列合在一起形成一个Matlab/Octave风格变量y其中第i列y(:,i)就是y(i)。 现在我们要计算梯度项和。对于梯度中的第一项就像过去在反向传播算法中所描述的那样对于每个训练样本(x,y)我们可以这样来计算 在这里表示对两个向量按对应元素相乘的运算译者注其结果还是一个向量。为了描述简单起见我们这里暂时忽略对参数b(l)的求导不过在你真正实现反向传播时还是需要计算关于它们的导数的。 假定我们已经实现了向量化的正向传播方法如前面那样计算了矩阵形式的变量z2, a2, z3和h那么反向传播的非向量化版本可如下实现 gradW1 zeros(size(W1));
gradW2 zeros(size(W2));
for i1:m,delta3 -(y(:,i) - h(:,i)) .* fprime(z3(:,i)); delta2 W2*delta3(:,i) .* fprime(z2(:,i));gradW2 gradW2 delta3*a2(:,i);gradW1 gradW1 delta2*a1(:,i);
end; 在这个实现中有一个for循环。而我们想要一个能同时处理所有样本、且去除这个for循环的向量化版本。 为做到这一点我们先把向量delta3和delta2替换为矩阵其中每列对应一个训练样本。我们还要实现一个函数fprime(z)该函数接受矩阵形式的输入z并且对矩阵的按元素分别执行。这样上面for循环中的4行Matlab代码中每行都可单独向量化以一行新的向量化的Matlab代码替换它不再需要外层的for循环。 在向量化练习中我们要求你自己去推导出这个算法的向量化版本。如果你已经能从上面的描述中了解如何去做那么我们强烈建议你去实践一下。虽然我们已经为你准备了反向传播的向量化实现提示但还是鼓励你在不看提示的情况下自己去推导一下。 稀疏自编码网络 稀疏自编码网络中包含一个额外的稀疏惩罚项目的是限制神经元的平均激活率使其接近某个预设的目标激活率ρ。其实在对单个训练样本上执行反向传播时我们已经考虑了如何计算这个稀疏惩罚项如下所示 在非向量化的实现中计算代码如下 % 稀疏惩罚Delta
sparsity_delta - rho ./ rho_hat (1 - rho) ./ (1 - rho_hat);
for i1:m,...delta2 (W2*delta3(:,i) beta*sparsity_delta).* fprime(z2(:,i)); ...
end; 但在上面的代码中仍旧含有一个需要在整个训练集上运行的for循环这里delta2是一个列向量。 作为对照回想一下在向量化的情况下delta2现在应该是一个有m列的矩阵分别对应着m个训练样本。还要注意稀疏惩罚项sparsity_delta对所有的训练样本一视同仁。这意味着要向量化实现上面的计算只需在构造delta2时往矩阵的每一列上分别加上相同的值即可。因此要向量化上面的代码我们只需简单的用repmat命令把sparsity_delta加到delta2的每一列上即可译者注这里原文描述得不是很清楚看似应加到上面代码中delta2行等号右边第一项即W2*delta3上。 中英文对照 向量化 vectorization正向传播 forward propagation反向传播 backpropagation训练样本 training examples激活函数 activation function稀疏自编码网络 sparse autoencoder稀疏惩罚 sparsity penalty平均激活率 average firing rate
