Давайте теперь действительно решим нашу задачу. Допустим, что мы взяли ε=0,1сек (если после того, как мы проделаем все вычисления, окажется, что этот интервал не достаточно мал, то необходимо повторить всё сначала с меньшим интервалом времени, например 0,01сек). Чему будет равно z(0,1), если в начальный момент времени z(0)=1?

z(t + ε)=z(t)+v_z(t)ε

=> z(0,1)= 1+0*0,1=1

скорость в начальный момент времени равна нулю. Таким образом, z(0,1)=1, ибо грузик ещё не начал двигаться. Но новая скорость:

v_z(t+ε)= v_z(t)+a_zε= v_z(t)−z(t)ε

=> v_z (0,1)=0 −1*0,1= −0,1

В момент 0,20 сек

z(t + ε)=z(t)+v_z(t)ε

=> z(0,2)= z(0,1)+v_z(0,1)*0,1=1−0,1*0,1=0,99;

И

v_z(t+ε)= v_z(t)+a_zε= v_z(t)−z(t)ε

=> v_z (0,2)= v_z(0,1)−z(0,1)*0,1 = −0,1 −1*0,1 = −0,2

Продолжая эту процедуру еще и еще, можно найти положение и скорость в любой момент времени, а это как раз то, что нам нужно. Однако практически мы используем нехитрый прием, который позволит увеличить точность вычислении. Если бы мы продолжали начатые нами расчеты, то они оказались бы довольно грубыми, поскольку интервал ε=0,1сек довольно большой. Пришлось бы уменьшить его, скажем, до 0,01сек. Но тогда, чтобы проследить движение за какой-то разумный отрезок времени, потребовалось бы сделать множество шагов. Мы же организуем процесс таким образом, что сможем увеличить точность, используя тот же интервал ε=0,1сек. Этого можно достичь, несколько изменив метод расчета.

Заметьте, что новое положение тела равно старому плюс интервал времени ε, умноженный на скорость. Но что это за скорость? В какой момент? В начале интервала одна скорость, а в конце она совсем другая. Прием состоит в том, чтобы брать скорость в середине интервала. Если известна скорость в настоящий момент и известно, что она меняется, как же можно надеяться получить удовлетворительный результат, считая, что тело все время движется с той же скоростью, что и в настоящий момент? Более разумно использовать какую-то среднюю скорость между началом и концом интервала. Те же рассуждения применимы к изменению самой скорости: для подсчета ее изменений нужно использовать ускорение в средней точке между двумя моментами времени, в которых необходимо найти скорость. Таким образом, реально мы будем пользоваться следующими уравнениями: положение в конце интервала равно положению в начале плюс интервал ε, умноженный на скорость в середине интервала. Эта скорость в свою очередь равна скорости в середине предыдущего интервала (т. е. на отрезок ε меньше) плюс ускорение в начале интервала, умноженное на ε.

z(t + ε)=z(t)+ ε v_z(t+ε/2).

v_z(t+ε/2)= v_z(t − ε/2)+a_z(t)ε= v_z(t − ε/2)−z(t)ε

a_z(t)=−z(t)

Таблица 9.1

РЕШЕНИЕ УРАВНЕНИЯ (dv_z/dt) =−z

Интервал ε = 0,10 сек

t	z	vz	аz
0,0	1,000	0,000	-1,000
		—0,050
0,1	0,995		-0,995
		—0,150
0,2	0,980		—0,980
		—0,248
0,3	0,955		-0,955
		-0,343
0,4	0,921		-0,921
		' —0,435
0,5	0,877	0523	—0,877
0,6	0,825	'	—0,825
		—0,605
0,7	0,765		—0,765
		—0,682
0,8	0,696		-0,693
		—0,752
0,9	0,621		—0,621
		—0,814
1,0	0,540		—0,540
		. __ п ЯР8
1,1	0,453	"~"1^,ООи	—0,453
		-0,913
1,2	0,362		-0,362
		—0,949
1,3	0,267		-0,267
		—0,976
1,4	0,169		-0,169
		—0,993
1,5	0,069		—0,070
		-1,000
1 (У	Г» АОГ\

Остается еще один небольшой вопрос: что такое v_z(ε/2)? Вначале у нас было v_z(0), а не v_z( — ε /2). Но теперь, чтобы начать наши вычисления, необходимо использовать дополнительное уравнение

v_z(ε/2)= v_z(0)+a_z(0)ε/2= v_z(0)−z(0)ε/2

Ну, а теперь все готово для расчетов. Для удобства можно их выполнить в виде таблицы, в столбцах которой стоят время, положение, скорость и ускорение, причем скорость пишется в промежутках между строками (табл. 9.1). Такая таблица есть, конечно, просто удобный способ записи результатов, полученных из уравнений, и фактически полностью заменяет их. Мы просто заполняем одно за другим свободные места в ней и получаем очень интересную картину движения: сначала грузик находится в покое, затем понемногу приобретает отрицательную скорость (вверх), а это приводит к уменьшению его расстояния от точки равновесия. При этом хотя ускорение и становится меньше, оно все еще «подгоняет» скорость. Однако по мере приближения к положению равновесия (z=0) ускорение становится все меньше и меньше, скорость нарастает все медленней и медленней, но все же еще нарастает вплоть до точки z=0, которая достигается примерно через 1,5сек. Скажем по секрету, что произойдет дальше. Грузик, конечно, не остановится в точке z=0, а пойдет дальше, но теперь все будет наоборот: его положение z станет отрицательным, а ускорение—положительным. Скорость начнет уменьшаться. Интересно сравнить полученные нами числа с функцией соs(t). Результат этого сравнения представлен на фиг. 9.4. Оказывается, что в пределах точности наших расчетов (три знака после запятой) совпадение полное! Позднее вы узнаете, что функция соs(t) — точное решение нашего уравнения, так что у вас теперь есть наглядное представление о мощи численного анализа: столь простой расчет дает столь точный результат.