预备知识5-自动微分学习笔记

1年前浏览1255

深度学习框架通过自动计算导数，即自动微分（automatic differentiation）来加快求导。实际中，根据设计好的模型，系统会构建一个计算图（computational graph），来跟踪计算是哪些数据通过哪些操作组合起来产生输出。自动微分使系统能够随后反向传播梯度。这里，反向传播（backpropagate）意味着跟踪整个计算图，填充关于每个参数的偏导数。计算图是深度学习框架实现自动微分的核心机制：

前向传播构建计算图，记录数据流动路径。
反向传播按图回溯，通过链式法则高效计算梯度。

一、一个简单的例子

作为一个例子，假设我们想对函数y = 2x⊤x关于列向量x求导。注意: 一个标量函数关于向量x的梯度是向量，并且与x具有相同的形状。

import torchx = torch.arange(4.0) # tensor([0., 1., 2., 3.])# requires_grad_(True) 是一个用于启用张量（Tensor）梯度跟踪的方法。它的作用是告诉# PyTorch 需要记录该张量参与的所有操作，以便后续通过自动微分（Autograd）计算梯度。x.requires_grad_(True) # # 等价于x=torch.arange(4.0,requires_grad=True)# 函数一：y = 2*dot(x, x)y = 2*torch.dot(x, x)print(y) # tensor(28., grad_fn=<MulBackward0>)# grad_fn=<DotBackward0>是张量（Tensor）的一个属性，表示该张量是通过某个操作(这里是# 矩阵乘法或向量点积)生成的，并且系统已记录该操作的梯度计算规则，用于后续的反向传播。# 通过调用反向传播函数(backward)来自动计算y关于x每个分量的梯度y.backward()print(x.grad) # tensor([ 0.,  4.,  8., 12.])# 函数二: y = sum(x)# 在默认情况下，PyTorch会积累梯度，我们需要清除之前的值x.grad.zero_()y = x.sum()y.backward()print(x.grad) # tensor([1., 1., 1., 1.])

二、非标量变量的反向传播

对非标量调用backward需要传入一个gradient参数，该参数指定微分函数关于self的梯度。

gradient 参数的作用

将非标量输出的梯度加权求和为标量，再反向传播。
数学意义：若 y = f(x) 是非标量，则 gradient 是 ∂L/∂y（损失函数 L 对 y 的梯度），PyTorch 据此计算 ∂L/∂x = (∂y/∂x)^T · (∂L/∂y)。

非标量直接反向传播（错误）

import torchx = torch.tensor([1.0, 2.0], requires_grad=True)y = x ** 2  # y = [1.0, 4.0]# 报错：RuntimeError: grad can be implicitly created only for scalar outputsy.backward()

正确用法：传入 `gradient` 参数

# 定义梯度权重（形状需与 y 一致）gradient = torch.tensor([0.1, 0.01])  # 假设 ∂L/∂y = [0.1, 0.01]y.backward(gradient) # 等价于先求和：L = 0.1*y1 + 0.01*y2，再反向传播
print(x.grad) # 输出: tensor([0.2000, 0.0400])# 计算过程: ∂L/∂x1 = 0.1*2*x1 = 0.2, ∂L/∂x2 = 0.01*2*x2 = 0.04

实际应用：Jacobian 矩阵计算

若需计算完整 Jacobian 矩阵（而非加权求和），可对每个输出分量单独反向传播：

x = torch.tensor([1.0, 2.0], requires_grad=True)y = x ** 2                 # y = [y1, y2]
# 计算 ∂y1/∂x 和 ∂y2/∂xjacobian = torch.zeros(2, 2)for i in range(2):    y.backward(torch.tensor([1.0 if j == i else 0.0 for j in range(2)]), retain_graph=True)    jacobian[i] = x.grad    x.grad.zero_()
print(jacobian)# 输出: tensor([[2., 0.], [0., 4.]]) （即 [[∂y1/∂x1, ∂y1/∂x2], [∂y2/∂x1, ∂y2/∂x2]]）

三、分离计算

在 PyTorch 中，分离计算（通过.detach() 或 torch.no_grad()）是一种关键机制，用于控制梯度传播和优化内存使用。

假设y是x的函数，而z是y和x的函数。现需计算z关于x的梯度，但希望将y视为一个常数，且只考虑到x在y被计算后发挥的作用。

