我正在使用scipy.integrate.solve_ivp在Python语言中运行一个长的ODE集成。是否可以在例程运行时访问集成的状态或检查该例程处于哪个集成步骤？我的集成花费的时间比预期的要长，我想知道是积分器卡在某个步骤上了，还是个别步骤花了很长时间。对于未来的任务；如果我将与solve_ivp的集成拆分为子间隔，以便在其间打印状态消息，这是否会影响某些求
浏览 22提问于2020-10-04得票数 0
回答已采纳
所以我在想:某些函数的积分之和不是应该等于函数之和的积分吗？在这里，我借助scipy的集成了三个任意函数from scipy.integrate import solve_ivp
def fun(t, y): return (0.5 * y) + (0.05 *
浏览 0提问于2020-09-16得票数 0
我正在尝试使用SciPy来解决一组ODE。分配给我的任务要求我在500个时间步长内求解微分方程。如何使用SciPy实现这一点？到目前为止，我已经尝试使用scipy.integrate.solve_ivp，它给了我一个正确的解决方案，但我无法控制它运行的时间步数。t_span参数允许我配置t的初始值和终值，但我实际上对此并不感兴趣--我只对我集成了多少次感兴趣。(例如，当我使
浏览 9提问于2021-05-13得票数 1
我尝试使用solve_ivp解决一个非常简单的问题import numpy as np
Traceback (most recent call last):
File "<ipython-input= 'LSODA'，尽
浏览 0提问于2019-02-21得票数 1
我正在使用Scipy.Integrate的solve_ivp方法来解决ivp，我希望能够在我为积分提供的时间步长上计算函数，但我不知道如何做到这一点。我可以回顾积分中的每个元素，但是除了已经花费在求解ivp上的时间之外，这将花费大量的时间，所以我更希望能够在实际方法计算积分过
浏览 18提问于2019-06-18得票数 0
我使用scipy.integrate.solve_ivp来解决一个ODE系统，因为它有事件函数。我需要这个函数的原因是，在积分过程中，有时我会得到一个奇异矩阵，每次发生这种情况时，我都需要完成积分并使用新参数重新启动它。我想知道在终端事件发生后，是否可以使用新参数重新启动scipy.integrate.solve_ivp，如果可能，我如何做到这一
浏览 20提问于2019-06-23得票数 1
2回答
我有一个函数，我想得到它的积分函数，如下所示：也就是说，我不需要在x点获得一个单一的集成值，而是需要在multiple points上获得值。在我的实际代码中，f(x)非常复杂，而且它是一个多重集成，因此运行时间太慢()。点Phi(20)处的积分值是来自Phi(19)的
浏览 1提问于2016-02-06得票数 6
我正在将我的代码从使用scipy的odeint转换为使用scipy的solve_ivp。当使用odeint时，我会使用一个while循环，如下所示： solver.integrate(t_final, step=True)# do other operations
这个方法允许我在每个时间步之后存储依赖于解决方案<e
浏览 3提问于2019-10-18得票数 1
1回答
在python中，我想为任意函数y(x)计算该函数在x的两个值之间的绝对积分。我知道我可以用以下函数计算出“正规”积分：如何计算函数的绝对积分？用‘绝对值’，我的意思是，y为负的区域也是，将加
浏览 1提问于2020-01-04得票数 0
回答已采纳
问题是，我希望能够立即对网格的每个点开始的微分方程进行积分，而不是在每个坐标的scipy积分器上循环。我认为代码的其余部分并不重要，因为它是与scipy集成的一部分，也是我使用meshgrid函数numpy的一部分。import solve_ivp
L = 1total_run_time = 10*3
grid_poin
浏览 0提问于2019-05-18得票数 0
回答已采纳
1回答
编辑t_list = np.linspace(0, 30, 100)[1. 1.端编辑
我正在尝试重新创建一个简单的matlab系统的微分方程到python使用scipy。我不知道为什么在执行过程中我会收到一个RuntimeWarning: invalid value encountered in double_scalars。我在odeint调用
浏览 10提问于2021-09-19得票数 1
回答已采纳
1回答
抛射运动的奇异值
、、
在给定初始条件的情况下，我绘制了抛射运动的曲线图。from scipy.integrate import odeint
import m
浏览 24提问于2019-10-28得票数 0
我用solve_ivp数值解微分方程，以便从初始条件积分位置和速度。这在Python中是直接的，而且很容易做到。但是，我也在使用numdifftools.Jacobian计算最终位置和速度的部分，与初始位置/速度和时间有关。这个我更难掌握了。让我解释一下。我有一组微分方程(参见rhs_CR3BP函数)，我使用Python solve_ivp在timespan [0,T]上进行数值积
浏览 1提问于2022-01-07得票数 1
回答已采纳
4回答
在枕中积分多维积分
、、、、
动机：i有一个多维积分，为了完整起见，我在下面复制了这个积分。当存在明显的各向异性时，它来自于第二维数系数的计算：这里W是所有变量的函数。它是一个已知的函数，我可以为它定义一个python函数。
编程问题：如何让scipy集成这个表达式？我想把两个三重四角()链接在一起，但我担心性能和准确性。在scipy中是否有一个
浏览 7提问于2012-12-28得票数 21
回答已采纳
1回答
我用的是二重积分，内积分有可变界。我写了一个函数，使用SciPy的四边形积分，允许我计算这个积分。然而，我想要的只是计算内部积分，所以我剩下的只是一个关于某个变量的，未计算的积分。然后我想画出这个“中途”计算的二重积分与该
浏览 14提问于2020-03-14得票数 0
回答已采纳
回答这一问题的过程已经在下面的链接问题中开始，但这一专题具体是关于整合一个功能，并得到了答复。所以我增加了一个新问题。
import matplotlib.pyplot as plt
ax2.plot(s.t[::2], s.y[1][::2], label="
浏览 3提问于2020-11-26得票数 0
回答已采纳
我是matplotlib/python的新手。我想用odeint或solve_ivp在某个区间内解决一些简单的常微分方程，例如segment0，10，其中“初始条件”是t_0= 5，而不是t_0 = 0。在我看来，odeint和solve_ivp都将初始条件作为值y_0，也就是说，y_0是函数在指定时间数组t开始时的值。
在我的积分间隔中间有一个条件，我如何</e
浏览 2提问于2020-05-17得票数 1
我期待着计算一个三维积分，在不规则网格上带有预定义值的函数。< coords = self._grid.weights / d我真的不知道如何获得combined的积分，假设coords
浏览 15提问于2022-04-04得票数 0
1回答
我试图使用python/scipy计算两个极限之间的积分。import scipy as sp return np.exp(-0.5*x**2)
def int_testand upper
浏览 2提问于2016-10-12得票数 0
回答已采纳
1回答
下面是对数正态分布的概率密度函数： from scipy.stats import lognorm
def f(x): return lognorm.pdf(x, s=0.2, loc=0, scale=np.exp(10)) 此函数的y值非常小(max ~ 1E-5)，并且分布在x值~1E5上。我们知道PDF的积分应该是1，但当使用以下代码直接计算积分时，答案是1E-66，因为计算</em
浏览 29提问于2020-09-08得票数 0
回答已采纳
云服务器
ICP备案
云直播
腾讯会议
实时音视频
腾讯云开发者
