numerical-integration 相关问题

使用近似技术在一个或多个维度上集成函数的算法,而不是使用符号代数和微积分的精确闭合形式解决方案。包括自适应求积法,蒙特卡罗方法,有限元分析,马尔可夫链等概念。

通过推广矩形方法进行 3 维积分

def RectanglesPointMilieu(f,a,b,n): 间隔 = 1.* (b-a) / n 总和 = 0 对于 np.arange(a,b,interval) 中的 i: sumfct += f(i + 间隔/2。) 返回间隔 * sumfct 我该怎么办

回答 1 投票 0

C++ nan 不断提出集成[已关闭]

我正在尝试运行一个集成程序,但在计算时我不断收到nan。我不知道我的代码有什么问题。 #包括 #包括 你...

回答 2 投票 0

在紧凑支持 [0,alpha], alpha 上集成类似 beta 函数的函数 < 1

我想集成以下名为 betalog 的函数 克<- function(x,a,b){ if (a < 0 | b < 0) stop() temp <- (a-1)*log(x) + (b-1)*log(1-x) return( exp(

回答 1 投票 0

maple 中的 EllipticE 与 matematica 和 python 相比

尝试在maple中实现Euler Elastica,其中Jacobi椭圆函数和椭圆积分至关重要。 Jacobi 椭圆函数:JacobiSN、JacobiCN、JacobiDN 和 JacobiAM 似乎有效...

回答 1 投票 0

计算不规则物体的精确体积

对于作业,我必须使用数值积分技术来计算圆柱表面的体积 R3 中的 Ω={(x,y,z),其中 (x−0.5)² +(y−0.5)² ≤ 0.5² 且 0 ≤ z ≤ |ln(x+y)|}。 我用过

回答 1 投票 0

Scipy tplquad 语法

我希望使用 scipy 中的 tplquad 函数在 python 中执行以下积分: 到目前为止,我已经使用了以下语法: def 函数(z,y,x): 返回 (x**(b1 + x1 - 1))*(y**(b2 + x2 - 1))*(z**(...

回答 2 投票 0

三重积分:ValueError:负数不能升到分数幂

我正在尝试对函数进行三重积分。 从 scipy.integrate 导入 tplquad S = 40 P1=0.37 P2 = 0.43 P3V = .05 UND = 1 - (P1+P2+P3V) b1 = S*P1 b2 = S*P2 b3 = S*P3 b4 = S*UND x1 = 48 x...

回答 1 投票 0

计算余弦积分和 Ci(x)

我有Java基础背景,所以我对不同结构和实现的了解有限。 我想知道如何在Java中计算余弦积分。 具体来说,我有 f(x...

回答 1 投票 0

编写代码以在 Python 中集成一维菲涅耳衍射

我的任务是对以下 1d 菲涅尔衍射方程进行积分(红色部分): 重点是,您正在将孔径傅里叶变换为屏幕上的图案,但现在只需关注...

回答 3 投票 0

比较 Haskell 中的大数以通过快速检查

我制作了一个函数,使用 n 个梯形的梯形规则计算某个函数 g 从 a 到 b 的定积分。 即定积分 (-4) 5 (\x -> x) 5 = 4.

回答 1 投票 0

如何在python中使用tplquad、dblquad进行积分?

如何使用 scipy 在 python 中执行这种积分: $$\int_{0}^{1} f(x) \, dx \int_{0}^{x} g(y) \, dy \int_{0}^{x} h(z) \,dz $$ 尝试使用 tplquad,但我认为事实上两个内部

回答 1 投票 0

python 中的积分: y(t,E)=-A+2A*Integral [g(t,tau)*f(Ea)]dEa

我需要进行图中写的积分: 我想针对给定的 E 值绘制 y(t,E)。在我的代码中,我收到以下消息:“在 double_scalars 中遇到除以零” 我不明白...

回答 1 投票 0

“数组索引必须是正整数或逻辑值”

问题详情及线索 %针对此问题编写一个名为 NC.m 的脚本文件,该文件实现 %任意函数 f(x) 的 Newton-Cotes 积分法。它 %应该将 fu...

回答 2 投票 0

尝试积分函数给出“ufunc 循环不支持 Symbol 类型的参数 0”

我的包裹: 从 pylab 导入 plt 将 pandas 导入为 pd 将 numpy 导入为 np 将 sympy 导入为 sy 将 matplotlib.pyplot 导入为 plt 将日期时间导入为 dt 从 pylab 导入 mpl 将 scipy.integrate 导入为 sci

回答 1 投票 0

如何在Python中与x^2进行(梯形)积分?

我的任务是首先对 f(x)=x^2 进行积分,然后用 Python 进行梯形积分 将 numpy 导入为 np 将 matplotlib.pyplot 导入为 plt x = np.arange(-10,10) y = x**2 l=plt.plot(x,...

回答 3 投票 0

Python 指数积分:exp 中遇到溢出

这个问题可能已经被问过一百万次了;但我无法找出针对我的具体问题的解决方案。这是我的代码: 将 numpy 导入为 np 将 matplotlib.pyplot 导入为 plt 从 sympy 导入 Symb...

回答 1 投票 0

Sympy:当源于莱布尼茨规则时,Integral(0, (R, b, r)) 没有简化为零

我正在积分下微分。由于被积函数并不明确依赖于变量,因此莱布尼茨规则中的相应项不会出现。 然而,sympy 并不

回答 1 投票 0

从一组加速度数据到位置的积分

我正在尝试一个项目来整合加速度数据以获得位置的近似值。我从一组真正简单的数据开始,并不断加速。 来自 scipy。

回答 1 投票 0

RK4 与欧拉方法比较

我选修了计算物理学课程,并接受了一项作业来求解简谐振子 x_dot_dot = -x 的方程。我用 Julia 编写了这两种方法的脚本 函数 R...

回答 1 投票 0

在Python中数值求解二重积分

我正在尝试使用矩形面积来数值求解二重积分,其中积分看起来只是两个变量(theta 和 phi)的函数。 到目前为止,这是我的尝试: 来自我...

回答 1 投票 0

© www.soinside.com 2019 - 2024. All rights reserved.