zoukankan      html  css  js  c++  java
  • 自编码算法与稀疏性

    目前为止,我们已经讨论神经网络在有监督学习中的应用。在有监督学习中,训练样本是有类别标签的。现在假设我们只有一个没有带标签的训练样本集合	extstyle {x^{(1)}, x^{(2)}, x^{(3)}, ldots} ,其中 	extstyle x^{(i)} in Re^{n} 。自编码神经网络是一种无监督学习算法,它使用了反向传播算法,并让目标值等于输入值,比如 	extstyle y^{(i)} = x^{(i)} 。下图是一个自编码神经网络的示例。

    Autoencoder636.png
     
         自编码神经网络尝试学习一个 	extstyle h_{W,b}(x) approx x 的函数。换句话说,它尝试逼近一个恒等函数,从而使得输出 	extstyle hat{x} 接近于输入 	extstyle x恒等函数虽然看上去不太有学习的意义,但是当我们为自编码神经网络加入某些限制,比如限定隐藏神经元的数量,我们就可以从输入数据中发现一些有趣的结构。举例来说,假设某个自编码神经网络的输入x是一张10×10图像(共100个像素)的像素灰度值,于是n=100,其隐藏层	extstyle L_2 中有50个隐藏神经元,注意,输出也是100维的 	extstyle y in Re^{100} 。由于只有50个隐藏神经元,我们迫使神经网络去学习输入数据的压缩表示,也就是说,它必须从50维的隐藏神经元激活度向量	extstyle a^{(2)} in Re^{50}重构出100维的像素灰度值输入 	extstyle x 。如果网络的输入数据是完全随机的,比如每一个输入 	extstyle x_i 都是一个跟其它特征完全无关的独立同分布高斯随机变量,那么这一压缩表示将会非常难学习。但是如果输入数据中隐含着一些特定的结构比如某些输入特征是彼此相关的,那么这一算法就可以发现输入数据中的这些相关性事实上,这一简单的自编码神经网络通常可以学习出一个跟主元分析(PCA)结果非常相似的输入数据的低维表示
         我们刚才论述的是基于隐藏元数量较小的假设。但是即使隐藏神经元的数量较大(可以比输入像素的个数还要多),我们仍然通过给自编码神经网络施加一些其他的限制条件来发现输入数据中的结构。具体来说,如果我们给隐藏神经元加入稀疏性限制,那么自编码神经网络即使在隐藏神经元数量较多的情况下仍然可以发现输入数据中一些有趣的结构
    稀疏性可以被简单地解释如下。如果当神经元的输出接近于1的时候我们认为它被激活,而输出接近于0的时候认为它被抑制,那么使得神经元大部分的时间都被抑制的限制则被称作稀疏性限制。这里我们假设的神经元的激活函数是sigmoid函数。如果你使用tanh作为激活函数的话,当神经元输出为-1的时候,我们认为神经元是被抑制的。
         注意	extstyle a^{(2)}_j表示隐藏神经元j的激活度,但是这一表示方法中并未明确指出哪一个输入x带来了这一激活度。所以我们将使用	extstyle a^{(2)}_j(x)来表示在给定输入为x情况下,自编码神经网络隐藏神经元j激活度。进一步,让
    egin{align}
hat
ho_j = frac{1}{m} sum_{i=1}^m left[ a^{(2)}_j(x^{(i)}) 
ight]
end{align}
    表示隐藏神经元j的平均活跃度(在训练集上取平均)。我们可以近似的加入一条限制
     
    egin{align}
hat
ho_j = 
ho,
end{align}
     
    其中,	extstyle 
