封面:https://www.pixiv.net/artworks/134297330
# 参数曲线从运动学的角度来看,曲线可以看作是一个质点在空间中的运动轨迹。参数曲线刻画了质点的位置随时间变化的规律。
# 曲线Definition. 映射 r : ( a , b ) → R n \pmb r:(a,b)\to \mathbb R^n r r : ( a , b ) → R n 称为 (正则)曲线 ,如果
r k ∈ C ∞ ( a , b ) r^k\in C^\infty (a,b) r k ∈ C ∞ ( a , b ) ;∣ r ′ ( t ) ∣ ≠ 0 |\pmb r'(t)|\ne 0 ∣ r r ′ ( t ) ∣ = 0 对任意 t ∈ ( a , b ) t\in (a,b) t ∈ ( a , b ) 成立。Remark. 正则曲线是光滑且切空间非退化的曲线。我们通常只考虑正则曲线。曲线的定义域只要求是开区间。
# 弧长与弧长参数Definition. 曲线 r : ( a , b ) → R n \pmb r:(a,b)\to \mathbb R^n r r : ( a , b ) → R n 的切向量 定义为
r ′ ( t ) = d r d t = ( d r 1 d t , d r 2 d t , ⋯ , d r n d t ) . \pmb r'(t)=\dfrac {\mathrm d\pmb r}{\mathrm dt}=\left(\dfrac{\mathrm dr^1}{\mathrm dt},\dfrac{\mathrm dr^2}{\mathrm dt},\cdots ,\dfrac{\mathrm dr^n}{\mathrm dt}\right). r r ′ ( t ) = d t d r r = ( d t d r 1 , d t d r 2 , ⋯ , d t d r n ) .
Definition. 曲线的弧长 定义为
s ( t ) = ∫ t 0 t ∣ r ′ ( u ) ∣ d u = ∫ t 0 t ⟨ r ′ ( u ) , r ′ ( u ) ⟩ 1 / 2 d u , s(t)=\int ^t_{t_0}|\pmb r'(u)|\mathrm du=\int^t_{t_0}\langle \pmb r'(u),\pmb r'(u)\rangle ^{1/2}\mathrm du, s ( t ) = ∫ t 0 t ∣ r r ′ ( u ) ∣ d u = ∫ t 0 t ⟨ r r ′ ( u ) , r r ′ ( u ) ⟩ 1 / 2 d u ,
其中 t 0 t_0 t 0 是定义域内的任一点。
Remark. 根据选点的不同,弧长函数族之间差一个常数,这个差异在求导时会被消除。此外,弧长是严格单增的,所以存在反函数,使得 t = t ( s ) t=t(s) t = t ( s ) 。
Lemma. s s s 称为曲线的弧长参数 。如果 s s s 本身作为参数时,记
r ˙ ( s ) = d r d s = d r d t d t d s \dot{\pmb r}(s)=\dfrac {\mathrm d\pmb r}{\mathrm ds}=\dfrac {\mathrm d\pmb r}{\mathrm dt}\dfrac {\mathrm dt}{\mathrm ds} r r ˙ ( s ) = d s d r r = d t d r r d s d t
则有 ∣ r ˙ ( s ) ∣ = 1 |\dot{\pmb r}(s)|=1 ∣ r r ˙ ( s ) ∣ = 1 。
Remark. 弧长参数表示下,曲线是匀速率运动的,切向量为单位,则加速度与速度是正交的。用 x ˙ \dot x x ˙ 来表示弧长参数下的导数。
Corollary. 在弧长参数下,下述由分量确定的切线方程就是 r ( s ) \pmb r(s) r r ( s ) 在 s = s 0 s=s_0 s = s 0 点的一阶近似 (基于 Taylor 展开的形式):
r ( s ) = r ( s 0 ) + r ˙ ( s 0 ) ( s − s 0 ) . \pmb r(s)=\pmb r(s_0)+\dot{\pmb r}(s_0)(s-s_0). r r ( s ) = r r ( s 0 ) + r r ˙ ( s 0 ) ( s − s 0 ) .
Remark. Taylor 展开指的是对分量上进行展开,那么总体形式上会保持 Taylor 展开形式。事实上,还有高阶近似,这里采取弧长参数展开:
r ( s ) = r ( s 0 ) + r ˙ ( s 0 ) ( s − s 0 ) + r ¨ ( s 0 ) 2 ! ( s − s 0 ) 2 + ⋯ \pmb r(s)=\pmb r(s_0)+\dot{\pmb r}(s_0)(s-s_0)+\dfrac {\ddot{\pmb r}(s_0)}{2!}(s-s_0)^2+\cdots r r ( s ) = r r ( s 0 ) + r r ˙ ( s 0 ) ( s − s 0 ) + 2 ! r r ¨ ( s 0 ) ( s − s 0 ) 2 + ⋯
# R2 上的曲线# Frenet 标架Definition. 曲线 r ( t ) \pmb r(t) r r ( t ) 的单位法向量 是 n ( t ) \pmb n(t) n n ( t ) ,满足 { t , n } \{\pmb t,\pmb n\} { t t , n n } 为正向正交基。即与标准正交基同定向(右手系)。
Notation. 记 t ( t ) \pmb t(t) t t ( t ) 是曲线 r ( t ) \pmb r(t) r r ( t ) 的单位切向量 ,即 t ( t ) = r ˙ \pmb t(t)=\dot {\pmb r} t t ( t ) = r r ˙ 。
Remark. 即 n ( t ) \pmb n(t) n n ( t ) 是 t ( t ) \pmb t(t) t t ( t ) 顺时针旋转 9 0 ∘ 90^\circ 9 0 ∘ 后的向量,写成坐标就是
n = ( 0 1 − 1 0 ) t . \pmb n=\begin{pmatrix}0&1\\ -1&0\end{pmatrix}\pmb t. n n = ( 0 − 1 1 0 ) t t .
Definition. 称 { r ; t , n } \{\pmb r;\pmb t,\pmb n\} { r r ; t t , n n } 为曲线 r ( t ) \pmb r(t) r r ( t ) 的 Frenet 标架 。
Theorem. 在弧长参数下,Frenet 标架满足 Frenet 方程 :
( t ˙ n ˙ ) = ( 0 κ − κ 0 ) ( t n ) \begin{pmatrix}\dot{\pmb t}\\ \dot{\pmb n}\end{pmatrix}=\begin{pmatrix}0&\kappa \\ -\kappa &0\end{pmatrix}\begin{pmatrix}\pmb t\\ \pmb n\end{pmatrix} ( t t ˙ n n ˙ ) = ( 0 − κ κ 0 ) ( t t n n )
其中 κ \kappa κ 称为曲线的曲率 ,由上式确定。
# 几何意义# 切线角度Remark. 曲率 κ \kappa κ 描述了曲线偏离切线的程度。因为
r ( s ) = r ( s 0 ) + t ( s 0 ) ( s − s 0 ) + κ ( s 0 ) n ( s 0 ) 2 ! ( s − s 0 ) 2 + ⋯ \pmb r(s)=\pmb r(s_0)+\pmb t(s_0)(s-s_0)+\dfrac {\kappa (s_0)\pmb n(s_0)}{2!}(s-s_0)^2+\cdots r r ( s ) = r r ( s 0 ) + t t ( s 0 ) ( s − s 0 ) + 2 ! κ ( s 0 ) n n ( s 0 ) ( s − s 0 ) 2 + ⋯
如果 κ ( s 0 ) > 0 \kappa(s_0)>0 κ ( s 0 ) > 0 ,则在该点处曲线向内(法方向)弯曲;κ ( s 0 ) < 0 \kappa(s_0)<0 κ ( s 0 ) < 0 时向外弯曲。
# 法线角度Remark. 从 Gauss 映射角度,曲率反映了法向量变化的快慢。因为
n ( s ) = n ( s 0 ) − κ ( s 0 ) t ( s 0 ) ( s − s 0 ) + ⋯ \pmb n(s)=\pmb n(s_0)-\kappa (s_0)\pmb t(s_0)(s-s_0)+\cdots n n ( s ) = n n ( s 0 ) − κ ( s 0 ) t t ( s 0 ) ( s − s 0 ) + ⋯
如果 κ ( s 0 ) > 0 \kappa(s_0)>0 κ ( s 0 ) > 0 ,则法向量向切线反方向(逆时针)转动;κ ( s 0 ) < 0 \kappa(s_0)<0 κ ( s 0 ) < 0 时向切线方向(顺时针)转动。
上面两种观点,当 κ ( s 0 ) = 0 \kappa(s_0)=0 κ ( s 0 ) = 0 时都不能确定曲线的弯曲性质。
# R3 上的曲线# Frenet 标架Definition. 曲线 r ( t ) \pmb r(t) r r ( t ) 的曲率向量 定义为 t ˙ ( s ) \dot{\pmb t}(s) t t ˙ ( s ) ,其模长称为曲线的曲率 κ ( s ) \kappa(s) κ ( s ) ;当 κ ( s ) > 0 \kappa(s)>0 κ ( s ) > 0 时,主法向量 和副法向量 定义为
n ( s ) = t ˙ ( s ) ∣ t ˙ ( s ) ∣ = t ˙ ( s ) κ ( s ) , b ( s ) = t ( s ) ∧ n ( s ) \pmb n(s)=\dfrac {\dot{\pmb t}(s)}{|\dot{\pmb t}(s)|}=\dfrac {\dot{\pmb t}(s)}{\kappa(s)},\quad \pmb b(s)=\pmb t(s)\wedge \pmb n(s) n n ( s ) = ∣ t t ˙ ( s ) ∣ t t ˙ ( s ) = κ ( s ) t t ˙ ( s ) , b b ( s ) = t t ( s ) ∧ n n ( s )
从而正交标架 { r ; t , n , b } \{\pmb r;\pmb t,\pmb n,\pmb b\} { r r ; t t , n n , b b } 称为空间曲线的 Frenet 标架 。
Remark. Frenet 标架的三根坐标轴分别是曲线的切线、主法线和副法线,三个坐标平面分别是切法平面、密切平面(以 b \pmb b b b 为法向量)和从切平面(以 n \pmb n n n 为法向量)。
Theorem. 在弧长参数下,空间曲线的 Frenet 标架满足 Frenet 方程 :
( t ˙ n ˙ b ˙ ) = ( 0 κ 0 − κ 0 τ 0 − τ 0 ) ( t n b ) \begin{pmatrix}\dot{\pmb t}\\ \dot{\pmb n}\\ \dot{\pmb b}\end{pmatrix}=\begin{pmatrix}0&\kappa &0\\ -\kappa &0&\tau \\ 0&-\tau &0\end{pmatrix}\begin{pmatrix}\pmb t\\ \pmb n\\ \pmb b\end{pmatrix} ⎝ ⎛ t t ˙ n n ˙ b b ˙ ⎠ ⎞ = ⎝ ⎛ 0 − κ 0 κ 0 − τ 0 τ 0 ⎠ ⎞ ⎝ ⎛ t t n n b b ⎠ ⎞
其中 τ \tau τ 称为曲线的挠率 ,由上式确定。
Proposition. 设空间曲线 r \pmb r r r 的曲率 κ > 0 \kappa>0 κ > 0 ,则 r \pmb r r r 是平面曲线,当且仅当 τ ≡ 0 \tau\equiv 0 τ ≡ 0 。
# 几何意义# 渐进展开角度Remark. 考虑在 s 0 = 0 s_0=0 s 0 = 0 处的 Taylor 展开:
r ( s ) = r ( s 0 ) + s r ˙ ( s 0 ) + s 2 2 ! r ¨ ( s 0 ) + s 3 3 ! r ¨ ˙ ( s 0 ) + o ( s 3 ) \pmb r(s)=\pmb r(s_0)+s\dot{\pmb r}(s_0)+\dfrac {s^2}{2!}\ddot{\pmb r}(s_0)+\dfrac {s^3}{3!}\dot{\ddot{\pmb r}}(s_0)+\pmb {o}(s^3) r r ( s ) = r r ( s 0 ) + s r r ˙ ( s 0 ) + 2 ! s 2 r r ¨ ( s 0 ) + 3 ! s 3 r r ¨ ˙ ( s 0 ) + o o ( s 3 )
所以
r ( s ) = r ( s 0 ) + ( s − κ 2 s 3 6 ) t + ( κ s 2 2 ! + κ ˙ s 3 3 ! ) n + κ τ s 3 3 ! b + o ( s 3 ) \pmb r(s)=\pmb r(s_0)+\left(s-\dfrac {\kappa^2s^3}6\right)\pmb t+\left(\dfrac {\kappa s^2}{2!}+\dfrac {\dot\kappa s^3}{3!}\right)\pmb n+\dfrac {\kappa\tau s^3}{3!}\pmb b+\pmb {o}(s^3) r r ( s ) = r r ( s 0 ) + ( s − 6 κ 2 s 3 ) t t + ( 2 ! κ s 2 + 3 ! κ ˙ s 3 ) n n + 3 ! κ τ s 3 b b + o o ( s 3 )
当 s → s 0 s\to s_0 s → s 0 时,曲线在该点处的三阶近似由切线、主法线和副法线的分量给出。由于 κ \kappa κ 的符号是确定的,所以 τ \tau τ 的符号决定了曲线穿越密切平面的方向。当 τ ( s 0 ) > 0 \tau(s_0)>0 τ ( s 0 ) > 0 时,曲线沿 b \pmb b b b 穿越 r ( s 0 ) \pmb r(s_0) r r ( s 0 ) 的密切平面;当 τ ( s 0 ) < 0 \tau(s_0)<0 τ ( s 0 ) < 0 时,曲线沿 b \pmb b b b 反向穿越密切平面。
# 投影角度Remark. 曲率 κ \kappa κ 描述了曲线偏离切线的程度。考虑上述的渐进展开,将曲线在密切平面上投影,得到的平面曲线的曲率即为 κ \kappa κ 。因为在密切平面的投影曲线满足
r ( s ) = r ( s 0 ) + s t ( s 0 ) + κ ( s 0 ) s 2 2 ! n ( s 0 ) + o ( s 2 ) \pmb r(s)=\pmb r(s_0)+s\pmb t(s_0)+\dfrac {\kappa (s_0)s^2}{2!}\pmb n(s_0)+\pmb o(s^2) r r ( s ) = r r ( s 0 ) + s t t ( s 0 ) + 2 ! κ ( s 0 ) s 2 n n ( s 0 ) + o o ( s 2 )
曲率 κ \kappa κ 描述了曲线在密切平面内的弯曲程度,和 b \pmb b b b 方向无关。
Remark. 挠率 τ \tau τ 描述了曲线偏离密切平面的程度。同样地,有在副法向量的投影:
⟨ r ( s ) − r ( s 0 ) , b ( s 0 ) ⟩ = κ ( s 0 ) τ ( s 0 ) s 3 3 ! + o ( s 3 ) \langle \pmb r(s)-\pmb r(s_0),\pmb b(s_0)\rangle =\dfrac {\kappa (s_0)\tau (s_0)s^3}{3!}+\pmb o(s^3) ⟨ r r ( s ) − r r ( s 0 ) , b b ( s 0 ) ⟩ = 3 ! κ ( s 0 ) τ ( s 0 ) s 3 + o o ( s 3 )
# 曲线论基本定理# ODE 存在唯一性定理Theorem. Picard-Lindelof/ Cauchy-Lipschitz 定理 指的是,考虑初值问题
d y d t = f ( t , y ) , y ( t 0 ) = y 0 \dfrac {\mathrm d\pmb y}{\mathrm dt}=f(t,\pmb y),\quad \pmb y(t_0)=\pmb y_0 d t d y y = f ( t , y y ) , y y ( t 0 ) = y y 0
其中 y ∈ R n , t ∈ R \pmb y\in\mathbb R^n,t\in\mathbb R y y ∈ R n , t ∈ R ,假设
f f f 在 R × R n \mathbb R\times \mathbb R^n R × R n 中 ( t 0 , y 0 ) (t_0,\pmb y_0) ( t 0 , y y 0 ) 某个闭域内连续;f f f 关于 y \pmb y y y 满足 Lipschitz 条件,即存在常数 L > 0 L>0 L > 0 ,使得对该闭域的任意 ( t , y 1 ) , ( t , y 2 ) (t,\pmb y_1),(t,\pmb y_2) ( t , y y 1 ) , ( t , y y 2 ) 有∥ f ( t , y 1 ) − f ( t , y 2 ) ∣ ≤ L ∣ y 1 − y 2 ∣ \|f(t,\pmb y_1)-f(t,\pmb y_2)|\le L|\pmb y_1-\pmb y_2| ∥ f ( t , y y 1 ) − f ( t , y y 2 ) ∣ ≤ L ∣ y y 1 − y y 2 ∣
则存在局部唯一解 y ( t ) \pmb y(t) y y ( t ) 。
Remark. 光滑函数是 Lipschitz 连续的。
# R3 上的曲线Proposition. 空间曲线的弧长、曲率和挠率在刚体运动下不变。
# 唯一性Theorem. 定义在相同弧长参数区间上的空间曲线,如果它们的曲率和挠率函数相同,则它们在刚体运动下是等价的。反之亦然。
# 存在性Theorem. 给定区间上的可微函数 κ ( s ) > 0 , τ ( s ) \kappa(s)>0,\tau(s) κ ( s ) > 0 , τ ( s ) ,则存在 R 3 \mathbb R^3 R 3 的弧长参数空间曲线 r ( s ) , s ∈ ( a , b ) \pmb r(s),s\in (a,b) r r ( s ) , s ∈ ( a , b ) ,它以 s s s 为弧长参数,以 κ ( s ) , τ ( s ) \kappa(s),\tau(s) κ ( s ) , τ ( s ) 分别为曲率和挠率。
证明 Proof. 考虑关于四个未知向量值函数 r ( s ) , t ( s ) , n ( s ) , b ( s ) \pmb r(s),\pmb t(s),\pmb n(s),\pmb b(s) r r ( s ) , t t ( s ) , n n ( s ) , b b ( s ) 的常微分方程组(Frenet 方程和运动方程),先证明常微分方程的唯一解 { t , n , b } \{\pmb t,\pmb n,\pmb b\} { t t , n n , b b } 是正交标架,再由 r ˙ = t \dot{\pmb r}=\pmb t r r ˙ = t t 得到 r \pmb r r r 。
Remark. 需要用到常微分方程解对初值的唯一性和存在性定理。
# 特殊曲线分类Proposition. 曲率或挠率为常数的曲线:
κ ≡ 0 \kappa\equiv 0 κ ≡ 0 ,当且仅当 r \pmb r r r 是直线;τ ≡ 0 \tau\equiv 0 τ ≡ 0 ,当且仅当 r \pmb r r r 是平面曲线;κ > 0 , τ = 0 \kappa>0,\tau=0 κ > 0 , τ = 0 为常数,当且仅当 r \pmb r r r 是半径为 1 / ∣ κ ∣ 1/|\kappa| 1 / ∣ κ ∣ 的圆周,或它的一部分;κ > 0 , τ ≠ 0 \kappa>0,\tau\neq 0 κ > 0 , τ = 0 为常数,当且仅当 r \pmb r r r 是圆柱螺旋线,或它的一部分。Remark. 圆柱螺旋线是绕某直线等速旋转并沿该直线等速平移的曲线,具体地
( a cos t , a sin t , b t ) , a > 0 : κ = a a 2 + b 2 , τ = b a 2 + b 2 (a\cos t,a\sin t,bt),a>0:\quad \kappa=\dfrac {a}{a^2+b^2},\quad \tau =\dfrac {b}{a^2+b^2} ( a cos t , a sin t , b t ) , a > 0 : κ = a 2 + b 2 a , τ = a 2 + b 2 b
如果 a < 0 a<0 a < 0 ,则为(反向)圆柱螺旋线。
# R2 上的曲线# 存在性和唯一性Theorem. 设 κ ( s ) \kappa(s) κ ( s ) 是连续可微函数,则
存在平面 R 2 \mathbb R^2 R 2 的曲线 r ( s ) \pmb r(s) r r ( s ) ,以 s s s 为弧长参数,以 κ ( s ) \kappa(s) κ ( s ) 为曲率; 定义在相同弧长参数区间上的平面曲线,如果它们的曲率相同,则它们在刚体运动下是等价的。反之亦然。 # 特殊曲线分类Proposition. 曲率为常数的曲线:
κ ≡ 0 \kappa\equiv 0 κ ≡ 0 ,当且仅当 r \pmb r r r 是直线;κ ≡ c ≠ 0 \kappa\equiv c\neq 0 κ ≡ c = 0 ,当且仅当 r \pmb r r r 是半径为 1 / ∣ c ∣ 1/|c| 1 / ∣ c ∣ 的圆周,或它的一部分。# 计算# 隐式曲线对于难以找到参数化的曲线,我们可以通过隐式方程来刻画曲线的局部性质。
# 基于隐函数定理Example. 设 F F F 是 R 2 \mathbb R^2 R 2 上的光滑函数,C = { ( x , y ) ∈ R 2 : F ( x , y ) = 0 } C=\{(x,y)\in\mathbb R^2:F(x,y)=0\} C = { ( x , y ) ∈ R 2 : F ( x , y ) = 0 } ,则 P = ( x 0 , y 0 ) ∈ C P=(x_0,y_0)\in C P = ( x 0 , y 0 ) ∈ C 如果满足 F y ( P ) ≠ 0 F_y(P)\neq 0 F y ( P ) = 0 ,则由隐函数定理知道 P P P 的一个邻域内,F ( x , y ) = 0 F(x,y)=0 F ( x , y ) = 0 有显式表示 y = f ( x ) y=f(x) y = f ( x ) ,所以 P P P 的某个邻域内 C C C 是一条曲线,有参数表示 ( t , f ( t ) ) (t,f(t)) ( t , f ( t ) ) ,满足
f ′ ( x ) = d y d x = − F x F y f'(x)=\dfrac {\mathrm d y}{\mathrm dx}=-\dfrac {F_x}{F_y} f ′ ( x ) = d x d y = − F y F x
# 基于水平集Example. 光滑函数 F ( x , y ) = 0 F(x,y)=0 F ( x , y ) = 0 的切线由方程
F x ( x 0 , y 0 ) ( x − x 0 ) + F y ( x 0 , y 0 ) ( y − y 0 ) = 0 F_x(x_0,y_0)(x-x_0)+F_y(x_0,y_0)(y-y_0)=0 F x ( x 0 , y 0 ) ( x − x 0 ) + F y ( x 0 , y 0 ) ( y − y 0 ) = 0
确定,其中 ( x 0 , y 0 ) (x_0,y_0) ( x 0 , y 0 ) 是曲线上的一点。这是因为考虑梯度
∇ F ( x 0 , y 0 ) = ( F x ( x 0 , y 0 ) , F y ( x 0 , y 0 ) ) \nabla F(x_0,y_0)=(F_x(x_0,y_0),F_y(x_0,y_0)) ∇ F ( x 0 , y 0 ) = ( F x ( x 0 , y 0 ) , F y ( x 0 , y 0 ) )
它垂直于水平集 F ( x , y ) = 0 F(x,y)=0 F ( x , y ) = 0 的切空间。
# 一般参数下曲率和挠率的计算Example. 曲线 r ( t ) = ( x ( t ) , y ( t ) ) \pmb r(t)=(x(t),y(t)) r r ( t ) = ( x ( t ) , y ( t ) ) 的曲率为
κ ( t ) = x ′ y ′ ′ − y ′ x ′ ′ ( x ′ 2 + y ′ 2 ) 3 / 2 \kappa(t)=\dfrac {x'y''-y'x''}{(x'^2+y'^2)^{3/2}} κ ( t ) = ( x ′ 2 + y ′ 2 ) 3 / 2 x ′ y ′ ′ − y ′ x ′ ′
Example. 曲线 r ( t ) = ( x ( t ) , y ( t ) , z ( t ) ) \pmb r(t)=(x(t),y(t),z(t)) r r ( t ) = ( x ( t ) , y ( t ) , z ( t ) ) 的曲率和挠率为
κ ( t ) = ∣ r ′ ∧ r ′ ′ ∣ ∣ r ′ ∣ 3 , τ ( t ) = ( r ′ , r ′ ′ , r ′ ′ ′ ) ∣ r ′ ∧ r ′ ′ ∣ 2 \kappa(t)=\dfrac {|\pmb r'\wedge \pmb r''|}{|\pmb r'|^3},\quad \tau(t)=\dfrac {(\pmb r', \pmb r'', \pmb r''')}{|\pmb r'\wedge \pmb r''|^2} κ ( t ) = ∣ r r ′ ∣ 3 ∣ r r ′ ∧ r r ′ ′ ∣ , τ ( t ) = ∣ r r ′ ∧ r r ′ ′ ∣ 2 ( r r ′ , r r ′ ′ , r r ′ ′ ′ )