微分几何笔记(2) —— 曲线的参数化
第二周讲完了Klingenberg的第一章Curves,做一点微小的笔记。
分成三个部分,本篇讲曲线的弧长参数;下一篇讲一般的Frenet标架及方程组;再下一篇讲二维三维空间曲线的curvature。
GTM51对入门者会难一些,因为直接从最一般的nnn维情况入手,再回头看二三维空间中的曲线,相比之下 Calculus and Analysis in Euclidean Space 这本UTM的Chapter 8 Parametrized Curves 算是相当新手友好的入门内容,由二三维曲线推广到一般情况。
借用宋老师的话,微分几何的本质是几何,主要通过微积分的手段,用数和函数来刻画几何图形的性质。可能在学习的时候会觉得明明很直观,几何画出来就是这样的,但是要用规范的严格语言说出来很困难。所以说也是一个打基础的课程,先学会怎样严格的表述,在这基础之上才能正确的使用微积分这个工具。
2.1 参数曲线(Parametrized Curves)
Definition 2.1.1 参数曲线 是指一个光滑映射(smooth mapping):c:I→Rnc: I \rightarrow \mathbb{R}^nc:I→Rn 这里n≥1n \geq 1n≥1, I⊆RnI\subseteq \mathbb{R}^nI⊆Rn 且非空。
更进一步,如果 ∀t∈I,c′(t)≠0\forall t\in I, c'(t)\neq0∀t∈I,c′(t)=0,称这样的曲线是正则的(regular).
也就是说我们平时说的“曲线”,其实是一个映射,平时说“点在曲线上”,是指这一点,如c(t)c(t)c(t)在这个映射的像集中。
我们一般都要求曲线是正则的,因为是光滑映射,所以导数不等于零⇔\Leftrightarrow⇔导数的范数不等于零,直观理解是这样的曲线是不会停下的,也因此我们可以用隐函数定理,反函数定理等。
如果III是不是开区间,定义的意思是,存在一个包含III的开区间I∗I^*I∗,使得c∗c^*c∗在I∗I^*I∗上是一个光滑映射,并且c∗∣I=cc^*|I=cc∗∣I=c,这样就解决了端点无定义的问题。
这里我们说的参数曲线都是“光滑”的,这个要求也太高了,不过现在考虑的都是很显然,一眼能够看明白的曲线,以后会逐步放宽。
Definition 2.1.2
i)切空间(Tangent Space):
对x0∈Rnx_0 \in \mathbb{R}^nx0∈Rn,切空间是指所有以x0x_0x0为起点的nnn维向量所构成的空间,记作Tx0RnT_{x_0} \mathbb{R}^nTx0Rn.
ii)沿着曲线ccc的向量场(Vector field along c:I→Rnc: I \rightarrow \mathbb{R}^nc:I→Rn):
指一个可微映射X:I→RnX: I\rightarrow \mathbb{R}^nX:I→Rn,∀t∈I\forall t\in I∀t∈I,X(t)∈Tc(t)RnX(t) \in T_{c(t)} \mathbb{R}^nX(t)∈Tc(t)Rn.
iii)切向量场(Tangent vector field of c:I→Rnc: I \rightarrow \mathbb{R}^nc:I→Rn):
指一个沿着ccc的向量场,其中在c(t)c(t)c(t)处的向量由切向量t↦c′(t)t \mapsto c'(t)t↦c′(t)给出。
这里主要是规定向量的起点位于什么地方。
切空间是相当于把x0x_0x0当成原点所建立的Rn\mathbb{R}^nRn空间,所有x0∈Rnx_0 \in \mathbb{R}^nx0∈Rn中向量起于x0x_0x0.
沿着曲线的向量场是指在曲线的每一点上都"生长"着一些向量,比如下图中所示那样:
图1:生长在曲线上的向量场
接下来一节讲讲弧长参数曲线(parametrization by arc length),用曲线的弧长做参数,这样做可以简化计算过程。
2.2 弧长参数曲线(Parametrization by Arc Length)
Definition 2.2.1 如果一个映射ϕ:I→I′\phi: I\rightarrow I'ϕ:I→I′是光滑的,而且其逆映射ϕ−1:I′→I\phi^{-1}: I'\rightarrow Iϕ−1:I′→I也是光滑的,称映射ϕ\phiϕ是微分同胚(Diffeomorphism)
Definition 2.2.2(Equivalence of Curves) 两条曲线 c:I→Rnc: I \rightarrow \mathbb{R}^nc:I→Rn 和 c~:I~→Rn\tilde{c} : \tilde{I} \rightarrow \mathbb{R}^nc~:I~→Rn 之间如果存在一个微分同胚 ϕ\phiϕ,使得 c~=c∘ϕ\tilde{c}=c\circ \phic~=c∘ϕ,则称ϕ\phiϕ是一个参数变换(parameter transformation).并且如果 ϕ′>0\phi'>0ϕ′>0,称这样的参数变换是保持定向的(orientation preserving) .
由参数变换所得的曲线构成一个所有 Rn\mathbb{R}^nRn中曲线的等价类(满足自反性(reflexive),对称性(symmetric)和传递性(transitive)),给这个等价类起个名字叫做 非参数化曲线(unparametrized curve).
其实我们可以把III看成是关于时间的集合,把c(t)c(t)c(t)看成一个空间中随时间运动的轨迹,这样c′(t)c'(t)c′(t)就是速度,而∣c′(t)∣|c'(t)|∣c′(t)∣则是速率。通过不同的参数变换,改变了质点沿曲线运动一定长度所需要的时间(区间III的长度),从而改变了质点运动的速率。我们特别关注的,是其中以单位速度运动的参数曲线,也就是∣c′(s)∣=1|c'(s)|=1∣c′(s)∣=1 的曲线。
Definition 2.2.3 (Arc length of a curve) 光滑曲线c:I→Rnc: I \rightarrow \mathbb{R}^nc:I→Rn上两点t0,t∈I,t0<tt_0,t\in I, t_0 <tt0,t∈I,t0<t之间的弧长为:
L(t0,t)=∫t0t∣c′(τ)∣dτL(t_0,t)=\int_{t_0}^{t}|c'(\tau)|d\tauL(t0,t)=∫t0t∣c′(τ)∣dτ.
事实上,这个积分对于C1C^1C1及更好的曲线都是well-define的,但是对连续的曲线不一定对,因为存在在闭区间上长度无穷的连续曲线。
比如 y=sin1xy=sin\frac{1}{x}y=sinx1 在 x=0x=0x=0 附近非道路连通;
比如把[0,1][0,1][0,1]映射到[0,1]×[0,1[0,1]\times[0,1[0,1]×[0,1的皮亚诺曲线;
再比如可以构造出[0,1]→[0,1]×[0,1][0,1]\rightarrow [0,1]\times[0,1][0,1]→[0,1]×[0,1]的双射:https://www.zhihu.com/question/301263376
把有限弧长的连续曲线,称为可求长曲线(rectifiable curve)
Definition 2.2.4 (Parametrization by arc length) 曲线c(s)c(s)c(s)称为弧长参数化的,如果L(s0,s)=s−s0L(s_0,s)=s-s_0L(s0,s)=s−s0,等价的说,∣c′(s)∣=1,∀s∈I|c'(s)|=1, \forall s\in I∣c′(s)∣=1,∀s∈I.
试想,现在我们通过改变曲线的参数,使得之前需要通过复杂的含绝对值积分才能得到的弧长,现在只需要做差就能够得到。通过这种方式,让计算变得如此简单,如果对后面计算曲率的过程有了解,就能感觉到这是一件非常令人振奋的事情,而且要是能把任意曲线都弧长参数化,那这事已经好的不能再好了。
直觉上来说这事肯定是对的,因为弧长参数化的直观理解就是以单位速度沿着曲线运动,这只要曲线本身的光滑程度够就行。
那既然是以单位速度沿着曲线运动,我只要以弧长为参数不就好了吗?这样用原来给定参数曲线的弧长来做参数,经过一个单位时间,等于经过一个单位弧长,等于沿弧长以单位速度运动。太妙了。
接下来是本篇笔记铺垫到现在,为了阐述的主要命题:
Proposition 2.2.5 每一个正则曲线(Regular Curve)都等价于一个弧长参数化曲线。
Proof .
存在性:假设c:I→Rn,t↦c(t)c: I \rightarrow \mathbb{R}^n,t\mapsto c(t)c:I→Rn,t↦c(t)为一光滑曲线,我们想找一个c~:I~→Rn,s↦c~(s)\tilde{c}: \tilde{I} \rightarrow \mathbb{R}^n, s\mapsto \tilde{c}(s)c~:I~→Rn,s↦c~(s),满足c(t)c(t)c(t)和c~(s)\tilde{c}(s)c~(s)的像集是一模一样的。
为此我们需要找一个映射ℓ:I→I~\ell: I\rightarrow\tilde{I}ℓ:I→I~. 之前我们考虑过,以ccc的弧长,来作为c~\tilde{c}c~的参数,也就是s=ℓ(t)=∫t0t∣c′(τ)∣dτs=\ell(t)=\int^{t}_{t_0}|c'(\tau)|d\taus=ℓ(t)=∫t0t∣c′(τ)∣dτ,写成交换图的形式:
因为ℓ′(t)=∣c′(t)∣>0\ell'(t)=|c'(t)|>0ℓ′(t)=∣c′(t)∣>0 (c is regular), 由反函数定理知道反函数存在且可导。这样,我们把交换图上边以ttt为参数的映射,转化成了下边以sss为参数的映射,接下来只需要验证∣c~′(s)∣=1|\tilde{c}'(s)|=1∣c~′(s)∣=1,由隐函数求导法则以及变上限函数的导数:
c′(t)=c~′(ℓ(t))=c~′(s)ℓ′(t)=c~′(s)∣c′(t)∣c'(t)=\tilde{c}'(\ell(t))=\tilde{c}'(s)\ell'(t)=\tilde{c}'(s)|c'(t)|c′(t)=c~′(ℓ(t))=c~′(s)ℓ′(t)=c~′(s)∣c′(t)∣
两边同时取绝对值,就得到了结论,存在性得证,并且知道c~(s)=c∘ℓ−1(s)\tilde{c}(s)=c \circ \ell^{-1} (s)c~(s)=c∘ℓ−1(s).
唯一性是显然的,假设有两条参数曲线都满足条件导数为1,因为是正则曲线,所以其导数可以去掉绝对值,同时积分。这样可以看出,参数起始于同一点的弧长参数化曲线是唯一的,证毕。
2.3 一般参数曲线化成弧长参数曲线的例子
怎么把一般参数曲线化成弧长参数曲线,其实在上面命题的证明中已经构造出来了。
数学证明很多时候是构造性的,怎么证明,就怎么使用;而有的证明只是用严格话的语言论述,说明条件满足某些定义,这种证明算是在教我们该怎么“说话”。二者都是需要的,但是对于推动数学进展来说,更重要的是前者,因为启发性更强,毕竟,好的证明使我们更加聪明。
这里求ℓ−1\ell^{-1}ℓ−1可能会是很困难的一步,因为ℓ\ellℓ本来是个积分函数,且变量在积分上限,这就意味着很多时候,我们只能在形式上表示出来,并不能真正求出ℓ−1\ell^{-1}ℓ−1(which mathematicians always say.)
举两个能求出来弧长参数化曲线的例子:
Example 2.3.1 c:R→R2c: \mathbb{R}\rightarrow \mathbb{R}^2c:R→R2, given by c(t)=(coset,sinet)c(t)=(\cos e^t, \sin e^t)c(t)=(coset,sinet).
这个形式一看就跟圆脱不了关系,只是问题在于他不是以单位速度沿圆周运动的,那么将他弧长参数化之后,应该会变成标准圆的参数方程才对。
Solution:
先确定自变量的范围,这里t∈(−∞,∞)t\in(- \infty , \infty)t∈(−∞,∞),
s=ℓ(t)=∫−∞t∣c′(τ)∣dτ=∫−∞teτdτ=et,s∈(0,∞)s=\ell(t)=\int^t_{-\infty}|c'(\tau)|d\tau=\int^t_{-\infty}e^{\tau}d\tau=e^t, s\in(0,\infty)s=ℓ(t)=∫−∞t∣c′(τ)∣dτ=∫−∞teτdτ=et,s∈(0,∞)
所以t=ℓ−1(s)=lns,s∈(0,∞)t=\ell^{-1}(s)=\ln s,s\in(0,\infty)t=ℓ−1(s)=lns,s∈(0,∞).
最终的弧长参数化曲线为:
c~(s)=c∘ℓ−1(s)=(coss,sins).\tilde{c}(s)=c \circ \ell^{-1}(s)=(\cos s, \sin s).c~(s)=c∘ℓ−1(s)=(coss,sins).
这也印证了我们的想法。
Example 2.3.2 c:[π2,3π2]→R2c:[\frac{\pi}{2},\frac{3\pi}{2}]\rightarrow \mathbb{R}^2c:[2π,23π]→R2, given by c(t)=(t−sint,1−cost)c(t)=(t-\sin t,1-\cos t)c(t)=(t−sint,1−cost).
这里的曲线称为摆线(cycloid),且θ\thetaθ取的是其中间的一段:
Solution:
s=ℓ(t)=∫π2t∣c′(τ)∣dτ=∫π2t(1−cosτ)2−sin2τdτ=∫π2t2sinτ2dτs=\ell(t)=\int^t_{\frac{\pi}{2}}|c'(\tau)|d\tau=\int^t_{\frac{\pi}{2}}\sqrt{(1-\cos \tau)^2-\sin^2 \tau}d\tau=\int^t_{\frac{\pi}{2}}2\sin \frac{\tau}{2}d\taus=ℓ(t)=∫2πt∣c′(τ)∣dτ=∫2πt(1−cosτ)2−sin2τdτ=∫2πt2sin2τdτ
⇒s=−4cost2+22,s∈[0,42],\Rightarrow s=-4\cos \frac{t}{2} +2\sqrt{2}, s\in[0,4\sqrt{2}],⇒s=−4cos2t+22,s∈[0,42],
因为在t2∈[π4,3π4]\frac{t}{2}\in[\frac{\pi}{4},\frac{3\pi}{4}]2t∈[4π,43π]上cos\coscos函数导数恒不为0,反函数存在,所以最终:
c~(s)=c∘ℓ−1(s)=(t−sint,1−cost)\tilde{c}(s)=c \circ \ell^{-1}(s)=(t-\sin t,1-\cos t)c~(s)=c∘ℓ−1(s)=(t−sint,1−cost)
⇒c~(s)=(2arccos22−s4−22−s8(22+s)2−2s2,2−2(22−s4)2).\Rightarrow \tilde{c}(s)=(2\arccos\frac{2\sqrt{2}-s}{4}-\frac{2\sqrt{2}-s}{8}\sqrt{(2\sqrt{2}+s)^2-2s^2}, 2-2(\frac{2\sqrt{2}-s}{4})^2).⇒c~(s)=(2arccos422−s−822−s(22+s)2−2s2,2−2(422−s)2).
我要知道结果会是这样,一开始我就不愿意求了。所以接下来的问题是有没有对一般非弧长参数适用的公式呢,这个问题将会在第四次笔记中讲到。
至此,曲线的参数化就讲完了,下一篇就初入微分几何的正题——Frenet标架。
参考:
[1]W. Klingenberg. A course in differential geometry. Translated from the German by David Hoffman. Graduate Texts in Mathematics, Vol. 51. Springer-Verlag, 1978.
[2]J. Shurman. Calculus and Analysis in Euclidean Space Undergraduate Texts in Mathematics,Springer-Verlag, 2016.
[3]厦门大学杨波老师的主页:http://math-faculty.xmu.edu.cn/display.aspx?tid=149
微分几何笔记(2) —— 曲线的参数化相关推荐
- 光滑曲线_微分几何笔记(2) —— 曲线的参数化
第二周讲完了Klingenberg的第一章Curves,做一点微小的笔记. 分成三个部分,本篇讲曲线的弧长参数:下一篇讲一般的Frenet标架及方程组:再下一篇讲二维三维空间曲线的curvature ...
- 绝对不能错过!计算机视觉Polygon Mesh Processing读书笔记——4微分几何中的曲线
流形 3D模型必须为流形.通俗地说,如果一个网格模型中存在多个(3个或以上)面共一条边,那么它就是非流形的(non-manifold),因为这个局部区域由于自相交而无法摊开展平为一个平面了.请看如图所 ...
- 光滑曲线_微分几何笔记(4) —— 二维三维空间中曲线的曲率以及环绕数
本篇文章我们从一般化的 空间回到我们生活的 空间,看看低维空间中的曲线有哪些性质,主要计算下在非弧长参数下的曲线,曲率挠率的一般表达式. 最后引入环绕数的概念,讲讲怎么数曲线转了多少圈. 4.1 ...
- 微分几何笔记(4) —— 二维三维空间中曲线的曲率以及环绕数
本篇文章我们从一般化的 Rn\mathbb{R}^nRn 空间回到我们生活的 R2,R3\mathbb{R}^2,\mathbb{R}^3R2,R3空间,看看低维空间中的曲线有哪些性质,主要计算下在非 ...
- soapUI学习笔记--用例字段参数化
字段参数化的简单操作 1.把Request新增一个TestCase 增加TestCase,下方会出现: 2.案例中,请求参数只有一个.先运行下请求,可以运行成功(保证接口是通的) 3.添加参数.见图中 ...
- 红橙Darren视频笔记 贝塞尔曲线实现消息拖拽粘性效果 画笔练习
效果图 要实现这样的效果 我们要先知道那段弧线是如何绘制的,实际上那段弧线就是贝塞尔曲线 一. 什么是贝塞尔曲线 https://www.jianshu.com/p/8f82db9556d2 上面有个 ...
- canvas学习笔记-贝塞尔曲线
3.4 贝塞尔曲线 canvas提供了两个绘制贝塞尔曲线api: ctx.quadraticCurveTo(cpx, cpy, x, y); 二次贝塞尔曲线,(cpx, cpy)控制点 (x, y)终 ...
- Unity笔记-贝塞尔曲线
定义 贝塞尔曲线(Bezier curve),又称贝兹曲线或贝济埃曲线,是应用于二维图形应用程序的数学曲线.一般的矢量图形软件通过它来精确画出曲线,贝兹曲线由线段与节点组成,节点是可拖动的支点,线段像 ...
- 学习笔记:曲线插值法
一.算法简介 1.曲线插值的方法是按照车辆在某些特定条件(安全.快速.高效)下,进行路径的曲线拟合,常见的有多项式曲线.双圆弧段曲线.正弦函数曲线.贝塞尔曲线.B样条曲线等. 二.算法思想 1.曲线插 ...
- 微分几何笔记(1):概念与定义
文章目录 概念与定义 微分流形 张量场与微分形式 联络与曲率 概念与定义 微分流形 拓扑流形: 设MMM是一个Hausdorff拓扑空间,若对每一点p∈Mp\in Mp∈M,都有ppp的一个邻域UUU ...
最新文章
- MySQL 高级 - 存储过程 - 语法 - while循环
- Study 1 —— HTML5概述
- python程序流程控制_python流程控制
- Actor并发模型入门
- php和js序列化,PHP中serialize和json序列化与反序列化的区别
- oracle脚本如何写,怎样写sql脚本
- 2020年接近尾声,我选择来鲲鹏开发者技术峰会学点干货!
- qt 串口助手_新手如何从零开发ROS小车 (ros串口通讯、PID运动控制、双轮差速模型解算...
- 小米盒子升级android tv后不能安装第三方软件下载,小米电视,小米盒子无法安装第三方应用,怎么办?...
- 30.STM32 DS18B20
- 七种常见的电子邮件安全协议简析
- 使用ffmpeg将视频切片并加密
- Android-PickerView系列之介绍与使用篇(一)
- TTL、RS232、485到底能传输多远距离?
- 嵌入式开发正在日薄西山
- 论文汇网站第三期改版完成
- Linux系统aboutyou,Linux字符设备驱动高级
- 红米3s进不了recovery_红米3s卡刷教程_红米3s用recovery刷第三方系统包
- win下搭建小程序服务器,win下搭建小程序服务器
- 2021CSP-J3