ho稀疏性参数,通常是一个接近于0的较小的值(比如	extstyle 
ho = 0.05换句话说,我们想要让隐藏神经元 	extstyle j平均活跃度接近0.05。为了满足这一条件,隐藏神经元的活跃度必须接近于0。

         我们将会在我们的优化目标函数中加入一个额外的惩罚因子,而这一惩罚因子将惩罚那些 	extstyle hat
ho_j 和 	extstyle 
ho 有显著不同的情况从而使得隐藏神经元的平均活跃度保持在较小范围内。惩罚因子的具体形式有很多种合理的选择,我们将会选择以下这一种:

    egin{align}
sum_{j=1}^{s_2} 
ho log frac{
ho}{hat
ho_j} + (1-
ho) log frac{1-
ho}{1-hat
ho_j}.
end{align}

    这里, 	extstyle s_2 是隐藏层中隐藏神经元的数量,而索引 	extstyle j依次代表隐藏层中的每一个神经元。如果你对相对熵(KL divergence)比较熟悉,这一惩罚因子实际上是基于它的。于是惩罚因子也可以被表示为

    egin{align}
sum_{j=1}^{s_2} {
m KL}(
ho || hat
ho_j),
end{align}

    其中 	extstyle {
m KL}(
ho || hat
ho_j)
 = 
ho log frac{
ho}{hat
ho_j} + (1-
ho) log frac{1-
ho}{1-hat
ho_j} 是一个以 	extstyle 
ho为均值和一个以 	extstyle hat
ho_j 为均值的两个伯努利随机变量之间的相对熵。相对熵是一种标准的用来测量两个分布之间差异的方法。

          这一惩罚因子有如下性质,当	extstyle hat
ho_j = 
ho	extstyle {
m KL}(
ho || hat
ho_j) = 0 ,并且随着	extstyle hat
ho_j 与 	extstyle 
ho 之间的差异增大而单调递增。举例来说,在下图中,我们设定 	extstyle 
ho = 0.2并且画出了相对熵值 	extstyle {
m KL}(
ho || hat
ho_j) 随着 	extstyle hat
ho_j 变化的变化。
    KLPenaltyExample.png

    我们可以看出,相对熵在 	extstyle hat
ho_j = 
ho 时达到它的最小值0,而当 	extstyle hat
ho_j 靠近0或者1的时候,相对熵则变得非常大(其实是趋向于	extstyle infty)。所以,最小化这一惩罚因子具有使得 	extstyle hat
ho_j 靠近 	extstyle 
ho 的效果。 现在,我们的总体代价函数可以表示为:

    egin{align}
J_{
m sparse}(W,b) = J(W,b) + eta sum_{j=1}^{s_2} {
m KL}(
ho || hat
ho_j),
end{align}

    其中 	extstyle J(W,b) 如之前所定义,而 	extstyle eta控制稀疏性惩罚因子的权重	extstyle hat
ho_j 项则也(间接地)取决于 	extstyle W,b ,因为它是隐藏神经元 	extstyle j 的平均激活度,而隐藏层神经元的激活度取决于 	extstyle W,b 。

          为了对相对熵进行导数计算,我们可以使用一个易于实现的技巧,这只需要在你的程序中稍作改动即可。具体来说,前面在后向传播算法中计算第二层( 	extstyle l=2更新的时候我们已经计算了
    egin{align}
delta^{(2)}_i = left( sum_{j=1}^{s_{2}} W^{(2)}_{ji} delta^{(3)}_j 
ight) f'(z^{(2)}_i),
end{align}

    现在我们将其换成

    egin{align}
delta^{(2)}_i =
  left( left( sum_{j=1}^{s_{2}} W^{(2)}_{ji} delta^{(3)}_j 
ight)
+ eta left( - frac{
ho}{hat
ho_i} + frac{1-
ho}{1-hat
ho_i} 
ight) 
ight) f'(z^{(2)}_i) .
end{align}

    就可以了。

         有一个需要注意的地方就是我们需要知道	extstyle hat
ho_i 来计算这一项更新。所以在计算任何神经元的后向传播之前,你需要对所有的训练样本计算一遍前向传播,从而获取平均激活度
    如果你的训练样本可以小到被整个存到内存之中(对于编程作业来说,通常如此),你可以方便地在你所有的样本上计算前向传播并将得到的激活度存入内存并且计算平均激活度 。
    然后你就可以使用事先计算好的激活度来对所有的训练样本进行后向传播的计算。
    如果你的数据量太大,无法全部存入内存,你就可以扫过你的训练样本并计算一次前向传播,然后将获得的结果累积起来并计算平均激活度 	extstyle hat
ho_i(当某一个前向传播的结果中的激活度 	extstyle a^{(2)}_i 被用于计算平均激活度 	extstyle hat
ho_i 之后就可以将此结果删除)。然后当你完成平均激活度 	extstyle hat
ho_i 的计算之后,你需要重新对每一个训练样本做一次前向传播从而可以对其进行后向传播的计算。对于后一种情况,你对每一个训练样本需要计算两次前向传播,所以在计算上的效率会稍低一些。
    中英文对照:
    自编码算法 Autoencoders
    稀疏性 Sparsity
    神经网络 neural networks
    监督学习 supervised learning 
    无监督学习 unsupervised learning 
    反向传播算法 backpropagation
    隐藏神经元 hidden units
    像素灰度值 IID
    主元分析 PCA
    抑制 inactive
    激活函数 activation function
    激活度 activation
    平均激活度 the average activation
    稀疏性参数 sparsity parameter
    惩罚因子 penalty term相对熵 KL divergence伯努利随机变量 Bernoulli random variable总体代价函数 overall cost function后向传播 backpropagation前向传播 forward pass梯度下降 gradient descent目标函数 the objective梯度验证方法 the derivative checking method
    转自:
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     

    比如 	extstyle y^{(i)} = x^{(i)} 。下图是一个自编码神经网络的示例。

  • 相关阅读:
    actionbar
    记:使用Xenocode加壳混淆后,无法“自杀覆盖”的自动更新
    用c# 实现设计器 DesignSurface
    将C#程序嵌入资源中(C# 调用嵌入资源的EXE文件方法)
    【转】给自己的软件制作注册码
    免费的Web Service网址
    MS SQL 2005/8 无法远程连接的解决方案【留】
    IIS 部署 ASP.NET 的一个注意点
    成语验证
    C#利用msscript.ocx控件执行JS脚本
  • 原文地址:https://www.cnblogs.com/rong86/p/3519094.html
Copyright © 2011-2022 走看看