arclength 函数计算一般空间曲线的长度。 如果曲线被假定为分段线性（单线性），这并不太难，但如果曲线是平滑插值的样条，则问题变得更加困难。
在三次情况下解决这个问题的一个好方法是将其表述为数值积分，因此可以从 quadgk 获得解决方案。 仔细编写的代码将问题分解为曲线上每对断点之间的一系列积分，以避免出现否则会发现的导数奇点。
2021-11-03 15:21:03
3KB
1
我必须计算地球子午线的弧长。 我注意到对此没有解析解，所以我写了这个非常简单的函数。
它计算以 (0,0) 为中心、半径为 a（沿 OX）和半径为 b（沿 OY）的椭圆的弧长
x(t) = a.cos(t) y(t) = b.sin(t)
t1 和 t2 之间的角度 t（以弧度为单位）。
通过将圆弧划分为多个小直线段，可以得到数值上的解。 另外，第二个输出是用近似公式（Ramanujan）补全的完整椭圆的长度。
吕克
2021-10-10 20:44:12
2KB
1
大地测量学基础日常作业，包含二维、三维坐标转换，大地主题正反算，大地坐标系与空间直角坐标系的转化，贝塞尔大地问题正反算，子午线弧长正反算，兰勃投影，高斯投影正反算，附加Eign矩阵运算
2019-12-21 20:54:55
2.82MB
1
三次样条插值，拉格朗日插值，牛顿插值，B样条插值，cardinal样条插值。使用Romberg积分求弧长，使用二分法在曲线上找特定弧长的点
2019-12-21 19:52:59
32.28MB
1