深度学习与智能故障诊断学习笔记（三）——RNN与LSTM推导详解_深度学习故障诊断-程序员宅基地

1.RNN

1.1网络结构

标准神经网络的输入输出在不同例子中可能有不同的长度，在学习中并不共享从不同位置上学到的特征。因为标准神经网络的训练集是稳定的，即所有的特征域表达的内容是同一性质的，一旦交换位置，就需要重新学习。故障诊断和健康管理属于带有时间序列的任务场景，在进行学习时参数量巨大，标准神经网络无法体现出时序上的前因后果，所以引入循环神经网络。如图所示为RNN循环神经网络的单元。

其中 $x^{<t>}$ 为当前输入， $a^{<t-1>}$ 为前一个状态，b为偏置项，tanh为激活函数，用于学习非线性部分。当前输入和前一个状态分别乘以对应权重并相加，在加上偏置项，乘激活函数得到当前状态 $a^{<t>}$ ，此状态在下一个神经元学习时又作为 $a^{<t-1>}$ 进行运算，由此实现时序关联。输出 $\widehat{y}^{t}$ 的激活函数根据任务类型来选择，若是多分类可以选择softmax，若二分类则可直接选择sigmod。

（注：输出并非每个神经元都必须有，RNN可以是多输入多输出，也可以是多输入单输出，仅在学习完成后输出）

1.2RNN网络特点

RNN网络为串联结构，可以体现出“前因后果”，后面结果的生成要参考前面的信息，且所有特征共享一套参数。这使得RNN在面对不同的输入（两个方面），可以学习到不同的相应结果，并极大的减少了训练参数量。

（RNN输入和输出数据在不同场景中可以有不同的长度）

1.3损失函数

单个时间步的损失函数可根据多分类和二分类进行自定义

整个序列的损失函数是将所有单步损失函数相加，如式。

1.4传播过程

前向传播如图一所示。

反向传播

（图源自吴恩达老师课件）

求解梯度即复合函数求导，按照链式法则进行求导。

反向传播具体过程需要按照损失函数来具体求解，但上式对所有RNN模型都适用。

1.5缺点

当序列太长时，容易产生梯度消失，参数更新只能捕捉到局部以来关系，没法再捕捉序列之间长期的关联或依赖关系。

如图为RNN连接，输入x，输出o（简单线性输出），权重w，s为生成状态。

根据前向传播可得：

假设使用平方误差作为损失函数，对单个时间点进行求梯度，假设再t=3时刻，损失函数为L3 = $\frac{1}{2}\left ( {Y}_3 - {O}_3 \right )^{2}$ 。然后根据网络参数Wx，Ws，Wo，b1，b2等求梯度。

Wo:

Wx（具体求解过程在下边）:

经整理可得：

具体求解过程：

首先，所求目标为L3对Wx的偏导，通过链式法则进行展开。对比前向传播公式图可知，O3中并不能直接对Wx求偏导，而是包含在S3中，所以要展开成如下形式。

但在S3中又包含S2，S2中包含Wx和S1，S1中又包含Wx，嵌套了很多层，为了方便表示，我们用 $\theta$ 3来表示S3括号中的内容。进一步简化可得：

由S3演变为S2，同理可递推求出 $\frac{\partial S_{2}}{\partial W_{x}}$ 和 $\frac{\partial S_{1}}{\partial W_{x}}$

再将所求出结果回代到公式中，可以得出 $\frac{\partial S_{3}}{\partial W_{x}}$

再回带至 $\frac{\partial L_{3}}{\partial W_{x}}$

由该式可以看出，梯度的更新同时依赖于x3，x2,x1包括其梯度值。将该式处理为

此为t=3时刻的梯度公式，推广至任意时刻的梯度公式为：

此式括号中的项为求导的连乘，此处求出的导数是介于0-1之间的，有一定的机率导致梯度消失（但非主要原因）。造成梯度消失和梯度爆炸的主要原因是最后一项：当Ws很小的时候，它的k-1的次方会无限接近于0，而当Ws大于1时，它的k-1次方会很大。

如下为t=20时梯度更新计算的结果：

从式中可以看出，t=3的节点由于连乘过多导致梯度消失，无法将信息传给t=20，因此t=20的更新无法引入t=3时的信息，认为t=20节点跟t=3的节点无关联。

对于梯度爆炸和梯度消失，可以通过梯度修剪来解决。相对于梯度爆炸，梯度消失更难解决。而LSTM很好的解决了这些问题。

2.LSTM

2.1设计思路

RNN是想把所有信息都记住，不管是有用的信息还是没用的信息。而LSTM设计了一个记忆细胞，具备选择性记忆功能，可以选择记忆重要信息，过滤掉噪声信息，减轻记忆负担。

2.2整体结构

如图为LSTM与RNN结构对比

LSTM

RNN

2.3单元结构

在LSTM每个时间步中，都有一个记忆细胞，这个东西给予了LSTM选择记忆功能，使得LSTM有能力自由选择每个时间步里面记忆的内容。

下图中Ct-1为上一个记忆细胞，ht-1为上一个时间点的状态，经过该单元，输出一个新的记忆细胞和一个新的状态。在单元中有三个 $\sigma$ ， $\sigma$ 被称为门单元，它的输出值介于0-1之间。ft为遗忘门，it为更新门，Ot为输出门。

门是一种选择性地让信息通过的方法。它们由sigmoid神经网络层和逐点乘法运算组成。Sigmoid层输出0到1之间的数字，描述每个组件应允许通过多少。值为零表示“不让任何内容通过”，而值为 1 表示“允许所有信息通过”

公式所示为前向传播，Ct与Ct-1，ht-1，xt等参数都有关，其中Wxf，Whf，Wxi分别代表相应权重。在单元结构图中可以看出ft与Ct-1进行×运算（对应元素相乘），gt与it进行×运算，两者相加为新生成的ct。

2.4 缓解梯度爆炸和梯度消失

此过程为公式推导（以求Wxf为例）。

通过调节Whf，Whi，Whg的值，可以灵活控制Ct对Ct-1的偏导值，当要从n时刻长期记忆某个东西到m时刻时，该路径上的

从而大大缓解了梯度消失和梯度爆炸。

（B站搜索老弓的学习日记，本篇博客为RNN与LSTM的学习笔记）

本文链接：https://blog.csdn.net/Vanderbiol/article/details/124256099

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

oracle 12c 集群安装后的检查_12c查看crs状态-程序员宅基地

文章浏览阅读1.6k次。安装配置gi、安装数据库软件、dbca建库见下：http://blog.csdn.net/kadwf123/article/details/784299611、检查集群节点及状态：[root@rac2 ~]# olsnodes -srac1 Activerac2 Activerac3 Activerac4 Active[root@rac2 ~]_12c查看crs状态

解决jupyter notebook无法找到虚拟环境的问题_jupyter没有pytorch环境-程序员宅基地

文章浏览阅读1.3w次，点赞45次，收藏99次。我个人用的是anaconda3的一个python集成环境，自带jupyter notebook，但在我打开jupyter notebook界面后，却找不到对应的虚拟环境，原来是jupyter notebook只是通用于下载anaconda时自带的环境，其他环境要想使用必须手动下载一些库：1.首先进入到自己创建的虚拟环境(pytorch是虚拟环境的名字)activate pytorch2.在该环境下下载这个库conda install ipykernelconda install nb__jupyter没有pytorch环境

国内安装scoop的保姆教程_scoop-cn-程序员宅基地

文章浏览阅读5.2k次，点赞19次，收藏28次。选择scoop纯属意外，也是无奈，因为电脑用户被锁了管理员权限，所有exe安装程序都无法安装，只可以用绿色软件，最后被我发现scoop，省去了到处下载XXX绿色版的烦恼，当然scoop里需要管理员权限的软件也跟我无缘了（譬如everything）。推荐添加dorado这个bucket镜像，里面很多中文软件，但是部分国外的软件下载地址在github，可能无法下载。以上两个是官方bucket的国内镜像，所有软件建议优先从这里下载。上面可以看到很多bucket以及软件数。如果官网登陆不了可以试一下以下方式。_scoop-cn

Element ui colorpicker在Vue中的使用_vue el-color-picker-程序员宅基地

