如何用Micrograd实现反向传播?完整数学推导与代码解析
Micrograd是一个轻量级的标量自动求导引擎,它以简洁的代码实现了神经网络的核心机制——反向传播。本文将通过数学原理与代码实现的双重解析,帮助你彻底理解反向传播的工作机制,并掌握如何用Micrograd构建自己的神经网络模型。## 反向传播的数学基础:从梯度到链式法则反向传播算法的核心是利用链式法则(Chain Rule)计算损失函数对各参数的梯度。假设我们有一个计算图 \( y =
如何用Micrograd实现反向传播?完整数学推导与代码解析
Micrograd是一个轻量级的标量自动求导引擎,它以简洁的代码实现了神经网络的核心机制——反向传播。本文将通过数学原理与代码实现的双重解析,帮助你彻底理解反向传播的工作机制,并掌握如何用Micrograd构建自己的神经网络模型。
反向传播的数学基础:从梯度到链式法则
反向传播算法的核心是利用链式法则(Chain Rule)计算损失函数对各参数的梯度。假设我们有一个计算图 ( y = f(g(x)) ),则梯度 ( \frac{dy}{dx} = \frac{dy}{dg} \cdot \frac{dg}{dx} )。在神经网络中,这个过程从输出层的损失值开始,逐层反向计算各权重的梯度。
基础运算的梯度推导
Micrograd的Value类实现了基本运算的反向传播逻辑:
- 加法运算:若 ( z = x + y ),则 ( \frac{\partial z}{\partial x} = 1 ),( \frac{\partial z}{\partial y} = 1 )
- 乘法运算:若 ( z = x \cdot y ),则 ( \frac{\partial z}{\partial x} = y ),( \frac{\partial z}{\partial y} = x )
- ReLU激活:若 ( z = \text{ReLU}(x) ),则 ( \frac{\partial z}{\partial x} = 1 )(当 ( x > 0 ))或 ( 0 )(当 ( x \leq 0 ))
Micrograd核心实现:Value类的设计与反向传播
Value类的核心结构
Micrograd的反向传播功能集中在micrograd/engine.py文件的Value类中。每个Value对象不仅存储数据值,还记录梯度和计算图关系:
class Value:
def __init__(self, data, _children=(), _op=''):
self.data = data # 存储标量值
self.grad = 0 # 梯度初始化为0
self._backward = lambda: None # 反向传播函数
self._prev = set(_children) # 子节点(计算图前驱)
self._op = _op # 产生当前节点的操作
反向传播的执行流程
- 拓扑排序:通过
build_topo函数生成计算图的拓扑序列,确保梯度计算顺序正确 - 梯度初始化:将输出节点的梯度设为1(损失对自身的梯度)
- 链式法则应用:按拓扑逆序执行每个节点的
_backward函数,累积梯度
核心代码实现:
def backward(self):
# 拓扑排序所有节点
topo = []
visited = set()
def build_topo(v):
if v not in visited:
visited.add(v)
for child in v._prev:
build_topo(child)
topo.append(v)
build_topo(self)
# 应用链式法则计算梯度
self.grad = 1 # 输出节点梯度初始化为1
for v in reversed(topo):
v._backward() # 执行每个节点的反向传播函数
基础运算的反向传播实现
以加法运算为例,其反向传播函数实现为:
def __add__(self, other):
other = other if isinstance(other, Value) else Value(other)
out = Value(self.data + other.data, (self, other), '+')
def _backward():
self.grad += out.grad # 梯度累加
other.grad += out.grad
out._backward = _backward
return out
乘法运算的反向传播则考虑了输入数据的影响:
def __mul__(self, other):
other = other if isinstance(other, Value) else Value(other)
out = Value(self.data * other.data, (self, other), '*')
def _backward():
self.grad += other.data * out.grad # 梯度 = 另一个输入的值 * 输出梯度
other.grad += self.data * out.grad
out._backward = _backward
return out
神经网络构建:从神经元到多层感知机
Micrograd在micrograd/nn.py中实现了神经网络组件,包括Neuron(神经元)、Layer(层)和MLP(多层感知机):
class Neuron(Module):
def __init__(self, nin):
self.w = [Value(random.uniform(-1,1)) for _ in range(nin)]
self.b = Value(random.uniform(-1,1))
def __call__(self, x):
# 计算 wx + b 并应用ReLU激活
act = sum((wi*xi for wi,xi in zip(self.w, x)), self.b)
return act.relu()
class Layer(Module):
def __init__(self, nin, nout):
self.neurons = [Neuron(nin) for _ in range(nout)]
def __call__(self, x):
return [n(x) for n in self.neurons]
class MLP(Module):
def __init__(self, nin, nouts):
sz = [nin] + nouts
self.layers = [Layer(sz[i], sz[i+1]) for i in range(len(nouts))]
def __call__(self, x):
for layer in self.layers:
x = layer(x)
return x if len(x) > 1 else x[0]
反向传播在神经网络中的应用
当我们用MLP完成前向传播后,只需调用输出值的backward()方法,即可自动计算所有参数的梯度:
# 示例:训练一个简单的MLP
model = MLP(3, [4, 4, 1]) # 输入维度3,隐藏层4-4,输出维度1
x = [Value(1.0), Value(2.0), Value(3.0)]
y_pred = model(x)
loss = (y_pred - y_true) ** 2 # 计算损失
loss.backward() # 自动计算所有参数梯度
实践案例:用Micrograd训练图像分类模型
虽然Micrograd是一个轻量级框架,但它可以用于训练简单的图像分类模型。下面是使用Micrograd构建图像分类器的基本步骤:
- 准备数据集:加载图像数据并转换为数值向量
- 构建模型:定义MLP架构,如
MLP(784, [128, 64, 10])(适用于MNIST) - 训练循环:
- 前向传播:计算预测值
- 计算损失:如均方误差或交叉熵
- 反向传播:调用
loss.backward()计算梯度 - 参数更新:按梯度下降更新权重
Micrograd神经网络训练流程示意图
总结:Micrograd反向传播的核心优势
Micrograd以不到200行代码实现了自动求导和神经网络功能,其设计哲学值得学习:
- 极简设计:用最简洁的代码实现核心功能,易于理解
- 透明性:计算图构建和反向传播过程完全可见,便于调试
- 教育价值:帮助理解深度学习底层原理,是学习反向传播的理想工具
通过本文的解析,你应该已经掌握了Micrograd反向传播的实现原理。想要深入学习,可以查看项目中的demo.ipynb和trace_graph.ipynb文件,通过交互式示例进一步探索自动求导的奥秘。
Micrograd证明,即使是复杂的深度学习算法,也可以通过清晰的数学原理和简洁的代码实现。这正是开源项目的魅力所在——让知识变得透明且可访问。
更多推荐


所有评论(0)