Matlab - 在任意點取一個矩陣的梯度

Question

我想在一個二維矩陣（x，y）上實現龍格庫塔的四階積分方法。矩陣包含每個點的高度值（浮點數）。這個想法是在矩陣中放置一個粒子，並觀察它的路徑，因爲它被方向向量u和v'推動'。

我已經實現了歐拉的積分方法，它已經使用了簡單的matlab中構建的'gradient'函數。但是，對於RK4積分，我需要在可變點上取4個斜率，而不是斜率給出的預定斜率。

我認爲Matlab漸變函數本身沒有這些類型的參數（從我從Matlab文檔中理解的）。我應該採取什麼樣的方法？

如果我要編寫自己的漸變函數來計算任意點處的斜率，我不太確定如何編寫它，因爲它不是直向函數。這種方法有什麼建議嗎？

謝謝。

Answer 1

如果我寫我自己的梯度 功能，在 任意點計算斜率，我不太肯定 如何寫它，因爲它不是 直盼着功能。有沒有 這種方法的建議？

我懷疑interp2可能是你要找的。你可以給它包含離散點上的梯度的矩陣，並且它會內插以在任何中間點找到梯度的值。