zoukankan      html  css  js  c++  java
  • Deep learning:二十八(使用BP算法思想求解Sparse coding中矩阵范数导数)

      前言:

      关于Sparse coding目标函数的优化会涉及到矩阵求数问题,因为里面有好多矩阵范数的导数,加上自己对矩阵运算不熟悉,推导前面博文Deep learning:二十六(Sparse coding简单理解)中关于拓扑(非拓扑的要简单很多)Sparse coding代价函数对特征变量s导数的公式时,在草稿纸上推导了大半天也没有正确结果。该公式表达式为:

        

      后面继续看UFLDL教程,发现这篇文章Deriving gradients using the backpropagation idea中已经给出了我想要的答案,作者是应用BP神经网络中求网络代价函数导数的思想,将上述代价函数演变成一个多层的神经网络,然后利用每层网络中节点的误差值来反向推导出每一层网络节点的导数。Andrew Ng真值得人佩服,给出的教程切中了我们的要害。

      在看怎样使用BP思想计算矩阵范数的导数时,先看下针对这种问题求解的BP算法形式(和以前经典的BP算法稍有不同,比如说最后一层网络的误差值计算方法,暂时还没弄明白这样更改的理由):

    1. 对网络(由代价函数转换成的网络)中输出层中节点的误差值,采用下面公式计算:

       

      2.  从网络的倒数第2层一直到第2层,依次计算网络每层的误差值:

       

      3.  计算网络中l层的网络参数的偏导(如果是第0层网络,则表示是求代价函数对输入数据作为参数的偏导):

       

      比如在上篇博文中Deep learning:二十七(Sparse coding中关于矩阵的范数求导),就使用过将矩阵范数转换成矩阵的迹形式,然后利用迹的求导公式得出结果,那时候是求sparse coding中非拓扑网络代价函数对权值矩阵A的偏导数,现在用BP思想来求对特征矩阵s的导数,代价函数为:

       

      将表达式中s当做网络的输入,依次将公式中各变量和转换关系变成下面的网络结构:

       

      列出每一层网络的权值,activation函数及其偏导数,误差值,每一层网络的输入,如下所示:

       

      求最后一层网络的误差值时按照前面BP算法的方法此处是:最后一层网络的输出值之和J对最后一层某个节点输入值的偏导,这里的J为:

                                                                                             

      因为此时J对Zi求导是只对其中关于Zi的那一项有效,所以它的偏导数为2*Zi。

      最终代价函数对输入X(这里是s)的偏导按照公式可以直接写出如下:

       

      下面继续来看那个我花了解决一天时间也没推倒出来的偏导数,即在拓扑sparse coding代价函数中关于特征矩阵s的偏导公式。也就是本文一开始给出的公式。

      用同样的方法将其转换成对应的网络结构如下所示:

       

      也同样的,列出它对应网络的参数:

       

       

      其中的输出函数J如下:

      

      最终那个神奇的答案为:

        

      看来这种方法得掌握,如果日后自己论文用到什么公式需要推导的话。

      参考资料:

         Deep learning:二十六(Sparse coding简单理解)

         Deriving gradients using the backpropagation idea

         Deep learning:二十七(Sparse coding中关于矩阵的范数求导)

    作者:tornadomeet 出处:http://www.cnblogs.com/tornadomeet 欢迎转载或分享,但请务必声明文章出处。 (新浪微博:tornadomeet,欢迎交流!)
  • 相关阅读:
    Nginx负载均衡+代理+ssl+压力测试
    Nginx配置文件详解
    HDU ACM 1690 Bus System (SPFA)
    HDU ACM 1224 Free DIY Tour (SPFA)
    HDU ACM 1869 六度分离(Floyd)
    HDU ACM 2066 一个人的旅行
    HDU ACM 3790 最短路径问题
    HDU ACM 1879 继续畅通工程
    HDU ACM 1856 More is better(并查集)
    HDU ACM 1325 / POJ 1308 Is It A Tree?
  • 原文地址:https://www.cnblogs.com/tornadomeet/p/3022361.html
Copyright © 2011-2022 走看看