-dc/dt = kc(t) - k'n(t)
Я вычисляю c (t) и n(t) из траектории моделирования с помощью CuPy. Затем я вычисляю k и k' с помощью Curve_fit с помощью следующего кода:
Код: Выделить всё
# compute dc/dt
dc_dt = np.gradient(c_t, times)
# define the model function
def model(t, k, k_prime):
return -k * np.array(c_t) + k_prime * np.array(n_t)
initial_guess = [0.5, 0.5]
bounds = ([0.0, 0.0], [np.inf, np.inf])
# fit the model to the data
params, covariance = curve_fit(model, times, dc_dt)
k1_optimized, k2_optimized = params
Я попробовал сохранить значения c(t) и n( t) результаты в лист Excel и вычисление k1 и k2 на другом компьютере; в этом случае k2 не было 1,0. Мне было интересно, что происходит не так?
Подробнее здесь: https://stackoverflow.com/questions/793 ... -computers