HELLO WORLD--一起加油(🍺)!|

kingwzun

园龄:3年6个月粉丝:111关注:0

机器学习相关作业/公式应用

work4

证明函数是否是凸函数

题目:要求证明两个性质:

  1. Logistic函数 f(x;β)=11+eβTx 对参数 β 是非凸的
  2. 对数似然函数 L(β)=y0βTx+ln(1+eβTx)β 是凸的,其中 y0 是常数。

基本知识补充

问题一:f(x;β) 什么含义?
可以理解为函数f(x;β)xβ两个变量。
但是 f(x;β) 通常表示一个参数化的函数,xβ 的先后顺序在符号上有不同的含义,通常,会把参数 β 写在后面,以强调 x 是输入特征,而 β 是控制模型行为的参数。具体来说:

  • x 是输入变量(特征向量),例如在机器学习模型中代表数据样本的特征。
  • β 是参数向量,表示模型的参数,需要通过训练来优化。

进一步说明,假设我们有一个线性回归模型:f(x;β)=βTxf(β;x)=βTx。从数学上来说,这两者给出的输出值是相同的,因为内积运算在交换变量时不受影响。但是

  • 前者强调:在固定 β 的情况下,x 的变化如何影响输出。
  • 后者强调:在固定 x 的情况下,β 的变化如何影响输出。

举例: Logistic 回归中的 f(x;β) 在逻辑回归中,f(x;β) 通常定义为:f(x;β)=σ(βTx)
其中:

  • βTx 是参数 β 和特征向量 x 的内积,结果是一个标量。
  • σ(z) 是 Sigmoid 函数,定义为:σ(z)=11+ez

这个函数的输出是一个介于 0 和 1 之间的数值,表示输入特征 x 属于某一类的概率。

解释含义:

  • 输入 x:代表特征或数据样本,它可以是多维的向量。
  • 参数 β:表示模型参数的集合,它们是训练过程中需要调整的,以使模型更好地拟合训练数据。
  • 输出 f(x;β):模型对输入 x 的输出结果,通常代表一种预测值或概率。

问题二:βT是什么
βT 表示向量 β 的转置。它是将一个列向量变成行向量,或者将行向量变成列向量的操作。

具体来说:

  • 如果 β 是一个列向量,例如:

    β=(β1β2βn)

    那么 βT 就是一个行向量: βT=(β1,β2,,βn)

  • 如果 β 是一个行向量:β=(β1,β2,,βn)
    那么 βT 就是一个列向量:

    βT=(β1β2βn)

问题三:什么是凸函数,如何证明?
凸函数(Convex Function)是定义在一个凸集上的实值函数,其满足以下性质:如果在该集合中任意两点之间画一条直线,该直线上的函数值不低于连接两点的函数值。直观地说,凸函数的图形总是“向上弯曲”的。

几何含义:在凸函数图形中,连接两点 (x1,f(x1))(x2,f(x2)) 的线段不会位于函数图形的下方。

使用二阶导数条件(Hessian 矩阵)可以证明是否为凸函数(下面证明都是使用该方法),即满足以下所有条件即可证明是凸函数:

  • 适用于两次可微的函数。
  • 计算函数的 Hessian 矩阵 2f(x)
  • 证明 2f(x) 是正半定矩阵,即对于任意非零向量 z,都有:zT2f(x)z0

举例:证明 f(x)=x2 是凸函数

  • 前提: 函数 f(x)=x2R 上可微且二次可微。

  • 步骤 1:计算一阶导数和二阶导数
    f(x)=2xf(x)=2

  • 步骤 2:检查二阶导数
    f(x)=2>0xR由于二阶导数恒为正,Hessian 矩阵在此情况下是 2,是正半定的。因此,f(x)=x2 是凸函数。

问题四:Hessian 矩阵是什么?二阶导数就是Hessian 矩阵吗?
Hessian 矩阵的定义

Hessian 矩阵是一个多变量函数的二阶导数矩阵。具体来说,如果 多变量函数( 即x 是一个 n 维向量) f(x)是一个可微两次的函数,其 Hessian 矩阵定义为:

2f(x)=(2fx122fx1x22fx1xn2fx2x12fx222fx2xn2fxnx12fxnx22fxn2)

其中:

  • 2fxi2f 对第 i 个变量的二阶偏导数。
  • 2fxixjf 对第 i 个变量和第 j 个变量的混合二阶偏导数。

