如何使用 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}$$
示例
fromscipy.linalgimport 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])
热门推荐
10 八一幼儿祝福语大全简短
11 公司乔迁食堂祝福语简短
12 婚礼结束聚餐祝福语简短
13 儿媳买车妈妈祝福语简短
14 毕业送礼老师祝福语简短
15 同事辞职正常祝福语简短
16 恭贺新婚文案祝福语简短
17 金店立秋祝福语简短英文
18 婆婆高寿祝福语大全简短