当前位置:职场发展 > 表面局部理论介绍-从表面的概念和基本形式到高斯曲率及其Pthyon计算

表面局部理论介绍-从表面的概念和基本形式到高斯曲率及其Pthyon计算

  • 发布:2023-09-30 10:55


前言

本文尽可能详细地介绍了曲面的局部理论,最后使用Python编程求出曲面的高斯曲率 。

解释:第1章到第6章是局部表面理论的介绍。这是表面入门的必读之作。这将是非常数学化的,需要耐心去研究。第七章是Python在曲面曲率计算中的应用,你将熟悉常见曲面的曲率。这篇文章虽然内容很多,但并不是不可读。我花了差不多一天的时间才完成了从头到尾全部内容的学习和理解。具体学习时间取决于读者自己的数学知识。


1.曲面的概念

说明:本文 R 3 \mathbb{R}^3 R3 表示具有内积结构的欧氏空间。

定义1.1 从平面区域 D D D R 3 \mathbb{R}^3 R3 映射 r ( u , v ) = ( x ( u , v ) , y ( u , v ) , z ( u , v ) ) , \bold{r}(u,v) = (x(u,v), y(u, v), z(u,v)), r(u,v)=(x (u,v),y (u,v), z(u,v)) , 满足 (1) 每个分量函数无限可微,(2) 向量 r u \bold{r}_u ru r v \bold{r }_v rv 线性独立,据说 r \bold{r} r R 3 \mathbb{R}^3 R 3 A 表面 , ( u , v ) (u,v) (u,v )称为坐标参数表面。

满足定义1.1的曲面在有些书中也称为正则曲面

设置 r ( u , v ) \bold{r}(u,v) r(u ,v)如上所定义,对于任意点 ( u 0 , v 0 ) ∈ D (u_0,v_0)\in D (u0,v0)εD,易验证有下式成立 r u × r v ∣ ( u 0 , v 0 ) = ( ∂ ( y , z ) ∂ ( u , v ) , ∂ ( z , x ) ∂ ( u , v ) , ∂ ( x , y ) ∂ ( u , v ) ) ∣ ( u 0 , v 0 ) \bold{r}_u \times \bold{r}_v |_{(u_0,v_0)}=\left.\left( \frac{\partial (y,z)}{\partial (u,v)}, \frac{\partial (z,x)}{\部分 (u,v)}, \frac{\部分 (x,y)}{\部分 (u,v)} \right) \right|_{(u_0,v_0)} ru ×rv(u0,v0)=((u) ,v)(y) ,z),∂ (u,v)(z,x), ∂(u,v) (x,y) ​) u 0v0 表面有两种表达形式, 的形状显示表达 z = f ( x , y ) z=f(x,y) z = f (x,y) ,即 r ( x , y ) = ( x , y , f ( x , y ) ) \bold{ r}(x,y)=(x,y,f(x,y)) r(x,y )= ( x,y)) ,由上式易知 r u × r v = ( − f x , − f y , 1 ) ≠ 0 , \bold{r}_u \times \bold{r}_v = (-f_x \;, -f_y \;, \;1 ) \neq \bold{0} , ru×rv= (fx ,fy ,1) =0 r v \bold{r}_v rv 线性无关。稍后我们将在向量中使用叉积符号 × \times × 楔积 符号 ∧\楔来替换以与后面的几何图形对齐。

隐式表达式看起来像 F ( x , y , z ) = 0 F(x,y,z)=0 F(x, y,z)=0,当 F z ( x 0 , y 0 , z 0 ) ≠ 0 F_z(x_0,y_0,z_0 ) \neq 0 Fz(x0,y0,z 当​)=0时,由隐函数定理 ,在 ( x 0 , y 0 ) (x_0,y_0) ( x0,y0) 当地社区内 D D D, F ( x 0 , y 0 , z 0 ) F(x_0,y_0,z_0) F(x0, y0,z0) 有一个显式表达式 z = f ( x , y ) , x , y ∈ D , z 0 = f ( x 0 , y 0 ) , z= f(x,y), \; x,y),x ,yεD,z0=f(x0, y0), 类似地,当 F x ( x 0 , y 0 , z 0 ) ≠ 0 F_x(x_0,y_0,z_0) \neq 0 Fx(x0,y0,z0)=0 F y ( x 0 , y 0 , z 0 ) 0 ​,y0​ 当,z时也可以这样做​)=0确定表面。所以当 ∇ F ( x 0 , y 0 , z 0 ) ≠ 0 \nabla F(x_0,y_0,z_0) \neq 0 F (x0,y0,z0 )=0, F ( x , y , z ) F(x,y,z) F(x,y ,z) ( x 0 , y 0 , z 0 ) (x_0,y_0,z_0) (x0, y0, z 在 0) 附近定义一个表面。

考虑表面的不同参数表示,给定 r ( u , v ) : D → R 3 \bold{r}(u,v): \; D \to \mathbb{R}^3 r(u,v ) :D R3参数变换 σ : ( u ˉ , v ˉ ) ∈ D ˉ → ( u , v ) ε D , \sigma: (\bar{u} ,\栏{v}) \in \bar{D} \ 到 (u,v) \in D, σ:(uˉ,v ˉ)εDˉ→ (u,v)εD, 其中 σ \sigma σ 是 bi投影和该变换的 雅可比行列式满足 ∣ ∂ ( u , v ) ∂ ( u ˉ , v ˉ ) ∣ = ∣ ∂ u ( u ˉ , v ˉ ) ∂ u ˉ ∂ v ( u ˉ , v ˉ ) ∂ u ˉ ∂ u ( u ˉ , v ˉ ) ∂ v ˉ ∂ v ( u , v ) ∂ v ˉ ∣ ≠ 0 。 \left|\frac{\partial(u,v)}{\partial(\bar{u},\bar{v})}\right| = \left |\begin{array}{cccc} \frac{ \partial u(\bar{u},\bar{v})}{\partial\bar{u}} &\frac{\partial v(\条{u},\条{v})}{\部分\条{u}} \\ \frac{\部分u(\条{u},\条{v})}{\部分\条{v }} &\frac{\partial v(\bar{u} ,\bar{v})}{\partial\bar{v}}\\ \end{ar射线}\右| \neq 0.(uˉ,vˉ)(u, v)=uˉ u(uˉ,vˉ)vˉu(uˉ,v ˉ)u ˉv(uˉ,vˉ)vˉv(uˉ ,vˉ)  =0.这样,我们得出曲面的新参数表示 r ( u ˉ , v ˉ ) = r ∘ σ ( u ˉ , v ˉ ) = r ( u ( u ˉ , v ˉ ) , v ( u ˉ , v ˉ ) ) \bold{r}(\bar{u},\ bar{v}) = \bold{r} \circ \sigma(\bar{u},\bar{v})=\bold{r}(u(\bar{u}, \bar{v}), v(\bar{u}, \bar{v})) r(uˉ, vˉ)=rσ(uˉ ​​,vˉ)=r (u(uˉ, vˉ),v(uˉ,vˉ)) 我们记下表面为 S S S,以下称为 r = r ( u , v ) \bold{r}=\bold{r}(u,v) r=r( u , v) 是曲面 参数表示。以下是表面参数表示的几个示例。