二阶导数和 Hessian 矩阵的关系

  • 单变量函数:对于单变量函数 f(x),Hessian 矩阵退化为一个 1×1 的矩阵,也就是二阶导数 f(x) 本身。它描述了函数的曲率。

    • 例如,若 f(x)=x2,则其二阶导数为 f(x)=2,对应的 Hessian 矩阵是 [2]
  • 多变量函数:对于多变量函数 f(x),Hessian 矩阵的维度为 n×n,它包含了所有变量之间的二阶变化率。

    • 例如,若 f(x,y)=x2+y2,则其 Hessian 矩阵为:

      2f(x,y)=(2fx22fxy2fyx2fy2)=(2002)

证明思路及其过程

证明思路

  1. Logistic函数 f(x;β)β 的非凸性:
    Logistic函数 f(x;β) 对参数 β 不是凸的,是因为其二阶导数在某些情况下会出现负值。
    通过计算 f(x;β)β 的Hessian矩阵(即二阶导数矩阵),并证明它不是正半定的,从而证明其非凸性。

  2. 对数似然函数 L(β) 的凸性:
    对数似然函数 L(β) 常用于Logistic回归中拟合参数 β
    证明凸性即证明 L(β)β 的二阶导数(即Hessian矩阵)是正半定的。

证明过程
image

证明 Logistic 函数 f(x;β)=11+eβTx 对参数 β 的非凸性

第一步:计算 Logistic 函数的一阶导数

Logistic函数可以重写为:

f(x;β)=11+eβTx

对参数 β 求一阶导数:

f(x;β)β=eβTxx(1+eβTx)2

第二步:计算 Logistic 函数的二阶导数(Hessian矩阵)

接下来,对 f(x;β)β 再对 β 求导,得到Hessian矩阵:

2f(x;β)β2=eβTxxxT(1eβTx)(1+eβTx)3

第三步:验证Hessian矩阵的正负定性

为了判断函数的凸性,我们需要验证Hessian矩阵的正负定性。可以观察到,当 eβTx 较大时(即 βTx 较小时),1eβTx 为负,因此Hessian矩阵的部分项可能为负值。

由于Hessian矩阵不是正半定的,说明Logistic函数对参数 β非凸的


证明 对数似然函数 L(β)=y0βTx+ln(1+eβTx) 对参数 β 的凸性

第一步:计算对数似然函数L(β)的一阶导数

L(β)β=y0x+eβTx1+eβTxx

第二步:计算对数似然函数的二阶导数(Hessian矩阵)
对一阶导数再求导,得到Hessian矩阵:

2L(β)β2=eβTx(1+eβTx)2xxT

第三步:验证Hessian矩阵的正定性

Hessian矩阵的结构中,eβTx(1+eβTx)2 始终为非负数,且 xxT 是正半定矩阵。
因此,Hessian矩阵是 正半定的,这表明对数似然函数 L(β) 对参数 β凸的

解惑

为什么二阶导数会出现xT

关键在于理解 Hessian 矩阵的计算原理及向量外积的来源。

背景知识:

  • 梯度 l(β):对多维参数 β 求导数,结果是一个向量。它代表的是损失函数相对于每个参数的变化率。
  • Hessian 矩阵 2l(β):对梯度 l(β) 再求一次导数,结果是一个矩阵。Hessian 矩阵捕捉的是损失函数相对于参数的二阶变化,即二阶偏导数。

以损失函数对数似然函数:l(β)=i=1m(yiβTxi+ln(1+eβTxi)) 为例

计算一阶导数 l(β)
对于第 i 项,l(β) 的推导过程如下:

  1. yiβTxi 求导

    (yiβTxi)=yixi

  2. ln(1+eβTxi) 求导

    • σ(βTxi)=eβTxi1+eβTxi,即 Sigmoid 函数。
    • (ln(1+eβTxi)) 的导数为:eβTxi1+eβTxixi=σ(βTxi)xi
  3. 综合一阶导数

    l(β)=i=1m(σ(βTxi)xiyixi)
    或者写作:
    l(β)=i=1m((σ(βTxi)yi)xi)

计算二阶导数 2l(β)(Hessian 矩阵):
计算 Hessian 矩阵,即对梯度向量 l(β)β 再求导数:

2l(β)=i=1m(σ(βTxi)xiyixi)

注意到,yixi 是一个关于 β 的线性项,它的导数为零,所以我们只需要对 σ(βTxi)xi 进行求导。

σ(βTxi)xi 求导:
这里是关键步骤:σ(βTxi) 是一个标量,而 xi 是一个向量。

为什么σ(βTxi) 是一个标量,而 xi 是一个向量?
因为

  • β 是一个向量,假设它的维度是 n,表示模型的参数:β=(β1,β2,,βn)T
  • xi 是一个向量,维度同样是 n,表示第 i 个样本的特征:xi=(xi1,xi2,,xin)T
  • 内积 βTxi 是一个标量,即:βTxi=β1xi1+β2xi2++βnxin

