数学杂谈：高维空间向量夹角小记

最新推荐文章于 2024-04-12 13:58:12 发布

Espresso Macchiato

最新推荐文章于 2024-04-12 13:58:12 发布

阅读量4.9k

点赞数 4

分类专栏：基础数学文章标签：算法几何学线性代数向量夹角正交

本文链接：https://blog.csdn.net/codename_cys/article/details/122799384

版权

基础数学专栏收录该内容

24 篇文章 5 订阅

订阅专栏

数学杂谈：高维空间向量夹角小记

1. 问题描述

故事起源于long long ago的时候看到的苏剑林某一篇博客当中提到了一个结论：

高维空间中两个随机向量大概率是相互正交的。

当时就对这个结论比较在意了，今天突然地又想到了这个问题，就想着趁着放假把这个结论验证一下。

显然，苏剑林在博客中提到的这个结论的表述是比较随意的，为了更好地定义我们的问题，我们将其细化一下：

对于n维中间的单位球面上的两个随机向量，他们之间的夹角 $\theta$ 在n取大值时，趋近于90度。

要证明这个结论感觉比较难，但是要演示说明一下这个结果确实相对容易的，我们的第一反应就是蒙特卡洛模拟一下，实际生成 $N$ 组 $n$ 维空间上的均匀单位向量，然后考察一下他们之间的角度分布。

但是，要做到这件事，我们首先需要生成一下n维空间当中的均匀分布的单位向量。

2. n维空间中的均匀向量

1. 2维以及3维空间中的特殊情况

首先，我们来考察一些简单的情况，即2维和3维空间当中的情况。

1. 2维空间中的均匀分布向量

二维空间中的均匀分布向量其实就是单位圆上面的均匀分布向量，因此，我们只需要给出一个 $0$ 到 $\pi$ 之间均匀分布的角度 $\phi$ 即可得到一个均匀分布的单位向量 $\vec{v} = (sin\theta, cos\theta)$ 。

3. 3维空间中的均匀分布向量

对于三维的情况，事实上我相信大部分读者只要熟悉坐标系变换的话就还是很轻易能够写出解答的。

我们极坐标系下的坐标如下：
$\left \{ \begin{aligned} x & = r\cdot sin\theta \cdot sin\phi \\ y & = r\cdot sin\theta \cdot cos\phi \\ z & = r\cdot cos\theta \end{aligned} \right.$

进而我们可以得到，单位体积元的表达公式为：
$\begin{aligned} \rho & = dxdydz \\ & = r^2sin\theta drd\theta d\phi \end{aligned}$

可以看到，对于单位面元而言，其具体的表达式为 $\rho = C \cdot sin\theta d\theta d\phi = C' \cdot dcos\theta \cdot d\phi$ 。因此，要生成一个均匀分布，我们只需要按照 $cos\theta$ 的分布生成一个 $\theta$ ，然后生成一个 $0$ 到 $2\pi$ 上面均匀分布的 $\phi$ 即可。

给出具体的python实现如下：

import numpy as np

def dummy():
    theta = np.arccos(np.random.uniform(-1, 1))
    phi = np.random.uniform() * 2 * np.pi

    x = np.sin(theta) * np.sin(phi)
    y = np.sin(theta) * np.cos(phi)
    z = np.cos(theta)

    return (x, y, z)

2. n维坐标系中的均匀向量

现在，我们来考察n维空间中的情况。

我们仿照3维空间的情况，只要先给出体积元的极坐标表达式，然后考察其中空间角的表达式即可。

给出n维空间下的极坐标转换如下：

$\left\{ \begin{aligned} & x_1 = r \cdot cos\theta_1 \\ & x_2 = r \cdot sin\theta_1 \cdot cos\theta_2 \\ & x_3 = r \cdot sin\theta_1 \cdot sin\theta_2 \cdot cos\theta_3 \\ & ... \\ & x_{n-1} = r \cdot sin\theta_1 \cdot sin\theta_2 \cdot ... \cdot sin\theta_{n-2} \cdot cos\theta_{n-1} \\ & x_n = r \cdot sin\theta_1 \cdot sin\theta_2 \cdot ... \cdot sin\theta_{n-2} \cdot sin\theta_{n-1} \end{aligned} \right.$

可以得到n维空间上的体积元：
$\begin{aligned} dx_1dx_2...dx_n & = \frac{\partial(x_1, x_2, ..., x_n)}{\partial(r, \theta_1, \theta_2, ..., \theta_{n-1}))} \cdot drd\theta_1d\theta_2...d\theta_{n-1} \\ \\ & = det \begin{vmatrix} \frac{\partial x_1}{\partial r} & \frac{\partial x_1}{\partial \theta_1} & ... & \frac{\partial x_1}{\partial \theta_{n-1}} \\ \frac{\partial x_2}{\partial r} & \frac{\partial x_2}{\partial \theta_1} & ... & \frac{\partial x_2}{\partial \theta_{n-1}} \\ ... \\ \frac{\partial x_n}{\partial r} & \frac{\partial x_n}{\partial \theta_1} & ... & \frac{\partial x_n}{\partial \theta_{n-1}} \end{vmatrix} \cdot drd\theta_1d\theta_2...d\theta_{n-1} \\ \\ & = r^{n-1}sin^{n-2}\theta_1sin^{n-3}\theta_2...sin^2\theta_{n-3}sin\theta_{n-2} \cdot drd\theta_1d\theta_2...d\theta_{n-1} \end{aligned}$

由此，我们只需要令任意的 $\theta_i$ 满足分布条件 $sin^{n-1-i} \theta_i d\theta_{i}$ 是均匀分布，即可得到全空间范围内的一个针对空间角均匀分布的随机向量。

当然，这件事并不容易实现就是了。

当然，如果可以令 $x_i$ 是 $(-\infty, \infty)$ 范围内的均匀分布，那么事实上随机生成的向量也是在空间角上均匀的，但是这显然也同样难以实现。

3. 正态分布的巧妙应用

这里，我们给出一个黑科技，即虽然我们无法在 $(-\infty, \infty)$ 范围内生成一个均匀随机分布，但是我们可以退而求其次，如果对于一个n维向量，他的每一个维度上的值都满足正态分布 $N (0, 1)$ ，那么这样随机生成的向量在任意的n维空间角上也是均匀分布的。

我们考察其在任意n维空间体积元内的概率密度如下：

$\begin{aligned} \rho & = \Pi_{i=1}^{n} \frac{1}{\sqrt{2\pi}} e^{-\frac{x_i^2}{2}} \\ & = (\frac{1}{\sqrt{2\pi}})^n \cdot e^{\sum_{i=1}^{n} x_i^2 / 2} \\ & = (\frac{1}{\sqrt{2\pi}})^n \cdot e^{r^2 / 2} \\ \end{aligned}$