python scipy integrate_使用scipy.integrate.odeint（）计算衰减

我正试图用 scipy.integrate.odeint() 来计算一些物种的衰变 . 其中一个衰减系数非常大，这似乎导致了一个问题，我不知道为什么或如何解决它 . 这是我的代码：

from scipy.integrate import odeint

import numpy as np

import matplotlib.pyplot as plt

def decay(init,t):

a,b,c,d = init

f0 = - 6*a

f1 = 3 * a - 1e12*b

f2 = 3 * a

f3 = 1e12 * b

return [f0,f1,f2,f3]

if __name__ == '__main__':

init = [1,0,0,0]

time = np.linspace(0, 1e12, 100001)

soln = odeint(decay, init ,time)

这给了我以下错误：

lsoda-- at t (=r1) and step size h (=r2), the

corrector convergence failed repeatedly

or with abs(h) = hmin ls

in above, r1 = 0.0000000000000D+00 r2 = 0.1552206225945D-09

/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/integrate/odepack.p

python scipy integrate_使用scipy.integrate.odeint（）计算衰减相关推荐

python寻峰算法_python/scipy的寻峰算法
我可以自己写一些东西,通过找到一阶导数的零交叉点或其他东西,但它似乎是一个足够通用的函数,可以包含在标准库中.有人知道吗? 我的特殊应用是一个二维数组,但通常它会用于在FFT等中查找峰值. 具体地说, ...
为python安装numpy和scipy(federo)
为了进行数值计算,例如积分等等,需要安装numpy和scipy,其中scipy是依赖于numpy的,所以先要装numpy. 1, 通过下载http://pypi.python.org/pypi/num ...
统计学（检验、分布）的 python（numpy/pandas/scipy）实现
scipy 中统计相关的 api:https://docs.scipy.org/doc/scipy/reference/stats.html https://zhuanlan.zhihu.com/p/ ...
（Python）使用Gdal+Scipy获得Dem的经纬度的高程值（双线性和三次样条内插）
(Python)使用Gdal+Scipy获得Dem的经纬度的高程值(双线性和三次样条内插) 前言基本原理双线性内插三次样条内插代码实现后记前言最近使用python进行一些的遥感影像处理, ...
linux 版本的scipy,linux安装scipy
1.1.Numba的约5分钟指南 Numba是Python的即时编译器,它最适用于使用NumPy数组和函数以及循环的代码.使用Numba的最常用方法是通过其装饰器集合,可以应用于您的函数来指示Numb ...
SciPy稀疏矩阵类 scipy.sparse（持续更新ing...）
诸神缄默不语-个人CSDN博文目录 scipy官方文档:SciPy documentation - SciPy v1.9.0 Manual scipy.sparse官方文档:Sparse matric ...
【SciPy库】scipy.optimize.fmin_l_bfgs_b进行L-BFGS优化
[时间]2020.01.07 [题目][SciPy库]scipy.optimize.fmin_l_bfgs_b进行L-BFGS优化具体用法参考官方文档:scipy.optimize.fmin_l_b ...
scipy中的scipy.optimize.curve_fit
scipy中的scipy.optimize.curve_fit 这里写目录标题 scipy中的scipy.optimize.curve_fit 参数 Return scipy.optimize.``c ...
Python使用numpy包编写自定义函数计算平均绝对误差(MAE、Mean Absolute Error)、评估回归模型和时间序列模型、解读MAE
Python使用numpy包编写自定义函数计算平均绝对误差(MAE.Mean Absolute Error).评估回归模型和时间序列模型.解读MAE 目录

python scipy integrate_使用scipy.integrate.odeint（）计算衰减

python scipy integrate_使用scipy.integrate.odeint（）计算衰减相关推荐

最新文章

热门文章