python - 与 Python / Numpy / Scipy 三重集成

Question

我是 python 新手，并通过遵循 Python“科学讲义 2013.1 版”教程进行学习。请帮我解决下面屏幕截图中的这个三重积分问题（第 70 页）。我已经介绍了该教程的先前内容。如果可能，请提供带有解释的逐步命令，因为作为航空航天工程师的编程概念对我来说是新的。

谢谢你。

练习：粗略的积分近似 编写一个函数 f(a, b, c)，它返回 a^b - c。形成一个 24x12x6 数组，其中包含参数范围 [0,1] x [0,1] x [0,1] 中的值。

Answer 1

这可能看起来令人生畏，但这个问题会一步一步地告诉你该怎么做。

编写一个f(a, b, c)返回的函数a^b-c。

def f(a, b, c):
    return a ** b - c

形成一个包含其参数范围的 24x12x6 数组[0,1] x [0,1] x [0,1]。给你ogrid提示。所以阅读文档我猜它看起来像：

x = np.ogrid[0:1:24j, 0:1:12j, 0:1:6j]

然后你可以做

f(x[0], x[1], x[2])

然后取平均值

np.mean(f(x[0], x[1], x[2]))

这给了我0.18884234602967925

Answer 2

你误解了这个问题。他们没有要求您计算积分。他们要求您计算该函数在一组点上的平均值。练习的动机是该平均值将是该积分值的近似值。

您需要做的就是形成他们请求的数组，包含 24x12x6 的值，将所有这些值相加，然后除以该数组中的元素数。