← Сtrl

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

fongostev

Да! Я про эту ссылку говорил. Вот мой вариант кода

from matplotlib import pyplot as plt
import numpy as np
from scipy.integrate import odeint
G=6.67*10**-8
M=2*10**33.
def motion(x_vec, t):
    x, y = x_vec  # tuple unpacking
    
    return [y, G*M/y**2.]
y0 = [10**13., 10**13.]
t_output = np.arange(0, 10**7, 10**3)
y_result = odeint(motion, y0, t_output)
plt.plot(y_result[:,0],y_result[:,1])
plt.show()

Я не понимаю, что происходит в теле функции, перед return.

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

Olmer

Там идет распаковка набора аргументов. Грубо говоря, там происходит следующее:

a, b, c = (a, b, c)

UPD: Т.е. принимается левая часть уравнений, а возвращается правая.

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

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

fongostev
x, y = x_vec

Вот здесь? “x, y = x_vec ” Немного не догоняю…

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

Olmer

Да, именно. Просто в интерпретаторе введи

t = (1,2)
a, b = t

И посмотри, что получится.

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

fongostev

Но x_vec, что такое? Что это за вектор?

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

Olmer

это твой набор независимых переменных, в твоем случае х и х с точкой.

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

fongostev

А ты можешь записать, как будет выглядеть тело функции, или у меня правильно написано?

fongostev · Апрель 21, 2014 11:02:16

Ты просто неаккуратно расписал систему. Замени строку и попробуй. Должно быть так:

return [y,  G*M/x**2]

Знаки только проверь

Отредактировано fongostev (Апрель 21, 2014 11:03:08)

Olmer · Апрель 21, 2014 11:08:26

fongostev
Ты просто неаккуратно расписал систему. Замени строку и попробуй. Должно быть так:

Он рисует, что-то похожее. Но если поменять знак на минут перед G*M, то получается страшное решение.

fongostev · Апрель 21, 2014 11:09:34

Ясное дело, ты же меняешь вектор действия потенциала на противоположный

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

Уведомления

#1 Апрель 21, 2014 10:37:20

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

#2 Апрель 21, 2014 10:38:42

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

#3 Апрель 21, 2014 10:47:11

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

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

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

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

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

#6 Апрель 21, 2014 10:54:25

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

#7 Апрель 21, 2014 10:57:22

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

#8 Апрель 21, 2014 11:02:16

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

#9 Апрель 21, 2014 11:08:26

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

#10 Апрель 21, 2014 11:09:34

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

Board footer