线性代数¶
矩阵分解¶
矩阵分解 是将一个矩阵分解为数个矩阵的乘积,是线性代数中的一个核心概念。
下面的表格总结了在 Julia 中实现的几种矩阵分解方式。具体的函数可以参考标准库文档的 Linear Algebra 章节。
Cholesky |
Cholesky 分解 |
CholeskyPivoted |
主元 Cholesky 分解 |
LU |
LU 分解 |
LUTridiagonal |
LU factorization for Tridiagonal matrices |
UmfpackLU |
LU factorization for sparse matrices (computed by UMFPack) |
QR |
QR factorization |
QRCompactWY |
Compact WY form of the QR factorization |
QRPivoted |
主元 QR 分解 |
Hessenberg |
Hessenberg 分解 |
Eigen |
特征分解 |
SVD |
奇异值分解 |
GeneralizedSVD |
广义奇异值分解 |
特殊矩阵¶
线性代数中经常碰到带有对称性结构的特殊矩阵,这些矩阵经常和矩阵分解联系到一起。 Julia 内置了非常丰富的特殊矩阵类型,可以快速地对特殊矩阵进行特定的操作.
下面的表格总结了 Julia 中特殊的矩阵类型,其中也包含了 LAPACK 中的一些已经优化过的运算。
Hermitian |
埃尔米特矩阵 |
Triangular |
上/下 三角矩阵 |
Tridiagonal |
三对角矩阵 |
SymTridiagonal |
对称三对角矩 |
Bidiagonal |
上/下 双对角矩阵 |
Diagonal |
对角矩阵 |
UniformScaling |
缩放矩阵 |
基本运算¶
| 矩阵类型 | + |
- |
* |
\ |
其它已优化的函数 |
|---|---|---|---|---|---|
Hermitian |
XY | inv, sqrtm, expm |
|||
Triangular |
XY | XY | inv, det |
||
SymTridiagonal |
X | X | XZ | XY | eigmax/min |
Tridiagonal |
X | X | XZ | XY | |
Bidiagonal |
X | X | XZ | XY | |
Diagnoal |
X | X | XY | XY | inv, det, logdet, / |
UniformScaling |
X | X | XYZ | XYZ | / |
图例:
| X | 已对矩阵-矩阵运算优化 |
| Y | 已对矩阵-向量运算优化 |
| Z | 已对矩阵-标量运算优化 |
矩阵分解¶
| 矩阵类型 | LAPACK | eig |
eigvals |
eigvecs |
svd |
svdvals |
|---|---|---|---|---|---|---|
Hermitian |
HE | ABC | ||||
Triangular |
TR | |||||
SymTridiagonal |
ST | A | ABC | AD | ||
Tridiagonal |
GT | |||||
Bidiagonal |
BD | A | A | |||
Diagonal |
DI | A |
图例:
| A | 已对寻找特征值和/或特征向量优化 | 例如 eigvals(M) |
| B | 已对寻找 ilth 到 ihth 特征值优化 |
eigvals(M, il, ih) |
| C | 已对寻找在 [vl, vh] 之间的特征值优化 |
eigvals(M, vl, vh) |
| D | 已对寻找特征值 x=[x1, x2,...] 所对应的特征向量优化 |
eigvecs(M, x) |
缩放运算¶
A UniformScaling operator represents a scalar times the identity operator, λ*I. The identity operator I is defined as a constant and is an instance of UniformScaling. The size of these operators are generic and match the other matrix in the binary operations +,``-,``* and \. For A+I and A-I this means that A must be square. Multiplication with the identity operator I is a noop (except for checking that the scaling factor is one) and therefore almost without overhead.