我试图了解如何使用中心差分方程计算曲率半径。因此,我使用圆半径上的三个点的简单示例。请在下图中找到示例。
我使用的方程式如下图所示。
我的代码如下:
z_1=0
z_2=50
z_3=0
x_1=0
x_2=50
x_3=100
Eq_1 = (1 + (dz_dx) ** 2) ** (3 / 2)
dz_dx = (z_3 - z_1) / (2 * (x_3 - x_1))
Eq_2 = (z_3 - 2 * z_2 + z_1) / (x_3 - x_1) ** 2
r = Eq_1 / Eq_2
我的问题是计算的半径是 100。但它应该是 50。
你看到我的计算有什么错误吗?
我会欢迎每一个帮助。
提前致谢!