Preface

查看损失函数的时候，突然想手工模拟梯度传递过程。

Script

损失函数是这个:\(l(x)=-||x| - |x-1||\)。程序脚本:

import mxnet as mx;import numpy as np;import mxnet.autograd as ag
x=mx.nd.zeros((,))+
x.attach_grad()
with ag.record():
    l=-mx.nd.abs( mx.nd.abs(x) - mx.nd.abs(x-1.)   )
    l.backward(mx.nd.ones_like(l))
print x,l,'x.grad:',x.grad

先回顾一下梯度传递的理论基础:如果\(y=f(z(x))\)，其中\(f,z\)都可看作不同layer代表的映射，那么需要求得对\(x\)的梯度，易于推广的方法是链式法则:
\[y_x = f_z(z)\times z_x(x)\]
得到普遍应用的原因是易于模块化，每个layer只需要把自己那部分设计好，就可以支持构建自定义系统。
明白这点后，模拟过程最好就是列表:
比如现在\(x=2\)

绝对值的导数:\(f(x)=|x| \Rightarrow f^,(x)=\frac{x}{f(x)}\)
填写Input/Output的过程类似于Forward过程；计算Grad类比于Backward，需要从表格底部算起。
x的节点被复用一次:\(grad_x=1+(-1)=0\)

Note

在\(f(0)=|0|\)处，MXNet使用的应该是置零操作。参考\(x=0\Rightarrow 1, x=0.5\Rightarrow 0, x=1\Rightarrow -1\)