本文简要介绍 python 语言中 scipy.linalg.solve_sylvester
的用法。
用法:
scipy.linalg.solve_sylvester(a, b, q)#
计算 Sylvester 方程 的解 (X)。
- a: (M, M) 数组
Sylvester 方程的前导矩阵
- b: (N, N) 数组
Sylvester 方程的尾随矩阵
- q: (M, N) 数组
右侧
- x: (M, N) ndarray
西尔维斯特方程的解。
- LinAlgError
如果没有找到解决方案
参数 ::
返回 ::
抛出 ::
注意:
通过 Bartels-Stewart 算法计算 Sylvester 矩阵方程的解。 A 和 B 矩阵首先进行 Schur 分解。生成的矩阵用于构造替代 Sylvester 方程 (
RY + YS^T = F
),其中 R 和 S 矩阵采用 quasi-triangular 形式(或者,当 R、S 或 F 为复数时,为三角形形式)。然后使用 LAPACK 中的*TRSYL
直接求解简化方程。例子:
给定 a、b 和 q 求解 x:
>>> import numpy as np >>> from scipy import linalg >>> a = np.array([[-3, -2, 0], [-1, -1, 3], [3, -5, -1]]) >>> b = np.array([[1]]) >>> q = np.array([[1],[2],[3]]) >>> x = linalg.solve_sylvester(a, b, q) >>> x array([[ 0.0625], [-0.5625], [ 0.6875]]) >>> np.allclose(a.dot(x) + x.dot(b), q) True
相关用法
- Python SciPy linalg.solve_circulant用法及代码示例
- Python SciPy linalg.solve_banded用法及代码示例
- Python SciPy linalg.solve_discrete_lyapunov用法及代码示例
- Python SciPy linalg.solve_toeplitz用法及代码示例
- Python SciPy linalg.solve_continuous_lyapunov用法及代码示例
- Python SciPy linalg.solve_continuous_are用法及代码示例
- Python SciPy linalg.solve_discrete_are用法及代码示例
- Python SciPy linalg.solve_triangular用法及代码示例
- Python SciPy linalg.solve用法及代码示例
- Python SciPy linalg.solveh_banded用法及代码示例
- Python SciPy linalg.svd用法及代码示例
- Python SciPy linalg.spsolve用法及代码示例
- Python SciPy linalg.spsolve_triangular用法及代码示例
- Python SciPy linalg.splu用法及代码示例
- Python SciPy linalg.spilu用法及代码示例
- Python SciPy linalg.svdvals用法及代码示例
- Python SciPy linalg.sqrtm用法及代码示例
- Python SciPy linalg.svds用法及代码示例
- Python SciPy linalg.sinm用法及代码示例
- Python SciPy linalg.schur用法及代码示例
- Python SciPy linalg.sinhm用法及代码示例
- Python SciPy linalg.signm用法及代码示例
- Python SciPy linalg.subspace_angles用法及代码示例
- Python SciPy linalg.eigvalsh_tridiagonal用法及代码示例
- Python SciPy linalg.cdf2rdf用法及代码示例
注:本文由纯净天空筛选整理自scipy.org大神的英文原创作品 scipy.linalg.solve_sylvester。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。