1143 字
6 分钟
类条件概率密度是高斯分布时的贝叶斯分类器
2026-06-16
无标签

考虑特征向量 x∈Rd\mathbf{x} \in \mathbb{R}^d,假设类条件概率密度是高斯分布:

p(x∣ωi)=1(2π)d/2∣Σi∣1/2exp⁡(−12(x−μi)TΣi−1(x−μi))∼N(μi,Σi)p(\mathbf{x}|\omega_i) = \frac{1}{(2\pi)^{d/2} |\Sigma_i|^{1/2}} \exp\left( -\frac{1}{2} (\mathbf{x} - \mu_i)^T \Sigma_i^{-1} (\mathbf{x} - \mu_i) \right)\sim \mathcal{N}(\mu_i, \Sigma_i)

此时的判别函数为:

TIP

贝叶斯判别准则就是最小错误率准则,选择后验概率最大的类别。

gi(x)=ln⁡p(x∣ωi)+ln⁡P(ωi)=−12(x−μi)TΣi−1(x−μi)−12ln⁡∣Σi∣+ln⁡P(ωi)−d2ln⁡(2π)i=arg⁡max⁡igi(x)\begin{aligned} g_i(\mathbf{x}) &= \ln p(\mathbf{x}|\omega_i) + \ln P(\omega_i) \\ &= -\frac{1}{2} (\mathbf{x} - \mu_i)^T \Sigma_i^{-1} (\mathbf{x} - \mu_i) - \frac{1}{2} \ln |\Sigma_i| + \ln P(\omega_i) - \cancel{\frac{d}{2} \ln (2\pi)}\\[1em] i &= \arg\max_i g_i(\mathbf{x}) \end{aligned}

情况1:Σi=σ2I,P(ωi)=1c\Sigma_i = \sigma^2 I,P(\omega_i)=\frac{1}{c}#

假设每一类数据在特征空间中都长成一个正圆球体(二维是圆,三维是球),并且所有类别的球体大小(半径)完全一样。

gi(x)=−12σ2∥x−μi∥2+constant∝−∥x−μi∥2i=arg⁡max⁡igi(x)  ⟹  i=arg⁡min⁡i∥x−μi∥2g_i(\mathbf{x}) = -\frac{1}{2\sigma^2} \|\mathbf{x} - \mu_i\|^2 + \text{constant}\propto -\|\mathbf{x} - \mu_i\|^2\\[1em] i = \arg\max_i g_i(\mathbf{x}) \implies i = \arg\min_i \|\mathbf{x} - \mu_i\|^2

此分类器称为距离分类器,判别函数可以用待识模式 x\mathbf{x} 与类别均值 μi\mu_i 之间的距离表示 gi(x)=−d(x,μi)g_i(\mathbf{x}) = -d(\mathbf{x}, \mu_i),最大化判别函数等于最小化距离,因此也称为最小距离分类器。

决策边界:两个球心连线的垂直平分线(超平面)。

情况2:Σi=Σ\Sigma_i = \Sigma#

