Мне дается следующий код в 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
@BlueRineS Да, и я знаю, как его использовать, как и много раз раньше, но на этот раз Ива получила неизвестную ошибку, и, кстати, я хотел бы точно перевести эти строки из Mathematica
в Python
У вас есть неизвестная ошибка? Можете ли вы поделиться им прямо в своем посте? включая всю трассировку стека и включая весь собственный код того, что вы пробовали до сих пор?
в течение 1 минуты я обновлю, подождите плз
Если подпись callable(t, y, ...)
, то должен быть установлен аргумент tfirst
True
?
@Джо, я не понял тебя
Посмотрите документ на odeint
, последнее ключевое слово.
@Joe Моя проблема в том, что когда я использую mathematica, я получаю ответ, но используя ODINT, я получаю сообщение об ошибке
Трудно сказать, что здесь не так. Для отладки начните с простого уравнения, которое работает, возможно, используйте пример из документации odeint и постепенно усложняйте. Тогда скажите нам, где это терпит неудачу.
@joe ODINT работает хорошо, на этот раз странно. Я обнаружил, что это противоречит (1-od)**
в уравнении. Mmathematica решает это легко, но ODEINT делает это невозможным. Я должен знать, почему.
Я получаю ValueError: negative number cannot be raised to a fractional power
, что более логично. То, что вы получаете, является последующей ошибкой внезапного изменения типа с реального на сложный. Какие версии python и scipy вы используете?
Кроме того, можете ли вы указать, что такое зависимые и независимые переменные? Возможно, вам нужно поставить tfirst=True
, чтобы вообще получить правильную интеграцию. При этом интеграция проходит без ошибок.
@LutzL Python 3.6.3, последняя версия Scipy и Numpy.
Затем были изменения, которые сделали его слишком умным для их же блага, автоматически превращая результаты рациональной силы в сложные всякий раз, когда реальных результатов не существовало. Внутри odeint
передает значения производных в код Fortran, работающий с действительными числами, требуя преобразования комплексного числа в действительное, что невозможно без числовой ошибки.
Из вызова 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, и я запускал его для каждого скопированного-вставленного шага, все было правильно, пока весь код не был скопирован и вставлен. Теперь я не получаю ошибку, пока я вообще ничего не меняю
это мой третий раз, когда в одной оболочке Python я получил ошибку, но в новой все правильно.
Итак, у вас есть что-то в старой рабочей области, которая меняет работу интерпретатора Python. Вам придется проанализировать всю историю команд этого сеанса, чтобы выяснить, что это может быть. Если в только что запущенной оболочке все хорошо, то проблема решена? Если проблема возникает в более крупном сценарии, максимально уменьшите его, сохранив при этом ошибку, либо вы найдете источник ошибки, либо зададите новый хороший вопрос.
Да, tfirst=True
сделал это правильно. Теперь это то, что делает NDSolve
математики. Я никогда раньше не использовал tfirst=True.
Это относительно новая опция, возможно, с тем же капитальным ремонтом, что и интерфейс solve_ivp
.
Вы читали документацию odeint?