如何使用 Python SciPy 解決厄米特正帶狀矩陣方程？

可以使用名為 scipy.linalg.solveh_banded 的線性函式來求解帶狀矩陣方程。在下面給出的示例中，我們將求解帶狀系統 Hx = b，其中 −

$$\mathrm{H} = \begin{bmatrix} 8 & 2-1j&0 &0 \ 2+1j & 5& 1j & -2-1j0\ 0 & -1j& 9& \ 0 & 0& -2+1j& 6 \end{bmatrix} \mathrm{b}=\begin{bmatrix} 1\ 1+1j\ 1-2j\ 0 \end{bmatrix}$$

對於下面的示例，我們將上對角線放入陣列 hb −

示例

from scipy.linalg import solveh_banded
hb = np.array([[0, 2-1j, 1j, -2-1j],  [8, 5, 9, 6 ]])
b = np.array([1, 1+1j, 1-2j, 0])
x = solveh_banded(hb, b) 
print(x)

輸出

array([ 0.07318536-0.02939412j, 0.11877624+0.17696461j, 0.10077984-0.23035393j, -0.00479904-0.09358128j])

Gaurav Kumar

更新時間： 2021 年 11 月 25 日

165 次瀏覽

開啟你的 職業生涯

完成課程認證

開始