1节1.1{x = r cos⁡(u)cos⁡(v)y = r cos⁡(u)sin⁡(v)z = rsin⁡(u)\ left \ { \begin{对齐} &x=R\cos(u)\cos(v) \\ &y=R\cos(u)\sin(v) \\ &z=R\sin(u)。 \ 结束 { 对齐 } \ 右 x=Rcos(u)cos (v)y=Rcosusin() v)z=Rsin (u).其中 D ˉ = { ( u , v ) : − π 2 < u < π 2 , 0 < v < 2 π } \bar{D} =\{(u,v): -\frac{\in}{2} Dˉ=2{(u, v):2π 2<u<2π ,0<v<2π}

示例 1.2 环面可以视为 x z xz xz 平面上的圆绕轴旋转 z z z 的曲面,以x轴上的点 R R R为圆心,取 r ( r < R ) r (r (r<R) 构造一个半径为 { x = R + r cos ⁡ ( u ) , z = r sin ⁡ ( u ) , \ left\{ \begin{aligned} &x=R+r\cos(u), \\ &z=r\sin(u), \end{aligned} \right。 {x =R+rcos(u),z=r(u ),然后把这个圆圈包起来 z z z 旋转轴得到圆环 { x ( u , v ) = ( R + r cos ⁡ ( u ) ) cos ⁡ ( v ) y ( u , v ) = ( R + r cos ⁡ ( u ) ) sin ⁡ ( v ) z ( u , v ) = R sin ⁡ ( u ) 。 \left\{ \begin{对齐} &x(u,v)=(R+r\cos(u))\cos(v) \\ &y(u,v)=( R+r\cos(u)) \sin(v)\\ &z(u,v)=R\sin(u)。 \end{对齐} \right。 x(u ,v) =(R) +rcos(u))cos(v)y(u, v)=(R+ rcos(u))sin (v)z(u ,v)=R (u). 一般来说,我们有旋转曲面的参数表示:

示例 1.3 x z xz xz 参数曲线 x = f ( u ) , z = g ( u ) x=f(u), \ ;\; z=g(u) x=f(u ),z=g( u ) 围绕 z z z 旋转获取旋转面的轴 r ( u , v ) = ( f ( u ) cos ⁡ (v),f(u)sin⁡(v),g(u))。 \bold{r}(u,v)=(f(u)\cos(v), f(u)\sin(v), g(u))。 r(u,v) =(f(u) cos (v),f( u)(v) ,g(u)).

示例1.4圆柱编号公式r(u,v)=(a cos⁡(u),sinsin⁡(u),b v)\ bold {r}(u,v)=(u,v)=( a\cos(u), a\sin(u), bv) r(u,v )= (acos(u),a 罪恶(u), bv) 其中 a > 0 , b a>0, b a >0,b城市不变。若定 − π < u < π -\piπ<u<π - ∞ < v < + ∞ -\infty<v<+,显示移除线圆柱面上 x = − 1 , y = 0 , z = b v x=-1, y=0, z=bv x= - , y=0,z=bv 收入实现部分。


2。切割平面给定方向

考虑弯曲 S S S ,参数表示为 r ( u , v ) = ( x ( u , v ), y ( u , v ), z ( u , v ) ) ) \bold{r}(u,v)=(x(u,v), y(u,v), z(u,v)) r(u . * D 。固定 u = a u=a u=a ,曲线 r ( a , v ) \bold{r}(a ,v) r(a,v) 在 v = b v=b v=切影兰是 r v ( a , b ) = d r d v ( a , b ) \bold{r}_v(a,b)=\frac{d\bold{r}}{dv}(a,b) rv(a, b)=r(a,b) v = b v=b b ,有 u = a u=a u=ar u ( a , b ) = d r d u ( a , b ) \bold{r}_u(a,b)=\frac{d\bold{r}}{du}(a,b) r u(a,b) =dudr(a,b)

表面确定知识, r u ( a , b ) \bold{r}_u(a,b) ru(a r u ( a , b ) \bold{r}_u(a,b) ru (a ,b)没有电线,就是 r u ( a , b ) \bold{r}_u(a,b) ru( a,b) 总和 r v ( a , b ) \bold{r}_v(a,b) rv (a, S S 相切优 P 0 = P ( a , b ) P_0 = P(a,b) P0=P(a,b) 点状平面,对称面 S S S 目前 P 0 P_0 P0切割平面,记录 T P 0 S T_{P_0}S TP ​S.

T T_{P_0}S
TP0S
垂直直线称为该点表面的 法线 r u ( a , b ) ∧ r v ( a , b ) \bold{r}_u(a,b)\wedge\bold{r}_v(a,b) ru(a,b)rv(a, b) 是曲面 S S S P 0 P_0 P0 法向量点 的 ,以及 { P 0 : r u , r v , r u ∧ r v } \{ P_0: \; \bold{r}_u \;,\bold{r}_v \;, \bold{r}_u \楔形\bold{r}_v\} {P0ru,rv,rurv} 成分 R 3 \mathbb {R}^3 R3 是一个 自然定向框架

下面给出曲面 S S S 的切平面与参数选择的关系。设 σ : ( u ˉ , v ˉ ) ∈ D ˉ → ( u , v ) ε D \sigma: (\bar{u}, \bar{v}) \in \bar{D} \to ( u,v) \in D σ:(uˉ ,vˉ)ε Dˉ(u,v )εD是参数变换,那么基变换矩阵为参数 变换后的 雅可比矩阵 ∂ ( u , v ) ∂ ( u ˉ , v ˉ ) \frac{\partial (u,v)}{\partial (\bar{ u},\bar {v})} (uˉ,vˉ)( u,v),即 ( r u ˉ r v ˉ ) = ( ∂ u ∂ u ˉ ∂ v ∂ u ˉ ∂ u ∂ v ˉ ∂ v ∂ v ˉ ) ( r u r v ) \left( \begin{aligned} \bold{r}_{\bar{u}} \\ \bold{r}_{\bar{v}} \end {aligned} \right) = \left( \begin{array}{cc} %矩阵共有3列,每列居中 \frac{ \partial u}{ \partial \bar{u}} & \frac { \partial v}{ \partial \bar{u}} \\ %第一行元素\frac{ \partial u}{ \partial \bar{v}} & \frac{ \partial v}{ \partial \bar{ v}} \\ %th 两行元素\end{array} \right) \left( \begin{Aligned} \ bold {r} _ {u} \\ \ bold {r} _ {v} \ end {aligned} \ right)ruˉrvˉ)= (uˉuvˉ∂ uuˉvvˉ∂ vr) urv) 特别地,在 P 0 P_0 P0 点, r u ˉ ∧ r v ˉ = ∂ ( u , v ) ∂ ( u ˉ , v ˉ ) r u ∧ r v ≠ 0 \bold{r}_{\bar{u}} \wedge \bold{r }_{\bar{v}}= \frac{\partial (u,v)}{\partial (\bar{u},\bar{v})} \bold{r}_{u} \wedge \粗体{r}_{v} \neq 0 ruˉr vˉ =(uˉ,v ˉ)(u,v)rurv=0 。于是我们就有了

性质2.1曲面的切平面和法线与曲面参数的稀疏相关。

很容易看到,当 ∂ ( u , v ) ∂ ( u ˉ , v ˉ ) > 0 \frac{\partial (u,v)}{\partial (\bar{u},\bar {v})}>0 (uˉ,v ˉ)() u,v)>0时, r u ˉ ∧ r v ˉ \bold{r}_{\bar{u}}\wedge \粗体{r}_{\bar{v}} ruˉr vˉ r u ∧ r v \bold{r}_u \wedge \bold{r}_v rurv 的方向相同;当 ∂ ( u , v ) ∂ ( u ˉ , v ˉ ) < 0 \frac{\partial (u,v)}{\partial (\bar{u},\bar {v})}<0 (uˉ,v ˉ)() u,v)<0时, r u ˉ ∧ r v ˉ \bold{r}_{\bar{u}}\wedge \粗体{r}_{\bar{v}} ruˉr vˉ r u ∧ r v \bold{r}_u \wedge \bold{r}_v rurv的可爱方向。

因此,当 ∂ ( u , v ) ∂ ( u ˉ , v ˉ ) > 0 \frac{\partial (u,v)}{\partial (\bar{u},\bar{v })}>0 (uˉ,v ˉ)(u)当,v)>0时,对应的参数变换称为 同方向参数变换;反之称为反向参数变换

我们再次检查微分 d r = r u ( u , v ) d u + r v ( u , v ) d v d\bold{r} = \bold{r}_u(u,v)du + \bold{ r}_v(u,v)dv dr=ru(u,v)d u+rv ​(u,v)dv 其中 d u , d v du, dv du ,d v 正是关于自然碱基 { r u , r v } \{ \bold{r}_u, \bold{r}_v \} { r u,rv} 的重量。

现在 d u , d v du, dv du, dv 被认为独立于 u , v u,v 对于另外两个自变量uv,微分公式 r u ( u , v ) d u + r v ( u , v ) d v \bold{r}_u( u ,v)du + \bold{r}_v(u,v)dv ru(u ,v )du+rv​ (u,v)d v 可以代表 T P S T_PS TPS 中的任何切向量。考虑以下定理 2.2:

定理2.2 V V V n n n 维向量空间,取一组基为 { α 1 , . 。 。 , α n } \{ \alpha_1, ..., \alpha_n \} {α1,... ,αn},其双空间 V ∗ V^* V*的底数为 f 1 , . 。 。 , f n f_1,...,f_n f1,..., fn,然后 α = Σ i = 1 n f i ( α ) α i , ∀ α ∈ V , f = Σ i = 1 n f ( α i ) f i , ∀ f ∈ V ∗ \begin{aligned} &\alpha = \sum_{ i=1}^{n} f_i(\alpha)\alpha_i, \quad \forall \alpha \in V, \\ &f = \sum_{i=1}^{n} f(\alpha_i)f_i , \quad \forall f \in V^* \end{对齐} α=i=1Σ nfi(α)α i,αV,f=i=1Σ nf(α i)fi, ∀fεV* 因为 T PS T_PS TPS 底数为 { r u , r v } \{\bold{r}_u, \bold{r}_v\} {ru,rv } ,我们假设 { f 1 , f 2 } \{ f_1, f_2 \} {f1,f2} 是对偶空间 T P ∗ S T_P ^*STP*S v ∈ T P S \bold{v}\in T_PS vεTPS ,根据定理 2.2 我们有 v = f 1 ( v ) r u + f 2 ( v ) r v \bold{ v} = f_1(\粗体{v})\粗体{r}_u + f_2(\粗体{v})\粗体{r}_v v=f1( v)ru+f2(v)r v 对比度微分 r u ( u , v ) d u + r v ( u , v ) d v \bold{r}_u(u,v)du + \bold{r} _v(u,v )dv ru(u, v)du +r v(u,v )dv 可以表示 v \bold{v } v d u , d v du, dv du,dv可以视为切空间 T P S T_PS TPS 上的线性函数,即 d u ( v ) = f 1 ( v ) , d v ( v ) = f 2 ( v ) 。 du(\bold{v})=f_1(\bold{v}), \\ dv(\bold{v})=f_2(\bold{v})。 du(v)= f1(v),dv(v)=f 2(v).因此, { d u , d v } \{du, dv\} {du,dv}成分 T P * S T_P^*S TP*S 空间基为 { r u , r v } \{ \bold {r}_u, \bold{r}_v \} {ru ​, rv} 双基础。在上面 v \bold{v} v 的展开中,尤其是 v = r u \bold{v}=\bold{r}_u v =ru,则我们得到 r u = d u ( r u ) r u + d v ( r u ) r v \bold{ r}_u = du(\bold{r} _u)\bold{r}_u+dv(\bold{r}_u)\bold{r}_v ru=du(ru) ru+d v(ru)rv以及取 v = r v \bold{v}= \bold{r}_v v=rv ,则 r v = d u ( r v ) r u + d v ( r v ) r v \bold {r}_v = du(\bold{r}_v)\bold{r}_u+dv(\bold{r}_v)\bold{r}_v rv=du(rv​ )ru+d v(rv) rv 从而结果 d u ( r u ) = d v ( r v ) = 1 du(\bold{r}_u)=dv(\bold{r}_v)=1 du(r u)=dv(r v)=1 d u ( r v ) = d v ( r u ) = 0 du(\bold{r}_v)=dv(\bold{r}_u)=0 du(r) v )=dv( ru)=0

事实上,任何向量 v \bold{v} v 都在切空间中 T P S T_PS TP S { r u , r v } \{ \bold{r}_u, \bold{r}_v \} {r 的基础u,rv} 下的坐标为 ( d u ( v ) , d v ( v ) ) (du(\bold{v}), dv (\bold {v})) (du(v) )dv() v)) ,其中 d u , d v du, dv du ,dv 是对偶空间 T P * S T_P^*S TP*S


3。曲面的第一种基本形式

假设 S S S21R3、 的表面 r = r ( u , v ) \bold{r}=\bold{r}(u, v) r=r(u ,v)为其参数表示,则任意切向量 S S S v \bold{v} v 都可以表示为以下形式 v = λ r u + μ r v \bold{v} = \lambda \bold{ r}_u + \mu \bold{r}_v v=λr u ​+μrv 那么 < v , v > = λ 2 < r u , r u > + 2 λ μ < r u , r v > + μ 2 < r v , r v > <\bold{v},\;\bold{v}>=\lambda^2 <\bold{r}_u, \; \bold{r}_u> + 2 \lambda\mu<\bold{r}_u, \; \bold{r}_v> + \mu^2<\bold{r}_v, \; \bold{r}_v> <v, v>= λ2<ru,ru> +2λμ<ru ,rv>+μ2< rvrv>

E = < r u , r u > E=<\bold{r}_u, \; \bold{r}_u> E=< ru,ru > F = < r u , r v > F = <\bold{r}_u, \; \bold{r}_v> F=<ru, rv> G = G=<\bold{r}_v, \; \bold{r}_v> G=<rv ,rv> ,考虑 S S S 上一曲线 r ( t ) = r ( u ( t ) , v ( t ) ) \bold{r}(t)=\bold{r}(u(t),v(t)) r(t )=r(u( t),v(t)) 其切向量为 d r ( t ) d t = r u u ′ ( t ) + r v v ′ ( t ) \frac{d\bold{r}(t)}{dt} = \bold{ r}_u u'(t) + \bold{r}_v v'(t) dtdr() t)=ruu(t)+rv v(t) 曲线在 a < t ′ < t aa <t<t间的弧长​​为 s = ∫ a t ∣ d r ( t ′ ) d t′ ∣ d t ′ = ∫ a t E ( d u d t ′ ) 2 + 2 F d u d t ′ d v d t ′ + G ( d v d t ′ ) 2 d t ′ \begin{aligned} s &= \int_{a}^{t} \left| \frac{d\bold{r}(t')}{dt'} \right| dt' \\ &=\int_{a}^{t} \sqrt{E\left( \frac{du}{dt'} \right)^2+2F\frac{du}{dt'}\frac{ dv}{dt'} + G\left( \frac{dv}{dt'} \right)^2} \;dt' \end{对齐} s=at dtd r(t)dt=atE(dt du)2+2Fd tdudtdv+G (dtd v)2 d t因此, I ≜ d s 2 = E d u ⋅ d u + 2 F d u ⋅ d v + G d v ⋅ d v I \triangleq ds^2 = Edu\cdot du+2Fdu \cdot dv +Gdv\cdot dv Ids2=Edu ⋅du+2Fdudv+Gdvdv称为曲面 S S S 的第一个基本形式。

定理3.1 曲面 S S S 的第一基本形式与参数选择无关。

结合 d r = r u d u + r v d v d\bold{r}=\bold{r}_u du+\bold{r}_v dv d r =2ru du+rv sdmv,最好的我= E d u d u + 2 F d u d v + G d v d v = < r u , r u > d u d v + 2 d u d v + < r v , r v > d v d v = \be gin {对齐} I &= Edudu+2Fdudv +Gdvdv \\ &= <\bold{r}_u, \; \bold{r}_u> dudv+ 2<\bold{r}_u, \; \bold{r}_v>dudv+<\bold{r}_v, \; \bold{r}_v>dvdv \\ &= \end{对齐} I =Edud u+Fdud v+Gdvdv=<ru,ru >dudv+2<ru,rv>du dv+<rv ,rv> ddr, dr> I I I 最昂贵的雪花之一

回想一下 合约变换 ,令 T T T 是一个正交矩阵,令 T ( P ) = P ⋅ T + P 0 \mathfrak{ T }(P) = P\cdot T+P_0 T(P)= PT+P0 ,则 T \mathfrak{T} T 是一个契约变换。常见的旋转、平移和反射变换都是契约变换。我们给出定理3.2,但没有证明。

定理3.2 曲面的第一基本形式在 R 3 \mathbb{R}^3 R的契约变换下保持不变。

这里有一些例子。

? (
u,v)= (u,v,c ) ( c c c系统常数), r u = ( 1 , 0 , 0 ) \bold {r}_u=(1,0,0) ru=(1,0, ) , r v = ( 0 , 1 , 0 ) \bold{r }_v=(0,1,0) rv=(0, 1,0),所以第一个基本形式为 I = d u d u + d v d v I=dudu+dvdv I=d udu+dvdv

? : (x(u) , y(u)) x y xy xy 平面上的一条正则参数曲线,设 C C C 沿 z z z 轴移动获得曲面的方向 r ( u , v ) = ( x ( u ) , y ( u ) , v ) \bold{r}(u,v) = (x(u), y(u) , v) r(u,v) =(x(u), y(u), v)称为圆柱体。由于 r u = ( x ′ , y ′ , 0 ) \bold{r}_u=(x',y',0) ru​ = (x,y, ) , r v = ( 0 , 0 , 1 ) \bold{r } _v=(0,0,1) rv=(0, 0,1),则圆柱A基本形式为 I = ( ( x ′ ) 2 + ( y ′ ) 2 ) d u d u + d v d v I = ((x')^2+(y')^2)dudu+dvdv I=((x)2+ () y)2)d udu+dvdv输入 u u u C C C ∥ ∥ r u ∥ = 1 \parallel\bold{r}_u \parallel=1 2r u ∥=1,故 I = dudu + dv dv I = dudu+dvdv I=dud u+dvdv 例3.5
(球❢)无穷小 R R R 无穷小 r (θ, phi)。 = ( R cos ⁡ ( θ ) cos ⁡ ( Φ ) , R cos ⁡ ( θ ) sin ⁡ ( Φ ) , R sin ⁡ ( θ ) ) \bold{r}(\theta, \phi) = (R\cos (\theta)\cos(\phi), R\cos(\theta)\sin(\phi), R\sin(\theta)) Rcos (θ)cos(phi),Rcosθsin(φ),Rsin(θ) ) r θ = ( − R sin ⁡ ( θ ) cos ⁡ ( Φ ), − R sin ⁡ ( θ ) ) sin ⁡ ( Φ ), R cos ⁡ ( . θ ) ) r Φ = ( − R cos ⁡ ( θ ) sin ⁡ ( Φ ), Rcos ⁡ ( θ ) cos ⁡ ( phi ) , 0 ) \begin{aligned} \bold{r}_{\theta} &=(-R\sin(\theta)\cos(\phi), -R\sin (\theta)\sin(\phi), R\cos(\theta)) \\\bold{r}_{\phi} &=(-R\cos(\theta)\sin(\phi), R \cos(\theta)\cos(\phi), 0) \end{对齐} rθ r = (R1(θ)cos(phi), Rθ罪 (phi),Rcos(θ))=(−) Rcos(θ) sin(phi),Rcos (θ)cos(Φ) ),0) I = R 2 ( d θ d θ + cos ⁡ 2 ( θ ) d Φ d Φ ) I = R^2 (d \theta d\ θ + \cos^2(\theta)d\phi d\phi) I=R2(d θdθ +cos2θ) d φdφ