假设每一类数据都长成椭球体,而且所有类别的椭球体形状、大小、朝向完全一致(比如都是东北-西南方向拉长的椭球),只是它们中心点(均值μ的位置不同。

gi(x)=−12(x−μi)TΣ−1(x−μi)+ln⁡P(ωi)+constant∝−12(x−μi)TΣ−1(x−μi)+ln⁡P(ωi)=−12xTΣ−1x+μiTΣ−1x−12μiTΣ−1μi+ln⁡P(ωi)\begin{aligned} g_i(\mathbf{x}) &= -\frac{1}{2} (\mathbf{x} - \mu_i)^T \Sigma^{-1} (\mathbf{x} - \mu_i) + \ln P(\omega_i) + \text{constant}\\ &\propto -\frac{1}{2} (\mathbf{x} - \mu_i)^T \Sigma^{-1} (\mathbf{x} - \mu_i) + \ln P(\omega_i)\\ &= -\frac{1}{2} \mathbf{x}^T \Sigma^{-1} \mathbf{x} + \mu_i^T \Sigma^{-1} \mathbf{x} - \frac{1}{2} \mu_i^T \Sigma^{-1} \mu_i + \ln P(\omega_i) \\ \end{aligned}

里面的纯二次 −12xTΣ−1x-\frac{1}{2} \mathbf{x}^T \Sigma^{-1} \mathbf{x} 对于所有类别相同,可以视为常数的一部分忽略。进一步展开:

gi(x)=μiTΣ−1x−12μiTΣ−1μi+ln⁡P(ωi)\begin{aligned} g_i(\mathbf{x}) &= \mu_i^T \Sigma^{-1} \mathbf{x} - \frac{1}{2} \mu_i^T \Sigma^{-1} \mu_i + \ln P(\omega_i) \end{aligned}

简化成线性形式,这就是 线性判别分析(LDA)。

gi(x)=wiTx+wi0g_i(\mathbf{x}) = \mathbf{w}_i^T \mathbf{x} + w_{i0}

决策边界:依旧是超平面,但不再是两个球心连线的垂直平分线了,而是根据 Σ\Sigma 的形状调整过的超平面。

情况3:Σi\Sigma_i 不同#

彻底放飞。每个类别的数据可以有自己的形状、大小和朝向。A类可以是细长的椭球,B类可以是扁平的圆盘,旋转的角度也可以完全不同。

和上面的情况2一样,纯二次项 −12xTΣi−1x-\frac{1}{2} \mathbf{x}^T \Sigma_i^{-1} \mathbf{x} 不能视为常数了,因为 Σi\Sigma_i 不同了。忽略常数展开后:

gi(x)=−12xTΣi−1x+μiTΣi−1x−12μiTΣi−1μi−12ln⁡∣Σi∣+ln⁡P(ωi)g_i(\mathbf{x}) = -\frac{1}{2} \mathbf{x}^T \Sigma_i^{-1} \mathbf{x} + \mu_i^T \Sigma_i^{-1} \mathbf{x} - \frac{1}{2} \mu_i^T \Sigma_i^{-1} \mu_i - \frac{1}{2} \ln |\Sigma_i| + \ln P(\omega_i)

决策边界:不再是超平面了,而是二次曲面(比如椭圆、双曲线等)。因此也称为二次判别分析(QDA)。

问题#

QDA 每类需要估计一个 d×dd \times d 的协方差矩阵 Σi\Sigma_i,即 d(d+1)2\frac{d(d+1)}{2} 个参数。当 dd 较大而样本数不足时,Σ^i\hat{\Sigma}_i 可能奇异或不稳定。

思路:减少需要估计的参数数量。一个极端做法是假设各维特征条件独立——协方差矩阵退化为对角矩阵,每类只需 dd 个方差参数。

朴素贝叶斯分类器#

TIP

参数估计里面的独立性假设是针对数据集中的特征向量之间而言的,而朴素贝叶斯分类器的独立性假设是针对具体特征向量的各个维度之间而言的。

假设特征向量的各个维度之间条件独立,即 p(x∣ωi)=∏j=1dp(xj∣ωi)p(\mathbf{x}|\omega_i) = \prod_{j=1}^d p(x_j|\omega_i) (第i类,第j维),每个特征单独建模为一维高斯分布:

p(xj∣ωi)=12πσijexp⁡(−(xj−μij)22σij2)p(x_j|\omega_i) = \frac{1}{\sqrt{2\pi} \sigma_{ij}} \exp\left( -\frac{(x_j - \mu_{ij})^2}{2\sigma_{ij}^2} \right)

特征向量的类条件概率密度函数为:

p(x∣ωi)=∏j=1dp(xj∣ωi)=∏j=1d12πσijexp⁡(−(xj−μij)22σij2)p(\mathbf{x}|\omega_i) = \prod_{j=1}^d p(x_j|\omega_i) = \prod_{j=1}^d \frac{1}{\sqrt{2\pi} \sigma_{ij}} \exp\left( -\frac{(x_j - \mu_{ij})^2}{2\sigma_{ij}^2} \right)

判别函数为:

gi(x)=ln⁡p(x∣ωi)+ln⁡P(ωi)=∑j=1d(−(xj−μij)22σij2−12ln⁡(2π)−ln⁡σij)+ln⁡P(ωi)=−12∑j=1d(xj−μij)2σij2−∑j=1dln⁡σij+ln⁡P(ωi)−d2ln⁡(2π)\begin{aligned} g_i(\mathbf{x}) &= \ln p(\mathbf{x}|\omega_i) + \ln P(\omega_i) \\ &= \sum_{j=1}^d \left( -\frac{(x_j - \mu_{ij})^2}{2\sigma_{ij}^2} - \frac{1}{2} \ln (2\pi) - \ln \sigma_{ij} \right) + \ln P(\omega_i) \\ &= -\frac{1}{2} \sum_{j=1}^d \frac{(x_j - \mu_{ij})^2}{\sigma_{ij}^2} - \sum_{j=1}^d \ln \sigma_{ij} + \ln P(\omega_i)- \cancel{\frac{d}{2} \ln (2\pi)} \end{aligned}

小结:协方差假设 → 分类器对照#

高斯判别分析的不同假设直接决定了决策面的形状和可训练参数的数量:

假设决策面每类参数数量分类器名称
Σi=σ2I\Sigma_i = \sigma^2 I线性(超平面)d+1d+1最小距离分类器
Σi=Σ\Sigma_i = \Sigma线性(超平面)d(d+3)2\frac{d(d+3)}{2}线性判别分析(LDA)
Σi\Sigma_i 任意二次(超二次曲面)d(d+3)2+d(d+1)2\frac{d(d+3)}{2} + \frac{d(d+1)}{2}二次判别分析(QDA)
Σi\Sigma_i 对角(Naive Bayes)二次(轴线与坐标轴平行)2d2d朴素贝叶斯

趋势:对协方差限制越强,参数越少、越不容易过拟合,但模型表达能力也越低。从 Σi=σ2I\Sigma_i = \sigma^2 I 一路放宽到任意 Σi\Sigma_i,就是在 偏差与方差之间做权衡。

类条件概率密度是高斯分布时的贝叶斯分类器
https://biscuit0613.github.io/posts/dl-basics/bayes-gauss/
作者
Biscuit
发布于
2026-06-16
许可协议
CC BY-NC-SA 4.0
概率密度估计:参数法与非参数法
Bayes判别准则