使用scipy.integrate的quad进行数值积分结果中的虚假凹槽
使用scipy.integrate的quad进行数值积分时，虚假凹槽是指在积分结果中出现的不合理的凹槽形状，即积分结果在某些区间内出现了负值或者不符合预期的形状。
虚假凹槽可能是由于积分函数在某些区间内具有非光滑或者不连续的特性导致的。这种情况下，quad函数可能无法准确地捕捉到函数的真实形状，从而导致虚假凹槽的出现。
为了解决虚假凹槽的问题，可以尝试以下方法：
- 调整积分方法：scipy.integrate的quad函数提供了多种积分方法，可以尝试使用不同的方法进行积分，例如使用"quad"方法的参数full_output=True，可以返回更详细的积分信息，包括积分误差估计值。
- 增加积分精度：可以通过调整quad函数的参数，例如设置参数epsabs和epsrel来增加积分的精度，从而减少虚假凹槽的出现。
- 分段积分：如果积分函数在某些区间内具有非光滑或者不连续的特性，可以将积分区间分成多个子区间，分别进行积分，然后将结果进行合并。
- 使用其他数值积分方法：除了scipy.integrate的quad函数，还可以尝试使用其他数值积分方法，例如trapz、simps等，看是否能够得到更准确的积分结果。
总之，虚假凹槽是数值积分中常见的问题，可以通过调整积分方法、增加积分精度、分段积分或者尝试其他数值积分方法来解决。具体的解决方法需要根据具体情况进行调试和优化。
腾讯云相关产品和产品介绍链接地址：
- 腾讯云计算服务：https://cloud.tencent.com/product/cvm
- 腾讯云数据库：https://cloud.tencent.com/product/cdb
- 腾讯云服务器运维：https://cloud.tencent.com/product/cvm
- 腾讯云音视频处理：https://cloud.tencent.com/product/vod
- 腾讯云人工智能：https://cloud.tencent.com/product/ai
- 腾讯云物联网：https://cloud.tencent.com/product/iot
- 腾讯云移动开发：https://cloud.tencent.com/product/mc
- 腾讯云存储：https://cloud.tencent.com/product/cos
- 腾讯云区块链：https://cloud.tencent.com/product/baas
- 腾讯云元宇宙：https://cloud.tencent.com/product/mu
相关·内容
我试图通过scipy.integrate.quad()计算这样一个积分(实际上是指数分布的cdf及其pdf)。
import numpy as np
from scipy.integrate import quad
def g(x):
return .5 * np.exp(-.5 * x)
print quad(g, a=0., b=np.inf)
print quad(g, a=0., b=10**6)
print quad(g, a=0., b=10**5)
print quad(g, a=0., b=10**4)
结果如下：
(1.0, 3.5807346295637055e
浏览 1提问于2016-09-15得票数 4
回答已采纳
1回答
如何正确显示图形上的积分？
、、、、
我计算了积分，我想在图上显示它，但我想知道它应该如何正确地放置在图上。在我看来，光靠plt.plot()是不够的，或者我错了，我想知道在图中显示这个结果的正确方法。
import matplotlib.pyplot as plt
import numpy as np
from scipy.integrate import quad
def integral(x, a, b):
return a * np.log(x + b)
a = 3
b = 2
I = quad(integral, 1, 5, args=(a, b))
print(I)
plt.plot()
plt.
浏览 3提问于2021-05-31得票数 1
回答已采纳
在依赖于两个变量(q,r)的函数中执行双定积分时，我遇到了一个问题，其中还有一个额外的积分。我想用高斯函数加权的函数是：
F(q,r)=f(q,r)+int_{0,r}(h(q,r')dr')
而in必须再次积分，才能与高斯加权：
I(q)=int_{0,inf}(F(q,r)^2*g(r)dr)
高斯g(r)以坐标R为中心。
您可以看到的主要问题是，我将数组与标量混合在一起。使用用于高斯(np.ogrid和轴上和)的相同方法可能是一个解决方案，但我不知道如何实现它。
import numpy as np
from scipy.integrate import quad
imp
浏览 2提问于2013-09-27得票数 0
我正在编写一个脚本，用于将峰值形状与Python中的光谱数据相匹配，使用Scipy、Numpy和Matplotlib。它可以同时适应多个峰值。峰值轮廓(目前)是伪Voigt，这是一个高斯(又名法线)和洛伦兹(又名柯西)分布的线性组合。
我有一个选项开关，可以让软件优化Gaussian和Lorentzian的贡献，或者将其设置为固定值(0=纯高斯和1=纯Lorentzian)。按其应有的工作原理，绘制合适的山峰看上去和预期的一样。当我尝试用scipy.integrate计算峰值的积分时，问题就开始了。
到目前为止，我尝试了scipy.integrate.quad、scipy.integrate.
浏览 1提问于2016-06-06得票数 1
回答已采纳
我试图计算Ecut，inf上1/((s - q02)*(s - q2))的整数(超过s)的主值，其中q02 < Ecut < q2。用手(或Mathematica)做主值可以得到一般结果。
ln((Q2-E)/(Ecut-q02))/ (q02 -q2)
在下面的具体示例中，这将给出结果-1.58637*10^-11。还可以通过将积分分成两部分，积分到q2 - eps，然后从q2 + eps开始，再加上这两个结果(偏差应该取消)，从而得到相同的结果。取eps越小越好，应能恢复上述结果。但是，当我使用quad在the中实现这一点时，我的结果会收敛到错误的结果6.04685e-11，如
浏览 0提问于2017-04-10得票数 1
回答已采纳
5回答
我试图计算一个函数的定积分，这个函数只在一个变量上有多个变量。这有点像我的代码-
from scipy.integrate import quad
import numpy as np
def integrand(x,y):
return x*np.exp(x/y)
quad(integrand, 1,2, args=())
它返回这个类型的错误：
TypeError: integrand() takes exactly 2 arguments (1 given)
不过，如果我把一个数字放进args里就行了。但我不想，因为我希望你和你一样，而不是一个数字。有人知道怎么做吗？
编辑:对
浏览 3提问于2014-08-14得票数 2
回答已采纳
我正在尝试用scipy.integrate中的NumPy和quad数值求解下面的积分。代码在某种程度上是有效的，但我在结果数据中得到了虚假的凹口：
有谁知道为什么会发生这种情况，以及如何获得正确的平滑结果？
以下是Python的原始代码：
#!/usr/bin/env python
import numpy as np
from scipy.integrate import quad
import matplotlib.pyplot as pl
numpts = 300
t_min = -4
t_max = 100
tt = np.linspace(t_min, t_max,
浏览 0提问于2017-02-17得票数 0
回答已采纳
我试图通过为epsabs指定参数scipy.integrate.quad来更精确地计算积分，比如我们将函数sin(x) / x^2从1e-16集成到1.0。
from scipy.integrate import quad
import numpy
integrand = lambda x: numpy.sin(x) / x ** 2
integral = quad(integrand, 1e-16, 1.0)
这给了我们
(36.760078801255595, 0.01091187908038005)
为了使结果更加精确，我们通过epsabs指定了绝对误差容限。
from scipy.i
浏览 2提问于2018-12-21得票数 3
回答已采纳
我在使用scipy.integrate.quad时发现了一种奇怪的行为。这种行为也出现在Octave的四重函数中，这使我相信它可能与QUADPACK本身有关。有趣的是，使用完全相同的Octave代码，这种行为不会出现在MATLAB中。
关于这个问题。我是数值积分一个对数正态分布在不同的边界上。对于F是对数正规的cdf，a是下界，b是上界，我发现在某些条件下，
当b是“很大的数”时，积分(F，a，b) =0
当b是np.inf时，积分(F，a，b) =正确的极限。(或者只是八达通的Inf。)
下面是一些示例代码来显示它的作用：
from __future__ import division
im
浏览 6提问于2014-06-02得票数 2
回答已采纳
我正在尝试计算两个概率分布之间的。为此，我需要执行此积分。
下面是我的简化代码，它目前失败了：
from scipy.integrate import quad
import numpy as np
def f(x):
return sum([ps[idx]*lambdas[idx]*np.exp(- lambdas[idx] * x) for idx in range(len(ps))])
def g(x):
return scipy.stats.weibull_min.pdf(x, c=c)
c = 0.9
ps = [1]
lambdas = [1]
eps
浏览 0提问于2019-10-20得票数 0
2回答
这里有一个有结果的例子：
我通过高斯分布(mu=800，sigma=1)与~+-2σppf和从-inf到+inf的相同积分进行积分。由于某些原因，第二次积分的结果为零，但在实践中，它应该更准确。
有人能解释一下，为什么会发生这样的异常，或者我在哪里犯了错误？
代码：
from scipy.integrate import quad
import numpy as np
from scipy.stats import norm
def integrand(x):
return x*norm.pdf(x, 800, 1)
print quad(integrand, norm.ppf(0.
浏览 3提问于2016-01-19得票数 4
回答已采纳
在Python集成中，我遇到了一个问题，用已知的分析解决方案返回积分的不正确值。所讨论的积分是
对于我使用的σ值(1e-15)，这个积分的解的值是~ 1.25e-45。然而，当我使用has集成包来计算时，我得到了零，我认为这与计算所需的精度有关。
#scipy method
import numpy as np
from scipy.integrate import quad
sigma = 1e-15
f = lambda x: (x**2) * np.exp(-x**2/(2*sigma**2))
#perform the integral and print the result
浏览 6提问于2022-08-25得票数 3
1回答
具有可变边界的超数组集成
、、、
我试图在一个点的列表上集成一个函数，并将整个数组传递给一个集成函数，以便将事情向量化。首先，打电话给scipy.integrate.quad太慢了，因为我有大约10000点要集成。使用scipy.integrate.romberg的速度要快得多，几乎是瞬间的，而四边形是慢的，因为你必须循环它或矢量化它。
我的函数很复杂，但是为了证明目的，假设我想从a到b集成x^2，但是x是一个计算x的标量数组。
import numpy as np
from scipy.integrate import quad, romberg
def integrand(x, y):
return x
浏览 7提问于2013-12-05得票数 1
我有一个复杂的函数，它的集成行为我无法理解：
pointsWhatever points确实指定了可能的奇异点的间隔，或者对上面的答案，如何解释2.1工作，而不是2？(还有1，5，.)请注意，奇点的唯一主题是0，但情况并非如此，因为没有/ w，函数在那里是0，即使删除了/ w，/ w仍然存在。
import numpy as np
from scipy.integrate import quad
def fn(w, mu=11.316582914572866):
return (np.exp(-1 / (1 - ((w - mu) * (np.abs(w - mu) < .
浏览 2提问于2020-11-24得票数 0
回答已采纳
3回答
有一个函数决定了圆孔的夫琅和费衍射图的强度。()
函数在距离x= -3.8317，3.8317中的积分必须约为83.8% (假设I0为100)，当将距离增加到-13.33，13.33时，约为95%。但是当我在python中使用积分时，答案是错误的。我不知道我的代码出了什么问题
from scipy.integrate import quad
from scipy import special as sp
I0=100.0
dist=3.8317
I= quad(lambda x:( I0*((2*sp.j1(x)/x)**2)) , -dist, dist)[0]
print I
积分的结果
浏览 3提问于2014-06-28得票数 4
回答已采纳
1回答
我试图用Python来解决积分问题，我做了以下代码
from sympy import integrate
import sympy as sy
from scipy.integrate import quad
import math
k_2=Symbol('k_2')
k_1=Symbol('k_1')
t_d=Symbol('t_d')
C_0=Symbol('C_0')
y=Symbol('y')
f=lambda y:((k_2+C_0*sy.exp(-y/t_d))/(k_2+C_0))**(k_1*t_
浏览 10提问于2021-11-03得票数 0
这个函数从0到inf是正数，为什么给大于1000会产生0，不合理。
import scipy.integrate as integrate
from math import *
integrate.quad(lambda x: x*exp(-x**2/2)/sqrt(2*pi), 0.0, 1000)
Out[52]:
(0.3989422804014328, 1.6471510195390376e-11)
integrate.quad(lambda x: x*exp(-x**2/2)/sqrt(2*pi), 0.0, 100000)
Out[54]:
(0.0, 0.0)
浏览 0提问于2017-11-09得票数 0
回答已采纳
我想使用scipy来集成以下函数
此函数w上的列表为w=linspace(0,10,300)。W是实验数据，所以我把它作为一个列表。我已经尝试了以下代码
from scipy.integrate import quad
import numpy as np
def integrand(t,R):
return np.exp(R*t)
R=0.04
w=np.linspace(0,10,300)
print(quad(integrand, 0, 4, args=(R)))
问题是，每当我试图将w与函数scipy.integrate相乘时，它都不起作用。有没有办法克服这个问题
浏览 42提问于2019-11-02得票数 0
回答已采纳
对于计算这种集成的任何帮助，F函数都是使用f函数定义的，这涉及到第一个集成，最后是集成F。
from scipy.integrate import quad
f = lambda x,a : a**2*x
def F(s,a):
return quad(f,0,s,args=(a,))
quad(F,0,5,args=(4,))
得到了错误：
2 def F(s,a):
3 return quad(f,0,s,args=(a,))
----> 4 quad(F,0,5,args=(4,))
5
446 if points
浏览 0提问于2019-07-25得票数 0
回答已采纳
1回答
我编写的代码看起来可能很愚蠢，因为它是一个派生函数的集成。因为它是我正在编写的关于声学分析的其他代码的基本基础。这个分析包含了乘法中不同的导数函数的积分。为此，我使用SciPy进行集成，使用渐近进行区分。但是它给出了一个显示TypeError的错误(“无法将表达式转换为浮点数”)。下面是我写的代码。希望能找到解决办法。
import sympy
from sympy import *
from scipy.integrate import quad
var('r')
def diff(r):
r=symbols('x')
Z = 64.25 *
浏览 3提问于2022-02-06得票数 0
如果我运行类似下面的代码，在椭圆的小圆弧和大圆弧中传递3.05和2.23，由圆弧形成的50度角，我如何能够将2.531432761012828的输出作为弧长，并将其传递回来以求解t？谢谢!
import math
from scipy.integrate import quad
import numpy as np
t = math.atan(3.05*math.tan(5*math.pi/18)/2.23)
Dx = lambda t: -3.05 * np.sin(t)
Dy = lambda t: 2.23 * np.cos(t)
quad(lambda t: np.sqrt(Dx
浏览 14提问于2017-07-07得票数 0
回答已采纳
我需要集成一个精度很高的函数，例如：
from decimal import *
from scipy.integrate import quad
getcontext().prec = 30
a = Decimal('1.1') + Decimal('1.2')
print type(a)
def f(X):
return quad(lambda x: X, 0, 1)[0]
b = f(a)
print "{0:.30f}".format(b), type(b)
我得到了
<class 'decimal.Deci
浏览 2提问于2015-12-10得票数 0
回答已采纳
我有一个函数，它应该计算一个积分，接受某个函数作为输入。我想让代码计算一个确定的积分为：<some function, in terms of x. e.g., 3*x or 3*x*(1-x), etc.> * np.sin(np.pi * x))。我用枕木来做这个：
import scipy.integrate as integrate
def calculate(a):
test = integrate.quad(a*np.sin(np.pi * x), 0, 1)
return test
a = lambda x: 3*x
calculate(a)
现在，
浏览 4提问于2020-11-20得票数 0
回答已采纳
1回答
我在计算积分
sin(x)/x , x = [0,inf]
我所做的工作如下：
import math
from scipy.integrate import quad
t = float("inf")
def integrand(x):
return (math.sin(x))/x
ans, err = quad(integrand, 0, t)
print(ans)
然而，我得到了错误的答案:它应该是Pi/2，并且来自Python的答案是2.247。
你知道有什么问题吗？
浏览 0提问于2017-10-01得票数 1
回答已采纳
我正在使用scipy.integrate中的quad从一个对象获取一个有限范围内的积分。假设目标对象在吹风中： ∫expm(A*X).expm(B*X)dx 其中A和B都是数值矩阵。 为了解决这个问题，我使用了blow代码： from scipy.integrate import quad
from scipy.linalg import expm
import numpy as np
def integrand(X, A, B):
return np.dot(expm(A*X),expm(B*X))
A = np.array([[1, 2], [3, 4]])
B = np.
浏览 63提问于2019-02-16得票数 0
回答已采纳
我一直在使用Python中的符号表达式，并且已经达到了一个我希望在一定的时间间隔内集成的表达式。表达式包含pi。
麻烦的是，我还没有弄清楚如何将这个表达式转换为一个可以作为参数输入到scipy.integrate.quad的函数。我的代码的相关部分如下：
from sympy import *
import numpy as np
import scipy.integrate as integrate
from sympy.utilities.lambdify import lambdify
# this defines the symbols that
# we will be using
浏览 0提问于2018-11-28得票数 0
回答已采纳
我试图用scipy.integrate中的四值函数来集成一个奇点函数，但我没有得到想要的答案。以下是代码：
from scipy.integrate import quad
import numpy as np
def fun(x):
return 1./(1-x**2)
quad(fun, -2, 2, points=[-1, 1])
这将导致IntegrationWarning并返回有关0.4的值。
这个函数的极点是-1,1。答案应该是大约1.09 (用钢笔和纸计算)。
浏览 0提问于2018-10-07得票数 5
回答已采纳
1回答
scipy.integrate的quad需要参数func，a，b，其中func是要集成的函数，a和b分别是下、上的积分极限。A和b必须是数字。
在这种情况下，我需要对数十万个不同的a，b函数的积分进行求值，并对结果进行求和。这需要很长时间才能循环通过。我试着给出a和b的四叉数组，希望quad能返回相应的数组，但这是行不通的。
这里有一段代码，说明了我正在尝试做什么，这两个Python循环都能工作，但速度非常慢，而且我试图将其矢量化，但不起作用。对于如何解决这个问题，有什么建议可以用一种快速的方式解决吗？
import numpy as np
from scipy.integrate impor
浏览 2提问于2015-03-13得票数 2
回答已采纳
我试图使用数值积分方法在我的程序中对任意(已知的代码时)函数进行数值集成。我正在使用Python2.5.2以及SciPy的数值集成包。为了对它有感觉，我决定尝试集成sin(x)并观察到这种行为-
>>> from math import pi
>>> from scipy.integrate import quad
>>> from math import sin
>>> def integrand(x):
... return sin(x)
...
>>> quad(integrand, -p
浏览 3提问于2009-02-24得票数 5
回答已采纳
我目前正在尝试执行高斯函数的定积分，当我确信情况并非如此时，我得到的答案是0。
这让我不禁要问，在执行定积分时，quad函数到底能做什么，有没有限制？我在正确的应用程序中使用quad了吗？quad到底是如何找到一个整数的呢？
import math
from scipy.integrate import quad
def g(λ,a,u,o):
return a*math.exp((λ-u)**2/(-2*o**2))
exc = quad(g, 4000, 8000, args=(1,6700,2.125))[0]
print(exc)
我绘制了这个高斯曲线，所以我知道它在我设
浏览 21提问于2021-10-06得票数 0
1回答
我需要用MATLAB集成一个表达式。这些功能之间有什么区别？
QUAD，QUAD8，QUADL，DBLQUAD
浏览 2提问于2011-02-23得票数 0
2回答
from scipy.integrate import quad
def integrand(a, b):
return a * x ** 2 + b
a = 2
b = 1
I = quad(integrand, 0, 1, args=(a,b))
I
这是我的程序。当我尝试运行它时，它显示错误：
integrand () takes 2 positional arguments but 3 were given ....
我不明白为什么当只有两个变量，即a和b时，它会要求3个参数。
有谁可以帮我？有人能澄清我的疑虑吗？
浏览 0提问于2018-08-22得票数 5
我使用Python3使用scipy的integrate.quad()函数对给定输入s的表达式(t^s / cosh^2(t))进行积分。我得到了OverflowError:数学范围错误。 我的代码如下所示： from math import inf
from math import cosh
from scipy.integrate import quad
import numpy as np
def Integral_Expression(s):
ans, err = quad(lambda t: (np.power(t,s))/(np.power(cosh(t),2)), 0,
浏览 113提问于2020-07-03得票数 0
回答已采纳
4回答
我想解决这个问题：
我不想用"int"，我想用"quad“家族(quad，dblquad，triplequad)，但我不能。你能帮我吗？
浏览 0提问于2010-07-27得票数 0
在这种情况下， from scipy.integrate import quad
import numpy as np
exponential_distribution = lambda x, lam: lam*np.exp(-lam*x)
result = quad(exponential_distribution, 0.25, 0.75, args=0.1)[0] 我正在想像这样的事情 from magic_library import integrate_distribution
result = integrate_distribution(distribution='e
浏览 14提问于2019-12-15得票数 2
回答已采纳
我现在正在使用scipy.integrate.quad来成功地集成一些真实的被积分器。现在出现了一个情况，我需要对一个复杂的被积函数进行积分。quad似乎不能像其他scipy.integrate例程那样做到这一点，所以我问:有没有办法用quad来积分一个复杂的被积函数，而不必分离实部和虚部的积分？
浏览 13提问于2011-05-11得票数 29
回答已采纳
我如何最小化一个函数(不受约束)，尊重a和a1？示例(这是一个简单的示例，用于I uderstand、numpy和py)：
import numpy as np
from scipy.integrate import *
from scipy.optimize import *
def function(a):
return(quad(lambda t: ((np.cos(a[0]))*(np.sin(a[1]))*t),0,3))
我试过：
l=np.array([0.1,0.2])
res=minimize(function,l, method='nelder-mead&#
浏览 2提问于2016-10-10得票数 0
回答已采纳
1回答
我刚开始使用Python进行集成，所以请容忍我。代码很短，所以我会发布全部内容：
h = 6.626E-34
c = 3.0E8
k = 1.3806E-23
import scipy.integrate
from scipy.integrate import quad
from math import exp, pi
def f(x, T):
return 4*pi*h*c**2/(x**5*exp(h*c/(x*k*T))-1)
res = quad(f, 0.0, 000000.3)
print(res)
我在这里要做的是在x上进行积分，并让python返回一个T的数值。显然
浏览 0提问于2020-04-29得票数 3
回答已采纳
2回答
我必须证明：
烦人的是c_i等于函数G的积分，这是我的尝试。
import numpy as np
from scipy.integrate import quad
def G(x,n):
P = (np.sqrt(735))*(np.sqrt(2))*np.sin(n*np.pi*x)*((x**3.0) - (11.0/7.0)*(x**2.0) + (4.0/7.0)*(x))
return P
def Sum(x, n):
i = 1
S = 0
I, err = quad(G, 0, 1, args=(i))
while
浏览 0提问于2013-10-24得票数 4
回答已采纳
我在试着找一个积分
b = ((8*WL* Requ)/(pie* WC *Eequ))^.5
z = quad( 'exp(sqrt(1-((x.^2)./b^2) ))',-b,b)
为此，我使用了内置命令的辛普森求积函数(quad)。我得到了正确的结果，但我的问题是，我希望quad命令从我的代码中编写的(b)公式中自动拾取(b)的变化值。但它没有做到这一点。我必须将(b)的每个变化值放入函数中。
浏览 0提问于2011-11-27得票数 0
回答已采纳
1回答
我正在尝试使用Scipy优化以下问题
目前，我已经定义了我的方程，并尝试使用四次函数和最小化来求解方程。
import numpy as np
from scipy.integrate import quad
from scipy.optimize import minimize
def g_x(x, cv):
cv2 = np.square(cv)
return (x+1-np.sqrt(np.square(x)+2*cv2*x+1))/(1-cv2)
def g_approx(x, c):
return (c*x)/(1+c*x)
def inte
浏览 2提问于2022-07-15得票数 1
回答已采纳
1回答
慢收敛积分的数值积分
、、、
基本上，这种类型的数值积分有一个问题：
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import quad
# Parameter
T = 94
Eps_inf = 3.05
Eps_s = 19.184544603857724
tau_0 = 1.27*10**-16
tau_CD = 0.34580390675331274
gamma = 0.49
C_pek = 1/Eps_inf - 1/Eps_s
hbar = 6.582119569*10**-16
# Functions
def
浏览 6提问于2022-02-22得票数 0
回答已采纳
1回答
我试图用渐近解出以下积分：
V是速度，C_i是时间步长t_0的浓度。到目前为止，这就是我所拥有的：
import sympy as smp
from scipy.integrate import quad
to,x = smp.symbols(('to','x'), real=True)
def f(to,x,c,v,t):
return c*smp.DiracDelta((x/v) - t + to)
c_arr = 0.5
v = 0.1
x = 10
tt = x/v
t_arr = np.arange(0,1000,1)
integ
浏览 3提问于2022-07-07得票数 1
回答已采纳
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，因为计算精度不够。 from scipy.integrate import quad
import pandas as pd
ans, err = quad(f, -np.inf, np.inf) 你能帮我正确计算这样的
浏览 29提问于2020-09-08得票数 0
回答已采纳
我使用scipy.integrate来计算定积分。我已经编写了以下代码片段来评估概率密度函数( PDF )和累积分布函数( CDF )，以评估PDF中的CDF，如下所示： inf = 10000
def f_normal(m, tau): #normal pdf with mean np.sqrt(m- 1/2) and var 1/pi
tmp = 1/np.sqrt(np.pi) * np.exp( - ( tau - np.sqrt(m- 1/2) )**2 )
return tmp
def F_normal(m, x): #integral of f_no
浏览 18提问于2020-08-31得票数 1
1回答
在python中集成无限限制
、、、、
我正面临着在(0，无限)限制内进行积分的困难。我正在尝试做这个Jupyter笔记本(不知道这是否是相关信息)。代码如下： import numpy as np
import matplotlib.pyplot as plt
import scipy as sc
import math
from scipy.integrate import quad
def integrand3(x,z,m,n):
return ((np.cosh(x))**m)*((np.sinh(x))**n)*(np.exp(-z*np.cosh(x)))
def IgK0(z,m,n):
IntK
浏览 32提问于2020-09-07得票数 1
回答已采纳
2回答
python中的振荡积分
、、、
我编写了下面的代码来绘制出光学元件的光强，这基本上是入射场上的球面Fourier积分，所以它有一个Bessel函数。其参数依赖于集成变量(x)和绘图变量(r)。
from sympy import *
import matplotlib.pyplot as plt
import numpy as np
from scipy.integrate import quad
from scipy.special import jn
#constants
mm = 1
um = 1e-3 * mm
nm = 1e-6 * mm
wavelength = 532*nm
klaser = 2*np.pi
浏览 0提问于2018-03-01得票数 1
回答已采纳
我在做简单的积分，唯一的事情是我想把'n‘作为一个变量。我如何才能在集成的同时做到这一点呢？
import numpy as np
import matplotlib as mpl
from matplotlib import pyplot as plt
import scipy.integrate as integrate
from scipy.integrate import quad
import math as m
y = lambda t: 3*t
T = 4 #period
n = 1
w = 2*np.pi*n/T
#for odd functions
def
浏览 1提问于2020-02-09得票数 0
2回答
递归积分
、、
这是我正在研究的实际问题的简化版本。
假设我们有一系列(数学)函数，f_0，f_1，f_2，.
我们修正了一个函数f_0，比如说，f_0(x)=x，在这个简单的例子中，f_n(x) =1/(n+1)！)x^{n+1}是否有一种优雅而高效的方法来递归地对此进行编程，从而使f_n返回函数给一个任意的函数f_0？
我一开始试图返回f_2，但这已经失败了：
from scipy.integrate import quad
f_0=lambda x: x
f_1=lambda x: quad(f_0,0,x)
f_2=lambda x: quad(f_1,0,x)
返回错误
error: Sup
浏览 12提问于2017-05-15得票数 3
回答已采纳
我已经设法使用scipy.integrate.quad为我的随机过程类编写了几行代码
我有标准布朗运动的马尔可夫转移函数
import numpy as np
def p(x,t):
return (1/np.sqrt(2*np.pi*t))*np.exp(-x**2/(2*t))
但我想计算一下，我将用不起作用的代码来写。我这样写，这样我们就可以在不使用latex的情况下理解这个问题。
from scipy.integrate import quad
integral = quad(quad(p(y-x),1,np.inf)*p(x,1),1,np.inf)
你可能注意到了，问
浏览 8提问于2018-07-30得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