这是一个住宿的好地方

定义4.1 令曲面 S S S 的参数表示为 r = r ( u , v ) \bold{r}=\bold{r }(u,v)r=ruv) ,则 n = r u ∧ r v ∣ r u ∧ r v ∣ \bold{n}=\frac{\bold{r}_u \wedge \bold{r}_v}{|\bold{r}_u \wedge \bold{r}_v|} n=ru∧ rvru rv S S S 单位法线向量,曲面 S S S 的第二种基本形式定义为 I I = − < d r , d n > II=- II =< dr,dn> by n \粗体{n} n 定义为 < r u , n > = 0 <\bold{r}_u, \;\bold{n}>=0 <ru, n>=0 , < r v , n > = 0 <\bold{r}_v, \;\bold{n}> =0 <rv,n >=0 ,这两个方程的导数给出 L ≜ < r u u , n > = − < r u , n u > M ≜ < r u v , n > = − < r u , n v > = − < r v , n u > N ≜ = − \begin{对齐} L &\triangleq <\bold{r}_{uu}, \;\bold{n}> = -<\bold{r}_u, \;\bold{n}_u> \\ M&\triangleq <\bold{r}_{uv}, \;\bold{n}> = -<\bold{r}_u, \;\bold{n}_v> =-<\bold{r}_v, \;\bold{n}_u> \\ N &\triangleq <\bold{r}_{vv}, \;\bold{n}> = -<\bold{r}_v, \;\bold{n}_v> \end{对齐} LMN<ruu,n >=<ru,nu>≜ <ruv,n >=<ru,nv>=<rv,nu><rvv,n>=<rv,n v>第二基本形式可表示为 I I = − < d r , d n > = − < r u d u + r v d v , n u d u + n v d v > = L d u d u + 2 M d u d v + N d v d v \begin{对齐} II &= - \\ &=-<\bold{r}_udu+\bold{r}_vdv, \;\bold{n}_udu+\bold{n}_vdv> \\ &= Ldudu+2Mdudv+Ndvdv \end{对齐} II= <dr,dn> =<rudu+rvd v,nu du+nvd v>=Ldud u+2Mdudv+Ndvdv ​ 第二种基本形式实际上反映了表面的形状。

