Как изменить математический код ODE на Python

Мне дается следующий код в Mathematica

solution = NDSolve[{\[CapitalOmega]\[CapitalLambda]'[x] == 
(1 - \[CapitalOmega]\[CapitalLambda][x]) \[CapitalOmega]\[CapitalLambda][x] ((1 - \[CapitalOmega]\[CapitalLambda][x])^((1 - \[Delta])/(4 -\2 \[Delta]))  \[CapitalOmega]\[CapitalLambda][x]^(1/(4 - 2 \[Delta]))), \[CapitalOmega]\[CapitalLambda][0] == \[CapitalOmega]\[CapitalLambda]f}, \[CapitalOmega]\\[CapitalLambda], {x, xi, xf}];

Я хочу написать это на Python, используя ODEINT но я действительно не понимаю, как это написать, потому что он отправляет `error1 комплексного числа.

def OD_H(z, od, delt):
    dMdt = od * (1 - od) * ((1 - od)**((1 - delt)/(4-2*delt)) * od**(1/(2 *(2-delt))))
    return dMdt                                                                                                                      

def ant(z, od0, delt):
    z1 = 0
    od = odeint(OD_H, od0, [z1, z], args=(delt,))[-1]       
    return od     

for z in np.arange(0,3.1,0.1):
    print(ant(z, 0.7, 1.1))

ошибка

od = odeint(OD_H, od0, [z1, z], args=(delt,))[-1] File "C:\Python36-32\lib\site-packages\scipy\integrate\odepack.py", line 244, in odeint int(bool(tfirst))) TypeError: can't convert complex to float

Вы читали документацию odeint?

WiseDev 04.07.2019 12:49

@BlueRineS Да, и я знаю, как его использовать, как и много раз раньше, но на этот раз Ива получила неизвестную ошибку, и, кстати, я хотел бы точно перевести эти строки из Mathematica в Python

Ma Y 04.07.2019 12:51

У вас есть неизвестная ошибка? Можете ли вы поделиться им прямо в своем посте? включая всю трассировку стека и включая весь собственный код того, что вы пробовали до сих пор?

WiseDev 04.07.2019 12:52

в течение 1 минуты я обновлю, подождите плз

Ma Y 04.07.2019 12:52

Если подпись callable(t, y, ...), то должен быть установлен аргумент tfirstTrue ?

Joe 04.07.2019 15:10

@Джо, я не понял тебя

Ma Y 04.07.2019 16:35

Посмотрите документ на odeint, последнее ключевое слово.

Joe 04.07.2019 19:08

@Joe Моя проблема в том, что когда я использую mathematica, я получаю ответ, но используя ODINT, я получаю сообщение об ошибке

Ma Y 04.07.2019 21:59

Трудно сказать, что здесь не так. Для отладки начните с простого уравнения, которое работает, возможно, используйте пример из документации odeint и постепенно усложняйте. Тогда скажите нам, где это терпит неудачу.

Joe 04.07.2019 22:04

@joe ODINT работает хорошо, на этот раз странно. Я обнаружил, что это противоречит (1-od)** в уравнении. Mmathematica решает это легко, но ODEINT делает это невозможным. Я должен знать, почему.

Ma Y 04.07.2019 23:16

Я получаю ValueError: negative number cannot be raised to a fractional power, что более логично. То, что вы получаете, является последующей ошибкой внезапного изменения типа с реального на сложный. Какие версии python и scipy вы используете?

Lutz Lehmann 04.07.2019 23:32

Кроме того, можете ли вы указать, что такое зависимые и независимые переменные? Возможно, вам нужно поставить tfirst=True, чтобы вообще получить правильную интеграцию. При этом интеграция проходит без ошибок.

Lutz Lehmann 04.07.2019 23:34

@LutzL Python 3.6.3, последняя версия Scipy и Numpy.

Ma Y 04.07.2019 23:54

Затем были изменения, которые сделали его слишком умным для их же блага, автоматически превращая результаты рациональной силы в сложные всякий раз, когда реальных результатов не существовало. Внутри odeint передает значения производных в код Fortran, работающий с действительными числами, требуя преобразования комплексного числа в действительное, что невозможно без числовой ошибки.

Lutz Lehmann 04.07.2019 23:58
Почему в Python есть оператор "pass"?
Почему в Python есть оператор "pass"?
Оператор pass в Python - это простая концепция, которую могут быстро освоить даже новички без опыта программирования.
Некоторые методы, о которых вы не знали, что они существуют в Python
Некоторые методы, о которых вы не знали, что они существуют в Python
Python - самый известный и самый простой в изучении язык в наши дни. Имея широкий спектр применения в области машинного обучения, Data Science,...
Основы Python Часть I
Основы Python Часть I
Вы когда-нибудь задумывались, почему в программах на Python вы видите приведенный ниже код?
LeetCode - 1579. Удаление максимального числа ребер для сохранения полной проходимости графа
LeetCode - 1579. Удаление максимального числа ребер для сохранения полной проходимости графа
Алиса и Боб имеют неориентированный граф из n узлов и трех типов ребер:
Оптимизация кода с помощью тернарного оператора Python
Оптимизация кода с помощью тернарного оператора Python
И последнее, что мы хотели бы показать вам, прежде чем двигаться дальше, это
Советы по эффективной веб-разработке с помощью Python
Советы по эффективной веб-разработке с помощью Python
Как веб-разработчик, Python может стать мощным инструментом для создания эффективных и масштабируемых веб-приложений.
0
14
200
1
Перейти к ответу Данный вопрос помечен как решенный

Ответы 1

Ответ принят как подходящий

Из вызова odeint мы узнаем, что вы используете z как переменную времени и od как переменную пространства/состояния. Поскольку ваша производная функция сначала имеет время, а по умолчанию в odeint сначала должно быть состояние, вам нужно установить опцию tfirst=True. С этим изменением ошибок нет. Вы получаете те же значения с меньшими вычислениями, используя всю мощь odeint, предоставляя ему список всех желаемых точек выборки.

z = np.arange(0,3.1,0.1)
od = odeint(OD_H, od0, z, args=(delt,), tfirst=True)

где тогда пары z[k], od[k] соответствуют парам аргумент-значение, которые вы вычисляете в цикле.

Это решено, вы знаете, я сделал это, но получил ту же ошибку. Я скопировал и вставил код шаг за шагом в новый файл Python, и я запускал его для каждого скопированного-вставленного шага, все было правильно, пока весь код не был скопирован и вставлен. Теперь я не получаю ошибку, пока я вообще ничего не меняю

Ma Y 05.07.2019 12:04

это мой третий раз, когда в одной оболочке Python я получил ошибку, но в новой все правильно.

Ma Y 05.07.2019 12:06

Итак, у вас есть что-то в старой рабочей области, которая меняет работу интерпретатора Python. Вам придется проанализировать всю историю команд этого сеанса, чтобы выяснить, что это может быть. Если в только что запущенной оболочке все хорошо, то проблема решена? Если проблема возникает в более крупном сценарии, максимально уменьшите его, сохранив при этом ошибку, либо вы найдете источник ошибки, либо зададите новый хороший вопрос.

Lutz Lehmann 05.07.2019 12:55

Да, tfirst=True сделал это правильно. Теперь это то, что делает NDSolve математики. Я никогда раньше не использовал tfirst=True.

Ma Y 05.07.2019 13:31

Это относительно новая опция, возможно, с тем же капитальным ремонтом, что и интерфейс solve_ivp.

Lutz Lehmann 05.07.2019 14:12

Другие вопросы по теме