预备知识

定理a(t)a(t)是一个非零连续可微函数,则

  1. 函数a(t)a(t)的长度是常数,当且仅当a(t)a(t)0a'(t)\cdotp a(t)\equiv 0;
  2. 函数a(t)a(t)的方向不变,当且仅当a(t)×a(t)0a'(t)\times a(t)\equiv 0;
  3. 函数a(t)a(t)与一个确定方向垂直,当且仅当$(a(t),a'(t),a''(t))\text{(混合积)}\equiv 0;

proof.proof.
1.
ddta(t)2=2a(t)a(t)a(t)2=Consta(t)a(t)=0\begin{aligned} &\because \frac{d}{dt}|a(t)|^2=2a'(t)\cdotp a(t)\\ &\therefore |a(t)|^2=Const\iff a'(t)\cdotp a(t)=0 \end{aligned}
2.
\exists单位常向量b,a(t)=f(t)bb,\quad a(t)=f(t)\cdotp b
其中f(t)f(t)是连续可微函数
a(t)=f(t)ba(t)×a(t)=0\begin{aligned} \therefore\quad &a'(t)=f'(t)\cdotp b\\ \therefore\quad&a'(t)\times a(t)=0 \end{aligned}
反之,设a(t)×a(t)=0a'(t)\times a(t)=0,命b(t)=a(t)/a(t)b(t)=a(t)/|a(t)|,则b(t)=1|b(t)|=1由(1)知,b(t)b(t)=0b'(t)\cdotp b(t)=0,即b(t)a(t)=0b'(t)\cdotp a(t)=0.由定义可知
a(t)=f(t)b(t)a(t)=f(t)\cdotp b(t)
其中f(t)=a(t)f(t)=|a(t)|,则
a(t)=f(t)b(t)+f(t)b(t)a'(t)=f'(t)\cdotp b(t)+f(t)\cdotp b'(t)
a(t)×a(t)=f(t)b(t)×a(t)=0a'(t)\times a(t)=f(t)\cdotp b'(t)\times a(t)=0
因此b(t)a(t)b'(t)\parallel a(t).由于b(t)a(t)=0b'(t)\cdotp a(t)=0,故b(t)=0b'(t)=0,即b(t)是常向量,即向量函数a(t)a(t)的方向不变。

曲线论

  1. 参数方程

    将曲线看做动点的运动轨迹。三维空间上的曲线是连续映射a:IR3a:\mathbf{I\to R^3},对于tIt\in \mathbf{I},我们用
    a(t)=(a1(t),a2(t),a3(t))a(t)=(a^1(t),a^2(t),a^3(t))
    作为aa的参数化法(parametrization)。

    直线的参数方程:p\mathbf p,q\mathbf q为直线上的两点,
    a(t)=p+t(qp)a(t)=\mathbf p + t(\mathbf q -\mathbf p)
    方向向量qp\mathbf q -\mathbf p是该时刻的速度向量,记为v\mathbf v

  2. 正则曲线

    • 如果曲线的参数表示式中的函数一阶连续可微函数,这称其为光滑曲线。
    • 对于光滑曲线r=r(t)r=r(t), 假设对于曲线r=r(t)r=r(t)t=t0t=t_0r(t0)̸=0r'(t_0)\not =0,则称这一点为曲线上的正常点。
    • 曲线上所有点都是正常点时,则称曲线为正则曲线。
  3. 自然参数

    • 定义:对于正则曲线Γ:r=r(t)\Gamma: r=r(t),称积分s(t)=t0tr(η)dηs(t)=\int_{t_0}^t |r'(\eta)|d\eta 为曲线Γ\Gamma从点t0t_0tt的弧长。
    • 自然参数:以弧长作为曲线的参数,记为ss。曲线表示为r=r(s)r=r(s)
    • 定理r=r(t),atbr=r(t),a\le t\le bE3E^3中的一条正则曲线,则t是它的弧长参数的充要条件是r(t)=1|r'(t)|=1
      proof.proof. ds/dt=r(t),ds=dtr(t)=1.\because ds/dt=|r'(t)|,\therefore ds=dt \iff |r'(t)|=1.
      上述命题的直观意义是:曲线以弧长为参数的充要条件是它的切向量场是单位切向量场。为了强调已经把弧长取为曲线的参数,我们通常用“·”表示对于弧长参数的微商,如r˙(s)\dot r(s)等等。
    • 自然参数的性质:
      • 切向量具有单位长度r˙(s)=dr(s)ds=1|\dot{r}(s)|=|\frac{dr(s)}{ds}|=1,即单位速率曲线。
      • r(s)r˙(s),r˙(s)r¨(s)r(s)\perp \dot{r}(s),\dot{r}(s)\perp \ddot{r}(s)
  4. 曲率、挠率

    • 单位切向量(unit tangent vector):T(s)=r˙(s),T=1T(s)=\dot{r}(s),|T|=1
      p7
    • 定理Δθ\Delta \theta表示向量α(s+Δs)\alpha(s+\Delta s)α(s)\alpha(s)之间的夹角,则
      T˙(s)=limΔs0ΔθΔs|\dot T(s)|=\lim_{\Delta s\to 0}|\frac{\Delta \theta}{\Delta s}|
      proof.proof. 我们把曲线上的单位切向量T(s)T(s)平行移动到原点OO,则它的端点便描出一条单位球面上的曲线,Δθ\Delta \theta就是T(s+Δs)|T(s+\Delta s)T(s)T(s)所张的角,而T(s+Δs)T(s)|T(s+\Delta s)-T(s)正好是它所对的弦长。所以:
      T˙(s)=limΔs0T(s+Δs)T(s)Δs=limΔs02sinΔθ2Δs=limΔs0ΔθΔs\begin{aligned} |\dot T(s)|=&\lim_{\Delta s \to 0}\frac{T(s+\Delta s)-T(s)|}{\Delta s|}\\ =&\lim_{\Delta s\to 0}\frac{2|\sin \frac{\Delta\theta}{2}|}{|\Delta s|}\\ =&\lim_{\Delta s\to 0}\frac{|\Delta \theta|}{|\Delta s|} \end{aligned}
    • 曲率(curvature): κ(s)=T˙(s)=r¨(s)κ(s)=limΔs0ΔφΔs\kappa (s)=|\dot T(s)|=|\ddot{r}(s)| \quad \kappa(s)=\lim\limits_{\Delta s\to 0}|\frac{\Delta\varphi}{\Delta s}|,称T˙(s)\dot T(s)为该曲线的曲率向量。
    • 主法向量(principal normal vector):N(s)=r¨(s)k(s),N=1N(s)=\frac{\ddot{r}(s)}{k(s)},|N|=1
    • 副法向量(binormal):B(s)=T(s)×N(s),B=1B(s)=T(s)\times N(s),|B|=1
      :单位切向量,主法向量,副法向量也作α(s),β(s),γ(s)\alpha(s),\beta(s),\gamma(s)
    • 挠率:τ(s)=B˙(s)N(s)\tau(s)=-\dot{B}(s)\cdotp N(s)
    • 定理β(s)β(s)是具有单位速度的曲线,则
      1. κ=0\kappa =0 当且仅当β(s)β(s)是直线。
      2. κ>0,τ=0\kappa >0,\tau =0当且仅当β(s)β(s)是平面曲线。
  5. Frenet标架

    • Frenet标架:{r(s):T(s),N(s),B(s)}\{r(s):T(s),N(s),B(s)\}
    • 曲率刻画了曲线的弯曲程度,刻画了曲线偏离切线的程度。
    • 挠率刻画了曲线偏离密切面的程度,是曲线非平面化的量
      • 若曲线在T正向的一方,且在密切面上面,即也在B正向的这一面,这时挠率τ\tau取正。
      • 若曲线在正向的一方,且在密切面的下面,挠率τ\tau取负。
        切面
    • 定义:过空间曲线上P点的切线和P点邻近一点Q可作一平面σ,当Q点沿曲线趋于P时,平面σ的极限位置π称为曲线在P点的密切平面
    • 定义:过曲线上一点P的主法线的正侧取线段PC,使PC的长为1/k。以C为圆心以1/k为半径在密切平面上确定一个圆,这个圆称为曲线在P点的密切圆或曲率圆,圆的中心叫曲率中心,圆的半径叫曲率半径。
      Δs=RΔακ=limΔs0ΔαΔs=1R\begin{aligned} &\because \Delta s=R\Delta\alpha \\ &\therefore \kappa=\lim\limits_{\Delta s\to 0}|\frac{\Delta \alpha}{\Delta s}|=\frac{1}{R} \end{aligned}
      在点P处的曲率圆与曲线有下列密切关系:
      1. 有公切线
      2. 凹向一致
      3. 曲率相同
        密切平面
    • Frenet公式,其系数阵为反对称方阵
      (T˙(s)N˙(s)B˙(s))=(0k(s)0k(s)0τ(s)0τ(s)0)(T(s)N(s)B(s))\begin{pmatrix} \dot{T}(s) \\ \dot{N}(s) \\ \dot{B}(s) \end{pmatrix} = \begin{pmatrix} 0 & k(s) & 0 \\ -k(s) & 0 & \tau (s)\\ 0 & -\tau(s) & 0 \end{pmatrix} \begin{pmatrix} T(s) \\ N(s) \\ B(s) \end{pmatrix}
    • 达布向量(Darboux vector): ω=τT+κB\omega =\tau T+\kappa B
  6. 一般参数曲线的曲率、挠率的计算

    • 空间曲线:r=r(t)r=r(t)
      • 曲率:κ=r(t)×r(t)r(t)3\kappa=\frac{|r'(t)\times r''(t)|}{|r'(t)|^3}
      • 挠率:τ=r(t)×r(t)r(t)r(t)×r(t)2(混合积)\tau =\frac{r'(t)\times r''(t)\cdotp r'''(t)}{|r'(t)\times r''(t)|^2}\text{(混合积)}
        推导:
        已经知道:T(s)=dr(s)ds,κ=T˙(s)=r¨(s)T(s)=\frac{dr(s)}{ds},\quad \kappa=|\dot T(s)|=|\ddot r(s)|
        如果κ0\kappa \ne 0,则
        N(s)=r¨(s)r¨(s),B(s)=T(s)×N(s)=r˙(s)×r¨(s)r¨(s)N(s)=\frac{\ddot r(s)}{|\ddot r(s)|},\quad B(s)=T(s)\times N(s)=\frac{\dot r(s)\times \ddot r(s)}{|\ddot r(s)|}
        如果曲线的方程是r=r(t),tr=r(t),t不是弧长参数,则
        dsdt=r(t)\frac{ds}{dt}=|r'(t)|

        T(t)=r(t)r(t),r(t)=r(t)T(t).T(t)=\frac{r'(t)}{|r'(t)|},\quad r'(t)=|r'(t)|\cdotp T(t).
        对后一式再微分得到
        r(t)=dr(t)dtT(t)+r(t)dTdsdsdt=dr(t)dtT(t)+r(t)2κN(t),\begin{aligned} r''(t)=&\frac{d|r'(t)|}{dt}T(t)+|r'(t)|\cdotp \frac{dT}{ds}\cdotp \frac{ds}{dt}\\ =&\frac{d|r'(t)|}{dt}T(t)+|r'(t)|^2\cdotp \kappa N(t), \end{aligned}
        所以
        r(t)×r(t)=r(t)3κ(T(t)×N(t))=r(t)3κB(t),\begin{aligned} r'(t)\times r''(t)=&|r'(t)|^3\cdotp \kappa(T(t)\times N(t))\\ =&|r'(t)|^3\cdotp \kappa B(t), \end{aligned}
        由此得到
        κ=r(t)×r(t)r(t)3,B(t)=r(t)×r(t)r(t)×r(t)\kappa=\frac{|r'(t)\times r''(t)|}{|r'(t)|^3},\quad B(t)=\frac{r'(t)\times r''(t)}{|r'(t)\times r''(t)|}
        N(t)=B(t)×T(t)=r(t)2r(t)(r(t)r(t))r(t)r(t)r(t)×r(t)\begin{aligned} N(t)=&B(t)\times T(t)\\ =&\frac{|r'(t)|^2r''(t)-(r'(t)\cdotp r''(t))r'(t)}{|r'(t)|\cdotp |r'(t)\times r''(t)|} \end{aligned}
    • 平面曲线 y=f(x)y=f(x)
      • y=f(x)y=f(x)曲率:κ=y(1+y2)3/2\kappa=\frac{|y''|}{(1+y'^2)^{3/2}}
      • x=x(t),y=y(t)x=x(t),y=y(t)曲率:κ=x˙y¨x¨y˙(x˙2+y˙2)3/2\kappa=\frac{|\dot x\ddot y-\ddot x\dot y|}{(\dot x^2+\dot y^2)^{3/2}}
      • 挠率:0
  7. 曲线论基本定理

    • 曲线的弧长、曲率、挠率这三个量是由曲线本身的形状决定的,与它的参数表示和坐标系的选取无关。
    • 定理r=r1(s),r=r2(s)r=r_1(s),r=r_2(s)E3E^3中两条以弧长s为参数的正则曲线,如果他们的曲率出处不为0,并且κ1(s)=κ2(s),τ1(s)=τ2(s)\kappa_1(s)=\kappa_2(s),\tau_1(s)=\tau_2(s),则必存在E3E^3的一个刚体运动σ\sigma,它把曲线r=r2(s)r=r_2(s)变为曲线r=r1(s)r=r_1(s)

曲面论

  1. 参数曲面片
    所谓参数曲面片是指从E2E^2的一个区域DD到空间E3E^3中的一个连续映射。若在E2,E3E^2,E^3中分别建立笛卡尔坐标系,用(u,v)(u,v)E2E^2中点的坐标,用(x,y,z)(x,y,z)E3E^3中的点的坐标,则一个参数曲面排列的方程可以表示为:
    {x=x(u,v)y=y(u,v)z=z(u,v)(u,v)D,\left\{ \begin{aligned} x=&x(u,v)\\ y=&y(u,v)\\ z=&z(u,v) \end{aligned}\right.\qquad (u,v)\in D,
    或者写成向量方程:
    r=r(u,v)=(x(u,v),y(u,v),z(u,v))r=r(u,v)=(x(u,v),y(u,v),z(u,v))
    我们首先要求函数x(u,v),y(u,v),z(u,v)x(u,v),y(u,v),z(u,v)都是三次以上连续可微的。变量u,vu,v称为该曲面片的参数。
    p11
    在曲面S上取一定点P0=r(u0,u0)P_0=r(u_0,u_0),让uu变化,而u=u0u=u_0,则动点描出一条落在曲面S上的曲线,这条曲线称为过点P0P_0的v-曲线,它的方程是
    u=u0u=u_0
    或者
    r=r(u0,v)r=r(u_0,v)
    同理,也有u-曲线,这样,经过每一点有一条u-曲线和一条v-曲线,它们构成曲面上的参数曲线网。通常把(u,v)(u,v)称为曲面上的曲纹坐标
  2. 正则曲面片
    如果ru(u0,v0),rv(u0,v0)r_u(u_0,v_0),r_v(u_0,v_0)是线性无关的,则称曲面片S在点P0P_0是正则的。三次以上连续可微的、处处是正则点的曲面片,即所谓的正则曲面片
  3. 曲面的Monge形式
    z=z(u,v)=z(f(x,y),g(x,y))z=z(u,v)=z(f(x,y),g(x,y))
    用Monge形式给出的曲面片总是正则的
  4. 曲面正侧
    规定,向量ru×rvr_u\times r_v所指的一侧为曲面正侧。参数u,v的次序决定了参数曲面的定向。
  5. 切平面、法线
    切平面的单位法向量
    n=ru×rvru×rvn=\frac{r_u\times r_v}{|r_u\times r_v|}
    曲面r=r(u,v)r=r(u,v)在点(u,v)(u,v)的切平面的参数方程是
    X(λ,μ)=r(u,v)+λru+μrvX(\lambda,\mu)=r(u,v)+\lambda r_u+\mu r_v
    其中λ,μ\lambda,\mu是切平面上动点的参数,法线的参数方程是
    X(t)=r(u,v)+tn(u,v)X(t)=r(u,v)+tn(u,v)
    曲面上的每一点,由其参数方程定义了一个标架{r;ru,rv,n}\{r;r_u,r_v,n\},称为曲面上的自然标架。

曲面的第一基本形式

设有正则参数曲面片S:r=r(u,v)S:r=r(u,v),曲面S在每一点的切空间是有切向量ru(u,v),rv(u,v)r_u(u,v),r_v(u,v)张成的二维向量空间,它是R3R^3的子空间,因此S的切向量作为R3R^3中的向量,可以求它们的长度及夹角。前面已经说过,曲面S在任意一点r(u,v)r(u,v)的切向量是
dr(u,v)=ru(u,v)du+rv(u,v)dvdr(u,v)=r_u(u,v)du+r_v(u,v)dv

(1)

其中(du,dv)(du,dv)是切向量dr(u,v)dr(u,v)在自然基底ru,rv{r_u,r_v}下的分量。
但是一般说来,ru,rv{r_u,r_v}不是单位正交基底,如果我们知道基底的度量系数,则表示成(1)式的切向量的内积就能够用他们在基底ru,rv{r_u,r_v}下的分量du,dvdu,dv来表达了。命
E(u,v)=ruruF(u,v)=rurvG(u,v)=rvrv\begin{aligned} E(u,v)=&r_u\cdotp r_u \\ F(u,v)=&r_u\cdotp r_v \\ G(u,v)=&r_v\cdotp r_v \end{aligned}
他们就是基底ru,rv{r_u,r_v}的度量系数,称为曲面S的第一类基本量
通常记作:
(EFFG)\begin{pmatrix} E & F\\ F & G \end{pmatrix}
显然这是一个正定矩阵
I=Edu2+2Fdudv+Gdv2=(du,dv)(EFFG)(dudv)\begin{aligned} I=&Edu^2+2Fdudv+Gdv^2\\ =&(du,dv) \begin{pmatrix} E & F\\ F & G \end{pmatrix} \begin{pmatrix} du \\ dv \end{pmatrix} \end{aligned}
二次微分式II在容许参数变换下是不变的,我们称II为曲面的第一基本形式
II的几何意义
I=drdrI=dr\cdotp dr
此外若命
δr=ruδu+rvδv\delta r=r_u\delta u+r_v\delta v
这切向量的内积是
drδr=(du,dv)(EFFG)(δuδv)dr\cdotp \delta r=(du,dv) \begin{pmatrix} E & F\\ F & G \end{pmatrix} \begin{pmatrix} \delta u \\ \delta v \end{pmatrix}

(2)

因此
dr=Edu2+2Fdudv+Gdv2cos(dr,δr)=drδrdrδr\begin{aligned} |dr|=&\sqrt{Edu^2+2Fdudv+Gdv^2}\\ \cos \angle(dr,\delta r)=&\frac{dr\cdotp \delta r}{|dr||\delta r|} \end{aligned}

(3)

特别的,drdrδr\delta r正交(1)=0\iff(1)=0

若在曲面上有一条曲线,它的方程是
u=u(t),v=v(t),atb.u=u(t),\quad v=v(t),\quad a\le t\le b.
由(3)知,它的弧长是
L=abr(t)dt=ab[E(u(t),v(t))(dudt)2+2F(u(t),v(t))dudtdvdt+G(u(t),v(t))(dvdt)2]1/2dt\begin{aligned} L=&\int_a^b|r'(t)|dt\\ =&\int_a^b[E(u(t),v(t))(\frac{du}{dt})^2+2F(u(t),v(t))\frac{du}{dt}\frac{dv}{dt}+G(u(t),v(t))(\frac{dv}{dt})^2]^{1/2}dt \end{aligned}

最后我们来讨论曲面上一个区域的面积计算。
p18
设S的方程是r=r(u,v),(u,v)DE2r=r(u,v),(u,v)\in D\subset E^2,考虑曲面上由参数曲线u=u0,u=u0+Δu,v=v0,v=v0+Δvu=u_0,u=u_0+\Delta u,v=v_0,v=v_0+\Delta v所围成的一小块,它的面积与切平面上由Δuru,Δvrv\Delta u\cdotp r_u,\Delta v\cdotp r_v所张成的平行是变形的面积在略去更高阶无穷小之后是相同的,而后者的面积是
(Δuru)×(Δvrv)=ru×rvΔuΔv=EGF2ΔuΔv\begin{aligned} |(\Delta u\cdotp r_u)\times(\Delta v\cdotp r_v)|=&|r_u\times r_v|\Delta u\Delta v\\ =&\sqrt{EG-F^2}\Delta u\Delta v \end{aligned}

dA=EGF2dudvdA=\sqrt{EG-F^2}dudv
称为曲面S的面积元素,那么S的面积是
A=DEGF2dudvA=\iint_D\sqrt{EG-F^2}dudv
根据重积分的替换法则以及第一类基本量的变换规律不难知道,A与曲面上参数的容许变换是无关的。

  1. 保长对应和保角对应
  2. 可展曲面

曲面的第二基本形式

参考资料

  1. 微分几何初步 陈维桓