记scipy中odeint函数用法 |
您所在的位置:网站首页 › obedient怎么记 › 记scipy中odeint函数用法 |
记scipy中odeint函数用法 小白艰难自学之路-(一) odeint介绍odeint()函数是scipy库中一个数值求解微分方程的函数 odeint()函数需要至少三个变量,第一个是微分方程函数,第二个是微分方程初值,第三个是微分的自变量。 一个一阶微分方程例子 import numpy as np import matplotlib.pyplot as plt from scipy.integrate import odeint def diff(y, x): return np.array(x) # 上面定义的函数在odeint里面体现的就是dy/dx = x x = np.linspace(0, 10, 100) # 给出x范围 y = odeint(diff, 0, x) # 设初值为0 此时y为一个数组,元素为不同x对应的y值 # 也可以直接y = odeint(lambda y, x: x, 0, x) plt.plot(x, y[:, 0]) # y数组(矩阵)的第一列,(因为维度相同,plt.plot(x, y)效果相同) plt.grid() plt.show()odeint()函数中第一个变量微分方程的函数中可以定义不止一个一阶微分方程,定义多个一阶微分方程就可以解高阶方程,下面是一个解 |
今日新闻 |
推荐新闻 |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |