方程的数值解.pdf

方程的数值解 现在我们来真正解前面的问题。假定取ε 0.100 秒。当我们做好这一切工 作后，如果发现这还不够小，我们可以再回过头来，以ε 0.010 秒重做一次。 从初值x 0 1.00 开始，x 0.1 是多少呢？它是原来的位置x 0 加上速度 ( ) ( ) ( ) (此时为 )乘 秒。于是x 0.1 仍是 ，因为它还没有开始运动。但在 0 0.10 ( ) 1.00 秒时的新速度就是原速度v 0 0 加 乘以加速度。加速度是 0.10 ( ) ε −x 0 =−1.00 。于是， ( ) v 0.1 0.01=−0.10 ×1.00 =−0.10 ( ) 现在，在0.20 秒时 x 0.20 x 0.1 =+εv 0.1 ( ) ( ) ( ) 1.00=−0.10×0.10 0.99 和 v 0.2 v 0.1 =+εa 0.1 ( ) ( ) ( ) =−0.10 −0.10 ×1.00 =−0.20 依此类推，一直做下去，就可计算出其余的运动，这正是我们要做的事。然而， 实际上，这里有点小小的技巧可用来提高准确度。假如我们继续已经开始的计算， 就会发现由于ε 0.100 秒是相当粗糙的，因而运动也是相当粗糙的，我们不 得不取一个更小的时间间隔，比如说ε 0.010 秒。于是，要对一段适当的总 时间间隔进行研究，就要作大量的重复计算。所以我们将在用同样粗糙的间隔 ε 0.100 秒的条件下，把要计算的工作组织一下以提高准确度。这一点在分 析技巧上略加改进就可以办到。 我们注意到，新的位置是老的位置加上时间间隔 乘以速度。但这是什么时ε 刻的速度呢？在时间间隔开始时是一个速度，在时间间隔结束时又是另外一个速 度。我们的改进就是利用二者之间的速度。同样的考虑也可用于速度：为了计算 速度变化，我们将使用要求出它的速度的那两个时刻中间的加速度。这样我们实 第 1 页，共 7 页 完整版，请访问完整版，请访问科大科院考研网，专注于中科大、中科院考研科大科院考研网，专注于中科大、中科院考研 际使用的方程就有如下述：后来的位置等于先前的位置加上ε乘以在间隔中间的 那个时刻的速度。类似地，间隔中间那个时刻的速度等于比它早ε时(它正处在 t 前一个时间间隔中间) 的速度加上ε乘以 时刻的加速度，也就是说，我们利用的 方程是 x t +ε x t =+εv t +ε 2 ( ) ( ) ( ) v t +ε 2 v t =−ε 2 +εa t (1) ( ) ( ) () a t =−x t