import torchx = torch.arange(4.0, requires_grad=True) # tensor([0., 1., 2., 3.])# 本例只想求偏导数的和，所以传递一个1的梯度是合适的y = x * x # tensor([0., 1., 4., 9.], grad_fn=<MulBackward0>)# 等价于y.backward(torch.ones(len(x)))y.sum().backward()print(x.grad) # tensor([0., 2., 4., 6.])
x.grad.zero_()y=x*xu=y.detach()z = u*xz.backward(torch.ones_like(x))print(x.grad) # tensor([0., 1., 4., 9.])
# 由于记录了y的计算结果，我们可以随后在y上调用反向传播，# 得到y=x*x关于的x的导数，即2*x。x.grad.zero_()y.sum().backward()pprint(x.grad == 2 * x) # tensor([True, True, True, True])

四、Python控制流的梯度计算

使用自动微分的一个好处是：即使构建函数的计算图需要通过Python控制流（例如，条件、循环或任意函数调用），我们仍然可以计算得到的变量的梯度。在下面的代码中，while循环的迭代次数和if语句的结果都取决于输入a的值。

import torchdef f(a):  b = a * 2  while b.norm() < 1000:    b = b * 2  if b.sum() > 0:    c = b  else:    c = 100 * b  return c
#a = torch.randn(size=(), requires_grad=True)a = torch.tensor(-1, dtype=torch.float32, requires_grad=True)print(a)d = f(a)d.backward()print(a.grad) # tensor(102400.)

来源：檐苔

线性神经网络4-softmax回归