内积的结果是一个实数(标量),因为它是将向量 βxi 对应元素相乘并求和得到的。

点击查看后续更详细解释

2. Sigmoid 函数 σ(z)

  • Sigmoid 函数的定义为:

    σ(z)=11+ez

    其中,z 是一个标量输入,σ(z) 的输出也是一个标量,取值范围在 (0,1) 之间。

  • 在这个问题中,z=βTxi,所以:

    σ(βTxi)=11+eβTxi

    这是一个标量,因为它是对内积结果(标量)进行非线性变换。

3. 向量 xi 的性质:

  • xi 是一个向量,维度为 n,表示第 i 个样本的特征向量。它在 Sigmoid 函数之外起作用,在计算梯度时需要考虑它的分量。

4. 结合梯度和 Hessian 矩阵:

  • 当我们计算 梯度 l(β) 时,会对每个 βjβ 的第 j 个分量)求导。因为 σ(βTxi) 是标量,所以对它求导时,向量 xi 会参与进来,表示为:

    βσ(βTxi)=σ(βTxi)(1σ(βTxi))xi

  • 当进一步计算 Hessian 矩阵 2l(β) 时,由于梯度项中有 xi,因此会产生向量的外积 xixiT,这个外积结果是一个 n×n 的矩阵。

总结:

  • σ(βTxi) 是一个标量,因为它是对内积 βTxi 的结果(一个数值)进行 Sigmoid 变换。
  • xi 是一个向量,因为它表示样本的特征向量(有多个维度)。
  • 在计算梯度和 Hessian 时,由于 xi 参与到对 β 的导数中,会导致出现 xixiT 这样的外积形式。

我们应用乘积法则来对这个项求导:

(σ(βTxi)xi)=σ(βTxi)xi+σ(βTxi)(xi)

  • 计算 σ(βTxi)xi

    • σ(βTxi)β 的导数为:

      σ(βTxi)=σ(βTxi)(1σ(βTxi))xi

      这里使用了 Sigmoid 函数的导数性质。

    • 将其代入,我们得到:

      σ(βTxi)xi=σ(βTxi)(1σ(βTxi))xixiT

      • 表示外积运算,xixiT 是外积的结果。
      • 由于 σ(βTxi) 是一个 n-维向量,外积的结果 xixiT 是一个 n×n 的矩阵。
  • 计算 σ(βTxi)(xi)

    • 这里 xi 是相对于 β 的常数向量,因此它的导数为零。这个部分对 Hessian 没有贡献。
  • 综合二阶导数

    最终得到 Hessian 矩阵:

    2l(β)=i=1mσ(βTxi)(1σ(βTxi))xixiT

为什么会出现 xixiT

  • 在对梯度求二阶导数时,由于 Sigmoid 函数的导数是一个标量,但我们同时需要对向量 xi 的分量进行导数运算,这就产生了向量之间的外积 xixiT
  • 外积 xixiT 的出现 是因为当我们在多维空间中对标量函数进行二次导数时,需要考虑所有分量之间的相互关系。

总结:

  • xixiT 出现在 Hessian 矩阵中是因为我们在求梯度的导数时,涉及到向量对自身的导数操作。
  • 导数结果中的 xixiT 表示向量的外积,它是一个矩阵,捕捉了 xi 的每个分量对参数 β 的二阶导数的变化。

实现对率回归

用python实现对率回归,对西瓜数据集(见表)进行划分。请附程序(含注释)
image

实现感知机算法

用 python 写一个感知机算法,分别实现 and 和 or。请附程序 (含注释)及运行结果

证明点到超平面的距离

image

计算梯度下降

image

work5

✨求解支持向量机(对偶形式)

image
或许参考:https://github.com/familyld/Machine_Learning/blob/master/06support_vector_machine.md

image

好的,以下是使用对偶方法求解支持向量机,划分三个数据点,并给出详细公式证明的过程。

1. 问题描述

给定三个数据点,其中正例为 x1=(3,3)x2=(3,2),反例为 x3=(1,1)。要求使用对偶方法求解支持向量机,并给出原问题、对偶形式和各自的约束条件,求解 (w,b),并给出支持向量。

2. 原问题

支持向量机的原问题可以表示为以下优化问题:

minw,b12||w||2s.t.yi(wTxi+b)1,i=1,2,3

其中,w 是权重向量,b 是偏置,yi 是数据点 xi 的标签(正例为 +1,反例为 -1)。

3. 对偶问题

为了求解原问题,我们可以将其转换为对偶问题。首先,我们引入拉格朗日乘子 αi0,构造拉格朗日函数:

