我在做一个有4个变量的多重积分,其中2个是函数的极限。但是错误出现在我的一个常量极限变量上。真的不明白为什么。非常感谢你的建议!在
from numpy import sqrt, sin, cos, pi, arcsin, maximum
from sympy.functions.special.delta_functions import Heaviside
from scipy.integrate import nquad
def bmax(x):
return 1.14*10**9/sin(x/2)**(9/7)
def thetal(x,y,z):
return arcsin(3.7*10**15*sqrt(cos(x/2)**2/10**6-1.23*10**10/z+0.003*sin(x/2)**2*(2.51*10**63/sin(x/2)**9/y**7-1))/(z*sin(x/2)**2*cos(x/2)*(2.51*10**63/sin(x/2)**9/y**7-1)))
def rt(x,y):
return 3.69*10**12/(2.5*10**63/sin(x/2)**7*y**7-sin(x/2)**2)
def rd(x,y):
return maximum(1.23*10**10,rt(x,y))
def rl(x,y):
return rd(x,y)*(sqrt(1+5.04*10**16/(rd(x,y)*cos(x/2)**2))-1)/2
def wbound():
return [1.23*10**10,3.1*10**16]
def zbound():
return [10**(-10),pi-10**(-10)]
def ybound(z):
return [0,bmax(z)-10**(-10)]
def xbound(z,y,w):
return [thetal(z,y,w),pi-thetal(z,y,w)]
def f(x,y,z,w):
return [5.77/10**30*sin(z)*sin(z/2)*y*sin(x)*Heaviside(w-rl(z,y))*Heaviside(w-rd(z,y))/w**2]
result = nquad(f, [xbound, ybound,zbound,wbound])
这个错误的原因是,尽管您不希望这些边界依赖于变量,
nquad
仍然将变量传递给您提供给它的函数。因此,绑定函数必须采用正确数量的变量:现在函数
zbound
和ybound
接受额外的变量,但忽略它们。在我不确定最后一个边界}翻转吗?根据^{} 的定义,假定正确的顺序是
^{pr2}$xbound(...)
:你想让变量y
和{编辑:正如kazemakase所指出的,函数
f
应该返回一个float
,而不是一个列表,因此应该删除return语句中的方括号[...]
。在nquad
需要一个bounds
序列作为其第二个参数,语法相当严格。在如果被积函数},其中每个边界可以是2元组,例如
f
依赖于x, y, z, w
,并且这是定义的顺序,bounds
中的术语必须依次为xb
、yb
、zb
和{xb = (xmin, xmax)
或返回2元组的函数。在关键是,这些函数的参数。。。当我们在}来计算
dx
中执行内部积分时,我们有y
、z
和{x
中的边界,因此它必须是def xb(y,z,w): return(..., ...)
-同样def yb(z,w): return (..., ...)
和def zb(w): return (..., ...)
。关于积分的最后一个变量的界限必须是常数。在
总结
相关问题 更多 >
编程相关推荐