Решение набора экспоненциальных уравнений в питоне

Пытаюсь решить набор экспоненциальных уравнений, но продолжаю получать ошибки, говорящие о том, что «невозможно преобразовать число с плавающей запятой» или «число с плавающей запятой не может быть вызвано».

Вы можете запустить код, просто замените переменные любыми значениями, которые вы хотите.

from sympy.abc import x, y
import numpy as np
import sympy as sp
import math 

A0 = 58
G0 = 44
Gmax = 117
tmax = 40

s=[((sp.log(x) - sp.log(y)) - (tmax*(x-y))), (Gmax- G0 - (A0 * x /(y-x))*((sp.exp((-x)*tmax)) - (sp.exp((-y)*tmax))))]
sp.solve(s, x, y)

python numpy math exponential equations
person Salchem    schedule 30.10.2017    source источник
comment
Возможный дубликат ошибка с вычислением интеграла python sympy для функции косинуса   -  person pacholik    schedule 30.10.2017
comment
@Salchem ​​Пожалуйста, предоставьте минимальный рабочий пример (MWE), чтобы опубликованный вами код мог работать в любой системе. Например, я не могу запустить all_training_df и solve_poly_system. Я не знаю, откуда берутся эти предметы.   -  person ThomasGuenet    schedule 30.10.2017
comment
Хорошо, см. выше, я заменил их   -  person Salchem    schedule 30.10.2017
comment
Я получаю следующую ошибку: не удалось решить -73*(y + LambertW(-40*yexp(-40*y))/40)*exp(40*y)*exp(-LambertW(-40 *yexp(-40*y))) - 29*(exp(40*y) - exp(-LambertW(-40*yexp(-40*y))))*LambertW (-40*yexp(-40*y))/20   -  person Salchem    schedule 30.10.2017

Ответы (1)


arrow_upward
1
arrow_downward

удалите математику и используйте sympy для журнала и опыта

from sympy.abc import x, y
import numpy as np
import sympy as sp
Gmax = np.amax(all_training_df.iloc[1])
tmax = np.argmax(all_training_df.iloc[1])
A0 = all_combined_df.iloc[1]['CHO (g)']
G0 = all_combined_df.iloc[1]['mg/dL']
s=[((sp.log(x) - sp.log(y)) - (tmax*(x-y))), (Gmax- G0 - (A0 * x /(y-x))*((sp.exp((-x)*tmax)) - (sp.exp((-y)*tmax))))]
sp.solve_poly_system(s, x, y)
person Smart Manoj    schedule 30.10.2017
comment
Да, пробовал, но не смог решить -73*(y + LambertW .......... См. выше - person Salchem; 30.10.2017