返回介绍

数学基础

统计学习

深度学习

工具

Scala

二、线性支持向量机

发布于 2023-07-17 23:38:26 字数 7761 浏览 0 评论 0 收藏 0

  1. 对于线性不可分训练数据,线性支持向量机不再适用,但可以想办法将它扩展到线性不可分问题。

2.1 原始问题

  1. 设训练集为 $ MathJax-Element-511 $ ,其中 $ MathJax-Element-329 $ 。

    假设训练数据集不是线性可分的,这意味着某些样本点 $ MathJax-Element-467 $ 不满足函数间隔大于等于 1 的约束条件。

    • 对每个样本点 $ MathJax-Element-467 $ 引进一个松弛变量 $ MathJax-Element-195 $ ,使得函数间隔加上松弛变量大于等于 1。

      即约束条件变成了: $ MathJax-Element-196 $ 。

    • 对每个松弛变量 $ MathJax-Element-382 $ ,支付一个代价 $ MathJax-Element-382 $ 。目标函数由原来的 $ MathJax-Element-199 $ 变成:

      $ \min \frac 12 ||\mathbf {\vec w}||^{2}_2+C\sum_{i=1}^{N}\xi_i $

      这里 $ MathJax-Element-200 $ 称作惩罚参数,一般由应用问题决定。

      • $ MathJax-Element-504 $ 值大时,对误分类的惩罚增大,此时误分类点凸显的更重要
      • $ MathJax-Element-504 $ 值较大时,对误分类的惩罚增加,此时误分类点比较重要。
      • $ MathJax-Element-504 $ 值较小时,对误分类的惩罚减小,此时误分类点相对不重要。
  2. 相对于硬间隔最大化, $ MathJax-Element-204 $ 称为软间隔最大化。

    于是线性不可分的线性支持向量机的学习问题变成了凸二次规划问题:

    $ \min_{\mathbf {\vec w},b,\vec\xi} \frac 12||\mathbf {\vec w}||^{2}_2+C\sum_{i=1}^{N}\xi_i\\ s.t. \quad \tilde y_i(\mathbf {\vec w}\cdot \mathbf {\vec x}_i+b) \ge 1-\xi_i,\quad i=1,2,\cdots,N\\ \xi_i \ge 0,\quad i=1,2,\cdots,N $
    • 这称为线性支持向量机的原始问题。

    • 因为这是个凸二次规划问题,因此解存在。

      可以证明 $ MathJax-Element-205 $ 的解是唯一的; $ MathJax-Element-493 $ 的解不是唯一的, $ MathJax-Element-493 $ 的解存在于一个区间。

  3. 对于给定的线性不可分的训练集数据,通过求解软间隔最大化问题得到的分离超平面为: $ MathJax-Element-238 $ ,

    以及相应的分类决策函数: $ MathJax-Element-209 $ ,称之为线性支持向量机。

  • 线性支持向量机包含线性可分支持向量机。
  • 现实应用中训练数据集往往是线性不可分的,线性支持向量机具有更广泛的适用性。

