python求解一阶常微分方程
一.用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求解一阶常微分方程相关推荐
- 一阶欧拉近似matlab,MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程.doc
MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程 姓名:樊元君 学号:2012200902 日期:2012.11.06 一.实验目的 掌握改进欧拉法与四阶龙格-库塔求解一阶常微分方程的初值问题 ...
- python求解一阶线性偏微分方程通解举例
python求解一阶线性偏微分方程的通解举例 Python求解偏微分方程也是其一个应用方面,下面举例说明. 一.问题: 求一阶线性偏微分方程 x ∂ f ( x , y ) ∂ x − y ∂ f ( ...
- MATLAB和Python求解非线性常微分方程
龙格-库塔方法 M阶Runge-Kutta方法用于求解方程, dydx=f(x,y),y(a)=y0(1)\frac{\boldsymbol{dy}}{\boldsymbol{dx}}=\boldsy ...
- 使用C++,用四阶Runge-Kutta的方法来求解一阶常微分方程
#include <iostream> using namespace std;/* dy/dx = y - 2x/y, 0< x <= 1 步长h = 0.2 */ cons ...
- python求解四阶微分方程_用Python求解二阶常微分方程组的RungeKutta四阶解
本文试图用龙格库塔四阶法数值求解两个常微分方程组. 初始系统: 要解决的系统: 我有非常奇怪的解图... 我有: 正确的图形: 我在我的龙格库塔找不到麻烦.请帮帮我.在 我的代码在这里:dt = 0. ...
- 经典龙格-库塔法(四阶龙格-库塔法)求解求一阶常微分方程相应的特解的Python程序
基本原理 例题 代码 #四阶龙格-库塔法 #求一阶常微分方程,相应的特解 #x变量的区间 a = 0 b = 1 #已知条件 X = [0] Y = [1] h = 0.2 #设置步长 n = (b- ...
- python解常微分方程龙格库_求解二阶常微分方程的RungeKutta四阶方法
我试着做一个简谐振子的例子,它将用龙格-库塔四阶法求解.要求解的二阶常微分方程(ODE)和初始条件为: y''+y=0 y(0)=0和y'(0)=1/pi 范围在0到1之间,共有100步.我用u作为辅 ...
- 数值计算一阶常微分方程求解实现
原创作品,出自 "晓风残月xj" 博客,欢迎转载,转载时请务必注明出处(http://blog.csdn.net/xiaofengcanyuexj). 由于各种原因,可能存在诸多不 ...
- Python求解常微分方程——sympy
[常微分方程简介] 方程中未知量是函数而不是变量,且未知量涉及未知函数的导数的方程称为微分方程. 常微分方程(ordinary differential equation, ODE)是一类特殊情况,未 ...
最新文章
- Linux软件包命令
- 阿里云高级工程师认证机会!
- 启动jar包并生成日志的linux脚本
- CTFshow 反序列化 web275
- 在sqlserver 中with(nolock)详解
- 浅谈对Ubuntu桌面系统的实验性理解
- TDSQL 全时态数据库系统--核心技术
- 论文浅尝 | Convolutional 2D knowledge graph embedding
- bool查询原理 es_es6.2.4,使用bool查询查出的结果,SearchResponse的Hits[]总是比t..._慕课猿问...
- android安全权限管理,Android 11 中的权限更新
- 使用nginx反向代理获取百度MP3的真实网址
- BZOJ 2594: [Wc2006]水管局长数据加强版( LCT )
- 51单片机初值计算方法
- 地图制作:Google Earth Pro的下载及功能介绍(详细介绍)(上)
- 传言阿里P10赵海平,被P11多隆打3.25后离职,如何评价赵海平对王垠的面试?
- 柴静十年看见了什么——一个央视记者的心灵史
- 暴风酷播云二期配置_暴风酷播云 一期-N3160版: 硬件折解及安装Proxmox VE-服务器虚拟化系统...
- uniapp 之 禁用手机物理返回键
- Python获取下周一日期
- 赛前采访里皮表示:希望队员用真心踢球。
热门文章
- HTML5 CSS控制Table内外边框、颜色、大小示例
- LCD不带显存,是如何内存映射屏幕。S5PV210SoC在内存中选一段内存存放颜色数据,通过配置将LCD控制器和这一段内存连接起来,构成映射关系,LCD控制器就自动从显存中读取像素数据传给LCD驱动器
- 顶尖C++程序员5分钟编程:逆战凉了,绝地求生黄了,唯有他站住脚!
- 针对刀具磨损的日志读取的曲线分析(一维高斯滤波及波形拟合)
- Linux中TTY是什么意思
- 斯人若彩虹,遇上方知有——dbGet(一)
- php能抓抖音短视频教程,初学者必须要掌握的抖音短视频拍摄的小技巧
- FPGA 结构分析 -IO 资源
- 华清远见上海中心培训感言
- vue使用iview中Upload上传组件