例4.2讨论平面和圆柱的第二种基本形式。

平面 r ( u , v ) = ( u , v , c ) \bold{r}(u,v) = (u,v,c) r() u,v)=(u,v,c) , n = ( 0 , 0 , 1 ) \bold{n}=(0,0,1) n=() 、01 , I I = − < d r , d n > = 0 II= -=0 II = −<dr,dn>=0

圆柱面 r ( u , v ) = ( x ( u ) , y ( u ) , v ) \bold{r}(u,v)=(x(u), y(u), v) r(u,v)=(x(u) ,y(u), ( x ( u ) , y ( u ) ) (x(u),y(u)) (x(u),y (u)) 为平面曲线, u u u 请检查曲线的弧长, r u u = ( x u u , y u u , 0 ) , r u v = ( 0 , 0 , 0 ) , r v v = ( 0 , 0 , 0 ) , n = ( y u , − x u , 0 ) \begin{对齐} &\bold{r}_{uu} =(x_{uu},y_{uu}, 0), \\ &\bold{r}_ {uv} =(0,0,0), \\ &\bold{r}_{vv} =(0,0,0), \\ &\bold {n} = (y_u, -x_u,0) \end{对齐} ruu=( xuu,yuu,0,ruv = (0,0,0 ),rvv=( ,0,0),n= (yu,− xu ,0) κ \kappa κ 是一个曲面 ( x ( u ) , y ( u ) ) (u)) 定义,在 κ = − < r u u , n > = − x u y u + x u y u u \kappa=-<\bold{r}_{uu}, \; \bold{n}>=-x_ {uu}y_u+x_uy_{uu} − uu ,n>=x x = − κ , M = N = 0. L. =<\bold{r}_{uu}, \;\bold{n}>=-\kappa , \quad M=N=0。 L=<ruu, n>=k, =N =0. 热电偶中随机分散的当量 I I = − k d u d u II=-\black kappa 0 =κdudu 尤其是 ( x ( u ) , y ( u ) ) ( x ( u ) , y(u) )半径 R R R κ = 1 R \kappa=\fra c{1}{R} κ =R1 ,圆柱体的第二种基本形式是 I I = − 1 R d u d u II = -\frac{1}{R}dudu II =R1dudu 以上结果表明,虽然平面和圆柱的第一基本形式是相同的,第二种基本形式则不同,表现为不同的形状。

示例 4.3 很容易验证半径为 R R R 的球的第二基本形状是 I I = R ( d θ d θ + cos ⁡ 2 ( θ ) d phi d phi ) II=R(d\theta d\theta+\cos^2(\theta)d\phi d\phi) II=Rdθdθ+cos2(θ)ddφ ) 其中, n \bold{n} n取为 ( − cos ⁡ ( θ ) cos ⁡ ( phi ) , − cos ⁡ ( θ ) sin ⁡ ( phi ) , − sin ⁡ ( phi ) ) (-\cos(\theta)\cos(\phi), -\cos(\theta)\sin(\phi), -\sin( \phi)) ( cos(θ) cos(phi),cos(θ)sin( φ) ,sin ()

性质 4.4 在本节末尾给出。

性质 4.4 曲面的第二基本形式 I I = L d u d u + 2 M d u d v + N d v d v II = Ldudu+2Mdudv+Ndvdv II=Ldudu+2M d udv+Ndvdv 正定点或负定点,即 L N − M 2 > 0 LN- M^2>0 L N− 在点M2>0 附近,曲面形状为凸面(或凹面,由法线选择决定)。在 L N − M 2 < 0 LN-M^2< 0 LNM <0 靠近点,表面呈马鞍形的。


5。法向曲率和 Weingarten 变换

我们在没有证据的情况下定义了表面的法向曲率。

曲面 S S S P 0 P_0 P0 ω = xi r u + η r v \ omega=\xi \bold{r}_u \eta \bold {r}_v ω=editru​​ + ηrv​​ 法曲率 定义为 κ n ( ω ) = I I ( ω , ω ) = L Σ 2 + 2 M Σ 2 + N 2 η 2 E ψ 2 + 2 F ψ η + G 2 η 2 \kappa_n(\omega) = \frac{II(\omega,\omega)}{I(\omega, \omega)} = \frac{L\xi^2+2M ) ω)=I(ω) )II())​​= E编辑 + 4 N2η2​​

