5 Результаты расчета
Результаты расчета занесены в результирующий файл, который имеет следующий вид:
Rechenie differencialnogo yravnenia metodom Runge-Kutta
Rachet vupolnil:Ivan Byivid gr.101071
Isxodnie dannie
m = 1.000 b = 0.200
c = 5.000 F = 10.000
Shag scheta h = 0.028
Shag pechati hp = 0.700
Vrema scheta tmax = 14.000
---------------------------------------------------
Tablica reshenia differencialnogo yravnenia
---------------------------------------------------
t | y | y` | y " |
---------------------------------------------------
0.00 0.000 0.000 0.000
0.70 1.770 4.183 2.512
1.43 5.013 4.177 -2.030
2.13 7.245 1.987 -3.855
2.83 7.668 -0.757 -3.694
3.53 6.343 -2.833 -2.026
4.23 4.065 -3.370 0.593
4.90 2.136 -2.085 3.109
5.60 1.552 0.486 3.674
6.30 2.640 2.351 1.406
7.00 4.417 2.444 -0.997
7.70 5.755 1.233 -2.241
8.40 6.047 -0.396 -2.216
9.10 5.295 -1.632 -1.179
9.80 3.987 -1.918 0.399
10.53 2.837 -1.072 1.793
11.23 2.567 0.318 1.925
11.93 3.187 1.321 0.809
12.63 4.196 1.404 -0.523
13.33 4.973 0.724 -1.294
14.03 5.147 -0.227 -1.298
---------------------------------------------------
По данным таблицы строим график результата решения (рис. 3).
Рисунок 3 – График решения дифференциального уравнения
6 Анализ результатов
Статическое значение y =
F/c = 10/5 = 2.
Это же значение, полученное решением
уравнения на компьютере, также равно
2. Период собственных колебаний равен
.
Из графика получаем Тм = 5,5.
Относительная ошибка равна
.
Отсюда следует, уравнение решено верно.
7 Литература
1. Молибошко л.А. Компьютерное моделирование автомобилей. – Мн.: ивц Минфина рб, 2007. – 280 с.