梯度下降法公式推导过程--再次补充:导数部分化简


前面一篇就是基础性的推导过程。从反馈的情况看,总体还是讲明白了。但是在导数的部分,仍有不少的存疑。
其实在数学方面,我也是学渣。所以尽我所能,希望再次的补充能讲的明白。若有谬误,期盼指正。

基础公式

所需基础公式抄录于下,不明白的请至上篇查看详解。

假设函数

$$ y' = h_θ(x) = sum_{i=0}^nθ_ix_i $$

#### 均方差损失函数

$$ J(θ) = frac1{2m}sum_{i=1}^m(h_θ(x^{(i)}) - y^{(i)})^2 $$

#### 梯度下降求解θ

$$ θ_j := θ_j - αfrac∂{∂θ_j}J(θ) $$ 摘出来上面公式步长α之后的部分: $$ egin{align} frac∂{∂θ_j}J(θ) & = frac∂{∂θ_j}frac1{2m}sum_{i=1}^m(h_θ(x^{(i)}) - y^{(i)})^2 \ & = frac1msum_{i=1}^m(h_θ(x^{(i)}) - y^{(i)})x_j^i end{align} $$

嗯,问题一般就是出在这里了,很多人尝试了化简,得不到上面的化简结果。

导数公式

化简上面的式子,需要微积分导数的一些知识,我抄录用到的部分于此,以方便对照查看:

导数

导数的目的是求得在给定点的切线方向,以保证梯度下降的下一步会向收敛方向(也即上面的损失函数最小化方向)迭代一个步长α。这个很多教程都讲过了,这里不再废话。

(偷懒从网上搜了张图,侵删。图中的W实际是我们公式中的θ,J(W)就是我们讲的J(θ))

首先公式(frac∂{∂θ_j})就是求导数的意思,别当做普通的分式,直接分子、分母把∂化简掉成为(frac1{θ_j})。当然大多数人不会这样做了,我只是见过这样的情况,说出来以防万一。

事实上,你把(frac∂{∂θ_j})换成常用的函数描述(f(θ_j))可能更贴切。

#### 对函数的和求导法则 为了描述起来方便,我们下面使用'符号来代表求导:

$$ (u + v)' = u' + v' $$ 在上面的公式中推广一下,Sigma求和不影响求导的传导,直接把Sigma符号提到前面就好: $$ (sum_{i=1}^mu^{(i)})' = sum_{i=1}^m(u^{(i)})' $$

#### 对函数的积求导法则

$$ (uv)' = u'v+uv' $$

#### 幂函数求导法则

$$ (x^u)' = ux^{(u-1)} $$

#### 对常数求导 这是我最爱的部分:

$$ (C)' = 0 $$

#### 链式法则 这是我最不喜欢的部分: 假设我们希望对变量z求导,而变量z依赖变量y,变量y又依赖变量x。例如:

$$ z = f(y) \ y = g(x) $$ 也即: $$ z = f(g(x)) $$ 那么对z求导就构成了链式法则: $$ (z)' = (f(g(x)))'·(g(x))' $$ 注意最后面乘上内部依赖函数求导的过程,简直是反人类的天外来客,经常会忘。但我等遵循自然界规则的凡人又能如何,死记而已。

推导

基本公式列完,开始推导过程:

$$ frac∂{∂θ_j}J(θ) = frac∂{∂θ_j}frac1{2m}sum_{i=1}^m(h_θ(x^{(i)}) - y^{(i)})^2 $$ 根据上面说的求和函数求导法则: $$ = frac1{2m}sum_{i=1}^m(frac∂{∂θ_j}(h_θ(x^{(i)}) - y^{(i)})^2) $$ 别急着对幂求导,考虑对中间的损失函数的依赖,实际要先处理链式法则: $$ = frac1{2m}sum_{i=1}^m(frac∂{∂θ_j}(h_θ(x^{(i)}) - y^{(i)})^2)·(frac∂{∂θ_j}(h_θ(x^{(i)}) - y^{(i)}) $$ 现在方程式前面的部分可以幂求导了,后面的部分把假设函数先展开: $$ = frac1{2m}sum_{i=1}^m2·(h_θ(x^{(i)}) - y^{(i)}))·(frac∂{∂θ_i}(sum_{i=0}^nθ_ix_i - y^{(i)})) $$ 因为展开的假设函数中使用i代表第i个权重,所以前面的求导也换成了(θ_i),不是指第i个批次的样本数据。这里原来没有打算展开讲,所以使用的符号名称有点容易混,但概念清楚的话不应当闹误会。

继续,式子前半部分的2跟1/2会抵消掉,这是前篇做均方差时候乘1/2的目的;后面的Sigma求导继续使用求和函数求导法则展开:

$$ = frac1{m}sum_{i=1}^m(h_θ(x^{(i)}) - y^{(i)}))·(sum_{i=0}^nfrac∂{∂θ_i}θ_ix_i - frac∂{∂θ_i}y^{(i)}) $$ 前半部分的化简已经完成,简单起见,我们只把后面部分摘出来: $$ sum_{i=0}^nfrac∂{∂θ_i}θ_ix_i - frac∂{∂θ_i}y^{(i)}\ = frac∂{∂θ_i}(θ_0x_0+θ_1x_1+...+θ_ix_i+...+θ_nx_n) - frac∂{∂θ_i}y^{(i)} $$ 根据求和函数求导法则展开,等于对其中每一项求导。而我们在对(θ_i)进行求导的时候,其余各项对我们来说,实际上就是一个常数,它们在求导这一刻是固定不能变的。嗯嗯,记得上一篇最后的提醒吗?θ在每个循环内固定不变,在计算完所有的θ之后,才一次代入,并在下个循环内保持不变。

而对常数求导,刚才说过了,那是我的最爱,因为结果是0。还有我们抄了好几行的(y^{(i)})求导,我忍得好辛苦,因为那也是样本集给出的常数,所以结果也是0: $$ = 0 + 0 + ... + frac∂{∂θ_i}θ_ix_i + ... + 0 - 0 $$ 现在需要对乘积函数求导展开了: $$ = frac∂{∂θ_i}θ_i·x_i + θ_i·frac∂{∂θ_i}x_i $$ 你看,这世界不总是那么残酷的,后面的(x_i)又双叒叕是一个常量,所以求导之后乘上(θ_i)仍然是0。 前面对(θ_i)的求导结果是1,原因很简单,你可以把(θ_i)看做1次幂。 $$ egin{align} & = frac∂{∂θ_i}(θ_i)^1·x_i + 0 \ & = 1·θ_i^{(1-1)}·x_i \ & = 1·1·x_i \ & = x_i \ end{align} $$ 只是瞬间,这个世界就清净了。原来对假设函数求导的最终结果,不过是(θ_i)的系数(x_i)。

前面我们两次把等式的局部摘出来化简,现在是把它们组合回去的时候了:

$$ egin{align} θ_j & = θ_j - αfrac∂{∂θ_j}J(θ) \ & = αfrac∂{∂θ_j}frac1{2m}sum_{i=1}^m(h_θ(x^{(i)}) - y^{(i)})^2 \ & = αfrac1msum_{i=1}^m(h_θ(x^{(i)}) - y^{(i)})x_j^i end{align} $$

希望不用再写补充的补充的补充了吧。
原文地址:https://www.cnblogs.com/andrewwang/p/11075395.html