如何解决如何在 Python 中计算这个数值积分?
x = [x1,x2,x3,x4]
b = array([[0,2],[0,3],[1,[2,3]])
最后我想计算如下:
Integral[ np.exp(-np.dot(np.ones(len(b)),np.cos(x[b[:,0]]-x[b[:,1]]))),{x1,Pi},{x2,Pi}...]
所以看起来当x
的大小很大时,最好象征性地定义x
from sympy import *
x = symarray('x',4)
并且由于x
包含将被积分的变量,有没有什么好的方法可以通过符号变量设置积分并在最后进行评估?
更新:
x =symarray('x',M0)
ij_args = x[b[:,1]]
ij_cosargs = [cos(ij_args[i]) for i in range(len(ij_args))]
integrand = exp(np.sum(ij_cosargs))
integrate(integrand,(x[0],np.pi),(x[1],(x[2],(x[3],np.pi))
上述代码的另一个问题是我必须分别为每个符号输入边界,例如 (x[0],np.pi)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。