一.用python求解一阶常微分方程

1.求解微分方程需要用到scipy库,pycharm中安装即可,同时需要导入numpy库和matplotlib两个库

2.使用scipy.integrate.odeint()来进行求解,一般指使用三个参数,默认其他参数

func : 导数函数f(y,t),即y在t处的导数,用函数形式表示

y0 :初始条件y0,用数组的形式表示

t : 求解函数值对应的时间点的序列。序列的第一个元素是与初始条件 y0 对应的初始时间 t0;时间序列必须是单调递增或单调递减的,允许重复值

scipy.integrate.odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, ml=None, mu=None, rtol=None, atol=None, tcrit=None, h0=0.0, hmax=0.0, hmin=0.0, ixpr=0, mxstep=0, mxhnil=0, mxordn=12, mxords=5, printmessg=0, tfirst=False)
3.举例

对于该函数而言

#导入库函数
from scipy.integrate import odeint
import numpy as np
import matplotlib.pyplot as plt
#定义函数的导数
def dy_dt(y,t):return np.sin(t**2)
#给定初始值和时间范围
y0=[1]
t = np.arange(-10,10,0.01)
#使用odeint()方法:
y=odeint(dy_dt,y0,t)
#绘图
plt.plot(t, y)
plt.title("picture")
plt.show()

python求解一阶常微分方程相关推荐

  1. 一阶欧拉近似matlab,MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程.doc

    MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程 姓名:樊元君 学号:2012200902 日期:2012.11.06 一.实验目的 掌握改进欧拉法与四阶龙格-库塔求解一阶常微分方程的初值问题 ...

  2. python求解一阶线性偏微分方程通解举例

    python求解一阶线性偏微分方程的通解举例 Python求解偏微分方程也是其一个应用方面,下面举例说明. 一.问题: 求一阶线性偏微分方程 x ∂ f ( x , y ) ∂ x − y ∂ f ( ...

  3. MATLAB和Python求解非线性常微分方程

    龙格-库塔方法 M阶Runge-Kutta方法用于求解方程, dydx=f(x,y),y(a)=y0(1)\frac{\boldsymbol{dy}}{\boldsymbol{dx}}=\boldsy ...

  4. 使用C++,用四阶Runge-Kutta的方法来求解一阶常微分方程

    #include <iostream> using namespace std;/* dy/dx = y - 2x/y, 0< x <= 1 步长h = 0.2 */ cons ...

  5. python求解四阶微分方程_用Python求解二阶常微分方程组的RungeKutta四阶解

    本文试图用龙格库塔四阶法数值求解两个常微分方程组. 初始系统: 要解决的系统: 我有非常奇怪的解图... 我有: 正确的图形: 我在我的龙格库塔找不到麻烦.请帮帮我.在 我的代码在这里:dt = 0. ...

  6. 经典龙格-库塔法(四阶龙格-库塔法)求解求一阶常微分方程相应的特解的Python程序

    基本原理 例题 代码 #四阶龙格-库塔法 #求一阶常微分方程,相应的特解 #x变量的区间 a = 0 b = 1 #已知条件 X = [0] Y = [1] h = 0.2 #设置步长 n = (b- ...

  7. python解常微分方程龙格库_求解二阶常微分方程的RungeKutta四阶方法

    我试着做一个简谐振子的例子,它将用龙格-库塔四阶法求解.要求解的二阶常微分方程(ODE)和初始条件为: y''+y=0 y(0)=0和y'(0)=1/pi 范围在0到1之间,共有100步.我用u作为辅 ...

  8. 数值计算一阶常微分方程求解实现

    原创作品,出自 "晓风残月xj" 博客,欢迎转载,转载时请务必注明出处(http://blog.csdn.net/xiaofengcanyuexj). 由于各种原因,可能存在诸多不 ...

  9. Python求解常微分方程——sympy

    [常微分方程简介] 方程中未知量是函数而不是变量,且未知量涉及未知函数的导数的方程称为微分方程. 常微分方程(ordinary differential equation, ODE)是一类特殊情况,未 ...

最新文章

  1. Linux软件包命令
  2. 阿里云高级工程师认证机会!
  3. 启动jar包并生成日志的linux脚本
  4. CTFshow 反序列化 web275
  5. 在sqlserver 中with(nolock)详解
  6. 浅谈对Ubuntu桌面系统的实验性理解
  7. TDSQL 全时态数据库系统--核心技术
  8. 论文浅尝 | Convolutional 2D knowledge graph embedding
  9. bool查询原理 es_es6.2.4,使用bool查询查出的结果,SearchResponse的Hits[]总是比t..._慕课猿问...
  10. android安全权限管理,Android 11 中的权限更新
  11. 使用nginx反向代理获取百度MP3的真实网址
  12. BZOJ 2594: [Wc2006]水管局长数据加强版( LCT )
  13. 51单片机初值计算方法
  14. 地图制作:Google Earth Pro的下载及功能介绍(详细介绍)(上)
  15. 传言阿里P10赵海平,被P11多隆打3.25后离职,如何评价赵海平对王垠的面试?
  16. 柴静十年看见了什么——一个央视记者的心灵史
  17. 暴风酷播云二期配置_暴风酷播云 一期-N3160版: 硬件折解及安装Proxmox VE-服务器虚拟化系统...
  18. uniapp 之 禁用手机物理返回键
  19. Python获取下周一日期
  20. 赛前采访里皮表示:希望队员用真心踢球。

热门文章

  1. HTML5 CSS控制Table内外边框、颜色、大小示例
  2. LCD不带显存,是如何内存映射屏幕。S5PV210SoC在内存中选一段内存存放颜色数据,通过配置将LCD控制器和这一段内存连接起来,构成映射关系,LCD控制器就自动从显存中读取像素数据传给LCD驱动器
  3. 顶尖C++程序员5分钟编程:逆战凉了,绝地求生黄了,唯有他站住脚!
  4. 针对刀具磨损的日志读取的曲线分析(一维高斯滤波及波形拟合)
  5. Linux中TTY是什么意思
  6. 斯人若彩虹,遇上方知有——dbGet(一)
  7. php能抓抖音短视频教程,初学者必须要掌握的抖音短视频拍摄的小技巧
  8. FPGA 结构分析 -IO 资源
  9. 华清远见上海中心培训感言
  10. vue使用iview中Upload上传组件