贡献者: JierPeter; Giacomo
三维欧几里得空间是我们最为熟悉的空间,其性质也非常好,易于研究。本节讨论的是三维空间中的 $C^2$ 参数曲线 $f: I = (a, b) \to \mathbb{R}^3$ 的性质,介绍了 Frenet 向量组以及该组中向量之间的关系,并引入了曲率、挠率等概念。
1. Frenet 向量组
图 1:Frenet 向量组示意图。途中三个蓝色向量就是红色曲线在某一点处的 Frenet 向量组,分别是单位切向量 $T$、单位主法向量 $N$ 和单位副法向量 $B$。
Frenet 向量组是三维欧几里得空间 $\mathbb{R}^3$ 中特有的概念,是刻画可微曲线性质的优良工具。
Frenet 向量组中的第一个向量,是单位切向量。关于 $t$ 对曲线 $f(t)$ 求导就能得到切向量,如果把 $t$ 看作时间,那么切向量 $\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }f(t)$ 也可以看成是曲线上一点的切速度。单位切向量就是方向与切速度相同,但长度为 $1$ 的向量。
为了得到单位切向量,我们可以用切速度除以切速度的大小来获得:$\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }f(t)/ \left\lvert \frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }f(t) \right\rvert $。但更简单的方法是,令切速度的大小恒为 $1$,或者说只讨论切速度大小恒为 $1$ 的情况,即弧长参数曲线定义 4 ;我们将用字母 $s$ 来表示弧长参数。
定义 1 单位切向量
使用弧长参数来描述可微曲线 $f(s)$,记 $T(s)=\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }f(s)$,称为曲线在参数为 $s$ 处的单位切向量(unit tangent vector)。
Frenet 向量组中的第二个向量,是单位主法向量。对于处处可微的曲线,切向量总是存在的,但是主法向量不一定存在。主法向量的定义如下:
定义 2 主法向量
使用弧长参数来描述可微曲线 $f(s)$,其单位切向量为 $T(s)$。如果 $\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }T(s)$ 不为零,则称 $\kappa(s) N(s)=\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }T(s)$ 为曲线在 $s$ 处的主法向量(normal vector)。其中 $N(s)$ 是单位向量,$\kappa$ 是一个正数,称作曲线在这一点处的曲率(curvature)。
由于使用弧长参数以后,单位切向量的大小不变,因此 $N(s)$ 恒与 $T(s)$ 垂直。$N(s)$ 就是那第二个成员,单位主法向量。
Frenet 向量组中的最后一个成员,是单位副法向量,是用前两个向量叉乘而来的。
定义 3 副法向量
条件设定如定义 1 和定义 2 ,称 $B(s)=T(s)\times N(s)$ 为曲线在 $s$ 处的单位副法向量(unit binormal vector)。
以上定义的副法向量必然是单位向量,因为 $T$ 和 $N$ 都是单位向量且正交。
Frenet 向量组可以用于构成所谓的Frenet 坐标架,也称活动坐标架。相比预先规定的自然坐标系,活动坐标架忽略了曲线在空间中的具体位置,更简洁和方便地表示了曲线本身的局部性质。
定义 4 密切平面
条件设定如定义 1 和定义 2 ,由 $T(s)$ 和 $N(s)$ 两向量在点 $f(s)$ 处所张成的二维空间,称为曲线在这一点处的密切平面(osculating plane)。
对于平面曲线,即整体都在一个平面上的曲线,其所在的平面就是每个点的密切平面。
2. 曲率和扭率
曲率已经在定义 2 中定义清楚了,是表征单位切向量旋转速率的数字。要说明的是,如果单位切向量在某点处不变化,那么主法向量的概念就不复存在,也就没有 Frenet 向量组的后两位成员了。在这种情况下,我们说曲线在这一点处的曲率为 $0$。更一致的曲率定义方法如下:
定义 5 曲率
使用弧长参数来描述可微曲线 $f(s)$,其单位切向量为 $T(s)$。称 $\kappa(s)= \left\lvert \frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }T(s) \right\rvert $ 为曲线在 $s$ 处的曲率(curvature)。
当 Frenet 向量组存在时,还有一个重要的概念叫扭率或挠率,在有的文献中也会被称为第二曲率。
在定义扭率之前,我们要先讨论副法向量的一个性质:$\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }B(s)$ 和 $N(s)$ 平行。首先,因为单位副法向量的长度不变,故有 $\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }B(s)$ 和 $B(s)$ 垂直。由于 $B=T\times N$,故
\begin{equation}
\begin{aligned}
\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }B(s)&=[\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }T(s)]\times N(s)+T(s)\times\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }B(s)\\&=T(s)\times\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }B(s)~,
\end{aligned}
\end{equation}
进而可知 $\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }B(s)$ 垂直于 $T(s)$。因此,$\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }B(s)$ 垂直于 $T$ 和 $B$,必然和 $N$ 在一条线上了。
定义 6 扭率
使用弧长参数来描述可微曲线 $f(s)$,令其 Frenet 向量组为 $T(s), N(s), B(S)$。令实数 $\tau$ 满足 $\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }B(s)=-\tau N(s)$,那么称实数 $\tau$ 为曲线在 $s$ 处的扭率(torsion)。
曲率描述了曲线单位切向量方向变化的速率,越快则曲线弯曲程度越强。扭率则描述了曲线偏离平面曲线的程度,也就是其扭曲程度;也可以理解为,扭率描述了密切平面旋转的速率。
3. Frenet-Serret 公式
三维欧几里得空间中的可微曲线 $f(s)$,如果处处有非零曲率 $\kappa(s)$,且其扭率为 $\tau(s)$(可以为零),则它一定有 Frenet 向量组。$T$,$N$ 和 $B$ 具有如下关系:
定理 1 Frenet-Serret 公式
对于曲线的 Frenet 向量组 $\{T(s), N(s), B(s)\}$,我们有如下关系:
\begin{equation}
\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }T(s)=\kappa N(s)~,
\end{equation}
\begin{equation}
\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }N(s)=-\kappa(s) T(s)+\tau(s) B(s)~,
\end{equation}
\begin{equation}
\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }B(s)=-\tau(s) N(s)~.
\end{equation}
用矩阵的形式可以更简洁地表示为
\begin{equation}
\begin{pmatrix}T'\\N'\\B'\end{pmatrix} = \begin{pmatrix}0&\kappa&0\\-\kappa&0&\tau\\0&-\tau&0\end{pmatrix} \begin{pmatrix}T\\N\\B\end{pmatrix} ~.
\end{equation}
矩阵表示也更好记忆。
证明:
- 根据 $N$ 的定义,直接可得式 2 。
- 根据扭率的定义,直接可得式 4 。
- 考虑到三个向量的方向关系,以及它们都是单位向量的事实,我们可以得到 $N=B\times T$。因此,$\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{s} }N=\frac{ \,\mathrm{d}{B} }{ \,\mathrm{d}{s} }\times T+B\times\frac{ \,\mathrm{d}{T} }{ \,\mathrm{d}{s} }=-\tau N\times T+B\times\kappa N=\tau B-\kappa T$。
证毕。
致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者
热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元,我们一周就能脱离亏损, 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。