本文简要介绍 python 语言中 scipy.linalg.expm_frechet
的用法。
用法:
scipy.linalg.expm_frechet(A, E, method=None, compute_expm=True, check_finite=True)#
A 的矩阵 index 在 E 方向上的 Frechet 导数。
- A: (N, N) 数组
其中矩阵取矩阵 index 。
- E: (N, N) 数组
采用 Frechet 导数的矩阵方向。
- method: str,可选
算法的选择。应该是其中之一
SPS(默认)
blockEnlarge
- compute_expm: 布尔型,可选
除了 expm_frechet_AE 之外,是否还计算 expm_A。默认为真。
- check_finite: 布尔型,可选
是否检查输入矩阵是否仅包含有限数。禁用可能会提高性能,但如果输入确实包含无穷大或 NaN,则可能会导致问题(崩溃、非终止)。
- expm_A: ndarray
A的矩阵 index 。
- expm_frechet_AE: ndarray
A 的矩阵 index 在 E 方向上的 Frechet 导数。
- 为了
compute_expm = False
, 只要expm_frechet_AE被退回。
参数 ::
返回 ::
注意:
本节介绍可以通过方法参数选择的可用实现。默认方法是 SPS。
方法 blockEnlarge 是一种简单的算法。
方法SPS是Scaling-Pade-Squaring[1].这是一个复杂的实现,只需要大约 3/8 的时间作为天真的实现。渐近线是相同的。
参考:
[1]Awad H. Al-Mohy 和 Nicholas J. Higham (2009) 计算矩阵 index 的 Frechet 导数,并应用于条件数估计。 SIAM 矩阵分析与应用杂志,30 (4)。第 1639-1657 页。 ISSN 1095-7162
例子:
>>> import numpy as np >>> from scipy import linalg >>> rng = np.random.default_rng()
>>> A = rng.standard_normal((3, 3)) >>> E = rng.standard_normal((3, 3)) >>> expm_A, expm_frechet_AE = linalg.expm_frechet(A, E) >>> expm_A.shape, expm_frechet_AE.shape ((3, 3), (3, 3))
创建一个包含 [[A, E], [0, A]] 的 6x6 矩阵:
>>> M = np.zeros((6, 6)) >>> M[:3, :3] = A >>> M[:3, 3:] = E >>> M[3:, 3:] = A
>>> expm_M = linalg.expm(M) >>> np.allclose(expm_A, expm_M[:3, :3]) True >>> np.allclose(expm_frechet_AE, expm_M[:3, 3:]) True
相关用法
- Python SciPy linalg.expm_multiply用法及代码示例
- Python SciPy linalg.expm_cond用法及代码示例
- Python SciPy linalg.expm用法及代码示例
- Python SciPy linalg.eigvalsh_tridiagonal用法及代码示例
- Python SciPy linalg.eig_banded用法及代码示例
- Python SciPy linalg.eigvalsh用法及代码示例
- Python SciPy linalg.eigh_tridiagonal用法及代码示例
- Python SciPy linalg.eig用法及代码示例
- Python SciPy linalg.eigh用法及代码示例
- Python SciPy linalg.eigs用法及代码示例
- Python SciPy linalg.eigvals用法及代码示例
- Python SciPy linalg.eigvals_banded用法及代码示例
- Python SciPy linalg.eigsh用法及代码示例
- Python SciPy linalg.cdf2rdf用法及代码示例
- Python SciPy linalg.LaplacianNd用法及代码示例
- Python SciPy linalg.solve_circulant用法及代码示例
- Python SciPy linalg.polar用法及代码示例
- Python SciPy linalg.clarkson_woodruff_transform用法及代码示例
- Python SciPy linalg.rsf2csf用法及代码示例
- Python SciPy linalg.hessenberg用法及代码示例
- Python SciPy linalg.tril用法及代码示例
- Python SciPy linalg.triu用法及代码示例
- Python SciPy linalg.svd用法及代码示例
- Python SciPy linalg.ishermitian用法及代码示例
- Python SciPy linalg.invhilbert用法及代码示例
注:本文由纯净天空筛选整理自scipy.org大神的英文原创作品 scipy.linalg.expm_frechet。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。