Редактирование: Моделирование движения системы двойной звезды
Внимание! Вы не авторизовались на сайте. Ваш IP-адрес будет публично видимым, если вы будете вносить любые правки. Если вы войдёте или создадите учётную запись, правки вместо этого будут связаны с вашим именем пользователя, а также у вас появятся другие преимущества.
Правка может быть отменена. Пожалуйста, просмотрите сравнение версий, чтобы убедиться, что это именно те изменения, которые вас интересуют, и нажмите «Записать страницу», чтобы изменения вступили в силу.
| Текущая версия | Ваш текст | ||
| Строка 24: | Строка 24: | ||
При написании программы я использовала метод численного интегрирования Верле для нахождения следующей координаты траектории, зная текущее и предыдущее значение: | При написании программы я использовала метод численного интегрирования Верле для нахождения следующей координаты траектории, зная текущее и предыдущее значение: | ||
| − | x(t+Δt)=2x(t)-x(t-Δt)+ax(t)Δt^2 | + | [x(t+Δt)=2x(t)-x(t-Δt)+ax(t)Δt^2 |
| − | y(t+Δt)=2y(t)-y(t-Δt)+ay(t)Δt^2 | + | y(t+Δt)=2y(t)-y(t-Δt)+ay(t)Δt^2 ] |
Ускорение считается из формулы второго закона Ньютона, а сила из закона гравитации: | Ускорение считается из формулы второго закона Ньютона, а сила из закона гравитации: | ||
| + | {{ | ||
F(t)=G*m1*m2/r^2; | F(t)=G*m1*m2/r^2; | ||
a(t)=F(t)/m; | a(t)=F(t)/m; | ||
| − | + | }} | |
Программа сохраняет координаты движения звезд в файл, графики строила в программе gnuplot. | Программа сохраняет координаты движения звезд в файл, графики строила в программе gnuplot. | ||