2.2 对偶问题

  1. 定义拉格朗日函数为:

    $ L(\mathbf {\vec w},b,\vec\xi,\vec\alpha,\vec\mu)=\frac 12||\mathbf {\vec w}||_2^{2}+C\sum_{i=1}^{N}\xi_i-\sum_{i]1}^{N}\alpha_i[\tilde y_i(\mathbf {\vec w}_i \cdot \mathbf {\vec x}_i+b)-1+\xi_i]-\sum_{i=1}^{N}\mu_i\xi_i \\ \alpha_i \ge 0, \mu_i \ge 0 $

    原始问题是拉格朗日函数的极小极大问题;对偶问题是拉格朗日函数的极大极小问题。

    • 先求 $ MathJax-Element-210 $ 对 $ MathJax-Element-211 $ 的极小。根据偏导数为0:

      $ \nabla_{\mathbf {\vec w}}L(\mathbf {\vec w},b,\vec\xi,\vec\alpha,\vec\mu)=\mathbf {\vec w}-\sum_{i=1}^{N}\alpha_i\tilde y_i\mathbf {\vec x}_i=\mathbf{\vec 0}\\ \nabla_b L(\mathbf {\vec w},b,\vec\xi,\vec\alpha,\vec\mu)=-\sum_{i=1}^{N}\alpha_i\tilde y_i=0\\ \nabla_{\xi_i} L(\mathbf {\vec w},b,\vec\xi,\vec\alpha,\vec\mu)=C-\alpha_i-\mu_i=0 $

      得到:

      $ \mathbf {\vec w}=\sum_{i=1}^{N}\alpha_i\tilde y_i\mathbf {\vec x}_i\\ \sum_{i=1}^{N}\alpha_i\tilde y_i=0\\ C-\alpha_i-\mu_i=0 $
    • 再求极大问题:将上面三个等式代入拉格朗日函数:

      $ \max_{\vec\alpha,\vec\mu}\min_{\mathbf {\vec w},b,\vec\xi}L(\mathbf {\vec w},b,\vec\xi,\vec\alpha,\vec\mu)=\max_{\vec\alpha,\vec\mu}\left[-\frac 12 \sum_{i=1}^{N}\sum_{j=1}^{N}\alpha_i\alpha_j\tilde y_i\tilde y_j(\mathbf {\vec x}_i \cdot \mathbf {\vec x}_j)+\sum_{i=1}^{N}\alpha_i\right] $

      于是得到对偶问题:

      $ \min_{\vec\alpha}\frac 12 \sum_{i=1}^{N}\sum_{j=1}^{N}\alpha_i\alpha_j\tilde y_i\tilde y_j(\mathbf {\vec x}_i \cdot \mathbf {\vec x}_j)-\sum_{i=1}^{N}\alpha_i\\ s.t. \quad \sum_{i=1}^{N}\alpha_i\tilde y_i=0\\ 0 \le \alpha_i \le C, i=1,2,\cdots,N $
  2. 根据 KKT 条件:

    $ \nabla_{\mathbf {\vec w}}L(\mathbf {\vec w}^*,b^*,\vec\xi^*,\vec\alpha^*,\vec\mu^*)=\mathbf {\vec w}^*-\sum_{i=1}^{N}\alpha_i^*\tilde y_i\mathbf {\vec x}_i=\mathbf{\vec 0}\\ \nabla_b L(\mathbf {\vec w}^*,b^*,\vec\xi^*,\vec\alpha^*,\vec\mu^*)=-\sum_{i=1}^{N}\alpha_i^*\tilde y_i=0\\ \nabla_{\xi_i} L(\mathbf {\vec w}^*,b^*,\vec\xi^*,\vec\alpha^*,\vec\mu^*)=C-\alpha_i^*-\mu_i^*=0\\ \alpha_i^*[\tilde y_i(\mathbf{\vec w}^*\cdot \mathbf{\vec x}_i+b^*)-1+\xi_i^*]=0\\ \mu_i^*\xi_i^*=0\\ \tilde y_i(\mathbf{\vec w}^*\cdot\mathbf{\vec x}_i+b^*)-1+\xi_i^*\ge 0\\ \xi_i^*\ge 0\\ C\ge \alpha_i^*\ge 0\\ \mu_i^*\ge0\\ i=1,2,\cdots,N $

    则有: $ MathJax-Element-212 $ 。

    • 若存在 $ MathJax-Element-335 $ 的某个分量 $ MathJax-Element-214 $ ,则有: $ MathJax-Element-215 $ 。

      若 $ MathJax-Element-218 $ 的所有分量都等于 0 ,则得出 $ MathJax-Element-217 $ 为零,没有任何意义。

      若 $ MathJax-Element-218 $ 的所有分量都等于 $ MathJax-Element-504 $ ,根据 $ MathJax-Element-220 $ ,则要求 $ MathJax-Element-221 $ 。这属于强加的约束,

      • 根据 $ MathJax-Element-245 $ ,有 $ MathJax-Element-223 $ 。
      • 考虑 $ MathJax-Element-224 $ ,则有: $ MathJax-Element-225 $
    • 分离超平面为: $ MathJax-Element-226 $ 。

    • 分类决策函数为: $ MathJax-Element-227 $ 。

  3. 线性支持向量机对偶算法:

    • 输入:训练数据集 $ MathJax-Element-511 $ ,其中 $ MathJax-Element-329 $

    • 输出:

      • 分离超平面
      • 分类决策函数
    • 算法步骤:

      • 选择惩罚参数 $ MathJax-Element-384 $ ,构造并且求解约束最优化问题:

        $ \min_{\vec\alpha} \frac 12 \sum_{i=1}^{N}\sum_{j=1}^{N}\alpha_i\alpha_j\tilde y_i\tilde y_j(\mathbf {\vec x}_i \cdot \mathbf {\vec x}_j) -\sum_{i=1}^{N} \alpha_i\\ s.t. \quad \sum_{i=1}^{N}\alpha_i\tilde y_i=0\\ C \ge \alpha_i \ge 0,i=1,2,\cdots,N $

        求得最优解 $ MathJax-Element-332 $ 。

      • 计算 : $ MathJax-Element-334 $ 。

      • 选择 $ MathJax-Element-335 $ 的一个合适的分量 $ MathJax-Element-336 $ ,计算: $ MathJax-Element-235 $ 。

        可能存在多个符合条件的 $ MathJax-Element-236 $ 。这是由于原始问题中,对 $ MathJax-Element-493 $ 的解不唯一。所以

        实际计算时可以取在所有符合条件的样本点上的平均值。

      • 由此得到分离超平面: $ MathJax-Element-238 $ ,以及分类决策函数: $ MathJax-Element-239 $ 。