例5.1 R R R 的球的法曲率至 1 R \frac{1}{R} R R .

这是因为 I I = I ⋅ 1 R II = I \cdot \frac{1}{R} II= I R1。这意味着球体的曲率在任何方向上都是相同的。

例5.2讨论平面和圆柱体的法向曲率。

平面法向量是常数向量,所以 d n = 0 d\bold{n}=0 dn= 0 , I I = 0 II =0 II=0

定义 r(u, v) = (R cos ⁡(u R), R sin ⁡(u R), v) \bold{r}(u,v)=\left(R \cos \left( \frac{u}{R} \right), R\sin\left( \frac{u}{R}\right), v \right) u,v) =(Rcos(R u2),Rsin(Ru2)2,v) I = d u d u + d v d v I = dudu+dvdv I=d udu+ d vdv I I = − 1 R d u d u II=-\frac{1}{ . R}嘟嘟II=R1d udu结合 r u \bold{r }_u ru r v \bold{r}_v r v 如果你有雪花 ω \omega m ω θ \theta θ 之一 ) ω \omega ω r u \bold{r}_u ru ω = cos ⁡ ( θ ) r u + sin ⁡ ( θ ) r v \omega=\cos(\theta)\粗体{r}_u+\sin(\theta)\bold{r}_vω=cos(θ )ru +(θ)rv κ n ( ω ) = L cos ⁡ 2 ( θ ) + 2 M cos ⁡ ( θ ) sin ⁡ ( θ ) + N sin ⁡ 2 ( θ ) E cos ⁡ 2 ( θ ) + 2 F cos ⁡ ( θ ) sin ⁡ ( θ ) + G sin ⁡ 2 ( θ ) = − 1 R cos ⁡ 2 ( θ ) cos ⁡ 2 ( θ ) + sin ⁡ 2 ( θ ) = − 1 R cos ⁡ 2 ( θ ) \begin{aligned} \kappa_n(\omega ) &= \frac{L\ cos ^2(\theta)+2M\cos(\theta)\sin(\theta)+N\sin^2(\theta)} {E\cos^2(\theta); +2F\cos(\theta)\sin(\theta)+G\sin^2(\theta)} \\ &=\frac{-\frac{1}{R}\cos^2(\theta)} {\cos^2(\theta)+\sin^2(\theta)} \\ &= -\frac{1}{R}\cos^2(\theta) \end{对齐} κn (ω) = Ecos2(θ) + Fcos(θ)sin (θ)+G2(θ)Lcos2(θ) )+2Mcos( θ )sin(θ)+Nsin2(θ)=cos2( θ)+sin2(θ) −R1cos2θ )=R1cos 2(θ) 以下定义高斯映射。

定义 5.3 令曲面 S S S 的参数表示为 r = r ( u , v ) \bold{r}=\bold{r }(u,v)r=ruv) ,单位法向量为 n ( u , v ) \bold{n}(u,v) n(u,v ) ,记录单位球面作为s 2 s^ 2S2,称为映射g:s→s 2,r(u,v)↦n(u,v),g:\ quad s \ to S^2 , \\ \bold {r}(u,v) \mapsto \bold{n}(u,v) , g:S→ S2,r(u,v )n(u,v) , 是表面 的 高斯图 SSS