一、分类问题标签表示方法：以不同的整数表示不同的类别，如以y∈{1,2,3}分别代表{狗,猫,鸡}；独热编码(one‐hotencoding)：独热编码是一个向量，它的分量和类别一样多。类别对应的分量设置为1，其他所有分量设置为0。如标签定义y为三维向量，其中(1,0,0)对应于“猫”、(0,1,0)对应于“鸡”、(0,0,1)对应于“狗”。二、网络架构仿射函数：数学和机器学习中的基础概念，它描述了输入与输出之间的线性关系，同时允许存在平移（截距项）。为了解决线性模型的分类问题，定义与输出一样多的仿射函数，即每个输出对应一个仿射函数。如有4个特征和3个可能的输出类别，利用仿射函数可表达为：可简化为向量形式表达式：o=Wx+b。与线性回归一样，softmax回归也是一个单层神经网络。由于计算每个输出o1、o2和o3取决于所有输入x1、x2、x3和x4，所以softmax回归的输出层也是全连接层。三、softmax运算softmax函数：对每个未规范化的预测求幂，再让每个求幂后的结果除以它们的总和，即：从以上定义可以看出，softmax函数能够将未规范化的预测变换为非负数并且总和为1，同时让模型保持可导的性质。softmax运算不会改变未规范化的预测o之间的大小次序，只会确定分配给每个类别的概率。尽管softmax是一个非线性函数，但softmax回归的输出仍然由输入特征的仿射变换决定。因此，softmax回归是一个线性模型(linearmodel)。softmax函数的导数及推导：四、交叉熵损失函数Cross-entropy（交叉熵损失函数)交叉熵是用来评估当前训练得到的概率分布与真实分布的差异情况。它刻画的是实际输出(概率)与期望输出(概率)的距离，也就是交叉熵的值越小，两个概率分布就越接近。常作为分类问题的损失函数。具体公式如下：其中，y_j是一个表示真实标签的长度为q的独热编码向量，即除了第j项为1其余项均为0；ŷ_j表示第j项的预测的概率值。交叉熵损失函数的导数及推导：五、分类问题采用交叉熵而不是均方差的为损失函数的解释交叉熵损失在分类任务中因其梯度高效、概率解释性强和优化稳定性，成为比MSE更优的选择。六、softmax函数上溢的解决方案在softmax计算之前，先从所有x_k中减去max(x_k):在减法和规范化步骤之后，可能有些x_i−max(x)具有较大的负值，exp(x_i−max(x))将有接近零的值而出现下溢(underflow)，进而导致ˆy_i为零，并且使得log(ˆy_i)的值为-inf。反向传播几步后，我们可能会发现自己面对一屏幕可怕的nan结果。尽管我们要计算指数函数，但我们最终在计算交叉熵损失时会取它们的对数。通过将softmax和交叉熵结合在一起，可以避免反向传播过程中可能会困扰我们的数值稳定性问题。七、softmax分类的实现importtorchfromtorch.utilsimportdataimporttorchvisionfromtorchvisionimporttransformsimportmatplotlib.pyplotaspltimporttypingfrommatplotlib_inlineimportbackend_inlinefromIPythonimportdisplayimportplotShow#定义一个累加器，实现对n个变量进行累加classAccumulator:'''在n个变量上累加'''def__init__(self,n:int):self.data=[0.0]*ndefadd(self,*args):self.data=[a+float(b)fora,binzip(self.data,args)]defreset(self):self.data=[0.0]*len(self.data)def__getitem__(self,idx):returnself.data[idx]classAnimator:'''在动画中绘制数据'''def__init__(self,xlabel=None,ylabel=None,legend=None,xlim=None,ylim=None,xscale='linear',yscale='linear',fmts=('-','m--','g-.','r:'),nrows=1,ncols=1,figsize=(3.5,2.5)):#增量地绘制多条线条iflegendisNone:legend=[]backend_inline.set_matplotlib_formats('svg')self.fig,self.axes=plt.subplots(nrows,ncols,figsize=figsize)ifnrows*ncols==1:self.axes=[self.axes,]#使用lambda函数捕获参数self.configAxes=lambda:plotShow.setAxes(self.axes[0],xlabel,ylabel,xlim,ylim,xscale,yscale,legend)self.X,self.Y,self.fmts=None,None,fmtsdefadd(self,x,y):#向图表中添加多个数据点ifnothasattr(y,'__len__'):y=[y]n=len(y)ifnothasattr(x,'__len__'):x=[x]*nifnotself.X:self.X=[[]for_inrange(n)]ifnotself.Y:self.Y=[[]for_inrange(n)]fori,(a,b)inenumerate(zip(x,y)):ifaisnotNoneandbisnotNone:self.X[i].append(a)self.Y[i].append(b)self.axes[0].cla()forx,y,fmtinzip(self.X,self.Y,self.fmts):self.axes[0].plot(x,y,fmt)self.configAxes()display.display(self.fig)plt.draw()plt.pause(0.001)display.clear_output(wait=True)defshow(self):display.display(self.fig)defgetFashionMnistLabels(labels:typing.Sequence):'''返回FashionMnist数据集的文本标签'''textLabels=['t-shirt','trouser','pullever','dress','coat','sandal','shirt','sneaker','bag','ankleboot']return[textLabels[int(i)]foriinlabels]defshowImages(imgs:list,numRows:int,numCols:int,titles:list=None,scale:float=1.5):figsize=(numCols*scale,numRows*scale)_,axes=plt.subplots(numRows,numCols,figsize=figsize)axes=axes.flatten()fori,(ax,img)inenumerate(zip(axes,imgs)):iftorch.is_tensor(img):#图片张量ax.imshow(img.numpy())else:#PIL图片ax.imshow(img)ax.axes.get_xaxis().set_visible(False)ax.axes.get_yaxis().set_visible(False)iftitles:ax.set_title(titles[i])returnaxesdefgetDataloaderWorkers()->int:'''使用4个进程来读取数据'''return4#读取数据集defloadDataFashionMnist(batchSize:int,resize=None)->tuple:'''下载Fashion-MNIST数据集，然后将其加载到内存中'''trans=[transforms.ToTensor()]ifresize:trans.insert(0,transforms.Resize(resize))trans=transforms.Compose(trans)mnistTrain=torchvision.datasets.FashionMNIST(root='../data',train=True,transform=trans,download=True)mnistTest=torchvision.datasets.FashionMNIST(root='../data',train=False,transform=trans,download=True)return(data.DataLoader(mnistTrain,batchSize,shuffle=True,num_workers=getDataloaderWorkers()),data.DataLoader(mnistTest,batchSize,shuffle=False,num_workers=getDataloaderWorkers()))#softmax函数defsoftmax(X:torch.Tensor)->torch.Tensor:'''1.对每个项求幂(使用exp);2.对每一行求和(小批量中每个样本是一行),得到每个样本的规范化常数；3.将每一行除以其规范化常数,确保结果的和为1。'''XExp:torch.Tensor=torch.exp(X)partition:torch.Tensor=XExp.sum(1,keepdim=True)returnXExp/partition#定义分类精度defaccuracy(yHat:torch.Tensor,y:torch.Tensor)->float:'''计算预测正确的数量'''iflen(yHat.shape)>1andyHat.shape[1]>1:yHat=yHat.argmax(axis=1)cmp=yHat.type(y.dtype)==yreturnfloat(cmp.type(y.dtype).sum())defevaluateAccuracy(net:torch.nn.Module,dataIter:typing.Tuple[torch.Tensor])->float:'''计算在指定数据集上模型的精度'''ifisinstance(net,torch.nn.Module):net.eval()#将模型设置为评估模式metric=Accumulator(2)#统计正确预测数、预测总数withtorch.no_grad():forX,yindataIter:metric.add(accuracy(net(X),y),y.numel())returnmetric[0]/metric[1]#定义训练模型deftrainEpochCh3(net:torch.nn.Module,trainIter:typing.Tuple[torch.Tensor],loss:typing.Callable[[torch.Tensor,torch.Tensor],torch.Tensor],updater:torch.optim.Optimizer)->typing.Tuple[float]:'''训练模型一个周期'''#将模型设置为训练模型ifisinstance(net,torch.nn.Module):net.train()#训练损失总和、训练准确度总和、样本数metric=Accumulator(3)forX,yintrainIter:yHat:torch.Tensor=net(X)l:torch.Tensor=loss(yHat,y)ifisinstance(updater,torch.optim.Optimizer):#使用PyTorch内置的优化器和损失函数updater.zero_grad()l.mean().backward()updater.step()else:#使用定制的优化器和损失函数l.sum().backward()updater(X.shape[0])metric.add(float(l.sum()),accuracy(yHat,y),y.numel())#返回训练损失和训练精度returnmetric[0]/metric[2],metric[1]/metric[2]deftrainCh3(net:torch.nn.Module,trainIter:data.DataLoader,testIter:data.DataLoader,loss:typing.Callable[[torch.Tensor,torch.Tensor],torch.Tensor],numEpochs:int,updater:torch.optim.Optimizer)->None:'''训练模型'''animator=Animator(xlabel='epoch',xlim=[1,numEpochs],ylim=[0.3,0.9],legend=['trainloss','trainacc','testacc'])forepochinrange(numEpochs):trainMetrics=trainEpochCh3(net,trainIter,loss,updater)testAcc=evaluateAccuracy(net,testIter)animator.add(epoch+1,trainMetrics+(testAcc,))trainLoss,trainAcc=trainMetricsasserttrainLoss<0.5,trainLossasserttrainAcc<=1andtrainAcc>0.7,trainAccasserttestAcc<=1andtestAcc>0.7,testAcc#定义优化算法defSGD(params:list[torch.Tensor],lr:float,batchSize:int):'''小批量随机梯度下降法params:模型参数集合lr:学习率，确定每一步更新的大小batchSize:批量样本大小'''withtorch.no_grad():forparaminparams:param-=lr*param.grad/batchSizeparam.grad.zero_()#预测defpredictCh3(net,testIter,n=6):'''预测标签'''forX,yintestIter:breaktrues=getFashionMnistLabels(y)preds=getFashionMnistLabels(net(X).argmax(axis=1))titles=[true+'\n'+predfortrue,predinzip(trues,preds)]showImages(X[0:n].reshape((n,28,28)),1,n,titles=titles[0:n])plt.show()#softmax回归的从零开始实现if__name__=='__main__':#1.读取数据集batchSize:int=256trainIter,testIter=loadDataFashionMnist(batchSize)#2.模型参数初始化numInputs:int=784#原始数据集每个样本是28*28的图像，因此以构造一个长度为784的向量来表示每一个像素点numOutPuts:int=10#数据集中有10个类别W=torch.normal(0.0,0.01,size=(numInputs,numOutPuts),requires_grad=True)#权重项b=torch.zeros(numOutPuts,requires_grad=True)#偏置项#3.定义softmax模型defnet(X:torch.Tensor)->torch.Tensor:returnsoftmax(torch.matmul(X.reshape((-1,W.shape[0])),W)+b)#4.定义损失函数defcrossEntropy(yHat:torch.Tensor,y:torch.Tensor)->torch.Tensor:'''yHat:预测概率y:标签向量'''return-torch.log(yHat[range(len(yHat)),y])#5.定义优化函数lr=0.1defupdater(batchSize:int):returnSGD([W,b],lr,batchSize)#6.模型训练numEpochs:int=10trainCh3(net,trainIter,testIter,crossEntropy,numEpochs,updater)#7.预测predictCh3(net,testIter,n=10)#softmax回归的简洁实现if__name__=='__main__':#1.读取数据集batchSize:int=256trainIter,testIter=loadDataFashionMnist(batchSize)#2.初始化模型参数#PyTorch不会隐式地调整输入的形状。因此，#我们在线性层前定义了展平层（flatten），来调整网络输入的形状net=torch.nn.Sequential(torch.nn.Flatten(),torch.nn.Linear(784,10))definitWeight(m):iftype(m)==torch.nn.Linear:torch.nn.init.normal_(m.weight,std=0.01)net.apply(initWeight)#3.定义损失函数loss=torch.nn.CrossEntropyLoss(reduction='none')#4.定义优化算法trainer=torch.optim.SGD(net.parameters(),lr=0.1)#5.训练numEpochs:int=10trainCh3(net,trainIter,testIter,loss,numEpochs,trainer)来源：檐苔