【论文笔记】Towards Certifying l-infinity robustness using neural networks with l-infinity-dist neurons
作者:互联网
本文发表于 2021 ICML,提出了一个新颖的神经网络计算方式:对于网络中的每个神经元,不采用传统的线性转换+非线性激活函数的方式,而是计算输入与参数之间的 \(\ell_{\infty}\)-distance,作者将其称为 \(\ell_{\infty}\)-dist net,网络中的神经元称为 \(\ell_{\infty}\)-dist neuron。作者理论证明了 \(\ell_{\infty}\)-dist net 具有很好地表达能力(expressivity)和泛化能力(generalization ability),还给出了 \(\ell_{\infty}\)-dist net 在训练中的优化策略。\(\ell_{\infty}\)-dist net 还可以作为特征提取器与其他模型结构(如卷积网络)结合使用,实验发现这样的设计在很多数据集上都能获得很好的 certified robustness。
1. \(\ell_{\infty}\)-dist net
1.1) Preliminaries
问题描述 考虑一个标准的分类任务:
- 样本 \(x \in \mathcal{X}\) 和相应标签 \(y\in\mathcal{Y}=\{1,...,M\}\) 的数据分布为 \(\mathcal{D}\),\(\mathcal{D}\) 通常是已知的。
- 训练集为 \(\tau=\{(x_1,y_1),...,(x_n,y_n)\}\),其中 \((x_i,y_i)\) i.i.d. 取自分布 \(\mathcal{D}\)。
- \(f:x\in\mathcal{X}\rightarrow \mathcal{Y}\) 是分类器。
- \(x'=x+\delta\) 是对抗样本。
- 扰动约束:\(\epsilon\)-bounded \(\ell_\infty\)-norm constraint(或称为 \(\ell_\infty\) 扰动),即 \(||\delta||_\infty \leq \epsilon\)。
任务目标: 使用训练集 \(\tau\) 学习一个模型,该模型可以抵抗带有任意\(\ell_\infty\) 扰动的样本 \((x,y)\),其中\((x,y)\sim \mathcal{D}\)。
这需要计算以 \(x\) 为中心的、不会改变 \(f\) 对它预测的 \(\ell_\infty\)-ball 的最大半径(称为 robust radius):
\[R(f;x,y)= \left\{\begin{matrix} \inf_{f(x')\neq f(x)} ||x'-x||_\infty & ,f(x)=y\\ 0& ,f(x)\neq y \end{matrix}\right. \tag{1} \]但对于标准的 DNNs,robust radius 很难计算,所以转为计算 \(R(f;x,y)\) 的下限 \(CR(f;x,y)\),称为 certified radius。对于任意的 \(f,x,y\) 有 \(CR(f;x,y)\leq R(f;x,y)\)。
1.2) \(\ell_{\infty}\)-dist neurons
图1的左图是传统的卷积神经元计算方式(线性变换+非线性激活函数),右图是 \(\ell_{\infty}\)-dist neuron 的计算方式:
\[u(x,\theta)=||x-w||_\infty + b \tag{2} \]其中 \(\theta = \{w,b\}\) 是参数集合。公式(2)本身就是非线性的,所以不需要像传统网络那样添加一个激活函数 \(\sigma\)。传统神经元使用点积计算来度量 \(x\) 与 \(w\) 之间的相似度(similarity),同样地,\(\ell_{\infty}\)-dist neuron 使用距离公式作为度量,且距离是非负数,值越小表示相似度越强。
1.3) MLP networks using \(\ell_{\infty}\)-dist neurons
考虑将 \(\ell_{\infty}\)-dist neurons 应用于最简单的模型结构 MLP:
定义一个 \(L\) 层的 \(\ell_{\infty}\)-dist net:假设第 \(l\) 个隐藏层有 \(d_l\) 个隐藏单元,网络的输入为 \(x^{(0)}\triangleq x \in \mathbb{R}^{d_{input}}\),第 \(l\) 个隐藏层中的第 \(k\) 个神经元的输出为:
\[x_k^{(l)}=u(x^{(l-1)},\theta^{(l,k)})=||x^{(l-1)}-w^{(l,k)}||_\infty+b^{(l,k)} \tag{3} \]其中 \(x^{(l)}=(x^{(l)}_{1},x^{(l)}_{2},...,x^{(l)}_{d_{l}})\) 为第 \(l\) 层的输出,\(1 \leq l\leq L\),\(1 \leq k \leq d_l\)。
对于 1.1 节描述的分类任务,输出维度 \(d_L\) 等于类别 \(M\)。取网络最后一层的输出的负数用于预测,即 \(g(x)=(-x^{(L)}_1,-x^{(L)}_2,...,-x^{(L)}_M)\),预测输出为 \(f(x) = \arg\max_{i\in [M]}g_i(x)\)。与标准网络一样,可以对\(\ell_{\infty}\)-dist net 使用任何标准的损失函数,比如交叉熵或 hinge loss。
1.4) 1-Lipschitz w.r.t. \(\ell_{\infty}\)-norm
这节首先证明 \(\ell_{\infty}\)-dist net 就是 1-Lipschitz w.r.t. \(\ell_{\infty}\)-norm;然后根据这一性质推导出模型的 certified robustness。
\[||g(z_1)-g(z_2)||_p\leq\lambda||z_1 - z_2||_p \]定义: 如果函数 \(g(z):\mathbb{R}^{m} \rightarrow \mathbb{R}^{n}\) 对任意 \(z_1,z_2\) 都满足下式,则 \(g(z)\) 被称为 \(\lambda\)-Lipschitz w.r.t. \(\ell_{p}\)-norm (\(||\cdot||_p\)):
Fact 1: \(\ell_{\infty}\)-dist net \(g(\cdot)\) 是 1-Lipschitz w.r.t. \(\ell_{\infty}\)-norm,即对任意 \(x_1,x_2 \in \mathbb{R}^{d_{input}}\;\;\) 都有 \(||g(x_1)-g(x_2)||_\infty\leq||x_1 - x_2||_\infty\).
Proof 1: 网络中每个神经元的计算(即公式3)是 1-Lipschitz → 从 \(x^{(l)}\) 到 \(x^{(l+1)}\) (层到层)的映射是 1-Lipschitz → 整个网络是 1-Lipschitz。
因此,当扰动很小时,输出的变化是有界的,直接约束了 certified radius。
Fact 2: 设 \({\rm margin}(x;g)\) 是输出向量 \(g(x)\) 中最大元素和第二大元素之间的差,那么对于任意 \(x'\) 满足 \(||x-x'||_\infty < {\rm margin}(x;g)/2\),有 \(f(x)=f(x')\)。即:
\[CR(f,x,y)\geq {\rm margin}(x;g)/2 \tag{4} \]Proof 2: \(g(x)\) 是 1-Lipschitz,因此当将输入从 \(x\) 变为 \(x'\)时,模型输出 \(g(x)\) 中的每个元素移动不超过 \({\rm margin}(x;g)/2\),那么最大的元素(即预测的类)是不变的。
使用 Fact 2 中的 bound,仅需一次正向传播就能够计算 \(\ell_{\infty}\)-dist net 的 certified robustness,计算成本小。
2. \(\ell_{\infty}\)-dist net 的训练
经验发现传统网络训练方法并不适用于 \(\ell_{\infty}\)-dist net,对此作者提出了一系列相应的优化策略。
2.1) Normalization
问题: 传统网络的线性层的输出是无偏的 unbiased(期望均值为0),而 \(\ell_{\infty}\)-dist neuron 的输出是有偏的 biased(假设没有bias项 \(b\),总是非负的)。这会导致每层的输出会随层数增加而线性增长。
解决方法:
考虑 Batch Normalization (BN),BN 使用 shift 和 scale 两个操作,但若直接在 \(\ell_{\infty}\)-dist net 中使用 BN 会导致 Lipschitz 常数发生变化(由于 scale 操作),从而无法保证模型的鲁棒性。
不过作者发现,只使用 shift 操作有助于优化,因此在所有的中间层计算完距离后添加 shift 操作,并移除了 bias 项 \(b\)(冗余了),但最后一层不做 normalization。与 BN 类似,在 inference 时使用 running mean。
2.2) Smoothed Approximated Gradients
问题: \(\ell_{\infty}\)-dist 的梯度向量(如 $\triangledown_w ||z-w||_\infty $ 和 \(\triangledown_z ||z-w||_\infty\))十分稀疏,通常只包含一个非零元素。通过实验观察到,如果直接使用 SGD/Adam 训练(随机初始化的) \(\ell_{\infty}\)-dist net,那么在每个 epoch 中,只有不到 1% 的参数在更新。
解决方法:
用 \(\ell_{p}\)-dist neuron 替换整个网络的神经元,取得近似的 & 非稀疏的参数梯度。在训练中,最开始时将 \(p\) 设置为一个很小的值,在接下来的每次迭代中不断增加 \(p\) 的值,直到逼近无穷。在最后几次 epochs 中,将 \(p\) 设置为无穷。
2.3) Parameter Initialization
问题: 深层模型的训练准确率比浅层模型的差
解决方法:
参考 ResNet ,可以在初始化 weights 和 biases 时直接构建恒等映射(identity mapping)。具体来说,对于输入-输出维度相同的 \(\ell_{\infty}\)-dist layer,首先用标准高斯分布随机初始化 weights,然后将对角元素(即公式3中的 \(w^{(l,j)}_j\),\(l\) 层第 \(j\) 个神经元与 \(l-1\) 层第 \(j\) 个神经元之间的weight)修改为一个很大的负数 \(C_0\)。在实验中作者设置 \(C_0 = -10\)。当应用了 mean shift normalization 后,不再需要添加偏差 biases,并且 running mean 会自动进行恒等映射。
2.4) Weight Decay
问题: 将 Weight Decay 应用到 \(\ell_{\infty}\)-dist net 会使模型的性能变差,可能是由于 Weight Decay 与 \(\ell_{\infty}\)-norm 不兼容。
解决方法:
传统网络的计算方式是点积,所以权重的 \(\ell_{2}\)-norm 可以控制输出的大小。而 \(||w||_2\) 与 \(\ell_{\infty}\)-dist layer 的输出大小无关。所以对于 Weight Decay Regularizer,可以用 \(||w||_{\infty}\) 取代 \(||w||_2\)。类似于 \(\left \langle x,w \right \rangle \leq ||x||_2||w||_2\),我们有 \(||x-w||_\infty \leq ||x||_\infty + ||w||_\infty\)。
对于训练中的 \(\ell_p\)-dist neurons,使用 \(\ell_p\)-norm regularization。通过对权重 \(w\) 的求导,有关权重的 weight decay 公式为:
\[\triangle_{w_{i}}=-\lambda \triangledown_{w_{i}}||w||^2_p=-\lambda\left ( \frac{|w_i|}{||w||_p}\right )^{p-2}w_i \tag{5} \]其中 \(\lambda\) 是 weight decay 的系数,当 \(p\rightarrow \infty\) 时,weight decay 往往只对绝对值最大的元素 \(w_i\) 产生影响。
3. 实验
3.1) 实验设置
四个基准数据集:MNIST, Fashion-MNIST, CIFAR-10, TinyImagenet
> 模型配置:
主要研究两种模型:1) 仅 \(\ell_\infty\)-dist net;2) \(\ell_\infty\)-dist net + MLP:\(\ell_\infty\)-dist net 作为特征提取器。
- 对于 MNIST 和 Fashion-MNIST,使用 5 层 \(\ell_\infty\)-dist net;对于 CIFAR-10 和TinyImageNet,使用 6 层 \(\ell_\infty\)-dist net。
- 每个隐藏层包含 5120 个单元。
- 每个中间层都应用 normalization。
- Top 层包含 10 个单元 (TinyImageNet 为 200 个单元)
- 对于 \(\ell_\infty\)-dist net + MLP,移除 top 层,并添加传统的全连接层,隐藏层包含 512 个单元,且添加了 tanh 激活函数。
> 训练配置:
- optimizer: 使用 Adam, 其中 \(\beta_1 = 0.9, \beta_2 = 0.99, \epsilon = 10^{-10}\).
- data augmentation: 对于 MNIST 和 Fashion-MNIST, 使用随机剪裁 random crop (padding = 1); 对于 CIFAR-10,使用随机剪裁 (padding = 4) 和 random horizontal flip; 对于 TinyImageNet,使用 random horizontal flip, 训练中将每张图剪裁到 $ 56 \times 56$ 像素,测试中使用中心剪裁 center crop.
- 损失函数: 对于 \(\ell_\infty\)-dist net, 使用 hinge loss, threshold 超参数 \(t\) 取决于 \(\epsilon\); 对于 \(\ell_\infty\)-dist net + MLP, 使用 IBP loss, 涉及两个超参数 \(\kappa, \epsilon_{train}\).
- 训练过程: 首先, 用 \(\ell_p\) 代替 \(\ell_\infty\), 并令 \(p=8\),此时训练 \(e_1\) 个 epochs; 然后逐渐将 \(p\) 从 \(8\) exponentially 增大到 \(1000\), 此时训练 \(e_2\) 个 epochs; 最后将 \(p\) 设为无穷, 此时训练 \(e_3\) 个 epochs. 其中 \(e_1,e_2,e_3\) 是超参数, 不同数据集设置的值不同.
- \(e_1\) epochs 中 \(lr = 0.02\), 在接下来的 \(e_2,e_3\) epochs 中使用 cosine annealing 降低学习率.
- weight decay: \(\lambda=0.005\). 对于 \(\ell_\infty\)-dist nets, 使用 \(\ell_p\)-norm weight decay; 对于 MLP, 使用 \(\ell_2\)-norm weight decay.
> Evaluation:
使用两种指标来评估模型的 robustness: 1) robust test accuracy: 使用 PGD 攻击, 攻击步长设为 \(20\); 2) certified radius: 计算每个样本的 CR, 并计算在 CR 内的测试样本的百分比. Note: 第二个指标始终 lower than 第一个指标.
> Baselines:
对比了先进的方法, 包括: 1) relaxation methods: CAP, PVT, DiffAI, IBP, CROWN-IBP, CROWN-IBP with loss function, COLT; 2) Lipschitz networks: GroupSort.
3.2) 实验结果:
"Test" 表示干净样本的测试准确率; "Robust" 表示PGD 样本的测试准确率; "Certified" 表示 certified robust 测试准确率. "FLOPs" 表示前向传播中所需的基本浮点运算的数量 (即传统网络中的点积和加法或 \(\ell_\infty\)-dist net 中的减法).
> General profermance of \(\ell_\infty\)-dist net
从表 1 中可以看到,单独使用 \(\ell_\infty\)-dist net 已经在所有数据集上获得了不错的 certified accuracy. 尤其是, 在 CIFAR10 数据集上达到了最好的 certified accuracy,且获得了比其他方法更高的标准准确率(干净样本).
Note: 只使用标准损失函数来训练 \(\ell_\infty\)-dist net, 无需任何对抗训练。
> General profermance of \(\ell_\infty\)-dist net + MLP
如表1和表2, 对于所有数据集,\(\ell_\infty\)-dist net + MLP 比单独的 \(\ell_\infty\)-dist net 获得了更好的 certified accuracy.
> Efficiency
如表4, 训练和 certification 都很快. 训练 \(\ell_\infty\)-dist net 的计算成本与训练相同大小的常规网络大致相同,并且 certification 过程只需要一次向前传播 forward pass.
> 与 GroupSort Network 对比
由于 GroupSort 也使用了1-Lipschitz 且使用标准损失函数训练即可(不需要对抗训练), 作者特别地将这两个模型进行了比较. 在 GroupSort Network 中,所有权重矩阵 \(W\) 都被限制为 bounded \(\ell_\infty\)-norm,即 \(||W||_\infty \leq 1\),这导致了耗时的 projection 操作, 带来了优化难度,进一步限制了网络结构的可扩展性. 作者将 \(\ell_\infty\)-dist net 在 MNIST 数据集上显着优于 GroupSort 的原因也解释为这一点.
> Ablation Studies
作者还实验观察前述的 smoothed approximated gradients, parameter initialization(使用 identity map 构建) 和 \(\ell_p\)-norm weight decay 的影响, 结果如表 3 所示, 可以看出:
- smoothed approximated gradients 对 \(\ell_\infty\)-dist net 的训练十分重要, 添加它后模型的 certified accuracy 可达到 32.56%.
- smoothed approximated gradients 和 parameter initialization 均对 \(\ell_\infty\)-dist net+MLP 训练十分重要, 结合它们一同使用后, 模型的 certified accuracy 达到 35.02%.
- \(\ell_\infty\)-norm weight decay 可以进一步提升结果,尽管效果可能很小(两个模型的 certified accuracy 分别提高了 0.59% 和 0.4%).
- 传统的 \(\ell_2\)-norm weight decay 会损害 \(\ell_\infty\)-dist net 的性能.
总之, 这几个训练策略都对模型的性能有帮助.
标签:infty,Towards,dist,训练,ell,neural,infinity,net,norm 来源: https://www.cnblogs.com/setdong/p/16456887.html