从方程的角度求解曲面上的测地线,是我们求测地线的常用方法。但是这种方法在求一些具体问题时计算量很大,同时也不利于人们对于曲面上测地线的直观想象。本文从测地曲率几何意义的角度出发,求出了椭圆抛物面、旋转抛物面和双曲抛物面上经过原点的所有平面曲线类型的测地线,求解过程简单,有助于人们直观理解,并利用数学软件Mathematica画出了测地线的图像。 Solving geodesics on surfaces from the perspective of equations is a common method we use to find geodesics. However, this method is very computationally intensive when solving some specific problems, and it is also not conducive to people’s intuitive imagination of geodesics on surfaces. In this paper, from the perspective of the geometric meaning of geodesic curvature, the geodesic of all plane curve types passing through the origin on elliptic paraboloid, rotating paraboloid and hyperbolic paraboloid are obtained, and the solution process is simple and helpful for people to intuitively understand, and use the mathematical software Mathematica to draw the geodesic image.
从方程的角度求解曲面上的测地线,是我们求测地线的常用方法。但是这种方法在求一些具体问题时计算量很大,同时也不利于人们对于曲面上测地线的直观想象。本文从测地曲率几何意义的角度出发,求出了椭圆抛物面、旋转抛物面和双曲抛物面上经过原点的所有平面曲线类型的测地线,求解过程简单,有助于人们直观理解,并利用数学软件Mathematica画出了测地线的图像。
测地曲率,抛物面,测地线
Yuying Wu, Tuya Bao*
College of Mathematics and Physics, Inner Mongolia Minzu University, Tongliao Inner Mongolia
Received: Oct. 23rd, 2022; accepted: Nov. 22nd, 2022; published: Nov. 30th, 2022
Solving geodesics on surfaces from the perspective of equations is a common method we use to find geodesics. However, this method is very computationally intensive when solving some specific problems, and it is also not conducive to people’s intuitive imagination of geodesics on surfaces. In this paper, from the perspective of the geometric meaning of geodesic curvature, the geodesic of all plane curve types passing through the origin on elliptic paraboloid, rotating paraboloid and hyperbolic paraboloid are obtained, and the solution process is simple and helpful for people to intuitively understand, and use the mathematical software Mathematica to draw the geodesic image.
Keywords:Geodesic Curvature, Parabolic Surface, Geodesic
Copyright © 2022 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
测地线又称为短程线,在数学上可以理解为平面上的直线在曲面上的推广,用于描述曲面上两点之间的最短曲线,是微分几何中的重要概念。它不仅在数学内部有着十分重要的作用,在物理学、工程学和大地测量学等方面也有广泛的应用。通过曲面上测地线的形状和性质,可以分析曲面的结构和性质。比如,球面上的测地线必在球面的大圆上 [
定义1 [
下面我们给出测地曲率的几何意义。
命题1 [
直线的曲率为0 [
引理1 [
π 1 : A 1 x + B 1 y + C 1 z + D 1 = 0 , π 2 : A 2 x + B 2 y + C 2 z + D 2 = 0
交于一条直线L,那么以直线L为轴的有轴平面束的方程是
l ( A 1 x + B 1 y + C 1 z + D 1 ) + m ( A 2 x + B 2 y + C 2 z + D 2 ) = 0
其中 l , m 是不全为零的任意实数。
引理2当两平面垂直时,其中一个平面上的曲线在另一个平面上的正投影曲线为直线,否则为曲线。
证:当两平面垂直时(见图1), π 1 上的任意点正投影到两平面的交线上,所以 π 1 上任意曲线在 π 2 上的正投影为直线。当两平面不垂直时(见图2), π 1 上非直线的任意曲线在 π 2 上的正投影都不是直线。
图1. 两平面垂直
图2. 两平面不垂直
曲面上如果存在直线,则此直线一定是测地线 [
定理1椭圆抛物面 r = { a u cos v , b u sin v , 1 2 u 2 } 上经过原点的非直线的平面曲线类型的测地线有且仅有两条,分别为 r 1 = { x , 0 , x 2 2 a 2 } 和 r 2 = { 0 , y , y 2 2 b 2 } 。
证:椭圆抛物面 r = { a u cos v , b u sin v , 1 2 u 2 } 的一般方程为 x 2 a 2 + y 2 b 2 = 2 z ( a > 0 , b > 0 ) 。
令 F ( x , y , z ) = x 2 a 2 + y 2 b 2 − 2 z ,则 F x ( x , y , z ) = 2 x a 2 , F y ( x , y , z ) = 2 y b 2 , F z ( x , y , z ) = − 2 。
所以椭圆抛物面在点 ( x 0 , y 0 , z 0 ) 处的切平面方程为
x 0 a 2 ( x − x 0 ) + y 0 b 2 ( y − y 0 ) − ( z − z 0 ) = 0 .
法线方程为
x − x 0 x 0 a 2 = y − y 0 y 0 b 2 = z − z 0 − 1 .
则法线的一般方程为
{ x − x 0 x 0 a 2 = z − z 0 − 1 y − y 0 y 0 b 2 = z − z 0 − 1 ,
即
{ a 2 x + x 0 z − a 2 x 0 − x 0 z 0 = 0 b 2 y + y 0 z − b 2 y 0 − y 0 z 0 = 0 .
所以椭圆抛物面一点处的法截面束为
l ( a 2 x + x 0 z − a 2 x 0 − x 0 z 0 ) + m ( b 2 y + y 0 z − b 2 y 0 − y 0 z 0 ) = 0 ,
即
a 2 l x + b 2 m y + ( x 0 l + y 0 m ) z − a 2 x 0 l − x 0 z 0 l − b 2 y 0 m − y 0 z 0 m = 0 .
其中 l , m 是不全为零的任意实数。
因为
a 2 l ⋅ x 0 a 2 + b 2 m ⋅ y 0 b 2 − ( x 0 l + y 0 m ) = 0 ,
所以椭圆抛物面每点处的法截面与该点处的切平面垂直。
当 ( x 0 , y 0 , z 0 ) 为原点 ( 0 , 0 , 0 ) 点时,法截面束与椭圆抛物面相交得到的法截线为
{ a 2 l x + b 2 m y = 0 x 2 a 2 + y 2 b 2 = 2 z ,
存在三种情况:
① l = 0 , m ≠ 0 时法截面方程为 y = 0 ,法截线的参数方程为
r 1 = { x , 0 , x 2 2 a 2 } . (1)
② l ≠ 0 , m = 0 时法截面方程为 x = 0 ,法截线的参数方程为
r 2 = { 0 , y , y 2 2 b 2 } . (2)
③ l ≠ 0 , m ≠ 0 时法截面方程为 a 2 l x + b 2 m y = 0 ,法截线的参数方程为
r 3 = { x , − a 2 l x b 2 m , ( a 6 l 2 + b 6 m 2 ) x 2 2 a 2 b 6 m 2 } . (3)
情况①中法截线(1)上除原点外任意一点 ( x 1 , 0 , x 1 2 2 a 2 ) 处的切平面方程为
x 1 a 2 ( x − x 1 ) − ( z − x 1 2 2 a 2 ) = 0 .
因为
x 1 a 2 × 0 + 0 × 1 + ( − 1 ) × 0 = 0 ,
所以通过(1)的法截面 y = 0 与(1)上任意一点处的切平面都垂直。
从而由引理2知法截线(1)在其上任意点处的切平面上的正投影曲线为直线,是测地线。
情况②中法截线(2)上除原点外任意一点 { 0 , y 1 , y 1 2 2 b 2 } 处的切平面方程为
y 1 b 2 ( y − y 1 ) − ( z − y 1 2 2 b 2 ) = 0 .
因为
0 × 1 + y 1 b 2 × 0 + ( − 1 ) × 0 = 0 ,
所以通过(2)的法截面 x = 0 与(2)上任意一点处的切平面都垂直。
从而由引理2知法截线(2)在其上任意点处的切平面上的正投影曲线为直线,是测地线。
情况③中法截线(3)上除原点外任意一点 { x 1 , − a 2 l x 1 b 2 m , ( a 6 l 2 + b 6 m 2 ) x 1 2 2 a 2 b 6 m 2 } 处的切平面方程为
x 1 a 2 ( x − x 1 ) − a 2 l x 1 b 4 m ( y + a 2 l x 1 b 2 m ) − ( z − ( a 6 l 2 + b 6 m 2 ) x 1 2 2 a 2 b 6 m 2 ) = 0 .
因为
x 1 a 2 ⋅ a 2 l + ( − a 2 x 1 l b 4 m ) ⋅ b 2 m + ( − 1 ) × 0 = x 1 l − a 2 x 1 l b 2 ≠ 0 ,
所以通过(3)的法截面 a 2 l x + b 2 m y = 0 与(3)上除原点外任意一点处的切平面都不垂直,由引理2知法截线(3)在其上任意点处的切平面上的正投影曲线都不是直线,则法截线(3)不是测地线。
从而椭圆抛物面 r = { a u cos v , b u sin v , 1 2 u 2 } 上经过原点的非直线的平面曲线类型的测地线方程有且仅有两条,分别为 r 1 = { x , 0 , x 2 2 a 2 } 和 r 2 = { 0 , y , y 2 2 b 2 } 。
下面给出椭圆抛物面上测地线向原点处切平面作正投影的图像(见图3),其中红色线曲线是测地线,绿色曲线是测地线在原点处切平面上的正投影曲线。
定理2旋转抛物面 r = { a u cos v , a u sin v , 1 2 u 2 } 上原点处的法截面束与旋转抛物面相交得到的法截线都是测地线。
图3. 椭圆抛物面上测地线在原点处作正投影(a) (b)
证:旋转抛物面 r = { a u cos v , a u sin v , 1 2 u 2 } 的一般方程为 x 2 a 2 + y 2 a 2 = 2 z ( a > 0 , b > 0 ) 。
令 F ( x , y , z ) = x 2 a 2 + y 2 a 2 − 2 z ,则 F x ( x , y , z ) = 2 x a 2 , F y ( x , y , z ) = 2 y a 2 , F z ( x , y , z ) = − 2 。
所以旋转抛物面在点 ( x 0 , y 0 , z 0 ) 处的切平面方程为
x 0 a 2 ( x − x 0 ) + y 0 a 2 ( y − y 0 ) − ( z − z 0 ) = 0 .
法线方程为
x − x 0 x 0 a 2 = y − y 0 y 0 a 2 = z − z 0 − 1 .
则法线的一般方程为
{ x − x 0 x 0 a 2 = z − z 0 − 1 y − y 0 y 0 a 2 = z − z 0 − 1 ,
即
{ a 2 x + x 0 z − a 2 x 0 − x 0 z 0 = 0 a 2 y + y 0 z − a 2 y 0 − y 0 z 0 = 0 .
所以旋转抛物面一点处的法截面束为
l ( a 2 x + x 0 z − a 2 x 0 − x 0 z 0 ) + m ( a 2 y + y 0 z − a 2 y 0 − y 0 z 0 ) = 0 ,
即
a 2 l x + a 2 m y + ( x 0 l + y 0 m ) z − a 2 x 0 l − x 0 z 0 l − a 2 y 0 m − y 0 z 0 m = 0 .
其中 l , m 是不全为零的任意实数。
因为
a 2 l ⋅ x 0 a 2 + a 2 m ⋅ y 0 a 2 − ( x 0 l + y 0 m ) = 0 ,
所以旋转抛物面每点处的法截面与该点处的切平面垂直。
当 ( x 0 , y 0 , z 0 ) 为原点 ( 0 , 0 , 0 ) 点时,法截面束与旋转抛物面相交得到的法截线为
{ a 2 l x + a 2 m y = 0 x 2 a 2 + y 2 a 2 = 2 z ,
存在三种情况:
① l = 0 , m ≠ 0 时法截面方程为 y = 0 ,法截线的参数方程为
r 4 = { x , 0 , x 2 2 a 2 } . (4)
② l ≠ 0 , m = 0 时法截面方程为 x = 0 ,法截线的参数方程为
r 5 = { 0 , y , y 2 2 a 2 } . (5)
③ l ≠ 0 , m ≠ 0 时法截面方程为 a 2 l x + a 2 m y = 0 ,法截线的参数方程为
r 6 = { x , − l x m , ( l 2 + m 2 ) x 2 2 a 2 m 2 } . (6)
情况①中法截线(4)上除原点外任意一点 ( x 1 , 0 , x 1 2 2 a 2 ) 处的切平面方程为
x 1 a 2 ( x − x 1 ) − ( z − x 1 2 2 a 2 ) = 0 .
因为
x 1 a 2 × 0 + 0 × 1 + ( − 1 ) × 0 = 0 ,
所以通过(4)的法截面 y = 0 与(4)上任意一点处的切平面都垂直。
从而由引理2知法截线(4)在其上任意点处的切平面上的正投影曲线为直线,是测地线。
情况②中法截线(5)上除原点外任意一点 { 0 , y 1 , y 1 2 2 a 2 } 处的切平面方程为
y 1 a 2 ( y − y 1 ) − ( z − y 1 2 2 a 2 ) = 0 .
因为
0 × 1 + y 1 a 2 × 0 + ( − 1 ) × 0 = 0 ,
所以通过(5)的法截面 x = 0 与(5)上任意一点处的切平面都垂直。
从而由引理2知法截线(5)在其上任意点处的切平面上的正投影曲线为直线,是测地线。
情况③中法截线(6)上除原点外任意一点 { x 1 , − l x 1 m , ( l 2 + m 2 ) x 1 2 2 a 2 m 2 } 处的切平面方程为
x 1 a 2 ( x − x 1 ) − l x 1 a 2 m ( y + l x 1 m ) − ( z − ( l 2 + m 2 ) x 1 2 2 a 2 m 2 ) = 0
因为
x 1 a 2 ⋅ a 2 l + ( − x 1 l a 2 m ) ⋅ a 2 m + ( − 1 ) × 0 = 0 ,
所以通过(6)的法截面 a 2 l x + a 2 m y = 0 与(6)上任意一点处的切平面都垂直。
从而由引理2知法截线(6)在其上任意点处的切平面上的正投影曲线为直线,是测地线。
综上旋转抛物面 r = { a u cos v , a u sin v , 1 2 u 2 } 上原点处的法截面束与旋转抛物面相交得到的法截线都是测地线。
下面给出旋转抛物面上部分测地线向原点处切平面作正投影的图像(见图4),其中红色线曲线是测地线,绿色曲线是测地线在原点处切平面上的正投影曲线。
图4. 旋转抛物面上测地线在原点处作正投影
定理3双曲抛物面 r = { a ( u + v ) , b ( u − v ) , 2 u v } 上经过原点的非直线的平面曲线类型的测地线有且仅有两条,分别为 r 7 = { x , 0 , x 2 2 a 2 } 和 r 8 = { 0 , y , − y 2 2 b 2 } 。
证:双曲抛物面 r = { a ( u + v ) , b ( u − v ) , 2 u v } 的一般方程为 x 2 a 2 − y 2 b 2 = 2 z ( a > 0 , b > 0 ) 。
令 F ( x , y , z ) = x 2 a 2 − y 2 b 2 − 2 z ,则 F x ( x , y , z ) = 2 x a 2 , F y ( x , y , z ) = − 2 y b 2 , F z ( x , y , z ) = − 2 ,
所以双曲抛物面在点 ( x 0 , y 0 , z 0 ) 处的切平面方程为
x 0 a 2 ( x − x 0 ) − y 0 b 2 ( y − y 0 ) − ( z − z 0 ) = 0 .
法线方程为
x − x 0 x 0 a 2 = y − y 0 − y 0 b 2 = z − z 0 − 1 .
法线的一般方程为
{ x − x 0 x 0 a 2 = z − z 0 − 1 y − y 0 − y 0 b 2 = z − z 0 − 1 ,
即
{ a 2 x + x 0 z − a 2 x 0 − x 0 z 0 = 0 b 2 y − y 0 z − b 2 y 0 + y 0 z 0 = 0 .
则双曲抛物面一点处的法截面束为
l ( a 2 x + x 0 z − a 2 x 0 − x 0 z 0 ) + m ( b 2 y − y 0 z − b 2 y 0 + y 0 z 0 ) = 0 ,
即
a 2 l x + b 2 m y + ( x 0 l − y 0 m ) z − a 2 x 0 l − x 0 z 0 l − b 2 y 0 m + y 0 z 0 m = 0 .
其中 l , m 是不全为零的任意实数。
因为
a 2 l ⋅ x 0 a 2 − y 0 b 2 ⋅ b 2 m − ( x 0 l − y 0 m ) = 0 ,
所以双曲抛物面每点处的法截面与该点处的切平面垂直。
当 ( x 0 , y 0 , z 0 ) 为原点 ( 0 , 0 , 0 ) 点时,法截面束与双曲抛物面相交得到的法截线为
{ a 2 l x + b 2 m y = 0 x 2 a 2 − y 2 b 2 = 2 z ,
存在三种情况:
① l = 0 , m ≠ 0 时法截面方程为 y = 0 ,法截线的参数方程为
r 7 = { x , 0 , x 2 2 a 2 } . (7)
② l ≠ 0 , m = 0 时法截面方程为 x = 0 ,法截线的参数方程为
r 8 = { 0 , y , − y 2 2 b 2 } . (8)
③ l ≠ 0 , m ≠ 0 时法截面方程为 a 2 l x + b 2 m y = 0 ,法截线的参数方程为
r 9 = { x , − a 2 l x b 2 m , ( b 6 m 2 − a 6 l 2 ) x 2 2 a 2 b 6 m 2 } . (9)
情况①中法截线(7)上除原点外任意一点 { x 1 , 0 , x 1 2 2 a 2 } 处的切平面方程为
x 1 a 2 ( x − x 1 ) − ( z − x 1 2 2 a 2 ) = 0 .
因为
x 1 a 2 × 0 + 0 × 1 + ( − 1 ) × 0 = 0 ,
所以通过(7)的法截面 y = 0 与(7)上任意一点处的切平面都垂直。
从而由引理2知法截线(7)在其上任意点处的切平面上的正投影曲线为直线,是测地线。
情况②中法截线(8)上除原点外任意一点 { 0 , y 1 , − y 1 2 2 b 2 } 处的切平面方程为
− y 1 b 2 ( y − y 1 ) − ( z + y 1 2 2 b 2 ) = 0 .
因为
0 × 1 + ( − y 1 b 2 ) × 0 + ( − 1 ) × 0 = 0 ,
所以通过(8)的法截面 x = 0 与(8)上任意一点处的切平面都垂直。
从而由引理2知法截线(8)在其上任意点处的切平面上的正投影曲线为直线,是测地线。
情况③中法截线(9)上除原点外任意一点 { x 1 , − a 2 l x 1 b 2 m , ( b 6 m 2 − a 6 l 2 ) x 1 2 2 a 2 b 6 m 2 } 处的切平面方程为
x 1 a 2 ( x − x 1 ) + a 2 l x 1 b 4 m ( y + a 2 l x 1 b 2 m ) − ( z − ( b 6 m 2 − a 6 l 2 ) x 1 2 2 a 2 b 6 m 2 ) = 0 .
因为
x 1 a 2 ⋅ a 2 l + a 2 x 1 l b 4 m ⋅ b 2 m + ( − 1 ) × 0 = x 1 l + a 2 x 1 l b 2 ≠ 0 ,
所以通过(9)的法截面 a 2 l x + b 2 m y = 0 与(9)上任意一点处的切平面都不垂直,由引理2知法截线(9)在其上任意点处的切平面上的正投影曲线都不是直线,则法截线(9)不是测地线。
从而双曲抛物面 r = { a ( u + v ) , b ( u − v ) , 2 u v } 上经过原点的非直线的平面曲线类型的测地线方程有且仅有两条,分别为 r 7 = { x , 0 , x 2 2 a 2 } 和 r 8 = { 0 , y , − y 2 2 b 2 } 。
下面给出双曲抛物面上测地线向原点处切平面作正投影的图像(见图5),其中红色线曲线是测地线,绿色曲线是测地线在原点处切平面上的正投影曲线。
图5. 双曲抛物面上测地线在原点处作正投影(a) (b)
本文利用测地曲率的几何意义求出了抛物面上经过原点的平面曲线类型的测地线,为测地线的求法提供了新思路。目前本文的工作仅针对于经过曲面原点的测地线,将本文的方法应用到求经过曲面任意一点的测地线的可行性也是值得研究的。
作者非常感谢相关文献对本文的启发以及审稿专家提出的宝贵意见。
国家自然科学基金(11661062);
内蒙古自然科学基金(2022MS01015);
内蒙古自治区高等学校青年科技英才支持计划(NJYT-19-A09)。
吴钰莹,包图雅. 有关抛物面上测地线的研究Research on Geodesics on Parabolic Surfaces[J]. 理论数学, 2022, 12(11): 2033-2043. https://doi.org/10.12677/PM.2022.1211220