高斯映射 g g g 沿曲线 ( u ( t ) , v ( t ) ) (u(t), v(t)) ( u(t),v( t)) 的导数为 d n ( t ) d t = n u d u d t + n v d v d t \frac{d\bold{n}(t)}{dt} = \bold{n}_u \frac{du}{dt} + \bold{n}_v \frac{dv}{dt } dtdn(t )=nudtdu+ nvdtdv < d n ( t ) d t , n > = 0 <\frac{d\bold{n}(t)}{dt} , \; \bold{n}> =0 <dtdn() t ),n>=0,我们得到 d n ( t ) d t \frac{d\bold{n}(t)} {dt} dtdn(t ) ​ S S 的切向量S,所以 n v \b old{n}_v nv 都是切向量。这表明 ( d u d t , d v d t ) → d n d t = n u d u d t + n v d v d t \left( \frac{du}{dt}, \frac{dv}{dt} \right) \to \frac{d\bold{n}}{dt} = \bold{n}_u \frac{du}{dt} +\bold{n}_v \frac{dv}{dt} (dtdu,dtdvd tdn=nudtdu+nvdt dv是切向量之间的对应关系。因此,可以定义从一个切平面到另一个切平面的线性变换 W : T P S → T P S , v = λ r u + μ r v ↦ W ( v ) = − ( λ n u + μ n v ) \begin{aligned } \mathfrak{W}: \quad &T_PS \to T_PS, \\ &\bold{v}=\lambda\bold{r}_u+\mu\bold{r}_v \mapsto \mathfrak{W}(\bold{ v })=-(\lambda\bold{n}_u+\mu\bold{n}_v) \end{对齐} W:T P STPS ,v=λr u+μrvW(v)=(λnu​❓ +μnv W \mathfrak{W} W 称为 Weingarten 变换

Weingarten 变换是将切向量 v \bold{v} v 以切平面 为基础进行变换{ r u , r v } \{ \bold{r}_u , \; \bold{r}_v \} {ru,r v} 转换为 { n u , n v } \{ \bold{n}_u, \; \bold{n}_v \} {nun v} ,系数变为原系数的相反数。

性质5.4 Weingarten变换与曲面的参数选择无关。

性质 5.5 到曲面的任何单位切向量 0 S S S 21 v \bold{v} v ,表面 S S S 沿 v \bold{v} v 方向的法向曲率是 κ n ( v ) = < W ( v ) , v > \kappa_n(\bold{v} ) = <\mathfrak{W}(\bold{v}), \; \bold{v}> κn(v )=<W(v) ,v> 证明 让我们假设 r u \bold {r}_u ru r v \bold{r}_v rv 是单位正交向量(可以通过施密特正交化得到),设 v = λ r u + μ r v \bold{v} = \lambda \bold{r}_u+\mu \bold{r}_v v=λru + μrv。那么 < W ( v ) , v > = − < λ n u + μ n v , λ r u + μ r v > = λ 2 L + 2 λ μ M + μ 2 N = κ n ( v ) \begin{aligned} <\mathfrak{W(\bold{v})}, \; \bold{v}> &= - \\ &= \ lambda^2L+2\lambda\mu M +\mu^2N \\ &= \kappa_n(\bold{v}) \end{对齐} <W(v) ),v>=<λnu+μnv,λ ru​v M+μ2N=κn​​(v) among them. . u G =<rv​​,rv​​>= , λ 2 + μ 2 = 1 \lambda^2+\mu^2=1 λ2+μ2=1

定理5.6 Weingarten transformation是曲面切电影的这些的自共业电影,即 ∀ v \forall \; \bold{v} v, w ∈ T P S \bold{w} \in T_PS wT_PS ​S , < W ( v ) , w > = < v , W ( w ) > = <\bold{v}, \; \mathfrak{W(\bold{w})}><W () w>=<v,W(w) >, Weingarten has two characteristics of transformation.值的实数。


6. 主曲率和高斯曲率

S ,令 κ \kappa κ 为 Weingarten 变换的特征值, v \bold{v} v 是相应的单位特征向量,由 给出 < W ( v ) , v > = = κ <\mathfrak{W(\bold{v})}, \; \bold{v}> =<\kappa \bold{v}, \; \bold{v}> = \kappa < W(v),v>=<κv,v>= κ属性 5.5已知, κ \ kappa κ 是曲面 v \bold{v} v 的法向曲率。我们将 P P P处的Weingarten变换的两个特征值称为曲面 S S S 位于 P P P 主曲率。特征值对应的两个实特征向量称为点处曲面的主方向 P P P 。

当两个主曲率不相等时( W \mathfrak{W} W有两个不同的特征值),对应的两个主方向完全确定且互为正值。支付。当两个主曲率相等时,不能唯一确定主方向。此时,曲面在该点的任意切线方向就是主方向(后面会解释)。

要计算曲面的主曲率,首先要求出坐标切向量下Weingarten变换的系数矩阵。令 S S S 的参数表示为 r ( u , v ) \bold{r}(u,v) r (u ,v) , 基于切平面 { r u , r v } \{\bold{r}_u , \; \bold{r}_v\} {ru,r v} 系数矩阵Weingarten 变换的值设置为 ( a b c d ) \left ( \begin{array}{cc} a & b \\ c & d \end{array} \right) (acb d然后我们有 W ( r u ) = − n u = a r u + br r v W ( r v ) = − n v = c r u + d r v \begin{对齐} &\mathfrak{W (\bold {r}_u)}=-\bold{n}_u=a \bold{r}_u+b \bold{r}_v \\ & \mathfrak{W(\bold{r}_v)}=-\bold {n}_v=c \bold{r}_u+d \bold{r}_v \end{对齐} W (ru)= −nu= aru+brv ​W(rv)=nv​ =cru+drv < − n u , r u > = a + b < r v , r u > <-\bold{n}_u, \; \bold{r}_u>=a<\bold{r}_u, \; \bold{r}_u>+b<\bold{r}_v, \; \bold{r}_u> nu,ru>=a<r1u,ru>+b<rv,ru>增益 L = a E + b F L= a E+ b F L=aE+bF于 < − n u , r v > = a < r u , r v > + b < r v , r v > <-\bold{n}_u, \; \bold{r}_v>=a<\bold{r}_u, \; \bold{r}_v>+b<\bold{r}_v, \; \bold{r}_v><nu,rv>=a< ru,rv ​>+b<r v,rv > 增益 M = a F + b G M=aF+bG M=aF+bG 余 < − n v , r u > = c < r u , r u > + d < r v , r u > <-\bold{n}_v, \; \bold{r}_u>=c <\bold{r}_u, \; \bold{r}_u> +d<\bold{r}_v, \; \bold{r}_u> <nv ​,ru>= c <ru, ru>+d<rv,ru> M = c E + d F M =cE+dF M=cE+d F < − n v , r v > = c < r u , r v > + d < r v , r v > <-\bold{n}_v, \; \bold{r}_v>=c<\bold{r}_u, \; \bold{r}_v>+d<\bold{r}_v, \; \bold{r}_v> <− nv,rv>=c<ru​ ,rv>+ d<rv, rv> N = c F + d G N=cF+dG N=cF+ dG联立上面两个方程组,求得 a = L G − M F E G − F 2 , b = ME − L F E G − F 2 c = MG − N F E G − F 2 , d = N E − M F E G − F 2 \begin{aligned} a=\frac{LG-MF} {EG-F^2},\quad & b=\frac{ME-LF}{EG-F^2} \\ c=\frac{MG-NF}{EG-F^2},\quad & d =\frac{NE-MF}{EG-F^2} \end{对齐} a=EGF2 LGMF,c=EG −F2MGNF,b=EGF2MELFd=EGF2NEMF因此 Weingarten 变换的矩阵为 ( a b c d ) = ( L M M N ) ( E F F G ) − 1 = 1 E G − F 2 ( LG − M F M E − L F M G − N F N E − M F ) \begin{aligned} \left( \begin{array}{cc} a & b \\ c & d \end{array} \right) &= \left( \begin{array}{cc} L & M \\ M & N \end{array} \right) \left( \begin{array}{cc} E & F \\ F & G \end{array} \right) ^{-1} \\ &= \frac{1}{EG-F^2} \left( \begin{array}{cc} LG-MF& ME-LF \\ MG-NF & NE-MF \end{array} \right) \end{对齐} (acbd)=(LMMNEFFG)1=EGF21​ (LGMFMGNFMELFNEMF )归因于主曲率 κ \kappa κ系上述矩阵的特征值,故有特征方程 ∣κ + M F − L G E G − F 2 L F − M E E G − F 2 N F − M G E G − F 2 κ + M F − N E E G − F 2 ∣ = 0 \begin{aligned} &\left| \begin{array}{cc} \kappa+\frac{MF-LG}{EG-F^2}& \frac{LF-ME}{EG-F^2} \\ \\ \frac{NF-MG} {EG-F^2} & \kappa + \frac{MF-NE}{EG-F^2} \end{array} \right| =0 \end{对齐} κ+EGF2MFLGEGF2NFMGEGF2LFMEκ+EGF2MFNE =0 展开即得 κ 2 − L G − 2 M F + N E E G − F 2 κ + L N − M 2 E G − F 2 = 0 \kappa^2-\frac{LG-2MF+ ne} {eg-f^2} \ kappa + \ frac {ln-m^2} {eg-f^2} = 0κ2-g--f 2LG2MF+NE ​κ+EGF2LNM2=0记断的两个主曲 κ 1 \kappa_1 κ1 κ 2 \kappa_2 κ2 ,则称为 H = 1 2 ( κ 1 + κ 2 ) H= \frac{1}{2}(\kappa_1+\kappa_2) H=21(κ1+κ2) 表面的平均曲率 , K = κ 1 ⋅ κ 2 \Kappa= \kappa_1 \ cdot \kappa_2 K= κ 1κ2被称为表面的高斯曲率。根据吠陀定理, H = 1 2 L G − 2 M F + N E E G − F 2 , K = L N − M 2 E G − F 2 = ∣ a b c d ∣ \begin{aligned} H &=\frac{1}{ 2}\frac{LG-2MF+NE}{EG-F^2}, \\ \Kappa &= \frac{LN-M^2}{EG-F^2} = \left| \begin{array} {cc} a & b \\ c & d \end{array} \right| \end{对齐} HK=21EG F2LG2MF+NE,= EGF2LNM2=acb d ∣ ​ 因为 n u = − a r u − b r v \bold{n} _u=-a \bold{r}_u-b\bold{r}_v nu=arub rv n v = − c r u − d r v \bold{n}_v=-c \bold{r}_u- d \bold{r}_v nv=cr udrv ,得 n u ∧ n v = a d r u ∧ r v + b c r v ∧ r u = ( a d − b c ) r u ∧ r v = K r u ∧ r v \begin{对齐} \bold{n}_u \楔形\bold{n}_v &= 广告\; \bold{r}_u \楔形 \bold{r}_v +bc\; \bold{r}_v \wedge \bold{r}_u \\ &=(ad-bc)\; \bold{r}_u \wedge \bold{r}_v \\ &=\Kappa \; \bold{r}_u \wedge \bold{r}_v \end{对齐} nun v=adrurv+bcrvru=(adbc)ru​ ∧rv=Krurv高斯曲率稀疏的曲面是可展曲曲面,平均曲率为零的曲面是,最小曲面

让我们考虑一下法曲率和主曲率之间的关系。设 P ∈ S P\in S PεS ,取 e 1 , e 2 e_1, e_2 e1, e2 是 中的曲面 S S S 在 P P P 该点的主方向,以及 { e 1 , e 2 } \{ e_1, e_2 \} {e1,e ​ } 成分 T P S T_PS T PS 的单位正交基(注意 W \mathfrak{W} W 的定义域是 T P S T_PS TPS),主曲率的定义为 W ( e i ) = κ i e i \mathfrak{W }(e_i)= \ kappa_ie_i W (ei) =κiei , i = 1 , 2 i=1,2 i=1,2,对于任何单位向量 v ∈ T P S \bold{ v}\in T_PS vεTPS ,设 v = cos ⁡ ( θ ) e 1 + sin ⁡ ( θ ) e 2 \ 粗体{v}=\cos(\theta)e_1+\sin(\theta)e_2 v=cos( θ)e1+θ) )e2, θ \theta θ v \bold{v} v 一 e 1 e_1 e1 娱乐然后 κ n ( v ) = = = κ 1 cos ⁡ 2 ( θ ) + κ 2 sin ⁡ 2 ( θ ) \begin{aligned} \kappa_ n(\bold{v}) &= <\mathfrak{W}(v), \; v> \\ &= \\ &=\kappa_1 \cos^2(\theta) + \kappa_2 \sin^2(\theta) \end{对齐} κn (v)=<W (v),v >=<cos(θ)κ e1+θ)κ2e2, cos(θ)e 1 +(θ)e2>=κ1cos2(θ)+κ2罪 (θ)娱乐更新简介 版本 6.1 输入 κ 1 \kappa_1 κ1 κ 2 \kappa_2 κ 2 P P P 版权所有, 和 1 and_1 e2 是对应的正交主方向​​,令 v ∈ T P S \bold{v} \in T_PS vεTPS 是单位向量, v \bold{v} v e 1 e_1 e1 角度为 θ \ theta θ ,则曲面沿 位于 P P P 点 v \b old{v} v 方向的法向曲率是 κ n ( v ) = κ 1 cos ⁡ 2 ( θ ) + κ 2 sin ⁡ 2 ( θ ) \kappa_n(\bold{v}) =\kappa_1 \cos^2(\theta) + \kappa_2 \sin^2(\theta ) κn(v) =κ1cos2(θ)+κ2 sin2( θ) 由上式可知当主曲率 κ 1 \kappa_1 当 κ1 等于 κ 2 \kappa_2 时κ2 速率和切线方向 v \bold{v} v 没关系。当主曲率 κ 1 \kappa_1 κ1 κ 2 \kappa_2 κ2 不相等时,根据第二个足够值条件已知的极值,当 θ = 0 \theta=0 θ=0 θ = π 2 \theta =\frac{\pi }{2} θ 当=2π时, κ n ( v ) \kappa_n(\bold{v} ) κn(v)同时取最大值和最小值(或同时取最小值和最大值)。到这里我们已经解释了之前的问题了。

