常微分方程解法（2）例题

用SciPyのscipy.integrate.odeint解联立微分方程式，以时间t（tの范围在0~几秒就可以，例：t=0~2秒）和解x(t), y(t), z(t)作图。

作为t=0的初期条件，x(t=0) = -10, y (t=0) = 0, z (t=0) = 35.0。还有，关于系数a,b,c，可以尝试a = 40, b = 5, c = 35 的情况和a = 40, b = 10, c = 35的情况。还有时间分割Δt要适当取小一点的值。另外，Δt如果取的过小，计算量会就变得很大。根据实际实行后的错误进行调整吧。

下面是解决的python程序：

import numpy as np
import scipy.integrate as sciin
import matplotlib.pyplot as plt#把与自变量t相关的导函数的函数名字放在F里
def f(F, t,params):    x,y,z = F       f_values = [a*(y-x),(c-a)*x-x*z+c*y,x*y-b*z]   #分别写x,y,z的导数等于的右边的式子 return f_values# 放系数
a = 40
b= 5
c = 35
#把上面三个系数放到parameters里
parameters = [a,b,c]#设定x,y,z的初始值
x0 = -10
y0 = 0.0
z0 = 35
#　把初始值放到Y0里
Y0 = [x0,y0,z0]# 开始点，结束点，间隔的设置
tStart = 0.0
tStop  = 2
tInc   = 0.01       #间隔
#把上面归纳到t里
t = np.arange(tStart, tStop, tInc)# sciin.odeint解ODE
solution = sciin.odeint(f, Y0, t, args=(parameters,))# 作图
plt.figure(figsize=(9.5, 6.5))
plt.plot(t, solution[:, 0], color='black')
plt.plot(t, solution[:, 1], color='green')
plt.plot(t, solution[:, 2], color='red')
plt.xlabel('time, t' , fontsize=14)
plt.ylabel('theta(t)', fontsize=14)
plt.show()

我简直是天才。

常微分方程解法（2）例题相关推荐

三元一次方程组例题_最新《三元一次方程组及其解法》例题与讲解
最新<三元一次方程组及其解法>例题与讲解 <三元一次方程组及其解法>例题与讲解 1.三元一次方程及三元一次方程组 (1)三元一次方程: 含有三个未知数,并且含未知数的项的次数都 ...
常微分方程及其基本例题
文章目录零基本概念微分方程向量场一可分离变量 1.1 变量分离 1.2 变量变换 1.3 齐次微分方程二一阶微分方程 2.1 齐次线性 2.2 非齐次线性三伯努利微分方程 3.1 ...
常微分方程解法：Euler‘s Method 欧拉方法
常微分方程的解法 (一): 常微分方程的离散化 :差商近似导数.数值积分方法.Taylor 多项式近似 Euler's Method 欧拉方法常微分方程的解法 (二): 欧拉(Euler)方法
最大全1子矩阵的两种解法（例题：City Game HDU - 1505）
以前牛客多校遇到过两道,都没做出来,这次来系统性的补习一下. 例题:City Game HDU - 1505 题意:给你一个矩阵,求最大全1子矩阵,最后结果乘以3... 全1矩阵可以参考下图这个框就 ...
协方差矩阵的解法和例题
看了一下百度出来的例题,全是晦涩难懂的公式,例题里面高赞票数普遍都是错误的一传十,十传百,于是写在这里供需要的人查阅,自己也可以加深印象.
lesson3 一阶线性常微分方程解法
微积分(三)-- 一阶常微分方程解法
一.积分公式二.积分因子三.温度-浓度模型(别名:传导-扩散模型) 四.利用积分因子求解积分方程五.实例
一阶线性常微分方程解法总结 Summary of Solving First Order Linear Ordinary Differential Equation (ODE)
文章目录 1. Separable Differentiable Equations 2. Linear Equations First-order linear equation Solution ...
matlab求解微分代数方程组,微分代数方程(DAE)的Matlab 解法.PDF
微分代数方程(DAE)的Matlab 解法微分代数方程(DAE)的Matlab解法所谓微分代数方程,是指在微分方程中,某些变量满足某些代数方程的约束.假设微分方程的更一般形式可以写成前面所介绍 ...

常微分方程解法（2）例题

常微分方程解法（2）例题相关推荐

最新文章

热门文章