用Numpy搭建神经网络第二期:梯度下降法的实现

栏目: Python · 发布时间: 5年前

内容简介:小伙伴们大家好呀~~用Numpy搭建神经网络,我们已经来到第二期了。第一期文摘菌教大家如何用Numpy搭建一个简单的神经网络,完成了前馈部分。这一期,为大家带来梯度下降相关的知识点,和上一期一样,依然用Numpy实现梯度下降。在代码开始之前,先来普及一下梯度下降的知识点吧。

用Numpy搭建神经网络第二期:梯度下降法的实现

大数据文摘出品

作者:蒋宝尚

小伙伴们大家好呀~~用Numpy搭建神经网络,我们已经来到第二期了。第一期文摘菌教大家如何用Numpy搭建一个简单的神经网络,完成了前馈部分。

这一期,为大家带来梯度下降相关的知识点,和上一期一样,依然用Numpy实现梯度下降。在代码开始之前,先来普及一下梯度下降的知识点吧。

梯度下降:迭代求解模型参数值

第一期文章中提到过,最简单的神经网络包含三个要素,输入层,隐藏层以及输出层。关于其工作机理其完全可以类比成一个元函数:Y=W*X+b。即输入数据X,得到输出Y。

如何评估一个函数的好坏,专业一点就是拟合度怎么样?最简单的方法是衡量真实值和输出值之间的差距,两者的差距约小代表函数的表达能力越强。

这个差距的衡量也叫损失函数。显然,损失函数取值越小,原函数表达能力越强。

那么参数取何值时函数有最小值?一般求导能够得到局部最小值(在极值点处取)。而梯度下降就是求函数有最小值的参数的一种方法。

梯度下降数学表达式

比如对于线性回归,假设函数表示为h θ (x 1, x 2 …x n )=θ 01 x 1 +..+θ n x n ,其中w i (i=0,1,2...n)为模型参数,x i (i=0,1,2...n)为每个样本的n个特征值。这个表示可以简化,我们增加一个特征x 0 =1,这样h(x o ,x 1 ,.…x n )=θ 0 x 01 x 1 +..+θ n x n 。同样是线性回归,对应于上面的假设函数,损失函数为(此处在损失函数之前加上1/2m,主要是为了修正SSE让计算公式结果更加美观,实际上损失函数取MSE或SSE均可,二者对于一个给定样本而言只相差一个固定数值):

用Numpy搭建神经网络第二期:梯度下降法的实现

算法相关参数初始化:主要是初始化θ 0 ,θ 1 ..,θ n ,我们比较倾向于将所有的初始化为0,将步长初始化为1。在调优的时候再进行优化。

对θ i 的梯度表达公式如下:

用Numpy搭建神经网络第二期:梯度下降法的实现  

用步长(学习率)乘以损失函数的梯度,得到当前位置下降的距离,即:

用Numpy搭建神经网络第二期:梯度下降法的实现

梯度下降法的矩阵方式描述

对应上面的线性函数,其矩阵表达式为:

用Numpy搭建神经网络第二期:梯度下降法的实现

损失函数表达式为:

用Numpy搭建神经网络第二期:梯度下降法的实现  

其中Y为样本的输出向量。

梯度表达公式为:

用Numpy搭建神经网络第二期:梯度下降法的实现

还是用线性回归的例子来描述具体的算法过程。损失函数对于向量的偏导数计算如下:

用Numpy搭建神经网络第二期:梯度下降法的实现

迭代:

用Numpy搭建神经网络第二期:梯度下降法的实现

两个矩阵求导公式为:

用Numpy搭建神经网络第二期:梯度下降法的实现

Python 实现梯度下降

<span style="">import pandas as pd</span>

<span>import numpy as np</span>

导入两个必要的包。

<span><span>def</span> <span>regularize(xMat):</span></span>

<span> <span>inMat</span>=<span>xMat. copy()</span></span>

<span> <span>inMeans</span>=<span>np. mean(inMat, axis=0)</span></span>

<span> <span>invar</span>=<span>np. std(inMat, axis=0)</span></span>

<span> <span>inMat</span>=<span>(inMat-inMeans)/invar </span></span>

<span> <span>return</span> <span>inMat</span></span>

定义标准化函数,不让过大或者过小的数值影响求解。

定义梯度下降函数:


 

def BGD_LR(data alpha=0.001, maxcycles=500):

xMat=np. mat(dataset)

yMat=np. mat(dataset).T

xMat=regularize(xMat)

m,n=xMat.shape

weights=np. zeros((n,1))

for i in range(maxcycles):

grad=xMat.T*(xMat * weights-yMat)/m

weights=weights -alpha* grad

return weights

其中,dataset代表输入的数据,alpha是学习率,maxcycles是最大的迭代次数。

即返回的权重就是说求值。np.zeros 是初始化函数。grad的求取是根据梯度下降的矩阵求解公式。

本文参考B站博主 菊安酱的机器学习。感兴趣的同学可以打开链接观看 视频哟~

https://www.bilibili.com/video/av35390140

好了,梯度下降这个小知识点就讲解完了,下一期,我们将第一期与第二期的知识点结合,用手写数字的数据完成一次神经网络的训练。

实习/全职编辑记者招聘ing

加入我们,亲身体验一家专业科技媒体采写的每个细节,在最有前景的行业,和一群遍布全球最优秀的人一起成长。坐标北京·清华东门,在大数据文摘主页对话页回复 “招聘” 了解详情。简历请直接发送至zz@bigdatadigest.cn

用Numpy搭建神经网络第二期:梯度下降法的实现

点「在看」的人都变好看了


以上就是本文的全部内容,希望本文的内容对大家的学习或者工作能带来一定的帮助,也希望大家多多支持 码农网

查看所有标签

猜你喜欢:

本站部分资源来源于网络,本站转载出于传递更多信息之目的,版权归原作者或者来源机构所有,如转载稿涉及版权问题,请联系我们

社群营销实战手册

社群营销实战手册

秋叶、邻三月、秦阳 / 人民邮电出版社 / 2018-1 / 69.00元

互联网正从“物以类聚”,走向“人以群分”的时代。秋叶等人的“社群营销”,并非单纯靠社群卖东西,而是建立一种中心化的、自行运转的生态,让“同好”们形成紧密的联系,创造出海量营销机会。 《社群营销实战手册 从社群运营到社群经济》共5章内容,从社群的定位、建立、扩张、变现、运营,到社群的生命周期延长、社群运营团队的打造和管理以及社群管理工具,大量干货秘笈一应俱全,并提供丰富的运营实战案例,全面解读社群的......一起来看看 《社群营销实战手册》 这本书的介绍吧!

图片转BASE64编码
图片转BASE64编码

在线图片转Base64编码工具

URL 编码/解码
URL 编码/解码

URL 编码/解码

UNIX 时间戳转换
UNIX 时间戳转换

UNIX 时间戳转换