Ctrl →

Olmer · Апрель 21, 2014 09:37:39

Привет.
Есть задача у меня, решить уравнение движения в обычном симметрично потенциале гравитационном. Я её решил методом Эйлера и всё хорошо. Но потом, пытаясь решить с помощью integrate.odeint эту систему, я очень сильно застрял.
Как это можно реализовать? Гугление дало некоторое представление, но не помогло. Сначала думаю, что всё ок и начало всё работать, но потом результат выходит такой, что сразу видно, как всё плохо. Помогите!

FishHook · Апрель 21, 2014 09:45:14

Olmer
Помогите!

В чем помощь то должна заключаться?

Olmer · Апрель 21, 2014 09:55:14

FishHook

Ой, да, точно. Есть Центр притяжения с массой M=2*10**33, постоянная тяготения G=6.67*10**-8, от этого центра на расстоянии x=10**13 мы запускаем пробную точку с начальной скоростью по оси y=3*10**6, наша задача построить кривую.
Ну и решить надо с помощью, integrate.odeint.

Отредактировано Olmer (Апрель 21, 2014 10:00:04)

fongostev · Апрель 21, 2014 10:08:08

Эх вы, решаем в Лагранжевом формализме аналитически, получаем интегралы движения и строим траектории численно, делов-то.

Olmer · Апрель 21, 2014 10:10:21

fongostev
Эх вы, решаем в Лагранжевом формализме аналитически, получаем интегралы движения и строим траектории численно, делов-то.

А если нам нужно построить 100 000 тысяч таких решений? А если потенциал не симметричный? Я спрашиваю не для того, чтобы мне задачу именно эту решить, а для того, чтобы понять, как вообще всё это работает.

fongostev · Апрель 21, 2014 10:13:40

Вот это видел? Все вроде понятно написано
http://stackoverflow.com/a/16008005

Olmer · Апрель 21, 2014 10:17:33

fongostev
Вот это видел? Все вроде понятно написаноhttp://stackoverflow.com/a/16008005

Видел. Но там ОДУ 1 порядка, а нам надо решить второго порядка.

Отредактировано Olmer (Апрель 21, 2014 10:20:02)

fongostev · Апрель 21, 2014 10:22:15

Olmer

Понизь степень уравнения, перейди к экв. системе и интегрируй.

Отредактировано fongostev (Апрель 21, 2014 10:24:23)

Olmer · Апрель 21, 2014 10:27:05

fongostev

Понизив, будем вынуждены перейти к системе из двух уравнений, и тут я тоже туплю. В гугле были примеры, когда решали уравнения колебаний и строили фигуру Лисажу(там ур-ие второго порядка), но попытки написать по той аналогии свою задачу у меня провалились.

Отредактировано Olmer (Апрель 21, 2014 10:34:19)

fongostev · Апрель 21, 2014 10:33:33

Olmer

Ты делал как написано тут?
http://nbviewer.ipython.org/gist/anonymous/d417c1ffbb76f13f678c

UPD: У тебя же уравнения Лагранжа сводятся к такому виду только так.

Отредактировано fongostev (Апрель 21, 2014 10:34:41)

Python-сообщество

Уведомления

#1 Апрель 21, 2014 09:37:39

Решение ОДУ с scipy.inegrate.odeint.

#2 Апрель 21, 2014 09:45:14

Решение ОДУ с scipy.inegrate.odeint.

#3 Апрель 21, 2014 09:55:14

Решение ОДУ с scipy.inegrate.odeint.

#4 Апрель 21, 2014 10:08:08

Решение ОДУ с scipy.inegrate.odeint.

#5 Апрель 21, 2014 10:10:21

Решение ОДУ с scipy.inegrate.odeint.

#6 Апрель 21, 2014 10:13:40

Решение ОДУ с scipy.inegrate.odeint.

#7 Апрель 21, 2014 10:17:33

Решение ОДУ с scipy.inegrate.odeint.

#8 Апрель 21, 2014 10:22:15

Решение ОДУ с scipy.inegrate.odeint.

#9 Апрель 21, 2014 10:27:05

Решение ОДУ с scipy.inegrate.odeint.

#10 Апрель 21, 2014 10:33:33

Решение ОДУ с scipy.inegrate.odeint.

Board footer