文章浏览阅读4.5k次，点赞2次，收藏3次。首先要有一个color-picker组件 <el-color-picker v-model="headcolor"></el-color-picker>在data里面data() { return {headcolor: ’ #278add ’ //这里可以选择一个默认的颜色} }然后在你想要改变颜色的地方用v-bind绑定就好了，例如：这里的:sty..._vue el-color-picker

迅为iTOP-4412精英版之烧写内核移植后的镜像_exynos 4412 刷机-程序员宅基地

文章浏览阅读640次。基于芯片日益增长的问题，所以内核开发者们引入了新的方法，就是在内核中只保留函数，而数据则不包含，由用户（应用程序员）自己把数据按照规定的格式编写，并放在约定的地方，为了不占用过多的内存，还要求数据以根精简的方式编写。boot启动时，传参给内核，告诉内核设备树文件和kernel的位置，内核启动时根据地址去找到设备树文件，再利用专用的编译器去反编译dtb文件，将dtb还原成数据结构，以供驱动的函数去调用。firmware是三星的一个固件的设备信息，因为找不到固件，所以内核启动不成功。_exynos 4412 刷机

Linux系统配置jdk_linux配置jdk-程序员宅基地

文章浏览阅读2w次，点赞24次，收藏42次。Linux系统配置jdkLinux学习教程，Linux入门教程（超详细）_linux配置jdk

随便推点

matlab(4)：特殊符号的输入_matlab微米怎么输入-程序员宅基地

文章浏览阅读3.3k次，点赞5次，收藏19次。xlabel('\delta');ylabel('AUC');具体符号的对照表参照下图：_matlab微米怎么输入

C语言程序设计-文件(打开与关闭、顺序、二进制读写)-程序员宅基地

文章浏览阅读119次。顺序读写指的是按照文件中数据的顺序进行读取或写入。对于文本文件，可以使用fgets、fputs、fscanf、fprintf等函数进行顺序读写。在C语言中，对文件的操作通常涉及文件的打开、读写以及关闭。文件的打开使用fopen函数，而关闭则使用fclose函数。在C语言中，可以使用fread和fwrite函数进行二进制读写。‍ Biaoge 于2024-03-09 23:51发布阅读量：7 ️文章类型：【 C语言程序设计】在C语言中，用于打开文件的函数是____，用于关闭文件的函数是____。

Touchdesigner自学笔记之三_touchdesigner怎么让一个模型跟着鼠标移动-程序员宅基地

文章浏览阅读3.4k次，点赞2次，收藏13次。跟随鼠标移动的粒子以grid（SOP）为partical（SOP）的资源模板，调整后连接【Geo组合+point spirit（MAT)】，在连接【feedback组合】适当调整。影响粒子动态的节点【metaball(SOP)+force(SOP)】添加mouse in（CHOP)鼠标位置到metaball的坐标，实现鼠标影响。..._touchdesigner怎么让一个模型跟着鼠标移动

【附源码】基于java的校园停车场管理系统的设计与实现61m0e9计算机毕设SSM_基于java技术的停车场管理系统实现与设计-程序员宅基地

文章浏览阅读178次。项目运行环境配置：Jdk1.8 + Tomcat7.0 + Mysql + HBuilderX（Webstorm也行）+ Eclispe（IntelliJ IDEA,Eclispe,MyEclispe,Sts都支持）。项目技术：Springboot + mybatis + Maven +mysql5.7或8.0+html+css+js等等组成，B/S模式 + Maven管理等等。环境需要1.运行环境：最好是java jdk 1.8，我们在这个平台上运行的。其他版本理论上也可以。_基于java技术的停车场管理系统实现与设计

Android系统播放器MediaPlayer源码分析_android多媒体播放源码分析时序图-程序员宅基地

文章浏览阅读3.5k次。前言对于MediaPlayer播放器的源码分析内容相对来说比较多，会从Java-&amp;gt;Jni-&amp;gt;C/C++慢慢分析，后面会慢慢更新。另外，博客只作为自己学习记录的一种方式，对于其他的不过多的评论。MediaPlayerDemopublic class MainActivity extends AppCompatActivity implements SurfaceHolder.Cal..._android多媒体播放源码分析时序图

java 数据结构与算法 ——快速排序法-程序员宅基地

文章浏览阅读2.4k次，点赞41次，收藏13次。java 数据结构与算法 ——快速排序法_快速排序法