立体角：定义与公式

立体角（solid angle）这个概念很明显就是平面角向三维的推广，在我长达二十年的学习生涯中我觉得我应该没有正式接触过这个概念，所以上学期在关于辐射的一个公式里看到才会觉得很懵，这一篇就简单总结一下我现在对立体角的认识吧。

立体角的概念

既然立体角是平面角向三维的推广，那么我们先回顾一下平面角的定义。在高一数学必修一以前，我们对角的认识应该还停留在两边夹一尖儿，用量角器能量角的大小的阶段，用量角器量出来的角单位是度。高一数学必修一对平面角的定义是单位圆上对应的一段弧长，以此为基础定义了三角函数，下面的图是维基百科上的：

这种定义可以简单记为，
$$\theta =某个圆上对应的弧长/这个的圆的半径$$

它的单位用$rad$表示，因为弧长与半径的量纲都是长度，所以平面角是一个无量纲的量。小学数学介绍的理解二维到三维最简单的方法就是把这张圆绕着它的某条直径旋转一周，扫过的部分就变成了一个三维的球，而这段弧就变成了球上的一个区域，圆的半径扫过的部分变成了圆，根据这个思想，平面角推广到立体角就很顺理成章了，用$\Omega$表示立体角，那么
$$\Omega =某个球面上对应的一块面积/球半径的平方$$

立体角的单位是球面度，记为$sr$，是单词steradian，因为面积除以半径的平方后单位约掉了，所以球面度也是一个无量纲的量，它就是角度$rad$的推广。

下面这张图也是维基百科上的，我们可以借此理解一下球面度的含义。假设我们站在半径为$r$的球的球心观察，不转动脖子可以看到球面上面积为$r^2$的范围，那么我们的视野范围就可以用球面度$\frac{r^2}{r^2}=1sr$来衡量；然后再把自己想象成一只猫头鹰，只要动一动脖子就能看清整个球面，那么我们现在的视野范围就是$\frac{4\pi r^2}{r^2}=4\pi sr$。

立体角的公式

按照立体角的定义，球坐标系中的立体角公式是最好求的。球坐标系中的面积微元是
$$dA=(r\sin \theta d\varphi )(rd\theta )=r^2\sin \theta d\theta d\varphi$$

所以
$$d\Omega =\frac{dA}{r^2}=\sin \theta d\theta d\varphi$$

当对应面积为$S$时，立体角等于
$$\Omega ={\iint }_{S}d\Omega ={\iint }_S\sin \theta d\theta d\varphi$$

下图也是维基百科的：

在一般情况下，我们需要根据立体角的定义式进行计算：
$$d\Omega =\frac{dA}{r^2}$$

这里的$A$是在球面上的投影面积。

例 毕奥-萨伐尔定理可以写成
$$\text{B}=\frac{I}{c}\nabla \Omega$$

证 加‘的${\text{x}}^{\prime }$代表磁场的源（比如小型环形电流）的位置，不加的$\text{x}$代表观察者的位置，加粗的量代表向量，毕奥-萨伐尔定理是
$$\text{B}=\frac{I}{c}{\oint }_C\frac{d{\text{l}}^{\prime }\times (\text{x}-{\text{x}}^{\prime })}{|\text{x}-{\text{x}}^{\prime }{|}^3}$$

所以我们只需要计算$\nabla \Omega$并把它与曲线积分比较即可。

根据定义
$$d\Omega =\frac{dA}{r^2}$$

这里$r$代表磁场的源与观察者之间的距离，
$$r^2=|\text{x}-{\text{x}}^{\prime }{|}^2$$

而$dA$代表（观察者向磁场的源的方向看过去的视野范围）在contour $C$围成的曲面上的投影面积，这个曲面的外法向为$\hat{n}$，因此
$$d\Omega =\frac{1}{|\text{x}-{\text{x}}^{\prime }{|}^2}\frac{\hat{n}\cdot (\text{x}-{\text{x}}^{\prime })}{|\text{x}-{\text{x}}^{\prime }|}dS$$

所以
$$\Omega =-{\int }_{S(C)}\frac{\hat{n}\cdot (\text{x}-{\text{x}}^{\prime })}{|\text{x}-{\text{x}}^{\prime }{|}^3}dS$$

$S(C)$是强调这是contour $C$围成的曲面，$-$是因为$\hat{n}$取的外法线。现在我们考虑一下$\hat{n}dS$，这是contour $C$围成的曲面，而$C$承载的是磁场的源，也就是电流，用${\text{l}}^{\prime }$表示电流方向，也就是$C$的方向，
$$\Omega =-{\oint }_C\frac{-(d\text{x}\times d{\text{l}}^{\prime })\cdot (\text{x}-{\text{x}}^{\prime })}{|\text{x}-{\text{x}}^{\prime }{|}^3}$$

根据引理1（放在最后证明），
$$\Omega =d\text{x}\cdot {\oint }_C\frac{d{\text{l}}^{\prime }\times (\text{x}-{\text{x}}^{\prime })}{|\text{x}-{\text{x}}^{\prime }{|}^3}$$

因此
$$\nabla \Omega ={\oint }_C\frac{d{\text{l}}^{\prime }\times (\text{x}-{\text{x}}^{\prime })}{|\text{x}-{\text{x}}^{\prime }{|}^3}$$

证毕。

引理1
$$d\text{x}\times d{\text{l}}^{\prime }\cdot (\text{x}-{\text{x}}^{\prime })=[d{\text{l}}^{\prime }\times (\text{x}-{\text{x}}^{\prime })]\cdot \text{x}$$

引入Eddington张量${ϵ}_{ijk}$，用Einstein求和约定改写左边的式子：
$$\begin{gathered} d\text{x}\times d{\text{l}}^{\prime }\cdot (\text{x}-{\text{x}}^{\prime })=({ϵ}_{ijk}d{x}_{j}d{l}_j^{\prime })(x_i-x_i^{\prime }) \\ =-{ϵ}_{ijk}d{l}_j^{\prime }(x_i-x_i^{\prime })d{x}_j \\ ={ϵ}_{jki}d{l}_j^{\prime }(x_i-x_i^{\prime })d{x}_j \\ =[d{\text{l}}^{\prime }\times (\text{x}-{\text{x}}^{\prime })]\cdot \text{x} \end{gathered}$$