使用近似技术在一个或多个维度上集成函数的算法，而不是使用符号代数和微积分的精确闭合形式解决方案。包括自适应求积法，蒙特卡罗方法，有限元分析，马尔可夫链等概念。
我目前正在尝试在Python中计算以下表达式 该表达式来自 t 分布，v 是表示自由度的标量，mu 和 sigma 是向量和矩阵
为什么积分 scipy.multivariate_normal 会给出错误的概率？
我正在尝试在方形区域上整合独立的二元正态分布。数值积分与蒙特卡罗模拟不匹配。这里出了什么问题？ 导入numpy...
我有一个问题。 假设我们有一个依赖 3 个变量的函数： f(x,y,z) = x+y+z 我想对 x 和 y 进行积分，并保持 z 固定。我们可以取 z = 0 或其他任何值。 我试过了
MATLAB 和 SciPy 贝塞尔函数数值积分结果之间的差异
我希望在 SciPy 中复制以下 MATLAB 代码的结果。 MATLAB版本： f = @(x, y) 贝塞尔(0, 2, x.^2 + y.^2); 积分2(f, -0.1, 0.1, -0.1, 0.1) MATLAB 结果为： 答案= ...
我正在尝试用 Verilog 编写一个集成模块，它可以改变集成周期的持续时间。它正在编译，没有错误，并且信号的数据寄存器实际上存储...
我遇到了一个乍一看似乎根本没有问题的问题。 我想编写一个函数，它将采用发行版名称和参数作为参数，并返回预期的...
我正在尝试在Python中集成一些矩阵条目。我想避免循环，因为我的任务包括 1 个 Mio 模拟。我正在寻找能够有效解决我的问题的规范...
说我想整合 f=@(x)[1,0]； 积分2( @(x,y) f(x)*f(y)' ,0,1 ); 但这给出了错误 被积函数输出大小与输入大小不匹配。 尽管 f 是一个向量，但我正在...
我在数值积分标准化后验分布方面遇到问题。我的代码是： a1 <- 0.012 sigma_a1 <- 0.001 likelihood <- function(a_true) { (1 / (sqrt(2 * pi) * sig...
我在数值积分标准化后验分布方面遇到问题。我的代码是： ＃ 数据 a1 <- 0.012 sigma_a1 <- 0.001 likelihood <- function(a_true) { (1 / (sqrt(2 ...
有人可以帮我在python库环境中安装Odespy库吗？ https://github.com/NisangaN/odespy.git 我尝试使用 pip install，但它不起作用。我对编程有点陌生...
我想在 Matlab 中使用 gpuArray 实现二重积分，并将结果绘制为 128x128 点的 2D 图像。我在 Mathworks 中找到了 Joss Knight 分享的示例，标签为：“How to use ...
我想知道是否有可能将自适应时间步进器并行化，例如 Runge-Kutta 23 方法（或更高阶）与 Dask？ 假设我有 10 个 Dask 工作人员（每个工作人员有 1 个核心/等级），我分配 1...
我想要一个方程式的数值结果来做一些测试。 但它由无限积分组成，没有明确的解决方案。我应该逐点求解积分以获得结果吗？ 我要……
如何解决scipy.integrate.quad中特定被积函数的舍入误差问题？
我正在做下面的数值计算。 将 numpy 导入为 np 从 scipy.integrate 导入四边形 def E_k（间隙，q，mu）： return np.sqrt(abs(gap) * abs(gap) + (q * q - mu) * (q * q - mu)) 四边形（...
我必须在 3D 空间上进行大量集成。我对所述空间进行了（网格）离散化（例如，在每个方向上等距 30 个点），我知道我想要的函数值
我想定义一个函数，它接受一定长度的数组输入，然后为每个值找到一个值的最佳参数。 它看起来像这样： 这里 a, b, c 和 measu...
我想通过求解运动方程来跟踪一个小时的粒子路径。对于时间步长，t= np。排列 (0, 3600,10)。这是正确的吗？
我想通过求解运动方程来跟踪一个小时的粒子路径。对于时间步长， t= np。排列 (0, 3600,10)。 它是否正确？ 持续时间是否取决于时间步长。如果我...
我正在尝试模拟流体（气体或液体，例如空气或水）中的可操纵车辆（飞机、气垫船、船）。 我的 JavaScript 实现是基于 3D 刚体飞机模拟...
我正在尝试加速以下计算代码： 我只需要为 x > y 从 0 到 1 计算此函数（但需要非常高的离散化，如 dt = 0.001）。我已经矢量化了...