本人python小白，以前都在matlab里做流程模拟之类的，现在越来越多的人建议我试试python （numpy，scipy，matplotlib）。
现在学习过程中有个问题，在python里如何用4阶龙格库塔法解微分方程组呢？（在matlab里是写一个函数文件，然后ode45调用并赋初值）
大佬可否根据实例在python写一段程序我参考学习一下 （最后还需要绘制曲线t-y）
2条回答 默认 最新
CSDN专家-黄老师 2021-05-25 19:28关注参考一下:https://www.pythonf.cn/read/103447,如有帮助，望采纳，谢谢
本回答被题主选为最佳回答 , 对您是否有帮助呢?解决 无用评论 打赏举报
分享
- 2022-06-09 16:52回答 1 已采纳 看这样行不行 from scipy.integrate import odeint import numpy as np import matplotlib.pyplot as plt def d
- 2021-06-16 12:07回答 2 已采纳 这个错误是使用了未定义的变量.
- 2020-11-29 20:52回答 1 已采纳 from scipy.integrate import odeint 这个可以 但是方程太多了不知道行不行 之前解过一个三个等式的可以用
- 2023-03-14 08:30微小冷的博客 `scipy`文档中将`odeint`函数和`ode, comples_ode`这两个类称为旧API，是scipy早期使用的微分方程求解器，但由于是Fortran实现的，速度没得说，所以有的时候还挺推荐使用的。
- 2021-07-16 16:21回答 1 已采纳 参考一下：https://blog.csdn.net/your_answer/article/details/79234275https://blog.csdn.net/linshen1213/art
- 2022-09-26 16:59回答 4 已采纳 可以用sympy 求解 from sympy import * x = Symbol('x') y = Symbol('y') z = Symbol('z') print(solve([5x+3y+(
- 2022-04-27 00:34回答 3 已采纳 Numerov方法是数值求解常微分方程： 均匀格点与对数格点上的Numerov方法 | Shigaro minyez's blog on l
- 2024-06-04 08:30作家小姐的博客 用Matlab写了这一点是因为这学期选了一个科学计数法的网课选修课，它里边讲Matlab，所以对它有一点接触，但不多，刚好提到可以用这个写作业，有点好奇，所以就在网上找了破解版下载了，目前只会求解一阶微分方程和...
- 2021-11-02 11:55回答 2 已采纳 math中的sin()函数有错误，用sympy的sin()函数；另外，solve(fx,x)你的x的位置写到solve外面去了；另外，我改好了以后发现程序也是执行不出来，拆解成两部分运行后，发现后半部
- 2022-05-09 19:14回答 1 已采纳 import turtle as tt tt.pensize(2) d = -45 for i in range(4): tt.seth(d) d += 90 tt.fd(2
- 2022-09-28 18:30回答 2 已采纳 你的函数是你自己定义的：f(x) = b * a^x + c拟合就是就是求a,b,c参数，popt[0]是参数a，popt[1]是参数b，popt[2]是参数c
- 2022-02-02 21:15伊滴小朋友的博客 本篇博客主要讲解了二阶常系数微分方程的解析解和数值解的解法以及一阶常系数微分方程组的通解解法，并对每一个步骤做了详细的说明（网上很多帖子讲的并不是很详细），并做了一些探究性的尝试
- 2020-12-24 14:04小野校长的博客 引入首先用高数知识求解非齐次常系数微分方程image.png再利用信号与系统中冲激响应求解验证image.png利用MATLAB求解验证y=dsolve('D2y+3*Dy+2*y=exp(-t)','y(0)=1','Dy(0)=2','t')得出结果：y =(t - 2 exp(-t) + 3) ...
- 2019-08-10 04:16在Python中，我们通常需要手动计算微分方程的解，然后使用matplotlib或其他可视化库来绘制。但通过plotdf，用户可以直接输入微分方程的定义，模块会自动进行数值求解和绘图。 首先，让我们了解一下如何使用plotdf。...
- 2020-12-20 18:48接下来，第二个例子展示了如何求解一个一阶非线性微分方程，并绘制其解的图形： \[ f'(x) + f(x) + f(x)^2 = 0 \] 同样，我们使用`dsolve()`函数求解： ```python from sympy import * f = symbols('f', cls=...
- 没有解决我的问题, 去提问
悬赏问题
- ¥15 宇视监控服务器无法登录
- ¥15 PADS Logic 原理图
- ¥15 PADS Logic 图标
- ¥15 电脑和power bi环境都是英文如何将日期层次结构转换成英文
- ¥20 气象站点数据求取中~
- ¥15 如何获取APP内弹出的网址链接
- ¥15 wifi 图标不见了 不知道怎么办 上不了网 变成小地球了
- ¥50 STM32单片机传感器读取错误
- ¥50 power BI 从Mysql服务器导入数据，但连接进去后显示表无数据
- ¥15 (关键词-阻抗匹配，HFSS，RFID标签天线)