人工智能基础
常用激活函数
Sigmoid函数
\(f(x)=\frac{1}{1+e^{-x}}:(-\infty,+\infty)\to(0,1)\).
\(f'(x)=f(x)(1-f(x))\).
tanh函数
\(f(x)=\frac{e^x-e^{-x}}{e^x+e^{-x}}:(-\infty,+\infty)\to(-1,1)\).
\(f'(x)=1-f^2(x)\).
ReLU函数
\(f(x)=\max(0,x):(-\infty,+\infty)\to(0,+\infty)\).
Leaky ReLU函数
\(f(x)=\max(\alpha x,x),0<\alpha<1:(-\infty,+\infty)\to(-\infty,+\infty)\).
Softmax函数
\(f(x_i)=\frac{e^{x_i}}{\sum_j e^{x_j}}\).
损失函数
Least Square
即\(\arg \min \sum_{i=1}^n(f(x_i)-y_i)^2=\arg \min(A\beta-Y\mid A\beta-Y)\).用最小二乘法取\(\hat \beta=(A^TA)^{-1}A^TY\).
Cross Entropy
用错误的分布\(q\)来表示真实分布\(p\)的样本,则平均编码长度应该是: \[ H(p,q)=\sum_i p(i)\log(\frac{1}{q(i)})=-\sum_i p(i)\log{q(i)} \] 此为交叉熵.
特别地,当最终样本只有两个的时候,例如Logistical Regression问题,可以写成: \[ H=-(y\log a+(1-y)\log(1-a))\\ \frac{\partial H}{\partial a}=-(\frac{y}{a}-\frac{1-y}{1-a}) \]
那如果有多个呢?考虑直接对归一化条件作偏导,先有: \[ H(p,q)=-\sum_i p_i(\log{q_i}-\log(\sum_j q_j))\\ \frac{\partial H(p,q)}{\partial q_k}=-\frac{p_k}{q_k}+\frac{\sum p_k}{\sum_j q_j}\\ =-\frac{p_k}{q_k}+1 \] 再乘以softmax那里的\(q_k\),得到\(-p_k+q_k\).
神经网络实现
通过若干隐藏层,假设最后的输出层为第\(L\)层,则:
- 对于第\(l\)层,取\(\vec z_l=(W_l)^t \vec a_{l-1}+\vec b_l\).这里对\(W_l\)作转置的目的是写代码的时候需要用行向量.
- 对于第\(l\)层,取\(\vec a_l=f(\vec z_l)\),这里意味着将每一个分量对\(f\)操作.
- 对于最终答案,取误差\(\mathcal L=\frac{1}{m}(\vec y-\vec a_L\mid \vec y-\vec a_L)\).
梯度下降法
换言之就是让\(w:=w-\alpha\frac{\partial \mathcal{L}}{\partial w}\),其中\(\alpha\)是一个选定的小常数,也可以采用类似模拟退火的方式动态决定.事实上可以把各个位置分开,写作\(w_j:=w_j-\alpha \frac{\partial\mathcal{L}}{\partial w_j}\).
另外,虽然是这么写,应当见到去掉下标\(k\)的记号仍然合理,无非是逐分量做此操作,因此下面如无特殊说明,运算均采用逐分量运算.例如可以定义\(\vec a\circ \vec b\)为两个向量逐分量相乘后得到的新向量,为表区分用\(\times\)表示正常的矩阵乘法.甚至采取\((\vec a-\vec y)^2\)表示其自点积.坦白而言,笔者对此符号相当无奈,可也想不出什么更好的写法了.但总之这种写法总是强于部分参考资料上所将下标放上面的写作\(a^L\)的做法.笔者所能维持的精神数院人的唯一做法也只能是在下面加上向量符号,藉此泄愤.
顺便一提,应当见到\(\times\)和\(\circ\)这两种运算比较随意,用线性映射来理解,你这个$W$任意作用在一个向量上就行.
误差反向传播
既然要用梯度下降法,就应该把每一层的偏导都求出来.然而\(\mathcal L\)是最后一层的结果,因此应该用链式法则一路求出前面的偏导.
更具体地,不妨设误差函数选的是\((\vec a-\vec y)^2\),激活函数选的是cross entropy有:
- \(\frac{\partial{\mathcal L}}{\partial \vec a_{L}}=(\vec a_{L}-\vec y)\).(如若选择不同的误差函数,这里作适当变化)
- \(\frac{\partial \vec a_{l}}{\partial \vec z_{l}}=f'( \vec z_{l})=\vec a_{l}\circ (1-\vec a_{l})\).(如若选取不同的激活函数,这里作适当变化)
- \(\frac{\partial \vec z_{l+1}}{\partial \vec z_{l}}=\frac{\partial \vec z_{l+1}}{\partial \vec a_{l}}\frac{\partial \vec a_{l}}{\partial \vec z_{l}}=( W_{l+1})^t\times(\vec a_l\circ (1-\vec a_l))\).
- \(\frac{\partial \vec z_{l}}{\partial W_{l}}=(\vec a_{l-1})\).结果理应是一个矩阵,其实就是这个列向量不断复制若干遍,或者写成\((\vec a_{l-1})^tM(1)\),其中\(M(1)\)是全\(1\)矩阵.
- \(\frac{\partial \vec z_{l}}{\partial \vec b_{l}}=1\).
我们应当见到:
不妨设\(\delta_l=\frac{\partial {\mathcal L}}{\partial \vec z_l}\).见到:
- \(\delta_L=\frac{\partial\mathcal L}{\partial \vec z_L}=\frac{\partial\mathcal L}{\partial \vec a_L}\frac{\partial\vec a_L}{\partial \vec z_L}=(\vec a_{L}-\vec y)\circ \vec a_{L}\circ (1-\vec a_{L})\).前者会因为误差函数的选取而改变,后者会因为激活函数的选取而改变.
- \(\delta_l=\frac{\partial {\mathcal L}}{\partial \vec z_l}=\delta_{l+1}\frac{\partial \vec z_{l+1}}{\partial \vec z_l}=\delta_{l+1}\circ (W_{l+1})^t\times (\vec a_{l}\circ (1-\vec a_{l}))\).
- \(\frac{\partial \mathcal L}{\partial W_{l}}=\frac{\partial \mathcal L}{\partial \vec z_{l}}\frac{\partial \vec z_l}{\partial W_{l}}=\delta_l\times a_{l-1}^t\).
- \(\frac{\partial \mathcal L}{\partial \vec b_l}=\delta_l\).
如此以上更新即可.
卷积神经网络(CNN)
神经网络受矩阵乘法的限制,导致对于真实的尺寸巨大的图像难以快速识别,因此产生了卷积神经网络的概念,大概有以下特征:
- 空间上权值共享:不同位置使用同一个卷积核(滤波器)
- 稀疏链接:每一层只链接前一层的感受野.
- 等变表示:卷积神经网络有某种平移不变性.
对于2D卷积,其公式如下: \[ S_{r,c}=(X*W)_{r,c}=\sum_i\sum_j X_{r+i,c+j}\times w_{i,j} \] 其中\(W\)是卷积核,\(X\)是输入图像,\(S\)是输出的结果.如果一个图像有多个通道(比如色彩层之类的),每个通道上都需要应用一个卷积核.
下面引入一些名词:
- input size:输入图像的尺寸.
- padding:填充的像素数.
- filter size:卷积核的尺寸.有时也写作两个变量:filter height和filter width.3D卷积还会有一个filter depth的变量.
- stride:步长.
- output size:卷积后输出的尺寸.有时也写作feature size.
- input channels:输入图像的通道数.
- n filters:卷积核的数量.
- dilation rate:膨胀率,用于空洞卷积.膨胀率为\(d\)的时候,卷积核中间会插入\(d-1\)个\(0\)间隔.
感受野计算
先看output size的计算,容易见到,其各个维度方面计算是独立的.只要对于单个维度算出卷积核在上面移动的次数,最后将不同维度相乘即可.
对于单个维度,这个维度的移动次数应该是: \[ \text{output\_size}=\lceil\frac{\text{input\_size}+2\times \text{padding}-\text{filter\_size}+1}{\text{stride}}\rceil\\ =\lfloor\frac{\text{input\_size}+2\times \text{padding}-\text{filter\_size}}{\text{stride}}\rfloor+1\\ \] 这个公式相当容易理解,原因是\(\text{stride}=1\)的时候,上面恰好是移动的次数,而\(\text{stride}\)变化的时候,当然要拿到一个上取整.
至于所谓的空洞卷积,只需在上面的基础上改\(\text{filter\_size}\)就好.
至于乘法操作,每得到一个\(\text{output}\)当然都会需要\(\text{filter\_size}\)次乘法操作.
池化(Pooling)
池化操作它没有一个可学习的参数,只是对输入数据进行固定的操作.简单来说就是降低输入的规模,以实现更好的鲁棒性以及提高效率.
常见的池化操作包括:
- MaxPooling:取区域内的最大值.
- MeanPooling:取区域内的平均值.
- PyramidPooling:多次进行尺度不同的池化.
常见卷积架构
AlexNet
首次引入ReLU激活函数,Dropout 技术,以及数据增强,提高了模型的训练效率和泛化能力.
采用了\(8\)层深的网络结构,证明了深度网络的潜力.
VGG
开始堆叠小尺寸的卷积核,获得与大卷积核相似的感受野的同时可以增加网络深度.
ResNet
引入残差的概念,直接将输入数据累加(跳跃连接)到最后的输出中,这样网络学习的实际上是输入和输出之间的残差,从而提高了网络学习能力.
SqueezeNet
SqueezeNet的基本构建单元是Fire模块.Fire模块由一个squeeze层和一个expand层组成.squeeze层使用\(1\times1\)卷积核减少通道数,而expand层则使用\(1\times 1\)和\(3\times 3\)卷积核增加通道数.这种设计有效地减少了参数数量和计算量.
MobileNet
- 深度卷积:在这个操作中,每个输入通道独立地进行卷积,这意味着在进行卷积时,不同通道之间没有交互.这样可以减少计算量和参数数量.
- 逐点卷积:逐点卷积使用\(1\times 1\)的卷积核,它作用在深度卷积的输出上,将不同通道的信息整合在一起.逐点卷积可以减少参数数量,同时保持较高的性能.
ShuffleNet
- 组卷积(Group Convolution):将通道分成几个组,并使用不同的卷积神经网络层执行标准卷积.
- 打乱层(Shuffle layer):通过对通道进行洗牌,将不同组的信息合并.
反卷积
也就是将较小的数据特征图扩大到较大的尺寸.有的时候也把这个操作说成上采样.
- 插值步骤(Interpolation Step):首先,在输入特征图的元素之间插入零,增加特征图的尺寸.
- 卷积步骤(Convolution Step):接下来,对扩大后的特征图应用一个标准的卷积操作.此步骤相当于在扩大的特征图上滑动卷积核,计算卷积输出.