L(w,b,α)=12||w||2i=13αi[yi(wTxi+b)1]

然后,对偶问题可以表示为:

maxαminw,bL(w,b,α)s.t.αi0,i=1,2,3

4. 对偶问题求解

为了求解对偶问题,我们首先需要求解 minw,bL(w,b,α)

L(w,b,α) 分别求 wb 的偏导数,并令其等于 0,得到:

wL(w,b,α)=wi=13αiyixi=0bL(w,b,α)=i=13αiyi=0

由第一个式子可以得到:

w=i=13αiyixi

w 代入拉格朗日函数,并利用第二个式子,得到:

L(w,b,α)=12||i=13αiyixi||2i=13αi[yi((j=13αjyjxj)Txi+b)1]=12i=13j=13αiαjyiyjxiTxji=13j=13αiαjyiyjxiTxjbi=13αiyi+i=13αi=12i=13j=13αiαjyiyjxiTxj+i=13αi

因此,对偶问题可以进一步表示为:

maxα12i=13j=13αiαjyiyjxiTxj+i=13αis.t.αi0,i=1,2,3i=13αiyi=0

这是一个二次规划问题,可以使用现成的优化工具求解。

5. 求解 (w, b)

将数据点代入对偶问题,并利用 KKT 条件,可以求解得到 α1=α2=18α3=14

αi 代入 w 的表达式,得到:

w=i=13αiyixi=181[33]+181[32]141[11]=[1234]

根据 KKT 条件,对于支持向量,有 yi(wTxi+b)=1。选择支持向量 x1=(3,3),代入求解 b,得到:

b=y1wTx1=1[1234][33]=134

因此,得到 (w,b)=([1234],134)

6. 支持向量

根据 KKT 条件,αi>0 对应的数据点为支持向量。因此,支持向量为 x1=(3,3)x2=(3,2)

总结

通过以上步骤,我们使用对偶方法求解了支持向量机,并给出了详细的公式证明。最终得到的分类超平面为 12x1+34x2134=0,支持向量为 x1=(3,3)x2=(3,2)

1. 问题描述

你有三个数据点:正例 x1=(3,3)x2=(3,2),反例 x3=(1,1)。目标是找到一个最优超平面来划分这些点。

2. 原始问题

支持向量机的原始问题可以表示为:

最小化: 12||w||2+Ci=1nξi

约束条件:

yi(wxi+b)1ξi,i=1,2,...,n

ξi0,i=1,2,...,n

其中:

  • w 是超平面的法向量
  • b 是超平面的截距
  • ξi 是松弛变量,允许一些点落在间隔内或错误的一侧
  • C 是惩罚参数,控制对错误分类的惩罚程度

3. 对偶问题

为了更容易求解,我们将原始问题转换为其对偶问题:

最大化: i=1nαi12i=1nj=1nαiαjyiyj(xixj)

约束条件:

i=1nαiyi=0

0αiC,i=1,2,...,n

其中 αi 是拉格朗日乘子。

4. 求解对偶问题

在本例中,我们有三个数据点,因此 n=3。让我们将数据点和标签代入对偶问题:

最大化: α1+α2+α312[α12(x1x1)+α22(x2x2)+α32(x3x3)+2α1α2(x1x2)2α1α3(x1x3)2α2α3(x2x3)]

约束条件:

α1+α2α3=0

0αiC,i=1,2,3

计算点积:

x1x1=(3,3)(3,3)=18
x2x2=(3,2)(3,2)=13
x3x3=(1,1)(1,1)=2
x1x2=(3,3)(3,2)=15
x1x3=(3,3)(1,1)=6
x2x3=(3,2)(1,1)=5

将这些值代入目标函数并简化,然后使用约束条件求解 α1α2α3。这通常需要使用二次规划求解器。

5. 计算 w 和 b

一旦你找到了最优的 αi,你可以使用以下公式计算 wb

w=i=1nαiyixi

b=yiwxi (对于任何 αi>0)

注意: 由于没有提供惩罚参数 C 的值,我无法提供完整的数值解。你需要选择一个合适的 C 值,然后使用二次规划求解器来找到最优的 αi,进而计算 wb

希望这些步骤能帮助你理解如何使用对偶方法解决这个支持向量机问题。

本文作者:kingwzun

本文链接:https://www.cnblogs.com/kingwz/p/18461719

版权声明:本作品采用知识共享署名-非商业性使用-禁止演绎 2.5 中国大陆许可协议进行许可。

posted @   kingwzun  阅读(51)  评论(0编辑  收藏  举报
点击右上角即可分享
微信分享提示
评论
收藏
关注
推荐
深色
回顶
收起