如何使用 Python SciPy 求解三角矩陣方程?

線性函式 scipy.linalg.solveh_triangular 用於求解帶狀矩陣方程。在下面給出的示例中，我們將求解三角系統 ax = b 其中 −

$$\mathrm{a} = \begin{bmatrix} 3 & 0 & 0 & 0\ 2 & 1 & 0 & 0\ 1 &0 &1 &0 \ 1& 1& 1& 1 \end{bmatrix};\; \mathrm{b} =\begin{bmatrix} 1\ 2\ 1\ 2 \end{bmatrix}$$

示例

from scipy.linalg import solve_triangular
import numpy as np
a = np.array([[3, 0, 0, 0], [2, 1, 0, 0], [1, 0, 1, 0], [1, 1, 1, 1]])
b = np.array([1, 2, 1, 2])
x = solve_triangular(a, b, lower=True)
print (x)

輸出

array([ 0.33333333, 1.33333333, 0.66666667, -0.33333333])

Gaurav Kumar

更新日期: 24-11-2021

376 次觀看

開啟你的 職業生涯

完成課程獲得認證

開始學習