散点图(散点矩阵)相关介绍

发布于:2025-07-29 ⋅ 阅读:(25) ⋅ 点赞:(0)

前言

提醒:
文章内容为方便作者自己后日复习与查阅而进行的书写与发布,其中引用内容都会使用链接表明出处(如有侵权问题,请及时联系)。
其中内容多为一次书写,缺少检查与订正,如有问题或其他拓展及意见建议,欢迎评论区讨论交流。

内容由AI辅助生成,仅经笔者审核整理,请甄别食用。


散点矩阵与高维函数可视化的数学原理

在高维数据分析中,散点矩阵(Scatter Matrix) 是连接数学理论与可视化实践的核心工具。它不仅是一个几何图形,更是多维统计分析的数学抽象。下面结合公式深入解析其原理。

一、散点矩阵的数学定义与性质

1. 数据矩阵与中心化

设高维函数采样数据为 X ∈ R n × p \mathbf{X} \in \mathbb{R}^{n \times p} XRn×p,其中:

  • n n n是样本数(如100个点)
  • p p p是维度(如5个变量)
  • x i j x_{ij} xij表示第 i i i个样本的第 j j j个特征值

中心化数据矩阵
X c = X − 1 n ⋅ X ˉ T \mathbf{X}_c = \mathbf{X} - \mathbf{1}_n \cdot \bar{\mathbf{X}}^T Xc=X1nXˉT
其中 X ˉ \bar{\mathbf{X}} Xˉ是各维度均值向量, 1 n \mathbf{1}_n 1n n × 1 n \times 1 n×1的全1向量。

2. 散点矩阵的数学表达

散点矩阵 S ∈ R p × p \mathbf{S} \in \mathbb{R}^{p \times p} SRp×p定义为:
S = X c T X c = ∑ i = 1 n ( x i − x ˉ ) ( x i − x ˉ ) T \mathbf{S} = \mathbf{X}_c^T \mathbf{X}_c = \sum_{i=1}^n (\mathbf{x}_i - \bar{\mathbf{x}})(\mathbf{x}_i - \bar{\mathbf{x}})^T S=XcTXc=i=1n(xixˉ)(xixˉ)T
其元素 s j k s_{jk} sjk表示维度 j j j k k k的交叉离差和:
s j k = ∑ i = 1 n ( x i j − x ˉ j ) ( x i k − x ˉ k ) s_{jk} = \sum_{i=1}^n (x_{ij} - \bar{x}_j)(x_{ik} - \bar{x}_k) sjk=i=1n(xijxˉj)(xikxˉk)

3. 散点矩阵的核心性质
  • 对称性 S = S T \mathbf{S} = \mathbf{S}^T S=ST
  • 半正定性:对任意非零向量 v \mathbf{v} v,有 v T S v ≥ 0 \mathbf{v}^T \mathbf{S} \mathbf{v} \geq 0 vTSv0
  • 与协方差矩阵的关系 Cov ( X ) = 1 n − 1 S \text{Cov}(\mathbf{X}) = \frac{1}{n-1} \mathbf{S} Cov(X)=n11S

二、散点矩阵在高维函数可视化中的应用

1. 可视化与矩阵元素的对应关系

散点矩阵图中的每个子图对应 S \mathbf{S} S的一个元素 s j k s_{jk} sjk

  • 对角线元素 s j j s_{jj} sjj:对应变量 j j j的直方图/密度图(反映方差)
  • 非对角线元素 s j k s_{jk} sjk:对应变量 j j j k k k的散点图(反映协方差)

例如,对于4维函数 f ( x 1 , x 2 , x 3 , x 4 ) f(x_1, x_2, x_3, x_4) f(x1,x2,x3,x4),其散点矩阵图为4×4网格,包含12个散点图和4个直方图。

2. 高维函数的局部线性近似

对于高维函数 y = f ( x ) y = f(\mathbf{x}) y=f(x),其在点 x 0 \mathbf{x}_0 x0附近的泰勒展开式为:
f ( x ) ≈ f ( x 0 ) + ∇ f ( x 0 ) T ( x − x 0 ) f(\mathbf{x}) \approx f(\mathbf{x}_0) + \nabla f(\mathbf{x}_0)^T (\mathbf{x} - \mathbf{x}_0) f(x)f(x0)+f(x0)T(xx0)
其中 ∇ f \nabla f f是梯度向量。通过散点矩阵分析各维度间的线性关系,可近似函数的局部行为。

3. 特征选择与降维

散点矩阵的特征值分解(EVD):
S = V Λ V T \mathbf{S} = \mathbf{V} \mathbf{\Lambda} \mathbf{V}^T S=VT
其中:

  • Λ = diag ( λ 1 , λ 2 , … , λ p ) \mathbf{\Lambda} = \text{diag}(\lambda_1, \lambda_2, \ldots, \lambda_p) Λ=diag(λ1,λ2,,λp)是特征值矩阵
  • V \mathbf{V} V的列向量是特征向量(对应主成分方向)

通过保留前 k k k个最大特征值对应的特征向量,可将 p p p维数据投影到 k k k维子空间(如PCA降维),并在散点矩阵图中可视化关键维度。

三、高维函数可视化的数学挑战

1. 维度灾难与信息损失

p p p很大时:

  • 散点矩阵图规模为 p × p p \times p p×p,信息过载
  • 直接可视化高维结构需依赖降维算法(如t-SNE、UMAP),但这些算法可能扭曲原始数据的几何性质
2. 非线性关系的捕捉

对于非线性函数 f ( x ) f(\mathbf{x}) f(x),散点矩阵只能反映线性相关性。需结合非线性变换(如核方法):
K i j = ϕ ( x i ) T ϕ ( x j ) \mathbf{K}_{ij} = \phi(\mathbf{x}_i)^T \phi(\mathbf{x}_j) Kij=ϕ(xi)Tϕ(xj)
其中 ϕ \phi ϕ是将数据映射到高维特征空间的非线性变换, K \mathbf{K} K是核矩阵。

四、总结

散点矩阵是高维函数可视化的数学基石,通过矩阵运算揭示变量间的线性关系和数据的几何结构。其核心价值在于:

  1. 将高维数据压缩为可解释的矩阵形式
  2. 为降维算法(PCA、t-SNE等)提供理论依据
  3. 通过可视化直观验证数学模型的合理性

理解散点矩阵的数学原理,能更深刻地把握高维数据的内在结构,避免被可视化结果误导。


网站公告

今日签到

点亮在社区的每一天
去签到