我们来考虑一下高斯图图像的面积和曲面的面积。

假设曲面 S S S 的面积元素为 d A dA dA,回忆以下混合产品 结论: a ⋅ ( b ∧ c ) = b ⋅ ( c ∧ a ) = c ⋅ ( a ∧ b ) a \cdot (b \wedge c) = b \cdot (c \wedge a) = c \ cdot (a \楔 b) a(b∧ c)=b (ca)=c(ab )拉格朗日公式 a ∧ ( b ∧ c ) = ( a ⋅ c ) b − ( a ⋅ b ) c a \楔形 (b \楔形 c) = (a \cdot c)b - (a\cdot b)c a(b∧ c)=(acb− (ab)c r u \bold{r}_u ru r v \ 粗体{r}_v rv 如上所述, d A = ∣ r u ∧ r v ∣ d u d v = ( r u ∧ r v ) ⋅ ( r u ∧ r v ) d u d v = r u ⋅ ( r v ∧ ( r u ∧ r v ) ) d u d v = [ ( r v ⋅ r v ) r u − ( r u ⋅ r v ) r v ] ⋅ r u d u d v = E G − F 2 d ud v \begin{对齐} dA &= \left| \bold{r}_u \楔形 \bold{r}_v \right| dudv \\ &=\sqrt{(\bold{r}_u \wedge\bold{r}_v)\cdot (\bold{r}_u \wedge\bold{r}_v)} \;dudv \\ &= \sqrt{\bold{r}_u \cdot(\bold{r}_v \wedge (\bold{r}_u \wedge\bold{r}_v))} \; dudv \\ &=\sqrt{[(\bold{r}_v \cdot \bold{r}_v)\bold{r}_u - (\bold{r}_u \cdot \bold{r}_v)\bold {r}_v]\cdot \bold{r}_u} \; dudv \\ &=\sqrt{EG-F^2} \; dudv \end{对齐} dA=r urv dudv =(rurv)(r urv)d udv=ru(rv(rur v))21dudv=[() rvrv)ru(rurvr v]rudu d v=EGF2 dudv S S S u u u u + d u u+du u+du v v v到达 v + d v v+dv v+dv 常吕小平行四足形状 Δ \Delta Δ下、 g ( Δ ) g(\Delta) g (Δ) d σ = ∣ ( n ( u + d u , v ) − n ( u , v ) ) ∧ ( n ( u , v + d v ) − n ( u , v ) ) ∣ ≈ ∣ n u ∧ n v ∣ d u d v \begin{aligned} d\sigma &= \left| (\bold{n}(u+du, v)-\bold{ n}(u,v)) \wedge (\bold{n}(u, v+dv)-\bold{n}(u,v) )) \右| \\ &\大约\左| \bold{n }_u \楔形 \bold{n}_v \right| dudv \end{对齐ned} dσ=( n(u+du ,v)n( u,v))() n(u,v+ dv)n( u,v))nunv dudv ​我们可以获取后面这个估计的面积,即 d σ = ∣ n u ∧ n v ∣ d u d v = K ∣ r u ∧ r v ∣ d u d v = K d A d\sigma = \left| \bold{n}_u \楔形 \bold{n}_v \right| dudv = \Kappa \left| \bold{r}_u \楔形 \bold{r}_v \right| dudv = \Kappa dA dσ=nunvdudv= Kru rvdud v=KdA其中 K \Kappa K 为高斯曲率。取 D D D 为 S S S 上的含 PP P, g ( D ) g(D) g() 的面积 D) D D D 就像高斯映射下的曲面,那么它的面积为 Area a ( g ( D ) ) = ∫ g ( D ) d σ = ∫ D K d A Area(g(D) ) = \int_{ g(D)} d\sigma = \int_{D} \Kappa dA Area( g (D )) =g(D )dσ=DKdA当时 D → P D \到 P 当 DP时, lim ⁡ D → P 区域 ( g ( D ) ) 区域 ( D ) = lim ⁡ D → P ∫ D K d A ∫ D d A = K ( P ) \lim_{D \to P} \frac{面积(g(D) )}{面积(D)} =\lim_{D \to P} \frac{\int_{D} \Kappa dA}{\int_{D} dA}= \Kappa(P) DPlimAre a (D)Area(g(D)) ​ =DPlimDdADKd A=K(P) 上一页式子的最后一个方程可以通过第一类的曲面积分中值定理直接推导出来。积分项 ∫ D K d A \int_{D} \Kappa dA DKdA还有更多come 出现在高斯-邦内定理

