在Python中使用dblquad对函数的正弦进行积分
在Python中，使用dblquad函数可以对函数的正弦进行积分。dblquad函数是scipy库中的一个函数，用于计算二重积分。
具体使用方法如下：
- 导入所需的库：
from scipy import integrate
import numpy as np- 定义要积分的函数：
def f(x, y):
return np.sin(x)- 调用dblquad函数进行积分：
result, error = integrate.dblquad(f, 0, np.pi, lambda x: 0, lambda x: np.pi)在上述代码中，f是要积分的函数，0和np.pi是积分的上下限，lambda x: 0和lambda x: np.pi是积分的边界函数。
- 输出积分结果：
print("积分结果：", result)
print("误差：", error)上述代码将输出积分结果和误差。
dblquad函数的优势是可以方便地进行二重积分计算，适用于各种需要对函数进行积分的场景。
腾讯云相关产品和产品介绍链接地址：
- 腾讯云产品：云服务器（https://cloud.tencent.com/product/cvm）
- 腾讯云产品：云数据库 MySQL 版（https://cloud.tencent.com/product/cdb_mysql）
- 腾讯云产品：云函数（https://cloud.tencent.com/product/scf）
- 腾讯云产品：人工智能（https://cloud.tencent.com/product/ai）
- 腾讯云产品：物联网（https://cloud.tencent.com/product/iotexplorer）
- 腾讯云产品：移动开发（https://cloud.tencent.com/product/mobdev）
- 腾讯云产品：对象存储（https://cloud.tencent.com/product/cos）
- 腾讯云产品：区块链（https://cloud.tencent.com/product/baas）
- 腾讯云产品：元宇宙（https://cloud.tencent.com/product/mu）
请注意，以上链接仅供参考，具体产品选择应根据实际需求进行评估。
相关·内容
我有一个关于加速度a( t )的函数，如果我对t积分两次，就可以找到位置x(t)。我正试图找到t= 10秒的位置。
因为a(t)不是多变量，所以我在使用Scipy dblquad函数来计算我需要的二重积分时遇到了问题。请看我到目前为止所拥有的：
def a(t):
return (2.5 / (1 + math.exp((t-8)/0.8)))
def upperbound():
return 10
def lowerbound():
return 0
x = dblquad(a,0,10,lowerbound,upperbound)
浏览 3提问于2017-10-26得票数 3
1回答
虚数的四次积分
、、、
我要计算包含虚数的隐函数的积分。
其中f(iz)类似于：
g(ix)类似于：
我想用数字来计算。Python scipy.quad不计算虚数的积分(在下面的代码1中解释)。Quadpy也没有效率，因为它传递整个numpy数组，而不是单个整值(在下面的代码2中解释)，因此需要额外的操作。因此，我正在考虑用如下方式来划分积分(其中Re是实部，I是虚部)：
并展开上述方程：
我能这么做吗？
下面是代码，我展示了两种方法来解决这个问题。
代码1：
第一步接近scipy.quad。根据：，我尝试将实部和虚部分开，并分别计算它们的积分值：
from scipy
浏览 4提问于2021-05-26得票数 0
回答已采纳
我想使用dblquad集成Int sin(f(x,y,a,b)) dx dy，其中a和b是两个参数。 下面显示了一个示例，其中x和y都从0到10，并且a= 0，b=1。 dblquad(f*np.sin(f),0,10,0,10,args=(0,1)) 但是，这会导致输出 Traceback (most recent call last):
File "./update.py", line 60, in <module>
dblquad(f*np.sin(f),0,10,0,10,args=(0,1))
AttributeError: 'fun
浏览 15提问于2019-01-03得票数 0
4回答
我想解决这个问题：
我不想用"int"，我想用"quad“家族(quad，dblquad，triplequad)，但我不能。你能帮我吗？
浏览 0提问于2010-07-27得票数 0
1回答
我想双倍集成一个函数。但是当使用dblquad而不是scipy.integrate和matlab时，我得到了不同的结果。用于双重集成的python实现如下所示：
###Python implementation##
import numpy as np
from scipy.integrate import dblquad
def InitialCondition(x_b, y_b, m10, m20, N0):
IC = np.zeros((len(x_b)-1,len(y_b)-1))
for i in xrange(len(x_b) - 1):
for
浏览 1提问于2016-01-28得票数 1
回答已采纳
4回答
我试图从高斯光束中求出光圈中的功率，在高斯光束中，光圈与光束中心相距。解决方案是以下方程() (对不起，这里没有LaTeX )：
Wz是一个常数，加上a和r，我不知道如何用MATLAB来做这样的事情。有人有什么建议吗？我知道有一个dblquad()函数，但它假设集成的极限是固定的，而不是相互依赖的。
浏览 9提问于2012-01-10得票数 2
我有一个巨大的功能要整合：
syms x y
f=(228155022448185.*(cos((2.*pi).*y)./exp(131738205584307./(35184372088832*x)) - 1)*(cos((8.*pi).*y)/exp(131738205584307./(8796093022208*x)) - 1)*(cos((8.*pi).*y)/exp(131738205584307./(8796093022208.*x)) + cos((18.*pi).*y)/exp(1185643850258763./(35184372088832.*x)) -2)/(18014
浏览 3提问于2013-05-28得票数 0
1回答
我试着整合这一点：
integrate(integrate(integrate(2*sin(z)*cos(atan((2*cos(y)-0.5+x)/(2*sin(y))))，y,0，pi/2)，x,0,1)，z,0，pi/2)；
Wolfram找到了解决方案，但我想控制它的准确性。我尝试使用tplquad，但是有一些错误。
def f(x,y,z):
return 2*sin(z)*cos(atan((2*cos(y)-0.5+x)/(2*sin(y))))
tplquad(f,0,1,0,pi/2,0,pi/2)
这些错误是：
“"/usr/lib/python2.7/dis
浏览 1提问于2016-11-04得票数 0
回答已采纳
1回答
我想在矩形区域(-pi/2, pi/2)x(-pi/2, pi/2)上计算二重积分。
我使用了int函数：
double(int(int(...)));
问题是运行时间太长了。但我不需要这么准确的结果。我如何以较低的精度计算这个积分，以便I可以减少运行时间？
浏览 3提问于2011-11-21得票数 1
我试图用n阶的高斯求积来计算一个双积分(在一个三角形上，节点在(0,0)，(0,1)，(1,0))。然而，运行
import scipy.integrate as integrate
f = lambda x,y: x+y
inside = lambda x: integrate.fixed_quad(f, 0, 1-x, args=(x,), n=5)[0]
outside = integrate.fixed_quad(inside, 0, 1, n=5)
给出
追溯(最近一次调用)：文件""，第1行，文件"/Users/username/anaconda/lib/
浏览 2提问于2017-02-16得票数 4
回答已采纳
2回答
下面的示例出现在Scipy的integration参考页面上。
from scipy import integrate
N = 5
def f(t, x):
return np.exp(-x*t) / t**N
integrate.nquad(f, [[1, np.inf],[0, np.inf]])
以下是从我的IPython笔记本(在cloud.sagemath.com上)得到的错误：
据我猜测，cloud.sagemath.com没有升级到Scipy的最新版本，因此缺乏nquad模块。但是，我所需要的只是对两个变量进行集成，因此希望使用已经在云中可用的dblquad。
因此，
浏览 3提问于2014-05-06得票数 6
回答已采纳
不言而喻--根据的说法，通用的集成方法是“四边形”。命名约定继续适用于双重集成(dblquad)、三重集成(tplquad)等。
据我所知，积分和术语"quad“之间没有数学上的联系。那么，为什么scipy会这样命名这些方法呢？
浏览 1提问于2020-04-09得票数 0
1回答
考虑一般2D高斯函数，中心是(0.5,0.5)，
A*exp(-a*(-0.5 + x)**2-b*(-0.5 + x)*(-0.5 + y)-c*(-0.5 + y)**2)
其中协方差矩阵可以用系数a、b和c作为
2a和B
B& 2c
逆时针旋转45度
(a-b+c)及(a-c)
(A)& (a+b+c)
但是，给定a=1.25、b=0和c=10000，并使用Python在单元平方上进行集成，
import numpy as np
import matplotlib.pyplot as plt
a=1.25
b=0
c=10000
d=(a-b+c)/2
e=a-c
浏览 1提问于2020-12-14得票数 0
回答已采纳
1回答
双积分与sipy
、、、、
我想问一下如何用python计算这个积分。
from scipy import integrate
Hight_above_ground = 2.5
surface_ground = 81
surface_collector = 65
f = (np.cos(tilt)*np.cos(tilt)) / (PI * s **2 )
(1/surface_ground) * integrate.dblquad(f, 0, surface_ground, 0, surface_collector)
注: A1和A2为2个比表面积，A1为surfa
浏览 3提问于2022-08-03得票数 0
import astropy.io.fits as ap
import numpy as np
import matplotlib
import matplotlib.pyplot as plt
from scipy.interpolate import interp1d
from scipy.ndimage import gaussian_filter
import math
from scipy.special import sph_harm
from scipy.special import legendre
浏览 13提问于2019-11-13得票数 1
2回答
我想将关键字参数传递给dblquad或nquad中的被积函数。这里有没有可能有一个关键字参数，或者我应该选择只有位置参数？
基本上，我试着把字典作为一个普通的参数传递。下面是我这样做的尝试：
def foo(A, B, **kwargs):
alpha = kwargs.get('alpha', 1.0)
beta = kwargs.get('beta', 1.0)
return A*alpha+B*beta
def integrator(**kwargs):
alpha = kwargs.get('alpha',
浏览 11提问于2019-03-25得票数 2
回答已采纳
我试图使用python生成一个双积分图，matplotlib。定义了X和Y域，Z是X和Y的双重积分。
import matplotlib.pyplot as plt
import numpy as np
import scipy
from scipy.integrate import dblquad
fig, ax = plt.subplots(subplot_kw={"projection": "3d"})
def function_I(x, y):
return np.sin(x)*np.cos(y/5)
#Make data
X = np.a
浏览 5提问于2022-03-11得票数 0
回答已采纳
1回答
我有以下的积分(我不知道为什么屏幕截图这么低，所以你需要滚动一点)。我尝试了以下代码：
from scipy import pi, sqrt
from scipy.integrate import dblquad
func = lambda x,y: ((1-sqrt(x**2+y**2))/(sqrt(x**2+y**2)))*(x**2)
x1,x2 = -1, 1
y1,y2 = lambda x: -sqrt(1-x**2), lambda x: sqrt(1-x**2)
print (dblquad(func, x1, x2, y1, y2))
但我一直
The occurrence
浏览 3提问于2020-12-01得票数 1
回答已采纳
import scipy as sp
from scipy import integrate
from sympy.abc import x, i
from scipy.integrate import dblquad
from sympy import *
eq = dblquad(lambda y, t: (5*exp(-y**2) - 5 * exp(-y**2) * exp(-10/(1+x**2))) * exp(i*x*t) * exp(-i*y*t) / (1 - exp(-10/1+t**2)), -oo,oo,-oo, oo)
print(eq)
我已经将所有可能的
浏览 9提问于2022-04-05得票数 -2
我有一个函数，返回一个需要集成的numpy数组。我想使用scipy.integrate.dblquad，但它需要函数才能返回浮点数。
我试着在dblquad上使用dblquad，它不起作用。经过一些思考，我意识到这种方法的问题是，矢量化的dblquad需要函数向量，而不是函数返回向量。
我想用任何方法来解决这个问题，除了拆分数组并逐个案例地集成(已经这样做了，代码就一团糟)。我想到的第一件事是在返回等效函数的数组中转换返回数组的函数，但我不知道这是否可能。
我会放些说明性的代码。下面，x和y是numpy数组，a和b是浮点数。
import numpy as np
from scipy.inte
浏览 1提问于2014-03-14得票数 3
我正在使用以下函数计算两个正态二元分布的重叠
function [ oa ] = bivariate_overlap_integral(mu_x1,mu_y1,mu_x2,mu_y2)
%calculating pdf. Using x as vector because of MATLAB requirements for integration
bpdf_vec1=@(x,y,mu_x,mu_y)(exp(-((x-mu_x).^2)./2.-((y-mu_y)^2)/2)./(2*pi));
%calcualting overlap of two distrib
浏览 0提问于2013-04-16得票数 2
回答已采纳
import numpy as np
import scipy.integrate as integrate
import scipy.special as special
import matplotlib.pyplot as plt
omega = 5
omega *= 2*np.pi
time = np.arange(0.0, 1, 0.0001)
amp = np.sin(omega * time)
在这里，我设定了upa时间，并为所需的正弦波创建了一个omega。在这种情况下，所需的正弦波是sin(5 *2* PI *时间过去了)。
这里我画出正弦波正弦波的正弦波正弦波，得到了
浏览 6提问于2022-05-09得票数 2
1回答
Matlab中数值积分的精度
、、、、
我正在尝试在矩形区域上集成一个解析函数( sqrt和trig函数的组合)。它在该区域没有奇异点，似乎是使用dblquad的完美候选者。我的问题是如何评估Matlab提供给我的数值的准确性。在不知道积分的确切值的情况下，我们如何证明有效数字的合理性？当需要给出一个精度为特定位数的值时，您应该能够对齐。考虑到这个值是用Matlab计算的，有可能实现吗？
浏览 4提问于2013-05-16得票数 4
2回答
我对做二维数值积分感兴趣。现在我正在使用scipy.integrate.dblquad，但是它非常慢。请看下面的代码。我需要用完全不同的参数来评估这个积分的100次。因此，我想使处理尽可能快和有效。守则是：
import numpy as np
from scipy import integrate
from scipy.special import erf
from scipy.special import j0
import time
q = np.linspace(0.03, 1.0, 1000)
start = time.time()
def f(q, z, t):
re
浏览 12提问于2020-04-04得票数 4
回答已采纳
我正在将一篇天体物理学论文转录成python (这篇文章没有定义x和y，所以现在我已经将它们设置为零)。但我在策划错误时遇到了困难：
TypeError:只有size-1数组可以转换为Python标量。
请注意，如果我将plt.plot(R,Vd)放在这里，就会产生错误：
ValueError：x和y必须具有相同的一维，但必须具有(100,)和(1,)形状。
import numpy as np
import matplotlib as mpl
from matplotlib import pyplot as plt
import scipy.integrate as inte
浏览 1提问于2020-02-02得票数 1
2回答
我试图在u和xx上集成我的函数，然后将这些值存储在一个矩阵中，以便用imshow或pcolormesh绘制它们。集成的界限是0 < u < inf和-inf < xx < inf。目前，我的极限只有10，直到我能弄清楚这一点。
import numpy as np
import pylab as pl
from scipy.integrate import dblquad
b = 50.0
x = np.linspace(-10, 10, 1000)
y = np.linspace(0, 10, 1000)
T = pl.zeros([len(x), len(y)
浏览 2提问于2013-10-17得票数 0
回答已采纳
1回答
多处理和枕木(dblquad)
、、、、
我试图加快python中的以下代码：
import numpy as np
import matplotlib.pyplot as plt
from scipy import interpolate
from scipy import integrate
import camb
from tqdm import tqdm
import os
#Reading a PS
dir = os.getcwd()
data = np.loadtxt(dir+"/ps1-peacock.txt")
kh = data[:,0]
p_lin = data[:,1]
p_nlin = dat
浏览 0提问于2020-08-13得票数 0
1回答
我想创建一个从频率f1开始到频率f2结束的正弦波。下面是我使用的代码： import matplotlib.pyplot as plt
import numpy as np
def freq_interp(dur,f1,f2,fs=44100):
num_samples = fs*dur
t = np.linspace(0,dur,num_samples)
a = np.linspace(0,1,num_samples)
f = (1-a)*f1+a*f2 # interpolate
samples = np.cos(2*np.pi*f*t)
浏览 21提问于2019-06-07得票数 0
回答已采纳
2回答
我是社会学的PhD学生，正在写我的论文。在一些数据分析过程中，我遇到了以下问题。
我有一张表，列出了一系列年来一些变量的测量值。值算作“在给定的一年中有多少特定类型的事件”？以下是它的外观示例：
year var
1983 22
1984 55
1985 34
1986 29
1987 15
1988 20
1989 41
因此，例如在一九八四年，全年有55宗这类事件发生。
在[1983,1990]实数域上表示这些数据的一种方法是使用分段函数f：
f(x) = var if floor(x) == year, for all x in [1983
浏览 3提问于2016-01-25得票数 0
回答已采纳
3回答
我很难找到如何从正弦波形成一个三角形波(而不是锯齿波)。
我知道如何为Square wave创建它：
if( sineValue >= 0 )
value = amp;
else
value = -amp;
但我不确定如何改变这一点，以适应三角波。
浏览 2提问于2011-09-13得票数 1
回答已采纳
1回答
我正试图用动画的线性度来解决一个挑战。假设我们有一个拔河技工。一支球队要强大得多，它以每秒1单位的速度从A点(0)移动到B点(10)。计算起来很容易，从0移动到10需要10秒。
现在的问题是，我想模拟两支球队都在互相争斗，试图赢得胜利，所以价值是波动的，非常混乱，但最终的结果是，10秒后，它就在B( 10 )点。
最初的想法是使用区间为<0,1>的随机放松函数，所以没有outBack、OutElastic等，因为我不想让值低于0或超过10。问题是放松函数对平滑的动画很好，但不适合模拟战斗。
我的第二个想法是使用多个随机宽松函数。例如，在<0,1>之间，我使用InQuad
浏览 0提问于2022-12-07得票数 0
回答已采纳
1回答
我正在尝试使用dblquad计算坐标网格上的fresnel integral。但它花了很长时间，最后没有任何结果。
下面是我的代码。在这段代码中，我只集成了一个10x10的网格，但我至少需要集成一个500x500的网格。
import time
st = time.time()
import pylab
import scipy.integrate as inte
import numpy as np
print 'imhere 0'
def sinIntegrand(y,x, X , Y):
a = 0.0001
R = 2e-3
z = 10e-3
浏览 1提问于2015-05-12得票数 1
2回答
我正在试着用scipy计算一个定的二重积分。被积函数有点复杂，因为它包含一些概率分布，以加权表示x和y的每个值的可能性(类似于混合模型)。以下代码的计算结果为负数，但它应该绑定为0,1。此外，它还需要大约半个小时的计算时间。
我有两个问题。
1)有没有更好的方法来计算这个积分？
2)这个负值是从哪里来的？对我来说，最大的问题是如何加快计算速度，因为我可以在我的代码中找到导致负值的错误。
from scipy import stats
from scipy.integrate import dblquad
import itertools
p= [list whose entries are
浏览 1提问于2010-10-28得票数 0
我想对以下内容进行数值积分：
哪里
a、b和β都是常量，为了简单起见，它们都可以设置为1。
无论是使用dblquad的Matlab，还是使用NIntegrate的Mathematica，都无法处理由分母产生的奇点。因为它是二重积分，所以我不能在Mathematica中指定奇点在哪里。
我确信它不是无穷大的，因为这个积分是基于微扰理论的，并且没有
已经被发现了(只是不是我发现的，所以我不知道它是怎么做到的)。
有什么想法吗？
浏览 0提问于2011-08-16得票数 8
回答已采纳
3回答
SciPy有三种对样本进行一维积分的方法(trapz、simps和romb)和一种对函数进行二维积分的方法(dblquad)，但它似乎没有对样本进行二维积分的方法--甚至在矩形网格上也没有。
我看到的最接近的是scipy.interpolate.RectBivariateSpline.integral --您可以从矩形网格上的数据创建一个RectBivariateSpline，然后集成它。然而，这并不是很快。
我想要比矩形方法更准确的方法(也就是把所有的东西加起来)。比方说，我可以使用二维辛普森规则，通过制作一个具有正确权重的数组，将其乘以我想要积分的数组，然后对结果求和。
然而，如果已经有更
浏览 2提问于2013-12-19得票数 23
回答已采纳
1回答
函数作为积分变量
、、
我想使用一个函数作为积分变量，例如，在计算积分时：
e^(cos(x))f(sin(x))dcos(x)
其中f(x)是任意函数。
如何更改集成变量，使我在dcos(x)上进行集成，而不是dx (默认情况下)？
编辑：
在计算积分(如a*x+b )时，可以通过以下方式对其进行评估：
def f(x,a,b):
return a*x+b
print quadrature(f,0,1,args=(2,3))
在这个例子中，积分变量是dx，但是现在我想将积分变量更改为函数，比如x**2 (我知道这可以解析地解决，但是我想应用它，这样一个函数是不可解析的。)，或者cos(x)。
浏览 2提问于2014-10-30得票数 1
回答已采纳
我需要执行一个2D积分(一维有一个无限界)。在MatLab中，我使用了integral2：
int_x = integral2(fun, 0, inf, 0, a, 'abstol', 0, 'reltol', 1e-6);
在Python中，我曾尝试过“‘ve”的dblquad：
int_x = scipy.integrate.dblquad(fun, 0, numpy.inf, lambda x: 0, lambda x: a, epsabs=0, epsrel=1e-6)
也尝试过使用嵌套的单四角体。不幸的是，这两种选择都比MatLab的要长80倍。
我的
浏览 0提问于2018-09-15得票数 1
2回答
有效计算双积分
、、
I使用scipy计算下列积分：
from scipy.stats import norm
def integrand(y, x):
# print "y: %s x: %s" % (y,x)
return (du(y)*measurment_outcome_belief(x, 3)(y))*fv_belief(item.mean, item.var)(x)
return dblquad(
integrand, norm.ppf(0.001, item.mean, item.var),
norm.ppf(0.999, i
浏览 4提问于2015-05-24得票数 2
回答已采纳
1回答
我试图使用matlab集成x，通过使用以下命令，这个任务很简单：
syms x
a=int(x)
问题是我不知道如何实现数值积分。我希望使用不同的技术使用一定数量的间隔来集成x。
有人能帮我处理数值积分的语法要求吗？MathWorks站点并不是很有帮助。
我也知道有一种叫做陷阱的方法，但是我正在寻找matlab中的其他方法，比如Riemann和近似。
更新
具体来说，我要找的是一个函数，它会把x分解成8块区域，然后把这8块加在一起。是否有除trapz之外的预定义函数来执行这样的事情？
好吧，我想我得出的结论是不存在这样的事情。你得自己做。
浏览 4提问于2013-05-20得票数 0
回答已采纳
我看过用PID控制器控制dc到DC buck变换器输出直流电压的例子。
我想实现直流电交转换器。
我试过用PID控制器。
由于我是比较输出交流信号(连续变化)和参考正弦信号，PID控制器没有提供预期的输出。
在DC到DC的情况下，两个参考信号都是恒定的。
在我的例子中，它是连续变化的正弦波。
有人能提出任何方法来获得准确的直流到交流转换吗？
PID控制器能在高频(100 high )连续变化的参考信号下工作吗？
浏览 3提问于2015-06-25得票数 0
1回答
曲线拟合，三角模型
、、、
我试图从表中建模给定x和y值的特定trig函数。我应该找到参数的函数是$A\cdot sin(x\cdot + ph) +d$。$A$是振幅，$x $是周期，$ph$是相移，$d$是垂直位移。这里，$t$是变量。
from pylab import *
from scipy.optimize import curve_fit
time = [1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23]
waterlevel = [125,199,223,176,101,73,109,181,219,188,122,85]
def h(t,x, A , p
浏览 3提问于2022-08-05得票数 0
回答已采纳
2回答
我想使用scipy.integrate中的dblquad重复计算二维复杂积分。由于评估的数量将是相当高的，我想增加我的代码的评估速度。
Dblquad似乎不能处理复杂的被积函数。因此，我将复数被积分为实数部分和虚数部分：
def integrand_real(x, y):
R1=sqrt(x**2 + (y-y0)**2 + z**2)
R2=sqrt(x**2 + y**2 + zxp**2)
return real(exp(1j*k*(R1-R2)) * (-1j*z/lam/R2/R1**2) * (1+1j/k/R1))
def integrand_imag(
浏览 0提问于2013-04-04得票数 11
回答已采纳
我正试图通过处理来画出正弦波。然而，它看起来像是输出的正弦波的一个很差的采样版本。我是否需要用一系列的线条替换形状，还是有其他的解决方案？
我尝试将变量转换为浮点数，并将c++更改为c += 1。
noFill();
stroke(255);
beginShape();
translate(0, 100);
for (int c = 0; c <= width; c += 1)
{
vertex(c, (float) 100 * sin(c / 50));
}
endShape();
我希望它能逐像素地穿越窗口，形成一个平滑的形状.我实际得到的似乎
浏览 0提问于2019-07-03得票数 2
回答已采纳
2回答
我有一个两个变量x和y的联合密度函数，我需要在函数f(X, Y) = y*e**(-y(x+1))中使用python中的quad计算X和Y中的边缘密度函数
from scipy.integrate import dblquad import numpy as np import math
def f(x,y):
return y*math.exp(-y(x+1)) # Joint Density Function
ans,err = dblquad(f,0,math.inf, lambda x: 0 , lambda x:math.inf)
ans
我在Jupyter noteb
浏览 0提问于2019-10-12得票数 0
2回答
我要根据三维离散的非均匀密度分布计算球体的质量。假设一组不同密度的3x3x3立方体是由一个球体所刻下的。使用Python对分区块进行总结的最快方法是什么？
我试着用scipy.integrate.dblquad计算一个球的数学方程下的体积: x^2 +y^2 +z^2 = R^2。但是，只有在边界小于球体半径的情况下，结果才是有效的，而对于50,000个每个含有27个立方体的球体，重复计算将非常慢。
另一方面，由于质量分布比较粗糙和离散，CoM计算的一般方程不能用在我看来。
浏览 0提问于2010-12-10得票数 2
4回答
标题是这样说的，这似乎很简单，但我遇到了困难。
我想创造一个正弦波，在我输入的功率下，它在频率上呈指数衰减。我尝试了很多不同的东西，我要么变得胡言乱语，要么我相信是一种声音的拍子？
胡言乱语
在这个例子中，数组乘法器从1到>4。当我绘制cos(70000000*t)和cos(4*70000000*t)时，这两个图看起来都很好，但是下面的代码中的图看起来就像噪音。
t = 1:.0000000004:1.0000004;
multiplier = linspace(1,2,1001).^-2;
reference_signal = cos(700000000*t.*multiplier);
浏览 5提问于2014-04-13得票数 0
回答已采纳
1回答
方波频率表示法
、、、
这是信号的新手，我对方波的频率表示感到困惑。如果我错了，请纠正我，周期方波是由奇次谐波正弦波组成的，它是基频的倍数，在频域中，这只表示基频的每个奇数倍的递减脉冲。正弦函数与方波有什么关系？方波的正确频谱表示是什么？谢谢!
浏览 3提问于2016-02-28得票数 0
2回答
我现在使用的是Python，我有一个函数需要对不同的常量进行乘法运算。
例如，如果我有f(x,y)= x^2y+a，其中'a‘是某个常量(可能是常量列表)。
如果'a‘是一个列表(大小未知，因为它取决于输入)，那么如果我们说a = [1,3,7]，我想做的操作是
(x^2y+1)*(x^2y+3)*(x^2y+7)
而是泛化为'a‘中的n个元素。有没有一个简单的方法可以在Python中做到这一点，因为我想不出解决这个问题的好办法？如果'a‘中的大小是固定的，那么它看起来会简单得多，因为我可以单独定义函数，然后将它们相乘在一个新的函数中，但由于大小不是固定的，这种
浏览 0提问于2019-12-16得票数 0
这个脚本计算(x0，y0，z0)的电场来自均匀带电的环空，使用SciPy的。
我的问题是如何使用gfun和hfun来定义内部积分极限对外部积分变量的函数依赖性。如果我选择在笛卡尔坐标中进行积分，这将是一个曲奇刀的形状，但是当我使用圆柱坐标时，函数返回常数浮动。
是否有一种方法可以消除这些简单返回常量的函数调用，以避免函数调用的时间损失？
该示例可能没有以其他方式进行优化，但它只是一个简单的示例，可以演示gfun和hfun的使用。
def Excalc(r, th):
x, y, z = r*np.cos(th), r*np.sin(th), 0.0
return (x0-x)
浏览 0提问于2018-03-23得票数 2
回答已采纳
我正在使用VS代码中Python的朱庇特笔记本扩展和Sympy模块做作业。我正在尝试用Latex字体呈现方程式，但它并不像我想象的那样有效。有谁能猜出问题会是什么吗？我已经添加了与问题相关的捕获图像。
我在用
v2020.12.414227025Python v2020.11.371526539扩展VS代码1.52
谢谢你，祝你星期天愉快！
当我键入像下面这样的随机方程时，它会以这种方式显示。
然而，当我输入我在网上找到的某些等式时，它就会以适当的Latex字体显示。
浏览 3提问于2020-12-13得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
