如何在 Python 中计算这个数值积分?

如何解决如何在 Python 中计算这个数值积分?

我有一个向量列表:

x = [x1,x2,x3,x4]

其中 x1,.. 是一些变量。 我有一个矩阵

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 举报,一经查实,本站将立刻删除。

相关推荐


Selenium Web驱动程序和Java。元素在(x,y)点处不可单击。其他元素将获得点击?
Python-如何使用点“。” 访问字典成员?
Java 字符串是不可变的。到底是什么意思?
Java中的“ final”关键字如何工作?(我仍然可以修改对象。)
“loop:”在Java代码中。这是什么,为什么要编译?
java.lang.ClassNotFoundException:sun.jdbc.odbc.JdbcOdbcDriver发生异常。为什么?
这是用Java进行XML解析的最佳库。
Java的PriorityQueue的内置迭代器不会以任何特定顺序遍历数据结构。为什么?
如何在Java中聆听按键时移动图像。
Java“Program to an interface”。这是什么意思?
Java在半透明框架/面板/组件上重新绘画。
Java“ Class.forName()”和“ Class.forName()。newInstance()”之间有什么区别?
在此环境中不提供编译器。也许是在JRE而不是JDK上运行?
Java用相同的方法在一个类中实现两个接口。哪种接口方法被覆盖?
Java 什么是Runtime.getRuntime()。totalMemory()和freeMemory()?
java.library.path中的java.lang.UnsatisfiedLinkError否*****。dll
JavaFX“位置是必需的。” 即使在同一包装中
Java 导入两个具有相同名称的类。怎么处理?
Java 是否应该在HttpServletResponse.getOutputStream()/。getWriter()上调用.close()?
Java RegEx元字符(。)和普通点?