这就是高斯曲率的几何意义。 高斯映射图像表面在某一点的曲率度正是该点处的表面高斯曲率。

至此,我们已经介绍了曲面的局部理论,接下来我们将进行曲率的Python编程计算。


7。 Python中的高斯曲率计算

说明:数学基础一旦打好,后续的编程往往是最容易实现的。

7.1 微分运算

先回顾一下如何实现求导

来自 sympy 导入符号差异#自己定义一个多元函数
def 功能(x, y):返回 y*  x**2 + 2*y  **2# 首先将需要的变量(x,y)符号化。 
x, y =符号('x y' ,真实=真实) print("求x的偏导数:")打印(差异(功能() x, y),2 ))
打印("求y的偏导数")
打印 (iDiff (功能 (x , y) , y))
打印("求x,y的二阶混合偏导数")
打印 (iDiff (功能 (x , y) , x,y))
打印("求x的二阶偏导数:")
打印(差异(功能() x, y),2 ,2))
打印("求y的二阶偏导数:")
打印(差异(功能() ,2))

输出以下结果

 求 x 的偏导数:
2*x*y
求 y 的偏导数x**2 + 4* y
求 x, y 的二阶混合偏导数
2*x
求 x  的二阶偏导数:
2*y
求 y  的二阶偏导数:
4

7.2 球体的高斯曲率

基于前面六章的内容,我们现在可以编写常见的曲面曲率计算的Python代码了。球面曲率相对容易实现。

#计算半径为R的球体曲率import sympy
来自 sympy 导入符号 diff,简化#球体的度量张量g11
defg11(theta, phi):R =符号   ('R',真实=) return R**2# 球体的公制张量g22
defg22(theta, phi):R =符号   ('R',真实=) 返回R**2*(sympy.theta) ) **2# 度量矩阵的行列式defg(theta,phir ==符号  ('R',真实=)返回R**4*( sympy.theta) **2#曲率(此曲率公式已简化,不能用于一般曲面情况)defK(theta, phi):返回( / (2*g(theta, *()  -差异g11theta) , phi), phi, 2)-差异(g22() theta, phi),theta  ,2))\-() g22(theta, phi)  /(4*g(theta, **2 *(差异(g11(theta) , phi),theta)*(-diff(g22) (theta, phi),  theta))-diff(g11) (theta, phi), phi)**2)\- (g11(theta, phi)/(4*g(theta) ,phi**2*diff theta, phi), phi) *(-diff(g11) (theta, phi), phi))-diff(  g22(theta, phi) ,theta)**2)#先将所求变量(theta,phi)符号化theta, phi =符号 ('theta phi',真实=True)#计算半径为R的球面的曲率
简化(K(theta, phi))

输出结果为 1 R 2 \frac{1}{R^2} R21 ,说明球面具有恒定常值的正曲率。

事实上,球体的法向曲率之前已经给出为 1 R \frac{1}{R} R1,因为它的法向曲率等于两个主曲率,因此其高斯曲率也应为 1 R ⋅ 1 R = 1 R 2 \frac{1}{R}\cdot \frac{1}{R}=\frac{1}{R ^2} R1R1=R21

7.3 双曲抛物线的高斯曲率

#双曲抛物面曲率计算 z=xyimport sympy
来自 sympy 导入符号 diff,简化#双曲抛物面g11的公制张量
defg11(u, v):返回1+  v **2#双曲抛物面的公制张量 g22
defg22(u, v):返回1+  u **2# 双曲抛物面的公制张量 g12
defg12(u, v):返回u*   v# 度量矩阵的行列式defg(u, ? g11 (u,v)*g22(u,v) -g12(u , v )**2)#曲率(这里的曲率公式为一般情况) defK(theta, phi):返回( / (2*g(theta, *()  *差异(g12(theta , phi),theta,phi )-差异(g11() theta, phi), phi  , 2)-diff(g22) (theta, phi),theta,2))\ -(g22(theta,)/(4*g() theta, phi)**2))*(diff(g11) (theta, phi),  theta)*(2*diff() g12(theta, phi), phi)-diff( g22(theta, phi) ,theta))-diff (g11(theta, phi),)**2)\ -(g11(theta, phi)/(4*g  (theta, phi)**2) )*(差异() g22(theta, phi), phi)*(2*diff(g12(theta, phi ),theta)-diff  (g11(theta, phi), phi))-diff(g22(theta, phi ),theta)**2)   \+(g12(theta,  phi)/(4*g (theta, phi)**2)*((diff() g11(theta, phi)  ,theta))*(diff) (g22(theta, phi) ),))- *(差异(g11()  theta, phi), phi) )*(差异(g22) (theta, phi),theta))+(2*  diff(g12(theta, phi), theta)-diff (g11(theta, phi),))*  (2*差异(g12(the ta, phi), phi) -差异g22theta) , phi),theta)  ) )#先将所求符号(theta,phi)符号化u, v =符号('u v',真实=)#计算双曲抛物面的曲率z=xy
sympy因子简化K) (u, v)))

输出结果为 − 1 ( u 2 + v 2 + 1 ) 2 -\frac{1}{\left(u^{2} + v^{2} + 1\right)^{ 2 }} (u2+v2+ )21 ,这是一个具有恒定但不恒定的鞍面恒定的负曲率。

7.4 锥面高斯曲率

#圆锥体曲率计算 z^2 = x^2 + y^2import sympy
来自 sympy 导入符号 diff, 简化# 双曲抛物面的度量张量 g11 v = z u=theta 
defg11u,vreturnv**  # 双曲抛物面的公制张量 g22
defg22(u, v):返回 2# 双曲抛物面度量张量g12defg12(u, v):返回 0# 度量矩阵行列式 
defg(u, ? g11 (u,v)*g22(u,v) -g12(u , v )**2)#曲率defK(theta, phi):返回( / (2*g(theta, *()  *差异(g12(theta , phi),theta,phi )-差异(g11() theta, phi), phi  , 2)-diff(g22) (theta, phi),theta,2))\ -(g22(theta,)/(4*g() theta, phi)**2))*(diff(g11) (theta, phi),  theta)*(2*diff() g12(theta, phi), phi)-diff( g22(theta, phi) ,theta))-diff (g11(theta, phi),)**2)\ -(g11(theta, phi)/(4*g  (theta, phi)**2) )*(差异() g22(theta, phi), phi)*(2*diff(g12(theta, phi ),theta)-diff  (g11(theta, phi), phi))-diff(g22(theta, phi ),theta)**2)   \+(g12(theta,  phi)/(4*g (theta, phi)**2)*((diff() g11(theta, phi)  ,theta))*(diff) (g22(theta, phi) ),))- *(差异(g11()  theta, phi), phi) )*(差异(g22) (theta, phi),theta))+(2*  diff(g12(theta, phi), theta)-diff (g11(theta, phi),))*  (2*差异(g12(the ta, phi), phi) -差异g22theta) , phi),theta)  ) )#先将所求符号(theta,phi)符号化u, v =符号('u v',真实=)#计算圆锥的曲率 z^2 = x^2 + y^2 
sympy因子简化K) (u, v)))

的输出结果为 0 0 0 ,即锥面为可展曲面,平坦如平面。

至于其他面,你可以自己重写代码,这里不再赘述。值得一提的是,莫比乌斯带是一个二维不可定向曲面。它的高斯曲率是不确定。其原因还需要进一步研究,本文不再解释。 。

相关文章