优质解答
如果旋转面的参数方程是
r(u,v) = (u*cos(v), u*sin(v), f(u))
则其测地线公式是(一般微分几何书上都有)
v(u) = v0 + \\int_{u0}^{u} c * \\sqrt(1+[f'(w)]^2) / [w * \\sqrt(w^2-c^2)] dw
其中 (u0, v0)是出发点位置, \\int_{u0}^{u}表示从u0到u的积分,
\\sqrt是开根号, f'(w)表示 f 的导数, c是一个常数,由
v(u1) = v1 决定
由Liouville公式,还可以推出
dv/ds = c/u^2
其中 s 是测地线弧长参数.这样
ds/dv = u^2/c
ds = (u^2/c)dv = (u^2/c) dv/du *du
= (u^2/c) * c * \\sqrt(1+[f'(u)]^2) / [u * \\sqrt(u^2-c^2)] du
= u * \\sqrt(1+[f'(u)]^2) / \\sqrt(u^2-c^2) du
所以弧长等于 \\int_{u0}^{u1} ds
= \\int_{u0}^{u1} u * \\sqrt(1+[f'(u)]^2) / \\sqrt(u^2-c^2) du
对于 x^2+y^2-z^2=1, 令f(u) = 正负\\sqrt(u^2-1) 即可.这样:
x = u*cos(v), y = u*sin(v), z = 正负 \\sqrt(u^2-1)
x^2+y^2-z^2 = u^2-(u^2-1) = 1
上面的积分不一定能够有积出来的表达式, 实际计算的时候, 可以数值计算.
以u0=1, v0=0, u1=2, v1=0为例, 这对应着空间两点:
p1 = ( 1, 0, 0 ), p2 = ( 2, 0, \\sqrt(3) )
显然最短路径就是沿着母线走, 通过这两点的母线方程是
x^2-z^2=1, 也就是 z = \\sqrt(x^2-1) (这里只取上半部分)
由求曲线弧长的微积分公式, 最短路径长度是
\\int_{1}^{2} \\sqrt(1+[dz/dx]^2) dx
= \\int_{1}^{2} \\sqrt[(2x^2-1)/(x^2-1)] dx
这个积分没有封闭表达式结果,只能数值计算.
如果旋转面的参数方程是
r(u,v) = (u*cos(v), u*sin(v), f(u))
则其测地线公式是(一般微分几何书上都有)
v(u) = v0 + \\int_{u0}^{u} c * \\sqrt(1+[f'(w)]^2) / [w * \\sqrt(w^2-c^2)] dw
其中 (u0, v0)是出发点位置, \\int_{u0}^{u}表示从u0到u的积分,
\\sqrt是开根号, f'(w)表示 f 的导数, c是一个常数,由
v(u1) = v1 决定
由Liouville公式,还可以推出
dv/ds = c/u^2
其中 s 是测地线弧长参数.这样
ds/dv = u^2/c
ds = (u^2/c)dv = (u^2/c) dv/du *du
= (u^2/c) * c * \\sqrt(1+[f'(u)]^2) / [u * \\sqrt(u^2-c^2)] du
= u * \\sqrt(1+[f'(u)]^2) / \\sqrt(u^2-c^2) du
所以弧长等于 \\int_{u0}^{u1} ds
= \\int_{u0}^{u1} u * \\sqrt(1+[f'(u)]^2) / \\sqrt(u^2-c^2) du
对于 x^2+y^2-z^2=1, 令f(u) = 正负\\sqrt(u^2-1) 即可.这样:
x = u*cos(v), y = u*sin(v), z = 正负 \\sqrt(u^2-1)
x^2+y^2-z^2 = u^2-(u^2-1) = 1
上面的积分不一定能够有积出来的表达式, 实际计算的时候, 可以数值计算.
以u0=1, v0=0, u1=2, v1=0为例, 这对应着空间两点:
p1 = ( 1, 0, 0 ), p2 = ( 2, 0, \\sqrt(3) )
显然最短路径就是沿着母线走, 通过这两点的母线方程是
x^2-z^2=1, 也就是 z = \\sqrt(x^2-1) (这里只取上半部分)
由求曲线弧长的微积分公式, 最短路径长度是
\\int_{1}^{2} \\sqrt(1+[dz/dx]^2) dx
= \\int_{1}^{2} \\sqrt[(2x^2-1)/(x^2-1)] dx
这个积分没有封闭表达式结果,只能数值计算.