如何找到Python中矩陣指數的積分

我有一個形式的矩陣，比如e^(Ax)，其中A是一個方形矩陣。我怎樣才能將它從一個給定的值a整合到另一個值b，以便輸出是一個相應的數組？如何找到Python中矩陣指數的積分

2013-02-15 user2076206

在這種情況下，「a」和「b」是什麼？ – tacaswell 2013-02-15 17:36:14

我懷疑你想使用'scipy' http://docs.scipy.org/doc/scipy/reference/integrate.html – tacaswell 2013-02-15 17:36:54

如果A具有正確的屬性，可以通過計算其特徵向量和特徵值將其轉換爲對角線形式A0。在對角線形式中，解決方案是sol = [exp(A0*b) - exp(A0*a)] * inv(A0)，其中A0是具有特徵值的對角矩陣，並且inv(A0)僅包含其對角線中的特徵值的倒數。最後，通過將解與左邊的特徵值和右邊的特徵值的轉置相乘來轉換回解決方案：transpose(eigvecs) * sol * eigvecs。

來源

2013-02-15 17:38:15

非常感謝Aradi對於問題的好方法和tcaswell的參考。但是，那麼，我怎麼能在Python中實現這個例如scipy.integrate.quad或類似的東西呢？我知道如何做到這一點，而不是一個矩陣。形成一個循環是唯一的解決方案？再次感謝。哦，a和b是標量，x是整合的限制。 – user2076206 2013-02-15 18:06:36

如何找到Python中矩陣指數的積分

回答

相關問題