2.3 支持向量

  1. 在线性不可分的情况下,对偶问题的解 $ MathJax-Element-332 $ 中,对应于 $ MathJax-Element-241 $ 的样本点 $ MathJax-Element-467 $ 的实例点 $ MathJax-Element-243 $ 称作支持向量,它是软间隔的支持向量。

  2. 线性不可分的支持向量比线性可分时的情况复杂一些:

    根据 $ MathJax-Element-244 $ ,以及 $ MathJax-Element-245 $ ,则:

    • 若 $ MathJax-Element-246 $ ,则 $ MathJax-Element-247 $ , 则松弛量 $ MathJax-Element-248 $ 。此时:支持向量恰好落在了间隔边界上。

    • 若 $ MathJax-Element-249 $ , 则 $ MathJax-Element-250 $ ,于是 $ MathJax-Element-382 $ 可能为任何正数:

      • 若 $ MathJax-Element-252 $ ,则支持向量落在间隔边界与分离超平面之间,分类正确。
      • 若 $ MathJax-Element-253 $ ,则支持向量落在分离超平面上。
      • 若 $ MathJax-Element-254 $ ,则支持向量落在分离超平面误分类一侧,分类错误。

2.4 合页损失函数

  1. 定义取正函数为:

    $ \text{plus}(z)= \begin{cases} z, & z \gt 0 \\ 0, & z \le 0 \end{cases} $

    定义合页损失函数为: $ MathJax-Element-255 $ ,其中 $ MathJax-Element-256 $ 为样本的标签值, $ MathJax-Element-257 $ 为样本的模型预测值。

    则线性支持向量机就是最小化目标函数:

    $ \sum_{i=1}^{N}\text{plus}(1-\tilde y_i(\mathbf {\vec w}\cdot\mathbf {\vec x}_i+b))+\lambda||\mathbf {\vec w}||^{2}_2,\quad \lambda \gt 0 $
  2. 合页损失函数的物理意义:

    • 当样本点 $ MathJax-Element-467 $ 被正确分类且函数间隔(确信度) $ MathJax-Element-261 $ 大于 1 时,损失为0
    • 当样本点 $ MathJax-Element-467 $ 被正确分类且函数间隔(确信度) $ MathJax-Element-261 $ 小于等于 1 时损失为 $ MathJax-Element-264 $
    • 当样本点 $ MathJax-Element-467 $ 未被正确分类时损失为 $ MathJax-Element-264 $
  3. 可以证明:线性支持向量机原始最优化问题等价于最优化问题:

    $ \min_{\mathbf {\vec w},b}\sum_{i=1}^{N}\text{plus}(1-\tilde y_i(\mathbf {\vec w}\cdot\mathbf {\vec x}_i+b))+\lambda||\mathbf {\vec w}||^{2}_2,\quad \lambda \gt 0 $
  4. 合页损失函数图形如下:

    • 感知机的损失函数为 $ MathJax-Element-265 $ ,相比之下合页损失函数不仅要分类正确,而且要确信度足够高(确信度为1)时,损失才是0。即合页损失函数对学习有更高的要求。

    • 0-1损失函数通常是二分类问题的真正的损失函数,合页损失函数是0-1损失函数的上界。

      • 因为0-1损失函数不是连续可导的,因此直接应用于优化问题中比较困难。
      • 通常都是用0-1损失函数的上界函数构成目标函数,这时的上界损失函数又称为代理损失函数。

    hinge_loss

  5. 理论上SVM 的目标函数可以使用梯度下降法来训练。但存在三个问题:

    • 合页损失函数部分不可导。这可以通过sub-gradient descent 来解决。
    • 收敛速度非常慢。
    • 无法得出支持向量和非支持向量的区别。

如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。

扫码二维码加入Web技术交流群

发布评论

需要 登录 才能够评论, 你可以免费 注册 一个本站的账号。
列表为空,暂无数据
    我们使用 Cookies 和其他技术来定制您的体验包括您的登录状态等。通过阅读我们的 隐私政策 了解更多相关信息。 单击 接受 或继续使用网站,即表示您同意使用 Cookies 和您的相关数据。
    原文