與Python/Numpy/Scipy的三重集成

Question

我是Python新手，學習遵循Python「科學講義 Rlease 2013.1」教程。 請幫我解決下面srcreenshot中的這個三重一體化問題（第70頁）。我已經介紹了該教程的以前的內容。 如果可能的話，請提供步驟明確的命令和解釋，因爲作爲航空航天工程師編程概念對我而言是新的。

謝謝。

練習：粗積分近似值 寫函數f（A，B，C）返回一個^ B - ℃。在參數 範圍[0,1] x [0,1] x [0,1]中形成一個包含其值的24x12x6陣列。

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的值，將所有這些值相加，然後除以該數組中的元素數。