Python：求解自治二阶 ODE
我有以下自治常微分方程系统，其中x和v是向量（初始条件x0和 >v0) dx/dt = v dv/dt = F(x, v) 等价地，可以将其编写为二阶 ODE 系统 d^2 x/ dt^2 = F(…
有没有办法对 NumPy 数组中具有相同值的所有行应用函数？
假设我们有一个矩阵 A，它具有以下值： In [2]: A Out[2]: array([[1, 1, 3], [1, 1, 5], [1, 1, 7], [1, 2, 3], [1, 2, 9], [2, 1, 5], [2, 2, 1], […
尝试将表达式与 lambda 集成时出现问题
我正在尝试集成一个具有实数和复数值的表达式，将其定义为 lambda 表达式。积分变量是 kx，积分的结果解将在 x 和 y 维度上求值，但是在我积分并尝试…
慢速收敛积分的数值积分
我基本上对这种类型的数值积分有一个问题： import numpy as np import matplotlib.pyplot as plt from scipy.integrate import quad # Parameter T =…
与 c++ 的数值积分在具有固定常数离散化的给定网格上
我有以下问题： 我的 C++ 代码可以 计算两个函数f1(i1,i2,i3,i4) f2(j1,j2) 为每组 {i1,i2,i3,i4} 我得到 f1 的一些值，并且对于每组 {j1,j2} 我都会…
集成到matlab中
我如何求解 ZV 的以下方程 ZV = integral_{0}^{ZD}(z.*A)dz ： ZD = 40; z = [1:2:12];%depth A = [112,118,123,143,188,200];%area …
Matlab 求解器中的相对和绝对公差定义
我试图理解以下RelTol和AbsTol参数的Matlab定义： RelTol — This tolerance is a measure of the error relative to the size of each solution comp…
在java中实现gamma不完整
我在 Java 中实现 MATLAB 的 gammainc 时遇到问题。 我试着一部分一部分地做。我使用了这些函数，这些函数使用了我在网上获得的 Lanczos 近似来求解 …
AdaptIntegrate 与集成的性能
我想在一维中执行数值积分，其中被积函数是向量值。 integrate() 只允许标量被积函数，因此我需要多次调用它。 cubature 包似乎很适合，但对于一维积…
Matlab：如何计算函数在多个极限上的定积分？
假设： z = [0.4 0.5 0.75]' function y = myfunct(x) y = quad(@sin, 0, x) 我想计算 sin(x) 从 0 到 0.4、0.5 和 0.75 的定积分，使用： myfunct(z)…
MATLAB 中 R^5 超立方体的蒙特卡罗积分
我需要编写 MATLAB 代码，使用 Monte Carlo 在 R^5 超立方体上进行积分。当我有通用函数时，我有一个基本算法。但我需要积分的函数是： ∫dA A 是 R^5…
这个积分可以在 Matlab 或 Mathematica 中进行数值计算吗？
我希望能够完全以数值方式进行下面的积分。 其中 、 和 、 和 是常量，为简单起见，可以全部设置为 1。 对 x 的积分可以通过手动或使用 Mathematica …
求解受限于给出非负解的时滞微分方程 (DDE) 系统
在 MATLAB 中，ode45 有一个名为非负这将解决方案限制为非负。 他们甚至写了一篇论文，介绍这种方法的工作原理以及它为何不是那么愚蠢的事情每当 y_i …
